<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>Dear Konrad Gruszka</div>
<div><br>
</div>
<div>This is a difficult question in fact.</div>
<div>There are to my knowledge three ways to extract Jij from DFT calculation:</div>
<div><br>
</div>
<div>i) Your approach that consist in doing several DFT calculation of various spin (usually collinear) configurations and mapping to Heinsenberg by inverting a series of linear equation.</div>
<div>ii) The explicit approach that uses Liechtenstein formula (not implemented in QE)</div>
<div>iii) The spin-spiral approach (based on the generalized Bloch theorem) that consists in calculation E(q) for various spin sipral vectors q and then fitting the E(q) to get the Jij.</div>
<div><br>
</div>
<div>All these approaches have their advantages and inconveniences but they all rely on the validity of Heisenberg Hamiltonian.</div>
<div><br>
</div>
<div>If in your case some magnetic configurations lead to solutions with a strong decrease of the magnetic moment that means that you are no longer in a pure Heisenberg picture.</div>
<div><br>
</div>
<div>Either you ignore these configurations or you generalize your Hamiltonian by adding a Landau-like term AM^2+BM^4 for example to reproduce the variation of magnetic amplitude...</div>
<div><br>
</div>
<div>There are some publications using this generalized Hamiltonian.</div>
<div><br>
</div>
<div>Cyrille<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="divRpF573343" style="direction: ltr;"><font size="2" face="Tahoma" color="#000000"><b>De :</b> users [users-bounces@lists.quantum-espresso.org] de la part de kgruszka [kgruszka@wip.pcz.pl]<br>
<b>Envoyé :</b> mercredi 5 juin 2019 13:52<br>
<b>À :</b> users@lists.quantum-espresso.org<br>
<b>Objet :</b> [QE-users] Technical questions about calculation of exchage coupling constant J in Heisenberg model<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">Dear users, <br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Recently I'm trying to do the so called mapping of DFT energies to the Heisenberg spin hamiltonian in order to calculate spin coupling constants J(i, j). To do so, I need to calculate Total energy for system with specific various spin configurations,
that I should force my system to be in. Then the procedure is to solve equations obtained in that way in respect to J(i, j). In those eq. spin i or j can have only 2 values: up (1) and down (-1).
<br>
</div>
<div dir="ltr">The main question is how to force specific spin configurations? For simplicity let consider BCC iron.
<br>
</div>
<div dir="ltr">I define in input supercell for BCC iron with different species (fe1, fe2 etc and nspin=2) so any site is defined separately and each different site can have diffetent spin. Then, I define starting_magnetisation for each site for eg. starting_magnetisation(1)=1,
starting_magnetisation(2)=-1 etc, for lets say antiferromagnetic configuration. The problem is, that during SCF cycles I can see spin flips, sometimes spin magnitude is going close to 0 (I conidered not only BCC iron) ruining my starting configuration,
so my system evolves (as it is supposed to) searching for the lowest energy state, that is not always my enforced state.
<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Therefore, as I know that this approach is quite bad, could you please suggest the correct way to do this mapping? Searching for local energy minima in metastable spin configurations seems to be some way, but I hope that there is another,
quicker/better/correct way; how to "lock" my spin configuration? <br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Also, is in case of obtaining J(i, j) noncollinear spin calculation a neccesity? My another doubt is that in case of low spin moment for example = - 0.2 uB can I treat it like spin down configuration in model Hamiltonian, or do I need exact
-1 uB, and what in case of magnitude greater/lower than 1/-1?. <br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">As you can see by my questions I am really confused and I need some basic understanding...
<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Hope I didn't offended any one by my very basic questions, but despite my efforts, I did not find an article describing the technicalities of this procedure, only some general considerations.<br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Regards, <br>
</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">Konrad Gruszka</div>
</div>
</div>
</div>
</body>
</html>