We have applied the AFM-CG method to parametrize uniaxial CG models of benzene and perylene in the liquid phase. Force and torque matching of both models gave reasonable errors, with the force error for the AFM-CG benzene model comparable with that obtained previously for a multi-site CG benzene model that considered many-body interactions. As the AFM-CG method uses the pair-additive approximation for the non-bonded potential, the AFM-CG potential was fitted to an analytical form of pair-wise potential to enhance simulation speed. A linear correction was added to the fitted CG potential to match the FG NPT simulation pressure for CG simulations in the NPT ensemble. The fitted AFM-CG benzene model reproduced the FG liquid structures well, especially compared with another recent single-site model of benzene. Reasonable agreement with the FG benzene model was also obtained for system density and dynamical properties. Although the CG model was not optimized for dynamics, the translational and rotational diffusion coefficients were within a factor-of-two of the FG values; however, the relative rate of rotational diffusion of the out-of-plane and in-plane molecular axes differed qualitatively from the FG model. The perylene CG model showed significant improvement over previous single-site CG models, in terms of both structural and thermodynamic properties. Even though there was slightly stronger molecular alignment in the AFM-CG model compared with the FG model, which is in part due to the more symmetric CG uniaxial shape, agreement between the models for the relative probability of orientations of molecular pairs was good. Higher accuracy of the AFM-CG perylene model could be obtained using biaxial CG particles.