Accuracy of approximate methods for the calculation of fluorescence-type linear spectra with a complex system-bath coupling

Much can be learned about molecular aggregates by modeling their fluorescence-type spectra. In this study, we systematically describe the accuracy of various methods for simulating fluorescence-type linear spectra in a dimer system with a complex system-environment interaction, which serves as a model for various molecular aggregates, including most photosynthetic light-harvesting complexes (LHCs). We consider the approximate full cumulant expansion (FCE), complex time-dependent Redfield (ctR), time-independent Redfield, and modified Redfield methods and calculate their accuracy as a function of the site energy gap and coupling, excitonic energy gap, and dipole factor (i.e., type of spectrum). We find that the FCE method is the most accurate method for couplings smaller than 300 cm-1 at 300 K, but this method fails for very strong couplings or low temperatures due to inaccurate modeling of the equilibrium initial state. The ctR method performs well for the calculation of fluorescence and linear anisotropy spectra but poorer for circularly polarized fluorescence spectra or for all spectra when the coupling is strong (∼100cm-1). The Redfield and modified Redfield methods generally perform much more poorly than the ctR and FCE methods-especially for small excitonic energy gaps and strong couplings. We show that accurate modeling of the Stokes shift is crucial and present a version of the ctR method that treats both the Stokes shift and initial state correctly for the parameter ranges in plant LHCs. Apart from the application to LHCs, our results will be useful for the spectral characterization and design of organic molecular aggregates.

Inspecting molecular aggregate quadratic vibronic coupling effects using squeezed coherent states

We present a systematic comparison of three quantum mechanical approaches describing excitation dynamics in molecular complexes using the time-dependent variational principle (TDVP) with increasing sophistication trial wavefunctions (ansatze): Davydov D₂, squeezed D₂ (sqD₂) and a numerically exact multiple D₂ (mD₂) ansatz in order to characterize validity of the sqD₂ ansatz. Numerical simulations of molecular aggregate absorption and fluorescence spectra with intra- and intermolecular vibrational modes, including quadratic electronic-vibrational (vibronic) coupling term, which is due to vibrational frequency shift upon pigment excitation are presented. Simulated absorption and fluorescence spectra of a J type molecular dimer with high frequency intramolecular vibrational modes obtained with D₂ and sqD₂ ansatze match the spectra of mD₂ ansatz only in the single pigment model without quadratic vibronic coupling. In general, the use of mD₂ ansatz is required to model an accurate dimer and larger aggregate's spectra. For a J dimer aggregate coupled to a low frequency intermolecular phonon bath, absorption and fluorescence spectra are qualitatively similar using all three ansatze. The quadratic vibronic coupling term in both absorption and fluorescence spectra manifests itself as a lineshape peak amplitude redistribution, static frequency shift and an additional shift, which is temperature dependent. Overall the squeezed D₂ model does not result in a considerable improvement of the simulation results compared to the simplest Davydov D₂ approach.

Extracting the excitonic Hamiltonian of a chlorophyll dimer from broadband two-dimensional electronic spectroscopy

We apply Frenkel exciton theory to model the entire Q-band of a tightly bound chlorophyll dimer inspired by the photosynthetic reaction center of photosystem II. The potential of broadband two-dimensional electronic spectroscopy experiment spanning the Q_x and Q_y regions to extract the parameters of the model dimer Hamiltonian is examined through theoretical simulations of the experiment. We find that the local nature of Q_x excitation enables identification of molecular properties of the delocalized Q_y excitons. Specifically, we demonstrate that the cross-peak region, where excitation energy is resonant with Q_y while detection is at Q_x, contains specific spectral signatures that can reveal the full real-space molecular Hamiltonian, a task that is impossible by considering the Q_y transitions alone. System-bath coupling and site energy disorder in realistic systems may limit the resolution of these spectral signatures due to spectral congestion.

Simulation of Ab Initio Optical Absorption Spectrum of β-Carotene with Fully Resolved S0 and S2 Vibrational Normal Modes

