An Ab Initio Description of the Excitonic Properties of LH2 and Their Temperature Dependence

Exciton Transfer Simulations in a Light-Harvesting 2 Complex Reveal the Transient Delocalization Mechanism

The striking efficiency of exciton transfer in light-harvesting (LH) complexes has remained a topic of debate since the revision of the long-held role of electronic coherences. To address this issue, we have developed a neural network for the pigments in the LH2 complex of Rhodospirillum molischianum that allows nonadiabatic molecular dynamic (NAMD) simulations of exciton transfer in a coupled quantum mechanical/molecular mechanics (QM/MM) embedding. The calculated exciton occupations are averaged over hundreds of trajectories, each lasting several picoseconds. We have obtained transitions within the B800 and B850 rings that agree well with the experimental results, indicating an incoherent hopping process in the B800 ring and a more delocalized transfer in the B850 subsystem. The reorganization energies and excitonic couplings are comparable to each other, indicating that the "transient delocalization" transport model is the underlying cause of the highly efficient exciton transport in the B850 ring. This phenomenon can be attributed to a localized exciton that shows occasional large delocalization events. Our results indicate that the reason for the striking efficiency is the unusual electronic property of bacteriochlorophyll, manifested in minimal inner and outer sphere reorganization energies.

Theory of photosynthetic membrane influence on B800-B850 energy transfer in the LH2 complex

Photosynthetic organisms rely on a network of light-harvesting protein-pigment complexes to efficiently absorb sunlight and transfer excitation energy to reaction center proteins where charge separation occurs. In photosynthetic purple bacteria, these complexes are embedded within the cell membrane, with lipid composition affecting complex clustering, thereby impacting inter-complex energy transfer. However, the impact of the lipid bilayer on intra-complex excitation dynamics is less understood. Recent experiments have addressed this question by comparing photo-excitation dynamics in detergent-isolated light-harvesting complex 2 (LH2) to LH2 complexes embedded in membrane discs mimicking the biological environment, revealing differences in spectra and energy-transfer rates. In this paper, we use available quantum chemical and spectroscopy data to develop a complementary theoretical study on the excitonic structure and intra-complex energy-transfer kinetics of the LH2 of photosynthetic purple bacteria Rhodoblastus (Rbl.) acidophilus (formerly Rhodopseudomonas acidophila) in two different conditions: the LH2 in a membrane environment and detergent-isolated LH2. We find that dark excitonic states, crucial for B800-B850 energy transfer within LH2, are more delocalized in the membrane model. Using nonperturbative and generalized Förster calculations, we show that such increased quantum delocalization results in a 30% faster B800 to B850 transfer rate in the membrane model, in agreement with experimental results. We identify the dominant energy-transfer pathways in each environment and demonstrate how differences in the B800 to B850 transfer rate arise from changes in LH2's electronic properties when embedded in the membrane. Furthermore, by accounting for the quasi-static variations of electronic excitation energies in the LH2, we show that the broadening of the distribution of the B800-B850 transfer rates is affected by the lipid composition. We argue that such variation in broadening could be a signature of a speed-accuracy trade-off, commonly seen in biological process.

Multiconfigurational Excitonic Couplings in Homo- and Heterodimer Stacks of Azobenzene-Derived Dyes

Molecular excitons play a major role within dye aggregates and hold significant potential for (opto)electronic and photovoltaic applications. Numerous studies have documented alterations in the spectral properties of dye homoaggregates, but only limited work has been reported for heteroaggregates. In this article, dimeric dye stacks were constructed from azobenzene-like dyes with identical or distinct structures, and their excitonic features were computationally investigated. Our results show that strong exciton coupling is not limited to identical chromophores, as often assumed, based on a recently made available Frenkel Exciton Hamiltonian and multiconfigurational plus second-order perturbation theory energetics methodology. Heteroaggregate stacks were found to exhibit different absorption features from the corresponding interacting monomers, indicating considerable coupling interactions between units. We analyzed how such coupling may vary according to various aspects, such as the relative positions of the interacting monomers or the differences in their energetics. Such qualitative and semiquantitative analyses allow the evaluation of the excitonic behavior of these dye aggregates to encourage further efforts toward a deeper understanding of the excitonic properties of tailored dye heteroaggregate systems.

