On Light-Induced Photoconversion of B800 Bacteriochlorophylls in the LH2 Antenna of the Purple Sulfur Bacterium Allochromatium vinosum

Efficiently Transplanting Potential Energy Interpolation Database between Two Systems: Bacteriochlorophyll Case with FMO and LH2 Complexes

Generating a reliable potential energy surface (PES) is an important issue for studying the dynamics of any system with computational simulations. Interpolation mechanics/molecular mechanics (IM/MM) based on a PES interpolation scheme is a useful tool in that regard as it provides an accuracy of a quantum chemistry (QC) level while maintaining its computational cost comparable to conventional MM force fields. Despite this benefit, constructing the database for interpolation itself is still challenging and time-consuming. Here, we present a method with which we can construct the IM database of one system based on a preexisting data set for another related system. We adopt the case of constructing bacteriochlorophyll PESs for the light-harvesting 2 (LH2) complex by utilizing already available IM database for the BChls from the Fenna-Matthews-Olson (FMO) complex. In this method, the IM database from FMO is first transplanted to LH2 by considering BChl displacement vectors that take into account the geometry differences induced by the protein scaffolds. From this transplanted primitive database entries, a relatively small number of effective ones are selected by a survival process based on a genetic algorithm such that the IM energies evaluated at geometries in a conveniently collected prediction set can closely match with the reference QC energies. The selection process is expedited by using two different levels of basis sets for the QC calculations. To demonstrate the utility of the PES thus constructed, we carry out 1 ns of IM/MM dynamics simulations with the finally optimized BChl database for LH2. Indeed, the energy profiles of the snapshots are found to be closely matching with the reference QC calculation data, with only ∼0.07 eV of errors in the ground- and excited-state energies and ∼0.008 eV of errors in the transition energies. We also show that properly selecting data points is actually quite important for generating an IM PES toward performing molecular dynamics simulations.

Dichotomous Disorder versus Excitonic Splitting of the B800 Band of Allochromatium vinosum

The LH2 antenna complex of the purple bacterium Allochromatium vinosum has a distinct double peak structure of the 800 nm band (B800). Several hypotheses were proposed to explain its origin. Recent 77 K two-dimensional electronic spectroscopy data suggested that excitonic coupling of dimerized bacteriochlorophylls (BChls) within the B800 ring is largely responsible for the B800 split [M. Schröter et al., J. Phys. Chem. Lett. 2018, 9, 1340]. Here we argue that the excitonic interactions between BChls in the B800 ring, though present, are weak and cannot explain the B800 band split. This conclusion is based on hole-burning data and modeling studies using an exciton model with dichotomous protein conformation disorder. Therefore, we uphold our earlier interpretation, first reported by Kell et al. [ J. Phys. Chem. B 2017, 121, 9999], that the two B800 sub-bands are due to different site-energies (most likely due to weakly and strongly hydrogen-bonded B800 BChls).

Two-dimensional electronic spectroscopy of anharmonic molecular potentials

Two-dimensional electronic spectroscopy (2DES) is a powerful tool in the study of coupled electron-phonon dynamics, yet very little is known about how nonlinearities in the electron-phonon coupling, arising from anharmonicities in the nuclear potentials, affect the spectra. These become especially relevant when the coupling is strong. From the linear spectroscopies, anharmonicities are known to give structure to the zero-phonon line and to break mirror-symmetry between absorption and emission, but the 2D analogues of these effects have not been identified. Using a simple two-level model where the electronic states are described by (displaced) harmonic oscillators with differing curvatures or displaced Morse oscillators, we find that the zero-phonon line shape is essentially transferred to the diagonal in 2DES spectra, and that anharmonicities break a horizontal mirror-symmetry in the infinite waiting time limit. We also identify anharmonic effects that are only present in 2DES spectra: twisting of cross-peaks stemming from stimulated emission signals; and oscillation period mismatch between ground state bleach and stimulated emission (for harmonic oscillators with differing curvatures), or inherently chaotic oscillations (for Morse oscillators). Our findings will facilitate an improved understanding of 2DES spectra and aid the interpretation of signals that are more realistic than those arising from simple models.