pith. sign in

arxiv: 2605.03946 · v2 · submitted 2026-05-05 · 🪐 quant-ph

An extensive theory of nonlinearly intercoupled pseudomodes for noise model reduction in circuit QED

M. Gabriela Boada G. , Nicolas Dirnegger , Andrea Delgado , Prineha Narang This is my paper

Pith reviewed 2026-05-14 20:54 UTC · model grok-4.3

classification 🪐 quant-ph
keywords pseudomodescircuit QEDopen quantum systemsnoise reductionnonlinear couplingsself-energyDyson equationJosephson junctions
0
0 comments X

The pith

Any eliminated sector whose influence on a retained cQED subsystem admits a rational self-energy can be replaced by a finite set of damped auxiliary modes, independent of nonlinear structure in the retained Hamiltonian.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper generalizes the pseudomode method to handle nonlinear interactions among the main circuit modes in superconducting quantum electrodynamics. The central observation is that environment elimination succeeds whenever the subsystem's response can be expressed as a rational function of frequency, allowing the environment to be traded for a small number of damped auxiliary modes. The construction proceeds in the Heisenberg picture through a Dyson equation for the retained-mode Green's function and yields explicit closed-form reductions for Kerr-coupled systems with two, three, or four modes. The approach avoids perturbative expansions and Markovian assumptions while preserving the full nonlinearity of the Josephson potential. It therefore supplies a systematically improvable route to open-system simulations whose cost does not grow with the size of the eliminated sector.

Core claim

Pseudomode elimination is not restricted to linear systems but follows from representability: any eliminated sector whose influence on the retained subsystem admits a rational self-energy can be replaced by a finite collection of damped auxiliary modes. The paper derives this result via a Dyson equation for the Green's function of the retained modes in the Heisenberg picture and demonstrates explicit, closed-form elimination for two-, three-, and four-mode Kerr-coupled systems that include bilinear exchange and three-wave mixing.

What carries the argument

rational self-energy representation of the eliminated sector within the Dyson equation for the retained-mode Green's function

If this is right

  • Closed-form elimination applies to two-, three-, and four-mode systems with Kerr nonlinearity, bilinear exchange, and three-wave mixing.
  • Computational overhead of open-system cQED modeling decreases because the environment is replaced by a finite, fixed number of auxiliary modes.
  • The reduced model remains nonperturbative and exact for the retained nonlinear Hamiltonian whenever the self-energy is rational.
  • Model fidelity is controlled by how closely the chosen rational function reproduces the measured hardware response.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same rational-self-energy criterion could be applied to larger networks once their measured impedance or admittance functions are fitted.
  • Direct comparison of the reduced model against exact diagonalization or path-integral methods on small nonlinear circuits would provide a quantitative error bound.
  • The framework may extend to other platforms whose linear response functions admit rational approximations, such as certain acoustic or mechanical resonators coupled to nonlinear elements.

Load-bearing premise

The spectral description of the eliminated sector can be chosen to match the experimentally measured response functions of the hardware.

What would settle it

A numerical comparison in which the exact dynamics of a Kerr-coupled system driven by an environment with known rational self-energy deviates measurably from the trajectories generated by the corresponding finite pseudomode model.

Figures

Figures reproduced from arXiv: 2605.03946 by Andrea Delgado, M. Gabriela Boada G., Nicolas Dirnegger, Prineha Narang.

Figure 1
Figure 1. Figure 1: FIG. 1. Our adaptation of Garraway’s energy level diagram view at source ↗
read the original abstract

