NPT CG simulations of 1000 CG perylene molecules were carried out with the optimized modified S-function potential at 1 atm pressure and the same temperatures as the FG simulations using a time step of 10 fs. CG simulations were also carried out of two existing biaxial ellipsoidal models of perylene developed by Babadi, Everaers, and Ejtehadi 25 and by Berardi, Fava, and Zannoni 86 under the same conditions and compared to the FG and AFM-CG models (we will refer to these models as the "Babadi CG model" and the "Berardi CG model"). These CG models of perylene were parametrized from atomistic FG models by matching the interaction energy between pairs of rigid perylene molecules in a variety of configurations, which neglects finite-temperature many-body effects that are taken into account in the AFM-CG method. It should be noted that the Babadi and Berardi CG models were parametrized using different FG models (MM3 87 and the Universal Force Field (UFF) , respectively) from the OPLS-AA FG model used here. So the comparison with the FG and AFM-CG models in this work is not completely equitable and for this reason we also compare with experiment below. For temperatures 550-670 K, which are approximately at or above the experimental melting point of perylene of 554 K, 89 NPT simulations of the models showed that both the FG and AFM-CG models were in the liquid phase, whereas the Babadi and Berardi CG models rapidly crystallized and remained solid at all temperatures, as illustrated in the simulation snapshots of the models at T = 570 K in Figure ,