2025-11-21T22:13:16.282978

Marcus Theory and The Condon Approximation Revisited I: E-SHAKE and Seam Sampling

Cofer-Shabica, DeRosa, Subotnik

Marcus theory is the workhorse of theoretical chemistry for predicting the rates of charge and energy transfer. Marcus theory overwhelmingly agrees with experiment -- both in terms of electron transfer and triplet energy transfer -- for the famous set of naphthalene-bridge-biphenyl and naphthalene-bridge-benzophenone systems studied by Piotrowiak, Miller, and Closs. That being said, the agreement is not perfect, and in this manuscript, we revisit one key point of disagreement: the molecule C-13-ae ([3,equatorial]-naphthalene-cyclohexane-[1,axial]-benzophenone). To better understand the theory-experiment disagreement, we introduce and employ a novel scheme to sample the seam between two diabatic electronic states (E-SHAKE) through which we reveal the breakdown of the Condon approximation and the presence of a conical intersection for the C-13-ae molecule; we also predict an isotopic effect on the rate of triplet-triplet energy transfer.

academic

Marcus Theory and The Condon Approximation Revisited I: E-SHAKE and Seam Sampling

基本信息

论文ID: 2510.11810
标题: Marcus Theory and The Condon Approximation Revisited I: E-SHAKE and Seam Sampling
作者: D. Vale Cofer-Shabica, Jennifer R. DeRosa, Joseph E. Subotnik
分类: physics.chem-ph
发表时间: 2025年10月13日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2510.11810v1

摘要

Marcus理论是理论化学中预测电荷和能量转移速率的重要工具。对于Piotrowiak、Miller和Closs研究的萘-桥-联苯和萘-桥-二苯甲酮系统，Marcus理论在电子转移和三重态能量转移方面与实验结果高度一致。然而，理论与实验之间仍存在分歧，本文重新审视了一个关键的不一致点：分子C-13-ae (3,赤道-萘-环己烷-1,轴向-二苯甲酮)。为了更好地理解理论-实验分歧，作者引入并采用了一种新颖的采样方案来采样两个非绝热电子态之间的接缝(E-SHAKE)，通过该方法揭示了Condon近似的失效以及C-13-ae分子中圆锥交叉的存在，并预测了同位素效应对三重态-三重态能量转移速率的影响。

研究背景与动机

问题背景

Marcus理论的重要性：Marcus理论是预测电子转移(ET)和三重态能量转移(TET)速率的核心理论工具，其表达式为： $k_{D→A} = \frac{2π}{\hbar}|H_{DA}|^2(4πλk_BT)^{-1/2}e^{-(λ+ΔG^0)^2/(4λk_BT)}$
理论假设：Marcus理论基于几个关键假设：
- 热化的谐振子浴
- 激活交叉极限
- 黄金规则适用性
- Condon近似：电子耦合被视为常数
实验观察的异常：在Closs分子系列中，'ee'(赤道-赤道)系统表现出完美的指数衰减关系，但'ae'(轴向-赤道)和'ea'(赤道-轴向)系统显示出更复杂的行为。

研究动机

理论-实验差异：特别是C-13-ae分子的实验速率与理论预测存在显著差异
Condon近似的有效性：需要验证在"柔性"分子中Condon近似是否仍然适用
分子刚性vs柔性：探究分子柔性如何影响能量转移过程

核心贡献

开发了E-SHAKE方法：一种新的电子约束分子动力学采样技术，能够采样非绝热电子态之间的接缝
揭示了Condon近似的失效：通过接缝采样发现C-13-ae分子中电子耦合的剧烈变化
发现了圆锥交叉：在C-13-ae分子中识别出能量可及的圆锥交叉点
预测同位素效应：预测氘代替换将对三重态能量转移速率产生显著影响
提供新的理论洞察：解释了为什么某些Closs分子的行为偏离Marcus理论预测

方法详解

局域化非绝热态构建

非绝热-绝热转换

通过旋转绝热态构建非绝热态： $|Ξ_i⟩ = \sum_{j=1}^{N_{states}} (U)_{ji} |Φ_j⟩$

Boys和BoysOV方法

Boys标准：最大化偶极矩差异 $f_{Boys}(U) = \sum_{i,j=1}^{N_{states}} |⟨Ξ_i|\hat{μ}|Ξ_i⟩ - ⟨Ξ_j|\hat{μ}|Ξ_j⟩|^2$
BoysOV方法：分别局域化占据和虚轨道 $f_{BoysOV}(U) = \sum_{i,j=1}^{N_{states}} |⟨Ξ_i|\hat{μ}_{occ}|Ξ_i⟩ - ⟨Ξ_j|\hat{μ}_{occ}|Ξ_j⟩|^2 + |⟨Ξ_i|\hat{μ}_{virt}|Ξ_i⟩ - ⟨Ξ_j|\hat{μ}_{virt}|Ξ_j⟩|^2$

E-SHAKE算法

核心思想

E-SHAKE是SHAKE/RATTLE算法的电子版本，用于在分子动力学中施加电子约束。

约束条件

位置约束： $σ(t) = E_D(R(t)) - E_A(R(t)) = 0$
速度约束： $\dot{σ}(t) = \dot{R}(t) · ∇σ|_{R(t)} = 0$

算法实现

使用速度Verlet积分器的约束版本： $R(t+h) = R(t) + hQ$ $\dot{R}(t+h) = Q + \frac{1}{2}M^{-1} · [hF(t+h) + kG(t+h)]$

其中通过数值方法求解Lagrange乘数以满足约束条件。

技术创新点

电子约束的实现：首次将SHAKE算法扩展到电子自由度
接缝空间采样：能够系统性地探索非绝热态交叉的整个接缝空间
实时耦合监测：在动力学过程中实时监测电子耦合的变化
通用性：可应用于任意电子约束，不仅限于能量差约束

实验设置

计算细节

电子结构方法：CIS/6-31G*基组，使用Q-Chem 6.0
约束精度：能量差约束精度为4×10^-5 a.u. (≈1 meV)
温控：使用随机修正Berendsen方法，时间常数τ = 10 fs
积分步长：标准分子动力学时间步长

分子体系

研究对象为Closs分子系列：

C-13-ee, C-13-ea, C-13-ae
C-14-ee, C-14-ea, C-14-ae
D-26, D-27 (十氢化萘桥联系统)

分析方法

接缝能量分析：以最低T1能量为零点，分析高能区域的耦合变化
原子贡献分解：通过梯度分析确定主导耦合变化的原子
附着-脱离密度分析：分析供体-受体波函数重叠的变化

实验结果

主要发现

1. 电子耦合的显著差异

C-13-ae vs C-13-ea：C-13-ae的电子耦合比C-13-ea小两个数量级
接缝采样结果：C-13-ae的耦合在接缝空间中跨越两个数量级的变化范围

2. Condon近似的失效

C-13-ee和C-13-ea：耦合相对稳定，符合Condon近似
C-13-ae：耦合围绕零值分布，严重违反Condon近似

3. 圆锥交叉的识别

当 $H_{DA} → 0$ 且 $E_D → E_A$ 时，系统从非绝热接缝转变为圆锥交叉： $\lim_{H_{DA}→0, E_D→E_A=E} \begin{pmatrix} E_D & H_{DA} \\ H_{DA} & E_A \end{pmatrix} = \begin{pmatrix} E & 0 \\ 0 & E \end{pmatrix}$