Tuning by Hydrogen Bonding in Photosynthesis

Hydrogen bonding plays a crucial role in stabilizing proteins throughout their folding process. In photosynthetic light-harvesting chromoproteins, enriched with pigment chromophores, hydrogen bonds also fine-tune optical absorption to align with the solar irradiation spectrum. Despite its significance for photosynthesis, the precise mechanism of spectral tuning through hydrogen bonding remains inadequately understood. This study investigates wild-type and genetically engineered LH2 and LH1 light-harvesting complexes from Rhodobacter sphaeroides using a unique set of advanced spectroscopic techniques combined with simple exciton modeling. Our findings reveal an intricate interplay between exciton and site energy shift mechanisms, challenging the prevailing belief that spectral changes observed in these complexes upon the modification of tertiary structure hydrogen bonds almost directly follow shifting site energies. These deeper insights into natural adaptation processes hold great promise for advancing sustainable solar energy conversion technologies.

Prominent Role of Charge Transfer in the Spectral Tuning of Photosynthetic Light-Harvesting I Complex

Purple bacteria possess two ring-shaped protein complexes, light-harvesting 1 (LH1) and 2 (LH2), both of which function as antennas for solar energy utilization for photosynthesis but exhibit distinct absorption properties. The two antennas have differing amounts of bacteriochlorophyll (BChl) a; however, their significance in spectral tuning remains elusive. Here, we report a high-precision evaluation of the physicochemical factors contributing to the variation in absorption maxima between LH1 and LH2, namely, BChl a structural distortion, protein electrostatic interaction, excitonic coupling, and charge transfer (CT) effects, as derived from detailed spectral calculations using an extended version of the exciton model, in the model purple bacterium Rhodospirillum rubrum. Spectral analysis confirmed that the electronic structure of the excited state in LH1 extended to the BChl a 16-mer. Further analysis revealed that the LH1-specific redshift (∼61% in energy) is predominantly accounted for by the CT effect resulting from the closer inter-BChl distance in LH1 than in LH2. Our analysis explains how LH1 and LH2, both with chemically identical BChl a chromophores, use distinct physicochemical effects to achieve a progressive redshift from LH2 to LH1, ensuring efficient energy transfer to the reaction center special pair.

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Identification of Residues Potentially Involved in Optical Shifts in the Water-Soluble Chlorophyll a-Binding Protein through Molecular Dynamics Simulations

Reversible light and thermally induced spectral shifts are universally observed in a wide variety of pigment-protein complexes at temperatures ranging from cryogenic to ambient. In this paper, we employed large-scale molecular dynamics (MD) simulations of a prototypical pigment-protein complex to better understand these shifts at a molecular scale. Although multiple mechanisms have been proposed over the years, no verification of these proposals via MD simulations has thus far been performed; our work represents the first step in this direction. From simulations of the water-soluble chlorophyll-binding protein complex, we determined that rearrangements of long hydrogen bonds were unlikely to be the origin of the multiwell landscape features necessary to explain observed spectral shifts. We also assessed small motions of amino acid residues and identified side chain rotations of some of these residues as likely candidates for the origin of relevant multiwell landscape features. The protein free-energy landscapes associated with side chain rotations feature energy barriers of around 1100-1600 cm-1, in agreement with optical spectroscopy results, with the most promising residue type associated with experimental signatures being serine, which possesses a symmetric triple-well landscape and moment of inertia of a relevant magnitude.

Excitation Energy Transfer between Higher Excited States of Photosynthetic Pigments: 2. Chlorophyll b is a B Band Excitation Trap

