何种从头算（ab initio）计算方法描述离子液体蒸气中的电荷转移？以EMIM-OTF为案例研究揭示前沿轨道与电离阈值

《RSC Advances》：Which ab initio calculation methods describe the charge transfer in ionic liquid vapors? Case study on EMIM-OTF reveals the frontier orbitals and the ionization threshold

　　

研究人员对1-乙基-3-甲基咪唑三氟甲磺酸盐（[EMIM][OTF]）蒸气进行了时间飞行质谱（TOF-MS）测量，以深入了解真空紫外（VUV）吸收阈值。实验数据通过杂化密度泛函理论（DFT）、时间依赖DFT（TD-DFT）以及多种多体微扰理论（MBPT）方法进行解读：包括二阶MBPT、Hedin的GW和Bethe–Salpeter方程（BSE）计算。据研究人员所知，这是对离子液体（IL）蒸气的DFT泛函进行的首次全面比较。研究人员评估了超过60种交换关联近似，包括许多范围分离杂化（RSH）泛函及其最优调谐变体。理论方法预测了广泛的能量范围；只有实验TOF-MS数据能够识别出真正具有预测性的方法。实验光谱揭示了一个光学带隙为6.75 eV的电荷转移（CT）态，以及一个在8.3 eV处的最大离子产额信号。结合实验与理论，研究人员得以重新审视离子对几何结构并重构[EMIM][OTF]的前沿轨道。结果表明，最高占据分子轨道（HOMO）和最低未占据分子轨道（LUMO）的能量分别为-9.05 eV和-2.30 eV；同时还确定了HOMO?1和LUMO+1的能量。尽管采用非杂化M11-L泛函进行的ΔSCF计算给出了合理的电离阈值，但正确描述光学带隙和紫外光电子能谱（UPS）需要≥50%的精确Hartree–Fock（HF）交换。RCAM-B3LYP泛函以及Bartlett的QTP泛函家族在预测TD-DFT计算中的CT激发方面表现出良好性能。没有方法能够预测所提出的LUMO能量，尽管M08-SO泛函在整体上最为接近。在本工作的局限性内，全局杂化PBE-2X泛函展现出最一致的整体性能。其范围分离杂化（RSH）类似物CAM-PBE50*在计算CT激发方面也表现良好。

**论文解读：离子液体蒸气中电荷转移的从头算方法评估——以EMIM-OTF为例**

**研究背景与问题**

离子液体（ILs）是熔融温度低于100°C的有机盐，因其极低蒸气压、宽液程、高热电化学稳定性及强溶解能力，在电沉积、杂环合成及多种电化学过程中得到广泛应用。准确掌握其前沿轨道能量和价带组成对于基础理解及器件设计至关重要。然而，许多ILs（包括[EMIM][OTF]）的气相电子结构仍知之甚少。例如，文献中报道的[EMIM][OTF]气相HOMO-LUMO（最高占据分子轨道-最低未占据分子轨道）带隙跨度几乎达6 eV，这种不一致主要源于两方面：一是缺乏可靠的气相实验数据，导致同时发表差异巨大的带隙值；二是实验数据的解释并非总是直接，例如液相[EMIM][TFSI]的紫外光电子能谱（UPS）中存在低至~2 eV的未解释低能尾。此外，采用杂化DFT计算虽能与大多数实验UPS谱定性吻合，但仍需系统性偏移，而预先知道特定杂化DFT泛函所需偏移量并不容易。更重要的是，ILs中激发态和电荷转移（CT）能量几乎未被探索。为此，研究人员对[EMIM][OTF]蒸气进行了VUV（真空紫外）光激发下的能量依赖TOF-MS（时间飞行质谱）测量，以获取其激发态和未占据电子态的信息。论文发表在《RSC Advances》。

**主要技术方法**

研究人员采用了多种从头算方法：杂化密度泛函理论（DFT）和时间依赖DFT（TD-DFT）用于基态和激发态计算；多体微扰理论（MBPT）包括G₀W₀（单次GW近似）和Bethe–Salpeter方程（BSE）用于准粒子能量和中性激发。计算使用MolGW和NWChem软件包，基组主要为aug-cc-pVDZ，部分采用aug-cc-pVTZ。特别评估了超过60种交换关联泛函，包括全局杂化、范围分离杂化（RSH）及其最优调谐（OT）变体。

