An official website of the United States government
We investigate the quantum dynamics of a spin coupling to a bath of independent spins via the dissipaton equation of motion (DEOM) approach. The bath, characterized by a continuous spectral density function, is composed of spins that are independent level systems described by the su(2) Lie algebra, representing an environment with a large magnitude of anharmonicity. Based on the previous work by Suarez and Silbey [J. Chem. Phys. 95, 9115 (1991)] and by Makri [J. Chem. Phys. 111, 6164 (1999)] that the spin bath can be mapped to a Gaussian environment under its linear response limit, we use the time-domain Prony fitting decomposition scheme to the bare–bath time correlation function (TCF) given by the bosonic fluctuation–dissipation theorem to generate the exponential decay basis (or pseudo modes) for DEOM construction. The accuracy and efficiency of this strategy have been explored by a variety of numerical results. We envision that this work provides new insights into extending the hierarchical equations of motion and DEOM approach to certain types of anharmonic environments with arbitrary TCF or spectral density.
more »
« less
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.
more »
« less
- PAR ID:
- 10516892
- Publisher / Repository:
- AIP
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 20
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima–Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.more » « less
-
Abstract We unravel the correlated quantum quench dynamics of a single impurity immersed in a bosonic environment confined in an one-dimensional double-well potential. A particular emphasis is placed on the structure of the time-evolved many-body (MB) wave function by relying on a Schmidt decomposition whose coefficients directly quantify the number of configurations that are macroscopically populated. For a non-interacting bosonic bath and weak postquench impurity-bath interactions, we observe the dynamical formation of a two-fold fragmented MB state which is related to intra-band excitation processes of the impurity and manifests as a two-body phase separation (clustering) between the two species for repulsive (attractive) interactions. Increasing the postquench impurity-bath coupling strength leads to the destruction of the two-fold fragmentation since the impurity undergoes additional inter-band excitation dynamics. By contrast, a weakly interacting bath suppresses excitations of the bath particles and consequently the system attains a weakly fragmented MB state. Our results explicate the interplay of intra- and inter-band impurity excitations for the dynamical generation of fragmented MB states in multi-well traps and for designing specific entangled impurity states.more » « less
-
We present a general approach to derive Lindblad master equations for a subsystem whose dynamics is coupled to dissipative bosonic modes. The derivation relies on a Schrieffer-Wolff transformation which allows to eliminate the bosonic degrees of freedom after self-consistently determining their state as a function of the coupled quantum system. We apply this formalism to the dissipative Dicke model and derive a Lindblad master equation for the atomic spins, which includes the coherent and dissipative interactions mediated by the bosonic mode. This master equation accurately predicts the Dicke phase transition and gives the correct steady state. In addition, we compare the dynamics using exact diagonalization and numerical integration of the master equation with the predictions of semiclassical trajectories. We finally test the performance of our formalism by studying the relaxation of a NOON state and show that the dynamics captures quantum metastability beyond the mean-field approximation.more » « less
-
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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0198567
