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

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