Laser-Limited Signatures of Quantum Coherence

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.

Multimode vibrational dynamics and orientational effects in fluorescence-encoded infrared spectroscopy. II. Analysis of early-time signals

Developing fluorescence-encoded infrared (FEIR) vibrational spectroscopy for single-molecule applications requires a detailed understanding of how the molecular response and external experimental parameters manifest in the detected signals. In Paper I [L. Whaley-Mayda, A. Guha, and A. Tokmakoff, J. Chem. Phys. 159, 194201 (2023)] we introduced a nonlinear response function theory to describe vibrational dynamics, vibronic coupling, and transition dipole orientation in FEIR experiments with ultrashort pulses. In this second paper, we apply the theory to investigate the role of intermode vibrational coherence, the orientation of vibrational and electronic transition dipoles, and the effects of finite pulse durations in experimental measurements. We focus on measurements at early encoding delays-where signal sizes are largest and therefore of most value for single-molecule experiments, but where many of these phenomena are most pronounced and can complicate the appearance of data. We compare experiments on coumarin dyes with finite-pulse response function simulations to explain the time-dependent behavior of FEIR spectra. The role of the orientational response is explored by analyzing polarization-dependent experiments and their ability to resolve relative dipole angles in the molecular frame. This work serves to demonstrate the molecular information content of FEIR experiments, and develop insight and guidelines for their interpretation.

The influence of a Hamiltonian vibration vs a bath vibration on the 2D electronic spectra of a homodimer

We elucidate the influence of the system–bath boundary placement within an open quantum system, with emphasis on the two-dimensional electronic spectra, through the application of the hierarchical equations of motion formalism for an exciton system. We apply two different models, the Hamiltonian vibration model (HVM) and bath vibration model (BVM), to a monomer and a homodimer. In the HVM, we specifically include the vibronic states in the Hamiltonian capturing vibronic quenching, whereas in the BVM, all vibrational details are contained within the bath and described by an underdamped spectral density. The resultant spectra are analyzed in terms of energetic peak position and thermodynamic broadening precision in order to evaluate the efficacy of the two models. The HVM produces 2D spectra with accurate peak positional information, while the BVM is well suited to modeling dynamic peak broadening. For the monomer, both models produce equivalent spectra in the limit where additional damping associated with the underdamped vibration in the BVM approaches zero. This is supported by analytical results. However, for the homodimer, the BVM spectra are redshifted with respect to the HVM due to an absence of vibronic quenching in the BVM. The computational efficiency of the two models is also discussed in order to inform us of the most appropriate use of each method.

Modeling and Simulating the Excited-State Dynamics of a System with Condensed Phases: A Machine Learning Approach

Simulating the irreversible quantum dynamics of exciton- and electron-transfer problems poses a nontrivial challenge. Because the irreversibility of the system dynamics is a result of quantum thermal activation and dissipation caused by the surrounding environment, it is necessary to include infinite environmental degrees of freedom in the simulation. Because the capabilities of full quantum dynamics simulations that include the surrounding molecular degrees of freedom are limited, employing a system-bath model is a practical approach. In such a model, the dynamics of excitons or electrons are described by a system Hamiltonian, while the other degrees of freedom that arise from the environmental molecules are described by a harmonic oscillator bath (HOB) and system-bath interaction parameters. By extending on a previous study of molecular liquids [ J. Chem. Theory Comput. 2020, 16, 2099], here, we construct a system-bath model for exciton- and electron-transfer problems by means of a machine learning approach. We determine both the system and system-bath interaction parameters, including the spectral distribution of the bath, using the electronic excitation energies obtained from a quantum mechanics/molecular mechanics (QM/MM) simulation that is conducted as a function of time. Using the analytical expressions of optical response functions, we calculate linear and two-dimensional electronic spectra (2DES) for indocarbocyanine dimers in methanol. From these results, we demonstrate the capability of our approach to elucidate the nonequilibrium exciton dynamics of a quantum system in a nonintuitive manner.

Signatures of exciton dynamics and interaction in coherently and fluorescence-detected four- and six-wave-mixing two-dimensional electronic spectroscopy

Two-dimensional electronic spectroscopy (2DES) can be realized in increasing nonlinear orders of interaction with the electric field, bringing new information about single- and multi-particle properties and dynamics. Furthermore, signals can be detected both coherently (C-2DES) and by fluorescence (F-2DES), with fundamental and practical differences. We directly compare the simultaneous measurements of four- and six-wave mixing C-2DES and F-2DES on an excitonic heterodimer of squaraine molecules. Spectral features are described in increasing orders of nonlinearity by an explicit excitonic model. We demonstrate that the four-wave-mixing spectra are sensitive to one-exciton energies, their delocalization and dynamics, while the six-wave-mixing spectra include information on bi-exciton and higher excited states including the state energies, electronic coupling, and exciton-exciton annihilation. We focus on the possibility to extract the dynamics arising from exciton-exciton interaction directly from the six-wave-mixing spectra. To this end, in analogy to previously demonstrated fifth-order coherently detected exciton-exciton-interaction 2DES (EEI2D spectroscopy), we introduce a sixth-order fluorescence-detected EEI2D spectroscopy variant.

