Transitions between liquid crystalline phases investigated by dielectric and infra-red spectroscopies

这是一份关于利用介电光谱和红外光谱研究液晶化合物 11OS5 相变行为的详细技术总结。

1. 研究背景与问题 (Problem)

液晶材料在不同温度下会经历多种介晶相（如向列相、近晶相）以及晶体相的转变。理解这些相变过程中的分子动力学变化（如分子重排、取向有序度变化）对于材料科学至关重要。

• 研究对象：化合物 4-戊基苯基 -4'-十一烷氧基硫代苯甲酸酯（简称 11OS5）。该化合物在冷却过程中表现出丰富的相序列：各向同性液体 (Iso) → 向列相 (N) → 近晶 A 相 (SmA) → 近晶 C 相 (SmC) → 近晶 X 相 (SmX，一种具有六边形层内有序的倾斜近晶相) → 晶体 (Cr)。
• 核心挑战：
  1. 某些相变（如 SmA → SmC）在介电谱中信号微弱，难以检测。
  2. 需要区分不同相态下分子振动模式的细微变化，以揭示分子间相互作用和构象变化。
  3. 需要确定晶体相中分子的具体排列构型（头对头 vs. 头对尾）。
  4. 需要评估不同密度泛函理论（DFT）计算水平在预测红外光谱时的精度与计算成本的平衡。

2. 方法论 (Methodology)

本研究采用了实验测量与理论计算相结合的综合方法：

• 宽带介电光谱 (BDS)：

  • 在 0.1-10^7 Hz 频率范围内，对 11OS5 样品进行降温测量（373 K 至 273 K）。
  • 使用 Cole-Cole 模型拟合复介电常数，提取弛豫时间、介电强度、离子电导率及活化能。
  • 重点分析分子短轴旋转（flip-flop 运动，即 s-过程）和麦克斯韦 - 瓦格纳 - 席拉尔斯 (MWS) 界面极化过程。

• 红外光谱 (IR)：

  • 实验：在室温下使用反射模式，在降温过程中使用透射模式（KBr 压片），收集 370-4000 cm⁻¹ 范围内的光谱。
  • 数据分析：
    1. 相关系数矩阵：计算不同温度光谱之间的皮尔逊相关系数，以识别相变引起的谱图相似性变化。
    2. K-means 聚类分析：将光谱数据分为不同簇，以自动检测相变边界。
    3. 峰位追踪：监测特定吸收带（如 C=O 伸缩振动）随温度的位移。

• 理论计算 (DFT)：

  • 使用 Gaussian 16 软件，采用多种基组（def2SVP, def2SVPP, def2TZVP, def2TZVPP, 631+Gd, 6311+Gdp）和泛函（BLYP-D3(BJ), B3LYP-D3(BJ)）计算孤立分子的红外光谱。
  • 通过缩放系数校正计算频率，并与实验值对比，评估均方根误差 (RMSE) 和计算时间。
  • 构建二聚体模型（头对头 head-to-head 和头对尾 head-to-tail），模拟晶体相中的分子间相互作用，解释实验观测到的频移。

3. 主要结果 (Key Results)

A. 介电光谱分析

• 弛豫过程：在 Iso、N、SmA 和 SmC 相中观察到两个弛豫过程：低频的 MWS 过程和高频的分子 s-过程（分子短轴翻转）。
• SmC → SmX 相变：
  • 在 SmC 到 SmX 相变时，s-过程的弛豫时间突然增加（频率降低一个数量级，比率约为 46 倍）。
  • s-过程的活化能从 SmA/SmC 相的 ~94 kJ/mol 激增至 SmX 相的 144 kJ/mol。
  • 介电强度显著下降。这表明 SmX 相中出现了长程位置有序，极大地阻碍了分子的翻转运动，尽管 SmX 相仍保留了一定的取向无序（s-过程依然存在）。
• SmA → SmC 相变：在介电谱中仅表现为极微小的台阶，难以通过介电数据清晰分辨。

B. 红外光谱与 DFT 计算