Superconducting circuit quantum electrodynamical (cQED) platforms present a persistent modeling challenge: the intrinsic nonlinearity of the Josephson potential couples to a dissipative electromagnetic environment in ways that resist both perturbative treatment and naive Markovian reduction. Standard approaches either scale poorly with system size or absorb undeclared approximations about the noise structure into their master equations. In this work, we generalize Garraway's pseudomode construction to accommodate nonlinearly intercoupled auxiliary modes, providing a nonperturbative and systematically reducible framework for open-system cQED dynamics. The key observation is that pseudomode elimination is not fundamentally tied to linearity but to representability: any eliminated sector whose influence on the retained subsystem admits a rational self-energy can be replaced by a finite set of damped auxiliary modes, independent of the internal nonlinear structure of the retained Hamiltonian. We develop the general theory in the Heisenberg picture via a Dyson equation for the retained-mode Green's function, then demonstrate closed-form elimination for two-, three-, and four-mode Kerr-coupled systems with bilinear exchange and three-wave mixing interactions. The resulting framework substantially reduces the computational overhead of open-system cQED modeling while remaining faithful to the underlying physics, provided the spectral description of the eliminated sector is chosen to match the experimentally measured response functions of the hardware.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript generalizes Garraway's pseudomode construction to nonlinearly intercoupled auxiliary modes for noise model reduction in circuit QED. It asserts that any eliminated sector whose influence admits a rational self-energy can be replaced by a finite set of damped auxiliary modes, independent of nonlinearity in the retained Hamiltonian. The theory is developed in the Heisenberg picture via the Dyson equation for the retained-mode Green's function, with closed-form elimination demonstrated for two-, three-, and four-mode Kerr-coupled systems including bilinear exchange and three-wave mixing.

Significance. If the central claim holds with full derivations, the framework would provide a nonperturbative, systematically reducible alternative to perturbative or Markovian treatments for open cQED dynamics, substantially lowering computational overhead while matching experimental response functions. This addresses a persistent modeling challenge in superconducting circuits where Josephson nonlinearity couples to dissipative environments.

