<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div>The level occupations are not negative as seen from the eigenvalues. </div><div id="AppleMailSignature">Some of the matrix elements are negative but this only tells the relative sign of the atomic orbitals components in the occupied states<br><br>stefano <div>(sent from my phone)</div></div><div><br>On 22 Mar 2016, at 07:54, Mostafa Youssef <<a href="mailto:myoussef@mit.edu">myoussef@mit.edu</a>> wrote:<br><br></div><blockquote type="cite"><div>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">



<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear all,<br>
<br>
I'm relaxing the structure of rutile TiO2 using ultrasoft pseudopotentials for Ti and O and PBEsol + U with U=5 eV on Ti.  For that purpose I use occupations = 'fixed'.  But I found that one of the Ti ions has negative occupations. To me this looks a bit strange
 since I'm not using any smearing.  The same result can be obtained by another set of PP. Improving the cutoffs or the k-point sampling does not resolve the problem. Any thoughts about the reason and possible remedies are appreciated.  I copy below the occupation
 matrices for the two Ti ions in the cell and the input file.<br>
<br>
<br>
Thank you,<br>
Mostafa Youssef<br>
MIT<br>
<br>
The matrices are:<br>
<br>
atom    1   Tr[ns(na)] =   2.70769<br>
    eigenvalues: <br>
  0.142  0.177  0.182  0.417  0.436<br>
    eigenvectors:<br>
  0.675  0.000  0.000  0.325  0.000<br>
  0.000  0.000  0.500  0.000  0.500<br>
  0.000  0.000  0.500  0.000  0.500<br>
  0.000  1.000  0.000  0.000  0.000<br>
  0.325  0.000  0.000  0.675  0.000<br>
    occupations:<br>
  0.231  0.000  0.000  0.000 -0.129<br>
  0.000  0.309 -0.127  0.000  0.000<br>
  0.000 -0.127  0.309  0.000  0.000<br>
  0.000  0.000  0.000  0.177  0.000<br>
 -0.129  0.000  0.000  0.000  0.328<br>
atom    2   Tr[ns(na)] =   2.70768<br>
    eigenvalues: <br>
  0.142  0.177  0.182  0.417  0.436<br>
    eigenvectors:<br>
  0.675  0.000  0.000  0.325  0.000<br>
  0.000  0.000  0.500  0.000  0.500<br>
  0.000  0.000  0.500  0.000  0.500<br>
  0.000  1.000  0.000  0.000  0.000<br>
  0.325  0.000  0.000  0.675  0.000<br>
    occupations:<br>
  0.231  0.000  0.000  0.000  0.129<br>
  0.000  0.309  0.127  0.000  0.000<br>
  0.000  0.127  0.309  0.000  0.000<br>
  0.000  0.000  0.000  0.177  0.000<br>
  0.129  0.000  0.000  0.000  0.328<br>
N of occupied +U levels =    5.415371<br>
<br>
The input file:<br>
&CONTROL<br>
   calculation =   'vc-relax'   ,<br>
   verbosity = 'high' ,<br>
   restart_mode = 'from_scratch' ,<br>
   outdir='./' ,<br>
   pseudo_dir = './' ,<br>
   tstress = .true.  ,<br>
   tprnfor= .true.  ,<br>
   nstep = 2000 ,<br>
   etot_conv_thr = 1.0D-7  ,<br>
   forc_conv_thr = 4.0D-5 ,<br>
 /<br>
 &SYSTEM<br>
   ibrav = 0 ,<br>
    nat  =  6 , <br>
   ntyp  =  2 , <br>
   nbnd  = 48 ,<br>
  ecutwfc  =  45 ,<br>
  ecutrho  =  360 , <br>
  nosym =.true.<br>
  occupations = 'fixed',<br>
  lda_plus_u = .true.,<br>
  Hubbard_U(1) = 5.0 ,<br>
  Hubbard_U(2) = 0.0 ,<br>
 /<br>
 &ELECTRONS<br>
                    diagonalization='david',<br>
                    mixing_mode = 'plain' ,<br>
                    mixing_beta = 0.7,<br>
                    startingwfc = 'random',<br>
                    conv_thr  =  1.0d-8  , <br>
  /<br>
&ions<br>
    ion_dynamics='bfgs'<br>
/<br>
&CELL<br>
   cell_dynamics = 'bfgs' ,<br>
   cell_dofree = 'shape',<br>
   press_conv_thr = 0.5 ,<br>
 /<br>
CELL_PARAMETERS  angstrom<br>
    4.6320252669841118    0.0000000000000000    0.0000000000000000<br>
    0.0000000000000000    4.6320252669841118    0.0000000000000000<br>
    0.0000000000000000    0.0000000000000000    2.9608930000414442<br>
ATOMIC_SPECIES<br>
  Ti   47.867   ti_pbesol_v1.4.uspp.F.UPF<br>
  O    15.9994  o_pbesol_v1.2.uspp.F.UPF<br>
ATOMIC_POSITIONS crystal <br>
Ti   0.0000000000000000  0.0000000000000000  0.0000000000000000<br>
Ti   0.5000000000000000  0.5000000000000000  0.5000000000000000<br>
O    0.3048741584772117  0.3048741584772117  0.0000000000000000<br>
O    0.6951258415227883  0.6951258415227883  0.0000000000000000<br>
O    0.8048741584772117  0.1951258415227883  0.5000000000000000<br>
O    0.1951258415227883  0.8048741584772117  0.5000000000000000<br>
K_POINTS automatic <br>
  4 4 6  1 1 1<br>
<br>
</div>


</div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>Pw_forum mailing list</span><br><span><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a></span><br><span><a href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a></span></div></blockquote></body></html>