<div dir="ltr"><div>Dear users and developers,</div><div><br></div><div>Are there any tricks for running an NPT ensemble with pw.x (vc-md) for high-temperatures? I do not have problems with 300-350-400K (reaching minimum Etot+Ekin fluctuations after a reasonable amount of steps) but cannot reach it for 600K. Below I show the basic input parameters used.</div><div><br></div><div>Thank you for any suggestions!</div><div><br></div><div><br></div><div><br></div><div>&CONTROL<br>  restart_mode = 'restart'<br>  calculation = 'vc-md'<br>  nstep = 1000<br>  dt = 60   ! 3fs, 1/10 of the highest vibration T<br>!  verbosity = 'high'<br>  wf_collect = .true.<br>  pseudo_dir = '.'<br>/<br>&SYSTEM</div><div>[all good]<br>  nosym = .True.<br>/<br>&ELECTRONS<br></div><div>[all defaut]<br></div><div>/<br>&IONS<br>  ion_velocities = 'default'<br>  tempw = 600<br>  ion_dynamics = 'beeman'<br>  pot_extrapolation = 'atomic'   ! 1st, 2nd order from prev step, no gain<br>  wfc_extrapolation = 'none'   ! 1st, 2nd order from prev step, no gain<br>! remove_rigid_rot = .True. trick for not big enough supercells<br>  ion_temperature = 'rescaling'<br>  tolp = 10<br>  nraise = 1<br>  refold_pos = .True.<br>/<br>&CELL<br>  cell_dynamics = 'pr'   ! Parrinello-Rahman, also Wentzcovitch 'w' is avail<br>!  wmass - fictitious cell mass, see default for 'pr' or 'w'<br>  cell_factor = 2<br>/<br><br>ATOMIC_SPECIES<br>...<br><br>ATOMIC_POSITIONS crystal<br>...<br><br>K_POINTS automatic<br>[converged grid tested to give accurate cell stresses]</div><div><br></div><div><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Best wishes,</div><div>Al.<span style="border-collapse:separate;line-height:normal;text-align:-webkit-auto;text-indent:0px;border-spacing:0px"><span style="text-align:-webkit-auto">,</span></span> phd candidate in chemical sciences<br>'I.. a universe of atoms, an atom in the universe' (Richard P. Feynman)<br><a href="https://cpms.kaust.edu.sa/" target="_blank">https://cpms.kaust.edu.sa/</a></div></div></div></div></div>

<br>
<div><hr></div><font face="Arial" size="1">This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.</font>