A polarizable force field for water using an artificial neural network

A force field for liquid water including polarization effects has been constructed using an artificial neural network (ANN). It is essential to include a many-body polarization effect explicitly into a potential energy function in order to treat liquid water which is dense and highly polar. The new potential energy function is a combination of empirical and nonempirical potentials. The TIP4P model was used for the empirical part of the potential. For the nonempirical part, an ANN with a back-propagation of error algorithm (BPNN) was introduced to reproduce the complicated many-body interaction energy surface from ab initio quantum mechanical calculations. BPNN, described in terms of a matrix, provides enough flexibility to describe the complex potential energy surface (PES). The structural and thermodynamic properties, calculated by isobaric-isothermal (constant-NPT) Monte Carlo simulations with the new polarizable force field for water, are compatible with experimental results. Thus, the simulation establishes the validity of using our estimated PES with a polarization effect for accurate predictions of liquid state properties. Applications of this approach are simple and systematic so that it can easily be applied to the development of other force fields besides the water-water system.     

Effect of interprotein polarization on protein–protein binding energy

Molecular dynamics simulation in explicit water for the binding of the benchmark barnase‐barstar complex was carried out to investigate the effect polarization of interprotein hydrogen bonds on its binding free energy. Our study is based on the AMBER force field but with polarized atomic charges derived from fragment quantum mechanical calculation for the protein complex. The quantum‐derived atomic charges include the effect of polarization of interprotein hydrogen bonds, which was absent in the standard force fields that were used in previous theoretical calculations of barnase‐barstar binding energy. This study shows that this polarization effect impacts both the static (electronic) and dynamic interprotein electrostatic interactions and significantly lowers the free energy of the barnase‐barstar complex. © 2012 Wiley Periodicals, Inc.     

Nano‐Ferroelectric for High Efficiency Overall Water Splitting under Ultrasonic Vibration

Piezocatalysis, converting mechanical vibration into chemical energy, has emerged as a promising candidate for water‐splitting technology. However, the efficiency of the hydrogen production is quite limited. We herein report well‐defined 10 nm BaTiO₃ nanoparticles (NPs) characterized by a large electro‐mechanical coefficient which induces a high piezoelectric effect. Atomic‐resolution high angle annular dark field scanning transmission electron microscopy (HAADF‐STEM) and scanning probe microscopy (SPM) suggests that piezoelectric BaTiO₃ NPs display a coexistence of multiple phases with low energy barriers and polarization anisotropy which results in a high electro‐mechanical coefficient. Landau free energy modeling also confirms that the greatly reduced polarization anisotropy facilitates polarization rotation. Employing the high piezoelectric properties of BaTiO₃ NPs, we demonstrate an overall water‐splitting process with the highest hydrogen production efficiency hitherto reported, with a H₂ production rate of 655 μmol g^?1 h^?1, which could rival excellent photocatalysis system. This study highlights the potential of piezoelectric catalysis for overall water splitting.     

Toward a consistent treatment of polarization in model QM/MM calculations

The concept of model chemistries within hybrid QM/MM calculations has been addressed through analysis of the polarization energy determined by two distinct approaches based on (i) induced charges and (ii) induced dipoles. The quantum mechanical polarization energy for four configurations of the water dimer has been determined for a range of basis sets using Morokuma energy decomposition analysis. This benchmark value has been compared to the fully classical polarization energy determined using the induced dipole approach, and the molecular mechanics polarization energy calculated using induced charges within the MM region of hybrid QM/MM calculations. From the water dimer calculations, it is concluded that the induced charge approach is consistent with medium sized basis set calculations whereas the induced dipole approach is consistent with large basis set calculations. This result is highly relevant to the concept of QM/MM model chemistries.     

Key role of the polarization anisotropy of water in modeling classical polarizable force fields

