Taming the third order cumulant approximation to linear optical spectroscopy

The second order cumulant method offers a promising pathway to predicting optical properties in condensed phase systems. It allows for the computation of linear absorption spectra from excitation energy fluctuations sampled along molecular dynamics (MD) trajectories, fully accounting for vibronic effects, direct solute-solvent interactions, and environmental polarization effects. However, the second order cumulant approximation only guarantees accurate line shapes for energy gap fluctuations obeying Gaussian statistics. A third order correction has recently been derived but often yields unphysical spectra or divergent line shapes for moderately non-Gaussian fluctuations due to the neglect of higher order terms in the cumulant expansion. In this work, we develop a corrected cumulant approach, where the collective effect of neglected higher order contributions is approximately accounted for through a dampening factor applied to the third order cumulant term. We show that this dampening factor can be expressed as a function of the skewness and kurtosis of energy gap fluctuations and can be parameterized from a large set of randomly sampled model Hamiltonians for which exact spectral line shapes are known. This approach is shown to systematically remove unphysical contributions in the form of negative absorbances from cumulant spectra in both model Hamiltonians and condensed phase systems sampled from MD and dramatically improves over the second order cumulant method in describing systems exhibiting Duschinsky mode mixing effects. We successfully apply the approach to the coumarin-153 dye in toluene, obtaining excellent agreement with experiment.

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.

Electronic-Vibrational Resonance Does Not Significantly Alter Steady-State Transport in Natural Light-Harvesting Systems

Oscillations in time-dependent two-dimensional electronic spectra appear as evidence of quantum coherence in light-harvesting systems related to electronic-vibrational resonant interactions. Nature, however, takes place in a non-equilibrium steady-state; therefore, the relevance of these arguments to the natural process is unclear. Here, we examine the role of intramolecular vibrations in the non-equilibrium steady-state of photosynthetic dimers in the natural scenario of incoherent light excitation. Specifically, we analyze the PEB dimer in the cryptophyte algae PE545 antenna protein. It is found that vibrations resonant with the energy difference between exciton states only marginally increase the quantum yield and the imaginary part of the intersite coherence that is relevant for transport compared to non-resonant vibrations in the natural non-equilibrium steady-state. That is, the electronic-vibrational resonance interaction does not significantly enhance energy transport under natural incoherent light excitation conditions.

Low-intensity pulsed ultrasound-generated singlet oxygen induces telomere damage leading to glioma stem cell awakening from quiescence

Cancer stem cells, quiescent and drug resistant, have become a therapeutic target. Unlike high-intensity focused ultrasound directly killing tumor, low-intensity pulsed ultrasound (LIPUS), a new noninvasive physical device, promotes pluripotent stem cell differentiation and is primarily applied in tissue engineering but rarely in oncotherapy. We explored the effect and mechanism of LIPUS on glioma stem cell (GSC) expulsion from quiescence. Here, we observed that LIPUS led to attenuated expression of GSC biomarkers, promoted GSC escape from G0 quiescence, and significantly weakened the Wnt and Hh pathways. Of note, LIPUS transferred sonomechanical energy into cytochrome c and B5 proteins, which converted oxygen molecules into singlet oxygen, triggering telomere crisis. The in vivo and in vitro results confirmed that LIPUS enhanced the GSC sensitivity to temozolomide. These results demonstrated that LIPUS "waked up" GSCs to improve their sensitivity to chemotherapy, and importantly, we confirmed the direct targeted proteins of LIPUS in GSCs.

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.

Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions

First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.

Transient synchronisation and quantum coherence in a bio-inspired vibronic dimer

