Coherent two-dimensional optical spectroscopy

Elucidating phonon dephasing mechanisms in layered perovskites with coherent Raman spectroscopies

Organic-inorganic hybrid perovskite quantum wells exhibit electronic structures with properties intermediate between those of inorganic semiconductors and molecular crystals. In these systems, periodic layers of organic spacer molecules occupy the interstitial spaces between perovskite sheets, thereby confining electronic excitations to two dimensions. Here, we investigate spectroscopic line broadening mechanisms for phonons coupled to excitons in lead-iodide layered perovskites with phenyl ethyl ammonium (PEA) and azobenzene ethyl ammonium (AzoEA) spacer cations. Using a modified Elliot line shape analysis for the absorbance and photoluminescence spectra, polaron binding energies of 11.2 and 17.5 meV are calculated for (PEA)2PbI4 and (AzoEA)2PbI4, respectively. To determine whether the polaron stabilization processes influence the dephasing mechanisms of coupled phonons, five-pulse coherent Raman spectroscopies are applied to the two systems under electronically resonant conditions. The prominence of inhomogeneous line broadening mechanisms detected in (AzoEA)2PbI4 suggests that thermal fluctuations involving the deformable organic phase broaden the distributions of phonon frequencies within the quantum wells. In addition, our data indicate that polaron stabilization primarily involves photoinduced reorganization of the organic phases for both systems, whereas the impulsively excited phonons represent less than 10% of the total polaron binding energy. The signal generation mechanisms associated with our fifth-order coherent Raman experiments are explored with a perturbative model in which cumulant expansions are used to account for time-coincident vibrational dephasing and polaron stabilization processes.

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.

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Ultrabroadband two-dimensional electronic spectroscopy in the pump-probe geometry using conventional optics

We report ultrabroadband two-dimensional electronic spectroscopy (2D ES) measurements obtained in the pump-probe geometry using conventional optics. A phase-stabilized Michelson interferometer provides the pump-pulse delay interval, τ1, necessary to obtain the excitation-frequency dimension. Spectral resolution of the probe beam provides the detection-frequency dimension, ω3. The interferometer incorporates active phase stabilization via a piezo stage and feedback from interference of a continuous-wave reference laser detected in quadrature. To demonstrate the method, we measured a well-characterized laser dye sample and obtained the known peak structure. The vibronic peaks are modulated as a function of the waiting time, τ2, by vibrational wave packets. The interferometer simplifies ultrabroadband 2D ES measurements and analysis.

Direct Probe of Vibrational Fingerprint and Combination Band Coupling

Vibrational fingerprints and combination bands are a direct measure of couplings that control molecular properties. However, most combination bands possess small transition dipoles. Here we use multiple, ultrafast coherent infrared pulses to resolve vibrational coupling between CH3CN fingerprint modes at 918 and 1039 cm-1 and combination bands in the 2750-6100 cm-1 region via doubly vibrationally enhanced (DOVE) coherent multidimensional spectroscopy (CMDS). This approach provides a direct probe of vibrational coupling between fingerprint modes and near-infrared combination bands of large and small transition dipoles in a molecular system over a large frequency range.

High-Resolution Frequency-Domain Spectroscopic and Modeling Studies of Photosystem I (PSI), PSI Mutants and PSI Supercomplexes

Photosystem I (PSI) is one of the two main pigment-protein complexes where the primary steps of oxygenic photosynthesis take place. This review describes low-temperature frequency-domain experiments (absorption, emission, circular dichroism, resonant and non-resonant hole-burned spectra) and modeling efforts reported for PSI in recent years. In particular, we focus on the spectral hole-burning studies, which are not as common in photosynthesis research as the time-domain spectroscopies. Experimental and modeling data obtained for trimeric cyanobacterial Photosystem I (PSI3), PSI3 mutants, and PSI3-IsiA18 supercomplexes are analyzed to provide a more comprehensive understanding of their excitonic structure and excitation energy transfer (EET) processes. Detailed information on the excitonic structure of photosynthetic complexes is essential to determine the structure-function relationship. We will focus on the so-called "red antenna states" of cyanobacterial PSI, as these states play an important role in photochemical processes and EET pathways. The high-resolution data and modeling studies presented here provide additional information on the energetics of the lowest energy states and their chlorophyll (Chl) compositions, as well as the EET pathways and how they are altered by mutations. We present evidence that the low-energy traps observed in PSI are excitonically coupled states with significant charge-transfer (CT) character. The analysis presented for various optical spectra of PSI3 and PSI3-IsiA18 supercomplexes allowed us to make inferences about EET from the IsiA18 ring to the PSI3 core and demonstrate that the number of entry points varies between sample preparations studied by different groups. In our most recent samples, there most likely are three entry points for EET from the IsiA18 ring per the PSI core monomer, with two of these entry points likely being located next to each other. Therefore, there are nine entry points from the IsiA18 ring to the PSI3 trimer. We anticipate that the data discussed below will stimulate further research in this area, providing even more insight into the structure-based models of these important cyanobacterial photosystems.

Strategies for phase matching control in a multidimensional Floquet state spectroscopy

