<div dir="ltr">Dear Prof. Blaha and other wien2k users:<br><br>(I posted a similar message yesterday, apologies in case this appears as a repeat; the first message has not appeared on the list, perhaps reflected due to included images.)<br>
<br>Regarding tests of the hyperfine fields in aluminum metal, we had thought about the issue of insufficient k-points, however we thought we had a handle on this issue. In a 9 T field, a rough calculation shows that the thin spin-polarized shell at Ef represents about 1/3000 of the BZ volume for fcc-Al. We ran a script gradually increasing the number of k-points, with a result (shown in <a href="http://people.physics.tamu.edu/jhchen/points.png">http://people.physics.tamu.edu/jhchen/points.png</a>) that the HFF settles down within about 20% of the expected value for 10,000 k-points in B=9T, with fluctuations dying down to the order of 10% and less in the range 30,000 - 80,000 k-points. We also ran a test for linearity in B at a setting of 10,000 k-points, and the results appeared to be quite linear up to 100 T (shown in <a href="http://people.physics.tamu.edu/jhchen/field.png">http://people.physics.tamu.edu/jhchen/field.png</a>).<br>
<br>We ran the test treating fcc-Al as simple cubic with 4 sites in order to be sure we understood how the field is applied in ORB, and expected if anything better convergence since the expanded cell gives a greater k-point density. However the results seem strange: with several k-point settings we found that in general, the HFF approached the expected value for fcc-Al after a relatively small number of iterations, yet without quite converging, and finally the HFF values diverged, with one or more going large and negative. We had not tried as many variations as for fcc since the results are much slower to obtain converged HFF.<br>
<br>Following the suggestion of Prof. Blaha after our last posting we tried increasing to very large field and k-point values, and did finally get convergence (more than 10 last iterations of HFF is the same) for a setting of 100000 k-points and 10000 T, yielding 4 reasonably close positive values as in the following:<br>
<br>------<br>:HFF001: 143.345 0.000 0.572 143.917 (KGAUSS)<br>:HFF002: 143.344 0.000 0.572 143.916 (KGAUSS)<br>:HFF003: 144.427 0.000 0.583 145.010 (KGAUSS)<br>
:HFF004: 143.344 0.000 0.572 143.916 (KGAUSS)<br>------<br><br>However we are concerned that the HFF values are still not identical, whereas at 10,000 T the spin-polarized shell at Ef represents a significant fraction of the BZ, and the spin energy is quite large. We expected this to be more than enough k-points for random sampling of the shell at Ef. For this reason, and in particular in light of the strange behavior in which the HFF values almost converge before diverging to widely separated values, is it possible that there might be some other issue that we are overlooking?<br>
<br>Any suggestions would be appreciated.</div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/10/7 Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">The hyperfine field for a metal is coming mainly from the contact term due to the induced spin-polarization by the magnetic field.<br>
<br>
You should notice, that a field of 9 T is (for theoretical calculations) an extremely small field, causing a very small spin-splitting of the states near EF, which causes the HFF.<br>
I suppose all you see is numerical noise.<br>
<br>
Since only the states at EF are of interest (the field can only reoccupy states within a few mRy (or less) around EF), you need to converge your calculation with respect to:<br>
<br>
a) the k-mesh (test MUCH larger meshes (10000, 50000 100000 k or more)<br>
b) the magnetic field (increase it and test fields up to 1000 T), You are not interested in the absolute number, but in ppm, i.e. the relative induced field.<br>
<br>
c) The angular momentum component of the hFF introduced by case.vorbup/dn is NOT correct. I would even suggest that you put l=0 to<br>
minimize the effect (or use -orbc with case.vorbup/dn , where all elements are set to zero.)<br>
<br>
d) In principle the orbital contribution should be obtainable from the NMR-module of wien2k_13. However, also there we observed for metals that it is very hard to converge with respect to k-mesh and the final results (sum of spin and orbital contribution) does not seem right, while spin-only has the correct magnitude (within 10% of the experiment). This is an unresolved issue for us so far.<br>
<br>
<br>
Am 07.10.2013 04:01, schrieb Jing-Han Chen:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Dear WIEN2k users and authors<br>
<br>
We are currently working on the hyperfine field calculation by using<br>
ORB package. In fcc aluminum case, we got 0.154 (KGAUSS) when the<br>
following case.inorb and case.indm are used<br>
<br>
case.inorb<br>
3 1 0 nmod, natorb, ipr<br>
PRATT, 1.0 mixmod, amix<br>
1 1 0 iatom nlorb, lorb<br>
9. Bext in T<br>
0. 0. 1. direction of Bext in terms of lattice vectors<br>
<br>
case.indm<br>
-9. Emin cutoff energy<br>
1 number of atoms for which density matrix is<br>
calculated<br>
1 1 0 index of 1st atom, number of L's, L1<br>
0 0 r-index, (l,s)index<br>
<br>
In order to confirm how the magnetic field is applied for the<br>
multiple sites crystal, we made aluminum as a simple cubic with 4<br>
inequivalent sites and we believe it should be physically identical to<br>
fcc. The following case.inorb and case.indm are used.<br>
<br>
case.inorb<br>
3 4 0 nmod, natorb, ipr<br>
PRATT, 1.0 mixmod, amix<br>
1 1 0 iatom nlorb, lorb<br>
2 1 0 iatom nlorb, lorb<br>
3 1 0 iatom nlorb, lorb<br>
4 1 0 iatom nlorb, lorb<br>
9. Bext in T<br>
0. 0. 1. direction of Bext in terms of lattice vectors<br>
<br>
case.indm<br>
-9. Emin cutoff energy<br>
4 number of atoms for which density matrix is<br>
calculated<br>
1 1 0 index of 1st atom, number of L's, L1<br>
2 1 0 index of 1st atom, number of L's, L1<br>
3 1 0 index of 1st atom, number of L's, L1<br>
4 1 0 index of 1st atom, number of L's, L1<br>
0 0 r-index, (l,s)index<br>
<br>
Both fcc and simple cubic are run by the same way (-orb -cc 0.00001).<br>
A complete different HFFs are obtained as the following<br>
<br>
:HFF001: 0.059 0.000 0.001<br>
0.060 (KGAUSS)<br>
:HFF002: -1.193 0.000 -0.010<br>
-1.204 (KGAUSS)<br>
:HFF003: 1.681 0.000 0.011<br>
1.692 (KGAUSS)<br>
:HFF004: 0.046 0.000 0.001<br>
0.047 (KGAUSS)<br>
<br>
We got four different HFFs which we thought they are supposed to be the<br>
same. Also all of them are very far from the fcc result (0.154 KGAUSS).<br>
Does anyone know why it happens?<br>
<br>
Any suggestion and comment are appreciated.<br>
<br>
--<br>
Jing-Han Chen<br>
Graduate Student<br>
Department of Physics<br>
Texas A&M University<br>
4242 TAMU<br>
College Station TX 77843-4242<br>
</div></div><a href="mailto:jhchen@tamu.edu" target="_blank">jhchen@tamu.edu</a> <mailto:<a href="mailto:jhchen@tamu.edu" target="_blank">jhchen@tamu.edu</a>> <<a href="mailto:jhchen@tamu.edu" target="_blank">jhchen@tamu.edu</a><br>
<mailto:<a href="mailto:jhchen@tamu.edu" target="_blank">jhchen@tamu.edu</a>>> / <a href="http://people.physics.tamu.edu/jhchen/" target="_blank">http://people.physics.tamu.<u></u>edu/jhchen/</a><br>
<br>
<br>
<br>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>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" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
<br>
</blockquote>
<br>
-- <br>
Peter Blaha<br>
Inst.Materials Chemistry<br>
TU Vienna<br>
Getreidemarkt 9<br>
A-1060 Vienna<br>
Austria<br>
<a href="tel:%2B43-1-5880115671" value="+4315880115671" target="_blank">+43-1-5880115671</a><br>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>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" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Jing-Han Chen<br>Graduate Student<br>Department of Physics<br>Texas A&M University<br>4242 TAMU<br>College Station TX 77843-4242<br><a href="mailto:jhchen@tamu.edu">jhchen@tamu.edu</a> <<a href="mailto:jhchen@tamu.edu">jhchen@tamu.edu</a>> / <a href="http://people.physics.tamu.edu/jhchen/">http://people.physics.tamu.edu/jhchen/</a><br>
</div>