Environment of Ca2+ ions in aqueous solvent

Hydration Structure of 102No2+: A Density Functional Theory-Molecular Dynamics Study

The hydration structure of No2+, the divalent cation of nobelium in water, was investigated by ab initio molecular dynamics (MD) simulations. First, a series of benchmark calculations were performed to validate the density functional theory (DFT) calculation methods for a molecule containing a No atom. The DFT-MD simulation of the hydration structure of No2+ was conducted after the MD method was validated by simulating the hydration structures of Ca2+ and Sr2+, whose behavior was previously reported to be similar to that of No2+. The model cluster containing M2+ (M = Ca, Sr, or No) and 32 water molecules was used for DFT-MD simulation. The results showed that the hydration distance of No2+ was intermediate between those of Ca2+ and Sr2+. This trend in the hydration distance is in good agreement with the elution position trend obtained in a previous radiochemical experiment. The calculated No-O bond lengths in the optimized structure of [No(H2O)8]2+ was 2.59 Å, while the average No-O bond length of [No(H2O)8]2+ in water by DFT-MD was 2.55 Å. This difference implies the importance of dynamic solvent effects, considering the second (and further) coordination sphere in the theoretical calculation of solution chemistry for superheavy elements.

On Local Structure Equilibration of Ca2+ in Solution by Ab Initio Molecular Dynamics

Analyzing the stable isotopic ratio of Ca offers valuable insights into a wide range of applications from climate reconstruction to bone cancer diagnosis and agricultural nutrient improvement. While the first hydration shell of Ca in solution is expected to play a major role in its fractionation properties, the coordination of Ca in water remains a subject of debate. In this work, Ca2+ in water has been modeled by means of ab initio molecular dynamics simulations using various exchange and correlation functionals and at different temperatures. Results show a significant effect of the functional on the average Ca2+ coordination, depending on its tendency to over- or understructure liquid water. The BLYP functional with Grimme-D2 correction was judged as the most accurate among those tested based on its accuracy to reproduce water structural and diffusion properties. Using this functional, the effect of temperature has been systematically investigated, focusing on means to limit the uncertainty in our assessments of the average coordination of Ca2+ ions by (1) estimating the number of water exchanges in the simulations and (2) implementing a statistical approach based on Markov chains. The findings indicate, especially, that our simulations at 300, 350, and 400 K do not yield converged results due to potential equilibration problems. These observations impose substantial constraints on the trustworthiness of numerous estimates in the existing literature that depend on trajectories with insufficient exchanges. We estimate Ca2+ coordination values of 6.8 ± 0.1, 6.8 ± 0.1, 6.7 ± 0.2, and 6.7 ± 0.2 at 600, 550, 500, and 450 K respectively. At lower temperatures (300, 350, and 400 K), while obtaining definitive values for Ca2+ coordination remains challenging, our research does indicate a potential temperature-related influence on coordination with an average Ca2+ coordination at 300 K as low as 6.2.

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Hydration and, in particular, the coordination number of a metal ion is of paramount importance as it defines many of its (bio)physicochemical properties. It is not only essential for understanding its behavior in aqueous solutions but also determines the metal ion reference state and its binding energy to (bio)molecules. In this paper, for divalent metal cations Ca2+, Cd2+, Cu2+, Fe2+, Hg2+, Mg2+, Ni2+, Pb2+, and Zn2+, we compare two approaches for predicting hydration numbers: (1) a mixed explicit/continuum DFT-D3//COSMO-RS solvation model and (2) density functional theory based ab initio molecular dynamics. The former approach is employed to calculate the Gibbs free energy change for the sequential hydration reactions, starting from [M(H2O)2]2+ aqua complexes to [M(H2O)9]2+, allowing explicit water molecules to bind in the first or second coordination sphere and determining the most stable [M(H2O)n]2+ structure. In the latter approach, the hydration number is obtained by integrating the ion-water radial distribution function. With a couple of exceptions, the metal ion hydration numbers predicted by the two approaches are in mutual agreement, as well as in agreement with the experimental data.

