<div>hi everyone:</div>
<div>my question is somewhat old ,i have seen a lot of related discussion in this forum .but my results does not like previous instance.</div>
<div>i only want to relax c of crystal lattice ,and fixed c,computes it,and then change it  ,compute it,and .....</div>
<div>i get the  E VS C tabular </div>
<div>so strange is that :</div>
<div> </div>
<div>c               E<br clear="all">2.655     -137.78747<br>2.725     -137.7912<br>2.76      -137.79219<br>2.77       -137.79233</div>
<div>2.78      -137.79248<br>2.795      -137.79266<br>2.8         -137.7927<br>2.81       -137.79275<br>2.82       -137.79275<br>2.83       -137.79278<br>2.84       -137.79276<br>2.85       -137.79279<br>(1)the energy is down  before 2,81,and c=
2.83 is like local minimum ,and then the energy is go on to down. the graph is smooth ,but   i can not get the parabola shape .so the best relaxed structure is unkown to me.</div>
<div>  and  i saw many discuss about the down energy ,maybe the  ecutrho issue.  but i am not sure .if ti is ,i should test the ecutrho?i post my input file here:</div>
<div>&CONTROL<br>calculation  = "relax",<br>  prefix       = "xyh",<br>  pseudo_dir   = "/home/xyh/intel_espresso-3.2/pseudo",<br>  outdir       = "/raid/xyh/pwscf/tmp/",<br>  etot_conv_thr=
1.D-5<br>  forc_conv_thr=1.D-4<br>  nstep=50<br>  tstress=.TRUE.<br>/<br>&SYSTEM<br>  ibrav     = 4,<br>  a=20,b=20,c=2.7950,cosab=-0.5,cosac=0,cosbc=0,<br>  nat       = 12,<br>  ntyp      = 2,<br>  ecutwfc   =29.9876
 ,<br>  ecutrho   =200<br>/<br>&ELECTRONS<br>  conv_thr    = 1.D-9,<br>  mixing_beta = 0.2D0,<br>/<br>&IONS<br>  pot_extrapolation = "second_order",<br>  wfc_extrapolation = "second_order",<br>
/<br> ATOMIC_SPECIES<br>H  1.0  H.pbe-van_bm.UPF<br>O 15.999 O.pbe-van_bm.UPF<br>ATOMIC_POSITIONS { crystal }<br>H  0.0373836515057430  0.9503896146857272  0.0102606085503574<br>H  0.1115317646909714  0.9930799089192295  
0.3344160938275307<br>H  0.0789318950169224  0.0718359337861800  0.9884868039481935<br>H  0.0724243415282819  0.1328116898677112  0.6645236821235255<br>H  0.9626484865634389  0.0496365043168296  0.0101367101334339<br>H  0.8884876867936907
  0.0069387333102180  0.3342097796899214<br>H  0.9210924986163674  0.9281832378776828  0.9882981504770443<br>H  0.9276176392628244  0.8672101853812219  0.6643346513014633<br>O  0.0926597580972745  0.9929556282114975  </div>

<div>......</div>
<div> </div>
<div> </div>
<div>..</div>
<div>(2)and another question is although i can not get the lowest erengy,but i can get the parabola of  the stress. maybe someone the stress is not consistent with energy,for i only want to    know  which  structure is best ,so the exactly stress is not so import for me, i only want to know the trend .in this point,whether i can rely on the stress to determine the best structure.
</div>
<div> </div>
<div>i need your help.</div>
<div>thanks <br> -- <br>Xu Yuehua<br>physics Department of Nanjing university<br>China </div>