Nitric acid trihydrate (NAT) in polar stratospheric clouds

Effect of Hydration on Electron Attachment to Methanesulfonic Acid Clusters

We report an experimental and computational study of the electron-induced chemistry of methanesulfonic acid (MSA, MeSO₃H) in clusters. We combine the mass spectra after the 70 eV electron ionization with the negative ion spectra after electron attachment (EA) at low electron energies of 0-15 eV of the MSA molecule, small MSA clusters, and microhydrated MSA clusters to reveal the solvation effects. The MSA/He coexpansion only generates small MSA clusters with up to four molecules, but adding water substantially hydrates the MSA clusters, resulting in clusters composed of 1-2 MSA molecules accompanied by quite a few water molecules. The clustering strongly suppresses the fragmentation of the MSA molecules upon both the positive ionization and EA. The electron-energy-dependent ion yield for different negative ions is measured. For the MSA molecule and pure MSA clusters, EA leads to an H-abstraction yielding MeSO₃^-. It proceeds efficiently at low electron energies below 2 eV with a shoulder at 3-4 eV and a broad, almost 2 orders of magnitude weaker, peak around 8 eV. The hydrated (H₂O)_nMeSO₃^- ions with n ≤ 3 exhibit only a broad peak around 7 eV similar to EA of pure water clusters. Thus, for the small clusters, the electron attachment and hydrogen abstraction from water occur. On the other hand, the larger clusters with n > 4 display a peak below 2 eV, which quickly dominates the spectrum with increasing n. This peak is related to the formation of the H₃O⁺·MeSO₃^- ion pair upon hydration and subsequent dipole-supported electron attachment followed by the hydronium neutralization and H₃O^• radical dissociation. The size-resolved experimental data indicate that the ionic dissociation of MSA starts to occur in the neutral MeSO₃H(H₂O)_N clusters with about four water molecules.

Isomer-specific cryogenic ion vibrational spectroscopy of the D2 tagged Cs+(HNO3)(H2O)n=0-2 complexes: ion-driven enhancement of the acidic H-bond to water

We report how the binary HNO3(H2O) interaction is modified upon complexation with a nearby Cs+ ion. Isomer-selective IR photodissociation spectra of the D2-tagged, ternary Cs+(HNO3)H2O cation confirms that two structural isomers are generated in the cryogenic ion source. In one of these, both HNO3 and H2O are directly coordinated to the ion, while in the other, the water molecule is attached to the OH group of the acid, which in turn binds to Cs+ with its -NO2 group. The acidic OH stretching fundamental in the latter isomer displays a ∼300 cm-1 red-shift relative to that in the neutral H-bonded van der Waals complex, HNO3(H2O). This behavior is analyzed with the aid of electronic structure calculations and discussed in the context of the increased effective acidity of HNO3 in the presence of the cation.

A systematic ab initio optimization of monohydrates of HCl•HNO3•H2SO4 aggregates

Hydrates of HCl, HNO₃ and H₂SO₄ involved in polar stratospheric clouds capture the attention of researchers due to the mixtures composed with them. The molecular aggregates generated with these strong acids show different behaviors, geometries and nucleation reactions at atmospheric temperatures. Here is presented a systematic ab initio optimization study of monohydrates of HCl•HNO₃•H₂SO₄ using the Density Functional Theory, by means of geometry optimizations carried out with B3LYP hybrid method and aug-cc-pVTZ basis set, a high level of theory, within Gaussian 09 program. This systematic optimization procedure consists to situate systematically the H₂O molecule around the cluster in study, on the favorable positions to develop higher quantity of hydrogen bonds as possible, in order to obtain major quantity of different electronic structures of these monohydrates. Applying this systematic optimization methodology over previously optimized complexes of HCl, HNO₃ and H₂SO₄, the present theoretical approach provides thirty-two different optimized electronic structures of monohydrates that were yielded from seven initial groups of (HCl•HNO₃•H₂SO₄)-complex, placing the H₂O in eight positions around them. Moreover, their Infrared spectra have been predicted for all (HCl•HNO₃•H₂SO₄)-monohydrates achieved. Likewise, It is shown the outcomes of the electronic energies, relative Gibbs free energies, Infrared spectra, the wavenumbers of hydrogen bonds, inter-monomeric parameters, electronic structures of (HCl•HNO₃•H₂SO₄)-monohydrates. These monohydrates could be considered precursors of the atmospheric heterogeneous nucleation reactions. These results can be useful to experimentalists of Catalysis, Astrophysics, Corrosion of metals and ceramics, aromatic compounds reactions, even environmental pollution and industrial smog.

