Pith. sign in

REVIEW 3 major objections 5 minor 12 references

Performance of the LABDOS01 spectrometer in dosimetric measurements

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A silicon diode spectrometer can gauge biological dose to within 15%.

desk verdict Solid instrument characterization; the H*(10) formula's low-LET term is unvalidated at ground level, so the 'better than 15%' claim is a bit strong. read the letter →

arxiv 2412.01515 v1 pith:SG4FZITB submitted 2024-12-02 physics.ins-det

Silvia Vernetto , Carlo Francesco Vigorito , Martin Kakona , Roman Dvorak , Satoshi Kodaira , Stefano Cirilli , Fabio Romanelli , Daniele Savietto
show 1 more author
Alba Zanini
This is my paper · ORCID
classification physics.ins-det PACS 29.40.Wk
keywords dosimetrycosmicrayssilicondiodespectrometerLiulin-typedetectorambientdoseequivalentLETspectroscopyCARI-7Amixedradiationfields
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

LABDOS01 is a small, low-cost silicon diode spectrometer that records the spectrum of energy deposited by cosmic-ray secondaries in a thin silicon layer. Along two airline routes covering rigidity cutoffs from 2.7 to 17.2 GV, its measured absorbed dose in silicon matches CARI-7A simulations with ratios between 1.00 and 1.16. By splitting the deposited energy at 1.5 MeV into low- and high-LET contributions, the paper derives a simple two-term formula for the ambient dose equivalent $H^{*}(10)$: $H^{+} = 1.19\,D_L + 10.3\,D_H$. Applied to an independent flight, the formula reproduces the simulated $H^{*}(10)$ within statistical errors, and the high-LET/neutron part matches Rem-counter measurements at three ground sites to better than 15%. The paper argues this makes LABDOS01 a practical monitor for long-term environmental dose at extreme-altitude and Antarctic sites.

What carries the argument

The load-bearing object is the deposited-energy spectrum measured by the LABDOS01 PIN diode (2 cm² area, nominal 300 µm thickness, effective depletion thickness 210 µm), sampled every ~12 s over roughly 60 keV to 7 MeV. From that spectrum, the absorbed dose in silicon is split at $E_{\mathrm{dep}} = 1.5$ MeV into a low-LET dose $D_L$ and a high-LET dose $D_H$; the paper treats the high-LET part as essentially neutron-driven, following earlier Liulin-type work. Two constant conversion ratios, $\bar{R}_L = 1.19$ and $\bar{R}_H = 10.3$, are derived by matching $D_L$ and $D_H$ to the corresponding low- and high-LET components of $H^{*}(10)$ computed with CARI-7A on four flight segments, and are then used linearly in Eq. (2) to convert any LABDOS01 spectrum into $H^{*}(10)$.

What would settle it

Take the same LABDOS01 unit to a location with a distinctly different radiation field, such as a polar flight with rigidity cutoff below 2.7 GV or a ground site with a strong local gamma background, measure $D_L$ and $D_H$, compute $H^{+} = 1.19\,D_L + 10.3\,D_H$, and compare it with a reference $H^{*}(10)$ instrument such as a tissue-equivalent proportional counter or Bonner-sphere spectrometer; a deviation beyond about 15% that tracks the neutron-to-photon ratio would falsify the constant-coefficient claim.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that a silicon diode of the Liulin type, calibrated only against $\alpha$ sources and accelerated ions, can act as an ambient-dose dosimeter in mixed cosmic-ray fields without a full response model. The instrument's measured absorbed dose in silicon agrees with the standard CARI-7A transport code to within 16% along flight routes that vary in altitude and geomagnetic cutoff. More specifically, the paper establishes an empirical identity: separating measured dose by a threshold at 1.5 MeV deposited energy, the ambient dose equivalent is $H^{*}(10) \approx 1.19\,D_L + 10.3\,D_H$, where $D_L$ and $D_H$ are the absorbed doses from particles depositing below and above that threshold. The coefficients are fixed by cross-correlating four calibration flight segments with CARI-7A, and the formula is then checked against an independent flight and against ground neutron doses, supporting a claimed accuracy better than 15% over the tested range.

Load-bearing premise

The two conversion factors are assumed to stay constant across altitudes, latitudes, and ground environments, but they are derived from only four calibration flight segments and then applied more widely.

Editorial extensions

