Chinese Journal of Chemical Physics  2018, Vol. 31 Issue (4): 503-509

The article information

Emily L. Tao, Jessica Y. Li, Samantha Soriano, Fu-Ming Tao
陶琳, 李宇杰, 萨曼莎-索里亚诺, 陶福明
Quantum Chemical Study of Potential Energy Surface in the Formation of Atmospheric Sulfuric Acid
大气中硫酸生成反应势能面的量子化学研究
Chinese Journal of Chemical Physics, 2018, 31(4): 503-509
化学物理学报, 2018, 31(4): 503-509
http://dx.doi.org/10.1063/1674-0068/31/cjcp1805126

Article history

Received on: May 30, 2018
Accepted on: June 23, 2018
Quantum Chemical Study of Potential Energy Surface in the Formation of Atmospheric Sulfuric Acid
Emily L. Tao, Jessica Y. Li, Samantha Soriano, Fu-Ming Tao     
Dated: Received on May 30, 2018; Accepted on June 23, 2018
Department of Chemistry and Biochemistry, California State University, Fullerton, 800 N. State College Blvd. Fullerton, CA 92834, USA
Author: Fu-Ming Tao is Professor of Theoretical/Computational Chemistry at California State University Fullerton (CSUF). He entered the undergraduate class of 1978 at the University of Science and Technology of China where he received his B.S. degree in Chemistry in 1982. He also earned his M.S. degree in Physical Chemistry at Suzhou University in 1985 and Ph.D. degree in Physical Chemistry at Boston College in 1991. From 1991 to 1995, he did his post-doctoral research at Brown University, where he studied physical properties of real gases and liquids using statistical mechanics, and at Harvard University, where he studied the structure and vibrational states of van der Waals complexes using ab initio theory. He started his tenure-track position at CSUF in 1995. Since then, his research has been focused on molecular properties and reaction mechanisms of atmospheric species using ab initio and density functional methods. His main contributions include the development of bond function basis sets for accurate ab initio calculations of intermolecular potentials and the use of explicit solvent (water) models in understanding several key problems in atmospheric and environmental reactions. He has published more than 140 research papers, mostly in Journal of Chemical Physics, Journal of Physical Chemistry, and Chemical Physics Letters. He has received over $500,000 in research grants from National Science Foundation (NSF), ACS-Petroleum Research Foundation, Dreyfus Foundation, and Research Corporation. He has been honored with numerous national awards, including Henry Dreyfus Teacher/Scholar Award in 2000 and Oversea Young Scientist Award (National Science Foundation of China) in 2005.
*Author to whom correspondence should be addressed. Fu-Ming Tao, E-mail: ftao@fullerton.edu
Part of the special issue for celebration of "the 60th Anniversary of University of Science and Technology of China and the 30th Anniversary of Chinese Journal of Chemical Physics"
Abstract: A new potential energy surface (PES) for the atmospheric formation of sulfuric acid from OH+SO2 is investigated using density functional theory and high-level ab initio molecular orbital theory. A pathway focused on the new PES assumes the reaction to take place between the radical complex SO3·HO2 and H2O. The unusual stability of SO3·HO2 is the principal basis of the new pathway, which has the same final outcome as the current reaction mechanism in the literature but it avoids the production and complete release of SO3. The entire reaction pathway is composed of three consecutive elementary steps:(1) HOSO2+O2→SO3·HO2, (2) SO3·HO2+H2O→SO3·H2O·HO2, (3) SO3·H2O·HO2→H2SO4+HO2. All three steps have small energy barriers, under 10 kcal/mol, and are exothermic, and the new pathway is therefore favorable both kinetically and thermodynamically. As a key step of the reactions, step (3), HO2 serves as a bridge molecule for low-barrier hydrogen transfer in the hydrolysis of SO3. Two significant atmospheric implications are expected from the present study. First, SO3 is not released from the oxidation of SO2 by OH radical in the atmosphere. Second, the conversion of SO2 into sulfuric acid is weakly dependent on the humidity of air.
Key words: Quantum chemistry    Atmospheric chemistry    Sulfur dioxide    Sulfuric acid    Computational chemistry    Density functional theory    Ab initio methods    
Ⅰ. INTRODUCTION