Persistent Water-Nitric Acid Condensate with Saturation Water Vapor Pressure Greater than That of Hexagonal Ice

A laboratory chilled mirror hygrometer (CMH), exposed to an airstream containing water vapor (H2O) and nitric acid (HNO3), has been used to demonstrate the existence of a persistent water-nitric acid condensate that has a saturation H2O vapor pressure greater than that of hexagonal ice (Ih). The condensate was routinely formed on the mirror by removing HNO3 from the airstream following the formation of an initial condensate on the mirror that resembled nitric acid trihydrate (NAT). Typical conditions for the formation of the persistent condensate were a H2O mixing ratio greater than 18 ppm, pressure of 128 hPa, and mirror temperature between 202 and 216 K. In steady-state operation, a CMH maintains a condensate of constant optical diffusivity on a mirror through control of only the mirror temperature. Maintaining the persistent condensate on the mirror required that the mirror temperature be below the H2O saturation temperature with respect to Ih by as much as 3 K, corresponding to up to 63% H2O supersaturation with respect to Ih. The condensate was observed to persist in steady state for up to 16 h. Compositional analysis of the condensate confirmed the co-condensation of H2O and HNO3 and thereby strongly supports the conclusion that the Ih supersaturation is due to residual HNO3 in the condensate. Although the exact structure or stoichiometry of the condensate could not be determined, other known stable phases of HNO3 and H2O are excluded as possible condensates. This persistent condensate, if it also forms in the upper tropical troposphere, might explain some of the high Ih supersaturations in cirrus and contrails that have been reported in the tropical tropopause region.

Metastable Nitric Acid Trihydrate in Ice Clouds

The composition of high‐altitude ice clouds is still a matter of intense discussion. The constituents in question are ice and nitric acid hydrates, but the exact phase composition of clouds and its formation mechanisms are still unknown. In this work, conclusive evidence for a long‐predicted phase, alpha‐nitric acid trihydrate (alpha‐NAT), is presented. This phase was characterized by a combination of X‐ray and neutron diffraction experiments, allowing a convincing structure solution. Furthermore, vibrational spectra (infrared and inelastic neutron scattering) were recorded and compared with theoretical calculations. A strong interaction between water ice and alpha‐NAT was found, which explains the experimental spectra and the phase‐transition kinetics. On the basis of these results, we propose a new three‐step mechanism for NAT formation in high‐altitude ice clouds.

Relative stabilities of HCl•H2SO4•HNO3 aggregates in polar stratospheric clouds

Strong acids such as HCl (C), HNO3 (N) and H2SO4 (S) acquire relevance in Polar Stratospheric Clouds (PSCs) and aerosols in which nucleation processes occur. Ab initio quantum chemical studies of aggregates were performed for these strong acids. Structures were calculated using DFT methods with the B3LYP hybrid functional and aug-cc-pVTZ basis set. As an initial constraint, an H2SO4 moiety was placed in all candidate structures. A total of 11 optimized structures was found: a global minimum (CSN-a) plus ten local minima on the Potential Energy Surface (PES). The global minimum aggregate gave four hydrogen bonds, yielding a hexagonal ring in its structure. HNO3 acts as proton donor in all clusters; nevertheless, using trans-H2SO4 as the proton donor yielded the most stable structures, whereas HCl acts mainly as a proton donor/acceptor. Real harmonic frequencies, IR spectra, and inter-monomeric parameters were obtained. CSN-a symmetric stretching modes were shifted to 2805.56 cm(-1) and 3520.00 cm(-1) for H-Cl modes, while O-H modes shifted to 3256.87 cm(-1) and 3362.47 cm(-1). On the other hand, relative stabilities improved for 5 of the 11 aggregates when the temperature decreased from 298 K to 210 K, 195 K and 188 K. The aggregate CSN-f remained unstable only at 210 K. Moreover, the relative Gibbs free energy, ΔG(0-298K) was -9.26 kcalmol(-1) with respect to CSN-a; relative reaction Gibbs free energy [Δ(ΔG)] values ranged from 0.0 at 298 K, to -6.9 kcalmol(-1) at 188 K. It seems that CSN aggregates remain slightly more stable than CNS aggregates with a HNO3 moiety when the temperature decreases from 298 to 188 K. Five structures remained relatively stable under both study conditions.

