<div dir="ltr"><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">Thanks Dr Gergard or the full explanation <br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">I want to deterine the exchange splitting which allows me to determine the Stoner parameter value I</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">If we take the example of Fe, in the paramagnetic state the number of up and the dn electrons equal both 4 and the DOS bands for up and Dn spins are fully symmetrized whereas in the ferromagnetic state a double shift occurs for the Up and Dn bands ( a shift of the Up bands in the left side and the number of the up electrons becomes 5.1 ( the number of the d elctrons is only 4.62 ) and a shift of the Dn bands in the right side and the number of dn electrons becomes 2.90 ( the number of d electrons is only 2.42 ).</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">I think that the exchange splitting energy as Pr Peter said is the sum of the energy shift of the Dn bands in the right side and the energy shift of the Up bands in the left side.</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">Is my conclusion correct or not ? and if is it correct how to determine this exchange splitting ?</div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:comic sans ms,sans-serif;font-size:small">Best regards<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 26 Apr 2020 at 09:21, Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Actually it is not clear what you want to do, and for what purpose !<br>
... and answer to "Is correct ?" : NO, the exchange splitting is not the energy difference between two maxima of the density of states.<br>
<br>
Check the band structure of Fe<br>
you see that there are in the majority case three bands (states) occupied at the Gamma point, <br>
these are from bottom (lowest energy) a1g (s), t2g (dxy, dxz, dyz) , and eg (dz^2, dx^2-y^2)<br>
in the minority case only the a1g and t2g are at Gamma below the Fermi energy.<br>
(the degeneracy is seen from the dispersion in the Sigma direction)<br>
<br>
You should easily find out that there is no unique splitting, but it depends on the state, <br>
and as I already mentioned in earlier mail you have different splittings (including not only exchange splitting) depending on the k-points and their symmetry.<br>
<br>
The energies of the bands, or better the ranges of the band energies are (roughly, depending on the k-mesh) given in the case.scf2up/dn files.<br>
<br>
To find the correct energies at Gamma you van either have a look in the spaghetti files or better use irrep.<br>
(Note their might be changes in the band series if you have a more complicated system,<br>
and then you have to use irrep do find the correct band character)<br>
The case.outputirruo/dn contains all informations, in a shortened way you will find for iron something like <br>
(note that I cut the values after 3 digits)<br>
<br>
for the states at Gamma (using irrep, eigval in Ry, all other information skipped):<br>
from up from down<br>
bnd ndg eigval eigval diff * 13.605 my info: band character<br>
4 1 -0.021 .......... =G1+ 0.005 .......... =G1+ ==> 0.22 eV A1g (s)<br>
5 3 0.436 .......... =G5+ 0.57 .......... =G5+ ==> 1.91 eV T2g (dxy, dxz, dyz)<br>
8 2 0.525 .......... =G3+ 0.728 .......... =G3+ ==> 2.75 eV Eg (dz^2, dx^2-y^2)<br>
<br>
you see that there is a very different splitting that depends on the character of the state !<br>
(Principally one could calculate the average (1.9eV all, 2.2eV only d, note the degeneracy !),<br>
but I think that this is an oversimplification and might be used for very rough comparison only.)<br>
<br>
To use single maxima of the density of states for determination of the splitting does not make any sense, <br>
just plot the density and shift either up or dn DOS by the value you found to see that it doesn't tell anything.<br>
(The idea of Peter would only work if the up/dn densities are identical with only a rigid shift, what is usually never the case.)<br>
<br>
Another way may be to try to find the average energy of the bands, as I mentioned before, but this will only give a very rough idea as it includes not only<br>
the exchange but also the crystal field and possible other splittings. Further one has to take care on suitable energy ranges (see case.scf2) and normalization.<br>
The difference in the band energies for majority and minority electrons, however, is still a good measure for the magnetic energy.<br>
<br>
You where also asking for " the exchange splitting between the paramagnetic and ferromagnetic states", it is not clear what you mean ?<br>
Do you mean the shift of the states between non-spinpolarized and spinpolarized up, and between non-spinpolarized and spinpolarized dn ?<br>
Again, there is no unique shift of the states but you may compare the energies at Gamma similar to what I did above.<br>
<br>
I suggest to read and understand the textbook of J Kübler "Theory of itinerant electron magnetism", he lists the splittings at different <br>
high symmetry points of the Brillouin zone (for different calculaton schemes compared to experiment) for iron and cobalt.<br>
<br>
Ciao<br>
Gerhard<br>
<br>
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>
"I think the problem, to be quite honest with you,<br>
is that you have never actually known what the question is."<br>
<br>
====================================<br>
Dr. Gerhard H. Fecher<br>
Institut of Inorganic and Analytical Chemistry<br>
Johannes Gutenberg - University<br>
55099 Mainz<br>
and<br>
Max Planck Institute for Chemical Physics of Solids<br>
01187 Dresden<br>
________________________________________<br>
Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>] im Auftrag von Abderrahmane Reggad [<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a>]<br>
Gesendet: Freitag, 24. April 2020 16:42<br>
An: A Mailing list for WIEN2k users<br>
Betreff: Re: [Wien] How to determine the exchange splitting in metals ?<br>
<br>
Hello again<br>
I have adopted another procedure as follows:<br>
<br>
I have extracted the energy values corresping to the maximum values of DOS in the up and down spin from the files case.dos1evup and case.dos1evdn and I have found the following values:<br>
<br>
E (max DOS up = 2.40) = - 0.95185 eV<br>
E (maxDOS dn = 2.80) = + 1.95977 eV<br>
<br>
and I have calculated the exchange splitting from this formula dE= + 1.95977 - ( - 0.95185) = 2.91162 eV<br>
<br>
Is correct ?<br>
<br>
<br>
<br>
On Fri, 24 Apr 2020 at 14:54, Abderrahmane Reggad <<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a><mailto:<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a>>> wrote:<br>
Thanks Gerhard for the explanation but I couldn't apply the inforation to get the exchange splitting neither from the DOS nor from the band structure<br>
<br>
Now i have some questions about my idea using the DOS picture and I want from both of you to answer me<br>
<br>
- is it possible from the files case.outputtup and case.outputtdn to get the energies corresponding to the integrated DOS values for spin up (5.1 e) and spin down (2.9 e) ? in the case o iron Fe<br>
- Is it possible to determine the exchange splitting as follows: dE = E (5.1 e) - E (2.9 e) in abslute value<br>
<br>
On Fri, 24 Apr 2020 at 14:52, Abderrahmane Reggad <<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a><mailto:<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a>>> wrote:<br>
Thanks Gerhard for the explanation but I couldn't apply the inforation to get the exchange splitting neither from the DOS nor from the band structure<br>
<br>
Now i have some questions about my idea using the DOS picture and I want from both of you to answer me<br>
<br>
- is it possible from the files case.outputtup and case.outputtdn to get the energies corresponding to the integrated DOS values for spin up (5.1 e) and spin down (2.9 e) ? in the case o iron Fe<br>
- Is it possible to determine the exchange splitting as follows: dE = E (5.1 e) - E (2.9 e) in abslute value<br>
<br>
I have joinded the case.outputt files or the paraagnetic and ferromagnetic state of iron Fe<br>
<br>
Best regards<br>
<br>
<br>
<br>
On Fri, 24 Apr 2020 at 12:26, Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de" target="_blank">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de" target="_blank">fecher@uni-mainz.de</a>>> wrote:<br>
maybe use irrep to see how much the bands at the Gamma point are splitted between up and down<br>
(single k-point no shift of BZ before calculating lapw1 -up, -dn; but be carefull which states at Gamma you compare)<br>
Note the splitting depends on k, what you easily see from the bands, therfore a comparison of the bands mmight not be very helpful.<br>
<br>
the mean or state resolved splittings may also be calculated by the difference in the band energies for up and down states, that is<br>
integral n(E)up * E dE - integral n(E)down * E dE<br>
where the integrals run over all occupied states of the valence bands, or you use only particular states, e.g. all d pr eg, or t2g only.<br>
(Note the sum is the overall band energy, that you may compare to that of the paramagnetic state, if you wish to do for whatever reason)<br>
<br>
Ciao<br>
Gerhard<br>
<br>
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>
"I think the problem, to be quite honest with you,<br>
is that you have never actually known what the question is."<br>
<br>
====================================<br>
Dr. Gerhard H. Fecher<br>
Institut of Inorganic and Analytical Chemistry<br>
Johannes Gutenberg - University<br>
55099 Mainz<br>
and<br>
Max Planck Institute for Chemical Physics of Solids<br>
01187 Dresden<br>
________________________________________<br>
Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>>] im Auftrag von Abderrahmane Reggad [<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a><mailto:<a href="mailto:abde.reggad@gmail.com" target="_blank">abde.reggad@gmail.com</a>>]<br>
Gesendet: Freitag, 24. April 2020 11:30<br>
An: A Mailing list for WIEN2k users<br>
Betreff: Re: [Wien] How to determine the exchange splitting in metals ?<br>
<br>
Thanks Pr Plaha for the explanation<br>
<br>
Now how to determine it through the band structure<br>
<br>
Best regards<br>
<br>
On Fri, 24 Apr 2020 at 09:14, Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>><mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a><mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>>>> wrote:<br>
I would do this with the band structure (because this could be<br>
k-dependent), but DOS is also fine.<br>
<br>
Just shift up and dn DOS in energy until they overlap as much as<br>
possible. This shift is your exchange splitting.<br>
<br>
Am 24.04.2020 um 00:07 schrieb Abderrahmane Reggad:<br>
> Hello wien2k users<br>
><br>
> I have calculated the DOS of the paramagnetic and ferromagnetic of 3d<br>
> transition metals Ni , Fe and I want to determine the exchange splitting<br>
> between the paramagnetic and ferromagnetic states.<br>
><br>
> How to do that ?<br>
><br>
> Best regards<br>
><br>
> --<br>
> Dr. Abderrahmane Reggad<br>
> Engineering Physics Laboratory<br>
> Faculty of Material Sciences, Ibn Khaldoun University, Tiaret, 14000,<br>
> Algeria<br>
> Tel: +213(0)561861963 - Algeria<br>
><br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>>><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
><br>
<br>
--<br>
Peter Blaha<br>
Inst.Materials Chemistry<br>
TU Vienna<br>
Getreidemarkt 9<br>
A-1060 Vienna<br>
Austria<br>
+43-1-5880115671<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>>><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br>
<br>
--<br>
Dr. Abderrahmane Reggad<br>
Engineering Physics Laboratory<br>
Faculty of Material Sciences, Ibn Khaldoun University, Tiaret, 14000, Algeria<br>
Tel: +213(0)561861963 - Algeria<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br>
<br>
--<br>
Dr. Abderrahmane Reggad<br>
Engineering Physics Laboratory<br>
Faculty of Material Sciences, Ibn Khaldoun University, Tiaret, 14000, Algeria<br>
Tel: +213(0)561861963 - Algeria<br>
<br>
<br>
--<br>
Dr. Abderrahmane Reggad<br>
Engineering Physics Laboratory<br>
Faculty of Material Sciences, Ibn Khaldoun University, Tiaret, 14000, Algeria<br>
Tel: +213(0)561861963 - Algeria<br>
<br>
<br>
--<br>
Dr. Abderrahmane Reggad<br>
Engineering Physics Laboratory<br>
Faculty of Material Sciences, Ibn Khaldoun University, Tiaret, 14000, Algeria<br>
Tel: +213(0)561861963 - Algeria<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</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><div><span style="color:rgb(92,92,92);font-family:Arial,Helvetica,"Lucida Sans Unicode","Microsoft Sans Serif","Segoe UI Symbol",STIXGeneral,"Cambria Math","Arial Unicode MS",sans-serif;font-size:13px;font-style:normal;font-weight:400;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);display:inline;float:none">Dr. Abderrahmane Reggad<br>Engineering Physics Laboratory</span></div><span style="color:rgb(92,92,92);font-family:Arial,Helvetica,"Lucida Sans Unicode","Microsoft Sans Serif","Segoe UI Symbol",STIXGeneral,"Cambria Math","Arial Unicode MS",sans-serif;font-size:13px;font-style:normal;font-weight:400;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);display:inline;float:none"><span></span>Faculty of Material Sciences, Ibn Khaldoun University, Tiaret, 14000, Algeria</span></div>Tel: +213(0)561861963 - Algeria<br></div></div></div></div></div></div>