Two-dimensional IR-spectroscopy: polarization anisotropy effects

Imaging the Photochemistry of the Hydrogen-Bonded Heptazine-Water Complex with Femtosecond Time-Resolved Spectroscopy: A Computational Study

Graphitic carbon nitride (g-CN) has attracted vast interest as a promising inexpensive metal-free photocatalyst for water splitting with solar photons. The heptazine (Hz) molecule is the building block of graphitic carbon nitride. The photochemistry of the Hz molecule and derivatives thereof in protic environments has been the subject of several recent experimental and computational studies. In the present work, the hydrogen-bonded Hz···H2O complex was adopted as a model system for the exploration of photoinduced electron and proton transfer processes in this complex with quasi-classical nonadiabatic trajectory simulations, using the ab initio ADC(2) electronic-structure method and a computationally efficient surface-hopping algorithm. The population of the optically excited bright 1ππ* state of the Hz chromophore relaxes through three 1nπ* states and a low-lying charge-transfer state, which drives proton transfer from H2O to Hz, to the long-lived optically dark S1(ππ*) state of Hz. The imaging of this ultrafast and complex dynamics with femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectroscopy and two-dimensional (2D) electronic spectroscopy (ES) was computationally explored in the framework of the quasi-classical doorway-window approximation. By comparison of the spectra of the Hz···H2O complex with those of the free Hz molecule, the effects of the hydrogen bond on the ultrafast internal conversion dynamics can be identified in the spectroscopic signals. Albeit the TA PP and 2D ES spectroscopies are primarily sensitive to electronic excited-state dynamics and less so to proton transfer dynamics, they nevertheless can provide mechanistic insights which can contribute to the acceleration of the optimization of photocatalysts for water splitting.

The first step of cyanine dye self-assembly: Dimerization

Self-assembling amphiphilic cyanine dyes, such as C8S3, are promising candidates for energy storage and optoelectronic applications due to their efficient energy transport properties. C8S3 is known to self-assemble in water into double-walled J-aggregates. Thus far, the molecular self-assembly steps remain shrouded in mystery. Here, we employ a multiscale approach to unravel the first self-assembly step: dimerization. Our multiscale approach combines molecular dynamics simulations with quantum chemistry calculations to obtain a Frenkel exciton Hamiltonian, which we then use in spectral calculations to determine the absorption and two-dimensional electronic spectra of C8S3 monomer and dimer systems. We model these systems solvated in both water and methanol, validating our model with experiments in methanol solution. Our theoretical results predict a measurable anisotropy decay upon dimerization, which is experimentally confirmed. Our approach provides a tool for the experimental probing of dimerization. Moreover, molecular dynamics simulations reveal that the dimer conformation is characterized by the interaction between the hydrophobic aliphatic tails rather than the π-π stacking previously reported for other cyanine dyes. Our results pave the way for future research into the mechanism of molecular self-assembly in similar light-harvesting complexes, offering valuable insights for understanding and optimizing self-assembly processes for various (nano)technological applications.

The Long and Short of Coupling and Uncoupling via 2D IR Spectroscopy

Determining dynamic structural changes along with the functional movements in biological systems has been a significant challenge for scientists for several decades. Utilizing vibrational coupling with the aid of 2D IR probe pairs has aided in uncovering structural dynamics and functional roles of chemical moieties involved in actions such as membrane peptide folding and transport, ion and water transport, and drug-protein interactions. Both native and non-native vibrational probe pairs have been developed for infrared studies, and their efficacy has been tested in various systems. With these probe pairs, 2D IR spectroscopy captures frozen snapshots of the structural events involved in biological function through vibrational coupling and correlated spectral diffusion. In this Perspective, different treatments of vibrational coupling and coupling models will be addressed, and a review of some of the specific vibrational probe pairs used to study these coupling mechanisms is presented. Overall, the intrinsic molecular dynamics detected on these ultrafast time scales will provide an atomic level view of how chosen structures traverse reaction paths. Thus, it is important to evaluate and assess the accuracy of the different vibrational coupling models and their consistency with the prediction of different molecular structures.

Two-Dimensional Infrared Spectroscopy Resolves the Vibrational Landscape in Donor-Bridge-Acceptor Complexes with Site-Specific Isotopic Labeling

Donor-bridge-acceptor complexes (D-B-A) are important model systems for understanding of light-induced processes. Here, we apply two-color two-dimensional infrared (2D-IR) spectroscopy to D-B-A complexes with a trans-Pt(II) acetylide bridge (D-C≡C-Pt-C≡C-A) to uncover the mechanism of vibrational energy redistribution (IVR). Site-selective 13C isotopic labeling of the bridge is used to decouple the acetylide modes positioned on either side of the Pt-center. Decoupling of the D-acetylide- from the A-acetylide- enables site-specific investigation of vibrational energy transfer (VET) rates, dynamic anharmonicities, and spectral diffusion. Surprisingly, the asymmetrically labeled D-B-A still undergoes intramolecular IVR between acetylide groups even though they are decoupled and positioned across a heavy atom usually perceived as a "vibrational bottleneck". Further, the rate of population transfer from the bridge to the acceptor was both site-specific and distance dependent. We show that vibrational excitation of the acetylide modes is transferred to ligand-centered modes on a subpicosecond time scale, followed by VET to solvent modes on the time scale of a few picoseconds. We also show that isotopic substitution does not affect the rate of spectral diffusion, indicating that changes in the vibrational dynamics are not a result of differences in local environment around the acetylides. Oscillations imprinted on the decay of the vibrationally excited acceptor-localized carbonyl modes show they enter a coherent superposition of states after excitation that dephases over 1-2 ps, and thus cannot be treated as independent in the 2D-IR spectra. These findings elucidate the vibrational landscape governing IR-mediated electron transfer and illustrate the power of isotopic labeling combined with multidimensional IR spectroscopy to disentangle vibrational energy propagation pathways in complex systems.

