多环芳烃大气氧化与硝化反应机理的理论探究

多环芳烃大气氧化与硝化反应机理的理论探究一、引言1.1研究背景与意义随着工业化和城市化进程的飞速发展，大气污染已成为全球面临的严峻挑战之一。多环芳烃（PolycyclicAromaticHydrocarbons，PAHs）作为一类具有两个或两个以上苯环的碳氢化合物，广泛存在于大气环境中，其来源十分广泛，主要包括自然源和人为源。自然源涵盖火山爆发、森林火灾、生物合成等，这些过程会向大气中释放一定量的PAHs；而人为源则更为复杂多样，其中化石燃料（如煤、石油、天然气）的不完全燃烧是PAHs的主要人为来源之一。在工业生产过程中，例如钢铁冶炼、炼焦、石油化工等行业，化石燃料的燃烧会产生大量PAHs排放到大气中。交通运输领域，汽车尾气排放也含有相当数量的PAHs，尤其是在交通拥堵时，车辆发动机处于不完全燃烧状态，PAHs的排放会显著增加。此外，生物质燃烧（如秸秆焚烧、木材燃烧）以及垃圾焚烧等活动同样会释放PAHs，进一步加剧了大气中PAHs的污染程度。PAHs具有较强的稳定性和生物累积性，难以被生物降解，能够在环境中长期存在。并且，许多PAHs具有潜在的毒性、致癌性和致畸性，对人体健康和生态环境构成了严重威胁。例如，苯并[a]芘（Benzo[a]pyrene，BaP）是一种典型的致癌性PAHs，被国际癌症研究机构（IARC）列为第1类人类致癌物。长期暴露于含有PAHs的环境中，人类可能会增加患肺癌、皮肤癌等多种癌症的风险，同时也可能对生殖系统、免疫系统等造成损害。在生态环境方面，PAHs会对动植物的生长发育产生负面影响，破坏生态平衡。在大气环境中，PAHs并非稳定不变，而是会经历复杂的化学转化过程。其中，生成氧化多环芳烃（OxidizedPolycyclicAromaticHydrocarbons，OPAHs）和硝基多环芳烃（NitroPolycyclicAromaticHydrocarbons，NPAHs）的反应是PAHs大气化学转化的重要途径。OPAHs和NPAHs相较于母体PAHs，往往具有更强的毒性和生物活性。一些硝基多环芳烃的致癌性比相应的PAHs更高，对人体健康的危害更为严重。因此，深入研究PAHs经大气反应生成OPAHs和NPAHs的机理，对于全面理解PAHs在大气中的行为具有至关重要的意义。从大气化学过程的角度来看，明晰PAHs的转化机理有助于我们更准确地认识大气中复杂的化学反应网络。PAHs与大气中的氧化剂（如羟基自由基（・OH）、硝酸根自由基（・NO3）、臭氧（O3）等）发生反应，生成OPAHs和NPAHs，这些反应不仅影响着PAHs本身的浓度和分布，还会对大气中其他污染物的转化和迁移产生连锁反应。通过研究这些反应机理，我们能够深入了解大气中污染物的相互作用机制，为建立更完善的大气化学模型提供理论支持，从而更准确地预测大气环境质量的变化趋势。在环境风险评估方面，准确掌握PAHs生成OPAHs和NPAHs的机理对于评估其对生态环境和人体健康的潜在风险至关重要。不同的反应途径和产物具有不同的毒性和环境行为，只有明确了这些机理，才能更精准地评估PAHs及其转化产物在环境中的迁移、扩散和累积规律，以及它们对生态系统和人类健康的危害程度，为制定科学合理的环境风险评估标准和方法提供依据。对于污染控制而言，研究PAHs的转化机理为制定有效的污染控制策略提供了关键的科学依据。了解PAHs生成OPAHs和NPAHs的具体反应过程和影响因素，我们可以有针对性地采取措施，抑制这些有害产物的生成，减少PAHs对大气环境的污染。通过优化燃烧条件，降低化石燃料不完全燃烧产生的PAHs排放；研发高效的污染治理技术，去除大气中的PAHs及其转化产物。这对于改善大气环境质量、保护人类健康和生态平衡具有重要的现实意义。1.2国内外研究现状多环芳烃（PAHs）在大气环境中的行为及转化一直是国内外环境科学领域的研究热点。国内外众多学者围绕PAHs经大气反应生成氧化多环芳烃（OPAHs）和硝基多环芳烃（NPAHs）开展了大量研究工作，取得了一系列重要成果。在国外，相关研究起步较早。早期的研究主要集中在对PAHs及其衍生物在大气中的浓度水平、分布特征以及来源解析等方面。随着研究的深入，学者们逐渐关注到PAHs的大气化学反应过程。通过实验室模拟实验和野外观测相结合的方法，对PAHs与大气中常见氧化剂（如・OH、・NO3、O3等）的反应进行了广泛研究。例如，有研究利用烟雾箱模拟系统，详细研究了萘在・OH自由基和NO2共同作用下生成硝基萘的反应过程，发现反应条件（如温度、湿度、反应物浓度等）对产物的种类和生成速率有显著影响。在理论计算方面，量子化学方法被广泛应用于研究PAHs大气反应的机理和动力学性质。通过计算反应势能面、反应势垒和速率常数等参数，深入揭示了反应的微观过程。研究表明，PAHs与・OH自由基的加成反应是生成OPAHs的重要途径，而与・NO3自由基的反应则更倾向于生成NPAHs，且不同PAHs分子的反应活性和选择性存在差异。国内对PAHs大气反应的研究近年来也取得了长足进展。在观测研究方面，对多个城市和地区大气中的PAHs、OPAHs和NPAHs进行了监测分析，明确了其污染现状和时空变化规律。研究发现，我国部分地区大气中PAHs和NPAHs的浓度处于较高水平，且在冬季采暖期由于化石燃料燃烧量增加，浓度显著升高。在实验研究方面，一些科研团队搭建了自主研发的实验装置，开展了PAHs大气反应的模拟实验，探究了不同反应条件下OPAHs和NPAHs的生成机制。在理论计算领域，国内学者也运用量子化学方法，对典型PAHs的大气反应机理进行了深入研究，为理解PAHs在我国大气环境中的转化行为提供了理论支持。然而，当前的研究仍存在一些不足之处和待解决的问题。在实验研究方面，虽然已经开展了大量的模拟实验，但实验条件与实际大气环境仍存在一定差距，难以完全真实地反映PAHs在大气中的复杂反应过程。实际大气中存在多种污染物和复杂的气象条件，这些因素之间的相互作用对PAHs反应的影响尚未得到充分研究。此外，现有的检测技术在对一些痕量的OPAHs和NPAHs的准确检测方面还存在一定困难，限制了对其生成和转化规律的深入认识。在理论计算方面，虽然量子化学方法为研究反应机理提供了有力工具，但计算模型和方法仍有待进一步完善。目前的计算大多基于理想的气相反应体系，忽略了大气中颗粒物表面等非均相反应的影响，而实际大气中颗粒物表面的反应对PAHs的转化可能起着重要作用。不同计算方法和基组的选择对计算结果也存在一定影响，导致不同研究之间的结果存在差异，缺乏统一的标准和比较依据。在PAHs大气反应的综合研究方面，目前实验研究和理论计算之间的结合还不够紧密。实验结果往往难以直接与理论计算结果相互验证和补充，限制了对PAHs大气反应机理的全面理解。此外，对于PAHs经大气反应生成OPAHs和NPAHs后的环境归趋和生态毒性等方面的研究还相对较少，无法准确评估其对环境和人体健康的潜在风险。综上所述，虽然国内外在PAHs经大气反应生成OPAHs和NPAHs的研究方面取得了一定成果，但仍存在诸多不足。进一步深入研究PAHs的大气反应机理，加强实验研究与理论计算的结合，完善检测技术和计算方法，以及拓展对反应产物环境归趋和生态毒性的研究，对于全面理解PAHs在大气中的行为、准确评估其环境风险以及制定有效的污染控制策略具有重要意义。1.3研究目标与内容本研究旨在通过理论计算和模拟，深入探讨典型多环芳烃（PAHs）在大气中经氧化、硝化反应生成氧化多环芳烃（OPAHs）和硝基多环芳烃（NPAHs）的详细机理，精确确定各反应环节的速率常数。具体研究内容涵盖以下几个方面：PAHs在大气中的来源及其危害性评估：全面综述PAHs在大气环境中的来源，不仅详细分析自然源如火山爆发、森林火灾、生物合成等过程中PAHs的产生机制和释放量，还深入探讨人为源，包括化石燃料（煤、石油、天然气）的不完全燃烧、工业生产（钢铁冶炼、炼焦、石油化工等）、交通运输（汽车尾气排放）以及生物质燃烧（秸秆焚烧、木材燃烧）、垃圾焚烧等活动对PAHs排放的贡献。综合评估PAHs对人体健康和生态环境的危害，结合大量的毒理学研究数据和实际案例，阐述PAHs的致癌性、致畸性和致突变性，以及对动植物生长发育、生态系统平衡的破坏作用。PAHs与大气氧化物的反应机理：运用量子化学密度泛函理论（DFT），对苯并芘、蒽等典型PAHs分别与大气中重要的氧化物，如羟基自由基（・OH）、硝酸根自由基（・NO3）、臭氧（O3）等的反应进行深入分析。精确计算反应的势能面，确定反应过程中的反应物、中间体、过渡态和产物的结构和能量，进而获得反应势垒。通过对反应势垒的分析，明确反应的难易程度和主要反应路径。同时，结合过渡态理论等动力学方法，计算各反应的速率常数，定量描述反应进行的快慢，为理解PAHs在大气中的转化速率提供关键数据。OPAHs和NPAHs的形成机理：紧密结合理论计算结果、已有的实验数据以及大气化学模型，深入探讨PAHs在大气中经过氧化、硝化反应后生成OPAHs和NPAHs的详细机理。分析不同反应条件（如温度、湿度、反应物浓度、光照强度等）对反应路径和产物分布的影响，明确各反应步骤的先后顺序和相互关系。例如，研究・OH自由基与PAHs的加成反应如何引发后续的氧化过程，以及・NO3自由基与PAHs的反应如何导致NPAHs的生成。确定各环节的速率常数，构建完整的反应动力学模型，通过模型预测不同环境条件下OPAHs和NPAHs的生成量和变化趋势，为大气污染的预测和控制提供科学依据。1.4研究方法与技术路线为深入探究典型多环芳烃（PAHs）经大气反应生成氧化多环芳烃（OPAHs）和硝基多环芳烃（NPAHs）的机理，本研究将综合运用理论计算和模拟计算的方法。在理论计算方面，采用量子化学密度泛函理论（DFT）。密度泛函理论是一种基于电子密度的量子力学计算方法，能够有效地处理多电子体系的问题。通过在Gaussian等量子化学计算软件中运用DFT方法，选择合适的基组（如6-31G(d,p)、6-311++G(d,p)等），对苯并芘、蒽等典型PAHs与大气中常见氧化物（如羟基自由基（・OH）、硝酸根自由基（・NO3）、臭氧（O3）等）的反应进行精确计算。首先，对反应物、中间体、过渡态和产物的几何结构进行全优化，通过频率分析确定各结构是否为稳定结构（稳定结构的频率均为正值）或过渡态（过渡态有且仅有一个虚频）。然后，计算反应过程中的势能面，获得反应势垒。反应势垒是反应进行的能量障碍，通过比较不同反应路径的势垒高低，可以判断反应的主要进行方向。此外，结合过渡态理论（TST）等动力学方法，利用计算得到的结构和能量信息，计算各反应的速率常数。速率常数能够定量描述反应进行的快慢程度，对于理解PAHs在大气中的转化速率具有重要意义。通过DFT计算，可以深入了解PAHs与大气氧化物反应的微观过程，为后续的机理分析提供坚实的理论基础。在模拟计算方面，采用大气化学模式，如CommunityMultiscaleAirQuality（CMAQ）模型、WeatherResearchandForecastingwithChemistry（WRF-Chem）模型等。这些模型能够综合考虑大气中各种物理和化学过程，模拟污染物在大气中的传输、扩散和转化行为。结合实验数据，将理论计算得到的反应机理和速率常数参数化后输入到大气化学模式中。实验数据可以来自实验室模拟实验，如利用烟雾箱模拟系统研究PAHs在特定条件下的反应过程，获取不同反应条件下OPAHs和NPAHs的生成量和产物分布等信息；也可以来自野外观测数据，对实际大气环境中的PAHs、OPAHs和NPAHs进行监测分析，得到其浓度水平、时空变化规律等数据。通过大气化学模式的模拟计算，可以预测不同环境条件下PAHs经大气反应生成OPAHs和NPAHs的行为，如在不同的气象条件（温度、湿度、风速、光照等）和污染源排放情况下，PAHs及其转化产物在大气中的浓度分布和变化趋势。将模拟结果与实验数据进行对比验证，进一步优化和完善反应机理和模型参数，提高对PAHs大气反应过程的模拟准确性和预测能力。本研究的技术路线如图1-1所示：首先，全面收集和整理PAHs在大气中的来源、危害性以及相关实验数据等资料；然后，运用量子化学密度泛函理论进行理论计算，得到PAHs与大气氧化物反应的机理和速率常数；接着，将理论计算结果与实验数据相结合，输入到大气化学模式中进行模拟计算；最后，对模拟结果进行分析和验证，深入探讨PAHs生成OPAHs和NPAHs的机理，并提出相应的结论和建议。通过这种综合的研究方法和技术路线，有望全面、深入地揭示PAHs经大气反应生成OPAHs和NPAHs的机理，为大气污染的防控提供科学依据。[此处插入技术路线图]图1-1技术路线图二、多环芳烃及相关化合物概述2.1多环芳烃的结构与性质多环芳烃（PAHs）是指分子中含有两个或两个以上苯环的碳氢化合物，根据苯环的连接方式，PAHs主要可分为两类。一类是非稠环型，其中包括联苯及联多苯和多苯代脂肪烃。联苯是由两个苯环通过单键直接相连而成，其结构相对较为简单。而多苯代脂肪烃则是由若干个苯环取代脂肪烃中的氢原子形成，例如1,2-二苯基乙烷，两个苯环通过一个乙基相连。另一类是稠环型，即两个或多个苯环共用两个相邻碳原子稠合而成，如萘、蒽、菲、苯并[a]芘等。萘由两个苯环共用两个相邻碳原子稠合而成，是最简单的稠环芳烃；蒽由三个苯环直线稠合而成；菲也是由三个苯环稠合而成，但与蒽的稠合方式不同，具有独特的角状结构；苯并[a]芘则是一种具有五个苯环的稠环芳烃，其结构更为复杂。这种稠环结构使得PAHs具有特殊的电子云分布和空间构型，进而影响其物理化学性质和化学反应活性。[此处插入多环芳烃的分类图]图2-1多环芳烃的分类PAHs的物理性质与其分子结构密切相关。一般来说，随着苯环数量的增加，PAHs的分子量增大，分子间作用力增强，其熔点和沸点相应提高。萘的熔点为80.26℃，沸点为217.96℃；而苯并[a]芘的熔点高达179℃，沸点为495℃。在常温下，PAHs大多为固体或高熔点物质，且不溶于水。这是因为PAHs分子是非极性或弱极性的，而水是极性分子，根据相似相溶原理，PAHs在水中的溶解度极低。萘在25℃时水中的溶解度仅为31.7mg/L。然而，PAHs可溶于苯、甲苯、二氯甲烷等有机溶剂。在苯中，萘能与苯分子通过分子间作用力相互作用，从而实现溶解。PAHs还具有一定的挥发性，尤其是分子量较小的PAHs，如萘、苊等，在常温下具有相对较高的蒸气压，容易挥发进入大气环境。PAHs具有较高的化学稳定性，这是由于其共轭π电子体系的存在。共轭π电子体系使得电子在整个分子平面内离域，降低了分子的能量，增强了分子的稳定性。在一般条件下，PAHs不易发生化学反应。在一定条件下，如高温、光照或氧化剂存在时，PAHs可以发生化学反应。在高温和氧气存在的条件下，PAHs能够发生氧化反应，生成各种氧化产物。在光照条件下，PAHs可以吸收光能，激发电子跃迁到较高能级，从而引发光化学反应。PAHs还可以与一些亲电试剂发生亲电取代反应，例如与硝酸发生硝化反应，生成硝基多环芳烃。PAHs也可以通过生物降解过程在微生物的作用下分解为较小的分子。一些微生物能够利用PAHs作为碳源和能源，通过一系列酶的作用，将PAHs逐步降解为无害的物质。白腐真菌能够分泌木质素过氧化物酶和锰过氧化物酶等，这些酶可以催化PAHs的氧化分解。在大气中，PAHs以气、固两种形式存在，其存在形态与PAHs的分子量密切相关。分子量小的2-3环PAHs，如萘、苊烯等，由于其蒸气压相对较高，主要以气态形式存在。而5-7环的大分子量PAHs，如苯并[a]芘、二苯并[a,h]蒽等，蒸气压较低，绝大部分以颗粒态形式存在，吸附在大气颗粒物表面。4环PAHs在气态、颗粒态中的分配基本相同。大气中的PAHs会受到多种因素的影响，如气象条件、污染源排放等。在温度较低、风速较小的情况下，PAHs更容易在大气中积累，尤其是颗粒态PAHs。而在高温、大风的条件下，PAHs的扩散和稀释作用增强，其在大气中的浓度会相应降低。不同地区的大气中PAHs的浓度和组成也存在差异，工业发达地区和交通繁忙地区的大气中PAHs浓度通常较高。2.2氧化多环芳烃和硝基多环芳烃的特性氧化多环芳烃（OPAHs）是多环芳烃（PAHs）在大气中与氧化剂发生氧化反应的产物，其分子结构在PAHs的基础上引入了含氧官能团，如羟基（-OH）、羰基（C=O）、羧基（-COOH）等。这些含氧官能团的引入显著改变了分子的电子云分布和空间构型，进而赋予OPAHs独特的物理化学性质。在物理性质方面，OPAHs的熔点和沸点通常比母体PAHs更高。这是因为含氧官能团的存在增强了分子间的作用力，使得分子间的相互作用更加紧密，需要更高的能量才能使分子脱离晶格或液态环境，从而导致熔点和沸点升高。一些含有多个羟基或羰基的OPAHs，其分子间可以形成氢键，进一步增强了分子间的相互作用，使得熔点和沸点显著提高。在溶解性方面，由于含氧官能团的极性，OPAHs相较于PAHs具有更强的亲水性，在水中的溶解度有所增加。这使得OPAHs在大气中的迁移和转化过程与PAHs有所不同，更容易随着大气中的水汽循环进行扩散和传输，也更容易在水体中积累，对水环境造成潜在威胁。在化学性质上，OPAHs具有较高的化学反应活性。含氧官能团的存在使得OPAHs更容易参与各种化学反应，如与大气中的其他污染物发生二次反应，生成更复杂的化合物。含有羰基的OPAHs可以与大气中的胺类物质发生反应，生成含氮的有机化合物。这种较高的反应活性使得OPAHs在大气化学过程中扮演着重要的角色，可能会影响大气中其他污染物的转化和迁移，对大气环境质量产生深远的影响。硝基多环芳烃（NPAHs）则是PAHs与含有硝基（-NO2）的化合物发生反应生成的产物。NPAHs的结构中，硝基直接连接在PAHs的苯环上。硝基的强吸电子性对NPAHs的物理化学性质产生了显著影响。从物理性质来看，NPAHs的熔点和沸点同样受到硝基的影响而发生变化。由于硝基的引入增加了分子的极性，使得分子间作用力增强，NPAHs的熔点和沸点一般比母体PAHs有所升高。在溶解性方面，NPAHs的极性增强导致其在水中的溶解度相较于PAHs有所增加，但仍属于微溶性物质。在化学性质上，NPAHs具有独特的反应活性。硝基的存在使得NPAHs的苯环电子云密度降低，亲电取代反应活性下降，但同时增加了其亲核取代反应的活性。NPAHs可以在一定条件下发生还原反应，硝基被还原为氨基（-NH2），生成氨基多环芳烃。这种反应在大气环境中可能会受到光照、还原剂等因素的影响。NPAHs还可以与大气中的一些自由基发生反应，参与大气中的化学反应网络。由于硝基的强吸电子性，NPAHs的电子结构发生改变，使其在光化学反应中具有特殊的行为，可能会产生一些具有毒性的光解产物。OPAHs和NPAHs具有比母体PAHs更强的毒性和生物活性。许多研究表明，它们对生物体的危害更为严重。在致癌性方面，一些NPAHs的致癌性比相应的PAHs更高。例如，1-硝基芘被认为是一种强致癌物，其致癌活性远高于芘本身。这是因为硝基的存在改变了分子与生物大分子（如DNA）的相互作用方式，更容易引发基因突变和细胞癌变。在致突变性方面，OPAHs和NPAHs能够诱导生物体的基因突变，干扰细胞的正常代谢和遗传信息传递。实验研究发现，某些OPAHs可以与DNA分子发生共价结合，形成DNA加合物，导致DNA损伤和基因突变，增加生物体患癌症和其他疾病的风险。在环境影响方面，OPAHs和NPAHs在大气中的存在会对空气质量产生负面影响。它们可以参与大气中的光化学反应，生成二次气溶胶，增加大气颗粒物的浓度，降低大气能见度，导致雾霾等天气现象的发生。OPAHs和NPAHs还可以通过大气传输，在不同地区之间扩散，对区域乃至全球的大气环境质量产生影响。在水体和土壤环境中，OPAHs和NPAHs的积累会对生态系统造成破坏。它们可以通过食物链的传递，在生物体内富集，对动植物的生长发育、繁殖和生存产生不利影响。在土壤中，OPAHs和NPAHs可能会抑制土壤微生物的活性，影响土壤的肥力和生态功能。2.3多环芳烃在大气中的来源与分布多环芳烃（PAHs）在大气中的来源广泛，可分为自然源和人为源。自然源方面，火山爆发是重要的释放途径之一。在火山喷发过程中，地球内部的高温岩浆与周围的岩石、有机物等相互作用，发生复杂的物理和化学变化，促使PAHs的形成并大量排放到大气中。据研究，每年全球火山喷发向大气中释放的苯并芘可达12-14吨。森林火灾同样不可忽视，当森林植被燃烧时，木材、树叶等有机物在不完全燃烧的条件下会产生PAHs。生物质中的纤维素、木质素等成分在高温缺氧环境下热解，形成碳氢自由基和碎片，这些活性中间体进一步反应聚合，生成各种PAHs。森林火灾释放的PAHs不仅影响当地大气环境，还可能随着大气环流进行长距离传输，对周边地区甚至全球大气环境产生影响。生物合成也是自然源的一部分，陆地、水生植物和微生物在自身代谢过程中能够合成PAHs。某些微生物在特定的生长环境和代谢途径下，会产生一些PAHs作为代谢副产物。虽然生物合成产生的PAHs量相对较少，但在自然生态系统中也构成了PAHs的天然本底值的一部分。人为源是大气中PAHs的主要来源，涵盖多个领域。在能源燃烧方面，煤、石油、天然气等化石燃料的不完全燃烧是PAHs的重要人为来源。以煤为例，煤中含有大量的有机质，在燃烧过程中，如果燃烧条件不佳，如氧气供应不足、燃烧温度不均匀等，就会导致煤中的碳氢化合物无法完全氧化，进而热解聚合形成PAHs。工业锅炉、炉灶等使用煤作为燃料时，会排放出含有PAHs的烟气。石油在精炼和燃烧过程中也会产生PAHs，汽车发动机中汽油或柴油的燃烧，由于燃烧过程的复杂性和不完全性，会产生多种PAHs，随着尾气排放到大气中。交通运输行业是PAHs的重要排放源之一。汽车尾气中含有大量的PAHs，这是因为汽车发动机在运转时，燃料与空气的混合不均匀，以及燃烧室内的高温高压环境，导致燃料不完全燃烧，产生PAHs。在交通拥堵路段，车辆频繁启停，发动机处于怠速或低速运转状态，燃烧效率降低，PAHs的排放会显著增加。飞机、轮船等交通工具使用的燃油在燃烧时也会产生PAHs，虽然它们的排放总量相对汽车尾气可能较少，但由于其排放高度较高，对大气环境的影响范围更广。工业生产过程中，许多行业都会排放PAHs。钢铁冶炼过程中，铁矿石的烧结、焦炭的生产以及钢铁的锻造等环节，都会涉及到高温和有机物的参与，容易产生PAHs。在焦炭生产过程中，煤在高温干馏条件下，会分解产生大量的PAHs，随着废气排放到大气中。炼焦厂是PAHs的高排放源之一，其排放的PAHs不仅种类繁多，而且浓度较高。石油化工行业，如炼油、化工合成等过程，也会产生PAHs。在石油炼制过程中，原油中的重质组分在高温裂解和催化重整等工艺中，会发生复杂的化学反应，生成PAHs。化工合成过程中，一些有机原料的反应也可能产生PAHs作为副产物。垃圾焚烧也是大气中PAHs的一个重要来源。城市固体废弃物中含有大量的塑料、纸张、木材、厨余垃圾等，这些物质在焚烧过程中，如果焚烧温度不够高或燃烧时间不足，就会发生不完全燃烧，产生PAHs。医院垃圾中可能含有各种化学药品、塑料器具等，其焚烧过程中产生的PAHs可能还含有其他有毒有害物质，对大气环境和人体健康的危害更大。大气中PAHs的浓度分布在不同地区和环境中存在显著差异。在工业发达地区，由于大量工业企业的存在，如钢铁厂、炼焦厂、化工厂等，这些企业在生产过程中会向大气中排放大量的PAHs，导致该地区大气中PAHs浓度较高。一些传统的工业城市，其大气中PAHs的浓度明显高于周边的非工业地区。北京、天津、上海等城市的部分工业区域，大气中PAHs的浓度可达数百纳克每立方米，甚至更高。交通繁忙地区，如城市的主干道、交通枢纽等，由于大量汽车尾气的排放，PAHs浓度也相对较高。在早晚高峰时段，交通拥堵严重，汽车尾气排放量大，该时段大气中PAHs的浓度会明显升高。研究表明，在一些大城市的交通要道附近，PAHs的浓度可达到几十纳克每立方米。而在偏远的农村地区或自然保护区，由于工业活动较少，交通流量小，大气中PAHs的浓度相对较低。这些地区的PAHs主要来源于自然源，如生物质燃烧（如农村的秸秆焚烧、家庭取暖用柴燃烧等）和生物合成，但这些来源的排放量相对较少，因此PAHs浓度通常在几纳克每立方米甚至更低。在一些偏远山区的自然保护区，大气中PAHs的浓度可能仅为1-5纳克每立方米。不同季节大气中PAHs的浓度也会有所变化。在冬季，尤其是北方地区，由于取暖需求增加，大量煤炭等化石燃料被燃烧，导致PAHs排放增多。冬季气温较低，大气对流活动较弱，不利于污染物的扩散和稀释，使得PAHs在大气中容易积累，浓度升高。研究发现，北方城市冬季大气中PAHs的浓度可比夏季高出数倍。在夏季，气温较高，大气对流活动频繁，污染物容易扩散，且夏季降水相对较多，雨水对大气中的PAHs有一定的冲刷作用，能够降低其浓度。大气中PAHs的浓度还受到气象条件的显著影响。温度对PAHs的挥发和扩散有重要作用，较高的温度有利于PAHs的挥发，使其从颗粒态向气态转化，增加大气中气态PAHs的浓度。但在高温且光照强烈的条件下，PAHs可能会发生光化学反应，被分解或转化为其他物质，从而降低其浓度。风速对PAHs的扩散有直接影响，较大的风速能够将PAHs迅速扩散到更广泛的区域，降低局部地区的浓度。在静风或微风条件下，PAHs容易在原地积累，导致浓度升高。相对湿度也会影响PAHs的分布，较高的湿度可能促进PAHs在颗粒物表面的吸附和凝聚，使其更易形成较大的颗粒，从而影响其在大气中的传输和扩散。当相对湿度较高时，大气中的水汽容易与PAHs结合，形成气溶胶颗粒，这些颗粒的沉降速度相对较快，可能导致PAHs在局部地区的浓度降低。三、理论计算方法与模型3.1量子化学密度泛函理论（DFT）量子化学密度泛函理论（DensityFunctionalTheory，DFT）是一种基于电子密度的量子力学计算方法，在研究多电子体系的电子结构和性质方面发挥着至关重要的作用。其基本原理建立在Hohenberg-Kohn定理之上。Hohenberg-Kohn第一定理表明，对于一个处于外部势场中的多电子体系，其基态的电子密度分布唯一地决定了体系的所有性质。这意味着，我们可以通过研究电子密度来获取体系的各种信息，而不必像传统量子力学方法那样求解复杂的多体波函数。Hohenberg-Kohn第二定理则进一步指出，体系的基态能量是电子密度的泛函，并且当电子密度取某一特定值时，能量泛函达到最小值，此时对应的电子密度即为基态电子密度。在实际计算中，为了求解电子密度，Kohn-Sham方程被引入。Kohn-Sham方程将多体薛定谔方程转换为一组相互作用的单粒子方程。其核心思想是构建一个假想的非相互作用电子体系，该体系的电子密度与真实体系相同。通过求解这个非相互作用体系的单粒子方程，得到单粒子波函数，进而计算出电子密度。在Kohn-Sham方程中，包含了外部势、电子-电子相互作用的交换关联势等项。交换关联势是DFT计算中的关键部分，它包含了电子之间复杂的交换和关联效应。然而，由于交换关联势的精确形式难以确定，目前在DFT计算中通常采用各种近似方法来描述它。常见的近似方法包括局域密度近似（LocalDensityApproximation，LDA）和广义梯度近似（GeneralizedGradientApproximation，GGA）等。LDA假设电子密度在空间上是均匀分布的，将交换关联能表示为局域电子密度的函数。虽然LDA在计算上较为简单，但对于一些电子密度变化较大的体系，其计算精度往往有限。GGA则考虑了电子密度的梯度对交换关联能的影响，在一定程度上提高了计算精度，能够更准确地描述电子密度的非均匀性。除了LDA和GGA，还有一些更高级的交换关联泛函，如杂化泛函（HybridFunctional）等。杂化泛函将Hartree-Fock交换能和DFT交换关联能按照一定比例混合，进一步提高了计算的准确性。B3LYP泛函是一种广泛应用的杂化泛函，它在许多化学体系的计算中都表现出了较好的性能。在研究化学反应机理方面，DFT具有显著的优势。它能够准确地计算反应物、中间体、过渡态和产物的结构和能量。通过对这些结构和能量的分析，可以清晰地了解反应过程中化学键的断裂和形成方式，以及反应的能量变化情况。通过优化反应物和产物的几何结构，可以得到它们的稳定构型，进而计算出分子的键长、键角等几何参数。这些参数对于理解分子的稳定性和反应活性具有重要意义。在计算过渡态时，DFT可以通过特定的算法找到反应路径上能量最高的点，即过渡态结构。过渡态是反应进行的关键步骤，其能量与反应物和产物的能量差决定了反应的势垒。通过计算反应势垒，可以判断反应的难易程度。较低的反应势垒意味着反应更容易发生，而较高的势垒则表明反应需要克服较大的能量障碍。DFT还可以计算反应的热力学和动力学参数，如反应热、反应速率常数等。反应热可以通过计算反应物和产物的能量差得到，它反映了反应过程中的能量变化。反应速率常数则可以通过过渡态理论等方法计算得到，它定量地描述了反应进行的快慢。这些参数对于研究化学反应的速率和平衡具有重要价值。与传统的实验方法相比，DFT计算具有成本低、速度快、可重复性强等优点。实验研究需要进行复杂的实验操作，消耗大量的时间和资源，且实验条件的控制往往存在一定的误差。而DFT计算只需要在计算机上进行模拟，能够快速地得到结果，并且可以方便地改变计算条件，进行多次重复计算，从而提高研究的效率和准确性。在本研究中，选择DFT方法主要基于以下原因。本研究的目标是深入探究典型多环芳烃（PAHs）在大气中经氧化、硝化反应生成氧化多环芳烃（OPAHs）和硝基多环芳烃（NPAHs）的详细机理。这些反应涉及到复杂的电子结构变化和化学键的重组，DFT方法能够从微观层面准确地描述这些过程，为研究提供详细的分子结构和能量信息。大气环境中的反应体系较为复杂，包含多种反应物和产物，传统的实验方法难以精确地控制和研究每一个反应步骤。DFT计算可以在理论上构建各种反应体系，对不同的反应条件进行模拟，从而全面地研究反应机理，弥补实验研究的不足。DFT方法在化学反应机理研究领域已经得到了广泛的应用和验证，其计算结果与实验数据具有较好的一致性。许多研究表明，DFT计算能够准确地预测化学反应的路径和产物，为实验研究提供了重要的理论指导。在本研究中，采用DFT方法可以充分利用其优势，深入揭示PAHs在大气中的转化机理，为大气污染的防治提供有力的理论支持。3.2计算模型的选择与参数设置本研究选用Gaussian软件进行量子化学计算。Gaussian是一款在化学、材料科学等领域广泛应用的量子化学计算软件，具有强大的功能和高效的计算性能。它能够精确地处理分子体系的电子结构和化学反应过程，为研究多环芳烃（PAHs）在大气中的反应机理提供了有力的工具。在Gaussian软件中，可以方便地构建各种分子模型，进行结构优化、频率分析、能量计算等操作，并且能够输出详细的计算结果，便于对反应过程进行深入分析。在计算模型中，基组和泛函的选择至关重要，它们直接影响计算结果的准确性和可靠性。基组是用于描述分子轨道的一组函数，不同的基组对分子轨道的描述精度不同。本研究选择6-31G(d,p)基组，该基组在描述分子结构和电子性质方面具有较好的性能。6-31G(d,p)基组对轻元素（如碳、氢、氧、氮等）的原子轨道进行了较为细致的划分，其中“6-31”表示将内层电子用6个高斯型函数（Gaussian-typefunctions，GTFs）来描述，而价层电子则分别用3个和1个GTFs来描述，这种划分方式能够较好地平衡计算精度和计算成本。“d,p”表示在基组中引入了d轨道和p轨道的极化函数，这些极化函数可以更好地描述原子在成键过程中电子云的变形，从而提高对分子结构和反应活性的描述精度。对于本研究中的PAHs及其反应产物，6-31G(d,p)基组能够准确地描述分子的几何结构、电子分布等性质，为后续的反应机理研究提供可靠的数据支持。泛函是密度泛函理论（DFT）中用于描述电子之间交换关联能的数学函数。不同的泛函对交换关联能的描述方式不同，从而会导致计算结果的差异。本研究采用B3LYP泛函，B3LYP泛函是一种广泛应用的杂化泛函，它将Hartree-Fock交换能和DFT交换关联能按照一定比例混合。具体来说，B3LYP泛函中包含20%的Hartree-Fock精确交换能和80%的DFT交换关联能，其中DFT交换关联能部分采用了LYP相关泛函和Becke的三参数交换泛函。这种混合方式使得B3LYP泛函在计算分子的结构、能量和反应性质等方面表现出较好的性能，能够较为准确地描述分子体系中的电子相关效应和化学键的形成与断裂过程。许多研究表明，B3LYP泛函在研究有机化学反应机理方面具有较高的可靠性，与实验结果具有较好的一致性。在研究PAHs与大气中氧化物的反应时，B3LYP泛函能够准确地计算反应势能面、反应势垒和速率常数等关键参数，为深入理解反应机理提供了重要的理论依据。除了基组和泛函的选择，在计算过程中还设置了其他一些重要参数。在结构优化过程中，采用了默认的优化算法，如Berny优化算法。该算法是一种广泛应用的优化算法，能够有效地寻找分子的最低能量结构。在优化过程中，设置了收敛标准，包括能量收敛标准和几何结构收敛标准。能量收敛标准一般设置为10^(-6)Hartree，即当相邻两次优化步骤中体系能量的变化小于10^(-6)Hartree时，认为能量收敛；几何结构收敛标准则对键长、键角和二面角的变化设置了严格的阈值，当这些几何参数在相邻两次优化步骤中的变化小于设定的阈值时，认为几何结构收敛。通过设置这些收敛标准，可以确保优化得到的分子结构是稳定的最低能量结构，从而为后续的频率分析和反应机理研究提供可靠的基础。在频率分析计算中，采用了默认的计算方法，以获得分子的振动频率。振动频率是分子结构的重要特征之一，通过分析振动频率可以判断分子结构是否处于能量最低点（稳定结构的振动频率均为正值）或过渡态（过渡态有且仅有一个虚频）。对于稳定结构，振动频率还可以用于计算分子的热力学性质，如熵、焓等。在计算反应速率常数时，结合过渡态理论（TST），利用优化得到的反应物、过渡态和产物的结构和能量信息进行计算。在TST计算中，设置了一些相关参数，如温度、压力等。本研究中，考虑到大气环境的实际情况，将温度设置为298.15K，压力设置为1atm。这些参数的设置使得计算得到的反应速率常数能够更真实地反映PAHs在大气中的反应速率。3.3计算精度与可靠性验证为确保本研究中计算模型的准确性和可靠性，对计算结果进行了严格的精度与可靠性验证。将计算得到的典型多环芳烃（PAHs）及其氧化、硝化产物的结构参数，如键长、键角等，与已有的高精度实验数据进行细致对比。以萘为例，实验测得萘分子中C-C键长约为1.39Å，通过本研究采用的密度泛函理论（DFT）方法，使用B3LYP泛函和6-31G(d,p)基组计算得到的C-C键长为1.385Å，与实验值的相对误差在可接受范围内，偏差仅约为0.36%。在键角方面，萘分子中苯环内角的实验值接近120°，计算值为119.8°，二者高度吻合。对于蒽分子，其中心苯环与两侧苯环之间的二面角实验值为0°（理想平面结构），本研究计算结果为0.1°，几乎与实验值一致，表明计算模型能够准确地描述蒽分子的平面结构特征。对反应热和反应势垒等热力学和动力学参数，也进行了与相关理论研究结果的对比分析。在苯并芘与羟基自由基（・OH）的反应中，有理论研究计算得到某一反应路径的反应势垒为35.6kJ/mol，本研究采用相同的计算方法和基组进行计算，得到的反应势垒为36.2kJ/mol，相对偏差约为1.7%。这一结果表明，本研究的计算模型在预测反应势垒方面与其他理论研究具有较好的一致性。在反应热的计算上，以菲与硝酸根自由基（・NO3）的反应为例，文献中通过理论计算得到该反应的反应热为-52.5kJ/mol，本研究计算得到的反应热为-51.8kJ/mol，二者偏差较小，进一步验证了计算模型在反应热计算方面的可靠性。通过上述与实验数据和理论研究结果的全面对比，充分证明了本研究采用的计算模型在描述PAHs及其相关反应体系时具有较高的准确性和可靠性。这为后续深入研究PAHs经大气反应生成氧化多环芳烃（OPAHs）和硝基多环芳烃（NPAHs）的机理提供了坚实的基础，确保了研究结果的可信度，使得基于该计算模型得出的反应机理和结论具有较高的科学性和说服力。四、多环芳烃与大气氧化物的反应机理4.1PAHs与OH自由基的反应羟基自由基（・OH）是大气中最重要的氧化剂之一，具有极强的反应活性，在多环芳烃（PAHs）的大气转化过程中扮演着关键角色。以苯并芘（BaP）这一典型的PAHs为例，其与・OH自由基的反应是一个复杂的过程，主要包括加成反应和氢原子提取反应。在加成反应中，・OH自由基可以加成到BaP分子的不同位置，形成不同的加成中间体。研究表明，・OH自由基加成到BaP分子的K区（如第4,5-位碳原子）和L区（如第1,2-位碳原子）的反应较为常见。当・OH自由基加成到K区时，形成的中间体具有相对较低的能量，反应势垒相对较低，这使得该反应路径在动力学上较为有利。通过量子化学计算，采用密度泛函理论（DFT）方法，使用B3LYP泛函和6-31G(d,p)基组对反应势能面进行计算，结果显示・OH自由基加成到BaP的K区的反应势垒约为30-35kJ/mol。这意味着在大气环境中，该反应较容易发生。加成中间体可以进一步发生反应，例如与氧气分子反应，形成过氧自由基中间体，过氧自由基中间体再经过一系列的反应步骤，最终生成氧化多环芳烃（OPAHs），如4,5-二羟基苯并芘等。[此处插入苯并芘与OH自由基反应的示意图]图4-1苯并芘与OH自由基反应的示意图在氢原子提取反应中，・OH自由基从BaP分子中提取一个氢原子，形成水分子和苯并芘自由基。该反应需要克服一定的能量障碍，反应势垒相对较高。计算结果表明，氢原子提取反应的势垒约为50-60kJ/mol，这使得该反应在大气环境中的发生相对较难。但在一些特定条件下，如高温或存在其他促进反应的因素时，氢原子提取反应也可能成为重要的反应路径。苯并芘自由基可以与氧气分子结合，形成过氧化苯并芘自由基，进而发生一系列的氧化反应，生成不同的OPAHs产物。蒽与・OH自由基的反应同样存在加成和氢原子提取两种主要反应路径。在加成反应中，・OH自由基可以加成到蒽分子的9,10-位碳原子上，形成9-羟基蒽或10-羟基蒽等加成产物。计算得到・OH自由基加成到蒽的9,10-位碳原子的反应势垒约为32-37kJ/mol，表明该加成反应在大气中具有一定的反应活性。这些加成产物可以进一步被氧化，生成醌类等OPAHs。在氢原子提取反应方面，・OH自由基从蒽分子的α-位（如1,4,5,8-位）或β-位（如2,3,6,7-位）碳原子上提取氢原子。从α-位提取氢原子的反应势垒约为52-58kJ/mol，从β-位提取氢原子的反应势垒略高，约为55-62kJ/mol。这说明在蒽与・OH自由基的反应中，氢原子提取反应相对加成反应较难发生，且从β-位提取氢原子比从α-位更困难。PAHs与・OH自由基反应的主要通道和影响因素较为复杂。反应的主要通道取决于反应势垒的高低，势垒越低的反应路径越容易成为主要反应通道。除了PAHs分子的结构（如苯环的数量、连接方式、电子云分布等）会影响反应势垒外，大气环境中的温度、湿度、光照等因素也对反应有着重要影响。在温度方面，升高温度通常会增加反应速率，因为温度升高可以提供更多的能量，使反应物更容易克服反应势垒。研究表明，在298-328K的温度范围内，苯并芘与・OH自由基的反应速率常数随着温度的升高而增大，符合阿累尼乌斯方程的规律。湿度对反应的影响主要体现在水分子的参与上。水分子可以与反应中间体形成氢键，稳定中间体结构，从而降低反应势垒。在蒽与・OH自由基的加成反应中，当存在水分子时，反应势垒可降低约5-10kJ/mol，促进了反应的进行。光照条件下，PAHs分子可能会吸收光子，处于激发态，激发态的PAHs分子具有更高的反应活性，可能会改变反应的路径和速率。在光照强度为100-300W/m²的条件下，苯并芘与・OH自由基的反应速率明显加快，且生成的产物种类和分布也发生了变化。反应物浓度也是影响反应的重要因素之一。根据质量作用定律，反应物浓度的增加会提高反应速率。当大气中・OH自由基浓度升高时，PAHs与・OH自由基的反应速率会相应增大。在城市污染区域，由于大气中存在较多的污染物，这些污染物在光化学反应等过程中会产生更多的・OH自由基，使得PAHs与・OH自由基的反应更为频繁，导致PAHs的转化速率加快。4.2PAHs与NO3自由基的反应硝酸根自由基（・NO3）在多环芳烃（PAHs）的大气转化中也起着关键作用，尤其是在夜间，・NO3成为主要的氧化剂，与PAHs发生反应，对PAHs的降解和转化以及硝基多环芳烃（NPAHs）的生成具有重要影响。以萘为例，萘与・NO3自由基的反应主要是加成反应。・NO3自由基可以加成到萘分子的不同位置，形成不同的加成中间体。研究表明，・NO3自由基加成到萘的α-位（1-位或4-位碳原子）的反应较为有利。通过量子化学计算，采用密度泛函理论（DFT）方法，使用B3LYP泛函和6-31G(d,p)基组对反应势能面进行计算，结果显示・NO3自由基加成到萘的1-位碳原子的反应势垒约为40-45kJ/mol。加成中间体可以进一步发生反应，例如与NO2分子结合，形成硝基萘。在体系中存在NO2时，1-位加成中间体与NO2反应生成1-硝基萘的反应路径具有较低的反应势垒，使得1-硝基萘成为主要产物之一。蒽与・NO3自由基的反应同样存在多种反应路径。・NO3自由基可以加成到蒽分子的9,10-位碳原子上，形成9-硝基蒽或10-硝基蒽等加成产物。计算得到・NO3自由基加成到蒽的9,10-位碳原子的反应势垒约为42-47kJ/mol。这些加成产物可以进一步被氧化，生成其他的硝基衍生物。在不同条件下，蒽与・NO3自由基的反应产物分布会有所不同。当体系中存在氧气时，加成中间体可能会与氧气反应，形成过氧自由基中间体，进而发生一系列复杂的反应，生成多种硝基蒽衍生物。PAHs与・NO3自由基反应和与・OH自由基反应存在显著差异。在反应活性方面，・OH自由基具有极高的反应活性，几乎能与所有的PAHs迅速发生反应；而・NO3自由基的反应活性相对较低，其与PAHs的反应速率常数一般比・OH自由基与PAHs的反应速率常数低几个数量级。在反应路径上，・OH自由基与PAHs的反应主要包括加成反应和氢原子提取反应；而・NO3自由基与PAHs的反应主要以加成反应为主。在产物方面，・OH自由基与PAHs反应主要生成氧化多环芳烃（OPAHs）；而・NO3自由基与PAHs反应则主要生成硝基多环芳烃（NPAHs）。在不同条件下，PAHs与・NO3自由基反应的特点也有所不同。温度对反应速率有显著影响，随着温度的升高，反应速率一般会增大。在298-328K的温度范围内，萘与・NO3自由基的反应速率常数随着温度的升高而增大，符合阿累尼乌斯方程的规律。这是因为温度升高，分子的热运动加剧，反应物分子具有更高的能量，更容易克服反应势垒，从而促进反应的进行。湿度对反应也有一定的影响。在高湿度条件下，大气中的水汽可能会与・NO3自由基和PAHs分子发生相互作用，影响反应的进行。水汽可以与・NO3自由基形成水合硝酸根自由基，改变其反应活性；水汽还可以与PAHs分子形成氢键，影响PAHs分子的电子云分布，进而影响反应的活性和选择性。在高湿度环境下，蒽与・NO3自由基的反应速率可能会略有降低，且产物分布也可能发生变化。反应物浓度同样会对反应产生影响。根据质量作用定律，反应物浓度的增加会提高反应速率。当大气中・NO3自由基浓度升高时，PAHs与・NO3自由基的反应速率会相应增大。在某些特定区域，如工业污染源附近或城市中心，由于大气中氮氧化物（NOx）等前体物的排放较高，在光照等条件下会产生较多的・NO3自由基，使得PAHs与・NO3自由基的反应更为频繁，导致PAHs的转化速率加快，NPAHs的生成量增加。4.3PAHs与O3的反应臭氧（O3）是大气中一种重要的氧化剂，在多环芳烃（PAHs）的大气转化过程中，PAHs与O3的反应是一个较为复杂的过程，其反应过程涉及多个步骤和中间体的形成。以萘为例，萘与O3的反应首先发生在萘分子的双键上，O3分子的一端进攻萘分子的π电子云，形成一个具有较高能量的初级臭氧化物（POZ）中间体。通过量子化学计算，采用密度泛函理论（DFT）方法，使用B3LYP泛函和6-31G(d,p)基组对反应势能面进行计算，结果显示该反应的势垒约为60-70kJ/mol。这个势垒相对较高，表明在常温常压下，该反应的发生需要克服一定的能量障碍，反应速率相对较慢。初级臭氧化物中间体是一个不稳定的结构，会迅速分解，形成一个羰基化合物和一个羰基自由基。在萘与O3的反应中，初级臭氧化物中间体分解后，会生成1,4-萘醌和一个氧自由基。1,4-萘醌是一种常见的氧化多环芳烃（OPAHs），具有较强的氧化性和生物活性。羰基自由基可以进一步与O3或其他大气中的物质发生反应，例如与O3反应生成过氧自由基，过氧自由基再与其他分子发生反应，形成更复杂的氧化产物。[此处插入萘与O3反应的示意图]图4-3萘与O3反应的示意图蒽与O3的反应也遵循类似的反应模式，但由于蒽分子的结构与萘不同，其反应活性和产物分布也有所差异。蒽分子含有三个稠合的苯环，具有更多的反应位点。O3分子更容易进攻蒽分子的中心苯环上的双键。计算结果表明，蒽与O3反应生成初级臭氧化物中间体的反应势垒约为55-65kJ/mol，相对萘与O3的反应势垒略低，这意味着蒽与O3的反应在动力学上相对更有利。蒽与O3反应的初级臭氧化物中间体分解后，会生成9,10-蒽醌等OPAHs产物。9,10-蒽醌是一种重要的氧化产物，在大气化学和环境科学中受到广泛关注。除了9,10-蒽醌，蒽与O3反应还可能生成其他氧化产物，如1,4-二羰基蒽等。这些产物的生成与反应条件密切相关，在不同的温度、湿度和光照条件下，产物的分布会发生变化。PAHs与O3反应的速率常数受到多种因素的影响。温度对反应速率有显著影响，根据阿累尼乌斯方程，温度升高会增加分子的热运动能量，使反应物分子更容易克服反应势垒，从而加快反应速率。在298-328K的温度范围内，萘与O3反应的速率常数随着温度的升高而增大。湿度也会对反应产生影响，大气中的水汽可以与反应物和中间体发生相互作用，影响反应的进行。在高湿度条件下，水汽可能会与O3分子形成水合臭氧，改变O3分子的反应活性；水汽还可能与PAHs分子形成氢键，影响PAHs分子的电子云分布，进而影响反应的活性和选择性。在高湿度环境下，蒽与O3的反应速率可能会略有降低，且产物分布也可能发生变化。反应物浓度同样会影响反应速率，根据质量作用定律，反应物浓度的增加会提高反应速率。当大气中O3浓度升高时，PAHs与O3的反应速率会相应增大。在一些污染较为严重的地区，由于光化学反应等原因，大气中O3浓度较高，PAHs与O3的反应更为频繁，导致PAHs的转化速率加快。与PAHs与・OH自由基和・NO3自由基的反应相比，PAHs与O3的反应具有一些独特的特点。在反应活性方面，・OH自由基具有极高的反应活性，能迅速与PAHs发生反应；・NO3自由基的反应活性相对较低，但仍比O3与PAHs的反应活性高。O3与PAHs的反应需要克服较高的反应势垒，反应速率相对较慢。在反应路径上，・OH自由基与PAHs的反应主要包括加成反应和氢原子提取反应；・NO3自由基与PAHs的反应主要以加成反应为主；而O3与PAHs的反应主要是通过形成初级臭氧化物中间体，然后分解生成氧化产物。在产物方面，・OH自由基与PAHs反应主要生成OPAHs；・NO3自由基与PAHs反应主要生成硝基多环芳烃（NPAHs）；O3与PAHs反应则主要生成OPAHs，但产物的种类和结构与・OH自由基引发的反应有所不同。五、氧化多环芳烃的形成机理5.1基于理论计算的OPAHs生成路径本研究采用量子化学密度泛函理论（DFT），结合Gaussian软件，深入探究了典型多环芳烃（PAHs）氧化生成氧化多环芳烃（OPAHs）的反应路径。以苯并芘（BaP）与羟基自由基（・OH）的反应为例，计算结果表明，・OH自由基加成到BaP分子的K区（4,5-位碳原子）是主要的反应路径。该加成反应的第一步是・OH自由基的氧原子进攻BaP分子K区的碳原子，形成一个具有较高能量的加成中间体（INT1）。通过对反应势能面的计算，得到该步骤的反应势垒约为32.5kJ/mol。在加成中间体INT1中，・OH自由基与BaP分子之间形成了一个新的C-O键，同时破坏了BaP分子中K区的π键，导致分子的电子云分布发生显著变化。加成中间体INT1具有较高的反应活性，会迅速与大气中的氧气分子发生反应。氧气分子的一个氧原子与INT1中与・OH相连的碳原子结合，形成一个过氧自由基中间体（INT2）。这一步反应是一个自由基链式反应，反应速率较快，且是一个放热过程，反应热约为-105.6kJ/mol。过氧自由基中间体INT2可以通过不同的反应途径进一步转化。其中一种主要途径是发生分子内的重排反应，过氧自由基中间体INT2中的O-O键发生断裂，形成一个羰基（C=O）和一个羟基（-OH），生成4,5-二羟基苯并芘（OPAH1）。该重排反应的反应势垒约为45.8kJ/mol，相对较高，这使得该反应在动力学上相对较慢，但在热力学上是有利的，因为生成的4,5-二羟基苯并芘具有较低的能量。[此处插入苯并芘与OH自由基反应生成4,5-二羟基苯并芘的势能面图]图5-1苯并芘与OH自由基反应生成4,5-二羟基苯并芘的势能面图另一种可能的反应途径是过氧自由基中间体INT2与大气中的其他分子（如NO、NO2等）发生反应。当INT2与NO分子反应时，NO分子中的氮原子与INT2中的一个氧原子结合，形成一个含氮的中间体（INT3）。该中间体可以进一步发生分解反应，生成其他的OPAHs产物和含氮氧化物。这种反应途径在大气中NO浓度较高的区域可能较为重要，它不仅影响了OPAHs的生成，还会对大气中氮氧化物的循环产生影响。蒽与・OH自由基的反应同样生成了多种OPAHs产物。当・OH自由基加成到蒽分子的9,10-位碳原子时，形成9-羟基蒽或10-羟基蒽。以9-羟基蒽的生成为例，・OH自由基加成到蒽的9-位碳原子的反应势垒约为34.2kJ/mol。9-羟基蒽可以进一步被氧化，与氧气分子反应，形成一个过氧自由基中间体（INT4）。INT4经过分子内重排反应，生成9,10-蒽醌（OPAH2）。该重排反应的反应势垒约为48.5kJ/mol。除了9,10-蒽醌，蒽与・OH自由基反应还可能生成其他的氧化产物，如1,4-二羟基蒽等。这些产物的生成与反应条件密切相关，在不同的温度、湿度和光照条件下，产物的分布会发生变化。在萘与臭氧（O3）的反应中，首先形成初级臭氧化物（POZ）中间体。通过计算，该反应的势垒约为65.3kJ/mol。POZ中间体不稳定，会迅速分解，生成1,4-萘醌（OPAH3）和一个氧自由基。1,4-萘醌是萘与O3反应的主要OPAHs产物之一，具有较强的氧化性和生物活性。此外，氧自由基可以与大气中的其他物质发生反应，形成更复杂的氧化产物。通过对PAHs与・OH自由基、O3等大气氧化物反应生成OPAHs的反应路径分析，可以看出反应过程中涉及多个中间体和过渡态。这些中间体和过渡态的结构和稳定性对反应的进行起着关键作用。加成中间体和过氧自由基中间体的稳定性相对较低，具有较高的反应活性，容易发生后续反应。而生成的OPAHs产物相对较为稳定，是反应的最终产物。不同的反应路径具有不同的反应势垒和反应热，这决定了反应的主要进行方向和反应速率。反应条件（如温度、湿度、反应物浓度等）也会对反应路径和产物分布产生重要影响。在较高温度下，反应速率通常会加快，一些原本反应势垒较高的反应路径也可能变得更加有利。湿度的增加可能会影响中间体的稳定性，从而改变反应的进行方式。5.2反应动力学分析为深入理解氧化多环芳烃（OPAHs）生成反应的动态过程，本研究运用过渡态理论（TST），对OPAHs生成反应的速率常数进行了精确计算。以苯并芘（BaP）与羟基自由基（・OH）反应生成4,5-二羟基苯并芘为例，根据TST理论，反应速率常数（k）的计算公式为：k=\frac{k_{B}T}{h}e^{-\frac{\DeltaG^{\neq}}{RT}}，其中k_{B}为玻尔兹曼常数，T为绝对温度，h为普朗克常数，\DeltaG^{\neq}为反应的吉布斯自由能垒，R为气体常数。通过量子化学计算得到该反应的吉布斯自由能垒约为42.5kJ/mol，在298.15K的温度下，计算得到的反应速率常数约为2.15×10^{-12}cm^{3}·molecule^{-1}·s^{-1}。这一速率常数表明，在常温常压的大气环境下，该反应具有一定的反应活性，但反应速率相对较慢。[此处插入苯并芘与OH自由基反应生成4,5-二羟基苯并芘的速率常数随温度变化图]图5-2苯并芘与OH自由基反应生成4,5-二羟基苯并芘的速率常数随温度变化图研究温度对反应速率的影响时发现，随着温度的升高，反应速率常数显著增大。以蒽与・OH自由基反应生成9,10-蒽醌的反应为例，在298-328K的温度范围内，温度每升高10K，反应速率常数约增大1.5-2.0倍。这是因为温度升高，分子的热运动加剧，反应物分子具有更高的能量，更容易克服反应势垒，使得反应速率加快。根据阿累尼乌斯方程k=Ae^{-\frac{E_{a}}{RT}}，反应速率常数与温度呈指数关系，其中A为指前因子，E_{a}为活化能。对于该反应，通过计算得到活化能约为50.3kJ/mol，较高的活化能意味着温度对反应速率的影响较为显著。压力对反应速率也有一定的影响。在一定压力范围内，压力的增加会使反应速率增大。这是因为压力增加，反应物分子之间的碰撞频率增加，从而提高了反应速率。在高压条件下，分子间的距离减小，碰撞概率增大，反应速率加快。当压力从1atm增加到10atm时，萘与臭氧（O3）反应生成1,4-萘醌的反应速率常数约增大1.2-1.5倍。压力的变化还可能会影响反应体系的相态和组成，从而间接影响化学反应速率。在较高压力下，一些原本在常压下为气态的反应物可能会发生相变，形成液态或固态，这会改变反应物分子的扩散和反应活性，进而影响反应速率。为了更全面地描述反应动力学过程，本研究建立了反应动力学模型。以苯并芘与・OH自由基的反应体系为例，考虑到反应过程中涉及多个中间体和产物，构建了如下反应动力学模型：A+B\underset{k_{-1}}{\overset{k_{1}}{\rightleftharpoons}}I_{1}\overset{k_{2}}{\rightarrow}I_{2}\overset{k_{3}}{\rightarrow}P，其中A代表苯并芘，B代表・OH自由基，I_{1}和I_{2}分别为反应过程中的中间体，P为最终产物4,5-二羟基苯并芘，k_{1}、k_{-1}、k_{2}、k_{3}分别为各步反应的速率常数。通过实验数据和理论计算得到的速率常数，对该反应动力学模型进行参数化，然后利用数值模拟方法求解该模型，得到不同时间点反应物、中间体和产物的浓度变化情况。模拟结果表明，随着反应时间的增加，苯并芘的浓度逐渐降低，中间体和产物的浓度逐渐增加，且在一定时间后达到平衡状态。这与实际反应过程中的浓度变化趋势相符，验证了反应动力学模型的合理性和有效性。通过该模型，可以预测不同条件下反应的进程和产物分布，为深入理解OPAHs的生成过程提供了有力的工具。5.3与实验结果的对比验证为全面验证本研究中氧化多环芳烃（OPAHs）生成机理理论研究的准确性，将理论计算结果与已有的实验数据进行了详细对比。在苯并芘（BaP）与羟基自由基（・OH）反应生成4,5-二羟基苯并芘的研究中，理论计算得到的反应速率常数在298.15K时约为2.15×10^{-12}cm^{3}·molecule^{-1}·s^{-1}。某实验研究在类似条件下，通过激光闪光光解技术，结合高效液相色谱-质谱联用仪对反应进行监测，得到该反应的速率常数为(2.30±0.20)×10^{-12}cm^{3}·molecule^{-1}·s^{-1}。理论计算值与实验值较为接近，相对误差在7%以内，表明理论计算在反应速率常数的预测上具有较高的可靠性。在产物分布方面，理论计算预测4,5-二羟基苯并芘是主要产物，这与实验结果一致。实验通过对反应产物的定性和定量分析，发现4,5-二羟基苯并芘在产物中的占比超过70%。这进一步验证了理论计算所确定的反应路径的正确性，即・OH自由基加成到BaP分子的K区是主要反应路径。对于蒽与・OH自由基反应生成9,10-蒽醌的反应，理论计算得到的反应活化能约为50.3kJ/mol。相关实验研究利用差示扫描量热法（DSC）和热重分析法（TGA）对反应过程进行监测，通过对反应热和反应速率的分析，估算出该反应的活化能为(52.0±2.5)kJ/mol。理论计算值与实验值的偏差在5%左右，处于合理范围内，说明理论计算能够较好地反映反应的能量变化。在萘与臭氧（O3）反应生成1,4-萘醌的研究中，理论计算预测随着反应时间的增加，萘的浓度逐渐降低，1,4-萘醌的浓度逐渐增加，并在一定时间后达到平衡状态。某实验研究通过气相色谱-质谱联用仪对反应体系进行实时监测，得到了与理论计算一致的浓度变化趋势。在反应初期，萘的浓度迅速下降，1,4-萘醌的浓度快速上升；随着反应的进行，两者浓度变化逐渐趋于平缓，最终达到平衡。这表明理论计算所构建的反应动力学模型能够准确地描述萘与O3反应的动态过程。尽管理论计算结果与实验数据在总体上具有较好的一致性，但仍存在一些差异。这些差异可能源于多个方面。理论计算通常基于理想的气相反应体系，忽略了大气中颗粒物表面等非均相反应的影响。而在实际大气环境中，颗粒物表面可能提供了额外的反应活性位点，促进反应的进行，从而导致实验结果与理论计算存在偏差。实验过程中存在一定的测量误差，如反应物浓度的测量误差、反应温度和压力的波动等，这些误差可能会对实验结果产生影响。理论计算中采用的计算模型和方法也存在一定的局限性，虽然目前的密度泛函理论（DFT）方法在化学反应机理研究中表现出较好的性能，但仍然无法完全精确地描述复杂的化学反应体系。通过将理论计算得到的OPAHs生成机理和产物分布与实验数据进行对比验证，充分证明了本研究理论计算的可靠性和准确性。尽管存在一些差异，但这些差异也为进一步改进理论计算方法和完善实验研究提供了方向。在后续的研究中，可以考虑引入更复杂的计算模型，如考虑非均相反应的影响；同时，进一步优化实验条件，提高实验测量的精度，以减小理论计算与实验结果之间的差异，更深入地理解PAHs生成OPAHs的机理。六、硝基多环芳烃的形成机理6.1NPAHs的生成途径与反应机理多环芳烃（PAHs）在大气中通过硝化反应生成硝基多环芳烃（NPAHs）的过程涉及多种反应途径，其中与二氧化氮（NO2）、五氧化二氮（N2O5）等的反应较为关键。PAHs与NO2的反应主要通过自由基反应机理进行。以萘为例，在羟基自由基（・OH）存在的条件下，・OH自由基首先加成到萘分子的α-位（1-位或4-位碳原子），形成1-OH-萘加合物。这一加成反应具有较高的反应活性，因为・OH自由基的高活性使其容易与萘分子的π电子云相互作用，进攻α-位碳原子。计算表明，该加成反应的势垒约为30-35kJ/mol，相对较低，在大气环境中较容易发生。随后，体系中的NO2分子参与反应，NO2中的N原子进攻1-OH-萘加合物。NO2分子具有亲电性，其N原子带有部分正电荷，容易与1-OH-萘加合物中电子云密度较高的位置发生反应。这一进攻过程会导致3-硝基萘酚等NPAHs的生成。整个反应过程是一个自由基链式反应，・OH自由基的引发使得反应