Synchronisation is a collective phenomenon widely investigated in classical oscillators and, more recently, in quantum systems. However it remains unclear what features distinguish synchronous behaviour in these two scenarios. Recent works have shown that investigating synchronisation dynamics in open quantum systems can give insight into this issue. Here we study transient synchronisation in a bio-inspired vibronic dimer, where electronic excitation dynamics is mediated by coherent interactions with intramolecular vibrational modes. We show that the synchronisation dynamics of local mode displacements exhibit a rich behaviour which arises directly from the distinct time-evolutions of different vibronic quantum coherences. Furthermore, our study shows that coherent energy transport in this bio-inspired system is concomitant with the emergence of positive synchronisation between mode displacements. Our work provides further understanding of the relations between quantum coherence and synchronisation in open quantum systems and suggests an interesting role for coherence in biomolecules, that of promoting synchronisation of vibrational motions driven out of thermal equilibrium.

Quantum dynamics of vibration-assisted excitation energy transfer in phycobiliprotein light-harvesting complex

Phycobiliprotein is a light-harvesting complex containing linear tetrapyrrole bilin pigments that are responsible for absorption and funneling the sun's energy in cryptophytes algae. In particular, the protein structure determines relative positions and orientations of the pigments and thus controls energy transfer pathways. The present research reveals the impact of molecular vibrations (in the 850-2700 cm^-1 region) on excitation energy transfer in phycobiliprotein. The analysis of the excitation energy transfer pathways indicates a possibility of the coherent mechanism of energy transfer (delocalization) in central dihydrobiliverdin pigments and incoherent vibration-assisted energy transfer to peripheral phycocyanobilin pigments at a sub-picosecond time scale. A computational approach that enables modeling the dynamics of the excitation energy transfer with the quantum master equation formalism employing Huang-Rhys factors to describe electronic-vibrational coupling has been developed. The computational methodology has been implemented in PyFREC software.

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.

Determination of the protonation preferences of bilin pigments in cryptophyte antenna complexes

The light-harvesting mechanisms of cryptophyte antenna complexes have attracted considerable attention due to their ability to exhibit maximal photosynthetic activity under very low-light conditions and to display several colors, as well as the observation of vibronic coherent features in their two-dimensional electronic spectra. However, detailed investigations on the interplay between the protein environment and their light-harvesting properties are hampered by the uncertainty related to the protonation state of the underlying bilin pigments. Here we study the protonation preferences of four types of bilin pigments including 15,16-dihydrobiliverdin (DBV), phycoerythrobilin (PEB), phycocyanobilin (PCB) and mesobiliverdin (MBV), which are found in phycoerythrin PE545 and phycocyanin PC577, PC612, PC630 and PC645 complexes. We apply quantum chemical calculations coupled to continuum solvation calculations to predict the intrinsic acidity of bilins in aqueous solution, and then combine molecular dynamics simulations with empirical pKa estimates to investigate the impact of the local protein environment on the acidity of the pigments. We also report measurements of the absorption spectra of the five complexes in a wide range of pH in order to validate our simulations and investigate possible changes in the light harvesting properties of the complexes in the range of physiological pH found in the lumen (pH ∼ 5-7). The results suggest a pKa > 7 for DBV and MBV pigments in the α polypeptide chains of PE545 and PC630/PC645 complexes, which are not coordinated to a negatively charged amino acid. For the other PEB, DBV and PCB pigments, which interact with a Glu or Asp side chain, higher pKa values (pKa > 8) are estimated. Overall, the results support a preferential population of the fully protonated state for bilins in cryptophyte complexes under physiological conditions regardless of the specific type of pigment and local protein environment.

Vibronic Wavepackets and Energy Transfer in Cryptophyte Light-Harvesting Complexes