Chlorophylls (Chls) are known for fast, subpicosecond internal conversion (IC) from ultraviolet/blue absorbing ("B" or "Soret" states) to the energetically lower, red light-absorbing Q states. Consequently, excitation energy transfer (EET) in photosynthetic pigment-protein complexes involving the B states has so far not been considered. We present, for the first time, a theoretical framework for the existence of B-B EET in tightly coupled Chl aggregates such as photosynthetic pigment-protein complexes. We show that according to a Förster resonance energy transport (FRET) scheme, unmodulated B-B EET has an unexpectedly high range. Unsuppressed, it could pose an existential threat-the damage potential of blue light for photochemical reaction centers (RCs) is well-known. This insight reveals so-far undescribed roles for carotenoids (Crts, cf. previous article in this series) and Chl b (this article) of possibly vital importance. Our model system is the photosynthetic antenna pigment-protein complex (CP29). The focus of the study is on the role of Chl b for EET in the Q and B bands. Further, the initial excited pigment distribution in the B band is computed for relevant solar irradiation and wavelength-centered laser pulses. It is found that both accessory pigment classes compete efficiently with Chl a absorption in the B band, leaving only 40% of B band excitations for Chl a. B state population is preferentially relocated to Chl b after excitation of any Chls, due to a near-perfect match of Chl b B band absorption with Chl a B state emission spectra. This results in an efficient depletion of the Chl a population (0.66 per IC/EET step, as compared to 0.21 in a Chl a-only system). Since Chl b only occurs in the peripheral antenna complexes of plants and algae, and RCs contain only Chl a, this would automatically trap potentially dangerous B state population in the antennae, preventing forwarding to the RCs.

A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates

We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

Quantum chemical elucidation of a sevenfold symmetric bacterial antenna complex

The light-harvesting complex 2 (LH2) of purple bacteria is one of the most studied photosynthetic antenna complexes. Its symmetric structure and ring-like bacteriochlorophyll arrangement make it an ideal system for theoreticians and spectroscopists. LH2 complexes from most bacterial species are thought to have eightfold or ninefold symmetry, but recently a sevenfold symmetric LH2 structure from the bacterium Mch. purpuratum was solved by Cryo-Electron microscopy. This LH2 also possesses unique near-infrared absorption and circular dichroism (CD) spectral properties. Here we use an atomistic strategy to elucidate the spectral properties of Mch. purpuratum LH2 and understand the differences with the most commonly studied LH2 from Rbl. acidophilus. Our strategy exploits a combination of molecular dynamics simulations, multiscale polarizable quantum mechanics/molecular mechanics calculations, and lineshape simulations. Our calculations reveal that the spectral properties of LH2 complexes are tuned by site energies and exciton couplings, which in turn depend on the structural fluctuations of the bacteriochlorophylls. Our strategy proves effective in reproducing the absorption and CD spectra of the two LH2 complexes, and in uncovering the origin of their differences. This work proves that it is possible to obtain insight into the spectral tuning strategies of purple bacteria by quantitatively simulating the spectral properties of their antenna complexes.

Recent progress in atomistic modeling of light-harvesting complexes: a mini review

In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.

Phasor Analysis Reveals Multicomponent Fluorescence Kinetics in the LH2 Complex from Marichromatium purpuratum

We investigated the fluorescence kinetics of LH2 complexes from Marichromatium purpuratum, the cryo-EM structure of which has been recently elucidated with 2.4 Å resolution. The experiments have been carried out as a function of the excitation density by varying both the excitation fluence and the repetition rate of the laser excitation. Instead of the usual multiexponential fitting procedure, we applied the less common phasor formalism for evaluating the transients because this allows for a model-free analysis of the data without a priori knowledge about the number of processes that contribute to a particular decay. For the various excitation conditions, this analysis reproduces consistently three lifetime components with decay times below 100 ps, 500 ps, and 730 ps, which were associated with the quenched state, singlet-triplet annihilation, and fluorescence decay, respectively. Moreover, it reveals that the number of decay components that contribute to the transients depends on whether the excitation wavelength is in resonance with the B800 BChl a molecules or with the carotenoids. Based on the mutual arrangement of the chromophores in their binding pockets, this leads us to conclude that the energy transfer pathways within the LH2 complex of this species differ significantly from each other for exciting either the B800 BChl molecules or the carotenoids. Finally, we speculate whether the illumination with strong laser light converts the LH2 complexes studied here into a quenched conformation that might be related to the development of the non-photochemical quenching mechanism that occurs in higher plants.

The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol

The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

