<div dir="ltr">Hello All!<div><br></div><div>I'm trying to calculate a system with a <b>known</b> AFM ground state. In the attached example, I provide the input data I'm using for CuMnAs -- a tetragonal anti-ferromagnet.</div><div>I've started on this project by first executing the examples, and particularly FeO. </div><div>I tested this example against all sorts of variations: different functionals, with/without U, and so on. Without any further tweaking, the system always converged to the AFM ground-state, provided the initial moments were also oriented in the AFM configuration.</div><div><br></div><div>Which makes my failure in this (CuMnAs) system even more puzzling. Firstly, <b>without</b> any constraints, the system does not converge to an AFM state. </div><div>1. Using 'constrained_magnetization=total' leads to completely wrong results, with a divergent "Magnetic field".</div><div>2. A more-or-less sensible result can be obtained with 'constrained_magnetization='atomic' (as shown), however, the resultant magnetization is not altogether anti-ferromagnetic. Note that the system is in general endowed with PT-symmetry. The resultant eigenvalues <b>DO NOT</b> show this and you can also see the disparity in the magnetic moments of the Mn atoms, as well as eigenvalue difference of more than 1meV for some bands and k-points.</div><div>3. This behavior is <i>weakly</i> dependent on lambda. I tried fiddling around with the values. A certain increase worsens the results, then seems to improve it, only to worsen again. What is a reason value for the constraint, in your estimation? I take it to be 5% of the unperturbed energy (i.e., energy without the constraint).</div><div>4. Testing the <b>exact</b> same system on different software (VASP, in this case), converged very well to the AFM state (i.e., PT symmetry was recovered to less than 1meV).</div><div><br></div><div>What am I doing wrong, therefore? </div><div>I would appreciate any advice.</div><div>Yours thankful,</div><div>Daniel Kaplan</div><div>Dept. of Condensed Matter Physics</div><div>Weizmann Institute of Science</div></div>