We have evaluated the extent to which classical polarizable force fields, based either on the chemical potential equalization principle or on distributed polarizabilities in the framework of the Sum of Interactions Between Fragments Ab initio computed (SIBFA), can reproduce the ab initio polarization energy and the dipole moment of three distinct water oligomers: bifurcated chains, transverse hydrogen-bonded chains, and longitudinal hydrogen-bonded chains of helical shape. To analyze the many-body polarization effect, chains of different size, i.e., from 2 to 12 water monomers, have been considered. Although the dipole moment is a well-defined quantity in both classical polarizable models and quantum mechanical methods, polarization energy can be defined unequivocally only in the former type of approaches. In this study we have used the Kitaura-Morokuma (KM) procedure. Although the KM approach is on the one hand known to overestimate the polarization energy for strongly interacting molecules, on the other hand it can account for the many-body polarization effectively, whereas some other procedures do not. Our data show that, if off-centered lone pair polarizabilities are explicitly represented, classical polarizable force fields can afford a close agreement with the ab initio results, both in terms of polarization energy and in terms of dipole moment.     

Water under temperature gradients: polarization effects and microscopic mechanisms of heat transfer

We report non-equilibrium molecular dynamics simulations (NEMD) of water under temperature gradients using a modified version of the central force model (MCFM). This model is very accurate in predicting the equation of state of water for a wide range of pressures and temperatures. We investigate the polarization response of water to thermal gradients, an effect that has been recently predicted using Non-Equilibrium Thermodynamics (NET) theory and computer simulations, as a function of the thermal gradient strength. We find that the polarization of the liquid varies linearly with the gradient strength, which indicates that the ratio of phenomenological coefficients regulating the coupling between the polarization response and the heat flux is independent of the gradient strength investigated. This notion supports the NET theoretical predictions. The coupling effect leading to the liquid polarization is fairly strong, leading to polarization fields of ~10(3-6) V m(-1) for gradients of ~10(5-8) K m(-1), hence confirming earlier estimates. Finally we employ our NEMD approach to investigate the microscopic mechanism of heat transfer in water. The image emerging from the computation and analysis of the internal energy fluxes is that the transfer of energy is dominated by intermolecular interactions. For the MCFM model, we find that the contribution from hydrogen and oxygen is different, with the hydrogen contribution being larger than that of oxygen.     

Solvent effects on the nπ* transition of pyrimidine in aqueous solution

A hybrid quantum mechanical and molecular mechanical potential is used in Monte Carlo simulations to examine the solvent effects on the electronic excitation energy for the n→π* transition of pyrimidine in aqueous solution. In the present study, the pyrimidine molecule is described by the semi-empirical AM1 model, while the solvent molecules are treated classically. Two sets of calculations are performed: the first involves the use of the pairwise three-point charge TIP3P model for water, and the second computation employs a polarizable many-body potential for the solvent. The latter calculation takes into account the effect of solvent polarization following the solute electronic excitation, and makes a correction to the energies determined using pairwise potentials, which neglects such fast polarization effects and overestimates the solute-solvent interactions on the Franck-Condon excited states. Our simulation studies of pyrimidine in water indicate that the solvent charge redistribution following the solute electronic excitation makes modest corrections (about −130␣cm⁻¹) to the energy predicted by using pairwise potentials. Specific hydrogen bonding interactions between pyrimidine and water are important for the prediction of solvatochromic shifts for pyrimidine. The computed n→π* blue shift is 2275±110 cm⁻¹, which may be compared with the experimental value (2700 cm⁻¹) from isooctane to water. Received: 14 January 1997 / Accepted: 21 February 1997     

Polarizable solute in polarizable and flexible solvents: simulation study of electron transfer reaction systems

A polarizable solute model, based on the empirical valence bond approach, is developed and applied to electron transfer (ET) reactions in polarizable and flexible water solvents. The polarization effect is investigated in comparison with a nonpolarizable solute and solvent model. With free energy curves constructed by a molecular dynamics simulation, the activation energy barrier and the reorganization energy related to ET processes are investigated. The present simulation results show that the activation energy barrier becomes larger in the polarizable model than in the nonpolarizable model and that this makes the ET rate slower than that with the nonpolarizable model. It is shown that the effect of the electronic energy difference of solute molecule on free energy profiles is remarkable and that, corresponding to this effect, the reorganization energy is significantly modified. These results indicate that the process of solvent polarization by the polarized solute to enhance the solute-solvent interaction is a key factor and that treating the polarization of both solute and solvent at the same time is essential. Also, the polarization effect on the diffusive motion of the solute molecule in the polarization solvent is studied. The polarized solute molecule shows slower diffusive motion compared with that in the nonpolarizable model.     

