Reaction of Chlorine Nitrate with Hydrogen Chloride and Water at Antarctic Stratospheric Temperatures

Molecular Insights into the Spontaneous Generation of Cl2O in the Reaction of ClONO2 and HOCl at the Air-Water Interface

Chemical processes involving chlorine nitrate (ClONO2) at the surface of stratospheric aerosols are crucial to ozone depletion. Herein, we show a reaction route for the formation of Cl2O, which is a source of stratospheric chlorine, in the ClONO2 + HOCl reaction at the air-water interface. Our ab initio molecular dynamics (AIMD) simulations show that the (ClONO2)Cl···O(HOCl) halogen bond plays a key role in the reaction and is the main interaction between ClONO2 and HOCl both at the air-water interface and in the bulk liquid water. Furthermore, metadynamics-based AIMD simulations reveal two pathways: (i) The OCl fragment of HOCl binds to the Cl atom in ClONO2, resulting in the formation of Cl2O and NO3-. Simultaneously, the remaining hydrogen atom is transferred to a water molecule to form H3O+. (ii) HOCl acts as a bridge for Cl atom transfer from ClONO2 to the O atom of a water molecule, and this water molecule transfers one of its H atoms to another water molecule, forming two HOCl molecules, NO3-, and H3O+. Free-energy calculations show that the former is the energetically more favorable process. More importantly, the free-energy barrier for Cl2O formation at the air-water interface is only ∼0.8 kcal/mol, and the reaction is exothermic. These findings provide insights into the importance of fundamental chlorine chemistry and the broader implications of the aerosol air-water interface for atmospheric chemistry.

Formulation of the cosmic ray-driven electron-induced reaction mechanism for quantitative understanding of global ozone depletion

This paper formulates the cosmic ray-driven electron-induced reaction as a universal mechanism to provide a quantitative understanding of global ozone depletion. Based on a proposed electrostatic bonding mechanism for charge-induced adsorption of molecules on surfaces and on the measured dissociative electron transfer (DET) cross sections of ozone-depleting substances (ODSs) adsorbed on ice, an analytical equation is derived to give atmospheric chlorine atom concentration: [Formula: see text] where Φe is the prehydrated electron (epre-) flux produced by cosmic ray ionization on atmospheric particle surfaces, [Formula: see text] is the surface coverage of an ODS, and ki is the ODS's effective DET coefficient that is the product of the DET cross section, the lifetimes of surface-trapped epre- and Cl-, and the particle surface area density. With concentrations of ODSs as the sole variable, our calculated results of time-series ozone depletion rates in global regions in the 1960s, 1980s, and 2000s show generally good agreement with observations, particularly with ground-based ozonesonde data and satellite-measured data over Antarctica and with satellite data in a narrow altitude band at 13 to 20 km of the tropics. Good agreements with satellite data in the Arctic and midlatitudes are also found. A previously unreported effect of denitrification on ozone loss is found and expressed quantitatively. But this equation overestimates tropospheric ozone loss at northern midlatitudes and the Arctic, likely due to increased ozone production by the halogen chemistry in polluted regions. The results render confidence in applying the equation to achieve a quantitative understanding of global ozone depletion.

Mechanistic Insights into the Reactive Uptake of Chlorine Nitrate at the Air-Water Interface

It is well-known that the aqueous-phase processing of chlorine nitrate (ClONO₂) plays a crucial role in ozone depletion. However, many of the physical and chemical properties of ClONO₂ at the air-water interface or in bulk water are unknown or not understood on a microscopic scale. Here, the solvation and hydrolysis of ClONO₂ at the air-water interface and in bulk water at 300 K were investigated by classical and ab initio molecular dynamics (AIMD) simulations combined with free energy methods. Our results revealed that ClONO₂ prefers to accumulate at the air-water interface rather than in the bulk phase. Specifically, halogen bonding interactions (ClONO₂)Cl···O(H₂O) were found to be the predominant interactions between ClONO₂ and H₂O. Moreover, metadynamics-biased AIMD simulations revealed that ClONO₂ hydrolysis is catalyzed at the air-water interface with an activation barrier of only ∼0.2 kcal/mol; additionally, the difference in free energy between the product and reactant is only ∼0.1 kcal/mol. Surprisingly, the near-barrierless reaction and the comparable free energies of the reactant and product suggested that the ClONO₂ hydrolysis at the air-water interface is reversible. When the temperature is lowered from 300 to 200 K, the activation barrier for the ClONO₂ hydrolysis at the air-water interface is increased to ∼5.4 kcal/mol. These findings have important implications for the interpretation of experiments.

