Symmetry Breaking and Hydration Structure of Carbonate and Nitrate in Aqueous Solutions: A Study by Ab Initio Quantum Mechanical Charge Field Molecular Dynamics

Unveiling the Structural and Dynamic Characteristics of Concentrated LiNO3 Aqueous Solutions through Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted a significant amount of attention for their potential applications in lithium-ion batteries. Nevertheless, a comprehensive understanding of the Li+ solvation structure and its migration within electrolyte solutions remains elusive. This study employs linear vibrational spectroscopy, ultrafast infrared spectroscopy, and molecular dynamics (MD) simulations to elucidate the structural dynamics in LiNO3 solutions by using intrinsic and extrinsic vibrational probes. The N-O stretching vibrations of NO3- exhibit a distinct spectral splitting, attributed to its asymmetric interaction with the surrounding solvation structure. Analysis of the vibrational relaxation dynamics of intrinsic and extrinsic probes, in combination with MD simulations, reveals cage-like networks formed through electrostatic interactions between Li+ and NO3-. This microscopic heterogeneity is reflected in the intertwined arrangement of ions and water molecules. Furthermore, both vehicular transport and structural diffusion assisted by solvent rearrangement for Li+ were analyzed, which are closely linked with the bulk concentration.

Uncovering the binding nature of thiocyanate in contact ion pairs with lithium ions

Ion pair formation is a fundamental molecular process that occurs in a wide variety of systems, including electrolytes, biological systems, and materials. In solution, the thiocyanate (SCN-) anion interacts with cations to form contact ion pairs (CIPs). Due to its ambidentate nature, thiocyanate can bind through either its sulfur or nitrogen atoms, depending on the solvent. This study focuses on the binding nature of thiocyanate with lithium ions as a function of the solvents using FTIR, 2D infrared spectroscopy (2DIR) spectroscopies, and theoretical calculations. The study reveals that the SCN- binding mode (S or N end) in CIPs can be identified through 2DIR spectroscopy but not by linear IR spectroscopy. Linear IR spectroscopy shows that the CN stretch frequencies are too close to one another to separate N- and S-bound CIPs. Moreover, the IR spectrum shows that the S-C stretch presents different frequencies for the salt in different solvents, but it is related to the anion speciation rather than to its binding mode. A similar trend is observed for the anion bend. 2DIR spectra show different dynamics for N-bound and S-bound thiocyanate. In particular, the frequency-frequency correlation function (FFCF) dynamics extracted from the 2DIR spectra have a single picosecond exponential decay for N-bound thiocyanate and a biexponential decay for S-bound thiocyanate, consistent with the binding mode of the anion. Finally, it is also observed that the binding mode also affects the line shape parameters, probably due to the different molecular mechanisms of the FFCF for N- and S-bound CIPs.

Asymmetry of the Electrostatic Environment as the Origin of the Symmetry Breaking Effect of the Nitrate Ion in Aqueous Solution

Elucidating the mechanism of how vibrational modes are affected by intermolecular interactions is important for a better understanding of the nature of the former as probes of the latter. Here, such an analysis is carried out for the N-O stretching modes of the nitrate ion interacting with water, with an emphasis on the symmetry breaking effect. On the basis of theoretical calculations on the structural, vibrational, and electrostatic properties of molecular clusters and spectral simulations for an aqueous solution, a transparent view is demonstrated on the mechanism that modulations of spatially local electrostatic environment give rise to structural and spectroscopic symmetry breaking effect. The electrostatic interaction model constructed here is a seven-parameter model; the use of a single electrostatic parameter, such as the electric field on a single atomic site, is found to be insufficient for quantitative evaluation. It is also shown that the frequency modulations of the N-O stretching modes in aqueous solution occur on a time scale much shorter than 0.1 ps.

On the existence of CO32- microsolvated clusters: a theoretical study