Quantum chemical study of atmospheric aggregates: HCl•HNO3•H2SO4

HCl, HNO3 and H2SO4 are implicated in atmospheric processes in areas such as polar stratospheric clouds in the stratosphere. Ternary complexes of HCl, HNO3 and H2SO4 were investigated by ab initio calculations at B3LYP level of theory with aug-cc-pVTZ and aug-cc-pVQZ basis sets, taking into account basis set superposition error (BSSE). The results were assessed in terms of structures (five hexagonal cyclic structures and two quasi-pentagonal cyclic structures), inter-monomeric parameters (all ternary complexes built three hydrogen bonds), energetics (seven minima obtained), infrared harmonic vibrational frequencies (red shifting of complexes from monomers), and relative stability of complexes, which were favorable when the temperature decreases under stratospheric conditions, from 298 K to 188 K, and in concrete, at 210 K, 195 K and 188 K.

Ionization of Nitric Acid on Crystalline Ice: The Role of Defects and Collective Proton Movement

Ionization of nitric acid (HNO3) on a model ice surface is studied using ab initio molecular dynamics at temperatures of 200 and 40 K with a surface slab model that consists of the ideal ice basal plane with locally optimized and annealed defects. Pico- and subpicosecond ionization of nitric acid can be achieved in the defect sites. Key features of the rapid ionization are (a) the efficient solvation of the polyatomic nitrate anion, by stealing hydrogen bonds from the weakened hydrogen bonds at defect sites, (b) formation of contact ion pairs to stable "presolvated" molecular species that are present at the defects,

The crystal structures of the low-temperature and high-pressure polymorphs of nitric acid

A new high-pressure phase of pure nitric acid (HNO(3)) has been characterised at 1.6 GPa at room temperature by high-pressure neutron powder and X-ray single-crystal diffraction techniques. This is the first crystalline phase obtained upon compression of liquid nitric acid at room temperature and appears to be the stable phase up to pressures of at least 4 GPa. The crystal structure of this new phase shows some similarities to that of the low-temperature phase of nitric acid at ambient pressure, which has been redetermined as part of this study. Both structures share a herringbone packing of hydrogen-bonded molecular catemers, although the presence of disorder within the hydrogen bonds within one of the catemers of the low-temperature phase makes its structure comparatively more complex.

Theoretical study on hydrogen-bond effects in IR spectra of high- and low-temperature phases of nitric acid dihydrate

The low- and high-temperature phases (alpha and beta, respectively) of solid nitric acid dihydrate (NAD) are studied in depth by DFT methods. Each phase contains two types of complex structures (H(3)O(+)) x (H(2)O), designated A and B, with different hydrogen-bonding (HB) characteristics. The theoretical study reveals that type A complexes are weakly bound and could be described as (H(3)O)(+) and H(2)O aggregates, with decoupled vibrational modes, whereas in type B structures the proton is situated close to the centre of the O...O bond and induces strong vibrational coupling. The proton-transfer mode is predicted at quite different wavenumbers in each complex, which provides an important differentiating spectral feature, together with splitting of some bands in beta-NAD. Theoretical spectra are estimated by using two GGA parameterizations, namely, PBE and BLYP. The potential-energy surface for each type of HB in NAD is also studied, as is the spectral influence of displacement of the shared H atom along the O-O bond. The results are compared to literature infrared spectra recorded by different techniques, namely, transmission and reflection-absorption, with both normal and tilted incident radiation. This work provides a thorough assignment of the observed spectra, and predictions for some spectra not yet available. The usefulness of high-level theoretical calculations as performed herein to discriminate between two phases of a solid crystal is thus evidenced.

