On the performance of molecular polarization methods. II. Water and carbon tetrachloride close to a cation

Our initial study on the performance of molecular polarization methods close to a positive point charge [M. Masia, M. Probst, and R. Rey, J. Chem. Phys. 121, 7362 (2004)] is extended to the case in which a molecule interacts with a real cation. Two different methods (point dipoles and shell model) are applied to both the ion and the molecule. The results are tested against high-level ab initio calculations for a molecule (water or carbon tetrachloride) close to Li+, Na+, Mg2+, and Ca2+. The monitored observable is in all cases the dimer electric dipole as a function of the ion-molecule distance for selected molecular orientations. The moderate disagreement previously obtained for point charges at intermediate distances, and attributed to the linearity of current polarization methods (as opposed to the nonlinear effects evident in ab initio calculations), is confirmed for real cations as well. More importantly, it is found that at short separations the phenomenological polarization methods studied here substantially overestimate the dipole moment induced if the ion is described quantum chemically as well, in contrast to the dipole moment induced by a point-charge ion, for which they show a better degree of accord with ab initio results. Such behavior can be understood in terms of a decrease of atomic polarizabilities due to the repulsion between electronic charge distributions at contact separations. It is shown that a reparametrization of the Thole method for damping of the electric field, used in conjunction with any polarization scheme, allows to satisfactorily reproduce the dimer dipole at short distances. In contrast with the original approach (developed for intramolecular interactions), the present reparametrization is ion and method dependent, and corresponding parameters are given for each case.     

Dipole preserving and polarization consistent charges

A method for estimating dipole preserving and polarization consistent (DPPC) charges is described, which reproduces exactly the molecular dipole moment as well as the local, atomic hybridization dipoles determined from the corresponding wave function and can yield accurate molecular polarization. The method is based on a model described by Thole and van Duijnen and a new feature is introduced to treat molecular polarization. Thus, the DPPC method offers a convenient procedure to describe molecular polarization in applications using semiempirical models and ab initio molecular orbital theory with relatively small basis functions such as 6‐31+G(d,p) or without inclusion of electron correlation; these methods tend to underestimate molecular polarizability. The trends of the DPPC partial atomic charges are found to be in good accord with those of the CM2 model, a class IV charge analysis method that has been used in a variety of applications. The DPPC method is illustrated to mimic the correct molecular polarizability in a water dimer test case and in water‐halide ion complexes using the explicit polarization (X‐Pol) potential with the Austin model 1 Hamiltonian. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Structure and dynamics of the aqueous liquid-vapor interface: a comprehensive particle-based simulation study

This research addresses a comprehensive particle-based simulation study of the structural, dynamic, and electronic properties of the liquid-vapor interface of water utilizing both ab initio (based on density functional theory) and empirical (fixed charge and polarizable) models. Numerous properties such as interfacial width, hydrogen bond populations, dipole moments, and correlation times will be characterized with identical schemes to draw useful conclusions on the strengths and weakness of the proposed models for interfacial water. Our findings indicate that all models considered in this study yield similar results for the radial distribution functions, hydrogen bond populations, and orientational relaxation times. Significant differences in the models appear when examining both the dipole moments and surface relaxation near the aqueous liquid-vapor interface. Here, the ab initio interaction potential predicts a significant decrease in the molecular dipole moment and expansion in the oxygen-oxygen distance as one approaches the interface in accordance with recent experiments. All classical polarizable interaction potentials show a less dramatic drop in the molecular dipole moment, and all empirical interaction potentials studied yield an oxygen-oxygen contraction as the interface is approached.     

Geometry-dependent atomic charges: Methodology and application to alkanes,aldehydes, ketones,and amides

A general methodology for deriving geometry-dependent atomic charges is presented. The main ingredient of the method is a model that describes the molecular dipole moment in terms of geometry-dependent point charges. The parameters of the model are determined from ab initio calculations of molecular dipole moments and their Cartesian derivatives at various molecular geometries. Transferability of the parameters is built into the model by fitting ab initio calculations for various molecules simultaneously. The results show that charge flux along the bonds is a major contributing factor to the geometry dependence of the atomic charges, with additional contributions from fluxes along valence angles and adjacent bonds. Torsion flux is found to be smaller in magnitude than the bond and valence angle fluxes but is not always unimportant. A set of electrostatic parameters is presented for alkanes, aldehydes, ketones, and amides. Transferability of these parameters for a host of molecules is established to within 3 ?5% error in the predicted dipole moments. A possible extension of the method to include atomic dipoles is outlined. With the inclusion of such atomic dipoles and with the set of transferable point charges and charge flux parameters, it is demonstrated that molecular electrostatic potentials as well as electrostatic forces on nuclei can be reproduced much better than is possible with other models (such as potential derived charges). © 1995 by John Wiley & Sons, Inc.     

