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

Comparative Investigation of the Microstructure of MgCl2 Aqueous Solutions Using Different X-ray Scattering Sources, Raman Spectroscopy, and Atomistic Simulations

Aqueous solutions of magnesium chloride (MgCl2(aq)) are often used to test advances in the theory of electrolyte solutions because they are considered an ideal strong 2:1 electrolyte. However, there is evidence that some ion association occurs in these solutions, even at low concentrations. Even a small ion-pairing constant can have a significant impact on the chemical speciation of ions, so it is important to determine whether ion pairing actually occurs. In this study, MgCl2(aq) with concentrations ranging from 1 to 35% was studied using three methods: X-ray scattering (XRS) with the Shanghai Synchrotron Radiation Facility (SSRF) and silver-anode laboratory sources, Raman spectroscopy, and molecular dynamics (MD) simulations with the COMPASS-II and Madrid force fields. XRS results were analyzed in the framework of PDF theory to obtain the reduced structure function F(Q) and the reduced pair distribution function G(r). The F(Q) values from synchrotron radiation and laboratory sources both showed that the tetrahedral hydrogen bonds in bulk water were destroyed with the increased MgCl2 concentration. The results of G(r) indicated that the main peaks centered at 2.05 and 2.80 Å can be ascribed to the interactions of Mg-O and O-O, respectively. The peak at 3.10 Å is attributed to the combined effect of O-O and Cl-O. By comparing the structural information on MgCl2 solution obtained from the two light sources, it was found that both SSRF and silver-anode laboratory sources can reflect the above-mentioned structural information on MgCl2 solution. The radial distribution function (RDF) obtained from MD simulations of MgCl2 solutions assigned the peaks at 2.0, 2.8, and 3.2 Å to the Mg-O, O-O, and Cl-O interatomic pairs, respectively. The decrease in the O-O coordination number confirms that the hydrogen-bonding network of water is disrupted by increasing MgCl2 observed by X-ray scattering. The proportion of Mg-Cl contact ion pairs gradually increases with MgCl2 concentration as does the coordination number. Raman spectroscopy results show that the bond type changes from double donor double acceptor (DDAA) to single donor-single acceptor (DA) with increasing concentration, providing explicit details of the hydrogen-bond evolution in the aqueous solution.

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.

Impact of Solution {Ba2+}:{SO42-} on Charge Evolution of Forming and Growing Barite (BaSO4) Crystals: A ζ-Potential Measurement Investigation

The impact of solution stoichiometry on formation of BaSO4 (barite) crystals and the development of surface charge was investigated at various predefined stoichiometries (raq = 0.01, 0.1, 1, 10, and 100, where raq = {Ba2+}:{SO42-}). Synthesis experiments and zeta potential (ζ-potential) measurements were conducted at a fixed initial degree of supersaturation (Ωbarite = 1000, where Ωbarite = {Ba2+}{SO42-}/Ksp), at circumneutral pH of ∼6, 0.02 M NaCl, and ambient temperature and pressure. Mixed-mode measurement-phase analysis light scattering (M3-PALS) showed that the particles stayed negative for raq < 1 during barite crystal formation and positive for raq > 1. At raq = 1, two populations with a positive or negative ζ-potential prevailed for ∼2.5 h before a population with a circumneutral ζ-potential (-10 to +10 mV) remained. We relate the observations of particle charge evolution to particle size and morphology evolution under the experimental conditions. Furthermore, we showed that the ζ-potential became more negative when the pH was increased for every raq. In addition, our results demonstrated that the type of monovalent background electrolyte did not influence the ζ-potential of barite crystals significantly, although NaCl showed slightly different behavior compared to KCl and NaNO3. Our results show the important role of surface charge (evolution) during ionic crystal formation under nonstoichiometric conditions. Moreover, our combined scanning electron microscopy and ζ-potential results imply that the surface charge during particle formation can be influenced by solution stoichiometry, besides the pH and ionic strength, and may aid in predicting the fate of barite in environmental settings and in understanding and improving industrial barite (surface chemistry) processes.

Effect of Divalent Cations and Other Ions on the Tetrahydrofuran Crystal Inhibition of Quaternary Ammonium Salts-Relevance to the Efficiency of Gas Hydrate Quaternary Anti-agglomerants

