Anomalous x-ray diffraction studies of Sr2+ hydration in aqueous solution

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.

Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2

In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.

Identification of strontium substitution mechanism in hematite via calcium solution

Nonradioactive strontium (Sr) are produced as a result of radioactive decay of heavier elements such as uranium and thorium. Nonradioactive Sr shares physicochemical similarities with Ca and can replace it during bone formation, which may cause bone cancer in humans. Hence, concerning the potential hazards associated with strontium, it is imperative to eliminate it. The present study aimed to investigate the removal mechanisms of hematite-adsorbed strontium by calcium solution. Strontium was adsorbed to hematite at pH 8 and 10 and washed with calcium solution. X-ray absorption near-edge structure (XANES), extended X-ray absorption fine structure (EXAFS), scanning electron microscopy, X-ray diffraction (XRD), and X-ray photoelectron spectroscopy (after Ca washing) were performed on the samples before and after washing. Analyses and fitting by XANES and EXAFS confirmed the formation of an inner-sphere complex of strontium at pH 10. The XRD spectra showed that SrCO3 and SrFe2O4 formed at pH 8 and 10, respectively. After washing with the calcium solution, strontium was directly substituted to form CaCO3 and CaFe2O4. The X-ray photoelectron spectroscopy results provided a systematic analysis of the proportions of hematite and strontium, confirming the substitution of strontium with calcium. This substitution could be attributed to the physicochemical similarities between calcium and strontium. This study confirms the substitution of Sr with Ca, highlighting the physicochemical similarity of the Sr and Ca that facilitates substitution reactions.

Hydration and Charge-Transfer Effects of Alkaline Earth Metal Ions Binding to a Carboxylate Anion, Phosphate Anion, and Guanine Nucleobase

To investigate the ability of alkaline earth metal ions to tune ion-mediated DNA adsorption, hydrated Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺ ions bound to a carboxylate anion, phosphate anion, and guanine nucleobase were modeled using density functional theory (DFT) and a combined explicit and continuum solvent model. The large first solvation shell of Ba²⁺ requires a larger solute cavity defined by a solvent-accessible surface, which is used to model all hydrated ions. Alkaline earth metal ions bind indirectly or directly to each binding site. DFT binding energies decrease with increasing ion size, which is likely due to ion size and hydration structure, rather than quantum effects such as charge transfer. However, charge transfer explains weaker ion binding to guanine compared to phosphate or carboxylate. Overall, carboxylate and phosphate anions are expected to compete equally for hydrated Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺ ions and larger alkaline earth metal ions may induce weaker ion-mediated adsorption. The ion size and hydration structure of alkaline earth metal ions may effectively tune ion-mediated adsorption processes, such as DNA adsorption to functionalized surfaces.

Hydration of Heavy Alkaline-Earth Cations Studied by Molecular Dynamics Simulations and X-ray Absorption Spectroscopy

The physicochemical properties of the three heaviest alkaline-earth cations, Sr²⁺, Ba²⁺, and Ra²⁺ in water have been studied by means of classical molecular dynamics (MD) simulations. A specific set of cation-water intermolecular potentials based on ab initio potential energy surfaces has been built on the basis of the hydrated ion concept. The polarizable and flexible model of water MCDHO2 was adopted. The theoretical-experimental comparison of structural, dynamical, energetic, and spectroscopical properties of Sr²⁺ and Ba²⁺ aqueous solutions is satisfactory, which supports the methodology developed. This good behavior allows a reasonable reliability for the predicted Ra²⁺ physicochemical data not experimentally determined yet. Simulated extended X-ray absorption fine-structure (EXAFS) and X-ray absorption near-edge spectroscopy spectra have been computed from the snapshots of the MD simulations and compared with the experimental information available for Sr²⁺ and Ba²⁺. For the Ra²⁺ case, the Ra L₃-edge EXAFS spectrum is proposed. Structural and dynamical properties of the aqua ions for the three cations have been obtained and analyzed. Along the [M(H₂O)_n]^m+ series, the M-O distance for the first-hydration shell is 2.57, 2.81, and 2.93 Å for Sr²⁺, Ba²⁺, and Ra²⁺, respectively. The hydration number also increases when one is going down along the group: 8.1, 9.4, and 9.8 for Sr²⁺, Ba²⁺, and Ra²⁺, respectively. Whereas [Sr(H₂O)₈]²⁺ is a typical aqua ion with a well-defined structure, the Ba²⁺ and Ra²⁺ hydration provides a picture exhibiting an average between the ennea- and the deca-hydration. These results show a similar chemical behavior of Ba²⁺ and Ra²⁺ aqueous solutions and support experimental studies on the removal of Ra-226 of aquifers by different techniques, where Ra²⁺ is replaced by Ba²⁺. A comparison of the heavy alkaline ions, Rb⁺ and Cs⁺, with the heavy alkaline-earth ions is made.

