Quantum mechanical study of solvent effects in a prototype SN2 reaction in solution: Cl- attack on CH3Cl

Accelerated Quantum Mechanics/Molecular Mechanics Simulations via Neural Networks Incorporated with Mechanical Embedding Scheme

A powerful tool to study the mechanism of reactions in solutions or enzymes is to perform the ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations. However, the computational cost is too high due to the explicit electronic structure calculations at every time step of the simulation. A neural network (NN) method can accelerate the QM/MM-MD simulations, but it has long been a problem to accurately describe the QM/MM electrostatic coupling by NN in the electrostatic embedding (EE) scheme. In this work, we developed a new method to accelerate QM/MM calculations in the mechanic embedding (ME) scheme. The potentials and partial point charges of QM atoms are first learned in vacuo by the embedded atom neural networks (EANN) approach. MD simulations are then performed on this EANN/MM potential energy surface (PES) to obtain free energy (FE) profiles for reactions, in which the QM/MM electrostatic coupling is treated in the mechanic embedding (ME) scheme. Finally, a weighted thermodynamic perturbation (wTP) corrects the FE profiles in the ME scheme to the EE scheme. For two reactions in water and one in methanol, our simulations reproduced the B3LYP/MM free energy profiles within 0.5 kcal/mol with a speed-up of 30-60-fold. The results show that the strategy of combining EANN potential in the ME scheme with the wTP correction is efficient and reliable for chemical reaction simulations in liquid. Another advantage of our method is that the QM PES is independent of the MM subsystem, so it can be applied to various MM environments as demonstrated by an S_N2 reaction studied in water and methanol individually, which used the same EANN PES. The free energy profiles are in excellent accordance with the results obtained from B3LYP/MM-MD simulations. In future, this method will be applied to the reactions of enzymes and their variants.

Tackling Solvent Effects by Coupling Electronic and Molecular Density Functional Theory

Solvation effects can have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike a previous work, both densities are minimized self-consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effects in quantum calculation. It is computationally more effective than a direct quantum mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods, it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data, but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway.

A reaction density functional theory study of the solvent effect in prototype SN2 reactions in aqueous solution

Reaction density functional theory (RxDFT), combining quantum DFT with classical DFT, has been employed to investigate the solvent effect and free energy profiles of S_N2 reactions in aqueous solution.

Solvent Effect on the Potential Energy Surfaces of the F- + CH3CH2Br Reaction

Although substantial work has been undertaken on reaction pathways involved in base-promoted elimination reactions and bimolecular nucleophilic substitution reaction of F^- on CH₃CH₂X (X = Cl, Br, I), the effect of solvents with varying dielectric constants on the stereochemistry of each of the reaction species involved across the reaction profile have not yet been clearly understood. The present investigation reports the effect of solvents on the potential energy surfaces (PES) and structures of the species appearing in the reaction pathway of F^- with bromoethane. The PESs in the gas phase have been computed at MP2 level and CCSD(T) level. The performance of several hybrid density functional, such as B3LYP, M06, M06L, BHandH, X3LYP, M05, M05-2X, and M06-2X have also been investigated toward describing the elimination and nucleophilic substitution reactions. With respect to MAE values and to make the computation cost-effective, we have explored the implicit continuum solvent model, CPCM in solvents like cyclohexane, methanol, acetonitrile, dimethyl sulfoxide and water. The reactant complexes proceed through the subsequent steps to produce fluoroethane as the substitution product and ethylene as one of the elimination products. For elimination reaction both syn and anti elimination have been explored. The calculated relatives energies values, which are negative in the gas phase, are found to be positive in polar solvents since the point charge in the separated reactants are more stabilized than the dispersed charge in the transient complex, which has also been analyzed through NBO analysis.

Nucleophilic Substitution in Solution: Activation Strain Analysis of Weak and Strong Solvent Effects

We have quantum chemically studied the effect of various polar and apolar solvents on the shape of the potential energy surface (PES) of a diverse collection of archetypal nucleophilic substitution reactions at carbon, silicon, phosphorus, and arsenic by using density functional theory at the OLYP/TZ2P level. In the gas phase, all our model SN 2 reactions have single-well PESs, except for the nucleophilic substitution reaction at carbon (SN 2@C), which has a double-well energy profile. The presence of the solvent can have a significant effect on the shape of the PES and, thus, on the nature of the SN 2 process. Solvation energies, charges on the nucleophile or leaving group, and structural features are compared for the various SN 2 reactions in a spectrum of solvents. We demonstrate how solvation can change the shape of the PES, depending not only on the polarity of the solvent, but also on how the charge is distributed over the interacting molecular moieties during different stages of the reaction. In the case of a nucleophilic substitution at three-coordinate phosphorus, the reaction can be made to proceed through a single-well [no transition state (TS)], bimodal barrier (two TSs), and then through a unimodal transition state (one TS) simply by increasing the polarity of the solvent.

Acceleration of Semiempirical QM/MM Methods through Message Passage Interface (MPI), Hybrid MPI/Open Multiprocessing, and Self-Consistent Field Accelerator Implementations

