Study of Electronic Structures and Pigment–Protein Interactions in the Reaction Center of Thermochromatium tepidum with a Dynamic Environment

Gaussian process regression for absorption spectra analysis of molecular dimers

A common task is the determination of system parameters from spectroscopy, where one compares the experimental spectrum with calculated spectra, that depend on the desired parameters. Here we discuss an approach based on a machine learning technique, where the parameters for the numerical calculations are chosen from Gaussian Process Regression (GPR). This approach does not only quickly converge to an optimal parameter set, but in addition provides information about the complete parameter space, which allows for example to identify extended parameter regions where numerical spectra are consistent with the experimental one. We consider as example dimers of organic molecules and aim at extracting in particular the interaction between the monomers, and their mutual orientation. We find that indeed the GPR gives reliable results which are in agreement with direct calculations of these parameters using quantum chemical methods.

Fully Quantum Modeling of Exciton Diffusion in Mesoscale Light Harvesting Systems

It has long been a challenge to accurately and efficiently simulate exciton-phonon dynamics in mesoscale photosynthetic systems with a fully quantum mechanical treatment due to extensive computational resources required. In this work, we tackle this seemingly intractable problem by combining the Dirac-Frenkel time-dependent variational method with Davydov trial states and implementing the algorithm in graphic processing units. The phonons are treated on the same footing as the exciton. Tested with toy models, which are nanoarrays of the B850 pigments from the light harvesting 2 complexes of purple bacteria, the methodology is adopted to describe exciton diffusion in huge systems containing more than 1600 molecules. The superradiance enhancement factor extracted from the simulations indicates an exciton delocalization over two to three pigments, in agreement with measurements of fluorescence quantum yield and lifetime in B850 systems. With fractal analysis of the exciton dynamics, it is found that exciton transfer in B850 nanoarrays exhibits a superdiffusion component for about 500 fs. Treating the B850 ring as an aggregate and modeling the inter-ring exciton transfer as incoherent hopping, we also apply the method of classical master equations to estimate exciton diffusion properties in one-dimensional (1D) and two-dimensional (2D) B850 nanoarrays using derived analytical expressions of time-dependent excitation probabilities. For both coherent and incoherent propagation, faster energy transfer is uncovered in 2D nanoarrays than 1D chains, owing to availability of more numerous propagating channels in the 2D arrangement.

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.

Coupled excitation energy and charge transfer dynamics in reaction centre inspired model systems

Functional operation conditions of reaction centre core complexes require the tight coupling of exciton states to concomitant charge separation. Rigorous theoretical treatment of such integrated excitation energy transfer (EET) and charge transfer (CT) dynamics is particularly challenging due to (i) appreciable system sizes, (ii) inter-site and system-bath couplings of similar magnitude that render the Born-Markov approximation invalid, (iii) substantial reorganization energies of CT states, and (iv) the presence of complex structured spectral densities due to vibrational modes of the surroundings. We present numerical simulations on bacterial reaction centre (bRC) inspired model systems that utilize the recently developed MACGIC-iQUAPI method [Richter et al., J. Chem. Phys., 2017, 146, 214101]. The simulations demonstrate that the method provides a rigorous framework for the investigation of such integrated EET-CT dynamics. First, the applicability of the MACGIC-iQUAPI method is explored for a transition from monotonically decaying to oscillatory system-bath influence coefficients, a behavior inherently imposed by structured bath spectral densities. Tightly coupled EET and CT dynamics is further addressed for an excitonic subsystem that resembles strong coupling of special pair states and serves as donor towards a generic bridge-acceptor system. By solving the dissipative quantum dynamics of such bRC inspired model systems, the quenching of excitonic coherence on the hundreds of femtoseconds timescale is explored via a variation of the bridge state energetics, resembling a continuous transition from sequential to superexchange mediated CT regimes. Further, the simulations explore the influence of resonant vibrational modes on the quenching of excitonic coherence via CT. The results reveal a moderate influence of vibrational mode on charge separation dynamics in regimes of biologically relevant EET and CT dynamics.

Spectrally Selective Spectroscopy of Native Ca-Containing and Ba-Substituted LH1-RC Core Complexes from Thermochromatium tepidum