Gas hydrate anti-agglomerants (AAs) are a class of low-dosage hydrate inhibitor that are used to prevent plugging of gas hydrates in oil and condensate upstream flow lines. Industrial AAs are mostly cationic surfactants which are "hydrate-philic", i.e., they are designed to interact with and modify gas hydrate crystal growth. Tetrahydrofuran (THF) hydrate crystal growth studies have been used for many years to determine useful functional groups to incorporate into AA surfactants. In particular, quaternary ammonium and phosphonium salts with optimized alkyl groups show good THF crystal growth inhibition, which is a key property for AAs. AAs are often screened and tested in model brines containing sodium chloride despite the produced water containing various divalent cations. Recent studies have shown that AAs performed better when tested in brines containing both sodium and calcium ions rather than just sodium ions. Here, we present THF hydrate crystal growth studies on quaternary ammonium and phosphonium salts and other related molecules including guanidinium salts and amine oxides. Tests were carried out with a variety of cations including sodium, calcium, magnesium, and lithium at identical pre-determined subcooling, in order to investigate the effect of the ion size and charge density on the crystal growth inhibition. We also investigate the effect of using the more polarizable iodide ions compared to chloride ions. Our results show that crystal growth inhibition in solutions with calcium ions is somewhat greater than that with sodium ions, in agreement with past studies on the effect of AA performance with mono- and divalent cations. However, the variation does not seem to be primarily related to the charge density and polarizing ability of the cations. This study therefore provides evidence that AAs should be tested in brines containing all the ions present in the produced water and not just sodium chloride brine.

New insights into the early stage nucleation of calcium carbonate gels by reactive molecular dynamics simulations

The precipitation of calcium carbonate (CaCO₃) is a key mechanism in carbon capture applications relying on mineralization. In that regard, Ca-rich cementitious binders offer a unique opportunity to act as a large-scale carbon sink by immobilizing CO₂ as calcium carbonate by mineralization. However, the atomistic mechanism of calcium carbonate formation is still not fully understood. Here, we study the atomic scale nucleation mechanism of an early stage amorphous CaCO₃ gel based on reactive molecular dynamics (MD) simulations. We observe that reactive MD offers a notably improved description of this reaction as compared to classical MD, which allows us to reveal new insights into the structure of amorphous calcium carbonate gels and formation kinetics thereof.

Calcite Kinks Grow via a Multistep Mechanism

The classical model of crystal growth assumes that kinks grow via a sequence of independent adsorption events where each solute transitions from the solution directly to the crystal lattice site. Here, we challenge this view by showing that some calcite kinks grow via a multistep mechanism where the solute adsorbs to an intermediate site and only transitions to the lattice site upon the adsorption of a second solute. We compute the free energy curves for Ca and CO₃ ions adsorbing to a large selection of kink types, and we identify kinks terminated both by Ca ions and by CO₃ ions that grow in this multistep way.

A Database of Solution Additives Promoting Mg2+ Dehydration and the Onset of MgCO3 Nucleation

Formed via aqueous carbonation of Mg²⁺ ions, the crystallization of magnesite (MgCO₃) is a promising route to carbon capture and reuse, albeit limited by the slow precipitation of MgCO₃. Although magnesite is naturally abundant, forming at low temperature conditions, its industrial production is an energy-intensive process due to the temperatures required to prevent the formation of hydrated phases. The principal difficulty in aqueous conditions arises from the very strong Mg²⁺···H₂O interaction, with high barriers to Mg²⁺ dehydration. Using atomistic simulations, we have investigated the influence of 30 additive anions (X ^n-, n = 1-3), ranging from simple halides to more complex molecules, on the first two steps of MgCO₃ aggregation from solution, as follows: Mg²⁺ dehydration and subsequent prenucleative Mg²⁺···CO₃ ^2- pairing. We have computed the thermodynamic stabilities of solvent shared ion pairs (Mg²⁺···H₂O···X ^n-) and contact ion pairs (Mg²⁺···X ^n-) to reveal the propensity of solution additives to inhibit or promote Mg²⁺···CO₃ ^2- formation. We have determined the stabilization of undercoordinated hydrated Mg²⁺ states with a vacant coordination site to which CO₃ ^2- can bind, subsequently initiating MgCO₃ nucleation or Mg²⁺ incorporation into the crystal lattice. Extensive molecular dynamics simulations of electrolyte solutions containing Na₂CO₃ with different sources of Mg²⁺ (i.e., MgCl₂, MgSO₄, and Mg(CH₃COO)₂) further show that the degree of dehydration of Mg²⁺ and the structure of prenucleation MgCO₃ clusters change depending on the counterion identity. Through a fundamental understanding of the role of solution additives in the mechanism of Mg²⁺ dehydration, our results help to rationalize previously reported experimental observation of the effect of solvation environments on the growth of magnesite. This understanding may contribute to identifying the solution composition and conditions that could promote the low-temperature CO₂ conversion into MgCO₃ at industrially relevant scales.

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na⁺, Ca²⁺, and Mg²⁺ cations and the CO₃^2- anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me-O-C (Me = metal) three-body valence angle parameters and Me-C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na⁺ predominantly forms carbonates due to weaker angular strain, while Ca²⁺ and Mg²⁺ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. The nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca²⁺ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na⁺ and Mg²⁺ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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.

