The FMO complex in a glycerol-water mixture

Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex

A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.

Multiscale QM/MM molecular dynamics simulations of the trimeric major light-harvesting complex II

Photosynthetic processes are driven by sunlight. Too little of it and the photosynthetic machinery cannot produce the reductive power to drive the anabolic pathways. Too much sunlight and the machinery can get damaged. In higher plants, the major Light-Harvesting Complex (LHCII) efficiently absorbs the light energy, but can also dissipate it when in excess (quenching). In order to study the dynamics related to the quenching process but also the exciton dynamics in general, one needs to accurately determine the so-called spectral density which describes the coupling between the relevant pigment modes and the environmental degrees of freedom. To this end, Born-Oppenheimer molecular dynamics simulations in a quantum mechanics/molecular mechanics (QM/MM) fashion utilizing the density functional based tight binding (DFTB) method have been performed for the ground state dynamics. Subsequently, the time-dependent extension of the long-range-corrected DFTB scheme has been employed for the excited state calculations of the individual chlorophyll-a molecules in the LHCII complex. The analysis of this data resulted in spectral densities showing an astonishing agreement with the experimental counterpart in this rather large system. This consistency with an experimental observable also supports the accuracy, robustness, and reliability of the present multi-scale scheme. To the best of our knowledge, this is the first theoretical attempt on this large complex system is ever made to accurately simulate the spectral density. In addition, the resulting spectral densities and site energies were used to determine the exciton transfer rate within a special pigment pair consisting of a chlorophyll-a and a carotenoid molecule which is assumed to play a role in the balance between the light harvesting and quenching modes.

DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex

Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum-classical description. To this end, Born-Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna-Matthews-Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.

Intramolecular vibrations enhance the quantum efficiency of excitonic energy transfer

We study the impact of underdamped intramolecular vibrational modes on the efficiency of the excitation energy transfer in a dimer in which each state is coupled to its own underdamped vibrational mode and, in addition, to a continuous background of environmental modes. For this, we use the numerically exact hierarchy equation of motion approach. We determine the quantum yield and the transfer time in dependence of the vibronic coupling strength, and in dependence of the damping of the incoherent background. Moreover, we tune the vibrational frequencies out of resonance with the excitonic energy gap. We show that the quantum yield is enhanced by up to 10% when the vibrational frequency of the donor is larger than at the acceptor. The vibronic energy eigenstates of the acceptor acquire then an increased density of states, which leads to a higher occupation probability of the acceptor in thermal equilibrium. We can conclude that an underdamped vibrational mode which is weakly coupled to the dimer fuels a faster transfer of excitation energy, illustrating that long-lived vibrations can, in principle, enhance energy transfer, without involving long-lived electronic coherence.

Effect of disorder and polarization sequences on two-dimensional spectra of light-harvesting complexes

Two-dimensional electronic spectra (2DES) provide unique ways to track the energy transfer dynamics in light-harvesting complexes. The interpretation of the peaks and structures found in experimentally recorded 2DES is often not straightforward, since several processes are imaged simultaneously. The choice of specific pulse polarization sequences helps to disentangle the sometimes convoluted spectra, but brings along other disturbances. We show by detailed theoretical calculations how 2DES of the Fenna-Matthews-Olson complex are affected by rotational and conformational disorder of the chromophores.

Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation

Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.

Relation between Vibrational Dephasing Time and Energy Gap Fluctuations

Dephasing processes are present in basically all applications in which quantum mechanics plays a role. These applications certainly include excitation energy and charge transfer in biological systems. In a previous study, we have analyzed the vibrational dephasing time as a function of energy gap fluctuation for a large set of molecular simulations. In that investigation, individual molecular subunits were the focus of the calculations. The set of studied molecules included bacteriochlorophylls in Fenna-Matthews-Olson and light-harvesting system 2 complexes as well as bilins in PE545 aggregates. The present work extends this study to entire complexes, including the respective intermolecular couplings. Again, it can be concluded that a universal and inverse proportionality exists between dephasing time and variance of the excitonic energy gap fluctuations, whereas the respective proportionality constants can be rationalized using the energy gap autocorrelation functions. Furthermore, these findings can be extended to the gaps between higher-lying neighboring excitonic states.