Determining the key features of high-efficiency photosynthetic energy transfer remains an ongoing task. Recently, there has been evidence for the role of vibronic coherence in linking donor and acceptor states to redistribute oscillator strength for enhanced energy transfer. To gain further insights into the interplay between vibronic wavepackets and energy-transfer dynamics, we systematically compare four structurally related phycobiliproteins from cryptophyte algae by broad-band pump-probe spectroscopy and extend a parametric model based on global analysis to include vibrational wavepacket characterization. The four phycobiliproteins isolated from cryptophyte algae are two "open" structures and two "closed" structures. The closed structures exhibit strong exciton coupling in the central dimer. The dominant energy-transfer pathway occurs on the subpicosecond timescale across the largest energy gap in each of the proteins, from central to peripheral chromophores. All proteins exhibit a strong 1585 cm^-1 coherent oscillation whose relative amplitude, a measure of vibronic intensity borrowing from resonance between donor and acceptor states, scales with both energy-transfer rates and damping rates. Central exciton splitting may aid in bringing the vibronically linked donor and acceptor states into better resonance resulting in the observed doubled rate in the closed structures. Several excited-state vibrational wavepackets persist on timescales relevant to energy transfer, highlighting the importance of further investigation of the interplay between electronic coupling and nuclear degrees of freedom in studies on high-efficiency photosynthesis.

Mechanistic Regimes of Vibronic Transport in a Heterodimer and the Design Principle of Incoherent Vibronic Transport in Phycobiliproteins

Following the observation of coherent oscillations in nonlinear spectra of photosynthetic pigment protein complexes, in particular, phycobilliproteins such as PC645, coherent vibronic transport has been suggested as a design principle for novel light-harvesting materials. Vibronic transport between energetically remote pigments is coherent when the presence of a vibration resonant with the electronic energy gap supports transient delocalization between the electronic excited states. We establish the mechanism of vibronic transport for a model heterodimer across a wide range of molecular parameter values. The resulting mechanistic map demonstrates that the molecular parameters of phycobiliproteins in fact support incoherent vibronic transport. This result points to an important design principle: Incoherent vibronic transport is more efficient than a coherent mechanism when energetic disorder exceeds the coupling between the donor and vibrationally excited acceptor states. Finally, our results suggest that the role of coherent vibronic transport in pigment protein complexes should be reevaluated.

Local protein solvation drives direct down-conversion in phycobiliprotein PC645 via incoherent vibronic transport

The mechanisms controlling excitation energy transport (EET) in light-harvesting complexes remain controversial. Following the observation of long-lived beats in 2D electronic spectroscopy of PC645, vibronic coherence, the delocalization of excited states between pigments supported by a resonant vibration, has been proposed to enable direct excitation transport from the highest-energy to the lowest-energy pigments, bypassing a collection of intermediate states. Here, we instead show that for phycobiliprotein PC645 an incoherent vibronic transport mechanism is at play. We quantify the solvation dynamics of individual pigments using ab initio quantum mechanics/molecular mechanics (QM/MM) nuclear dynamics. Our atomistic spectral densities reproduce experimental observations ranging from absorption and fluorescence spectra to the timescales and selectivity of down-conversion observed in transient absorption measurements. We construct a general model for vibronic dimers and establish the parameter regimes of coherent and incoherent vibronic transport. We demonstrate that direct down-conversion in PC645 proceeds incoherently, enhanced by large reorganization energies and a broad collection of high-frequency vibrations. We suggest that a similar incoherent mechanism is appropriate across phycobiliproteins and represents a potential design principle for nanoscale control of EET.

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.

Exciton transport in the PE545 complex: insight from atomistic QM/MM-based quantum master equations and elastic network models

In this paper, we work out a parameterization of environmental noise within the Haken-Strobl-Reinenker (HSR) model for the PE545 light-harvesting complex, based on atomic-level quantum mechanics/molecular mechanics (QM/MM) simulations. We use this approach to investigate the role of various auto- and cross-correlations in the HSR noise tensor, confirming that site-energy autocorrelations (pure dephasing) terms dominate the noise-induced exciton mobility enhancement, followed by site energy-coupling cross-correlations for specific triplets of pigments. Interestingly, several cross-correlations of the latter kind, together with coupling-coupling cross-correlations, display clear low-frequency signatures in their spectral densities in the 30-70 [Formula: see text] region. These slow components lie at the limits of validity of the HSR approach, which requires that environmental fluctuations be faster than typical exciton transfer time scales. We show that a simple coarse-grained elastic-network-model (ENM) analysis of the PE545 protein naturally spotlights collective normal modes in this frequency range that represent specific concerted motions of the subnetwork of cysteines covalenty linked to the pigments. This analysis strongly suggests that protein scaffolds in light-harvesting complexes are able to express specific collective, low-frequency normal modes providing a fold-rooted blueprint of exciton transport pathways. We speculate that ENM-based mixed quantum classical methods, such as Ehrenfest dynamics, might be promising tools to disentangle the fundamental designing principles of these dynamical processes in natural and artificial light-harvesting structures.