Bridging Carotenoid-to-Bacteriochlorophyll Energy Transfer of Purple Bacteria LH2 With Temperature Variations: Insights From Conformational Changes

Photosynthesis is a key process for converting light energy into chemical energy and providing food for lives on Earth. Understanding the mechanism for the energy transfers could provide insights into regulating energy transfers in photosynthesis and designing artificial photosynthesis systems. Many efforts have been devoted to exploring the mechanism of temperature variations affecting the excitonic properties of LH2. In this study, we performed all-atom molecular dynamics (MD) simulations and quantum mechanics calculations for LH2 complex from purple bacteria along with its membrane environment under three typical temperatures: 270, 300, and 330 K. The structural analysis from validated MD simulations showed that the higher temperature impaired interactions at N-terminus of both α and β polypeptide helices and led to the dissociation of this hetero polypeptide dimer. Rhodopin-β-D-glucosides (RG1) moved centripetally with α polypeptide helices when temperature increased and enlarged their distances with bacteriochlorophylls molecules that have the absorption peak at 850 nm (B850), which resulted in reducing the coupling strengths between RG1 and B850 molecules. The present study reported a cascading mechanism for temperature regulating the energy transfers in LH2 of purple bacteria.

Spectral Tuning Mechanism of Photosynthetic Light-Harvesting Complex II Revealed by Ab Initio Dimer Exciton Model

Excited states of two kinds of bacteriochlorophyll (BChl) aggregates, B850 and B800, in photosynthetic light-harvesting complex II (LH2) are theoretically investigated by developing and using an extended exciton model considering efficiently evaluated excitonic coupling. Our exciton model based on dimer fragmentation is shown to reproduce the experimental absorption spectrum of LH2 with good accuracy, entailing their different redshifts originating from aggregations of B850 and B800. The systematic analysis has been performed on the spectra by quantitatively decomposing their spectral shift energies into the contributions of various effects: structural distortion, electrostatic, excitonic coupling, and charge-transfer (CT) effects. Our results show that the spectral redshift of B800 is mainly attributed to its electrostatic interaction with the protein environment, while that of B850 arises from the marked effect of the excitonic coupling between BChl units. The interchromophore CT excitation also plays a key role in the spectral redshift of B850. This CT effect can be effectively described using our dimer model. This suited characterization reveals that the pronounced CT effect originates from the characteristics of B850 that has closely spaced BChls as dimers. We highlight the importance of the refinement of the crystal structure with the use of quantum chemical methods for prediction of the spectrum.

Photosynthetic Light-Harvesting (Antenna) Complexes-Structures and Functions

Chlorophylls and bacteriochlorophylls, together with carotenoids, serve, noncovalently bound to specific apoproteins, as principal light-harvesting and energy-transforming pigments in photosynthetic organisms. In recent years, enormous progress has been achieved in the elucidation of structures and functions of light-harvesting (antenna) complexes, photosynthetic reaction centers and even entire photosystems. It is becoming increasingly clear that light-harvesting complexes not only serve to enlarge the absorption cross sections of the respective reaction centers but are vitally important in short- and long-term adaptation of the photosynthetic apparatus and regulation of the energy-transforming processes in response to external and internal conditions. Thus, the wide variety of structural diversity in photosynthetic antenna “designs” becomes conceivable. It is, however, common for LHCs to form trimeric (or multiples thereof) structures. We propose a simple, tentative explanation of the trimer issue, based on the 2D world created by photosynthetic membrane systems.

Hybrid theoretical models for molecular nanoplasmonics

The multidisciplinary nature of the research in molecular nanoplasmonics, i.e., the use of plasmonic nanostructures to enhance, control, or suppress properties of molecules interacting with light, led to contributions from different theory communities over the years, with the aim of understanding, interpreting, and predicting the physical and chemical phenomena occurring at molecular- and nano-scale in the presence of light. Multiscale hybrid techniques, using a different level of description for the molecule and the plasmonic nanosystems, permit a reliable representation of the atomistic details and of collective features, such as plasmons, in such complex systems. Here, we focus on a selected set of topics of current interest in molecular plasmonics (control of electronic excitations in light-harvesting systems, polaritonic chemistry, hot-carrier generation, and plasmon-enhanced catalysis). We discuss how their description may benefit from a hybrid modeling approach and what are the main challenges for the application of such models. In doing so, we also provide an introduction to such models and to the selected topics, as well as general discussions on their theoretical descriptions.