Structure and electronic properties of a benzene-water solution

Electronic properties of benzene in water were investigated by a sequential quantum mechanical/molecular dynamics approach. Emphasis was placed on the analysis of the structure, polarization effects, and ionization spectrum. By adopting a polarizable model for both benzene and water the structure of the benzene-water solution is in good agreement with data from first principles molecular dynamics. Further, strong evidence that water molecules acquire enhanced orientational order near the benzene molecule is found. Upon hydration, the quadrupole moment of benzene is not significantly changed in comparison with the gas-phase value. We are also reporting results for the dynamic polarizability of benzene in water. Our results indicate that the low energy behaviour of the dynamic polarizability of gas-phase and hydrated benzene is quite similar. Outer valence Green's function calculations for benzene in liquid water show a splitting of the gas-phase energy levels associated with the 1e(1g)(π), 2e(2g), and 2e(1u) orbitals upon hydration. Lifting of the orbitals degeneracy and redshift of the outer valence bands is related to symmetry breaking of the benzene structure in solution and polarization effects from the surrounding water molecules.     

Energy components of aqueous solution: Insight from hybrid QM/MM simulations using a polarizable solvent model

An energy decomposition method is present in statistical Monte Carlo simulations of aqueous solutions of a series of organic solutes, making use of a hybrid quantum mechanical and polarizable molecular mechanical (QM/MM-PIPF) approach. In the hybrid QM/MM-PIPF method, the mutual solute–solvent polarization effect is specifically considered through a coupled iterative procedure that ensures the convergence of solvent induced dipoles and the solute wave function. It should be noted that the method is an approximate approach without specifically considering the electronic correlation effect between solute and solvent electrons, and energetic results have not been verified by free energy calculations. Nevertheless, the energy decomposition analysis provides insight into the details of the molecular polarization effect. Qualitative trends of the energy components from such analyses provide guidance in the understanding of the nature of intermolecular interactions in biomolecular systems, whereas quantitative results on specific terms may be utilized to derive empirical, yet computationally more efficient, force fields. Polarization effects are found to be significant, which contribute 10% to 23% to the total solute–solvent interaction energies. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 :1061–1071, 1997     

Charge-transfer excitons in DNA

There have been a number of theoretical treatments of excitons in DNA, most neglecting both the intrachain and interchain wavefunction overlaps of the electron and hole, treating them as Frenkel excitons. Recently, the importance of the intrachain and interchain coupling has been highlighted. Experiments have shown that in (dA)n oligomers and in duplex (dA)n.(dT)n, to be abbreviated (A/T), where A is adenine and T is thymine, the exciton wavefunction is delocalized over several bases. In duplexes it is possible to have charge-transfer (CT) excitons. Theoretical calculations have suggested that CT excitons in DNA may have lower energy than single chain excitons. In all the calculations of excitons in DNA, the polarization of the surrounding water has been neglected. Calculations have shown, however, that polarization of the water by an excess electron or a hole in DNA lowers its energy by approximately 1/2 eV, causing it to become a polaron. It is therefore to be expected that polarization charge induced in the surrounding water has a significant effect on the properties of the exciton. In what follows, we present calculations of some properties CT excitons would have in an A/T duplex taking into account the wavefunction overlaps, the effect of the surrounding water, which results in the electron and hole becoming polarons, and the ions in the water. As expected, the CT exciton has lowest energy when the electron and hole polarons are directly opposite each other. By appropriate choice of the dielectric constant, we can obtain a CT exciton delocalized over the number of sites found in photoinduced absorption experiments. The absorption threshold that we then calculate for CT exciton creation in A/T is in reasonable agreement with the lowest singlet absorption deduced from available data.     

