Coupled jump rotational dynamics in aqueous nitrate solutions

Tracking nitrite's deviation from Stokes-Einstein predictions with pulsed field gradient 15N NMR spectroscopy

Predicting the behavior of oxyanions in radioactive waste stored at the Department of Energy legacy nuclear sites requires the development of novel analytical methods. This work demonstrates 15N pulsed field gradient nuclear magnetic resonance spectroscopy to quantify the diffusivity of nitrite. Experimental results, supported by molecular dynamics simulations, indicate that the diffusivity of free hydrated nitrite exceeds that of free hydrated sodium despite the greater hydrodynamic radius of nitrite. Investigations are underway to understand how the compositional and dynamical heterogeneities of the ion networks at high concentrations affect rheological and transport properties.

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.

Ab Initio Molecular Dynamics Study of Aqueous Solutions of Magnesium and Calcium Nitrates: Hydration Shell Structure, Dynamics and Vibrational Echo Spectroscopy

Ab initio molecular dynamics simulations are performed to study the hydration shell structure, dynamics, and vibrational echo spectroscopy of aqueous Mg(NO₃)₂ and Ca(NO₃)₂ solutions. The hydration shell structure is probed through calculations of various ion-ion and ion-water radial and spatial distribution functions. On the dynamical side, calculations have been made for the hydrogen bond dynamics of hydration shells and also residence dynamics and lifetimes of water in different solvation environments. Subsequently, we looked at the dynamics of frequency fluctuations of OD modes of heavy water in different hydration environments. Specifically, the temporal decay of spectral observables of two-dimensional infrared (2DIR) spectroscopy, three pulse echo peak shift (3PEPS) measurements and also of time correlations of frequency fluctuations are calculated to investigate the dynamics of vibrational spectral diffusion of water in different hydration environments in these solutions. The OD stretch frequencies of water molecules in the vicinity of both divalent cations are found to be red-shifted and also fluctuating at a slower rate than other water molecules present in the solutions. The Mg²⁺ ions are found to be strongly hydrated which can be linked to their lower tendency to form contact ion-pairs and essentially no water exchange between the cationic hydration shells and bulk during the time scale of the current simulations. The stronger hydration of Mg²⁺ ions make their hydration shells structurally and dynamically more rigid and make the dynamics of hydrogen bonds and vibrational spectral diffusion, as revealed through spectral observables of 2DIR and 3PEPS slower than that for the Ca²⁺ ions. The structural and spectral dynamics of water molecules outside the cationic solvation shells in the Mg(NO₃)₂ solution are also found to be relatively slower than that of the Ca(NO₃)₂ solution and pure water which show the effects of stronger electric fields of Mg²⁺ ions extending beyond their first hydration shells. Also, water molecules in the hydration shells of the NO₃^- ions are found to relax at a slower rate in the Mg(NO₃)₂ solution which manifests the effect countercations have on anionic hydration shells for divalent metal nitrate solutions.

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.

Structural and dynamic properties of some aqueous salt solutions

Aqueous salt solutions are utilized and encountered in wide-ranging technological applications and natural settings. Towards improved understanding of the effect of salts on the dynamic properties of such systems, dilute aqueous salt solutions (up to 1 molar concentration) are investigated here, via experiments and molecular simulations. Four salts are considered: sodium chloride, for which published results are readily available for comparison, ammonium acetate, barium acetate and barium nitrate, for which published data are scarce. In the present work, molecular dynamics (MD) simulations are conducted to quantify viscosity and water self-diffusion coefficients, together with rheometry and Pulsed Field Gradient Spin Echo (PFGSE)-NMR experiments for validation. Simulation predictions are consistent with experimental observations in terms of trend and magnitude of salt-specific effects. Combining insights from the approaches considered, an interpretation of the results is proposed whereby the capacity of salts to influence bulk dynamics arises from their molecular interfacial area and strength of interaction with first hydration-shell water molecules. For the concentration range investigated, the interpretation could be useful in formulating aqueous systems for applications including the manufacturing of advanced catalysts.

Microscopic structural features of water in aqueous-reline mixtures of varying compositions

