Photoinduced Intra- and Intermolecular Energy Transfer in Chlorophyll a Dimer

Coupled Nuclear and Electron Dynamics in Chlorophyll Unraveled by XMS-CASPT2 X-ray Absorption Spectra

Attosecond spectroscopy, especially time-resolved X-ray absorption spectra (XAS), enables direct observation of ultrafast molecular dynamics. The complementary and even preceding development of theoretical simulations can offer the necessary guidance and stimulate new experiments. In this work, we simulated high-level XAS for the magnesium and nitrogen K-edge of chlorophyll a. In our previous work on the ultrafast relaxation process in the Q-band, our quantum dynamics simulations found the Qx and Qy states to be energetically close and therefore strongly coupled. Here, we analyze the strong coupling between Qx and Qy via XAS, indicating promising possibilities for experimental observation. The excited-state energies, potential energy surfaces, and XAS are computed at the XMS-CASPT2 level of theory to capture the complex multireference character of chlorophyll excitations. In our simulated spectra, we could follow the ultrafast population transfer between Qx and Qy and thus draw conclusions about the strong vibrational coupling between them.

Complementing Adiabatic and Nonadiabatic Methods To Understand Internal Conversion Dynamics in Porphyrin Derivatives

We analyze the internal conversion dynamics within the Qy and Qx excited states of both bare and functionalized porphyrins, which are known to exhibit significantly different time constants experimentally. Through the integration of two complementary approaches, static calculation of per-mode reorganization energies and nonadiabatic molecular dynamics, we achieve a comprehensive understanding of the factors determining the different behavior of the two molecules. We identify the key normal and essential modes responsible for the population transfer between excited states and discuss the efficacy of different statistical and nonstatistical analyses in providing a full physics-based description of the phenomenon.

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair PD1-PD2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state PD1--PD2+ stabilizes energetically below the excitonic state. This underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

Excitation Energy Transfer between Higher Excited States of Photosynthetic Pigments: 2. Chlorophyll b is a B Band Excitation Trap

Chlorophylls (Chls) are known for fast, subpicosecond internal conversion (IC) from ultraviolet/blue absorbing ("B" or "Soret" states) to the energetically lower, red light-absorbing Q states. Consequently, excitation energy transfer (EET) in photosynthetic pigment-protein complexes involving the B states has so far not been considered. We present, for the first time, a theoretical framework for the existence of B-B EET in tightly coupled Chl aggregates such as photosynthetic pigment-protein complexes. We show that according to a Förster resonance energy transport (FRET) scheme, unmodulated B-B EET has an unexpectedly high range. Unsuppressed, it could pose an existential threat-the damage potential of blue light for photochemical reaction centers (RCs) is well-known. This insight reveals so-far undescribed roles for carotenoids (Crts, cf. previous article in this series) and Chl b (this article) of possibly vital importance. Our model system is the photosynthetic antenna pigment-protein complex (CP29). The focus of the study is on the role of Chl b for EET in the Q and B bands. Further, the initial excited pigment distribution in the B band is computed for relevant solar irradiation and wavelength-centered laser pulses. It is found that both accessory pigment classes compete efficiently with Chl a absorption in the B band, leaving only 40% of B band excitations for Chl a. B state population is preferentially relocated to Chl b after excitation of any Chls, due to a near-perfect match of Chl b B band absorption with Chl a B state emission spectra. This results in an efficient depletion of the Chl a population (0.66 per IC/EET step, as compared to 0.21 in a Chl a-only system). Since Chl b only occurs in the peripheral antenna complexes of plants and algae, and RCs contain only Chl a, this would automatically trap potentially dangerous B state population in the antennae, preventing forwarding to the RCs.

Excitation Energy Transfer between Higher Excited States of Photosynthetic Pigments: 1. Carotenoids Intercept and Remove B Band Excitations

