I am currently trying to implement a new respa integrator, which is
region based. i.e. I defined three regions integrated at different
time steps. The work flow seems to be correct.
The temperature remains however very low (~7K), despite a Langavin
thermostat set at 300 K. Does anyone know where/when/how is the
temperature computed in lammps? And which parameter could affect this