聚乙烯微塑料中弱键的多尺度识别：键离解能与等转化率动力学的桥接

《Microplastics》：Multiscale Identification of Weak Links in Polyethylene Microplastics: Bridging Bond Dissociation Energy and Isoconversional Kinetics

　　

聚乙烯（PE）微塑料是环境持久性污染物，其渐进式碎裂表明降解可能并非由均匀的主链稳定性所控制。本研究将密度泛函理论（DFT）与非等温热重分析（TGA）相结合，从分子和动力学互补视角评估PE降解。在M06-2X/def2-TZVP水平下优化了C₉₀H₁₈₂聚乙烯齐聚物，并沿链计算了位置分辨的C–H和C–C键离解能（BDE）。C–H键表现出相对较高且均匀的稳定性，而C–C主链则显示出较低的离解能，并在中心区域出现局部能量凹陷。最小C–C BDE出现在C₄₄位置，数值为85.73 kcal·mol^?1，从而在有限的全反式衍生齐聚物中识别出一个模型特异性的低BDE区域，该区域在所评估的计算条件下可能有利于主链断裂。在氮气气氛下、升温速率为5、10和20 °C min^?1的热重分析显示出一个主要降解事件，微分热重（DTG）峰温随升温速率增加从462.6 °C移至494.6 °C。Flynn–Wall–Ozawa（FWO）和Kissinger–Akahira–Sunose（KAS）分析揭示，表观活化能在低转化率时约为170–175 kJ·mol^?1，在高转化率时逐渐增加至280–285 kJ·mol^?1。尽管BDE与表观活化能并不直接等价，但二者的联合解释支持一种异质降解模型，即PE碎裂优先在局部的低BDE C–C环境中引发，随后推进至常规主链断裂和二次降解反应。

聚乙烯（PE）微塑料是环境中广泛存在的持久性污染物。PE因低成本、耐化学性、易加工和长期力学耐久性而被大量生产，但其废弃后经光氧化、热老化、机械磨损和化学风化等过程，可逐步转化为次生微塑料和纳米塑料。传统降解机理常将PE内部C–C键视为能量等价，但局部构象弛豫、链端效应以及自由基碎片稳定化差异可能导致键强度随位置变化。以往计算多采用C3–C10短链烷烃模型，平均内部C–C键离解能（BDE，即均裂特定化学键所需的内在电子能量）约为87–89 kcal·mol^?1，但此类模型难以评估链内局域能量变化。因此，研究人员提出一个更具体的机制问题：PE碎裂的引发是由统计等价的骨架键控制，还是由具有较低内在键稳定性的局域分子环境控制。为检验该假设，研究人员将密度泛函理论（DFT）与非等温热重分析（TGA）相结合，从分子和动力学两个互补尺度分析PE降解，并在此基础上提出异质主链降解模型。该研究突破短链模型的局限，为理解微塑料碎裂的早期分子起源提供了新视角。论文发表在《Microplastics》。

研究人员开展的研究主要包括：构建含90个碳原子的线性C₉₀H₁₈₂齐聚物，在M06-2X/def2-TZVP水平优化并计算位置分辨的C–H和C–C BDE；对采集自哥伦比亚某聚乙烯生产厂的原生高密度聚乙烯（HDPE）微塑料样品进行拉曼光谱和差示扫描量热（DSC）表征；在氮气气氛下以5、10、20 °C min^?1进行非等温TGA，采用Flynn–Wall–Ozawa（FWO）和Kissinger–Akahira–Sunose（KAS）等转化率方法求取转化率依赖性表观活化能E_α；并对BDE数据进行分组统计与机制关联分析。

研究结果分述如下：

3.1 Optimized Polyethylene Oligomer and Structural Basis for BDE Mapping（优化聚乙烯齐聚物及BDE作图的结构基础）：通过M06-2X/def2-TZVP优化，C₉₀H₁₈₂齐聚物保持扩展的全反式锯齿构象，为位置分辨BDE比较提供了受控分子框架。该模型未再现实际PE的折叠、缠结、支化等形态复杂性，因此仅用于检测有限类PE主链内的能量异质性。

