Structural and Dynamical Nature of Hydration Shells of the Carbonate Ion in Water: An Ab Initio Molecular Dynamics Study

Agglomeration Drives the Reversed Fractionation of Aqueous Carbonate and Bicarbonate at the Air-Water Interface

In the course of our investigations of the adsorption of ions to the air-water interface, we previously reported the surprising result that doubly charged carbonate anions exhibit a stronger surface affinity than singly charged bicarbonate anions. In contrast to monovalent, weakly hydrated anions, which generally show enhanced concentrations in the interfacial region, multivalent (and strongly hydrated) anions are expected to show a much weaker surface propensity. In the present work, we use resonantly enhanced deep-UV second-harmonic generation spectroscopy to measure the Gibbs free energy of adsorption of both carbonate (CO32-) and bicarbonate (HCO3-) anions to the air-water interface. Contrasting the predictions of classical electrostatic theory and in support of our previous findings from X-ray photoelectron spectroscopy, we find that carbonate anions do indeed exhibit much stronger surface affinity than do the bicarbonate anions. Extensive computer simulations reveal that strong ion pairing of CO32- with the Na+ countercation in the interfacial region results in the formation of near-neutral agglomerate clusters, consistent with a theory of interfacial ion adsorption based on hydration free energy and capillary waves. Simulated X-ray photoelectron spectra predict a 1 eV shift in the carbonate spectra compared to that of bicarbonate, further confirming our experiments. These findings not only advance our fundamental understanding of ion adsorption chemistry but also impact important practical processes such as ocean acidification, sea-spray aerosol chemistry, and mammalian respiration physiology.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg²⁺ and SO₄ ^2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO₄ ^2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO₄ ^2- ions around 1100 cm^-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO₄ and for Na₂SO₄ in the presence of Mg²⁺ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na₂SO₄ solutions in the absence of Mg²⁺ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO₄ ^2--Mg²⁺ ion pairs.

Time-resolved terahertz-Raman spectroscopy reveals that cations and anions distinctly modify intermolecular interactions of water

The solvation of ions changes the physical, chemical and thermodynamic properties of water, and the microscopic origin of this behaviour is believed to be ion-induced perturbation of water's hydrogen-bonding network. Here we provide microscopic insights into this process by monitoring the dissipation of energy in salt solutions using time-resolved terahertz-Raman spectroscopy. We resonantly drive the low-frequency rotational dynamics of water molecules using intense terahertz pulses and probe the Raman response of their intermolecular translational motions. We find that the intermolecular rotational-to-translational energy transfer is enhanced by highly charged cations and is drastically reduced by highly charged anions, scaling with the ion surface charge density and ion concentration. Our molecular dynamics simulations reveal that the water-water hydrogen-bond strength between the first and second solvation shells of cations increases, while it decreases around anions. The opposite effects of cations and anions on the intermolecular interactions of water resemble the effects of ions on the stabilization and denaturation of proteins.

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.

Kosmotropic Electrolyte (Na2CO3, NaF) Perturbs the Air/Water Interface through Anion Hydration Shell without Forming a Well-Defined Electric Double Layer

The ion-driven electric double layer (EDL) and the structural transformation of interfacial water are implicated in unusual reaction kinetics at the air/water interface. By combining heterodyne-detected vibrational sum frequency generation (HD-VSFG) with differential spectroscopy involving simultaneous curve fitting (DS-SCF) analysis, we retrieve electrolyte (Na₂CO₃ and NaF)-correlated OH-stretch spectra of water at the air/water interface. Vibrational mapping of the perturbed interfacial water with the hydration shell spectra (obtained by DS-SCF analysis of Raman spectra) of the corresponding anion discloses that the kosmotropic electrolytes do not form well-defined EDL at the air/water interface. Instead, the interfacial water forms a stronger hydrogen-bond with the surface-expelled anions (CO₃^2- and F^-) and becomes more inhomogeneous than the pristine air/water interface. Together, the results reveal that the perturbation of interfacial water by the kosmotropic electrolyte is a "local phenomenon" confined within the hydration shell of the surface-expelled anion.

Insignificant Effect of Temperature on the Structure and Angular Jumps of Water near a Hydrophobic Cation