Room temperature quantum coherence vs. electron transfer in a rhodanine derivative chromophore

Understanding electron transfer in organic molecules is of great interest in quantum materials for light harvesting, energy conversion and integration of molecules into solar cells. This, however, poses the challenge of designing specific optimal molecular structure for which the processes of ultrafast quantum coherence and electron transport are not so well understood. In this work, we investigate subpicosecond time scale quantum dynamics and electron transfer in an efficient electron acceptor rhodanine chromophoric complex. We consider an open quantum system approach to model the complex-solvent interaction, and compute the crossover from weak to strong dissipation on the reduced system dynamics for both a polar (methanol) and a non polar solvent (toluene). We show that the electron transfer rates are enhanced in the strong chromophore-solvent coupling regime, being the highest transfer rates those found at room temperature. Even though the computed dynamics are highly non-Markovian, and they may exhibit a quantum character up to hundreds of femtoseconds, we show that quantum coherence does not necessarily optimise the electron transfer in the chromophore.

Correction of spectral distortion in two-dimensional electronic spectroscopy arising from the wedge-based delay line

Unlike the probe wavelength, which is spectrally resolved by monochromator, the excitation wavelength in two-dimensional electronic spectroscopy is retrieved by means of Fourier transform of the interference signal introduced by the coherence delay time between the first and second excitation laser pulses. Hence, the calibration of delay lines would determine its accuracy. In this work, we showed that an inaccurate calibration factor of wedge-based delay line would result in a global peak shift and asymmetric spectral twists along the excitation axis. Both theoretical analysis and experiments have shown that such spectral distortions can be corrected by an accurately predetermined calibration factor. The relative accuracy of calibration factor reaches 3 × 10^-5 in our setup. The dispersion effect of wedges also has been considered for the broadband excitation.

Characterization of Vibrational Coherence in Monomeric Bacteriochlorophyll a by Two-Dimensional Electronic Spectroscopy

Bacteriochlorophyll a (BChla) is the most abundant pigment found in the Bacterial Reaction Center (BRC) and light-harvesting proteins of photosynthetic purple and green bacteria. Recent two-dimensional electronic spectroscopy (2DES) studies of photosynthetic pigment-protein complexes including the BRC and the Fenna-Matthews-Olson (FMO) complex have shown oscillatory signals, or coherences, whose physical origin has been hotly debated. To better understand the observations of coherence in larger photosynthetic systems, it is important to carefully characterize the spectroscopic signatures of the monomeric pigments. Prior spectroscopic studies of BChla have differed significantly in their observations, with some studies reporting little to no coherence. Here we present evidence of strong coherences in monomeric BChla in isopropanol using 2DES at 77 K. We resolve many modes with frequencies that correspond well with known vibrational modes. We confirm their vibrational origin by comparing the 2D spectroscopic signatures with expectations based on a purely vibrational model.

Strongly coupled bacteriochlorin dyad studied using phase-modulated fluorescence-detected two-dimensional electronic spectroscopy

Fluorescence-detected two-dimensional electronic spectroscopy (F-2DES) projects the third-order non-linear polarization in a system as an excited electronic state population which is incoherently detected as fluorescence. Multiple variants of F-2DES have been developed. Here, we report phase-modulated F-2DES measurements on a strongly coupled symmetric bacteriochlorin dyad, a relevant 'toy' model for photosynthetic energy and charge transfer. Coherence map analysis shows that the strongest frequency observed in the dyad is well-separated from the excited state electronic energy gap, and is consistent with a vibrational frequency readily observed in bacteriochlorin monomers. Kinetic rate maps show a picosecond relaxation timescale between the excited states of the dyad. To our knowledge this is the first demonstration of coherence and kinetic analysis using the phase-modulation approach to F-2DES.

Spectral Filtering as a Tool for Two-Dimensional Spectroscopy: A Theoretical Model