Hydration Mimicry by Membrane Ion Channels

Ions transiting biomembranes might pass readily from water through ion-specific membrane proteins if these protein channels provide environments similar to the aqueous solution hydration environment. Indeed, bulk aqueous solution is an important reference condition for the ion permeation process. Assessment of this hydration mimicry concept depends on understanding the hydration structure and free energies of metal ions in water in order to provide a comparison for the membrane channel environment. To refine these considerations, we review local hydration structures of ions in bulk water and the molecular quasi-chemical theory that provides hydration free energies. In doing so, we note some current views of ion binding to membrane channels and suggest new physical chemical calculations and experiments that might further clarify the hydration mimicry concept.

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

Classical molecular dynamics simulations reveal size-dependent trends of alkaline earth metal ions binding to DNA are due to ion size and hydration behavior.

Ab initio molecular dynamics studies of hydroxide coordination of alkaline earth metals and uranyl

Ab initio molecular dynamics (AIMD) simulations of the Mg2+, Ca2+, Sr2+ and UO22+ ions in either a pure aqueous environment or an environment containing two hydroxide ions have been carried out at the density functional level of theory, employing the generalised gradient approximation via the PBE exchange-correlation functional. Calculated mean M-O bond lengths in the first solvation shell of the aquo systems compared very well to existing experimental and computational literature, with bond lengths well within values measured previously and coordination numbers in line with previously calculated values. When applied to systems containing additional hydroxide ions, the methodology revealed increased bond lengths in all systems. Proton transfer events (PTEs) were recorded and were found to be most prevalent in the strontium hydroxide systems, likely due to the low charge density of the ion and the consequent lack of hydroxide coordination. For all alkaline earths, intrashell PTEs which occurred outside of the first solvation shell were most prevalent. Only three PTEs were identified in the entire simulation data of the uranium dihydroxide system, indicating the clear impact of the increased charge density of the hexavalent uranium ion on the strength of metal-oxygen bonds in aqueous solution. Broadly, systems containing more charge dense ions were found to exhibit fewer PTEs than those containing ions of lower charge density.

Ion Hydration and Ion Pairing as Probed by THz Spectroscopy

Ion hydration is of pivotal importance for many fundamental processes. Various spectroscopic methods are used to study the retardation of the hydration bond dynamics in the vicinity of anions and cations. Here we introduce THz-FTIR spectroscopy as a powerful method to answer the open questions. We show through dissection of THz spectra that we can pinpoint characteristic absorption features that can be attributed to the rattling modes of strongly hydrating ions within their hydration cages as well as vibrationally induced charge fluctuations in the case of weakly hydrating ions. Further analysis yields information on anion-cation cooperativity, the size of the dynamic hydration shell, as well as the lifetimes of these collective ion-hydration water modes and their connecting thermal bath states. Our study provides evidence for a non-additive behavior, thus questioning the simplified Hofmeister model. THz spectroscopy enables ion pairing to be observed and quantified at a high salt concentration.

Strontium and barium in aqueous solution and a potassium channel binding site

Ion hydration structure and free energy establish criteria for understanding selective ion binding in potassium (K⁺) ion channels and may be significant to understanding blocking mechanisms as well. Recently, we investigated the hydration properties of Ba²⁺, the most potent blocker of K⁺ channels among the simple metal ions. Here, we use a similar method of combining ab initio molecular dynamics simulations, statistical mechanical theory, and electronic structure calculations to probe the fundamental hydration properties of Sr²⁺, which does not block bacterial K⁺ channels. The radial distribution of water around Sr²⁺ suggests a stable 8-fold geometry in the local hydration environment, similar to Ba²⁺. While the predicted hydration free energy of -331.8 kcal/mol is comparable with the experimental result of -334 kcal/mol, the value is significantly more favorable than the -305 kcal/mol hydration free energy of Ba²⁺. When placed in the innermost K⁺ channel blocking site, the solvation free energies and lowest energy structures of both Sr²⁺ and Ba²⁺ are nearly unchanged compared with their respective hydration properties. This result suggests that the block is not attributable to ion trapping due to +2 charge, and differences in blocking behavior arise due to free energies associated with the exchange of water ligands for channel ligands instead of free energies of transfer from water to the binding site.

A Kirkwood-Buff derived force field for alkaline earth halide salts