If this is right

  • If the 15% accuracy holds, LABDOS01 can serve as a low-cost continuous dosimeter at high-altitude and polar sites such as Chacaltaya and Concordia, where it is already installed.
  • The method gives a way to estimate biologically relevant dose from a single silicon spectrum, without needing a separate neutron detector or tissue-equivalent counter.
  • Because the device records full spectra, it can also track dose changes during space-weather events, such as ground-level enhancements, rather than only reporting time-integrated dose.
  • The empirical formula is validated for cruise altitudes and rigidity cutoffs from 2.7 to 17.2 GV; applying it at other altitudes, ground sites, or rigidity cutoffs would require new validation.
  • The high-LET/neutron branch of the formula agrees with Wendi-2 Rem-counter measurements at three ground sites within 15%, so the neutron-dose channel can be monitored with the same detector.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the 1.19 and 10.3 coefficients are universal, archival data from older Liulin-type flight dosimeters could be reprocessed into $H^{*}(10)$ estimates even where no reference dosimeter was flown; this is a testable extension the paper does not attempt.
  • The effective diode thickness is about 210 µm, while CARI-7A assumes 300 µm; the paper estimates a few percent bias, but a dedicated simulation would be needed to know whether the conversion coefficients shift for other detector thicknesses.
  • A decisive extension would be an intercalibration flight during a ground-level enhancement or on a polar route with rigidity cutoff below 2.7 GV, where the neutron spectrum may soften or harden enough to break the 1.5 MeV LET split.
  • The paper's ground validation covers only the high-LET/neutron branch; the low-LET coefficient 1.19 has not yet been compared against a reference instrument for gamma- and electron-dominated ground fields.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper reports on the LABDOS01, a silicon diode-based spectrometer (Liulin-type) intended for environmental dosimetry in mixed radiation fields. The instrument was calibrated using radioisotope alpha sources and HIMAC ion beams, then flown on two long-range routes (Milan–Christchurch and Milan–La Paz) where its absorbed dose in silicon (DSi) was compared with CARI-7A simulations. The measured DSi rates follow the simulated curves, with integrated segment ratios between 1.00 and 1.16. By splitting the deposited-energy spectrum into low-LET (Edep < 1.5 MeV) and high-LET (1.5–7 MeV) parts, the authors derive an empirical formula, H+ = 1.19 DL + 10.3 DH (Eq. 2), to estimate the ambient dose equivalent H*(10). The coefficients are calibrated on the Milan–Christchurch flight segments using CARI-7A as reference, then applied to the Milan–La Paz flight and to three ground-level sites (Turin, Concordia, Chacaltaya). The ground check validates only the high-LET term against neutron Rem-counter measurements. The paper concludes that H+(10) can be estimated with an accuracy better than 15% over a wide range of altitudes and rigidity cutoffs.

Significance. If the central claims are supported, the paper provides a useful low-cost, open-source detector for monitoring cosmic-ray dose at remote high-altitude and polar sites, and a practical empirical route from silicon absorbed dose to ambient dose equivalent. The strengths include the use of two flight campaigns over a wide rigidity-cutoff range, a clear calibration procedure with independent accelerator data, and a direct ground comparison of the neutron component against Wendi-2 Rem-counters. The paper is also careful to state several limitations (e.g., missing tracking data, different-day tracking for one segment, and the need for ad hoc simulations for full ground validation). However, the headline accuracy claim ('better than 15%') is not yet fully established for the total H*(10) at the ground-deployments where the instrument is actually operated, because the low-LET conversion is calibrated only in cruise-altitude cosmic-ray fields and the ground validation does not independently test that component.

major comments (3)
  1. [Sec. 4.1, Eq. (2), Table 7] The conversion coefficients RL and RH in Eq. (2) are calibrated only on four cruise-altitude flight segments (Tables 4 and 5), at altitudes of 10–12.5 km and rigidity cutoffs of 2.7–17.2 GV. The formula is then applied in Table 7 to the ground-level sites (Turin 240 m, Concordia 3233 m, Chacaltaya 5240 m), where the radiation field differs materially: the low-LET component includes local radioactivity and muons, and the neutron spectrum is softer. The ground validation compares only the high-LET product RH·DH against Wendi-2 neutron H*(10); the low-LET product RL·DL is not measured against any independent reference at ground. The paper itself states that ground validation of the formula 'would require ad hoc simulations... but this is beyond the scope of this work.' Consequently, the abstract's and Section 5's claim that H*(10) is obtained 'with an accuracy better than 15%' at the deployed sites is not supported for the total H+; only the high-LET contribution is checked.
  2. [Sec. 2, Tables 2–7] The paper quotes a systematic uncertainty of ±8% in the estimated active mass of the PIN diode, but this uncertainty is not propagated into any of the reported integrated doses, the ratios RD, or the H+ values; only statistical errors are shown in Tables 2–7. As a result, statements such as 'the measured doses are slightly larger than the simulation results by a factor of less than 10%' (Sec. 3.1) and 'within the quoted statistical and systematic errors' (Sec. 5) lack a quantitative basis. The authors should either propagate the ±8% uncertainty through all derived numbers, or explicitly show which conclusions are insensitive to it. One important subtlety is that the mass uncertainty cancels in the ratios RL and RH if the same detector is used for both calibration and measurement, but not when the formula is applied to the second unit (LABDOS01-EF) in Table 6, whose active mass may differ.
  3. [Sec. 4.1, Table 6] The out-of-sample check of Eq. (2) is limited to a single route (Milan–La Paz). One segment, VVI-LPZ, shows a relative deviation of 27% between H+ (0.70±0.16 µSv) and the CARI-7A H*(10) (0.55 µSv), which is large compared with the claimed 15% accuracy. Although the deviation is within about 1.9σ given the statistical error, the segment is only 0.26 h long, and the simulated values for the preceding MAD-VVI leg rely on tracking data from a different day. The statement that the H+ values are in 'excellent agreement' with the simulated H*(10) is therefore too strong; the data demonstrate consistency within statistical uncertainties, but do not firmly corroborate a 15% accuracy bound. Additional out-of-sample segments or a longer independent dataset would be needed to substantiate the claimed accuracy.
