Hi, I would like calculate hyperfine properties of metals, is possible
this with QE 4.1...??  I read a mail ago, but, in the manual about
GIPAW no appears the job (job='hyperfine') and i'm  don't know the
inputs...!!<br>
<br>
Best.<br>
<br>
Arles V. Gil Rebaza<br>
Intituto de Física de La Plata<br>
La Plata - Argentina.<br><br>---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Davide Ceresoli</b> <span dir="ltr"><<a href="mailto:ceresoli@mit.edu">ceresoli@mit.edu</a>></span><br>Date: 2009/9/29<br>
Subject: Re: [Pw_forum] Hyperfine interaction in Li-Mn compounds?<br>To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br><br><br><div class="im">Gregor Mali wrote:<br>
> Dear PWscf/GIPAW users.<br>
><br>
> I am studying lithium-manganese oxides and lithium-manganese silicates.<br>
> Isotropic shifts of 6Li MAS NMR signals in these compounds are<br>
> predominantly determined by the contact hyperfine interaction (Fermi<br>
> shifts), i.e. by the interaction between 6Li nuclei and unpaired<br>
> electronic spins.<br>
><br>
> Is it possible to calculate these hyperfine shifts or the magnitude of<br>
> hyperfine interaction by the PWsfc/GIPAW module using job = 'hyperfine'.<br>
> I also wonder if the corresponding GIPAW input file requires any<br>
> additional input parameters. The option 'hyperfine' is, namely, not<br>
> documented in the INPUT_GIPAW.html file.<br>
><br>
<br>
</div>Dear Gregor,<br>
     the extra information to enter are the nuclear g-factors<br>
for the isotopes. Here is an example for MnO (55Mn and 18O).<br>
<br>
&inputgipaw<br>
         job = 'hyperfine'<br>
         prefix = 'mno'<br>
         tmp_dir = './scratch/'<br>
         spline_ps = .true.<br>
         hfi_output_unit = 'MHz'<br>
         hfi_nuclear_g_factor(1) =  1.387487      ! 55Mn<br>
         hfi_nuclear_g_factor(2) = -0.757520      ! 18O<br>
/<br>
<br>
I've got the numbers from <a href="http://webelements.com/" target="_blank">webelements.com</a>, under "NMR properties".<br>
There, the gyromangnetic ratio \gamma is listed in 10^7 rad T^{-1} s^{-1}.<br>
To obtain the g-factor, multiply \hbar and divide by the nuclear<br>
magneton (\mu_N):<br>
<br>
E.g. using the command 'units':<br>
You have: 6.6452546e7 hbar / nuclearmagneton T s<br>
You want:<br>
         Definition: 1.3874872<br>
<br>
If the nuclear spin is N/2, N>1, you need to divide the output<br>
quantities by N.<br>
<br>
For Li, you should get good results for the Fermi contact, because<br>
the unpaired electron are in the 2s orbital. But if the magnetization<br>
comes form p and d orbitals, I've found that the core electrons (those<br>
not included in the pseudopotential!) polarize opposite to the<br>
valence, giving rise to large cancellations in the Fermi contact.<br>
<br>
I've implemented the core relaxation in an experimental version of<br>
the code. I haven't put it in the CVS, because it's not perfect, and<br>
one should carefully check if the results make sense. If you are<br>
interested, send me an e-mail.<br>
<br>
> Best regards.<br>
><br>
> Gregor Mali<br>
Best,<br>
<font color="#888888">     Davide</font><br clear="all"><br>-- <br>###--------->   Arles V.   <---------###<br>