热变质煤的化学结构与电子性质：多技术分析与建模的见解

《ACS Omega》：Chemical Structure and Electronic Properties of Thermally Altered Coals: Insights from Multi-Technique Analysis and Modeling

　　

随着浅层煤炭资源的逐渐枯竭，热变质煤的研究日益重要，特别是理解其化学结构与性质。在本论文中，基于元素分析（ultimate analysis）、¹³C 核磁共振（¹³C NMR）、X射线光电子能谱（XPS）和高分辨率透射电子显微镜（HRTEM）的结果，研究人员建立了两种热变质煤的分子模型。根据实验结果，构建了SL-A和SL-B的分子模型，其分子化学式分别为C_{207H118N3O5和C242H130N4O5。模拟的¹³C NMR谱图与实验数据吻合良好，验证了两个模型的可靠性。通过分子力学（MM）和分子动力学（MD）模拟，进一步获得了能量最小构象的三维分子结构。聚集模型显示SL-A和SL-B的晶胞尺寸分别为3.31 × 3.31 × 3.31 nm³和3.54 × 3.54 × 3.54 nm³。利用量子化学方法，计算了静电势（ESP）和前沿分子轨道（FMO）以评估模型的反应活性。结果表明，两种煤模型中的负电势主要与含氧官能团（如?C＝O）和芳香环相关，而氮和氢原子表现出强正电势。FMO分析揭示，热变质煤中高的变质程度与增加的反应活性相关。因此，研究热变质煤的化学结构和电子性质为其进一步加工利用（包括气化、燃烧和煤石墨化）提供了有价值的见解。}

**论文解读文章：热变质煤的化学结构与电子性质——基于多技术分析与建模的研究**

**研究背景与问题**

随着浅层煤炭资源的日益枯竭，热变质煤（thermal altered coal）的研究愈发重要。岩浆侵入（magmatic intrusion）煤层会创造高温高压条件，加速煤层的热演化，导致其物理性质与化学成分改变。前人已在岩石学、光学特性及化学组成方面对热变质煤开展了广泛研究，但对其化学结构与电子性质的深入理解仍显不足。特别是，不同变质程度的热变质煤中芳香环的缩合程度、杂原子赋存形态及其反应活性位点尚未系统揭示。这些结构信息对预测煤在后续加工利用（如气化、燃烧和石墨化）中的行为至关重要。尽管FTIR、Raman、XRD、XPS、HRTEM和¹³C NMR等实验手段已用于表征热变质煤，但基于实验数据构建合理的分子模型并结合量子化学计算评价其反应活性的研究仍较少。因此，开展热变质煤的分子建模与电子性质研究，对于从微观层面理解其化学反应机理具有重要理论和实际意义。

**研究内容与结论**

研究人员选取了中国安徽水俚煤田5号煤层中受岩浆岩影响的两份热变质煤样品（SL-A和SL-B，分别距岩床0.8 m和0.3 m），通过工业分析、元素分析、密度测定、¹³C NMR、XPS和HRTEM等实验手段表征其基本性质。基于实验结果，构建了SL-A（C<sub>207H118N3O5）和SL-B（C242H130N4O5）的二维分子模型，并通过分子力学（MM）和分子动力学（MD）模拟获得能量最低的三维构象。利用量子化学方法计算静电势（ESP）和前沿分子轨道（FMO），评价模型的反应活性。研究论文发表在《ACS Omega》。

**主要关键技术方法**

研究涉及以下关键方法：（1）¹³C核磁共振（¹³C NMR）确定碳骨架结构参数；（2）X射线光电子能谱（XPS）分析氧和氮官能团类型与分布；（3）高分辨率透射电子显微镜（HRTEM）定量统计芳香层尺寸与类型；（4）元素分析与原子数比例计算分子化学式；（5）分子力学（MM）与分子动力学（MD）模拟优化三维构象并计算密度；（6）基于密度泛函理论（DFT）的量子化学计算（GGA-PBE泛函）获取ESP和FMO。样品来源：安徽水俚煤田5号煤层，SL-A距岩床0.8 m，SL-B距岩床0.3 m。

**研究结果**

**3.1 煤常规特征**
通过镜质组最大反射率（Ro,max）、工业分析与元素分析发现：随距岩床距离减小，Ro,max、C和N含量增加，S、H和O含量降低。SL-B的Ro,max高达6.78%，属于超级无烟煤（mate-anthracite）。

