<div>Dear all,</div><div> </div><div>I am running vc-md of a cuprite structure with 2 by 2 by 2 supercell. When I set the tempw=300, it does converge but the equilibrium temperature seems to be only 190K. It starts from 280 and all the way down to 190K and converges there. Actually I did not find many controls on it, below is what my current setups for &ions and &cells.</div>
<div> </div><div>I kept trying different setups, but has little effect, for example, I tried different 'dt' from 20 to 50, different 'ion_temperature' choices, 'tolp' or 'nraise'. I also used a primitive cell, the actual temperature can be only 108K or so. Therefore, I am not sure whether it has something to do with the cell size, or the temperature calibration issue, or some setup would substantially influence the temperature control.</div>
<div> </div><div> &ions<br> tempw=300,<br> ion_temperature= 'rescaling',<br> tolp=20,<br> nraise=1,<br> pot_extrapolation='second-order',<br> wfc_extrapolation='second-order',<br>
/<br> &cell<br> cell_dynamics = 'pr' ,<br> press = 0.00<br> /<br clear="all"></div><div>Besides of it, I have a stupid question because I did not find any information about the atom velocities. I only find updated coordinates and forces in .out file. In .save folders, I did not find anything either. My control icons for print are as follows:</div>
<div> </div><div> disk_io='high'<br>! wf_collect=.TRUE.<br> pseudo_dir = '.',<br> outdir='\temp'<br> tstress =.true.,<br> tprnfor =.true.,<br> </div><div>Thank you in advance if you have any idea about these.</div>
<div>Best Regards,</div><div>Tian<br>-- <br>Lan, Tian<br>
Ph.D. Candidate, Department of Applied Physics and Materials Science<br>
California Institute of Technology,<br>
Caltech M/C 138-78, Pasadena, CA, 91125<br>
</div>