minor comments (5)
  1. [Sec. 4.1, last paragraph] The phrase 'cruise altitudes and Rc lower than 2.7 GV' appears to be a typo; the data cover Rc down to 2.7 GV, not lower than 2.7 GV. Please correct to 'Rc down to 2.7 GV'.
  2. [Sec. 3.2, Fig. 7 and Table 3] The airport code for La Paz is written inconsistently as 'VVI-LPX' in the main text and 'VVI-LPZ' in Fig. 7 and Table 3. Please standardize (the common code is VVI-LPZ).
  3. [Sec. 5, first paragraph] There is a typo: 'the capability of the device to to disentangle' should read 'the capability of the device to disentangle'.
  4. [Throughout] The instrument is variably called 'LABDOS01' and 'LABDOS1' (e.g., in Section 5); please use one consistent name.
  5. [Fig. 4 caption] The caption states 'h=11000 km'; this should be 'h=11000 m'.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Eq. (2) is an explicitly empirical calibration against an external Monte Carlo model, and the Milan-La Paz flight is a genuine out-of-sample transfer test.

full rationale

The paper's central derivation is an empirical calibration chain, not a self-referential one. The LABDOS01 absorbed dose DSi is measured from the silicon deposited-energy spectrum (Eq. 1), and compared to CARI-7A, an external FAA-standard code (Copeland 2017). The H*(10) formula H+ = 1.19 DL + 10.3 DH (Eq. 2) is obtained by fitting the two ratios RL = H*_L/DL and RH = H*_H/DH to CARI-7A reference values over four Milan-Christchurch flight segments (Tables 4 and 5). This is openly stated: 'We derive here a simplified empirical approach ... by cross-correlating the data of the four flight segments from Milan to Christchurch with the corresponding CARI-7A calculations used as reference.' The subsequent application to the Milan-La Paz flight is explicitly out-of-sample: 'As a check, we apply this method to the LABDOS01-2 data recorded on the Milan-La Paz flight, not used for calibration purposes.' The coefficients are not refit to this route, so the agreement in Table 6 is not forced by construction. The ground-level comparison validates the high-LET (neutron) term against independent Thermo Wendi-2 Rem-counter measurements (Table 7), and the paper candidly states that full ground validation of the low-LET term 'would require ad hoc simulations ... beyond the scope of this work.' No load-bearing self-citation is used: the cited Ploc, Dachev, Spurny, and Zhang works are external, and the authors' own SAMADHA citation is only a project description. There is no uniqueness theorem, no ansatz smuggled through citation, and no renaming of a known result. The main caveats—constant conversion ratios across altitudes and ground/aircraft environments, and the 27% deviation in the short VVI-LPZ segment—are limitations on the accuracy claim, not circularity. The derivation is therefore self-contained against external benchmarks, and no circular step can be exhibited from the paper's equations or citations.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central empirical formula rests on two fitted coefficients and four domain assumptions. RL and RH are calibrated to CARI-7A simulated H*(10) values for the Milan-Christchurch flight, so the reference model is assumed correct. The method also assumes the 1.5 MeV threshold cleanly separates low and high LET signals and that the conversion ratios are transferable outside the calibration data. No new entities are introduced.

