The solvation of Li and Na in acetonitrile from ab initio-derived many-body ion–solvent potentials

Several Li⁺- and Na⁺-acetonitrile models were derived from ab initio calculations at the counterpoise-corrected MP2/TZV++(d,p) level for distorted ion-(MeCN)_n clusters with n=1, 4 and 6. Two different many-body ion-acetonitrile models were constructed: an effective three-body potential for use with the six-site effective pair model of Böhm et al., and an effective polarizable many-body model. The polarizable acetonitrile model used in the latter model is a new empirical model which was also derived in the present paper. Mainly for comparative purposes, two ion-acetonitrile pair potentials were also constructed from the ab initio cluster calculations: one pure pair potential and one effective pair potential. Using all these potential models, MD simulations in the NPT ensemble were performed for the pure acetonitrile liquid and for Li⁺(MeCN) and Na⁺(MeCN) solutions with 1 ion in 512 solvent molecules and with a simulation time of at least 120 ps per system. Thermodynamic properties, solvation-shell structure and the self-diffusion coefficient of the ions and of the solvent molecules were calculated and compared between the different models and with experimental data, where available. The Li⁺ ion is found to be four-coordinated when the new many-body potentials are used, in contrast to the six-coordinated structure obtained for the pure pair and effective pair potentials. The coordination number of Na⁺ is close to six for all the models derived here, although the coordination number becomes slightly smaller with the many-body potentials. For both ions, the solvent molecules in the first shell point their nitrogen ends towards the cation, while in the second shell the opposite orientation is the most common.     

Computer simulation of the supercritical carbon dioxide fluid (II) internal energy and structure of carbon dioxide system containing one benzene molecule

Using potential models based on ab initio quantum chemical calculations, we study a supercritical CO₂ fluid containing one benzene molecule using Monte Carlo simulations. First, molecular average internal energy is calculated for the whole system and for the first solvation shell of the benzene molecule. This analysis shows that the CO₂ molecules in the first solvation shell have a large energetic stabilization owing to the shape of the solute. In addition to the stabilization, the solute-solvent interactions in the first solvation shell show large fluctuations for both the in-plane and out-of-plane parts. Secondly, an orientational distribution function is defined to investigate the CO₂ fluid structure. This function indicates that the CO₂---CO₂ intermolecular configuration has a large dependence on the temperature of the system for both the whole system, and for the first solvation shell of the solute. Moreover, the benzene molecule is confirmed to control the mutual arrangement between neighboring CO₂ molecules.     

Theoretical study of the solvation of fluorine and chlorine anions by water

The solvation of fluoride and chloride anions (F(-) and Cl(-), respectively) by water has been studied using effective fragment potentials (EFPs) for the water molecules and ab initio quantum mechanics for the anions. In particular, the number of water molecules required to fully surround each anion has been investigated. Monte Carlo calculations have been used in an attempt to find the solvated system X(-)(H(2)O)(n) (X = F, Cl) with the lowest energy for each value of n. It is predicted that 18 water molecules are required to form a complete solvation shell around a Cl(-) anion, where "complete solvation" is interpreted as an ion that is completely surrounded by solvent molecules. Although fewer water molecules may fully solvate the Cl(-) anion, such structures are higher in energy than partially solvated molecules, up to n > or = 18. Calculations on the F(-) anion suggest that 15 water molecules are required for a complete solvation shell. The EFP predictions are in good agreement with the relative energies predicted by ab initio energy calculations at the EFP geometries.     

Monte Carlo simulation of a 1,4,7,10-tetraazacyclododecane/lithium complex in aqueous solution

Monte Carlo simulations have been carried out for the system consisting of a 1,4,7,10-tetraazacyclododecane (cyclen)-lithium complex in 201 water molecules. The volume of the periodic cube was calculated using the experimental density of pure water at 298 K and 1 atm of 1 g.cm(-)(3), plus additional space occupied by the complex. The geometry of the complex is the alternated form, where the ion is located at the center of the cyclen. The complex-water interaction was represented by the cyclen-water and lithium-water pair potentials, both of which were developed on the basis of ab initio calculations. The results show two layers of solvation shells consisting of 2 and 6.9 water molecules. Two water molecules in the first solvation shell (O(1) and O(2)) bind directly to the ion in which the ion-oxygen distance is 2.38 A, the dipole vector points to the ion, and rotation takes place around the ion-oxygen axis. In the next layer, 4 water molecules coordinate simultaneously to the first 2 water molecules in the first shell and the NH functional groups of cyclen. The remaining 2.9 water molecules in the second layer are also coordinated to be in the first half-hydration shell of O(1) and O(2).     

Probing the transition from hydrophilic to hydrophobic solvation with atomic scale resolution