Atmospheric reactions of sulfur compounds, particularly for the formation of sulfuric acid, have been a topic of considerable study in recent years [1, 2]. Sulfuric acid, a major component of acid rain, is also known to influence atmospheric nucleation processes for the formation of atmospheric aerosols [3-6]. Because virtually all sulfur in the atmosphere is either emitted as or eventually converted into sulfur dioxide (SO$_2$), many studies focus on the oxidation of SO$_2$ as an initial step for sulfuric acid formation.The general reaction process of sulfur dioxide oxidation is fairly well-established, beginning with the reaction of SO$_2$ with hydroxyl radical (OH) [1, 2, 7]. Sulfur compounds emitted to the atmosphere are oxidized to form SO$_2$, which reacts first with OH to form the HOSO$_2$ radical. The HOSO$_2$ radical then forms SO$_3$ and HO$_2$ through collisions with O$_2$. Finally SO$_3$ goes on to react with water via a water/acid catalyzed reaction to form sulfuric acid:

(1) SO$_2$ + OH + M $\rightarrow$ HOSO$_2$ + M

(2) HOSO$_2$ + O$_2$ $\rightarrow$ SO$_3$ + HO$_2$

(3) SO$_3$ + H$_2$O $\rightarrow$ H$_2$SO$_4$

The first step of this mechanism has been more or less confirmed by theoretical and experimental evidence; its product, the HOSO$_2$ radical, has been detected directly in the atmosphere [1]. Step (2) has been proposed as a dominant channel for the reaction of the HOSO$_2$ radical in the atmosphere [7]. While both the thermochemical estimation [8] and theoretical calculation [9-11] showed that the reaction was endothermic with an energy change of 8 kcal/mol or larger, a recent study revealed it to be slightly exothermic with an enthalpy of -2.3 kcal/mol [12]. Some studies propose an alternate step involving the formation of a HOSO$_4$ radical instead of SO$_3$ and HO$_2$ [13]. However, this step, while a likely atmospheric process, is of little importance in the overall reaction mechanism because the HOSO$_4$ radical may form SO$_3$ and HO$_2$ on its own through a fast hydrogen transfer step [2].

Reaction (3) also has been a subject of several theoretical studies [14, 15]. Sulfuric acid is unlikely to form directly from an elementary reaction of SO$_3$ and H$_2$O due to the high energy barrier involved (about 30 kcal/mol) [16]. Although various studies [15, 17] show that reaction of SO$_3$ in a cluster of two, three, or four waters significantly lowers the barrier, making this step appear more favorable, clusters of multiple waters are in fact highly unstable in the gas phase under normal atmospheric conditions. Therefore, it is still uncertain that such a process actually occurs in the atmosphere, and an alternate reaction pathway should be considered.

This study attempts to find an alternate pathway to better describe the sulfuric acid formation mechanism. One possible pathway for consideration is the direct reaction of a water molecule with the intermediate product SO$_3$$\cdot$HO$_2$ from reaction (2), rather than releasing the separate intermediate products SO$_3$ and HO$_2$. The reaction steps can be revised as follows.

(4) HOSO$_2$ + O$_2$ $\rightarrow$ SO$_3$$\cdot$HO$_2$

(5) SO$_3$$\cdot$HO$_2$ + H$_2$O $\rightarrow$ H$_2$SO$_4$ + HO$_2$

The SO$_3$$\cdot$HO$_2$ complex is known to be highly stable. A binding energy of 12.5 kcal/mol was found for SO$_3$$\cdot$HO$_2$ with respect to the separate SO$_3$ and HO$_2$ species [11]. Water, abundant under normal atmospheric conditions, is widely available to facilitate the hydrolysis of SO$_3$$\cdot$HO$_2$ to form sulfuric acid and HO$_2$. The final outcome of the new mechanism is identical to that of the original mechanism in the literature, except that the production and complete release of SO$_3$ is avoided in the new mechanisms. Other possible pathways could also be considered. For example, a radical complex HOSO$_2$$\cdot$H$_2$O could be formed prior to the reaction with O$_2$. The radical complex HOSO$_2$$\cdot$H$_2$O is also very stable with a binding energy of 10.6 kcal/mol [18].

