📄 化学

High-level G4 and G4(MP2) theoretical study into the gas-phase thermodynamic and electronic properties of atmospherically relevant terrestrial vegetation derived C10H16 monoterpenes and C15H24 sesquiterpenes

本研究采用高水平的 G4 和 G4(MP2) 计算方法，旨在确定 74 种具有大气相关性的 C10 和 C15 倍半萜类的气相热力学性质和绝热电离能，揭示了紧凑且高度桥联的结构最为稳定，且计算值与实验数据吻合良好，并受特定的电子和结构因素所支配。

原作者： Sierra Rayne, Kaya Forest

发布于 2026-10-23

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Sierra Rayne, Kaya Forest

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

我们呼吸的空气并非空无一物；它是一个由生物界塑造的、繁忙的化学景观。树木和植物不断释放出大量的不可见气体，即挥发性有机化合物，这些化合物漂浮到大气中，并与阳光及其他化学物质发生相互作用。其中最显著的是萜烯（terpenes），这是一类赋予森林特征性气味的大型碳氢化合物。虽然其中一些分子结构简单，但另一些则具有复杂的结构，其碳原子排列成环或链状。一旦释放，这些化合物并不会就此消失；它们会经历一系列化学转化，进而导致大气中微小颗粒的形成，而这些颗粒反过来又会影响云的形成和地球的气候。为了理解这些分子的行为方式，科学家需要了解它们的基本能量特性：其化学键内储存了多少能量、它们彼此之间的稳定性如何，以及它们失去电子变成带电状态的难易程度。这些细节对于预测它们在大气中的反应速度至关重要，但对于许多这类复杂的植物衍生分子而言，此类数据一直以来都难以获取或难以直接测量。

最近的一项研究通过使用强大的计算机模拟，绘制了陆地植被中七十四种不同萜烯分子的能量图谱。研究人员重点关注了两大类分子：单萜（monoterpenes，含有十个碳原子）和倍半萜（sesquiterpenes，含有十五个碳原子）。这些分子具有多种不同的形状，从开放的链状结构到紧密堆积的多环结构不等。由于这些化合物在自然界中通常不稳定或仅以极微量存在，在实验室中测量其性质往往是不可能的。相反，该团队采用了高水平的理论方法，本质上是在超级计算机上进行精确的数学计算，以确定形成每种分子所需的精确能量、一种形状转变为另一种形状时释放或吸收的能量，以及剥离分子中一个电子所需的特定能量。这种方法使他们能够生成一套完整且一致的数据，涵盖各种结构，为大气科学家提供了一个可靠的参考。

计算结果揭示了分子的形状如何决定其稳定性的清晰规律。研究发现，具有开放链状结构的分子通常是最不稳定的，与它们的环状对应物相比，承受着显著的能量惩罚。当碳原子形成环时，分子会变得更加稳定，但这种稳定程度很大程度上取决于环是如何融合的以及双键位于何处。研究发现，最稳定的分子是那些具有紧凑、高度桥联框架且能有效分散结构应力的分子。例如，一个将双键隐藏在环系统内部的分子，始终比一个双键从环外突出的分子更稳定。这种对内部、更拥挤的双键的偏好遵循了一个众所周知的化学原理，但该研究在广泛的复杂结构中证实了这一点。研究人员还发现，最稳定的倍半萜是那些能够最小化内部应力的刚性、融合环系统，而最不稳定的则是相同化学式下的柔性开放链版本。

除了稳定性之外，研究还绘制了这些分子释放电子的难易程度，这一属性被称为电离能。这个数值至%%于预测分子与大气清除剂（如羟基自由基）发生反应的速度。结果表明，失去电子的难易程度受分子能否有效分散产生的正电荷所支配。具有扩展共轭双键系统的分子（即电子可以在结构中自由移动的分子），其电离最为容易。同样，紧邻双键的应变三元环也提供了一种稳定带电状态的独特方式，从而降低了电离所需的能量。相比之下，那些被周围原子保护不足的孤立双键分子，则需要更多的能量来进行电离。研究发现，最稳定、最刚性的分子通常具有最高的电离能，而最活跃、最柔性的分子则具有最低的电离能。

研究人员将计算值与现有的少量此类化合物的实验数据进行了对比，发现两者高度吻之一，证实了其模拟的准确性。他们还测试了结果与其它理论方法的对比，发现一些较旧的、不够先进的方法会产生错误的稳定性排序，有时甚至会暗示一个稳定性较低的分子实际上是最稳定的。这凸显了使用高水平计算工具来确保物理机制正确的必要性。该研究还澄清了关于特定化合物（如蒎烯和柠檬烯）电离能的一些历史性混淆，通过提供一致且高精度的基准，调和了文献中相互矛盾的报告。通过建立这些精确的能量值，这项工作为改进大气化学模型奠定了坚实的基础。这些模型依赖于准确的数据来预测生物源排放物如何转化为二次有机气溶胶，而这些气溶胶在气候调节中起着重要作用。有了这个全新的数据集，科学家现在可以更好地理解这些植物衍生气体的去向及其对全球环境的影响。

技术摘要：萜类化合物的高级 G4 与 G4(MP2) 研究