Floquet state spectroscopy is an optical analogue of multiple quantum coherence nuclear magnetic resonance (MQC-NMR). Tunable ultrafast excitation pulses resonantly excite multiple states in a sample to form the Floquet state. The Floquet state emits multiple coherent beams at frequencies and in directions that conserve energy and momenta. The different output beams differ in the time ordering and coherences created by the excitation beams. They correspond to the different methodologies in the NMR family. Isolating a specific beam and monitoring the output intensity as a function of excitation frequencies creates multidimensional spectra containing cross-peaks between coupled states. The frequency range of the multidimensional spectra is limited by phase matching constraints. This paper presents a new, to the best of our knowledge, active phase matching strategy that increases the versatility of multidimensional Floquet state spectroscopy through both longer sample path lengths and larger spectral ranges.

Efficient formulation of multitime generalized quantum master equations: Taming the cost of simulating 2D spectra

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light-matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori-Nakajima-Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method's capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O(103) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies.

Water inside the Selectivity Filter of a K+ Ion Channel: Structural Heterogeneity, Picosecond Dynamics, and Hydrogen Bonding

Water inside biological ion channels regulates the key properties of these proteins, such as selectivity, ion conductance, and gating. In this article, we measure the picosecond spectral diffusion of amide I vibrations of an isotope-labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100-2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope-labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D line shapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent or nonadjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations was observed on a picosecond timescale. These dynamics are in stark contrast with liquid water, which remains highly dynamic even in the presence of ions at high concentrations.

Homogeneous Dephasing in Photosynthetic Bacterial Reaction Centers: Time Correlation Function Approach

A new tractable linear electronic transition dipole moment time correlation function (ETDMTCF) that accurately accounts for electronic dephasing, asymmetry, and width of 1-phonon profile, which the zero-phonon line (ZPL) contributes to it, in Rhodopseudomonas viridis bacterial reaction center is derived. This time correlation function proves to be superior to other frequency-domain expressions in case of strong electron-phonon coupling (which is often the case in bacterial RCs and pigment-protein complexes), many vibrational modes involved, and high temperature, whereby more vibronic and electronic (sequence) transitions would arise. The Fourier transform of this ETDMTCF leads to asymmetric multiphonon profiles composed of Lorentzian distribution and Gaussian distribution on the high- and low-energy sides, respectively, whereby the overtone widths fold themselves with that of the one-phonon profile. This ETDMTCF also features expedient computation in large systems using asymmetric phonon profiles to account correctly for dephasing and pigment-protein interaction (electron-phonon coupling). The derived ETDMTCF allows computing all nonlinear optical signals in both time and frequency domains, through the nonlinear dipole moment time correlation functions (as guided by nonlinear optical response theory) in line with the eight Liouville space pathways. The linear transition dipole moment time correlation function is of a central value as the nonlinear transition dipole moment time correlation function is expressed in terms of the linear transition dipole moment time correlation function, derived herein. One of the great advantages of presenting this ETDMTCF is its applicability to nonlinear transition dipole moment time correlation functions in line with the eight Liouville space pathways needed in computing nonlinear signals. As such, there is more to the utility and applicability of the presented ETDMTCF besides computational expediency and efficiency. Results show good agreement with the reported literature. The intimate connection between a one-phonon profile and the corresponding bath spectral density in photosynthetic complexes is discussed.

Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.

Theoretical Study of the Two-Dimensional Vibrational Sum Frequency Generation Spectroscopy of the Air-Water Interface at Varying Temperature and Its Connections to the Interfacial Structure and Dynamics

We performed a theoretical study of the temperature variation of two-dimensional vibrational sum frequency generation (2D-VSFG) spectra of the OH stretch modes at air-water interfaces in the mid-IR region. The calculations are performed at four different temperatures from 250 to 325 K by using a combination of techniques involving response function formalism of nonlinear spectroscopy, electronic structure calculations, and molecular dynamics simulations. Also, the calculations are performed for isotopically dilute solutions so that the intra- and intermolecular coupling between the vibrational modes of interest can be ignored. We have established the connections of temperature variation of various frequency- and time-dependent features of the calculated spectra to the changes in the underlying structure and dynamics of the interfaces. The results reveal that interfacial water is dynamically more heterogeneous than bulk water, with three dominant dynamical processes exhibiting their corresponding time-dependent features in the 2D-VSFG spectrum. These are the spectral diffusion of hydrogen-bonded OH groups at the interface, conversion of an initially hydrogen-bonded OH group to a dangling OH which is a stable state for surface water, unlike the bulk water, and the third one, which involves the conversion of an initially free or dangling OH group to its hydrogen-bonded state at the interface. The temporal appearance of the cross peaks corresponding to interconversion of the hydrogen-bonded state to the dangling state or vice versa of an interfacial OH group is found to take place at a slower rate than the dynamics of spectral diffusion of hydrogen-bonded molecules at the interface, which, in turn, is slower than the corresponding spectral diffusion of bulk water molecules. The temperature variation of these dynamic processes can be linked to the decay of appropriate hydrogen-bond and non-hydrogen-bond time correlation functions of interfacial water molecules for the different air-water systems studied in this work.

Water inside the selectivity filter of a K+ ion channel: structural heterogeneity, picosecond dynamics, and hydrogen-bonding

Water inside biological ion channels regulates the key properties of these proteins such as selectivity, ion conductance, and gating. In this Article we measure the picosecond spectral diffusion of amide I vibrations of an isotope labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100 - 2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D lineshapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent, or non-adjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations were observed on a picosecond timescale. These dynamics are in stark contrast with liquid water that remains highly dynamic even in the presence of ions at high concentrations.

