dynamical matrix problem

chwang45 · Post by **chwang45** » Fri Jul 29, 2005 6:00 am

When I generate a dynamical matrix, I have a bunch of 0's after several lines to the end. Would anyone give me the possible reasons? Is the end part of the matrix file related to the displaced ions in the end part of POSCAR?

Post by **graeme** » Sun Jul 31, 2005 5:05 pm

I think the issue here is probably that there are fewer force calls in the OUTCAR file than degrees of freedom in the DISPLACECAR file. The NSW variable should equal to the number of degrees of freedom plus one. If vasp terminates early, or if this was miscounted, there will be fewer force calls in the OUTCAR than are needed to general the full dynamical matrix. In this case, the final lines of the matrix will contain only zeros.

andri · Post by **andri** » Mon Aug 01, 2005 3:02 am

It can also be that the NSW can be higher than the number of degrees of freedom. When all lines in the DISPLACECAR have been exhausted then the system doesn't move and the forces stay the same after that. So when the matrix is set up with finite difference of the force, zero values are obtained.