A Combined QM/MM Molecular Dynamics Simulations Study of Nitrate Anion (NO3-) in Aqueous Solution

The isosbestic point in the Raman spectra of the hydration shell

Isosbestic point is often observed in a series of spectra, but their interpretation is still controversial, such as whether the continuum model can produce an isosbestic point. In order to answer this question, the Raman spectra of hydration shell with continuous distribution structure in different ionic aqueous solutions were separated by Raman ratio spectra, and an isosbestic point was successfully observed. Our experimental results show that the continuum model can indeed produce the isosbestic point. In order to deepen the understanding of the isosbestic point, we calculate the first moment of the Raman spectra and conduct molecular dynamics (MD) simulations. Both experimental and theoretical findings indicate that elevated temperatures lead to increased disorder among water molecules within the hydration shell.

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.

Liquid Structure and Hydrogen Bonding in Aqueous Hydroxylammonium Nitrate

Hydroxylammonium nitrate (HAN) has emerged as a promising component in ionic liquid-based spacecraft propellants. However, the physicochemical and structural properties of aqueous HAN have been largely overlooked. The purpose of this study is to investigate the hydrogen bonding in aqueous HAN and understand its implications on these properties and the proton transfer mechanism as a function of concentration. Classical polarizable molecular dynamics simulations have been employed with the APPLE&P force field to analyze the geometry of individual hydrogen bonds and the overall hydrogen-bonding network in various concentrations of aqueous HAN. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) indicate the structural arrangement of the species and their hydrogen bonds. Projections of water density and the orientation of its electric dipole moment near the ions provide insight into the hydrogen-bonding network. The incorporation of water into the hydrogen-bonding network at high ion concentrations occurs via interstitial accommodation around the ions immediately outside the first solvation shell. While ion pairs are observed at all concentrations considered, the frequency of Ha···On hydrogen bonds increases substantially with the ion concentration. The findings contribute to a better fundamental understanding of HAN and the precursors of reactivity, crucial to the development of "green" spacecraft propellants.

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.

Unveiling the structure of aqueous magnesium nitrate solutions by combining X-ray diffraction and theoretical calculations

The structure of aqueous magnesium nitrate solution is gaining significant interest among researchers, especially whether contact ion pairs exist in concentrated solutions. Here, combining X-ray diffraction experiments, quantum chemical calculations and ab initio molecular dynamics simulations, we report that the [Mg(NO₃)₂] molecular structure in solution from the coexistence of a free [Mg(H₂O)₆]²⁺ octahedral supramolecular structure with a free [NO₃(H₂O)_n]^- (n = 11-13) supramolecular structure to an [Mg²⁺(H₂O)_n(NO₃^-)_m] (n = 3, 4, 5; m = 3, 2, 1) associated structure with increasing concentration. Interestingly, two hydration modes of NO₃^--the nearest neighbor hydration with a hydration distance less than 3.9 Å and the next nearest neighbor hydration with hydration distance ranging from 3.9 to 4.3 Å-were distinguished. With an increase in the solution concentration, the hydrated NO₃^- ions lost outer layer water molecules, and the hexagonal octahedral hydration structure of [Mg(H₂O)₆²⁺] was destroyed, resulting in direct contact between Mg²⁺ and NO₃^- ions in a monodentate way. As the concentration of the solution further increased, NO₃^- ions replaced water molecules in the hydration layer of Mg²⁺ to form three-ion clusters and even more complex chains or linear ion clusters.

Structure of Aqueous KNO3 Solutions by Wide-Angle X-ray Scattering and Density Functional Theory

The structure of aqueous KNO₃ solutions was studied by wide-angle X-ray scattering (WAXS) and density functional theory (DFT). The interference functions were subjected to empirical potential structure refinement (EPSR) modeling to extract the detailed hydration structure information. In aqueous KNO₃ solutions with a concentration from 0.50 to 2.72 mol·dm^-3, water molecules coordinate K⁺ with a mean coordination number (CN) from 6.6 ± 0.9 to 5.8 ± 1.2, respectively, and a K-O_W(H₂O) distance of 2.82 Å. To further describe the hydration properties of K⁺, a hydration factor (f_h) was defined based on the orientational angle between the water O-H vector and the ion-oxygen vector. The f_h value obtained for K⁺ is 0.792, 0.787, 0.766, and 0.765, and the corresponding average hydration numbers (HN) are 5.2, 5.1, 4.8, and 4.5. The reduced HN is compensated by the direct binding of oxygen atoms of NO₃^- with the average CN from 0.3 ± 0.7 to 2.6 ± 1.7, respectively, and the K-O_N(NO₃^-) distance of 2.82 Å. The average number of water molecules around NO₃^- slightly reduces from 10.5 ± 1.5 to 9.6 ± 1.7 with r_N-OW = 3.63 Å. K⁺-NO₃^- ion association is characterized by K-O_N and K-N pair correlation functions (PCFs). A K-N peak is observed at 3.81 Å as the main peak with a shoulder at 3.42 Å in g_K-N(r). This finding indicates that two occupancy sites are available for K⁺, i.e., the edge-shared bidentate (BCIP) and the corner-shared monodentate (MCIP) contact ion pairs. The structure and stability of the KNO₃(H₂O)₁₀ cluster were investigated by a DFT method and cross-checked with the results from WAXS.

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.

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.

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.

