REVIEW 5 major objections 5 minor 26 references
Simulation and measurement of Black Body Radiation background in a Transition Edge Sensor
T0 review · 5 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Black-body radiation from the room-temperature laboratory is the primary background source for a fiber-coupled transition-edge sensor, and a modular simulation reproduces the measured spectrum within uncertainties.
desk verdict A genuinely useful BBR simulation framework for fiber-coupled TES detectors, but the claimed 'within uncertainties' agreement with the 72-hour background measurement is overstated and needs a quantitative test. 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 central object is a modular radiative-transfer simulation built on the black-body spectrum $B(E,T)=2E^2/(h^3c^2(\exp(E/kT)-1))$ for a body at temperature $T$. For each optical path, the flux entering the fiber is integrated over the core area and the acceptance cone set by the numerical aperture, then multiplied by the wavelength-dependent transmissions of the fiber, two fiber-curling sections, and the transition-edge sensor optical stack; the result is convolved with a Gaussian energy response of width $\sigma_{1064\,\mathrm{nm}}$ and supplemented by a pile-up term with minimum resolvable time separation $\Delta t_{\min}=0.75\,\mu\mathrm{s}$. The same machinery produces the predictions at 11.3% and 5.3% energy resolution, and its uncertainty bands are dominated by the assumed curling diameters and by a 2 K temperature variation of the warm fiber end.
What would settle it
Measure a 72-hour extrinsic background spectrum with the fiber's path inside the cryostat mapped to better than a centimeter and a thermometer at the fiber's covered warm end; if the simulated black-body spectrum with those inputs does not reproduce the measured shape between 0.95 eV and 1.2 eV, or if the unexplained excess around 1.165 eV is independent of the measured fiber-end temperature, the central claim is falsified.
Extended reading notes
Core claim
The paper's central claim is that black-body radiation from the room-temperature laboratory, coupled through the optical fiber that guides light to the sensor, is the primary background source in the fiber-coupled setup, and that this source can be modeled quantitatively. The load-bearing mechanism is a black-body source spectrum $B(E,T)$, multiplied along every optical path by wavelength-dependent transmissions—fiber, fiber curling, and the sensor's band-pass absorbance—and finally convolved with the sensor's Gaussian energy resolution and a pile-up term. The simulated spectrum is consistent with the measured extrinsic background between roughly 0.95 eV and 1.2 eV, and the simulation correctly predicts the size of the improvement obtained when the analysis energy resolution is improved from 11.3% to 5.3%: the background rate in the 1064 nm signal window falls by an order of magnitude, to about $10^{-4}$ counts per second. The paper also notes an excess of events in the one-$\sigma$ region around 1.165 eV that the black-body model does not explain, which it attributes to non-black-body events surviving the analysis cuts.
Load-bearing premise
The analysis assumes that the pulse-shape cuts—fitted rise time, decay time, and reduced chi-squared, all calibrated on 1064 nm laser pulses—reject non-photon events without biasing the reconstructed photon energies; the paper itself reports an excess of events in the one-sigma signal region, indicating that this assumption is at least partially violated exactly where the signal would be.
Editorial extensions
If this is right
- If black-body radiation is indeed the dominant background, then reducing the temperature of the fiber's warm end or shielding it inside the cryostat attacks the dominant background source; the simulation estimates a factor-of-five rate reduction when the lab temperature drops from 295 K to 283 K.
- The factor-of-two improvement in energy resolution, from 11.3% to 5.3%, lowers the measured background rate in every one-sigma window around 1064 nm by an order of magnitude, in agreement with the simulation.
- At 5.3% energy resolution the simulated upper-limit black-body rate in the [0, 3-sigma] window reaches $5\times10^{-6}$ counts per second, which would meet the experiment's target, while the measured value of $6.9\times10^{-5}$ counts per second shows the gap that remains.
- Fiber curling suppresses both the low-energy tail of the black-body spectrum and pile-up pairs that could mimic 1064 nm photons, making the in-cryostat fiber geometry a controllable background-rejection handle.
- The residual excess in the one-sigma signal region implies that even a perfect black-body background model will not suffice; those non-black-body events must be identified and suppressed separately.
Reading between the lines
- A consequence the authors leave implicit: the unexplained events in the one-sigma signal region should be nearly independent of the fiber-end temperature, unlike the black-body component, so a temperature-sweep measurement would separate the two backgrounds without changing the analysis.
- Beyond this specific detector, the same radiative-transfer machinery applies to other fiber-coupled cryogenic photon counters; the uncertainty analysis here suggests that precisely documenting the fiber path inside a cryostat is the cheapest way to sharpen any such background prediction.
- Because the black-body rate at 1.165 eV falls steeply with temperature, actively cooling just the fiber feedthrough could push the black-body rate well below the factor-of-five gain the paper quotes, an engineering route the paper leaves open.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports on a tungsten Transition Edge Sensor (TES) being developed for single-photon detection in the ALPS II light-shining-through-a-wall experiment. The main technical contribution is a simulation framework that computes the expected black-body radiation (BBR) background propagating through an optical fiber to the TES, including fiber transmission, fiber curling losses, TES absorbance, energy resolution, and pile-up. The framework is exercised for the specific geometry of an extrinsics background measurement (warm fiber end covered, fiber curled inside the cryostat, 72 h of data). The authors compare the simulated BBR spectrum with measured spectra from the extrinsics run and claim that the background is consistent with BBR as the primary contributor and that the simulation reproduces the spectral distribution within uncertainties. They also show that improving the energy resolution from 11.3% to 5.3% reduces the predicted background in the 1064 nm signal region by about an order of magnitude, although the measured rates in that region still exceed the ALPS II requirement of 7.7e-6 cps.
Significance. If the BBR simulation framework is validated, it would be a useful tool for designing background-reduction strategies for the ALPS II single-photon detector and for similar cryogenic single-photon experiments. The framework is modular, and most input parameters are taken from independent measurements (manufacturer fiber-loss data, spectrometer measurements of curling transmission, reflectance data from a thesis, and calibration-derived energy resolution), which is a strength. The paper also documents a concrete reduction in measured extrinsic background achieved through analysis improvements. However, the central claim of quantitative agreement between simulation and measurement is currently not supported by the quantitative comparison in the paper: in the most important 1064 nm signal region, the measured rates exceed the simulated upper limits by factors of 3 to 14, and the paper itself attributes part of the excess to non-BBR events. The manuscript therefore requires substantial revision of its claims and a more careful uncertainty treatment before the validation claim can be accepted.
major comments (5)
- [Abstract; Section V, Tables IV and V] The abstract and Section V state that the simulation reproduces the observed background 'within uncertainties', but the numbers in Tables IV and V do not support this for the signal region. For the [-1,1] sigma window, the measured pulse-height rate is 1.2e-4 cps while the simulated upper limit is 4.6e-5 cps; for the [0,3] sigma window the measured rate is 6.9e-5 cps versus a simulated upper limit of 5e-6 cps, a factor of about 14. The paper should quantify the level of agreement (or disagreement) explicitly, state which energy ranges the 'within uncertainties' claim applies to, and restrict the validation claim accordingly.
- [Section IV, Fig. 8; Section V, Figs. 10 and 11] The simulation uncertainty band is constructed from two ad hoc variations (a +/-2 K temperature change and an assumed extra half-loop with an estimated diameter), rather than from a measured or propagated uncertainty on the fiber path inside the cryostat. In addition, the measured data are presented without error bars and no chi-square, Kolmogorov-Smirnov, or similar quantitative comparison is reported. The authors should either provide a more principled uncertainty estimate (e.g., varying the curling parameters within independently measured ranges, including the 10 cm vs 11.2 cm approximation) or explicitly frame Fig. 11 as a qualitative comparison.
- [Section V, cuts on pulse-shape parameters] The pulse-shape cuts (rise time, decay time, reduced chi-squared) are defined from 1064 nm calibration pulses and are assumed to be energy independent, but this assumption is load-bearing for interpreting the measured spectrum as a photon energy spectrum. The paper itself notes an excess of events in the 1-sigma signal region that is attributed to non-BBR events, indicating that the cuts do not fully reject non-photon events in exactly the region relevant for ALPS II. The authors should test or justify the energy independence of the cuts (for example, by injecting known-energy sources across the measured range or by comparing the measured spectrum before and after cuts) and discuss how residual non-photon events affect the comparison with the BBR simulation.
- [Section V, Table V and Section IV, Tables III and IV] The simulated rates in Tables III and IV use an analysis efficiency that only includes the sigma-window acceptance, while the measured rates in Table V include an additional 98.4% cut acceptance and 90% system detection efficiency. The comparison between Tables IV and V is therefore not apples-to-apples. The authors should state explicitly whether the simulated rates should be multiplied by the same acceptance and efficiency factors before comparison, or provide a direct corrected comparison in the text.
- [Section III D, energy-resolution assumption] Equation 6 and the surrounding text assume a constant energy resolution sigma(E) = sigma_1064nm over the range [0 eV, 1.5 eV], citing previous measurements and other groups. Given that the low-energy part of the spectrum is not well described by the simulation and that the TES response may be nonlinear at low energies, the authors should provide quantitative evidence for this assumption in their own setup or explicitly discuss its impact on the comparison in Fig. 11.
minor comments (5)
- [Fig. 9 caption] The caption appears to have inconsistent labeling: subfigures (b), (c), and (d) are referenced, but the text seems to refer to (c) twice and (e) twice; please correct the panel labels and their descriptions.
- [Section IV, Fig. 8] The figure legend states that the upper limit is produced by varying the temperature by 2 K and the lower limit by adding a non-accounted bending, but the text in Section IV also mentions a sharper quarter-loop; please clarify in the legend which exact variations correspond to the band edges.
- [Figs. 10 and 11] The measured spectral points are plotted without error bars, which is especially problematic because the rates in Table V are derived from a single 72-hour run; adding Poisson error bars would help the reader judge the significance of the observed discrepancies.
- [Section III B, fiber transmission] The extrapolation of fiber loss outside the displayed wavelength range is described as 'linear', but the linearity is stated without further justification; a brief comment on the validity range would be useful.
- [Section IV, Eq. 11] The text says 'The accuracy of Eq. 11 is affected mainly by the fiber curling component and the temperature', but the actual contributions of the fiber-loss data, the TES reflectance approximation, and the numerical aperture are not quantified; a short sensitivity discussion would strengthen the uncertainty treatment.
Circularity Check
No significant circularity: the BBR simulation is a forward model built from independent measured and literature inputs, and it is not fitted to the extrinsics spectrum.
full rationale
The paper's central derivation chain is not circular. The simulated BBR spectrum is computed from Planck's law (Eq. 1), the fiber numerical aperture and core area (Eq. 10), manufacturer fiber loss data [18], a spectrometer-measured fiber-curling transmission fitted to an error-function form (Eq. 4 and Table II), a TES reflectance curve approximated from the independent thesis data of [17], and an energy resolution taken from prior data analyses [13,25]. None of these parameters is adjusted to match the 72-hour extrinsics histogram of Section V, so the simulation is a genuine forward prediction rather than a fit renamed as a prediction. The qualitative comparison in Fig. 11 and the abstract's 'reproduces, within uncertainties' claim are quantitatively under-supported (no chi-square or KS statistic is reported), and the authors themselves concede that the low-energy part is not described and that the excess in the 1-sigma signal region 'indicates that these events are not due to BBR.' However, an overstated agreement claim is a validation and presentation concern, not circularity. The energy-resolution reduction result is also not a fitted-input-called-prediction: reducing sigma_1064nm in Eq. 6 narrows the Gaussian detector response, and because the BBR spectrum rises steeply toward low energies, the rate in sigma-scaled windows around 1064 nm decreases by roughly an order of magnitude. This is a consequence of the model, not an input relabeled as an output. The self-citations used in the paper, e.g., [13] and [25] for energy resolution and system detection efficiency, supply independently measured detector parameters that do not assume the BBR hypothesis and are not used to enforce agreement with the background data. No equation in the paper reduces to another by construction, and no load-bearing claim is justified solely by a self-citation chain. Therefore no circularity is identified.
Assumptions & free parameters
free parameters (7)
- TES reflectance model parameters (lambda1, b1, a1, lambda2, b2, a2) =
700 nm, 0.55, 0.01/nm, 2100 nm, 0.2, 0.002/nm
- Fiber one-loop curling fit parameters (a, lambda0 per diameter) =
d=12.0 mm: a=9.71e-3/nm, lambda0=1231 nm; d=112 mm: a=1.75e-2/nm, lambda0=1768 nm (Table II)
- Room temperature T =
295 K with +2 K upper variation
- Fiber length inside cryostat =
4 m
- Fiber bending geometry inside cryostat =
4 loops of 11.2 cm diameter plus 1/4 loop of 2.5 cm; lower bound adds 1/2 loop of 4.2 cm and 1/4 loop of 1.2 cm
- Energy resolution sigma_1064nm =
0.13 eV (11.3%) and 0.061 eV (5.3%)
- Pile-up time window Delta_t_min =
0.75 microseconds
assumptions (8)
- domain assumption Black-body emission follows Planck's law and Lambert's cosine law (Eq. 1)
- domain assumption Fiber at room temperature re-emits absorbed BBR inside the fiber, equivalent to BBR production at the cryostat feedthrough
- domain assumption Fiber accepts all photons with incidence angle below theta_max set by NA and rejects all others; component transmission is angle-independent
- domain assumption TES energy response is a Gaussian with constant sigma(E) = sigma_1064nm over 0 to 1.5 eV
- domain assumption Pulse-shape parameters (rise time, decay time, reduced chi-squared) are independent of photon energy, so 1064 nm calibration cuts apply to all background events
- domain assumption TES response is linear: pulse integral and pulse height are proportional to photon energy
- domain assumption BBR coupling at fiber junctions inside the cryostat is negligible because of the lower temperature
- domain assumption Measured extrinsics environment is a black body at 295 K with the fiber end covered by black cloth
Cite this review
Pith. "Pith review of Simulation and measurement of Black Body Radiation background in a Transition Edge Sensor." pith.science (2026). https://pith.science/paper/TQS454BJ
@misc{pith2026250508555,
author = {Pith},
title = {Pith review of: Simulation and measurement of Black Body Radiation background in a Transition Edge Sensor},
year = {2026},
howpublished = {\url{https://pith.science/paper/TQS454BJ}},
note = {Machine review of arXiv:2505.08555}
}
abstract
The Any Light Particle Search II (ALPS II) experiment at DESY, Hamburg, is a Light-Shining-through-a-Wall (LSW) experiment aiming to probe the existence of axions and axion-like particles (ALPs), which are candidates for dark matter. Data collection in ALPS II is underway utilizing a heterodyne-based detection scheme. A complementary run for confirmation or as an alternative method is planned using single photon detection, requiring a sensor capable of measuring low-energy photons ($1064\,\mathrm{nm}$, $1.165\,\mathrm{eV}$) with high efficiency (higher than $50\,\%$) and a low background rate (below $7.7\cdot10^{-6}\,\mathrm{cps}$). To meet these requirements, we are investigating a tungsten Transition Edge Sensor (TES) provided by NIST, which operates in its superconducting transition region at millikelvin temperatures. This sensor exploits the drastic change in resistance caused by the absorption of a single photon. We find that the background observed in the setup with a fiber-coupled TES is consistent with Black Body Radiation (BBR) as the primary background contributor. A framework was developed to simulate BBR propagation to the TES under realistic conditions. The framework not only allows the exploration of background reduction strategies, such as improving the TES energy resolution, but also reproduces, within uncertainties, the spectral distribution of the observed background. These simulations have been validated with experimental data, in agreement with the modeled background distribution, and show that the improved energy resolution reduces the background rate in the $1064\,\mathrm{nm}$ signal region by one order of magnitude, to approximately $10^{-4}\,\mathrm{cps}$. However, this rate must be reduced further to meet the ALPS II requirements.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
Ringwald, Axions and axion-like particles, (2014), arXiv:1407.0546 [hep-ph]
A. Ringwald, Axions and axion-like particles, (2014), arXiv:1407.0546 [hep-ph]
arXiv 2014
-
[2]
A. A. Anselm, Arion ↔ Photon Oscillations in a Steady Magnetic Field. (In Russian), Yad. Fiz. 42, 1480 (1985)
work page 1985
-
[3]
K. Van Bibber, N. R. Dagdeviren, S. E. Koonin, A. K. Kerman, and H. N. Nelson, Proposed experiment to pro- duce and detect light pseudoscalars, Phys. Rev. Lett. 59, 759 (1987)
work page 1987
-
[4]
R. B¨ ahre, B. D¨ obrich, J. Dreyling-Eschweiler, S. Ghaz- aryan, R. Hodajerdi, D. Horns, F. Januschek, E. A. Kn- abbe, A. Lindner, D. Notz, A. Ringwald, J. E. v. Seggern, R. Stromhagen, D. Trines, and B. Willke, Any light par- ticle search II — technical design report, Journal of In- strumentation 8 (09), T09001–T09001
-
[5]
I. G. Irastorza and J. Redondo, New experimental ap- proaches in the search for axion-like particles, Progress in Particle and Nuclear Physics 102, 89–159 (2018)
work page 2018
-
[6]
R. Ballou et al. (OSQAR), New exclusion limits on scalar and pseudoscalar axionlike particles from light shining through a wall, Phys. Rev. D 92, 092002 (2015), arXiv:1506.08082 [hep-ex]
arXiv 2015
-
[7]
316720, Proceedings, 11th Patras Workshop on Axions, WIMPs and WISPs (Axion-WIMP 2015), 11th Patras Workshop on Axions, WIMPs and WISPs, Zaragoza (Spain), 22 Jun 2015 - 26 Jun 2015 (Verlag Deutsches Elektronen-Synchrotron, Hamburg, 2015) pp. 1–251
work page 2015
-
[8]
M. D. Ortiz et al., Design of the ALPS II optical system, Phys. Dark Univ. 35, 100968 (2022), arXiv:2009.14294 [physics.optics]
arXiv 2022
Show all 26 references
-
[9]
Hallal, G
A. Hallal, G. Messineo, M. D. Ortiz, J. Gleason, H. Hollis, D. Tanner, G. Mueller, and A. Spector, The heterodyne sensing system for the alps ii search for sub-ev weakly interacting particles, Physics of the Dark Universe 35, 100914 (2022)
2022
-
[10]
Shah, K.-S
R. Shah, K.-S. Isleif, F. Januschek, A. Lindner, and M. Schott, TES detector for ALPS II, Pro- ceedings of The European Physical Society Confer- ence on High Energy Physics — PoS(EPS-HEP2021) 10.22323/1.398.0801 (2022)
2022 doi
-
[11]
Irwin and G
K. Irwin and G. Hilton, Transition-edge sensors, in Cryo- genic Particle Detection (Springer Berlin Heidelberg, Berlin, Heidelberg, 2005) pp. 63–150
2005
-
[12]
Shah, K.-S
R. Shah, K.-S. Isleif, F. Januschek, A. Lindner, and M. Schott, Characterising a Single-Photon Detector for ALPS II, Journal of Low Temperature Physics 209, 355 (2022)
2022
-
[13]
J. A. Rubiera Gimeno, K.-S. Isleif, F. Januschek, A. Lind- ner, M. Meyer, G. Othman, M. Schott, R. Shah, and L. Sohl, The TES detector of the ALPS II experiment, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associate...
2023
-
[14]
A. J. Miller, A. Lita, D. Rosenberg, S. Gruber, and S. Nam, Superconducting photon number resolving de- tectors: Performance and promise, in Proc. 8th Int. Conf. Quantum Communication, Measurement and Computing (QCMC’06) (2007) pp. 445–450
2007
-
[15]
Greiner, L
W. Greiner, L. Neise, and H. St¨ ocker, Thermodynam- ics and statistical mechanics, Classical theoretical physics (Springer, 1995)
1995
-
[16]
A. E. Lita, B. Calkins, L. A. Pellouchoud, A. J. Miller, and S. Nam, Superconducting transition-edge sensors op- timized for high-efficiency photon-number resolving de- tectors, in Advanced Photon Counting Techniques IV , Vol. 7681, edited by M. A. Itzler and J. C. Campbell, I...
-
[17]
J. Dreyling-Eschweiler, A superconducting mi- crocalorimeter for low-flux detection of near-infrared single photons , Dr., University of Hamburg, Hamburg (2014), university of Hamburg, Diss., 2014
2014
-
[18]
M. A. Arshad, A. Hartung, and M. J¨ ager, A stimulated stokes raman scattering-based approach for continuous wave supercontinuum generation in optical fibers, Laser Physics Letters 16, 035108 (2019)
2019
-
[19]
Marcuse, Influence of curvature on the losses of doubly clad fibers, Appl
D. Marcuse, Influence of curvature on the losses of doubly clad fibers, Appl. Opt. 21, 4208 (1982)
1982
-
[20]
Smirnov, Y
K. Smirnov, Y. Vachtomin, A. Divochiy, A. Antipov, and G. Goltsman, Dependence of dark count rates in super- conducting single photon detectors on the filtering effect of standard single mode optical fibers, Applied Physics Express 8, 022501 (2015)
2015
-
[21]
Cabrera, R
B. Cabrera, R. M. Clarke, P. Colling, A. J. Miller, S. Nam, and R. W. Romani, Detection of single infrared, optical, and ultraviolet photons using superconducting transition edge sensors, Applied Physics Letters 73, 735 (1998)
1998
-
[22]
A. Lita, A. Miller, and S. Nam, Energy collection effi- ciency of tungsten transition-edge sensors in the near- infrared, Journal of Low Temperature Physics 151, 125 (2008)
2008
-
[23]
P. Li, J. Zhong, W. Zhang, Z. Wang, K. Zhou, W. Miao, Y. Ren, J. Li, Q. Yao, and S. Shi, Multi-color photon detection with a single superconducting transition-edge sensor, Nuclear Instruments and Methods in Physics Re- search Section A: Accelerators, Spectrometers, Detectors a...
2023
-
[24]
Eckart and F
C. Eckart and F. R. Shonka, Accidental coincidences in 13 counter circuits, Phys. Rev. 53, 752 (1938)
1938
-
[25]
J. A. Rubiera Gimeno, F. Januschek, K.-S. Isleif, A. Lind- ner, M. Meyer, G. Othman, C. Schwemmbauer, and R. Shah, A TES system for ALPS II - Status and Prospects, PoS EPS-HEP2023, 567 (2023)
2023
-
[26]
Meyer, G
M. Meyer, G. Othman, K. Isleif, F. Januschek, A. Lindner, J. A. Rubiera Gimeno, C. Schwemmbauer, M. Schott, and R. Shah, A first application of ma- chine and deep learning for background rejection in the ALPS II TES detector, 17th Workshop on Axions, WIMPs and WISPs 10.1002/an...
2023 arXiv
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.