Revisiting Ultrafast Dynamics in Carbonate-Based Electrolytes for Li-Ion Batteries: Clarifying 2D-IR Cross-Peak Interpretation

Understanding chemical exchange in carbonate-based electrolytes employed in Li-ion batteries (LIBs) is crucial for elucidating ion transport mechanisms. Ultrafast two-dimensional (2D) IR spectroscopy has been widely used to investigate the solvation structure and dynamics of Li-ions in organic carbonate-based electrolytes. However, the interpretation of cross-peaks observed in picosecond carbonyl stretch 2D-IR spectra has remained contentious. These cross-peaks could arise from various phenomena, including vibrational couplings between neighboring carbonyl groups in the first solvation shell around Li-ions, vibrational excitation transfers between carbonyl groups in distinct solvation environments, and local heating effects. Therefore, it is imperative to resolve the interpretation of 2D-IR cross-peaks to avoid misinterpretations regarding ultrafast dynamics found in LIB carbonate-based electrolytes. In this study, we have taken a comprehensive investigation of carbonate-based electrolytes utilizing 2D-IR spectroscopy and molecular dynamics (MD) simulations. Through meticulous analyses and interpretations, we have identified that the cross-peaks observed in the picosecond 2D-IR spectra of LIB electrolytes predominantly arise from intermolecular vibrational excitation transfer processes between the carbonyl groups of Li-bound and free carbonate molecules. We further discuss the limitations of employing a picosecond 2D-IR spectroscopic technique to study chemical exchange and intermolecular vibrational excitation transfer processes, particularly when the effects of the molecular photothermal process cannot be ignored. Our findings shed light on the dynamics of LIB electrolytes and resolve the controversy related to 2D-IR cross-peaks. By discerning the origin of these features, we could provide valuable insights for the design and optimization of next-generation Li-ion batteries.

Two-Dimensional Infrared Spectroscopy of Isolated Molecular Ions

Two-dimensional infrared (2D IR) spectroscopy of mass-selected, cryogenically cooled molecular ions is presented. Nonlinear response pathways, encoded in the time-domain photodissociation action response of weakly bound N2 messenger tags, were isolated using pulse shaping techniques following excitation with four collinear ultrafast IR pulses. 2D IR spectra of Re(CO)3(CH3CN)3+ ions capture off-diagonal cross-peak bleach signals between the asymmetric and symmetric carbonyl stretching transitions. These cross peaks display intensity variations as a function of pump-probe delay time due to coherent coupling between the vibrational modes. Well-resolved 2D IR features in the congested fingerprint region of protonated caffeine (C8H10N4O2H+) are also reported. Importantly, intense cross-peak signals were observed at 3 ps waiting time, indicating that tag-loss dynamics are not competing with the measured nonlinear signals. These demonstrations pave the way for more precise studies of molecular interactions and dynamics that are not easily obtainable with current condensed-phase methodologies.

Structural Dependence of Extended Amide III Vibrations in Two-Dimensional Infrared Spectra

Two-dimensional infrared (2D-IR) spectroscopy is a powerful experimental method for probing the structure and dynamics of proteins in aqueous solution. So far, most experimental studies have focused on the amide I vibrations, for which empirical vibrational exciton models provide a means of interpreting such experiments. However, such models are largely lacking for other regions of the vibrational spectrum. To close this gap, we employ an efficient quantum-chemical methodology for the calculation of 2D-IR spectra, which is based on anharmonic theoretical vibrational spectroscopy with localized modes. We apply this approach to explore the potential of 2D-IR spectroscopy in the extended amide III region. Using calculations for a dipeptide model as well as alanine polypeptides, we show that distinct 2D-IR cross-peaks in the extended amide III region can potentially be used to distinguish α-helix and β-strand structures. We propose that the extended amide III region could be a promising target for future 2D-IR experiments.

Machine Learning Approach Based on a Range-Corrected Deep Potential Model for Efficient Vibrational Frequency Computation

As an ensemble average result, vibrational spectrum simulation can be time-consuming with high accuracy methods. We present a machine learning approach based on the range-corrected deep potential (DPRc) model to improve the computing efficiency. The DPRc method divides the system into "probe region" and "solvent region"; "solvent-solvent" interactions are not counted in the neural network. We applied the approach to two systems: formic acid C═O stretching and MeCN C≡N stretching vibrational frequency shifts in water. All data sets were prepared using the quantum vibration perturbation approach. Effects of different region divisions, one-body correction, cut range, and training data size were tested. The model with a single-molecule "probe region" showed stable accuracy; it ran roughly 10 times faster than regular deep potential and reduced the training time by about four. The approach is efficient, easy to apply, and extendable to calculating various spectra.

Counteracting Effects of Trimethylamine N-Oxide against Urea in Aqueous Solutions: Insights from Theoretical Two-Dimensional Infrared Spectroscopy

