Efficient Simulation of Finite-Temperature Open Quantum Systems

Effective modeling of open quantum systems by low-rank discretization of structured environments

The accurate description of the interaction of a quantum system with its environment is a challenging problem ubiquitous across all areas of physics and lies at the foundation of quantum mechanics theory. Here, we pioneer a new strategy to create discrete low-rank models of the system-environment interaction, by exploiting the frequency and time domain information encoded in the fluctuation-dissipation relation connecting the system-bath correlation function and the spectral density. We demonstrate the effectiveness of our methodology by combining it with tensor-network methodologies and simulating the quantum dynamics of complex excitonic systems in a highly structured bosonic environment. The new modeling framework sets the basis for a leap in the analysis of open quantum systems, providing controlled accuracy at significantly reduced computational costs, with benefits in all connected research areas.

Efficient simulation of open quantum systems coupled to a reservoir through multiple channels

It is challenging to simulate open quantum systems that are connected to a reservoir through multiple channels. For example, vibrations may induce fluctuations in both energy gaps and electronic couplings, which represent two independent channels of system-bath couplings. Systems of this kind are ubiquitous in the processes of excited state radiationless decay. Combined with density matrix renormalization group (DMRG) and matrix product states (MPS) methods, we develop an interaction-picture chain mapping strategy for vibrational reservoirs to simulate the dynamics of these open systems, resulting in time-dependent spatially local system-bath couplings in the chain-mapped Hamiltonian. This transformation causes the entanglement generated by the system-bath interactions to be restricted within a narrow frequency window of vibrational modes, enabling efficient DMRG/MPS dynamical simulations. We demonstrate the utility of this approach by simulating singlet fission dynamics using a generalized spin-boson Hamiltonian with both diagonal and off-diagonal system-bath couplings. This approach generalizes an earlier interaction-picture chain mapping scheme, allowing for efficient and exact simulation of systems with multi-channel system-bath couplings using matrix product states, which may further our understanding of nonlocal exciton-phonon couplings in exciton transport and the non-Condon effect in energy and electron transfer.

Extended system-bath entanglement theorem with multiple baths in the presence of external fields

The system-bath entanglement theorem (SBET) was established in terms of linear response functions [Du et al., J. Chem. Phys. 152, 034102 (2020)] and generalized to correlation functions [Su et al., J. Chem. Phys. 160, 084104 (2024)] in our previous studies. This theorem connects the entangled system-bath properties to the local system and bare-bath ones. In this work, we extend the SBET to field-dressed conditions with multiple baths at different temperatures. As in reality, the external fields may interact with not only the system but also environments. The extended SBET facilitates, for example, photo-acoustic, photo-thermal, pump-probe related studies. The theorem under the field-free condition (multiple baths) and its counterpart in the classical limit is also presented.

Quantum neural network approach to Markovian dissipative dynamics of many-body open quantum systems

Numerous variational methods have been proposed for solving quantum many-body systems, but they often face exponentially increasing computational complexity as the Hilbert space dimension grows. To address this, we introduce a novel approach using quantum neural networks to simulate the dissipative dynamics of many-body open quantum systems. This method combines neural-network quantum state representation with the time-dependent variational principle, both implemented via quantum algorithms. This results in accurate open quantum dynamics described by the Lindblad quantum master equation, exemplified by the spin-boson and transverse field Ising models. Our approach avoids the computational expense of classical algorithms and demonstrates the potential advantages of quantum computing for many-body simulations. To reduce measurement errors, we introduce a projection reset procedure, which could benefit other quantum simulations. In addition, our approach can be extended to simulate non-Markovian quantum dynamics.

MPSDynamics.jl: Tensor network simulations for finite-temperature (non-Markovian) open quantum system dynamics

