Mortaza Aghtar

From atomistic modeling to excitation transfer and two-dimensional spectra of the FMO light-harvesting complex

The experimental observation of long-lived quantum coherences in the Fenna-Matthews-Olson (FMO) light-harvesting complex at low temperatures has challenged general intuition in the field of complex molecular systems and provoked considerable theoretical effort in search of explanations. Here we report on room-temperature calculations of the excited-state dynamics in FMO using a combination of molecular dynamics simulations and electronic structure calculations. Thus we obtain trajectories for the Hamiltonian of this system which contains time-dependent vertical excitation energies of the individual bacteriochlorophyll molecules and their mutual electronic couplings. The distribution of energies and couplings is analyzed together with possible spatial correlations. It is found that in contrast to frequent assumptions the site energy distribution is non-Gaussian. In a subsequent step, averaged wave packet dynamics is used to determine the exciton dynamics in the system. Finally, with the time-dependent Hamiltonian, linear and two-dimensional spectra are determined. The thus-obtained linear absorption line shape agrees well with experimental observation and is largely determined by the non-Gaussian site energy distribution. The two-dimensional spectra are in line with what one would expect by extrapolation of the experimental observations at lower temperatures and indicate almost total loss of long-lived coherences.

Influence of Force Fields and Quantum Chemistry Approach on Spectral Densities of BChl a in Solution and in FMO Proteins

Studies on light-harvesting (LH) systems have attracted much attention after the finding of long-lived quantum coherences in the exciton dynamics of the Fenna-Matthews-Olson (FMO) complex. In this complex, excitation energy transfer occurs between the bacteriochlorophyll a (BChl a) pigments. Two quantum mechanics/molecular mechanics (QM/MM) studies, each with a different force-field and quantum chemistry approach, reported different excitation energy distributions for the FMO complex. To understand the reasons for these differences in the predicted excitation energies, we have carried out a comparative study between the simulations using the CHARMM and AMBER force field and the Zerner intermediate neglect of differential orbital (ZINDO)/S and time-dependent density functional theory (TDDFT) quantum chemistry methods. The calculations using the CHARMM force field together with ZINDO/S or TDDFT always show a wider spread in the energy distribution compared to those using the AMBER force field. High- or low-energy tails in these energy distributions result in larger values for the spectral density at low frequencies. A detailed study on individual BChl a molecules in solution shows that without the environment, the density of states is the same for both force field sets. Including the environmental point charges, however, the excitation energy distribution gets broader and, depending on the applied methods, also asymmetric. The excitation energy distribution predicted using TDDFT together with the AMBER force field shows a symmetric, Gaussian-like distribution.

Different Types of Vibrations Interacting with Electronic Excitations in Phycoerythrin 545 and Fenna–Matthews–Olson Antenna Systems

The interest in the phycoerythrin 545 (PE545) photosynthetic antenna system of marine algae and the Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria has drastically increased since long-lived quantum coherences were reported for these complexes. For the PE545 complex, this phenomenon is clearly visible even at ambient temperatures, while for the FMO system it is more prominent at lower temperatures. The key to elucidate the role of the environment in these long-lived quantum effects is the spectral density. Here, we employ molecular dynamics simulations combined with quantum chemistry calculations to study the coupling between the biological environment and the vertical excitation energies of the bilin pigment molecules in PE545 and compare them to prior calculations on the FMO complex. It is found that the overall strength of the resulting spectral densities for the PE545 system is similar to the experiment-based counterpart but also to those in the FMO complex. Molecular analysis, however, reveals that the origin for the spectral densities in the low frequency range, which is most important for excitonic transitions, is entirely different. In the case of FMO, this part of the spectral density is due to environmental fluctuations, while, in case of PE545, it is essentially only due to internal modes of the bilin molecules. This finding sheds new light on possible explanations of the long-lived quantum coherences and that the reasons might actually be different in dissimilar systems.

The FMO Complex in a Glycerol–Water Mixture

Experimental findings of long-lived quantum coherence in the Fenna-Matthews-Olson (FMO) complex and other photosynthetic complexes have led to theoretical studies searching for an explanation of this unexpected phenomenon. Extending in this regard our own earlier calculations, we performed simulations of the FMO complex in a glycerol-water mixture at 310 K as well as 77 K, matching the conditions of earlier 2D spectroscopic experiments by Engel et al. The calculations, based on an improved quantum procedure employed by us already, yielded spectral densities of each individual pigment of FMO, in water and glycerol-water solvents at ambient temperature that compare well to prior experimental estimates. Due to the slow solvent dynamics at 77 K, the present results strongly indicate the presence of static disorder, i.e., disorder on a time scale beyond that relevant for the construction of spectral densities.

Juxtaposing density matrix and classical path-based wave packet dynamics

In many physical, chemical, and biological systems energy and charge transfer processes are of utmost importance. To determine the influence of the environment on these transport processes, equilibrium molecular dynamics simulations become more and more popular. From these simulations, one usually determines the thermal fluctuations of certain energy gaps, which are then either used to perform ensemble-averaged wave packet simulations, also called Ehrenfest dynamics, or to employ a density matrix approach via spectral densities. These two approaches are analyzed through energy gap fluctuations that are generated to correspond to a predetermined spectral density. Subsequently, density matrix and wave packet simulations are compared through population dynamics and absorption spectra for different parameter regimes. Furthermore, a previously proposed approach to enforce the correct long-time behavior in the wave packet simulations is probed and an improvement is proposed.

Impact of Electronic Fluctuations and Their Description on the Exciton Dynamics in the Light-Harvesting Complex PE545

Temperature-dependent fluctuations of both site energies and electronic couplings are known to affect the excitation energy transfer in light-harvesting complexes. Environment effects on such fluctuations as well as possible spatial correlations among them are investigated here in the PE545 complex from cryptophyte algae using ensemble-averaged wave packet dynamics to extract the exciton dynamics. This strategy directly uses the time-dependent fluctuations of the system Hamiltonian, as described by quantum mechanics/molecular mechanics calculations performed along a classical MD trajectory. Neither the fluctuations in the couplings nor spatial correlations including cross-correlations between site energies and couplings are found to be important in the exciton dynamics of the complex. This finding does not change if a polarizable embedding is used instead of its electrostatic counterpart. The impact of variations in spectral densities and screening of excitonic couplings based on the electrostatic and polarizable embeddings are discussed as well.

Relation between Dephasing Time and Energy Gap Fluctuations in Biomolecular Systems.

Excitation energy and charge transfer are fundamental processes in biological systems. Because of their quantum nature, the effect of dephasing on these processes is of interest especially when trying to understand their efficiency. Moreover, recent experiments have shown quantum coherences in such systems. As a first step toward a better understanding, we studied the relationship between dephasing time and energy gap fluctuations of the individual molecular subunits. A larger set of molecular simulations has been investigated to shed light on this dependence. This set includes bacterio-chlorophylls in Fenna-Matthews-Olson complexes, the PE545 aggregate, the LH2 complexes, DNA, photolyase, and cryptochromes. For the individual molecular subunits of these aggregates it has been confirmed quantitatively that an inverse proportionality exists between dephasing time and average gap energy fluctuation. However, for entire complexes including the respective intermolecular couplings, such a relation still needs to be verified.