<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">First of all thanks a lot for your suggestions. <br>As you mentioned that the super cell atomic positions are not correct then why the code is not giving error? <br><br>Coming to the magnetic moment querry, as you said the d electrons get filled and thus loosing the magnetic moment. But then one can expect that the moment may be trasffered to some other atom in the cell. But this also in not true in my result. Can this error be because i am doing the scf calculations for gamma point only and then the DOS calculations for 2*2*2 k points mesh? This is the only thing coming to my mind where i can be wrong. Is the difference in k point mesh give such a big difference in the magnetization value, total magnetization is > 1 and local moment is zero? <br><br>Please suggest.<br>Regards<br><br>Padmaja Patnaik<br>
Research Scholar<br>
Dept of Physics<br>
IIT Bombay<br>
Mumbai, India<br><br>-<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div class="plainMail"><br>------------------------------<br><br>Message: 3<br>Date: Mon, 05 Jul 2010 12:43:52 +0200<br>From: Gabriele Sclauzero <<a ymailto="mailto:sclauzer@sissa.it" href="/mc/compose?to=sclauzer@sissa.it">sclauzer@sissa.it</a>><br>Subject: Re: [Pw_forum] Total magnetization and local magnetic moment<br>To: PWSCF Forum <<a ymailto="mailto:pw_forum@pwscf.org" href="/mc/compose?to=pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br>Message-ID: <<a ymailto="mailto:4C31B768.80500@sissa.it" href="/mc/compose?to=4C31B768.80500@sissa.it">4C31B768.80500@sissa.it</a>><br>Content-Type: text/plain; charset="iso-8859-1"<br><br>On 07/05/2010 12:05 PM, Padmaja Patnaik wrote:<br>> I have proceeded with the supercell which is giving correct DOS. Here <br>> is the value of total magnetizatin from the scf.out file
without <br>> relaxation.<br>><br>> ! total energy = -1494.09037638 Ry<br>> Harris-Foulkes estimate = -1494.09037639 Ry<br>> estimated scf accuracy < 6.7E-09 Ry<br>><br>> The total energy is the sum of the following terms:<br>><br>> one-electron contribution = 257.19918198 Ry<br>> hartree contribution = 202.08158140 Ry<br>> xc contribution = -548.49776710 Ry<br>> ewald contribution = -1404.85946098 Ry<br>> smearing contrib. (-TS) = -0.01391169 Ry<br>><br>> total
magnetization = 1.02 Bohr mag/cell<br>> absolute magnetization = 1.12 Bohr mag/cell<br>><br>> convergence has been achieved in 17 iterations<br>> ------------------------------------------------------------------------<br>> Here is the portion of the file *.pdos.out without relaxation<br>><br>> Atom # 1: total charge = 14.3653, s = 2.4622, p = 6.8403, d = <br>> 5.0628,<br>> spin up = 7.7466, s = 1.2349, p = 3.4229, d <br>> = 3.0889,<br>> spin down = 6.6187, s = 1.2274, p = 3.4174, d <br>> =
1.9739,<br>> polarization = 1.1279, s = 0.0075, p = 0.0054, d <br>> = 1.1150,<br>> Atom # 2: total charge = 2.9244, s = 0.8587, p = 2.0656, d <br>> = 0.0000,<br>> spin up = 1.4798, s = 0.4351, p = 1.0446, d <br>> = 0.0000,<br>> spin down = 1.4446, s = 0.4236, p = 1.0210, d <br>> = 0.0000,<br>> polarization = 0.0352, s = 0.0115, p = 0.0236, d <br>> = 0.0000,<br>> Atom # 3: total charge
= 2.9244, s = 0.8587, p = 2.0656, d <br>> = 0.0000,<br>> spin up = 1.4798, s = 0.4351, p = 1.0446, d <br>> = 0.0000,<br>> spin down = 1.4446, s = 0.4236, p = 1.0210, d <br>> = 0.0000,<br>> polarization = 0.0352, s = 0.0115, p = 0.0236, d <br>> = 0.0000,<br>> -----------------------------------------------------------------------------------------------------------------<br>> As you can see atom number 1(this is the impurity) is showing a local <br>> magnetic moment. Now, below i add portions of the same two files from <br>> with relaxation calculations.<br>> !
total energy = -1494.15610974 Ry<br>> Harris-Foulkes estimate = -1494.15610974 Ry<br>> estimated scf accuracy < 1.9E-10 Ry<br>><br>> The total energy is the sum of the following terms:<br>><br>> one-electron contribution = 270.52901167 Ry<br>> hartree contribution = 195.52169853 Ry<br>> xc contribution = -548.43270406 Ry<br>> ewald contribution = -1411.76028316 Ry<br>> smearing contrib. (-TS) = -0.01383273 Ry<br>><br>> total magnetization =
1.13 Bohr mag/cell<br>> absolute magnetization = 1.26 Bohr mag/cell<br>><br>> convergence has been achieved in 9 iterations<br>><br>><br><br><br>It looks like your impurity (Ge, I guess) acquires a lot of charge upon <br>relaxing the structure, so that the d electrons get filled and you lose <br>the local moment. This may be a real description or some spurious effect <br>due to orbitals overlap (though it's unlikely if these are very <br>localized 3d electrons).<br><br>About the total magnetization, well I don't know. It may come from <br>neighboring atoms if they have donated some charge to the impurity.<br><br>HTH<br><br>GS<br><br>> ------------------------------------------------------------------------------------------------------------<br>> Atom # 1: total charge = 19.3947, s = 2.5253, p
= 7.6414, d = <br>> 9.2280,<br>> spin up = 9.7222, s = 1.2688, p = 3.8289, d <br>> = 4.6245,<br>> spin down = 9.6725, s = 1.2566, p = 3.8125, d <br>> = 4.6035,<br>> polarization = 0.0496, s = 0.0122, p = 0.0164, d <br>> = 0.0210,<br>> Atom # 2: total charge = 3.0599, s = 0.9045, p = 2.1555, d <br>> = 0.0000,<br>> spin up = 1.5315, s = 0.4527, p = 1.0787, d <br>> = 0.0000,<br>>
spin down = 1.5285, s = 0.4517, p = 1.0767, d <br>> = 0.0000,<br>> polarization = 0.0030, s = 0.0010, p = 0.0020, d <br>> = 0.0000,<br>> Atom # 3: total charge = 3.0599, s = 0.9045, p = 2.1555, d <br>> = 0.0000,<br>> spin up = 1.5315, s = 0.4527, p = 1.0787, d <br>> = 0.0000,<br>> spin down = 1.5285, s = 0.4517, p = 1.0767, d <br>> = 0.0000,<br>> polarization = 0.0030, s
= 0.0010, p = 0.0020, d <br>> = 0.0000,<br>> ----------------------------------------------------------------------------------------------------------------------<br>> Here as you can see the total magnetization increases but the atom <br>> number 1 (impurity) is showing zero local magnetic moment.<br>> I hope i made the problem more clear this time.<br>><br>> Regards<br>><br><br><br>-- <br><br>Gabriele Sclauzero, EPFL SB ITP CSEA<br>PH H2 462, Station 3, CH-1015 Lausanne<br><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://www.democritos.it/pipermail/pw_forum/attachments/20100705/f9d2e088/attachment.htm" target="_blank">http://www.democritos.it/pipermail/pw_forum/attachments/20100705/f9d2e088/attachment.htm</a> <br><br>------------------------------<br><br>_______________________________________________<br>Pw_forum mailing list<br><a
ymailto="mailto:Pw_forum@pwscf.org" href="/mc/compose?to=Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br><br><br>End of Pw_forum Digest, Vol 37, Issue 11<br>****************************************<br></div></blockquote></td></tr></table><br>