丙醇异构体中的介电弛豫与平衡热力学描述符：射频测量与openCOSMO-RS分析

《Molecules》：Dielectric Relaxation and Equilibrium Thermodynamic Descriptors in Propanol Isomers: Radio-Frequency Measurements and openCOSMO-RS Analysis

　　

研究人员通过射频测量（频率范围0.3–32 MHz，共七个固定频率）在+20至?160 °C温度区间内考察了正丙醇与异丙醇的介电弛豫行为，并结合openCOSMO-RS计算进行分析。研究人员根据介电损耗峰的位置确定了表观特征时间，并采用Arrhenius与Eyring形式化方法进行处理。每种样品均获得九组温度扫描数据，这些数据来自三个相互独立的测量电容器，每个电容器重复扫描三次。基于电容器的Arrhenius拟合给出的表观活化能分别为：正丙醇21.30 ± 0.16 kJ·mol^?1（标准差SD；平均值95%置信区间20.91–21.69），异丙醇21.99 ± 0.16 kJ·mol^?1（95%置信区间21.60–22.38）。对三个独立电容器估计值进行的保守Welch比较给出ΔE_a = 0.687 kJ·mol^?1（95%置信区间0.332–1.042，p = 0.0058）。尽管这一亚千焦（sub-kJ）量级的差异在电容器层级分析中具有统计学可分辨性，但由于该比较未纳入系统温度不确定性和未解析的光谱重叠，研究人员不认为其具有机理性意义。Eyring活化焓分别为19.60 ± 0.15和20.23 ± 0.15 kJ·mol^?1（SD）。熵和吉布斯能值仅在κ = 1条件下作为表观量报告。openCOSMO-RS给出的液体至正己烷无限稀释转移焓为30.5和28.0 kJ·mol^?1。上述转移量包含氢键以及非特异性溶剂化、色散和堆积贡献，因此研究人员未将其等同于氢键焓。文献量热数据给出的纯醇氢键焓量级约为17 kJ·mol^?1。Kirkwood因子大于1表明存在强烈的静态取向关联。由于七频率窗口无法独立分辨Debye过程与结构α过程，动力学参数被解释为主导介电损耗过程的有效描述符，而非特定模式的势垒。

论文解读：丙醇异构体介电弛豫与平衡热力学描述符的联合分析

1. 研究背景与问题

单羟基醇的介电弛豫同时反映局部分子运动和与瞬态氢键网络相关的集体极化关联。宽带介电研究表明，慢Debye型贡献与较快的结构α弛豫可以共存，氢键交换、链重排和长程偶极互相关未必共享单一微观时间尺度。因此，从介电弛豫时间导出的活化能本质上是一种有效动力学描述符，而非某一特定O–H···O键的解离焓。正丙醇与异丙醇具有相同分子式和羟基官能团，但分子几何与堆积方式不同，为研究支化对氢键组织和弛豫行为的影响提供了理想对比体系。然而，异构体之间的细微差异必须经过定量不确定度分析，才能判断其是否具有机理性意义。已有介电和结构研究表明，支化可以改变氢键组织和弛豫行为，但小差异的机理解释需要严格的统计支撑。

现有实验数据在较小频率窗口内仅能获得主导损耗过程的有效参数，无法独立分解Debye和α贡献。同时，动力学量（如Arrhenius活化能）与平衡热力学量（如转移焓、氢键焓）在物理意义上存在本质区别，但数值上的相似性往往被误读为机制等价。因此，本研究旨在不预设等价性的前提下，比较丙醇异构体的表观活化能尺度、Kirkwood静态取向关联因子以及openCOSMO-RS平衡分子表面描述符，并特别强调有限频率窗口带来的局限、异构体间活化能小差异的不确定度，以及动力学与平衡量之间不同的物理意义。

2. 研究内容与主要技术方法

研究人员对高纯度正丙醇和异丙醇（≥99.9%，分析纯，Sigma-Aldrich）经精馏干燥后，采用Keysight E4990A精密阻抗分析仪，在0.3、0.8、2.0、5.0、9.0、20和32 MHz七个固定频率下，于+20至?160 °C温度范围内进行射频介电测量。每种样品使用三个独立测量电容器，每个电容器重复三次扫描，共获得九组温度扫描数据。在每一固定频率下，由介电损耗ε″的温度极大值确定损耗峰温度T₀，并利用关系式2πf₀τ_app ≈ 1赋予表观特征时间。基于电容器的Arrhenius和Eyring拟合用于获取活化参数。量子化学计算采用RDKit距离几何方法生成分子结构，GFN2-xTB预优化后，在DFT/BP86/def2-TZVP(-f)/COSMO水平优化，再于BP86/def2-TZVPD/CPCM水平计算屏蔽表面电荷分布；热力学性质采用开源openCOSMO-RS 3.11（24a参数化）计算σ-剖面、σ-电势、无限稀释活度系数及转移焓。

3. 研究结果

3.1 介电弛豫的温度–频率依赖关系

