Browsing by Subject "Polarizability"

Now showing 1 - 2 of 2
  • Thumbnail Image
    Item
    Combining Semiempirical QM Methods with Atom Dipole Interaction Model for Accurate and Efficient Polarizability Calculations
    (2022-12) Young, Ryan; Pu, Jingzhi; Long, Eric C.; Naumann, Christoph; Sardar, Rajesh
    Molecular polarizability plays a significant role in chemistry, biology, and medicine. Classical prediction of polarizability often relies on atomic-type specific polarizability optimized for training set molecules, which limits the calculations to systems of similar chemical environment. Although ab initio (AI) quantum mechanical (QM) methods are more transferable in predicting molecular polarizability, their high computational costs especially when used with large basis sets for obtaining quantitatively reliable results make them less practical. To obtain accurate QM polarizability in an efficient manner, we have developed a dual-level approach, where the polarizability (α) obtained from the efficient semiempirical QM (SE) method is corrected using a set of element-base atomic polarizabilities derived from the atomic dipole interaction model (ADIM) to reproduce the density functional theory (DFT) results. We have optimized the atomic polarizability correction parameters for CHON-containing systems using a small training set of molecules and tested the resulting SE-ADIM model on the neutral drug-like molecules in the QM7B database. SE-ADIM corrected AM1 showed substantial improvement with its relative percent error (RPE) compared to B3LYP reduced from 33.81% to 3.35%. To further test its robustness for larger molecules in broader chemical bonding situations, we applied this method to a collection of drug molecules from the e-Drug3D database. For the 1004 molecules tested, our SE-ADIM model, which only contains four empirical parameters, greatly reduces the RPE in AM1 polarizability relative to B3LYP from 26.8% to 2.9%. Error decomposition shows consistent improvements across molecules with diverse bond saturations, molecular sizes, and charge states. In addition, we have applied AlphaML, a promising machine learning (ML) technique for predicting molecular polarizability, to the e-Drug3D dataset to compare its performance with our SE-ADIM correction of AM1. We found SE-ADIM performs competitively with AlphaML bolstering our confidence in the value of our method. Errors distinct to AlphaML were also discovered. We found four molecules for which AlphaML predicts negative molecular polarizabilities, all of which were peroxides. In contrast, SE-ADIM has no such issue with these molecules or this chemical type. Finally, to improve performance of SE-ADIM when correcting AM1 molecular polarizability calculations for charged molecules, we introduce a charge dependent polarizability (CDP) enabled SE-ADIM. Training the CDP enabled SE-ADIM with a single additional parameter, B, we were able to reduce error in AM1 molecular polarizability calculations of charged molecules relative to B3LYP from 29.57% to 5.16%. By contrast, SE-ADIM without CDP corrected AM1 relative to B3LYP had an RPE of 8.56%. The most benefit of CDP was evident within negatively charged molecules where AM1 error relative to B3LYP fell from 32.20% to 3.77% while SE-ADIM without CDP enabled error for these same negative molecules was 10.06%.
  • Thumbnail Image
    Item
    Doubly Polarized QM/MM with Machine Learning Chaperone Polarizability
    (American Chemical Society, 2021) Kim, Bryant; Shao, Yihan; Pu, Jingzhi; Chemistry and Chemical Biology, School of Science
    A major shortcoming of semiempirical (SE) molecular orbital methods is their severe underestimation of molecular polarizability compared with experimental and ab initio (AI) benchmark data. In a combined quantum mechanical and molecular mechanical (QM/MM) treatment of solution-phase reactions, solute described by SE methods therefore tends to generate inadequate electronic polarization response to solvent electric fields, which often leads to large errors in free energy profiles. To address this problem, here we present a hybrid framework that improves the response property of SE/MM methods through high-level molecular-polarizability fitting. Specifically, we place on QM atoms a set of corrective polarizabilities (referred to as chaperone polarizabilities), whose magnitudes are determined from machine learning (ML) to reproduce the condensed-phase AI molecular polarizability along the minimum free energy path. These chaperone polarizabilities are then used in a machinery similar to a polarizable force field calculation to compensate for the missing polarization energy in the conventional SE/MM simulations. Because QM atoms in this treatment host SE wave functions as well as classical polarizabilities, both polarized by MM electric fields, we name this method doubly polarized QM/MM (dp-QM/MM). We demonstrate the new method on the free energy simulations of the Menshutkin reaction in water. Using AM1/MM as a base method, we show that ML chaperones greatly reduce the error in the solute molecular polarizability from 6.78 to 0.03 Å3 with respect to the density functional theory benchmark. The chaperone correction leads to ~10 kcal/mol of additional polarization energy in the product region, bringing the simulated free energy profiles to closer agreement with the experimental results. Furthermore, the solute-solvent radial distribution functions show that the chaperone polarizabilities modify the free energy profiles through enhanced solvation corrections when the system evolves from the charge-neutral reactant state to the charge-separated transition and product states. These results suggest that the dp-QM/MM method, enabled by ML chaperone polarizabilities, provides a very physical remedy for the underpolarization problem in SE/MM-based free energy simulations.