Development of first-principles interaction model potentials. An application to the study of the bromide hydration

Probing the geometry reorganization from solution to gas-phase in putrescine derivatives by IRMPD, 1H-NMR and theoretical calculations

Geometry reorganization of ESI formed ions are demonstrated and explicit calculations of the solution phase are shown to be relevant.

Force field development for actinyl ions via quantum mechanical calculations: an approach to account for many body solvation effects

Advances in computational algorithms and methodologies make it possible to use highly accurate quantum mechanical calculations to develop force fields (pair-wise additive intermolecular potentials) for condensed phase simulations. Despite these advances, this approach faces numerous hurdles for the case of actinyl ions, AcO2(n+) (high-oxidation-state actinide dioxo cations), mainly due to the complex electronic structure resulting from an interplay of s, p, d, and f valence orbitals. Traditional methods use a pair of molecules (“dimer”) to generate a potential energy surface (PES) for force field parametrization based on the assumption that many body polarization effects are negligible. We show that this is a poor approximation for aqueous phase uranyl ions and present an alternative approach for the development of actinyl ion force fields that includes important many body solvation effects. Force fields are developed for the UO2(2+) ion with the SPC/Fw, TIP3P, TIP4P, and TIP5P water models and are validated by carrying out detailed molecular simulations on the uranyl aqua ion, one of the most characterized actinide systems. It is shown that the force fields faithfully reproduce available experimental structural data and hydration free energies. Failure to account for solvation effects when generating PES leads to overbinding between UO2(2+) and water, resulting in incorrect hydration free energies and coordination numbers. A detailed analysis of arrangement of water molecules in the first and second solvation shell of UO2(2+) is presented. The use of a simple functional form involving the sum of Lennard-Jones + Coulomb potentials makes the new force field compatible with a large number of available molecular simulation engines and common force fields.

Theoretical study of the aqueous solvation of HgCl2: Monte Carlo simulations using second-order Moller-Plesset-derived flexible polarizable interaction potentials

A study of the solvation of HgCl(2) including ab initio aggregates of up to 24 water molecules and the results of extensive Monte Carlo simulations for the liquid phase using MP2-derived interaction potentials is presented. The interaction potentials are flexible, polarizable, and include non-additive effects. We conclude that a cluster description of the solvation mechanism is limited when compared to the condensed phase. The molecular image derived from the MC simulations is peculiar. It resembles that of a hydrophobic solute, which explains the rather easy passage of this neutral molecule through the cell membrane; however, it also shows an intermittent binding of one, two, or three water molecules to HgCl(2) in the fashion of a hydrophilic solute.

A refined potential for hydroxylamine clusters and the liquid phase

A detailed study including ab initio calculations and classic Monte-Carlo simulations of hydroxylamine in the gas and liquid phases is presented. A classical interaction potential for hydroxylamine, which includes polarizability, many-body effects, and intramolecular relaxation, was constructed. The results of the simulation were compared to the available experimental data in order to validate the model. We conclude that liquid hydroxylamine has a multitude of hydrogen bonds leading to a large density where the existence of cis conformers and clusters of these conformers is possible. This explains the occurrence of the classical [R. Nast and I. Z. Foppl, Z. Anorg. Allg. Chem. 263, 310 (1950)] scheme for the molecule's decomposition at room temperature and its large exothermicity and instability.