The MPSDynamics.jl package provides an easy-to-use interface for performing open quantum systems simulations at zero and finite temperatures. The package has been developed with the aim of studying non-Markovian open system dynamics using the state-of-the-art numerically exact Thermalized-Time Evolving Density operator with Orthonormal Polynomials Algorithm based on environment chain mapping. The simulations rely on a tensor network representation of the quantum states as matrix product states (MPS) and tree tensor network states. Written in the Julia programming language, MPSDynamics.jl is a versatile open-source package providing a choice of several variants of the Time-Dependent Variational Principle method for time evolution (including novel bond-adaptive one-site algorithms). The package also provides strong support for the measurement of single and multi-site observables, as well as the storing and logging of data, which makes it a useful tool for the study of many-body physics. It currently handles long-range interactions, time-dependent Hamiltonians, multiple environments, bosonic and fermionic environments, and joint system-environment observables.

ACE: A general-purpose non-Markovian open quantum systems simulation toolkit based on process tensors

We describe a general-purpose computational toolkit for simulating open quantum systems, which provides numerically exact solutions for composites of zero-dimensional quantum systems that may be strongly coupled to multiple, quite general non-Markovian environments. It is based on process tensor matrix product operators (PT-MPOs), which efficiently encapsulate environment influences. The code features implementations of several PT-MPO algorithms, in particular Automated Compression of Environments for general environments comprised of independent modes as well as schemes for generalized spin boson models. The latter includes a divide-and-conquer scheme for periodic PT-MPOs, which enable million time step simulations for realistic models. PT-MPOs can be precalculated and reused for efficiently probing different time-dependent system Hamiltonians. They can also be stacked together and combined to provide numerically complete solutions of small networks of open quantum systems. The code is written in C++ and is fully controllable by configuration files, for which we have developed a versatile and compact human-readable format.

Kylin-V: An open-source package calculating the dynamic and spectroscopic properties of large systems

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system's complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Managing temperature in open quantum systems strongly coupled with structured environments

In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer.

General framework for quantifying dissipation pathways in open quantum systems. II. Numerical validation and the role of non-Markovianity

In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1-214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics.

Bexcitonics: Quasiparticle approach to open quantum dynamics

We develop a quasiparticle approach to capture the dynamics of open quantum systems coupled to bosonic thermal baths of arbitrary complexity based on the Hierarchical Equations of Motion (HEOM). This is done by generalizing the HEOM dynamics and mapping it into that of the system in interaction with a few bosonic fictitious quasiparticles that we call bexcitons. Bexcitons arise from a decomposition of the bath correlation function into discrete features. Specifically, bexciton creation and annihilation couple the auxiliary density matrices in the HEOM. The approach provides a systematic strategy to construct exact quantum master equations that include the system-bath coupling to all orders even for non-Markovian environments. Specifically, by introducing different metrics and representations for the bexcitons it is possible to straightforwardly generate different variants of the HEOM, demonstrating that all these variants share a common underlying quasiparticle picture. Bexcitonic properties, while unphysical, offer a coarse-grained view of the correlated system-bath dynamics and its numerical convergence. For instance, we use it to analyze the instability of the HEOM when the bath is composed of underdamped oscillators and show that it leads to the creation of highly excited bexcitons. The bexcitonic picture can also be used to develop more efficient approaches to propagate the HEOM. As an example, we use the particle-like nature of the bexcitons to introduce mode-combination of bexcitons in both number and coordinate representation that uses the multi-configuration time-dependent Hartree to efficiently propagate the HEOM dynamics.

Environmentally Driven Symmetry Breaking Quenches Dual Fluorescence in Proflavine

Nonadiabatic couplings between several electronic excited states are ubiquitous in many organic chromophores and can significantly influence optical properties. A recent experimental study demonstrated that the proflavine molecule exhibits surprising dual fluorescence in the gas phase, which is suppressed in polar solvent environments. Here, we uncover the origin of this phenomenon by parametrizing a linear-vibronic coupling Hamiltonian from spectral densities of system-bath coupling constructed along molecular dynamics trajectories, fully accounting for interactions with the condensed-phase environment. The finite-temperature absorption, steady-state emission, and time-resolved emission spectra are then computed using powerful, numerically exact tensor network approaches. We find that the dual fluorescence in vacuum is driven by a single well-defined coupling mode but is quenched in solution due to dynamic solvent-driven symmetry breaking that mixes the two low-lying electronic states. We expect the computational framework developed here to be widely applicable to the study of non-Condon effects in complex condensed-phase environments.

