[Pw_forum] temperature control in PWscf

Eric Schwegler schwegler at llnl.gov
Thu Jan 12 18:23:53 CET 2006


Alain,
I've run into this as well. It seems to me that in version 3.0, 
dynamics.f90 does not save the
velocities between MD timesteps. So after the initial velocity scaling 
they are basically set
to zero. Simply adding the vel array to the unit=4 read and write 
statements in dynamics.f90
seems to fix the problem for me. Is this a known bug?
-Eric Schwegler



> Dear all,
> I am trying to perform a temperature rescaling MD with PWscf
>
> &IONS
>     ion_dynamics = 'verlet'
>     ion_temperature = 'rescaling'
>     tempw = 500.0
>     nraise = 1
>     tolp = 1.d-8
>     wfc_extrapolation = 'first_order'
> /
>
> The starting temperature is correct but it decreases very rapidly
>
> 358:     Starting temperature  =   500.00 K
> 405:     temperature           =   496.70493449 K
> 693:     temperature           =     0.24166052 K
> 903:     temperature           =     0.21332563 K
> 1113:     temperature           =     0.18863928 K
> 1323:     temperature           =     0.16711531 K
> 1533:     temperature           =     0.14835294 K
> 1743:     temperature           =     0.13194319 K
> 1953:     temperature           =     0.11759977 K




More information about the Pw_forum mailing list