The ambiguity in the behavior of water molecules around hydrophobic solutes is a matter of interest for many studies. Motivated by the earlier results on the dynamics of water molecules around tetramethylammonium (TMA) cation, we present the effect of temperature on the structure and angular jumps of water due to hydrophobicity using first principles molecular dynamics simulations. The average intermolecular distance between the central oxygen and four nearest neighbors is found to be the highest for water molecules in the solvation shell of TMA at 400 K, followed by the same at 330 K. The hydrogen bond (HB) donor-acceptor count, HB per water molecule, and tetrahedral order parameter suggests the loss of tetrahedrality in the solvation shell. Elevated temperature affects the tetrahedral parameter in local regions. The HB jump mechanism is studied for methyl hydrogen and water molecules in the solvation shell. Observations hint at the presence of dangling water molecules in the vicinity of the hydrophobic cation, and no evidence is found for the enhanced structural ordering of nearby water molecules.

On the interference of urea with CO2/CO32- chemistry of cellulose model solutions in NaOH(aq)

The CO₂/CO₃^2- chemistry of the cellulose/NaOH(aq) solutions has been recently reported to comprise a CO₂ incorporation through formation of a transient cellulose carbonate intermediate along with cellulose - CO₃^2- interactions. This work explores on molecular interactions arising when this chemistry is brought together with urea, the most common stabiliser of these solutions. ¹H, ¹³C and steady-state heteronuclear Overhauser effect NMR studies on the cellulose analogues (methyl-β-glucopyranoside (β-MeO-Glcp) and microcrystalline cellulose), combined with pH and ATR-FTIR measurements, reveal concurrent interactions of urea with both CO₂ and CO₃^2-- leading to increased uptake of CO₂ and a buffering effect. Yet, regardless of the presence of urea, the route of conversion from CO₂ to CO₃^2-, whether going through reaction with the carbohydrate alkoxides or OH^-, is likely to determine the chemical environment of the formed CO₃^2-. These findings shed a new light on rather overlooked, albeit prominent, interactions in these solutions with the readily absorbed air CO₂, essential for further development and implementation, whether towards regenerated and modified cellulose or CO₂-capturing concepts.

Transport of hydrated nitrate and nitrite ions through graphene nanopores in aqueous medium

Nitrate ( NO 3 - ) and nitrite ( NO 2 - ) ions are naturally occurring inorganic ions that are part of the nitrogen cycle. High doses of these ions in drinking water impose a potential risk to public health. In this work, molecular dynamics simulations are carried out to study the passage of nitrate and nitrite ions from water through graphene nanosheets (GNS) with hydrogen-functionalized narrow pores in presence of an external electric field. The passage of ions through the pores is investigated through calculations of ion flux, and the results are analyzed through calculations of various structural and thermodynamic properties such as the density of ions and water, ion-water radial distribution functions, two-dimensional density distribution functions, and the potentials of mean force of the ions. Current simulations show that the nitrite ions can pass more in numbers than the nitrate ions in a given time through GNS hydrogen-functionalized pore of different geometry. It is found that the nitrite ions can permeate faster than the nitrate ions despite the former having higher hydration energy in the bulk. This can be explained in terms of the competition between the number density of the ions along the pore axis and the free energy barrier calculated from the potential of mean force. Also, an externally applied electric field is found to be important for faster permeation of the nitrite over the nitrate ions. The current study suggests that graphene nanosheets with carefully created pores can be effective in achieving selective passage of ions from aqueous solutions.

Ion Pairing and Multiple Ion Binding in Calcium Carbonate Solutions Based on a Polarizable AMOEBA Force Field and Ab Initio Molecular Dynamics

The speciation of calcium carbonate in water is important to the geochemistry of the world's oceans and has ignited significant debate regarding the mechanism by which nucleation occurs. Here, it is vital to be able to quantify the thermodynamics of ion pairing versus higher order association processes in order to distinguish between possible pathways. Given that it is experimentally challenging to quantify such species, here we determine the thermodynamics for ion pairing and multiple binding of calcium carbonate species using bias-enhanced molecular dynamics. In order to examine the uncertainties underlying these results, we derived a new polarizable force field for both calcium carbonate and bicarbonate in water based on the AMOEBA model to compare against our earlier rigid ion model, both of which are further benchmarked against ab initio molecular dynamics for the ion pair. Both force fields consistently indicate that the association of calcium carbonate ion pairs to form larger species is stable, though with an equilibrium constant that is lower than for ion pairing itself.