Characterizing the Mechanisms of Ca and Mg Carbonate Ion-Pair Formation with Multi-Level Molecular Dynamics/Quantum Mechanics Simulations

The carbonate minerals of Ca and Mg are abundant throughout the lithosphere and have recently garnered significant research interest as possible long-term carbon sinks in the sequestration of atmospheric carbon dioxide. Nonetheless, an understanding of the atomic-level processes comprising their mineralization remains limited. Here, we characterize and contrast the mechanisms of contact ion-pair formation in aqueous Ca and Mg carbonate systems, which represents the most fundamental step leading to the formation of their mineral solids. Utilizing multilevel embedded correlated wavefunction-based ab initio molecular dynamics/quantum mechanics simulations, we characterize not only the dynamics of these processes but also factors arising from the electronic structure of the involved species, revealing further details of the fundamentally different mechanisms for the interconversion between the contact ion-pairs and solvent-shared ion-pairs of Ca versus Mg carbonate.

Physiology of intracellular calcium buffering

Calcium signaling underlies much of physiology. Almost all the Ca2+ in the cytoplasm is bound to buffers, with typically only ∼1% being freely ionized at resting levels in most cells. Physiological Ca2+ buffers include small molecules and proteins, and experimentally Ca2+ indicators will also buffer calcium. The chemistry of interactions between Ca2+ and buffers determines the extent and speed of Ca2+ binding. The physiological effects of Ca2+ buffers are determined by the kinetics with which they bind Ca2+ and their mobility within the cell. The degree of buffering depends on factors such as the affinity for Ca2+, the Ca2+ concentration, and whether Ca2+ ions bind cooperatively. Buffering affects both the amplitude and time course of cytoplasmic Ca2+ signals as well as changes of Ca2+ concentration in organelles. It can also facilitate Ca2+ diffusion inside the cell. Ca2+ buffering affects synaptic transmission, muscle contraction, Ca2+ transport across epithelia, and the killing of bacteria. Saturation of buffers leads to synaptic facilitation and tetanic contraction in skeletal muscle and may play a role in inotropy in the heart. This review focuses on the link between buffer chemistry and function and how Ca2+ buffering affects normal physiology and the consequences of changes in disease. As well as summarizing what is known, we point out the many areas where further work is required.

Probing pH-Dependent Dehydration Dynamics of Mg and Ca Cations in Aqueous Solutions with Multi-Level Quantum Mechanics/Molecular Dynamics Simulations

The dehydration of aqueous calcium and magnesium cations is the most fundamental process controlling their reactivity in chemical and biological phenomena, such as the formation of ionic solids or passing through ion channels. It holds particular relevance in light of recent advancements in the development of carbon capture techniques that rely on mineralization for long-term carbon storage. Specifically, dehydration of Ca2+ and Mg2+ is a key step in proposed carbon capture processes aiming to exploit the relatively high concentration of dissolved carbon dioxide in seawater via the formation of carbonate minerals from solvated Ca2+ and Mg2+ cations for sequestration and storage. Nevertheless, atomic-scale understanding of the dehydration of aqueous Ca2+ and Mg2+ cations remains limited. Here, we utilize rare event sampling via density functional theory molecular dynamics and embedded wavefunction theory calculations to elucidate the dehydration dynamics of aqueous Ca2+ and Mg2+. Emphasis is placed on the investigation of the effect pH has on the stability of the different coordination environments. Our results reveal significant differences in the dehydration dynamics of the two cations and provide insight into how they may be modulated by pH changes.

Solvation structures of calcium and magnesium ions in water with the presence of hydroxide: a study by deep potential molecular dynamics