free parameters (4)
  • RL (low-LET conversion coefficient) = 1.19
    Fit to the Milan-Christchurch data by comparing measured DL with CARI-7A H*_L for four flight segments (Table 4).
  • RH (high-LET conversion coefficient) = 10.3
    Fit to the Milan-Christchurch data by comparing measured DH with CARI-7A H*_H for four flight segments (Table 5).
  • Low/high LET threshold = 1.5 MeV
    Chosen by hand, acknowledged as having no precise physical meaning (Section 4.1), separates low and high LET deposits.
  • Upper energy bound = 7 MeV
    Chosen to exclude saturated ADC channels (Section 4.1).
assumptions (5)
  • domain assumption CARI-7A provides accurate reference values for absorbed dose in silicon and H*(10) in the flight conditions.
    Used as ground truth for both DSi validation (Section 3) and H+ calibration (Section 4.1).
  • domain assumption The high-LET component of H*(10) is almost exclusively due to neutrons.
    Section 4.1, citing Dachev et al. (2020); used to set H*_H = H*_n.
  • ad hoc to paper The conversion ratios H*_L/DL and H*_H/DH are constant across the application range.
    Needed to generalize Eq. (2) beyond the four calibration segments; not derived from first principles.
  • ad hoc to paper The two-bin split at 1.5 MeV adequately separates low and high LET contributions.
    Assumed in Section 4.1; the threshold has no physical basis.
  • domain assumption Energy calibration from alpha sources and the HIMAC ion beam is valid for all particle types and deposited energies.
    Section 2 calibrates using 239Pu, 241Am, and 572.8 MeV He ions; assumes linear response over the full range.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Performance of the LABDOS01 spectrometer in dosimetric measurements." pith.science (2026). https://pith.science/paper/SG4FZITB

@misc{pith2026241201515,
  author       = {Pith},
  title        = {Pith review of: Performance of the LABDOS01 spectrometer in dosimetric measurements},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SG4FZITB}},
  note         = {Machine review of arXiv:2412.01515}
}
abstract

This paper describes the performance of the LABDOS01, a silicon diode-based spectrometer suitable for dose measurements in mixed radiation fields. The instrument is currently being used in two high-altitude environmental dose monitoring projects: SAMADHA (South Atlantic Magnetic Anomaly Dosimetry at High Altitude) at Chacaltaya (Bolivia, 5240 m a.s.l.) and CORDIAL (COsmic Rays Dosimetry In Antarctic Latitudes) at the Concordia station (Antarctica, 3233 m a.s.l.). Before installing two of these devices at the measurement sites, the detectors were tested on flight routes covering a wide range of geomagnetic latitudes. The collected dosimetric data were compared with the expectations derived by the CARI-7A software, which provides the absorbed dose rate in silicon due to cosmic ray secondaries at a given position on the Earth. The measured dose rates along the flights at variable altitude and rigidity cutoff agree well with the simulated ones. By analyzing the spectrum of the energy deposited in the silicon layer, we derive an empirical method to approximately evaluate the ambient dose equivalent $H^{*}(10)$, a quantity directly related to the biological damage caused by environmental radiation.

Figures

Figures reproduced from arXiv: 2412.01515 by the authors.