The study of small osmolytes in their aqueous solutions has gained significant attention because of their relevance to structure and thermodynamics of proteins in aqueous media. Special attention has been given to the binary and ternary aqueous solutions of urea and trimethylamine N-oxide (TMAO). Urea is a well-known protein denaturant, while TMAO protects proteins in their native states. Interestingly, TMAO counteracts urea's ability to denature proteins when present in solutions with approximately half of the concentration of urea. Vibrational spectroscopy can improve our understanding of the molecular origin of this counteracting effect because of its sensitivity to local structure and dynamics. We present results of theoretical linear vibrational and two-dimensional infrared (2DIR) spectroscopy of water in the binary and ternary aqueous solutions of TMAO and urea. The 2DIR spectra are calculated using the electronic structure/molecular dynamics approach. The non-Condon effects in spectral transitions are incorporated in the theoretical calculations of 2DIR spectra. It is found that TMAO disrupts the local structure of water, while urea leaves it essentially unaffected. The 2DIR results show that both TMAO and urea slow down the dynamics of spectral diffusion of water. The extent of slowing down is found to be particularly significant for both hydration and bulk water in the presence of TMAO which can be attributed to strong hydrogen bonds between the water and TMAO molecules. The water molecules present in the hydration layer of the solutes in the ternary solutions are found to relax at even slower rates compared to that in their binary solutions in water. The hydrogen bonds between TMAO and urea are found to be not stable. Thus, the counteracting effect of TMAO against urea is seen to take place mainly through water-mediated interactions. Such TMAO-induced effects giving rise to more structured and slower hydrogen-bonded network are successfully captured through 2DIR spectroscopic calculations.

Progress and Prospects in Optical Ultrafast Microscopy in the Visible Spectral Region: Transient Absorption and Two-Dimensional Microscopy

Ultrafast optical microscopy, generally employed by incorporating ultrafast laser pulses into microscopes, can provide spatially resolved mechanistic insight into scientific problems ranging from hot carrier dynamics to biological imaging. This Review discusses the progress in different ultrafast microscopy techniques, with a focus on transient absorption and two-dimensional microscopy. We review the underlying principles of these techniques and discuss their respective advantages and applicability to different scientific questions. We also examine in detail how instrument parameters such as sensitivity, laser power, and temporal and spatial resolution must be addressed. Finally, we comment on future developments and emerging opportunities in the field of ultrafast microscopy.

Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure

Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. We show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.

Photoelectron spectroscopy with entangled photons; enhanced spectrotemporal resolution

In this theoretical study, we show how photoelectron signals generated by time-energy entangled photon pairs can monitor ultrafast excited state dynamics of molecules with high joint spectral and temporal resolutions, not limited by the Fourier uncertainty of classical light. This technique scales linearly, rather than quadratically, with the pump intensity, allowing the study of fragile biological samples with low photon fluxes. Since the spectral resolution is achieved by electron detection and the temporal resolution by a variable phase delay, this technique does not require scanning the pump frequency and the entanglement times, which significantly simplifies the experimental setup, making it feasible with current instrumentation. Application is made to the photodissociation dynamics of pyrrole calculated by exact nonadiabatic wave packet simulations in a reduced two nuclear coordinate space. This study demonstrates the unique advantages of ultrafast quantum light spectroscopy.

Theoretical Two Dimensional Infrared Spectroscopy of Aqueous Solutions of tert-Butyl Alcohol: Variation of the Dynamics of Spectral Diffusion along the Percolation Transition

Binary mixtures of water and tert-butyl alcohol (TBA) are known to exhibit the so-called percolation transition where small clusters of TBA molecules span into large aggregates beyond a threshold concentration of the alcohol. In the present study, we have investigated the linear and two-dimensional infrared spectral features of aqueous solutions of TBA for varying concentration of the alcohol along the percolation transition. The percolation transition is characterized through calculations of intermolecular radial distribution functions and average size of the largest cluster of TBA molecules. It is found that, with variation of alcohol concentration, the radial distribution functions of the central carbon atoms of TBA molecules show a nonmonotonic change in the height of the first peak and also the size of the largest cluster of TBA molecules show a jump in the increase of its size for TBA mole fraction between 0.04 and 0.06 corresponding to a transition from smaller clusters to larger spanning aggregates. However, it is found that the linear infrared spectrum of water does not exhibit any noticeable changes on variation of TBA concentration along the percolation transition. Subsequently, two-dimensional infrared (2DIR) spectra and vibrational frequency time correlation function of water are calculated for all the TBA-water solutions considered in this study. The spectral diffusion of water calculated from 2DIR is found to slow down with increase of the TBA concentration. The time scales of spectral diffusion of water, as characterized by the relaxation of frequency time correlation function, 2DIR metric of central line slope, and also the hydrogen bond time correlation functions, are found to exhibit a noticeable jump along the percolation transition. The hydrophilic group of TBA is found to retard the water dynamics more effectively than the hydrophobic groups. Also, the jump in the dynamical slowdown along the percolation transition is found to be more significant for water molecules at the hydrophilic sites.

Pulse overlap artifacts and double quantum coherence spectroscopy

The double quantum coherence (DQC) signal in nonlinear spectroscopy gives information about the many-body correlation effects not easily available by other methods. The signal is short-lived, consequently, a significant part of it is generated during the pulse overlap. Since the signal is at two times the laser frequency, one may intuitively expect that the pulse overlap-related artifacts are filtered out by the Fourier transform. Here, we show that this is not the case. We perform explicit calculations of phase-modulated two-pulse experiments of a two-level system where the DQC is impossible. Still, we obtain a significant signal at the modulation frequency, which corresponds to the DQC, while the Fourier transform over the pulse delay shows a double frequency. We repeat the calculations with a three-level system where the true DQC signal occurs. We conclude that with realistic dephasing times, the pulse-overlap artifact can be significantly stronger than the DQC signal. Our results call for great care when analyzing such experiments. As a rule of thumb, we recommend that only delays larger than 1.5 times the pulse length should be used.

