
==== Front
Proc Natl Acad Sci U S A
Proc Natl Acad Sci U S A
PNAS
Proceedings of the National Academy of Sciences of the United States of America
0027-8424
1091-6490
National Academy of Sciences

38285944
202314819
10.1073/pnas.2314819121
research-articleResearch ArticlechemChemistry410
Physical Sciences
Chemistry
Probing the dynamics and bottleneck of the key atmospheric SO2 oxidation reaction by the hydroxyl radical
Yuan Dao-Fu a b 1 https://orcid.org/0000-0001-8461-6889

Liu Yang c d 1 https://orcid.org/0000-0002-5939-0846

Trabelsi Tarek e f 1 https://orcid.org/0000-0001-6258-7191

Zhang Yue-Rou b https://orcid.org/0000-0002-4405-2967

Li Jun jli15@cqu.edu.cn
c 2 https://orcid.org/0000-0003-2392-8322

Francisco Joseph S. frjoseph@sas.upenn.edu
e f 2 https://orcid.org/0000-0002-5461-1486

Guo Hua d https://orcid.org/0000-0001-9901-053X

Wang Lai-Sheng Lai-Sheng_Wang@brown.edu
b 2 https://orcid.org/0000-0003-1816-5738

aHefei National Research Center for Physical Science at Microscale, University of Science and Technology of China, Hefei 230026, China
bDepartment of Chemistry, Brown University, Providence, RI 02912
cSchool of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China
dDepartment of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, NM 87131
eDepartment of Earth and Environmental Sciences, University of Pennsylvania, Philadelphia, PA 19104
fDepartment of Chemistry, University of Pennsylvania, Philadelphia, PA 19104
2To whom correspondence may be addressed. Email: jli15@cqu.edu.cn, frjoseph@sas.upenn.edu, or Lai-Sheng_Wang@brown.edu.
Contributed by Joseph S. Francisco; received August 29, 2023; accepted December 29, 2023; reviewed by Martin Head-Gordon and Mark A. Johnson

1D.-F.Y., Y.L., and T.T. contributed equally to this work.

29 1 2024
6 2 2024
29 7 2024
121 6 e231481912129 8 2023
29 12 2023
Copyright © 2024 the Author(s). Published by PNAS.
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ This article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).

Significance

The oxidation of SO2 (Sulfur dioxide) in the atmosphere by the hydroxyl radical to form the HOSO2 radical intermediate is the first and a rate-limiting step to the formation of acid rain and sulfate aerosols. Despite extensive kinetic studies, the underlying dynamics of this key reaction step is poorly understood. The current study provides a dynamic insight into the HO + SO2 → HOSO2 reaction using a combination of a state-of-the-art experimental technique and theoretical calculations. The implications of this work extend to atmospheric chemistry, providing fresh insights into the mechanisms underlying acid rain formation. The current study also makes significant advances to probe the dynamics of complex bimolecular chemical reactions.

SO2 (Sulfur dioxide) is the major precursor to the production of sulfuric acid (H2SO4), contributing to acid rain and atmospheric aerosols. Sulfuric acid formed from SO2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth’s atmosphere to increase the planetary albedo to counteract global warming. SO2 oxidation in the atmosphere by the hydroxyl radical HO to form HOSO2 is a key rate-limiting step in the mechanism for forming acid rain. However, the dynamics of the HO + SO2 → HOSO2 reaction and its slow rate in the atmosphere are poorly understood to date. Herein, we use photoelectron spectroscopy of cryogenically cooled HOSO2– anion to access the neutral HOSO2 radical near the transition state of the HO + SO2 reaction. Spectroscopic and dynamic calculations are conducted on the first ab initio-based full-dimensional potential energy surface to interpret the photoelectron spectra of HOSO2– and to probe the dynamics of the HO + SO2 reaction. In addition to the finding of a unique pre-reaction complex (HO⋯SO2) directly connected to the transition state, dynamic calculations reveal that the accessible phase space for the HO + SO2 → HOSO2 reaction is extremely narrow, forming a key reaction bottleneck and slowing the reaction rate in the atmosphere, despite the low reaction barrier. This study underlines the importance of understanding the full multidimensional potential energy surface to elucidate the dynamics of complex bimolecular reactions involving polyatomic reactants.

