arxiv: 2405.17192 · v3 · submitted 2024-05-27 · ⚛️ physics.ins-det · astro-ph.IM· hep-ex· quant-ph

Recognition: unknown

Quantum Parity Detectors: a qubit based particle detection scheme with meV thresholds for rare-event searches

Karthik Ramanathan , Brandon J. Sandoval , John E. Parker , Lalit M. Joshi , Andrew D. Beyer , Pierre M. Echternach , Serge Rosenblum , Sunil R. Golwala

Authors on Pith no claims yet

Pith reviewed 2026-05-06 19:00 UTC · model claude-opus-4-7

classification ⚛️ physics.ins-det astro-ph.IMhep-exquant-ph PACS 85.25.-j95.35.+d74.40.-n

keywords quantum parity detectorsuperconducting qubit sensorquasiparticle tunnelingrare-event searchdark matter detectioncoherent neutrino scatteringcryogenic phonon detectormultiplexed resonator readout

0 comments

The pith

A particle detector that reads each event as a binary flip in a superconducting qubit's charge parity, with a designed path to sub-eV thresholds.

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

The paper introduces quantum parity detectors, a particle-detection concept built on the same physics that makes charge parity a nuisance in superconducting qubits: a single quasiparticle tunneling across a weak link flips the device's parity, and that flip is observable. Patterned on a crystalline substrate, such a qubit becomes a phonon sensor: energy deposited in the crystal generates phonons that break Cooper pairs in the qubit's superconducting film, producing quasiparticles whose tunneling events register as parity flips. The authors provide an analytic formulation that turns the time series of parity flips into a reconstructed energy, and describe resonator-coupled readout schemes compatible with frequency-multiplexed arrays. They estimate sensitivities reaching the meV scale and outline an R&D path toward sub-eV thresholds, the regime relevant for low-mass dark matter searches and coherent neutrino scattering.

Core claim

The paper proposes quantum parity detectors (QPDs): superconducting qubit elements patterned on a crystalline absorber that record particle interactions as binary flips of the device's charge parity. A phonon shower from an energy deposit breaks Cooper pairs in the qubit film and produces a single quasiparticle tunneling event across a coherent weak link; counting such parity flips, with the right readout, encodes the deposited energy. The authors lay out resonator-coupled readout architectures, give an analytic prescription for reconstructing energies from parity-flip statistics, and describe how to multiplex many sensors. They argue this scheme can push particle-detection thresholds 100–10

What carries the argument

A surface-patterned superconducting qubit whose charge parity is flipped, one quasiparticle at a time, by phonons from particle interactions in the underlying crystal; combined with an analytic energy-reconstruction map from parity-flip statistics to deposited energy, and resonator-multiplexed readout for large arrays.

If this is right

Rare-event searches for sub-GeV dark matter and coherent neutrino-nucleus scattering gain a sensor class with thresholds 100–1000× below current calorimeters.
Existing superconducting-qubit fabrication and dispersive readout infrastructure can be repurposed for particle detection rather than developed from scratch.
Frequency-multiplexed resonator readout makes kilopixel-scale QPD arrays plausible, enabling target masses competitive with cryogenic bolometer arrays.
An analytic energy-reconstruction map means each detected event yields a calibrated energy, not just a hit/no-hit tag, supporting spectroscopy at meV scales.
The same devices could probe quasiparticle poisoning physics in qubits as a side benefit, since the detection signal IS the poisoning event.

Where Pith is reading between the lines

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

The achievable threshold will likely be set not by readout sensitivity but by the intrinsic non-equilibrium quasiparticle background — stray infrared, cosmic-ray-induced phonon bursts, and substrate stress-driven pair breaking — which already produce parity flips at tens of Hz to kHz in current transmons.
Phonon collection efficiency from a macroscopic crystal into a small qubit-film footprint is the natural place where the meV-threshold promise is won or lost; phonon-focusing or interface engineering may be needed to make the geometry work at scale.
Because each quasiparticle tunneling event is binary, the energy resolution is fundamentally counting-statistics limited: doubling resolution requires four times the parity-flip yield per unit deposited energy, which constrains film gap engineering.
If the parity-flip rate can be made cleanly Poissonian above background, QPDs could double as a calibration tool for the long-standing quasiparticle-poisoning problem in quantum computing — turning a noise source into a measurand.