The solvation structures of calcium (Ca²⁺) and magnesium (Mg²⁺) ions with the presence of hydroxide (OH^-) ion in water are essential for understanding their roles in biological and chemical processes but have not been fully explored. Ab initio molecular dynamics (AIMD) is an important tool to address this issue, but two challenges exist. First, an accurate description of OH^- from AIMD needs an appropriate exchange-correlation functional. Second, a long trajectory is needed to reach an equilibrium state for the Ca²⁺-OH^- and Mg²⁺-OH^- ion pairs in aqueous solutions. Herein, we adopt a deep potential molecular dynamics (DPMD) method to simulate 1 ns trajectories for the Ca²⁺-OH^- and Mg²⁺-OH^- ion pairs in water; the DPMD method provides efficient machine-learning-based models that have the accuracy of the SCAN exchange-correlation functional within the framework of density functional theory. The solvation structures of the cations and the OH^- in terms of three different species have been systematically investigated. On the one hand, we find that OH^- have more significant effects on the solvation structure of Ca²⁺ than that of Mg²⁺. We observe that the OH^- substantially affects the orientation angles of water molecules surrounding the cation. Through the time correlation functions, we conclude that the water molecules in the first solvation shell of Ca²⁺ change their preferred orientation faster than those of Mg²⁺. On the other hand, with the presence of the cation in the first solvation shell of OH^-, we find that the hydrogen bonds of OH^- are severely altered, and the adjacent water molecules of OH^- are squeezed. The two cations have substantially different effects on the solvation structure of OH^-. Our work provides new insight into the solvation structures of Ca²⁺ and Mg²⁺ in water with the presence of OH^-.

Water-in-Salt Electrolyte-Based Extended Voltage Range, Safe, and Long-Cycle-Life Aqueous Calcium-Ion Cells

The narrow electrochemical stability window (1.23 V) of an aqueous electrolyte hinders the practical realization of calcium-ion chemistries of high-energy-density and long-cycle-life batteries. Furthermore, developing an aqueous electrolyte that is low cost, is environmentally friendly, and has a wide voltage window is essential to designing safe, high-energy-density, and sustainable calcium-ion batteries. A calcium-based water-in-salt (WISE) aqueous electrolyte surpasses the narrow stability window by offering a 2.12 V wide window by suppressing the hydrogen evolution at the anode and minimizing the overall water activity at the cathode. A comprehensive theoretical study predicts the preferential reduction of salt aggregates over water to form a passivation layer at the electrode-electrolyte interface and enhance the electrolyte stability window. Additionally, Raman spectroscopy reveals that the calcium ion coordination number, which is the number of nitrate ions surrounding the calcium ions in the aqueous electrolyte, gradually increases with an increase in the electrolyte concentration, leading to a gradual decrease in the hydration number of the calcium ions. A full cell in WISE was demonstrated to exhibit an excellent rate capability and cycling stability with negligible capacity loss (0.01 per cycle), maintaining 80% capacity retention over 1800 cycles with ∼99.99% Coulombic efficiency. The full cell provides an energy density of 232 Wh kg^-1 at a power density of 69 W kg^-1 and a current rate of 0.15 A g^-1. Even at a higher current rate of 5 A g^-1, the battery delivers an energy density of 182 Wh kg^-1 (based on the active mass of the anode). This is one of the best performances to date of all previously reported full-cell aqueous calcium-ion batteries. A fundamental understanding of the storage mechanism and a electrode degradation study was achieved. This work suggests and expands new avenues for the practical realization of low-cost, safe, eco-friendly, and high-performance aqueous calcium-ion batteries for future large storage applications.

Ground-State Structures of Hydrated Calcium Ion Clusters From Comprehensive Genetic Algorithm Search