Reline, a mixture of urea and choline chloride in a 2 : 1 molar ratio, is one of the most frequently used deep eutectic solvents. Pure reline and its aqueous solution have large scale industrial use. Owing to the presence of active hydrogen bond formation sites, urea and choline cations can disrupt the hydrogen-bonded network in water. However, a quantitative understanding of the microscopic structural features of water in the presence of reline is still lacking. We carry out extensive all-atom molecular dynamics simulations to elucidate the effect of the gradual addition of co-solvents on the microscopic arrangements of water molecules. We consider four aqueous solutions of reline, between 26.3 and 91.4 wt%. A disruption of the local hydrogen-bonded structure in water is observed upon inclusion of urea and choline chloride. The extent of deviation of the water structure from tetrahedrality is quantified using the tetrahedral order parameter (q_tet). Our analyses show a monotonic increase in the structural disorder as the co-solvents are added. Increase in the q_tet values are observed when highly electro-negative hetero-atoms like nitrogen, oxygen of urea and choline cations are counted as partners of the central water molecules. Further insights are drawn from the characterization of the hydrogen-bonded network in water and we observe the gradual rupturing of water-water hydrogen bonds and their subsequent replacement by the water-urea hydrogen bonds. A negligible contribution from the hydrogen bonds between water and bulky choline cations has also been found. Considering all the constituents as the hydrogen bond partners we calculate the possibility of a successful hydrogen bond formation with a central water molecule. This gives a clear picture of the underlying mechanism of water replacement by urea.

Role of local order in anomalous ion diffusion: Interrogation through tetrahedral entropy of aqueous solvation shells

Small rigid ions perturb the water structure around them significantly. At constant viscosity, alkali cations (Li⁺, Na⁺, and so on) exhibit an anomalous non-monotonic dependence of diffusivity on ion-size, in stark violation of the Stokes-Einstein expression. Although this is a well-known problem, we find that an entropic view of the problem can be developed, which provides valuable insight. The local entropy experienced by the solute ion is relevant here, which leads to the connection with local viscosity, discussed earlier by many. Due to the strong interactions with ions, the translational and rotational entropy of solvation water decreases sharply; however, an opposite effect comes from the disruption of the tetrahedral network structure of water near the charges. We compute the tetrahedral order of water molecules (q_tet) around the ion and suitably defined tetrahedral entropy [S(q_tet)] that is a contribution to the excess entropy of the system. Our results reveal that although the structural properties of the second shell become nearly identical to the bulk, S(q_tet) of the second shell is found to play an important role in giving rise to the non-monotonic ion-size dependence. The detailed study of the static and dynamic fluctuations in q_tet and the number of hydration water molecules provides interesting insights into correlation between the structure and dynamics; the smallest static fluctuation of q_tet for the first hydration shell water molecules of Li⁺ is indicative of the iceberg picture. The study of fluctuation properties of q_tet and the coordination number also reveals the role of the second hydration layer and could explain the anomalous behavior of the Rb⁺ ion.

Two-Dimensional Infrared Spectroscopy of Aqueous Solutions of Metal Nitrates: Slowdown of Spectral Diffusion in the Presence of Divalent Cations

The hydrogen-bonded network of water can be affected both structurally and dynamically by the presence of ions. In the present study, we have considered three aqueous solutions of metal nitrates to investigate the effects of divalent cations (Mg²⁺ and Ca²⁺), compared to that of monovalent Na⁺ ions, on hydrogen-bond fluctuations and vibrational spectral diffusion through calculations of linear and two-dimensional infrared spectra of these solutions at room temperature. We have employed the methods of molecular dynamics simulations using effective polarizable models of ions combined with quantum mechanical calculations of transition variables and statistical mechanical calculations of spectral response functions of vibrational spectroscopy. Divalent cations are found to have much stronger and longer-ranged effects on the structure and dynamics of the hydrogen-bonded network than that induced by the monovalent sodium ions. The blue shifts in the calculated linear spectra are found to follow the Hofmeister trend for the cations. The 2D-IR spectral lineshape and intensity corresponding to three-pulse echo peak shift (3PEPS) experiments are calculated. The timescales of these nonlinear spectral responses and also frequency-time correlations show significant slowing down of spectral diffusion for solutions containing divalent Mg²⁺ and Ca²⁺ ions compared to the corresponding dynamics of the solution containing Na⁺ ions. Unlike NaNO₃ solution, the relaxation of frequency and dipole orientational fluctuations of anion-bound water in Mg(NO₃)₂ and Ca(NO₃)₂ solutions are found to be somewhat slower than bulk water, which can be attributed to the presence of divalent cations whose effects go beyond their first solvation shells. This is also seen in the dynamics of bulk water in these solutions which is found to be notably slower for the solutions containing divalent cations than that in the NaNO₃ solution. Unlike Mg²⁺ and Ca²⁺ ions, no specific cationic effect is observed for the Na⁺ ions.

