Origin of Unexpectedly Simple Oscillatory Responses in the Excited-State Dynamics of Disordered Molecular Aggregates

2D Electronic Spectroscopy Uncovers 2D Materials: Theoretical Study of Nanocavity-Integrated Monolayer Semiconductors

Transition metal dichalcogenides (TMDs) have emerged as promising 2D semiconductors due to their strong excitonic effects, spin-valley coupling, and tunable light-matter interactions. Here, we employ a fully quantum, numerically "exact" multi-Davydov Ansatz approach to simulate two-dimensional electronic spectroscopy signals in hBN-encapsulated WSe2 monolayers integrated with a tunable nanocavity. By incorporating both momentum-bright and momentum-dark excitons alongside detailed phonon dispersion, our model captures vibrational resonances and exciton-polariton behaviors, enabling the evaluation of beating maps (3D spectra) that disentangle ground-state bleach and stimulated emission pathways. The results highlight the essential role of vibronic coherence in TMD monolayers and offer quantitative guidance for the design of next-generation optoelectronic devices based on cavity-coupled 2D materials.

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

TD-DMRG Study of Exciton Dynamics with both Thermal and Static Disorders for Fenna-Matthews-Olson Complex

Photosynthesis is a fundamental process that converts solar energy into chemical energy. Understanding the microscopic mechanisms of energy transfer in photosynthetic systems is crucial for the development of novel optoelectronic materials. Simulating these processes poses significant challenges due to the intricate interactions between electrons and phonons, compounded by static disorder. In this work, we present a numerically nearly exact study using the time-dependent density matrix renormalization group (TD-DMRG) method to simulate the quantum dynamics of the Fenna-Matthews-Olson (FMO) complex considering an eight-site model with both thermal and static disorders. We employ the thermo-field dynamics formalism for temperature effects. We merge all electronic interactions into one large matrix product state (MPS) site, boosting accuracy efficiently without increasing complexity. Previous combined experimental and computational studies indicated that the static disorders range from 30 to 90 cm-1 for different FMO sites. We employ a Gaussian distribution and the auxiliary bosonic operator approach to consider the static disorder in our TD-DMRG algorithm. We investigate the impact of different initial excitation sites, temperatures, and degrees of static disorder on the exciton dynamics and temporal coherence. It is found that under the influence of the experimentally determined static disorder strength, the exciton population evolution shows a non-negligible difference at zero temperature, while it is hardly affected at room temperature.

Finite temperature dynamics in a polarized sub-Ohmic heat bath: A hierarchical equations of motion-tensor train study

The dynamics of the sub-Ohmic spin-boson model under polarized initial conditions at finite temperatures is investigated by employing both analytical tools and the numerically accurate hierarchical equations of motion-tensor train method. By analyzing the features of nonequilibrium dynamics, we discovered a bifurcation phenomenon, which separates two regimes of the dynamics. It is found that before the bifurcation time, increasing temperature slows down the population dynamics, while the opposite effect occurs after the bifurcation time. The dynamics is highly sensitive to both initial preparation of the bath and thermal effects.

Finite temperature dynamics of the Holstein-Tavis-Cummings model

By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.

Two-dimensional fluorescence excitation spectroscopy: A novel technique for monitoring excited-state photophysics of molecular species with high time and frequency resolution

We propose a novel UV/Vis femtosecond spectroscopic technique, two-dimensional fluorescence-excitation (2D-FLEX) spectroscopy, which combines spectral resolution during the excitation process with exclusive monitoring of the excited-state system dynamics at high time and frequency resolution. We discuss the experimental feasibility and realizability of 2D-FLEX, develop the necessary theoretical framework, and demonstrate the high information content of this technique by simulating the 2D-FLEX spectra of a model four-level system and the Fenna-Matthews-Olson antenna complex. We show that the evolution of 2D-FLEX spectra with population time directly monitors energy transfer dynamics and can thus yield direct qualitative insight into the investigated system. This makes 2D-FLEX a highly efficient instrument for real-time monitoring of photophysical processes in polyatomic molecules and molecular aggregates.

Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore

