where R c is the cutoff distance for the non-bonded CG interactions and b is a parameter that was chosen so that the pressure in the CG simulation matched that of the corresponding FG simulation under the parametrization conditions. This pressure correction was relatively small compared with the fitted modified S-function potential for the systems studied, as shown in the SM in Figures and. More rigorous methods of pressure-matching for CG parametrization in the NPT ensemble involving thermodynamic state-dependent corrections are possible, along the lines of those previously applied to isotropic CG models, and will be the focus of future study.