photoelectron
potential energy surface
polyatomic reactants
bimolecular reaction
spectroscopy
DOD | USAF | AMC | Air Force Office of Scientific Research (AFOSR) 100000181 FA9550-22-1-0350 Hua Guo U.S. Department of Energy (DOE) 100000015 DE-SC0018679 Lai-Sheng Wang
==== Body
pmcSulfur dioxide (SO2) is a major atmospheric pollutant. The main sources of SO2 emissions are the burning of fossil fuels and the smelting of sulfur-rich ores. SO2 is accepted as a major precursor to the production of sulfuric acid (H2SO4), contributing to acid rain (1–3) and atmospheric aerosols (4–7). Consequently, atmospheric SO2 significantly impacts the environment, human health, and climate change. Furthermore, atmospheric sulfuric acid formed from SO2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth’s atmosphere to increase the planetary albedo to counteract global warming (8–11).

The atmospheric conversion of SO2 to sulfuric acid starts via SO2 oxidation by the hydroxyl radical HO (12–16) or the Criegee intermediate in the presence of biogenic volatile organic compounds (17). The atmospheric formation of sulfuric acid via SO2 oxidation by the ubiquitous hydroxyl radical occurs through three steps (12–16, 18, 19). The first step (R1) involves gas-phase oxidation of SO2 by HO to form the HOSO2 radical intermediate, which is further oxidized by O2 to form SO3 and HO2 (R2) (12, 17). The final step (R3) corresponds to the hydration of SO3 to produce sulfuric acid, which is one of the major components of acid rain:

[R1] HO + SO2 → HOSO2,

[R2] HOSO2 + O2 → HO2 + SO3,

[R3] SO3 + nH2O → H2SO4 + (n-1)H2O.

The first step to form the key HOSO2 radical intermediate (R1) has been extensively studied both experimentally (13, 20–25) and theoretically (26–34). The existence of the HOSO2 radical was first detected by mass spectrometry (35). The radical was trapped in low-temperature matrices and studied by infrared spectroscopy (15, 16, 36). R1 and R3 are significantly slower than R2 and the production of H2SO4 requires a high concentration of water in the atmosphere (37). The reaction rate for the first step (R1) was comprehensively studied in the 1980s (21–23), and the rate coefficients were determined to be a function of pressure and temperature (20). Given the reaction rate for R1 and the atmospheric HO radical and SO2 abundance profile, the atmospheric lifetime for SO2 was estimated to be 13 d (38, 39). Recently, highly accurate ab initio theories have been used to construct the HO + SO2 reaction pathways to predict the temperature- and pressure-dependent rate coefficients (26, 27).

Although the reaction rates and kinetics for reaction R1 have been investigated extensively, the dynamics of the HO + SO2 → HOSO2 reaction remain poorly understood. For example, despite the fact that the potential barrier of this reaction is extremely low, the rate of the reaction is rather slow, leading to the relatively long lifetime of SO2 in the atmosphere (38, 39). A key issue in chemical reaction dynamics is to understand how the transition state controls the reactivity (40). As shown in Fig. 1, the reaction path for R1 features several shallow pre-reaction wells (RC1, RC2, and RC3) and an extremely low barrier (TS1) (0.77 kJ/mol, zero-point corrected) for the formation of the HOSO2 product (33, 34). A natural question is, how is this low barrier related to a slow reaction rate?

Fig. 1. Energetics for the photodetachment of HOSO2− and reaction pathway for the HO + SO2 system. The energies in kJ/mol are relative to the reactant asymptote HO + SO2 at the UCCSD(T)-F12a/aug-cc-pV(T+d)Z level and are ZPE corrected values. RC: reaction complex. TS: transition state.

The key to understanding this question is the dynamics near the transition state. Experimentally, photoelectron spectroscopy of suitable precursor anions is a promising technique for probing the spectroscopy and dynamics of the transition state region (41, 42), especially when aided by reaction dynamics analyses on accurate potential energy surfaces. This approach has proven to be powerful for probing the dynamics of elementary chemical reactions, primarily involving hydrogen transfer reactions or reactions that can be effectively described by a one-dimensional reaction coordinate (43–52). However, R1 corresponds to an association reaction involving heavy reactants and potentially multi-dimensional reaction coordinates. To further advance our understanding of fundamental polyatomic reaction dynamics, such as R1, it is essential to go beyond one-dimensional systems. In particular, it is critical to understand the multidimensional potential energy surface of the HO + SO2 reaction to gain dynamic insight and unravel the mechanisms of the slow reaction rate.

Herein, we report a photoelectron spectroscopy study of cryogenically cooled HOSO2– anion and ab initio calculations of the structures and multidimensional potential energy surfaces of the HOSO2– anion and the HOSO2 radical. The experimental and theoretical investigations uncover the importance of the shape of the multidimensional potential energy surface near the transition state in controlling the dynamics of the HO + SO2 reaction. This finding results from the analyses of the photoelectron spectra of the cold HOSO2– anion and the quantum mechanical and quasi-classical trajectory (QCT) calculations on a full-dimensional potential energy surface that is ab initio based.

