Computational Investigation of the Solvation of Nitric Acid: Formation of the NO3- and H3O+ Ion Pair

Understanding the Formation and Growth of New Atmospheric Particles at the Molecular Level through Laboratory Molecular Beam Experiments

Atmospheric new particle formation (NPF), which exerts comprehensive implications for climate, air quality and human health, has received extensive attention. From molecule to cluster is the initial and most important stage of the nucleation process of atmospheric new particles. However, due to the complexity of the nucleation process and limitations of experimental characterization techniques, there is still a great uncertainty in understanding the nucleation mechanism at the molecular level. Laboratory-based molecular beam methods can experimentally implement the generation and growth of typical atmospheric gas-phase nucleation precursors to nanoscale clusters, characterize the key physical and chemical properties of clusters such as structure and composition, and obtain a series of their physicochemical parameters, including association rate coefficients, electron binding energy, pickup cross section and pickup probability and so on. These parameters can quantitatively illustrate the physicochemical properties of the cluster, and evaluate the effect of different gas phase nucleation precursors on the formation and growth of atmospheric new particles. We review the present literatures on atmospheric cluster formation and reaction employing the experimental method of laboratory molecular beam. The experimental apparatuses were classified and summarized from three aspects of cluster generation, growth and detection processes. Focus of this review is on the properties of nucleation clusters involving different precursor molecules of water, sulfuric acid, nitric acid and NxOy, respectively. We hope this review will provide a deep insight for effects of cluster physicochemical properties on nucleation, and reveal the formation and growth mechanism of atmospheric new particle at the molecular level.

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.

Does HNO3 dissociate on gas-phase ice nanoparticles?

We investigated the dissociation of nitric acid on large water clusters (H2O)N, N̄ ≈ 30-500, i.e., ice nanoparticles with diameters of 1-3 nm, in a molecular beam. The (H2O)N clusters were doped with single HNO3 molecules in a pickup cell and probed by mass spectrometry after a low-energy (1.5-15 eV) electron attachment. The negative ion mass spectra provided direct evidence for HNO3 dissociation with the formation of NO3-⋯H3O+ ion pairs, but over half of the observed cluster ions originated from non-dissociated HNO3 molecules. This behavior is in contrast with the complete dissociation of nitric acid on amorphous ice surfaces above 100 K. Thus, the proton transfer is significantly suppressed on nanometer-sized particles compared to macroscopic ice surfaces. This can have considerable implications for heterogeneous processes on atmospheric ice particles.

DFT Study of Microsolvated [NO3·(H2O)n]- (n = 1-12) Clusters and Molecular Dynamics Simulation of Nitrate Solution

Investigation of the process of the NO₃^- anion solvation is central to understanding the chemical and physical properties of its aqueous solutions. The importance of this topic can be seen in atmospheric chemistry, as well as in nuclear waste processing research. In this work, we used a particle swarm optimization technique driven by density functional theory to sample the potential energy surface of various microsolvated [NO₃·(H₂O)_n]^- (n = 1-12) clusters. We found that the charge transfer plays a crucial role in the stabilization of the investigated species. Moreover, by conducting ab initio molecular dynamics simulations, we showed that at low concentrations (∼0.2 M) the NO₃^- species tend to be located on the surface of water solution. We also observed that the contact ion pair K⁺-NO₃^- undergoes a fast dissociation and each of the ions is solvated separately. As a result, from our calculations, we expect that at low concentration there could be oppositely signed concentration gradients for NO₃^- and K⁺ ions in a thin water film.

DFT Analysis of Organotin Catalytic Mechanisms in Dehydration Esterification Reactions for Terephthalic Acid and 2,2,4,4-Tetramethyl-1,3-cyclobutanediol

