Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach

A stochastic Schrödinger equation and matrix product state approach to carrier transport in organic semiconductors with nonlocal electron-phonon interaction

Evaluation of the charge transport property of organic semiconductors requires exact quantum dynamics simulation of large systems. We present a numerically nearly exact approach to investigate carrier transport dynamics in organic semiconductors by extending the non-Markovian stochastic Schrödinger equation with complex frequency modes to a forward-backward scheme and by solving it using the matrix product state (MPS) approach. By utilizing the forward-backward formalism for noise generation, the bath correlation function can be effectively treated as a temperature-independent imaginary part, enabling a more accurate decomposition with fewer complex frequency modes. Using this approach, we study the carrier transport and mobility in the one-dimensional Peierls model, where the nonlocal electron-phonon interaction is taken into account. The reliability of this approach was validated by comparing carrier diffusion motion with those obtained from the hierarchical equations of motion method across various parameter regimes of the phonon bath. The efficiency was demonstrated by the modest virtual bond dimensions of MPS and the low scaling of the computational time with the system size.

Simulation of the Pump-Probe Spectra and Excitation Energy Relaxation of the B850 Band of the LH2 Complex in Purple Bacteria

Ultrafast spectroscopic techniques have been vital in studying excitation energy transfer (EET) in photosynthetic light harvesting complexes. In this paper, we simulate the pump-probe spectra of the B850 band of the light harvesting complex 2 (LH2) of purple bacteria, by using the hierarchical equation of motion method and the optical response function approach. The ground state bleach, stimulated emission, and excited state absorption components of the pump-probe spectra are analyzed in detail. The laser pulse-induced population dynamics are also simulated to help understand the main features of the pump-probe spectra and the EET process. It is shown that the excitation energy relaxation is an ultrafast process with multiple time scales. The first 40 fs of the pump-probe spectra is dominated by the relaxation of the k = ±1 states to both the k = 0 and higher energy states. Dynamics on a longer time scale around 200 fs reflects the relaxation of higher energy states to the k = 0 state.

MesoHOPS: Size-invariant scaling calculations of multi-excitation open quantum systems

The photoexcitation dynamics of molecular materials on the 10-100 nm length scale depend on complex interactions between electronic and vibrational degrees of freedom, rendering exact calculations difficult or intractable. The adaptive Hierarchy of Pure States (adHOPS) is a formally exact method that leverages the locality imposed by interactions between thermal environments and electronic excitations to achieve size-invariant scaling calculations for single-excitation processes in systems described by a Frenkel-Holstein Hamiltonian. Here, we extend adHOPS to account for arbitrary couplings between thermal environments and vertical excitation energies, enabling formally exact, size-invariant calculations that involve multiple excitations or states with shared thermal environments. In addition, we introduce a low-temperature correction and an effective integration of the noise to reduce the computational expense of including ultrafast vibrational relaxation in Hierarchy of Pure States (HOPS) simulations. We present these advances in the latest version of the open-source MesoHOPS library and use MesoHOPS to characterize charge separation at a one-dimensional organic heterojunction when both the electron and hole are mobile.

Dynamics of a strongly coupled quantum heat engine-Computing bath observables from the hierarchy of pure states

We present a fully quantum dynamical treatment of a quantum heat engine and its baths based on the Hierarchy of Pure States (HOPS), an exact and general method for open quantum system dynamics. We show how the change of the bath energy and the interaction energy can be determined within HOPS for arbitrary coupling strength and smooth time dependence of the modulation protocol. The dynamics of all energetic contributions during the operation can be carefully examined both in its initial transient phase and, also later, in its periodic steady state. A quantum Otto engine with a qubit as an inherently nonlinear work medium is studied in a regime where the energy associated with the interaction Hamiltonian plays an important role for the global energy balance and, thus, must not be neglected when calculating its power and efficiency. We confirm that the work required to drive the coupling with the baths sensitively depends on the speed of the modulation protocol. Remarkably, departing from the conventional scheme of well-separated phases by allowing for temporal overlap, we discover that one can even gain energy from the modulation of bath interactions. We visualize these various work contributions using the analog of state change diagrams of thermodynamic cycles. We offer a concise, full presentation of HOPS with its extension to bath observables, as it serves as a universal tool for the numerically exact description of general quantum dynamical (thermodynamic) scenarios far from the weak-coupling limit.

Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach

In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.