Coherent Vibronic Wavepackets Show Structure-Directed Charge Flow in Host-Guest Donor-Acceptor Complexes

Designing and controlling charge transfer (CT) pathways in organic semiconductors are important for solar energy applications. To be useful, a photogenerated, Coulombically bound CT exciton must further separate into free charge carriers; direct observations of the detailed CT relaxation pathways, however, are lacking. Here, photoinduced CT and relaxation dynamics in three host-guest complexes, where a perylene (Per) electron donor guest is incorporated into two symmetric and one asymmetric extended viologen cyclophane acceptor hosts, are presented. The central ring in the extended viologen is either p-phenylene (ExV²⁺) or electron-rich 2,5-dimethoxy-p-phenylene (ExMeOV²⁺), resulting in two symmetric cyclophanes with unsubstituted or methoxy-substituted central rings, ExBox⁴⁺ and ExMeOBox⁴⁺, respectively, and an asymmetric cyclophane with one of the central viologen rings being methoxylated ExMeOVBox⁴⁺. Upon photoexcitation, the asymmetric host-guest ExMeOVBox⁴⁺ ⊃ Per complex exhibits directional CT toward the energetically unfavorable methoxylated side due to structural restrictions that facilitate strong interactions between the Per donor and the ExMeOV²⁺ side. The CT state relaxation pathways are probed using ultrafast optical spectroscopy by focusing on coherent vibronic wavepackets, which are used to identify CT relaxations along charge localization and vibronic decoherence coordinates. Specific low- and high-frequency nuclear motions are direct indicators of a delocalized CT state and the degree of CT character. Our results show that the CT pathway can be controlled by subtle chemical modifications of the acceptor host in addition to illustrating how coherent vibronic wavepackets can be used to probe the nature and time evolution of the CT states.

Temperature Dependence of Non-Condon Effects in Two-Dimensional Vibrational Spectroscopy of Water

Non-Condon effects in vibrational spectroscopy refers to the dependence of a molecule's vibrational transition dipole and polarizability on the coordinates of the surrounding environment. Earlier studies have shown that such effects can be pronounced for hydrogen-bonded systems like liquid water. Here, we present a theoretical study of two-dimensional vibrational spectroscopy under the non-Condon and Condon approximations at varying temperatures. We have performed calculations of both two-dimensional infrared and two-dimensional vibrational Raman spectra to gain insights into the temperature dependence of non-Condon effects in nonlinear vibrational spectroscopy. The two-dimensional spectra are calculated for the OH vibration of interest in the isotopic dilution limit where the coupling between the oscillators is ignored. Generally, both the infrared and Raman line shapes undergo red shifts with decrease in temperature due to strengthening of hydrogen bonds and decrease in the fraction of OH modes with weaker or no hydrogen bonds. The infrared line shape is further red-shifted under the non-Condon effects at a given temperature, while the Raman line shape does not show any such red shift due to non-Condon effects. The spectral dynamics becomes slower on decrease of temperature due to slower hydrogen bond relaxation and, for a given temperature, the spectral diffusion occurs at a faster rate upon inclusion of non-Condon effects. The time scales of spectral diffusion extracted from different metrics agree well with each other and also with experiments. The changes in the spectrum due to non-Condon effects are found to be more significant at lower temperatures.

Narrow homogeneous linewidths and slow cooling dynamics across infrared intra-band transitions in n-doped HgSe colloidal quantum dots

Intra-band transitions in colloidal quantum dots (QDs) are promising for opto-electronic applications in the mid-IR spectral region. However, such intra-band transitions are typically very broad and spectrally overlapping, making the study of individual excited states and their ultrafast dynamics very challenging. Here, we present the first full spectrum two-dimensional continuum infrared (2D CIR) spectroscopy study of intrinsically n-doped HgSe QDs, which exhibit mid-infrared intra-band transitions in their ground state. The obtained 2D CIR spectra reveal that underneath the broad absorption line shape of ∼500 cm^-1, the transitions exhibit surprisingly narrow intrinsic linewidths with a homogeneous broadening of 175-250 cm^-1. Furthermore, the 2D IR spectra are remarkably invariant, with no sign of spectral diffusion dynamics at waiting times up to 50 ps. Accordingly, we attribute the large static inhomogeneous broadening to the distribution of size and doping level of the QDs. In addition, the two higher-lying P-states of the QDs can be clearly identified in the 2D IR spectra along the diagonal with a cross-peak. However, there is no indication of cross-peak dynamics indicating that, despite the strong spin-orbit coupling in HgSe, transitions between the P-states must be longer than our maximum waiting time of 50 ps. This study illustrates a new frontier of 2D IR spectroscopy enabling the study of intra-band carrier dynamics in nanocrystalline materials across the entire mid-infrared spectrum.

Forced Interactions: Ionic Polymers at Charged Surfactant Interfaces