**研究结果**

1. **GW结果**：G₀W₀和自洽GW（scGW）对大多数杂化泛函给出相似结果，但基组依赖性显著。G₀W₀@aug-cc-pVDZ对HOMO能量给出优异结果，但难以推广至其他ILs。

2. **非杂化DFT**：HLE16、HLE17等功能优化用于带隙，表现优于HCTH类，但所有非杂化泛函（如TPSS、SCAN、M11-L）均严重低估吸收阈值，且ΔSCF电离能与HOMO能量偏差大（约3 eV）。HCTH-407+产生虚假低能CT激发，表明非杂化泛函无法准确描述ILs。

3. **Meta-GGA泛函**：非杂化meta-GGA（如SCAN）仍严重低估HOMO能量，但杂化版本SCAN0和R2SCAN0的ΔSCF电离能与实验HOMO能量一致。M08-SO（含56% HF交换）预测正确的HOMO能量，且对LUMO能量表现最佳。

4. **全局杂化泛函**：标准全局杂化（如B3LYP、PBE0）产生虚假CT跃迁和错误带隙。HF交换分数是关键：≥50% HF交换（如PBE-2X、BHHLYP）可抑制虚假CT态，但过高HF会损害热化学精度。PBE-2X是唯一能定量再现所有关键实验能量的泛函。

5. **RSH泛函**：Bartlett的QTP泛函避免虚假CT激发，其中QTP17在TD-DFT激发能和UPS描述上最佳，但高估A_1,2双峰强度。α=0的RSH（如LC-BLYP、LC-BLYP-EA）需BSE计算描述中性激发，LC-BLYP-EA在基态DOS和UPS上表现优异。

6. **最优调谐（OT）**：OT显著改善RSH性能。LC-BLYP的默认μ=0.33a₀^-1恰好是[EMIM][OTF]的最优值。CAM-LDA0虽非最优调谐，但表现良好。LRC-PBE0经OT后（μ=0.23a₀^-1）能量改善超1.5 eV。CAM-PBE50*（μ=0.09）预测与实验良好吻合。

7. **与文献基准比较**：文献中CAM-B3LYP、M06-2X、ωB97X-D等常用泛函在IL离子对上表现并非最佳，因HF交换要求更高。QTP17和PBE-2X等高HF泛函更优。

8. **LUMO**：几乎所有方法高估所提出的LUMO能量（-2.3 eV）。meta-GGA泛函（如M08、SCAN）对LUMO描述较好，而高HF泛函（如PBE-2X、ωB97X）给出正LUMO能量，不适用于LUMO计算。

**结论与讨论**

研究人员强调，结论受限于单一构象、主要使用aug-cc-pVDZ基组及实验PIY（部分离子产额）数据分辨率较粗。研究重新评估了蒸发温度180°C下的IL构象，发现预测正确构象极具挑战。[EMIM][OTF]的CT过程可理解为从阴离子到阳离子的电子转移，观测到的CT能量与估算的离子对静电稳定能吻合。通过结合从头算与实验数据，推导出一组准粒子前沿轨道，并实现了对多种计算方法的严格评估。杂化泛函对准确描述离子对电子结构至关重要，所有结果强烈依赖于HF交换含量。高HF泛函有效消除虚假低能CT态，但多数TD-DFT仍高估A_1,2双峰强度。最优调谐RSH泛函优于标准版本。BSE计算有助于揭示激发态，但系统性地低估CT能量并引入虚假态。ΔSCF计算（如M11-L）或GW方法可成功再现实验垂直电离能。PBE-2X泛函是唯一在62种方法中定量满足所有关键实验能量（垂直电离能9.05 eV，第一CT激发6.75 eV，主CT跃迁8.3 eV）的泛函，且无虚假CT态。然而，没有任何方法能同时提供正确的HOMO和LUMO能量，建议在需要LUMO能量时使用meta-GGA泛函。