Coordination Numbers of Bivalent Ions in Organic Solvents

Abstract

Molecular dynamic models are created for properties of bivalent ions in organic solvents. It is shown that molecules of the considered solvents bound to ions via oxygen atoms. A theoretical model is developed that describes the ion coordination number. The coordination number in this model is determined by the ratio between the sizes of the ion and the atom organic molecule bound to it. It is shown that the coordination number depends weakly on the solvent and strongly on the type of ion. A value of 0.13 nm is obtained for the effective size of an oxygen atom bound to a bivalent ion. The constructed theoretical model agrees with the results from molecular dynamic calculations and the available experimental data.

Strontium increasing calcium accessibility from calcium citrate

Strontium chloride added to aqueous suspensions of metastable calcium citrate tetrahydrate increased calcium ion activity measured electrochemically without transition of metastable tetrahydrate to stable calcium citrate hexahydrate as shown by DSC. Calcium activity increase was explained by lower solubility of strontium citrate pentahydrate formed (8.9 × 10^-4 M at 25 °C) increasing with temperature compared to calcium citrate tetrahydrate (1.6 × 10^-3 M) decreasing with temperature. Strontium binding to citrate was found endothermic, ΔH⁰ = 45 kJ∙mol^-1 at 25 °C, while calcium binding shows variation from ΔH⁰ = 94 kJ∙mol^-1 at 10 °C becoming exothermic above physiological temperature with ΔH⁰ = -9 kJ∙mol^-1 at 45 °C as determined from temperature and concentration variation in electric conductivity. These differences in solution thermodynamics and pH effect on complex formation between calcium and strontium citrate are discussed in relation to biomineralization.

First Steps towards Understanding the Non-Linear Impact of Mg on Calcite Solubility: A Molecular Dynamics Study

Magnesium (Mg2+) is one of the most common impurities in calcite and is known to have a non-linear impact on the solubility of magnesian calcites. Using molecular dynamics (MD), we observed that Mg2+ impacts overall surface energies, local free energy profiles, interfacial water density, structure and dynamics and, at higher concentrations, it also causes crystal surface deformation. Low Mg concentrations did not alter the overall crystal structure, but stabilised Ca2+ locally and tended to increase the etch pit nucleation energy. As a result, Ca-extraction energies over a wide range of 39 kJ/mol were observed. Calcite surfaces with an island were less stable compared to flat surfaces, and the incorporation of Mg2+ destabilised the island surface further, increasing the surface energy and the calcium extraction energies. In general, Ca2+ is less stable in islands of high Mg2+ concentrations. The local variation in free energies depends on the amount and distance to nearest Mg in addition to local disruption of interfacial water and the flexibility of surface carbonate ions to rotate. The result is a complex interplay of these characteristics that cause variability in local dissolution energies. Taken together, these results illustrate molecular scale processes behind the non-linear impact of Mg2+ concentration on the solubility of magnesium-bearing calcites.

Slow lactate gluconate exchange in calcium complexes during precipitation from supersaturated aqueous solutions

