where ρ is the number density of molecules. The ARDFs at 300 K of the FG and AFM-CG benzene models are compared in Figure , along with that of a recent single-site benzene CG model by Bowen et al. , which we will refer to as the "Bowen CG model". The Bowen CG model was parametrized from the same FG OPLS-AA force field but by matching the interaction energy between pairs of rigid benzene molecules in a variety of configurations, and so its parametrization does not account for finite-temperature many-body effects that are taken into account in the AFM-CG method. Figure shows that the ARDF for the Bowen CG model is relatively isotropic, except at around 5 Å where the edge-to-face configuration is most favorable, whereas the distributions vary anisotropically for the FG and AFM-CG models over a wider range of intermolecular separations. The ARDFs for the FG and AFM-CG models are in good agreement, except for some slight differences at a distance of around 4.5 Å and 6.5 Å. These deviations can be explained by minor differences between the fitted S-function potential compared with the AFM-CG basis expansion potential shown in Figure and in Figure in the SM. in the SM. The diffusion coefficients calculated from the fitted slopes of these plots are plotted as a function of temperature along with experimental data in Figure . By linear interpolation and extrapolation of the line-of-best-fit to the experimental points, the FG diffusion coefficients are found to overestimate the experimental values by approximately 31%, 21%, 7% and 10% at 280, 300, 320 and 330 K, respectively. Finite-size effects lower the translational diffusion coefficient in simulated systems with periodic boundaries relative to the infinite system-size limit, and so cannot explain the discrepancy between the FG model and experiment. The most likely explanation for the overestimated FG diffusion coefficients is the underestimated densities shown in Figure . The AFM-CG model shows faster translational diffusion than the FG model at all temperatures studied, with the highest relative error of 32% at 280 K. Faster dynamics compared with the underlying FG model is a general feature of CG models with only conservative forces and can be explained by the coarse-graining process integrating out degrees of freedom of the FG system that act as sources of dissipation.