Polyesters synthesized from 2,2,4,4-tetramethyl-1,3-cyclobutanediol (TMCD) and terephthalic acid (TPA) are improved alternatives to toxic polycarbonates based on bisphenol A. In this work, we use ωB97X-D/LANL2DZdp calculations, in the presence of a benzaldehyde polarizable continuum model solvent, to show that esterification of TMCD and TPA will reduce and subsequently dehydrate a dimethyl tin oxide catalyst, becoming ligands on the now four-coordinate complex. This reaction then proceeds most plausibly by an intramolecular acyl-transfer mechanism from the tin complex, aided by a coordinated proton donor such as hydronium. These findings are a key first step in understanding polyester synthesis and avoiding undesirable side reactions during production.

Predictive models for the phase behaviour and solution properties of weak electrolytes: nitric, sulphuric, and carbonic acids

The distribution of ionic species in electrolyte systems is important in many fields of science and engineering, ranging from the study of degradation mechanisms to the design of systems for electrochemical energy storage. Often, other phenomena closely related to ionic speciation, such as ion pairing, clustering and hydrogen bonding, which are difficult to investigate experimentally, are also of interest. Here, we develop an accurate molecular approach, accounting for reactions as well as association and ion pairing, to deliver a predictive framework that helps validate experiment and guides future modelling of speciation phenomena of weak electrolytes. We extend the SAFT-VRE Mie equation of state [D. K. Eriksen et al., Mol. Phys., 2016, 114, 2724-2749] to study aqueous solutions of nitric, sulphuric, and carbonic acids, considering complete and partially dissociated models. In order to incorporate the dissociation equilibria, correlations to experimental data for the relevant thermodynamic equilibrium constants of the dissociation reactions are taken from the literature and are imposed as a boundary condition in the calculations. The models for water, the hydronium ion, and carbon dioxide are treated as transferable and are taken from our previous work. We present new molecular models for nitric acid, and the nitrate, bisulfate, sulfate, and bicarbonate anions. The resulting framework is used to predict a range of phase behaviour and solution properties of the aqueous acids over wide ranges of concentration and temperature, including the degree of dissociation, as well as the activity coefficients of the ionic species, and the activity of water and osmotic coefficient, density, and vapour pressure of the solutions. The SAFT-VRE Mie models obtained in this manner provide a means of elucidating the mechanisms of association and ion pairing in the systems studied, complementing the experimental observations reported in the literature.

Mass spectrometry of aerosol particle analogues in molecular beam experiments

Nanometer-size particles such as ultrafine aerosol particles, ice nanoparticles, water nanodroplets, etc, play an important, however, not yet fully understood role in the atmospheric chemistry and physics. These species are often composed of water with admixture of other atmospherically relevant molecules. To mimic and investigate such particles in laboratory experiments, mixed water clusters with atmospherically relevant molecules can be generated in molecular beams and studied by various mass spectrometric methods. The present review demonstrates that such experiments can provide unprecedented details of reaction mechanisms, and detailed insight into the photon-, electron-, and ion-induced processes relevant to the atmospheric chemistry. After a brief outline of the molecular beam preparation, cluster properties, and ionization methods, we focus on the mixed clusters with various atmospheric molecules, such as hydrated sulfuric acid and nitric acid clusters, N_x O_y and halogen-containing molecules with water. A special attention is paid to their reactivity and solvent effects of water molecules on the observed processes.

A detailed hydrogen bonding analysis on the compositions of H2SO4/HNO3/H2O ternary systems: A computational study

Hydrogen bonding properties of H₂SO₄/HNO₃/H₂O ternary molecular clusters have been studied by means of structural, energetic, topological, and spectroscopic perspectives. The roles of the hydrogen bonds in the formation of these clusters are considered according to the molecule positions (proton donor or proton acceptor) in the clusters. 33 stable conformers were identified on the potential energy surface. The global minimum one was obtained when HNO₃ donates a proton to H₂SO₄, however, the cooperativity effect contribution was found to be significant for the cluster where HNO₃ acts as a proton donor to an H₂O molecule.

Hydrogen-bonding behavior of various conformations of the HNO3…(CH3OH)2 ternary system

