Abstract
We present an ab initio polarizable representation of classical molecular mechanics (MM) atoms by employing an angular momentum-based expansion scheme of the point charges into partial wave orbitals. The charge density represented by these orbitals can be fully polarized, and for hybrid quantum-mechanical-molecular-mechanical (QM/MM) calculations, mutual polarization within the QM/MM Hamiltonian can be obtained. We present the mathematical formulation and the analytical expressions for the energy and forces pertaining to the method. We further develop a variational scheme to appropriately determine the expansion coefficients and then validate the method by considering polarizations of ions by the QM system employing the hybrid GROMACS-CPMD QM/MM program. Finally, we present a simpler prescription for adding isotropic polarizability to MM atoms in a QM/MM simulation. Employing this simpler scheme, we present QM/MM energy minimization results for the classic case of a water dimer and a hydrogen sulfide dimer. Also, we present single-point QM/MM results with and without the polarization to study the change in the ionization potential of tetrahydrobiopterin (BH(4)) in water and the change in the interaction energy of solvated BH(4) (described by MM) with the P(450) heme described by QM. The model can be employed for the development of an extensive classical polarizable force-field.