Structure of hydration shell of calcium cation by NMR relaxation, Car-Parrinello molecular dynamics and quantum-chemical calculations

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.

Computational and solubility equilibrium experimental insight into Ca2+-fluoride complexation and their dissociation behaviors in aqueous solutions: implication for the association constant measured using fluoride ion selective electrodes

Although the Ca²⁺-F^- association is of great importance for aqueous environments and industrial systems containing F^-, as well as for defluorination processes, many details of the association solvation structures and behavior remain unclear. Herein, a combination of classical/ab initio molecular dynamics simulations and density functional theory calculations was used to investigate the structure and hydration of CaF_x^2-x (x = 1, 2) and the association/dissociation behavior of Ca²⁺-F^- in aqueous CaF₂ solutions. The primary shell of Ca²⁺ is found to be very flexible in the association of Ca²⁺-F^-, with coordination numbers dynamically oscillating in the range of 6-9, with 6 and 7 being the most favorable. The calculations show that for CaF(H₂O)₁₄⁺, the contact ion pair (CIP) is more favorable and occurs with no energy barrier, whereas the formation of CaF₂(aq.) must overcome a ∼3.6 kJ mol^-1 energy barrier; moreover, the CIP and solvent shared ion pair (SSIP) dynamically coexist for CaF₂(H₂O)₁₄ in aqueous CaF₂ solutions. Calculations for the dissociation process of CaF(H₂O)₆⁺ show a dramatic energy increase going from SSIP to free Ca²⁺ and F^-, ascribed to the surprisingly long-range electrostatic attraction between Ca²⁺ and F^- rather than to special F⋯H interactions. The energy increase results in the estimated association constant of CaF⁺ being larger than that previously measured using fluoride ion selective electrodes. This is attributed to the fact that the latter value might correspond to the ligand reaction of free Ca²⁺ and F^- to form the Ca²⁺-F^- SSIP. The combination of these results with CaF_2(s) solubility measurements suggests that the higher-order Ca²⁺-F^- complexes are absent in aqueous CaF₂ solutions.

Influence of Divalent Metal Ions on the Precipitation of the Plasma Protein Fibrinogen

Fibrinogen nanofibers are very attractive biomaterials to mimic the native blood clot architecture. Previously, we reported the self-assembly of fibrinogen nanofibers in the presence of monovalent salts and have now studied how divalent salts influence fibrinogen precipitation. Although the secondary fibrinogen structure was significantly altered with divalent metal ions, morphological analysis revealed exclusively smooth fibrinogen precipitates. In situ monitoring of the surface roughness facilitated predicting the tendency of various salts to form fibrinogen fibers or smooth films. Analysis of the chemical composition revealed that divalent salts were removed from smooth fibrinogen films upon rinsing while monovalent Na⁺ species were still present in fibrinogen fibers. Therefore, we assume that the decisive factor controlling the morphology of fibrinogen precipitates is direct ion-protein contact, which requires disruption of the ion-surrounding hydration shells. We conclude that in fibrinogen aggregates, this mechanism is effective only for monovalent ions, whereas divalent ions are limited to indirect fibrinogen adsorption.

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.

Multivalent Ions as Reactive Crosslinkers for Biopolymers-A Review

Many biopolymers exhibit a strong complexing ability for multivalent ions. Often such ions form ionic bridges between the polymer chains. This leads to the formation of ionic cross linked networks and supermolecular structures, thus promoting the modification of the behavior of solid and gel polymer networks. Sorption of biopolymers on fiber surfaces and interfaces increases substantially in the case of multivalent ions, e.g., calcium being available for ionic crosslinking. Through controlled adsorption and ionic crosslinking surface modification of textile fibers with biopolymers can be achieved, thus altering the characteristics at the interface between fiber and surrounding matrices. A brief introduction on the differences deriving from the biopolymers, as their interaction with other compounds, is given. Functional models are presented and specified by several examples from previous and recent studies. The relevance of ionic crosslinks in biopolymers is discussed by means of selected examples of wider use.

Enhanced Formation of Solvent-Shared Ion Pairs in Aqueous Calcium Perchlorate Solution toward Saturated Concentration or Deep Supercooling Temperature and Its Effects on the Water Structure

As a candidate of Martian salts, calcium perchlorate [Ca(ClO₄)₂] has the potential to stabilize liquid water on the Martian surface because of its hygroscopicity and low freezing temperature when forming aqueous solution. These two properties of electrolytes in general have been suggested to result from the specific cation-anion-water interaction (ion pairing) that interrupts the structure of solvent water. To investigate how this concentration-dependent and temperature-dependent ion pairing process in aqueous Ca(ClO₄)₂ solution leads to its high hygroscopic property and the extreme low eutectic temperature, we have conducted two sets of experiments. First, the effects of concentration on aqueous calcium perchlorate from 3 to 7.86 m on ion pairing were investigated using Raman spectroscopy. Deconvolution of the Raman symmetric stretching band (ν₁) of ClO₄^- showed the enhanced formation of solvent-shared ion pairs upon increasing salt concentration at room temperature. We have confirmed that the low tendency of forming contact ion pairs in concentrated solution contributes to the high hygroscopicity of the salt. Second, the near eutectic samples were studied as a function of temperature by both combined differential scanning calorimetry-Raman spectroscopic experiments and in situ X-ray diffraction. The number of solvent-shared ion pairs was found to increase with decreasing temperature when cooled below the temperature of maximum density of the solution, driven by a change in water toward an ice-like structure in the supercooled regime. The massive presence of solvent-shared ion pairs in turn limits the development of the long-range order in the tetrahedral networks of water molecules, which is responsible for the extremely low eutectic point and deep supercooling effects observed in the Ca(ClO₄)₂-H₂O system.

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.