Arrhenius图中表观特征时间τ_app的?lg值对100/T作图近似线性。三个电容器独立拟合给出的活化能分别为：正丙醇21.146、21.301、21.458 kJ·mol^?1，异丙醇21.832、21.988、22.146 kJ·mol^?1。平均值为正丙醇21.30 ± 0.16 kJ·mol^?1（95% CI 20.91–21.69，平均R2 = 0.9870），异丙醇21.99 ± 0.16 kJ·mol^?1（95% CI 21.60–22.38，平均R2 = 0.9938）。Welch检验给出ΔE_a = 0.687 kJ·mol^?1（95% CI 0.332–1.042，p = 0.0058），但该亚千焦差异未被视为机制性证据。Cole–Cole图近似半圆形，但高频侧存在系统偏差，表明存在重叠的更快色散过程，单一时间常数描述仅为近似。拟合的静态介电常数ε₀和高频极限参数ε_∞显示，ε_∞超过光学极限n2，证明测量窗口高频边缘仍有较快介电色散。Kirkwood–Fr?hlich关系计算的g_K在正丙醇中约为3.4–3.7，异丙醇中约为3.3–3.8，远大于1，表明强烈的静态取向关联。

3.2 Eyring参数化

表观弛豫速率k_rel = 1/τ_app的Eyring分析给出正丙醇ΔH_act = 19.60 ± 0.15 kJ·mol^?1（95% CI 19.22–19.97），异丙醇20.23 ± 0.15 kJ·mol^?1（95% CI 19.85–20.60）。在κ = 1条件下，表观活化熵分别为?3.55 ± 0.35和?4.36 ± 0.34 J·mol^?1·K^?1。200 K时ΔG_act,app分别为20.31 ± 0.08和21.10 ± 0.08 kJ·mol^?1。研究人员指出，透射系数κ的不确定性会影响熵和吉布斯能的绝对解释，因此这些量仅作为表观值。

3.3 σ-剖面与分子极性

两种异构体的σ-剖面均呈现单羟基醇的特征三区域分布：中央烷基区域（σ ≈ 0）、负σ供体区（羟基氢）和正σ受体区（氧孤对电子环境）。供体和受体区域在两种异构体间非常相似，表明羟基本征供受体特征相近；而中央区域差异显著，正丙醇的烷基区最大值较紧凑，异丙醇因支化呈现更结构化的分布。该结果支持异构体主要差异源于分子几何和堆积的结论。

3.4 无限稀释转移焓与实验热力学背景

openCOSMO-RS计算的液体至正己烷无限稀释转移焓ΔH_tr^∞为正丙醇30.5 kJ·mol^?1、异丙醇28.0 kJ·mol^?1。实验量热文献中，Solomonov等报道的脂肪醇纯液体氢键焓量级为16.9–17.7 kJ·mol^?1，液体至环己烷转移焓约23–24 kJ·mol^?1。这些数据显示“缔合能”的数值强烈依赖于热力学参考过程。介电活化能（21.30–21.99 kJ·mol^?1）、实验转移值和氢键焓均处于数十kJ·mol^?1的宽范围内，但该量级相似不具机制诊断意义。非氢键液体中同样存在活化转动、结构扩散和粘滞过程，因此不能将介电势垒定量拆分为氢键与范德华/堆积贡献。

3.5 σ-电势与静态介电常数的温度演化

在220、250和298 K下，两种醇的σ-电势均显示供体侧特征。冷却时供体分支变得更负，正丙醇从约?7变为?15 kJ·mol^?1，表明在COSMO-RS框架内供体–受体相互作用越来越有利。实验上ε₀也随冷却增加，两者趋势定性相容，但均不决定介电损耗过程的微观时间依赖。

3.6 超额焓与热力学一致性检验

openCOSMO-RS计算的醇–正己烷混合物在298 K的摩尔超额焓H^E均为正值，正丙醇最大约0.95 kJ·mol^?1，异丙醇约1.01 kJ·mol^?1，与稀释时有利醇–醇相互作用被吸热破坏一致。该曲线仅作为定性一致性检验，未声称定量验证指标。

4. 讨论总结与结论翻译

研究的主要局限包括：七个固定频率无法独立拟合Debye和α贡献，因此报道的τ_app、E_a和Eyring参数为有效量而非模式分辨参数；峰温重复性已从九次扫描量化，但未包含绝对温度校准和谱峰重叠的完整系统不确定度预算；g_K的形式化传播不确定度无法计算；残留水和冷却速率未定量记录，可能影响最深过冷区。openCOSMO-RS作为平衡热力学模型，其转移焓包含特异性和非特异性相互作用，无法与介电实验共同定量分离氢键、色散、堆积和集体动力学贡献。

研究结论：本研究报道了正丙醇与异丙醇的射频介电测量及openCOSMO-RS平衡描述符。基于三个独立电容器（每电容器三次重复扫描）的电容器级Arrhenius分析给出E_a = 21.30 ± 0.16和21.99 ± 0.16 kJ·mol^?1（SD），95%置信区间分别为20.91–21.69和21.60–22.38 kJ·mol^?1。电容器级差异为0.687 kJ·mol^?1（Welch 95% CI 0.332–1.042，p = 0.0058），但该亚千焦分离未被赋予机制意义，因为系统温度不确定度和未解析光谱重叠未包含在统计模型中。Eyring分析给出表观活化焓19.60 ± 0.15和20.23 ± 0.15 kJ·mol^?1（SD）。Kirkwood因子远大于1支持强静态取向关联，但其形式化传播不确定度无法从现有输入重建。计算的液体至正己烷转移焓（30.5和28.0 kJ·mol^?1）不同于实验纯醇氢键焓（约17 kJ·mol^?1），且包含大量非特异性相互作用贡献，因此这些能量尺度的数值相似性不能作为共同微观势垒的证据。由于七频率测量未分离Debye和α过程，提取的动力学参数应视为主导射频损耗过程的有效描述符。实验与计算结果的联合最有价值之处在于为缔合和介电动力学提供互补约束，更具体的机制归属需要宽带模式分辨光谱。