Transient temperature dynamics of reservoirs connected through an open quantum system

The dynamics of open quantum systems connected with several reservoirs attract great attention due to their importance in quantum optics, biology, quantum thermodynamics, transport phenomena, etc. In many problems, the Born approximation is applicable, which implies that the influence of the open quantum system on the reservoirs can be neglected. However, in the case of long-time dynamics or mesoscopic reservoirs, the reverse influence can be crucial. In this paper, we investigate the transient dynamics of several bosonic reservoirs connected through an open quantum system. We use an adiabatic approach to study the temporal dynamics of temperatures of the reservoirs during relaxation to thermodynamic equilibrium. We show that there are various types of temperature dynamics that strongly depend on the values of dissipative rates and initial temperatures. We demonstrate that temperatures of the reservoirs, including the hottest and coldest ones, can exhibit nonmonotonic behavior. Moreover, there are moments of time during which the reservoir with an initially intermediate temperature becomes the hottest or coldest reservoir. The obtained results pave the way for managing energy flows in mesoscale and nanoscale systems.

Microscopic theory of exciton-polariton model involving multiple molecules: Macroscopic quantum electrodynamics formulation and essence of direct intermolecular interactions

Cavity quantum electrodynamics (CQED) and its extensions are widely used for the description of exciton-polariton systems. However, the exciton-polariton models based on CQED vary greatly within different contexts. One of the most significant discrepancies among these CQED models is whether one should include direct intermolecular interactions in the CQED Hamiltonian. To answer this question, in this article, we derive an effective dissipative CQED model including free-space dipole-dipole interactions (CQED-DDI) from a microscopic Hamiltonian based on macroscopic quantum electrodynamics. Dissipative CQED-DDI successfully captures the nature of vacuum fluctuations in dielectric media and separates them into free-space effects and dielectric-induced effects. The former include spontaneous emissions, dephasings, and dipole-dipole interactions in free space; the latter include exciton-polariton interactions and photonic losses due to dielectric media. We apply dissipative CQED-DDI to investigate the exciton-polariton dynamics (the population dynamics of molecules above a plasmonic surface) and compare the results with those based on the methods proposed by several previous studies. We find that direct intermolecular interactions are a crucial element when employing CQED-like models to study exciton-polariton systems involving multiple molecules.

Systematic Coarse Graining of Environments for the Nonperturbative Simulation of Open Quantum Systems

Conducting precise electronic-vibrational dynamics simulations of molecular systems poses significant challenges when dealing with realistic environments composed of numerous vibrational modes. Here, we introduce a technique for the construction of effective phonon spectral densities that capture accurately open-system dynamics over a finite time interval of interest. When combined with existing nonperturbative simulation tools, our approach can reduce significantly the computational costs associated with many-body open-system dynamics.

Generalized system-bath entanglement theorem for Gaussian environments

The entanglement between system and bath often plays a pivotal role in complex systems spanning multiple orders of magnitude. A system-bath entanglement theorem was previously established for Gaussian environments in J. Chem. Phys. 152, 034102 (2020) regarding linear response functions. This theorem connects the entangled responses to the local system and bare bath properties. In this work, we generalize it to correlation functions. Key steps in derivations involve using the generalized Langevin dynamics for hybridizing bath modes and the Bogoliubov transformation that maps the original finite-temperature reservoir to an effective zero-temperature vacuum by employing an auxiliary bath. The generalized theorem allows us to evaluate the system-bath entangled correlations and the bath mode correlations in the total composite space, as long as we know the bare-bath statistical properties and obtain the reduced system correlations. To demonstrate the cross-scale entanglements, we utilize the generalized theorem to calculate the solvation free energy of an electron transfer system with intramolecular vibrational modes.

Anomalous Transport of Small Polarons Arises from Transient Lattice Relaxation or Immovable Boundaries