We searched the lowest-energy structures of hydrated calcium ion clusters Ca2+(H2O)n (n = 10-18) in the whole potential energy surface by the comprehensive genetic algorithm (CGA). The lowest-energy structures of Ca2+(H2O)10-12 clusters show that Ca2+ is always surrounded by six H2O molecules in the first shell. The number of first-shell water molecules changes from six to eight at n = 12. In the range of n = 12-18, the number of first-shell water molecules fluctuates between seven and eight, meaning that the cluster could pack the water molecules in the outer shell even though the inner shell is not full. Meanwhile, the number of water molecules in the second shell and the total hydrogen bonds increase with an increase in the cluster size. The distance between Ca2+ and the adjacent water molecules increases, while the average adjacent O-O distance decreases as the cluster size increases, indicating that the interaction between Ca2+ and the adjacent water molecules becomes weaker and the interaction between water molecules becomes stronger. The interaction energy and natural bond orbital results show that the interaction between Ca2+ and the water molecules is mainly derived from the interaction between Ca2+ and the adjacent water molecules. The charge transfer from the lone pair electron orbital of adjacent oxygen atoms to the empty orbital of Ca2+ plays a leading role in the interaction between Ca2+ and water molecules.

Comparison of the Microsolvation of CaX2 (X = F, Cl, Br, I) in Water: Size-Selected Anion Photoelectron Spectroscopy and Theoretical Calculations

To understand the microsolvation of alkaline-earth dihalides in water and provide information about the dependence of solvation processes on different halides, we investigated CaBr₂(H₂O)_n^-, CaI₂(H₂O)_n^-, and CaF₂(H₂O)_n^- (n = 0-6) clusters using size-selected anion photoelectron spectroscopy and conducted theoretical calculations on these clusters and their neutrals. The results are compared with those of CaCl₂(H₂O)_n^-/0 clusters reported previously. It is found that the vertical detachment energies (VDEs) of CaCl₂(H₂O)_n^-, CaBr₂(H₂O)_n^-, and CaI₂(H₂O)_n^- show a similar trend with increasing cluster size, while the VDEs of CaF₂(H₂O)_n^- show a different trend. The VDEs of CaF₂(H₂O)_n^- are much lower than those of CaCl₂(H₂O)_n^-, CaBr₂(H₂O)_n^-, and CaI₂(H₂O)_n^-. A detailed probing of the structures shows that a significant increase of the Ca-X distance (separation of Ca²⁺-X^- ion pair) in CaCl₂(H₂O)_n^-/0, CaBr₂(H₂O)_n^-/0, and CaI₂(H₂O)_n^-/0 clusters occurred at about n = 5. However, for CaF₂(H₂O)_n^-/0, no abrupt change of the Ca-F distance with the increasing cluster size has been observed. In CaCl₂(H₂O)₆^-/0, CaBr₂(H₂O)₆^-/0, and CaI₂(H₂O)₆^-/0, the Ca atom coordinates directly with 5 H₂O molecules. However, in CaF₂(H₂O)_n^-/0, the Ca atom coordinates directly with only 2 or 3 H₂O molecules. The similarity or differences in the structures and coordination numbers are consistent with the fact that CaCl₂, CaBr₂, and CaI₂ have similar solubility, while CaF₂ has much lower solubility.

The determination of the thermodynamic constants of MgUO2(CO3)32- complex in NaClO4 and NaCl media by time-resolved luminescence spectroscopy, and applications in different geochemical contexts

The formation constants and specific ion interaction coefficients of MgUO₂(CO₃)₃^2- complex were determined in 0.1 to 1.0 mol kg_w^-1 NaCl and 0.10 to 2.21 mol kg_w^-1 NaClO₄ media in the framework of the specific ion interaction theory (SIT), by time-resolved laser-induced luminescence spectroscopy. The upper limits of ionic strength were chosen in order to limit luminescence quenching effects generated by high concentrations of Cl^- and ClO₄^- already observed during our earlier studies on Ca_nUO₂(CO₃)₃^(4-2n)- complexes (Shang & Reiller, Dalton Trans., 49, 466; Shang et al., Dalton Trans., 49, 15443). The cumulative formation constant determined is , and the specific ion interaction coefficients are ε(MgUO₂(CO₃)₃^2-, Na⁺) = 0.19 ± 0.11 kg_w mol^-1 in NaClO₄ and ε(MgUO₂(CO₃)₃^2-, Na⁺) = 0.09 ± 0.16 kg_w mol^-1 in NaCl. Two gratings of 300 and 1800 lines per mm were used to acquire MgUO₂(CO₃)₃^2- luminescence spectra, where the high-resolution 1800 lines per mm grating detected slight spectral shifts for the principal luminescent bands relative to Ca_nUO₂(CO₃)₃^(4-2n)-. The applications of the consistent set of thermodynamic constants and ε values for M_nUO₂(CO₃)₃^(4-2n)- (M = Mg and Ca) were examined in different geochemical contexts, where Mg over Ca concentration ratio varies to help defining the relative importance of these species.

