High accuracy exponential decomposition of bath correlation functions for arbitrary and structured spectral densities: Emerging methodologies and new approaches.

This study investigates the decomposition of bath correlation functions (BCFs) in terms of complex exponential functions, with an eye on the realistic modeling of open quantum systems based on the hierarchical equations of motion. We introduce the theoretical background of various decomposition schemes in both time and frequency domains and assess their efficiency and accuracy by demonstrating the decomposition of various BCFs. We further develop a new procedure for the decomposition of BCFs originating from highly structured spectral densities with a high accuracy and compare it with existing fitting techniques. Advantages and disadvantages of each methodology are discussed in detail with special attention to their application to the corresponding quantum dynamical problem. This work provides fundamental tools for choosing and using a variety of decomposition techniques of BCFs for the study of open quantum systems in structured environments.

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.

Thermo-Field Dynamics Approach to Photo-induced Electronic Transitions Driven by Incoherent Thermal Radiation.

The effects of thermal light-matter interaction on the dynamics of photo-induced electronic transitions in molecules are investigated using a novel first principles approach based on the thermo-field dynamics description of both the molecular vibrational modes and of the radiation field. The developed approach permits numerically accurate simulations of quantum dynamics of electronic/excitonic systems coupled to nuclear and photonic baths kept at different temperatures. The baths can be described by arbitrary spectral densities and can have any system-bath coupling strengths. In agreement with the results obtained previously by less rigorous methods, we show that the excitation process obtained by the continuous interaction with the suddenly turned-on thermal radiation field creates a mixed ensemble having a nonnegligible component consisting of a superposition of vibronic eigenstates which can sustain coherent oscillations for relatively long times. The results become especially relevant for the dynamics of electronic transitions upon sunlight excitation. Analytical results based on time-dependent perturbation theory support the numerical simulations and provide a simple interpretation of the time evolution of quantum observables.

Current-induced bond rupture in single-molecule junctions: Effects of multiple electronic states and vibrational modes.

Current-induced bond rupture is a fundamental process in nanoelectronic architectures, such as molecular junctions, and scanning tunneling microscopy measurements of molecules at surfaces. The understanding of the underlying mechanisms is important for the design of molecular junctions that are stable at higher bias voltages and is a prerequisite for further developments in the field of current-induced chemistry. In this work, we analyze the mechanisms of current-induced bond rupture employing a recently developed method, which combines the hierarchical equations of motion approach in twin space with the matrix product state formalism and allows accurate, fully quantum mechanical simulations of the complex bond rupture dynamics. Extending previous work [Ke et al. J. Chem. Phys. 154, 234702 (2021)], we consider specifically the effect of multiple electronic states and multiple vibrational modes. The results obtained for a series of models of increasing complexity show the importance of vibronic coupling between different electronic states of the charged molecule, which can enhance the dissociation rate at low bias voltages profoundly.

Effect of mixed Frenkel and charge transfer states in time-gated fluorescence spectra of perylene bisimides H-aggregates: Hierarchical equations of motion approach.

We theoretically investigated the effect of mixed Frenkel (F) and charge transfer (CT) states on the spectral properties of perylene bisimide (PBI) derivatives, focusing on the role of strong electron-phonon interactions. The model consists of a four-level system described by the Holstein Hamiltonian coupled to independent local heat-baths on each site, described by Brownian spectral distribution functions. We employ the reduced hierarchical equations of motion (HEOM) approach to calculate the time evolution of the system and compare it to the pure F exciton cases. We compute the absorption and time-gated fluorescence (TGF) spectra for different exciton transfer integrals and F-CT bandgap conditions. The coherence length of excitons (Ncoh) is evaluated employing two different definitions. We observe the presence of an excited hot state peak whose intensity is associated with the delocalization of the excited species and ultrafast dynamics that are solely dependent on the frequency of the local bath. The results indicate that the inclusion of CT states promotes localization of the excitons, which is manifested in a decrease in the intensity of the hot state peak and the 0-1 peak and an increase in the intensity of the 0-0 emission peak in the TGF spectrum, leading to a decrease of Ncoh.

Hierarchical equations of motion approach to hybrid fermionic and bosonic environments: Matrix product state formulation in twin space.

We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation [R. Borrelli, J. Chem. Phys. 150, 234102 (2019)] to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation.

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.

Transparent and Colorless Dye-Sensitized Solar Cells Exceeding 75% Average Visible Transmittance.