Elucidating transport mechanisms is crucial for advancing material design, yet state-of-the-art theory is restricted to exact simulations of small lattices with severe finite-size effects or approximate ones that assume the nature of transport. We leverage algorithmic advances to tame finite-size effects and exactly simulate small polaron formation and transport in the Holstein model. We further analyze the applicability of the ubiquitously used equilibrium-based Green-Kubo relations and nonequilibrium methods to predict charge mobility. We find that these methods can converge to different values and track this disparity to finite-size dependence and the sensitivity of Green-Kubo relations to the system's topology. Contrary to standard perturbative calculations, our results demonstrate that small polarons exhibit anomalous transport that manifests transiently due to nonequilibrium lattice relaxation or permanently as a signature of immovable boundaries. These findings can offer new interpretations of transport experiments on polymers and transition metal oxides.

Modeling the dynamics of quantum systems coupled to large-dimensional baths using effective energy states

The quantum dynamics of a low-dimensional system in contact with a large but finite harmonic bath is theoretically investigated by coarse-graining the bath into a reduced set of effective energy states. In this model, the couplings between the system and the bath are obtained from statistically averaging over the discrete, degenerate effective states. Our model is aimed at intermediate bath sizes in which non-Markovian processes and energy transfer between the bath and the main system are important. The method is applied to a model system of a Morse oscillator coupled to 40 harmonic modes. The results are found to be in excellent agreement with the direct quantum dynamics simulations presented in the work of Bouakline et al. [J. Phys. Chem. A 116, 11118-11127 (2012)], but at a much lower computational cost. Extension to larger baths is discussed in comparison to the time-convolutionless method. We also extend this study to the case of a microcanonical bath with finite initial internal energies. The computational efficiency and convergence properties of the effective bath states model with respect to relevant parameters are also discussed.

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.

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.

Unraveling the contributions to the spectral shape of flexible dyes in solution: insights on the absorption spectrum of an oxyluciferin analogue

We present a computational investigation of the absorption spectrum in water of 5,5-spirocyclopropyl-oxyluciferin (5,5-CprOxyLH), an analogue of the emitter compound responsible for the bioluminescence in fireflies. Several factors participate in determining the 5,5-CprOxyLH's spectral shape: (i) the contribution of the four close-energy excited states, which show significant non-adiabatic couplings, (ii) the flexible molecular structure and (iii) the specific interactions established with the surrounding environment, which strongly couple the protic solvent dynamics with the dye's spectral response. To tackle the challenge to capture and dissect the role of all these effects we preliminarily investigate the role of non-adiabatic couplings with quantum dynamics simulations and a linear vibronic coupling model in the gas phase. Then, we account for both the molecular flexibility and solvent interactions by resorting to a mixed quantum classical protocol, named Adiabatic Molecular Dynamics generalized Vertical Gradient (Ad-MD|gVG), which is built on a method recently proposed by some of us. It is rooted in the partition between stiff degrees of freedom of the dye, accounted for at the vibronic level within the harmonic approximation, and flexible degrees of freedom of the solute (and of the solvent), described classically through a sampling based on Molecular Dynamics (MD). Ad-MD|gVG avoids spurious effects arising in the excited state Hessians due to non-adiabatic couplings, and can therefore be applied to account for the contributions of the first four excited states to the 5,5-CprOxyLH absorption spectrum. The final simulated spectrum is in very good agreement with the experiment, especially when the MD is driven by a refined quantum-mechanically derived force-field. More importantly, the origin of each separate contribution to the spectral shape is appropriately accounted for, paving the way to future applications of the method to more complex systems or alternative spectroscopies, as emission or circular dichroism.

Compact and complete description of non-Markovian dynamics

Generalized master equations provide a theoretically rigorous framework to capture the dynamics of processes ranging from energy harvesting in plants and photovoltaic devices to qubit decoherence in quantum technologies and even protein folding. At their center is the concept of memory. The explicit time-nonlocal description of memory is both protracted and elaborate. When physical intuition is at a premium, one would desire a more compact, yet complete, description. Here, we demonstrate how and when the time-convolutionless formalism constitutes such a description. In particular, by focusing on the dissipative dynamics of the spin-boson and Frenkel exciton models, we show how to: easily construct the time-local generator from reference reduced dynamics, elucidate the dependence of its existence on the system parameters and the choice of reduced observables, identify the physical origin of its apparent divergences, and offer analysis tools to diagnose their severity and circumvent their deleterious effects. We demonstrate that, when applicable, the time-local approach requires as little information as the more commonly used time-nonlocal scheme, with the important advantages of providing a more compact description, greater algorithmic simplicity, and physical interpretability. We conclude by introducing the discrete-time analog and a straightforward protocol to employ it in cases where the reference dynamics have limited resolution. The insights we present here offer the potential for extending the reach of dynamical methods, reducing both their cost and conceptual complexity.