Nine minima were found on the intermolecular potential energy surface for the ternary system HNO₃(CH₃OH)₂ at the MP2/aug-cc-pVDZ level of theory. The cooperative effect, which is a measure of the hydrogen-bonding strength, was probed in these nine conformations of HNO₃…(CH₃OH)₂. The results are discussed here in terms of structures, energetics, infrared vibrational frequencies, and topological parameters. The cooperative effect was observed to be an important contributor to the total interaction energies of the cyclic conformers of HNO₃…(CH₃OH)₂, meaning that it cannot be neglected in simulations in which the pair-additive potential is applied. Graphical abstract The H-bonding behavior of various conformations of the HNO₃(CH₃OH)₂ trimer was investigated.

Electron-triggered chemistry in HNO3/H2O complexes

Polar stratospheric clouds, which consist mainly of nitric acid containing ice particles, play a pivotal role in stratospheric chemistry. We investigate mixed nitric acid-water clusters (HNO3)m(H2O)n, m ≈ 1-6, n ≈ 1-15, in a laboratory molecular beam experiment using electron attachment and mass spectrometry and interpret our experiments using DFT calculations. The reactions are triggered by the attachment of free electrons (0-14 eV) which leads to subsequent intracluster ion-molecule reactions. In these reactions, the nitrate anion NO3- turns out to play the central role. This contradicts the electron attachment to the gas-phase HNO3 molecule, which leads almost exclusively to NO2-. The nitrate containing clusters are formed through at least three different reaction pathways and represent terminal product ions in the reaction cascade initiated by the electron attachment. Besides, the complex reaction pathways represent a new hitherto unrecognized source of atmospherically important OH and HONO molecules.

Microhydration of a benzoic acid molecule and its dissociation

C₆H₅COOH requires at least eight H₂O molecules for its dissociation. A linear correlation is observed between pK_a and n_H2O; an unknown pK_a can thus be predicted.

Proton transfer and autoionization in HNO3·HCl·(H2O)n particles

The structure and spectroscopic properties of clusters of HNO(3)·HCl·(H(2)O)(n), with n = 1 to 6, have been calculated at the MP2/aug-cc-pVDZ level of theory. Altogether 22 different clusters have been found as stable structures, with minima in their potential energy surfaces. The clusters can be grouped in families with the same number of water molecules, and with close aggregation energies within each family. The addition of each new water molecule increments the aggregation energy of the clusters by a nearly constant value of 76.2 ± 0.1 Hartree. The proton transfer parameter and the coordination number of HNO(3) and HCl in each cluster have been evaluated, and the wavenumber shifts for the X(-)-H(+) vibration from the corresponding mode in the isolated molecules have also been predicted. These values allow classification of the acidic species in the clusters into three types, characterized by the strength of the hydrogen bond and the degree of ionization. A correspondence is found between the coordination number of HNO(3) and the magnitude of the X(-)-H(+) vibrational shift.

Reaction Paths of the Water-Assisted Solvolysis of N,N-Dimethylformamide

Density functional theory calculations were conducted on the title reactions with explicit inclusion of a variety of water molecules, H-CO-NMe2+MeOH+(H2O)n-->H-CO-OMe+HNMe2+(H2O)n. Geometries of transition states, reactant-like complexes and product-like ones were determined by the use of RB3LYP/6-31G(d) SCRF=dipole. Concerted paths were examined with n=0-3. Their Gibbs activation energies are larger than the experimental value. Stepwise paths were also investigated with n=2-4. The n=4 model has the energy close to the experimental value. However, when the catalytic water molecules were added to the n=4 one, the stepwise path was switched to the concerted one. A systematic comparison of the concerted path with n=2+1, 2+2, 2+3, 2+4, 2+5, 2+4+4, and 2+5+5 models was made, and the water-dimer based reaction path was found to be most favorable. The contrast between the concerted path of the amide solvolysis (and hydrolysis) and the stepwise one of the ester hydrolysis was discussed in terms of the frontier-orbital theory.