Most photovoltaic (PV) technologies are opaque to maximize visible light absorption. However, see-through solar cells open additional perspectives for PV integration. Looking beyond maximizing visible light harvesting, this work considers the human eye photopic response to optimize a selective near-infrared sensitizer based on a polymethine cyanine structure (VG20-C x ) to render dye-sensitized solar cells (DSSCs) fully transparent and colorless. This peculiarity was achieved by conferring to the dye the ability to strongly and sharply absorb beyond 800 nm (S0-S1 transition) while rejecting the upper S0-S n contributions far in the blue where the human retina is poorly sensitive. When associated with an aggregation-free anatase TiO2 photoanode, the selective NIR-DSSC can display 3.1% power conversion efficiency, up to 76% average visible transmittance (AVT), a value approaching the 78% AVT value of a standard double glazing window while reaching a color rendering index (CRI) of 92.1%. The ultrafast and fast charge transfer processes are herein discussed, clarifying the different relaxation channels from the dye monomer excited states and highlighting the limiting steps to provide future directions to enhance the performances of this nonintrusive NIR-DSSC technology.

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.

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.

Expanding the Range of Hierarchical Equations of Motion by Tensor-Train Implementation.

The non-equilibrium thermo-field dynamics formulation of the hierarchical equations of motion combined with the tensor-train representation of the density matrix is discussed, and a new numerical integration scheme is introduced. The numerical methodology is based on an adaptive low-rank Galerkin reduction scheme and can preserve linear invariants (such as the trace of the density matrix). The method is applied to the study of the charge transfer dynamics in model pentacene molecular aggregates. The combined effect of a discrete set of molecular vibrational modes and a thermal bath is investigated, paying special attention to the coherent-incoherent transition of the charge transport. The new computational framework is shown to be a very promising methodology for the study of the quantum dynamics of complex molecular systems in the condensed phase.

Probing photoinduced proton coupled electron transfer process by means of two-dimensional resonant electronic-vibrational spectroscopy.

We develop a detailed theoretical model of photo-induced proton-coupled electron transfer (PPCET) processes, which are at the basis of solar energy harvesting in biological systems and photovoltaic materials. Our model enables us to analyze the dynamics and the efficiency of a PPCET reaction under the influence of a thermal environment by disentangling the contribution of the fundamental electron transfer and proton transfer steps. In order to study quantum dynamics of the PPCET process under an interaction with the non-Markovian environment, we employ the hierarchical equations of motion. We calculate transient absorption spectroscopy (TAS) and a newly defined two-dimensional resonant electronic-vibrational spectroscopy (2DREVS) signals in order to study the nonequilibrium reaction dynamics. Our results show that different transition pathways can be separated by TAS and 2DREVS.

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.

Reliable Predictions of Benzophenone Singlet-Triplet Transition Rates: A Second-Order Cumulant Approach.

Fermi golden rule and second-order cumulant expansion of the time-dependent density matrix have been used to compute from first principles the rate of intersystem crossing in benzophenone, using minimum-energy geometries and normal modes of vibrations computed at the TDDFT/CAM-B3LYP level. Both approaches yield reliable values of the S1 decay rate, the latter being almost in quantitative agreement with the results of time-dependent spectroscopic measurements (0.154 ps-1 observed vs 0.25 ps-1 predicted). The Fermi golden rule slightly overestimates the decay rate of S1 state (kd = 0.45 ps-1) but provides better insights into the chemico-physical parameters, which govern the transition from a thermally equilibrated population of S1, showing that the indirect mechanism is much faster than the direct one because of the vanishingly small Franck-Condon weighted density of states at ΔE of transition.

Proton tunneling in a two-dimensional potential energy surface with a non-linear system-bath interaction: Thermal suppression of reaction rate.

We consider a proton transfer (PT) system described by a proton transfer reaction (PTR) coordinate and a rate promoting vibrational (RPV) coordinate interacting with a non-Markovian heat bath. While dynamics of PT processes has been widely discussed using two-dimensional potential energy surfaces, the role of the heat bath, in particular, in a realistic form of the system-bath interaction has not been well explored. Previous studies are largely based on a one-dimensional model and linear-linear system-bath interaction. In the present study, we introduce an exponential-linear (EL) system-bath interaction, which is derived from the analysis of a PTR-RPV system in a realistic situation. This interaction mainly causes vibrational dephasing in the PTR mode and population relaxation in the RPV mode. Numerical simulations were carried out using the hierarchical equations of motion approach. We analyze the role of the heat bath interaction in the chemical reaction rate as a function of the system-bath coupling strength at different temperatures and for different values of the bath correlation time. A prominent feature of the present result is that while the reaction rate predicted from classical and quantum Kramers theory increases as the temperature increases, the present EL interaction model exhibits opposite temperature dependence. The Kramers turn-over profile of the reaction rate as a function of the system-bath coupling is also suppressed in the present EL model, turning into a plateau-like curve for larger system-bath interaction strength. Such features arise from the interplay of the vibrational dephasing process in the PTR mode and the population relaxation process in the RPV mode.