Modulation of the calcium oxalate dihydrate to calcium oxalate monohydrate phase transition with citrate and zinc ions

Higher concentrations of Ca2+ and Ox2− can form COD which then transforms to COM. Citrate forms a protective layer to inhibit COD transition; whereas Zn2+ substitutes Ca2+ sites to generate a stable COD structure that retards COM formation.

Deciphering the Solvent Effect for the Solvation Structure of Ca2+ in Polar Molecular Liquids

Although the crystal structures for many inorganic compounds are readily available, researchers are still working hard to understand the relations between the structures and chemical properties of solutions because most of the chemical reactions take place in solutions. A huge amount of effort has been put toward modeling the ion solvation structure from the perspectives of both experiments and theories. In this study, the solvation structures of Ca²⁺ ions in aqueous and alcoholic solutions at different concentrations were carefully evaluated by Ca K-edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) analyses. Density functional theory (DFT) calculations were also performed to correlate the results with the experimental data and then further extended to other similar systems. It was found that the number of coordinating solvent molecules decreases with increasing Ca²⁺ concentration and increasing solvent molecule sizes. From the EXAFS data, it was observed that the first solvation shell of Ca²⁺ splits into two Ca-O distances in a methanol solution and the counter ion Cl^- might also be within the first shell at high concentrations. For the first time, the effects of solvents with different polarities and sizes on the ion solvation environment were systematically evaluated.

Voltage issue of aqueous rechargeable metal-ion batteries

Over the past two decades, a series of aqueous rechargeable metal-ion batteries (ARMBs) have been developed, aiming at improving safety, environmental friendliness and cost-efficiency in fields of consumer electronics, electric vehicles and grid-scale energy storage. However, the notable gap between ARMBs and their organic counterparts in energy density directly hinders their practical applications, making it difficult to replace current widely-used organic lithium-ion batteries. Basically, this huge gap in energy density originates from cell voltage, as the narrow electrochemical stability window of aqueous electrolytes substantially confines the choice of electrode materials. This review highlights various ARMBs with focuses on their voltage characteristics and strategies that can effectively raise battery voltage. It begins with the discussion on the fundamental factor that limits the voltage of ARMBs, i.e., electrochemical stability window of aqueous electrolytes, which decides the maximum-allowed potential difference between cathode and anode. The following section introduces various ARMB systems and compares their voltage characteristics in midpoint voltage and plateau voltage, in relation to respective electrode materials. Subsequently, various strategies paving the way to high-voltage ARMBs are summarized, with corresponding advancements highlighted. The final section presents potential directions for further improvements and future perspectives of this thriving field.

Aqueous solvation of Mg(ii) and Ca(ii): A Born-Oppenheimer molecular dynamics study of microhydrated gas phase clusters