Effect of sodium chloride adsorption on the surface premelting of ice

We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na⁺/Cl^-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na⁺/Cl^- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na⁺/Cl^- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.

Observation of large and all-season ozone losses over the tropics

This paper reveals a large and all-season ozone hole in the lower stratosphere over the tropics (30°N–30°S) existing since the 1980s, where an O3 hole is defined as an area of O3 loss larger than 25% compared with the undisturbed atmosphere. The depth of this tropical O3 hole is comparable to that of the well-known springtime Antarctic O3 hole, whereas its area is about seven times that of the latter. Similar to the Antarctic O3 hole, approximately 80% of the normal O3 value is depleted at the center of the tropical O3 hole. The results strongly indicate that both Antarctic and tropical O3 holes must arise from an identical physical mechanism, for which the cosmic-ray-driven electron reaction model shows good agreement with observations. The whole-year large tropical O3 hole could cause a great global concern as it can lead to increases in ground-level ultraviolet radiation and affect 50% of the Earth’s surface area, which is home to approximately 50% of the world’s population. Moreover, the presence of the tropical and polar O3 holes is equivalent to the formation of three “temperature holes” observed in the stratosphere. These findings will have significances in understanding planetary physics, ozone depletion, climate change, and human health.

Photo-generated hydroxyl radicals contribute to the formation of halogen radicals leading to ozone depletion on and within polar stratospheric clouds surface

Polar stratospheric clouds (PSCs), of which the surface is a dynamic liquid water layer and might consist of aqueous HNO₃ and H₂O₂, is a well-known key meteorological condition contributing to the ozone hole in the polar stratosphere. PSCs has been considered to provide abundant surface for the heterogeneous reactions causing the formation of the Cl₂ and HOCl, which are further photolyzed into Cl and ClO radicals leading to the ozone destruction. Here we demonstrated that the sunlight drives the massive and stable production of OH radicals in aqueous HNO₃ and its main photo-induced byproduct HNO₂. We also found that the photo-generated OH radicals in aqueous HNO₃, HNO₂ and H₂O₂ have the remarkable capability to react with the dissolved HCl, Cl^- and Br^- to form halogen radicals. In addition, we observed that the H₂O₂ can react with dissolved HCl and Br^- in darkness to form and release Cl₂ and Br₂ gases, which could further be photolyzed into reactive halogen radicals whenever sunlight is available. All these findings suggest that, except for the well-known heterogeneous reactions, photochemical reactions involving the aqueous HNO₃ and H₂O₂ on and within PSCs surface might constitute another important halogen activation pathway for ozone destruction. This study may shed deeper insights into the mechanism of halogen radicals resulting in ozone depletion in polar stratosphere.

The HOX⋯SO2 (X=F, Cl, Br, I) Binary Complexes: Implications for Atmospheric Chemistry

Sulfur dioxide and hypohalous acids (HOX, X=F, Cl, Br, I) are ubiquitous molecules in the atmosphere that are central to important processes like seasonal ozone depletion, acid rain, and cloud nucleation. We present the first theoretical examination of the HOX⋯SO₂ binary complexes and the associated trends due to halogen substitution. Reliable geometries were optimized at the CCSD(T)/aug-cc-pV(T+d)Z level of theory for HOF and HOCl complexes. The HOBr and HOI complexes were optimized at the CCSD(T)/aug-cc-pV(D+d)Z level of theory with the exception of the Br and I atoms which were modeled with an aug-cc-pwCVDZ-PP pseudopotential. 27 HOX⋯SO₂ complexes were characterized and the focal point method was employed to produce CCSDT(Q)/CBS interaction energies. Natural Bond Orbital analysis and Symmetry Adapted Perturbation Theory were used to classify the nature of each principle interaction. The interaction energies of all HOX⋯SO₂ complexes in this study ranged from 1.35 to 3.81 kcal mol^-1 . The single-interaction hydrogen bonded complexes spanned a range of 2.62 to 3.07 kcal mol^-1 , while the single-interaction halogen bonded complexes were far more sensitive to halogen substitution ranging from 1.35 to 3.06 kcal mol^-1 , indicating that the two types of interactions are extremely competitive for heavier halogens. Our results provide insight into the interactions between HOX and SO₂ which may guide further research of related systems.

