An official website of the United States government
Two-dimensional electronic spectroscopy (2DES) is a powerful experimental technique, as it directly probes the nonlinear (third-order) response function of the system, providing key insights into ultrafast energy transfer and relaxation processes. However, 2DES experiments are generally difficult to interpret, often relying on simulations in order to associate observed spectral features with specific underlying system dynamics. For this reason, the development of robust, computationally inexpensive theoretical methods for modeling these experiments remains an active area of research. We have recently derived such an approach for computing the exact finite-temperature nonlinear response function for harmonic Hamiltonians within the Condon approximation, assuming that the transition dipole moment is independent of nuclear coordinates. In this work, we extend our formalism to exactly account for non-Condon/Herzberg−Teller (HT) type contributions to the nonlinear response function, which are known to be crucial for accurately describing linear optical spectra in a wide range of molecular systems. We highlight the key insights that can be gained from our new method, named FC2DES+HT, by simulating the 2DES signals of two molecules with known non-Condon behavior, the phenolate anion and free-base porphyrin. The results demonstrate that Herzberg−Teller couplings substantially impact energy relaxation dynamics in these systems.
more »
« less
This content will become publicly available on June 10, 2026
FC2DES: Modeling 2D Electronic Spectroscopy for Harmonic Hamiltonians
Two-dimensional electronic spectroscopy (2DES) can provide detailed insight into the energy transfer and relaxation dynamics of chromophores by directly measuring the nonlinear response function of the system. However, experiments are often difficult to interpret, and the development of computationally affordable approaches to simulate experimental signals is desirable. For linear spectroscopy, optical spectra of small to medium-sized molecules can be efficiently calculated in the Franck-Condon approach. Approximating the nuclear degrees of freedom as harmonic around the ground- and excited-state minima, closed-form expressions for the exact finite-temperature linear response function can be derived using known solutions for the propagation operator between normal mode coordinate sets, fully accounting for Duschinsky mode-mixing effects. In the present work, we demonstrate that a similar approach can be utilized to yield analogous closed-form expression for the finite-temperature nonlinear (third-order) response function of harmonic nuclear Hamiltonians. The resulting approach, named FC2DES, is implemented on graphics processing units (GPUs), allowing efficient computations of 2DES signals for medium-sized molecules containing hundreds of normal modes. Benchmark comparisons against the widely used cumulant method for computing 2DES signals are performed on small model systems, as well as the nile red molecule. We highlight the advantages of the FC2DES approach, especially in systems with moderate Duschinsky mode mixing and for long delay times in the nonlinear response function, where low-order cumulant approximations are shown to fail.
more »
« less
- Award ID(s):
- 2441876
- PAR ID:
- 10650858
- Publisher / Repository:
- American Chemical Society
- Date Published:
- Journal Name:
- Journal of Chemical Theory and Computation
- Volume:
- 21
- Issue:
- 11
- ISSN:
- 1549-9618
- Page Range / eLocation ID:
- 5625 to 5641
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.more » « less
-
In an ultrafast nonlinear optical interaction, the electric field of the emitted nonlinear signal provides direct access to the induced nonlinear transient polarization or transient currents and thus carries signatures of ultrafast dynamics in a medium. Measurement of the electric field of such signals offers sensitive observables to track ultrafast electron dynamics in various systems. In this work, we resolve the real-time phase of the electric field of a femtosecond third-order nonlinear optical signal in the molecular frame. The electric field emitted from impulsively pre-aligned gas-phase molecules at room temperature, in a degenerate four-wave mixing scheme, is measured using a spectral interferometry technique. The nonlinear signal is measured around a rotational revival to extract its molecular-frame angle dependence from pump-probe time-delay scans. By comparing these measurements for two linear molecules, carbon dioxide and nitrogen, we show that the measured second-order phase parameter (temporal chirp) of the signal is sensitive to the valence electronic symmetry of the molecules, whereas the amplitude of the signal does not show such sensitivity. We compare measurements to theoretical calculations of the chirp observable in the molecular frame. This work is an important step towards using electric field measurements in nonlinear optical spectroscopy to study ultrafast dynamics of electronically excited molecules in the molecular frame.more » « less
-
For Rayleigh-Bénard-Poiseuille flows, thermal stratification resulting from a wall-normal temperature gradient together with an opposing gravitational field can lead to buoyancy-driven instability. Moreover, for sufficiently large Reynolds numbers, viscosity-driven instability can occur. Two higher-order-accurate methods based on the full and linearized Navier-Stokes equations were developed for investigating the temporal stability of such flows. The new methods employ a spectral discretization in the homogeneous directions. In the wall-normal direction, the convective and viscous terms are discretized with fifth-order-accurate biased and fourth-order-accurate central compact finite differences. A fourth-order-accurate explicit Runge-Kutta method is employed for time integration. To validate the methods, the primary instability was investigated for different combinations of the Reynolds and Rayleigh number. The results from these primary stability investigations are consistent with linear stability theory results from the literature with respect to both the onset of the instability and the dependence of the temporal growth rate on the wave angle. For the cases with buoyancy-driven instability, strong linear growth is observed for a broad range of spanwise wavenumbers. The largest growth rates are obtained for a wave angle of 90deg. For the cases with viscosity-driven instability, the linear growth rates are lower and the first mode to experience nonlinear growth is a higher harmonic with half the wavelength of the fundamental.more » « less
-
Network tomography aims at estimating source-destination traffic rates from link traffic measurements. This inverse problem was formulated by Vardi in 1996 for independent Poisson traffic over networks operating under deterministic as well as random routing regimes. Vardi used a second-order moment matching approach to estimate the rates where a solution for the resulting linear matrix equation was obtained using an iterative minimum I-divergence procedure. Vardi’s second-order moment matching approach was recently extended to higher order cumulant matching approach with the goal of improving the rank of the system of linear equations. In this paper we go one step further and develop a moment generating function matching approach for rate estimation, and seek a least squares as well as an iterative minimum I-divergence solution of the resulting linear equations. We also specialize this approach to a characteristic function matching approach which exhibits some advantages. These follow from the fact that the characteristic function matching approach results in fewer conflicting equations involving the empirical estimates. We demonstrate that the new approach outperforms the cumulant matching approach while being conceptually simpler.more » « less