Self-interaction correction, electrostatic, and structural influences on time-dependent density functional theory excitations of bacteriochlorophylls from the light-harvesting complex 2

First-principles calculations offer the chance to obtain a microscopic understanding of light-harvesting processes. Time-dependent density functional theory can have the computational efficiency to allow for such calculations. However, the (semi-)local exchange-correlation approximations that are computationally most efficient fail to describe charge-transfer excitations reliably. We here investigate whether the inexpensive average density self-interaction correction (ADSIC) remedies the problem. For the systems that we study, ADSIC is even more prone to the charge-transfer problem than the local density approximation. We further explore the recently reported finding that the electrostatic potential associated with the chromophores' protein environment in the light-harvesting complex 2 beneficially shifts spurious excitations. We find a great sensitivity on the chromophores' atomistic structure in this problem. Geometries obtained from classical molecular dynamics are more strongly affected by the spurious charge-transfer problem than the ones obtained from crystallography or density functional theory. For crystal structure geometries and density-functional theory optimized ones, our calculations confirm that the electrostatic potential shifts the spurious excitations out of the energetic range that is most relevant for electronic coupling.

Absorption and Circular Dichroism Spectra of Molecular Aggregates With the Full Cumulant Expansion

The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic techniques such as linear absorption and circular dichroism. To compare calculated Hamiltonians to experiments, a lineshape theory is needed, which takes into account the coupling of the excitons with inter- and intramolecular vibrations. This coupling is normally introduced in a perturbative way through the cumulant expansion formalism and further approximated by assuming a Markovian exciton dynamics, for example with the modified Redfield theory. Here, we present the implementation of the full cumulant expansion (FCE) formalism (J. Chem. Phys.142, 2015, 09410625747060) to efficiently compute absorption and circular dichroism spectra of molecular aggregates beyond the Markov approximation, without restrictions on the form of exciton–phonon coupling. By employing the LH2 system of purple bacteria as a challenging test case, we compare the FCE lineshapes with the Markovian lineshapes obtained with the modified Redfield theory, showing that the latter presents a less satisfying agreement with experiments. The FCE approach instead accurately describes the lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the FCE approach will become a valuable tool for accurately comparing model exciton Hamiltonians with optical spectroscopy experiments.

Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates

Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure–property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties.

Multiscale modeling resolves the molecular structure of a synthetic light-harvesting complex, unraveling the microscopic origin of its photophysical properties.

Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models

In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.

A comparative look at structural variation among RC-LH1 'Core' complexes present in anoxygenic phototrophic bacteria

All purple photosynthetic bacteria contain RC-LH1 'Core' complexes. The structure of this complex from Rhodobacter sphaeroides, Rhodopseudomonas palustris and Thermochromatium tepidum has been solved using X-ray crystallography. Recently, the application of single particle cryo-EM has revolutionised structural biology and the structure of the RC-LH1 'Core' complex from Blastochloris viridis has been solved using this technique, as well as the complex from the non-purple Chloroflexi species, Roseiflexus castenholzii. It is apparent that these structures are variations on a theme, although with a greater degree of structural diversity within them than previously thought. Furthermore, it has recently been discovered that the only phototrophic representative from the phylum Gemmatimonadetes, Gemmatimonas phototrophica, also contains a RC-LH1 'Core' complex. At present only a low-resolution EM-projection map exists but this shows that the Gemmatimonas phototrophica complex contains a double LH1 ring. This short review compares these different structures and looks at the functional significance of these variations from two main standpoints: energy transfer and quinone exchange.

Dynamic Disorder Drives Exciton Transfer in Tubular Chlorosomal Assemblies

