Miralem Sinanović,Alessandro Ciani,Shai Machnes,Frank K. Wilhelm
Abstract:Open quantum systems are ubiquitous in nature and central to quantum technologies. A common description of their dynamics is given by the celebrated Lindblad master equation, which can be generalized to the non-Markovian scenario. In this work, we introduce the Non-Markovian Ensemble Propagation (NMEP) method, which extends the Monte Carlo Wave-Function (MCWF) method to the non-Markovian case in a simple and general manner. We demonstrate its accuracy and effectiveness in a selection of examples, and compare the results with either analytic expressions or direct numerical integration of the master equation.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to effectively simulate the dynamic evolution of the system in non - Markovian open quantum systems. Specifically, the paper introduces the Non - Markovian Ensemble Propagation (NMEP) method, which is a new technique that extends the Monte Carlo Wave Function (MCWF) method and can handle dynamic problems under non - Markovian conditions, especially when the jump rates can take negative values.
### Background and Problem Description of the Paper
Open quantum systems are ubiquitous in nature and play a central role in quantum technology. Usually, the dynamics of these systems can be described by the well - known Lindblad master equation, which can be generalized to non - Markovian situations. However, the traditional MCWF method can only handle problems under Markovian conditions, that is, all jump rates must be non - negative. For many practical applications, especially when the environment is not a Markov process, this assumption is often not valid. For example, superconducting qubits need to use a non - Markovian description when affected by 1/f noise or other pink noise.
### Main Contributions of the Paper
The main contribution of the paper is the proposal of the NMEP method, which can effectively simulate the dynamic evolution of the system under non - Markovian conditions by introducing the concept of negative jump rates. Specifically, the NMEP method:
1. **Extends the MCWF method**: The NMEP method can handle not only positive jump rates but also negative jump rates, thus being applicable to a wider range of physical scenarios.
2. **Maintains computational efficiency**: Despite the introduction of negative jump rates, the NMEP method still maintains high computational efficiency, especially when dealing with large - scale Hilbert spaces.
3. **Provides a theoretical basis**: The paper elaborately derives the mathematical basis of the NMEP method and proves that it satisfies the non - Markovian master equation on average.
### Specific Implementation
The specific implementation steps of the NMEP method are as follows:
1. **Deterministic Trajectories and Jump Paths**: Starting from each initial state \( |\psi_\alpha \rangle \), calculate the deterministic trajectory \( |\psi^{(\alpha)}_{\beta_0} \rangle \) and the jump path \( |\psi^{(\alpha)}_{\beta_l} \rangle \):
\[
|\psi^{(\alpha)}_{\beta_0} \rangle = \frac{(1 - i H_{\text{eff}}(t) \delta t) |\psi_\alpha \rangle}{\| (1 - i H_{\text{eff}}(t) \delta t) |\psi_\alpha \rangle \|},
\]
\[
|\psi^{(\alpha)}_{\beta_l} \rangle = \frac{A_l(t) |\psi_\alpha \rangle}{\| A_l(t) |\psi_\alpha \rangle \|}.
\]
2. **Random Sampling**: Define the probability \( P^{(\alpha)}_l \) of each state \( |\psi_\alpha \rangle \) to jump:
\[
P^{(\alpha)}_l = \delta t | \gamma_l(t) | \| A_l(t) |\psi_\alpha \rangle \|^2.
\]
Then, randomly sample the occupation number \( N^{(\alpha)}_{\beta_l} \) of each jump path:
\[
N^{(\alpha)}_{\beta_l} \sim \text{sgn}(N_\alpha \gamma_l(t))