Dynamics of Water in the Solvation Shell of an Iodate Ion: A Born-Oppenheimer Molecular Dynamics Study

The iodate ion has an anisotropic structure and charge distribution. It has a pyramidal shape with the iodine atom located at the peak of the pyramid. The water molecules interact differently with the positively charged iodine and the negatively charged oxygen atoms of this anion, giving rise to two distinct solvation shells. In the present study, we have performed ab initio Born-Oppenheimer molecular dynamics simulations to investigate the dynamics of water molecules in the iodine and oxygen solvation shells of the iodate ion and compared the behavior with those of the bulk. The dynamics of water is calculated for both the BLYP and the dispersion-corrected BLYP-D3 functionals at room temperature. The dynamics of water in the solvation shells at higher temperatures of 353 and 330 K has also been investigated for the BLYP and BLYP-D3 functionals, respectively. The hydrogen bond dynamics, vibrational spectral diffusion, orientational and translational diffusion, and residence dynamics of water molecules in the two solvation shells are looked at in the current study. The ion-water hydrogen bond dynamics is found to be somewhat faster than that for water-water hydrogen bonds in the bulk, which can be attributed to a ring-like electron distribution on the iodate oxygens. The dynamical trends are connected to the water structure making/breaking properties of the positively charged iodine and negatively charged oxygen sites of the anion. Furthermore, orientational jumps of the iodate ion and also those of surrounding water molecules which are hydrogen bonded to the oxygen atoms of the iodate ion are also investigated. It is found that the nature of these orientational jumps can be different from those reported earlier for planar polyoxyanions such as the nitrate ion.

Distinctive behavior and two-dimensional vibrational dynamics of water molecules inside glycine solvation shell

We present a first principles molecular dynamics study of a deuterated aqueous solution of a single glycine moiety to explore the structure, dynamics, and two-dimensional infrared spectra of water molecules found in the solvation shell of glycine.

Resonance-driven ion transport and selectivity in prokaryotic ion channels

Ion channels exhibit a remarkably high accuracy in selecting uniquely its associated type of ion. The mechanisms behind ion selectivity are not well understood. Current explanations build mainly on molecular biology and bioinformatics. Here we propose a simple physical model for ion selectivity based on the driven damped harmonic oscillator (DDHO). The driving force for this oscillator is provided by self-organizing harmonic turbulent structures in the dehydrating ionic flow through the ion channel, namely, oscillating pressure waves in one dimension, and toroidal vortices in two and three dimensions. Density fluctuations caused by these turbulences efficiently transmit their energy to aqua ions that resonate with the driving frequency. Consequently, these release their hydration shell and exit the ion channel as free ions. Existing modeling frameworks do not express the required complex spatiotemporal dynamics, hence we introduce a macroscopic continuum model for ionic dehydration and transport, based on the hydrodynamics of a dissipative ionic flow through an ion channel, subject to electrostatic and amphiphilic interactions. This model combines three classical physical fields: Navier-Stokes equations from hydrodynamics, Gauss's law from Maxwell theory, and the convection-diffusion equation from continuum physics. Numerical experiments with mixtures of chemical species of ions in various degrees of hydration indeed reveal the emergence of strong oscillations in the ionic flow that are instrumental in the efficient dehydration and cause a strong ionic jet into the cell. As such, they provide a powerful engine for the DDHO mechanism. Theoretical predictions of the modeling framework match significantly with empirical patch-clamp data. The DDHO standard response curve defines a unique resonance frequency that depends on the mass and charge of the ion. In this way, the driving oscillations act as a selection mechanism that filters out one specific ion. Application of the DDHO model to real ion data shows that this mechanism indeed clearly distinguishes between chemical species and between aqua and bare ions with a large Mahalanobis distance and high oscillator quality. The DDHO framework helps to understand how SNP mutations can cause severe genetic pathologies as they destroy the geometry of the channel and so alter the resonance peaks of the required ion type.

Dynamics of linear molecules in water: Translation-rotation coupling in jump motion driven diffusion