Figure 2
Figure 2. The ADC spectrum recorded in the Chacaltaya Laboratory in 82 days. The peak around channel 475 is an artifact due to the amplifier saturation. The energy calibration of the LABDOS01 spectrometer was performed by NPI CAS using two reference radionuclide sources of α particle, 239Pu and 241Am, and He ions accelerated to an energy of 572.836 MeV at the HIMAC accelerator (NIRS Japan (Yamada, 1995)). The SRIM program (Zi… view at source ↗
Figure 4
Figure 4. Dose rates as a function of rigid￾ity cutoff at the typical cruise altitude h=11000 m. As the altitude changes from h = 2000 m to 12000 m, both dose rates DSi and H∗ (10) increase by a factor of ∼ 30-50, depending on the rigidity cutoff. Variations in the rigidity cutoff have smaller effects than altitude on both indices, as shown in [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Altitude and rigidity cutoff during the three flights from Milan to Christchurch. The grey area marks the part without tracking information. The rigidity cutoff Rc increases from 4.8 GV in Milan to 17.2 GV at the geomagnetic equator, to decrease up to 2.7 GV at Christchurch. Tracking may be interrupted due to lack of radar coverage or a receiver on the surface. This is the case for the second flight, when the plane … view at source ↗
Figures from the paper (4 more)
Figure 6
Figure 6. Figure 6: Upper panel: absorbed dose rate in silicon measured by LABDOS01-B8 during [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 7
Figure 7. Figure 7: Altitude and rigidity cutoff for the three flight segments from Milan to La Paz: [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]
Figure 8
Figure 8. Figure 8: Upper panel: absorbed dose rate in silicon measured by LABDOS01-EF during [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]
Figure 9
Figure 9. Figure 9: Spectra of the deposited energy in silicon measured by LABDOS01 devices in [PITH_FULL_IMAGE:figures/full_fig_p015_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references · 11 canonical work pages

  1. [10]

    Dachev, T., Tomov, B., Matviichuk, Y., Dimitrov, P., Semkova, J., Jor- danova, M., Bankov, N., Ploc, O., Litvak, M., Tretyakov, V., Shurshakov, V., Benghin, V.,

    1016/j.asr.2009.08.007. Dachev, T., Tomov, B., Matviichuk, Y., Dimitrov, P., Semkova, J., Jor- danova, M., Bankov, N., Ploc, O., Litvak, M., Tretyakov, V., Shurshakov, V., Benghin, V.,

    work page 2009

  2. [11]

    Neutron induced nuclear counter effect in Hamamatsu silicon PIN diodes and APDs, in: IEEE Nuclear Science Symposium & Medical Imaging Conference, pp. 601–607. doi:10.1109/ NSSMIC.2010.5873830. Ziegler, J.,

    work page arXiv 2010

  3. [12]

    Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 268, 1818–1823

    SRIM – the stopping and range of ions in matter (2010). Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 268, 1818–1823. doi:10.1016/j.nimb.2010.02.091. 22

    doi:10.1016/j.nimb.2010.02.091 2010

  4. [1315]

    Kákona, M., Šlegl, J., Kyselová, D., Sommer, M., Kákona, J., Lužová, M., Štěpán, V., Ploc, O., Kodaira, S., Chroust, J.,

    doi:10.1140/epjp/s13360-022-03439-3. Kákona, M., Šlegl, J., Kyselová, D., Sommer, M., Kákona, J., Lužová, M., Štěpán, V., Ploc, O., Kodaira, S., Chroust, J.,

    work page doi:10.1140/epjp/s13360-022-03439-3

  5. [2005]

    Radiation Measurements 39, 219–223

    Response on a Si-diode-based device to fast neutrons. Radiation Measurements 39, 219–223. doi:10.1016/j.radmeas.2004.05

    work page doi:10.1016/j.radmeas.2004.05 2004

  6. [2010]

    Use of energy deposi- tion spectrometric Liulin for individual monitor of aircrew. Radiat. Prot. Dosimetry 144, 1–4. doi:10.1093/rpd/ncq505. Sommer, M., Štěpánová, D., Kákona, M., Velychko, O., Ambrožová, I., Ploc, O.,

    work page doi:10.1093/rpd/ncq505

  7. [2017]

    Radiation Protection Dosimetry 175, 419–431

    CARI-7A: development and validations. Radiation Protection Dosimetry 175, 419–431. doi:10.1093/rpd/ncw369. Dachev, T.,

    work page doi:10.1093/rpd/ncw369

  8. [2018]

    Space Weather 16, 1269–1276

    First steps toward the verification of mod- els for the assessment of the radiation exposure at aviation altitudes during quiet space weather conditions. Space Weather 16, 1269–1276. doi:10.1029/2018SW001984. Ploc, O., Brabcová, K., Spurný, F., Dachev, T.,

    work page doi:10.1029/2018sw001984

Show all 12 references
  1. [2020]

    (Eds.), 20 Twelfth Workshop on Solar Influences on the Magnetosphere, Ionosphere and Atmosphere, pp

    Registration of neutrons with a single detector Liulin type spectrometer, in: Georgieva, K., Kirov, B., Danov, D. (Eds.), 20 Twelfth Workshop on Solar Influences on the Magnetosphere, Ionosphere and Atmosphere, pp. 109–116. doi:10.31401/WS.2020.proc. Fontanilla, A., Calamida, ...

    doi

  2. [2021]

    Journal of instrumentation 16, T03006

    AIRDOS — open-source PIN diode airborne dosimeter. Journal of instrumentation 16, T03006. doi:10.1088/1748-0221/16/03/T03006. Matthias M. Meier, M.M., Copeland, K., Matthiä, D.and Mertens, C., Schennetten, K.,

    doi

  3. [2022]

    Radiat Prot Dosimetry 198, 597–603

    Calibration of silicon detectors Liulin and AIRDOS using cosmic rays and timepix for use at flight altitudes. Radiat Prot Dosimetry 198, 597–603. doi:10.1093/rpd/ncac104. Spurny, F.,

    doi

  4. [2023]

    Yamada, S.,

    doi:10.22323/1.444.1254. Yamada, S.,

    doi

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.