<div dir="ltr">If you want to model antiferromagnetic structure, you should define two types of atoms, Gd1 and Gd2. <div style>Change ntyp to 2. Use, e.g.,</div><div style> <span style="font-family:arial,sans-serif;font-size:12.727272033691406px">starting_magnetization(1)=1.0</span></div>
<div style><span style="font-family:arial,sans-serif;font-size:12.727272033691406px">starting_magnetization(2)=-1.0</span><span style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></span></div><div style>
<span style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></span></div><div style><span style="font-family:arial,sans-serif;font-size:12.727272033691406px">ATOMIC_SPECIES</span><br style="font-family:arial,sans-serif;font-size:12.727272033691406px">
<span style="font-family:arial,sans-serif;font-size:12.727272033691406px"> Gd1  1.0  Gd.GGA-PBE-paw.UPF</span><span style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></span></div><div style><span style="font-family:arial,sans-serif;font-size:12.727272033691406px"> Gd2  1.0  Gd.GGA-PBE-paw.UPF</span><span style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br>
</span></div><div style><span style="font-family:arial,sans-serif;font-size:12.727272033691406px">ATOMIC_POSITIONS {crystal}</span><br style="font-family:arial,sans-serif;font-size:12.727272033691406px"><span style="font-family:arial,sans-serif;font-size:12.727272033691406px">Gd1  0.33333333  0.66666667  0.25</span><br style="font-family:arial,sans-serif;font-size:12.727272033691406px">
<span style="font-family:arial,sans-serif;font-size:12.727272033691406px">Gd2 0.66666667  0.33333333  0.75</span><span style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></span></div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Wed, Mar 13, 2013 at 11:05 AM, 琨陶 <span dir="ltr"><<a href="mailto:taokun76@gmail.com" target="_blank">taokun76@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
Yes, the bulk Gd without Hubbard U is antiferromagnetic. If I set the<br>
two atoms in the unit cell to be antiferromagnetic, the "total<br>
magnetization" should be equal to zero, and the "absolute<br>
magnetization" is what I have gotten. Now, however, the "total<br>
magnetization" is not zero. Could you give me some suggestions?<br>
Thanks!<br>
<br>
With regards,<br>
Kun Tao<br>
<div class="HOEnZb"><div class="h5"><br>
On 13 March 2013 10:48, Stefano de Gironcoli <<a href="mailto:degironc@sissa.it">degironc@sissa.it</a>> wrote:<br>
> are you sure you are defining the correct magnetic structure ?<br>
> is Gd maybe antiferromagnetic ?<br>
> stefano<br>
><br>
> On 03/13/2013 03:19 PM, 琨陶 wrote:<br>
>> Dear everyone,<br>
>><br>
>> I want to calculate the magnetic properties of bulk Gd whose magnetic<br>
>> moment is about 8 Bohr magneto/atom. Usually, in the output file<br>
>> values of "total magnetization" and "absolute magnetization" should be<br>
>> nearly the same. However, in my output file, I found a big difference<br>
>> between them as following:<br>
>><br>
>>       total magnetization       =     8.00 Bohr mag/cell<br>
>>       absolute magnetization    =    15.72 Bohr mag/cell<br>
>><br>
>> The good news is that the "absolute magnetization" is close to the<br>
>> experimental results (7.86 Bohr mag/atom), the bad news is that the<br>
>> "total magnetization" is wrong. Moreover, when I try to calculate the<br>
>> magnetic moment on each atom with projwfc.x program, I got about 4.1<br>
>> Bohr mag/atom.<br>
>><br>
>> Could you give me some suggestions? Any advice are appreciated, thanks!<br>
>><br>
>> Regards,<br>
>> Kun Tao<br>
>><br>
>><br>
>> Attached below is the input file:<br>
>><br>
>>   &control<br>
>>      calculation='scf'<br>
>>      restart_mode='from_scratch',<br>
>>      prefix='Gd',<br>
>>      pseudo_dir = './'<br>
>>      outdir='./',<br>
>>   /<br>
>>   &system<br>
>>      ibrav= 4, celldm(1)= 6.96, celldm(3)= 1.59,<br>
>>      nat= 2, ntyp=1,<br>
>>      nspin = 2,<br>
>>      starting_magnetization(1)=1.0,<br>
>>      ecutwfc=55.0,<br>
>>      ecutrho=220.0,<br>
>>      occupations ='smearing',<br>
>>      smearing ='gauss',<br>
>>      degauss = 0.01,<br>
>>   /<br>
>>   &electrons<br>
>>      diagonalization='david',<br>
>>      mixing_beta = 0.3,<br>
>>      conv_thr = 1.0d-8,<br>
>>   /<br>
>> ATOMIC_SPECIES<br>
>>   Gd  1.0  Gd.GGA-PBE-paw.UPF<br>
>> ATOMIC_POSITIONS {crystal}<br>
>> Gd  0.33333333  0.66666667  0.25<br>
>> Gd  0.66666667  0.33333333  0.75<br>
>> K_POINTS {automatic}<br>
>> 4 4 4 0 0 0<br>
>> _______________________________________________<br>
>> Pw_forum mailing list<br>
>> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
>> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></div></div></blockquote></div><br></div>