Search

We combine on-the-fly trajectory surface hopping simulations and the doorway–window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We formulate a contraction theorem that maps quantum dynamics of a multilevel degenerate system (DS) driven by a time-dependent external field to the dynamics of the corresponding contracted non-degenerate system (CNS) of lower dimension, provided transitions between each pair of degenerate levels in the DS have identical transition dipole moments. The theorem is valid for an external field of any strength and shape, with and without rotating wave approximation in the system–field interaction. It establishes explicit relations between DS and CNS observables, significantly simplifies numerical calculations, and clarifies physical origins of the field-induced DS dynamics. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac–Frenkel time-dependent variational approach with the multiple Davydov D1Ansatz in the interaction picture. The validity of our approach is carefully checked by comparing the results with those of the hierarchy equations of motion method. By analyzing the features of nonequilibrium dynamics, we identify the phase diagrams for different bath initial conditions. We find that for the spectral exponent s < sc, there exists a transition from coherent to quasicoherent dynamics with increasing coupling strengths. For sc < s ≤ 1, the coherent to incoherent crossover occurs at a certain coupling strength and the quasicoherent dynamics emerges at much larger couplings. The initial preparation of the bath has a considerable influence on the dynamics. [ABSTRACT FROM AUTHOR]

Landau‐Zener (LZ) transitions in the three‐level bow‐tie model (3L‐BTM) in a dissipative environment are investigated by using the numerically accurate method of multiple Davydov D2$\mathrm{D}_{2}$ Ansätze. The 3L‐BTM coupled to a single harmonic mode is first considered, and evolutions of transition probabilities for selected values of model parameters are studied and interpreted with the aid of the energy diagram method. Then the 3L‐BTM coupled to a boson bath of a continuum spectral density is explored. The simulations demonstrate that sub‐Ohmic, Ohmic and super‐Ohmic boson baths have substantially different influences on the 3L‐BTM dynamics, which cannot be grasped by the standard phenomenological Markovian single‐rate descriptions. Novel bath‐induced phenomena are also described, which are absent in two‐level LZ systems. [ABSTRACT FROM AUTHOR]

By employing the time-dependent variational principle and the versatile multi-D2 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. [ABSTRACT FROM AUTHOR]

A long 0.9 ps lifetime of the upper excited singlet state in perylene is resolved by femtosecond pump–probe measurements under ultraviolet (4.96 eV) excitation and further validated by theoretical simulations of transient absorption kinetics. This finding prompts exploration and development of novel perylene-based materials for upper excited state photochemistry applications. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We have extended the multiconfigurational Ehrenfest approach to the simulation of four-wave-mixing signals of systems involving multiple electronic and vibrational degrees of freedom. As an illustration, we calculate signals of three widely used spectroscopic techniques, time- and frequency-resolved fluorescence spectroscopy, transient absorption spectroscopy, and two-dimensional (2D) electronic spectroscopy, for a two-electronic-state, twenty-four vibrational-mode conical intersection model. It has been shown that all these three spectroscopic signals characterize fast population transfer from the higher excited electronic state to the lower excited electronic state. While the time- and frequency-resolved spectrum maps the wave packet propagation exclusively on the electronically excited states, the transient absorption and 2D electronic spectra reflect the wave packet dynamics on both electronically excited states and the electronic ground state. Combining trajectory-guided Gaussian basis functions and the nonlinear response function formalism, the present approach provides a promising general technique for the applications of various Gaussian basis methods to the calculations of four-wave-mixing spectra of polyatomic molecules. [ABSTRACT FROM AUTHOR]