A novel quantum mechanical/molecular mechanical approach to the free energy calculation for isomerization of glycine in aqueous solution

The free energy change associated with the isomerization reaction of glycine in water solution has been studied by a hybrid quantum mechanical/molecular mechanical (QM/MM) approach combined with the theory of energy representation (QM/MM-ER) recently developed. The solvation free energies for both neutral and zwitterionic form of glycine have been determined by means of the QM/MM-ER simulation. The contributions of the electronic polarization and the fluctuation of the QM solute to the solvation free energy have been investigated. It has been found that the contribution of the density fluctuation of the zwitterionic solute is estimated as -4.2 kcal/mol in the total solvation free energy of -46.1 kcal/mol, while that of the neutral form is computed as -3.0 kcal/mol in the solvation free energy of -15.6 kcal/mol. The resultant free energy change associated with the isomerization of glycine in water has been obtained as -7.8 kcal/mol, in excellent agreement with the experimental data of -7.3 or -7.7 kcal/mol, implying the accuracy of the QM/MM-ER approach. The results have also been compared with those computed by other methodologies such as the polarizable continuum model and the classical molecular simulation. The efficiency and advantage of the QM/MM-ER method has been discussed.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

Interactions between phosphate and water in solution: a natural bond orbital based analysis in a QM/MM framework

The natural bond orbital (NBO) and natural energy decomposition analysis (NEDA) calculations are used to analyze the interaction between mono-methyl phosphate-ester (MMP) and its solvation environment in a combined quantum mechanical/molecular mechanical (QM/MM) framework. The solute-solvent configurations are generated using a specific parametrization of the self-consistent-charge density functional tight-binding (SCC-DFTB) model for the MMP and TIP3P for water. The NBO and NEDA calculations are done with several QM/MM partitioning schemes with HF/6-31+G** as the QM level. Regardless of the size of the QM region, a notable amount of charge transfer is observed between MMP and the neighboring water molecules and the charge-transfer interactions are, in the NEDA framework, as important as the electric (electrostatic and polarization) components. This work illustrates that NBO based analyses are effective tools for probing intermolecular interactions in condensed phase systems.     

A multilayered representation,quantum mechanical and molecular mechanics study of the CH3F + OH− reaction in water

The bimolecular nucleophilic substitution (S_N2) reaction of CH₃F + OH^? in aqueous solution was investigated using a combined quantum mechanical and molecular mechanics approach. Reactant complex, transition state, and product complex along the reaction pathway were analyzed in water. The potentials of mean force were calculated using a multilayered representation with the DFT and CCSD(T) level of theory for the reactive region. The obtained free energy activation barrier for this reaction at the CCSD(T)/MM representation is 18.3 kcal/mol which agrees well with the experimental value at ～21.6 kcal/mol. Both the solvation effect and solute polarization effect play key roles on raising the activation barrier height in aqueous solution, with the former raising the barrier height by 3.1 kcal/mol, the latter 1.5 kcal/mol. © 2013 Wiley Periodicals, Inc.     

Variational calculation of quantum mechanical/molecular mechanical free energy with electronic polarization of solvent

Quantum mechanical/molecular mechanical (QM/MM) free energy calculation presents a significant challenge due to an excessive number of QM calculations. A useful approach for reducing the computational cost is that based on the mean field approximation to the QM subsystem. Here, we describe such a mean-field QM/MM theory for electronically polarizable systems by starting from the Hartree product ansatz for the total system and invoking a variational principle of free energy. The MM part is then recast to a classical polarizable model by introducing the charge response kernel. Numerical test shows that the potential of mean force (PMF) thus obtained agrees quantitatively with that obtained from a direct QM/MM calculation, indicating the utility of self-consistent mean-field approximation. Next, we apply the obtained method to prototypical reactions in several qualitatively different solvents and make a systematic comparison of polarization effects. The results show that in aqueous solution the PMF does not depend very much on the water models employed, while in nonaqueous solutions the PMF is significantly affected by explicit polarization. For example, the free energy barrier for a phosphoryl dissociation reaction in acetone and cyclohexane is found to increase by more than 10 kcal/mol when switching the solvent model from an empirical to explicitly polarizable one. The reason for this is discussed based on the parametrization of empirical nonpolarizable models.     