Saturated solutions of calcium l-lactate in water or in deuterium oxide continuously dissolve calcium l-lactate by addition of solid sodium d-gluconate and become strongly supersaturated in calcium d-gluconate due to no or slow precipitation. The quantification of total dissolved calcium allied with the calcium complexes equilibrium constants allowed an ion speciation, which shows an initial non-thermal and spontaneous supersaturation of more than a factor of 50 at 25 °C only slowly decreasing after initiation of precipitation of calcium d-gluconate after a lag phase of several hours. A mathematical model is proposed, based on numerical solution of coupled differential equations of dynamics of l-lactate and d-gluconate exchange during the lag phase for precipitation and during precipitation. A slow exchange of l-lactate coordinated to calcium with d-gluconate is indicated with a time constant of 0.20 h^-1 in water and of 0.15 h^-1 in deuterium oxide and a kinetic deuterium/hydrogen isotope effect of 1.25. Such spontaneous non-thermal supersaturation and slow ligand exchange with a pseudo first order equilibration process with a half-life of 3.5 h in water for calcium hydroxycarboxylates can help to understand the higher calcium bioavailability from calcium hydroxycarboxylates compared to simple salts.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

ReaxFF molecular dynamics simulations of electrolyte-water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal-chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li⁺ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na⁺ and K⁺ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

Understanding Calcium-Dependent Conformational Changes in S100A1 Protein: A Combination of Molecular Dynamics and Gene Expression Study in Skeletal Muscle

The S100A1 protein, involved in various physiological activities through the binding of calcium ions (Ca²⁺), participates in several protein-protein interaction (PPI) events after Ca²⁺-dependent activation. The present work investigates Ca²⁺-dependent conformational changes in the helix-EF hand-helix using the molecular dynamics (MD) simulation approach that facilitates the understanding of Ca²⁺-dependent structural and dynamic distinctions between the apo and holo forms of the protein. Furthermore, the process of ion binding by inserting Ca²⁺ into the bulk of the apo structure was simulated by molecular dynamics. Expectations of the simulation were demonstrated using cluster analysis and a variety of structural metrics, such as interhelical angle estimation, solvent accessible surface area, hydrogen bond analysis, and contact analysis. Ca²⁺ triggered a rise in the interhelical angles of S100A1 on the binding site and solvent accessible surface area. Significant configurational regulations were observed in the holo protein. The findings would contribute to understanding the molecular basis of the association of Ca²⁺ with the S100A1 protein, which may be an appropriate study to understand the Ca²⁺-mediated conformational changes in the protein target. In addition, we investigated the expression profile of S100A1 in myoblast differentiation and muscle regeneration. These data showed that S100A1 is expressed in skeletal muscles. However, the expression decreases with time during the process of myoblast differentiation.

Reconsidering Calcium Dehydration as the Rate-Determining Step in Calcium Mineral Growth

The dehydration of cations is generally accepted as the rate-limiting step in many processes. Molecular dynamics (MD) can be used to investigate the dynamics of water molecules around cations, and two different methods exist to obtain trajectory-based water dehydration frequencies. Here, these two different post-processing methods (direct method versus survival function) have been implemented to obtain calcium dehydration frequencies from a series of trajectories obtained using a range of accepted force fields. None of the method combinations reproduced the commonly accepted experimental water exchange frequency of 10-8.2 s-1. Instead, our results suggest much faster water dynamics, comparable with more accurate ab initio MD simulations and with experimental values obtained using neutron scattering techniques. We obtained the best agreement using the survival function method to characterize the water dynamics, and we show that different method combinations significantly affect the outcome. Our work strongly suggests that the fast water exchange kinetics around the calcium ions is not rate-limiting for reactions involving dissolved/solvated calcium. Our results further suggest that, for alkali and most of the earth alkali metals, mechanistic rate laws for growth, dissolution, and adsorption, which are based on the principle of rate-limiting cation dehydration, need careful reconsideration.

Do water's electrons care about electrolytes?

