[논문 리뷰] Excitation energy transfer rates: comparison of approximate methods to the exact solution

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

Spectroscopy and Quantum Chemical Studies참고 문헌 1인용 수 6

한 줄 요약

이 논문은 광합성 복합체에서의 자극 에너지 전달(EET) 속도를 정확한 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) 속도에 대해 널리 사용되는 근사 방법의 정확도를 평가하기 위해.
• 포르스터, 레드필드, 변분 마스터 방정식과 같은 근사 방법이 유효한 매개변수 영역을 규명하기 위해.
• 변분 마스터 방정식이 중간 재조정 에너지 영역에서의 성능을 평가하기 위해.
• 다양한 체계 매개변수에 걸쳐 어느 근사 방법이 일관되게 EET 속도를 과소 또는 과대평가하는지 확인하기 위해.
• 다이머 및 트리머 체계에 대해 HEOM 방법을 통해 정확한 EET 속도를 확립하여 향후 이론적 및 수치적 연구의 기준을 제공하기 위해.

제안 방법

• 시간 적분된 속도 커널 $\overline{K}$를 풀어 정확한 EET 속도를 계산하기 위해 계층적 운동 방정식(HEOM)을 사용하며, 이는 모든 양자 효과를 포괄한다.
• 계산 가능성을 확보하기 위해 적응형 절단 기법을 적용하여 최대 7개의 색소를 포함하는 체계에 대해 HEOM를 적용하고, 큰 희박 행렬의 영공간을 효율적으로 계산한다.
• 통일된 투사 연산자 형식을 사용하여 포르스터, 레드필드, 수정된 레드필드, 변분 마스터 방정식 방법의 근사 속도 커널 $\overline{K}$를 유도한다.
• 모든 모델 간 일관된 비교를 가능하게 하기 위해 진동 환경의 드루드-로렌츠 스펙트럼 밀도를 사용한다.
• 양자 마코프 핵심을 임의의 기저에서 고전적 속도 커널 $\overline{K}$로 매핑하여 실험 결과 및 근사 모델과 직접 비교할 수 있도록 한다.
• 다양한 재조정 에너지 $\lambda$와 전자 결합 상수 $H_{\rm el}$에 대해 정확한 HEOM 결과와 근사 방법 간의 수치적 비교를 수행한다.

실험 결과

연구 질문

• RQ1 포르스터, 레드필드, 수정된 레드필드 근사가 정확한 HEOM 해와 비교해 정확한 EET 속도를 제공하는 매개변수 영역는 어디인가?
• RQ2 변분 마스터 방정식은 임의의 $\lambda$를 대상으로 설계되었음에도 불구하고 중간 재조정 에너지 영역에서 정확한 EET 속도를 제공하지 못하는 이유는 무엇인가?
• RQ3 트리머 체계에서 특정 기저(예: 엑시톤 기저 또는 사이트 기저)에서 EET 속도가 음수가 되는 현상을 최적의 중간 기저를 선택함으로써 피할 수 있는가? 이는 EET 동역학의 물리적 해석에 어떤 함의를 갖는가?
• RQ4 근사 방법이 일관되게 EET 속도를 과소 또는 과대평가하는 정도는 어느 정도이며, 이 행동은 체계 크기나 결합 강도에 따라 달라지는가?
• RQ5 변분 폴라론 변환의 한계는 무엇인가? 특히 중간 결합 영역에서의 퍼티르베이티브 처리를 위한 상호작용 해밀토니안을 단순화하는 데서 발생하는가?

주요 결과

• 포르스터, 레드필드, 수정된 레드필드 방법은 각각의 페르투베이티브 한계 내에서만 정확하며, 이는 실제 광합성 복합체의 물리적으로 관련 있는 매개변수 범위를 벗어나는 경우가 많다.
• 트리머 체계에서 근사 방법은 $H_{\rm el}$과 $\lambda$의 값에 따라 EET 속도를 과소 또는 과대평가할 수 있으며, 이는 매개변수 공간 전반에 걸쳐 일관된 편향이 없음을 시사한다.
• 변분 마스터 방정식은 자유 에너지 기능 $A_0$의 국소 최소값 간의 전환으로 인해 최적의 이동 함수 $F_m(\omega)$에 급격한 변화가 생기면서 예측된 속도에 불연속적인 점프를 보인다.
• 트리머 체계에서 중간 $\lambda$ 영역에서는 사이트 기저와 엑시톤 기저 양쪽에서 EET 속도가 음수가 되며, 이는 표준 기저에서 동역학을 고전적 점프 과정으로 해석할 수 없음을 나타낸다.
• 사이트와 엑시톤 표현 사이를 보간하는 기저 $\ket{M}$ 이 존재할 경우 전체 매개변수 영역에서 양수의 속도를 얻을 수 있으며, 이는 EET 동역학을 더 견고하게 기술하는 잠재적 길을 시사한다.
• 변분 마스터 방정식이 중간 $\lambda$ 영역에서 실패하는 것은 현재 변분 안사드의 한계를 시사하며, 이는 폴라론 변환의 일반성 부족이나 $A_0$에 대한 최적화 기준의 비최적성 때문일 수 있다.