The electronic absorption spectrum of β-carotene (β-Car) is studied using quantum chemistry and quantum dynamics simulations. Vibrational normal modes were computed in optimized geometries of the electronic ground state S₀ and the optically bright excited S₂ state using the time-dependent density functional theory. By expressing the S₂-state normal modes in terms of the ground-state modes, we find that no one-to-one correspondence between the ground- and excited-state vibrational modes exists. Using the ab initio results, we simulated the β-Car absorption spectrum with all 282 vibrational modes in a model solvent at 300 K using the time-dependent Dirac-Frenkel variational principle and are able to qualitatively reproduce the full absorption line shape. By comparing the 282-mode model with the prominent 2-mode model, widely used to interpret carotenoid experiments, we find that the full 282-mode model better describes the high-frequency progression of carotenoid absorption spectra; hence, vibrational modes become highly mixed during the S₀ → S₂ optical excitation. The obtained results suggest that electronic energy dissipation is mediated by numerous vibrational modes.

Modeling Molecular J and H Aggregates using Multiple-Davydov D2 Ansatz

The linear absorption spectrum of J and H molecular aggregates is studied using the time-dependent Dirac-Frenkel variational principle (TDVP) with the multi-Davydov D2 (mD2) trial wavefunction (Ansatz). Both the electronic...

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.

Quantum-Classical Approach for Calculations of Absorption and Fluorescence: Principles and Applications

Absorption and fluorescence spectroscopy techniques provide a wealth of information on molecular systems. The simulations of such experiments remain challenging, however, despite the efforts put into developing the underlying theory. An attractive method of simulating the behavior of molecular systems is provided by the quantum-classical theory─it enables one to keep track of the state of the bath explicitly, which is needed for accurate calculations of fluorescence spectra. Unfortunately, until now there have been relatively few works that apply quantum-classical methods for modeling spectroscopic data. In this work, we seek to provide a framework for the calculations of absorption and fluorescence lineshapes of molecular systems using the methods based on the quantum-classical Liouville equation, namely, the forward-backward trajectory solution (FBTS) and the non-Hamiltonian variant of the Poisson bracket mapping equation (PBME-nH). We perform calculations on a molecular dimer and the photosynthetic Fenna-Matthews-Olson complex. We find that in the case of absorption, the FBTS outperforms PBME-nH, consistently yielding highly accurate results. We next demonstrate that for fluorescence calculations, the method of choice is a hybrid approach, which we call PBME-nH-Jeff, that utilizes the effective coupling theory [Gelzinis, A.; J. Chem. Phys. 2020, 152, 051103] to estimate the excited state equilibrium density operator. Thus, we find that FBTS and PBME-nH-Jeff are excellent candidates for simulating, respectively, absorption and fluorescence spectra of real molecular systems.

Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.

Vibration-mediated energy transport in bacterial reaction center: Simulation study

Exciton energy relaxation in a bacterial Reaction Center (bRC) pigment-protein aggregate presumably involves emission of high energy vibrational quanta to cover wide energy gaps between excitons. Here, we assess this hypothesis utilizing vibronic two-particle theory in modeling of the excitation relaxation process in bRC. Specific high frequency molecular vibrational modes are included explicitly one at a time in order to check which high frequency vibrations are involved in the excitation relaxation process. The low frequency bath modes are treated perturbatively within Redfield relaxation theory. The analysis of the population relaxation rate data indicates energy flow pathways in bRC and suggests that specific vibrations may be responsible for the excitation relaxation process.

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.

Linear and non-linear infrared response of one-dimensional vibrational Holstein polarons in the anti-adiabatic limit: Optical and acoustical phonon models

The theory of linear and non-linear infrared response of vibrational Holstein polarons in one-dimensional lattices is presented in order to identify the spectral signatures of self-trapping phenomena. Using a canonical transformation, the optical response is computed from the small polaron point of view which is valid in the anti-adiabatic limit. Two types of phonon baths are considered: optical phonons and acoustical phonons, and simple expressions are derived for the infrared response. It is shown that for the case of optical phonons, the linear response can directly probe the polaron density of states. The model is used to interpret the experimental spectrum of crystalline acetanilide in the C=O range. For the case of acoustical phonons, it is shown that two bound states can be observed in the two-dimensional infrared spectrum at low temperature. At high temperature, analysis of the time-dependence of the two-dimensional infrared spectrum indicates that bath mediated correlations slow down spectral diffusion. The model is used to interpret the experimental linear-spectroscopy of model α-helix and β-sheet polypeptides. This work shows that the Davydov Hamiltonian cannot explain the observations in the NH stretching range.