Fingerprint and Universal Markovian Closure of Structured Bosonic Environments

We exploit the properties of chain mapping transformations of bosonic environments to identify a finite collection of modes able to capture the characteristic features, or fingerprint, of the environment. Moreover we show that the countable infinity of residual bath modes can be replaced by a universal Markovian closure, namely, a small collection of damped modes undergoing a Lindblad-type dynamics whose parametrization is independent of the spectral density under consideration. We show that the Markovian closure provides a quadratic speedup with respect to standard chain mapping techniques and makes the memory requirement independent of the simulation time, while preserving all the information on the fingerprint modes. We illustrate the application of the Markovian closure to the computation of linear spectra but also to nonlinear spectral response, a relevant experimentally accessible many body coherence witness for which efficient numerically exact calculations in realistic environments are currently lacking.

Few-mode field quantization for multiple emitters

The control of the interaction between quantum emitters using nanophotonic structures holds great promise for quantum technology applications, while its theoretical description for complex nanostructures is a highly demanding task as the electromagnetic (EM) modes form a high-dimensional continuum. We here introduce an approach that permits a quantized description of the full EM field through a small number of discrete modes. This extends the previous work in ref. (I. Medina, F. J. García-Vidal, A. I. Fernández-Domínguez, and J. Feist, "Few-mode field quantization of arbitrary electromagnetic spectral densities," Phys. Rev. Lett., vol. 126, p. 093601, 2021) to the case of an arbitrary number of emitters, without any restrictions on the emitter level structure or dipole operators. The low computational demand of this method makes it suitable for studying dynamics for a wide range of parameters. We illustrate the power of our approach for a system of three emitters placed within a hybrid metallodielectric photonic structure and show that excitation transfer is highly sensitive to the properties of the hybrid photonic-plasmonic modes.

Excitations of Quantum Many-Body Systems via Purified Ensembles: A Unitary-Coupled-Cluster-Based Approach

State-average calculations based on a mixture of states are increasingly being exploited across chemistry and physics as versatile procedures for addressing excitations of quantum many-body systems. If not too many states should need to be addressed, calculations performed on individual states are also a common option. Here we show how the two approaches can be merged into one method, dealing with a generalized yet single pure state. Implications in electronic structure calculations are discussed and for quantum computations are pointed out.

A statistical quasi-particles thermofield theory with Gaussian environments: System-bath entanglement theorem for nonequilibrium correlation functions

For open quantum systems, environmental dissipative effect can be represented by statistical quasi-particles, namely dissipatons. We exploit this fact to establish the dissipaton thermofield theory. The resulting generalized Langevin dynamics of absorptive and emissive thermofield operators are effectively noise-resolved. The system-bath entanglement theorem is then readily followed between a important class of nonequilibrium steady-state correlation functions. All these relations are validated numerically. A simple corollary is the transport current expression, which exactly recovers the result obtained from the nonequilibrium Green's function formalism.

Exact simulation of pigment-protein complexes unveils vibronic renormalization of electronic parameters in ultrafast spectroscopy

The primary steps of photosynthesis rely on the generation, transport, and trapping of excitons in pigment-protein complexes (PPCs). Generically, PPCs possess highly structured vibrational spectra, combining many discrete intra-pigment modes and a quasi-continuous of protein modes, with vibrational and electronic couplings of comparable strength. The intricacy of the resulting vibronic dynamics poses significant challenges in establishing a quantitative connection between spectroscopic data and underlying microscopic models. Here we show how to address this challenge using numerically exact simulation methods by considering two model systems, namely the water-soluble chlorophyll-binding protein of cauliflower and the special pair of bacterial reaction centers. We demonstrate that the inclusion of the full multi-mode vibronic dynamics in numerical calculations of linear spectra leads to systematic and quantitatively significant corrections to electronic parameter estimation. These multi-mode vibronic effects are shown to be relevant in the longstanding discussion regarding the origin of long-lived oscillations in multidimensional nonlinear spectra.