ABEEMσπ/MM力场模型下的多肽构象

在ABEEM/MM蛋白质浮动电荷力场模型的基础上,加入孤对电子和 电荷位点,从而能够体现多肽和蛋白质分子中一些重要的各向异性极化性质,允许非化学键方向的电子转移和极化。利用从头计算数据拟合模型相关参数。计算得到的小分子团簇结合能、偶极矩、氢键键长等性质与从头计算结果符合很好。该经典极化模型力场能够重复量子场下丙氨酸二肽、丙氨酸四肽、甘氨酸三肽的各稳定构象,其稳定性顺序与精密从头计算结果相一致,其结构和能量性质较以往模型有一定提高,并优于其他力场模型。     

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.     

Measurement of the electronic transition dipole moment by Autler-Townes splitting: Comparison of three- and four-level excitation schemes for the Na2 A 1Sigma(u)+ - X 1Sigma(g)+ system

We present a fundamentally new approach for measuring the transition dipole moment of molecular transitions, which combines the benefits of quantum interference effects, such as the Autler-Townes splitting, with the familiar R-centroid approximation. This method is superior to other experimental methods for determining the absolute value of the R-dependent electronic transition dipole moment function mu(e)(R), since it requires only an accurate measurement of the coupling laser electric field amplitude and the determination of the Rabi frequency from an Autler-Townes split fluorescence spectral line. We illustrate this method by measuring the transition dipole moment matrix element for the Na2 A 1Sigma(u)+ (v' = 25, J' = 20e)-X 1Sigma(g)+ (v" = 38, J" = 21e) rovibronic transition and compare our experimental results with our ab initio calculations. We have compared the three-level (cascade) and four-level (extended Lambda) excitation schemes and found that the latter is preferable in this case for two reasons. First, this excitation scheme takes advantage of the fact that the coupling field lower level is outside the thermal population range. As a result vibrational levels with larger wave function amplitudes at the outer turning point of vibration lead to larger transition dipole moment matrix elements and Rabi frequencies than those accessible from the equilibrium internuclear distance of the thermal population distribution. Second, the coupling laser can be "tuned" to different rovibronic transitions in order to determine the internuclear distance dependence of the electronic transition dipole moment function in the region of the R-centroid of each coupling laser transition. Thus the internuclear distance dependence of the transition moment function mu(e)(R) can be determined at several very different values of the R centroid. The measured transition dipole moment matrix element for the Na2 A 1Sigma(u)+ (v' = 25, J' = 20e)-X 1Sigma(g)+ (v" = 38, J" = 21e) transition is 5.5+/-0.2 D compared to our ab initio value of 5.9 D. By using the R-centroid approximation for this transition the corresponding experimental electronic transition dipole moment is 9.72 D at Rc = 4.81 A, in good agreement with our ab initio value of 10.55 D.     

Accuracy of distributed multipoles and polarizabilities: comparison between the LoProp and MpProp models

Localized multipole moments up to the fifth moment as well as localized dipole polarizabilities are calculated with the MpProp and the newly developed LoProp methods for a total of 20 molecules, predominantly derived from amino acids. A comparison of electrostatic potentials calculated from the multipole expansion obtained by the two methods with ab initio results shows that both methods reproduce the electrostatic interaction with an elementary charge with a mean absolute error of approximately 1.5 kJ/mol at contact distance and less than 0.1 kJ/mol at distances 2 A further out when terms up to the octupole moments are included. The polarizabilities are tested with homogenous electric fields and are found to have similar accuracy. The MpProp method gives better multipole moments unless diffuse basis sets are used, whereas LoProp gives better polarizabilities.     

Ab initio based polarizable force field parametrization

Experimental and simulation studies of anion-water systems have pointed out the importance of molecular polarization for many phenomena ranging from hydrogen-bond dynamics to water interfaces structure. The study of such systems at molecular level is usually made with classical molecular dynamics simulations. Structural and dynamical features are deeply influenced by molecular and ionic polarizability, which parametrization in classical force field has been an object of long-standing efforts. Although when classical models are compared to ab initio calculations at condensed phase, it is found that the water dipole moments are underestimated by approximately 30%, while the anion shows an overpolarization at short distances. A model for chloride-water polarizable interaction is parametrized here, making use of Car-Parrinello simulations at condensed phase. The results hint to an innovative approach in polarizable force fields development, based on ab initio simulations, which do not suffer for the mentioned drawbacks. The method is general and can be applied to the modeling of different systems ranging from biomolecular to solid state simulations.     

