Unraveling the atmospheric oxidation mechanism and kinetics of naphthalene: Insights from theoretical exploration

Abstract

Naphthalene, the most abundant polycyclic aromatic hydrocarbon in the atmosphere, significantly influences OH consumption and secondary organic aerosol (SOA) formation. Naphthoquinone (NQ) is a significant contributor to ring-retaining SOA from naphthalene degradation, impacting the redox properties and toxicity of ambient particles. However, inconsistencies persist regarding concentrations of its isomers, 1,2-NQ and 1,4-NQ. In present work, our theoretical investigation into naphthalene's reaction with OH and subsequent oxygenation unveils their role in SOA formation. The reaction kinetics of initial OH and subsequent O2 oxidation was extensively studied using high-level quantum chemical methods (DLPNO-CCSD(T)/aug-ccpVQZ//M052x-D3/6–311++G(d,p)) combined with RRKM/master equation simulations. The reactions mainly proceed through electrophilic addition and abstraction from the aromatic ring. The total rate coefficient of naphthalene + OH at 300 K and 1 atm from our calculation (7.2 × 10−12 cm3 molecule−1 s−1) agrees well with previous measurements (∼1 × 10−11 cm3 molecule−1 s−1). The computed branching ratios facilitate accurate product yield determination. The largest yield of 1-hydroxynaphthalen-1-yl radical (add1) producing the major precursor of RO2 is computed to be 93.8 % in the ambient environment. Our calculated total rate coefficient (5.2 × 10−16 cm3 molecule−1 s−1) for add1 + O2 closely matches that of limited experimental data (8.0 × 10−16 cm3 molecule−1 s−1). Peroxy radicals (RO2) generated from add1 + O2 include 4-cis/trans-(1-hydroxynaphthalen-1-yl)-peroxy radical (add1-4OOadd-cis/trans, 66.0 %/17.5 %), 2-cis/trans-(1-hydroxynaphthalen-1-yl)-peroxy radical (add1-2OOadd-cis/trans, 10.3 %/6.3 %). Regarding the debated predominance of 1,4-NQ (corresponding to the parent RO2, i.e., add1-4OOadd-cis/trans) and 1,2-NQ (corresponding to the parent RO2, i.e., add1-2OOadd-cis/trans) in the atmosphere, our findings substantiate the dominance of 1,4-NQ. This study also indicates potential weakening of 1,4-NQ's dominance due to competition from decomposition reactions of add1-4OOadd-cis/trans and add1-2OOadd-cis/trans. Precise reaction kinetics data are essential for characterizing SOA transformation derived from naphthalene and assessing their climatic impacts within modeling frameworks.