Temperature effects on singlet fission dynamics mediated by a conical intersection

Cavity-Tuned Exciton Dynamics in Transition Metal Dichalcogenides Monolayers

A fully quantum, numerically accurate methodology is presented for the simulation of the exciton dynamics and time-resolved fluorescence of cavity-tuned two-dimensional (2D) materials at finite temperatures. This approach was specifically applied to a monolayer WSe2 system. Our methodology enabled us to identify the dynamical and spectroscopic signatures of polaronic and polaritonic effects and to elucidate their characteristic timescales across a range of exciton-cavity couplings. The approach employed can be extended to simulation of various cavity-tuned 2D materials, specifically for exploring finite temperature nonlinear spectroscopic signals.

Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study

This study leverages two-pulse femtosecond stimulated Raman spectroscopy (2FSRS) to characterize molecular systems with avoided crossings (ACs) and conical intersections (CIs) in their low-lying excited electronic states. By simulating 2FSRS spectra of microscopically inspired ACs and CIs models, we demonstrate that 2FSRS not only delivers valuable information on the molecular parameters characterizing ACs and CIs but also helps distinguish between these two systems.

Finite-Temperature Hole-Magnon Dynamics in an Antiferromagnet

Employing the numerically accurate multiple Davydov Ansatz in combination with the thermo-field dynamics approach, we delve into the interplay of the finite-temperature dynamics of holes and magnons in an antiferromagnet, which allows for scrutinizing previous predictions from the self-consistent Born approximation while offering, for the first time, accurate finite-temperature computation of detailed magnon dynamics as a response and a facilitator to the hole motion. The study also uncovers a pronounced temperature dependence of the magnon and hole populations, pointing to the feasibility of potential thermal manipulation and control of hole dynamics. Our methodology can be applied not only to the calculation of steady-state angular-resolved photoemission spectra but also to the simulation of femtosecond terahertz pump-probe and other nonlinear signals for the characterization of antiferromagnetic materials.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

Exciton Dynamics and Time-Resolved Fluorescence in Nanocavity-Integrated Monolayers of Transition-Metal Dichalcogenides

We have developed an ab initio-based, fully quantum, numerically accurate methodology for the simulation of the exciton dynamics and time- and frequency-resolved fluorescence spectra of the cavity-controlled two-dimensional materials at finite temperatures and applied this methodology to the single-layer WSe₂ system. Specifically, the multiple Davydov D₂ Ansatz has been employed in combination with the method of thermofield dynamics for the finite-temperature extension of accurate time-dependent variation. This allowed us to establish dynamical and spectroscopic signatures of the polaronic and polaritonic effects as well as uncover their characteristic time scales in the relevant range of temperatures. Our study reveals the pivotal role of multidimensional conical intersections in controlling the many-body dynamics of highly intertwined excitonic, phononic, and photonic modes.

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.

Plenty of Room on the Top: Pathways and Spectroscopic Signatures of Singlet Fission from Upper Singlet States

We investigate dynamic signatures of the singlet fission (SF) process triggered by the excitation of a molecular system to an upper singlet state S_N (N > 1) and develop a computational methodology for the simulation of nonlinear spectroscopic signals revealing the S_N → TT₁ SF in real time. We demonstrate that SF can proceed directly from the upper state S_N, bypassing the lowest excited state, S₁. We determine the main S_N → TT₁ reaction pathways and show by computer simulation and spectroscopic measurements that the S_N-initiated SF can be faster and more efficient than the traditionally studied S₁ → TT₁ SF. We claim that the S_N → TT₁ SF offers novel promising opportunities for engineering SF systems and enhancing SF yields.

Singlet Fission, Polaron Generation and Intersystem Crossing in Hexaphenyl Film

The ultrafast dynamics of triplet excitons and polarons in hexaphenyl film was investigated by time-resolved fluorescence and femtosecond transient absorption techniques under various excitation photon energies. Two distinct pathways of triplet formation were clearly observed. Long-lived triplet states are populated within 4.5 ps via singlet fission-intersystem crossing, while the short-lived triplet states (1.5 ns) are generated via singlet fission from vibrational electronic states. In the meantime, polarons were formed from hot excitons on a timescale of <30 fs and recombined in ultrafast lifetime (0.37 ps). In addition, the characterization of hexaphenyl film suggests the morphologies of crystal and aggregate to wide applications in organic electronic devices. The present study provides a universally applicable film fabrication in hexaphenyl system towards future singlet fission-based solar cells.

Accurate Simulation of Spectroscopic Signatures of Cavity-Assisted, Conical-Intersection-Controlled Singlet Fission Processes

A numerically accurate, fully quantum methodology has been developed for the simulation of the dynamics and nonlinear spectroscopic signals of cavity-assisted, conical-intersection-controlled singlet fission systems. The methodology is capable of handling several molecular systems strongly coupled to the photonic mode of the cavity and treats the intrinsic conical intersection and cavity-induced polaritonic conical intersections in a numerically exact manner. Contributions of higher-lying molecular electronic states are accounted for comprehensively. The intriguing process of cavity-modified fission dynamics, including all of its electronic, vibrational, and photonic degrees of freedom, together with its two-dimensional spectroscopic manifestation, is simulated for two rubrene dimers strongly coupled to the cavity mode.

