<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div>Oops, not new to it at this point just misread the docs ðŸ˜¬ unfortunately SOC is necessary for what im trying to do.</div>
<br>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nicola Marzari <nicola.marzari@epfl.ch><br>
<b>Sent:</b> Thursday, February 13, 2025 7:44 PM<br>
<b>To:</b> Quantum ESPRESSO users Forum <users@lists.quantum-espresso.org>; Johnson, Miles R. <mjohnso7@caltech.edu><br>
<b>Subject:</b> Re: [QE-users] Absolute Magnetization diverging</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText"><br>
<br>
Dear Miles,<br>
<br>
<br>
a few comments:<br>
<br>
1) there are no (or not yet, to be precise) fully relativisitc <br>
pseudopotentials in the SSSP and the Materials Cloud - so nothing that <br>
has been tested and recommended.<br>
<br>
2) you mention "increasing the mixing beta to 0.9"; but that's exactly <br>
the opposite of what you should do, i.e. use very small mixing betas, <br>
even 0.1 or 0.05, and play around with the different Thomas-Fermi <br>
approaches. If you are new to these kind of calculations, you might need <br>
to consider spending some time learning the basics (a very minimal intro <br>
is here: <br>
<a href="https://www.materialscloud.org/learn/sections/spVa3j/fireside-chats-for-lockdown-times)">https://www.materialscloud.org/learn/sections/spVa3j/fireside-chats-for-lockdown-times)</a><br>
<br>
3) maybe try collinear calculations first, to make sure you have a <br>
handle on running calculations, and double check your inputs against the <br>
QE input generator (<a href="https://www.materialscloud.org/work/tools/options">https://www.materialscloud.org/work/tools/options</a>)
<br>
to make sure you are in the right ballpark.<br>
<br>
HTH<br>
<br>
nicola<br>
<br>
<br>
On 14/02/2025 04:09, Johnson, Miles R. wrote:<br>
> Hello,<br>
> <br>
> I'm trying to use the PSlib pseudopotentials for an scf calculation on <br>
> an antiferromagnetic material, including spin-orbit, as the materials <br>
> cloud database and quantum espresso website seem to imply these are the <br>
> best. Unfortunately, the absolute magnetization has immediately diverged <br>
> in every computation I've tried, and I'm not sure how to stop it.<br>
> <br>
> Using the pseudodojo pseudopotentials it converges to around 6 bohr/unit <br>
> cell. With the PSlib pseudopotentials even by the first or second <br>
> iteration it is at like 1000 bohr/unit cell, and never recovers from <br>
> there. I've tried setting startingwfc to 'atomic', and this makes the <br>
> first iteration have a reasonable absolute magnetization, but then in <br>
> the second iteration it already increases to like 100 bohr/unit cell and <br>
> diverges from there. Increasing the mixing beta to 0.9, it still <br>
> diverges very quickly.<br>
> <br>
> The total energy also lowers by hundreds of Ry as the magnetization <br>
> diverges. This happens with both the fully reletavistic PBE and PBEsol <br>
> pseudopotentials. I've also tried modifying input_dft a bit, but the <br>
> functionals that I've tried either aren't compatible with noncolinear <br>
> spin, or haven't changed the divergence.<br>
> <br>
> Does anyone have tips for why this might be diverging so badly, and how <br>
> to get it to converge to a reasonable ground state? I realize I can use <br>
> the pseudodojo potentials, but I'm thinking my bands might not be <br>
> correct so I want to compare how they look for different pseudopotentials.<br>
> <br>
> Best,<br>
> Miles Johnson<br>
> PhD Candidate in Applied Physics<br>
> California Institute of Technology<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> _______________________________________________________________________________<br>
> The Quantum ESPRESSO Foundation stands in solidarity with all civilians worldwide who are victims of terrorism, military aggression, and indiscriminate warfare.<br>
> --------------------------------------------------------------------------------<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)<br>
> users mailing list users@lists.quantum-espresso.org<br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
<br>
-- <br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, SNSF<br>
Laboratory Head, PSI Center for Scientific Computing, Theory, and Data<br>
Contact info and websites at <a href="http://theossrv1.epfl.ch/Main/Contact">http://theossrv1.epfl.ch/Main/Contact</a><br>
</div>
</span></font></div>
</body>
</html>