with interactions truncated at a system-dependent cut-off distance R c (given in the sections below). After optimizing the pressure-correction potential (Eq. ( )) in the NVT ensemble at the same temperature used for the AFM-CG parametrization, CG simulations were undertaken in the NPT ensemble. CG simulations in the NVT ensemble with densities set to the FG average values were also conducted to compare with the FG simulations without the confounding effect of differing densities. Temperature and pressure were controlled with a Nosé-Hoover thermostat and barostat and particle positions were updated using the velocity Verlet 66 time integrator in both the FG and CG simulations. The equations of motion in the CG simulations were integrated using the ASPHERE package 67 in LAMMPS, in which particle orientations are stored as quaternions and updated using Richardson iterations.