Theoretical studies on photoelectron and IR spectral properties of Br2∙−(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.

Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field

A polarizable potential function for the hydration of alkali and halide ions is developed on the basis of the recent SWM4-DP water model [Lamoureux, G.; MacKerell, A. D., Jr.; Roux, B. J. Chem. Phys. 2003, 119, 5185]. Induced polarization is incorporated using classical Drude oscillators that are treated as auxiliary dynamical degrees of freedom. The ions are represented as polarizable Lennard-Jones centers, whose parameters are optimized to reproduce the binding energies of gas-phase monohydrates and the hydration free energies in the bulk liquid. Systematic exploration of the parameters shows that the monohydrate binding energies can be consistent with a unique hydration free energy scale if the computed hydration free energies incorporate the contribution from the air/water interfacial electrostatic potential (-540 mV for SWM4-DP). The final model, which can satisfyingly reproduce both gas and bulk-phase properties, corresponds to an absolute scale in which the intrinsic hydration free energy of the proton is -247 kcal/mol.

Structure, energy, and IR spectra of I2∙−.nH2O clusters (n=1–8): A theoretical study

The authors report theoretical results on structure, bonding, energy, and infrared spectra of iodine dimer radical anion hydrated clusters, I(2) (-).nH(2)O (n=1-8), based on a systematic study following density functional theory. Several initial guess structures are considered for each size cluster to locate minimum energy conformers with a Gaussian 6-311++G(d,p) split valence basis function (triple split valence 6-311 basis set is applied for iodine). It is observed that three different types of hydrogen bonded structures, namely, symmetrical double hydrogen bonding, single hydrogen bonding, and interwater hydrogen bonding structures, are possible in these hydrated clusters. But conformers having interwater hydrogen bonding arrangements are more stable compared to those of double or single hydrogen bonded structures. It is also noticed that up to four solvent H(2)O units can reside around the solute in interwater hydrogen bonding network. At the maximum six H(2)O units are independently linked to the dimer anion having four double hydrogen bonding and two single hydrogen bonding, suggesting the hydration number of I(2) (-) to be 6. However, conformers having H(2)O units independently linked to the iodine dimer anion are not the most stable structures. In all these hydrated clusters, the odd electron is found to be localized over two I atoms and the two atoms are bound by a three-electron hemi bond. The solvation, interaction, and vertical detachment energies are calculated for all I(2) (-).nH(2)O clusters. Energy of interaction and vertical detachment energy profiles show stepwise saturation, indicating geometrical shell closing in the hydrated clusters, but solvation energy profile fails to show such behavior. A linear correlation is observed between the calculated energy of interaction and vertical detachment energy. It is observed that formation of I(2) (-)-water cluster induces significant shifts from the normal O-H stretching modes of isolated H(2)O. However, bending mode of H(2)O remains insensitive to the successive addition of solvent H(2)O units. Weighted average energy profiles and IR spectra are reported for all the hydrated clusters based on the statistical population of individual conformers at room temperature.

On the performance of molecular polarization methods. I. Water and carbon tetrachloride close to a point charge

The three main methods to implement molecular polarization (point dipoles, fluctuating charges, and shell model) are tested against high level ab initio calculations for a molecule (water, carbon tetrachloride) close to a point charge (at the distance of a lithium or magnesium ion). The goal is to check whether an approximation (linear polarization) strictly valid at large intermolecular distances is sufficiently accurate for liquid state molecular dynamics simulations, where strong polarization effects are to be expected at short separations. The monitored observable is the molecular dipole moment as a function of the charge-molecule distance for selected molecular orientations. Analytic formulas are derived for the components of the molecular polarization tensor, facilitating the optimization of the performance for each polarization method as a function of its underlying parameters. Overall, the methods studied provide a remarkably good representation of the induced dipole, with no divergences appearing even at the shortest distances. For water close to a monovalent point charge the point dipole model, implemented with one or three dipoles, accurately reproduces the water dipole moment at all distances. Deficiencies appear as the molecular polarizability and/or charge increase: basically, the ab initio induced moments grow faster at intermediate distances than the linear increase characteristic of the phenomenological polarization methods, suggesting that nonlinear effects (hyperpolarizability) cannot be neglected in these cases. Regarding the capabilities of each method, the point dipole method is the one that performs best overall, with the shell model achieving acceptable results in most instances. The fluctuating charge method shows some noticeable limitations for implementations of comparable complexity (in terms of the number of sites required).

Microhydration shell structure in Cl2∙−∙nH2O clusters: A theoretical study

We present the results of a detailed study on structure and electronic properties of hydrated cluster Cl2*-.nH2O (n = 1-7) based on a nonlocal density functional, namely, Becke's [J. Chem. Phys. 98, 1372 (1993)] half and half hybrid exchange-correlation functional with a split valence 6-311++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with various possible initial guess structures without any symmetry restriction. Several minimum energy structures (conformers) are predicted with a small difference in total energy. There is a competition between the binding of solvent H2O units with Cl2*- dimer radical anion directly through ion-molecule interaction and forming interwater hydrogen-bonding network in Cl2*-.nH2O (n > or = 2) hydrated cluster. Structure having interwater H-bonded network is more stable over the structure where H2O units are connected to the solute dimer radical anion Cl2*- rather independently either by single or double H bonding in a particular size (n) of hydrated cluster Cl2*-.nH2O. At the maximum four solvent H2O units reside in interwater H-bonding network present in these hydrated clusters. It is observed that up to six H2O units are independently linked to the anion having four double H bondings and two single H bondings suggesting the primary hydration number of Cl2*- to be 6. In all these clusters, the odd electron is found to be mostly localized over the two Cl atoms and these two atoms are bound by a three-electron hemibond. Calculated interaction (between solute and different water clusters) and vertical detachment energy profiles show saturation at n = 6 in the hydrated cluster Cl2*-.nH2O (n = 1-7). However, calculated solvation energy increases with the increase in number of solvent H2O molecules in the cluster. Interaction energy varies linearly with vertical detachment energy for the hydrated clusters Cl2*-.nH2O (n < or = 6). Calculation of the vibration frequencies show that the formation of Cl2*(-)-water clusters induces significant shifts from the normal stretching modes of isolated water. A clear difference in the pattern of IR spectra is observed in the O-H stretching region of water from hexa- to heptahydrated cluster.

Density functional study of ion hydration for the alkali metal ions (Li+,Na+,K+) and the halide ions (F−,Br−,Cl−)

We performed first principles density functional calculations to study the effect of monovalent ions M+ (M = Li,Na,K) and A- (A = F,Cl,Br) in water with the aim of characterizing the local molecular properties of hydration. For this reason, several ion-water clusters, up to five or six water molecules were considered; such structures were optimized, and the Wannier analysis was then applied to determine the average molecular dipole moment of water. We found that with an increasing number of water molecules, the molecular polarization is determined by the water-water interaction rather than the water-ion interaction, as one would intuitively expect. These results are consistent with those obtained in previous density functional calculations and with other results obtained by employing classical polarizable water models. The main message of this work is that as one increases the number of molecules the average dipole moment of all water molecules and the ones in the first shell tends to the same value as the average of a similar sized cluster of pure water. This supports the use of nonpolarizable classical models of water in classical atomistic simulations.

Structure and Dynamics of Br- Ion in Liquid Methanol

A Car-Parrinello molecular dynamics simulation has been performed on a solution of Br- in liquid methanol analyzing with particular attention charge transfer and polarization effects. The first solvation shell has been characterized in terms of H-bonds, and it has been found that the high polarization of the bromide gives rise to a stable solvent cage. The differences in the coordination number with the chloride can be ascribed to the ionic radius and to the stronger perturbations brought by the solvent to the bromide ion.

A comparative ab initio study of Br2•− and Br2 water clusters

The work presents ab initio results on structure and electronic properties of Br2*-.nH2O(n=1-10) and Br2.nH2O(n=1-8) hydrated clusters to study the effects of an excess electron on the microhydration of the halide dimer. A nonlocal density functional, namely, Becke's half-and-half hybrid exchange-correlation functional is found to perform well on the present systems with a split valence 6-31++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with several initial guess structures and without imposing any symmetry restriction. Br2*-.nH2O clusters prefer to have symmetrical double hydrogen-bonding structures. Results on Br2.nH2O(n>or=2) cluster show that the O atom of one H2O is oriented towards one Br atom and the H atom of another H2O is directed to other Br atom making Br2 to exist as Br+-Br- entity in the cluster. The binding and solvation energies are calculated for the Br2*-.nH2O and Br2.nH2O clusters. Calculations of the vibrational frequencies show that the formation of Br2*- and Br2 water clusters induces significant shifts from the normal modes of isolated water. Excited-state calculations are carried out on Br2*-.nH2O clusters following configuration interaction with single electron excitation procedure and UV-VIS absorption profiles are simulated. There is an excellent agreement between the present theoretical UV-VIS spectra of Br2*-.10H2O cluster and the reported transient optical spectra for Br2*- in aqueous solution.

Dependence of ion hydration on the sign of the ion’s charge

The solvation of simple ions in water is studied using molecular dynamics simulations with a polarizable force field. Previous simulations using this potential demonstrated that anions are more favorably solvated in water than cations. The present work is an attempt to explain this result by examining the effects of ions on the surrounding water structure, with particular focus on the first solvation shell and its interactions with the surrounding water. We conclude that while the first solvation shell surrounding cations is frustrated by competition between ion-water and water-water interactions, solvation of anions is compatible with good water-water interactions.

Study of the stabilization energies of halide-water clusters: An application of first-principles interaction potentials based on a polarizable and flexible model

The aim of this work is to compute the stabilization energy E(stab)(n) of [X(H(2)O)(n)](-) (X identical with F, Br, and I for n=1-60) clusters from Monte Carlo simulations using first-principles ab initio potentials. Stabilization energy of [X(H(2)O)(n)](-) clusters is defined as the difference between the vertical photodeachment energy of the cluster and the electron affinity of the isolated halide. On one hand, a study about the relation between cluster structure and the E(stab)(n) value, as well as the dependence of the latter with temperature is performed, on the other hand, a test on the reliability of our recently developed first-principles halide ion-water interaction potentials is carried out. Two different approximations were applied: (1) the Koopmans' theorem and (2) calculation of the difference between the interaction energy of [X(H(2)O)(n)](-) and [X(H(2)O)(n)] clusters using the same ab initio interaction potentials. The developed methodology allows for using the same interaction potentials in the case of the ionic and neutral clusters with the proviso that the charge of the halide anion was switched off in the latter. That is, no specific parametrization of the interaction potentials to fit the magnitude under study was done. The good agreement between our predicted E(stab)(n) and experimental data allows us to validate the first-principles interaction potentials developed elsewhere and used in this study, and supports the fact that this magnitude is mainly determined by electrostatic factors, which can be described by our interaction potentials. No relation between the value of E(stab)(n) and the structure of clusters has been found. The diversity of E(stab)(n) values found for different clusters with similar interaction energy indicates the need for statistical information to properly estimate the stabilization energy of the halide anions. The effect of temperature in the prediction of the E(stab)(n) is not significant as long as it was high enough to avoid cluster trapping into local equilibrium configurations which guarantees an appropriate sampling of the configurational space. Parallel tempering method was applied in particular cases to guarantee satisfactory sampling of clusters at low temperature.

Many-body potentials for aqueous Li+, Na+, Mg2+, and Al3+: Comparison of effective three-body potentials and polarizable models

Many-body potentials for the aqueous Li(+), Na(+), Mg(2+), and Al(3+) ions have been constructed from ab initio cluster calculations. Pure pair, effective pair, effective three-body, and effective polarizable models were created and used in subsequent molecular dynamics simulations. The structures of the first and second solvation shells were studied using radial distribution functions and angular-radial distribution functions. The effective three-body and polarizable potentials yield similar first-shell structures, while the contraction of the O-O distances between the first and second solvation shells is more pronounced with the polarizable potentials. The definition of the tilt angle of the water molecules around the ions is discussed. When a proper definition is used, it is found that for Li(+), Mg(2+), and Al(3+) the water molecules prefer a trigonal orientation, but for Na(+) a tetrahedral orientation (ion in lone-pair direction) is preferred. The self-diffusion coefficients for the water molecules and the ions were calculated; the ionic values follow the order obtained from experiment, although the simulated absolute values are smaller than experiment for Mg(2+) and Al(3+).