Structure and dynamics of dissociated and undissociated forms of nitric acid and their implications in interfacial mass transfer: insights from molecular dynamics simulations

Nitric acid (HNO3) is widely used in various chemical and nuclear industries. Therefore, it is important to develop an understanding of the different forms of nitric acid for its practical applications. Molecular dynamics (MD) simulation is one of the best tools to investigate the behavior of concentrated nitric acid in aqueous solution with various forms together with pure nitric acid to identify a suitable model of nitric acid for use in simulations of biphasic systems for interfacial mass transfer. The Mulliken partial charge embedded OPLS-AA force field was used to model the neutral nitric acid, hydronium ion and nitrate ion, and it was found that the Mulliken partial charge embedded force field works quite well. The computed density of the dissociated and mixed-form acid was in good agreement with the experimental values. In water, the HNO3 molecule was seen to be coordinated with three water molecules in the first sphere of coordination. The distribution of water surrounding the HNO3 molecule and nitrate ion was corroborated by the DFT-optimized hydrated cluster. The calculated diffusivity values of the neutral acid and ions were significantly higher in the mixed form of nitric acid, which is an important dynamic quantity controlling the kinetics of the liquid-liquid interfacial extraction. The structural analysis revealed that the local aggregation is minimized when both forms of acid are present together in the solution. The water-ion and water-neutral acid interactions were predicted to be enhanced, as confirmed by H-bond studies. The shear viscosity of the mixed acid exhibited excellent agreement with the experimental values, which again confirms the consideration of the mixed form of nitric acid. The simulated value of surface tension for the mixed form of acid also appeared to be quite accurate based on the surface tension of water. The mixed form of nitric acid comprising both forms of acid is the best representation for nitric acid to be considered for MD simulations of biphasic systems. The mixed form of nitric acid established that the concentrated nitric acid may not be present either in the fully dissociated form or fully undissociated form in the solution.

Interaction of Nitrite Ions with Hydrated Portlandite Surfaces: Atomistic Computer Simulation Study

The nitrite admixtures in cement and concrete are used as corrosion inhibitors for steel reinforcement and also as anti-freezing agents. The characterization of the protective properties should account for the decrease in the concentration of free NO₂^- ions in the pores of cement concretes due to their adsorption. Here we applied the classical molecular dynamics computer simulation approach to quantitatively study the molecular scale mechanisms of nitrite adsorption from NaNO₂ aqueous solution on a portlandite surface. We used a new parameterization to model the hydrated NO₂^- ions in combination with the recently upgraded ClayFF force field (ClayFF-MOH) for the structure of portlandite. The new NO₂^- parameterization makes it possible to reproduce the properties of hydrated NO₂^- ions in good agreement with experimental data. In addition, the ClayFF-MOH model improves the description of the portlandite structure by explicitly taking into account the bending of Ca-O-H angles in the crystal and on its surface. The simulations showed that despite the formation of a well-structured water layer on the portlandite (001) crystal surface, NO₂^- ions can be strongly adsorbed. The nitrite adsorption is primarily due to the formation of hydrogen bonds between the structural hydroxyls on the portlandite surface and both the nitrogen and oxygen atoms of the NO₂^- ions. Due to that, the ions do not form surface adsorption complexes with a single well-defined structure but can assume various local coordinations. However, in all cases, the adsorbed ions did not show significant surface diffusional mobility. Moreover, we demonstrated that the nitrite ions can be adsorbed both near the previously-adsorbed hydrated Na⁺ ions as surface ion pairs, but also separately from the cations.

Hydrogen Bonding in Liquid Ammonia

The nature of hydrogen bonding in condensed ammonia phases, liquid and crystalline ammonia has been a topic of much investigation. Here, we use quantum molecular dynamics simulations to investigate hydrogen bond structure and lifetimes in two ammonia phases: liquid ammonia and crystalline ammonia-I. Unlike liquid water, which has two covalently bonded hydrogen and two hydrogen bonds per oxygen atom, each nitrogen atom in liquid ammonia is found to have only one hydrogen bond at 2.24 Å. The computed lifetime of the hydrogen bond is t ≅ 0.1 ps. In contrast to crystalline water-ice, we find that hydrogen bonding is practically nonexistent in crystalline ammonia-I.

Interpretation of the X-Ray Emission Spectra of Liquid Water through Temperature and Isotope Dependence

The interpretation of x-ray emission spectroscopy (XES) spectra in terms of their sensitivity to the hydrogen bonding and the consequent microheterogeneity in liquid water has been debated over a decade. To shed a light on this problem, we report the theoretical reproduction of the debated 1b_{1} peaks observed in the XES spectra of liquid water using semiclassical Kramers-Heisenberg formula. The essence of the temperature and isotope dependence of the 1b_{1} double peaks is explained by molecular dynamics simulations including full vibrational (O─H stretching, bending, and) modes, rotational combined with the density functional theory and core-hole induced dynamics. Some inconsistencies exist with the experimental XES profile, which illustrates the need to employ a more precise theoretical calculations for both geometry sampling and electronic structure using a more sophisticated procedure.