Signatures of coherent vibrational dynamics in ethylene carbonate

Despite having practical applications in battery technology and serving as a model system for Fermi resonance coupling, ethylene carbonate (EC) receives little direct attention as a vibrational probe in nonlinear vibrational spectroscopy experiments. EC contains a Fermi resonance that is well-characterized in the linear spectrum, and the environmental sensitivity of its Fermi resonance peaks could make it a good molecular probe for two-dimensional infrared spectroscopy (2DIR) experiments. As a model system, we investigate the linear and 2DIR vibrational spectrum of the carbonyl stretching region of ethylene carbonate in tetrahydrofuran. The 2DIR spectrum reveals peak dynamics that evolve coherently. We characterize these dynamics in the context of Redfield theory and find evidence that EC dynamics proceed through coherent pathways, including singular coherence transfer pathways that have not been widely observed in other studies. We find that coherent contributions play a significant role in the observed dynamics of cross-peaks in the 2DIR spectrum, which must be accounted for to extract accurate measurements of early waiting time dynamics.

On-the-Fly Simulation of Two-Dimensional Fluorescence-Excitation Spectra

Two-dimensional (2D) fluorescence-excitation (2D-FLEX) spectroscopy is a recently proposed nonlinear femtosecond technique for the detection of photoinduced dynamics. The method records a time-resolved fluorescence signal in its excitation- and detection-frequency dependence and hence combines the exclusive detection of excited state dynamics (fluorescence) with signals resolved in both excitation and emission frequencies (2D electronic spectroscopy). In this work, we develop an on-the-fly protocol for the simulation of 2D-FLEX spectra of molecular systems, which is based on interfacing the classical doorway-window representation of spectroscopic responses with trajectory surface hopping simulations. Applying this methodology to gas-phase pyrazine, we show that femtosecond 2D-FLEX spectra can deliver detailed information that is otherwise obtainable via attosecond spectroscopy.

Suppressing sidechain modes and improving structural resolution for 2D IR spectroscopy via vibrational lifetimes

Vibrational spectroscopy of protein structure often utilizes 13C18O-labeling of backbone carbonyls to further increase structural resolution. However, sidechains such as arginine, aspartate, and glutamate absorb within the same spectral region, complicating the analysis of isotope-labeled peaks. In this study, we report that the waiting time between pump and probe pulses in two-dimensional infrared spectroscopy can be used to suppress sidechain modes in favor of backbone amide I' modes based on differences in vibrational lifetimes. Furthermore, differences in the lifetimes of 13C18O-amide I' modes can aid in the assignment of secondary structure for labeled residues. Using model disordered and β-sheet peptides, it was determined that while β-sheets exhibit a longer lifetime than disordered structures, amide I' modes in both secondary structures exhibit longer lifetimes than sidechain modes. Overall, this work demonstrates that collecting 2D IR data at delayed waiting times, based on differences in vibrational lifetime between modes, can be used to effectively suppress interfering sidechain modes and further identify secondary structures.

Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe

Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.

Coarse-Grained Approach to Simulate Signatures of Excitation Energy Transfer in Two-Dimensional Electronic Spectroscopy of Large Molecular Systems

Two-dimensional electronic spectroscopy (2DES) has proven to be a highly effective technique in studying the properties of excited states and the process of excitation energy transfer in complex molecular assemblies, particularly in biological light-harvesting systems. However, the accurate simulation of 2DES for large systems still poses a challenge because of the heavy computational demands it entails. In an effort to overcome this limitation, we devised a coarse-grained 2DES method. This method encompasses the treatment of the entire system by dividing it into distinct weakly coupled segments, which are assumed to communicate predominantly through incoherent exciton transfer. We first demonstrate the efficiency of this method through simulation on a model dimer system, which demonstrates a marked improvement in calculation efficiency, with results that exhibit good concordance with reference spectra calculated with less approximate methods. Additionally, the application of this method to the light-harvesting antenna 2 (LH2) complex of purple bacteria showcases its advantages, accuracy, and limitations. Furthermore, simulating the anisotropy decay in LH2 induced by energy transfer and its comparison with experiments confirm that the method is capable of accurately describing dynamical processes in a biologically relevant system. This method presented lends itself to an extension that accounts for the effect of intrasegment relaxation processes on the 2DES spectra, which for computational efficiency are ignored in the implementation reported here. It is envisioned that the method will be employed in the future to accurately and efficiently calculate 2D spectra of more extensive systems, such as photosynthetic supercomplexes.

Theoretical model of femtosecond coherence spectroscopy of vibronic excitons in molecular aggregates