Incorporation of a QM/MM buffer zone in the variational double self-consistent field method

The explicit polarization (X-Pol) potential is an electronic-structure-based polarization force field, designed for molecular dynamics simulations and modeling of biopolymers. In this approach, molecular polarization and charge transfer effects are explicitly treated by a combined quantum mechanical and molecular mechanical (QM/MM) scheme, and the wave function of the entire system is variationally optimized by a double self-consistent field (DSCF) method. In the present article, we introduce a QM buffer zone for a smooth transition from a QM region to an MM region. Instead of using the Mulliken charge approximation for all QM/MM interactions, the Coulombic interactions between the adjacent fragments are determined directly by electronic structure theory. The present method is designed to accelerate the speed of convergence of the total energy and charge density of the system.     

Simple and exact approach to the electronic polarization effect on the solvation free energy: formulation for quantum-mechanical/molecular-mechanical system and its applications to aqueous solutions

We present a simple and exact numerical approach to compute the free energy contribution δμ in solvation due to the electron density polarization and fluctuation of a quantum-mechanical solute in the quantum-mechanical/molecular-mechanical (QM/MM) simulation combined with the theory of the energy representation (QM/MM-ER). Since the electron density fluctuation is responsible for the many-body QM-MM interactions, the standard version of the energy representation method cannot be applied directly. Instead of decomposing the QM-MM polarization energy into the pairwise additive and non-additive contributions, we take sum of the polarization energies in the QM-MM interaction and adopt it as a new energy coordinate for the method of energy representation. Then, it is demonstrated that the free energy δμ can be exactly formulated in terms of the energy distribution functions for the solution and reference systems with respect to this energy coordinate. The benchmark tests were performed to examine the numerical efficiency of the method with respect to the changes in the individual properties of the solvent and the solute. Explicitly, we computed the solvation free energy of a QM water molecule in ambient and supercritical water, and also the free-energy change associated with the isomerization reaction of glycine from neutral to zwitterionic structure in aqueous solution. In all the systems examined, it was demonstrated that the computed free energy δμ agrees with the experimental value, irrespective of the choice of the reference electron density of the QM solute. The present method was also applied to a prototype reaction of adenosine 5'-triphosphate hydrolysis where the effect of the electron density fluctuation is substantial due to the excess charge. It was demonstrated that the experimental free energy of the reaction has been accurately reproduced with the present approach.     

Novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: free energy calculation for the Beckmann rearrangement promoted by proton transfers in the supercritical water

The Beckmann rearrangement of acetone oxime promoted by proton transfers in the supercritical water has been investigated by means of the hybrid quantum mechanical/molecular mechanical approach combined with the theory of energy representation (QM/MM-ER) recently developed. The transition state (TS) structures have been explored by ab initio calculations for the reaction of hydrated acetone oxime on the assumption that the reaction is catalyzed by proton transfers along the hydrogen bonds connecting the solute and the solvent water molecules. Up to two water molecules have been considered as reactants that take part in the proton transfers. As a result of the density functional theory calculations with B3LYP functional and aug-cc-pVDZ basis set, it has been found that participation of two water molecules in the reaction reduces the activation free energy by -12.3 kcal/mol. Furthermore, the QM/MM-ER simulations have revealed that the TS is more stabilized than the reactant state in the supercritical water by 2.7 kcal/mol when two water molecules are involved in the reaction. Solvation free energies of the reactant and the TS have been decomposed into terms due to the electronic polarization of the solute, electron density fluctuation, and others to elucidate the origin of the stabilization of the TS as compared with the reactant. It has been revealed that the promotion of the chemical reaction due to the hydration mainly originates from the interaction between the nonpolarized solute and the solvent water molecules at the supercritical state.