Extraction of Gd3+ and UO22+ Ions Using Polystyrene Grafted Dibenzo Crown Ether (DB18C6) with Octanol and Nitrobenzene: A Molecular Dynamics Study

Atomistic molecular dynamics (MD) simulations are performed in order to derive thermodynamic properties important to understand the extraction of gadolinium (Gd³⁺) and uranium dioxide (UO₂) with dibenzo crown ether (DBCE) in nitrobenzene (NB) and octanol (OCT) solvents. The effect of polystyrene graft length, on DBCE, on the binding behavior of Gd³⁺ and UO₂²⁺ is investigated for the first time. Our simulation results demonstrate that the binding of Gd³⁺ and UO₂²⁺ onto the oxygens of crown ethers is favorable for polystyrene grafted crown ether in the organic solvents OCT and NB. The metal ion binding free energy (ΔG_Binding) in different solvent environments is calculated using the thermodynamic integration (TI) method. ΔG_Binding becomes more favorable in both solvents, NB and OCT, with an increase in the polystyrene monomer length. The metal ion transferability from an aqueous phase to an organic phase is estimated by calculating transfer free-energy calculations (ΔG_Transfer). ΔG_Transfer is significantly favorable for both Gd³⁺ and UO₂²⁺ for the transfer from the aqueous phase to the organic phase (i.e., NB and OCT) via ion-complexation to DBCE with an increase in polystyrene length. The partition coefficient (log P) values for Gd³⁺ and UO₂²⁺ show a 5-fold increase in separation capacity with polystyrene grafted DBCE. We corroborate the observed behavior by further analyzing the structural and dynamical properties of the ions in different phases.

Ab initio investigation of the aqueous solvation of the nitrate ion

The turning point between surface and interior solvation of NO₃⁻ is predicted to lie around a cluster size of (H₂O)₆₄.

DFT study on X⁻·(H₂O)(n=1-10) (X=OH, NO₂, NO₃, CO₃) anionic water cluster

A theoretical study to understand the interaction between anion and the water molecules through the hydration (X(-)·(H2O)n (X=OH, NO2, NO3, CO3), where n=1-10), using the density functional theory method with B3LYP functional and 6-311++G(d,p) basis set has been carried out systematically. In these hydrated clusters we notice three different cases of bond arrangements, namely, symmetrical double hydrogen bond, single hydrogen bond and inter-water hydrogen bond. All the complexes are dominated by the O-H⋯O hydrogen bond, in which the anion act as a proton acceptor, while the water molecule act as a proton donor. A linear correlation is obtained between the solvent stabilization energy and the size (n) of the hydrated cluster for all the anions. The weighted average interaction energy values, shows that the water molecules strongly bind with the OH(-) anion. Besides, the solvation of the OH(-) anion requires less number of water molecules when compared with the other anions. Energy decomposition analysis (EDA) shows the strong dominance of the electrostatic energy component within the interaction energy. The total NPA charges on the anions indicate an increase in the solvation due to hydration. From AIM analysis, excellent linear inverse correlation is observed for both the electron density and Laplacian of the electron density with respect to the hydrogen bond length. Natural bonding orbital analysis (NBO) predicts large charge transfer between the OH(-) anion and the water molecules.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Deep ultraviolet resonance Raman excitation enables explosives detection

We measured the 229 nm absolute ultraviolet (UV) Raman cross-sections of the explosives trinitrotoluene (TNT), pentaerythritol tetranitrate (PETN), cyclotrimethylene-trinitramine (RDX), the chemically related nitroamine explosive HMX, and ammonium nitrate in solution. The 229 nm Raman cross-sections are 1000-fold greater than those excited in the near-infrared and visible spectral regions. Deep UV resonance Raman spectroscopy enables detection of explosives at parts-per-billion (ppb) concentrations and may prove useful for stand-off spectroscopic detection of explosives.

Solution structure of NaNO3 in water: diffraction and molecular dynamics simulation study

The structure of a series of aqueous sodium nitrate solutions (1.9-7.6 M) was studied using a combination of experimental and theoretical methods. The results obtained from diffraction (X-ray, neutron) and molecular dynamics simulation have been compared and the capabilities and limitations of the methods in describing solution structure are discussed. For the solutions studied, diffraction methods were found to perform very well in description of hydration spheres of the sodium ion but do not yield detailed structural information on the anion's hydration structure. Molecular dynamics simulations proved to be a suitable tool in the detailed interpretation of the hydration sphere of ions, ion pair formation, and bulk structure of solutions.