Engineering Cavity Singlet Fission in Rubrene

By employing the numerically exact multiple Davydov D₂ ansatz, we study cavity-manipulated singlet fission that is mediated by polaritonic conical intersections for both one- and two-molecule systems. The population evolution of the TT state and the cavity photons is carefully examined in search for a high fission efficiency via cavity engineering. Several interesting mechanisms have been uncovered, such as photon-assisted singlet fission, system localization via a displaced photon state, and collective enhancement of the fission efficiency for the two-molecule system. It is also found that the system localization process in the two-molecule system differs substantially from that in the one-molecule system because of the appearance of a novel central polaritonic conical intersection in the two-molecule system. It has been demonstrated that the cavity-controlled singlet fission process can be switched on and off by controlling the average pumping photon number.

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.

Variational Squeezed Davydov Ansatz for Realistic Chemical Systems with Nonlinear Vibronic Coupling

Chemical systems normally possess strong nonlinear vibronic couplings at both zero and finite temperature. For the lowest-order quadratic couplings, here, we introduce a squeezing operator into a variational coherent-state-based method, Davydov ansatz, to simulate the quantum dynamics and the respective spectroscopy. Two molecular systems, pyrazine and the 2-pyridone dimer, are taken as calculated model systems, both of which involve nontrivial quadratic vibronic couplings in high- and low-frequency regions, respectively. Upon a comparison with the benchmarks, the method manifests its advantage for nonlinear couplings. The squeezed bases are also proven to be applicable for the finite temperature by adapting with the thermofield dynamics.

Dynamics of disordered Tavis-Cummings and Holstein-Tavis-Cummings models

By employing the time-dependent variational principle and the versatile multi-D₂ Davydov trial states, in combination with the Green's function method, we study the dynamics of the Tavis-Cummings model and the Holstein-Tavis-Cummings model in the presence of diagonal disorder and cavity-qubit coupling disorder. For the Tavis-Cummings model, time evolution of the photon population, the optical absorption spectra, and the hetero-entanglement between the qubits and the cavity mode are calculated by using the Green's function method to corroborate numerically exact results of Davydov's Ansätze. For the Holstein-Tavis-Cummings model, only the latter is utilized to simulate effects of disorder on the photon population dynamics and the absorption spectra. We have demonstrated that the complementary employment of analytical and numerical methods permits uncovering a fairly comprehensive picture of a variety of complex behaviors in disordered multidimensional polaritonic cavity quantum electrodynamics systems.

Quantum dynamics simulation of intramolecular singlet fission in covalently linked tetracene dimer

In this work, we study singlet fission in tetracene para-dimers, covalently linked by a phenyl group. In contrast to most previous studies, we account for the full quantum dynamics of the combined excitonic and vibrational system. For our simulations, we choose a numerically unbiased representation of the molecule's wave function, enabling us to compare with experiments, exhibiting good agreement. Having access to the full wave function allows us to study in detail the post-quench dynamics of the excitons. Here, one of our main findings is the identification of a time scale t₀ ≈ 35 fs dominated by coherent dynamics. It is within this time scale that the larger fraction of the singlet fission yield is generated. We also report on a reduced number of phononic modes that play a crucial role in the energy transfer between excitonic and vibrational systems. Notably, the oscillation frequency of these modes coincides with the observed electronic coherence time t₀. We extend our investigations by also studying the dependency of the dynamics on the excitonic energy levels that, for instance, can be experimentally tuned by means of the solvent polarity. Here, our findings indicate that the singlet fission yield can be doubled, while the electronic coherence time t₀ is mainly unaffected.

Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.

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

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

Simulation of Time- and Frequency-Resolved Four-Wave-Mixing Signals at Finite Temperatures: A Thermo-Field Dynamics Approach

We propose a new approach to simulate four-wave-mixing signals of molecular systems at finite temperatures by combining the multiconfigurational Ehrenfest method with the thermo-field dynamics theory. In our approach, the four-time correlation functions at finite temperatures are mapped onto those at zero temperature in an enlarged Hilbert space with twice the vibrational degrees of freedom. As an illustration, we have simulated three multidimensional spectroscopic signals, time- and frequency-resolved fluorescence spectra, transient-absorption pump-probe spectra, and electronic two-dimensional (2D) spectra at finite temperatures, for a conical intersection-mediated singlet fission model of a rubrene crystal. It is shown that a detailed dynamical picture of the singlet fission process can be extracted from the three spectroscopic signals. An increasing temperature leads to lower intensities of the signals and broadened vibrational peaks, which can be attributed to faster singlet-triplet population transfer and stronger bath-induced electronic dephasing at higher temperatures.

Singlet fission dynamics and optical spectra of pentacene and its derivatives

A multimode Brownian oscillator model is employed to investigate the absorption spectra of pentacene and its derivatives in solution and thin films. Excellent agreement has been obtained between simulated and measured absorption spectra. Furthermore, using parameters obtained from fitting the absorption spectra of these pentacene derivatives, the singlet fission dynamics and two-dimensional electronic spectra of an ab initio Hamiltonian are investigated by Dirac-Frenkel time-dependent variation with multiple Davydov trial states. It is found that the periodic wave packet motion induced in the displaced excited state, and the accompanying vibrational relaxation, can be visualized by two-dimensional electronic spectra at short times.