A numerical method is proposed for simulation of composite open quantum systems. It is based on Lindblad master equations and adiabatic elimination. Each subsystem is assumed to converge exponentially towards a stationary subspace, slightly impacted by some decoherence channels and weakly coupled to the other subsystems. This numerical method is based on a perturbation analysis with an asymptotic expansion. It exploits the formulation of the slow dynamics with reduced dimension. It relies on the invariant operators of the local and nominal dissipative dynamics attached to each subsystem. Second-order expansion can be computed only with local numerical calculations. It avoids computations on the tensor-product Hilbert space attached to the full system. This numerical method is particularly well suited for autonomous quantum error correction schemes. Simulations of such reduced models agree with complete full model simulations for typical gates acting on one and two cat qubits (z, zz, and cnot) when the mean photon number of each cat qubit is less than eight. For larger mean photon numbers and gates with three cat qubits (zzz and ccnot), full model simulations are almost impossible whereas reduced model simulations remain accessible. In particular, one observes numerically the simultaneous capture of the dominant phase-flip error rate and of the very small bit-flip error rate with its exponential suppression versus the mean photon number.
Consider an open quantum system governed by a Gorini-Kossakowski-Sudarshan-Lindblad master equation with two times-scales: a fast one, exponentially converging towards a linear subspace of quasi-equilibria; a slow one resulting from small decoherence and Hamiltonian dynamics. Usually adiabatic elimination is performed in the Schrödinger picture. We propose here a Heisenberg formulation where the invariant operators attached to the fast decay dynamics towards the quasi-equilibria subspace play a key role. Based on geometric singular perturbations, asymptotic expansions of the Heisenberg slow dynamics and of the fast invariant linear subspaces are proposed. They exploit Carr's approximation lemma from center-manifold and bifurcation theory. Second-order expansions are detailed and shown to ensure preservation, up to second-order terms, of the complete positivity for the slow propagator on a slow timescale. Such expansions can be exploited numerically to derive reduced-order dynamical models.
This paper discusses quantum adiabatic elimination, which is a model reduction technique for a composite Lindblad system consisting of a fast decaying subsystem coupled to another subsystem with a much slower timescale. Such a system features an invariant manifold that is close to the slow subsystem. This invariant manifold is reached subsequent to the decay of the fast degrees of freedom, after which the slow dynamics follow on it. By parametrizing the invariant manifold, the slow dynamics can be simulated via a reduced model. To find the evolution of the reduced state, we perform an asymptotic expansion with respect to the timescale separation. So far, the second-order expansion has mostly been considered. It has then been revealed that the second-order expansion of the reduced dynamics is generally given by a Lindblad equation, which ensures complete positivity of the time evolution. In this paper, we present two examples where complete positivity of the reduced dynamics is violated with higher-order contributions. In the first example, the violation is detected for the evolution of the partial trace without truncation of the asymptotic expansion. The partial trace is not the only way to parametrize the slow dynamics. Concerning this nonuniqueness, it was conjectured in [Quantum Sci. Technol. 2, 044011 (2017)] that there exists a parameter choice ensuring complete positivity. With the second example, however, we refute this conjecture by showing that complete positivity cannot be restored in any choice of parametrization. We discuss these results in terms of the invariant slow manifold consisting of quantum correlated states.
Adiabatic elimination methods allow the reduction of the space dimension needed to describe systems whose dynamics exhibit separation of timescale. For open quantum systems, it consists in eliminating the fast part assuming it has almost instantaneously reached its steady state and obtaining an approximation of the evolution of the slow part. These methods can be applied to eliminate a linear subspace within the system Hilbert space or, alternatively, to eliminate a fast subsystem in a bipartite quantum system. In this work, we extend an adiabatic elimination method used for removing fast degrees of freedom within a system [Phys. Rev. A 101, 042102 (2020)] to eliminate a subsystem from an open bipartite quantum system. As an illustration, we apply our technique to a dispersively coupled two-qubit system and in the case of the open Rabi model.
We develop a mathematically rigorous framework for simulating \emph{multiscale physical systems} using quantum computational resources, by translating the \emph{language of asymptotic-preserving (AP) schemes} into the formalism of quantum channels and Lindbladian dynamics. For stiff open quantum systems governed by singularly perturbed generators $\cL_\eps = \eps^{-1}\cL_{\mathrm{fast}} + \cL_{\mathrm{slow}}$ with $\eps \to 0$, we prove that layered quantum protocols, which implement fast-scale relaxation via native analog evolution or analytic manifold projection, converge uniformly in the diamond norm to consistent discretizations of the limiting slow dynamics, with explicit error bound $\mathcal{O}(\epsΔt + Δt^2)$ independent of stiffness. We establish precise resource-complexity bounds showing that superlinear gate-count savings
Ω(κ⋅(dtot/dslow)c) arise if and only if fast dynamics are resolved via (i) hardware-native analog evolution, or (ii) analytic adiabatic elimination reducing effective Hilbert space dimension. The framework is illustrated through cavity QED in the bad-cavity limit and a quantum-inspired AP discretization of kinetic equations converging to fluid limits, with quantified error propagation in trace and diamond norms. This work provides a principled mathematical bridge between classical multiscale numerical analysis and quantum simulation algorithms.
This thesis develops model reduction methods for open quantum system dynamics, yielding efficient models for bipartite systems, with direct application to circuit QED and broader applicability beyond this platform. First, we consider dissipative dynamics characterized by timescale separation, which translates to a gap on the real axis of the spectrum of the Liouvillian generating it. This allows to compute the reduced evolution for a target subsystem by adiabatically eliminating fast-decaying degrees of freedom. A Sylvester equation integral representation yields closed forms for the reduced generators and explicit reconstruction maps to arbitrary order without a preliminary interaction-picture averaging step, providing a linear center manifold description that separates long-lived from fast-decaying modes by the real parts of the spectrum. The reduced dynamics takes place on the center manifold, that is, the spectral subspace associated with eigenvalues whose real parts are small (long-lived modes), while placing no restriction on imaginary parts, thus allowing for fast unitary dynamics on the target subsystem. We further demonstrate the equivalence of such a geometric approach to adiabatic elimination with the time-convolutionless master equation framework, by showing that the maps defining the reduced model satisfy the same invariance equation in both frameworks. The time-convolutionless formulation of adiabatic elimination extends to cases where the target system undergoes fast unitary dynamics and retains explicit time dependence in the reduced model, hence capturing initial transients that are neglected in time-independent geometric reductions. Second, we apply a unitary chain mapping of microscopic system–environment Hamiltonians to a nearest-neighbor form, enabling efficient time evolution exploiting tensor network time evolution algorithms based on the time dependent variational principle. Accuracy of the tensor network representation, and thus the size of the reduced model, is controlled by entanglement entropy growth, which we quantify numerically. While the methods are general and apply beyond the specific platform considered, we focus on dispersive readout in circuit QED (a qubit dipolarly coupled to a driven, lossy cavity). We quantify the dependence of the relaxation time
T1 of the qubit on intracavity photon number
nˉ (and thus on the drive amplitude) and on the spectral structure of the environment, and compute the emission spectrum directly from environmental observables. Simulations including a Purcell notch filter predict a decrease of
T1 with increasing drive power, a qualitatively different behavior with respect to Born–Markov–secular GKSL predictions but consistent with experimental trends, indicating that explicit filter modeling is required in this case. Together, the results specify when reduced Liouvillian models under a dissipative gap are reliable and when microscopic, non-perturbative simulations are required.
Dissipative cat qubits are known for the exponential suppression of the bit-flip rate. However, there is significant discrepancy between experimental measurements and analytical predictions of the strength of the bit-flip suppression. In this paper we resolve this discrepancy for some of the perturbations, by developing a second-order perturbation theory on top of a nonlinear dissipative Lindbladian. Following this scheme, we derive an analytical expression for the exponentially small bit-flip rate due to single-photon loss, which shows good agreement with numerical simulations. We also apply our scheme to other perturbations, such as frequency detuning and the
Z gate, and find the corresponding bit-flip rates, which also show good agreement with the numerical simulation.