Chlorosomes stand out for their highly efficient excitation energy transfer (EET) in extreme low light conditions. Yet, little is known about the EET when a chlorosome is excited to a pure state that is an eigenstate of the exciton Hamiltonian. In this work, we consider the dynamic disorder in the intermolecular electronic coupling explicitly by calculating the electronic coupling terms in the Hamiltonian using nuclear coordinates that are taken from molecular dynamics simulation trajectories. We show that this dynamic disorder is capable of driving the evolution of the exciton, being a stationary state of the initial Hamiltonian. In particular, long-distance excitation energy transfer between domains of high exciton population and oscillatory behavior of the population in the site basis are observed, in line with two-dimensional electronic spectroscopy studies. We also found that in the high exciton population domains, their population variation is correlated with their overall coupling strength. Analysis in a reference state basis shows that such dynamic disorder, originating from thermal energy, creates a fluctuating landscape for the exciton and promotes the EET process. We propose such dynamic disorder as an important microscopic origin for the high efficient EET widely observed in different types of chlorosomes, bioinspired tubular aggregates, or other light-harvesting complexes.

Molecular dynamics simulations in photosynthesis

Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

Benchmark and performance of long-range corrected time-dependent density functional tight binding (LC-TD-DFTB) on rhodopsins and light-harvesting complexes

The chromophores of rhodopsins (Rh) and light-harvesting (LH) complexes still represent a major challenge for a quantum chemical description due to their size and complex electronic structure. Since gradient corrected and hybrid density functional approaches have been shown to fail for these systems, only range-separated functionals seem to be a promising alternative to the more time consuming post-Hartree-Fock approaches. For extended sampling of optical properties, however, even more approximate approaches are required. Recently, a long-range corrected (LC) functional has been implemented into the efficient density functional tight binding (DFTB) method, allowing to sample the excited states properties of chromophores embedded into proteins using quantum mechanical/molecular mechanical (QM/MM) with the time-dependent (TD) DFTB approach. In the present study, we assess the accuracy of LC-TD-DFT and LC-TD-DFTB for rhodopsins (bacteriorhodopsin (bR) and pharaonis phoborhodopsin (ppR)) and LH complexes (light-harvesting complex II (LH2) and Fenna-Matthews-Olson (FMO) complex). This benchmark study shows the improved description of the color tuning parameters compared to standard DFT functionals. In general, LC-TD-DFTB can exhibit a similar performance as the corresponding LC functionals, allowing a reliable description of excited states properties at significantly reduced cost. The two chromophores investigated here pose complementary challenges: while huge sensitivity to external field perturbation (color tuning) and charge transfer excitations are characteristic for the retinal chromophore, the multi-chromophoric character of the LH complexes emphasizes a correct description of inter-chromophore couplings, giving less importance to color tuning. None of the investigated functionals masters both systems simultaneously with satisfactory accuracy. LC-TD-DFTB, at the current stage, although showing a systematic improvement compared to TD-DFTB cannot be recommended for studying color tuning in retinal proteins, similar to some of the LC-DFT functionals, because the response to external fields is still too weak. For sampling of LH-spectra, however, LC-TD-DFTB is a viable tool, allowing to efficiently sample absorption energies, as shown for three different LH complexes. As the calculations indicate, geometry optimization may overestimate the importance of local minima, which may be averaged over when using trajectories. Fast quantum chemical approaches therefore may allow for a direct sampling of spectra in the near future.

The molecular mechanisms of light adaption in light-harvesting complexes of purple bacteria revealed by a multiscale modeling

The light-harvesting in photosynthetic purple bacteria can be tuned in response to the light conditions during cell growth. One of the used strategies is to change the energy of the excitons in the major fight-harvesting complex, commonly known as LH2. In the present study we report the first systematic investigation of the microscopic origin of the exciton tuning using three complexes, namely the common (high-light) and the low-light forms of LH2 from Rps. acidophila plus a third complex analogous to the PucD complex from Rps. palustris. The study is based on the combination of classical molecular dynamics of each complex in a lipid membrane and excitonic calculations based on a multiscale quantum mechanics/molecular mechanics approach including a polarizable embedding. From the comparative analysis, it comes out that the mechanisms that govern the adaptation of the complex to different light conditions use the different H-bonding environment around the bacteriochlorophyll pigments to dynamically control both internal and inter-pigment degrees of freedom. While the former have a large effect on the site energies, the latter significantly change the electronic couplings, but only the combination of the two effects can fully reproduce the tuning of the final excitons and explain the observed spectroscopic differences.

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.