In this study, the atmospheric reaction mechanism for the production of sulfuric acid from SO$_2$ is investigated with quantum chemical calculations using density functional theory and high-level ab initio molecular orbital theory. The new reaction pathway outlined above is the main focus of this study. All reactant complexes, transition states, reaction intermediates, product complexes, and separate products are characterized for their molecular and thermodynamic properties. A new reaction mechanism will be established based on the results, and the atmospheric implications will be discussed.

Ⅱ. THEORETICAL METHODS

Equilibrium geometries and energies of the reactant complexes, transition states, and reaction intermediates involved in reactions (1) to (5) were calculated using density functional theory (DFT) and ab initio molecular orbital theory. These geometries and energies represent the stationary points on a new potential energy surface in the atmospheric formation of H$_2$SO$_4$ from OH+SO$_2$. The DFT method in this study employed Becke's three-parameter functional with nonlocal correlation provided by Lee, Yang, and Parr (B3LYP) [19, 20], along with three different basis sets, 6-31+G$^*$ [21], aug-cc-pVDZ [22], and aug-cc-pVTZ [22]. Two ab initio methods for different levels of electron correlation involving valence electrons only were employed: second-order Møller-Plesset perturbation approximation (MP2) [23] and coupled-cluster method with single, double, and non-iterative triple excitations (CCSD(T)) [24-27], both using the aug-cc-pVDZ basis set. The CCSD(T) calculations were carried out using single-point energy calculations on the MP2/aug-cc-pVDZ geometries. Frozen-core approximation was used in both MP2 and CCSD(T) calculations. Level of consistency in the results from the diverse theoretical methods and basis sets may indicate the accuracy and reliability in the results of calculation. Transition state geometries were found using the synchronous transit-guided quasi-Newton (STQN) method [28, 29]. Transition states were verified by harmonic frequency calculations and by Intrinsic Reaction Coordinate (IRC) calculations [30, 31].

Relative electronic energy ($\Delta E$), enthalpy ($\Delta H$), and Gibbs free energy ($\Delta G$) were determined as the difference between the energy of the given molecular complex (or transition state) and the sum of the energies of the respective reactant monomers HOSO$_2$, H$_2$O, and O$_2$ (in the triplet ground state), that is,

$ \Delta E=E-(E_{\textrm{HOSO}_2}+E_{\textrm{H}_2\textrm{O}}+E_{\textrm{O}_2}) $ (1)
$ \Delta H=H-(H_{\textrm{HOSO}_2}+H_{\textrm{H}_2\textrm{O}}+H_{\textrm{O}_2}) $ (2)
$ \Delta G=G-(G_{\textrm{HOSO}_2}+G_{\textrm{H}_2\textrm{O}}+G_{\textrm{O}_2}) $ (3)

where $E$, $H$, and $G$ are the total electronic energy, enthalpy, and Gibbs free energy of the molecular complex (or transition state), and $E_m$, $H_m$, and $G_m$ are the corresponding energies of the reactant monomer (HOSO$_2$, H$_2$O, or O$_2$) which the molecular complex is composed of. The relative enthalpy ($\Delta H$) and Gibbs free energy ($\Delta$$G$), reported at 298.15 K and 1 atm, were calculated using harmonic vibrational frequencies. All calculations in this study were carried out on the Gaussian 09 computational chemistry program [32] running on a UNIX computer cluster.

Included in our calculations are all species involved in reactions (1), (2), and (3), particularly the reactant complex HOSO$_2$$\cdot$O$_2$, product complex SO$_3$$\cdot$HO$_2$, and the corresponding transition state TS0. Instead of the dissociation channel of SO$_3$$\cdot$HO$_2$ into SO$_3$ and the HO$_2$ radical, our focus is on the continuing reaction of SO$_3$$\cdot$HO$_2$ with water which is abundant in the atmosphere. As a new reaction pathway on the potential energy surface, a water molecule is introduced to the complex SO$_3$$\cdot$HO$_2$, resulting in the trimolecular complex SO$_3$$\cdot$HO$_2$$\cdot$H$_2$O as involved in reaction (5). It should be noted that several conformations are expected for the complex SO$_3$$\cdot$HO$_2$$\cdot$H$_2$O and the present study primarily focuses on the most stable conformation as well as the one which would directly facilitate the hydrolysis reaction to form sulfuric acid in the following step.

