Mechanistic Insights into Sulfuric Acid Formation in the Atmosphere via First-Principles Simulations
Abstract
Sulfuric acid (H2SO4) is essential in the formation of atmospheric nucleation and cloudy condensation nuclei (CCN). This study uses density functional theory (DFT) to investigate the formation of H2SO4 from sulfur dioxide (SO2) via three key reactions: SO2 oxidation by hydroxyl radicals (R1), reaction of HOSO2 with O2 (R2), and hydrolysis of SO3 with water (R3). Reaction and activation energies were computed using the PBE, r2SCAN, DC-r2SCAN, CAM-B3LYP, and PBE0 flavors of density functional theory. The key findings reveal that reaction R1 has a calculated reaction energy (ΔG) of −23.845 kcal/mol and an activation energy (ΔG*) of −0.628 kcal/mol by using the DC-r2SCAN functional. Reaction R2, which was previously assumed to be barrierless, showed a small but significant activation barrier of 1.225 kcal/mol at the CAM-B3LYP/6-31G** level. Reaction R3 led to a calculated chemical reaction energy of −23.218 kcal/mol, with an activation energy of 5.648 kcal/mol using r2SCAN/TZV2P. This study demonstrates that r2SCAN and DC-r2SCAN provide a computationally efficient alternative to high-level methods, achieving comparable accuracy in the description of sulfuric acid formation.
Copyright and License
© 2025 American Chemical Society.
Acknowledgement
The authors acknowledge financial support from the Deanship of Research at Sultan Qaboos University (Grant No. RF/SCI/PHYS/25/01). The authors also thank Dr. Yang Liu (University of New Mexico) for sharing the Cartesian coordinates of the RC1 complex. We are grateful to the anonymous reviewers for their insightful comments and suggestions that significantly improved the manuscript.
Contributions
M.A.-K. and A.K. performed the simulations and data analysis. A.K. supervised the project. All authors contributed to the interpretation of the results and the writing of the manuscript.
Supplemental Material
Computational details, list of functionals and basis sets used for different reactions (Table 1); molecules involved in each reaction step (Table 2); vibrational frequencies (in cm–1) of various modes in Reaction R2 molecules (Table 3); infrared (IR) spectra for HOSO2 obtained at the r2SCAN/aug-TZV2P level (Figure S1); list of molecules (Table S4); LSD presence, and multiplicity for reactions R1–R3, molecule’s computational time for different functionals (Table 5) (PDF)
Files
jz5c02165_si_001.pdf
Files
(257.9 kB)
| Name | Size | Download all |
|---|---|---|
|
md5:655d9d18841f30e10ff980df7c0b9f42
|
257.9 kB | Preview Download |
Additional details
Identifiers
- PMID
- 40947615
Related works
- Describes
- Journal Article: 40947615 (PMID)
Funding
- Sultan Qaboos University
- RF/SCI/PHYS/25/01
Dates
- Accepted
-
2025-09-08
- Available
-
2025-09-15Published online