<div dir="ltr">Dear all,<div><br></div><div>In recent calculations I have been attempting to calculate uniaxial magnetocrystalline anisotropies in transition-metal magnetic compounds using DFT+U method. </div><div><br></div><div>As I understand it, these calculations require (a) non-collinear spins `noncolin = .true.` and (b) DFT+U with both U and J parameter `lda_plus_u_kind = 1`, as the simpler  `lda_plus_u_kind = 0` is incompatible with non-collinear spin calculations. </div><div>(This is occasionally inconvenient, as DFT+U+J itself is not supported for structural relaxations.)</div><div><br></div><div>1. Are there any known issues (or theoretical concerns) with using DFT+U+J in calculations of magnetic anisotropy in QE?</div><div><br></div><div>2. For reference, I calculate E_MCA = E_IP – E_OP, i.e. the difference in (SCF) energy between configurations with in-plane and out-of-plane magnetization respectively, rather than the perhaps easier force theorem approach. While I don't expect this to change much, does this affect the considerations from (1), if any?</div><div><br></div><div>Any information (or corrections) regarding the above points would be greatly appreciated.</div><div><br></div><div>Peter Stoeckl</div><div>PhD Candidate, Department of Physics</div><div>University of Minnesota</div></div>