Dielectric Constant of Liquid Water Determined with Neural Network Quantum Molecular Dynamics

The static dielectric constant ϵ_{0} and its temperature dependence for liquid water is investigated using neural network quantum molecular dynamics (NNQMD). We compute the exact dielectric constant in canonical ensemble from NNQMD trajectories using fluctuations in macroscopic polarization computed from maximally localized Wannier functions (MLWF). Two deep neural networks are constructed. The first, NNQMD, is trained on QMD configurations for liquid water under a variety of temperature and density conditions to learn potential energy surface and forces and then perform molecular dynamics simulations. The second network, NNMLWF, is trained to predict locations of MLWF of individual molecules using the atomic configurations from NNQMD. Training data for both the neural networks is produced using a highly accurate quantum-mechanical method, DFT-SCAN that yields an excellent description of liquid water. We produce 280×10^{6} configurations of water at 7 temperatures using NNQMD and predict MLWF centers using NNMLWF to compute the polarization fluctuations. The length of trajectories needed for a converged value of the dielectric constant at 0°C is found to be 20 ns (40×10^{6} configurations with 0.5 fs time step). The computed dielectric constants for 0, 15, 30, 45, 60, 75, and 90°C are in good agreement with experiments. Our scalable scheme to compute dielectric constants with quantum accuracy is also applicable to other polar molecular liquids.

Hydration Properties of HnPO4n-3 (n = 0-3) From Ab Initio Molecular Dynamics Simulations

For a comprehensive and detailed microscopic understanding of the hydration properties of primary aqueous phosphorus species of valence states V (viz., H₃PO₄, H₂PO₄^-, HPO₄^2-, and PO₄^3-), a series of extensive ab initio molecular dynamics simulations is conducted at ambient temperature. In each of these cases, the spatially resolved, three-dimensional hydration shells are computed, allowing for a direct microscopic visual understanding of the hydration shells around the species. Since these species are excellent agents for the formation of hydrogen bonds (H-bonds) in water, which determine a wide range of their structural, dynamic, and spectroscopic features, a detailed analysis of the qualitative and quantitative aspects of the H-bonds, including their lifetime calculations, is performed. Vibrational density of states (VDOS) is calculated for each of the species in solute phases, resolved for each H-bonding site, and compared against the gas-phase normal modes of H₃PO₄ for the purpose of understanding the signatures of the peaks in VDOS plots and, in particular, the effects of solvation and H-bonding mechanisms. The results are well in line with available experimental data and other recent computer-aided studies in the literature.

Effect of the Side-Chain Length in Perfluorinated Sulfonic and Phosphoric Acid-Based Membranes on Nanophase Segregation and Transport: A Molecular Dynamics Simulation Approach

The effect of side-chain length on the nanophase-segregated structure and transport in perfluorinated sulfonic acid (PFSA)-based and perfluorinated phosphoric acid (PFPA)-based membranes is investigated at 20 and 5 wt % water content conditions using a molecular dynamics simulation method. It is found using the pair correlation analysis that the longer side chain leads to more developed local water structures in the water phase at 20 wt % water content, observable in both membrane chemistries albeit more distinct in PFPA-based membranes. It is also confirmed from the structure factor analysis that large-scale nanophase segregation is enhanced with increasing side-chain length for PFPA membranes, whereas no significant change is observed at these scales for PFSA membranes. Next, it is revealed that the proton transport is increased by 0.004 S/cm in PFSA-based membranes with increasing side-chain length due to the enhanced vehicular and hopping mechanisms, whereas the proton transport in PFPA-based membranes is decreased by 0.002 S/cm despite improved nanophase segregation. As confirmed by the pair correlation function analysis, the diminished proton transport in PFPA-based membranes is attributed to the molecular association of phosphate groups with hydronium ions via hydrogen bond in the longer side-chain case, which is namely a hydronium-mediated bridge configuration. Such bridge configurations and correspondingly similar trends in proton transport are also observed at 5 wt % water content condition to a lesser extent. Our simulation study demonstrates that the proton transport is affected by the hydrogen-bonding network as well as by the nanophase segregation.

How does hydrogen bond network analysis reveal the golden ratio of water-glycerol mixtures?

