<div dir="ltr"><div><div><div><div><div>Dear QE users<br></div>I am trying to model amorphous Si using heat-and-quench by BO molecular dynamics (based on example3). My starting temperature is 1600K. But, In my output file, after every step, I see a different instantaneous temperature printed. Why is that so? Am I missing any input parameter? My input file is as given below:<br>

<br> &control<br>    calculation='md'<br>    restart_mode='from_scratch',<br>    prefix = 'Si'<br>    pseudo_dir = '/home/deva/espresso-5.0.1/pseudo/',<br>    outdir='./',<br>
    dt=20,<br>    nstep=60,<br>    disk_io='high'<br> /<br> &system<br>    ibrav= 1, celldm(1)=10.18, nat=  8, ntyp= 1,<br>    ecutwfc = 20.0, nosym=.true.<br> /<br> &electrons<br>    conv_thr =  1.0d-8,<br>

    mixing_beta = 0.7<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>ATOMIC_SPECIES<br> Si  28.086  Si.pz-vbc.UPF<br>ATOMIC_POSITIONS<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></div>Thanks in advance<br></div>Devan<br></div><div>PhD student<br>
</div>UMD<br></div>College Park<br></div>