Characterizing electrostatic interactions at heterogeneous interfaces is critical for developing a fundamental description of the dynamic processes at charged interfaces. Water-in-oil reverse micelles (RMs) offer a high degree of tunability across composition, polarity, and temperature, making them ideal systems for studying interactions at heterogeneous liquid-liquid interfaces. In the present study, we use a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics (MD) simulations to determine the picosecond interfacial dynamics in RMs containing binary compositions of sorbitan monostearate and anionic or cationic cosurfactants, which are used to tune the ratio of charged to nonionic surfactants at the interface. The positively charged polyethylenimine (PEI) polymer is encapsulated within the RMs, and the carbonyl stretching mode of sorbitan monostearate reports on the interfacial hydrogen-bond populations and dynamics. The results show that hydrogen-bond populations are altered through the inclusion of both negatively and positively charged cosurfactants. Charged surfactants increase interfacial water penetration into the surfactant layer, and the surface localization of polymers decreases water penetration. Local hydrogen-bond dynamics undergo a slowdown with the inclusion of charged surfactants, and the encapsulation of polymers results in similar effects, irrespective of the charge.

An accurate and efficient Ehrenfest dynamics approach for calculating linear and nonlinear electronic spectra

Linear and nonlinear electronic spectra provide an important tool to probe the absorption and transfer of electronic energy. Here, we introduce a pure state Ehrenfest approach to obtain accurate linear and nonlinear spectra that is applicable to systems with large numbers of excited states and complex chemical environments. We achieve this by representing the initial conditions as sums of pure states and unfolding multi-time correlation functions into the Schrödinger picture. By doing this, we show that one can obtain significant improvements in accuracy over the previously used projected Ehrenfest approach and that these benefits are particularly pronounced in cases where the initial condition is a coherence between excited states. While such initial conditions do not arise when calculating linear electronic spectra, they play a vital role in capturing multidimensional spectroscopies. We demonstrate the performance of our method by showing that it is able to quantitatively capture the exact linear, 2D electronic spectroscopy, and pump-probe spectra for a Frenkel exciton model in slow bath regimes and is even able to reproduce the main spectral features in fast bath regimes.

Hierarchical Mapping for Efficient Simulation of Strong System-Environment Interactions

Quantum dynamics (QD) simulation is a powerful tool for interpreting ultrafast spectroscopy experiments and unraveling their microscopic mechanism in out-of-equilibrium excited state behaviors in various chemical, biological, and material systems. Although state-of-the-art numerical QD approaches such as the time-dependent density matrix renormalization group (TD-DMRG) already greatly extended the solvable system size of general linearly coupled exciton-phonon models with up to a few hundred phonon modes, the accurate simulation of larger system sizes or strong system-environment interactions is still computationally highly challenging. Based on quantum information theory (QIT), in this work, we realize that only a small number of effective phonon modes couple to the excitonic system directly regardless of a large or even infinite number of modes in the condensed phase environment. On top of the identified small number of direct effective modes, we propose a hierarchical mapping (HM) approach through performing block Lanczos transformations on the remaining indirect modes, which transforms the Hamiltonian matrix to a nearly block-tridiagonal form and eliminates the long-range interactions. Numerical tests on model spin-boson systems and realistic singlet fission models in a rubrene crystal environment with up to 7000 modes and strong system-environment interactions indicate HM can reduce the system size by 1-2 orders of magnitude and accelerate the calculation by ∼80% without losing accuracy.

Infrared Diffusion-Ordered Spectroscopy Reveals Molecular Size and Structure

Inspired by ideas from NMR, we have developed Infrared Diffusion-Ordered Spectroscopy (IR-DOSY), which simultaneously characterizes molecular structure and size. We rely on the fact that the diffusion coefficient of a molecule is determined by its size through the Stokes-Einstein relation, and achieve sensitivity to the diffusion coefficient by creating a concentration gradient and tracking its equilibration in an IR-frequency resolved manner. Analogous to NMR-DOSY, a two-dimensional IR-DOSY spectrum has IR frequency along one axis and diffusion coefficient (or equivalently, size) along the other, so the chemical structure and the size of a compound are characterized simultaneously. In an IR-DOSY spectrum of a mixture, molecules with different sizes are nicely separated into distinct sets of IR peaks. Extending this idea to higher dimensions, we also perform 3D-IR-DOSY, in which we combine the conformation sensitivity of femtosecond multi-dimensional IR spectroscopy with size sensitivity.

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.

In Situ Identification of Secondary Structures in Unpurified Bombyx mori Silk Fibrils Using Polarized Two-Dimensional Infrared Spectroscopy

The mechanical properties of biomaterials are dictated by the interactions and conformations of their building blocks, typically proteins. Although the macroscopic behavior of biomaterials is widely studied, our understanding of the underlying molecular properties is generally limited. Among the noninvasive and label-free methods to investigate molecular structures, infrared spectroscopy is one of the most commonly used tools because the absorption bands of amide groups strongly depend on protein secondary structure. However, spectral congestion usually complicates the analysis of the amide spectrum. Here, we apply polarized two-dimensional (2D) infrared spectroscopy (IR) to directly identify the protein secondary structures in native silk films cast from Bombyx mori silk feedstock. Without any additional peak fitting, we find that the initial effect of hydration is an increase of the random coil content at the expense of the helical content, while the β-sheet content is unchanged and only increases at a later stage. This paper demonstrates that 2D-IR can be a valuable tool for characterizing biomaterials.

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.

Ultrafast Two-Dimensional Infrared Spectroscopy Resolved a Structured Lysine 159 on the Cytoplasmic Surface of the Microbial Photoreceptor Bacteriorhodopsin