When used as pump pulses in transient absorption spectroscopy measurements, femtosecond laser pulses can produce oscillatory signals known as quantum beats. The quantum beats arise from coherent superpositions of the states of the sample and are best studied in the Fourier domain using Femtosecond Coherence Spectroscopy (FCS), which consists of one-dimensional amplitude and phase plots of a specified oscillation frequency as a function of the detection frequency. Prior works have shown ubiquitous amplitude nodes and π phase shifts in FCS from excited-state vibrational wavepackets in monomer samples. However, the FCS arising from vibronic-exciton states in molecular aggregates have not been studied theoretically. Here, we use a model of vibronic-exciton states in molecular dimers based on displaced harmonic oscillators to simulate FCS for dimers in two important cases. Simulations reveal distinct spectral signatures of excited-state vibronic-exciton coherences in molecular dimers that may be used to distinguish them from monomer vibrational coherences. A salient result is that, for certain relative orientations of the transition dipoles, the key resonance condition between the electronic coupling and the frequency of the vibrational mode may yield strong enhancement of the quantum-beat amplitude and, perhaps, also cause a significant decrease of the oscillation frequency to a value far lower than the vibrational frequency. Future studies using these results will lead to new insights into the excited-state coherences generated in photosynthetic pigment-protein complexes.

On-the-fly simulation of time-resolved fluorescence spectra and anisotropy

We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.

Observation of Coherent Symmetry-Breaking Vibration by Polarization-Dependent Femtosecond Spectroscopy

Understanding photoinduced chemical reactions beyond the Born-Oppenheimer paradigm requires a comprehensive examination of vibronic interactions. Although femtosecond studies have unveiled the influence of vibrational modes strongly coupled to ultrafast intramolecular reactions in the excited state, they often lack direct observations of how vibrations modulate electronic properties due to the rapid disappearance of reactants. To address this gap, our research investigates the dynamics of photoexcited molecules that do not react. Specifically, we focus on the coherent librational motion of molecular transition dipole moments, discovering that the coherent libration primarily originates from symmetry-breaking components in vibronically excited vibrational modes. Symmetry breaking motion can significantly impact the excited-state dynamics of highly symmetric molecules, potentially leading to nonadiabatic transitions. In essence, the data analysis framework introduced in this study can be harnessed to uncover potential reactivity in photoexcited molecules, further enhancing our understanding of the mechanisms governing these reactions.

Cross-phase modulation in two-dimensional spectroscopy

Developing from transient absorption (TA) spectroscopy, two-dimensional (2D) spectroscopy with pump-probe geometry has emerged as a versatile approach for alleviating the difficulty in implementing 2D spectroscopy with other geometries. However, the presence of cross-phase modulation (XPM) in TA spectroscopy introduces significant spectral distortions, particularly when the pump and probe pulses overlap. We demonstrate that this phenomenon is extended to the 2D spectroscopy with pump-probe geometry and the XPM is induced by the interference of the two pump pulses. We present the oscillatory behavior of XPM in the 2D spectrum and its displacement with respect to the waiting time delay through both experimental measurements and numerical simulations. Additionally, we explore the influence of probe pulse chirp on XPM and discover that by compressing the chirp, the impact of XPM on the desired signal can be reduced.

On selection rules in two-dimensional terahertz-infrared-visible spectroscopy

Two-dimensional terahertz-infrared-visible (2D TIRV) spectroscopy directly measures the coupling between quantum high-frequency vibrations and classical low-frequency modes of molecular motion. In addition to coupling strength, the signal intensity in 2D TIRV spectroscopy can also depend on the selection rules of the excited transitions. Here, we explore the selection rules in 2D TIRV spectroscopy by studying the coupling between the high-frequency CH3 stretching and low-frequency vibrations of liquid dimethyl sulfoxide (DMSO). Different excitation pathways are addressed using variations in laser pulse timing and different polarizations of exciting pulses and detected signals. The DMSO signals generated via different excitation pathways can be readily distinguished in the spectrum. The intensities of different excitation pathways vary unequally with changes in polarization. We explain how this difference stems from the intensities of polarized and depolarized Raman and hyper-Raman spectra of high-frequency modes. These results apply to various systems and will help design and interpret new 2D TIRV spectroscopy experiments.

Multimode vibrational dynamics and orientational effects in fluorescence-encoded infrared spectroscopy. I. Response function theory

Fluorescence-encoded infrared (FEIR) spectroscopy is an emerging technique for performing vibrational spectroscopy in solution with detection sensitivity down to single molecules. FEIR experiments use ultrashort pulses to excite a fluorescent molecule's vibrational and electronic transitions in a sequential, time-resolved manner, and are therefore sensitive to intervening vibrational dynamics on the ground state, vibronic coupling, and the relative orientation of vibrational and electronic transition dipole moments. This series of papers presents a theoretical treatment of FEIR spectroscopy that describes these phenomena and examines their manifestation in experimental data. This first paper develops a nonlinear response function description of Fourier-transform FEIR experiments for a two-level electronic system coupled to multiple vibrations, which is then applied to interpret experimental measurements in the second paper [L. Whaley-Mayda et al., J. Chem. Phys. 159, 194202 (2023)]. Vibrational coherence between pairs of modes produce oscillatory features that interfere with the vibrations' population response in a manner dependent on the relative signs of their respective Franck-Condon wavefunction overlaps, leading to time-dependent distortions in FEIR spectra. The orientational response of population and coherence contributions are analyzed and the ability of polarization-dependent experiments to extract relative transition dipole angles is discussed. Overall, this work presents a framework for understanding the full spectroscopic information content of FEIR measurements to aid data interpretation and inform optimal experimental design.

