Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer

Intermolecular Protein-Water Coupling Impedes the Coupling Between the Amide A and Amide I Mode in Interfacial Proteins

The coupling between different vibrational modes in proteins is essential for chemical dynamics and biological functions and is linked to the propagation of conformational changes and pathways of allosteric communication. However, little is known about the influence of intermolecular protein-H2O coupling on the vibrational coupling between amide A (NH) and amide I (C═O) bands. Here, we investigate the NH/CO coupling strength in various peptides with different secondary structures at the lipid cell membrane/H2O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy (SFG-VS) in which a femtosecond infrared pump is used to excite the amide A band, and SFG-VS is used to probe transient spectral evolution in the amide A and amide I bands. Our results reveal that the NH/CO coupling strength strongly depends on the bandwidth of the amide I mode and the coupling of proteins with water molecules. A large extent of protein-water coupling significantly reduces the delocalization of the amide I mode along the peptide chain and impedes the NH/CO coupling strength. A large NH/CO coupling strength is found to show a strong correlation with the high energy transfer rate found in the light-harvesting proteins of green sulfur bacteria, which may understand the mechanism of energy transfer through a molecular system and assist in controlling vibrational energy transfer by engineering the molecular structures to achieve high energy transfer efficiency.

First-principles simulation of excitation energy transfer and transient absorption spectroscopy in the CP29 light-harvesting complex

The dynamics of delocalized excitons in light-harvesting complexes (LHCs) can be investigated using different experimental techniques, and transient absorption (TA) spectroscopy is one of the most valuable methods for this purpose. A careful interpretation of TA spectra is essential for the clarification of excitation energy transfer (EET) processes occurring during light-harvesting. However, even in the simplest LHCs, a physical model is needed to interpret transient spectra as the number of EET processes occurring at the same time is very large to be disentangled from measurements alone. Physical EET models are commonly built by fittings of the microscopic exciton Hamiltonians and exciton-vibrational parameters, an approach that can lead to biases. Here, we present a first-principles strategy to simulate EET and transient absorption spectra in LHCs, combining molecular dynamics and accurate multiscale quantum chemical calculations to obtain an independent estimate of the excitonic structure of the complex. The microscopic parameters thus obtained are then used in EET simulations to obtain the population dynamics and the related spectroscopic signature. We apply this approach to the CP29 minor antenna complex of plants for which we follow the EET dynamics and transient spectra after excitation in the chlorophyll b region. Our calculations reproduce all the main features observed in the transient absorption spectra and provide independent insight on the excited-state dynamics of CP29. The approach presented here lays the groundwork for the accurate simulation of EET and unbiased interpretation of transient spectra in multichromophoric systems.

Impact of Spatial Inhomogeneity on Excitation Energy Transport in the Fenna-Matthews-Olson Complex

The dynamics of the excitation energy transfer (EET) in photosynthetic complexes is an interesting question both from the perspective of fundamental understanding and the research in artificial photosynthesis. Over the past decade, very accurate spectral densities have been developed to capture spatial inhomogeneities in the Fenna-Matthews-Olson (FMO) complex. However, challenges persist in numerically simulating these systems, both in terms of parameterizing them and following their dynamics over long periods of time because of long non-Markovian memories. We investigate the dynamics of FMO with the exact treatment of various theoretical spectral densities using the new tensor network path integral-based methods, which are uniquely capable of addressing the difficulty of long memory length and incoherent Förster theory. It is also important to be able to analyze the pathway of EET flow, which can be difficult to identify given the non-trivial structure of connections between bacteriochlorophyll molecules in FMO. We use the recently introduced ideas of relating coherence to population derivatives to analyze the transport process and reveal some new routes of transport. The combination of exact and approximate methods sheds light on the role of coherences in affecting the fine details of the transport and promises to be a powerful toolbox for future exploration of other open systems with quantum transport.

Time-resolved circular dichroism of excitonic systems: theory and experiment on an exemplary squaraine polymer

Experimental and theoretical foundations for femtosecond time-resolved circular dichroism (TRCD) spectroscopy of excitonic systems are presented. In this method, the system is pumped with linearly polarized light and the signal is defined as the difference between the transient absorption spectrum probed with left and with right circularly polarized light. We present a new experimental setup with a polarization grating as key element to generate circularly polarized pulses. Herein the positive (negative) first order of the diffracted light is left-(right-)circularly polarized and serves as a probe pulse in a TRCD experiment. The grating is capable of transferring ultrashort broadband pulses ranging from 470 nm to 720 nm into two separate beams with opposite ellipticity. By applying a specific chopping scheme we can switch between left and right circular polarizations and detect transient absorption (TA) and TRCD spectra on a shot-to-shot basis simultaneously. We perform experiments on a squaraine polymer, investigating excitonic dynamics, and we develop a general theory for TRCD experiments of excitonically coupled systems that we then apply to describe the experimental data in this particular example. At a magic angle of 54.7° between the pump-pulse polarization and the propagation direction of the probe pulse, the TRCD and TA signals become particularly simple to analyze, since the orientational average over random orientations of complexes factorizes into that of the interaction with the pump and the probe pulse, and the intrinsic electric quadrupole contributions to the TRCD signal average to zero for isotropic samples. Application of exciton theory to linear absorption and to linear circular dichroism spectra of squaraine polymers reveals the presence of two fractions of polymer conformations, a dominant helical conformation with close interpigment distances that are suggested to lead to short-range contributions to site energy shifts and excitonic couplings of the squaraine molecules, and a fraction of unfolded random coils. Theory demonstrates that TRCD spectra of selectively excited helices can resolve state populations that are practically invisible in TA spectroscopy due to the small dipole strength of these states. A qualitative interpretation of TRCD and TA spectra in the spectral window investigated experimentally is offered. The 1 ps time component found in these spectra is related to the slow part of exciton relaxation obtained between states of the helix in the low-energy half of the exciton manifold. The dominant 140 ps time constant reflects the decay of excited states to the electronic ground state.

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.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0–0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang–Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang–Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0–0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

The importance of reaction field and screening effects on the excitonic couplings is demonstrated, and from quantum-chemical calculations a single scaling factor is derived that can be used to improve simple models based on the Poisson equation.

Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes

Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.

Semiclassical Modified Redfield and Generalized Förster Theories of Exciton Relaxation/Transfer in Light-Harvesting Complexes: The Quest for the Principle of Detailed Balance

A conceptual problem of transfer theories that use a semiclassical description of the electron-vibrational coupling is the neglect of the correlation between momenta and coordinates of nuclei. In the Redfield theory of exciton relaxation, this neglect leads to a violation of the principle of detailed balance; equal “uphill” and “downhill” transfer rate constants are obtained. Here, we investigate how this result depends on nuclear reorganization effects, neglected in Redfield but taken into account in the modified Redfield theory. These reorganization effects, resulting from a partial localization of excited states, are found to promote a preferential “downhill” relaxation of excitation energy. However, for realistic spectral densities of light-harvesting antennae in photosynthesis, the reorganization effects are too small to compensate for the missing coordinate–momentum uncertainty. For weaker excitonic couplings as they occur between domains of strongly coupled pigments, we find the principle of detailed balance to be fulfilled in a semiclassical variant of the generalized Förster theory. A qualitatively correct description of the transfer is obtained with this theory at a significantly lower computational cost as with the quantum generalized Förster theory. Larger deviations between the two theories are expected for large energy gaps as they occur in complexes with chemically different pigments.