Bacteriorhodopsin (bR) is a light-driven microbial receptor, and lysine 159 (K159) is a charged residue on the cytoplasmic (CP) side of its E-F loop. However, its conformation and function remain unknown due to fast surface dynamics. By utilizing a ¹³C, ¹⁵N-labeled lysine (K) as an isotope probe, we created a network of site-specific amide-I vibrational signatures (backbone carbonyl stretch) to identify the frequency contribution of the labeled residues to the amide-I excitonic band structure. Thus, the red-shifted amide-I frequency in the ¹³C, ¹⁵N-lysine-labeled bR (uK-bR) to the unlabeled bR (WT-bR) could be differentiated and examined by ultrafast two-dimensional vibrational echo infrared (2D IR) spectroscopy. Our results showed that the backbone carbonyl of K159 is located at a high frequency of ca. 1693 cm^-1 and has a vibrational excited-state relaxation time shorter than the bulk helical amide-I mode at the same frequency, suggesting that K159 may possess a hydrogen-bonded γ-turn structure with E161, one of the carboxylate residues on the CP surface of bR. The 2D solid-state NMR study of uK-bR also revealed conformational dependent lysine residues, from which K159 was found to involve the turn motif. This γ-turn structure maintained by K159 may help to stabilize the E-F loop and support E161 in attracting protons from the bulk during the late stage of the bR photocycle. The combined spectroscopic approach illustrated in this work may be applied to map residue-specific local structures and dynamics of other receptors and large proteins.

Broadband Multidimensional Spectroscopy Identifies the Amide II Vibrations in Silkworm Films

We used two-dimensional infrared spectroscopy to disentangle the broad infrared band in the amide II vibrational regions of Bombyx mori native silk films, identifying the single amide II modes and correlating them to specific secondary structure. Amide I and amide II modes have a strong vibrational coupling, which manifests as cross-peaks in 2D infrared spectra with frequencies determined by both the amide I and amide II frequencies of the same secondary structure. By cross referencing with well-known amide I assignments, we determined that the amide II (N-H) absorbs at around 1552 and at 1530 cm-1 for helical and β-sheet structures, respectively. We also observed a peak at 1517 cm-1 that could not be easily assigned to an amide II mode, and instead we tentatively assigned it to a Tyrosine sidechain. These results stand in contrast with previous findings from linear infrared spectroscopy, highlighting the ability of multidimensional spectroscopy for untangling convoluted spectra, and suggesting the need for caution when assigning silk amide II spectra.

Calculating non-linear response functions for multi-dimensional electronic spectroscopy using dyadic non-Markovian quantum state diffusion

We present a methodology for simulating multi-dimensional electronic spectra of molecular aggregates with coupling of electronic excitation to a structured environment using the stochastic non-Markovian quantum state diffusion (NMQSD) method in combination with perturbation theory for the response functions. A crucial aspect of our approach is that we propagate the NMQSD equation in a doubled system Hilbert space, but with the same noise. We demonstrate that our approach shows fast convergence with respect to the number of stochastic trajectories, providing a promising technique for numerical calculation of two-dimensional electronic spectra of large molecular aggregates.

Mapping Simulated Two-Dimensional Spectra to Molecular Models Using Machine Learning

Two-dimensional (2D) spectroscopy encodes molecular properties and dynamics into expansive spectral data sets. Translating these data into meaningful chemical insights is challenging because of the many ways chemical properties can influence the spectra. To address the task of extracting chemical information from 2D spectroscopy, we study the capacity of simple feedforward neural networks (NNs) to map simulated 2D electronic spectra to underlying physical Hamiltonians. We examined hundreds of simulated 2D spectra corresponding to monomers and dimers with varied Franck-Condon active vibrations and monomer-monomer electronic couplings. We find the NNs are able to correctly characterize most Hamiltonian parameters in this study with an accuracy above 90%. Our results demonstrate that NNs can aid in interpreting 2D spectra, leading from spectroscopic features to underlying effective Hamiltonians.

Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations

Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.

Measuring Protein Conformation at Aqueous Interfaces with 2D Infrared Spectroscopy of Emulsions

Determining the secondary and tertiary structures of proteins at aqueous interfaces is crucial for understanding their function, but measuring these structures selectively at the interface is challenging. Here we demonstrate that two-dimensional infrared (2D-IR) spectroscopy of protein stabilized emulsions offers a new route to measuring interfacial protein structure with high levels of detail. We prepared hexadecane/water oil-in-water emulsions stabilized by model LK peptides that are known to fold into either α-helix or β-sheet conformations at hydrophobic interfaces and measured 2D-IR spectra in a transmission geometry. We saw clear spectral signatures of the peptides folding at the interface, with no detectable residue from remaining bulk peptides. Using 2D spectroscopy gives us access to correlation and dynamics data, which enables structural assignment in cases where linear spectroscopy fails. Using the emulsions allows one to study interfacial spectra with standard transmission geometry spectrometers, bringing the richness of 2D-IR to the interface with no additional optical complexity.

Stochastic exciton-scattering theory of optical lineshapes: Renormalized many-body contributions