Probing exciton dynamics with spectral selectivity through the use of quantum entangled photons

Quantum light is increasingly recognized as a promising resource for developing optical measurement techniques. Particular attention has been paid to enhancing the precision of the measurements beyond classical techniques by using nonclassical correlations between quantum entangled photons. Recent advances in the quantum optics technology have made it possible to manipulate spectral and temporal properties of entangled photons, and photon correlations can facilitate the extraction of matter information with relatively simple optical systems compared to conventional schemes. In these respects, the applications of entangled photons to time-resolved spectroscopy can open new avenues for unambiguously extracting information on dynamical processes in complex molecular and materials systems. Here, we propose time-resolved spectroscopy in which specific signal contributions are selectively enhanced by harnessing nonclassical correlations of entangled photons. The entanglement time characterizes the mutual delay between an entangled twin and determines the spectral distribution of photon correlations. The entanglement time plays a dual role as the knob for controlling the accessible time region of dynamical processes and the degrees of spectral selectivity. In this sense, the role of the entanglement time is substantially equivalent to the temporal width of the classical laser pulse. The results demonstrate that the application of quantum entangled photons to time-resolved spectroscopy leads to monitoring dynamical processes in complex molecular and materials systems by selectively extracting desired signal contributions from congested spectra. We anticipate that more elaborately engineered photon states would broaden the availability of quantum light spectroscopy.

Two-dimensional fluorescence excitation spectroscopy: A novel technique for monitoring excited-state photophysics of molecular species with high time and frequency resolution

We propose a novel UV/Vis femtosecond spectroscopic technique, two-dimensional fluorescence-excitation (2D-FLEX) spectroscopy, which combines spectral resolution during the excitation process with exclusive monitoring of the excited-state system dynamics at high time and frequency resolution. We discuss the experimental feasibility and realizability of 2D-FLEX, develop the necessary theoretical framework, and demonstrate the high information content of this technique by simulating the 2D-FLEX spectra of a model four-level system and the Fenna-Matthews-Olson antenna complex. We show that the evolution of 2D-FLEX spectra with population time directly monitors energy transfer dynamics and can thus yield direct qualitative insight into the investigated system. This makes 2D-FLEX a highly efficient instrument for real-time monitoring of photophysical processes in polyatomic molecules and molecular aggregates.

Equivalence of quantum and classical third order response for weakly anharmonic coupled oscillators

Two-dimensional (2D) infrared (IR) spectra are commonly interpreted using a quantum diagrammatic expansion that describes the changes to the density matrix of quantum systems in response to light-matter interactions. Although classical response functions (based on Newtonian dynamics) have shown promise in computational 2D IR modeling studies, a simple diagrammatic description has so far been lacking. Recently, we introduced a diagrammatic representation for the 2D IR response functions of a single, weakly anharmonic oscillator and showed that the classical and quantum 2D IR response functions for this system are identical. Here, we extend this result to systems with an arbitrary number of bilinearly coupled, weakly anharmonic oscillators. As in the single-oscillator case, quantum and classical response functions are found to be identical in the weakly anharmonic limit or, in experimental terms, when the anharmonicity is small relative to the optical linewidth. The final form of the weakly anharmonic response function is surprisingly simple and offers potential computational advantages for application to large, multi-oscillator systems.

Molecular Vibrational Polariton Dynamics: What Can Polaritons Do?

ConspectusWhen molecular vibrational modes strongly couple to virtual states of photonic modes, new molecular vibrational polariton states are formed, along with a large population of dark reservoir modes. The polaritons are much like the bonding and antibonding molecular orbitals when atomic orbitals form molecular bonds, while the dark modes are like nonbonding orbitals. Because the polariton states are half-matter and half-light, whose energy is shifted from the parental states, polaritons are predicted to modify chemistry under thermally activated conditions, leading to an exciting and emerging field known as polariton chemistry that could potentially shift paradigms in chemistry. Despite several published results supporting this concept, the chemical physics and mechanism of polariton chemistry remain elusive. One reason for this challenge is that previous works cannot differentiate polaritons from dark modes. This limitation makes delineating the contributions to chemistry from polaritons and dark states difficult. However, this level of insight is critical for developing a solid mechanism for polariton chemistry to design and predict the outcome of strong coupling with any given reaction. My group addressed the challenge of differentiating the dynamics of polaritons and dark modes by ultrafast two-dimensional infrared (2D IR) spectroscopy. Specifically, (1) we found that polaritons can facilitate intra- and intermolecular vibrational energy transfer, opening a pathway to control vibrational energy flow in liquid-phase molecular systems, and (2) by studying a single-step isomerization event, we verified that indeed polaritons can modify chemical dynamics under strong coupling conditions, but in contrast, the dark modes behave like uncoupled molecules and do not change the dynamics. This finding confirmed the central concept of polariton chemistry: polaritons modify the potential energy landscape of reactions. The result also clarified the role of dark modes, which lays a critical foundation for designing cavities for future polariton chemistry. Aside from using 2D IR spectroscopy to study polariton chemistry, we also used the same technique to develop molecular polaritons into a potential quantum simulation platform. We demonstrated that polaritons have Rabi oscillations, and using a checkerboard cavity design, we showed that polaritons could have large nonlinearity across space. We further used the checkerboard polaritons to simulate coherence transfer and visualize it. A unidirectional coherence transfer was observed, indicating non-Hermitian dynamics. The highlighted efforts in this Account provide a solid understanding of the capability of polaritons for chemistry and quantum information science. I conclude this Account by discussing a few challenges for moving polariton chemistry toward being predictable and making the polariton quantum platform a complement to existing systems.

