<div dir="ltr"><div>Dear all<br></div><div>The program mentioned below is for heating and quenching of cooling, i have some question regarding this, please tell me what mistake am i doing.<br><br></div><div>1This system doesn't quench and show higher temperature at the end,say initial temp is 1600K and, End temperature is around 6000K, when is using 'reduced-T' ,and around 2500K when berendsen is used, except these two, if i use any different way i.e rescaling etc. then system would not be converge, <b>why?<br>
<br></b></div><div>2.If i choose deltaT and nraise less than 400 and 20 respectively then system doesn't converge ????<br><br><br></div><div>Thanks in advance<a name="id3121444"></a><a name="ion_temperature"><pre><br>
</pre><pre><br></pre></a></div>&control<br><div><div><div>calculation='md'<br>restart_mode='from_scratch',<br>prefix = 'Si'<br>pseudo_dir = '/home/metal/Desktop/espresso-5.0.2/pseudo',<br>
outdir='./',<br>dt=10,<br>nstep=500,<br>disk_io='high'<br>/<br>&system<br>ibrav= 1, celldm(1)=10.18, nat= 8, ntyp= 1,<br>ecutwfc = 35.0, nosym=.true.<br>/<br>&electrons<br>conv_thr = 1.0d-8,<br>mixing_beta = 0.7<br>
<br>/<br>&ions<br>ion_dynamics='verlet'<br>pot_extrapolation='second-order'<br>wfc_extrapolation='second-order'<br>ion_temperature='reduce-T'<br>tempw=1600.D0<br>delta_t=-400.D0<br>nraise=20<br>
<br>/<br>ATOMIC_SPECIES<br>Si 28.086 Si.pz-vbc.UPF<br><br>Si -0.123 -0.123 -0.123<br>Si 0.377 0.377 -0.123<br>Si 0.377 -0.123 0.377<br>Si -0.123 0.377 0.377<br>Si 0.123 0.123 0.123<br>Si 0.623 0.623 0.123<br>Si 0.623 0.123 0.623<br>
Si 0.123 0.623 0.623<br>K_POINTS {automatic}<br>1 1 1 0 0 0<br><br><br><br clear="all"><div><div><div dir="ltr"><div>Anik shrivastava</div>
<div>Senior Research Fellow </div>
<div>Naval Materials Research Lab,DRDO</div>
<div>Mumbai-400085</div>
<div>India</div></div></div>
</div></div></div></div></div>