<div dir="ltr"><div>Dear Peter,</div><div><br></div><div>I don't know how to answer your question 1). As for 2) maybe this paper can help you:</div><div><br></div><div>PHYSICAL REVIEW B 96, 165132 (2017); DOI: 10.1103/PhysRevB.96.165132</div><div dir="ltr"><br></div><div>Best regards,</div><div><br></div><div>Matteo</div><div><br></div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno ven 18 feb 2022 alle ore 03:30 Peter Stoeckl via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Matteo Cococcioni<br></div><div>Department of Physics<br></div>University of Pavia<br>Via Bassi 6, I-27100 Pavia, Italy<br>tel +39-0382-987485<br><div>e-mail <a href="mailto:matteo.cococcioni@unipv.it" target="_blank">matteo.cococcioni@unipv.it</a></div></div></div></div></div></div></div></div>