<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <div class="moz-cite-prefix">are you sure the code is not correct?<br>
      how are the oxygen octahedra oriented  in your case?<br>
      a rotation of pi/4 around z would do..<br>
      stefano<br>
      <br>
      On 09/09/2012 04:35 AM, Peng Chen wrote:<br>
    </div>
    <blockquote
cite="mid:CAP_5U0d+Kqbf-JOfpJw_QN-i3eFH1T_FPdNVsFfaianPcyW-+Q@mail.gmail.com"
      type="cite">
      <pre wrap="">Hi All,

When I tried to do GGA/GGA+U computation on Ni3V2O8, the occupation number
of Ni looks a little weird. In this compound,  d orbitals of Ni  split into
t2g (dxy, dxz,dyz) and eg (dx^2-y^2, dz^2) because of NiO6 octahedral
coordination. So the down spin will occupy t2g states. But both GGA and
GGA+U results show the occupation of dxz, dyz and dx^2-y^2. I tried to
write occupation matrix using  starting_ns_eigenvalue so that the
occupation can be dxy, dxz and dyz, but the scf calculation can not reach
convergence. Any suggestion is appreciated!
 

GGA after the first iteration:
    eigenvalues:
  0.059  0.066  0.075  0.146  0.161
    eigenvectors:
  0.009  0.001  0.000  0.000  0.990
  0.000  0.000  0.997  0.003  0.000
  0.076  0.924  0.000  0.000  0.000
  0.916  0.075  0.000  0.000  0.009
  0.000  0.000  0.003  0.997  0.000
    occupations:
  0.160  0.001  0.000 -0.010  0.000
  0.001  0.075  0.000  0.000 -0.004
  0.000  0.000  0.066 -0.002  0.001
 -0.010  0.000 -0.002  0.060  0.000
  0.000 -0.004  0.001  0.000  0.146

GGA after scf convergence:
    eigenvalues:
  0.336  0.364  0.987  0.988  0.991
    eigenvectors:
  0.030  0.968  0.002  0.000  0.000
  0.001  0.000  0.026  0.822  0.151
  0.000  0.000  0.063  0.175  0.761
  0.000  0.001  0.909  0.002  0.088
  0.968  0.030  0.000  0.001  0.000
    occupations:
  0.364 -0.003 -0.014  0.022  0.005
 -0.003  0.988 -0.001  0.000 -0.024
 -0.014 -0.001  0.990  0.001  0.001
  0.022  0.000  0.001  0.987 -0.001
  0.005 -0.024  0.001 -0.001  0.338


GGA+U after the first iteration:
   eigenvalues:
  0.071  0.079  0.095  0.157  0.172
    eigenvectors:
  0.008  0.000  0.000  0.000  0.992
  0.024  0.363  0.607  0.006  0.000
  0.000  0.630  0.369  0.001  0.000
  0.967  0.006  0.018  0.000  0.008
  0.000  0.000  0.006  0.993  0.000
    occupations:
  0.171  0.001  0.001 -0.009  0.000
  0.001  0.089 -0.008 -0.003 -0.005
  0.001 -0.008  0.085  0.001  0.002
 -0.009 -0.003  0.001  0.072  0.000
  0.000 -0.005  0.002  0.000  0.157

GGA+U after scf convergence:
    eigenvalues:
  0.265  0.306  0.995  0.996  0.996
    eigenvectors:
  0.010  0.988  0.001  0.000  0.000
  0.001  0.000  0.003  0.919  0.077
  0.000  0.000  0.098  0.060  0.841
  0.000  0.001  0.897  0.020  0.082
  0.989  0.010  0.000  0.001  0.000
    occupations:
  0.307 -0.002 -0.014  0.023  0.004
 -0.002  0.995  0.000  0.000 -0.024
 -0.014  0.000  0.996  0.001  0.000
  0.023  0.000  0.001  0.994 -0.001
  0.004 -0.024  0.000 -0.001  0.266


</pre>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
    </blockquote>
    <br>
  </body>
</html>