Estimation of damped oscillation associated spectra from ultrafast transient absorption spectra

When exciting a complex molecular system with a short optical pulse, all chromophores present in the system can be excited. The resulting superposition of electronically and vibrationally excited states evolves in time, which is monitored with transient absorption spectroscopy. We present a methodology to resolve simultaneously the contributions of the different electronically and vibrationally excited states from the complete data. The evolution of the excited states is described with a superposition of damped oscillations. The amplitude of a damped oscillation cos(ω_nt)exp(-γ_nt) as a function of the detection wavelength constitutes a damped oscillation associated spectrum DOAS_n(λ) with an accompanying phase characteristic φ_n(λ). In a case study, the cryptophyte photosynthetic antenna complex PC612 which contains eight bilin chromophores was excited by a broadband optical pulse. Difference absorption spectra from 525 to 715 nm were measured until 1 ns. The population dynamics is described by four lifetimes, with interchromophore equilibration in 0.8 and 7.5 ps. We have resolved 24 DOAS with frequencies between 130 and 1649 cm^-1 and with damping rates between 0.9 and 12 ps^-1. In addition, 11 more DOAS with faster damping rates were necessary to describe the "coherent artefact." The DOAS contains both ground and excited state features. Their interpretation is aided by DOAS analysis of simulated transient absorption signals resulting from stimulated emission and ground state bleach.

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.

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.

Chebyshev Expansion Applied to Dissipative Quantum Systems

To determine the dynamics of a molecular aggregate under the influence of a strongly time-dependent perturbation within a dissipative environment is still, in general, a challenge. The time-dependent perturbation might be, for example, due to external fields or explicitly treated fluctuations within the environment. Methods to calculate the dynamics in these cases do exist though some of these approaches assume that the corresponding correlation functions can be written as a weighted sum of exponentials. One such theory is the hierarchical equations of motion approach. If the environment, however, is described by a complex spectral density or if its temperature is low, these approaches become very inefficient. Therefore, we propose a scheme based on a Chebyshev decomposition of the bath correlation functions and detail the respective quantum master equations within second-order perturbation theory in the environmental coupling. Similar approaches have recently been proposed for systems coupled to Fermionic reservoirs. The proposed scheme is tested for a simple two-level system and compared to existing results. Furthermore, the advantages and disadvantages of the present Chebyshev approach are discussed.

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.

Vibronic resonances facilitate excited-state coherence in light-harvesting proteins at room temperature

Until recently it was believed that photosynthesis, a fundamental process for life on earth, could be fully understood with semiclassical models. However, puzzling quantum phenomena have been observed in several photosynthetic pigment-protein complexes, prompting questions regarding the nature and role of these effects. Recent attention has focused on discrete vibrational modes that are resonant or quasi-resonant with excitonic energy splittings and strongly coupled to these excitonic states. Here we unambiguously identify excited state coherent superpositions in photosynthetic light-harvesting complexes using a new experimental approach. Decoherence on the time scale of the excited state lifetime allows low energy (56 cm(-1)) oscillations on the signal intensity to be observed. In conjunction with an appropriate model, these oscillations provide clear and direct experimental evidence that the persistent coherences observed originate from quantum superpositions among vibronic excited states.

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.