Spectral line-shapes provide a window into the local environment coupled to a quantum transition in the condensed phase. In this paper, we build upon a stochastic model to account for non-stationary background processes produced by broad-band pulsed laser stimulation. In particular, we consider the contribution of pair-fluctuations arising from the full bosonic many-body Hamiltonian within a mean-field approximation, treating the coupling to the system as a stochastic noise term. Using the It{\^o} transformation, we consider two limiting cases for our model which lead to a connection between the observed spectral fluctuations and the spectral density of the environment. In the first case, we consider a Brownian environment and show that this produces spectral dynamics that relax to form dressed excitonic states and recover an Anderson-Kubo-like form for the spectral correlations. In the second case, we assume that the spectrum is Anderson-Kubo like, and invert to determine the corresponding background. Using the Jensen inequality, we obtain an upper limit for the spectral density of the background. The results presented here provide the technical tools for applying the stochastic model to a broad range of problems.

π-Stacking-Dependent Vibronic Couplings Drive Excited-State Dynamics in Perylenediimide Assemblies

Vibronic coupling, the interplay of electronic and nuclear vibrational motion, is considered a critical mechanism in photoinduced reactions such as energy transfer, charge transfer, and singlet fission. However, our understanding of how particular vibronic couplings impact excited-state dynamics is lacking due to the limited number of experimental studies of model molecular systems. Herein, we use two-dimensional electronic spectroscopy (2DES) to launch and interrogate a range of vibronic coherences in two distinct types of perylenediimide slip stacks─along the short and long molecular axes, which form either an excimer or a mixed state between the Frenkel exciton (FE) and charge transfer states. We explore the functionality of these vibronic coherences using quantum beatmaps, which display the Fourier amplitude signal oscillations as a function of pump and probe frequencies, along with knowledge of the characteristic signatures of the FE, ionic, and excimer species. We find that a low-frequency vibrational mode of the short-axis slip stack appears concomitantly with the formation of the excimer state, survives 2-fold longer than in the FE state in the reference monomer, and shows a phase shift compared to other modes. For the long-axis slip stacks, a pair of low-frequency modes coupled to a high-frequency coordinate of the FE state were found to play a critical role in mixed-state generation. Our findings thus experimentally reveal the complex and varying roles of vibronic couplings in tightly packed multimers undergoing a range of photoinduced processes.

Optical Cavity Manipulation and Nonlinear UV Molecular Spectroscopy of Conical Intersections in Pyrazine

Optical cavities provide a versatile platform for manipulating the excited-state dynamics of molecules via strong light-matter coupling. We employ optical absorption and two-multidimensional electronic spectroscopy simulations to investigate the effect of optical cavity coupling in the nonadiabatic dynamics of photoexcited pyrazine. We observe the emergence of a novel polaritonic conical intersection (PCI) between the electronic dark state and photonic surfaces as the cavity frequency is tuned. The PCI could significantly change the nonadiabatic dynamics of pyrazine by doubling the decay rate constant of the S₂ state population. Moreover, the absorption spectrum and excited-state dynamics could be systematically manipulated by tuning the strong light-matter interaction, e.g., the cavity frequency and cavity coupling strength. We propose that a tunable optical cavity-molecule system may provide promising approaches for manipulating the photophysical properties of molecules.

Dynamic effect of polymers at the surfactant-water interface: an ultrafast study

Interfaces play a role in controlling the rates and outcomes of chemical processes. Characterizing the interactions at heterogeneous interfaces is critical to developing a comprehensive model of the role of interfaces and confinement in modulating chemical reactions. Reverse micelles are an ideal model system for exploring the effect of encapsulated species on interfacial environments. Here, we use a combination of ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics (MD) simulations to characterize the picosecond interfacial dynamics in reverse micelles (RMs) containing acrylamide monomers and polyacrylamide polymers within the aqueous phase. The ester carbonyl vibrations of the sorbitan monostearate surfactants are examined to extract interfacial hydrogen-bonding populations and dynamics. Hydrogen bond populations at the ester carbonyl positions remain unchanged with the inclusion of either polymer or monomer species. Hydrogen-bond dynamics are not altered with the addition of monomer but are slowed down twofold in the presence of encapsulated polyacrylamide polymer species as a result of polymer chains partially localizing to the interface. These findings imply that kinetics of reactions that occur at interfaces or in confined environments could be modulated by interfacial localization of the different components.

Transient measurement of phononic states with covariance-based stochastic spectroscopy

We present a novel approach to transient Raman spectroscopy, which combines stochastic probe pulses and a covariance-based detection to measure stimulated Raman signals in alpha-quartz. A coherent broadband pump is used to simultaneously impulsively excite a range of different phonon modes, and the phase, amplitude, and energy of each mode are independently recovered as a function of the pump-probe delay by a noisy-probe and covariance-based analysis. Our experimental results and the associated theoretical description demonstrate the feasibility of 2D-Raman experiments based on the stochastic-probe schemes, with new capabilities not available in equivalent mean-value-based 2D-Raman techniques. This work unlocks the gate for nonlinear spectroscopies to capitalize on the information hidden within the noise and overlooked by a mean-value analysis.

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.

Yb-Doped Fiber Chirped Pulse Amplification System Delivering 1 mJ, 231 fs at 1 kHz Repetition Rate

In this paper a single-channel chirped pulse amplification laser system based on Yb-doped photonic crystal fiber was constructed, which achieved a pulse energy output of 1 mJ with a beam quality close to the diffraction limit. Pulsed synchronous pumping was used to suppress amplified spontaneous emission at a repetition rate of 1 kHz. The de-chirped pulse width of 231 fs was achieved by precise systematic dispersion control, and the corresponding peak power reached 3.85 GW.