Ions have a profound effect on the geometrical structure of liquid water and an aqueous environment is known to change the electronic structure of ions. Here we combine photoelectron spectroscopy measurements from liquid microjets with molecular dynamical and quantum chemical calculations to address the reverse question, to what extent do ions affect the electronic structure of liquid water? We study aqueous solutions of sodium iodide (NaI) over a wide concentration range, from nearly pure water to 8 M solutions, recording spectra in the 5 to 60 eV binding energy range to include all water valence and the solute Na+ 2p, I- 4d, and I- 5p orbital ionization peaks. We observe that the electron binding energies of the solute ions change only slightly as a function of electrolyte concentration, less than 150 ± 60 meV over an ∼8 M range. Furthermore, the photoelectron spectrum of liquid water is surprisingly mildly affected as we transform the sample from a dilute aqueous salt solution to a viscous, crystalline-like phase. The most noticeable spectral changes are a negative binding energy shift of the water 1b2 ionizing transition (up to -370 ± 60 meV) and a narrowing of the flat-top shape water 3a1 ionization feature (up to 450 ± 90 meV). A novel computationally efficient technique is introduced to calculate liquid-state photoemission spectra using small clusters from molecular dynamics (MD) simulations embedded in dielectric continuum. This theoretical treatment captured the characteristic positions and structures of the aqueous photoemission peaks, reproducing the experimentally observed narrowing of the water 3a1 feature and weak sensitivity of the water binding energies to electrolyte concentration. The calculations allowed us to attribute the small binding energy shifts to ion-induced disruptions of intermolecular electronic interactions. Furthermore, they demonstrate the importance of considering concentration-dependent screening lengths for a correct description of the electronic structure of solvated systems. Accounting for electronic screening, the calculations highlight the minimal effect of electrolyte concentration on the 1b1 binding energy reference, in accord with the experiments. This leads us to a key finding that the isolated, lowest-binding-energy, 1b1, photoemission feature of liquid water is a robust energetic reference for aqueous liquid microjet photoemission studies.

Adhesive forces between two cleaved calcite surfaces in NaCl solutions: The importance of ionic strength and normal loading

The mechanical strength of calcite bearing rocks is influenced by pore fluid chemistry due to the variation in nano-scale surface forces acting at the grain contacts or close to the fracture tips. The adhesion of two contacting surfaces, which affects the macroscopic strength of the material, is not only influenced by the fluid chemistry but also by the surface topography. In this paper, we use Atomic Force Microscope (AFM) to measure the interfacial forces between two freshly cleaved calcite surfaces in CaCO₃-saturated solutions with varying NaCl concentration. We show that calcite contacts become stronger with increasing NaCl concentration (>100 mM), as a result of progressively weaker secondary hydration and increasing attraction due to instantaneous ion-ion correlation. Moreover, we discuss the effect of normal applied force (F_n) and surface roughness on the measured adhesion forces (F_ad). We show that the measured pull-off force (adhesion) is linearly correlated with the magnitude of F_n, where an increase in applied force results in increased adhesion. This is attributed to a larger number of contacting surface asperities and thus increase in real contact area and the contact-bond strength. We discuss that the possible variation in local topography at contacts, together with strong dependence on ionic strength of the solution, can explain the inconsistent behavior of calcite rocks in NaCl solutions.

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.

On the selectivity of the NaChBac channel: an integrated computational and experimental analysis of sodium and calcium permeation

Ion channel selectivity is essential for their function, yet the molecular basis of a channel's ability to select between ions is still rather controversial. In this work, using a combination of molecular dynamics simulations and electrophysiological current measurements we analyze the ability of the NaChBac channel to discriminate between calcium and sodium. Our simulations show that a single calcium ion can access the Selectivity Filter (SF) interacting so strongly with the glutamate ring so as to remain blocked inside. This is consistent with the tiny calcium currents recorded in our patch-clamp experiments. Two reasons explain this scenario. The first is the higher free energy of ion/SF binding of Ca²⁺ with respect to Na⁺. The second is the strong electrostatic repulsion exerted by the resident ion that turns back a second potentially incoming Ca²⁺, preventing the knock-on permeation mechanism. Finally, we analyzed the possibility of the Anomalous Mole Fraction Effect (AMFE), i.e. the ability of micromolar Ca²⁺ concentrations to block Na⁺ currents. Current measurements in Na⁺/Ca²⁺ mixed solutions excluded the AMFE, in agreement with metadynamics simulations showing the ability of a sodium ion to by-pass and partially displace the resident calcium. Our work supports a new scenario for Na⁺/Ca²⁺ selectivity in the bacterial sodium channel, challenging the traditional notion of an exclusion mechanism strictly confining Ca²⁺ ions outside the channel.

Molecular insight into the nanoconfined calcite-solution interface