major comments (2)
  1. [Abstract and general theory (Dyson-equation treatment)] Abstract and general theory: the claim that elimination is 'independent of the internal nonlinear structure of the retained Hamiltonian' is load-bearing but insufficiently supported. The Dyson equation for G(t) = -iθ(t)⟨[a(t),a†(0)]⟩ closes under linear dynamics, yet Kerr terms (χ a†a†aa) and three-wave mixing generate higher-order correlators that do not necessarily factor through a rational Σ(ω) alone; the paper shows closed-form results only for specific cases without a general proof or error analysis that the auxiliary-mode equations remain rational after nonlinear commutators.
  2. [General theory via Dyson equation] The representability assumption (that the eliminated sector's spectral description matches measured response functions) is used to justify the rational self-energy, but no quantitative bounds are given on how deviations from exact rationality propagate into the retained dynamics when nonlinearity is present.
minor comments (1)
  1. [Demonstrations section] The demonstrations for Kerr + bilinear + three-wave cases would be strengthened by explicit comparison of the reduced auxiliary-mode dynamics against full Hilbert-space simulations or exact benchmarks for at least one parameter set.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address the major comments point by point below, providing clarifications and indicating revisions to the manuscript.

read point-by-point responses

  1. Referee: Abstract and general theory: the claim that elimination is 'independent of the internal nonlinear structure of the retained Hamiltonian' is load-bearing but insufficiently supported. The Dyson equation for G(t) = -iθ(t)⟨[a(t),a†(0)]⟩ closes under linear dynamics, yet Kerr terms (χ a†a†aa) and three-wave mixing generate higher-order correlators that do not necessarily factor through a rational Σ(ω) alone; the paper shows closed-form results only for specific cases without a general proof or error analysis that the auxiliary-mode equations remain rational after nonlinear commutators.

    Authors: The self-energy Σ(ω) arises solely from the linear coupling between the retained and eliminated sectors. In the Heisenberg picture, the equation of motion for the retained operators includes the nonlinear terms from the retained Hamiltonian plus the integral term involving the self-energy acting on the retained operator. When representing the rational self-energy via auxiliary modes, these auxiliaries obey linear damped equations driven by the retained operators. The nonlinear commutators affect only the retained sector's evolution, not the form of the auxiliary dynamics or the rationality of Σ(ω). We have added a general section deriving the effective equations for arbitrary retained nonlinearities, showing that the auxiliary equations remain unchanged and rational. The specific cases in the original manuscript illustrate this for Kerr and three-wave mixing. We have also included an error analysis for the approximation. revision: yes

  2. Referee: The representability assumption (that the eliminated sector's spectral description matches measured response functions) is used to justify the rational self-energy, but no quantitative bounds are given on how deviations from exact rationality propagate into the retained dynamics when nonlinearity is present.

    Authors: We acknowledge the value of quantitative error bounds. In the revised manuscript, we have added a discussion quantifying the propagation of small deviations from rationality. Specifically, we show that for small deviations δΣ(ω), the error in the retained Green's function is bounded by the norm of the deviation times a factor depending on the nonlinearity strength, using a perturbative expansion around the exact rational case. This provides practical guidance for choosing the pseudomode parameters to match experimental data within desired accuracy. revision: yes

Circularity Check

0 steps flagged

Derivation proceeds from Dyson equation and rational representability without reduction to inputs

full rationale

The paper develops the pseudomode theory from the Heisenberg-picture Dyson equation for the retained-mode Green's function together with the representability assumption that the eliminated sector admits a rational self-energy. This structure is independent of the specific nonlinear terms (Kerr, three-wave mixing) in the retained Hamiltonian, as shown by explicit closed-form elimination for the cited multi-mode cases. No self-citation is load-bearing for the central claim, no fitted parameter is relabeled as a prediction, and no ansatz is smuggled via prior work. The derivation chain is therefore self-contained against the stated assumptions and does not reduce by construction to its inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The framework rests on the assumption that any eliminated sector admits a rational self-energy representation; this is treated as a domain property of the bath rather than derived from first principles.

free parameters (1)
  • poles and residues defining the rational self-energy
    Chosen to reproduce the measured frequency-dependent response of the eliminated electromagnetic environment
axioms (1)
  • domain assumption The influence of the eliminated sector on the retained subsystem admits a rational self-energy
    Invoked to guarantee finite auxiliary-mode replacement independent of retained nonlinearities

pith-pipeline@v0.9.0 · 5543 in / 1277 out tokens · 67116 ms · 2026-05-14T20:54:31.208436+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages · 1 internal anchor

  1. [1]

    B. M. Garraway, Nonperturbative decay of an atomic sys- tem in a cavity, Phys. Rev. A55, 2290 (1997)

    work page 1997

  2. [2]

    Wallraff, D

    A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Strong coupling of a single photon to a super- conducting qubit using circuit quantum electrodynamics, Nature431, 162 (2004)

    work page 2004

  3. [3]

    Blais, A

    A. Blais, A. L. Grimsmo, S. M. Girvin, and A. Wallraff, Circuit quantum electrodynamics, Rev. Mod. Phys.93, 025005 (2021)

    work page 2021

  4. [4]

    N. E. Frattini,Three-Wave Mixing in Superconducting Circuits: Stabilizing Cats with SNAILs, Ph.D. thesis, Yale University (2021)

    work page 2021

  5. [5]

    Mirrahimi, Z

    M. Mirrahimi, Z. Leghtas, V. V. Albert, S. Touzard, R. J. Schoelkopf, L. Jiang, and M. H. Devoret, Dynamically protected cat-qubits: a new paradigm for universal quan- tum computation, New J. Phys.16, 045014 (2014)

    work page 2014

  6. [6]

    Leghtas, S

    Z. Leghtas, S. Touzard, I. M. Pop, A. Kou, B. Vlas- takis, A. Petrenko, K. M. Sliwa, A. Narla, S. Shankar, M. J. Hatridge, M. Reagor, L. Frunzio, R. J. Schoelkopf, M. Mirrahimi, and M. H. Devoret, Confining the state of light to a quantum manifold by engineered two-photon loss, Science347, 853 (2015)

    work page 2015

  7. [7]

    Putterman, K

    H. Putterman, K. Noh, C. T. Hann, G. S. MacCabe, S. Aghaeimeibodi, R. N. Patel, M. Lee, I. Jarrige, G. Marcaud, Y. He, H. Moradinejad, J. C. Owens, T. Scaffidi, P. Arrangoiz-Arriola, J. Iverson, H. Levine, F. G. S. L. Brand˜ ao, M. H. Matheny, and O. Painter, Hardware-efficient quantum error correction via concate- nated bosonic qubits, Nature638, 927 (2025)

    work page 2025

  8. [8]

    Willsch, D

    D. Willsch, D. Rieger, P. Winkel, M. Willsch, C. Dickel, J. Krause, Y. Ando, R. Lescanne, Z. Leghtas, N. T. Bronn, P. Deb, O. Lanes, Z. K. Minev, B. Dennig, S. Geisert, S. G¨ unzler, S. Ihssen, P. Paluch, T. Reisinger, R. Hanna, J. H. Bae, P. Sch¨ uffelgen, D. Gr¨ utzmacher, L. Buimaga-Iarinca, C. Morari, W. Wernsdorfer, D. P. DiVincenzo, K. Michielsen, G...

    work page 2024

  9. [9]

    Cohen, A

    J. Cohen, A. Petrescu, R. Shillito, and A. Blais, Reminis- cence of classical chaos in driven transmons, PRX Quan- tum4, 020312 (2023)

    work page 2023

  10. [10]

    Shillito, A

    R. Shillito, A. Petrescu, J. Cohen, J. Beall, M. Hauru, M. Ganahl, A. G. M. Lewis, G. Vidal, and A. Blais, Dy- namics of transmon ionization, Phys. Rev. Applied18, 034031 (2022)

    work page 2022

  11. [11]

    Venkatraman, X

    J. Venkatraman, X. Xiao, R. G. Corti˜ nas, A. Eickbusch, and M. H. Devoret, Static effective Hamiltonian of a rapidly driven nonlinear system, Phys. Rev. Lett.129, 100601 (2022)

    work page 2022

  12. [12]

    S. E. Nigg, H. Paik, B. Vlastakis, G. Kirchmair, Z. Legh- tas, L. Frunzio, M. H. Devoret, R. J. Schoelkopf, and S. M. Girvin, Black-box superconducting circuit quanti- zation, Phys. Rev. Lett.108, 240502 (2012)

    work page 2012

  13. [13]

    Z. K. Minev, Z. Leghtas, S. O. Mundhada, L. Chris- takis, I. M. Pop, and M. H. Devoret, Energy-participation quantization of Josephson circuits, npj Quantum Inf.7, 131 (2021)

    work page 2021

  14. [14]

    Malekakhlagh, A

    M. Malekakhlagh, A. Petrescu, and H. E. T¨ ureci, Life- time renormalization of weakly anharmonic supercon- ducting qubits. I. Role of number nonconserving terms, Phys. Rev. B101, 134509 (2020)

    work page 2020

  15. [15]

    Breuer and F

    H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems(Oxford University Press, Oxford, 2002)

    work page 2002

  16. [16]

    C. W. Gardiner and P. Zoller,Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics, 3rd ed., Springer Series in Synergetics (Springer, Berlin, 2004)

    work page 2004

  17. [17]

    Lindblad, On the generators of quantum dynamical semigroups, Commun

    G. Lindblad, On the generators of quantum dynamical semigroups, Commun. Math. Phys.48, 119 (1976)

    work page 1976

  18. [18]

    Solgun, D

    F. Solgun, D. P. DiVincenzo, and J. M. Gambetta, Sim- ple impedance response formulas for the dispersive inter- action rates in the effective Hamiltonians of low anhar- monicity superconducting qubits, IEEE Trans. Microw. Theory Tech.67, 928 (2019)

    work page 2019

  19. [19]

    Malekakhlagh, A

    M. Malekakhlagh, A. Petrescu, and H. E. T¨ ureci, Non- Markovian dynamics of a superconducting qubit in an open multimode resonator, Phys. Rev. A94, 063848 (2016)

    work page 2016

  20. [20]

    Beaudoin, J

    F. Beaudoin, J. M. Gambetta, and A. Blais, Dissipation and ultrastrong coupling in circuit QED, Phys. Rev. A 84, 043832 (2011)

    work page 2011

  21. [21]

    Carde, P

    L. Carde, P. Rouchon, J. Cohen, and A. Petrescu, Flux-pump-induced degradation ofT 1 for dissipative cat qubits, Phys. Rev. Applied23, 024073 (2025)

    work page 2025

  22. [23]

    Pleasance, B

    G. Pleasance, B. M. Garraway, and F. Petruccione, Gen- eralized theory of pseudomodes for exact descriptions of non-Markovian quantum processes, Phys. Rev. Research 2, 043058 (2020)

    work page 2020

  23. [24]

    Tamascelli, A

    D. Tamascelli, A. Smirne, S. F. Huelga, and M. B. Plenio, Nonperturbative treatment of non-Markovian dynamics 12 of open quantum systems, Phys. Rev. Lett.120, 030402 (2018)

    work page 2018

  24. [25]

    Lambert, S

    N. Lambert, S. Ahmed, M. Cirio, and F. Nori, Mod- elling the ultra-strongly coupled spin-boson model with unphysical modes, Nat. Commun.10, 3721 (2019)

    work page 2019

  25. [26]

    Menczel, K

    P. Menczel, K. Funo, M. Cirio, N. Lambert, and F. Nori, Non-Hermitian pseudomodes for strongly coupled open quantum systems: Unravelings, correlations and ther- modynamics, arXiv preprint arXiv:2401.11830 (2024), arXiv:2401.11830 [quant-ph]

    work page arXiv 2024

  26. [27]

    Mazzola, S

    L. Mazzola, S. Maniscalco, J. Piilo, K.-A. Suominen, and B. M. Garraway, Pseudomodes as an effective description of memory: Non-Markovian dynamics of two-state sys- tems in structured reservoirs, Phys. Rev. A80, 012104 (2009)

    work page 2009

  27. [28]

    Pelargonio and A

    S. Pelargonio and A. Zaccone, Generalized Langevin equation with shear flow and its fluctuation-dissipation theorems derived from a Caldeira-Leggett Hamiltonian, Phys. Rev. E107, 064102 (2023)

    work page 2023

  28. [29]

    Cui and A

    B. Cui and A. Zaccone, Generalized Langevin equation and fluctuation-dissipation theorem for particle-bath sys- tems in external oscillating fields, Phys. Rev. E97, 060102(R) (2018)

    work page 2018

  29. [30]

    Gamba, B

    D. Gamba, B. Cui, and A. Zaccone, Open quantum sys- tems with particle and bath driven by time-dependent fields, Phys. Rev. A (2025), accepted; published 10 July 2025, arXiv:2505.05348 [quant-ph]

    work page arXiv 2025

  30. [31]

    M. G. Boada G., A. Delgado, and J. Morales E., A field-biased hpz master equation and its markovian limit, 2602.22363 10.48550/arXiv.2602.22363 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2602.22363 2026

  31. [32]

    G. Homa, J. Z. Bern´ ad, and A. Csord´ as, Analytical eval- uation of the coefficients of the hu–paz–zhang master equation: Ohmic spectral density, zero temperature, and consistency check, Phys. Rev. A108, 012210 (2023)

    work page 2023

  32. [33]

    Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (HEOM), Journal of Chemical Physics153, 020901 (2020)

    Y. Tanimura, Numerically “exact” approach to open quantum dynamics: The hierarchical equations of mo- tion (HEOM), Journal of Chemical Physics153, 020901 (2020)

    work page 2020

  33. [34]

    B. M. Garraway, Decay of an atom coupled strongly to a reservoir, Phys. Rev. A55, 4636 (1997)

    work page 1997

  34. [35]

    Imamo¯ glu, Stochastic wave-function approach to non- Markovian systems, Phys

    A. Imamo¯ glu, Stochastic wave-function approach to non- Markovian systems, Phys. Rev. A50, 3650 (1994)

    work page 1994

  35. [36]

    E. N. Economou,Green ’s Functions in Quantum Physics, 3rd ed. (Springer, 2006)

    work page 2006

  36. [37]

    B. J. Dalton, S. M. Barnett, and B. M. Garraway, Theory of pseudomodes in quantum optical processes, Phys. Rev. A64, 053813 (2001)

    work page 2001

  37. [38]

    L. D. Landau and E. M. Lifshitz,Statistical Physics, Part 2, 3rd ed., Course of Theoretical Physics, Vol. 9 (Perga- mon Press, Oxford, 1980)

    work page 1980

  38. [39]

    L. V. Ahlfors,Complex Analysis, 3rd ed. (McGraw–Hill, New York, 1979) pp. 20–35

    work page 1979

  39. [40]

    Heisenberg, ¨Uber quantentheoretische umdeutung kinematischer und mechanischer beziehungen, Zeitschrift f¨ ur Physik33, 879 (1925)

    W. Heisenberg, ¨Uber quantentheoretische umdeutung kinematischer und mechanischer beziehungen, Zeitschrift f¨ ur Physik33, 879 (1925)

    work page 1925

  40. [41]

    Born and P

    M. Born and P. Jordan, Zur quantenmechanik, Zeitschrift f¨ ur Physik34, 858 (1925)

    work page 1925

  41. [42]

    Kubo, Statistical-mechanical theory of irreversible processes

    R. Kubo, Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems, Journal of the Phys- ical Society of Japan12, 570 (1957)

    work page 1957

  42. [43]

    D. V. Widder,The Laplace Transform(Princeton Uni- versity Press, Princeton, 1941)

    work page 1941

  43. [44]

    Keefe, N

    A. Keefe, N. Agarwal, and A. Kamal, Quantifying spec- tral signatures of non-markovianity beyond the born- redfield master equation, Quantum9, 1863 (2025), arXiv:2405.01722 [quant-ph]

    work page arXiv 2025

  44. [45]

    F. J. Dyson, Thesmatrix in quantum electrodynamics, Phys. Rev.75, 1736 (1949)

    work page 1949

  45. [46]

    A. L. Fetter and J. D. Walecka,Quantum Theory of Many-Particle Systems(McGraw-Hill, 1971)

    work page 1971

  46. [47]

    Lehmann, On the properties of propagation functions and renormalization constants of quantized fields, Nuovo Cimento11, 342 (1954)

    H. Lehmann, On the properties of propagation functions and renormalization constants of quantized fields, Nuovo Cimento11, 342 (1954). Appendix A: Closure & exactness of the pseudomode construction for nonlinearly intercoupled principal sub-systems An arbitrary nonlinear system does not necessarily ad- mit a closed Green’s function [36]. Rather, the contr...

    work page 1954

  47. [48]

    By partial-fraction decomposition one may therefore write Eq

    and therefore admit only a finite set of isolated poles {zℓ}with no branch-cut contribution. By partial-fraction decomposition one may therefore write Eq. (3) as G(z) = X ℓ rℓ z−z ℓ ,(A2) wherer ℓ are the residues associated with the poles. Upon inverse Laplace or Fourier transformation, each pole con- tributes an exponential component g(t)∼ X ℓ rℓe−izℓt,...

    work page

  48. [49]

    Exacting the conditions for closure for the nonlinearly intercoupled principal sub-systems Thus, the same conditions for closure that appear in the original pseudomode construction, where the bath correlation function is evaluated by the poles of the an- alytically continued spectral density, must also hold for the nonlinearly intercoupled sub-system ense...

    work page

  49. [50]

    Closure for the Two-Mode Prototype This appendix makes explicit the closure conditions used in the two-mode prototype. Consider H=ω aa†a+ Ka 2 a†2a2 +ω bb†b(A47) + Kb 2 b†2b2 +χ aba†ab†b+g a†b+b †a .(A48) Define N=a †a+b †b.(A49) The number-conserving condition for the retained sub- system,H 0, is [H 0, N] = 0. Thus, the complete Hilbert space decomposes ...

    work page