《Journal of Chemical Theory and Computation》：Coarse-Grained Nonadiabatic Dynamics in Disordered Condensed Phases
编辑推荐：
在无序凝聚相中，临界能隙涨落与非绝热动力学背后的机制仍不清楚，由于这些体系固有的复杂性，其精确模拟构成挑战。研究人员在本研究中首次尝试利用高效粗粒化（CG）模型考察这些基本机制。采用严谨的自底向上参数化方法，研究人员为两个原型体系构建了CG模型：四氢呋喃中的胡
在无序凝聚相中，临界能隙涨落与非绝热动力学背后的机制仍不清楚，由于这些体系固有的复杂性，其精确模拟构成挑战。研究人员在本研究中首次尝试利用高效粗粒化（CG）模型考察这些基本机制。采用严谨的自底向上参数化方法，研究人员为两个原型体系构建了CG模型：四氢呋喃中的胡萝卜素-卟啉-富勒烯三元体，以及乙腈中的苝二酰亚胺二聚体。利用非绝热半经典映射动力学，研究人员评估了从全原子显式分子模拟、中间CG分辨率到高度约化的多态谐波（MSH）模型的一系列表征，以考察其对能隙涨落和非绝热性质的影响。这些显式分子模拟的一个显著优势是能产生随时间变化的径向分布函数，揭示光诱导电荷转移下局域溶剂重排的详细分子图像。结果表明，通过对溶剂与溶质进行适当的粗粒化，所得CG模型忠实保留了全原子基准的结构响应、成对重组能以及电子布居动力学。相反，激进的2-重原子溶质映射严重扭曲了局域非键环境。重要的是，从结构有效的CG模拟参数化的MSH模型在丢弃数万自由度的情况下成功复现了显式分子动力学。本工作为将非绝热动力学模拟扩展至更大更复杂凝聚相体系提供了高效的新理论基石。
研究背景与动机
在太阳能转换与有机光伏（OPV）材料中，光诱导电荷转移（CT）是无序凝聚相（如体异质结共混物）中的核心过程。凝聚相环境并非静态背景，静态无序（源于局域堆积形态差异）与动态无序（源于热运动与电声耦合）共同调制非绝热机理。全原子（AA）非绝热模拟虽准确但受限于计算成本，仅能处理数纳米空间尺度；而有效模型如多态谐波（MSH）模型虽能压缩自由度，却丢失显式分子运动信息。传统粗粒化（CG）力场多基于自由能匹配，过度平滑短程非键相互作用，破坏驱动电子转移的能隙时间关联函数（TCF）。因此，亟需一种自底向上CG框架，在保留局域非键能景与能隙TCF的前提下降低自由度。该研究发表于《Journal of Chemical Theory and Computation》，旨在桥接空间尺度，为大规模光化学模拟提供新工具。
关键技术方法
研究人员采用自底向上CG策略：键合项由直接Boltzmann反演（DBI）获取，静电项由组成原子电荷标量和赋值；范德华项通过Lennard-Jones静态势匹配（LJSPM）结合有效体积校正（EVC）参数化，以复现AA非键能景。非绝热传播统一采用基于Meyer-Miller-Stock-Thoss（MMST）哈密顿的半经典映射动力学，包括Ehrenfest平均场（MF）、线性化半经典（LSC）与对称准经典（SQC）等变体。显式分子动力学在AA与CG层级采样能隙TCF，据此构建MSH模型（200个简正模离散谱密度）。体系为四氢呋喃中胡萝卜素-卟啉-富勒烯（CPC60）三元体（三态）与乙腈中苝二酰亚胺（PDI）二聚体（四态），电子耦合取Condon近似常数，态特异性能量修正由量子化学激发能扣除AA力场内能给出。
研究结果
- 2.
理论
2.1 全原子多态哈密顿：在绝热表示下总哈密顿含核动能与对角/非对角势能矩阵，力场采用GAFF（键合+Lennard-Jones+静电），态特性能量修正Wα锚定垂直能隙至高精度电子结构。
2.2 自底向上CG多态哈密顿：CG势能同构AA，DBI定键合、电荷求和定静电、LJSPM-EVC定范德华；EVC以零势等值面约束珠半径避免密度畸变；WαCG同理修正。
2.3 多态谐波模型：将显式环境映射为共享简谐振子浴，由能隙TCF经F(F-1)/2个关联提取重组能Er(αβ)与谱密度J(αβ)(ω)，离散为200模，平衡位移Sj(d,α)严格满足所有成对重组能约束。
2.4 非绝热半经典映射动力学：MMST将电子态映射为经典变量，零级能（ZPE）参数γ区分MF(0)/SQC(1/3)/LSC(1/2)；约化密度矩阵由Wigner变换轨迹平均获得。
- 3.
模拟细节
CPC60/THF盒含1三元体+2700 THF（~70 ?），PDI/ACN盒含1二聚体+1440 ACN（~50 ?）；NPT 2 ns + NVT 1 ns平衡，20 ns基态采样取初值；NAMD在NVE下传播（三元体15 ps，二聚体5 ps）。CG层级涵盖AA、联合原子（UA）、2-重原子（2h）溶质与3珠（3b）/2珠（2b）溶剂。
- 4.
结果
结构响应：CG（AA/UA、AA/3b等）径向分布函数（RDF）与AA一致捕捉CT后局域溶剂重排；2h溶质因抹平非键环境使RDF失真。能隙涨落：适当CG保留能隙均值、方差及TCF衰减；激进CG破坏高频分量。重组能与谱密度：成对Er(αβ)在合理CG下与AA吻合，J(αβ)(ω)低频区保真。布居动力学：LSC/SQC/MF在CG（非2h）与AA基准一致；MSH由CG采样参数化后复现显式分子NAMD，自由度自数万降至~600。2h映射因局部摩擦异常导致布居转移速率偏差。
讨论与结论翻译
讨论指出，CG用于非绝热动力学必须超越结构复现，须保留驱动电子跃迁的能景与动力学特征；LJSPM-EVC自底向上传递非键能景是关键。MSH自有效CG模拟参数化可继承显式分子TCF，实现尺度跨越。结论：适当溶剂-溶质粗粒化CG模型忠实保留结构响应、成对重组能与电子布居动力学；激进2-重原子溶质映射严重扭曲局域非键环境；自结构有效CG模拟构建的MSH模型在丢弃数万自由度下成功复现显式分子动力学。本工作为将非绝热动力学模拟扩展至更大更复杂凝聚相体系提供了高效新理论基石。