Solvation Shell of the Nitrite Ion in Water: An Ab Initio Molecular Dynamics Study

We performed ab initio molecular dynamics simulation of a nitrite ion in water to investigate the structural and dynamical properties of its hydration shell. The nitrite ion is found to exhibit strong asymmetry toward hydrogen bonding due to its two different types of hydrogen bond acceptor sites. This difference is better captured through further partitioning of the hydration shell into its proximal and distal regions. The frequency shifts of the stretch modes of hydration shell water reveal that the nitrogen site forms a stronger hydrogen bond than its oxygen sites with the latter forming hydrogen bonds, which are similar in strength to that between a pair of water molecules. The escape dynamics of water from the hydration shell is found to be rather slow, which seems to classify the nitrite ion as a structure-maker. However, the dynamics of orientational and hydrogen bond relaxation reveal a faster mobility of water molecules in the hydration shell than bulk water in spite of strong ion-water interactions. It is found that the nitrite ion can hold water molecules in its solvation shell and still make them rotate fast in its vicinity through switching of their hydrogen bonds between its nitrogen and oxygen acceptor sites. The dipole moment of the solute in water is also calculated in the present study.

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.

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.

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.

Rotational dynamics of polyatomic ions in aqueous solutions: From continuum model to mode-coupling theory, aided by computer simulations

Due to the presence of the rotational mode and the distributed surface charges, the dynamical behavior of polyatomic ions in water differs considerably from those of the monatomic ions. However, their fascinating dynamical properties have drawn scant attention. We carry out theoretical and computational studies of a series of well-known polyatomic ions, namely, sulfate, nitrate, and acetate ions. All three ions exhibit different rotational diffusivity, with that of the nitrate ion being considerably larger than the other two. They all defy the hydrodynamic laws of size dependence. Study of the local structure around the ions provides valuable insight into the origin of these differences. We carry out a detailed study of the rotational diffusion of these ions by extensive computer simulation and by using the theoretical approaches of the dielectric friction developed by Fatuzzo-Mason (FM) and Nee-Zwanzig (NZ), and subsequently generalized by Alavi and Waldeck. A critical element of the FM-NZ theory is the decomposition of the total rotational friction, ζ_Rot, into Stokes and dielectric parts. The study shows a dominant role of dielectric friction in the sense that if the ions are made neutral, the nature of diffusion changes and the values become much larger. Our analyses further reveal that the decomposition of total friction into the Stokes and dielectric friction breaks down for sulfate ions but remains semi-quantitatively valid for nitrate and acetate ions. We discuss the relationship between translational and rotational dielectric friction on rigid spherical ions. We develop a self-consistent mode-coupling theory (SC-MCT) formalism that could provide a unified view of rotational friction of polyatomic ions in polar medium. Our SC-MCT shows that the breakdown can be attributed to the change in the microscopic structural features. The mode-coupling theory helps in elucidating the role of coupling between translational and rotational motion of these ions. In fact, these two motions self-consistently determine the value of each other. The reference interaction site model-based MCT suggests an interesting relation between the torque-torque and the force-force time correlation function with the proportionality constant being determined by the geometry and the charge distribution of the polyatomic molecule. We point out several parallelisms between the theories of translational and rotation friction calculations of ions in polar liquids.