Ⅲ. RESULTS AND DISCUSSION

FIG. 1 shows the equilibrium geometries of all complexes and transition states at the stationary points along the new reaction pathway on the potential energy surface. It should be noted that selected bond lengths in the figure were obtained from MP2/aug-cc-pVDZ calculations. The equilibrium geometries from B3LYP calculations with the three basis sets, 6-31+G$^*$, aug-cc-pVDZ, and aug-cc-pVQZ respectively, are similar to those from MP2/aug-cc-pVDZ calculations. The reaction pathway is composed of three reaction steps: (ⅰ) the radical complex HOSO$_2$ reacts with O$_2$ to form the radical complex SO$_3$$\cdot$HO$_2$, (ⅱ) SO$_3$$\cdot$HO$_2$ combines H$_2$O to form the radical complex SO$_3$$\cdot$HO$_2$$\cdot$H$_2$O which is rearranged to SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$, and (ⅲ) SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$ reacts to form the product complex H$_2$SO$_4$$\cdot$HO$_2$. Each of the reaction steps involves a transition state. For reaction step (ⅰ), the reactant complex, transition state, and the product complex are denoted HOSO$_2$$\cdot$O$_2$, TS0, and SO$_3$$\cdot$HO$_2$, respectively. Similarly, the corresponding species for step (ⅱ) are SO$_3$$\cdot$HO$_2$$\cdot$H$_2$O, TS1, and SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$, and those for step (ⅲ) are SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$, TS2, and H$_2$SO$_4$$\cdot$HO$_2$.

FIG. 1 Optimized geometries from MP2/aug-cc-pVDZ calculations. Selected bond lengths are in Å

The total electronic energies of the species of all species at the stationary points along the reaction pathways on the potential energy surface, including those of separate reactant and product molecules, are available from the authors directly. These energies are from B3LYP/6-31+G$^*$, B3LYP/aug-cc-pVDZ, B3LYP/ aug-cc-pVTZ, MP2/aug-cc-pVDZ, and CCSD(T)/aug-cc-pVDZ calculations, respectively. Table Ⅰ gives the relative electronic energies ($\Delta E$), relative enthalpies ($\Delta$$H$), and relative Gibbs free energies ($\Delta G$) with respect to the separate reactants, HOSO$_2$, O$_2$, and H$_2$O. Note that only electronic energies ($\Delta E$) are available from CCSD(T)/aug-cc-pVDZ calculations. FIG. 2 shows a schematic of the corresponding potential energy profile for the reactions of HOSO$_2$, O$_2$, and H$_2$O to form H$_2$SO$_4$ and HO$_2$.

Table Ⅰ Relative electronic energies ($\Delta E$), enthalpies ($\Delta H$), and Gibbs free energies ($\Delta G$) (all in kcal/mol) from B3LYP, MP2, and CCSD(T) calculations with different basis sets.
FIG. 2 Schematic of the potential energy profile for the reaction HOSO$_2$+O$_2$+H$_2$O$\rightarrow$H$_2$SO$_4$+HO$_2$.

As shown in Table Ⅰ and FIG. 2, the reaction HOSO$_2$+O$_2$$\rightarrow$SO$_3$+HO$_2$ has a very small energy barrier from the reactant complex HOSO$_2$$\cdot$O$_2$: about 4.0, 0.1, and 1.4 kcal/mol at the B3LYP, MP2, and CCSD(T) levels, respectively. The corresponding activation energies are either slightly negative or near zero, which would suggest the reaction is fast and kinetically favorable. However, the separate products SO$_3$ and HO$_2$ are about 10 kcal/mol higher in energy than the reactants HOSO$_2$ and O$_2$. The result is consistent with the thermal estimation [8] and the reported theoretical calculations [9-11]. It indicates that the reaction HOSO$_2$+O$_2$$\rightarrow$SO$_3$+HO$_2$ is not thermodynamically favorable. The similar, positive values of $\Delta G$ at the various levels confirm the thermodynamic un-favorability of the reaction. The strong intermolecular energy between SO$_3$ and HO$_2$, as in the product complex SO$_3$$\cdot$HO$_2$, is mainly responsible for this result. It can thus be concluded that the reaction HOSO$_2$+O$_2$ may take place, but it only produces the SO$_3$$\cdot$HO$_2$ complex.