Picosecond and femtosecond X-ray absorption spectroscopy is used to probe the changes of the solvent shell structure upon electron abstraction of aqueous iodide using an ultrashort laser pulse. The transient L(1,3) edge EXAFS at 50 ps time delay points to the formation of an expanded water cavity around the iodine atom, in good agreement with classical and quantum mechanical/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. These also show that while the hydrogen atoms pointed toward iodide, they predominantly point toward the bulk solvent in the case of iodine, suggesting a hydrophobic behavior. This is further confirmed by quantum chemical (QC) calculations of I(-)/I(0)(H(2)O)(n=1-4) clusters. The L(1) edge sub-picosecond spectra point to the existence of a transient species that is not present at 50 ps. The QC calculations and the QM/MM MD simulations identify this transient species as an I(0)(OH(2)) complex inside the cavity. The simulations show that upon electron abstraction most of the water molecules move away from iodine, while one comes closer to form the complex that lives for 3-4 ps. This time is governed by the reorganization of the main solvation shell, basically the time it takes for the water molecules to reform an H-bond network. Only then is the interaction with the solvation shell strong enough to pull the water molecule of the complex toward the bulk solvent. Overall, much of the behavior at early times is determined by the reorientational dynamics of water molecules and the formation of a complete network of hydrogen bonded molecules in the first solvation shell.     

Hydration of the Calcium Dication: Direct Evidence for Second Shell Formation from Infrared Spectroscopy

Infrared laser action spectroscopy in a Fourier‐transform ion cyclotron resonance mass spectrometer is used in conjunction with ab initio calculations to investigate doubly charged, hydrated clusters of calcium formed by electrospray ionization. Six water molecules coordinate directly to the calcium dication, whereas the seventh water molecule is incorporated into a second solvation shell. Spectral features indicate the presence of multiple structures of Ca(H₂O)₇²⁺ in which outer‐shell water molecules accept either one (single acceptor) or two (double acceptor) hydrogen bonds from inner‐shell water molecules. Double‐acceptor water molecules are predominately observed in the second solvent shells of clusters containing eight or nine water molecules. Increased hydration results in spectroscopic signatures consistent with additional second‐shell water molecules, particularly the appearance of inner‐shell water molecules that donate two hydrogen bonds (double donor) to the second solvent shell. This is the first reported use of infrared spectroscopy to investigate shell structure of a hydrated multiply charged cation in the gas phase and illustrates the effectiveness of this method to probe the structures of hydrated ions.     

MD studies of electrolyte solution|liquid mercury interfaces

We report molecular dynamics (MD) computer simulations of a single lithium or iodide ion near a water|liquid mercury interface. The ion–mercury and the water–mercury potentials are derived from ab initio calculations of an ion or a water molecule and a mercury cluster consisting of seven, nine or 10 atoms. The flexible BJH water model and a mercury–mercury potential derived from pseudopotential theory are employed. The ion–water potentials are also based on ab initio calculations. The structural properties at the interfaces are described in terms of various density profiles and the ion–mercury radial distribution functions (RDF). An analysis of the induced rearrangements of the mercury atoms at and below the surface is also performed. Finally, the spectral densities of the hindered translational motions of the ions parallel and perpendicular to the mercury surface are reported. We conclude that, while the I⁻-ion is contact adsorbed on the mercury surface, the Li⁺-ion is not.     

Redox entropy of plastocyanin: developing a microscopic view of mesoscopic polar solvation

We report applications of analytical formalisms and molecular dynamics (MD) simulations to the calculation of redox entropy of plastocyanin metalloprotein in aqueous solution. The goal of our analysis is to establish critical components of the theory required to describe polar solvation at the mesoscopic scale. The analytical techniques include a microscopic formalism based on structure factors of the solvent dipolar orientations and density and continuum dielectric theories. The microscopic theory employs the atomistic structure of the protein with force-field atomic charges and solvent structure factors obtained from separate MD simulations of the homogeneous solvent. The MD simulations provide linear response solvation free energies and reorganization energies of electron transfer in the temperature range of 280-310 K. We found that continuum models universally underestimate solvation entropies, and a more favorable agreement is reported between the microscopic calculations and MD simulations. The analysis of simulations also suggests that difficulties of extending standard formalisms to protein solvation are related to the inhomogeneous structure of the solvation shell at the protein-water interface combining islands of highly structured water around ionized residues along with partial dewetting of hydrophobic patches. Quantitative theories of electrostatic protein hydration need to incorporate realistic density profile of water at the protein-water interface.     

Semiempirical and ab initio calculations on the alkaline hydrolysis of the β-lactam ring Influence of the solvent