Molecular Interaction and Orientation of HOCl on Aqueous and Ice Surfaces

The interaction and orientation of hypochlorous acid (HOCl) on the ice surface has been of great interest as it has important implications to ozone depletion. As HOCl interacts with the ice surface, previous classical molecular dynamics simulations suggest its OH moiety orients to the outside of the ice surface, whereas the quantum calculations performed at 0 K indicate its Cl atom is exposed. To resolve this contradiction, herein, Born-Oppenheimer molecular dynamics simulations are adopted, and the results suggest that at ambient temperature, the interaction between HOCl with interfacial water is dominated by the robust H-bond of (HOCl)H-O(H₂O). As a result, the HOCl mainly acts as the proton donor to the water surface, which thus can participate in proton transfer reactions via the promotion of interfacial water. Moreover, the Cl atom of HOCl is found to be exposed to the outside of the water surface. Therefore, during the heterogeneous reactions of HOCl on the water surface, the Cl atom becomes the reactive site and is easily attacked by other species.

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid-base chemistry depending on temperature in small water assemblies.

Ion Pairing in HCl-Water Clusters: From Electronic Structure Investigations to Multiconfigurational Force-Field Development

In the bulk, condensed-phase HCl exists as a dissociated Cl^- ion and a proton that is delocalized over solvating water molecules. However, in the gas phase, HCl is covalent, and even on the introduction of hydrating water molecules, the HCl covalent state dominates small clusters and is relevant at larger clusters including 21 water molecules. Electronic structure calculations (at the MP2 level) and ab initio metadynamics simulations (at the DFT level) have been carried out on HCl-(H₂O)_n clusters with n = 2-22 to investigate distinct solvation environments in clusters from covalent HCl structure, to contact ion pairs and solvent-separated ion pairs. The data were further used to train and validate a multiconfigurational force-field for HCl-water clusters that incorporates covalent HCl states into the MS-EVB3.2 formalism. Additionally, the many-body interaction of the Cl^- ion with water and the excess proton was modeled by the introduction of two geometric three-body terms that incorporates the dominant many-body interaction in an efficient noniterative manner.

100 Years of Progress in Gas-Phase Atmospheric Chemistry Research

Remarkable progress has occurred over the last 100 years in our understanding of atmospheric chemical composition, stratospheric and tropospheric chemistry, urban air pollution, acid rain, and the formation of airborne particles from gas-phase chemistry. Much of this progress was associated with the developing understanding of the formation and role of ozone and of the oxides of nitrogen, NO and NO2, in the stratosphere and troposphere. The chemistry of the stratosphere, emerging from the pioneering work of Chapman in 1931, was followed by the discovery of catalytic ozone cycles, ozone destruction by chlorofluorocarbons, and the polar ozone holes, work honored by the 1995 Nobel Prize in Chemistry awarded to Crutzen, Rowland, and Molina. Foundations for the modern understanding of tropospheric chemistry were laid in the 1950s and 1960s, stimulated by the eye-stinging smog in Los Angeles. The importance of the hydroxyl (OH) radical and its relationship to the oxides of nitrogen (NO and NO2) emerged. The chemical processes leading to acid rain were elucidated. The atmosphere contains an immense number of gas-phase organic compounds, a result of emissions from plants and animals, natural and anthropogenic combustion processes, emissions from oceans, and from the atmospheric oxidation of organics emitted into the atmosphere. Organic atmospheric particulate matter arises largely as gas-phase organic compounds undergo oxidation to yield low-volatility products that condense into the particle phase. A hundred years ago, quantitative theories of chemical reaction rates were nonexistent. Today, comprehensive computer codes are available for performing detailed calculations of chemical reaction rates and mechanisms for atmospheric reactions. Understanding the future role of atmospheric chemistry in climate change and, in turn, the impact of climate change on atmospheric chemistry, will be critical to developing effective policies to protect the planet.

Effect of Hcl Addition and Temperature on the Heterogeneous Chemistryand Photochemistry of Clono2 Adsorbed on ice Crystals