Load-bearing premise

That the chain of phonon production, transport into the qubit film, Cooper-pair breaking, and quasiparticle tunneling is efficient and clean enough at meV energies, and that the rate of "fake" parity flips from stray quasiparticles is low enough to leave a real signal visible.

What would settle it

Build a prototype QPD on a crystalline substrate, expose it to a calibrated low-energy source (e.g., optical photons or a low-energy X-ray fluorescence line attenuated to meV-scale phonon depositions), and measure whether the parity-flip rate tracks deposited energy as the analytic model predicts, with intrinsic flip backgrounds low enough to resolve sub-eV depositions above the poisoning floor seen in transmons.

Figures

Figures reproduced from arXiv: 2405.17192 by Andrew D. Beyer, Brandon J. Sandoval, John E. Parker, Karthik Ramanathan, Lalit M. Joshi, Pierre M. Echternach, Serge Rosenblum, Sunil R. Golwala.

**Figure 1.** Figure 1: FIG. 1. Schematic of the operating principle of both QPD styles. view at source ↗

**Figure 2.** Figure 2: FIG. 2. Stylized mask layouts of QPDs. Items in red are made of a high gap material (e.g. Nb or Ta), used for the feedline and resonator. Blue view at source ↗

**Figure 3.** Figure 3: FIG. 3 view at source ↗

**Figure 4.** Figure 4: FIG. 4. Stylized chart showing the string of efficiencies that mod view at source ↗

**Figure 5.** Figure 5: FIG. 5. Example energy level structure from solving the CPB Hamil view at source ↗

**Figure 6.** Figure 6: FIG. 6 view at source ↗

**Figure 7.** Figure 7: FIG. 7 view at source ↗