Do NAD and NAT Form in Liquid Stratospheric Aerosols by Pseudoheterogeneous Nucleation?

Laboratory data of the freezing of nitric acid hydrates (NAD, NAT) from HNO(3)/H(2)O and HNO(3)/H(2)SO(4)/H(2)O solution droplets have been evaluated with respect to a "pseudoheterogeneous" (surface-induced) nucleation mechanism of NAD and NAT, which has been argued to possibly lead to the formation of polar stratospheric clouds (PSCs). In addition, a parametrization of pseudoheterogeneous nucleation of NAD and NAT suggested recently (Tabazadeh et al. J. Phys. Chem. A 2002, 106, 10238-10246) has been analyzed, showing that this parametrization should not be used in stratospheric modeling studies. The analysis of several laboratory data sets yields an upper limit of the pseudoheterogeneous nucleation rate coefficient of NAD of 2.2 x 10(-5) cm(-2) s(-1). In contrast, the upper limit of the pseudoheterogeneous nucleation rate coefficient of NAT could not be constrained satisfactorily, since formation of NAT has not been observed at stratospheric conditions in laboratory experiments applying small droplets. Maximum NAD production rates of 9.6 x 10(-9) cm(-3) (air) h(-1) in the stratosphere have been estimated assuming a pseudoheterogeneous nucleation mechanism that is constrained by the experimental observations. If maximum NAD supersaturation persisted for 4 weeks in the polar stratosphere the corresponding NAD particle number densities are estimated to be about 6 x 10(-6) cm(-3). These particle number densities are 3 orders of magnitude lower than particle number densities recently observed in the stratosphere. In conclusion, on the basis of laboratory data it is found that a pseudoheterogeneous nucleation mechanism is not sufficient to explain recent observations of large nitric acid containing particles in the polar stratosphere.

Orientation Effects on Nitric Acid Dihydrate Films

An investigation of orientation effects in films of nitric acid dihydrate (NAD) is presented, based on a systematic study of transmission and reflection-absorption infrared (RAIR) spectra of samples of varying thickness. The samples are prepared by vapor deposition on Ge (for transmission spectroscopy) and on Al substrates (for RAIR spectroscopy) at 175 K to produce crystalline alpha-NAD films. Transmission spectra were recorded at normal incidence, and RAIR spectra were recorded at a grazing angle of 75 degrees, with polarized radiation. The observed spectra are compared with predictions of a classical Fresnel model, to test the available optical indices of NAD, which are of great importance for the accurate interpretation of data from remote sensing measurements. Whereas the procedure yields satisfactory results for transmission and s-polarized RAIR spectra, it is found that the agreement is not acceptable for p-polarized RAIR spectra. An explanation is suggested in terms of a preferential alignment of the films, with the (10-1) crystallographic plane of the crystal situated parallel to the substrate. The infrared activity of a band at approximately 1170 cm(-1) is explained in terms of a preferential orientation of the crystal domains in the film.

Theoretical Calculations of Refractive Indices and Optical Effects in Spectra of Nitric Acid and Nitric Acid Monohydrate Crystals

The theoretical infrared refractive indices of two systems related to atmospheric research, nitric acid (NA) and nitric acid monohydrate (NAM) crystals, have been computed using a methodology based on first-principles. The effects of lack of coherence in the infrared beam in RAIR and transmission spectra have also been treated using a model based on classical optics. The optical constants of NA crystals are presented for the first time; the results on NAM are compared to empirical values previously published with good general agreement. With the optical constants of NA, polarized reflection-absorption infrared spectra are predicted and compared to experimental spectra recorded also for the first time, for a set of varying film thickness. The global agreement is satisfactory. The effects of a number of experimental factors in transmission spectra of NAM are assessed, in an attempt to explain observed differences among experimental spectra. It is concluded that the spectral disparities are probably due to differences in the nature of the samples.

Influence of Monolayer Amounts of HNO3 on the Evaporation Rate of H2O over Ice in the Range 179 to 208 K: A Quartz Crystal Microbalance Study