Multimode vibronic mixing in model photosynthetic systems revealed by numerically exact simulations is shown to strongly modify linear and non-linear optical responses and facilitate the persistence of coherent dynamics.

Correlations, Information Backflow, and Objectivity in a Class of Pure Dephasing Models

We critically examine the role that correlations established between a system and fragments of its environment play in characterising the ensuing dynamics. We employ a dephasing model with different initial conditions, where the state of the initial environment represents a tunable degree of freedom that qualitatively and quantitatively affects the correlation profiles, but nevertheless results in the same reduced dynamics for the system. We apply recently developed tools for the characterisation of non-Markovianity to carefully assess the role that correlations, as quantified by the (quantum) Jensen-Shannon divergence and relative entropy, as well as changes in the environmental state, play in whether the conditions for classical objectivity within the quantum Darwinism paradigm are met. We demonstrate that for precisely the same non-Markovian reduced dynamics of the system arising from different microscopic models, some exhibit quantum Darwinistic features, while others show that no meaningful notion of classical objectivity is present. Furthermore, our results highlight that the non-Markovian nature of an environment does not a priori prevent a system from redundantly proliferating relevant information, but rather it is the system's ability to establish the requisite correlations that is the crucial factor in the manifestation of classical objectivity.

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S₁ and the dark S₂ states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Accurate Truncations of Chain Mapping Models for Open Quantum Systems

The dynamics of open quantum systems are of great interest in many research fields, such as for the interaction of a quantum emitter with the electromagnetic modes of a nanophotonic structure. A powerful approach for treating such setups in the non-Markovian limit is given by the chain mapping where an arbitrary environment can be transformed to a chain of modes with only nearest-neighbor coupling. However, when long propagation times are desired, the required long chain lengths limit the utility of this approach. We study various approaches for truncating the chains at manageable lengths while still preserving an accurate description of the dynamics. We achieve this by introducing losses to the chain modes in such a way that the effective environment acting on the system remains unchanged, using a number of different strategies. Furthermore, we demonstrate that extending the chain mapping to allow next-nearest neighbor coupling permits the reproduction of an arbitrary environment, and adding longer-range interactions does not further increase the effective number of degrees of freedom in the environment.

Formally exact simulations of mesoscale exciton dynamics in molecular materials

Excited state carriers, such as excitons, can diffuse on the 100 nm to micron length scale in molecular materials but only delocalize over short length scales due to coupling between electronic and vibrational degrees-of-freedom. Here, we leverage the locality of excitons to adaptively solve the hierarchy of pure states equations (HOPS). We demonstrate that our adaptive HOPS (adHOPS) methodology provides a formally exact and size-invariant (i.e., ) scaling algorithm for simulating mesoscale quantum dynamics. Finally, we provide proof-of-principle calculations for exciton diffusion on linear chains containing up to 1000 molecules.

The adaptive hierarchy of pure states (adHOPS) algorithm leverages the locality of excitons in molecular materials to perform formally-exact simulations with size-invariant (i.e., ) scaling, enabling efficient simulations of mesoscale exciton dynamics.

Efficient Exploration of Hamiltonian Parameter Space for Optimal Control of Non-Markovian Open Quantum Systems

We present a general method to efficiently design optimal control sequences for non-Markovian open quantum systems, and illustrate it by optimizing the shape of a laser pulse to prepare a quantum dot in a specific state. The optimization of control procedures for quantum systems with strong coupling to structured environments-where time-local descriptions fail-is a computationally challenging task. We modify the numerically exact time evolving matrix product operator (TEMPO) method, such that it allows the repeated computation of the time evolution of the reduced system density matrix for various sets of control parameters at very low computational cost. This method is potentially useful for studying numerous optimal control problems, in particular in solid state quantum devices where the coupling to vibrational modes is typically strong.

Insights into Photosynthetic Energy Transfer Gained from Free-Energy Structure: Coherent Transport, Incoherent Hopping, and Vibrational Assistance Revisited

