{"id":"dce7e66f-2055-43f2-ab39-6cb9b610a0b8","arxiv_id":"2412.01515","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A low-cost silicon diode spectrometer is validated against the CARI-7A model in flight, and a two-parameter empirical conversion to ambient dose equivalent H*(10) is shown to hold within about 15 percent in the tested range.","lead":"LABDOS01, a low-cost silicon diode spectrometer, measured cosmic-ray dose rates on two intercontinental flights and matched the CARI-7A model within about 10 percent. The paper also fits a two-coefficient formula that converts the measured deposited-energy spectrum into the ambient dose equivalent H*(10), checked on a separate flight and against ground neutron counters.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The two conversion ratios in Eq. (2) are calibrated only from cruise-altitude flights but are applied to ground-level spectra; the ground check validates only the neutron term, so the claimed better-than-15% accuracy for H*(10) at the deployed sites is not established.","rationale":"The paper's strongest empirical support is the DSi comparison: measured absorbed dose rates track CARI-7A within roughly 10–16% across eight flight segments (Tables 2 and 3), and the instrument design, calibration, and dead-time corrections are clearly described. The load-bearing risk is not the instrument but the generalization of the two-coefficient formula. The four calibration segments span a wide rigidity-cutoff range but are all at cruise altitude; Figure 3 shows dose rate changes by a factor of 30–50 between 2 km and 12 km, so extrapolating flight-derived constants to ground level is a major extrapolation that the current data cannot bound. The ground check only validates the neutron-dominated high-LET term, and at Chacaltaya the H_H estimate (107 nSv/h) is already 12% below the Wendi-2 value (121.3 nSv/h), close to the claimed 15% accuracy limit. The one out-of-sample flight (Milan–La Paz) contains a segment with a 27% relative difference, and the substituted tracking data for MAD-VVI introduces an additional uncontrolled systematic. These limitations are consistent with the reader's conditional verdict: the core method is plausible and the DSi measurements are solid, but the accuracy claim is broader than the evidence supports. A decisive test—computing the conversion ratios over a wider altitude/rigidity grid or obtaining one independent total-H* ground measurement—would settle whether the constants are universal. Since the paper already warrants conditional acceptance pending such a check, the verdict remains unchanged rather than moving to reject.","tokens_in":11270,"tokens_out":8111,"duration_ms":71788,"concrete_test":"Run CARI-7A (or a Monte Carlo code such as MCNPX/FLUKA) over a grid of altitudes (0, 1, 3, 5, 8, 12 km) and rigidity cutoffs (0, 5, 10, 15 GV) that brackets the deployed environments, and compute the separate low-LET and high-LET ratios H*_L/DL and H*_H/DH using the same 1.5 MeV silicon-deposit partition as in the paper. If any grid point—particularly at ground level (0 km)—gives a ratio that deviates from 1.19 or 10.3 by more than 15%, the constant coefficients in Eq. (2) do not hold over the claimed domain.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (2) — H+ = 1.19 DL + 10.3 DH — assumes the ratios H*_L/DL and H*_H/DH are constant over the intended application domain. These coefficients are derived from only four flight segments at cruise altitudes (10–12.5 km, Rc 2.7–17.2 GV; Tables 4 and 5) and then applied to ground-level sites (Turin, 240 m; Concordia, 3233 m; Chacaltaya, 5240 m; Table 7). At ground level the radiation field differs materially: the low-LET component is enhanced by local radioactivity, muons, and a softer neutron spectrum, and the 1.5 MeV partition in silicon may not track the changes in the H*(10)-to-silicon-dose conversion. The ground validation in Table 7 only checks the high-LET term (RH·DH against Wendi-2 neutron H*), not the low-LET term (RL·DL); moreover, the Turin Wendi-2 data are from three years earlier. The out-of-sample flight test (Table 6) is limited to one route, one segment of which uses tracking data from a different day, and one segment (VVI-LPZ) shows a 27% relative difference (0.70 vs 0.55 µSv), albeit within statistical errors. The paper itself states that ground-level validation of the formula would require ad hoc simulations 'beyond the scope of this work,' so the claimed 'accuracy better than 15%' is not supported at the deployed sites.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11631,"tokens_out":6600,"duration_ms":57616,"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":[{"comment":"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.","section":"Sec. 4.1, Eq. (2), Table 7"},{"comment":"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.","section":"Sec. 2, Tables 2–7"},{"comment":"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.","section":"Sec. 4.1, Table 6"}],"minor_comments":[{"comment":"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'.","section":"Sec. 4.1, last paragraph"},{"comment":"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).","section":"Sec. 3.2, Fig. 7 and Table 3"},{"comment":"There is a typo: 'the capability of the device to to disentangle' should read 'the capability of the device to disentangle'.","section":"Sec. 5, first paragraph"},{"comment":"The instrument is variably called 'LABDOS01' and 'LABDOS1' (e.g., in Section 5); please use one consistent name.","section":"Throughout"},{"comment":"The caption states 'h=11000 km'; this should be 'h=11000 m'.","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid detector-characterization study with a useful application, but the central accuracy claim for H*(10) is broader than the evidence. The calibration of Eq. (2) is entirely against CARI-7A on one route, so the 'better than 15%' statement is essentially a statement about how well the formula reproduces that model in a limited domain. The only absolute check is the neutron comparison at three ground sites, which covers only the high-LET term. I recommend the editor ask for a revision that either limits the claim to the cruise-altitude regime (with a clear statement that ground H+ is an extrapolation), or provides additional validation of the low-LET term, e.g., by comparing total H+ with a directly calibrated reference instrument or by performing the ad hoc simulations mentioned in the text. The paper's overall contribution is valuable and the revision is within scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent, small-scope instrument paper. The LABDOS01 reproduces CARI-7A absorbed dose estimates within ~10% on most flight segments, and the two-coefficient H*(10) formula passes an out-of-sample flight test and a ground neutron check. The soft spot is that the conversion coefficients are calibrated at cruise altitudes and applied at ground level, where only the high-LET (neutron) term is tested; the low-LET term is untested on the ground. The 'better than 15%' claim is somewhat stronger than the tables support.\n\nWhat's new: the specific LABDOS01 detector, its open-hardware design, energy calibration with HIMAC and alpha sources, and the flight validation over a wide rigidity range (2.7-17.2 GV). The empirical H*(10) formula is an adaptation of Ploc et al. (2010), but the coefficients (1.19 and 10.3) are new for this detector, and the ground data from Concordia and Chacaltaya are valuable.\n\nWhat it does well: the DSi comparison is careful, with dead-time corrections, tracking interpolation, and a sensible discussion of the silicon thickness effect. The out-of-sample Milan-La Paz flight and the Wendi-2 comparison provide independent checks of the method.\n\nWhere it's soft: the two coefficients are fitted to CARI-7A on the same route used for the DSi validation, so the calibration and reference model are entangled, though the out-of-sample checks mitigate this. More importantly, the formula is applied to ground-level spectra without validating the low-LET term; the ground check only compares H_H+ to Wendi-2 neutron doses. The VVI-LPZ segment shows a 27% difference in H+ (0.70 vs 0.55 µSv), technically within statistical error but not reassuring. The 8% systematic uncertainty in the diode mass is not propagated into any of the quoted ratios. The 'better than 15%' claim in the conclusions is not fully supported by the tables, especially for ground sites. The substituted tracking data for the Madrid-Santa Cruz flight is a necessary but acknowledged limitation.\n\nWho it's for: applied dosimetry groups, especially those using Liulin-type detectors for environmental monitoring. It doesn't break new ground but is a solid characterization useful for the SAMADHA and CORDIAL projects.\n\nRecommendation: send to peer review. The core finding (the detector works, and the H* formula is plausible within its calibration domain) is sound. The authors should be asked to soften the accuracy claim, discuss the domain of applicability, and ideally provide some estimate of the low-LET ground validation or at least acknowledge the limitation more prominently.","headline":"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.","tokens_in":12198,"tokens_out":3368,"would_cite":true,"duration_ms":27025,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["29.40.Wk"],"model":"deepseek-v4-flash","headline":"A silicon diode spectrometer can gauge biological dose to within 15%.","keywords":["dosimetry","cosmic rays","silicon diode spectrometer","Liulin-type detector","ambient dose equivalent","LET spectroscopy","CARI-7A","mixed radiation fields"],"falsifier":"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.","tokens_in":11052,"feed_emoji":"☢️","tokens_out":7758,"duration_ms":63800,"temperature":0.7,"pith_summary":"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.","feed_headline":"Silicon diode turns cosmic-ray dose into biological dose to 15%","feed_subtitle":"A cheap, open-source silicon chip now tracks the dose living tissue absorbs, verified on flights and against neutron counters.","key_machinery":"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)$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Proposes the energy-threshold LET-splitting method that this paper adapts to estimate $H^{*}(10)$ from a Liulin-type silicon spectrum.","marker":"Ploc et al., 2010"},{"why":"Documents the CARI-7A code that supplies the reference absorbed-dose and $H^{*}(10)$ values used for calibration and comparison.","marker":"Copeland, 2017"},{"why":"Independently verifies the model chain behind CARI-7A, supporting its use as a reference benchmark in this work.","marker":"Matthias M. Meier et al., 2018"},{"why":"Characterizes the Liulin-type spectrometer class that LABDOS01 belongs to, establishing its suitability for mixed radiation fields.","marker":"Dachev, 2009"},{"why":"Supports the paper's assumption that the high-LET contribution to $H^{*}(10)$ is almost exclusively due to neutrons.","marker":"Dachev et al., 2020"},{"why":"Shows that silicon diode detectors respond to fast neutrons above roughly 1 MeV, justifying the high-LET interpretation of the spectrum.","marker":"Spurny, 2005"},{"why":"Demonstrates neutron-induced nuclear effects in Hamamatsu silicon PIN diodes, which relates to the high-energy deposits used to derive $D_H$.","marker":"Zhang et al., 2010"},{"why":"Describes the open-source AIRDOS predecessor whose design LABDOS01 follows, placing the instrument in an established hardware lineage.","marker":"Kákona et al., 2021"}],"fun_headline_variants":["Silicon diode turns cosmic-ray dose into tissue dose within 15%","Cosmic-ray dose to tissue: silicon diode accuracy 15%","Silicon chip maps cosmic rays to biological dose within 15%","Silicon sensor converts cosmic-ray dose to biological dose, 15% accurate","Diode dosimeter fits cosmic-ray dose to tissue within 15%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Silicon diode turns cosmic-ray dose into tissue dose within 15%","Cosmic-ray dose to tissue: silicon diode accuracy 15%","Silicon chip maps cosmic rays to biological dose within 15%","Silicon sensor converts cosmic-ray dose to biological dose, 15% accurate","Diode dosimeter fits cosmic-ray dose to tissue within 15%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002298,"raw_usage":{"total_tokens":8877,"prompt_tokens":966,"completion_tokens":7911,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":7816}},"tokens_in":582,"tokens_out":7911,"duration_ms":43556,"temperature":1.0,"reasoning_tokens":7816,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:17:20.138994+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}