The heterogeneous chemistry and photochemistry of chlorine nitrate adsorbed on HCl-doped crystals was studied mainly at 181 K and at 190 K. Under our experimental conditions the main gaseous products found at 181 K were Cl2O and Cl2, while at 190 K, Cl2O was mainly observed. At both temperatures a net enhancement in the rate of gaseous product formation was observed when light of wavelength longer than 350 nm entered the reactor. The obtained results in conditions where the ClONO2 concentration is higher than HCl concentration, and both are high enough to maintain saturation of the surface, showed that the rate of gaseous products formation at 181 K depends on the reactive species adsorbed while at 190 K, it is independence of them. The implications of these findings to the polar stratospheric chemistry are briefly discussed.

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.

In Situ Raman Spectroscopic Study of Gypsum (CaSO4·2H2O) and Epsomite (MgSO4·7H2O) Dehydration Utilizing an Ultrasonic Levitator

We present an original apparatus combining an acoustic levitator and a pressure-compatible process chamber. To characterize in situ the chemical and physical modifications of a levitated, single particle while heated to well-defined temperatures using a carbon dioxide laser, the chamber is interfaced to a Raman spectroscopic probe. As a proof-of-concept study, by gradually increasing the heating temperature, we observed the variations in the Raman spectra as 150 μg of crystals of gypsum and epsomite were dehydrated in anhydrous nitrogen gas. We display spectra showing the decreasing intensities of the ν1 symmetric and ν3 asymmetric stretching modes of water with time and the simultaneous shift of the ν1(SO4(2-)) symmetric stretch mode to higher wavenumbers. Our results demonstrate that the new apparatus is well suited to study the dehydration of levitated species such as minerals and offers potential advantages compared with previous experiments on bulk samples.

Organic constituents on the surfaces of aerosol particles from Southern Finland, Amazonia, and California studied by vibrational sum frequency generation

This article summarizes and compares the analysis of the surfaces of natural aerosol particles from three different forest environments by vibrational sum frequency generation. The experiments were carried out directly on filter and impactor substrates, without the need for sample preconcentration, manipulation, or destruction. We discuss the important first steps leading to secondary organic aerosol (SOA) particle nucleation and growth from terpene oxidation by showing that, as viewed by coherent vibrational spectroscopy, the chemical composition of the surface region of aerosol particles having sizes of 1 μm and lower appears to be close to size-invariant. We also discuss the concept of molecular chirality as a chemical marker that could be useful for quantifying how chemical constituents in the SOA gas phase and the SOA particle phase are related in time. Finally, we describe how the combination of multiple disciplines, such as aerosol science, advanced vibrational spectroscopy, meteorology, and chemistry can be highly informative when studying particles collected during atmospheric chemistry field campaigns, such as those carried out during HUMPPA-COPEC-2010, AMAZE-08, or BEARPEX-2009, and when they are compared to results from synthetic model systems such as particles from the Harvard Environmental Chamber (HEC). Discussions regarding the future of SOA chemical analysis approaches are given in the context of providing a path toward detailed spectroscopic assignments of SOA particle precursors and constituents and to fast-forward, in terms of mechanistic studies, through the SOA particle formation process.

Free Radicals Within the Antarctic Vortex: The Role of CFCs in Antarctic Ozone Loss

How strong is the case linking global release of chlorofluorocarbons to episodic disappearance of ozone from the Antarctic stratosphere each austral spring? Three lines of evidence defining a link are (i) observed containment in the vortex of ClO concentrations two orders of magnitude greater than normal levels; (ii) in situ observations obtained during ten high-altitude aircraft flights into the vortex as the ozone hole was forming that show a decrease in ozone concentrations as ClO concentrations increased; and (iii) a comparison between observed ozone loss rates and those predicted with the use of absolute concentrations of ClO and BrO, the rate-limiting radicals in an array of proposed catalytic cycles. Recent advances in our understanding of the kinetics, photochemistry, and structural details of key intermediates in these catalytic cycles as well as an improved absolute calibration for ClO and BrO concentrations at the temperatures and pressures encountered in the lower antarctic stratosphere have been essential for defining the link.

Direct ab initio dynamics calculations for rates and the kinetic isotope effects of multiproton transfer in ClONO2 + HCl --> HNO3 + Cl2 reactions with water clusters: breakdown of the rule of the geometric mean