Little is known about the influence of nanoconfinement on calcium carbonate mineralization. Here, colloidal probe atomic force microscopy is used to confine the calcite-solution interface with a silica microsphere and to measure Derjaguin-Landau-Verwey-Overbeek (DLVO) and non-DLVO forces as a function of the calcium concentration, also after charge reversal of both surfaces occurs. Through the statistical analysis of the oscillatory component of a strong hydration force, the subnanometer interfacial structure of the confined atomically flat calcite is resolved in aqueous solution. By applying a mechanical work, both water and hydrated counterions are squeezed out from the nanoconfined solution, leaving the calcite surface more negatively charged than the analogous unconfined surfaces. Layer size and applied work allow a distinction between the hydration states of the counterions in the Stern layer; we propose counterions to be inner- and outer-sphere calcium ions, with a population of inner-sphere calcium ions larger than on unconfined calcite surfaces. It is also shown that the composition of the nanoconfined solution can be tuned by varying calcium concentration. This is a fundamental study of DLVO and hydration forces, and of their connection, on atomically flat calcite. More broadly, our work scrutinizes the greatly unexplored relation between surface science and confined mineralization, with implications for diverse areas of inquiry, such as nanoconfined biomineralization, CO₂ sequestration in porous aquifers, and pressure solution and crystallization in confined hydrosystems.

Simulations reveal the role of composition into the atomic-level flexibility of bioactive glass cements

Bioactive glass ionomer cements (GICs), the reaction product of a fluoro-alumino-silicate glass and polyacrylic acid, have been in effective use in dentistry for over 40 years and more recently in orthopaedics and medical implantation. Their desirable properties have affirmed GIC's place in the medical materials community, yet are limited to non-load bearing applications due to the brittle nature of the hardened composite cement, thought to arise from the glass component and the interfaces it forms. Towards helping resolve the fundamental bases of the mechanical shortcomings of GICs, we report the 1st ever computational models of a GIC-relevant component. Ab initio molecular dynamics simulations were employed to generate and characterise three fluoro-alumino-silicate glasses of differing compositions with focus on resolving the atomic scale structural and dynamic contributions of aluminium, phosphorous and fluorine. Analyses of the glasses revealed rising F-content leading to the expansion of the glass network, compression of Al-F bonding, angular constraint at Al-pivots, localisation of alumino-phosphates and increased fluorine diffusion. Together, these changes to the structure, speciation and dynamics with raised fluorine content impart an overall rigidifying effect on the glass network, and suggest a predisposition to atomic-level inflexibility, which could manifest in the ionomer cements they form.

A DFT-based comparative equilibrium study of thermal dehydration and hydrolysis of CaCl2 hydrates and MgCl2 hydrates for seasonal heat storage

An equilibrium based comparative study of CaCl₂ and MgCl₂ hydrates reveals that CaCl₂ hydrates have better resistance against hydrolysis compared to MgCl₂ hydrates; therefore it can improve the durability of chloride based salt hydrates.

Ab Initio Studies of Calcium Carbonate Hydration

Ab initio simulations of large hydrated calcium carbonate clusters are challenging due to the existence of multiple local energy minima. Extensive conformational searches around hydrated calcium carbonate clusters (CaCO3·nH2O for n = 1-18) were performed to find low-energy hydration structures using an efficient combination of Monte Carlo searches, density-functional tight binding (DFTB+) method, and density-functional theory (DFT) at the B3LYP level, or Møller-Plesset perturbation theory at the MP2 level. This multilevel optimization yields several low-energy structures for hydrated calcium carbonate. Structural and energetics analysis of the hydration of these clusters revealed a first hydration shell composed of 12 water molecules. Bond-length and charge densities were also determined for different cluster sizes. The solvation of calcium carbonate in bulk water was investigated by placing the explicitly solvated CaCO3·nH2O clusters in a polarizable continuum model (PCM). The findings of this study provide new insights into the energetics and structure of hydrated calcium carbonate and contribute to the understanding of mechanisms where calcium carbonate formation or dissolution is of relevance.

Monitoring the intramolecular charge transfer process in the Z907 solar cell sensitizer: a transient Vis and IR spectroscopy and ab initio investigation

The excited state dynamics of Z907 in solution and on semiconductor substrates has been studied with ultrafast UV/Vis and IR spectroscopy and DFT/TDDFT calculations.