In the new reaction pathway, a water molecule is added to the SO$_3$$\cdot$HO$_2$ complex formed from the reaction between HOSO$_2$ and O$_2$. At least two conformations were found for the SO$_3$$\cdot$HO$_2$$\cdot$H$_2$O complex. In the most stable conformation, H$_2$O accepts a hydrogen bond from HO$_2$ while donating a hydrogen bond to an oxygen atom of SO$_3$. In the second conformation, SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$, the water is inserted between the O atom of HO$_2$ and the S atom of SO$_3$. In other words, HO$_2$ and H$_2$O exchange positions on SO$_3$. The SO$_3$$\cdot$H$_2$O unit is an electron donor-acceptor complex involved as a pre-reaction complex for reaction (3). There is a transition state, TS1, for the conversion between the two complexes, SO$_3$$\cdot$HO$_2$$\cdot$H$_2$O and SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$. The fairly low value of the relative energy at TS1, -1.6 kcal/mol at the MP2 level or -2.7 kcal/mol at CCSD(T), indicates that the conversion should be kinetically favorable. The further next step of the reaction is the hydrolysis of SO$_3$ in the complex SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$. The transition state, TS2, involves the transfer of an H atom from H$_2$O to HO$_2$ as the latter returns its own H atom to SO$_3$. The relative energy at TS2, -6.8 kcal/mol at the MP2 level or -4.7 kcal/mol at CCSD(T), is again very low, and so the reaction is also expected to be fast. The product complex H$_2$SO$_4$$\cdot$HO$_2$ is highly stable with a low relative energy of -13.4 kcal/mol at the MP2 level or -17.3 kcal/mol at CCSD(T). Finally, the H$_2$SO$_4$$\cdot$HO$_2$ complex is dissociated into the separate products H$_2$SO$_4$ and HO$_2$. The overall reaction, HOSO$_2$+O$_2$+H$_2$O$\rightarrow$SO$_3$+HO$_2$, has a value of $\Delta E$=-5.0 kcal/mol at the CCSD(T) level. The reaction is expected to be slightly endothermic or exothermic as shown by small values of $\Delta H$, -8.1 and 1.3 kcal/mol from B3LYP/aug-cc-pVTZ and MP2/aug-cc-pVDZ calculations, respectively.

The hydrolysis of SO$_3$, with the transition state TS2, is a critical step in the reaction pathway. The relative energy of the transition state is surprisingly low. The corresponding energy barrier, about 6 kcal/mol from both MP2 and CCSD(T) methods, is much smaller than the energy barriers for the hydrolysis of SO$_3$ in small water clusters SO$_3$$\cdot$(H$_2$O)$_n$, ca. 21 and 14 kcal/mol with one and two waters, respectively [14, 15]. The barrier height is comparable to the one found for the cluster of three water molecules [15]. Clearly, HO$_2$ plays a more efficient role in transferring the proton/hydrogen atom in the hydrolysis reaction than a single H$_2$O molecule or a water dimer (H$_2$O)$_2$. At the transition state TS2, the H$\cdots$O-O$\cdots$H chain appears to be equivalent to a stable molecule of hydrogen peroxide (H$_2$O$_2$), which might be responsible for the stability of TS2 and the low energy barrier for the hydrolysis of SO$_3$ in the complex SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$. Note that the progressively large unit of the hydrated proton was responsible for the step-wise decrease in the energy barrier in the clusters SO$_3$$\cdot$(H$_2$O)$_n$ [15].

The results should be reliable as shown by generally consistent values of the relative energies ($\Delta E$, $ \Delta H$, and $\Delta G$) from different methods and basis sets in Table Ⅰ. There are some noticeable fluctuations in these values at the B3LYP level among the three basis sets, 6-31+G$^*$, aug-cc-pVDZ, and aug-cc-pVTZ. It appears that the proposed new reactions are more favored by the use of the largest basis set (aug-cc-pVTZ). On the other hand, MP2 and B3LYP calculations using the same basis set (aug-cc-pVDZ) give surprisingly consistent results despite of the two very different theoretical methods. Furthermore, CCSD(T) calculations using the same basis set have nearly reproduced the results from MP2 calculations. It is expected that the energetic values for the proposed reactions would become more favorable if CCSD(T) calculations were performed using a larger basis set beyond aug-cc-pVDZ. Such calculations were attempted but were not complete because of limitations in our computer resources.