Machine learning exciton dynamics

Obtaining the exciton dynamics of large photosynthetic complexes by using mixed quantum mechanics/molecular mechanics (QM/MM) is computationally demanding. We propose a machine learning technique, multi-layer perceptrons, as a tool to reduce the time required to compute excited state energies. With this approach we predict time-dependent density functional theory (TDDFT) excited state energies of bacteriochlorophylls in the Fenna-Matthews-Olson (FMO) complex. Additionally we compute spectral densities and exciton populations from the predictions. Different methods to determine multi-layer perceptron training sets are introduced, leading to several initial data selections. In addition, we compute spectral densities and exciton populations. Once multi-layer perceptrons are trained, predicting excited state energies was found to be significantly faster than the corresponding QM/MM calculations. We showed that multi-layer perceptrons can successfully reproduce the energies of QM/MM calculations to a high degree of accuracy with prediction errors contained within 0.01 eV (0.5%). Spectral densities and exciton dynamics are also in agreement with the TDDFT results. The acceleration and accurate prediction of dynamics strongly encourage the combination of machine learning techniques with ab initio methods.

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.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Combining classical molecular dynamics and quantum mechanical methods for the description of electronic excitations: The case of carotenoids

Carotenoids are important actors both in light-harvesting (LH) and in photoprotection functions of photosynthetic pigment-protein complexes. A deep theoretical investigation of this multiple role is still missing owing to the difficulty of describing the delicate interplay between electronic and nuclear degrees of freedom. A possible strategy is to combine accurate quantum mechanical (QM) methods with classical molecular dynamics. To do this, however, accurate force-fields (FF) are necessary. This article presents a new FF for the different carotenoids present in LH complexes of plants. The results show that all the important structural properties described by the new FF are in very good agreement with QM reference values. This increased accuracy in the simulation of the structural fluctuations is also reflected in the description of excited states. Both the energy order and the different nature of the lowest singlet states are preserved during the dynamics when the new FF is used, whereas an unphysical mixing is found when a standard FF is used.

Open quantum system parameters for light harvesting complexes from molecular dynamics

We elucidate the difference between various parameter extraction methods and demonstrate sensitivity to molecular dynamics equilibration.

Limits and potentials of quantum chemical methods in modelling photosynthetic antennae

A critical overview of quantum chemical approaches to simulate the light-harvesting process in photosynthetic antennae is presented together with a perspective on the developments that need to be introduced to reach a quantitative predictive power.

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.

Molecular basis of the exciton–phonon interactions in the PE545 light-harvesting complex

A fully polarizable QM/MM approach is used in combination with classical MD simulations to predict the pigment-dependent spectral densities of the PE545 antenna complex and account for their effects on the exciton dynamics.

Quantification of non-Markovian effects in the Fenna-Matthews-Olson complex

The excitation energy transfer dynamics in the Fenna-Matthews-Olson complex is quantified in terms of a non-Markovianity measure based on the time evolution of the trace distance of two quantum states. We use a system description derived from experiments and different environmental fluctuation spectral functions, which are obtained either from experimental data or from molecular dynamics simulations. These exhibit, in all cases, a nontrivial structure with several peaks attributed to vibrational modes of the pigment-protein complex. Such a structured environmental spectrum can, in principle, give rise to strong non-Markovian effects. We present numerically exact real-time path-integral calculations for the transfer dynamics and find, in all cases, a monotonic decrease of the trace distance with increasing time which renders a Markovian description valid.

Disentangling electronic and vibronic coherences in two-dimensional echo spectra

The prevalence of long-lasting oscillatory signals in two-dimensional (2D) echo spectroscopy of light-harvesting complexes has led to a search for possible mechanisms. We investigate how two causes of oscillatory signals are intertwined: (i) electronic coherences supporting delocalized wavelike motion and (ii) narrow bands in the vibronic spectral density. To disentangle the vibronic and electronic contributions, we introduce a time-windowed Fourier transform of the signal amplitude. We find that 2D spectra can be dominated by excitations of pathways which are absent in excitonic energy transport. This leads to an underestimation of the lifetime of electronic coherences by 2D spectra.