An official website of the United States government
Accurately modeling absorption and fluorescence spectra for molecules in solution poses a challenge due to the need to incorporate both vibronic and environmental effects, as well as the necessity of accurate excited state electronic structure calculations. Nuclear ensemble approaches capture explicit environmental effects, Franck–Condon methods capture vibronic effects, and recently introduced ensemble-Franck–Condon approaches combine the advantages of both methods. In this study, we present and analyze simulated absorption and fluorescence spectra generated with combined ensemble-Franck–Condon approaches for three chromophore–solvent systems and compare them to standard ensemble and Franck–Condon spectra, as well as to the experiment. Employing configurations obtained from ground and excited state ab initio molecular dynamics, three combined ensemble-Franck–Condon approaches are directly compared to each other to assess the accuracy and relative computational time. We find that the approach employing an average finite-temperature Franck–Condon line shape generates spectra nearly identical to the direct summation of an ensemble of Franck–Condon spectra at one-fourth of the computational cost. We analyze how the spectral simulation method, as well as the level of electronic structure theory, affects spectral line shapes and associated Stokes shifts for 7-nitrobenz-2-oxa-1,3-diazol-4-yl and Nile red in dimethyl sulfoxide and 7-methoxy coumarin-4-acetic acid in methanol. For the first time, our studies show the capability of combined ensemble-Franck–Condon methods for both absorption and fluorescence spectroscopy and provide a powerful tool for simulating linear optical spectra.
more »
« less
This content will become publicly available on January 28, 2026
Time-resolved vibronic spectra with nuclear–electronic orbital time-dependent configuration interaction
Time-resolved spectroscopy is an important tool for probing photochemically induced nonequilibrium dynamics and energy transfer. Herein, a method is developed for the ab initio simulation of vibronic spectra and dynamical processes. This framework utilizes the recently developed nuclear–electronic orbital time-dependent configuration interaction (NEO-TDCI) approach, which treats all electrons and specified nuclei quantum mechanically on the same footing. A strategy is presented for calculating time-resolved vibrational and electronic absorption spectra from any initial condition. Although this strategy is general for any TDCI implementation, utilizing the NEO framework allows for the explicit inclusion of quantized nuclei, as illustrated through the calculation of vibrationally hot spectra. Time-resolved spectra produced by either vibrational or electronic Rabi oscillations capture ground-state absorption, stimulated emission, and excited-state absorption between vibronic states. This methodology provides the foundation for fully ab initio simulations of multidimensional spectroscopic experiments.
more »
« less
- PAR ID:
- 10608831
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 4
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Bacteriochlorophyll a (Bchl a) and chlorophyll a (Chl a) play important roles as light absorbers in photosynthetic antennae and participate in the initial charge-separation steps in photosynthetic reaction centers. Despite decades of study, questions remain about the interplay of electronic and vibrational states within the Q-band and its effect on the photoexcited dynamics. Here we report results of polarized two-dimensional electronic spectroscopic measurements, performed on penta-coordinated Bchl a and Chl a and their interpretation based on state-of-the-art time-dependent density functional theory calculations and vibrational mode analysis for spectral shapes. We find that the Q-band of Bchl a is comprised of two independent bands, that are assigned following the Gouterman model to Q x and Q y states with orthogonal transition dipole moments. However, we measure the angle to be ∼75°, a finding that is confirmed by ab initio calculations. The internal conversion rate constant from Q x to Q y is found to be 11 ps −1 . Unlike Bchl a, the Q-band of Chl a contains three distinct peaks with different polarizations. Ab initio calculations trace these features back to a spectral overlap between two electronic transitions and their vibrational replicas. The smaller energy gap and the mixing of vibronic states result in faster internal conversion rate constants of 38–50 ps −1 . We analyze the spectra of penta-coordinated Bchl a and Chl a to highlight the interplay between low-lying vibronic states and their relationship to photoinduced relaxation. Our findings shed new light on the photoexcited dynamics in photosynthetic systems where these chromophores are primary pigments.more » « less
-
Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born–Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.more » « less
-
The four lowest Ω substates (X2Π3/2,g, X2Π1/2,g, A2Π3/2,u and A2Π1/2,u) of the cation have been studied by high-precision ab initio calculations in comparison with experimental high-resolution absorption spectra. The potential energy curves were calculated using the multi-reference configuration interaction (MRCI) method and Dirac method, respectively. Rovibrational levels of these electronic states were derived by solving the radial Schrödinger rovibrational equation. Molecular constants were obtained in fitting energy levels to a spectroscopic model. Using the fit spectroscopic constants and newly calculated transition dipole moment matrix elements, line strengths of vibronic bands in the A2Π3/2,u- X2Π3/2,g system, as well as Einstein A coefficients for 45 of these bands with ν′ = 11–19 and ν′′ = 1–5, have been derived. The Einstein A coefficients were used to compute radiative lifetimes of the ν′ = 11–19 vibrational levels of the A2Π3/2,u state. Enhancement factors for detecting the variation of the fine-structure constant (α) and the proton-to-electron mass ratio(µ) using transitions between nearly degenerate rovibronic levels of these low-lying states have been calculated.more » « less
-
Attosecond extreme ultraviolet (XUV) and soft x-ray sources provide powerful new tools for studying ultrafast molecular dynamics with atomic, state, and charge specificity. In this report, we employ attosecond transient absorption spectroscopy (ATAS) to follow strong-field-initiated dynamics in vinyl bromide. Probing the Br M edge allows one to assess the competing processes in neutral and ionized molecular species. Using ab initio non-adiabatic molecular dynamics, we simulate the neutral and cationic dynamics resulting from the interaction of the molecule with the strong field. Based on the dynamics results, the corresponding time-dependent XUV transient absorption spectra are calculated by applying high-level multi-reference methods. The state-resolved analysis obtained through the simulated dynamics and related spectral contributions enables a detailed and quantitative comparison with the experimental data. The main outcome of the interaction with the strong field is unambiguously the population of the first three cationic states, D1, D2, and D3. The first two show exclusively vibrational dynamics while the D3 state is characterized by an ultrafast dissociation of the molecule via C–Br bond rupture within 100 fs in 50% of the analyzed trajectories. The combination of the three simulated ionic transient absorption spectra is in excellent agreement with the experimental results. This work establishes ATAS in combination with high-level multi-reference simulations as a spectroscopic technique capable of resolving coupled non-adiabatic electronic-nuclear dynamics in photoexcited molecules with sub-femtosecond resolution.more » « less
