ated every innermost sub-timestep
bonded 3- and 4-body forces are computed every 1st parameter sub-timesteps
short-range nonbond pairwise forces (LJ, Coulombic) are computed every 
  (2nd parameter * 1st parameter) sub-timesteps
long-range (Ewald, PPPM) forces are computed every 
  (3rd parameter * 2nd parameter * 1st parameter) sub-timesteps
the timestepping for all 3 inner loops (bond, 3/4-body, nonbond) is performed
  as sub-cycling within the long-range timestepping loop
the fastest (innermost) timestep