Development of an anisotropic polarizable model for the all-atom AMOEBA force field

For planar and rigid π-conjugated molecular systems, electrostatic and inductive interactions are pivotal in governing molecular packing structures and electron polarization energies. These electrostatic interactions typically exhibit an anisotropic nature within π-conjugated systems. In this study, we utilize the atoms in molecules (AIM) theory in conjunction with linear response theory to decompose molecular polarizability into distributed atomic polarizability tensors. On the basis of atomic polarizability tensors, we extended an anisotropic polarizable model into the AMOEBA polarizable force field. Both anisotropic and isotropic polarizable models in combination with various density functional theory (DFT)-derived atomic multipoles were applied to optimize the experimental crystals of naphthalene and anthracene. Furthermore, these two types of electrostatic models, coupled with the evolutionary algorithm USPEX program, are utilized to predict the crystal structures of oligoacenes. Our findings demonstrate that the anisotropic polarizable model exhibits superior performance in crystal refinement and crystal structure prediction. This enriched anisotropic polarizable model is seamlessly integrated into the AMOEBA polarizable force field and readily applicable within our modified Tinker program.