Characterization of dynamics and reactivities of solvated ions byab initio simulations

Hydration of iron-porphyrins: ab initio quantum mechanical charge field molecular dynamics simulation study

The ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) simulation approach was successfully applied to Fe²⁺-P and Fe³⁺-P in water to evaluate their structural, dynamical and energetic properties. Based on the structural data, it was found that Fe²⁺-P accommodates one water molecule in the first coordination sphere of the Fe²⁺ ion including the four nitrogen atoms of the porphyrin system coordinating with central metal species. On the other hand, two water molecules were coordinated to Fe³⁺-P, thus forming a hexa-coordinated species. Comparison of dynamical properties such as the vibrational power spectrum and ligand mean residence times to other metal-free porphyrin systems demonstrate the ions' influence on the hydration structure, enabling a characterisation of the strong interaction of the ions which greatly reduces the hydrogen bonding potential of the complex. The association of water molecules with the metal ions in both solutes was quantified by computing the free energy of binding obtained via the potential of mean force. This further confirmed the strong association of water to the metal ions which was conversely weak as inferred from the energetic data for the Fe²⁺-P system.

Zinc- and copper-porphyrins in aqueous solution - two similar complexes with strongly contrasting hydration

We present detailed analysis of the hydration behavior of zinc and copper bound porphyrins treated via ab initio quantum mechanical charge field molecular dynamics which agrees well with available experimental data. The computed metal-water coordination in the case of zinc bound porphyrin demonstrates a strong association of water with zinc compared to the copper-water interaction which correlates well with the calculated free energy of binding. The H-bond dynamics in these hydrated systems yield weaker H-bond interactions as compared to that observed in the case of metal-free porphyrin; nevertheless, the effect of metal association with porphyrin resulted in shifts in the vibrational frequencies. These characteristic data suggest a contrasting behavior between these metalloporphyrins in solution which could also serve to correlate with the properties of biological systems.

Palladium(ii) in liquid ammonia: an investigation of structural and dynamical properties by applying quantum mechanical charge field molecular dynamics (QMCF-MD)

The structural and dynamical properties of Pd²⁺ in liquid ammonia have been investigated via quantum mechanical charge field molecular dynamics.

The hydration structure of the heavy-alkalines Rb+ and Cs+ through molecular dynamics and X-ray absorption spectroscopy: surface clusters and eccentricity

Hydration shells around Rb⁺ and Cs⁺ are not symmetric; the cation and the 1st-shell water mass center are separated by ∼0.4 Å, and this is supported by agreement between the theoretical and experimental EXAFS spectrum.

Nuclear quantum effects in water exchange around lithium and fluoride ions

We employ classical and ring polymer molecular dynamics simulations to study the effect of nuclear quantum fluctuations on the structure and the water exchange dynamics of aqueous solutions of lithium and fluoride ions. While we obtain reasonably good agreement with experimental data for solutions of lithium by augmenting the Coulombic interactions between the ion and the water molecules with a standard Lennard-Jones ion-oxygen potential, the same is not true for solutions of fluoride, for which we find that a potential with a softer repulsive wall gives much better agreement. A small degree of destabilization of the first hydration shell is found in quantum simulations of both ions when compared with classical simulations, with the shell becoming less sharply defined and the mean residence time of the water molecules in the shell decreasing. In line with these modest differences, we find that the mechanisms of the exchange processes are unaffected by quantization, so a classical description of these reactions gives qualitatively correct and quantitatively reasonable results. We also find that the quantum effects in solutions of lithium are larger than in solutions of fluoride. This is partly due to the stronger interaction of lithium with water molecules, partly due to the lighter mass of lithium and partly due to competing quantum effects in the hydration of fluoride, which are absent in the hydration of lithium.

A Comparative Study of [CaEDTA](2-) and [MgEDTA](2-): Structural and Dynamical Insights from Quantum Mechanical Charge Field Molecular Dynamics