Simulation of two-dimensional infrared Raman spectroscopy with application to proteins

Two-dimensional infrared Raman spectroscopy is a powerful technique for studying the structure and interaction in molecular and biological systems. Here, we present a new implementation of the simulation of the two-dimensional infrared Raman signals. The implementation builds on the numerical integration of the Schrödinger equation approach. It combines the prediction of dynamics from molecular dynamics with a map-based approach for obtaining Hamiltonian trajectories and response function calculations. The new implementation is tested on the amide-I region for two proteins, where one is dominated by α-helices and the other by β-sheets. We find that the predicted spectra agree well with experimental observations. We further find that the two-dimensional infrared Raman spectra at least of the studied proteins are much less sensitive to the laser polarization used compared to conventional two-dimensional infrared experiments. The present implementation and findings pave the way for future applications for the interpretation of two-dimensional infrared Raman spectra.

Revealing Intermolecular Electronic and Vibronic Coherence with Polarization-Dependent Two-Dimensional Beating Maps

Two-dimensional electronic spectroscopy (2DES) has been widely employed as an efficient tool to reveal the impact of intermolecular electronic and/or vibronic quantum coherence on excitation energy transfer in light-harvesting complexes. However, intramolecular vibrational coherence would also contribute to oscillating signals in 2D spectra, along with the intermolecular coherence signals that are directly related to energy transfer. In this work, the possibility of screening the vibrational coherence signals is explored through polarization-dependent 2DES. The all-parallel (AP) and double-crossed (DC) polarization-dependent two-dimensional rephasing spectra (2DRS) are simulated for a minimalist heterodimer model with vibrational coupling. By combining the DC-2DRS and the 2D beating maps, we demonstrate that the population and vibrational coherence signals can be largely suppressed, resulting in highlighted intermolecular electronic and vibronic coherence signals. Moreover, the AP- and DC-2DBMs show rather different patterns at the vibrational frequency, indicating a possible way to identify pure vibrational coherence.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

Controlling Rotationally Resolved Two-Dimensional Infrared Spectra with Polarization

Recent advancements in infrared frequency combs will enable facile recording of coherent two-dimensional infrared spectra of gas-phase molecules with rotational resolution (RR2DIR). Using time-dependent density-matrix perturbation theory and angular momentum algebra techniques, we derive new polarization conditions unique to freely rotating molecules and absent in the condensed phase. These polarization conditions can be used to suppress parts of 2DIR rovibrational response, clarifying complicated RR2DIR spectra. With the polarization control methods described here, RR2DIR spectroscopy can be a powerful tool for studying complex gas mixtures of polyatomic molecules.

Resolving the Impact of Hydrogen Bonding on the Phylloquinone Cofactor through Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy (2DIR) was applied to phylloquinone (PhQ), an important biological cofactor, to elucidate the impact of hydrogen bonding on the ultrafast dynamics and energetics of the carbonyl stretching modes. 2DIR measurements were performed on PhQ dissolved in hexanol, which served as the hydrogen bonding solvent, and hexane, which served as a non-hydrogen bonding control. Molecular dynamics simulations and quantum chemical calculations were performed to aid in spectral assignment and interpretation. From the position of the peaks in the 2DIR spectra, we extracted the transition frequencies for the fundamental, overtone, and combination bands of hydrogen bonded and non-hydrogen bonded carbonyl groups of PhQ in the 1635-1680 cm^-1 region. We find that hydrogen bonding to a single carbonyl group acts to decouple the two carbonyl units of PhQ. Through analysis of the time-resolved 2DIR data, we find that hydrogen bonding leads to faster vibrational relaxation as well as an increase in the inhomogeneous broadening of the carbonyl groups. Overall, this work demonstrates how hydrogen bonding to the carbonyl groups of PhQ presents in the 2DIR spectra, laying the groundwork to use PhQ as a 2DIR probe to characterize the ultrafast fluctuations in the local environment of natural photosynthetic complexes.

Cavity-enabled enhancement of ultrafast intramolecular vibrational redistribution over pseudorotation

Vibrational strong coupling (VSC) between molecular vibrations and microcavity photons yields a few polaritons (light-matter modes) and many dark modes (with negligible photonic character). Although VSC is reported to alter thermally activated chemical reactions, its mechanisms remain opaque. To elucidate this problem, we followed ultrafast dynamics of a simple unimolecular vibrational energy exchange in iron pentacarbonyl [Fe(CO) 5 ] under VSC, which showed two competing channels: pseudorotation and intramolecular vibrational-energy redistribution (IVR). We found that under polariton excitation, energy exchange was overall accelerated, with IVR becoming faster and pseudorotation being slowed down. However, dark-mode excitation revealed unchanged dynamics compared with those outside of the cavity, with pseudorotation dominating. Thus, despite controversies around thermally activated VSC modified chemistry, our work shows that VSC can indeed alter chemistry through a nonequilibrium preparation of polaritons.

Rotationally-Resolved Two-Dimensional Infrared Spectroscopy of CO2(g): Rotational Wavepackets and Angular Momentum Transfer

