An Estimation of Hybrid Quantum Mechanical Molecular Mechanical Polarization Energies for Small Molecules Using Polarizable ForceField Approaches
Abstract
Here in this work, we report two polarizable molecular mechanics (polMM) force field models for estimating the polarization energy in hybrid quantum mechanical molecular mechanical (QM/MM) calculations. These two models, named the potential of atomic charges (PAC) and potential of atomic dipoles (PAD), are formulated from the ab initio quantum mechanical (QM) response kernels for the prediction of the QM density response to an external molecular mechanical (MM) environment (as described by external point charges). The PAC model is similar to fluctuating charge (FQ) models because the energy depends on external electrostatic potential values at QM atomic sites; the PAD energy depends on external electrostatic field values at QM atomic sites, resembling induced dipole (ID) models. To demonstrate their uses, we apply the PAC and PAD models to 12 small molecules, which are solvated by TIP3P water. The PAC model reproduces the QM/MM polarization energy with a R ^{2} value of 0.71 for aniline (in 10,000 TIP3P water configurations) and 0.87 or higher for other eleven solute molecules, while the PAD model has a much better performance with R ^{2} values of 0.98 or higher. The PAC model reproduces reference QM/MM hydration free energies for 12 solute molecules with amore »
 Authors:
 Publication Date:
 Research Org.:
 Brookhaven National Laboratory (BNL), Upton, NY (United States)
 Sponsoring Org.:
 USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC22)
 OSTI Identifier:
 1342645
 Report Number(s):
 BNL1134602017JA
Journal ID: ISSN 15499618; R&D Project: 16068; KC0403020
 Grant/Contract Number:
 SC0012704; SC0011297; GM09667802; GM072558; GM051501
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Journal of Chemical Theory and Computation
 Additional Journal Information:
 Journal Volume: 13; Journal Issue: 2; Journal ID: ISSN 15499618
 Publisher:
 American Chemical Society
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
Citation Formats
Huang, Jing, Mei, Ye, König, Gerhard, Simmonett, Andrew C., Pickard, Frank C., Wu, Qin, Wang, LeePing, MacKerell, Alexander D., Brooks, Bernard R., and Shao, Yihan. An Estimation of Hybrid Quantum Mechanical Molecular Mechanical Polarization Energies for Small Molecules Using Polarizable ForceField Approaches. United States: N. p., 2017.
Web. doi:10.1021/acs.jctc.6b01125.
Huang, Jing, Mei, Ye, König, Gerhard, Simmonett, Andrew C., Pickard, Frank C., Wu, Qin, Wang, LeePing, MacKerell, Alexander D., Brooks, Bernard R., & Shao, Yihan. An Estimation of Hybrid Quantum Mechanical Molecular Mechanical Polarization Energies for Small Molecules Using Polarizable ForceField Approaches. United States. doi:10.1021/acs.jctc.6b01125.
Huang, Jing, Mei, Ye, König, Gerhard, Simmonett, Andrew C., Pickard, Frank C., Wu, Qin, Wang, LeePing, MacKerell, Alexander D., Brooks, Bernard R., and Shao, Yihan. Tue .
"An Estimation of Hybrid Quantum Mechanical Molecular Mechanical Polarization Energies for Small Molecules Using Polarizable ForceField Approaches". United States.
doi:10.1021/acs.jctc.6b01125. https://www.osti.gov/servlets/purl/1342645.
@article{osti_1342645,
title = {An Estimation of Hybrid Quantum Mechanical Molecular Mechanical Polarization Energies for Small Molecules Using Polarizable ForceField Approaches},
author = {Huang, Jing and Mei, Ye and König, Gerhard and Simmonett, Andrew C. and Pickard, Frank C. and Wu, Qin and Wang, LeePing and MacKerell, Alexander D. and Brooks, Bernard R. and Shao, Yihan},
abstractNote = {Here in this work, we report two polarizable molecular mechanics (polMM) force field models for estimating the polarization energy in hybrid quantum mechanical molecular mechanical (QM/MM) calculations. These two models, named the potential of atomic charges (PAC) and potential of atomic dipoles (PAD), are formulated from the ab initio quantum mechanical (QM) response kernels for the prediction of the QM density response to an external molecular mechanical (MM) environment (as described by external point charges). The PAC model is similar to fluctuating charge (FQ) models because the energy depends on external electrostatic potential values at QM atomic sites; the PAD energy depends on external electrostatic field values at QM atomic sites, resembling induced dipole (ID) models. To demonstrate their uses, we apply the PAC and PAD models to 12 small molecules, which are solvated by TIP3P water. The PAC model reproduces the QM/MM polarization energy with a R2 value of 0.71 for aniline (in 10,000 TIP3P water configurations) and 0.87 or higher for other eleven solute molecules, while the PAD model has a much better performance with R2 values of 0.98 or higher. The PAC model reproduces reference QM/MM hydration free energies for 12 solute molecules with a RMSD of 0.59 kcal/mol. The PAD model is even more accurate, with a much smaller RMSD of 0.12 kcal/mol, with respect to the reference. Lastly, this suggests that polarization effects, including both local charge distortion and intramolecular charge transfer, can be well captured by induced dipole type models with proper parametrization.},
doi = {10.1021/acs.jctc.6b01125},
journal = {Journal of Chemical Theory and Computation},
number = 2,
volume = 13,
place = {United States},
year = {Tue Jan 24 00:00:00 EST 2017},
month = {Tue Jan 24 00:00:00 EST 2017}
}
Web of Science

