We provide a critical assessment of typical phonon spectral densities, J(ω), used to describe linear and nonlinear optical spectra in photosynthetic complexes. Evaluation is based on a more careful comparison to experiment than has been provided in the past. J(ω) describes the frequency-dependent coupling of the system to the bath and is an important component in calculations of excitation energy transfer times. On the basis of the shape of experimental J(ω) obtained for several photosynthetic complexes, we argue that the shape of J(ω) strongly depends on the pigment-protein complex. We show that many densities (especially the Drude-Lorentz/constant damping Brownian oscillator) display qualitatively wrong behavior when compared to experiment. Because of divergence of J(ω) at zero frequency, the Brownian oscillator cannot fit a single-site spectrum correctly. It is proposed that a log-normal distribution can be used to fit experimental data and exhibits desired attributes for a physically meaningful phonon J(ω), in contrast to several commonly used spectral densities which exhibit low-frequency behavior in qualitative disagreement with experiment. We anticipate that the log-normal J(ω) function proposed in this work will be further tested in theoretical modeling of both time- and frequency-domain data.

The Fenna-Matthews-Olson (FMO) trimer (composed of identical subunits) from the green sulfur bacterium Chlorobaculum tepidum is an important protein model system to study exciton dynamics and excitation energy transfer (EET) in photosynthetic complexes. In addition, FMO is a popular model for excitonic calculations, with many theoretical parameter sets reported describing different linear and nonlinear optical spectra. Due to fast exciton relaxation within each subunit, intermonomer EET results predominantly from the lowest energy exciton states (contributed to by BChl a 3 and 4). Using experimentally determined shapes for the spectral densities, simulated optical spectra are obtained for the entire FMO trimer. Simultaneous fits of low-temperature absorption, fluorescence, and hole-burned spectra place constraints on the determined pigment site energies, providing a new Hamiltonian that should be further tested to improve modeling of 2D electronic spectroscopy data and our understanding of coherent and dissipation effects in this important protein complex.

The vibrational spectral density is an important physical parameter needed to describe both linear and non-linear spectra of multi-chromophore systems such as photosynthetic complexes. Low-temperature techniques such as hole burning (HB) and fluorescence line narrowing are commonly used to extract the spectral density for a given electronic transition from experimental data. We report here that the lineshape function formula reported by Hayes et al. [J. Phys. Chem. 98, 7337 (1994)] in the mean-phonon approximation and frequently applied to analyzing HB data contains inconsistencies in notation, leading to essentially incorrect expressions in cases of moderate and strong electron-phonon (el-ph) coupling strengths. A corrected lineshape function L(ω) is given that retains the computational and intuitive advantages of the expression of Hayes et al. [J. Phys. Chem. 98, 7337 (1994)]. Although the corrected lineshape function could be used in modeling studies of various optical spectra, we suggest that it is better to calculate the lineshape function numerically, without introducing the mean-phonon approximation. New theoretical fits of the P870 and P960 absorption bands and frequency-dependent resonant HB spectra of Rb. sphaeroides and Rps. viridis reaction centers are provided as examples to demonstrate the importance of correct lineshape expressions. Comparison with the previously determined el-ph coupling parameters [Johnson et al., J. Phys. Chem. 94, 5849 (1990); Lyle et al., ibid. 97, 6924 (1993); Reddy et al., ibid. 97, 6934 (1993)] is also provided. The new fits lead to modified el-ph coupling strengths and different frequencies of the special pair marker mode, ωsp, for Rb. sphaeroides that could be used in the future for more advanced calculations of absorption and HB spectra obtained for various bacterial reaction centers.

Persistent/transient spectral hole burning (HB) and computer simulations are used to provide new insight into the excitonic structure and excitation energy transfer of the widely studied bacterial reaction center (bRC) of Rhodobacter (Rb.) sphaeroides. We focus on site energies of its cofactors and electrochromic shifts induced in the chemically oxidized (P(+)) and charge-separated (P(+)QM(-)) states. Theoretical models lead to two alternative interpretations of the H-band. On the basis of our experimental and simulation data, we suggest that the bleach near 813-825 nm in transient HB spectra in the P(+)QM(-) state, often assigned to the upper exciton component of the special pair, is mostly due to different electrochromic shifts of the BL/M cofactors. From the exciton compositions in the charge-neutral (CN) bRC, the weak fourth excitonic band near 780 nm can be denoted PY+, that is, the upper excitonic band of the special pair, which in the CN bRC behaves as a delocalized state over PM and PL pigments that weakly mixes with accessory BChls. Thus, the shoulder in the absorption of Rb. sphaeroides near 813-815 nm does not contain the PY+ exciton band.

No abstract