<div dir="ltr">Dear all,<div><br></div><div>I am trying to calculate the magnetic anisotropy energy of a molecule (VOPc). But I am facing a weird behavior that makes the total magnetization huge and strongly oscillating; I already tried to use "lambda" with "atomic direction" to align the magnetic moment but nothing seems to help. Usually, within the first 2 iterations the magnetic moment becomes huge</div><div><br></div><div>     total magnetization       =     1.41    17.01   150.39 Bohr mag/cell<br>     absolute magnetization    =   201.24 Bohr mag/cell<br>     lambda                    =      10.00 Ry<br></div><div><br></div><div>and then refuses to reduce (the system should be spin-1/2). I already tried to reduce the mixing  mixing_beta=0.01, tried variations of smearing, cutoffs - nothing seems to improve this behavior. What I cannot understand is the huge magnetization after setting it initially to </div><div><br></div><div>starting_magnetization(1)=0.1<br>starting_magnetization(2)=0.0<br>starting_magnetization(3)=0.0<br>starting_magnetization(4)=0.0<br>starting_magnetization(5)=0.0<br>angle1(1)=90<br>angle2(1)=0.0<br>angle1(2)=0.0<br>angle2(2)=0.0<br>angle1(3)=0.0<br>angle2(3)=0.0<br>angle1(4)=0.0<br>angle2(4)=0.0<br>angle1(5)=0.0<br>angle2(5)=0.0<br>noncolin=.true.<br>lspinorb=.true.<br></div><div><br></div><div>with atom types</div><div>ATOMIC_SPECIES<br>V     50.942    V.rel-pbe-spnl-kjpaw_psl.1.0.0.UPF<br>C     12.0107   C.pbe-n-kjpaw_psl.1.0.0.UPF<br>N     14.0067   N.pbe-n-radius_5.UPF<br>O     15.995    O.pbe-n-kjpaw_psl.1.0.0.UPF<br>H      1.000    H.pbe-kjpaw_psl.1.0.0.UPF<br></div><div><br></div><div>Any hint?</div><div><br></div><div>Thanks in advance for your help!</div><div>Chris </div><div><br></div><div><br></div><div><br></div><div><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div>