Microsolvated clusters of multiply charged anions play a crucial role in atmospheric chemistry and some of them were previously registered experimentally. At the same time, there are no experimental observations of [CO₃·(H₂O)_n]^2-. The reasons for this may be related to the thermodynamical or kinetical instability of microsolvated CO₃^2- toward autoionization or autoprotonation processes. In this study we theoretically investigate the potential stability of the [CO₃·(H₂O)_n]^2- microsolvated clusters from both perspectives - thermodynamic and kinetic - and we claim they are stable toward autoionization and kinetically semi-stable toward autoprotonation. In addition, the behaviour of CO₃^2- anions in bulk water solvent was analysed to highlight important precautions for synthetic purposes.

Molecular modeling and simulation of aqueous solutions of alkali nitrates

A set of molecular models for the alkali nitrates (LiNO3, NaNO3, KNO3, RbNO3, and CsNO3) in aqueous solutions is presented and used for predicting the thermophysical properties of these solutions with molecular dynamics simulations. The set of models is obtained from a combination of a model for the nitrate anion from the literature with a set of models for the alkali cations developed in previous works of our group. The water model is SPC/E and the Lorentz–Berthelot combining rules are used for describing the unlike interactions. This combination is shown to yield fair predictions of thermophysical and structural properties of the studied aqueous solutions, namely the density, the water activity and the mean ionic activity coefficient, the self-diffusion coefficients of the ions, and radial distribution functions, which were studied at 298 K and 1 bar; except for the density of the solutions of all five nitrates and the activity properties of solutions of NaNO3, which were also studied at 333 K. For calculating the water the activity and the mean ionic activity coefficient, the OPAS ( osmotic pressure for the activity of selvents) method was applied. The new models extend an ion model family for the alkali halides developed in previous works of our group in a consistent way.

Solubilities in aqueous nitrate solutions that appear to reverse the law of mass action

Non-ideal aqueous electrolyte solutions have been studied since the start of the application of thermodynamics to chemistry in the late 19th century. The present study examines some of the most extreme non-ideal behavior ever observed: solubilities of alkali and NH₄⁺ nitrate salts in water that appear to behave the opposite of how the Law of Mass Action would predict. A literature review discovered that the solubilities of NH₄NO₃ and many alkali nitrate salts increases when another nitrate-bearing electrolyte is added to solution. These occurrences were in concentrated solutions with insufficient water to provide all ions their preferred hydration number without sharing waters between ions. This water deficit results in the formation of contact ion-pairs as well as larger ion-clusters. These ion-clusters may be favored when there is more than one type of monovalent cation present.

The unexpected effect of aqueous ion pairs on the forbidden n →π* transition in nitrate

Aqueous nitrate is ubiquitous in the environment, found for example in stratospheric clouds, tropospheric particulate matter, rain and snow, fertilized fields, rivers and the ocean. Its photolysis is initiated by absorption into the strongly forbidden n →π* transition. Photolysis reactivates deposited nitrate, releasing nitrogen oxides, and UV light is commonly used to break down nitrate pollution. The transition is doubly forbidden unless its symmetry is broken, giving a powerful means of probing the interactions of nitrate with its environment and of using experiment to validate the results of theory. In this study we demonstrate the remarkably different effects of the addition of a series of mono- and di-valent metal chlorides on the nitrate UV transition. While they all shift the transition to shorter wavelengths, the shift changes significantly from one to another. For the monovalent series Li⁺, Na⁺, K⁺, the blue shift decreases down the column being strongest for Li⁺ and weakest for K⁺. For the divalent series Mg²⁺, Ca²⁺, Ba²⁺, the opposite effect is observed with the energy shift of Ba²⁺ being an order of magnitude larger than for Mg²⁺. The absorption intensity also changes; the addition of Na⁺ and K⁺ decrease intensity whereas Li⁺ increases intensity. For the divalent cations an increase is seen for all three members of the series Mg²⁺, Ca²⁺ and Ba²⁺. Paradoxically, the effect of addition of CaCl₂ to the solution is to decrease the environmental photolysis rate of nitrate; despite the increase in intensity, Ca²⁺ blue shifts the peak position above the tropospheric photolysis threshold around 300 nm. Using computational chemistry we conclude that the effects are due to the microscopic interactions of the nitrate anion and not continuum effects. Two microscopic mechanisms are investigated in detail, the formation of a nitrate monohydrate cluster and a contact ion pair. The contact ion pair shows the potential for significant impact on the energy and intensity of the transition.