Ultrafast energy transfer within the photosystem II core complex

We report 2D electronic spectroscopy on the photosystem II core complex (PSII CC) at 77 K under different polarization conditions. A global analysis of the high time-resolution 2D data shows rapid, sub-100 fs energy transfer within the PSII CC. It also reveals the 2D spectral signatures of slower energy equilibration processes occurring on several to hundreds of picosecond time scales that are consistent with previous work. Using a recent structure-based model of the PSII CC [Y. Shibata, S. Nishi, K. Kawakami, J. R. Shen and T. Renger, J. Am. Chem. Soc., 2013, 135, 6903], we simulate the energy transfer in the PSII CC by calculating auxiliary time-resolved fluorescence spectra. We obtain the observed sub-100 fs evolution, even though the calculated electronic energy shows almost no dynamics at early times. On the other hand, the electronic-vibrational interaction energy increases considerably over the same time period. We conclude that interactions with vibrational degrees of freedom not only induce population transfer between the excitonic states in the PSII CC, but also reshape the energy landscape of the system. We suggest that the experimentally observed ultrafast energy transfer is a signature of excitonic-polaron formation.

Dynamics of coherence, localization and excitation transfer in disordered nanorings

Self-assembled supramolecular aggregates are excellent candidates for the design of efficient excitation transport devices. Both artificially prepared and natural photosynthetic aggregates in plants and bacteria present an important degree of disorder that is supposed to hinder excitation transport. Besides, molecular excitations couple to nuclear motion affecting excitation transport in a variety of ways. We present an exhaustive study of exciton dynamics in disordered nanorings with long-range interactions under the influence of a phonon bath taking the LH2 system of purple bacteria as a model. Nuclear motion is explicitly taken into account by employing the Davydov ansatz description of the polaron and quantum dynamics are obtained using a time-dependent variational method. We reveal an optimal exciton-phonon coupling that suppresses disorder-induced localization and facilitate excitation de-trapping. This excitation transfer enhancement, mediated by environmental phonons, is attributed to energy relaxation toward extended, low-energy excitons provided by the precise LH2 geometry with anti-parallel dipoles and long-range interactions. An analysis of localization and spectral statistics is followed by dynamic measures of coherence and localization, transfer efficiency and superradiance. Linear absorption, 2D photon-echo spectra and diffusion measures of the exciton are examined to monitor the diffusive behavior as a function of the strengths of disorder and exciton-phonon coupling.

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.

Tracing feed-back driven exciton dynamics in molecular aggregates

Excitation, exciton transport, dephasing and energy relaxation, and finally detection processes shift molecular systems into a specific superposition of quantum states causing localization, local heating and finally excitonic polaronic effects.

Simulation of Quantum Dynamics of Excitonic Systems at Finite Temperature: an efficient method based on Thermo Field Dynamics

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on Thermo Field Dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. The solution of Thermo Field Dynamics equations with a novel technique for the propagation of Tensor Trains (Matrix Product States) is implemented and discussed. The methodology is applied to the study of the exciton dynamics in the Fenna-Mathews-Olsen complex using a realistic structured spectral density to model the electron-phonon interaction. The results of the simulations highlight the effect of specific vibrational modes on the exciton dynamics and energy transfer process, as well as call for careful modeling of electron-phonon couplings.

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.

Quantum electron-vibrational dynamics at finite temperature: Thermo field dynamics approach

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on the thermo field dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. A comparison with the theoretically equivalent density matrix formulation shows the key numerical advantages of the present approach. The solution of thermo field dynamics equations with a novel technique for the propagation of tensor trains (matrix product states) is discussed. Numerical applications to model spin-boson systems show that the present approach is a promising tool for the description of quantum dynamics of complex molecular systems at finite temperature.