Cavity-Modified Exciton Dynamics in Photosynthetic Units

Recently, exciton-photon strong coupling has been proposed as a means to control and enhance energy transfer in ensembles of organic molecules. Here, we demonstrate that the exciton dynamics in an archetypal purple bacterial photosynthetic unit, composed of six LH2 antennas surrounding a single LH1 complex, is greatly modified by its interaction with an optical cavity. We develop a Bloch-Redfield master equation approach that accounts for the interplay between the B800 and B850 bacteriochlorophyll molecules within each LH2 antenna, as well as their interactions with the central LH1 complex. Using a realistic parametrization of both the photosynthetic unit and optical cavity, we investigate the formation of polaritons in the system, revealing that these can be tuned to accelerate its exciton dynamics by 3 orders of magnitude. This yields a significant occupation of the LH1 complex, the stage immediately prior to the reaction center, with only a few-femtosecond delay after the initial excitation of the LH2 B800 pigments. Our theoretical findings unveil polaritonic phenomena as a promising route for the characterization, tailoring, and optimization of light-harvesting mechanisms in natural and artificial photosynthetic processes.

Coupling to Charge Transfer States is the Key to Modulate the Optical Bands for Efficient Light Harvesting in Purple Bacteria

The photosynthetic apparatus of purple bacteria uses exciton delocalization and static disorder to modulate the position and broadening of its absorption bands, leading to efficient light harvesting. Its main antenna complex, LH2, contains two rings of identical bacteriochlorophyll pigments, B800 and B850, absorbing at 800 and 850 nm, respectively. It has been an unsolved problem why static disorder of the strongly coupled B850 ring is several times larger than that of the B800 ring. Here we show that mixing between excitons and charge transfer states in the B850 ring is responsible for the effect. The linear absorption spectrum of the LH2 system is simulated by using a multiscale approach with an exciton Hamiltonian generalized to include the charge transfer states that involve adjacent pigment pairs, with static disorder modeled microscopically by molecular dynamics simulations. Our results show that sufficient inhomogeneous broadening of the B850 band, needed for efficient light harvesting, is only obtained by utilizing static disorder in the coupling between local excited and interpigment charge transfer states.

Shaping excitons in light-harvesting proteins through nanoplasmonics

Nanoplasmonics has been used to enhance molecular spectroscopic signals, with exquisite spatial resolution down to the sub-molecular scale. By means of a rigorous, state-of-the-art multiscale model based on a quantum chemical description, here we show that optimally tuned tip-shaped metal nanoparticles can selectively excite localized regions of typically coherent systems, eventually narrowing down to probing one single pigment. The well-known major light-harvesting complex LH2 of purple bacteria has been investigated because of its unique properties, as it presents both high and weak delocalization among subclusters of pigments. This finding opens the way to the direct spectroscopic investigation of quantum-based processes, such as the quantum diffusion of the excitation among the chromophores, and their external manipulation.

The role of charge-transfer states in the spectral tuning of antenna complexes of purple bacteria

The LH2 antenna complexes of purple bacteria occur, depending on light conditions, in various different spectroscopic forms, with a similar structure but different absorption spectra. The differences are related to point changes in the primary amino acid sequence, but the molecular-level relationship between these changes and the resulting spectrum is still not well understood. We undertook a systematic quantum chemical analysis of all the main factors that contribute to the exciton structure, looking at how the environment modulates site energies and couplings in the B800-850 and B800-820 spectroscopic forms of LH2. A multiscale approach combining quantum chemistry and an atomistic classical embedding has been used where mutual polarization effects between the two parts are taken into account. We find that the loss of hydrogen bonds following amino acid changes can only explain a part of the observed blue-shift in the B850 band. The coupling of excitonic states to charge-transfer states, which is different in the two forms, contributes with a similar amount to the overall blue-shift.

Structure and dynamics of propylammonium nitrate-acetonitrile mixtures: An intricate multi-scale system probed with experimental and theoretical techniques