Ion-ion interactions enhance aluminum solubility in alkaline suspensions of nano-gibbsite (α-Al(OH)3) with sodium nitrite/nitrate

Despite widespread industrial importance, predicting metal solubilities in highly concentrated, multicomponent aqueous solutions is difficult due to poorly understood ion-ion and ion-solvent interactions. Aluminum hydroxide solid phase solubility in concentrated sodium hydroxide (NaOH) solutions is one such case, with major implications for ore refining, as well as processing of radioactive waste stored at U.S. Department of Energy legacy sites, such as the Hanford Site, Washington State. The solubility of gibbsite (α-Al(OH)₃) is often not well predicted because other ions affect the activity of hydroxide (OH^-) and aluminate (Al(OH)₄^-) anions. In the present study, we systematically examined the influence of key anions, nitrite (NO₂^-) and nitrate (NO₃^-), as sodium salts on the solubility of α-Al(OH)₃ in NaOH solutions taking care to establish equilibrium from both under- and oversaturation. Rapid equilibration was enabled by use of a highly pure and crystalline synthetic nano-gibbsite of well-defined particle size and shape. Measured dissolved aluminum concentrations were compared with those predicted by an α-Al(OH)₃ solubility model derived for simple Al(OH)₄^-/OH^- systems. Specific anion effects were expressed as an enhancement factor (Al_enhc) conveying the excess of dissolved aluminum. At 45 °C, NaNO₂ and NaNO₃-containing systems exhibited Al_enhc values of 2.70 and 1.88, respectively, indicating significant enhancement. The solutions were examined by Raman and high-field ²⁷Al NMR spectroscopy, indicating specific interactions including Al(OH)₄^--Na⁺ contact ion pairing and Al(OH)₄^--NO₂^-/NO₃^- ion-ion interactions. Dynamic evolution of the α-Al(OH)₃ particles including growth and agglomeration was observed revealing the importance of dissolution/reprecipitation in establishing equilibrium. These studies indicate that incomplete ion hydration, as a result of the low water activity in these concentrated electrolytes, results in: (i) enhanced reactivity of the hydroxide ion with respect to α-Al(OH)₃; (ii) increased concentrations of Al(OH)₄^- in solution; and (iii) stronger ion-ion interactions that act to stabilize the supersaturated solutions. This information on the mechanisms by which α-Al(OH)₃ becomes supersaturated is essential for more energy-efficient aluminum processing technologies, including the treatment of millions of gallons of Al(OH)₄^--rich high-level radioactive waste.

Salt Solubilities in Aqueous Solutions of NaNO3, NaNO2, NaCl, and NaOH: A Hofmeister-like Series for Understanding Alkaline Nuclear Waste

Nonelectrolyte solubility in electrolyte solutions follow the Hofmeister series, but the applicability of the series to salt solubility has been less appreciated. This study, using solubility data for thirteen sodium-bearing salts, shows that salts are consistently salted out by electrolytes important to alkaline nuclear waste in the order NaOH > NaCl > NaNO₂ > NaNO₃ at 298.15 K, which is the same order as the Hofmeister series. Graphical presentation allowed for easy separation of the common ion effect (caused by the addition of Na⁺) from the salting-out effect (caused by the presence of anions) because there is a large difference between the solubility of a given salt in different background electrolytes at a common Na⁺ molality. The trend persists even in very high electrolyte concentrations where essentially all of the water molecules must be in the coordination sphere of an ion, which means that the effect of electrolytes on "bulk water" is not the cause of the trend. These specific interactions more likely result from the sharing of water molecules between ions, augmented by differences in ion-pairing of the electrolytes. The Hofmeister series has practical application to the management of alkaline high-level radioactive waste created at nuclear fuel reprocessing facilities, where a predictive understanding of salt solubility is essential for blending wastes of disparate compositions prior to treatment.