**3.2 ¹³C NMR分析**
¹³C NMR谱图显示芳香碳峰强度远大于脂肪碳和羰基碳峰。通过分峰拟合获得12个碳结构参数。SL-A和SL-B的芳香碳（fa'）分别为88.61%和92.76%，脂肪碳（fal）分别为9.42%和6.07%，羰基/羧基碳（faC）分别为1.97%和1.17%。芳香桥头碳比例（XBP）分别为0.23和0.21，表明芳香环数需结合HRTEM进一步确定。

**3.3 XPS分析**
氧元素主要以醚键（C–O–C）形式存在，少量羰基（?C＝O）和羧基/酯基（?COO?）。氮元素在两种煤中主要存在于吡咯氮（pyrrolic nitrogen），其次是氧化氮；SL-A中吡啶氮（pyridinic nitrogen）和季氮（quaternary nitrogen）含量分别为18.09%和7.74%，SL-B中为12.04%和16.01%，表明变质程度影响氮形态转换。

**3.4 HRTEM分析**
芳香层尺寸统计表明：SL-A平均芳香环数为4.16，其中萘（naphthalene）占51.07%，2×2和3×3芳香条纹分别占25.63%和19.02%；SL-B平均环数为5.40，萘占46.69%，2×2和3×3分别占26.16%和19.78%。靠近岩床的煤样中更大芳香环的比例增加。

**3.5 分子结构的构建与校正**
基于上述实验参数，在ACD/Chemsketch中构建了初始平面分子模型，并通过gNMR软件模拟¹³C NMR谱图与实验谱图对比，反复调整原子比例直至吻合。最终分子式：SL-A为C207H118N3O5，SL-B为C242H130N4O5，元素比例和芳香度与实验值高度一致。通过MM和MD优化获得三维构象，SL-B因强π–π相互作用呈现平行取向，而SL-A因刚性键未出现明显卷曲。总势能显著降低，主要归因于键伸缩能（EB）和范德华能（Evan）的减小。密度模拟得到SL-A和SL-B的适宜密度分别为1.38和1.44 g/cm³，略低于实验值（1.79和1.71 g/cm³），可能是实验密度包含无机矿物和孔隙填充小分子所致。

**3.6 量子化学分析**
**3.6.1 静电势（ESP）分析**：两种煤模型中负电势分布于芳香环和含氧官能团（如?C＝O）附近，而氢和氮原子区域显示正电势。SL-A最大负电势为?8.723×10^?2 au，最大正电势为6.384×10^?2 au；SL-B相应值分别为?7.155×10^?2 au和6.189×10^?2 au。
**3.6.2 前沿分子轨道（FMO）分析**：HOMO和LUMO主要定域于芳香环，呈π特征；SL-A中的羰基也显示较高活性。SL-A和SL-B的HOMO–LUMO能隙（ΔEgap）分别为0.216 eV和0.110 eV，暗示靠近侵入岩的煤样（SL-B）由于变质程度更高，分子结构反应活性更强，更易参与亲电和亲核反应。

**讨论与结论**

研究讨论部分指出，所建分子模型成功再现了实验¹³C NMR谱图，密度结果在可接受范围内。ESP分析揭示了煤分子中亲电和亲核反应位点；FMO分析表明高变质程度增强煤的反应活性。这些结果为热变质煤在气化、燃烧和石墨化等加工利用中的行为提供了微观理论支持。研究结论翻译如下：

本研究使用¹³C NMR、XPS、HRTEM、氦密度测试等实验技术对两种热变质煤样品进行了表征，成功构建了其分子结构模型，并采用量子化学计算评估了化学反应活性。主要结论总结如下：

（1）SL-A和SL-B模型的分子式分别为C207H118N3O5和C242H130N4O5，分子量分别为2727.17和3173.64。经分子动力学优化后，模拟的¹³C NMR谱图与实验结果高度吻合，验证了分子模型的准确性。两种模型中，C原子主要以芳香环形式存在，O原子主要以醚键交联芳香环，N原子以吡咯、吡啶或季氮基团形式存在。

（2）分子动力学优化后，热变质煤分子模型呈现增强的三维构象。桥键和脂肪键发生显著扭曲。特别是SL-B模型中，芳香环间的π–π相互作用导致这些结构几乎平行取向。SL-A和SL-B经分子动力学处理后的势能分别为1043.410和1315.874 kcal/mol，主要归因于键伸缩能和范德华能的降低。

（3）ESP分析表明，两种热变质煤模型中的负电势位于芳香环和含氧官能团（如?C＝O）附近，而正电势位于N和H原子附近。HOMO–LUMO能隙（ΔEgap）表明，随着变质程度增加，热变质煤对亲电和亲核反应的反应活性增强。</sub>