The evaporation flux J(ev) of H2O from thin H2O ice films containing between 0.5 and 7 monolayers of HNO3 has been measured in the range 179 to 208 K under both molecular and stirred flow conditions in isothermal experiments. FTIR absorption of the HNO3/H2O condensate revealed the formation of metastable alpha-NAT (HNO(3).3H2O) converting to stable beta-NAT at 205 K. After deposition of HNO3 for 16-80 s on a 1 mum thick pure ice film at a deposition rate in the range (6-60) x 10(12) molecules s(-1) the initial evaporative flux J(ev)(H2O) was always that of pure ice. J(ev)(H2O) gradually decreased with the evaporation of H2O and the concomitant increase of the average mole fraction of HNO3, chi(HNO3), indicating the presence of an amorphous mixture of H2O/HNO3 that is called complexed or (c)-ice whose vapor pressure is that of pure ice. The final value of J(ev) was smaller by factors varying from 2.7 to 65 relative to pure ice. Depending on the doping conditions and temperature of the ice film the pure ice thickness d(D) of the ice film for which J(ev) < 0.85J(ev)(pure ice) varied between 130 and 700 nm compared to the 1000 nm thick original ice film at 208 and 191 K, respectively, in what seems to be an inverse temperature dependence. There exist three different types of H2O molecules under the present experimental conditions, namely (a) free H2O corresponding to pure ice, (b) complexed H2O or c-ice, and (c) H2O molecules originating from the breakup of NAT or amorphous H2O/HNO3 mixtures. The significant decrease of J(ev)(H2O) with increasing chi(HNO3) leads to an increase of the evaporative lifetime of atmospheric ice particles in the presence of HNO3 and may help explain the occurrence of persistent and/or large contaminated ice particles at certain atmospheric conditions.

Determination of polar stratospheric cloud particle refractive indices by use of in situ optical measurements and T-matrix calculations

A new algorithm to infer structural parameters such as refractive index and asphericity of cloud particles has been developed by use of in situ observations taken by a laser backscattersonde and an optical particle counter during balloon stratospheric flights. All three main particles, liquid, ice, and a no-ice solid (NAT, nitric acid trihydrate) of polar stratospheric clouds, were observed during two winter flights performed from Kiruna, Sweden. The technique is based on use of the T-matrix code developed for aspherical particles to calculate the backscattering coefficient and particle depolarizing properties on the basis of size distribution and concentration measurements. The results of the calculations are compared with observations to estimated refractive indices and particle asphericity. The method has also been used in cases when the liquid and solid phases coexist with comparable influence on the optical behavior of the cloud to estimate refractive indices. The main results prove that the index of refraction for NAT particles is in the range of 1.37-1.45 at 532 nm. Such particles would be slightly prolate spheroids. The calculated refractive indices for liquid and ice particles are 1.51-1.55 and 1.31-1.33, respectively. The results for solid particles confirm previous measurements taken in Antarctica during 1992 and obtained by a comparison of lidar and optical particle counter data.

Optical and physical properties of stratospheric aerosols from balloon measurements in the visible and near-infrared domains. II. Comparison of extinction, reflectance, polarization, and counting measurements

The physical properties of stratospheric aerosols can be retrieved from optical measurements involving extinction, radiance, polarization, and counting. We present here the results of measurements from the balloonborne instruments AMON, SALOMON, and RADIBAL, and from the French Laboratoire de Météorologie Dynamique and the University of Wyoming balloonborne particle counters. A cross comparison of the measurements was made for observations of background aerosols conducted during the polar winters of February 1997 and January-February 2000 for various altitudes from 13 to 19 km. On the one band, the effective radius and the total amount of background aerosols derived from the various sets of data are similar and are in agreement with pre-Pinatubo values. On the other hand, strong discrepancies occur in the shapes of the bimodal size distributions obtained from analysis of the raw measurement of the various instruments. It seems then that the log-normal assumption cannot fully reproduce the size distribution of background aerosols. The effect ofthe presence of particular aerosols on the measurements is discussed, and a new strategy for observations is proposed.