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

Nonadiabatic Dynamics via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity

The symmetrical quasi-classical (SQC) method recently proposed by Miller and Cotton allows one to simulate nonadiabatic dynamics based on an algorithm with classical-like scaling with respect to system size. This is made possible by casting the electronic degrees of freedom in terms of mapping variables that can be propagated in a classical-like manner. While SQC was shown to be rather accurate when applied to benchmark models with harmonic electronic potential energy surfaces, it was also found to become inaccurate and to suffer numerical instabilities when applied to anharmonic systems. In this paper, we propose an extended SQC (E-SQC) methodology for overcoming those discrepancies by describing the anharmonic nuclear modes, which are coupled to the electronic degrees of freedom, in terms of classical-like mapping variables. The accuracy of E-SQC relative to standard SQC is demonstrated on benchmark models with quartic and Morse potential energy surfaces.

The modeling of the absorption lineshape for embedded molecules through a polarizable QM/MM approach

We present a polarizable QM/MM strategy to simulate the absorption line shape of chromophores embedded in complex matrices, including both homogeneous and inhomogeneous broadening.

How fine-tuned for energy transfer is the environmental noise produced by proteins around biological chromophores?

Analysis of intermolecular motions of pigment–protein complexes shows no significant difference in influence of local environment despite different biological functions.

Absence of Selection for Quantum Coherence in the Fenna-Matthews-Olson Complex: A Combined Evolutionary and Excitonic Study

We present a study on the evolution of the Fenna-Matthews-Olson bacterial photosynthetic pigment-protein complex. This protein complex functions as an antenna. It transports absorbed photons-excitons-to a reaction center where photosynthetic reactions initiate. The efficiency of exciton transport is therefore fundamental for the photosynthetic bacterium's survival. We have reconstructed an ancestor of the complex to establish whether coherence in the exciton transport was selected for or optimized over time. We have also investigated the role of optimizing free energy variation upon folding in evolution. We studied whether mutations which connect the ancestor to current day species were stabilizing or destabilizing from a thermodynamic viewpoint. From this study, we established that most of these mutations were thermodynamically neutral. Furthermore, we did not see a large change in exciton transport efficiency or coherence, and thus our results predict that exciton coherence was not specifically selected for.

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.

Observation of Electronic Excitation Transfer Through Light Harvesting Complex II Using Two-Dimensional Electronic-Vibrational Spectroscopy

Light-harvesting complex II (LHCII) serves a central role in light harvesting for oxygenic photosynthesis and is arguably the most important photosynthetic antenna complex. In this work, we present two-dimensional electronic-vibrational (2DEV) spectra of LHCII isolated from spinach, demonstrating the possibility of using this technique to track the transfer of electronic excitation energy between specific pigments within the complex. We assign the spectral bands via comparison with the 2DEV spectra of the isolated chromophores, chlorophyll a and b, and present evidence that excitation energy between the pigments of the complex are observed in these spectra. Finally, we analyze the essential components of the 2DEV spectra using singular value decomposition, which makes it possible to reveal the relaxation pathways within this complex.

Modeling Electronic-Nuclear Interactions for Excitation Energy Transfer Processes in Light-Harvesting Complexes

An accurate approach for computing intermolecular and intrachromophore contributions to spectral densities to describe the electronic-nuclear interactions relevant for modeling excitation energy transfer processes in light harvesting systems is presented. The approach is based on molecular dynamics (MD) calculations of classical correlation functions of long-range contributions to excitation energy fluctuations and a separate harmonic analysis and single-point gradient quantum calculations for electron-intrachromophore vibrational couplings. A simple model is also presented that enables detailed analysis of the shortcomings of standard MD-based excitation energy fluctuation correlation function approaches. The method introduced here avoids these problems, and its reliability is demonstrated in accurate predictions for bacteriochlorophyll molecules in the Fenna-Matthews-Olson pigment-protein complex, where excellent agreement with experimental spectral densities is found. This efficient approach can provide instantaneous spectral densities for treating the influence of fluctuations in environmental dissipation on fast electronic relaxation.

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.

Challenges facing an understanding of the nature of low-energy excited states in photosynthesis

While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.

On the Conflicting Estimations of Pigment Site Energies in Photosynthetic Complexes: A Case Study of the CP47 Complex

We focus on problems with elucidation of site energies (E0n) for photosynthetic complexes (PSCs) in order to raise some genuine concern regarding the conflicting estimations propagating in the literature. As an example, we provide a stern assessment of the site energies extracted from fits to optical spectra of the widely studied CP47 antenna complex of photosystem II from spinach, though many general comments apply to other PSCs as well. Correct values of E0n for chlorophyll (Chl) a in CP47 are essential for understanding its excitonic structure, population dynamics, and excitation energy pathway(s). To demonstrate this, we present a case study where simultaneous fits of multiple spectra (absorption, emission, circular dichroism, and nonresonant hole-burned spectra) show that several sets of parameters can fit the spectra very well. Importantly, we show that variable emission maxima (690–695 nm) and sample-dependent bleaching in nonresonant hole-burning spectra reported in literature could be explained, assuming that many previously studied CP47 samples were a mixture of intact and destabilized proteins. It appears that the destabilized subpopulation of CP47 complexes could feature a weakened hydrogen bond between the 13¹-keto group of Chl29 and the PsbH protein subunit, though other possibilities cannot be entirely excluded, as discussed in this work. Possible implications of our findings are briefly discussed.

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.

Nonadiabatic Dynamics in Atomistic Environments: Harnessing Quantum-Classical Theory with Generalized Quantum Master Equations

The development of methods that can efficiently and accurately treat nonadiabatic dynamics in quantum systems coupled to arbitrary atomistic environments remains a significant challenge in problems ranging from exciton transport in photovoltaic materials to electron and proton transfer in catalysis. Here we show that our recently introduced MF-GQME approach, which combines Ehrenfest mean field theory with the generalized quantum master equation framework, is able to yield quantitative accuracy over a wide range of charge-transfer regimes in fully atomistic environments. This is accompanied by computational speed-ups of up to 3 orders of magnitude over a direct application of Ehrenfest theory. This development offers the opportunity to efficiently investigate the atomistic details of nonadiabatic quantum relaxation processes in regimes where obtaining accurate results has previously been elusive.

Hybrid QM/MM study of FMO complex with polarized protein-specific charge

The Fenna-Matthews-Olson (FMO) light-harvesting complex is now one of the primary model systems for the study of excitation energy transfer (EET). However, the mechanism of the EET in this system is still controversial. In this work, molecular dynamics simulations and the electrostatic-embedding quantum-mechanics/molecular-mechanics single-point calculations have been employed to predict the energy transfer pathways utilizing the polarized protein-specific charge (PPC), which provides a more realistic description of Coulomb interaction potential in the protein than conventional mean-field charge scheme. The recently discovered eighth pigment has also been included in this study. Comparing with the conventional mean-field charges, more stable structures of FMO complex were found under PPC scheme during molecular dynamic simulation. Based on the electronic structure calculations, an exciton model was constructed to consider the couplings during excitation. The results show that pigments 3 and 4 dominate the lowest exciton levels whereas the highest exciton level are mainly constituted of pigments 1 and 6. This observation agrees well with the assumption based on the spatial distribution of the pigments. Moreover, the obtained spectral density in this study gives a reliable description of the diverse local environment embedding each pigment.

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.