<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.23487"></HEAD>
<BODY>
<DIV dir=ltr align=left><SPAN class=296151018-04062013><FONT color=#0000ff 
size=2 face=Arial>Hi Franco,</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=296151018-04062013><FONT color=#0000ff 
size=2 face=Arial>A few questions. Is it an antiferro system? Why do you think 
you need spin polarization for relaxation? Why do you shift the K points?  
Is it electronic structure or the atomic positions that are not converged? Are 
there any signs of convergence from the output?</FONT></SPAN></DIV>
<DIV dir=ltr align=left><SPAN class=296151018-04062013><FONT color=#0000ff 
size=2 face=Arial>fzm</FONT></SPAN></DIV><BR>
<DIV dir=ltr lang=fr class=OutlookMessageHeader align=left>
<HR tabIndex=-1>
<FONT size=2 face=Tahoma><B>De :</B> pw_forum-bounces@pwscf.org 
[mailto:pw_forum-bounces@pwscf.org] <B>De la part de</B> Franco 
Bonafé<BR><B>Envoyé :</B> 4 juin 2013 13:45<BR><B>À :</B> PWSCF 
Forum<BR><B>Objet :</B> [Pw_forum] SCF doesn't converge in charged 
supercell<BR></FONT><BR></DIV>
<DIV></DIV>
<DIV dir=ltr>
<DIV>Hi to all,</DIV>
<DIV>I am trying to perform a relaxation in a charged cell and it reaches the 
maximun number of steps and doesn't achieve convergence.</DIV>
<DIV>What am I doing wrong?</DIV>
<DIV>Thanks in advance.</DIV>
<DIV>Franco.</DIV>
<DIV><BR></DIV>
<DIV>Input script:</DIV>
<DIV><BR></DIV>
<DIV> &CONTROL</DIV>
<DIV>                 calculation = 
'relax' ,</DIV>
<DIV>                restart_mode = 
'from_scratch' ,</DIV>
<DIV>                    
  outdir = '/home/tmp' ,</DIV>
<DIV>                  pseudo_dir = 
'/home/pseudo/' ,</DIV>
<DIV>                    
  prefix = 'anath2oLi ,</DIV>
<DIV>               etot_conv_thr = 
1.0D-3 ,</DIV>
<DIV>               forc_conv_thr = 
1.0D-2 ,</DIV>
<DIV> /</DIV>
<DIV> &SYSTEM</DIV>
<DIV>                    
   ibrav = 0,</DIV>
<DIV>                  
 celldm(1) = 19.3262367 ,</DIV>
<DIV>                    
     nat = 73 ,</DIV>
<DIV>                    
    ntyp = 4 ,</DIV>
<DIV>                    
 ecutwfc = 30. ,</DIV>
<DIV>                    
 ecutrho = 300. ,</DIV>
<DIV>                  tot_charge = 
1.000000,</DIV>
<DIV>                 occupations = 
'smearing' ,</DIV>
<DIV>                    
 degauss = 0.03 ,</DIV>
<DIV>                    
smearing = 'methfessel-paxton' ,</DIV>
<DIV>                    
   nspin = 2 ,</DIV>
<DIV>           tot_magnetization = 0,</DIV>
<DIV> /</DIV>
<DIV> &ELECTRONS</DIV>
<DIV>                    
conv_thr = 1.0d-6 ,</DIV>
<DIV>                 mixing_mode = 
'local-TF' ,</DIV>
<DIV>                 mixing_beta = 
0.7 ,</DIV>
<DIV>             diagonalization = 'david' 
,</DIV>
<DIV> /</DIV>
<DIV> &IONS</DIV>
<DIV>                ion_dynamics = 
'bfgs' ,</DIV>
<DIV> /</DIV>
<DIV>CELL_PARAMETERS alat </DIV>
<DIV>     1.010387266    0.000000000   
 0.000000000 </DIV>
<DIV>     0.000000000    0.709834502   
 0.000000000 </DIV>
<DIV>     0.000000000    0.000000000   
 2.151168486 </DIV>
<DIV>ATOMIC_SPECIES</DIV>
<DIV>    O   15.99940  O.pbe-van_ak.UPF </DIV>
<DIV>   Ti   47.86700  Ti.pbe-sp-van_ak.UPF </DIV>
<DIV>    H    1.00794  H.pbe-van_ak.UPF </DIV>
<DIV>   Li    6.94100  Li.pbe-s-van_ak.UPF </DIV>
<DIV>ATOMIC_POSITIONS angstrom </DIV>
<DIV>...</DIV>
<DIV>K_POINTS automatic </DIV>
<DIV>  1 1 1   1 1 1 </DIV><BR clear=all>
<DIV><BR></DIV>-- <BR>
<DIV dir=ltr><SPAN style="BORDER-COLLAPSE: collapse"><FONT 
face="georgia, serif">
<DIV style="FONT-SIZE: 13px"><SPAN 
style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)">Franco P. 
Bonafé</SPAN></DIV>
<DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><FONT 
size=1><I>Chemistry Student</I></FONT></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><I 
style="FONT-SIZE: x-small">Undergrad Research Assistant</I></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><I><FONT 
size=1>Department of Math & Physics</FONT></I></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><I><FONT 
size=1>National University of Córdoba</FONT></I></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><FONT 
size=1><I>Argentina</I></FONT></DIV></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><FONT 
size=1><I><BR></I></FONT></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><FONT 
size=1><I>+57 9 0351 15 5472791</I></FONT></DIV>
<DIV style="BACKGROUND-COLOR: rgb(255,255,255); COLOR: rgb(34,34,34)"><FONT 
size=1><I><A href="mailto:fbonafe@fcq.unc.edu.ar" 
target=_blank>fbonafe@fcq.unc.edu.ar</A></I></FONT></DIV></FONT></SPAN></DIV></DIV></BODY></HTML>