Results

Photoelectron Spectroscopy.

The experiment is conducted using an electrospray ionization (ESI) photoelectron imaging apparatus (53). Further details about the experimental methods and procedures are given in SI Appendix and the Materials and Methods section. The photoelectron spectra of cryogenically-cooled HOSO2– anion obtained at photon energies ranging from 3.5399 to 4.1993 eV (341.55 to 405.17 kJ/mol) (Fig. 2A) exhibit well-resolved peaks (a to z) corresponding to detachment transitions from the anion ground state to low-lying vibrational levels of the neutral HOSO2 radical. There is another nearly isoenergetic isomer (54, 55), HSO3–, in which the H atom is directly bonded to the S atom (SI Appendix, Figs. S1 and S2). In the current experiment, the contribution from the HSO3– isomer can be excluded because of its high detachment threshold of 460.14 kJ/mol, according to our ab initio calculations (SI Appendix, Fig. S1). The S atom in HSO3– is in its highest oxidation state of +VI, which is similar to that in HOSO3–; HOSO3– has a known detachment threshold of 4.75 eV (458 kJ/mol) (56). The positions of peaks a-z are given in SI Appendix, Table S1. The observed peak a with the lowest binding energy of 3.367 eV (324.9 kJ/mol) defines the detachment threshold of HOSO2– or the electron affinity of HOSO2, which is in good agreement with the theoretical value of 320.04 kJ/mol (SI Appendix, Table S2). Peaks b–z correspond to excited vibrational levels of HOSO2.

Fig. 2. (A) Experimental photoelectron spectra of HOSO2– at 10 photon energies ranging from 3.5399 to 4.1993 eV (341.55 to 405.17 kJ/mol). The observed features are labeled as a–z. (B) Experimental photoelectron spectrum of HOSO2– at 3.6042 eV (347.75 kJ/mol) (Expt.), compared with the calculated Franck–Condon factors (FCFs). (C) Photoelectron spectra of HOSO2– at 4.661 eV (449.7 kJ/mol) (red curve) compared with the calculated Franck–Condon (blue curve) and Wigner (black curve) profiles. The black arrow indicates the dissociation threshold energy of the HOSO2– anion (to HO + SO2 + e–). TS denotes features near the transition state of the HO + SO2 reaction.

The low-energy wing of the photoelectron spectra is dominated by low-lying vibrational levels of HOSO2. Our computed anharmonic frequencies for HOSO2 are summarized in SI Appendix, Table S3 relative to previous theoretical values and experimental frequencies measured in low-temperature matrices (15, 16, 57). The atomic displacement vectors for the vibrational modes of HOSO2 and HOSO2− are presented in SI Appendix, Fig. S3. The simulated spectrum for the first few vibrational levels of HOSO2 is compared with the measured photoelectron spectrum at 3.6042 eV (347.75 kJ/mol) in Fig. 2B. The calculated spectrum is obtained by a Franck–Condon model using the normal modes of HOSO2− and HOSO2 on the anion and neutral potential energy surfaces with the experimentally determined vibrational frequencies. The active modes, including the HO pseudo-rotation (ν1, 265 cm−1), SO3 umbrella (ν4, 515 cm−1), and HO-SO2 stretch (ν5, 774 cm−1), are consistent with the differences between the neutral and anion geometries (SI Appendix, Fig. S2).

At high binding energies, the photoelectron spectrum approaches the region near TS1. The reactants HO + SO2 are 430.71 kJ/mol (4.4640 eV) above the ground state of HOSO2– (Fig. 1). To access the transition state region of the HO + SO2 reaction, we obtain the photoelectron spectrum of HOSO2– at 4.661 eV (449.7 kJ/mol) (Fig. 2C), revealing the full Franck–Condon region for the detachment transition from HOSO2– to the HOSO2 neutral final state. The vertical detachment energy of HOSO2– is determined to be 3.90 eV (376 kJ/mol) from the maximum point of the 4.661 eV (449.7 kJ/mol) spectrum (Fig. 2C), which is in good agreement with the theoretical value of 384.40 kJ/mol (3.984 eV) (SI Appendix, Table S2). By using the full potential energy surface, built from machine learning [PIP-NN (permutational invariant polynomial-neural network)] with high-level ab initio data, we simulate the photoelectron spectrum with two methods (see SI Appendix, Theoretical Methods for further details). The overall spectral profile at 4.661 eV (449.7 kJ/mol) is not well reproduced by the Franck–Condon model on the high-energy side presumably due to anharmonic effects; however, it can be reproduced using Wigner phase space sampling (Fig. 2C) (58). The Wigner method takes into account anharmonic effects in the final states and is more suitable for electronic transitions involving large geometry changes between the final state and the initial state (see SI Appendix, section II.D. for details). The latter spectrum is slightly broader than the experimental spectrum, mainly due to anharmonic effects.