We used semiempirical and ab initio calculations to investigate the nucleophilic attack of the hydroxyl ion on the β-lactam carbonyl group. Both allowed us to detect reaction intermediates pertaining to proton-transfer reactions. We also used ab initio calculations and the PM3 semiempirical method to investigate the influence of the solvent on the process. The AMSOL method predicts the occurrence of a potential energy barrier of 20.7 kcal mol⁻¹ due to the desolvation of the hydroxyl ion in approaching the β-lactam carbonyl group. Using the supermolecular approach and a water solvation sphere of 20 molecules around the solute, the potential energy barrier is lowered to 17.5 kcal mol⁻¹. Ab initio calculations using the SCRF method predict a potential energy barrier of 13.6 kcal mol⁻¹. These three values, especially the last two, are very close to the experimental value of 16.7 kcal mol⁻¹.     

The coupling constant polarizability and hyperpolarizabilty of 1J(NH) in N-methylacetamide, and its application for the multipole spin-spin coupling constant polarizability/reaction field approach to solvation

We present a benchmark study of a combined multipole spin-spin coupling constant (SSCC) polarizability/reaction field (MJP/RF) approach to the calculation of both specific and bulk solvation effects on SSCCs of solvated molecules. The MJP/RF scheme is defined by an expansion of the SSCCs of the solvated molecule in terms of coupling constant dipole and quadrupole polarizabilities and hyperpolarizabilities derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics (MD) simulations thereby accounting for solute-solvent dynamical effects. The MJP/RF method is benchmarked against polarizable QM/MM calculations for the one-bond N-H coupling constant in N-methylacetamide. The best agreement between the MJP/RF and QM/MM approaches is found by truncating the electric field expansion in the MJP/RF approach at the linear electric field level. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the coupling constant (hyper-)polarizabilities and find that they are affected by the basis set in a way similar to the coupling constants themselves.     

Be2+在水、甲醇和乙醇中结构性质的从头算分子动力学模拟

采用Car-Parrinello分子动力学(CPMD)方法分别研究了Be²⁺在水、甲醇和乙醇中的溶剂结构性质, 并对Be²⁺的第一溶剂壳结构的实验及理论结果进行了比较. 所得第一溶剂壳结构与已报道的实验和理论结果较为一致. 对径向分布函数、配位数以及角度分布等进行了详细的分析. 结果表明: 在水、甲醇和乙醇中, Be²⁺第一溶剂壳为稳定理想的四面体结构. 在本文的模拟时间尺度内,没有观察到第一溶剂壳中的分子与第二溶剂壳中的分子进行交换, 进一步证明Be²⁺第一溶剂壳为稳定的四配位结构. 根据计算得到的空间分布函数, Be²⁺在溶剂分子的等高面上主要集中分布在溶剂分子接受氢键的方向. 根据氧原子在Be²⁺周围的分布, 壳层分子主要集中分布在Be²⁺周围的四个区域, 进一步证实溶剂壳为四面体对称.     

IR, Raman and theoretical ab initio RHF study of aminoglutethimide — an anticancer drug

A comparative analysis of the IR and Raman spectra of aminoglutethimide (AG) dissolved in CCl₄, CHCl₃ and CH₃CN was performed. Most of the absorption bands were assigned to characteristic group vibrations with the use of the IR and Raman spectra of deuterated AG, glutethimide, N-methyl glutethimide and glutarimide. The AG samples very weakly interacting with the environment were studied with the use of the Ar matrix isolation IR spectra. For comparison, the IR and Raman spectra of the crystalline samples formed by hydrogen-bonded AG molecules were recorded. The spectra were analyzed also in terms of normal modes and the harmonic approximation with the use of the ab initio restricted Hartree-Fock theory. It was found that increasing the solute concentration in CCl₄ and CHCl₃ leads to formation of the autoassociates. In CH₃CN the solute–solvent AG–CH₃CN dimers occur. Possible structures of the associates were theoretically studied on the model systems: the centrosymmetric glutarimide dimer and the linear AG–CH₃CN dimer. By a comparison of the theoretical and experimental spectra we were able to identify several peaks originating from the solute–solvent interactions.     

First principles and classical molecular dynamics simulations of solvated benzene

We have performed extensive ab initio and classical molecular dynamics (MD) simulations of benzene in water in order to examine the unique solvation structures that are formed. Qualitative differences between classical and ab initio MD simulations are found and the importance of various technical simulation parameters is examined. Our comparison indicates that nonpolarizable classical models are not capable of describing the solute-water interface correctly if local interactions become energetically comparable to water hydrogen bonds. In addition, a comparison is made between a rigid water model and fully flexible water within ab initio MD simulations which shows that both models agree qualitatively for this challenging system.     

Ab initio QM/MM simulation of Ag+ in 18.6% aqueous ammonia solution: structure and dynamics investigations