Structure and dynamics of [MgEDTA](2-) and [CaEDTA](2-) complexes in aqueous solution have been investigated via quantum mechanical/molecular mechanical (QM/MM) simulations. While for the first a 6-fold octahedral complex has been observed, the presence of an additional coordinating water ligand has been observed in the latter case. Because of rapidly exchanging water molecules, this 7-fold coordination complex was found to form pentagonal bipyramidal as well as capped trigonal prismatic configurations along the simulation interchanging on the picosecond time scale. Also in the case of [MgEDTA](2-) a trigonal prismatic configuration has been observed for a very short time period of approximately 1 ps. This work reports for the first time the presence of trigonal prismatic structures observed in the coordination sphere of [MgEDTA](2-) and [CaEDTA](2-) complexes in aqueous solution. In addition to the detailed characterization of structure and dynamics of the systems, the prediction of the associated infrared spectra indicates that the ion-water vibrational mode found at approximately 250 cm(-1) provides a distinctive measure to experimentally detect the presence of the coordinating water molecule via low-frequency IR setups.

Hydration properties determining the reactivity of nitrite in aqueous solution

The knowledge of the hydration properties of the nitrite ion is key to understanding its reaction mechanism controlled by solvent effects. Here, ab initio quantum mechanical charge field molecular dynamics was performed to obtain the structural and dynamical properties of the hydration shell in an aqueous solution of nitrite ions, elucidated by data analysis using a molecular approach and an extended quantitative analysis of all superimposed trajectories with three-dimensional alignment (density map). The pattern of the power spectra corresponded to the experimental data, indicating the suitability of the Hartree-Fock method coupled with double-ζ plus polarization and diffuse functional basis sets to study this system. The density maps revealed the structure of the hydration shell, that presented a higher density in the N-O bond direction than in the axis vertical to the molecular plane, whereas the atomic and molecular radial distribution functions provided vague information. The number of actual contacts indicated 4.6 water molecules interacting with a nitrite ion, and 1.5 extra water molecules located in the molecular hydration shell, forming a H-bonding network with the bulk water. The mean residence times for the water ligands designated the strength of the hydration spheres for the oxygen sites, whilst the results for the nitrogen sites over-estimated the number of water molecules from other sites and indicated a weak structure. These results show the influence of the water molecules surrounding the nitrite ion creating an anisotropic hydration shell, suggesting that the reactive sites are situated above and below the molecular plane with a lower water density.

Modelling the effects of salt solutions on the hydration of calcium ions

Classical molecular dynamics simulations of several aqueous alkali halide salt solutions have been used to determine the effect of electrolytes on the structure of water and the hydration properties of calcium ions. Compared with the simulations of Ca(2+) ions in pure liquid water, the frequency of water exchange in the first hydration shell of calcium, which is a fundamental process in controlling the reactivity of calcium(ii) aqua-ions, is drastically reduced in the presence of other electrolytes in solution. The strong stabilization of the hydration shell of Ca(2+) occurs not only when the halide anions are directly coordinated to calcium, but also when the alkali and halide ions are placed at or outside the second coordination shell of Ca(2+), suggesting that the reactivity of the first solvation shell of the calcium ion can be influenced by the specific affinity of other ions in solution for the water molecules coordinated to Ca(2+). Analysis of the hydrogen-bonded structure of water in the vicinity of the calcium ion shows that the average number of hydrogen bonds per water molecules, which is 1.8 in pure liquid water, decreases as the concentration of alkali-halide salts in solution increases, and that the temporal fluctuations of hydrogen bonds are significantly larger than those obtained for Ca(2+) in pure liquid water. This effect has been explained in terms of the dynamics of reorganization of the O-H···X(-) (X = F, Cl and Br) hydrogen bond. This work shows the importance of solution composition in determining the hydrogen-bonding network and ligand-exchange dynamics around metal ions, both in solution and at the mineral-water interfaces, which in turn has implications for interactions occurring at the mineral-water interface, ultimately controlling the mobilization of ions in the environment as well as in industrial processes.

Structure and water exchange dynamics of hydrated oxo halo ions in aqueous solution using QMCF MD simulation, large angle X-ray scattering and EXAFS