Ab initio rotation-vibration spectra of HCN and HNC

We have calculated an ab initio HCN/HNC linelist for all transitions up to J= 25 and 18000 cm(-1) above the zero point energy. This linelist contains more than 200 million lines each with frequencies and transition dipoles. The linelist has been calculated using our semi-global HCN/HNC VQZANO + PES and dipole moment surface, which were reported in van Mourik et al. (J. Chem. Phys. 115 (2001) 3706). With this linelist we synthesise absorption spectra of HCN and HNC at 298 K and we present the band centre and band transition dipoles for the bands which are major features in these spectra. Several of the HCN bands and many of the HNC bands have not been previously studied. Our line intensities reproduce via fully ab initio methods the unusual intensity structure of the HCN CN stretch fundamental (00(0)1) for the first time and also the forbidden (02(2)0) HCN bending overtone. We also compare the J = 1-->0 pure rotational transition dipole in the HCN/HNC ground and vibrationally excited states with experimental and existing ab initio results.     

Study of peptide conformation in terms of the ABEEM/MM method

The ABEEM/MM model (atom-bond electronegativity equalization method fused into molecular mechanics) is applied to study of the polypeptide conformations. The Lennard-Jones and torsional parameters were optimized to be consistent with the ABEEM/MM fluctuating charge electrostatic potential. The hydrogen bond was specially treated with an electrostatic fitting function. Molecular dipole moments, dimerization energies, and hydrogen bond lengths of complexes are reasonably achieved by our model, compared to ab initio results. The ABEEM/MM fluctuating charge model reproduces both the peptide conformational energies and structures with satisfactory accuracy with low computer cost. The transferability is tested by applying the parameters of our model to the tetrapeptide of alanine and another four dipeptides. The overall RMS deviations in conformational energies and key dihedral angles for four di- or tetrapeptide, is 0.39 kcal/mol and 7.7 degrees . The current results agree well with those by the accurate ab initio method, and are comparable to those from the best existing force fields. The results make us believe that our fluctuating charge model can obtain more promising results in protein and macromolecular modeling with good accuracy but less computer cost.     

Gaussian induced dipole polarization model

A new induced dipole polarization model based on interacting Gaussian charge densities is presented. In contrast to the original induced point dipole model, the Gaussian polarization model is capable of finite interactions at short distances. Aspects of convergence related to the Gaussian model will be explored. The Gaussian polarization model is compared with the damped Thole-induced dipole model and the point dipole model. It will be shown that the Gaussian polarization model performs slightly better than the Thole model in terms of fitting to molecular polarizability tensors. An advantage of the model based on Gaussian charge distribution is that it can be easily generalized to other multipole moments and provide effective damping for both permanent electrostatic and polarization models. Finally, a method of parameterizing polarizabilities is presented. This method is based on probing a molecule with point charges and fitting polarizabilities to electrostatic potential. In contrast to the generic atom type polarizabilities fit to molecular polarizability tensors, probed polarizabilities are significantly more accurate in terms of reproducing molecular polarizability tensors and electrostatic potential, while retaining conformational transferability.     

An interpretation of the enhancement of the water dipole moment due to the presence of other water molecules

The dipole moment of the gas phase water monomer is 1.85 D. When solvated in bulk water, the dipole moment of an individual water molecule is observed to be enhanced to the much larger value of 2.9 +/- 0.6 D. To understand the origin of this dipole moment enhancement, the effective fragment potential (EFP) method is used to solvate an ab initio water molecule to predict the dipole moments for various cluster sizes. The dipole moment as a function of cluster size, nH 2O, is investigated [for n = 6-20 (even n), 26, 32, 41, and 50]. Localized charge distributions are used in conjunction with localized molecular orbitals to interpret the dipole moment enhancement. These calculations suggest that the enhancement of the dipole moment originates from the decrease of the angle between the dipole vectors of the lone pairs on oxygen as the number of hydrogen bonds to that oxygen increases. Thus, the decreased angle, and the consequent increase in water dipole moment, is most likely to occur in environments with a larger number of hydrogen bonds, such as the center of a cluster of water molecules.     

Electrostatic calculation of linear and non-linear optical properties of ice Ih,II, IX and VIII