We performed high-level quantum mechanical calculations and direct ab initio reaction dynamics calculations for multiple proton transfers in ClONO(2)+HCl-->HNO(3)+Cl(2) with water clusters containing one to two water molecules, which can be used as a model of the reactions occurring on ice surface in stratospheric clouds. The energy barriers of these reactions depend on the number of water molecules involved. Two and three protons in these reactions with one and two water molecules, respectively, were transferred concertedly and asynchronously. The potential energy barrier at the MP2/6-311++(3df,3pd)//MP2/6-31G(d,p) level was 4.8 kcal/mol for the triple proton transfer involving two water molecules with a rate constant of 1.6x10(3) s(-1) at 197 K. The potential energy curve near the saddle points was very flat and the tunneling effect on the proton transfer was negligible. The primary HH/DH kinetic isotope effect for the double proton transfer involving one water molecule was lower than unity due to the enhanced force constant at the transition state. The rule of the geometric mean for the concerted proton transfer does not hold in these reactions because the zero-point energy changes of each proton in flight at the transition state are not the same in the highly asynchronous processes.

Interaction in the ternary complexes of HCl-methanol-X, X = H2O or NH3: Ab initio calculations and on-the-fly molecular dynamics

Dynamics, structures, energetics, and vibrational spectra of the ternary complexes of hydrogen chloride with either methanol and water or methanol and ammonia were investigated by on-the-fly molecular dynamics and ab initio and density functional theory (DFT) with aug-cc-pvDZ basis sets. Addition of CH3OH to the HCl-NH3 system catalyzes the proton transfer from HCl to NH3. However, the dynamics of the system show that the proton is not localized on NH3; rather, it is shared between N and Cl.

The spectroscopy and chemical dynamics of microparticles explored using an ultrasonic trap

Microsized particles play an important role in many diverse areas of science and technology, for example, surface reactions of micron-sized particles play a key role in astrochemistry, plasma reactors and atmospheric chemistry. To date much of our knowledge of such surface chemistry is derived from 'traditional' surface science-based research. However, the large surface area and morphology of surface material commonly used in such surface science techniques may not necessarily mimic that on the surface of micron/nano scale particles. Hence, a new generation of experiments in which the spectroscopy (e.g., albedo) and chemical reactivity of micron-sized particles can be studied directly must be developed. One, as yet underexploited, non-invasive technique is the use of ultrasonic levitation. In this article, we describe the operation of an 'ultrasonic trap' to store and study the physical and chemical properties of microparticles.

Quantum state-resolved CO2 collisions at the gas-liquid interface: surface temperature-dependent scattering dynamics

Energy transfer dynamics at the gas-liquid interface are investigated as a function of surface temperature both by experimental studies of CO2 + perfluorinated polyether (PFPE) and by molecular dynamics simulations of CO2 + fluorinated self-assembled monolayers (F-SAMs). Using a normal incident molecular beam, the experimental studies probe scattered CO2 internal-state and translational distributions with high resolution infrared spectroscopy. At low incident energies [Einc = 1.6(1) kcal/mol], CO2 J-state populations and transverse Doppler velocity distributions are characteristic of the surface temperature (Trot approximately Ttrans approximately TS) over the range from 232 to 323 K. In contrast, the rotational and translational distributions at high incident energies [Einc = 10.6(8) kcal/mol] show evidence for both trapping-desorption (TD) and impulsive scattering (IS) events. Specifically, the populations are surprisingly well-characterized by a sum of Boltzmann distributions where the two components include one (TD) that equilibrates with the surface (TTD approximately TS) and a second (IS) that is much hotter than the surface temperature (TIS > TS). Support for the superthermal, yet Boltzmann, nature of the IS channel is provided by molecular dynamics (MD) simulations of CO2 + F-SAMs [Einc = 10.6 kcal/mol], which reveal two-temperature distributions, sticking probabilities, and angular distributions in near quantitative agreement with the experimental PFPE results. Finally, experiments as a function of surface temperature reveal an increase in both sticking probability and rotational/translational temperature of the IS component. Such a trend is consistent with increased surface roughness at higher surface temperature, which increases the overall probability of trapping, yet preferentially leads to impulsive scattering of more highly internally excited CO2 from the surface.