The hydration features of [Mg(H₂O)_n]²⁺ and [Ca(H₂O)_n]²⁺ clusters with n = 3-6, 8, 18, and 27 were studied by means of Born-Oppenheimer molecular dynamics simulations at the B3LYP/6-31+G** level of theory. For both ions, it is energetically more favorable to have all water molecules in the first hydration shell when n ≤ 6, but stable lower coordination average structures with one water molecule not directly interacting with the ion were found for Mg²⁺ at room temperature, showing signatures of proton transfer events for the smaller cation but not for the larger one. A more rigid octahedral-type structure for Mg²⁺ than for Ca²⁺ was observed in all simulations, with no exchange of water molecules to the second hydration shell. Significant thermal effects on the average structure of clusters were found: while static optimizations lead to compact, spherically symmetric hydration geometries, the effects introduced by finite-temperature dynamics yield more prolate configurations. The calculated vibrational spectra are in agreement with infrared spectroscopy results. Previous studies proposed an increase in the coordination number (CN) from six to eight water molecules for [Ca(H₂O)_n]²⁺ clusters when n ≥ 12; however, in agreement with recent measurements of binding energies, no transition to a larger CN was found when n > 8. Moreover, the excellent agreement found between the calculated extended X-ray absorption fine structure spectroscopy spectra for the larger cluster and the experimental data of the aqueous solution supports a CN of six for Ca²⁺.

Calcium isotope fractionation between aqueous compounds relevant to low-temperature geochemistry, biology and medicine

Stable Ca isotopes are fractionated between bones, urine and blood of animals and between soils, roots and leaves of plants by >1000 ppm for the 44Ca/40Ca ratio. These isotopic variations have important implications to understand Ca transport and fluxes in living organisms; however, the mechanisms of isotopic fractionation are unclear. Here we present ab initio calculations for the isotopic fractionation between various aqueous species of Ca and show that this fractionation can be up to 3000 ppm. We show that the Ca isotopic fractionation between soil solutions and plant roots can be explained by the difference of isotopic fractionation between the different first shell hydration degree of Ca2+ and that the isotopic fractionation between roots and leaves is controlled by the precipitation of Ca-oxalates. The isotopic fractionation between blood and urine is due to the complexation of heavy Ca with citrate and oxalates in urine. Calculations are presented for additional Ca species that may be useful to interpret future Ca isotopic measurements.

Cation solvation with quantum chemical effects modeled by a size-consistent multi-partitioning quantum mechanics/molecular mechanics method

In the condensed phase, quantum chemical properties such as many-body effects and intermolecular charge fluctuations are critical determinants of the solvation structure and dynamics.

Molecular dynamics simulation for the impact of an electrostatic field and impurity Mg2+ions on hard water

A detailed analysis of the structural parameters and dynamic parameters of hard water solutions under an external electrostatic field was performed by molecular dynamics (MD) simulations with the presence of impurity Mg²⁺ions.

Prediction of the pKa's of aqueous metal ion +2 complexes

Aqueous metal ions play an important role in many areas of chemistry. The acidities of [Be(H2O)4](2+), [M(H2O)6](2+), M = Mg(2+), Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), Zn(2+), Cd(2+), and Hg(2+), and [M(H2O)n](2+), M = Ca(2+) and Sr(2+), n = 7 and 8, complexes have been predicted using density functional theory, second-order Møller-Plesset perturbation theory (MP2), and coupled cluster CCSD(T) theory in the gas phase. pKa's in aqueous solution were predicted by using self-consistent reaction field (SCRF) calculations with different solvation models. The most common binding motif of the majority of the metal +2 complexes is coordination number (CN) 6, with each hexaaquo cluster having reasonably high symmetry for the best arrangement of the water molecules in the first solvation shell. Be(2+) is tetracoordinated, but a second solvation shell of 8 waters is needed to predict the pKa. The Ca(2+) and Sr(2+) aquo clusters have a coordination number of 7 or 8 as found in terms of the energy of the reaction M(H2O)7(2+) + H2O → M(H2O)8(2+) and the pKa values. The calculated geometries are in reasonable agreement with experiment. The SCRF calculations with the conductor-like screening model (COSMO), and the conductor polarized continuum model (CPCM) using COSMO-RS radii, consistently agree best with experiment at the MP2/aug-cc-pVDZ and CCSD(T)/aug-cc-pVDZ levels of theory. The CCSD(T) level provides the most accurate pKa's, and the MP2 level also provides reliable predictions. Our predictions were used to elucidate the properties of metal +2 ion complexes. The pKa predictions provide confirmation of the size of the first solvation shell sizes. The calculations show that it is still difficult to predict pKa's using this cluster/implicit solvent approach to better than 1 pKa unit.

