聚苯乙烯降解的分子能量学与非等温动力学：一项集成低聚体密度泛函理论–热重分析研究

《Microplastics》：Molecular Energetics and Non-Isothermal Kinetics of Polystyrene Degradation: An Integrated Oligomeric DFT–TGA Study

　　

聚苯乙烯（PS）热降解涉及局域分子键断裂事件，与宏观热分析获得的表观动力学参数并不直接等价。研究人员采用有限氢封端PS低聚体模型，在M06-2X/LANL2DZ水平进行密度泛函理论（DFT）计算，同时独立地在氮气下对该PS样品进行非等温热重分析（TGA）。计算分析考察了前沿分子轨道分布及与均裂C–C断裂和自由基介导β-断裂相关的位点特异性热力学描述符。计算所得HOMO–LUMO能隙为742.62 kJ mol?1，表明所选低聚体模型内轨道能分离较大，而前沿轨道在芳香区与苄位区的定域化揭示了空间异质电子分布。均裂C–C断裂的键解离能介于414.09至481.24 kJ mol?1，显示自由基生成的能垒需求依赖于局部分子环境。所选自由基β-断裂反应的吉布斯自由能变介于55.44至189.41 kJ mol?1，这些量属模型依赖的反应热力学，因未计算过渡态故不可解释为活化能垒。热重分析显示Tonset与Tmax随升温速率升高而系统性升高，符合非等温下动力学延迟与热滞后效应。Kissinger法给出全局表观活化能186.61 kJ mol?1，而残基质量校正后的Flynn–Wall–Ozawa（FWO）与Kissinger–Akahira–Sunose（KAS）法在α=0.05–0.95区间平均表观活化能分别为180.81与178.49 kJ mol?1；同区间内FWO表观活化能从143.10升至221.71 kJ mol?1，KAS从140.10升至220.23 kJ mol?1，表明宏观降解响应演化且高转化区不确定性增大。计算与实验数据被解释为互补但非等价的尺度依赖描述：DFT比较有限孤立低聚体内选定分子反应相对热力学，TGA表征凝聚态聚合物样品全局表观动力学行为。分子反应能与TGA表观活化能间未建立直接数值对应，亦未将单一断裂反应指派至特定转化区间。将结果外推至高相对分子质量、多分散、含添加剂、交联或环境老化PS微塑料须谨慎。

研究背景方面，聚苯乙烯（Polystyrene，PS）因产量大、用途广且天然降解抗性强，废弃后在环境中持续累积为微纳塑料污染物，同时其热解可回收苯乙烯富集液与芳烃具有资源价值。目前PS热降解普遍被描述为自由基介导的多步过程（引发、传播、解聚、β-断裂、氢转移及二次反应），但宏观热重分析（Thermogravimetric Analysis，TGA）给出的表观活化能（Apparent Activation Energy，Ea）是自由基引发、链传播、挥发、传热与二次反应的综合描述符，并不等价于单键均裂能垒；而分子尺度的密度泛函理论（Density Functional Theory，DFT）可给出位点特异性键解离能（Bond Dissociation Energy，BDE）与吉布斯自由能变（ΔG），却常脱离凝聚态多步网络。两尺度间长期缺乏严谨的非数值等价连接，单一常数能垒无法刻画低、中、高转化区主导机制切换，这正是本研究切入点。研究人员在《Microplastics》发表该文，旨在以互补但非等价的方式集成分子热力学与块体非等温动力学，明确PS降解不能由统一能参描述。

关键技术方法上，研究人员构建含十个苯乙烯重复单元的氢封端等规PS十聚体（C₈₀H₈₂），采用Gaussian 16以M06-2X/LANL2DZ在气相873 K、1 atm下做几何优化与频率校正，获得均裂C–C断裂BDE及自由基β-断裂反应焓/ΔG，并作HOMO–LUMO分析；实验侧取机械研磨过200目筛（<75 μm）的未使用反应器级 virgin PS微塑料馏分，真空65℃干燥24 h，以Discovery TGA 55在氮气流60 mL min^?1下以5、10、15、20、25 ℃ min^?1升温至600℃，由TGA/DTG提取T_onset、T_max与转化度α，再用Kissinger、Flynn–Wall–Ozawa（FWO）、Kissinger–Akahira–Sunose（KAS）三种非等温法求全局及α=0.05–0.95依赖的表观活化能，全程区分DFT热力学量与TGA表观动力学量。

