A remarkable amount of theoretical research has been carried out to elucidate the physical origins of the recently observed long-lived quantum coherence in the electronic energy transfer process in biological photosynthetic systems. Although successful in many respects, several widely used descriptions only include an effective treatment of the protein-chromophore interactions. In this work, by combining an all-atom molecular dynamics simulation, time-dependent density functional theory, and open quantum system approaches, we successfully simulate the dynamics of the electronic energy transfer of the Fenna-Matthews-Olson pigment-protein complex. The resulting characteristic beating of populations and quantum coherences is in good agreement with the experimental results and the hierarchy equation of motion approach. The experimental absorption, linear, and circular dichroism spectra and dephasing rates are recovered at two different temperatures. In addition, we provide an extension of our method to include zero-point fluctuations of the vibrational environment. This work thus presents, to our knowledge, one of the first steps to explain the role of excitonic quantum coherence in photosynthetic light-harvesting complexes based on their atomistic and molecular description.

Organic molecules store the energy of absorbed light in the form of charge-neutral molecular excitationsFrenkel excitons. Usually, in amorphous organic materials, excitons are viewed as quasiparticles, localized on single molecules, which diffuse randomly through the structure. However, the picture of incoherent hopping is not applicable to some classes of molecular aggregates -assemblies of molecules that have strong near-field interaction between electronic excitations in the individual subunits. Molecular aggregates can be found in nature, in photosynthetic complexes of plants and bacteria, and they can also be produced artificially in various forms including quasi-one dimensional chains, two-dimensional films, tubes, etc. In these structures light is absorbed collectively by many mole cules and the following dynamics of molecular excitation possesses coherent properties. This energy transfer mechanism, mediated by the coherent exciton dynamics, resembles the propagation of electromagnetic waves through a structured medium on the nanometer scale. The absorbed energy can be transferred resonantly over distances of hundreds of nanometers before exciton relaxation occurs. Furthermore, the spatial and energetic landscape of molecular aggregates can enable the funneling of the exciton energy to a small number of molecules either within or outside the aggregate. In this review we establish a bridge between the fields of photo nics and excitonics by describing the present understanding of exciton dynamics in molecular aggregates.

We investigate on the procedure of extracting a "spectral density" from mixed QM/MM calculations and employing it in open quantum systems models. In particular, we study the connection between the energy gap correlation function extracted from ground state QM/MM and the bath spectral density used as input in open quantum system approaches. We introduce a simple model which can give intuition on when the ground state QM/MM propagation will give the correct energy gap. We also discuss the role of higher order correlators of the energy-gap fluctuations which can provide useful information on the bath. Further, various semiclassical corrections to the spectral density, are applied and investigated. Finally, we apply our considerations to the photosynthetic Fenna-Matthews-Olson complex. For this system, our results suggest the use of the Harmonic prefactor for the spectral density rather than the Standard one, which was employed in the simulations of the system carried out to date.

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.

Machine learning ground state QM/MM for accelerated computation of exciton dynamics.