Theoretical ab initio quantum mechanical charge field molecular dynamics (QMCF MD) has been applied in conjunction with experimental large angle X-ray scattering (LAXS) and EXAFS measurements to study structure and dynamics of the hydrated oxo chloro anions chlorite, ClO2(-), chlorate, ClO3(-), and perchlorate, ClO4(-). In addition, the structures of the hydrated hypochlorite, ClO(-), bromate, BrO3(-), iodate, IO3(-) and metaperiodate, IO4(-), ions have been determined in aqueous solution by means of LAXS. The structures of the bromate, metaperiodate, and orthoperiodate, H2IO6(3-), ions have been determined by EXAFS as solid sodium salts and in aqueous solution as well. The results show clearly that the only form of periodate present in aqueous solution is metaperiodate. The Cl-O bond distances in the hydrated oxo chloro anions as determined by LAXS and obtained in the QMCF MD simulations are in excellent agreement, being 0.01-0.02 Å longer than in solid anhydrous salts due to hydration through hydrogen bonding to water molecules. The oxo halo anions, all with unit negative charge, have low charge density making them typical structure breakers, thus the hydrogen bonds formed to the hydrating water molecules are weaker and more short-lived than those between water molecules in pure water. The water exchange mechanism of the oxo chloro anions resembles those of the oxo sulfur anions with a direct exchange at the oxygen atoms for perchlorate and sulfate. The water exchange rate for the perchlorate ion is significantly faster, τ0.5 = 1.4 ps, compared to the hydrated sulfate ion and pure water, τ0.5 = 2.6 and 1.7 ps, respectively. The angular radial distribution functions show that the chlorate and sulfite ions have a more complex water exchange mechanism. As the chlorite and chlorate ions are more weakly hydrated than the sulfite ion the spatial occupancy is less well-defined and it is not possible to follow any well-defined migration pattern as it is difficult to distinguish between hydrating water molecules and bulk water in the region close to the ions.

QM/MM MD simulations of La(iii)–phosphopeptide complexes

Hybrid quantum mechanical/molecular mechanical simulations have been used to study the structural and dynamical properties of a La(iii)–phosphopeptide complex.

Structure and dynamics of chromatographically relevant Fe(III)-chelates

Immobilized metal ion affinity chromatography (IMAC) is an important chromatographic technique for biomolecules. In order to get a detailed understanding of the hydration of immobilized Fe(III), complexes of Fe(III) with methyl substituted iminodiacetate ([Fe(MSIDA)(H2O)3](+)) as well as with methyl substituted nitrilotriacetate ([Fe(MSNTA)(H2O)2]) were simulated in aqueous solutions with the quantum mechanical charge field molecular dynamics (QMCF MD) approach. The simulations were carried out at the Hartree-Fock (HF) level of theory, since cluster calculations at the HF, MP2, and B3LYP levels of theory showed that this method results in a good compromise between computational effort and accuracy. None of the coordinating water molecules were exchanged during the simulation period of 15 ps. The Fe-OH2O bond distances as well as the Fe-OH2O stretching motions differed among the coordinating water molecules, indicating different bond strengths. For the water molecules in the second hydration layer, mean residence times of 2.7 and 1.9 ps were obtained for [Fe(MSIDA)(H2O)3](+) and [Fe(MSNTA)(H2O)2], respectively. Furthermore, infrared measurements were carried out to characterize the most prominent bond features of aqueous Fe(III)-NTA and to discuss these results in conjunction with the computationally derived frequencies.

Structure and water exchange of the hydrated thiosulfate ion in aqueous solution using QMCF MD simulation and large angle X-ray scattering

Experimental and simulation data of the thiosulfate ion show large similarities in hydration structure and mechanism with the sulfate ion but with weaker hydration of the terminal sulfur atom in thiosulfate.

Dynamics of actinyl ions in water: a molecular dynamics simulation study

The dynamics of actinyl ions (AnO₂ⁿ⁺) in aqueous solutions is important not only for the design of advanced separation processes but also for understanding the fate of actinides in the environment.