[論文レビュー] 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

ひとこと要約

本稿は、光合成複合体における励起エネルギー移動（EET）率の評価において、Förster法、Redfield法、修正Redfield法、および変分マスター方程式といった近似手法を、階層的運動方程式（HEOM）法による正確な率と比較する。その結果、電子-格子結合が中程度の強さにある領域では、どの近似手法も信頼できる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率を確立することで、今後の理論的・計算的研究のベンチマークを提供すること。

提案手法

• 時間積分された率カーネル $\overline{K}$ を解くことにより、量子効果をすべて捉える階層的運動方程式（HEOM）を用いて正確なEET率を計算する。
• 適応的切断スキームを導入することで、最大7色素の系に対してもHEOMの計算を可能にし、大規模なスパース行列の核空間を効率的に計算できるようにする。
• 統一された射影作用素形式を用いて、Förster法、Redfield法、修正Redfield法、および変分マスター方程式の近似率カーネル $\overline{K}$ を導出する。
• すべてのモデル間で一貫した比較が可能となるよう、振動環境にドゥードル・ローレンツスペクトル密度を適用する。
• 任意の基底における量子的マコフ型カーネルを古典的率カーネル $\overline{K}$ に写像することで、実験的および近似モデルとの直接比較を可能にする。
• 再構成エネルギー $\lambda$ および電子的結合項 $H_{\rm el}$ を変化させたパラメータ領域において、正確なHEOM結果と近似手法との数値的比較を実施する。

実験結果

リサーチクエスチョン

• RQ1 Förster法、Redfield法、および修正Redfield法が、正確なHEOM解と比較して、どのパラメータ領域で正確な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$ の最適化基準の非最適性が関与している可能性がある。