3.2 Spatial Distribution of C–H and C–C Bond Dissociation Energies（C–H和C–C键离解能的空间分布）：位置分辨计算显示，内部C–H键BDE大多集中在约101.4 kcal·mol^?1，仅链端升至约104.3 kcal·mol^?1；C–C键则大部分接近89.2–89.5 kcal·mol^?1的骨架基线，但在链中央C₄₁、C₄₄、C₅₀出现局域能量凹陷，其中C₄₄最低，为85.73 kcal·mol^?1，比规则C–C骨架均值低约3.7 kcal·mol^?1。这表明所选齐聚物并非完全能量均匀，C–C裂解比C–H夺取更可能成为主链断裂的直接分子事件。

3.3 Statistical Evaluation of the Localized Low-BDE C–C Region（局部低BDE C–C区域的统计评估）：将C–C位点分为规则组与低BDE组，规则C–C均值为89.44 kcal·mol^?1，低BDE组均值为86.88 kcal·mol^?1，平均差为2.55 kcal·mol^?1，差异具有统计学意义（p=0.048）。但低BDE组仅含3个位点，该统计被解释为对空间BDE趋势的定量支持，而非普遍弱键行为的决定性证据。

3.4 Benchmarking and Mechanistic Relevance of the DFT-Derived BDE Values（DFT所得BDE值的基准比较与机制意义）：所得C–C和C–H BDE值与文献中PE相关烃类模型的基准范围一致，表明计算方案化学合理。与短链模型相比，扩展齐聚物的主要贡献在于揭示了低BDE域的局域位置，而非仅提供平均BDE值。

3.5 Physicochemical Characterization of the Virgin HDPE Microplastic Sample（原生HDPE微塑料样品的理化表征）：拉曼光谱在约1130、1440、2867 cm^?1处呈现PE特征峰，分别对应骨架C–C伸缩、CH₂变形/剪式振动和亚甲基C–H伸缩；DSC测得熔融峰T_m=132.6 °C、结晶峰T_c=116.3 °C，确认样品为半结晶原生HDPE。

3.6 Thermogravimetric Degradation Behavior of PE Microplastics（PE微塑料的热重降解行为）：TGA/DTG显示PE在氮气下以单一主要质量损失事件降解，DTG峰温随升温速率升高从462.6 °C移至494.6 °C。该单峰不代表单一基元反应，而是自由基引发、主链断裂、氢转移、β-断裂、重组、二次裂化和挥发等多个过程的宏观叠加。

3.7 Conversion-Dependent Activation Energy from Isoconversional TGA（等转化率TGA得到的转化率依赖性活化能）：FWO和KAS方法在α=0.1–0.9范围内均显示E_α随转化率单调上升，从低转化率约170–175 kJ·mol^?1升至高转化率280–285 kJ·mol^?1。两种方法结果一致，表明PE降解不受恒定能垒控制。

3.8 Mechanistic Relationship Between BDE and Apparent Activation Energy（BDE与表观活化能之间的机制关系）：BDE与E_α并非数值等价，前者是分子模型中断键的内在电子能量，后者是块体样品热失重的表观动力学壁垒。按降解阶段关联：低转化率阶段对应低BDE C–C环境参与引发；中转化率阶段对应规则内部C–C键和自由基传播；高转化率阶段对应二次降解、残留碎片分解和重组等过程。

综合讨论指出，该研究支持异质主链降解模型：PE碎片化优先在局部低BDE C–C域引发，随后推进至规则骨架断裂和二次反应。但C₉₀H₁₈₂模型仅为单一全反式优化构象，不含支化、折叠、缠结、氧化缺陷等，且低BDE组仅含三个位点，故C₄₄不应被外推为真实PE中的普适断裂位点。研究结论可翻译为：本研究基于PE主链不应被视为能量均匀C–C键序列的假设，通过位置分辨DFT键离解能与转化率依赖性等转化率动力学结合，建立了联系分子键稳定性与宏观热降解行为的多尺度框架。C–H键离解能较高且均匀，C–C键在链中央存在局域能量凹陷，最低值出现在C₄₄。热重分析显示单一主要失重事件，等转化率分析显示E_α随转化率升高而增大。BDE与E_α的机制关联表明，PE微塑料降解符合异质模型：碎片化优先在弱C–C键域引发，再扩展至常规骨架断裂和二次自由基过程。该发现修正了均匀主链的传统认识，为微塑料早期碎裂的分子机制提供了依据。