• 理论计算评估：
  • def2TZVPP/B3LYP-D3(BJ) 组合给出了最小的 RMSE (11.6 cm⁻¹)，但计算时间最长 (156.9 小时)。
  • def2TZVP/B3LYP-D3(BJ) 的 RMSE 略高 (11.9 cm⁻¹)，但计算时间减半 (85.8 小时)，性价比最高。
  • 631+Gd/B3LYP-D3(BJ) 提供了极佳的平衡，RMSE 为 12.5 cm⁻¹，计算时间仅需 17.6 小时。
• 相变检测：
  • 相关系数矩阵：能区分 Iso/N、SmA/SmC（部分）、SmX 和晶体相，但未能清晰区分 SmA 和 SmC。
  • K-means 聚类：全谱分析能识别 N→SmA 和 SmX→Cr 转变，但无法区分 SmC 和 SmX。分波段聚类分析发现，不同振动模式对不同相变的敏感度不同（例如，C-H 伸缩振动对 Iso→N 敏感，而芳香环变形对 N→SmA 敏感）。
• 关键频移现象：
  • C=O 伸缩振动 (~1663-1674 cm⁻¹)：随温度降低发生红移（向低频移动），特别是在 SmX→Cr 结晶过程中红移最大（~8 cm⁻¹）。这归因于弱氢键的形成。
  • 芳香环面内变形/C-S 伸缩 (~901-909 cm⁻¹)：随温度降低发生蓝移（向高频移动）。
  • SmA → SmC：由于倾斜角仅 11°且层间距变化极小，红外光谱无法区分这两个相。
  • SmC → SmX：频移很小（<1 cm⁻¹），表明尽管 SmX 具有三维位置有序，但分子取向无序度仍较高。

C. 晶体相分子构型推断

• 通过对比孤立分子、头对头二聚体和头对尾二聚体的理论光谱：
  • 头对头构型：C=O 基团靠近相邻分子的 S 原子，导致 C=O 振动蓝移。
  • 头对尾构型：C=O 基团与相邻分子的烷氧链 CH2 基团相互作用（O...H 距离 2.6 Å），导致 C=O 振动红移。
• 结论：实验观测到 C=O 带随降温发生红移，且蓝移不明显，因此推断 11OS5 在晶体相中主要采取头对尾 (head-to-tail) 或交错排列的构型，而非头对头构型。

4. 主要贡献 (Key Contributions)

1. 多模态相变表征：成功结合 BDS 和 IR 技术，揭示了 11OS5 从介晶相到晶体相的完整动力学和结构演变，特别是捕捉到了 SmC → SmX 相变中分子动力学的剧烈变化（活化能激增）。
2. 计算方法优化：系统评估了多种 DFT 基组和泛函在预测复杂有机液晶分子红外光谱时的性能，为未来类似研究提供了关于精度与计算成本平衡的实用指南（推荐 631+Gd/B3LYP-D3(BJ) 或 def2TZVP/B3LYP-D3(BJ)）。
3. 分子构型解析：利用实验频移趋势与二聚体 DFT 计算的对比，有力证明了晶体相中分子采取头对尾排列，揭示了 C=O 与烷氧链间的弱相互作用机制。
4. 数据分析方法应用：展示了相关系数矩阵和 K-means 聚类在复杂相变序列检测中的潜力与局限性，指出在缺乏先验知识时，需结合特定波段分析才能准确识别细微相变（如 SmA/SmC）。

5. 意义与影响 (Significance)

• 理论意义：深化了对具有复杂相序列（特别是包含 SmX 等晶体类近晶相）的硫代酯类液晶分子动力学行为的理解。揭示了位置有序（SmX）对分子旋转动力学的显著抑制作用。
• 方法论价值：证明了将机器学习（聚类分析）与传统光谱学结合，可以有效辅助相变检测，但也强调了在液态晶体这种渐变有序体系中，先验知识的重要性。
• 应用前景：对分子构型和相互作用的理解有助于设计具有特定相变温度和动力学特性的新型液晶材料，特别是在需要快速响应或特定介电性能的应用场景中。

综上所述，该研究通过实验与理论的深度交叉验证，不仅精确描绘了 11OS5 的相变图谱，还阐明了微观分子排列与宏观物理性质（介电、光谱）之间的内在联系。