Multisite ion model in concentrated solutions of divalent cations (MgCl2 and CaCl2): osmotic pressure calculations

Accurate force field parameters for ions are essential for meaningful simulation studies of proteins and nucleic acids. Currently accepted models of ions, especially for divalent ions, do not necessarily reproduce the right physiological behavior of Ca(2+) and Mg(2+) ions. Saxena and Sept (J. Chem. Theor. Comput. 2013, 9, 3538-3542) described a model, called the multisite-ion model, where instead of treating the ions as an isolated sphere, the charge was split into multiple sites with partial charge. This model provided accurate inner shell coordination of the ion with biomolecules and predicted better free energies for proteins and nucleic acids. Here, we expand and refine the multisite model to describe the behavior of divalent ions in concentrated MgCl2 and CaCl2 electrolyte solutions, eliminating the unusual ion-ion pairing and clustering of ions which occurred in the original model. We calibrate and improve the parameters of the multisite model by matching the osmotic pressure of concentrated solutions of MgCl2 to the experimental values and then use these parameters to test the behavior of CaCl2 solutions. We find that the concentrated solutions of both divalent ions exhibit the experimentally observed behavior with correct osmotic pressure, the presence of solvent separated ion pairs instead of direct ion pairs, and no aggregation of ions. The improved multisite model for (Mg(2+) and Ca(2+)) can be used in classical simulations of biomolecules at physiologically relevant salt concentrations.

Direct measurement of the Mn(II) hydration state in metal complexes and metalloproteins through 17O NMR line widths

Here we describe a simple method to estimate the inner-sphere hydration state of the Mn(II) ion in coordination complexes and metalloproteins. The line width of bulk H2(17)O is measured in the presence and absence of Mn(II) as a function of temperature, and transverse (17)O relaxivities are calculated. It is demonstrated that the maximum (17)O relaxivity is directly proportional to the number of inner-sphere water ligands (q). Using a combination of literature data and experimental data for 12 Mn(II) complexes, we show that this method provides accurate estimates of q with an uncertainty of ±0.2 water molecules. The method can be implemented on commercial NMR spectrometers working at fields of 7 T and higher. The hydration number can be obtained for micromolar Mn(II) concentrations. We show that the technique can be extended to metalloproteins or complex:protein interactions. For example, Mn(II) binds to the multimetal binding site A on human serum albumin with two inner-sphere water ligands that undergo rapid exchange (1.06 × 10(8) s(-1) at 37 °C). The possibility of extending this technique to other metal ions such as Gd(III) is discussed.

The aqueous Ca2+ system, in comparison with Zn2+, Fe3+, and Al3+: an ab initio molecular dynamics study

Herein, we report on the structure and dynamics of the aqueous Ca(2+) system studied by using ab initio molecular dynamics (AIMD) simulations. Our detailed study revealed the formation of well-formed hydration shells with characteristics that were significantly different to those of bulk water. To facilitate a robust comparison with state-of-the-art X-ray absorption fine structure (XAFS) data, we employ a 1st principles MD-XAFS procedure and directly compare simulated and experimental XAFS spectra. A comparison of the data for the aqueous Ca(2+) system with those of the recently reported Zn(2+), Fe(3+), and Al(3+) species showed that many of their structural characteristics correlated well with charge density on the cation. Some very important exceptions were found, which indicated a strong sensitivity of the solvent structure towards the cation's valence electronic structure. Average dipole moments for the 2nd shell of all cations were suppressed relative to bulk water.