Structure and dynamics investigations of Ag(+) in 18.6% aqueous ammonia solution have been carried out by means of the ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation method. The most important region, the first solvation shell, was treated by ab initio quantum mechanics at the Restricted Hartree-Fock (RHF) level using double-zeta plus polarization basis sets for ammonia and plus ECP for Ag(+). For the remaining region in the system, newly constructed three-body corrected potential functions were used. The average composition of the first solvation shell was found to be [Ag(NH(3))(2)(H(2)O)(2.8)](+). No ammonia exchange process was observed for the first solvation shell, whereas ligand exchange processes occurred with a very short mean residence time of 1.1 ps for the water ligands. No distinct second solvation shell was observed in this simulation.     

Modeling water exchange on an aluminum polyoxocation

For the first time, water exchange on a polymeric complex has been modeled using a combination of gas-phase ab initio calculations and molecular dynamics (MD) simulations. The GaO4Al12(OH)24(H2O)12(7+)aq ion (GaAl12) was chosen because high-quality experimental data exist, including an activation enthalpy (+63 +/- 7 kJ/mol) and an activation volume (+3 +/- 1 cm3/mol). We took a two-step approach. First, the local solvent structure and the initial states for reaction were inferred from the molecular dynamics simulations. Second, we used this information to evaluate initial-state structures in the ab initio calculations. The energy differences between the initial and transition states from the ab initio calculations varied from +59 kJ/mol to +53 kJ/mol depending upon details, closely approximating the activation enthalpy.     

Electronic structure and reactivity of Mg+(H2O)n cluster ions

Electronically excited states of magnesium-water cluster ions, Mg⁺(H₂O) _n ,n=1–5, are studied by photodissociation after mass selection. The observed photodissociation spectra are assigned to the²P–²S type transitions localized on the Mg⁺ ion with the aid of ab initio CI calculations. In addition to evaporation of water molecules, photoinduced intracluster reaction to produce MgOH⁺(H₂O) _n is found to occur efficiently, with pronounced size dependence. The intriguing features observed in the mass spectrum of nascent cluster ions are discussed in relation to the stepwise solvation of this reaction.     

Structure and dynamics of hydrated NH(4) (+): an ab initio QM/MM molecular dynamics simulation

A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to investigate solvation structure and dynamics of NH(4) (+) in water. The most interesting region, the sphere includes an ammonium ion and its first hydration shell, was treated at the Hartree-Fock level using DZV basis set, while the rest of the system was described by classical pair potentials. On the basis of detailed QM/MM simulation results, the solvation structure of NH(4) (+) is rather flexible, in which many water molecules are cooperatively involved in the solvation shell of the ion. Of particular interest, the QM/MM results show fast translation and rotation of NH(4) (+) in water. This phenomenon has resulted from multiple coordination, which drives the NH(4) (+) to translate and rotate quite freely within its surrounding water molecules. In addition, a "structure-breaking" behavior of the NH(4) (+) is well reflected by the detailed analysis on the water exchange process and the mean residence times of water molecules surrounding the ion.     

Theoretical study of C24N4 molecule

Both ab initio 6-31G, 3-21G and STO-3G basis sets and semiempirical PM3 and AM1 molecular orbital calculations are carried out on the C₂₄N₄ molecule of the T_d symmetry group. Results on the fully optimized structure which constrained T_d symmetry, molecular orbitals and vibrational frequency were obtained by both ab initio and semiempirical methods. The binding energy and various thermodynamic properties were also calculated via the PM3 and AM1 semiempirical methods. All the evidence of this work proves that the C₂₄N₄ molecule is stable and that its four six-membered rings with a remarkable delocalized C…C bond are similar to the related rings in the C₆₀ buckminsterfullerene structure.     

Cu(II) in liquid ammonia: an approach by hybrid quantum-mechanical/molecular-mechanical molecular dynamics simulation.

To investigate the solvation structure of the Cu(II) ion in liquid ammonia, ab initio quantum-mechanical/molecular-mechanical (QM/MM) molecular dynamics (MD) simulations were carried out at Hartree Fock (HF) and hybrid density functional theory (B3 LYP) levels. A sixfold-coordinated species was found to be predominant in the HF case whereas five- and sixfold-coordinated complexes were obtained in a ratio 2:1 from the B3 LYP simulation. In contrast to hydrated Cu(II), which exhibits a typical Jahn-Teller distortion, the geometrical arrangement of ligand molecules in the case of ammonia can be described as a [2 + 4] ([2 + 3]) configuration with 4 (3) elongated copper-nitrogen bonds. First shell solvent exchange reactions at picosecond rate took place in both HF and B3 LYP simulations, again in contrast to the more stable sixfold-coordinated hydrate. NH3 ligands apparently lead to strongly accelerated dynamics of the Cu(II) solvate due to the "inverse" [2 + 4] structure with its larger number of elongated copper-ligand bonds. Several dynamical properties, such as mean ligand residence times or ion-ligand stretching frequencies, prove the high lability of the solvated complex.