A velocity Verlet integrator 40 with a timestep of 1 fs for the first two equilibration runs and 2 fs for the rest of the runs was used. In order to avoid membrane and water systems from moving laterally and possibly drifting in opposite directions, the motions of the center of mass of the bilayer and the solvent were reset separately through the comm-grps keyword.