Probing the Structure of Salt Water under Confinement with First-Principles Molecular Dynamics and Theoretical X-ray Absorption Spectroscopy

We investigated the structure of liquid water around cations (Na(+)) and anions (Cl(-)) confined inside of a (19,0) carbon nanotube with first-principles molecular dynamics and theoretical X-ray absorption spectroscopy (XAS). We found that the ions preferentially reside near the interface between the nanotube and the liquid. Upon confinement, the XAS signal of water molecules surrounding Na(+) exhibits enhanced pre-edge and reduced post-edge features with respect to that of pure water, at variance with the solvation shell of Na(+) in bulk water. Conversely, the first solvation shell of confined Cl(-) has a main-edge intensity comparable to that of bulk solvated Cl(-), likely as a result of a high number of acceptor hydrogen bonds in the first solvation shell. Confined nonsolvating water molecules exhibit bulk-like or water-monomer-like properties, depending on whether they belong to core or interfacial layers, respectively.

Aqueous solutions of divalent chlorides: ions hydration shell and water structure

By combining neutron diffraction and Monte Carlo simulations, we have determined the microscopic structure of the hydration ions shell in aqueous solutions of MgCl(2) and CaCl(2), along with the radial distribution functions of the solvent. In particular the hydration shell of the cations, show cation specific symmetry, due to the strong and directional interaction of ions and water oxygens. The ions and their hydration shells likely form molecular moieties and bring clear signatures in the water-water radial distribution functions. Apart from these signatures, the influence of divalent salts on the microscopic structure of water is similar to that of previously investigated monovalent solutes, and it is visible as a shift of the second peak of the oxygen-oxygen radial distribution function, caused by distortion of the hydrogen bond network of water.

Hydration of calcium oxide surface predicted by reactive force field molecular dynamics

In this work, we present the parametrization of Ca-O/H interactions within the reactive force field ReaxFF, and its application to study the hydration of calcium oxide surface. The force field has been fitted using density functional theory calculations on gas phase calcium-water clusters, calcium oxide bulk and surface properties, calcium hydroxide, bcc and fcc Ca, and proton transfer reactions in the presence of calcium. Then, the reactive force field has been used to study the hydration of the calcium oxide {001} surface with different water contents. Calcium oxide is used as a catalyzer in many applications such as CO(2) sequestration and biodiesel production, and the degree of surface hydroxylation is a key factor in its catalytic performance. The results show that the water dissociates very fast on CaO {001} bare surfaces without any defect or vacancy. The surface structure is maintained up to a certain amount of water, after which the surface undergoes a structural rearrangement, becoming a disordered calcium hydroxyl layer. This transformation is the most probable reason for the CaO catalytic activity decrease.

Structures, stability, vibration entropy and IR spectra of hydrated calcium ion clusters [Ca(H(2)O)(n)](2+) (n = 1-20, 27): a systematic investigation by density functional theory

The low-lying candidates of hydrated calcium ion clusters, [Ca(H(2)O)(n)](2+) with n = 1-20 and 27, have been extensively sought by using density functional theory (DFT) at BLYP/6-311+G(d,p) level. The results showed that the first hydration shell around the calcium ion was fully occupied by six water molecules, whereas the second hydration shell might be fully occupied with different numbers of water molecules. This just corresponds to different growth patterns of the hydrated calcium ion clusters. Furthermore, we revealed that the vibration entropy contributed to the free energy of an isomer significantly. As a result, the stability of some low-lying candidates at zero-temperature was not maintained at finite temperatures. Therefore, we suggested that, at finite temperatures, the realistic products of [Ca(H(2)O)(n)](2+) should be a mixture of the best candidate and some of metastable isomers for a given cluster size. For a cluster having second and/or third shell of water molecules, we found structural transitions between a low-lying structure and the lowest-energy structure undergoing much lower energy barriers. In addition, the IR spectra of the best candidates were predicted, in which the evolution of hydrogen-bond configurations with the cluster size was revealed.