I wonder if anyone knows whether the thermostatted SLLOD equations for shear flow are implemented with the velocity Verlet algorithm or whether it acn be used with a multiple time-step algorithm? According the manual it is implemented with a velocity verlet algorithm but the formulation is a bit cryptic.
"As implemented in LAMMPS, they are coupled to a Nose/Hoover chain
thermostat in a velocity Verlet formulation, closely following the
implementation used for the fix nvt command."