**Figure 8.** Figure 8: FIG. 8. Simulated charge parity state time-series (blue telegraph view at source ↗

**Figure 9.** Figure 9: FIG. 9 view at source ↗

**Figure 10.** Figure 10: FIG. 10 view at source ↗

**Figure 11.** Figure 11: FIG. 11 view at source ↗

**Figure 12.** Figure 12: FIG. 12. Fraction of trapped quasiparticles as computed from Eq. view at source ↗

**Figure 13.** Figure 13: FIG. 13 view at source ↗

**Figure 14.** Figure 14: FIG. 14. Trap quasiparticle density as a function of transmission view at source ↗

read the original abstract

The next generation of rare-event searches, such as those aimed at determining the nature of particle dark matter or in measuring fundamental neutrino properties, will benefit from particle detectors with thresholds at the meV scale, 100-1000x lower than currently available. Quantum parity detectors (QPDs) are a class of proposed quantum devices, extending recent work on superconducting qubit sensors, that exploit the fingerprints of single quasiparticle tunneling across a coherent weak-link as their detection concept. As envisioned, phonons generated by particle interactions within a crystalline substrate cause an eventual quasiparticle cascade within a surface-patterned superconducting qubit element. This process alters the fundamental charge parity of the device in a binary manner, which can be used to deduce the initial properties of the energy deposition. This work lays out multiple resonator coupled readout schemes depending on qubit architecture, provides an analytic formulation for reconstructing sensor energies, and details strategies for multiplexing large arrays of sensors. We further compute the sensitivity of QPDs and detail an R&D pathway to demonstrating sub-eV energy deposit thresholds.

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.

Desk Editor's Note private letter to a colleague

Solid concept paper for a meV-threshold qubit-based particle detector; the projected sensitivity hinges on quasiparticle poisoning rates the abstract doesn't quantify, but that's a paper to read, not dismiss.

read the letter

Quick take: this is the Caltech/JPL group writing down, as a coherent detector class, what the qubit community has been circling for a few years — using charge-parity flips from single quasiparticle tunneling as the readout for athermal phonons in a substrate. The framing as "QPDs," with explicit readout architectures, an analytic energy-reconstruction formalism, and a multiplexing scheme, is the actual contribution. The underlying physics (parity jumps in transmons, phonon-mediated QP creation) is borrowed from Serniak, Wilen, Cardani et al., and they're not pretending otherwise.

What's genuinely new, as far as I can tell from the abstract: the integrated detector packaging, the analytic mapping from parity-flip statistics back to deposited energy, and a concrete R&D roadmap toward sub-eV thresholds. That's a useful piece of work even before any hardware. The dark-matter and CEvNS communities have been hungry for a credible meV-threshold path, and this slots in alongside TES, KID, and SNSPD-based proposals as a fourth lane worth tracking.

Soft spot, and it's the obvious one: a binary parity observable means energy reconstruction is essentially counting flips across an array, and the background that matters is the intrinsic non-equilibrium QP poisoning rate, which in current best transmons sits at 10²–10⁴ Hz per qubit. The stress-test note has this right — if you multiplex N qubits, that floor scales linearly unless coincidence/timing cuts are very tight. Whether the projected meV sensitivity assumes a poisoning floor below what's been demonstrated is the load-bearing question, and it can't be answered from the abstract. I'd want to see, in the body, (a) the assumed intrinsic flip-rate floor, (b) how the analytic reconstruction handles Poisson background flips, and (c) whether the sub-eV pathway requires hardware improvements that haven't been shown yet. My guess from knowing this group is that they address at least (a) and (b) explicitly; the reader's "abstract-only" caution is appropriate but shouldn't be read as skepticism about the paper itself.

Audience: rare-event-search detector people, qubit-sensor people, and anyone planning low-mass DM or coherent neutrino-scattering experiments in the next decade. It's a concept/architecture paper, not a measurement paper, and should be read as such.

Recommendation: send it out. This deserves a referee who knows both transmon QP physics and cryogenic calorimetry — there aren't many, but they exist. I'd cite it.

Referee Report

4 major / 4 minor

Summary. The manuscript proposes "Quantum Parity Detectors" (QPDs): rare-event particle detectors in which athermal phonons from an energy deposition in a crystalline substrate are absorbed in a surface-patterned superconducting qubit element, producing quasiparticles whose tunneling across the qubit's weak link flips the charge parity of the device. Counting parity flips serves as the energy estimator. The authors describe several resonator-coupled readout architectures suited to different qubit topologies, give an analytic formulation for reconstructing deposited energy from parity-flip statistics, outline multiplexing strategies for large arrays, compute projected sensitivities, and lay out an R&D pathway toward sub-eV thresholds — claiming a 100–1000× improvement over current technology in the meV regime relevant to light dark matter and coherent-neutrino searches.

Significance. If the projected meV-scale thresholds and multiplexing scheme hold up, this would be a genuinely impactful contribution to the rare-event detection community: it identifies a path that leverages the rapid maturation of superconducting-qubit fabrication and multiplexed dispersive readout for a problem (sub-eV calorimetry over kg-scale targets) that currently has no clean solution. The paper's strengths, as advertised, are (i) an analytic, rather than purely Monte-Carlo, energy-reconstruction formulation that should make sensitivity projections checkable; (ii) explicit consideration of multiplexing, which is the practical bottleneck for any large detector array; and (iii) a falsifiable R&D roadmap with numerical sensitivity targets. These are concrete, testable claims rather than aspirational ones, which is the right form for a proposal paper of this kind.

major comments (4)

[Sensitivity / R&D pathway] The load-bearing number for a meV-threshold claim is the intrinsic charge-parity-flip background rate from sources unrelated to the substrate phonon cascade — non-equilibrium quasiparticles from stray pair-breaking IR, cosmic-ray-induced phonon bursts in the chip itself, substrate stress relaxation, and ambient QP density. Published transmon parity-flip rates in shielded dilution-fridge setups are typically 10^2–10^4 Hz (Serniak et al.; Wilen et al.; correlated-burst work by Cardani, Vepsäläinen, McEwen et al.). The abstract asserts meV sensitivity but does not state the assumed intrinsic flip-rate floor. The manuscript should make explicit, in the sensitivity section, what background flip rate is assumed for each curve, and how that compares with state-of-the-art measurements; otherwise the projected thresholds cannot be evaluated.
[Detection concept (phonon → QP → tunneling chain)] The chain phonon generation in substrate → phonon transport to qubit film → Cooper-pair breaking → quasiparticle tunneling across the weak link → observable parity flip has several efficiency factors (phonon collection efficiency into the patterned film, downconversion above 2Δ, quasiparticle trapping vs. tunneling, parity-flip detection efficiency in the chosen readout cycle). The energy-reconstruction formulation depends on these. The manuscript should specify which of these factors are taken from measurement, which are assumed, and propagate their uncertainties into the sensitivity curves; a single-number 'efficiency' that absorbs the whole chain is not adequate for a 100–1000× threshold claim.
[Energy reconstruction / readout duty cycle] An analytic energy reconstruction from parity-flip counts has to confront the finite parity-measurement bandwidth: the observable is integrated flips per readout window, which saturates when the per-event flip rate becomes comparable to the readout repetition rate, and is Poisson-limited at the low-energy end. The manuscript should state the readout cycle time assumed for the sensitivity projections, the dynamic range that follows, and where on the response curve the quoted meV threshold sits. A related point: if the multiplexed readout reduces effective per-qubit sampling rate, the trade-off between channel count and threshold should be made explicit.
[Multiplexing / array scaling] The claim that the scheme supports large arrays via resonator-coupled multiplexed readout is central to the kg-scale-target use case but is in tension with crosstalk and correlated-poisoning effects observed in current multi-qubit chips (correlated parity bursts from cosmic-ray impacts on the chip itself, e.g., McEwen et al. 2022). The manuscript should address whether such correlated events constitute a veto-able background or a non-recoverable dead-time/livetime cost in the proposed architecture, since this directly affects the achievable exposure at the quoted thresholds.

minor comments (4)

[Abstract / terminology] The phrase 'fingerprints of single quasiparticle tunneling across a coherent weak-link' is suggestive but not precise; specify in the abstract or introduction what measurable signature (parity-dependent transition frequency, dispersive shift, etc.) is meant for each architecture considered.
[Abstract] '100–1000× lower than currently available' would benefit from a citation anchor in the introduction (e.g., current best TES/SNSPD/KID thresholds for the same target mass class), so the comparison is unambiguous.
[Sensitivity figures] Sensitivity curves should be presented with the assumed background rate, integration time, and target mass clearly labeled on the figure (not only in the caption or text), to allow direct comparison with other proposed technologies.
[References] Ensure the relevant prior superconducting-qubit-as-sensor literature is cited (Serniak, Vepsäläinen, Wilen, McEwen, Cardani, and the Aalto/PSI/Fermilab QP-poisoning programs), including those that quantify the parity-flip background floor that this proposal must beat.

Simulated Author's Rebuttal

4 responses · 2 unresolved

We thank the referee for a careful and constructive report that engages directly with the load-bearing technical claims of the manuscript rather than its framing. The four major comments — (1) the assumed intrinsic parity-flip background floor, (2) the efficiency budget along the phonon→QP→tunneling→readout chain, (3) the readout duty-cycle / dynamic-range implications of an analytic reconstruction, and (4) correlated parity-burst backgrounds in multiplexed arrays — are all legitimate, and we agree that the sensitivity curves cannot be properly evaluated without each being made explicit. We will revise accordingly. Our responses below indicate where existing material in the manuscript already addresses a point (and how we will make it more visible), and where new text, figures, or tabulated assumptions are required. We note that the paper is intended as a proposal/roadmap, so several of these numbers are necessarily projected from current best-measured values; we will be more transparent about which is which. None of the referee's points require us to walk back a quantitative claim, but together they require a substantially more explicit "assumptions and uncertainties" treatment than the current draft contains, and we will produce that in the revision.

read point-by-point responses

Referee: Sensitivity / R&D pathway: the meV-threshold claim depends on the intrinsic (non-signal) charge-parity flip-rate floor, which is not stated. State-of-the-art shielded transmons show ~10^2–10^4 Hz. The sensitivity section must declare the assumed background flip rate per curve and compare to measurements.

Authors: Agreed. The sensitivity projections in the current draft do assume an intrinsic background flip rate, but it is buried in the reconstruction formalism rather than displayed alongside each curve. In the revision we will (i) add an explicit table listing, for each sensitivity curve, the assumed intrinsic parity-flip rate, its origin (residual nonequilibrium QP density, IR pair-breaking, ambient phonon hits on the qubit film), and the corresponding literature benchmark (Serniak et al., Wilen et al., Vepsäläinen et al., McEwen et al.); (ii) plot the sensitivity as a band rather than a line, bracketed by today's best demonstrated rates (~few × 10^2 Hz in well-shielded devices) and the target rates needed to reach the quoted meV thresholds; and (iii) state clearly that the most aggressive (sub-eV) curves require an additional ~order-of-magnitude reduction in the intrinsic flip rate beyond present demonstrations, which is the explicit subject of the R&D roadmap. This makes the threshold claim falsifiable in exactly the way the referee asks. revision: yes
Referee: Detection concept: the chain phonon generation → transport → pair-breaking → QP tunneling → parity-flip detection has multiple efficiency factors. A single lumped efficiency is inadequate for a 100–1000× threshold claim; uncertainties should be propagated.

Authors: We accept this criticism. The current manuscript folds several efficiencies into a combined factor and only partially separates them. In the revised version we will break the chain into its individual links — (a) phonon collection into the patterned superconducting film (set by film coverage and acoustic mismatch, taken from TES/KID literature), (b) downconversion and pair-breaking yield above 2Δ (from QP-poisoning and kinetic-inductance studies), (c) QP tunneling vs. trapping/recombination at the junction (from recent transmon QP-dynamics measurements), and (d) parity-flip detection efficiency per readout cycle (set by readout fidelity and duty cycle; see our response to point 3). For each we will indicate whether the value is measured, extrapolated, or assumed, attach a stated uncertainty, and propagate these into the sensitivity bands. We will also add a sentence flagging the parameters that most strongly drive the threshold so future experimental work can target them. revision: yes
Referee: Energy reconstruction / readout duty cycle: the analytic estimator saturates at high flip rate (relative to readout repetition) and is Poisson-limited at low energy. The manuscript should give the assumed readout cycle time, the resulting dynamic range, and where on the response curve the meV threshold sits. Multiplexing reductions in per-qubit sampling rate must be folded in.

Authors: This is a fair and important point. The analytic reconstruction in Sec. (energy-reconstruction) is derived under the implicit assumption of a parity-measurement repetition rate fast compared to the per-event tunneling rate; this assumption is stated but its consequences for dynamic range are not displayed. In revision we will (i) state the assumed readout cycle time explicitly (we have been using ~1–10 μs per parity measurement, consistent with current dispersive-readout demonstrations), (ii) plot the response curve (mean detected flips vs. deposited energy) and mark where the saturation knee and the Poisson-limited toe lie, (iii) identify the meV threshold's location on that curve — it sits in the Poisson-limited regime, where the threshold is set by single-flip discrimination against the intrinsic floor of point 1 — and (iv) add a subsection quantifying the channel-count vs. per-channel-bandwidth trade-off in the multiplexed scheme, showing how the threshold degrades as the per-qubit sampling rate is reduced. This addresses the implicit coupling between points 3 and 4. revision: yes
Referee: Multiplexing / array scaling: correlated parity bursts from cosmic-ray impacts on the readout chip itself (e.g., McEwen et al. 2022) threaten the kg-scale claim. The manuscript should specify whether such events are veto-able background or non-recoverable livetime loss.

Authors: We agree this needs explicit treatment; the present draft mentions multiplexing crosstalk but does not address correlated bursts originating in the chip substrate as opposed to the target crystal. In the revision we will add a dedicated paragraph distinguishing two classes of correlated events: (a) bursts originating in the target crystal, which are the signal class and are reconstructed by the spatial flip pattern across the array — these are not a background; and (b) bursts originating in the readout chip substrate or housing, which produce array-wide correlated parity flips on a timescale set by chip-phonon ballistic transport. Class (b) is vetoable in principle by its characteristic multiplicity and timing signature (essentially a coincidence cut across many channels with no localization in the target), and we will quote an estimated livetime cost based on the cosmic-ray and ambient-radiogenic rates relevant to a deep-site deployment, drawing on the McEwen et al. burst rate and on the Cardani/Vepsäläinen correlated-event measurements. We will be candid that this is an estimate, not a measurement on a QPD prototype, and we will list it as an R&D milestone. revision: yes

standing simulated objections not resolved

We cannot, in this paper, demonstrate experimentally that the intrinsic parity-flip floor required for the most aggressive sub-eV curves is achievable; that demonstration is the explicit purpose of the R&D roadmap and will require dedicated prototype measurements. The revision will make this limitation, and the gap to current best-measured rates, fully explicit rather than try to argue it away.
The correlated-burst livetime estimate in our response to point 4 is based on extrapolation from existing multi-qubit chip measurements rather than on a QPD-architecture measurement; we will flag it as such and not present it as a validated number.

Circularity Check

0 steps flagged

No significant circularity visible from the abstract; the QPD sensitivity claim is an engineering projection benchmarked against external particle-physics thresholds, not a self-referential derivation.

full rationale

Only the abstract is available, so a full chain audit of equations is not possible. On the abstract's own terms, the derivation chain is: (i) phonons from particle interactions, (ii) quasiparticle cascade in a surface-patterned qubit, (iii) charge-parity flips read out via resonator coupling, (iv) an analytic energy-reconstruction formula, (v) projected sensitivity and an R&D pathway to sub-eV thresholds. None of these steps, as stated, redefine the output in terms of the input or rename a fitted parameter as a prediction. The "100–1000x lower threshold" comparison is to currently available particle detectors — an external benchmark, not an internal fit. The "extending recent work on superconducting qubit sensors" phrasing is ordinary literature grounding, not a load-bearing self-citation that imports a uniqueness theorem or smuggles in an ansatz. The reader's skeptic concern (that the projected sensitivity may assume a parity-flip background floor below what current transmons achieve) is a correctness/assumption-validity risk, not circularity: it would only become circular if the paper fitted the background rate from the same data used to claim the sensitivity, which the abstract gives no evidence of. Without the body text, I cannot inspect Eq.-level structure of the analytic energy reconstruction to check whether any "predicted" quantity reduces to a fitted input by construction; if such a reduction exists it would live in the sensitivity computation, not the abstract. Honest finding: no circularity demonstrable from the available text. Score 1 reflects only the residual unverifiable possibility, not any identified circular step.

Axiom & Free-Parameter Ledger

3 free parameters · 3 axioms · 0 invented entities

The proposal layers a particle-detector concept on top of established qubit-physics primitives (parity-sensitive transmon readout, phonon-mediated quasiparticle generation). The load-bearing inputs are efficiency parameters in a chain rather than new postulated entities — there are no new particles, forces, or fundamental constants introduced. Free parameters live in the noise/efficiency budget for the sensitivity projection. The main domain assumptions are about the cleanliness of the phonon→QP→parity-flip chain and the achievable suppression of background quasiparticle poisoning. Abstract-only assessment limits how exhaustively these can be enumerated.

free parameters (3)

Phonon-to-quasiparticle conversion efficiency
Not given in abstract; sets the achievable threshold in the sensitivity projection.
Quasiparticle poisoning / background parity-flip rate
Determines the signal-to-noise floor; abstract does not state the assumed value.
Resonator readout SNR / integration time per parity measurement
Affects multiplexing limits and reconstruction fidelity.

axioms (3)

domain assumption Phonons from substrate interactions reach the qubit film and produce quasiparticles with sufficient yield at meV deposits.
Standard cryo-detector physics; quantitative validity at meV is the open question.
domain assumption Charge-parity flips from single QP tunneling are cleanly readable through the coupled resonator scheme.
Demonstrated in QC literature; reused here as the detection primitive.
ad hoc to paper Background QP-poisoning rates can be suppressed below the signal-flip rate at meV thresholds.
Not quantified in the abstract; the sub-eV roadmap depends on this.

pith-pipeline@v0.9.0 · 9869 in / 5721 out tokens · 86713 ms · 2026-05-06T19:00:12.543612+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Characterizing charge-parity detection based on an offset-charge-tunable transmon qubit via randomized benchmarking
quant-ph 2026-04 unverdicted novelty 4.0

Offset-charge-tunable transmon qubit achieves 99.37% fidelity in charge-parity mapping and over 93.4% in continuous monitoring at 4 μs intervals via randomized benchmarking.