Experimental measurements of electron transport properties of molecular junctions are often performed in solvents. Solvent–molecule coupling and physical properties of the solvent can be used as the external stimulus to control the electric current through a molecule. In this paper, we propose a model that includes dynamical effects of solvent–molecule interaction in non-equilibrium Green's function calculations of the electric current. The solvent is considered as a macroscopic dipole moment that reorients stochastically and interacts with the electrons tunneling through the molecular junction. The Keldysh–Kadanoff–Baym equations for electronic Green's functions are solved in the time domain with subsequent averaging over random realizations of rotational variables using the Furutsu–Novikov method for the exact closure of infinite hierarchy of stochastic correlation functions. The developed theory requires the use of wideband approximation as well as classical treatment of solvent degrees of freedom. The theory is applied to a model molecular junction. It is demonstrated that not only electrostatic interaction between molecular junction and solvent but also solvent viscosity can be used to control electrical properties of the junction. Alignment of the rotating dipole moment breaks the particle–hole symmetry of the transmission favoring either hole or electron transport channels depending upon the aligning potential. [ABSTRACT FROM AUTHOR]

Finite-temperature dynamics of singlet fission in crystalline rubrene is investigated by utilizing the Dirac–Frenkel time-dependent variational method in combination with multiple Davydov D2 trial states. To probe temperature effects on the singlet fission process mediated by a conical intersection, the variational method is extended to include number state propagation with thermally averaged Boltzmann distribution as initialization. This allows us to simulate two-dimensional electronic spectroscopic signals of two-mode and three-mode models of crystalline rubrene in the temperature range from 0 K to 300 K. It is demonstrated that an elevated temperature facilitates excitonic population transfer and accelerates the singlet fission process. In addition, increasing temperature leads to dramatic changes in two-dimensional spectra, thanks to temperature-dependent electronic dephasing and to an increased number of system eigenstates amenable to spectroscopic probing. [ABSTRACT FROM AUTHOR]

We studied spectroscopic signatures of the nonadiabatic dynamics at conical intersections formed by the lowest excited singlet states in pyrazine. We considered two ab initio models of conical intersections in the excited states of pyrazine developed by Sala et al. [Phys. Chem. Chem. Phys. 16, 15957 (2014)]: a two-state (B2u and B3u), five-mode model and a three-state (B2u, B3u, and Au), nine-mode model. We simulated the signals of three widely used techniques: time- and frequency-resolved fluorescence spectroscopy, transient absorption pump–probe spectroscopy, and electronic two-dimensional spectroscopy. The signals were calculated through third-order response functions, which, in turn, were evaluated with the numerically accurate multiple Davydov ansatz. We establish spectroscopic signatures of the optically dark Au state and demonstrate that the key features of the photoinduced dynamics, such as electronic/nuclear populations, electronic/nuclear coherences, and electronic/nuclear energy transfer processes, are imprinted in the spectroscopic signals. We show that a fairly complete picture of the nonadiabatic dynamics at conical intersections can be obtained when several spectroscopic techniques are combined. Provided that the time resolution is sufficient, time- and frequency-resolved fluorescence may provide the best visualization of the nonadiabatic dynamics near conical intersections. [ABSTRACT FROM AUTHOR]

We have extended the multiconfigurational Ehrenfest (MCE) approach to investigate the dynamics of a one-dimensional Holstein molecular crystal model. It has been shown that the extended MCE approach yields results in perfect agreement with benchmark calculations by the hierarchy equations of motion method. The accuracies of the MCE approach in describing the dynamical properties of the Holstein polaron over a wide range of exciton transfer integrals and exciton-phonon couplings are carefully examined by a detailed comparison with the fully variational multiple Davydov D2 ansatz. It is found that while the MCE approach and the multi-D2 ansatz produce almost exactly the same results for a small transfer integral, the results obtained by the multi-D2 ansatz start to deviate from those by the MCE approach at longer times for a large transfer integral. A large number of coherent state basis functions are required to characterize the delocalized features of the phonon wavefunction in the case of large transfer integral, which becomes computationally too demanding for the multi-D2 ansatz. The MCE approach, on the other hand, uses hundreds to thousands of trajectory guided basis functions and converges very well, thus providing an effective tool for accurate and efficient simulations of polaron dynamics. [ABSTRACT FROM AUTHOR]