Tripathy, Agarwal, and Kumar 34 used the total instantaneous force between pairs of molecules extracted from condensed-phase all-atom simulations to fit the parameters in an anisotropic CG pair potential for several different polyaromatic hydrocarbons. Tanis et al. used a similar force-matching approach, but instead used average pair forces and torques from condensed-phase all-atom simulations as a function of the pair configuration to fit analytical CG pair potentials. Although demonstrated to accurately reproduce properties of the FG all-atom models for specific systems, these methods do not rigorously account for many-body effects on the CG PMF, even though the molecule pairs used in the parametrization were extracted from condensed-phase simulations. They also rely on CG potentials that can be decomposed into pair interactions, limiting their ability to be generalized to more complex effective interactions.