First-principles study of dispersive readout in circuit QED
Pith reviewed 2026-05-10 14:55 UTC · model grok-4.3
The pith
Microscopic simulation of the transmission line bath explains why qubit T1 falls at high readout drive amplitudes in circuit QED.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By coupling the circuit QED Hamiltonian to a microscopic model of the measurement transmission line that can have arbitrary spectrum including filters, the full unitary dynamics of dispersive readout can be simulated from first principles. Access to the bath degrees of freedom shows that qubit T1 decreases with increasing readout drive amplitude specifically when a Purcell notch filter is placed at the qubit frequency. In contrast, the Lindblad master equation exhibits general qualitative defects and fails to capture this behavior.
What carries the argument
Microscopic model of the measurement transmission line bath coupled to the circuit QED system, enabling simulation of unitary dynamics and direct access to emission spectra and bath evolution.
If this is right
- The qubit relaxation time T1 becomes drive-amplitude dependent based on bath details like filters.
- A Purcell notch at qubit frequency causes T1 to drop as drive increases.
- The Lindblad master equation has qualitative shortcomings for this system.
- The emission spectrum of the system can be studied versus drive power.
- This explains observed saturation in readout fidelity at high amplitudes.
Where Pith is reading between the lines
- Designing transmission line filters to avoid notches at qubit frequencies could allow higher drive powers without T1 loss.
- Similar first-principles methods might apply to other decoherence channels or readout schemes in quantum circuits.
- Experimental measurement of the predicted power-dependent emission spectrum would test the model's accuracy.
- Accounting for bath spectrum details may be necessary for accurate modeling of high-power operations in superconducting quantum processors.
Load-bearing premise
The model of the measurement transmission line correctly represents the actual bath spectrum, couplings, and filter effects present in the physical device.
What would settle it
Observe whether qubit T1 decreases with readout drive amplitude only when a notch filter is tuned to the qubit frequency, and check if Lindblad predictions match or deviate from the observed dependence.
Figures
read the original abstract
The speed and fidelity of dispersive readout of superconducting qubits should improve by increasing the amplitude of the measurement drive. Experiments show, however, that beyond some drive amplitude there is always a saturation or drop in fidelity, often associated with a decrease in qubit energy relaxation time $T_1$. A simple Lindblad master equation does not capture the latter effect. More involved approaches based on effective master equations rely on strong assumptions about the spectra of the system and the bath and only partially agree with observations. Here, we perform a first-principles simulation of the full unitary dynamics of dispersive readout by considering the circuit QED Hamiltonian coupled to a microscopic model for the measurement transmission line, allowing for its arbitrary spectrum, including filters. Our access to the dynamics of the bath degrees of freedom allows us to investigate the emission spectrum of the system as a function of drive power. We show how the dependence of qubit $T_1$ on readout drive amplitude is sensitive to the details of the bath spectrum. In particular, we find that $T_1$ drops with increasing drive amplitude when a Purcell notch filter is placed at the qubit frequency, and that the Lindblad master equation shows general qualitative defects compared to the first-principles model.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a first-principles simulation of dispersive readout by evolving the full unitary circuit-QED Hamiltonian coupled to a microscopic model of the measurement transmission line that permits arbitrary bath spectra, including Purcell filters. It reports that qubit T1 depends sensitively on readout drive amplitude and on the details of the bath spectrum; specifically, T1 decreases with increasing drive when a notch filter is placed at the qubit frequency. The work also finds qualitative discrepancies between these dynamics and predictions from the standard Lindblad master equation, and uses access to bath degrees of freedom to examine the emitted spectrum versus drive power.
Significance. If the microscopic bath model is representative of real devices, the results would clarify drive-induced decoherence mechanisms that limit readout fidelity and would demonstrate concrete limitations of Lindblad treatments for strong-drive dispersive readout. The approach of retaining the full unitary dynamics plus an explicit bath is a strength, as is the ability to vary the bath spectrum arbitrarily and to inspect emission spectra directly.
major comments (2)
- [Bath model and numerical results sections] The central claims—that T1 drops with drive amplitude under a Purcell notch at the qubit frequency and that Lindblad shows qualitative defects—rest on the specific form chosen for the microscopic transmission-line bath spectral density and system-bath coupling. The manuscript does not demonstrate that this bath reproduces the known weak-drive Purcell relaxation rate or matches measured transmission data, leaving open the possibility that the reported high-drive T1 degradation and Lindblad discrepancies are artifacts of the bath choice rather than generic features. (Bath model and numerical results sections)
- [Comparison to Lindblad subsection] The comparison between the first-principles evolution and the Lindblad master equation is described as showing 'general qualitative defects,' yet the manuscript provides only qualitative statements without quantitative metrics (e.g., relative error in predicted T1 versus drive amplitude or fidelity) that would allow a reader to assess the practical importance of the discrepancies. (Comparison to Lindblad subsection)
minor comments (1)
- [Methods] Notation for the bath spectral density and the precise definition of the notch-filter transmission function should be stated explicitly in an equation rather than only in prose, to facilitate reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review. The comments highlight important points for strengthening the manuscript, and we address each major comment below. We will incorporate revisions to address the concerns.
read point-by-point responses
-
Referee: [Bath model and numerical results sections] The central claims—that T1 drops with drive amplitude under a Purcell notch at the qubit frequency and that Lindblad shows qualitative defects—rest on the specific form chosen for the microscopic transmission-line bath spectral density and system-bath coupling. The manuscript does not demonstrate that this bath reproduces the known weak-drive Purcell relaxation rate or matches measured transmission data, leaving open the possibility that the reported high-drive T1 degradation and Lindblad discrepancies are artifacts of the bath choice rather than generic features. (Bath model and numerical results sections)
Authors: We agree that explicit verification of the weak-drive limit is essential to establish that the observed high-drive effects are not artifacts. In the revised manuscript, we will add a new subsection (or extended discussion) in the Bath model section, including a figure or calculation, demonstrating that for low readout drive amplitudes the qubit T1 extracted from the first-principles unitary evolution matches the analytical Purcell relaxation rate computed from the bath spectral density evaluated at the qubit frequency. This will confirm consistency with known linear-response physics. Regarding matching measured transmission data, our microscopic transmission-line model follows the standard derivation used in circuit QED literature and reproduces the expected linear transmission spectrum (including notch-filter features) when the system is in the ground state; we will add a short clarifying paragraph with a reference to this property. These additions should address the concern that the high-drive T1 degradation and Lindblad discrepancies could be bath-specific artifacts. revision: yes
-
Referee: [Comparison to Lindblad subsection] The comparison between the first-principles evolution and the Lindblad master equation is described as showing 'general qualitative defects,' yet the manuscript provides only qualitative statements without quantitative metrics (e.g., relative error in predicted T1 versus drive amplitude or fidelity) that would allow a reader to assess the practical importance of the discrepancies. (Comparison to Lindblad subsection)
Authors: We concur that quantitative metrics would make the comparison more rigorous and allow readers to judge the practical relevance of the discrepancies. In the revised manuscript, we will augment the Comparison to Lindblad subsection with quantitative measures, including plots or tables of the relative error in T1 between the first-principles simulation and the Lindblad prediction as a function of drive amplitude. We will also include additional quantitative comparisons, such as differences in the emitted spectrum or estimated readout fidelity, to better quantify the impact of the qualitative defects. revision: yes
Circularity Check
No significant circularity in first-principles unitary simulation
full rationale
The paper performs a direct numerical simulation of the full unitary dynamics starting from the standard circuit-QED Hamiltonian coupled to an explicit microscopic model of the transmission-line bath whose spectrum can be chosen arbitrarily (including filters). The reported drive-amplitude dependence of T1 and the qualitative differences from the Lindblad equation are obtained by evolving this model under different bath spectra; they are not obtained by fitting parameters to a subset of data and then relabeling the fit as a prediction, nor by any self-definitional closure or load-bearing self-citation chain. The derivation chain therefore remains self-contained and independent of the target observables.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
inside of a multi-level system representing the qubit, we consider atwo-level system coupled to a driven harmonic readout resonator [3, 4], and we extract the drive power dependence of the relaxation rate of the qubit, for arbitrary spectral distributions of the bath. This captures behaviors that are beyond the reach of previous models in the presence of ...
- [2]
-
[4]
A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. A69, 062320 (2004)
work page 2004
- [5]
-
[6]
R. Acharya, D. A. Abanin, L. Aghababaie-Beni, oth- ers, Google Quantum AI, and Collaborators, Nature638, 920–926 (2025)
work page 2025
-
[7]
S. Krinner, N. Lacroix, A. Remm, A. Di Paolo, E. Genois, C. Leroux, C. Hellings, S. Lazar, F. Swiadek, J. Her- rmann, G. J. Norris, C. K. Andersen, M. M¨ uller, A. Blais, C. Eichler, and A. Wallraff, Nature605, 669 (2022)
work page 2022
- [8]
-
[9]
V. Champain, S. Zihlmann, A. Chessari, B. Bertrand, H. Niebojewski, ´E. Dumur, X. Jehl, V. Schmitt, B. Brun, C. Winkelmann, Y. M. Niquet, M. Filippone, S. De Franceschi, and R. Maurand, Phys. Rev. Appl.23, 034067 (2025)
work page 2025
-
[10]
G. Burkard, T. D. Ladd, A. Pan, J. M. Nichol, and J. R. Petta, Rev. Mod. Phys.95, 025003 (2023)
work page 2023
-
[11]
V. B. Braginsky, Y. I. Vorontsov, and K. S. Thorne, Sci- ence209, 547 (1980)
work page 1980
-
[12]
D. Sank, Z. Chen, M. Khezri, J. Kelly, R. Barends, B. Campbell, Y. Chen, B. Chiaro, A. Dunsworth, A. Fowler, and et al., Physical Review Letters117, 10.1103/physrevlett.117.190503 (2016)
- [13]
-
[14]
W. Dai, S. Hazra, D. K. Weiss, P. D. Kurilovich, T. Con- nolly, H. K. Babla, S. Singh, V. R. Joshi, A. Z. Ding, P. D. Parakh, J. Venkatraman, X. Xiao, L. Frunzio, and 6 M. H. Devoret, arXiv 10.48550/arXiv.2506.24070 (2025), 2506.24070
-
[15]
D. H. Slichter, R. Vijay, S. J. Weber, S. Boutin, M. Bois- sonneault, J. M. Gambetta, A. Blais, and I. Siddiqi, Phys. Rev. Lett.109, 153601 (2012)
work page 2012
-
[16]
S. O. Mundhada, Bulletin of the American Physical So- ciety61, Number 2(2016)
work page 2016
-
[17]
D. Gusenkova, M. Spiecker, R. Gebauer, M. Willsch, F. Valenti, N. Karcher, L. Gr¨ unhaupt, I. Takmakov, P. Winkel, D. Rieger, A. V. Ustinov, N. Roch, W. Werns- dorfer, K. Michielsen, O. Sander, and I. M. Pop, arXiv 10.1103/PhysRevApplied.15.064030 (2020), 2009.14785
- [18]
-
[19]
A. Bista, M. Thibodeau, K. Nie, K. Chow, B. K. Clark, and A. Kou, arXiv 10.48550/arXiv.2501.17807 (2025), 2501.17807
-
[20]
A. S. May, L. Sutevski, J. Solard, G. Cardoso, L. Carde, L. Pallegoix, R. Lescanne, D. Vion, P. Bertet, and E. Flurin, arXiv 10.48550/arXiv.2502.14804 (2025), 2502.14804
-
[21]
M. Boissonneault, J. M. Gambetta, and A. Blais, Phys. Rev. A79, 013819 (2009)
work page 2009
-
[22]
E. A. Sete, J. M. Gambetta, and A. N. Korotkov, Phys. Rev. B89, 104516 (2014)
work page 2014
-
[23]
A. Petrescu, M. Malekakhlagh, and H. E. T¨ ureci, Phys. Rev. B101, 134510 (2020)
work page 2020
- [24]
-
[25]
R. Lescanne, L. Verney, Q. Ficheux, M. H. Devoret, B. Huard, M. Mirrahimi, and Z. Leghtas, Phys. Rev. Ap- plied11, 014030 (2019)
work page 2019
- [26]
-
[27]
R. Shillito, A. Petrescu, J. Cohen, J. Beall, M. Hauru, M. Ganahl, A. G. Lewis, G. Vidal, and A. Blais, Phys. Rev. Appl.18, 034031 (2022)
work page 2022
- [28]
-
[29]
M. F. Dumas, B. Groleau-Par´ e, A. McDonald, M. H. Mu˜ noz-Arias, C. Lled´ o, B. D’Anjou, and A. Blais, Phys. Rev. X14, 041023 (2024)
work page 2024
-
[30]
T. Thorbeck, Z. Xiao, A. Kamal, and L. C. G. Govia, Phys. Rev. Lett.132, 090602 (2024)
work page 2024
-
[31]
H. P. Breuer and F. Petruccione,The theory of open quantum systems(Oxford University Press, Great Clarendon Street, 2002)
work page 2002
- [32]
- [33]
-
[34]
H. P. Breuer and M. Holthaus, Zeitschrift f¨ ur Physik D Atoms, Molecules and Clusters11, 1–14 (1989)
work page 1989
-
[35]
A. O. Caldeira and A. J. Leggett, Ann. Phys.149, 374 (1983)
work page 1983
-
[36]
J. Haegeman, J. I. Cirac, T. J. Osborne, I. Piˇ zorn, H. Ver- schelde, and F. Verstraete, Phys. Rev. Lett.107, 070601 (2011)
work page 2011
-
[37]
A. W. Chin, A. Rivas, S. F. Huelga, and M. B. Plenio, Journal of Mathematical Physics51, 092109 (2010)
work page 2010
-
[38]
B. Foxen, J. Y. Mutus, E. Lucero, R. Graff, A. Megrant, Y. Chen, C. Quintana, B. Burkett, J. Kelly, E. Jeffrey, Y. Yang, A. Yu, K. Arya, R. Barends, Z. Chen, B. Chiaro, A. Dunsworth, A. Fowler, C. Gidney, M. Giustina, T. Huang, P. Klimov, M. Neeley, C. Neill, P. Roushan, D. Sank, A. Vainsencher, J. Wenner, T. C. White, and J. M. Martinis, Quantum Science ...
work page 2017
-
[39]
M. Mirhosseini, E. Kim, V. S. Ferreira, M. Kalaee, A. Sipahigil, A. J. Keller, and O. Painter, Nature Com- munications9, 3706 (2018)
work page 2018
-
[40]
H. Putterman, K. Noh, R. N. Patel, G. A. Peairs, G. S. MacCabe, M. Lee, S. Aghaeimeibodi, C. T. Hann, I. Jar- rige, G. Marcaud, Y. He, H. Moradinejad, J. C. Owens, T. Scaffidi, P. Arrangoiz-Arriola, J. Iverson, H. Levine, F. G. Brand˜ ao, M. H. Matheny, and O. Painter, Physical Review X15, 10.1103/physrevx.15.011070 (2025)
-
[41]
H. Putterman, K. Noh, C. T. Hann, G. S. MacCabe, et al., Nature638, 927–934 (2025)
work page 2025
- [42]
-
[43]
Y.-X. Xiao, D. Feng, X.-Y. Gu, G.-H. Liang, M.-C. Wang, Z.-Y. Peng, B.-J. Chen, Y. Yan, Z.-Y. Mei, S.-L. Zhao, Y.-Z. Bu, C.-L. Deng, K. Yang, Y. Tian, X. Song, D. Zheng, Y.-X. Zhang, Y.-H. Shi, Z. Xiang, K. Xu, and H. Fan, Flexible readout and unconditional reset for su- perconducting multi-qubit processors with tunable pur- cell filters (2025), arXiv:250...
-
[44]
M. D. Reed, B. R. Johnson, A. A. Houck, L. Di- Carlo, J. M. Chow, D. I. Schuster, L. Frunzio, and R. J. Schoelkopf, Applied Physics Letters96, 10.1063/1.3435463 (2010)
-
[45]
H. Yan, X. Wu, A. Lingenfelter, Y. J. Joshi, G. Anders- son, C. R. Conner, M.-H. Chou, J. Grebel, J. M. Miller, R. G. Povey, H. Qiao, A. A. Clerk, and A. N. Cleland, Applied Physics Letters123, 10.1063/5.0161893 (2023)
-
[46]
M. Cattaneo and G. S. Paraoanu, Advanced Quantum Technologies4, 2100054 (2021)
work page 2021
-
[47]
Schollw¨ ock, Annals of Physics326, 96 (2011)
U. Schollw¨ ock, Annals of Physics326, 96 (2011)
work page 2011
- [48]
-
[49]
A. Sharafiev, M. L. Juan, O. Gargiulo, M. Zanner, S. W¨ ogerer, J. J. Garc´ ıa-Ripoll, and G. Kirchmair, Quan- tum5, 474 (2021)
work page 2021
-
[50]
Z. Huang, J. S.-C. Hung, M. Raha, M.-H. Chou, H. Levine, A. Retzker, C. T. Hann, D. Hover, F. G. S. L. Brand˜ ao, A. A. Clerk, A. Haim, and O. Painter, arXiv 10.48550/arXiv.2603.13857 (2026), 2603.13857
-
[51]
For consistency across the different simulations, we la- bel thex-axis with ¯naccording to the expression under Eq. (5), where the drive frequencyω d is set according to the calibrated value ¯ωabove, andω a is replaced byω a1, the resonator frequency as pulled by the qubit excited state
-
[52]
´Emile Cochin, J. Keeling, B. W. Lovett, and A. W. Chin, Efficient construction of time-invariant process ten- sors for simulating high-dimensional non-markovian open quantum systems (2026), arXiv:2603.06840 [quant-ph]
-
[53]
T. Lacroix, B. Le D´ e, A. Riva, A. J. Dunnett, and A. W. Chin, The Journal of Chemical Physics161, 084116 (2024)
work page 2024
-
[54]
A. Dunnett, T. Lacroix, B. Le D´ e, and A. Riva, sharelo- qs/MPSDynamics (2021)
work page 2021
- [55]
-
[56]
A. Riva, D. Tamascelli, A. J. Dunnett, and A. W. Chin, Phys. Rev. B108, 195138 (2023). 7
work page 2023
-
[57]
F. Beaudoin, J. M. Gambetta, and A. Blais, Physical Review A84, 10.1103/physreva.84.043832 (2011). 8 SUPPLEMENTAL MATERIAL This Supplemental Material is organized as follows. Sec. A outlines how the initial model (Eq. (1)) is transformed into the equations actually simulated using the chain mapping [36], including the methodology for reconstructing the sp...
-
[58]
(2), obtaining the eigenvalues {ωk}and corresponding eigenvectors{v (k)}
Diagonalize theN×NJacobi matrix of the chain coefficients, the second line of Eq. (2), obtaining the eigenvalues {ωk}and corresponding eigenvectors{v (k)}
-
[59]
Givenk 2 0 defined in Eq. (9), we compute the weights (withg k ≡ √wk) as wk =k 2 0 v(k) 1 2 ,(11) wherev (k) 1 is the first element of thek−th eigenvector
-
[60]
To obtain a smoother function, we convolve the discretized spectral density of Eq. (10) with a Gaussian or Lorentzian kernelf η, whereηgives the broadening, Jsmooth(ω) = (J∗f η)(ω) = X k wkfη(ω−ω k).(12) The broadening parameterηhas to be chosen to be slightly larger than the interval between the sampled frequenciesω k to ensure a smooth profile without w...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.