The activity and function of many macromolecules in cellular environments are coupled with the binding of divalent ions such as calcium or magnesium. In principle, computer simulations can be used to understand the molecular level aspects of how many important macromolecules interact with ions. However, most of the force fields currently available often fail to accurately reproduce the properties of divalent ions in aqueous environments. Here we develop classical non-polarizable force fields for the aqueous alkaline earth metal halides (MX₂), where M = Mg²⁺, Ca²⁺, Sr²⁺, Ba²⁺ and X = Cl^-, Br^-, I^-, which can be used in bimolecular simulations and which are compatible with the Simple Point Charge/Extended (SPC/E) water model. The force field parameters are specifically developed to reproduce the experimental Kirkwood-Buff integrals for aqueous solutions and thereby the experimental activity derivatives, partial molar volumes, and excess coordination numbers. This ensures that a reasonable balance between ion-ion, ion-water, and water-water distributions is obtained. However, this requires a scaling of the cation to water oxygen interaction strength in order to accurately reproduce the integrals. The scaling factors developed for chloride salts are successfully transferable to the bromide and iodide salts. Use of these new models leads to reasonable diffusion constants and dielectric decrements. However, the performance of the models decreases with increasing salt concentration (>4m), and simulations of the pure crystals exhibited unstable behavior.

A transferable ab initio based force field for aqueous ions

We present a new polarizable force field for aqueous ions (Li(+), Na(+), K(+), Rb(+), Cs(+), Mg(2 +), Ca(2 +), Sr(2 +), and Cl(-)) derived from condensed phase ab initio calculations. We use maximally localized Wannier functions together with a generalized force and dipole-matching procedure to determine the whole set of parameters. Experimental data are then used only for validation purposes and a good agreement is obtained for structural, dynamic, and thermodynamic properties. The same procedure applied to crystalline phases allows to parametrize the interaction between cations and the chloride anion. Finally, we illustrate the good transferability of the force field to other thermodynamic conditions by investigating concentrated solutions.

Sr(II) in Water: A Labile Hydrate with a Highly Mobile Structure

Despite the large number of experimental as well as theoretical investigations available in the literature, some properties of the hydration structure of Sr(II), for example, the coordination number, are still ambiguous. The presented molecular dynamics study based on a most suitable ab initio QM/MM protocol allowed a detailed investigation of structural and dynamical properties of this hydrate, which shows a considerable degree of internal flexebility as well as ligand mobility within the first shell. Despite the high computational effort an exceptionally long QM/MM simulation had to be carried out to obtain sufficient information to investigate first shell ligand exchange reactions.

An anomalous x-ray diffraction study of the hydration structures of Cs+ and I− in concentrated solutions

Anomalous x-ray diffraction experiments were carried out on concentrated aqueous solutions of sodium iodide (6 molal) and cesium iodide (3 molal). Data were gathered at two energies below the absorption edges of the Cs+ and I- ions in order to avoid contributions from fluorescence. The statistics and quality of the raw data were improved by the use of a focusing analyzer crystal. Differences were taken between the data sets and used to calculate the hydration structures of Cs+ and I-. The structures found are more complex than anticipated for such large ions with relatively low charge densities and show evidence of ion-pair formation in both solutions. A two-Gaussian fit to the Cs+ data gives information about the Cs+-O and Cs+-I- correlations. The central position of the Gaussian representing the Cs+-O was fixed at 3.00 A, that is, the maximum of this contribution. The other parameters were allowed to vary freely, giving a Cs+-I- distance of 3.84+/-0.05 A and coordination numbers of 7.9 and 2.7, respectively, for the Cs+-O and Cs+-I- correlations. The results on the structure of I- in the 6 molal NaI aqueous solution were also fitted to a model based on Gaussians; this gives correlations for I- -O and I- -Na+ at 3.17+/-0.06 and 3.76+/-0.06 A with respective coordination numbers of 8.8 and 1.6. The structure of I- in the 3 molal CsI solution shows overlapping contributions due to I- -H, I- -O, and I- -Cs+. The best Gaussian fit gives two peaks centered at 3.00+/-0.08 and 3.82+/-0.04 A and shows that the latter two correlations are unresolved. The hydration structures are compared with those of other alkali and halide ions. The results are also found to be in good agreement with those obtained from standard x-ray diffraction and computer simulation.

Free Energy of Adsorption of Water and Metal Ions on the {101̄4} Calcite Surface

We calculated the free energy profiles of water and three metal ions (magnesium, calcium, and strontium) adsorbing on the [1014] calcite surface in aqueous solution. The approach uses molecular dynamics with parametrized equations to describe the interatomic forces. The potential model is able to reproduce the interactions between water and the metal ions regardless of whether they are at the mineral surface or in bulk water. The simulations predict that the free energy of adsorption of water is relatively small compared to the enthalpy of adsorption calculated in previous papers. This suggests a large change in entropy associated with the water adsorption on the surface. We also demonstrate that the free energy profile of a metal ion adsorbing on the surface correlates with the solvent density and that the rate of formation of an innersphere complex depends on overcoming a large free energy barrier, which is mainly electrostatic in nature. Furthermore, comparison among the rates of desorption of magnesium, calcium, and strontium from the calcite surface suggests that magnesium has a much lower rate of desorption due to its strong interactions with both water and the surface.