Flexible scheme to truncate the hierarchy of pure states

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.

Formally Exact Simulations of Mesoscale Exciton Diffusion in a Light-Harvesting 2 Antenna Nanoarray

The photosynthetic apparatus of plants and bacteria combine atomically precise pigment-protein complexes with dynamic membrane architectures to control energy transfer on the 10-100 nm length scales. Recently, synthetic materials have integrated photosynthetic antenna proteins to enhance exciton transport, though the influence of artificial packing on the excited-state dynamics in these biohybrid materials is not fully understood. Here, we use the adaptive hierarchy of pure states (adHOPS) to perform a formally exact simulation of excitation energy transfer within artificial aggregates of light-harvesting complex 2 (LH2) with a range of packing densities. We find that LH2 aggregates support a remarkable exciton diffusion length ranging from 100 nm at a biological packing density to 300 nm at the densest packing previously suggested in an artificial aggregate. The unprecedented scale of these formally exact calculations also underscores the efficiency with which adHOPS simulates excited-state processes in molecular materials.

Calculating non-linear response functions for multi-dimensional electronic spectroscopy using dyadic non-Markovian quantum state diffusion

We present a methodology for simulating multi-dimensional electronic spectra of molecular aggregates with coupling of electronic excitation to a structured environment using the stochastic non-Markovian quantum state diffusion (NMQSD) method in combination with perturbation theory for the response functions. A crucial aspect of our approach is that we propagate the NMQSD equation in a doubled system Hilbert space, but with the same noise. We demonstrate that our approach shows fast convergence with respect to the number of stochastic trajectories, providing a promising technique for numerical calculation of two-dimensional electronic spectra of large molecular aggregates.

Simulation of absorption spectra of molecular aggregates: A hierarchy of stochastic pure state approach

Simulation of spectroscopic observables for molecular aggregates with strong and structured coupling of electronic excitation to vibrational degrees of freedom is an important but challenging task. The Hierarchy of Pure States (HOPS) provides a formally exact solution based on local, stochastic trajectories. Exploiting the localization of HOPS for the simulation of absorption spectra in large aggregates requires a formulation in terms of normalized trajectories. Here, we provide a normalized dyadic equation where the ket- and bra-states are propagated in different electronic Hilbert spaces. This work opens the door to applying adaptive HOPS methods for the simulation of absorption spectra.