We study by computer simulations, and by theory, the coupled rotational and translational dynamics of three important linear diatomic molecules, namely, carbon monoxide (CO), nitric oxide (NO), and cyanide ion (CN^-) in water. Translational diffusion of these molecules is found to be strongly coupled to their own rotational dynamics which, in turn, are coupled to similar motions of the surrounding water. In particular, we find that coupled orientational jump motions play an important role in all three cases. While CO and NO show similar features, CN^- exhibits certain differences. Our results agree well with the known experimental values of the diffusion coefficient. We examined the validity of hydrodynamic predictions and found them to be inadequate, particularly for rotational diffusion. A mode coupling theory approach is developed and applied to understand the complexity of translation-rotation coupling.

Molecular insight into the wetting behavior and amphiphilic character of cellulose nanocrystals

The study of nanocellulose is a field of growing interest due to its many applications and its use in the development of biocompatible and eco-friendly materials. In spite of the vast number of studies in the field, many questions about the role of the molecular structure in the properties of cellulose are still subject of debate. One of these fundamental questions is the possible amphiphilic nature of cellulose and the relative role of hydrogen bonding and hydrophobic effect on the interactions of cellulose. In this work we present an extensive molecular dynamics simulation study of this question by analyzing the wetting of cellulose with water and organic solvent, its interaction with hydrophilic and hydrophobic ions and its interaction with a protein (human epidermal growth factor, hEGF). We consider two characteristic cellulose crystal planes of Iβ cellulose with very different roughness, different hydrogen bonding capability and different exposure of cellulose hydrophobic groups (the (010) plane which has exposed -OH groups and the (100) plane with buried -OH groups). Our results show that both surfaces are simultaneously hydrophilic and lipophilic, with both surfaces having very similar contact angles. In spite of the global similarity of wetting of both surfaces, the molecular details of wetting are very different and substantial local wetting heterogeneities (which strongly depend on the surface) appear for both solvents. We also observe a weak interaction of both surfaces with hydrophobic and hydrophilic solutes. These weak interactions are attributed to the simultaneous lipophilic and hydrophilic character of both (100) and (010) cellulose surfaces. Interestingly, we found a substantial interaction of both cellulose planes with polar and apolar residues of the hEGF protein.

Computational Study of the pH-Dependent Competition between Carbonate and Thymine Addition to the Guanine Radical

When oligonucleotides are oxidized by carbonate radical, thymine and carbonate can add to guanine radical, yielding either a guanine-thymine cross-link product (G∧T) or 8-oxo-7,8-dehydroguanine (8oxoG) and its further oxidation products such as spiroiminodihydantoin (Sp) and guanidinohydantoin (Gh). The ratio of thymine addition to carbonate addition depends strongly on the pH. Details of the mechanism have been explored by density functional calculations using the ωB97XD/6-31+G(d,p) level of theory with the SMD implicit solvation method, augmented with a few explicit waters. Free energies of intermediates and transition states in aqueous solution have been calculated along the pathways for addition of thymine, CO₃^2-/HCO₃^- and carbonate radical to guanine radical. The pH dependence was examined by using appropriate explicit proton donors/acceptors as computational models for buffers at pH 2.5, 7, and 10. Deprotonation of thymine is required for nucleophilic addition at C8 of guanine radical, and thus is favored at higher pH. The barrier for carbonate radical addition is lower than for bicarbonate or carbonate dianion addition; however, for low concentrations of carbonate radical, the reaction may proceed by addition of bicarbonate/carbonate dianion to guanine radical. Thymine and bicarbonate/carbonate dianion addition are followed by oxidation by O₂, loss of a proton from C8 and decarboxylation of the carbonate adduct. At pH 2.5, guanine radical cation can be formed by oxidization with sulfate radical. Water addition to guanine radical cation is the preferred path for forming 8oxoG at pH 2.5.

Ab initio calculations and reduced density gradient analyses of the structure and energetics of hydrated calcium fluoride and calcium carbonate

We studied microhydrated calcium fluoride, calcium carbonate and their ions at the MP2/6-311++G** level of theory and found that water–water non-covalent interactions destabilize the solvation shell, and are compensated by cooperative hydrogen bonds.

Freeman C, Harding JH. Interaction of stable aggregates drives the precipitation of calcium phosphate in supersaturated solutions

Using molecular dynamics simulations, we show for the first time that calcium phosphate nanoparticles of eight formula units are thermodynamically stable and could be key in the nucleation of amorphous calcium phosphate.