Decoding Oxyanion Aqueous Solvation Structure: A Potassium Nitrate Example at Saturation

The ability to probe the structure of a salt solution at the atomic scale is fundamentally important for our understanding of many chemical reactions and their mechanisms. The capability of neutron diffraction to "see" hydrogen (or deuterium) and other light isotopes is exceptional for resolving the structural complexity around the dissolved solutes in aqueous electrolytes. We have made measurements using oxygen isotopes on aqueous nitrate to reveal a small hydrogen-bonded water coordination number (3.9 ± 1.2) around a nitrate oxyanion. This is compared to estimates made using the existing method of nitrogen isotope substitution and those of computational simulations (>5-6 water molecules). The low water coordination number, combined with a comparison to classical molecular dynamics simulations, suggests that ion-pair formation is significant. This insight demonstrates the utility of experimental diffraction data for benchmarking atomistic computer simulations, enabling the development of more accurate intermolecular potentials.

Structure and dynamics of propylammonium nitrate-acetonitrile mixtures: An intricate multi-scale system probed with experimental and theoretical techniques

In this article, we report the study of structural and dynamical properties for a series of acetonitrile/propylammonium nitrate mixtures as a function of their composition. These systems display an unusual increase in intensity in their X-ray diffraction patterns in the low-q regime, and their ¹H-NMR diffusion-ordered NMR spectroscopy (DOSY) spectra display unusual diffusivities. However, the magnitude of both phenomena for mixtures of propylammonium nitrate is smaller than those observed for ethylammonium nitrate mixtures with the same cosolvent, suggesting that the cation alkyl tail plays an important role in these observations. The experimental X-ray scattering data are compared with the results of molecular dynamics simulations, including both ab initio studies used to interpret short-range interactions and classical simulations to describe longer range interactions. The higher level calculations highlight the presence of a strong hydrogen bond network within the ionic liquid, only slightly perturbed even at high acetonitrile concentration. These strong interactions lead to the symmetry breaking of the NO₃^- vibrations, with a splitting of about 88 cm^-1 in the ν₃ antisymmetric stretch. The classical force field simulations use a greater number of ion pairs, but are not capable of fully describing the longest range interactions, although they do successfully account for the observed concentration trend, and the analysis of the models confirms the nano-inhomogeneity of these kinds of samples.

Carbonate and carbonate anion radicals in aqueous solutions exist as CO3(H2O)62− and CO3(H2O)6˙− respectively: the crucial role of the inner hydration sphere of anions in explaining their properties

An effort to reproduce the physical properties of CO₃²⁻ and CO₃˙⁻ in water proves that one has to include an inner hydration sphere of six water molecules for both anions.

Properties of aqueous nitrate and nitrite from x-ray absorption spectroscopy

Nitrate and nitrite ions are of considerable interest, both for their widespread use in commercial and research contexts and because of their central role in the global nitrogen cycle. The chemistry of atmospheric aerosols, wherein nitrate is abundant, has been found to depend on the interfacial behavior of ionic species. The interfacial behavior of ions is determined largely by their hydration properties; consequently, the study of the hydration and interfacial behavior of nitrate and nitrite comprises a significant field of study. In this work, we describe the study of aqueous solutions of sodium nitrate and nitrite via X-ray absorption spectroscopy (XAS), interpreted in light of first-principles density functional theory electronic structure calculations. Experimental and calculated spectra of the nitrogen K-edge XA spectra of bulk solutions exhibit a large 3.7 eV shift between the XA spectra of nitrate and nitrite resulting from greater stabilization of the nitrogen 1s energy level in nitrate. A similar shift is not observed in the oxygen K-edge XA spectra of NO3 (-) and NO2 (-). The hydration properties of nitrate and nitrite are found to be similar, with both anions exhibiting a similar propensity towards ion pairing.

Effects of ions on water structure: a low-field ¹H T₁ NMR relaxometry approach

