Model order reduction encompasses mathematical techniques aimed at reducing the complexity of mathematical models in simulations while retaining the essential characteristics and behaviors of the original model. This is particularly useful in the context of large-scale dynamical systems, which can be computationally expensive to analyze and simulate. Here, we present a model order reduction technique to reduce the time complexity of open quantum systems, grounded in the principle of measurement-adapted coarse graining. This method, governed by a coarse-graining timescale τ and the spectral band center ω 0 of the measurement channel, organizes corrections to the lowest-order model which aligns with the rotating-wave approximation Hamiltonian in certain limits and rigorously justifies the resulting effective quantum master equation (EQME). The focus on calculating to a high degree of accuracy what can only be resolved by the measurement introduces a principled regularization procedure to address singularities and generates low-stiffness models suitable for efficient long-time integration. The closed-form EQME parameters greatly enhance interpretability, while the predicted EQME structure provides a foundation for deriving stochastic master equations under continuous measurement and for constructing efficient gray-box models for system identification and control. As a demonstration, we derive the fourth-order EQME for a challenging problem related to the dynamics of a superconducting qubit under high-power dispersive readout in the presence of a continuum of dissipative waveguide modes. This derivation shows that the lowest-order terms align with previous results, while higher-order corrections suggest new phenomena.
Model order reduction encompasses mathematical techniques aimed at reducing the complexity of mathematical models in simulations while retaining the essential characteristics and behaviors of the original model. This is particularly useful in the context of large-scale dynamical systems, which can be computationally expensive to analyze and simulate. Here, we present a model order reduction technique to reduce the time complexity of open quantum systems, grounded in the principle of measurement-adapted coarse graining. This method, governed by a coarse-graining timescale τ and the spectral band center ω 0 of the measurement channel, organizes corrections to the lowest-order model which aligns with the rotating-wave approximation Hamiltonian in certain limits and rigorously justifies the resulting effective quantum master equation (EQME). The focus on calculating to a high degree of accuracy what can only be resolved by the measurement introduces a principled regularization procedure to address singularities and generates low-stiffness models suitable for efficient long-time integration. The closed-form EQME parameters greatly enhance interpretability, while the predicted EQME structure provides a foundation for deriving stochastic master equations under continuous measurement and for constructing efficient gray-box models for system identification and control. As a demonstration, we derive the fourth-order EQME for a challenging problem related to the dynamics of a superconducting qubit under high-power dispersive readout in the presence of a continuum of dissipative waveguide modes. This derivation shows that the lowest-order terms align with previous results, while higher-order corrections suggest new phenomena.
Real-world experiments on quantum systems are always performed with measurement apparatus whose interaction times with the systems are finite. This restricts the observable quantum states to the space of time-coarse-grained density matrices ${\overline{\ensuremath{\rho}}(t)}$, providing motivation for the time-coarse-graining (TCG) approach that solves for $\overline{\ensuremath{\rho}}(t)$ without ever referring to the unobservable $\ensuremath{\rho}(t)$ of infinite time resolution. Phenomenologically, this implies that coherent transitions far outside the bandwidth would be filtered out, leaving only their effective impacts on the ``slow'' dynamics resolvable by the finite time resolution of the measurements. Therefore, the TCG framework provides rigorous justification for many existing effective Hamiltonian methods in the literature that aim at capturing the unitary part of the long-time dynamics. However, since time-coarse-graining is fundamentally irreversible, the resulting effective model allows for secular loss of information and dissipation of energy in general, which cannot be captured by any unitary effective models and has to be treated with explicit time-coarse-graining. Such incoherent effects are particularly prominent in driven nonlinear quantum systems where exchange of information and energy with the drive gives rise to incoherent effective dynamics at all finite time resolutions in general. While rigorous in principle, existing TCG methods in the literature can be applied only to simple systems by one iteratively solving superoperator equations at low orders in the coupling strengths. The complexity of such methods prevents systematic study of the time-coarse-grained dynamics and limits analytical results to the IR (low-resolution) limit in most cases. We address these limitations in this paper, presenting a systematic time-coarse-graining method that overcomes the complexities and restrictions of current techniques, offering a comprehensive and accurate modeling framework for driven nonlinear quantum systems. We derive closed-form formulas as well as diagrammatic representations for both unitary and nonunitary contributions, in the form of an effective Hamiltonian and nonunitary dissipators at arbitrary order in the coupling strengths, and complement this with a computer-algebra software package. We demonstrate the effectiveness of the method using several typical models of driven nonlinear systems in superconducting circuits, and show that it generalizes and improves on existing methods by providing more accurate results and explaining phenomena that have not been accounted for.
Many recent advancements in quantum computing leverage strong drives on nonlinear systems for state preparation, signal amplification, or gate operation. However, the interplay within such strongly driven system introduces multi-scale dynamics that affects the long-time behavior of the system in non-trivial ways that are very difficult to model. Therefore, the analysis of these systems often relies on effective Hamiltonian models that introduce additional nonlinear processes which approximate the long-time dynamics so that highly oscillatory terms may be ignored. However, the removal of such high frequency transitions can only be performed rigorously within a systematic framework of time-coarse graining, which is a fundamentally irreversible operation. This implies that standard approaches with unitary effective models cannot accurately capture the long-time behavior of strongly driven nonlinear quantum systems in general. We introduce a systematic perturbation theory for obtaining the complete non-unitary effective model of the time-coarse grained (TCG) dynamics of a driven quantum system to any order in the coupling strengths. We derive a closed-form analytical formula for both unitary and non-unitary contributions, in the form of an effective Hamiltonian and non-unitary (pseudo-)dissipators. Remarkably, even though the effective theory presumes unitary time evolution at the microscopic level, the time-coarse grained dynamics is found to follow a non-unitary time evolution in general. This occurs even when there is no open heat reservoir for the system to become entangled with or dissipate into. We demonstrate the effectiveness of the new method using several typical models of driven nonlinear systems in superconducting circuits, and show that it generalizes and improves on existing methods by providing more accurate results and explaining phenomena that have not been accounted for.
Driven-dissipative quantum devices are routinely operated with time-dependent drives and exposed to temporally correlated classical noise. When these perturbations vary on timescales comparable to the system's intrinsic dynamics, rotating-wave and Markovian approximations can fail, making effective static descriptions unreliable. As a result, accurately modeling systems in this regime often relies on simulating the full time-dependent master equation. This is particularly demanding for processes requiring long integration times, such as rare switching between metastable states, or for stochastic perturbations requiring averaging over many noise realizations. To address these challenges, we introduce Liouvillian spectral snapshots (LSS), a reduced-basis method that enables efficient simulation of such systems by projecting the master equation onto a low-dimensional subspace capturing their coarse-grained dynamics. LSS constructs this subspace by computing the slowly decaying, low-frequency eigenoperators of the Liouvillian at a few adaptively sampled perturbation values and compressing them into a common basis. This basis is built once and can be reused across noise realizations and perturbation parameters, reducing the cost of ensemble averages and parameter sweeps.
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.
Several approximation procedures, such as the full or partial rotating-wave, time-averaging, and geometric-arithmetic approximations, have been proposed to derive Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) generators from the Born-Markov quantum master equation (e.g., the Redfield equation). Establishing rigorous error bounds for these approximations is of fundamental and practical importance. However, existing bounds face two major limitations: they are highly specific to individual methods, and, more critically, they diverge in the long-time limit, ensuring the accuracy of the derived GKSL generator only in short-time regimes. In this Letter, we resolve both issues by deriving a unified, rigorous error bound for a general class of approximation methods -- termed temporal coarse graining -- that encompasses all aforementioned schemes. Crucially, our error bound is time-uniform. This guarantees that GKSL generators obtained via temporal coarse graining remain accurate for arbitrarily long times, provided the dissipation timescale is significantly longer than the bath correlation timescale.
Understanding how classical physics emerges from quantum mechanics remains a central problem in the foundations of physics. Here we derive a classical limit from finite-resolution measurements, modeled by continuous coarse-grained POVMs. When the resolved phase-space area is large compared with Planck's constant, the accessible statistics of any quantum state admit an effective classical description: coarse-grained observables become approximately jointly measurable and the induced probability density is positive. We derive the exact evolution equation for this density and show that, in the strong coarse-graining regime, its non-Liouville corrections are suppressed up to an Ehrenfest time, resulting in classical Hamiltonian flow generated by a Hamiltonian smoothed over the measurement cell. When the Hamiltonian varies negligibly across such a cell, the smoothed Hamiltonian reduces to the classical Hamiltonian whose quantization produced the quantum dynamics, thereby closing the quantization--classical-limit loop. Repeated finite-resolution measurements then generate stochastic records confined, with high probability, to tubes around classical trajectories. Our results provide a unified operational framework for the quantum-to-classical transition in microscopic and macroscopic systems, and establish the consistency of the quantization--classical-limit cycle.