Finally, it is important to point out that the spin-contamination problem is very mild throughout the calculations and it does not severely affect the reliability of the results in the present study. Typically unrestricted DFT methods such as B3LYP are subject to less spin contamination than molecular orbital methods such as UHF and UMP2. Coupled-cluster theory is even more effective in reducing the spin-contamination problem and CCSD(T) nearly completely eliminates the two major contaminants (S+1 and S+2). This is another reason why three different methods, B3LYP, MP2, and CCSD(T), were performed in our calculations. Almost all molecular species (reactant complexes, transition states, and product complexes) in the study are open-shell with the exact expectation eigenvalue of 0.75 for the spin-squared expectation operator $\langle S^2\rangle$. The unrestricted B3LYP, MP2, and CCSD(T) procedures were used for the calculations of these open-shell species. The B3LYP calculations with the different basis sets give the $\langle S^2\rangle$ values between 0.7550 and 0.7563 before annihilation (0.7500 after annihilation). The MP2 calculations give the values between 0.7639 and 0.7816 before annihilation (0.7501-0.7507 after annihilation). As expected, no spin-contamination problems were found in the CCSD(T) calculations. It is clear that the results in the present study are not affected by the spin-contamination problem.

Ⅳ. CONCLUSION AND ATMOSPHERIC IMPLICATIONS

Density functional theory and high-level ab initio molecular orbital theory have been used to explore the new potential energy surface for the atmospheric formation of sulfuric acid from OH+SO$_2$ by simultaneous inclusion of O$_2$ and H$_2$O molecules. Instead of producing the separate products of SO$_3$ and HO$_2$, as for the reaction HOSO$_2$+O$_2$ discussed in the literature, the intermediate complex SO$_3$$\cdot$HO$_2$ directly reacts with a water molecule and produces the final products H$_2$SO$_4$ and HO$_2$. The entire reaction pathway is composed of three consecutive elementary steps: (1) HOSO$_2$+O$_2$$\rightarrow$SO$_3$$\cdot$HO$_2$, (2) SO$_3$$\cdot$HO$_2$+H$_2$O$\rightarrow$SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$, (3) SO$_3$$\cdot$H$_2$O$\cdot$HO$_2$$\rightarrow$ H$_2$SO$_4$+HO$_2$. Molecular geometries at the stationary points of PES were calculated using B3LYP and MP2 methods with several basis sets, corresponding to the reactant complex, transition state, and product complex involved in each reaction step. The results consistently show that all three elementary steps have very low activation energies and are exothermic or nearly exothermic. It implies that the reactions are favorable both kinetically and thermodynamically. The results have been further confirmed by single-point energy calculations with high-level CCSD(T) method using MP2 geometries.

The new reaction pathway for the formation of sulfuric acid from SO$_2$ is more direct and favorable than existing mechanisms in the literature which assume the release of SO$_3$ from the initial steps of reaction. No additional species are required in the new pathway to initiate or complete the reaction other than OH, O$_2$, and H$_2$O, which are abundant in the atmosphere and known to be the principal ingredients in the conversion of SO$_2$ into sulfuric acid. There are several major advantages of the new pathway over existing mechanisms in the literature. First, the hydrolysis of SO$_3$ takes place before SO$_3$ is released from the SO$_3$$\cdot$HO$_2$ complex. In other words, no release of SO$_3$ is necessary in the new reaction pathway. There is very strong intermolecular binding in the SO$_3$$\cdot$HO$_2$ complex, about 10 kcal/mol, which is thermodynamically unfavorable to dissociate into SO$_3$ and HO$_2$. Second, the radical HO$_2$ in the SO$_3$$\cdot$HO$_2$ complex assists the hydrolysis by exchanging its hydrogen between H$_2$O and an O atom of SO$_3$ and lowering the energy barrier of reaction; without HO$_2$, additional water molecules, at least two, must be added to serve the same purpose of lowering the energy barrier. Third, the hydrogen of HO$_2$ generated in the new mechanism originates from the water molecule was directly used for the hydrolysis of SO$_3$, instead of coming from HOSO$_2$ (originally the OH radical). Note that the production of HO$_2$ is important to the overall chemistry of SO$_2$ oxidation [1]. In the presence of NO, the HO$_2$ reacts to regenerate OH by the reaction HO$_2$+NO$\rightarrow$OH+NO$_2$. As a result, the oxidation of SO$_2$ is a chain reaction without a net loss of OH. The production of HO$_2$ is predicted in the new mechanism as well as in existing mechanisms in the literature.