Chlorophylls (Chls) are known for fast, subpicosecond internal conversion (IC) from ultraviolet/blue-absorbing ("B" or "Soret" states) to the energetically lower, red light-absorbing Q states. Consequently, excitation energy transfer (EET) in photosynthetic pigment-protein complexes involving the B states has so far not been considered. We present, for the first time, a theoretical framework for the existence of B-B EET in tightly coupled Chl aggregates such as photosynthetic pigment-protein complexes. We show that according to a Förster resonance energy transport (FRET) scheme, unmodulated B-B EET has an unexpectedly high range. Unsuppressed, it could pose an existential threat: the damage potential of blue light for photochemical reaction centers (RCs) is well-known. This insight reveals so far undescribed roles for carotenoids (Crts, this article) and Chl b (next article in this series) of possibly vital importance. Our model system is the photosynthetic antenna pigment-protein complex (CP29). Here, we show that the B → Q IC is assisted by the optically allowed Crt state (S2): The sequence is B → S2 (Crt, unrelaxed) → S2 (Crt, relaxed) → Q. This sequence has the advantage of preventing ∼39% of Chl-Chl B-B EET since the Crt S2 state is a highly efficient FRET acceptor. The B-B EET range and thus the likelihood of CP29 to forward potentially harmful B excitations toward the RC are thus reduced. In contrast to the B band of Chls, most Crt energy donation is energetically located near the Q band, which allows for 74/80% backdonation (from lutein/violaxanthin) to Chls. Neoxanthin, on the other hand, likely donates in the B band region of Chl b, with 76% efficiency. Crts thus act not only in their currently proposed photoprotective roles but also as a crucial building block for any system that could otherwise deliver harmful "blue" excitations to the RCs.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Vibrational Funnels for Energy Transfer in Organic Chromophores

Photoinduced intramolecular energy transfers in multichromophoric molecules involve nonadiabatic vibronic channels that act as energy transfer funnels. They commonly take place through specific directions of motion dictated by the nonadiabatic coupling vectors. Vibrational funnels may support persistent coherences between electronic states and sometimes delineate the presence of minor alternative energy transfer pathways. The ultimate confirmation of their role on the interchromophoric energy transfer can be achieved by performing nonadiabatic excited-state molecular dynamics simulations by selectively freezing the nuclear motions in question. Our results point out this strategy as a useful tool to identify and evaluate the impact of these vibrational funnels on the energy transfer processes and guide the in silico design of materials with tunable properties and enhanced functionalities. Our work encourages applications of this methodology to different chemical and biochemical processes such as reactive scattering and protein conformational changes, to name a few.

Signatures of intramolecular vibrational and vibronic Q[Formula: see text]-Q[Formula: see text] coupling effects in absorption and CD spectra of chlorophyll dimers

An electron-vibrational coupling model that includes the vibronic (non-adiabatic) coupling between the Q[Formula: see text] and Q[Formula: see text] transitions of chlorophyll (Chl), created by Reimers and coworkers (Scientific Rep. 3, 2761, 2013) is extended here to chlorophyll dimers with interchlorophyll excitonic coupling. The model is applied to a Chl a dimer of the water-soluble chlorophyll binding protein (WSCP). As for isolated chlorophyll, the vibronic coupling is found to have a strong influence on the high-frequency vibrational sideband in the absorption spectrum, giving rise to a band splitting. In contrast, in the CD spectrum the interplay of vibronic coupling and static disorder leads to a strong suppression of the vibrational sideband in excellent agreement with the experimental data. The conservative nature of the CD spectrum in the low-energy region is found to be caused by a delicate balance of the intermonomer excitonic coupling between the purely electronic Q[Formula: see text] transition and the Q[Formula: see text] transition involving intramolecular vibrational excitations on one hand and the coupling to higher-energy electronic transitions on the other hand.

Q-Band relaxation in chlorophyll: new insights from multireference quantum dynamics

The ultrafast relaxation within the Q-bands of chlorophyll plays a crucial role in photosynthetic light-harvesting. We investigate this process via nuclear and electronic quantum dynamics on multireference potential energy surfaces.

Back-and-Forth Energy Transfer during Electronic Relaxation in a Chlorin-Perylene Dyad

Donor-acceptor dyads represent a practical approach to tuning the photophysical properties of linear conjugated polymers in materials chemistry. Depending on the absorption wavelength, the acceptor and donor roles can be interchanged, and as such, the directionality of the energy transfer can be controlled. Herein, nonadiabatic excited state molecular dynamics simulations have been performed in an arylethylene-linked perylene-chlorin dyad. After an initial photoexcitation at the Soret band of chlorin, we observe an ultrafast sequential electronic relaxation to the lowest excited state. This process is accomplished through an efficient round-trip chlorin-to-perylene-to-chlorin energy transfer. It is characterized by successive intermittent localized and delocalized vibronic dynamics. Nonradiative relaxation takes place mainly through energy transfer events with perylene acting as a "heat sink" through which the nonradiative relaxation is efficiently funneled, and the excess energy is dispersed in a larger space of vibrational degrees of freedom. Thus, our findings suggest the use of donor-acceptor dyads as a useful strategy when one needs to deactivate an electronic excitation.

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S₂/S₃(D_2d) → S₁(D_8h) → S₀(D_4h) → S₀(D_2d) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S₁ plays a crucial role in the photochemistry. The obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations

The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Q_x and Q_y states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Q_x-Q_y crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.

Photoexcitation dynamics in perylene diimide dimers

We utilize first-principles theory to investigate photo-induced excited-state dynamics of functionalized perylene diimide. This class of materials is highly suitable for solar energy conversion because of the strong optical absorbance, efficient energy transfer, and chemical tunability. We couple time-dependent density functional theory to a recently developed time-resolved non-adiabatic dynamics approach based on a semi-empirical description. By studying the monomer and dimer, we focus on the role stacking plays on the time-scales associated with excited-state non-radiative relaxation from a high excitonic state to the lowest energy exciton. We predict that the time-scale for energy conversion in the dimer is significantly faster than that in the monomer when equivalent excited states are accounted for. Additionally, for the dimer, the decay from the second to the nearly degenerate lowest energy excited-state involves two time-scales: a rapid decay on the order of ∼10 fs followed by a slower decay of ∼100 fs. Analysis of the spatial localization of the electronic transition density during the internal conversion process points out the existence of localized states on individual monomers, indicating that the strength of thermal fluctuations exceeds electronic couplings between the states such that the exciton hops between localized states throughout the simulation.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Dimerization Behavior of Methyl Chlorophyllide a as the Model of Chlorophyll a in the Presence of Water Molecules-Theoretical Study

A dimerization of methyl chlorophyllide a molecules and a role of water in stabilization and properties of methyl chlorophyllide a dimers were studied by means of symmetry-adapted perturbation theory (SAPT), functional-group SAPT (F-SAPT), density-functional theory (DFT), and time-dependent DFT approaches. The quantification of various types of interactions, such as π-π stacking, coordinative, and hydrogen bonding by applying the F-SAPT energy decomposition scheme shows the major role of the magnesium atom and the pheophytin macrocycle in the stability of the complex. The examination of interaction energy components with respect to a mutual orientation of monomers and in the presence or absence of water molecules reveals that the dispersion energy is the main binding factor of the interaction, while water molecules tend to weaken the attraction between methyl chlorophyllide a species. The dimerization can be seen in computed UV-vis spectra, and results in a doubling of the lowest peaks, as compared to the monomer spectrum, and in an intensity rise of the lowest 1.8 and 2.4 eV peaks at a cost of the 3.5 eV peaks for the majority of dimer configurations. The complexation of water has little effect on the peaks' position; however, it affects the overall shape of simulated spectra through changes in peak intensities, which is strongly dependent on the structure of the complex. The VCD spectra for the dimers show several characteristic features attributed to the interaction of substituting groups and/or water ligand attached to macrocycle groups belonging to different monomers. VCD is sensitive to the type of the formed dimer, but not to the number of water molecules it contains. This and several other features, as well as the differential UV-vis spectra, may serve as the indicator of the presence of a given dimer structure in the experiment.

Multi-stimuli responsive properties switch by intra- and inter-molecular charge transfer constructed from triphenylamine derivative

The compound TPA-BI exhibited multi-responsive fluorescence behaviors caused by inter-molecular charge transfer (CT) and intra-molecular CT formation.

Exciton localization in excited-state dynamics of a tetracene trimer: a surface hopping LC-TDDFTB study

Excitons in the molecular aggregates of chromophores are key participants in important processes such as photosynthesis or the functioning of organic photovoltaic devices. Therefore, the exploration of exciton dynamics is crucial. Here we report on exciton localization during excited-state dynamics of the recently synthesized tetracene trimer [Liu et al., Org. Lett., 2017, 19, 580]. We employ the surface hopping approach to nonadiabatic molecular dynamics in conjunction with the long-range corrected time-dependent density functional tight binding (LC-TDDFTB) method [Humeniuk and Mitrić, Comput. Phys. Commun., 2017, 221, 174]. Utilizing a set of descriptors based on the transition density matrix, we perform comprehensive analysis of exciton dynamics. The obtained results reveal an ultrafast exciton localization to a single tetracene unit of the trimer during excited-state dynamics, along with exciton transfer between units.

An ab initio multiple cloning approach for the simulation of photoinduced dynamics in conjugated molecules

Multidimensional wave function: a superposition of Gaussian coherent states guided by Ehrenfest trajectories suited to clone and swap their electronic amplitudes.

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.