The strategy and implementation of scalable and efficient semiempirical (SE) QM/MM methods in CHARMM are described. The serial version of the code was first profiled to identify routines that required parallelization. Afterward, the code was parallelized and accelerated with three approaches. The first approach was the parallelization of the entire QM/MM routines, including the Fock matrix diagonalization routines, using the CHARMM message passage interface (MPI) machinery. In the second approach, two different self-consistent field (SCF) energy convergence accelerators were implemented using density and Fock matrices as targets for their extrapolations in the SCF procedure. In the third approach, the entire QM/MM and MM energy routines were accelerated by implementing the hybrid MPI/open multiprocessing (OpenMP) model in which both the task- and loop-level parallelization strategies were adopted to balance loads between different OpenMP threads. The present implementation was tested on two solvated enzyme systems (including <100 QM atoms) and an S_N2 symmetric reaction in water. The MPI version exceeded existing SE QM methods in CHARMM, which include the SCC-DFTB and SQUANTUM methods, by at least 4-fold. The use of SCF convergence accelerators further accelerated the code by ∼12-35% depending on the size of the QM region and the number of CPU cores used. Although the MPI version displayed good scalability, the performance was diminished for large numbers of MPI processes due to the overhead associated with MPI communications between nodes. This issue was partially overcome by the hybrid MPI/OpenMP approach which displayed a better scalability for a larger number of CPU cores (up to 64 CPUs in the tested systems).

Revisiting the Dielectric Constant Effect on the Nucleophile and Leaving Group of Prototypical Backside SN2 Reactions: A Reaction Force and Atomic Contribution Analysis

The solvent effect on the nucleophile and leaving group atoms of the prototypical F^- + CH₃Cl → CH₃F + Cl^- backside bimolecular nucleophilic substitution reaction (S_N2) is analyzed employing the reaction force and the atomic contributions methods on the intrinsic reaction coordinate (IRC). Solvent effects were accounted for using the polarizable continuum solvent model. Calculations were performed employing 11 dielectric constants, ε, ranging from 1.0 to 78.5, to cover a wide spectrum of solvents. The reaction force data reveal that the solvent mainly influences the region of the IRC preceding the energy barrier, where the structural rearrangement to reach the transition state occurs. A detailed analysis of the atomic role in the reaction as a function of ε reveals that the nucleophile and the carbon atom are the ones that contribute the most to the energy barrier. In addition, we investigated the effect of the choice of nucleophile and leaving group on the ΔE₀ and ΔE^‡ of Y^- + CH₃X → YCH₃ + X^- (X, Y = F, Cl, Br, I) in aqueous solution. Our analysis allowed us to find relationships between the atomic contributions to the activation energy and leaving group ability and nucleophilicity.

Charge-dependent many-body exchange and dispersion interactions in combined QM/MM simulations

Accurate modeling of the molecular environment is critical in condensed phase simulations of chemical reactions. Conventional quantum mechanical/molecular mechanical (QM/MM) simulations traditionally model non-electrostatic non-bonded interactions through an empirical Lennard-Jones (LJ) potential which, in violation of intuitive chemical principles, is bereft of any explicit coupling to an atom's local electronic structure. This oversight results in a model whereby short-ranged exchange-repulsion and long-ranged dispersion interactions are invariant to changes in the local atomic charge, leading to accuracy limitations for chemical reactions where significant atomic charge transfer can occur along the reaction coordinate. The present work presents a variational, charge-dependent exchange-repulsion and dispersion model, referred to as the charge-dependent exchange and dispersion (QXD) model, for hybrid QM/MM simulations. Analytic expressions for the energy and gradients are provided, as well as a description of the integration of the model into existing QM/MM frameworks, allowing QXD to replace traditional LJ interactions in simulations of reactive condensed phase systems. After initial validation against QM data, the method is demonstrated by capturing the solvation free energies of a series of small, chlorine-containing compounds that have varying charge on the chlorine atom. The model is further tested on the SN2 attack of a chloride anion on methylchloride. Results suggest that the QXD model, unlike the traditional LJ model, is able to simultaneously obtain accurate solvation free energies for a range of compounds while at the same time closely reproducing the experimental reaction free energy barrier. The QXD interaction model allows explicit coupling of atomic charge with many-body exchange and dispersion interactions that are related to atomic size and provides a more accurate and robust representation of non-electrostatic non-bonded QM/MM interactions.

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

The effect of the environment on the methyl transfer reaction mechanism between trimethylsulfonium and phenolate

The methyl transfer reaction mechanism in different molecular environments were studied by electronic structure methods and QM/MM molecular dynamics simulations.

Why many semiempirical molecular orbital theories fail for liquid water and how to fix them

Water is an extremely important liquid for chemistry and the search for more accurate force fields for liquid water continues unabated. Neglect of diatomic differential overlap (NDDO) molecular orbital methods provide and intriguing generalization of classical force fields in this regard because they can account both for bond breaking and electronic polarization of molecules. However, we show that most standard NDDO methods fail for water because they give an incorrect description of hydrogen bonding, water's key structural feature. Using force matching, we design a reparameterized NDDO model and find that it qualitatively reproduces the experimental radial distribution function of water, as well as various monomer, dimer, and bulk properties that PM6 does not. This suggests that the apparent limitations of NDDO models are primarily due to poor parameterization and not to the NDDO approximations themselves. Finally, we identify the physical parameters that most influence the condensed phase properties. These results help to elucidate the chemistry that a semiempirical molecular orbital picture of water must capture. We conclude that properly parameterized NDDO models could be useful for simulations that require electronically detailed explicit solvent, including the calculation of redox potentials and simulation of charge transfer and photochemistry.