Properties of water-glycerol mixtures depend closely on the water/glycerol ratio. Around the 30 mol% glycerol concentration mark, the so-called golden ratio of water-glycerol mixtures, several of the mixture's properties have observed maxima or minima, without a clear fundamental explanation. In this work, a series of molecular dynamics simulations have been performed over a wide range of water-glycerol concentrations to analyze the intermolecular hydrogen bond (H-bond) network. The collected values from simulations are justified from both a probabilistic model of H-bonding and from observing the dynamic behavior of each type of H-bonds. The populations of H-bonds that exist at a given concentration of glycerol are largely governed by the probability of one oxygen atom randomly associating with another oxygen atom. However, the H-bonds that glycerol oxygen can form are dependent on the H-bonds that are formed by the other intramolecular glycerol oxygen. Based on the dynamic analysis of each type of H-bonds, there are deviations from randomly associating with another oxygen. Water preferentially donates a hydrogen to a glycerol than to another water molecule. Yet, glycerol has a near-equal likelihood for donating a hydrogen to either another glycerol or a water. This has an effect of increasing the number of H-bonds between water and glycerol molecules and decreasing H-bonds between two water molecules. A maximum contribution of H-bonds between water and glycerol occurs around 30 mol% glycerol which is a concentration where several of the mixture's properties have an observed maxima or minima.

Dynamic Properties of Water Confined in Graphene-Based Membrane: A Classical Molecular Dynamics Simulation Study

We performed molecular dynamics simulations of water molecules inside a hydrophobic membrane composed of stacked graphene sheets. By decreasing the density of water molecules inside the membrane, we observed that water molecules form a droplet through a hydrogen bond with each other in the hydrophobic environment that stacked graphene sheets create. We found that the water droplet translates as a whole body rather than a dissipate. The translational diffusion coefficient along the graphene surface increases as the number of water molecules in the droplet decreases, because the bigger water droplet has a stronger van der Waals interaction with the graphene surface that hampers the translational motion. We also observed a longer hydrogen bond lifetime as the density of water decreased, because the hydrophobic environment limits the libration motion of the water molecules. We also calculated the reorientational correlation time of the water molecules, and we found that the rotational motion of confined water inside the membrane is anisotropic and the reorientational correlation time of confined water is slower than that of bulk water. In addition, we employed steered molecular dynamics simulations for guiding the target molecule, and measured the free energy profile of water and ion penetration through the interstice between graphene sheets. The free energy profile of penetration revealed that the optimum interlayer distance for desalination is ~10 Å, where the minimum distance for water penetration is 7 Å. With a 7 Å interlayer distance between the graphene sheets, water molecules are stabilized inside the interlayer space because of the van der Waals interaction with the graphene sheets where sodium and chloride ions suffer from a 3–8 kcal/mol energy barrier for penetration. We believe that our simulation results would be a significant contribution for designing a new graphene-based membrane for desalination.

Molecular Dynamics Simulation Studies on Structure, Dynamics, and Thermodynamics of Uranyl Nitrate Solution at Various Acid Concentrations

The structural and dynamical characteristics of uranyl ions in an aqueous acidic environment are of immense importance in the field of nuclear fuel reprocessing. In view of that, the structural and dynamical behavior of the uranyl ion in water has been investigated by performing molecular dynamics (MD) simulations using different force fields. All the force fields have depicted similar structural and dynamical properties except the free energy of hydration where the Guilbaud-Wipff (GW) model performs well over the others. The calculated density using MD simulations is found to be in excellent agreement with the measured experimental density, which ensures the accuracy of the adopted GW force field. The calculated surface tension and shear viscosity are seen to be increased with uranyl nitrate concentrations. At a higher concentration of about 4.0 mol/L, the supersaturation effect has been captured by an inflection in the plot of surface tension and shear viscosity against concentration because of the solution heterogeneity, which was correlated by an inflection in the scattering intensity observed by performing the dynamic light scattering experiment. The binding mode of nitrate ions with the uranyl ion is found to be concentration-dependent, and at higher concentration, it is predominantly monodentate.

Molecular mechanism of stabilization of thin films for improved water evaporation protection

All-atom molecular dynamics simulations and experimental characterization have been used to examine the structure and dynamics of novel evaporation-suppressing films where the addition of a water-soluble polymer to an ethylene glycol monooctadecyl ether monolayer leads to improved water evaporation resistance. Simulations and Langmuir trough experiments demonstrate the surface activity of poly(vinyl pyrrolidone) (PVP). Subsequent MD simulations performed on the thin films supported by the PVP sublayer show that, at low surface pressures, the polymer tends to concentrate at the film/water interface. The simulated atomic concentration profiles, hydrogen bonding patterns, and mobility analyses of the water-polymer-monolayer interfaces reveal that the presence of PVP increases the atomic density near the monolayer film, improves the film stability, and reduces the mobility of interfacial waters. These observations explain the molecular basis of the improved efficacy of these monolayer/polymer systems for evaporation protection of water and can be used to guide future development of organic thin films for other applications.