The LH1-RC core complex from the thermophilic photosynthetic purple sulfur bacterium Thermochromatium tepidum has recently attracted interest of many researchers because of its several unique properties, such as increased robustness against environmental hardships and the much red-shifted near-infrared absorption spectrum of the LH1 antenna exciton polarons. The known near-atomic-resolution crystal structure of the complex well supported this attention. Yet several mechanistic aspects of the complex prominence remained to be understood. In this work, samples of the native, Ca²⁺-containing core complexes were investigated along with those destabilized by Ba²⁺ substitution, using various spectrally selective steady-state and picosecond time-resolved spectroscopic techniques at physiological and cryogenic temperatures. As a result, the current interpretation of exciton spectra of the complex was significantly clarified. Specifically, by evaluating the homogeneous and inhomogeneous compositions of the spectra, we showed that there is little to no effect of cation substitution on the dynamic or kinetic properties of antenna excitons. Reasons of the extra red shift of absorption/fluorescence spectra observed in the Ca-LH1-RC and not in the Ba-LH1-RC complex should thus be searched in subtle structural differences following the inclusion of different cations into the core complex scaffold.

Exciton-Coupled Electron Transfer Process Controlled by Non-Markovian Environments

We theoretically investigate an exciton-coupled electron transfer (XCET) process that is conversion of an exciton into a charge transfer state. This conversion happens in an exciton transfer (XT) process, and the electron moves away in an electron transfer (ET) process in multiple environments (baths). This XCET process plays an essential role in the harvesting of solar energy in biological and photovoltaic materials. We develop a practical theoretical model to study the efficiency of the XCET process that occurs either in consecutive or concerted processes under the influence of non-Markovian baths. The role of quantum coherence in the XT-ET system and the baths is investigated using reduced hierarchal equations of motion (HEOM). This model includes independent baths for each XT and ET state, in addition to a XCET bath for the conversion process. We found that, while quantum system-bath coherence is important in the XT and ET processes, coherence between the XT and ET processes must be suppressed in order to realize that an efficient irreversible XCET process through the weak off-diagonal interaction between the XT and ET bridge sites arises from an XCET bath.

Relation between Vibrational Dephasing Time and Energy Gap Fluctuations

Dephasing processes are present in basically all applications in which quantum mechanics plays a role. These applications certainly include excitation energy and charge transfer in biological systems. In a previous study, we have analyzed the vibrational dephasing time as a function of energy gap fluctuation for a large set of molecular simulations. In that investigation, individual molecular subunits were the focus of the calculations. The set of studied molecules included bacteriochlorophylls in Fenna-Matthews-Olson and light-harvesting system 2 complexes as well as bilins in PE545 aggregates. The present work extends this study to entire complexes, including the respective intermolecular couplings. Again, it can be concluded that a universal and inverse proportionality exists between dephasing time and variance of the excitonic energy gap fluctuations, whereas the respective proportionality constants can be rationalized using the energy gap autocorrelation functions. Furthermore, these findings can be extended to the gaps between higher-lying neighboring excitonic states.

Intermolecular Modes between LH2 Bacteriochlorophylls and Protein Residues: The Effect on the Excitation Energies

Light-harvesting system 2 (LH2) executes the primary processes of photosynthesis in purple bacteria; photon absorption, and energy transportation to the reaction center. A detailed mechanistic insight into these operations is obscured by the complexity of the light-harvesting systems, particularly by the chromophore-environment interaction. In this work, we focus on the effects of the protein residues that are ligated to the bacteriochlorophylls (BChls) and construct potential energy surfaces of the ground and first optically excited state for the various BChl-residue systems where we in each case consider two degrees of freedom in the intermolecular region. We find that the excitation energies are only slightly affected by the considered modes. In addition, we see that axial ligands and hydrogen-bonded residues have opposite effects on both excitation energies and oscillator strengths by comparing to the isolated BChls. Our results indicate that only a small part of the chromophore-environment interaction can be associated with the intermolecular region between a BChl and an adjacent residue, but that it may be possible to selectively raise or lower the excitation energy at the axial and planar residue positions, respectively.