Potential Energy Surface Calculations.

Fig. 1 shows that a pre-reaction complex is formed before the submerged transition state (TS1) along the HO + SO2 → HOSO2 reaction path. Previous studies (27, 34) have shown that the pre-reaction complex is a hydrogen-bonded complex between the H atom of HO and one of the O atoms of SO2 (RC2, Fig. 1). RC2 in this work is at an energy of −7.72 kJ/mol with zero-point energy (ZPE) correction, which is consistent with the previous value of −7.2 kJ/mol (27). However, the intrinsic reaction coordinate (IRC) analysis on the current potential energy surface, built from the PIP-NN machine-learning method with high-level ab initio data (see SI Appendix, Fig. S5 and Theoretical Methods for further details), reveals a unique pre-reaction complex RC1 with a ZPE-corrected energy of −8.88 kJ/mol, in which the O atom of HO is pointed toward the S atom of SO2. Importantly, the structure of RC1 is close to TS1 (Fig. 1), and it is easier to reach TS1 from RC1 than from RC2 or RC3. A cyclic complex (RC3, Fig. 1) is found at a ZPE-corrected energy of −8.10 kJ/mol. The optimized geometry of TS1 on the potential energy surface is similar to the recent high-level calculation at the CCSD(T)-F12a/cc-pVTZ-F12 level (26). The ZPE-corrected energy of TS1 is 0.77 kJ/mol, which is consistent with the previous values of 1.2 kJ/mol and 1.1 kJ/mol at CCSD(T)/CBS (27) and CCSDT(Q)/CBS (26), respectively. The fitting error (Efit − Etarget) of the potential energy surface as a function of the ab initio energy is given in SI Appendix, Fig. S4. The potential energies along the minimum energy paths of TS1, TS2, TS2′, and TS3 (SI Appendix, Fig. S1) on the potential energy surface are depicted in SI Appendix, Fig. S5. There is good agreement between the potential energy surface and the direct ab initio calculation at the UCCSD(T)-F12a/aug-cc-pV(T+d)Z level. Importantly, the intrinsic IRC analysis presented in SI Appendix, Fig. S5, shows that TS1 connects with the pre-reaction complex RC1 rather than the previously reported hydrogen-bonded RC2 (27, 34). Our calculated dissociation energy from HOSO2 to OH + SO2 is 111.83 kJ/mol (SI Appendix, Table S4), which agrees well with the previous experimental reaction enthalpy of −113.3 ± 6 kJ/mol (20) and the benchmark calculation at the CCSDT(Q)/CBS limit (−111.54 kJ/mol) (26) for the HO + SO2 → HOSO2 reaction.

Interestingly, the photoelectron spectrum at 4.661 eV (449.7 kJ/mol) exhibits reproducible oscillatory structures in the high binding energy side at approximately 4.5 eV (434 kJ/mol) (labeled as TS in Fig. 2C). These oscillatory structures occur in the same energy range predicted for the transition state of the HO + SO2 reaction (calculated to be 431.48 kJ/mol or 4.472 eV above the HOSO2– ground state, Fig. 1), strongly suggesting the existence of resonances. These resonances imply the existence of a narrow phase space bottleneck for the HO + SO2 reaction. However, a definitive identification of these resonances requires quantum dynamic calculations beyond the state-of-the-art computational methods. Nevertheless, the full-dimensional global potential energy surface obtained from the PIP-NN method can provide critical insight into the dynamics of this reaction.

Discussion

The contour plot of the potential energy surface along the reaction pathway of R1 is shown in Fig. 3, displaying the positions and energies of RC1, TS1, and HOSO2. In addition, the conversion between HOSO2 and its mirror image HOSO2′ via TS2 or TS2′ is given. Despite the extremely low-energy barrier, the entrance channel for R1 is very tight, forming the key bottleneck responsible for its slow reaction rate. The tight entrance channel is clearly illustrated in Fig. 4A, where the potential contour is displayed as functions of the angle θ and the distance ROS with other coordinates being relaxed (the definitions of θ and ROS are shown in the Inset). The reaction proceeds through a narrow entrance with a small range for the acceptance angle θ. The tight reactive entrance channel is revealed by the potential contour along the ROS and other angles, as depicted in SI Appendix, Fig. S6, resulting in the extremely low reactive encounters as shown in the trajectory calculations (vide infra).

