<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Times New Roman;color: #000000;font-size: 10pt;">
<div>Hi Marcos</div>
<div><br>
</div>
<div>The implementation of the Force Theorem has been described in the following paper:<br>
</div>
<div><a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.90.205409" target="_blank">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.90.205409</a></div>
<div><br>
</div>
<div>The procedure is the following:</div>
<div>first perform a scf calculation with scalar relativistic pseudo</div>
<div>then perform  nscf calculation with fully relativistic pseudo (option lforcet=.true., nosym=.true') starting from previous scf charge  (startingpot='file') with various spin orientations (theta=0,90 for example)</div>
<div>Finally perform a projwfc calculation with  lforcet=.true. and the value of the Fermi level from the nscf calculation (same ef_0 for all calculations).<br>
</div>
<div><br>
</div>
<div>Then you get a file with the energy decomposed over the various atoms and orbitals of the system..</div>
<div>The anisotropy is obtained by difference between the different spin orientations.<br>
</div>
<div><br>
</div>
<div>hope it helps..</div>
<div><br>
</div>
<div>Cyrille<br>
</div>
<div><br>
</div>
<div><br>
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px"><font size="2">========================<br>
</font>
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px"><font size="2"></font>
<div style="font-family:Tahoma; font-size:13px"><font size="2"></font>
<div style="font-family:Tahoma; font-size:13px"><font size="2"></font>
<div style="font-family:Tahoma; font-size:13px"><b>Cyrille Barreteau</b><br>
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px"><font size="2"><font face="Times New Roman">CEA Saclay<span style="background-color:white"></span></font><font face="Times New Roman">,
</font><font face="Times New Roman"><font face="Times New Roman"><font face="Times New Roman">IRAMIS,</font> SPEC</font></font><font face="Times New Roman"><font face="Times New Roman"><span style="background-color:white"> Bat. 771<br>
</span></font>91191 Gif sur Yvette Cedex<span style="background-color:white"></span>, FRANCE 
</font><br>
<font face="Times New Roman"><font face="Times New Roman">~~~~~~~~~~~~~~~~~~~~~~~~~~~~</font><br>
</font><font face="Times New Roman"><font face="Times New Roman"><font face="Times New Roman">+33 1 69 08 38 56 /</font></font></font><font face="Times New Roman"><font face="Times New Roman"><font face="Times New Roman"><font face="Times New Roman"><font face="Times New Roman">+33 
 6 47 53 66 52  (mobile)  </font></font></font> </font></font></font><br>
<font face="Times New Roman">email:     cyrille.barreteau@cea.fr </font><font face="Times New Roman"><br>
Web:     <a href="http://iramis.cea.fr/Pisp/cyrille.barreteau/" tabindex="0">http://iramis.cea.fr/Pisp/cyrille.barreteau/</a></font><font size="2"><br>
</font><font size="2">========================</font></div>
</div>
</div>
</div>
<span style="background-color:white; font-family:Times New Roman"><font size="2" color="black"><span dir="ltr" style="font-size:10pt"><font size="2"><span style="font-size:16px"></span></font></span></font></span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF200782" style="direction: ltr;"><font size="2" face="Tahoma" color="#000000"><b>De :</b> users [users-bounces@lists.quantum-espresso.org] de la part de Paolo Giannozzi [p.giannozzi@gmail.com]<br>
<b>Envoyé :</b> jeudi 12 juillet 2018 10:04<br>
<b>À :</b> Quantum Espresso users Forum<br>
<b>Objet :</b> Re: [QE-users] Magnetic anisotropy energy in QE 6.3<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div>here? PP/examples/ForceTheorem_example/<br>
<br>
</div>
P.<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Jul 12, 2018 at 10:01 AM, Marcos Veríssimo Alves
<span dir="ltr"><<a href="mailto:marcos.verissimo.alves@gmail.com" target="_blank" rel="noopener noreferrer">marcos.verissimo.alves@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
<div dir="auto">Hi all,
<div dir="auto"><br>
</div>
<div dir="auto">Browsing QE 6.3's documentation, I saw that MAE can be calculated as a post-processing step to a pw.x scf calculation. What is the exact procedure followed? I.e., in the post-processing calculation is the spin density rotated, and then SOC is
 included? Is there any reference that details the procedure used when the MFT is applied in Espresso?</div>
<div dir="auto"><br>
</div>
<div dir="auto">Best regards,</div>
<div dir="auto"><br>
</div>
<div dir="auto">Marcos</div>
</div>
<br>
______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noopener noreferrer">users@lists.quantum-espresso.<wbr>org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-<wbr>espresso.org/mailman/listinfo/<wbr>users</a><br>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
<div class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>