Hello, I see that in the first step the thermostat is not used. This is a problem when one wants to equilibrate a system where molecules overlap. Using a Berendsen thermostat with a time constant of timestep/10 would allow to remove overlaps for these systems. Thanks, Ramon