<div dir="ltr"><div>Dear all,</div><div><br></div><div>We are running NPT simulations 
on an 80 atom cell. The simulations have been set up with 'nose' 
thermostat for ions and the parrinello rahman 'pr' barostat. The 
simulations were running fine, but there was a slow increase in the 
fictitious kinetic energy 'ekincw' of the electrons. Therefore, the 
electronic thermostat was used to control this by setting 
`electron_temperature = 'nose'`. But when this was done, the constant of
 motion 'econt' which was constant started increasing. The econt values 
which had been constant (changes only occurring in the range of 10e-4 
a.u) until the application of electronic thermostat started increasing 
at a rate greater than 0.05 a.u per ps.</div><div><br></div><div>Is the 
constant of motion not a constant for electronic thermostat? or does 
this imply a discrepancy in the calculation? If so, what would cause 
this behaviour?</div><div><br></div><div>Any information regarding this will be helpful. The relevant inputs are given below:</div><div><br></div><div>---------------------------------------------------</div><div>&electrons<br>    electron_dynamics = 'verlet',<br>    emass = 200.d0,<br>    emass_cutoff = 1.d0,<br>    electron_velocities = 'default',<br>    electron_maxstep = 2000,<br>    conv_thr = 1.0D-10,<br>    electron_temperature = 'nose',<br>    ekincw = 0.0125<br>    fnosee = 150<br>/<br>&ions<br>    ion_dynamics = 'verlet',<br>    ion_temperature = 'nose',<br>    tempw= 300.0,<br>    fnosep = 55.0,<br>    ion_velocities = 'default',<br>    ion_positions = 'default',<br>/<br>&cell<br>    cell_dynamics = 'pr',<br>    cell_parameters = 'default'<br>    cell_velocities = 'default'<br>    press = 300.0,<br>/</div><div>--------------------------------------------------------------</div><div><br></div><div>Rashid Rafeek</div><div>Junior Research Fellow</div><div>IISER Bhopal, India</div></div>