My Community

Hi

As I posted the other day, I am running noncollinear calculations and I'm using the tetrahedron method (ISMEAR = -5). BUT, I'm experiencing crashes almost straight away with several occurrences of the message

magnetization density of overlapping atoms calculated
LAPACK: Routine ZPOTRF failed! 1

I've now tried the suggestion from a previous post of using the package lapack_double.o but I get the same error.

Any suggestions?

Thanks

Ben

PS recall my INCAR file is:

SYSTEM = Co monatomic chain (Noncollinear)
LSORBIT = .TRUE.
ISYM = 0
GGA_COMPAT = .FALSE.
ISMEAR = -5
MAGMOM = 4.0 0.05 -0.05
ICHARG = 1
PREC = High
NELM = 300
NELMIN = 4
NBANDS = 20
EDIFF = 1.0e-5
EDIFFG = 1.0e-4
ISPIN = 2
IBRION = 1
ISIF = 2
NSW = 100
LORBIT = 10
RWIGS = 1.302

You are NOT running noncolinear calculations with this INCAR file.
If you want to do so, you have to set LNONCOLLINEAR = .TRUE. -tag ...
read for example: http://cms.mpi.univie.ac.at/vasp/vasp/node142.html

I disagree. Check out this:

http://cms.mpi.univie.ac.at/vasp/vasp/node144.html

in particular the very first line which reads:

LSORBIT = .TRUE. switches on spin-orbit coupling and automatically sets LNONCOLLINEAR= .TRUE..

I rest my case

Ben

Sorry, I overlooked the LSORBIT tag ....

please note that
1) the appropriate k-mesh for a chain of atoms (not interacting with neighboring chains) probably is conflicting with ISMEAR=-5
2) depending on the geometry input / non-collinear setup: with your definition of MAGMOM, the magnetic moment for exactly one atom is specified. If there are more than 1 atoms in the chain, this may lead to a starting density which is not reasonable and causes the ZPOTRF error.