Coherent nuclear and electronic dynamics in primary charge separation in photosynthetic reaction centers: a Redfield theory approach

Primary charge separation dynamics is modeled in the pheophytin-modified Rhodobacter sphaeroides R-26 reaction center (RC). To explain the observed spectral evolution, it is assumed that the process is coupled to coherent nuclear motion. A density matrix equation with the Redfield relaxation superoperator is used for simulation of the electron-vibrational dynamics and its spectral signatures. The model includes two diabatic states, i.e., an excited state P * of the primary donor (i.e., special pair, P), and a charge-transfer state (P+B-, which is the primary photoproduct in the pheophytin-modified RC). The strong coupling of these states with two collective nuclear modes is supposed. The mixing of diabatic states (with different displacements along each of the two nuclear coordinates) results in a complicated potential surface that determines the dynamics of the excited-state wave packet. The coupled nuclear and charge-transfer dynamics is calculated in the basis of vibronic eigenstates obtained by numerical diagonalization of the electron-vibrational Hamiltonian. The third-order nonlinear response associated with excited-state dynamics is calculated, including the P * f P stimulated emission (SE) and the P+B- f P+(B-) * excited-state absorption (ESA). The model allowed us to obtain a quantitative fit of the experimental kinetics of the SE near 900-950 nm and the ESA in the 1020-nm region of the pheophytin-modified Rhodobacter sphaeroides R-26 RC (Yakovlev, A. G.; Shkuropatov, A. Ya.; Shuvalov, V. A. FEBS Lett. 2000, 466, 209). By use of the parameters adjusted from the fit, we hav