Angular momentum transfer and wavepacket dynamics of CO₂(g) were measured on the picosecond time scale using polarization-resolved two-dimensional infrared (2D-IR) spectroscopy. The dynamics of rotational levels up to J_max ≈ 50 are observed simultaneously at room temperature. Rotational wavepackets launched by the pump pulses cause oscillations in the intensity of individual peaks and beating patterns in the 2D-IR spectra. The structure of the rotationally resolved 2D-IR spectrum is explained using nonlinear response function theory. Spectral diffusion of the rotationally resolved 2D-IR peaks reveals information about angular momentum transfer. We demonstrate the ability to directly measure inelastic angular momentum dynamics simultaneously across the ∼50 thermally excited rotational levels over several hundred picoseconds.

Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM

We present a series of new implementations that we recently introduced in COBRAMM, the open-source academic software developed in our group. The goal of these implementations is to offer an automatized workflow and interface to simulate time-resolved transient absorption (TA) spectra of medium-to-big chromophore embedded in a complex environment. Therefore, the excited states absorption and the stimulated emission are simulated along nonadiabatic dynamics performed with trajectory surface hopping. The possibility of treating systems from medium to big size is given by the use of time-dependent density functional theory (TD-DFT) and the presence of the environment is taken into account employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. The full implementation includes a series of auxiliary scripts to properly setup the QM/MM system, the calculation of the wavefunction overlap along the dynamics for the propagation, the evaluation of the transition dipole moment at linear response TD-DFT level, and scripts to setup, run and analyze the TA from an ensemble of trajectories. Altogether, we believe that our implementation will open the door to the easily simulate the time-resolved TA of systems so far computationally inaccessible.

Phase stable, shot-to-shot measurement of third- and fifth-order two-quantum correlation spectra using a pulse shaper in the pump-probe geometry

We demonstrate the first phase stable measurement of a third-order 2Q spectrum using a pulse shaper in the pump-probe geometry. This measurement was achieved by permuting the time-ordering of the pump pulses, thus rearranging the signal pathways that are emitted in the probe direction. The third-order 2Q spectrum is self-heterodyned by the probe pulse. Using this method, one can interconvert between a 1Q experiment and a 2Q experiment by simply reprogramming a pulse shaper or delay stage. We also measure a fifth-order absorptive 2Q spectrum in the pump-probe geometry, which contains similar information as a third-order experiment but does not suffer from dispersive line shapes. To do so, we introduce methods to minimize saturation-induced artifacts of the pulse shaper, improving fifth-order signals. These techniques add new capabilities for 2D spectrometers that use pulse shapers in the pump-probe beam geometry.

The polarization dependence of 2D IR cross-peaks distinguishes parallel-stranded and antiparallel-stranded DNA G-quadruplexes

Guanine-rich nucleic acid sequences have a tendency to form four-stranded non-canonical motifs known as G-quadruplexes. These motifs may adopt a wide range of structures characterized by size, strand orientation, guanine base conformation, and fold topology. Using three K⁺-bound model systems, we show that vibrational coupling between guanine C6 = O and ring modes varies between parallel-stranded and antiparallel-stranded G-quadruplexes, and that such structures can be distinguished by comparison of the polarization dependences of cross-peaks in their two-dimensional infrared (2D IR) spectra. Combined with previously defined vibrational frequency trends, this analysis reveals key features of a 30-nucleotide unimolecular variant of the Bcl-2 proximal promoter that are consistent with its reported structure. This study shows that 2D IR spectroscopy is a convenient method for analyzing G-quadruplex structures that can be applied to complex sequences where traditional high-resolution methods are limited by solubility and disorder.

Nonresonant coherent two-dimensional spectroscopy

Coherent electronic two-dimensional spectroscopy is nowadays a matured experimental technique that monitors the time evolution of the studied sample after its resonant optical excitation. However, the experimental experience shows that even nonresonant interactions can provide detectable spectral contributions. These are often present as a weak parasitic signals originating in the solvent and/or cuvette walls underlying the resonant spectrum of the actual sample and as such they are usually discarded from the analysis. In this work, we adapt the formalism of double-sided Feynman diagrams for the needs of coherent two-dimensional spectroscopy in the nonresonant regime. We analytically calculate the third-order polarization of a two-level and several variants of three-level systems. As a result, we demonstrate the typical appearance of the optical Kerr-effect, cross-phase modulation, excited-state coherence, two-photon absorption and stimulated Raman scattering in the 2D spectrum. This provides a framework for studying these effects by means of coherent two-dimensional spectroscopy.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex

Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.

Theory of two-dimensional spectroscopy with intense laser fields

Two-dimensional vibrational and electronic spectroscopic observables of isotropically oriented molecular samples in solution are sensitive to laser field intensities and polarization. The third-order response function formalism predicts a signal that grows linearly with the field strength of each laser pulse, thus lacking a way of accounting for non-trivial intensity-dependent effects, such as saturation and finite bleaching. An analytical expression to describe the orientational part of the molecular response, which, in the weak-field limit, becomes equivalent to a four-point correlation function, is presented. Such an expression is evaluated for Liouville-space pathways accounting for diagonal and cross peaks for all-parallel and cross-polarized pulse sequences, in both the weak- and strong-field conditions, via truncation of a Taylor series expansion at different orders. The results obtained in the strong-field conditions suggest how a careful analysis of two-dimensional spectroscopic experimental data should include laser pulse intensity considerations when determining molecular internal coordinates.

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

Experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang-Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic-Vibrational (2D EV) spectroscopy and two-dimensional Vibrational-Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

Anisotropy in fifth-order exciton-exciton-interaction two-dimensional spectroscopy

Exciton-exciton-interaction two-dimensional (EEI2D) spectroscopy is a fifth-order variant of 2D electronic spectroscopy. It can be used to probe biexciton dynamics in molecular systems and to observe exciton diffusion in extended systems such as polymers or light-harvesting complexes. The exciton transport strongly depends on the geometrical and energetic landscape and its perturbations. These can be of both local character, such as molecular orientation and energetic disorder, and long-range character, such as polymer kinks and structural domains. In the present theoretical work, we investigate the anisotropy in EEI2D spectroscopy. We introduce a general approach for how to calculate the anisotropy by using the response-function formalism in an efficient way. In numerical simulations, using a Frenkel exciton model with Redfield-theory dynamics, we demonstrate how the measurement of anisotropy in EEI2D spectroscopy can be used to identify various geometrical effects on exciton transport in dimers and polymers. Investigating a molecular heterodimer as an example, we demonstrate the utility of anisotropy in EEI2D spectroscopy for disentangling dynamic localization and annihilation. We further calculate the annihilation in extended systems such as conjugated polymers. In a polymer, a change in the anisotropy provides a unique signature for exciton transport between differently oriented sections. We analyze three types of geometry variations in polymers: a kink, varying geometric and energetic disorder, and different geometric domains. Our findings underline that employing anisotropy in EEI2D spectroscopy provides a way to distinguish between different geometries and can be used to obtain a better understanding of long-range exciton transport.

Observation of Ultrafast Coherence Transfer and Degenerate States with Polarization-Controlled Two-Dimensional Electronic Spectroscopy

Optical spectroscopy is a powerful tool to interrogate quantum states of matter. We present simulation results for the cross-polarized two-dimensional electronic spectra of the light-harvesting system LH2 of purple bacteria. We identify a spectral feature on the diagonal, which we assign to ultrafast coherence transfer between degenerate states. The implication for the interpretation of previous experiments on different systems and the potential use of this feature are discussed. In particular, we foresee that this kind of feature will be useful for identifying mixed degenerate states and for identifying the origin of symmetry breaking disorder in systems like LH2. Furthermore, this may help identify both vibrational and electronic states in biological systems such as proteins and solid-state materials such as hybrid perovskites.

Effect of disorder and polarization sequences on two-dimensional spectra of light-harvesting complexes

Two-dimensional electronic spectra (2DES) provide unique ways to track the energy transfer dynamics in light-harvesting complexes. The interpretation of the peaks and structures found in experimentally recorded 2DES is often not straightforward, since several processes are imaged simultaneously. The choice of specific pulse polarization sequences helps to disentangle the sometimes convoluted spectra, but brings along other disturbances. We show by detailed theoretical calculations how 2DES of the Fenna-Matthews-Olson complex are affected by rotational and conformational disorder of the chromophores.

A Proposed Method to Obtain Surface Specificity with Pump-Probe and 2D Spectroscopies

Surfaces and interfaces are ubiquitous in nature. From cell membranes, to photovoltaic thin films, surfaces have important function in both biological and materials systems. Spectroscopic techniques have been developed to probe systems like these, such as sum frequency generation (SFG) spectroscopies. The advantage of SFG spectroscopy, a second-order spectroscopy, is that it can distinguish between signals produced from molecules in the bulk versus on the surface. We propose a polarization scheme for third-order spectroscopy experiments, such as pump-probe and 2D spectroscopy, to select for surface signals and not bulk signals. This proposed polarization condition uses one pulse perpendicular compared to the other three to isolate cross-peaks arising from molecules with polar and uniaxial (i.e., biaxial) order at a surface, while removing the signal from bulk isotropic molecules. In this work, we focus on two of these cases: XXXY and YYYX, which differ by the sign of the cross-peak they create. We compare this technique to SFG spectroscopy and vibrational circular dichroism to provide insight to the behavior of the cross-peak signal. We propose that these singularly cross-polarized schemes provide odd-ordered spectroscopies the surface-specificity typically associated with even-ordered techniques.

Intramolecular hydrogen bonding protects the hydroxyl group from attack by fluctuating solvent forces

Ultrafast spectroscopy of molecular systems involving hydrogen- (H-) bonding has been at the forefront of fundamental chemical and physical research for several decades. Among the spectroscopic observables of the ultrafast dynamics is the pure dephasing of vibrationally excited molecules. Using third-order nonlinear vibrational spectroscopy, including polarization-selective transient grating measurements of vibrational lifetime and orientational diffusion as well as two-dimensional infrared spectroscopy, we determined different individual line shape components of hydroxyl stretching (νOH) excitations in a homologous series of chlorophenols and obtained the corresponding pure dephasing rates. The pure dephasing rates are correlated with vibrational anharmonicity of the νOH mode, which is tuned remotely from the hydroxyl site by changing the position of the chlorine substituents on the phenol ring. We found that in molecules where the hydroxyl group is in its free form, the pure dephasing rates scale linearly with the mode's anharmonicity such that assuming it is dominated by the third-order diagonal term, the ultrafast dynamics follow the prediction of the Kubo-Oxtoby theory. However, in the intramolecularly H-bonded ortho-chlorophenols, this trend is reversed, and the pure dephasing slows down by ∼50% for an increase in anharmonicity of only a few wavenumbers. Because the νOH mode's anharmonicity is known to reflect the H-bonding strength, our results suggest that intramolecular H-bonding can serve as a mechanism of protection from fluctuating forces exerted by the solvent. Such an effect can be relevant for ultrafast dynamics in biomolecules, where H-bonding plays a central role.

