热加工技术有潜力将木质素转化为化学品、材料和燃料等高值化产品,热解是其中最基础的过程,深入认识木质素热解反应机理对发展木质素高效利用工艺非常重要。鉴于木质素结构和热解反应的复杂性,已有的实验方法和量子化学计算方法还难以获得木质素热解过程的全貌和动态演化细节,本论文利用大规模ReaxFFMD模拟新方法,采用基于GPU的高性能计算程序GMD-Reax直接模拟了 11800-77200原子规模的木质素模型热解过程,并利用化学反应与可视化分析工具VARxMD深入挖掘了热解过程全貌及其背后反应机理,系统地研究了软木、硬木和硫酸盐木质素热解过程中的连接键、中间体、官能团及主要产物的演化行为,揭示了不同种类木质素热解行为的差异与共性规律。本论文取得的主要结果如下:(1)基于木质素不同的植物来源、不同的分离方法、不同的单体聚合观点构建了软木、硬木、硫酸盐三类木质素的7种分子模型:Adler和Freudenberg软木木质素支化高分子模型、Ralph-SL软木木质素线型高分子模型、Nimz硬木木质素支化高分子模型、Ralph-HL硬木木质素线型高分子模型、以及Marton和Crestini硫酸盐木质素支化高分子模型。这些3D模型规模均大于10,000原子,结构覆盖了 7种含氧官能团和20种连接键。这些大规模木质素模型可以更好地描述木质素高分子的支化或线型结构的多样性,有助于考察木质素结构在热解反应过程中的相互影响。(2)提出了采用木质素单体间连接键和其所连接的两个苯环的组合结构作为识别特定连接键反应子结构的策略,获得了木质素中常见1 1种连接键整体结构的初始主导转化路径:Cα/β-O醚键断裂反应在α-O-4、β-O-4、α-O-4&β-5、α-O-4&β-O-4&5-5和β-1&α-O-α连接键整体结构的初始转化路径中占主导地位(>85%);Cα-O醚键断裂反应和苯环开环反应在β-β&γ-O-α连接键整体结构的初始转化路径中是两类势均力敌的主导反应路径(各约50%);苯环开环反应是4-O-5、5-5和β-N连接键整体结构的初始主导转化路径(>95%)。此外,相邻连接键整体结构之间的反应活性会相互影响:α-O-4和β-O-4等反应活性较高的连接键发生醚键断裂后生成的苯氧基会促进其相邻的、反应活性较低的4-O-5等连接键整体结构发生键断裂反应或苯环开环反应。(3)模拟获得的苯环可以发生开环反应的结论与前人的量子力学计算和冲击波实验得到的结果一致,本论文进一步揭示了木质素热解过程中导致稳定的芳环发生开环反应的重要中间体及主要路径:①苯环在苯氧基的作用下被插入一个O原子后可转化为含氧7元杂环;②苯环在二环[3.1.0]己基结构的作用下发生3元环上的碳碳键断裂反应后可转化为5元脂环结构;③苯环在二环[3.1.0]己基结构的作用下发生5元环上的碳碳键断裂反应后可转化为3元脂环结构;④苯环在苄基结构的作用下被插入一个C原子后可转化为7元脂环结构。(4)获得了软木、硬木和硫酸盐木质素模型热解行为的差异与共性规律。差异体现为:①热重曲线的演化快慢与C-O醚键的含量密切相关,C-O醚键含量高的硬木木质素热反应性最高,其次为软木木质素,硫酸盐木质素的热反应性最低。②提出了依据C原子数区分的产物类型及其演化特征划分木质素热解阶段的策略,发现软木和硬木木质素热解过程具有显著的三阶段特征,而硫酸盐木质素热解过程仅呈现出后两个阶段。③不同种类木质素模型中β-O-4连接键结构上所含有的羰基、甲氧基及邻接α-O-4连接键等活性取代基的不同结构环境会影响其转化速率:活性取代基含量最高的Nimz硬木木质素转化速率最快;Adler软木和Freudenberg软木木质素次之;含量最低的Marton硫酸盐木质素转化速率最慢。类似的活性取代基效应也影响β-β&γ-O-α连接键整体结构的转化反应:含羰基(Cα=O)的Marton硫酸盐木质素及含醚键(Cα-O-4和-OCH3)的Nimz硬木木质素中其转化更早发生。当木质素中含有较多非活性Ph-(C)n-Ph键和较少活性PhO-C键时,会使碳碳型连接键整体结构的转化变慢。不同种类木质素热解的共性体现为:木质素单体的芳环结构在不同种类木质素热解反应过程中均可以在苯氧基、苯上含氧官能团或苄基的作用下转化生成含氧七元杂环、五元脂环或七元脂环等环结构中间体且各环结构的演化趋势相同。(5)考察了同一类木质素中仍然存在争议的支化和线型分子结构对木质素热解行为可能产生的影响,结果表明:不同木质素模型热解模拟获得的6种主要小分子产物的开始生成温度和最大生成量与木质素中相关官能团或连接键结构的含量密切相关,与其支化或线型高分子结构的关联性不大。(6)采用Solid-Py/SR-VUV-PI-TOF-MS实验装置进行了木质素热解实验研究,实验结果与模拟结果对比表明:实验检测到的16种C0-C3产物在模拟中也能获得;不同木质素模型热解模拟获得的苯、苯酚及其衍生物的种类少于实验结果。更大模型的模拟可以更好地描述苯、苯酚及其衍生物等热解产物的演化行为,大规模硫酸盐木质素模型的热解模拟获得的乙烯、甲醛、甲醇、硫化氢、甲硫醇、苯、苯酚、邻甲基苯酚、邻苯二酚和愈创木酚等产物随温度的演化趋势和实验定性一致,且热解产物在不同温度下随时间的演化规律均呈现出生成速率随温度升高而增大、温度越高其产量越快达到最大值的趋势。尽管木质素的Solid-Py/SR-VUV-PI-TOF-MS实验和ReaxFFMD模拟的时间尺度和温度区间存在较大差距,但实验和模拟获得的热解产物随温度的演化趋势以及在不同温度下随时间的演化趋势定性一致。本论文的工作表明:大规模ReaxFF MD模拟方法不仅可以定性预测木质素热解产物随温度和时间的演化概貌,并且可通过所获得的反应路径对预测结果加以解释,从而可为木质素热解的高值化应用提供深入的分子反应机理的认识,为发展木质素高效利用的工艺提供理论支持。