研究结果部分，第三节“优化结构与前沿分子轨道分布”中，优化十聚体呈非线性构象，HOMO能?659.21 kJ mol^?1、LUMO能83.41 kJ mol^?1，HOMO–LUMO能隙742.62 kJ mol^?1，轨道定域于芳香/苄位区，说明PS具全局电子稳定但存在电子异质弱联区，支撑位点非均匀断裂而非直接替代BDE。第三节“位点分辨均裂与β-断裂热力学”中，均裂BDE为414.09–481.24 kJ mol^?1（ΔG_diss 279.12–504.80 kJ mol^?1），β-断裂ΔH_rxn 288.32–469.99 kJ mol^?1、ΔG_rxn 55.44–189.41 kJ mol^?1，均裂属高能耗引发、β-断裂属较低能耗传播，但ΔG皆正，仅为热力学非动力学能垒。第三节“计算裂解能量学描述性比较”以均值/极差汇总，均裂BDE均值446.23 kJ mol^?1、β-断裂ΔG_rxn均值远低于均裂ΔG_diss，强调C1–C10为结构溯源标号而非反应坐标，不作统计推断。第三节“选定低聚体内相对裂解敏感性”指出C2均裂BDE最低、C4/C6次之，β-断裂中C1焓最低、C10的ΔG最低，弱联图受苄位稳定与苯环电子效应共同控制。第三节“TGA降解窗口与多尺度分析作用”显示T_onset由334.0℃（5℃ min^?1）移至369.0℃（25℃ min^?1）、T_max由426.5℃移至461.5℃，R²=0.9685，属非等温动力学延迟而非本征稳定性变化，块体单峰质量损失包纳自由基网络而非单基元反应。第三节“PS降解表观活化能等转化演化”中Kissinger全局Ea=186.61 kJ mol^?1（R²=0.9977），FWO与KAS在α=0.05–0.95均值180.81与178.49 kJ mol^?1，FWO从143.10升至221.71 kJ mol^?1、KAS从140.10升至220.23 kJ mol^?1，低α对应弱联引发、中α对应β-断裂/解聚传播、高α对应芳残寡聚体抗降解。第三节“转化依赖动力学与低聚体热力学概念关联”定性映射三阶段，但拒绝将C1–C10指派至具体α区间。

讨论与结论部分，研究人员明确DFT量与TGA量物理尺度不同：前者是有限孤立气相低聚体选定反应的模型依赖热力学描述符（BDE、ΔH、ΔG），后者是凝聚态块体非等温质量损失导出的全局表观动力学参数，转化度α仅指可测失重比例而非特定键或片段；二者互补但无直接数值对应，亦未将单个裂解位置匹配特定温度/转化区间。结论重申十聚体均裂BDE 414.09–481.24 kJ mol^?1（C2最低）、β-断裂ΔG_rxn 55.44–189.41 kJ mol^?1，前沿能隙742.62 kJ mol^?1且轨道定域芳香/苄位；TGA侧Kissinger Ea=186.61 kJ mol^?1，残基校正后FWO/KAS均值180.81/178.49 kJ mol^?1，等转化Ea随α上升。研究意义在于示范“表观单峰失重+转化依赖Ea+分子非等价热力学”可共存，为PS化学回收中优先降低引发阶高能垒提供分子依据，但外推至高分子量、多分散、含添加剂、交联或环境老化微塑料须谨慎；局限含十聚体长度/构象/端基/气相假设及virgin样品未表征Mn/Mw/氧化指数等。全文未断言顺序拉链式机理，未给产率或催化性能，仅作多尺度机制框架。