<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<html>

<head>

  <meta content="text/html; charset=ISO-8859-1"

 http-equiv="Content-Type">

</head>

<body text="#000000" bgcolor="#ffffff">

On 07/05/2010 12:05 PM, Padmaja Patnaik wrote:

<blockquote cite="mid:322092.98137.qm@web26101.mail.ukl.yahoo.com"

 type="cite">

  <table border="0" cellpadding="0" cellspacing="0">

    <tbody>

      <tr>

        <td

 style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit; -x-system-font: none;"

 valign="top">I have proceeded with the supercell which is giving

correct DOS. Here is the value of total magnetizatin from the scf.out

file without relaxation.<br>

        <br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1494.09037638 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; = -1494.09037639 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6.7E-09 Ry<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; The total energy is the sum of the following terms:<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; one-electron contribution =&nbsp;&nbsp; 257.19918198 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; hartree contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 202.08158140 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; xc contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -548.49776710 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; ewald contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1404.85946098 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; smearing contrib. (-TS)&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -0.01391169 Ry<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; total magnetization&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 1.02 Bohr mag/cell<br>

&nbsp;&nbsp;&nbsp;&nbsp; absolute magnetization&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 1.12 Bohr mag/cell<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; convergence has been achieved in&nbsp; 17 iterations<br>

------------------------------------------------------------------------<br>

Here is the portion of the file *.pdos.out&nbsp; without relaxation<br>

        <br>

&nbsp;Atom #&nbsp;&nbsp; 1: total charge =&nbsp; 14.3653, s =&nbsp; 2.4622, p =&nbsp; 6.8403, d =&nbsp;

5.0628,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 7.7466, s =&nbsp; 1.2349, p =&nbsp; 3.4229, d

=&nbsp; 3.0889,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 6.6187, s =&nbsp; 1.2274, p =&nbsp; 3.4174, d

=&nbsp; 1.9739,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp; 1.1279, s =&nbsp; 0.0075, p =&nbsp; 0.0054, d

=&nbsp; 1.1150,<br>

&nbsp;&nbsp;&nbsp;&nbsp; Atom #&nbsp;&nbsp; 2: total charge =&nbsp;&nbsp; 2.9244, s =&nbsp; 0.8587, p =&nbsp; 2.0656, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.4798, s =&nbsp; 0.4351, p =&nbsp; 1.0446, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.4446, s =&nbsp; 0.4236, p =&nbsp; 1.0210, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp; 0.0352, s =&nbsp; 0.0115, p =&nbsp; 0.0236, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp; Atom #&nbsp;&nbsp; 3: total charge =&nbsp;&nbsp; 2.9244, s =&nbsp; 0.8587, p =&nbsp; 2.0656, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.4798, s =&nbsp; 0.4351, p =&nbsp; 1.0446, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.4446, s =&nbsp; 0.4236, p =&nbsp; 1.0210, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp; 0.0352, s =&nbsp; 0.0115, p =&nbsp; 0.0236, d

=&nbsp; 0.0000,<br>

-----------------------------------------------------------------------------------------------------------------<br>

As you can see atom number 1(this is the impurity) is showing a local

magnetic moment. Now, below i add portions of the same two files from

with relaxation calculations.<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1494.15610974 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; = -1494.15610974 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.9E-10 Ry<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; The total energy is the sum of the following terms:<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; one-electron contribution =&nbsp;&nbsp; 270.52901167 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; hartree contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 195.52169853 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; xc contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -548.43270406 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; ewald contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1411.76028316 Ry<br>

&nbsp;&nbsp;&nbsp;&nbsp; smearing contrib. (-TS)&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -0.01383273 Ry<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; total magnetization&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 1.13 Bohr mag/cell<br>

&nbsp;&nbsp;&nbsp;&nbsp; absolute magnetization&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 1.26 Bohr mag/cell<br>

        <br>

&nbsp;&nbsp;&nbsp;&nbsp; convergence has been achieved in&nbsp;&nbsp; 9 iterations<br>

        <br>

        <br>

        </td>

      </tr>

    </tbody>

  </table>

</blockquote>

<br>

<br>

It looks like your impurity (Ge, I guess) acquires a lot of charge upon

relaxing the structure, so that the d electrons get filled and you lose

the local moment. This may be a real description or some spurious

effect due to orbitals overlap (though it's unlikely if these are very

localized 3d electrons). <br>

<br>

About the total magnetization, well I don't know. It may come from

neighboring atoms if they have donated some charge to the impurity.<br>

<br>

HTH<br>

<br>

GS<br>

<br>

<blockquote cite="mid:322092.98137.qm@web26101.mail.ukl.yahoo.com"

 type="cite">

  <table border="0" cellpadding="0" cellspacing="0">

    <tbody>

      <tr>

        <td

 style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit; -x-system-font: none;"

 valign="top">------------------------------------------------------------------------------------------------------------<br>

Atom #&nbsp;&nbsp; 1: total charge =&nbsp; 19.3947, s =&nbsp; 2.5253, p =&nbsp; 7.6414, d =&nbsp;

9.2280,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 9.7222, s =&nbsp; 1.2688, p =&nbsp; 3.8289, d

=&nbsp; 4.6245,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 9.6725, s =&nbsp; 1.2566, p =&nbsp; 3.8125, d

=&nbsp; 4.6035,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp; 0.0496, s =&nbsp; 0.0122, p =&nbsp; 0.0164, d

=&nbsp; 0.0210,<br>

&nbsp;&nbsp;&nbsp;&nbsp; Atom #&nbsp;&nbsp; 2: total charge =&nbsp;&nbsp; 3.0599, s =&nbsp; 0.9045, p =&nbsp; 2.1555, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.5315, s =&nbsp; 0.4527, p =&nbsp; 1.0787, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.5285, s =&nbsp; 0.4517, p =&nbsp; 1.0767, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp; 0.0030, s =&nbsp; 0.0010, p =&nbsp; 0.0020, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp; Atom #&nbsp;&nbsp; 3: total charge =&nbsp;&nbsp; 3.0599, s =&nbsp; 0.9045, p =&nbsp; 2.1555, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.5315, s =&nbsp; 0.4527, p =&nbsp; 1.0787, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 1.5285, s =&nbsp; 0.4517, p =&nbsp; 1.0767, d

=&nbsp; 0.0000,<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp; 0.0030, s =&nbsp; 0.0010, p =&nbsp; 0.0020, d

=&nbsp; 0.0000,<br>

----------------------------------------------------------------------------------------------------------------------<br>

Here as you can see the total magnetization increases but the atom

number 1 (impurity) is showing zero local magnetic moment.<br>

I hope i made the problem more clear this time.<br>

        <br>

Regards<br>

        </td>

      </tr>

    </tbody>

  </table>

</blockquote>

<br>

<br>

<pre class="moz-signature" cols="72">-- 

Gabriele Sclauzero, EPFL SB ITP CSEA

PH H2 462, Station 3, CH-1015 Lausanne

</pre>

</body>

</html>