利用静电嵌入的透热化方法研究有机分子晶体中的光物理

《Journal of Chemical Theory and Computation》：Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

　　

高发射性有机分子晶体在多个领域具有应用，例如有机电子学、太阳能电池和传感器。理解这些应用背后的激发态机制对于优化和有效控制它们至关重要。耦合非绝热动力学（特别是量子动力学）的激子模型为分子晶体中的光化学和光物理过程提供了关键见解。然而，目前仍缺乏通用的工具和自动化工作流来促进此类模拟。在本文中，研究人员提出了一种计算策略来研究有机分子晶体的光激活动力学，桥接了传统上用于分子系统和材料科学的方法，特别关注局部激发和电荷转移（CT）过程之间的相互作用。研究人员在fromage和Overdia程序之间实现了一个接口，能够在激发态ONIOM(QM:QM′)框架内构建分子晶体的振动耦合哈密顿量，并通过基于RESP的Ewald求和纳入长程静电效应。片段基透热化（FrD）提供了弱到中等耦合区间的量子动力学模拟途径。该方法应用于二苯并[g,p]屈（DBC）晶体的光物理，使用时间依赖密度泛函理论（TD-DFT）。fromage/Overdia接口用于计算晶体中18个独特DBC二聚体的局部激发和CT态的耦合，并量化静电嵌入的影响，发现其影响适中（10–20%）。对π-堆积二聚体的模拟再现了实验观察到的从溶液到晶体的微小红移，这归因于固定单体之间的电子相互作用而非晶体静电效应。量子动力学模拟揭示了从明亮的局部激发到CT态的超快布居转移。该方法建立了一个连接分子和固态激发态动力学的稳健框架，具有研究分子晶体中激发、缺陷和杂质的潜在应用。

**论文解读：利用静电嵌入的透热化方法研究有机分子晶体中的光物理**

**研究背景与问题**

有机分子晶体因其高发射特性在有机电子学、太阳能电池和传感器等领域具有广泛应用。理解这些应用背后的激发态机制，尤其是局部激发（LE）与电荷转移（CT）过程之间的相互作用，对于优化和调控器件性能至关重要。激子模型结合非绝热动力学（特别是量子动力学，QD）为分子晶体中的光化学和光物理提供了关键见解，但现有方法和工具缺乏通用性和自动化工作流，难以模拟晶体环境中长程静电耦合和振动耦合效应。传统方法多聚焦于分子系统或单独的材料科学，未能有效整合两者。为此，研究人员开发了一种结合静电嵌入（EE）的透热化（Diabatization）方法，旨在构建一个统一的框架，以研究弱到中等耦合区间内分子晶体中的激发态动力学。

**主要技术方法**

研究人员采用从age和Overdia程序的接口（fro_overdia.py），基于ONIOM(QM:QM′)方法构建分子晶体的振动哈密顿量。技术核心包括：1）通过RESP（Restrained Electrostatic Potential）电荷的Ewald求和（Ewald summation）实现长程静电嵌入（EE），模拟晶体环境；2）利用片段基透热化（FrD）将超分子复合物的绝热态投影到准透热基组上，定义局部激发（LE）和电荷转移（CT）态；3）结合线性振动耦合（LVC）模型，通过数值微分计算透热态的线性耦合参数，构建FrD-LVC(EE)哈密顿量；4）采用多层多组态时间相关Hartree（ML-MCTDH）方法进行量子动力学（QD）模拟，包含90个振动模式。样本来源：二苯并[g,p]屈（DBC）晶体的实验晶体结构（CCDC: 1481424），经周期性DFT（PBE-D3）弛豫后用于簇模型构建。

**研究结果**

**4.1 18个DBC晶体二聚体的片段透热化**
通过从age自动筛选晶体中18个独特二聚体，计算了激子耦合（LE-LE）和LE-CT耦合。利用FrD方法与其他方法（VEET、VCoul、VTrESP）对比，发现π-堆积二聚体（Dimer 1）具有最大的Frenkel耦合（约27 meV），其LE-CT耦合也最强（63 meV）。静电嵌入（EE）对耦合的调节幅度为10–20%，对LE-CT耦合影响更显著。分析表明，激子耦合强度不仅取决于单体间距离，还受堆积几何和前线轨道形状影响。

**4.2 π-堆积结构**
聚焦Dimer 1（π-堆积二聚体），通过ONIOM(QM:QM′)优化S0和S1态，发现晶体环境中单体、二聚体和三聚体的吸收和发射均出现红移（二聚体S1垂直跃迁红移0.09 eV，发射红移0.06 eV）。在S1最小能量处，激发主要局域在单个单体上，但少量电子密度转移至相邻单体，导致红移和振子强度下降。实验观察到从溶液到晶体的红移（约0.21 eV）被部分再现，并归因于单体间电子相互作用而非晶体静电效应。

**4.3 π-二聚体的FrD与FrD(EE)分析**
在晶体、S0-ONIOM和S1-ONIOM三种几何结构下，计算了FrD和FrD(EE)模型。静电嵌入使较大耦合（≥5 meV）变化10–15%，特定CT耦合变化超过50%。在S1-ONIOM几何中，涉及单体a的透热态能量降低，耦合增强，而涉及单体b的状态不变。S0-ONIOM的FrD(EE)结果与π-三聚体结果一致，验证了二聚体模型的有效性。

**4.4 ML-MCTDH量子动力学模拟**
基于FrD-LVC(EE)哈密顿量，模拟了光激发后的布居演化。激发至|LE_a¹?或|LE_b¹?后，单体间布居转移缓慢（250 fs后约30–40%），CT态贡献极小。激发至明亮|LE_a^S?或|LE_b^S?后，超快（约100 fs）布居转移至|CT(b→a)?（约40%），同时伴随向LE₁态的衰减。模拟还发现，初始激发为二聚体透热态（S2和S6）时，产生超快退相干，布居快速分配至两个LE₁态。静电嵌入对动力学影响较小，主要使CT态布居略微降低。

**结论与讨论**
研究人员开发了从age与Overdia接口，实现了基于ONIOM(QM:QM′)和Ewald静电嵌入的FrD-LVC(EE)模型，用于分子晶体中激子动力学研究。该方法在弱到中等耦合区间内有效，适用于激发局域化但受晶体环境耦合的场景。以DBC晶体为例，通过二聚体筛选、几何优化和量子动力学模拟，揭示了π-堆积二聚体中激发态行为：LE-LE耦合弱，导致局域化；激发至明亮LE态时，超快CT态参与，形成主要弛豫路径。静电嵌入对耦合和动力学影响适中（10–20%），但可能对极性晶体更显著。该工作为连接分子和固态激发态动力学提供了可扩展工具，有望用于研究晶体中的缺陷、杂质和聚集诱导发光等。未来工作将扩展至更宽耦合区间和更大聚集模型。