In this article, we report the study of structural and dynamical properties for a series of acetonitrile/propylammonium nitrate mixtures as a function of their composition. These systems display an unusual increase in intensity in their X-ray diffraction patterns in the low-q regime, and their ¹H-NMR diffusion-ordered NMR spectroscopy (DOSY) spectra display unusual diffusivities. However, the magnitude of both phenomena for mixtures of propylammonium nitrate is smaller than those observed for ethylammonium nitrate mixtures with the same cosolvent, suggesting that the cation alkyl tail plays an important role in these observations. The experimental X-ray scattering data are compared with the results of molecular dynamics simulations, including both ab initio studies used to interpret short-range interactions and classical simulations to describe longer range interactions. The higher level calculations highlight the presence of a strong hydrogen bond network within the ionic liquid, only slightly perturbed even at high acetonitrile concentration. These strong interactions lead to the symmetry breaking of the NO₃^- vibrations, with a splitting of about 88 cm^-1 in the ν₃ antisymmetric stretch. The classical force field simulations use a greater number of ion pairs, but are not capable of fully describing the longest range interactions, although they do successfully account for the observed concentration trend, and the analysis of the models confirms the nano-inhomogeneity of these kinds of samples.

Computational Modeling of Exciton-Bath Hamiltonians for Light Harvesting 2 and Light Harvesting 3 Complexes of Purple Photosynthetic Bacteria at Room Temperature

Light harvesting 2 (LH2) complex is the primary component of the photosynthetic unit of purple bacteria that is responsible for harvesting and relaying excitons. The electronic absorption line shape of LH2 contains two major bands at 800 and 850 nm wavelength regions. Under low light conditions, some species of purple bacteria replace LH2 with light harvesting 3 (LH3), a variant form with almost the same structure as the former but with distinctively different spectral features. The major difference between the absorption line shapes of LH2 and LH3 is the shift of the 850 nm band of the former to a new 820 nm region. The microscopic origin of this difference has been the subject of some theoretical/computational investigations. However, the genuine molecular level source of such a difference is not clearly understood yet. This work reports a comprehensive computational study of LH2 and LH3 complexes so as to clarify different molecular level features of LH2 and LH3 complexes and to construct simple exciton-bath models with a common form. All-atomistic molecular dynamics simulations of both LH2 and LH3 complexes provide detailed molecular level structural differences of bacteriochlorophylls (BChls) in the two complexes, in particular, in their patterns of hydrogen bonding (HB) and torsional angles of the acetyl group. Time-dependent density functional theory calculation of the excitation energies of BChls for structures sampled from the MD simulations suggests that the observed differences in the HB and torsional angles cannot fully account for the experimentally observed spectral shift of LH3. Potential sources that can explain the actual spectral shift of LH3 are discussed, and their magnitudes are assessed through fitting of experimental line shapes. These results demonstrate the feasibility of developing simple exciton-bath models for both LH2 and LH3, which can be employed for large-scale exciton quantum dynamics in their aggregates.

Excited state energy fluctuations in the Fenna-Matthews-Olson complex from molecular dynamics simulations with interpolated chromophore potentials

We analyze the environment-induced fluctuation of pigment excitation energies in the Fenna-Matthews-Olson (FMO) complex from various perspectives, by employing an interpolation-based all-atom potential energy model for describing realistic pigment vibrations. We conduct molecular dynamics simulations on a 100 ns timescale, which is an extent that can enclose the effect of static disorder, and demonstrate its timescale separation from fast dynamic disorder. We extract the spectral densities of the complex by considering both the site and the exciton bases. We show that exciton delocalization reduces the effective environmental fluctuation and rationalize this aspect based on a model of fluctuating molecular aggregates. We also obtained the spectral density of the lowest exciton state under low temperature conditions and show that it reasonably well reproduces the experimental result. Finally, by additionally performing non-equilibrium excited state trajectory simulations, we show that the system lies well within the linear response regime after photo-absorption and that the pigments do not visit anharmonic regions of the potential surface to a significant extent. This indicates that methodologies based on harmonic bath models are indeed reasonable approaches for describing the excited state dynamics of the FMO complex.