大气中有机过氧自由基的反应动力学与机理研究综述

大气中有机过氧自由基的反应动力学与机理研究综述有机过氧自由基（RO₂）是大气光化学过程中关键的活性中间体，其反应动力学与反应机理直接决定了对流层臭氧、二次有机气溶胶（SOA）等污染物的生成效率，同时在全球碳循环、自由基循环过程中扮演着不可替代的角色。工业源排放的挥发性有机物（VOCs）、植物释放的异戊二烯、单萜烯等天然VOCs，在大气中经羟基自由基（OH）氧化、臭氧分解或光解过程生成不同结构的RO₂，其后续反应路径的差异直接影响区域空气质量与气候变化效应。近年来，随着同步辐射光电离质谱、激光诱导荧光等原位检测技术的发展，以及量子化学计算与动力学模拟方法的进步，学界对RO₂的反应机制有了诸多新认知，修正了传统模式中对部分关键反应参数的假设，为完善大气化学模式、准确评估污染物转化规律提供了重要支撑。一、大气中有机过氧自由基的来源与结构分类大气中RO₂的生成路径具有显著的源排放依赖性，不同前体物对应的RO₂结构差异直接决定其反应活性。对于烷烃类VOCs，其与OH自由基的反应以夺氢为主，生成的烷基自由基快速与氧气结合生成饱和RO₂，例如丙烷氧化生成的正丙基过氧自由基、异丙基过氧自由基，这类RO₂的α位无取代基团，反应活性相对稳定，是城市大气中RO₂的重要组成部分。烯烃类VOCs的反应路径更为复杂，OH既可以发生夺氢反应，也可以与碳碳双键发生加成反应，生成的β-羟基烷基自由基与氧气结合后得到羟基取代的RO₂，这类RO₂分子内存在氢键作用，其自反应与交叉反应速率常数与饱和RO₂存在数量级差异。天然源VOCs贡献了全球约90%的VOCs排放总量，其中异戊二烯每年排放量超过5亿吨，其氧化生成的RO₂结构最为多样。异戊二烯与OH加成后可生成4种异构的羟烷基自由基，与氧气结合后分别得到对应位置的过氧自由基，其中1,2-加成生成的RO₂在低NOₓ环境下易发生分子内氢迁移（即自氧化过程），生成高氧化态多官能团化合物（HOMs），是森林地区SOA生成的关键前体物。单萜烯如α-蒎烯的氧化过程更为复杂，臭氧与双键反应生成的Criegee中间体既可以直接分解生成RO₂，也可以发生环化反应形成双环过氧自由基，这类自由基的分子量更大，反应过程更容易生成低挥发性产物。除了VOCs氧化过程，大气中其他活性物种的反应也会生成特殊结构的RO₂。例如卤代烃在海洋边界层中经卤素自由基氧化生成含氯、溴取代的RO₂，这类自由基的反应活性更高，其与NO反应生成的卤代硝酸盐会影响大气卤素循环。芳香烃的氧化过程中，OH与苯环加成生成的环己二烯型自由基与氧气结合后生成芳香类RO₂，这类自由基容易发生开环反应，生成二羰基类化合物，是城市地区光化学烟雾的重要前体物。近年来研究发现，大气中的RO₂还可以通过非均相反应生成，例如吸附在气溶胶表面的VOCs被臭氧氧化后，过氧自由基可以从颗粒物表面脱附进入气相，这部分源在高污染地区的贡献可达10%以上。二、有机过氧自由基的主要反应通道与动力学参数RO₂的大气反应路径主要分为三类：与NO的反应、与HO₂的反应、与其他RO₂的自反应及交叉反应，不同反应通道的竞争关系由环境中NOₓ浓度、自由基浓度、温度、压力等因素共同决定。在城市等NOₓ浓度较高的区域，RO₂与NO的反应是主导路径，反应速率常数通常在10⁻¹¹~10⁻¹²cm³·molecule⁻¹·s⁻¹量级。该反应主要有两条分支：一条是氧化NO为NO₂，同时生成烷氧基自由基（RO），这是对流层臭氧生成的关键步骤，每生成一个NO₂在光照下即可分解产生一个氧原子，进而与氧气结合生成臭氧；另一条是生成有机硝酸酯（RONO₂），属于自由基的汇过程，会降低臭氧生成效率，同时生成的有机硝酸酯具有一定的低挥发性，是SOA的组成部分。不同结构RO₂与NO反应的硝酸酯分支比存在显著差异，这是当前模式模拟的重要不确定性来源。对于饱和直链烷烃生成的RO₂，碳链长度超过4个碳原子时，硝酸酯分支比随碳数增加线性上升，例如正辛基过氧自由基的硝酸酯分支比可达15%，而短链的乙基过氧自由基分支比仅为3%左右。羟基取代的RO₂由于分子内氢键的作用，硝酸酯分支比明显高于饱和RO₂，异戊二烯氧化生成的β-羟基RO₂的硝酸酯分支比最高可达25%，远高于模式中默认的10%假设值。近年来的实验研究还发现，对于含有多个官能团的RO₂，其与NO反应过程中还可能发生分子内的环化反应，生成环氧类化合物，该分支的占比最高可达10%，这一过程在之前的模式中完全被忽略，会导致臭氧生成模拟值偏高。在偏远地区如森林、海洋上空，NOₓ浓度低于100ppt时，RO₂与HO₂的反应成为主导路径，反应速率常数通常在10⁻¹¹~10⁻¹⁰cm³·molecule⁻¹·s⁻¹量级。该反应的主要产物是有机氢过氧化物（ROOH），这类化合物具有较低的挥发性，容易通过干沉降去除，或者在光照下光解重新释放RO₂和OH自由基，是大气中自由基的临时储库。对于α位含有羰基、羟基等吸电子基团的RO₂，与HO₂反应还会生成其他产物，例如乙醛基过氧自由基（CH₃C(O)O₂）与HO₂反应有20%的分支生成乙酸，这是大气中有机酸的重要来源之一。最新的研究发现，在高湿度条件下，RO₂与HO₂的反应过程会受到水分子的催化作用，反应速率可以提升30%以上，这一效应在热带雨林等高温高湿地区尤为显著，会显著改变当地的自由基收支平衡。RO₂之间的自反应和交叉反应是低NOₓ环境下另一类重要反应，饱和RO₂的自反应速率常数通常在10⁻¹⁴~10⁻¹⁶cm³·molecule⁻¹·s⁻¹量级，远低于与NO或HO₂的反应速率，但在RO₂浓度较高的污染边界层顶部、夜间化学过程中，这类反应的贡献不可忽视。RO₂自反应的产物非常复杂，通常包括醇、酮、醚类化合物，同时会生成部分RO自由基，继续参与后续氧化过程。对于酰基类RO₂，其自反应速率常数可以达到10⁻¹¹cm³·molecule⁻¹·s⁻¹，接近与HO₂的反应速率，因此在乙醛、丙酮等羰基化合物浓度较高的地区，这类自反应的贡献需要重点考虑。不同结构RO₂之间的交叉反应速率常数通常介于两种自由基自反应速率常数的几何平均值，但是近年来的实验发现，对于含有羟基的RO₂与酰基RO₂的交叉反应，其速率常数比几何平均值高2~3倍，这会导致模式中对这类反应的产物模拟出现偏差。除了上述气相反应，RO₂还可以发生分子内反应，即自氧化过程。当RO₂分子内存在容易被夺取的氢原子时，过氧基团可以通过六元环过渡态夺取相邻碳上的氢，生成氢过氧化物取代的烷基自由基，随后再与氧气结合生成新的过氧基团，这一过程可以连续发生多次，最终生成含有多个过氧羟基、羰基的高氧化态化合物，其O/C比通常大于0.8，饱和蒸气压低于10⁻⁵Pa，一旦生成即可几乎全部凝结进入颗粒相，是SOA生成的重要路径。异戊二烯氧化生成的RO₂在298K下的分子内氢迁移速率约为0.01s⁻¹，而单萜烯生成的双环RO₂的氢迁移速率可达1s⁻¹以上，即使在NOₓ浓度为1ppb的环境中，自氧化路径的占比仍可超过50%。温度对自氧化速率的影响非常显著，低温下反应速率会快速下降，因此这一过程在热带森林地区的贡献远高于寒带地区。三、有机过氧自由基反应机理的最新研究进展传统的大气化学模式中，通常假设RO₂的反应只发生在气相，且结构相似的RO₂具有相同的反应参数，但是近年来的研究发现，RO₂的非均相反应和构效关系对其反应机理有重要影响。RO₂在气溶胶表面的uptake过程是当前研究的热点，实验测定的RO₂在无机气溶胶表面的沉降系数约为10⁻³~10⁻²，在含有机质的气溶胶表面沉降系数更高。RO₂吸附在颗粒物表面后，既可以与颗粒相中的抗坏血酸、酚类等还原性物质发生反应，被还原为ROOH或醇类化合物，也可以与颗粒相中的其他自由基发生反应，生成低聚物。在高污染地区，气溶胶表面积浓度可达1000μm²·cm⁻³以上，RO₂的非均相反应寿命可低至几分钟，与气相反应寿命相当，这会显著降低气相中的RO₂浓度，进而抑制臭氧的生成。RO₂的构效关系研究为模式参数化提供了重要基础。通过量子化学计算和实验测定，研究者发现RO₂的反应速率常数与其最高占据分子轨道（HOMO）能量、电离能等参数存在显著的线性关系。对于与NO的反应，RO₂的HOMO能量越高，反应速率常数越大，因为反应过程中RO₂提供电子给NO的反键轨道。对于与HO₂的反应，反应速率常数与RO₂的电子亲和能相关，吸电子取代基团会降低RO₂的电子云密度，更容易与HO₂发生反应。基于构效关系建立的预测模型，可以估算不同结构RO₂的反应速率常数，误差通常在2倍以内，解决了大量复杂结构RO₂参数缺失的问题。近年来，同步辐射真空紫外光电离质谱技术的发展为RO₂反应的产物检测提供了有力手段。传统的检测方法如激光诱导荧光、化学电离质谱等只能检测特定种类的RO₂，而同步辐射光电离质谱可以实现对反应过程中所有物种的在线检测，通过调谐光子能量可以区分同分异构体。研究者利用该技术发现了RO₂反应过程中许多之前未被观测到的产物，例如异戊二烯RO₂自氧化过程中生成的含有4个过氧基团的HOMs，以及RO₂与NO反应生成的环氧硝酸酯类化合物。这些新发现的产物解释了传统模式中SOA生成模拟值偏低的问题，将异戊二烯SOA生成模拟值提高了2~3倍。夜间化学过程中RO₂的反应机理也取得了重要进展。夜间大气中OH自由基浓度很低，NO₃自由基成为主要的氧化剂，其与VOCs反应生成的RO₂与白天生成的RO₂结构存在差异，例如单萜烯与NO₃反应生成的是硝酸酯取代的RO₂。夜间环境中NO浓度通常很低，RO₂主要与HO₂、其他RO₂以及NO₃发生反应，其中RO₂与NO₃的反应速率常数约为10⁻¹²cm³·molecule⁻¹·s⁻¹量级，是夜间NO₃自由基的重要汇。该反应的产物主要是RO、NO₂和氧气，生成的RO可以继续被氧气氧化生成HO₂，维持夜间的自由基循环。最新的研究发现，夜间生成的含硝酸酯基团的RO₂更容易发生自氧化过程，生成的HOMs含有更多的含氮官能团，其挥发性更低，对夜间SOA的贡献可达40%以上。四、反应动力学参数在大气化学模式中的应用与不确定性RO₂的反应动力学参数是大气化学模式的核心输入参数，其准确性直接决定了臭氧、SOA等污染物的模拟效果。当前主流的区域空气质量模式如CMAQ、WRF-Chem中，RO₂的反应参数主要基于少量实验测定的结果，对于复杂结构的RO₂通常采用同系物近似的方法，导致模拟结果存在较大的不确定性。以有机硝酸酯的生成为例，早期模式中默认所有RO₂与NO反应的硝酸酯分支比为10%，而实际上不同结构RO₂的分支比在1%~30%之间，这导致模式对城市地区臭氧生成的模拟误差可达20%以上，对有机硝酸酯浓度的模拟误差甚至超过100%。RO₂自氧化过程的参数化是当前模式发展的重点。由于自氧化过程涉及多个连续的反应步骤，不同结构RO₂的氢迁移速率差异可达几个数量级，直接将所有反应加入模式会导致计算量显著增加。目前的参数化方案通常将RO₂分为几个大类，分别给出不同的自氧化速率和产率，例如将异戊二烯RO₂的自氧化分支比设为0.01s⁻¹，单萜烯RO₂的自氧化分支比设为0.5s⁻¹。这类参数化方案可以显著提升SOA的模拟效果，但是在不同温度、NOₓ浓度条件下的适用性仍需要进一步验证。非均相反应的参数化是另一个难点。目前模式中对于RO₂的非均相反应通常采用简化的沉降系数方法，但是实际气溶胶的成分、酸度、氧化态等性质都会影响RO₂的uptake系数，不同地区的气溶胶性质差异很大，采用统一的沉降系数会导致模拟结果出现偏差。例如在酸性气溶胶表面，RO₂更容易发生质子化反应，生成的阳离子自由基会发生不同的反应路径，产物分布与中性条件下存在显著差异。未来需要建立基于气溶胶成分的RO₂非均相反应参数化方案，进一步提高模式的模拟精度。不同反应通道的竞争关系对气象条