Tunable Redox Potential, Optical Properties, and Enhanced Stability of Modified Ferrocene-Based Complexes.

We report a series of ferrocene-based derivatives and their corresponding oxidized forms in which the introduction of simple electron donating groups like methyl or tert-butyl units on cyclopentadienyl-rings afford great tunability of Fe+III/Fe+II redox potentials from +0.403 V down to -0.096 V versus saturated calomel electrode. The spin forbidden d-d transitions of ferrocene derivatives shift slightly toward the blue region with an increasing number of electron-donating groups on the cyclopentadienyl-rings with very little change in absorptivity values, whereas the ligand-to-metal transitions of the corresponding ferricinium salts move significantly to the near-IR region. The electron-donating groups also contribute in the strengthening of electron density of Fe+III d-orbitals, which therefore improves the chemical stability against the oxygen reaction. Further, density functional theory calculations show a reducing trend in outer shell reorganization energy with an increasing number of the electron donating units.

Density matrix dynamics in twin-formulation: An efficient methodology based on tensor-train representation of reduced equations of motion.

The twin-formulation of quantum statistical mechanics is employed to describe a new methodology for the solution of the equations of motion of the reduced density matrix in their hierarchical formulation. It is shown that the introduction of tilde operators and of their algebra in the dual space greatly simplifies the application of numerical techniques for the propagation of the density matrix. The application of tensor-train representation of a vector to solve complex quantum dynamical problems within the framework of the twin-formulation is discussed. Next, applications of the hierarchical equations of motion to a dissipative polaron model are presented showing the validity and accuracy of the new approach.

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

Unraveling the many facets of coherent and incoherent exciton motion in an ensemble of chromophores is an inherently complex quantum mechanical problem that has triggered a vivid debate on the role of quantum effects in molecular materials and biophysical systems. Here the dynamics of a statistical ensemble of molecular aggregates consisting of identical chromophores is investigated within a new theoretical framework. Taking account of intrinsic properties of the system, the Hamiltonian of the aggregate is partitioned into two mutually commuting vibrational and vibronic operators. This representation paves the way for an analysis that reveals the role of static disorder in ensembles of aggregates. Using analytical methods, it is demonstrated that after a critical time τD ≃ 2π/σ (σ being the dispersion of the disorder) any dynamic variable of the aggregate exhibits purely vibrational dynamics. This result is confirmed by exact numerical calculations of the time-dependent site populations of the aggregate. These findings may be useful for the interpretation of recent femtosecond spectroscopic experiments on molecular aggregates.

Transient and Enduring Electronic Resonances Drive Coherent Long Distance Charge Transport in Molecular Wires.

It is shown that the yields of oxidative damage observed in double-stranded DNA oligomers consisting of two guanines separated by adenine-thymine (A:T) n bridges of various lengths are reliably accounted for by a multistep mechanism, in which transient and nontransient electronic resonances induce charge transport and solvent relaxation stabilizes the hole transfer products. The proposed multistep mechanism leads to results in excellent agreement with the observed yield ratios for both the short and the long distance regime; the almost distance independence of yield ratios for longer bridges ( n ≥ 3) is the consequence of the significant energy decrease of the electronic levels of the bridge, which, as the bridge length increases, become quasi-degenerate with those of the acceptor and donor groups (enduring resonance). These results provide significant guidelines for the design of novel DNA sequences to be employed in organic electronics.

Hole Hopping Rates in Organic Semiconductors: A Second-Order Cumulant Approach.

Second-order cumulant expansion of the time dependent reduced density matrix has been employed to evaluate hole hopping rates in pentacene, tetracene, picene, and rubrene homodimers. The cumulant expansion is a full quantum mechanical approach, which enables the use of the whole set of nuclear coordinates in computations and the inclusion of both the effects of the equilibrium position displacements and of normal mode mixing upon hole transfer. The time dependent populations predicted by cumulant approach are in good agreement with those obtained by numerical solution of time dependent Schrödinger equation, even for ultrafast processes, where the Fermi Golden Rule fails.