Giant strides in ultrashort laser pulse technology have enabled real-time observation of dynamical processes in complex molecular systems. Specifically, the discovery of oscillatory transients in the two-dimensional electronic spectra of photosynthetic systems stimulated a number of theoretical investigations exploring the possible physical mechanisms of the remarkable quantum efficiency of light harvesting processes. In this work, we revisit the elementary aspects of environment-induced fluctuations in the involved electronic energies and present a simple way to understand energy flow with the intuitive picture of relaxation in a funnel-type free-energy landscape. The presented free-energy description of energy transfer reveals that typical photosynthetic systems operate in an almost barrierless regime. The approach also provides insights into the distinction between coherent and incoherent energy transfer and the criteria by which the necessity of the vibrational assistance is considered.

Simulating Quantum Vibronic Dynamics at Finite Temperatures With Many Body Wave Functions at 0 K

For complex molecules, nuclear degrees of freedom can act as an environment for the electronic “system” variables, allowing the theory and concepts of open quantum systems to be applied. However, when molecular system-environment interactions are non-perturbative and non-Markovian, numerical simulations of the complete system-environment wave function become necessary. These many body dynamics can be very expensive to simulate, and extracting finite-temperature results—which require running and averaging over many such simulations—becomes especially challenging. Here, we present numerical simulations that exploit a recent theoretical result that allows dissipative environmental effects at finite temperature to be extracted efficiently from a single, zero-temperature wave function simulation. Using numerically exact time-dependent variational matrix product states, we verify that this approach can be applied to vibronic tunneling systems and provide insight into the practical problems lurking behind the elegance of the theory, such as the rapidly growing numerical demands that can appear for high temperatures over the length of computations.

Matrix Product State Simulations of Non-Equilibrium Steady States and Transient Heat Flows in the Two-Bath Spin-Boson Model at Finite Temperatures

Simulating the non-perturbative and non-Markovian dynamics of open quantum systems is a very challenging many body problem, due to the need to evolve both the system and its environments on an equal footing. Tensor network and matrix product states (MPS) have emerged as powerful tools for open system models, but the numerical resources required to treat finite-temperature environments grow extremely rapidly and limit their applications. In this study we use time-dependent variational evolution of MPS to explore the striking theory of Tamascelli et al. (Phys. Rev. Lett. 2019, 123, 090402.) that shows how finite-temperature open dynamics can be obtained from zero temperature, i.e., pure wave function, simulations. Using this approach, we produce a benchmark dataset for the dynamics of the Ohmic spin-boson model across a wide range of coupling strengths and temperatures, and also present a detailed analysis of the numerical costs of simulating non-equilibrium steady states, such as those emerging from the non-perturbative coupling of a qubit to baths at different temperatures. Despite ever-growing resource requirements, we find that converged non-perturbative results can be obtained, and we discuss a number of recent ideas and numerical techniques that should allow wide application of MPS to complex open quantum systems.

Analytic expressions for the steady-state current with finite extended reservoirs

Open-system simulations of quantum transport provide a platform for the study of true steady states, Floquet states, and the role of temperature, time dynamics, and fluctuations, among other physical processes. They are rapidly gaining traction, especially techniques that revolve around "extended reservoirs," a collection of a finite number of degrees of freedom with relaxation that maintains a bias or temperature gradient, and have appeared under various guises (e.g., the extended or mesoscopic reservoir, auxiliary master equation, and driven Liouville-von Neumann approaches). Yet, there are still a number of open questions regarding the behavior and convergence of these techniques. Here, we derive general analytical solutions, and associated asymptotic analyses, for the steady-state current driven by finite reservoirs with proportional coupling to the system/junction. In doing so, we present a simplified and unified derivation of the non-interacting and many-body steady-state currents through arbitrary junctions, including outside of proportional coupling. We conjecture that the analytic solution for proportional coupling is the most general of its form for isomodal relaxation (i.e., relaxing proportional coupling will remove the ability to find compact, general analytical expressions for finite reservoirs). These results should be of broad utility in diagnosing the behavior and implementation of extended reservoir and related approaches, including the convergence to the Landauer limit (for non-interacting systems) and the Meir-Wingreen formula (for many-body systems).