Determining the Orientation and Vibronic Couplings between Electronic and Vibrational Coordinates with Polarization-Selective Two-Dimensional Vibrational-Electronic Spectroscopy

We experimentally demonstrate polarization-selective two-dimensional (2D) vibrational-electronic (VE) spectroscopy on a transition-metal mixed-valence complex where the cyanide stretching vibrations are coupled to the metal-to-metal charge-transfer transition. A simultaneous fitting of the parallel and crossed polarized 2D VE spectra quantifies the relative vibronic coupling strengths and angles between the charge-transfer transition and three coupled cyanide stretching vibrations in a mode-specific manner. In particular, we find that the bridging vibration, which modulates the distance between the transition-metal centers, is oriented nearly parallel to the charge-transfer axis and is 9 times more strongly coupled to the electronic transition than the radial vibration, which is oriented almost perpendicular to the charge-transfer axis. The results from this experiment allow us to map the spectroscopically observed vibronic coordinates onto the molecular frame providing a general method to spatially resolve vibronic energy transfer on a femtosecond time scale.

Vibronic structure of photosynthetic pigments probed by polarized two-dimensional electronic spectroscopy and ab initio calculations

Bacteriochlorophyll a (Bchl a) and chlorophyll a (Chl a) play important roles as light absorbers in photosynthetic antennae and participate in the initial charge-separation steps in photosynthetic reaction centers. Despite decades of study, questions remain about the interplay of electronic and vibrational states within the Q-band and its effect on the photoexcited dynamics. Here we report results of polarized two-dimensional electronic spectroscopic measurements, performed on penta-coordinated Bchl a and Chl a and their interpretation based on state-of-the-art time-dependent density functional theory calculations and vibrational mode analysis for spectral shapes. We find that the Q-band of Bchl a is comprised of two independent bands, that are assigned following the Gouterman model to Q x and Q y states with orthogonal transition dipole moments. However, we measure the angle to be ∼75°, a finding that is confirmed by ab initio calculations. The internal conversion rate constant from Q x to Q y is found to be 11 ps-1. Unlike Bchl a, the Q-band of Chl a contains three distinct peaks with different polarizations. Ab initio calculations trace these features back to a spectral overlap between two electronic transitions and their vibrational replicas. The smaller energy gap and the mixing of vibronic states result in faster internal conversion rate constants of 38-50 ps-1. We analyze the spectra of penta-coordinated Bchl a and Chl a to highlight the interplay between low-lying vibronic states and their relationship to photoinduced relaxation. Our findings shed new light on the photoexcited dynamics in photosynthetic systems where these chromophores are primary pigments.

Advances and opportunities in ultrafast X-ray crystallography and ultrafast structural optical crystallography of nuclear and electronic protein dynamics

Both nuclear and electronic dynamics contribute to protein function and need multiple and complementary techniques to reveal their ultrafast structural dynamics response. Real-space information obtained from the measurement of electron density dynamics by X-ray crystallography provides aspects of both, while the molecular physics of coherence parameters and frequency-frequency correlation needs spectroscopy methods. Ultrafast pump-probe applications of protein dynamics in crystals provide real-space information through direct X-ray crystallographic structure analysis or through structural optical crystallographic analysis. A discussion of methods of analysis using ultrafast macromolecular X-ray crystallography and ultrafast nonlinear structural optical crystallography is presented. The current and future high repetition rate capabilities provided by X-ray free electron lasers for ultrafast diffraction studies provide opportunities for optical control and optical selection of nuclear coherence which may develop to access higher frequency dynamics through improvements of sensitivity and time resolution to reveal coherence directly. Specific selection of electronic coherence requires optical probes, which can provide real-space structural information through photoselection of oriented samples and specifically in birefringent crystals. Ultrafast structural optical crystallography of photosynthetic energy transfer has been demonstrated, and the theory of two-dimensional structural optical crystallography has shown a method for accessing the structural selection of electronic coherence.

Insights into the mechanisms and dynamics of energy transfer in plant light-harvesting complexes from two-dimensional electronic spectroscopy

During the past two decades, two-dimensional electronic spectroscopy (2DES) and related techniques have emerged as a potent experimental toolset to study the ultrafast elementary steps of photosynthesis. Apart from the highly engaging albeit controversial analysis of the role of quantum coherences in the photosynthetic processes, 2DES has been applied to resolve the dynamics and pathways of energy and electron transport in various light-harvesting antenna systems and reaction centres, providing unsurpassed level of detail. In this paper we discuss the main technical approaches and their applicability for solving specific problems in photosynthesis. We then recount applications of 2DES to study the exciton dynamics in plant and photosynthetic light-harvesting complexes, especially light-harvesting complex II (LHCII) and the fucoxanthin-chlorophyll proteins of diatoms, with emphasis on the types of unique information about such systems that 2DES is capable to deliver. This article is part of a Special Issue entitled Light harvesting, edited by Dr. Roberta Croce.