Macroscopic first- and third-order susceptibilities of ice Ih, ice II, ice IX and ice VIII are calculated using static and frequency-dependent electronic and static vibrational molecular (hyper)polarizabilities at the MP2 level. The molecular properties are in good agreement with experiment and with high-level ab initio calculations. Intermolecular electrostatic and polarization effects due to induced dipoles are taken into account using a rigorous local-field theory. The electric field due to permanent dipoles is used to calculate effective in-crystal (hyper)polarizabilities. The polarizability depends only weakly on the permanent field, but the dipole moment and the hyperpolarizabilities are strongly affected. The calculated linear susceptibility is in good agreement with available experimental data for ice Ih, and the third-order susceptibility for a third harmonic generation experiment is in reasonable agreement with experimental values for liquid water. The molecular vibrational contributions have a small effect on the susceptibilities. The electric properties of a water tetramer are calculated and used to estimate the effect of non-dipolar interactions on the susceptibilities of ice Ih, which are found to be small.     

Change in a protein's electronic structure induced by an explicit solvent: an ab initio fragment molecular orbital study of ubiquitin

The effect of solvation on the electronic structure of the ubiquitin protein was analyzed using the ab initio fragment molecular orbital (FMO) method. FMO calculations were performed for the protein in vacuo, and the protein was immersed in an explicit solvent shell as thick as 12 A at the HF or MP2 level by using the 6-31G* basis set. The protein's physical properties examined were the net charge, the dipole moment, the internal energy, and the solvent interaction energy. Comparison of the computational results revealed the following changes in the protein upon solvation. First, the positively charged amino acid residues on the protein surface drew electrons from the solvent, while the negatively charged ones transfer electrons to the solvent. Second, the dipole moment of the protein was enhanced as a result of the polarization. Third, the internal energy of the protein was destabilized, but the destabilization was more than compensated for by the generation of a favorable protein-solvent interaction. Finally, the energetic changes were elicited both by the electron correlation effect of the first solvent shell and by the electrostatic effect of more distant solvent molecules. These findings were consistent with the picture of the solvated protein being a polarizable molecule dissolved in a dielectric media.     

Reassessment of large dipole moment enhancements in crystals: a detailed experimental and theoretical charge density analysis of 2-methyl-4-nitroaniline

The molecular dipole moment of MNA in the crystal has been critically reexamined, to test the conclusion from an earlier experimental charge density analysis that it was substantially enhanced due to a combination of strong intermolecular interactions and crystal field effects. X-ray and neutron diffraction data have been carefully measured at 100 K and supplemented with ab initio crystal Hartree-Fock calculations. Considerable care taken in the measurement and reduction of the experimental data excluded most systematic errors, and sources of error and their effects on the experimental electron density have been carefully investigated. The electron density derived from a fit to theoretical structure factors assisted in the determination of the scale and thermal motion model. The dipole moment enhancement for MNA in the crystal is much less than that reported previously and only on the order of 30-40% (approximately 2.5 D). In addition to the dipole moment, experimental deformation electron density maps, bond critical point data, electric field gradients at hydrogen nuclei, and atomic and group charges all agree well with theoretical results and trends. Anisotropic modeling of the motion of hydrogen atoms, integral use of periodic ab initio calculations, and improved data quality are all aspects of this study that represent a considerable advance over previous work.     

Water liquid-vapor equilibria predicted by refined ab initio derived potentials

Coexistence properties for water near the critical point using several ab initio models were calculated using grand canonical Monte Carlo simulations with multiple histogram reweighting techniques. These models, that have proved to yield a good reproduction of the water properties at ambient conditions, perform rather well, improving the performance of a previous ab initio model. It is also shown that bulk geometry and dipole values, predicted by the simulation, can be used and a good approximation obtained with a polarizable rigid water model but not when polarization is excluded.     

Behavior of polarizable models in presence of strong electric fields. I. Origin of nonlinear effects in water point-charge systems

In the current opinion, the inclusion of polarization response in classical computer simulations is considered as one of the most important and urgent improvements to be implemented in modern empirical potential models. In this work we focus on the capability of polarizable models, based on the pairwise Coulomb interactions, to model systems where strong electric fields enter into play. As shown by Masia, Probst, and Rey (MPR) [in J. Chem. Phys. 121, 7362 (2004)], when a molecule interacts with point charges, polarizable models show underpolarization with respect to ab initio methods. We prove that this underpolarization, clearly related to nonlinear polarization effects, cannot be simply ascribed to the lack of hyperpolarization in the polarizable models, as suggested by MPR. Analysis of the electron-density rearrangement induced on a water molecule by a point charge reveals a twofold level of polarization response. One level involves intramolecular charge transfer on the whole molecular volume, with the related polarization exhibiting a seemingly linear behavior with the external electric field. The other nonlinear polarization level occurs only at strong electric fields and is found to be strictly correlated to the quantum-mechanical nature of the water molecule. The latter type of polarization has a local character, being limited to the space region of the water lone pairs.