Fig. 3. Contour plot for the potential energy surfaces (energy in kJ/mol) along the HO + SO → HOSO2 reaction pathway as functions of θ and ROS, with the other coordinates being relaxed. ROS is the distance between O2 and S3, and θ is the dihedral angle of H1-O2-S3-O4.

Fig. 4. (A) Potential energy contours (in kJ/mol) for the HO + SO2 system as functions of θ and ROS, with all other coordinates being relaxed. (B) Semilogarithmic plot of QCT survival probabilities as a function of time and the least-square fits at 298 K.

By using the PIP-NN potential energy surface, we further simulate the bimolecular reaction rate for R1 at room temperature (within the zero-pressure limit) using a QCT method. As shown in SI Appendix, Fig. S1, the formation of HSO3 is not possible under atmospheric conditions due to the large barriers (TS4 and TS5); nevertheless, the HOSO2 intermediate can be readily formed. The calculated formation rate for HOSO2 is 1.99 × 10−12 cm3 mole−1 s−1 using 9 × 105 trajectories, which is in reasonable agreement with the experimental values (24, 59) in the range of 7.10 × 10−13 to 3.61 × 10−12 cm3 mole−1 s−1. The trajectory calculations reveal that approximately 0.24% of the trajectories reach HOSO2, due to the very narrow bottleneck for R1, illustrating how the shape of the entrance channel restricts the number of reactive encounters. Furthermore, our QCT calculations suggest that nascent HOSO2 can survive a sufficiently long lifetime to experience collisional stabilization. As shown in Fig. 4B, the lifetime distribution has a bi-exponential behavior. Most of the trajectories have a relatively short lifetime (1.4 ps), but many other trajectories survive on the order of 57.0 ps. These long-lived trajectories are a direct result of the reaction bottleneck, which is conducive to resonances implicated in the photoelectron spectrum at 4.661 eV (Fig. 2C).

There is limited understanding of the dynamics and causes of the slow rate of SO2 oxidation by HO in the atmosphere. By employing photoelectron spectroscopy of cryogenically cooled HOSO2– anion, we have successfully accessed the neutral HOSO2 radical near the transition state of the HO + SO2 reaction. By comparing these observations with the first ab initio-based full-dimensional potential energy surface, we have made a significant breakthrough in understanding the dynamics of the HO + SO2 reaction. We find that the accessible phase space for the HO + SO2 → HOSO2 reaction is exceptionally narrow, resulting in a major bottleneck and contributing to the sluggish reaction rate in the atmosphere, despite the low energy barrier. Our findings demonstrate that the rate of the HO + SO2 reaction is influenced by the multidimensional reaction dynamics rather than by only the reaction energetics.

The implications of our study extend to atmospheric chemistry, providing fresh insights into the mechanisms underlying acid rain formation. The findings of this work on elucidating the dynamics of the HO + SO2 reaction underscore the crucial role of a full-dimensional potential energy surface in unravelling the intricate details of bimolecular reactions involving polyatomic reactants. It is highly probable that such dynamics-driven phenomena are ubiquitous in complex bimolecular reactions, emphasizing the importance of explicitly considering dynamics when characterizing and interpreting complicated reaction mechanisms.

Materials and Methods

Experimental Method.

Experiments were carried out using an ESI photoelectron spectroscopy apparatus equipped with a cryogenically cooled three-dimensional Paul trap and high-resolution photoelectron imaging. The photoelectron spectra of cryogenically-cooled HOSO2– anions were taken at photon energies ranging from 3.5399 to 4.1993 eV (341.55 to 405.17 kJ/mol), and at 4.661 eV (449.7 kJ/mol). Further details are available in SI Appendix.

Theoretical Method.

The geometries, energies, and harmonic frequencies of all stationary points were determined at the explicitly correlated unrestricted coupled-cluster singles, doubles, and perturbative triples level [UCCSD(T)-F12a] with the augmented correlation-consistent polarized split-valence set aug-cc-pV(T+d)Z using the 2020 version of the MOLPRO package. The aug-cc-pV(T+d)Z basis set is the improved version of the standard aug-cc-pVTZ set and can provide more accurate results for systems containing sulfur atoms. Full-dimensional potential energy surfaces were built from the PIP-NN machine-leaning method with high-level ab initio data. The FCclasses3 code with time-independent strategies was employed to calculate the photoelectron spectrum within the normal mode approximation. The photoelectron spectrum was also calculated using the Wigner phase space sampling approach. The QCT calculations were performed to obtain the rate coefficient of the OH + SO2 reaction at room temperature. Further details are available in SI Appendix.

