[论文解读] Excitation energy transfer rates: comparison of approximate methods to the exact solution

Simon Jesenko, Marko Žnidarič|arXiv (Cornell University)|May 16, 2014

Spectroscopy and Quantum Chemical Studies参考文献 1被引用 6

一句话总结

本论文通过将近似方法——Förster、Redfield、改进的Redfield以及变分主方程——与通过层级运动方程（HEOM）方法计算出的精确速率进行比较，评估了光合作用复合物中的激发能量转移（EET）速率。研究发现，在电子-声子耦合的中间 regime 中，没有任何近似方法能可靠预测EET速率，此时精度对系统参数极为敏感；即使变分主方程在此关键区域也失效。

ABSTRACT

We have evaluated excitation energy transfer rates in photosyntetic complexes using the exact HEOM method and various approximate methods frequently used in the literature, namely, the Förster method, the Redfield method, the modified Redfield method and the variational master equation. The rates are evaluated for the case of a simple dimer and a trimer photosynthetic complex, with vibrational environment characterized by the Drude-Lorentz spectral density. Comparing approximate rates to the exact ones, we have confirmed the validity of approximate methods within appropriate limits, however these limits are often well outside parameter ranges that are relevant for the dynamics in real photosynthetic complexes.

研究动机与目标

• 评估广泛使用的近似方法在光合作用复合物中激发能量转移（EET）速率方面的准确性。
• 识别Förster、Redfield和变分主方程等近似方法仍有效的参数区域。
• 评估变分主方程在中间重组能区域的表现。
• 确定是否存在某种近似方法在不同系统参数下始终高估或低估EET速率。
• 通过HEOM方法为二聚体和三聚体系统建立精确的EET速率基准，以供未来理论和计算研究参考。

提出的方法

• 使用层级运动方程（HEOM）求解时间积分速率核 $\overline{K}$，以计算精确的EET速率，该核捕捉了所有量子效应。
• 采用自适应截断方案，使HEOM在最多七个色素的系统中计算可行，从而高效计算大规模稀疏矩阵的零空间。
• 利用统一的投影算符形式，推导出Förster、Redfield、改进的Redfield及变分主方程方法的近似速率核 $\overline{K}$。
• 使用Drude-Lorentz谱密度描述振动环境，确保各模型间的一致性比较。
• 将量子马尔可夫核映射到任意基底下的经典速率核 $\overline{K}$，以实现与实验结果及近似模型的直接比较。
• 在不同重组能 $\lambda$ 和电子耦合 $H_{\rm el}$ 下，对精确HEOM结果与近似方法进行数值比较。

实验结果

研究问题

• RQ1 与精确HEOM解相比，Förster、Redfield和改进的Redfield近似在何种参数区域内能给出准确的EET速率？
• RQ2 尽管变分主方程被设计用于任意 $\lambda$，为何其在中间重组能区域仍无法提供准确的EET速率？
• RQ3 在三聚体系统中，能否通过选择最优中间基底避免某些基底（如激发子基或位点基）出现负的EET速率？这对EET动力学的物理解释有何含义？
• RQ4 近似方法在多大程度上系统性地高估或低估EET速率？这种行为是否依赖于系统尺寸或耦合强度？
• RQ5 变分极化子变换在中间耦合区域中，对简化相互作用哈密顿量以进行微扰处理的局限性是什么？

主要发现

• Förster、Redfield和改进的Redfield方法仅在其各自的微扰极限内准确，而这些极限通常不在真实光合作用复合物的物理相关参数范围内。
• 对于三聚体系统，近似方法可能因 $H_{\rm el}$ 和 $\lambda$ 的取值不同而高估或低估EET速率，表明在参数空间中不存在一致的偏差。
• 变分主方程预测的速率出现不连续跳跃，这是由于自由能泛函 $A_0$ 的局部极小值之间发生切换，导致最优位移函数 $F_m(\omega)$ 突变所致。
• 在三聚体系统中，当 $\lambda$ 处于中间值时，EET速率在位点基和激发子基中均可能出现负值，表明在标准基底下无法将动力学解释为经典的跃迁过程。
• 存在一个介于位点表示与激发子表示之间的基底 $\ket{M}$，可在整个参数区域内产生正值速率，提示这可能是实现更稳健的EET动力学描述的潜在路径。
• 变分主方程在中间 $\lambda$ 区域的失效表明当前变分试探波函数存在局限，可能源于极化子变换的通用性不足或 $A_0$ 最小化准则不够优化。