We present a hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics study of 18crown6 (18c6) in both polar hydrogen bonding and nonpolar solvents. The QM/MM method described here is based on our previous QM/MM study of K{sup +}/18c6 which employed the semiempirical AM1 method for 18c6 and the SPC/e model for H{sub 2}O. We improve our previous QM/MM method by including MM torsion terms to better describe the energetics of OCCO and COCC rotamers. The torsion terms are parametrized against the ab initio calculations of Jaffe et al. for 1,2dimethoxyethane. The resulting torsionmodified AM1 Hamiltonian is used with our previous QM/MMmore »

Multiple environment single system quantum mechanical/molecular mechanical (MESSQM/MM) calculations. 1. Estimation of polarization energies
In combined quantum mechanical/molecular mechanical (QM/MM) free energy calculations, it is often advantageous to have a frozen geometry for the quantum mechanical (QM) region. For such multipleenvironment singlesystem (MESS) cases, two schemes are proposed here for estimating the polarization energy: the first scheme, termed MESSE, involves a Roothaan step extrapolation of the selfconsistent field (SCF) energy; whereas the other scheme, termed MESSH, employs a Newton–Raphson correction using an approximate inverse electronic Hessian of the QM region (which is constructed only once). Both schemes are extremely efficient, because the expensive Fock updates and SCF iterations in standard QM/MM calculations are completelymore »Cited by 10 
Molecular dynamics study of response of liquid N,Ndimethylformamide to externally applied electric field using a polarizable force field
The behavior of Liquid N,Ndimethylformamide subjected to a wide range of externally applied electric fields (from 0.001 V/nm to 1 V/nm) has been investigated through molecular dynamics simulation. To approach the objective the AMOEBA polarizable force field was extended to include the interaction of the external electric field with atomic partial charges and the contribution to the atomic polarization. The simulation results were evaluated with quantum mechanical calculations. The results from the present force field for the liquid at normal conditions were compared with the experimental and molecular dynamics results with nonpolarizable and other polarizable force fields. The uniform externalmore » 
Hybrid Quantum Mechanical and Molecular Mechanics Study of the SN2 Reaction of CCl4 + OH in Aqueous Solution: The Potential of Mean Force, Reaction Energetics, and Rate Constants
The bimolecular nucleophilic substitution reaction of CCl{sub 4} and OH{sup } in aqueous solution was investigated on the basis of a combined quantum mechanical and molecular mechanics method. A multilayered representation approach is employed to achieve high accuracy results at the CCSD(T) level of theory. The potential of mean force calculations at the DFT level and CCSD(T) level of theory yield reaction barrier heights of 22.7 and 27.9 kcal/mol, respectively. Both the solvation effects and the solventinduced polarization effect have significant contributions to the reaction energetics, for example, the solvation effect raises the saddle point by 10.6 kcal/mol. The calculatedmore » 
Determining polarizable force fields with electrostatic potentials from quantum mechanical linear response theory
We developed a new method to calculate the atomic polarizabilities by fitting to the electrostatic potentials (ESPs) obtained from quantum mechanical (QM) calculations within the linear response theory. This parallels the conventional approach of fitting atomic charges based on electrostatic potentials from the electron density. Our ESP fitting is combined with the induced dipole model under the perturbation of uniform external electric fields of all orientations. QM calculations for the linear response to the external electric fields are used as input, fully consistent with the induced dipole model, which itself is a linear response model. The orientation of the uniformmore »