Aqueous salt solutions play an important role in nature because of their effects on environmental biogeochemical processes and on structural properties of biomolecules. Upon dissolution, salts split in ions that are solvated. Water in hydration shells is subjected to molecular motions that can be monitored by (1)H T1 NMR relaxometry. This technique allowed the evaluation of the nature of the interactions between water and ions via variable temperature experiments. Examination of relaxometry properties of aqueous solutions at variable salt concentrations allowed acknowledgement of the role played by ions in either structuring or destructuring water aggregates. A mathematical model has been applied on six environmentally relevant salts: NaCl, KCl, CaCl2, CaCO3, NaNO3, and NH4NO3. It was linear only for the concentration dependence of KCl-R1. This model accorded with the one reported in literature where it has been considered valid only for diluted solutions. However, in the present study, the range of linearity for KCl was extended up to the saturation point. The model was modified for NaCl, CaCl2, and CaCO3 by using it as an exponential form in order to account for the nonlinearity of the R1-versus-concentration curves. Nonlinearity was explained by the nonnegligible ion-ion interactions occurring as concentration was increased. Finally, further modification was needed to account for the asymmetric distribution of water around nitrate (in NaNO3 and NH4NO3) and ammonium (in NH4NO3). This study is preliminary to the comprehension of the diffusion mechanisms of ions in water solutions at the equilibrium condition with solid surfaces such as soils and biochar-amended soils.

Hydration properties determining the reactivity of nitrite in aqueous solution

The knowledge of the hydration properties of the nitrite ion is key to understanding its reaction mechanism controlled by solvent effects. Here, ab initio quantum mechanical charge field molecular dynamics was performed to obtain the structural and dynamical properties of the hydration shell in an aqueous solution of nitrite ions, elucidated by data analysis using a molecular approach and an extended quantitative analysis of all superimposed trajectories with three-dimensional alignment (density map). The pattern of the power spectra corresponded to the experimental data, indicating the suitability of the Hartree-Fock method coupled with double-ζ plus polarization and diffuse functional basis sets to study this system. The density maps revealed the structure of the hydration shell, that presented a higher density in the N-O bond direction than in the axis vertical to the molecular plane, whereas the atomic and molecular radial distribution functions provided vague information. The number of actual contacts indicated 4.6 water molecules interacting with a nitrite ion, and 1.5 extra water molecules located in the molecular hydration shell, forming a H-bonding network with the bulk water. The mean residence times for the water ligands designated the strength of the hydration spheres for the oxygen sites, whilst the results for the nitrogen sites over-estimated the number of water molecules from other sites and indicated a weak structure. These results show the influence of the water molecules surrounding the nitrite ion creating an anisotropic hydration shell, suggesting that the reactive sites are situated above and below the molecular plane with a lower water density.

Water in the formation of biogenic minerals: peeling away the hydration layers

Minerals of biogenic origin form and crystallize from aqueous environments at ambient temperatures and pressures. The in vivo environment either intracellular or intercellular, contains many components that modulate both the activity of the ions which associate to form the mineral, as well as the activity and structure of the crowded water. Most of the studies about the mechanism of mineralization, that is, the detailed pathways by which the mineral ions proceed from solution to crystal state, have been carried out in relatively dilute solutions and clean solutions. These studies have considered both thermodynamic and kinetic controls. Most have not considered the water itself. Is the water a passive bystander, or is it intimately a participant in the mineral ion densification reaction? A wide range of experiments show that the mineralization pathways proceed through a series of densification stages with intermediates, such as a "dense liquid" phase and the prenucleation clusters that form within it. This is in contrast to the idea of a single step phase transition, but consistent with the Gibbs concept of discontinuous phase transitions from supersaturated mother liquor to crystal. Further changes in the water structure at every surface and interface during densification guides the free energy trajectory leading to the crystalline state. In vertebrates, mineralization takes place in a hydrated collagen matrix, thus water must be considered as a direct participant. Although different in detail, the crystallization of calcium phosphates, as apatite, and calcium carbonates, as calcite, are mechanistically identical from the viewpoint of water.