Two-dimensional optical spectroscopy is a powerful technique for the probing of coherent quantum superpositions. Recently, the finite width of the laser spectrum has been employed to selectively tune experiments for the study of particular coherences. This involves the exclusion of certain transition frequencies, which results in the elimination of specific Liouville pathways. The rigorous analysis of such experiments requires the use of ever more sophisticated theoretical models for the optical spectroscopy of electronic and vibronic systems. Here we develop a nonimpulsive and non-Markovian model, which combines an explicit definition of the laser spectrum, via the equation of motion-phase matching approach (EOM-PMA), with the hierarchical equations of motion (HEOM). This theoretical framework is capable of simulating the 2D spectroscopy of vibronic systems with low frequency modes, coupled to environments of intermediate and slower time scales. In order to demonstrate the spectral filtering of vibronic coherences, we examine the elimination of lower energy peaks from the 2D spectra of a zinc porphyrin monomer upon blue-shifting the laser spectrum. The filtering of Liouville pathways is revealed through the disappearance of peaks from the amplitude spectra for a coupled vibrational mode.

Spectral line shapes in linear absorption and two-dimensional spectroscopy with skewed frequency distributions

The effect of Gaussian dynamics on the line shapes in linear absorption and two-dimensional correlation spectroscopy is well understood as the second-order cumulant expansion provides exact spectra. Gaussian solvent dynamics can be well analyzed using slope line analysis of two-dimensional correlation spectra as a function of the waiting time between pump and probe fields. Non-Gaussian effects are not as well understood, even though these effects are common in nature. The interpretation of the spectra, thus far, relies on complex case to case analysis. We investigate spectra resulting from two physical mechanisms for non-Gaussian dynamics, one relying on the anharmonicity of the bath and the other on non-linear couplings between bath coordinates. These results are compared with outcomes from a simpler log-normal dynamics model. We find that the skewed spectral line shapes in all cases can be analyzed in terms of the log-normal model, with a minimal number of free parameters. The effect of log-normal dynamics on the spectral line shapes is analyzed in terms of frequency correlation functions, maxline slope analysis, and anti-diagonal linewidths. A triangular line shape is a telltale signature of the skewness induced by log-normal dynamics. We find that maxline slope analysis, as for Gaussian dynamics, is a good measure of the solvent dynamics for log-normal dynamics.

Two-Dimensional Spectroscopy Is Being Used to Address Core Scientific Questions in Biology and Materials Science

Two-dimensional spectroscopy is a powerful tool for extracting structural and dynamic information from a wide range of chemical systems. We provide a brief overview of the ways in which two-dimensional visible and infrared spectroscopies are being applied to elucidate fundamental details of important processes in biological and materials science. The topics covered include amyloid proteins, photosynthetic complexes, ion channels, photovoltaics, batteries, as well as a variety of promising new methods in two-dimensional spectroscopy.

The origin of absorptive features in the two-dimensional electronic spectra of rhodopsin

A three-state three-mode model Hamiltonian reveals the origin of the absorptive features in the two-dimensional electronic spectra of rhodopsin.

The Primary Photochemistry of Vision Occurs at the Molecular Speed Limit

Ultrafast photochemical reactions are initiated by vibronic transitions from the reactant ground state to the excited potential energy surface, directly populating excited-state vibrational modes. The primary photochemical reaction of vision, the isomerization of retinal in the protein rhodopsin, is known to be a vibrationally coherent reaction, but the Franck-Condon factors responsible for initiating the process have been difficult to resolve with conventional time-resolved spectroscopies. Here we employ experimental and theoretical 2D photon echo spectroscopy to directly resolve for the first time the Franck-Condon factors that initiate isomerization on the excited potential energy surface and track the reaction dynamics. The spectral dynamics reveal vibrationally coherent isomerization occurring on the fastest possible time scale, that of a single period of the local torsional reaction coordinate. We successfully model this process as coherent wavepacket motion through a conical intersection on a ∼30 fs time scale, confirming the reaction coordinate as a local torsional coordinate with a frequency of ∼570 cm^-1. As a result of spectral features being spread out along two frequency coordinates, we unambiguously assign reactant and product states following passage through the conical intersection, which reveal the key vibronic transitions that initiate the vibrationally coherent photochemistry of vision.

Resolving Vibrational from Electronic Coherences in Two-Dimensional Electronic Spectroscopy: The Role of the Laser Spectrum

The observation of coherent quantum effects in photosynthetic light-harvesting complexes prompted the question whether quantum coherence could be exploited to improve the efficiency in new energy materials. The detailed characterization of coherent effects relies on sensitive methods such as two-dimensional electronic spectroscopy (2D-ES). However, the interpretation of the results produced by 2D-ES is challenging due to the many possible couplings present in complex molecular structures. In this work, we demonstrate how the laser spectral profile can induce electronic coherencelike signals in monomeric chromophores, potentially leading to data misinterpretation. We argue that the laser spectrum acts as a filter for certain coherence pathways and thus propose a general method to differentiate vibrational from electronic coherences.