问题陈述
生物源挥发性有机化合物（BVOCs），特别是单萜（C10H16）和倍半萜（C15H24），是驱动对流层化学的关键因素，影响着低层大气的氧化能力以及二次有机气溶胶（SOA）的形成。尽管其具有重要的环境意义，但由于缺乏高质量的基础热力学和电子数据，对其降解路径的全面理解受到了阻碍。这些化合物表现出巨大的结构多样性，包括无环、单环、双环和三环框架，其中许多结构不稳定或难以通过直接实验测量进行分离。因此，可靠的气相标准状态生成焓 ( $\Delta_fH^\circ(g)$ )、异构化焓 ( $\Delta_{isom}H^\circ(g)$ ) 和绝热电离能（AIEs）往往缺失或不一致，限制了大气动力学模型和结构-活性关系（SAR）预测的准确性。

研究方法
为了填补这些数据空白，本研究采用了高水平复合从头算量子化学方案，具体为 Gaussian-4 (G4) 和 Gaussian-4 改进的二阶 Møller-Plesset 微扰理论 (G4(MP2)) 方法。研究重点针对 74 种具有大气相关性的烃类：包括 20 种单萜异构体和 54 种源自陆地植被的倍半萜异构体。

计算工作流包括：

结构生成： 从 PubChem 数据库检索初始 3D 结构，并通过系统性的转子搜索使用 MMFF94 力场进行优化，以确定全局最小能量构象。
量子化学计算： 使用 Gaussian 16 软件包在加拿大数字研究联盟（Digital Research Alliance of Canada）的计算网格上进行高水平计算。
热力学性质： 使用原子化法计算气相标准状态生成焓 ( $\Delta_fH^\circ(g)$ )，利用氢和碳的实验原子生成焓。异构化焓 ( $\Delta_{isom}H^\circ(g)$ ) 相对于每一类中生成焓最低的异构体进行推导。
电子性质： 计算了气相标准状态（298.15 K）和 0 K 的绝热电离能（AIEs）。
验证： 将计算值与 NIST 化学数据库（NIST Chemistry WebBook）及既往文献中的可用实验数据进行对比，并与较低水平的方法（如 DFT、基团贡献法）进行对比，以评估准确性和一致性。

主要结果
本研究为所有 74 种化合物生成了一个自洽的热力学和电子性质数据集。主要发现包括：

方法论一致性： G4 和 G4(MP2) 方法产生了高度自洽的结果。G4(MP2) 计算的 $\Delta_fH^\circ(g)$ 和 $\Delta_{isom}H^\circ(g)$ 值系统性地比 G4 低几个 kJ mol⁻¹（范围为 -0.2 至 -6.6 kJ mol⁻¹）。同样，G4(MP2) 的 AIEs 也系统性地比 G4 低 0.01–0.05 eV（平均偏差 ≈0.03 eV 或 2.9 kJ mol⁻¹）。
热力学准确性： 计算的 $\Delta_fH^\circ(g)$ 值与有限的可用实验数据（例如针对金鸡菊烯 camphene 和 $\alpha$ -蒎烯 $\alpha$ -carene）表现出良好的一致性，表明其热化学准确度在 $\pm 4.2$ kJ mol⁻¹ 以内。
稳定性趋势：
- 单萜： 无环三烯被确定为热力学最不稳定的异构体。稳定性随环化程度增加而提高，紧凑且高度桥联的多环框架（如三环烯 tricyclene、金鸡菊烯 camphene、 $\alpha$ -芬烯 $\alpha$ -fenchene）代表了最稳定的构型。稳定性受环张力、双键取代（遵循扎伊采夫规则）和超共轭效应控制。内环双键的稳定性一致高于外环双键。
- 倍半萜： 观察到类似的趋势，无环异构体（ $\alpha$ -和 $\beta$ -法尼烯 $\alpha$ - and $\beta$ -farnesene）是最不稳定的。全局最小值被确定为 $\beta$ -补丁烯 $\beta$ -patchoulene，这是一种刚性的稠环三环结构。大环和单环框架（如锗曼烯 germacrenes）表现出中间稳定性。
电离能 (AIEs)： AIEs 范围约为 7.2 eV 至 8.6 eV。其数值主要受 $\pi$ $π$ -共轭、烷基取代和特定轨道相互作用的影响。
- 扩展的 $\pi$ -共轭（如无环三烯和共轭二烯中的共轭）显著降低了 AIEs。
- 烯烃与张力环之间的相互作用（特别是烯烃–环丙烷 Walsh 轨道相互作用）提供了额外的稳定作用，降低了如杜精烯 thujenes 和莎比宁 sabinene 等化合物的 AIEs。
- 高取代的内环双键通常比孤立或外环双键具有更低的 AIEs。
与其他方法的比较： 研究强调了高水平 G4/G4(MP2) 结果与较低水平密度泛函理论 (DFT) 或基团贡献法 (GCMs) 之间的显著差异。在若干案例中，较低水平的方法产生了错误的稳定性排序或错误的量级估计，凸显了对于这些复杂体系使用高水平复合方法的必要性。

意义与主张
作者声称，这项工作提供了一个全面的、内部一致的参考数据集，解决了萜类化合物热力学中长期存在的歧义。通过填补结构复杂异构体的数据空白，本研究为以下方面提供了精确的结构和能量输入：

精炼大气化学机制和氧化路径。
改进用于预测氧化剂攻击速率（如 $\cdot$ OH, O $_3$ , $\cdot$ NO $_3$ ）的结构-活性关系。
参数化二次有机气溶胶产量。
支持对对流层氧化能力和气溶胶驱动气候调节的稳健建模。

论文强调，G4 与 G4(MP2) 在刚性及柔性框架中的相互一致性，验证了使用这些方案来扩展此类萜类异构体数据集的可行性。此外，观察到的与常见较低水平方法之间的差异，进一步强调了这些高水平复合基准在应用到大气科学中的价值。

类似论文