Supplementary Material

Appendix 01 (PDF)

Click here for additional data file.

This work was supported by the US Department of Energy, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division under Grant DE-SC0018679 (to L.-S.W.), National Natural Science Foundation of China (21973009 to J.L.) and the Air Force Office of Scientific Research (FA9550-22-1-0350 to H.G.).

Author contributions

J.S.F., H.G., and L.-S.W. designed research; D.-F.Y., Y.L., T.T., Y.-R.Z., and J.L. performed research; D.-F.Y., Y.L., T.T., Y.-R.Z., J.L., J.S.F., H.G., and L.-S.W. analyzed data; and D.-F.Y., J.L., J.S.F., H.G., and L.-S.W. wrote the paper.

Competing interests

The authors declare no competing interest.

Data, Materials, and Software Availability

All study data are included in the article and/or SI Appendix.

Supporting Information

Reviewers: M.H.-G., University of California Berkeley; and M.A.J., Yale University.
==== Refs
1 B. J. Finlayson-Pitts, J. N. Pitts, Atmospheric Chemistry: Fundamental and Experimental Techniques (John Wiley and Sons, New York, 1986).
2 J. G. Calvert , Chemical mechanisms of acid generation in the troposphere. Nature 317 , 27–35 (1985).
3 R. P. Wayne, Chemistry of Atmospheres (Oxford University Press, Oxford, UK, ed. 3, 2000).
4 R. J. Weber , Measurements of enhanced H2SO4 and 3–4 nm particles near a frontal cloud during the first aerosol characterization experiment (ACE 1). J. Geophys. Res.: Atmos. 106 , 24107–24117 (2001).
5 M. Sipilä , The role of sulfuric acid in atmospheric nucleation. Science 327 , 1243–1246 (2010).20203046
6 C. Kuang, P. H. McMurry, A. V. McCormick, F. L. Eisele, Dependence of nucleation rates on sulfuric acid vapor concentration in diverse atmospheric locations. J. Geophys. Res.: Atmos. 113 , D10209 (2008).
7 R. Zhang, A. Khalizov, L. Wang, M. Hu, W. Xu, Nucleation and growth of nanoparticles in the atmosphere. Chem. Rev. 112 , 1957–2011 (2012).22044487
8 M. Hulme, Climate change: Climate engineering through stratospheric aerosol injection. Prog. Phys. Geogr.: Earth Environ. 36 , 694–705 (2012).
9 D. G. Macmartin , The climate response to stratospheric aerosol geoengineering can be tailored using multiple injection locations. J. Geophys. Res. Atmos. 122 , 12574–12590 (2017).
10 National Research Council, Climate Intervention: Reflecting Sunlight to Cool Earth (National Academies Press, Washington, D.C., 2015). 10.17226/18988.
11 D. G. Kravitz , Comparing surface and stratospheric impacts of geoengineering with different SO2 injection strategies. J. Geophys. Res.: Atmos. 124 , 7900–7918 (2019).
12 W. R. Stockwell, J. G. Calvert, The mechanism of the HO-SO2 reaction. Atmos. Environ. 17 , 2231–2235 (1983).
13 Y. Y. Lee, W. C. Kao, Y. P. Lee, Kinetics of the reaction hydroxyl + sulfur dioxide in helium, nitrogen, and oxygen at low pressure. J. Phys. Chem. 94 , 4535–4540 (1990).
14 D. Martin, J. L. Jourdain, G. Le Bras, Discharge flow measurements of the rate constants for the reaction OH + SO2 + He and HOSO2 + O2 in relation with the atmospheric oxidation of sulfur dioxide. J. Phys. Chem. 90 , 4143–4147 (1986).
15 S. Nagase, S. Hashimoto, H. Akimoto, Bisulfite and peroxybisulfate radicals studied by ab initio calculation and matrix isolation technique. J. Phys. Chem. 92 , 641–644 (1988).
16 S. Hashimoto, G. Iioue, H. Akimoto, Infrared spectroscopic detection of the HOSO2 radical in argon matrix at 11 K. Chem. Phys. Lett. 107 , 198–202 (1984).
17 R. L. Mauldin III , A new atmospherically relevant oxidant of sulphur dioxide. Nature 488 , 193–196 (2012).22874964
18 G. D. Greenblatt, C. J. Howard, Oxygen atom exchange in the interaction of hydroxyl-(18OH) with several small molecules. J. Phys. Chem. 93 , 1035–1042 (1989).
19 D. Majumdar , Ab initio investigations on the HOSO2 + O2 → SO3 + HO2 reaction. J. Chem. Phys. 112 , 723–730 (2000).
20 M. A. Blitz, K. J. Hughes, M. J. Pilling, Determination of the high-pressure limiting rate coefficient and the enthalpy of reaction for OH + SO2. J. Phys. Chem. A 107 , 1971–1978 (2003).
21 P. H. Wine , Kinetics of the reaction OH + SO2 + M → HOSO2 + M. Temperature and pressure dependence in the fall-off region. J. Phys. Chem. 88 , 2095–2104 (1984).
22 G. Paraskevopoulos, D. L. Singleton, R. S. Irwin, Rates of OH radical reactions. The reaction OH + SO2 + N2. Chem. Phys. Lett. 100 , 83–87 (1983).
23 M. T. Leu, Rate constants for the reaction of hydroxyl with sulfur dioxide at low pressure. J. Phys. Chem. 86 , 4558–4562 (1982).
24 D. Fulle, H. F. Hamann, H. Hippler, The pressure and temperature dependence of the recombination reaction HO + SO2 + M → HOSO2 + M. Phys. Chem. Chem. Phys. 1 , 2695–2702 (1999).
25 M. A. Blitz, R. J. Salter, D. E. Heard, P. W. Seakins, An experimental study of the kinetics of OH/OD (v = 1, 2, 3) + SO2: The limiting high-pressure rate coefficients as a function of temperature. J. Phys. Chem. A 121 , 3175–3183 (2017).28363245
26 B. Long, J. L. Bao, D. G. Truhlar, Reaction of SO2 with OH in the atmosphere. Phys. Chem. Chem. Phys. 19 , 8091–8100 (2017).28265640
27 D. J. Medeiros, M. A. Blitz, P. W. Seakins, Exploring the features on the OH + SO2 potential energy surface using theory and testing its accuracy by comparison to experimental data. Phys. Chem. Chem. Phys. 20 , 8984–8990 (2018).29557461
28 S. Sitha, L. L. Jewell, S. J. Piketh, G. Fourie, A quantum chemical calculation of the potential energy surface in the formation of HOSO2 from OH + SO2. Atmos. Environ. 45 , 745–754 (2011).
29 D. R. Glowacki, S. K. Reed, M. J. Pilling, D. V. Shalashilin, E. Martínez-Núñez, Classical, quantum and statistical simulations of vibrationally excited HOSO2: IVR, dissociation, and implications for OH + SO2 kinetics at high pressures. Phys. Chem. Chem. Phys. 11 , 963–974 (2009).19177215
30 H. Somnitz, Quantum chemical and dynamical characterisation of the reaction OH + SO2 ⇌ HOSO2 over an extended range of temperature and pressure. Phys. Chem. Chem. Phys. 6 , 3844–3851 (2004).
31 P. L. Moore-Plummer, R. Wu, J. J. Flenner, Ab initio investigation of the gas phase reaction SO2 + H2O → SO2OH radical. Mol. Phys. 100 , 1847–1853 (2002).
32 F.-Y. Leung, A. J. Colussi, M. R. Hoffmann, Sulfur isotopic fractionation in the gas-phase oxidation of sulfur dioxide initiated by hydroxyl radicals. J. Phys. Chem. A 105 , 8073–8076 (2001).
33 W.-K. Li, M. L. McKee, Theoretical study of OH and H2O addition to SO2. J. Phys. Chem. A 101 , 9778–9782 (1997).
34 V. M. Miriyala, P. Bhasi, Z. P. Nhlabatsi, S. Sitha, Formation of a pre-reaction hydrogen-bonded complex and its significance in the potential energy surface of the OH + SO2 → HOSO2 reaction: A computational study. J. Theor. Comput. Chem. 16 , 1750046 (2017).
35 H. Egsgaard, L. Carlsen, H. Florencio, T. Drewello, H. Schwarz, Experimental evidence for the gaseous HSO3 radical. The key intermediate in the oxidation of SO2 in the atmosphere. Chem. Phys. Lett. 148 , 537–540 (1988).
36 Y.-P. Kuo, B.-M. Cheng, Y.-P. Lee, Production and trapping of HOSO2 from the gaseous reaction OH + SO2: The infrared absorption of HOSO2 in solid argon. Chem. Phys. Lett. 177 , 195–199 (1991).
37 M. Torrent-Sucarrat, J. S. Francisco, J. M. Anglada, Sulfuric acid as autocatalyst in the formation of sulfuric acid. J. Am. Chem. Soc. 134 , 20632–20644 (2012).23198746
38 J. R. Green , Rates of wintertime atmospheric SO2 oxidation based on aircraft observations during clear-sky conditions over the eastern United States. J. Geophys. Res.: Atmos. 124 , 6630–6649 (2019).
39 C. Lee , SO2 emissions and lifetimes: Estimates from inverse modeling using in situ and global, space-based (SCIAMACHY and OMI) observations. J. Geophys. Res. 116 , D06304 (2011).
40 J. C. Polanyi, A. H. Zewail, Direct observation of the transition state. Acc. Chem. Res. 28 , 119–132 (1995).
41 D. M. Neumark, Transition state spectroscopy of bimolecular chemical reactions. Annu. Rev. Phys. Chem. 43 , 153–176 (1992).
42 D. M. Neumark, Probing the transition state with negative ion photodetachment: Experiment and theory. Phys. Chem. Chem. Phys. 7 , 433–442 (2005).19785127
43 D. E. Manolopoulos , The transition state of the F + H2 reaction. Science 262 , 1852–1855 (1993).17829631
44 P. G. Wenthold, D. A. Hrovat, W. T. Borden, W. C. Lineberger, Transition-state spectroscopy of cyclooctatetraene. Science 272 , 1456–1459 (1996).8662467
45 D. H. Zhang, M. Yang, M. A. Collins, S.-Y. Lee, Probing the transition state via photoelectron and photodetachment spectroscopy of H3O−. Proc. Natl. Acad. Sci. U.S.A. 99 , 11579–11582 (2002).12186977
46 E. Garand, J. Zhou, D. E. Manolopoulos, M. H. Alexander, D. M. Neumark, Nonadiabatic interactions in the Cl + H2 reaction probed by ClH2− and ClD2− photoelectron imaging. Science 319 , 72–75 (2008).18174436
47 R. Otto , Imaging dynamics on the F + H2O → HF + OH potential energy surfaces from wells to barriers. Science 343 , 396–399 (2014).24407479
48 J. B. Kim , Spectroscopic observation of resonances in the F + H2 reaction. Science 349 , 510–513 (2015).26228142
49 M. L. Weichman , Feshbach resonances in the exit channel of the F + CH3OH → HF + CH3O reaction observed using transition-state spectroscopy. Nat. Chem. 9 , 950–955 (2017).28937670
50 J. A. DeVine , Encoding of vinylidene isomerization in its anion photoelectron spectrum. Science 358 , 336–339 (2017).29051373
51 Y. Benitez, T. L. Nguyen, A. J. Parsons, J. F. Stanton, R. E. Continetti, Probing the exit channel of the OH + CH3OH → H2O + CH3O reaction by photodetachment of CH3O–(H2O). J. Phys. Chem. Lett. 13 , 142–148 (2022).34962408
52 M. C. Babin , Observation of resonances in the transition state region of the F + NH3 reaction using anion photoelectron spectroscopy. Nat. Chem. 15 , 194–199 (2022).36509851
53 L. S. Wang, Electrospray photoelectron spectroscopy: From multiply-charged anions to ultracold anions. J. Chem. Phys. 143 , 040901 (2015).26233095
54 D. A. Horner, R. E. Connick, Equilibrium quotient for the isomerization of bisulfite ion from HSO3– to SO3H–. Inorg. Chem. 25 , 2414–2417 (1986).
55 R. R. Squires, Gas-phase thermochemical properties of the bicarbonate and bisulfite ions. Int. J. Mass Spectrom. Ion Processes 117 , 565–600 (1992).
56 X. B. Wang, J. B. Nicholas, L. S. Wang, Photoelectron spectroscopy and theoretical calculations of SO4- and HSO4-: Confirmation of high electron affinities of SO4 and HSO4. J. Phys. Chem. A 104 , 504–508 (2000).
57 M. Wierzejewska, A. Olbert-Majkut, Photolysis of matrix isolated HONO/SO2 system. Identification and infrared spectra of nitrososulfonic acid NO(NO)SO2 and hydroxysulfonyl HOSO2 radical. J. Phys. Chem. A 107 , 10944–10952 (2003).
58 E. J. Heller, Wigner phase space method: Analysis for semiclassical applications. J. Chem. Phys. 65 , 1289–1298 (1976).
59 R. Atkinson , Evaluated kinetic and photochemical data for atmospheric chemistry. Atmos. Chem. Phys. 4 , 1461–1738 (2004).