Excitation Dynamics in Chain-Mapped Environments

The chain mapping of structured environments is a most powerful tool for the simulation of open quantum system dynamics. Once the environmental bosonic or fermionic degrees of freedom are unitarily rearranged into a one dimensional structure, the full power of Density Matrix Renormalization Group (DMRG) can be exploited. Beside resulting in efficient and numerically exact simulations of open quantum systems dynamics, chain mapping provides an unique perspective on the environment: the interaction between the system and the environment creates perturbations that travel along the one dimensional environment at a finite speed, thus providing a natural notion of light-, or causal-, cone. In this work we investigate the transport of excitations in a chain-mapped bosonic environment. In particular, we explore the relation between the environmental spectral density shape, parameters and temperature, and the dynamics of excitations along the corresponding linear chains of quantum harmonic oscillators. Our analysis unveils fundamental features of the environment evolution, such as localization, percolation and the onset of stationary currents.

Evolution Equations for Quantum Semi-Markov Dynamics

Using a newly introduced connection between the local and non-local description of open quantum system dynamics, we investigate the relationship between these two characterisations in the case of quantum semi-Markov processes. This class of quantum evolutions, which is a direct generalisation of the corresponding classical concept, guarantees mathematically well-defined master equations, while accounting for a wide range of phenomena, possibly in the non-Markovian regime. In particular, we analyse the emergence of a dephasing term when moving from one type of master equation to the other, by means of several examples. We also investigate the corresponding Redfield-like approximated dynamics, which are obtained after a coarse graining in time. Relying on general properties of the associated classical random process, we conclude that such an approximation always leads to a Markovian evolution for the considered class of dynamics.

Steady-state quantum transport through an anharmonic oscillator strongly coupled to two heat reservoirs

We investigate the transport properties of an anharmonic oscillator, modeled by a single-site Bose-Hubbard model, coupled to two different thermal baths using the numerically exact thermofield based chain-mapping matrix product states (TCMPS) approach. We compare the effectiveness of TCMPS to probe the nonequilibrium dynamics of strongly interacting system irrespective of the system-bath coupling against the global master equation approach in Gorini-Kossakowski-Sudarshan-Lindblad form. We discuss the effect of on-site interactions, temperature bias as well as the system-bath couplings on the steady-state transport properties. Last, we also show evidence of non-Markovian dynamics by studying the nonmonotonicity of the time evolution of the trace distance between two different initial states.

Rate Operator Unraveling for Open Quantum System Dynamics

Stochastic methods with quantum jumps are often used to solve open quantum system dynamics. Moreover, they provide insight into fundamental topics, such as the role of measurements in quantum mechanics and the description of non-Markovian memory effects. However, there is no unified framework to use quantum jumps to describe open-system dynamics in any regime. We solve this issue by developing the rate operator quantum jump (ROQJ) approach. The method not only applies to both Markovian and non-Markovian evolutions, but also allows us to unravel master equations for which previous methods do not work. In addition, ROQJ yields a rigorous measurement-scheme interpretation for a wide class of dynamics, including a set of master equations with negative decay rates, and sheds light on different types of memory effects which arise when using stochastic quantum jump methods.

Modelling the ultra-strongly coupled spin-boson model with unphysical modes

A quantum system weakly coupled to a zero-temperature environment will relax, via spontaneous emission, to its ground-state. However, when the coupling to the environment is ultra-strong the ground-state is expected to become dressed with virtual excitations. This regime is difficult to capture with some traditional methods because of the explosion in the number of Matsubara frequencies, i.e., exponential terms in the free-bath correlation function. To access this regime we generalize both the hierarchical equations of motion and pseudomode methods, taking into account this explosion using only a biexponential fitting function. We compare these methods to the reaction coordinate mapping, which helps show how these sometimes neglected Matsubara terms are important to regulate detailed balance and prevent the unphysical emission of virtual excitations. For the pseudomode method, we present a general proof of validity for the use of superficially unphysical Matsubara-modes, which mirror the mathematical essence of the Matsubara frequencies.