REVIEW 4 major objections 4 minor 22 references
Computed models of natural radiation backgrounds in qubits and superconducting detectors
T0 review · 4 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Natural radiation backgrounds in superconducting qubit substrates are summarized by a single scaling law that gives three event rates to roughly ±25 percent accuracy.
desk verdict Useful compact background-rate formulas, but the elevation scale height λ is never given, so Eq. 1 cannot be evaluated at nonzero altitude as published. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the three-rate summary and the closed-form scaling law of Eq. (1). The rates are: $R$, the rate of any event depositing energy $E>0$; $M$, the rate of events depositing $E>1$ MeV, chosen because it marks the transition from terrestrial-gamma and charged-particle events to proton and neutron events and corresponds to roughly one event per hour in the nominal substrate; and $P$, the total power deposited. Equation (1) is separable: an area factor $A/(100\,\mathrm{mm}^2)$, a sum over six background sources of terms $c_s+g_s\tau$, $p_s\tau^{\beta_s}$, and $m_s\tau^{\alpha_s}$, and corrections $\kappa_c\kappa_{sh}\kappa_\rho$ for ceiling, shape, and density. This factorization is what carries the argument, because it reduces a six-dimensional simulation space to a table of constants that any user can evaluate directly.
What would settle it
Measure the three rates $R$, $M$, and $P$ in a 10 mm by 10 mm, 500 µm silicon substrate at sea level under roughly 20 cm of concrete, using a detector that records the full deposited-energy spectrum; if the measured values differ from Eq. (1) with Table II by more than the claimed ±25 percent, or if the thickness dependence of $M$ departs from the quoted power laws, the central claim would be undercut.
Extended reading notes
Core claim
On the paper's own terms, the central claim is that all simulated background behavior in superconducting-circuit substrates can be summarized by Eq. (1), which expresses the three summary rates as products of the wafer area, a sum over six source terms (the $^{40}$K, $^{232}$Th, and $^{238}$U decay chains, with the uranium and thorium chains each split into two halves around radon, plus cosmic rays), and three multiplicative corrections for ceiling, shape, and density. Each source term is proportional to its relative activity $\tilde a_s$; the gamma contributions scale with thickness $\tau \equiv t/500\,\mu\text{m}$ as $c_s + g_s\tau$ for $R$ and as power laws $\tau^{\beta_s}$ and $\tau^{\alpha_s}$ for $P$ and $M$, while the cosmic-ray contribution grows with elevation as $\exp(H/\lambda)$ with different scale heights for different particle species. The paper states that the absolute rates are accurate to about $\pm25\%$ over a range of realistic conditions, and that the model reproduces the measured event-rate spectrum of the authors' thermal kinetic-inductance detector.
Load-bearing premise
The fragile premise is that the simulated backgrounds, built from standard cosmic-ray spectra and typical European concrete, match real labs closely enough that the scaling formulas hold to the claimed 25 percent for any substrate and elevation, even though only one measurement is used for comparison.
Editorial extensions
If this is right
- Any lab can compute approximate background rates for its own substrate geometry, material, shielding, and elevation by evaluating Eq. (1) with Table II, without specialized simulation expertise.
- Thinning the substrate is an effective lever: the high-energy event rate scales as $M_\gamma \propto \tau^5$ for gamma rays and as $\tau^{1.8}$ for cosmic rays, so a thinner wafer suppresses the rare damaging events far more than it suppresses the total rate.
- Adding concrete ceiling reduces cosmic-ray $R$ and $P$ by only about 2% per 10 cm, so overhead shielding is not a strong mitigation against cosmic-ray backgrounds.
- The $M$ rate in the nominal 500 µm silicon substrate is about once per hour, the level flagged as a concern for quantum error correction.
- Denser and gallium-containing substrates have higher rates; for gallium the gamma-ray rates scale as if the density were increased by a source-dependent amount, and high-energy events scale roughly as $\tilde\rho^{2.7}$.
Reading between the lines
- Extension: if Eq. (1) generalizes as claimed, a lab can identify its dominant background term without simulation---at sea level terrestrial gamma rays set the low-energy rate, while the once-per-hour MeV events are mostly cosmic-ray protons and neutrons---so the relative benefit of underground siting versus additional local shielding follows directly from the two scale heights.
- Extension: the gamma-ray $M_\gamma \propto \tau^5$ exponent is estimated from very few simulated events; a testable consequence is that measuring MeV events in substrates of thickness 30, 100, and 500 µm should reveal this steep dependence, and if it does not, the exponent needs revision.
- Extension: because the nominal concrete activities are European averages, labs using unusual aggregates or with elevated radon should expect systematic offsets; assaying local concrete and rescaling via $\tilde a_s$ is a natural extension of the formula.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses Geant4/TOPAS with PARMA cosmic-ray spectra to simulate energy deposition from natural radiation (terrestrial gamma rays and cosmic rays) in silicon and other substrates typical of superconducting qubits and detectors. It distills the simulation output into three rates—R (any energy event), M (events >1 MeV), and P (total deposited power)—and packages them in Eq. (1) with parameters in Table II, intended to let experimenters estimate background rates as functions of substrate area, thickness, material, elevation, and ceiling thickness. The authors report consistency with their earlier TKID measurement at 1640 m altitude and claim roughly ±25% accuracy over realistic conditions.
Significance. If Eq. (1) and Table II are reliable, the paper provides a genuinely useful practical tool: experimenters could obtain approximate background rates for superconducting qubit and sensor substrates without running dedicated Monte Carlo simulations. The two-step re-aiming Monte Carlo scheme and the explicit three-rate parametrization are sensible and clearly described. However, the central formula as printed is incomplete because the cosmic-ray elevation scale height λ is never given, and the claimed accuracy is not backed by quantitative uncertainties or a quantitative validation. The paper's practical value therefore depends on addressing these gaps.
major comments (4)
- [Section VI, Eq. (1), Table II] The cosmic-ray elevation factor is defined as ãCR ≡ exp(H/λ) with "an appropriate scale height λ," but no value of λ is given anywhere in the paper: Table II has no λ entry, and Section V reports only species-dependent scale heights (about 5 km for muons, 1 km for nuclear particles, and intermediate values for e±, γ, and the overall rate). This is a load-bearing omission because Eq. (1) cannot be evaluated at any nonzero elevation without λ, including the 1640 m altitude of the paper's own validation measurement. Moreover, the Section V observation that scale heights differ by particle species directly undermines the single-exponential form for all three rates; at 1640 m, exp(1.64/1) ≈ 5.1 versus exp(1.64/3) ≈ 1.7, a factor-of-three spread that dwarfs the claimed ±25% accuracy. The authors must either supply the λ value used for each rate or reformulate the elevation dependence to reflect the species-dependent composition changes.
- [Table II and Section VII] Table II lists fitted parameters cs, gs, ps, βs, ms, αs, and ρGa,s without any statistical or systematic uncertainties, yet Section VII claims the absolute rates are accurate to about ±25%. The uncertainty claim is therefore unsupported: the reader cannot tell whether 25% reflects Monte Carlo statistical errors, PARMA model uncertainty, concrete activity variation, or the re-aiming approximation. I ask the authors to provide uncertainties for the fitted parameters (at least for the dominant terms) or to weaken the accuracy claim to a qualitative estimate until such uncertainties are quantified.
- [Section IV, Mγ scaling] The >1 MeV terrestrial gamma rate Mγ is stated to scale approximately as t^5 (Section IV), but the paper acknowledges that this is "difficult to estimate from the few simulated events above that energy." No event counts, confidence intervals, or goodness-of-fit measures are given, and Section VI similarly states that too few MeV-scale events were generated to characterize any shape correction to M. A power-law exponent of 5 inferred from a handful of events is not a reliable interpolation or extrapolation basis, and it directly affects Eq. (1)'s M row. The authors should provide the number of simulated events above 1 MeV for the relevant thicknesses and a quantitative measure of the fit quality, or explicitly mark the Mγ scaling as preliminary.
- [Section VI, validation against Ref. [8]] The only validation is a qualitative consistency statement against the authors' own TKID measurement at 1640 m altitude: "A TKID-based spectroscopic measurement [8] is consistent with the models' results under nominal conditions." No quantitative comparison (e.g., fitted normalization, chi-square, ratio of modeled to measured rates) is shown, and the validation measurement is at an altitude where the model cannot currently be evaluated because λ is missing. Without a quantitative validation, the ±25% accuracy claim in Section VII cannot be checked. I recommend adding a direct comparison plot or table with the measured and modeled R, P, and energy spectra, including the elevation scaling with a specified λ.
minor comments (4)
- [Table II] The table's header row and last row (Cosmic rays) are hard to parse because the columns ns, cs, gs, ps, βs, ms, αs, ρGa,s are not all applicable to cosmic rays; consider using separate sub-tables or explicit placeholders (e.g., em dashes) with a footnote explaining which coefficients are used for the cosmic-ray term in Eq. (1).
- [Section VI] The sentence "Each term in the event rate R depends linearly on τ, with both constants depending on the source" is confusing because Eq. (1) shows R depending on τ through gs τ only, with cs independent of τ; consider rewording to "the gamma-ray contribution to R scales linearly with τ, while the charged-particle contribution cs is thickness-independent."
- [Section V, Figure 3] The y-axis of Figure 3 is labeled "Event rate (s−1)" but the axis values appear to span 10−3 to 10−1; please clarify the units (per substrate? per cm2?) and ensure the axis label matches the text's description of rates in the nominal substrate.
- [Section IV, Figure 2] Panel (b) is labeled "P/t (keV s−1 µm−1)", but the text discusses P growing as t^1.12; please define whether P is the total deposited power or power per unit thickness, and keep the notation consistent with Eq. (1).
Circularity Check
No significant circularity: Eq. 1 is a fit to Geant4/PARMA simulations, not to the validation measurement; the only self-citation is a consistency check.
full rationale
The paper's central equations (Eq. 1 and Table II) are obtained by fitting analytic scaling forms to Geant4 simulation outputs driven by PARMA cosmic-ray spectra and nominal European concrete activities. None of the constants in Table II is fitted to the TKID measurement of Ref. [8]; that measurement appears only as a consistency check ('consistent with our earlier measurement'), so it is not an input to the model. The acknowledgment that co-authors of Ref. [8] helped 'set the direction' does not amount to tuning parameters. The elevation factor aCR = exp(H/λ) is an ansatz, and λ is not tabulated, which is an incompleteness/robustness problem rather than a circular one. Likewise, the sparse Mγ proportional to t^5 extrapolation is a statistical weakness. No load-bearing step in the derivation is equivalent by construction to its inputs, and no uniqueness or ansatz is imported from a self-citation. Score 0.
Assumptions & free parameters
free parameters (9)
- cs, charged-particle event-rate coefficients for 40K, 232Th-a, 232Th-b, 238U-a, 238U-b, cosmic rays =
2.2, 0.6, 1.5, 0.02, 1.9, 40 (10^-3 s^-1)
- gs, gamma-ray event-rate coefficients for the six sources =
6.8, 4.9, 6.6, 0.6, 11.7, 1.4 (10^-3 s^-1)
- ps, deposited-power coefficients for the six sources =
1.4, 0.5, 0.9, 0.03, 1.4, 8.0 (keV/s)
- beta_s, thickness power-law exponents for P =
1.12 for gamma sources; 1.0 for cosmic rays
- ms, >1 MeV event-rate coefficients for the six sources =
15, 2, 20, 0, 13, 180 (10^-6 s^-1)
- alpha_s, thickness power-law exponents for M =
5.0 for 40K, 4 for 232Th chains and 238U-b, 15 for 238U-a, 1.8 for cosmic rays
- rho_Ga,s, effective density adjustments for gallium substrates =
2, 4, 4, 15, 4 g/cm^3 for the five gamma sources
- lambda, cosmic-ray elevation scale height =
Not specified precisely; approximately 5 km for muons, 1 km for nuclear particles, and between these for the overall…
- kappa_c, kappa_sh, kappa_rho, correction factors for ceiling, shape, and density =
Approximately 2% per 10 cm concrete, up to 20% increase in R and 3% decrease in P for a 10x1 mm wafer, and M…
assumptions (6)
- domain assumption Secular equilibrium within decay chains, with radon-gap half-chains treated as distinct sources
- domain assumption Nominal specific activities of European building materials represent typical laboratory concrete
- domain assumption PARMA accurately models ground-level cosmic-ray spectra by species and elevation
- ad hoc to paper The two-step re-aimed simulation preserves the energy and direction distributions relevant to a small substrate
- domain assumption Thin superconducting films can be ignored, and the substrate dominates energy deposition
- domain assumption Energy deposited in the substrate is a sufficient proxy for qubit decoherence impact
Cite this review
Pith. "Pith review of Computed models of natural radiation backgrounds in qubits and superconducting detectors." pith.science (2026). https://pith.science/paper/FCF5BD2J
@misc{pith2026241116974,
author = {Pith},
title = {Pith review of: Computed models of natural radiation backgrounds in qubits and superconducting detectors},
year = {2026},
howpublished = {\url{https://pith.science/paper/FCF5BD2J}},
note = {Machine review of arXiv:2411.16974}
}
read the original abstract
Naturally occurring radiation backgrounds cause correlated decoherence events in superconducting qubits. These backgrounds include both gamma rays produced by terrestrial radioisotopes and cosmic rays. We use the particle-transport code Geant4 and the PARMA summary of the cosmic-ray spectrum to model both sources of natural radiation and to study their effects in the typical substrates used in superconducting electronics. We focus especially on three rates that summarize radiation's effect on substrates. We give analytic expressions for these rates, and how they depend upon parameters including laboratory elevation, substrate material, ceiling thickness, and wafer area and thickness. The modeled rates and the distribution of event energies are consistent with our earlier measurement of radiation backgrounds using a silicon thermal kinetic-inductance detector.
Figures
Reference graph
Works this paper leans on
-
[8]
J. W. Fowler, P. Szypryt, R. Bunker, E. R. Edwards, I. F. Florang, J. Gao, A. Giachero, S. F. Hoogerheide, B. Loer, H. P. Mumm, N. Nakamura, G. C. O’Neil, J. L. Orrell, E. M. Scott, J. Stevens, D. S. Swetz, B. A. VanDevender, M. Vissers, and J. N. Ullom, “Spectroscopic measurements and models of energy deposition in the substrate of quantum circuits by na...
-
[1]
The future of quantum computing with superconducting qubits,
S. Bravyi, O. Dial, J. M. Gambetta, D. Gil, and Z. Nazario, “The future of quantum computing with superconducting qubits,” Journal of Applied Physics , vol. 132, no. 16, p. 160902, 10 2022. [Online]. Available: https://doi.org/10.1063/5.0082975
-
[2]
Superconducting qubits: current state of play,
M. Kjaergaard, M. E. Schwartz, J. Braum ¨uller, P. Krantz, J. I.- J. Wang, S. Gustavsson, and W. D. Oliver, “Superconducting qubits: current state of play,” Annual Review of Condensed Matter Physics, vol. 11, no. 1, pp. 369–395, 2020. [Online]. Available: https://doi.org/10.1146/annurev-conmatphys-031119-050605
-
[3]
Superconducting microresonators: Physics and applica- tions,
J. Zmuidzinas, “Superconducting microresonators: Physics and applica- tions,” Annual Review of Condensed Matter Physics , vol. 3, no. 1, pp. 169–214, Mar. 2012
work page 2012
-
[4]
Reducing the impact of radioactivity on quantum cir- cuits in a deep-underground facility,
L. Cardani et al., “Reducing the impact of radioactivity on quantum cir- cuits in a deep-underground facility,” Nature Communications, vol. 12, no. 1, p. 2733, 2021
work page 2021
-
[5]
Correlated charge noise and relaxation errors in superconducting qubits,
C. D. Wilen, S. Abdullah, N. A. Kurinsky, C. Stanford, L. Cardani, G. D’Imperio, C. Tomei, L. Faoro, L. B. Ioffe, C. H. Liu, A. Opremcak, B. G. Christensen, J. L. DuBois, and R. McDermott, “Correlated charge noise and relaxation errors in superconducting qubits,” Nature, vol. 594, no. 7863, pp. 369–373, 2021
work page 2021
-
[6]
Resolving catastrophic error bursts from cosmic rays in large arrays of superconducting qubits,
M. McEwen et al. , “Resolving catastrophic error bursts from cosmic rays in large arrays of superconducting qubits,” Nature Physics, vol. 18, no. 1, pp. 107–111, 2022
work page 2022
-
[7]
Two- level-system dynamics in a superconducting qubit due to background ionizing radiation,
T. Thorbeck, A. Eddins, I. Lauer, D. T. McClure, and M. Carroll, “Two- level-system dynamics in a superconducting qubit due to background ionizing radiation,” PRX Quantum , vol. 4, p. 020356, 6 2023
work page 2023
Show all 22 references
-
[9]
The radioactivity of building materials,
C. Papastefanou, S. Stoulos, and M. Manolopoulou, “The radioactivity of building materials,” Journal of Radioanalytical and Nuclear Chemistry , vol. 266, pp. 367–372, 11 2005
2005
-
[10]
The effects of different types of concrete on population doses,
A. Suzuki, T. Iida, J. Moriizumi, and Y . Sakuma, “The effects of different types of concrete on population doses,” Radiation Protection Dosimetry, vol. 90, pp. 437–443, 9 2000. [Online]. Available: https://academic. oup.com/rpd/article-lookup/doi/10.1093/oxfordjournals.rpd.a033171
-
[11]
Natural radioactivity in building materials in the european union: a database and an estimate of radiological significance,
R. Trevisi, S. Risica, M. D’Alessandro, D. Paradiso, and C. Nuccetelli, “Natural radioactivity in building materials in the european union: a database and an estimate of radiological significance,” Journal of Environmental Radioactivity, vol. 105, pp. 11–20, 2 2012
2012
-
[12]
Kovler, Radioactive materials
K. Kovler, Radioactive materials . Elsevier, 2012, pp. 196–240
2012
-
[13]
An extremely energetic cosmic ray observed by a surface detector array,
The Telescope Array Collaboration, “An extremely energetic cosmic ray observed by a surface detector array,” Science, vol. 382, pp. 903–907, 11 2023
2023
-
[14]
Geant4—a simulation toolkit,
S. Agostinelli et al. , “Geant4—a simulation toolkit,” Nuclear Instru- ments and Methods in Physics Research Section A: Accelerators, Spec- trometers, Detectors and Associated Equipment , vol. 506, pp. 250–303, 7 2003
2003
-
[15]
Geant4 developments and applications,
J. Allison et al. , “Geant4 developments and applications,” IEEE Trans- actions on Nuclear Science , vol. 53, pp. 270–278, 2 2006
2006
-
[16]
Recent developments in Geant4,
——, “Recent developments in Geant4,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment , vol. 835, pp. 186–225, 11 2016
2016
-
[17]
TOPAS: An innovative proton Monte Carlo platform for research and clinical applications,
J. Perl, J. Shin, J. Sch ¨umann, B. Faddegon, and H. Paganetti, “TOPAS: An innovative proton Monte Carlo platform for research and clinical applications,” Medical Physics , vol. 39, pp. 6818–6837, 11 2012
2012
-
[18]
The TOPAS tool for particle simulation, a Monte Carlo simulation tool for physics, biology and clinical research,
B. Faddegon, J. Ramos-M ´endez, J. Schuemann, A. McNamara, J. Shin, J. Perl, and H. Paganetti, “The TOPAS tool for particle simulation, a Monte Carlo simulation tool for physics, biology and clinical research,” Physica Medica , vol. 72, pp. 114–121, 4 2020
2020
-
[19]
Quantum error correction below the surface code threshold,
Google Quantum AI, “Quantum error correction below the surface code threshold,” arXiv, 8 2024. [Online]. Available: https://arxiv.org/ abs/2408.13687
2024 arXiv
-
[20]
Xcom: Photon cross section database (version 1.5), nist standard reference database 8,
M. J. Berger, J. Hubbell, S. Seltzer, J. Chang, J. Coursey, R. Sukumar, D. Zucker, and K. Olsen, “Xcom: Photon cross section database (version 1.5), nist standard reference database 8,” U.S. National Institute of Standards and Technology, Tech. Rep., 2010. [Online]. Available:...
2010
-
[21]
Particle and heavy ion transport code system, PHITS, version 2.52,
T. Sato, K. Niita, N. Matsuda, S. Hashimoto, Y . Iwamoto, S. Noda, T. Ogawa, H. Iwase, H. Nakashima, T. Fukahori, K. Okumura, T. Kai, S. Chiba, T. Furuta, and L. Sihver, “Particle and heavy ion transport code system, PHITS, version 2.52,” Journal of Nuclear Science and Technol...
2013
-
[22]
Analytical model for estimating the zenith angle dependence of terrestrial cosmic ray fluxes,
T. Sato, “Analytical model for estimating the zenith angle dependence of terrestrial cosmic ray fluxes,” PLOS ONE , vol. 11, p. e0160390, 8 2016
2016
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.