By employing the doorway-window (DW) on-the-fly simulation protocol, we performed ab initio simulations of peak evolutions and beating maps of electronic two-dimensional (2D) spectra of a polyatomic molecule in the gas phase. As the system under study, we chose pyrazine, which is a paradigmatic example of photodynamics dominated by conical intersections (CIs). From the technical perspective, we demonstrate that the DW protocol is a numerically efficient methodology suitable for simulations of 2D spectra for a wide range of excitation/detection frequencies and population times. From the information content perspective, we show that peak evolutions and beating maps not only reveal timescales of transitions through CIs but also pinpoint the most relevant coupling and tuning modes active at these CIs.

Dynamics of diverse coherences in primary charge separation of bacterial reaction center at 77 K revealed by wavelet analysis

To uncover the mechanism behind the high photo-electronic conversion efficiency in natural photosynthetic complexes it is essential to trace the dynamics of electronic and vibrational quantum coherences. Here we apply wavelet analysis to two-dimensional electronic spectroscopy data for three purple bacterial reaction centers with mutations that produce drastically different rates of primary charge separation. From the frequency distribution and dynamic evolution features of the quantum beating, electronic coherence with a dephasing lifetime of ~50 fs, vibronic coherence with a lifetime of ~150 fs and vibrational/vibronic coherences with a lifetime of 450 fs are distinguished. We find that they are responsible for, or couple to, different specific steps during the primary charge separation process, i.e., intradimer charge transfer inside the special bacteriochlorophyll pair followed by its relaxation and stabilization of the charge-transfer state. The results enlighten our understanding of how quantum coherences participate in, and contribute to, a biological electron transfer reaction.

Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning

Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.

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.

Hierarchical Equations-of-Motion Method for Momentum System-Bath Coupling

For a broad class of quantum models of practical interest, we demonstrate that the Hamiltonian of the system nonlinearly coupled to a harmonic bath through the system and bath coordinates can be equivalently mapped into the Hamiltonian of the system bilinearly coupled to the bath through the system and bath momenta. We show that the Hamiltonian with bilinear system-bath momentum coupling can be treated by the hierarchical equations-of-motion (HEOM) method and present the corresponding proof-of-principle simulations. The developed methodology creates the opportunity to scrutinize a new family of nonlinear quantum systems by the numerically accurate HEOM method.

Electronic coherence lifetimes of the Fenna-Matthews-Olson complex and light harvesting complex II

The study of coherence between excitonic states in naturally occurring photosynthetic systems offers tantalizing prospects for uncovering mechanisms of efficient energy transport.

The study of coherence between excitonic states in naturally occurring photosynthetic systems offers tantalizing prospects of uncovering mechanisms of efficient energy transport. However, experimental evidence of functionally relevant coherences in wild-type proteins has been tentative, leading to uncertainty in their importance at physiological conditions. Here, we extract the electronic coherence lifetime and frequency using a signal subtraction procedure in two model pigment-protein-complexes (PPCs), light harvesting complex II (LH2) and the Fenna–Matthews–Olson complex (FMO), and find that the coherence lifetimes occur at the same timescale (<100 fs) as energy transport between states at the energy level difference equal to the coherence energy. The pigment monomer bacteriochlorophyll a (BChla) shows no electronic coherences, supporting our methodology of removing long-lived vibrational coherences that have obfuscated previous assignments. This correlation of timescales and energy between coherences and energy transport reestablishes the time and energy scales that quantum processes may play a role in energy transport.

Orientational relaxation of a quantum linear rotor in a dissipative environment: Simulations with the hierarchical equations-of-motion method

We study the effect of a dissipative environment on the orientational relaxation of a three-dimensional quantum linear rotor. We provide a derivation of the Hamiltonian of a linear rotor coupled to a harmonic bath from first principles, confirming earlier conjectures. The dynamics generated by this Hamiltonian is investigated by the hierarchical equations-of-motion method assuming a Drude spectral density of the bath. We perform numerically accurate simulations and analyze the behavior of orientational correlation functions and the rotational structures of infrared absorption and Raman scattering spectra. We explore the features of orientational correlation functions and their spectra for a wide range of system-bath couplings, bath memory times, and temperatures. We discuss the signatures of the orientational relaxation in the underdamped regime, the strongly damped regime, and the librational regime. We show that the behavior of orientational correlation functions and their spectra can conveniently be analyzed in terms of three characteristic times, which are explicitly expressed in terms of the parameters of the Hamiltonian.