Clearly, there are several significant atmospheric implications with the new mechanism. First, no SO$_3$ is released or can be observed in the conversion of SO$_2$ into sulfuric acid. As an intermediate, SO$_3$ is predicted in the existing mechanism. To our best knowledge, however, no direct evidence has been reported in support of the generation of SO$_3$ in the oxidation of SO$_2$ by the OH radical in the atmosphere. Second, the rate of conversion in the new mechanism is expected to weakly depend on the humidity of air because only one water molecule is required for each SO$_2$. In contrast, the rate of conversion in the existing mechanism would be strongly dependent on the humidity because multiple water molecules are required for the conversion process of each SO$_2$.

Ⅴ. ACKNOWLEDGMENTS

This work was partially funded by National Science Foundation of the United States (No.1012994) and by California State University, Fullerton.

Reference
[1] B. J. Finlayson-Pitts, and J. N. Pitts Jr, Chemistry of the Upper and Lower Atmosphere. San Diego: Academic Press (2000).
[2] K. P. Wayne, Chemistry of Atmospheres. 3rd Edn. Oxford: Oxford University Press (2000).
[3] J. H Seifeld, and S. N. Pandis, Atmospheric Chemistry and Physics. New York: Wiley (1998).
[4] Acid Deposition, Atmospheric Process in Eastern North America National Research Council, Washington. DC: National Academy Press (1983).
[5] J. G. Calvert, A. Lazrus, G. L. Kok, B. G. Heikes, J. G. Walega, J. Lind, and C. A. Cantrell, Nature 317 , 27 (1985). DOI:10.1038/317027a0
[6] H. Berrwsheim, P. H. Wine, and D. D. Davis, in Composition, Chemistry and Climate of the Atmosphere, H. B. Singh Ed., Washington, DC: Van Nostrand Reinhold, (1994).
[7] W. R. Stockwell, and J. G. Calvert, Atmos. Environ. 17 , 2231 (1983). DOI:10.1016/0004-6981(83)90220-2
[8] S. W. Benson, Chem. Rev. 78 , 23 (1978). DOI:10.1021/cr60311a003
[9] S. Nagase, S. Hashimoto, and H. Akimoto, J. Phys. Chem. 92 , 641 (1988). DOI:10.1021/j100314a015
[10] W. K. Li, and M. L. McKee, J. Phys. Chem. 101 , 9778 (1997). DOI:10.1021/jp972389r
[11] D. Majumdar, G. S. Kim, J. Kim, K. S. Oh, and J. Y. Lee, J. Chem. Phys. 112 , 723 (2000). DOI:10.1063/1.480605
[12] W. Klopper, D. P. Tew, N. Gonzalez, - Garcia, and M. Olzmann, J. Chem. Phys. 129 , 114308 (2008). DOI:10.1063/1.2973637
[13] N. Gonzalez, - Garcia, W. Klopper, and M. Olzmann, Chem. Phys. Lett. 470 , 59 (2009). DOI:10.1016/j.cplett.2009.01.045
[14] K. Morokuma, and C. Muguruma, J. Am. Chem. Soc. 116 , 10316 (1994). DOI:10.1021/ja00101a068
[15] L. J. Larson, M. Kuno, and F. M. Tao, J. Chem. Phys. 112 , 8830 (2000). DOI:10.1063/1.481532
[16] C. E. Kolb, J. T. Jayne, D. R. Worsnop, M. J. Molina, R. J. Meads, and A. A. Viggiano, J. Am. Chem. Soc. 116 , 10314 (1994). DOI:10.1021/ja00101a067
[17] E. V. Akhmatskaya, C. J. Apps, I. H. Hillier, A. J. Masters, N. E. Wat, and J. C. Whitehead, Chem. Commun. , 707 (1997).
[18] E. T. Aaltonen, and J. S. Francisco, J. Phys. Chem. A 107 , 1216 (2003). DOI:10.1021/jp021977o
[19] A. D. Becke, J. Chem. Phys. 98 , 5648 (1993). DOI:10.1063/1.464913
[20] C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37 , 785 (1988). DOI:10.1103/PhysRevB.37.785
[21] G. A. Petersson, A. Bennett, T. G. Tensfeldt, M. A. AlLaham, W. A. Shirley, and J. A. Mantzaris, J. Chem. Phys. 89 , 2193 (1988). DOI:10.1063/1.455064
[22] R. A. Kendall, T. H. Dunning Jr, and R. J. Harrison, J. Chem. Phys. 96 , 6796 (1992). DOI:10.1063/1.462569
[23] C. Moller, and M. S. Plesset, Phys. Rev. 46 , 618 (1934). DOI:10.1103/PhysRev.46.618
[24] G. D. Purvis, and R. J. Bartlett, J. Chem. Phys. 76 , 1910 (1982). DOI:10.1063/1.443164
[25] J. F. Stanton, J. Gauss, J. D. Watts, and R. J. Bartlett, J. Chem. Phys. 94 , 4334 (1991). DOI:10.1063/1.460620
[26] K. Raghavachari, G. W. Truck, and J. A. Pople, M. Head-Gordon, Chem. Phys. Lett. 157 , 479 (1989). DOI:10.1016/S0009-2614(89)87395-6
[27] J. D. Watts, J. Gauss, and R. J. Bartlett, J. Chem. Phys. 98 , 8718 (1993). DOI:10.1063/1.464480
[28] C. Peng, and H. B. Schlegel, Israel J. Chem. 33 , 449 (1993). DOI:10.1002/ijch.v33.4
[29] C. Peng, P. Y. Ayala, H. B. Schlegel, and M. J. Frisch, J. Comp. Chem. 17 , 49 (1996). DOI:10.1002/(ISSN)1096-987X
[30] K. Fukui, Acc. Chem. Res. 14 , 363 (1981). DOI:10.1021/ar00072a001
[31] H. P. Hratchian and H. B. Schlegel, in Theory and Applications of Computational Chemistry: The First 40 Years, C. E. Dykstra, G. Frenking, K. S. Kim, and G. Scuseria Eds., Amsterdam: Elsevier, 195(2005).
[32] M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M.Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery Jr, J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, Ö. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and D. J. Fox, Gaussian 09, Revision B.01. Wallingford CT: Gaussian Inc. (2009).
大气中硫酸生成反应势能面的量子化学研究
陶琳, 李宇杰, 萨曼莎-索里亚诺, 陶福明     
美国加州州立大学化学和生物化学系, 加利福尼亚州, 富勒顿 92834
摘要: 本文运用密度泛函理论和高水平分子轨道从头计算理论,研究大气中由OH+SO2反应生成硫酸的一个新的势能面.新的势能面集中考虑自由基络合物SO3·HO2与H2O直接反应的一个通道,其主要原因是SO3·HO2的高度稳定性.其后续反应的最终结果与文献中流行的反应机理的结果相一致,但避免了产生与释放SO2.整个反应通道由三个连续的基元反应组成,(1)HOSO2+O2→SO3·HO2,(2)SO3·HO2+H2O→SO3·H2O·HO2,(3)SO3·H2O·HO2→H2SO4+HO2.这三个基元反应都只有很小的能垒,小于10 kcal/mol,并且都是放热反应.因此,这一新的反应通道在动力学及热力学上都是可行的.反应(3)是关键的一步,其中HO2作为一个桥梁分子促成了低能垒氢迁移过程而使SO3水解成为硫酸.这一研究结果主要有两个方面的大气化学意义.第一,大气中OH自由基氧化SO2并不产生SO3.第二,大气中SO2转化为硫酸受空气湿度的影响较小.
关键词: 量子化学    大气化学    二氧化硫    硫酸    计算化学    密度泛函理论    从头算方法