Establishing an independent measurement traceability for 60-Co Air Kerma

Cristina Garc\'ia-Mulas; N\'estor Cornejo-D\'iaz; Paz Avi\'es-Lucas

REVIEW 3 major objections 4 minor 16 references

Establishing an independent measurement traceability for 60-Co Air Kerma

T0 review · 3 major / 4 minor · reviewed 2026-07-31 · deepseek-v4-flash

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

Pith's one-line read A single calculated correction factor, k_Co,Cs = 0.988 ± 0.008, extends cesium-137 air-kerma traceability to cobalt-60 radiation-protection beams.

desk verdict A credible Monte-Carlo-based k_Co,Cs for a 1000 cm3 transfer chamber, with an internally consistent uncertainty budget, but the paper's headline external validation is asserted rather than shown. read the letter →

arxiv 2607.27875 v1 pith:VGYHDVO3 submitted 2026-07-30 physics.ins-det physics.comp-phphysics.med-ph

Néstor Cornejo-Díaz , Cristina García-Mulas , Paz Aviés-Lucas This is my paper · ORCID

classification physics.ins-detphysics.comp-phphysics.med-ph

keywords 60Coairkermabeam-qualitycorrectionfactorMonteCarlosimulationionizationchambermeasurementtraceabilityradiationprotectiondosimetryprimarystandardk_Q

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

The paper establishes an independent traceability chain for 60Co air-kerma measurements at radiation-protection levels without constructing a dedicated 60Co primary standard. It derives a beam-quality correction factor k_Co,Cs by Monte Carlo simulation for a large-volume ionization chamber with a flat energy response, then uses that factor to convert a cesium-137 primary-standard calibration into a cobalt-60 calibration. The computed factor is 0.988 with an expanded uncertainty of 0.8%, and the authors show it agrees with four historical chamber calibrations and with an international supplementary comparison. If correct, any laboratory with a well-characterized chamber and a cesium-137 primary standard can supply cobalt-60 protection-level traceability at comparable uncertainty, without depending on external primary calibrations.

What carries the argument

The beam-quality correction factor k_Co,Cs, defined as the ratio of the air-kerma calibration coefficients for the two radiation qualities, is the central object. The calculation cancels the chamber air mass and W/e, so the exact cavity volume is not required. Monte Carlo simulation supplies the air kerma per unit fluence and the absorbed dose in the chamber cavity for each quality, and sensitivity studies add uncertainty components from photon spectra, interaction cross sections, axial beam non-uniformity, and the chamber's geometric model.

What would settle it

Directly calibrate the same chamber in both beam qualities against two independent primary standards for cesium-137 and cobalt-60; if the measured ratio N_K,Co / N_K,Cs differs from 0.988 by more than the combined expanded uncertainties of the two calibrations and the simulation, the Monte Carlo model is biased. A simpler falsifier is to measure the chamber's response across a range of monoenergetic photon energies and compare that curve with the simulated energy response.

Watch

Extended reading notes

Core claim

The central claim is that the ratio of a chamber's air-kerma calibration coefficients between cobalt-60 and cesium-137 can be computed from first principles via Monte Carlo simulation, rather than measured against a primary standard for each beam. Because the air mass in the cavity and the mean energy per ion pair cancel in the ratio, the final factor does not depend on the cavity volume or on W/e. The paper obtains k_Co,Cs = 0.988, with a combined standard uncertainty of 0.41% and an expanded uncertainty of 0.8%, and validates it against historical calibration data collected over several years and against an international comparison in which the laboratory took part as a primary laboratory.

Load-bearing premise

The Monte Carlo model of the ionization chamber accurately represents how the real chamber's response changes between cesium-137 and cobalt-60, even though the model is built from manufacturer drawings held under a confidentiality agreement, supplemented only by X-ray images and conservative sensitivity variations.

Editorial extensions

If this is right

A cesium-137 primary standard can now provide cobalt-60 protection-level traceability with an expanded uncertainty near 0.8%.
Other laboratories with similarly flat-response, large-volume chambers and their own cesium-137 standards could adopt the same ratio-based method.
The uncertainty budget identifies chamber modelling as the dominant term, showing where effort would most reduce the final uncertainty.
The cancellation of cavity volume and W/e relaxes the need for exact dimensional knowledge of the transfer chamber.
The method removes the requirement for external primary calibration of cobalt-60 beams, simplifying long-term stability monitoring and comparison participation.

Reading between the lines

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

Because chamber model uncertainty dominates the budget, a direct measurement of the same chamber's response ratio against two independent primary standards would be a much stronger validation than the historical ratios, whose uncertainties are roughly twice as large.
The same ratio-symmetry approach could be extended to other beam qualities (for example, 241Am or lower-energy X-rays), but only where the chamber's energy response is flat enough that the Monte Carlo model can be trusted to the needed accuracy.
A testable extension is to compare the simulated energy-response curve of the chamber with measurements in a set of monoenergetic photon beams; this would directly expose any systematic bias in the geometric or material model.
If the chamber model is later revised from direct dimensional measurements, the k_Co,Cs value could shift by more than the statistical component, which suggests the 0.41% combined uncertainty may be sensitive to modelling choices.

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.

Desk Editor's Note

read the letter

Colleague,

The useful thing in this paper is a new number: k_Co,Cs = 0.988 ± 0.008 (k=2) for a specific PTW 32002 chamber, computed with EGSnrc, with a decomposed uncertainty budget. That is genuinely new. The trick of letting air mass and W/e cancel (eqs. 4-6) is clean, and the Monte Carlo N_K values line up with the historical certificates within 0.1% for 137Cs and 0.3% for 60Co. The arithmetic checks out: table 5 ratio is 0.9882, table 6 quadrature is 0.41%. For a metrology lab that wants to wean itself off external 60Co calibrations using a 137Cs primary, this is a workable template.

The soft spots are real but not fatal. The biggest is the S19 comparison: the paper says CIEMAT participated as a primary laboratory and showed 'excellent agreement' with NPL, but the degrees of equivalence are not reported. That is the single strongest external anchor, and the reader cannot check it. The authors should print the DoE values and the comparison reference value, or at least a table of results. Similarly, the analytic-vs-MC air kerma validation in §2.2.1 is asserted but not tabulated; a one-line table would close that gap. The chamber-model uncertainty (0.30%, the largest component) is judged from a sensitivity study that simplifies the central electrode support and varies wall mass thickness by ±10%. That is a reasonable first cut, but it is not a direct measurement of energy response, and the geometry is NDA-protected, so an independent check is not possible from the text. A correlated geometric bias at the 0.2-0.3% level cannot be fully excluded. The historical ratios (0.984-0.987) bracket the MC value with ±1.3-1.5% expanded uncertainties, which is reassuring but too coarse to resolve a 0.3% bias.

None of this amounts to a demonstrated error. If the S19 degrees of equivalence are indeed good, the central claim stands. As written, the paper's conclusion outruns the evidence shown in the text — 'conclusively validates' is too strong for an unquantified comparison. That is a revision, not a rejection.

This is a specialist paper for people in ionizing-radiation metrology, not a general audience. A serious referee should see it, mainly to pressure the authors into reporting the S19 numbers and the analytic comparison. I would send it to review, with the expectation that the revision closes the reporting gaps.

Referee Report

3 major / 4 minor

Summary. The paper reports a method for establishing an independent air-kerma traceability chain for 60Co radiation-protection level beams at CIEMAT. Because the existing 137Cs primary standard has insufficient signal in low-dose-rate 60Co beams, the authors calibrate a 1000 cm3 PTW 32002 secondary-standard chamber against the 137Cs primary standard and multiply the resulting N_K,Cs by a Monte Carlo-derived beam-quality correction factor k_Co,Cs. The factor is obtained from EGSnrc simulations of air kerma and cavity absorbed dose, exploiting the cancellation of the air mass and W/e in the ratio. The result is k_Co,Cs = 0.988 ± 0.008 (k = 2) with a combined standard uncertainty of 0.41%, supported by an uncertainty budget (Table 6), historical calibration certificates (Table 7), and an asserted validation through EURAMET.RI(I)-S19.

Significance. If the result is correct, the paper offers a practical independent route to 60Co protection-level traceability without a dedicated 60Co primary standard, with potential applicability to other secondary-standard chambers. The cancellation of the air mass and W/e is a clean and clearly explained formal step, and the MC calculation contains no free parameters fitted to the target quantity. The arithmetic is internally consistent: Table 5 gives N_K,Co/N_K,Cs = 2.4700/2.4996 = 0.9882, and Table 6 combines to 0.41%. The absolute MC-derived values agree with historical PSDL-linked certificates within about 0.3%, which is creditable. The principal weaknesses are that the largest uncertainty component (chamber model, 0.30%) is justified by a qualitative sensitivity study rather than a direct measurement of chamber energy response, and that the strongest external validation (EURAMET.RI(I)-S19) is reported only descriptively, without degrees of equivalence or numerical comparison data. These gaps leave the central claim not fully verifiable from the manuscript as submitted.

major comments (3)

[§3.2.2] The claimed international validation is not quantified. The text states that CIEMAT's S19 results 'demonstrate excellent agreement' and 'conclusively validates the accuracy of the calculated k_Co,Cs factor', but it gives no degrees of equivalence, comparison reference value, or associated uncertainties. A supplementary comparison normally reports DoE values for each participating laboratory; please include the actual numbers, either in a table or in the text, with their uncertainties and the comparison reference. Without them, the reader cannot assess the strength of the external validation or check its consistency with k_Co,Cs = 0.988 ± 0.008.
[§3.1 / Table 6 / §2.2.2] The dominant uncertainty component, the PTW 32002 chamber model (0.30%), is not documented sufficiently. The text says only that the evaluation used 'a simplified model—excluding some details of the central electrode support—and applying maximum variations (±10%) to the outer wall mass thickness'. This does not specify how many model variants were simulated, which parameters were varied, how the observed spread was converted to a standard uncertainty, or whether correlated variations (graphite coating thickness, internal gaps, electrode-support geometry) were considered. Because this component controls the 0.41% combined standard uncertainty, the derivation of the 0.30% value must be reproducible from the text, or the stated expanded uncertainty is not verifiable.
[§3.2.1] The historical validation is bracketing but not resolving. The expanded uncertainties on the historical k_Co,Cs values in Table 7 are 1.3–1.5%, roughly three to four times the MC combined standard uncertainty of 0.41%. The agreement of all four historical ratios with 0.988 is reassuring, but at this precision a 0.3–0.4% systematic bias in the chamber-model component would not be detectable. The paper should state this limitation explicitly and, if possible, provide an additional experimental cross-check with a smaller uncertainty.

minor comments (4)

[§2.2] Typographical error: 'the subscript Q my be omitted' should read 'may be omitted'.
[Table 3] The column header 'Energía / keV' is in Spanish; the rest of the paper is in English. Please change to 'Energy / keV'.
[§3.2.2 / Ref. [16]] There is an inconsistency in the comparison name: 'EURAMET.RI(I)-S19' appears in the text, while the reference list and URL use 'EURAMET.RI(I).S-19' and 's19'. Please standardize and, if possible, provide the full CIPM MRA comparison report identifier and date.
[§1 / Table 1] The notation 'S-Cs' and 'S-Co' is used in Table 1 and throughout without explicit definition in the text. Define these beam-quality labels at first use, or state that they follow ISO 4037-3.

Circularity Check

0 steps flagged · score 1.0 of 10

No material circularity: k_Co,Cs is simulated from first principles; self-citations are external comparisons, not load-bearing.

full rationale

The derivation of k_Co,Cs uses equations (2)-(6): the Monte Carlo ratio of K_air to chamber cavity dose. No parameter is fitted to the final k value; chamber air mass and W/e cancel by construction (eqs. 5-6), and the reported value 0.9882(3) follows from simulated K_air and D values in Table 5. The largest uncertainty component, the PTW 32002 chamber model (0.30%, Table 6), is bounded by a conservative sensitivity study (simplified support geometry and ±10% wall thickness), not by adjustment to the result. External validation does not rely on the same calculation: historical k ratios (Table 7) come from PSDL-calibrated PTW 32005 transfer calibrations, and EURAMET.RI(I)-S19 [16] benchmarks CIEMAT against the NPL primary standard; BIPM.RI(I)-K5 [2] is a BIPM-led comparison. Although refs [1], [2], and [16] include the first author, this is normal authorship overlap and does not reduce the argument to self-citation. The S19 agreement is reported as 'excellent agreement' without numerical degrees of equivalence, which is an unquantified validation claim and a verifiability gap, but it is not a circular reduction. No equation or fitted parameter forces the target result from its input.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The central result has no fitted free parameters: k_Co,Cs is the simulation output; W/e and chamber air mass cancel in the ratio; ECUT choices (1.173/1.844 MeV) are physically determined by the maximum electron energies; photon-splitting and ESTEPE are numerical settings. The load-bearing inputs are domain assumptions about the Monte Carlo physics, spectral models, chamber-model fidelity, and the correctness of the underlying 137Cs standard. No new physical entities are introduced.

assumptions (6)

domain assumption (W/e)_Q is independent of photon/electron energy across 0.66-1.33 MeV, permitting its cancellation in the k ratio
Invoked in §2.2 (eq. 4 through eqs. 5-6); standard ICRU-90 practice, but a residual ~0.1-0.2% energy dependence would partly propagate into the individual NK values (though largely cancelling in the ratio).
domain assumption EGSnrc cross-section libraries (xcom, MCDF-XCOM, MCDF-EPDL) describe photon interactions in air and the chamber materials (POM, graphite) at these energies
Used throughout §2.2.1-2.2.2; uncertainty bounded at 0.03% via library cross-checks; a tooling assumption, not independently verified in this paper.
domain assumption Fano-condition photon regeneration with ECUT set to the maximum electron energy yields true air kerma
§2.2.1; the analytical cross-check (eq. 7) would support it, but the comparison is asserted, not tabulated.
domain assumption The NDA-based chamber model, supplemented by X-ray imaging, captures the true cavity-dose ratio within the 0.30% sensitivity bound
§2.2.2 and table 6; the largest uncertainty term, bounded by a ±10% wall-mass variation and a simplified support model rather than by a direct response measurement.
domain assumption The three photon spectra used (pure emission lines, laboratory irradiator spectra, and the Mora et al. therapy spectrum) span the true S-Cs and S-Co beams
§2.2.2, refs [11-14]; 0.20% uncertainty component; a therapy-source spectrum is not obviously representative of a filtered protection-level beam, though the 0.20% allowance is meant to cover this.
domain assumption The 137Cs primary standard (ref [1]) and its BIPM.RI(I)-K5 equivalence (ref [2]) are correct
The whole chain is anchored to this standard; the externally administered key comparison gives independent support, but the standard itself is a self-cited prior result (ref [1], first author overlaps with the present first author).

how reviews work

0 comments

Cite this review

Pith. "Pith review of Establishing an independent measurement traceability for 60-Co Air Kerma." pith.science (2026). https://pith.science/paper/VGYHDVO3

@misc{pith2026260727875,
  author       = {Pith},
  title        = {Pith review of: Establishing an independent measurement traceability for 60-Co Air Kerma},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VGYHDVO3}},
  note         = {Machine review of arXiv:2607.27875}
}

abstract

An independent air kerma traceability chain for $^{60}$Co radiation protection levels has been successfully established at the CIEMAT Ionizing Radiation Metrology Laboratory (LMRI-CIEMAT), based on the reference provided by the $^{137}$Cs primary standard. To achieve this, a secondary standard ionization chamber with an appropriate energy response was characterized, and its beam-quality correction factor, $k_Q$, was accurately determined via Monte Carlo simulations using the EGSnrc code. The method's accuracy was validated through a comparison with long-term historical calibration data and peer-confirmed through CIEMAT's successful participation in the EURAMET.RI(I)-S19 supplementary comparison, recently officially published in the BIPM KCDB. The excellent agreement achieved with international reference values demonstrates the high robustness and traceability independence of this newly implemented methodology.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references

[1]

Cornejo Díaz, New LMRI-CIEMAT primary standard for 137Cs air-kerma

N. Cornejo Díaz, New LMRI-CIEMAT primary standard for 137Cs air-kerma. Measurement 188 (2022) 110374

2022
[2]

Kessler, P

C. Kessler, P. Roger, N. Cornejo Díaz, Key comparison BIPM.RI(I) -K5 of the air -kerma standards of the CIEMAT, Spain and the BIPM in 137Cs gamma radiation. Metrologia 61 - 1A (2024) 06007

2024
[3]

ISO 4037-3 (2019)

International Standard Organization, Radiological protection – X and gamma reference radiation for calibrating dosemeters and dose rate meters and for determining their response as a function of photon energy – Part 3: Calibration of area and personal dosemeters and the measurement of their response as a function of energy and angle of incidence. ISO 4037...

2019
[4]

Cornejo Díaz, Maximum likelihood estimation using expectation maximization applied to ambient dose equivalent measurements

N. Cornejo Díaz, Maximum likelihood estimation using expectation maximization applied to ambient dose equivalent measurements. Radiation Protection Dosimetry (2018), Vol. 182, No. 2, pp. 285–293

2018
[5]

PTW Freiburg GmbH D466.131.00/04 EN (2009)

Physikalisch-Technische Werkstätten, User manual 1L Spherical chamber Type 32002 . PTW Freiburg GmbH D466.131.00/04 EN (2009)

2009
[6]

TECHNICAL REPORTS SERIES No

International Atomic Energy Agency, Absorbed Dose Determination in External Beam Radiotherapy . TECHNICAL REPORTS SERIES No . 398 Rev. 1 . VIENNA (2024). https://www- pub.iaea.org/MTCD/Publications/PDF/p15048-DOC-010-398-Rev1_web.pdf

2024
[7]

Kawrakow, E

I. Kawrakow, E. Mainegra -Hing, D. W. O. Rogers, F. Tessier, B. Walters, The EGSnrc Code System: Monte Carlo simulation of electron and photon transport . Technical Report PIRS -701, N ational Research Council Canada (2025). https://nrc-cnrc.github.io/EGSnrc/doc/pirs701-egsnrc.pdf

2025
[8]

ICRU Report No

International Commission on Radiation Units and Measurements, Key data for ionizing – radiation dosimetry: Measurement standards and applications. ICRU Report No. 90. ICRU Vol. 14 No.1 (2016)

2016

Show all 16 references

[9]

Freiburg GmbH

PHYSIKALISCH-TECHNISCHE WERKSTÄTTEN, Non-Disclosure Agreement between PTW and CIEMAT. Freiburg GmbH. LS-01 T32002. 2023-06-28 TS/Ir (2023)

2023
[10]

Gijón Tiradas, Radiografías a escala real de la cámara PTW 32002 N/S: 00345

D. Gijón Tiradas, Radiografías a escala real de la cámara PTW 32002 N/S: 00345. División de Materiales de Interés Energético. Unidad de Tecnología. CIEMAT. Personal communication (2024)

2024
[11]

Nuclear Data – Table –LNE-LNHB/CEA (lnhb.fr)

Laboratoire National Henri Becquerel, Table de radionucléides. Nuclear Data – Table –LNE-LNHB/CEA (lnhb.fr). http://www.lnhb.fr/accueil/donnees-nucleaires/donnees-nucleaires-tableau/. Last accessed January 2026 (2026)

2026
[12]

Cornejo Díaz, Uncertainties of the calibration of dosemeters above ground and deep underground compared

N. Cornejo Díaz, Uncertainties of the calibration of dosemeters above ground and deep underground compared. Sources of uncertainty during calibrations at the CIEMAT metrology laboratory. Deliverable D 3.1.7. Metrology for radiological early warning networks in Europe. MetroERM...

2016
[13]

Cornejo Díaz, E

N. Cornejo Díaz, E. Guantes Die z, Caracterización del haz de fotones del irradiador NI -645 mediante simulación estadística. LMRI-CIEMAT. DT-LMRI-1804 (2018)

2018
[14]

G. M. Mora, A. Maio, D. W. O. Rogers, Monte Carlo simulation of a typical 60Co therapy source. Med. Phys. 26 (11), 2494 – 2502 (1999)

1999
[15]

Kawrakow, E

I. Kawrakow, E. Mainegra-Hing, F. Tessier, R. Townson, B. Walters B., Geometry module. EGSnrc C++ class library. ReportPIRS-898 (2021). https://nrc-cnrc.github.io/EGSnrc/doc/pirs898/modules.html. – 10 –

2021
[16]

Persson, E.L

L. Persson, E.L. Hansen, P.O. Hetland, N. Cornejo Díaz, M. Kelly, J. Lillhök, R. Nylund, J. Huikari, Air kerma in 241Am an d 60Co radiation protection beams . CIPM MRA Comparison R eports (2026) 06003. EURAMET.RI(I).S-19. https://www.bipm.org/en/d/-euramet-ri-i-s19

2026

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.

[1] [1]

Cornejo Díaz, New LMRI-CIEMAT primary standard for 137Cs air-kerma

N. Cornejo Díaz, New LMRI-CIEMAT primary standard for 137Cs air-kerma. Measurement 188 (2022) 110374

2022

[2] [2]

Kessler, P

C. Kessler, P. Roger, N. Cornejo Díaz, Key comparison BIPM.RI(I) -K5 of the air -kerma standards of the CIEMAT, Spain and the BIPM in 137Cs gamma radiation. Metrologia 61 - 1A (2024) 06007

2024

[3] [3]

ISO 4037-3 (2019)

International Standard Organization, Radiological protection – X and gamma reference radiation for calibrating dosemeters and dose rate meters and for determining their response as a function of photon energy – Part 3: Calibration of area and personal dosemeters and the measurement of their response as a function of energy and angle of incidence. ISO 4037...

2019

[4] [4]

Cornejo Díaz, Maximum likelihood estimation using expectation maximization applied to ambient dose equivalent measurements

N. Cornejo Díaz, Maximum likelihood estimation using expectation maximization applied to ambient dose equivalent measurements. Radiation Protection Dosimetry (2018), Vol. 182, No. 2, pp. 285–293

2018

[5] [5]

PTW Freiburg GmbH D466.131.00/04 EN (2009)

Physikalisch-Technische Werkstätten, User manual 1L Spherical chamber Type 32002 . PTW Freiburg GmbH D466.131.00/04 EN (2009)

2009

[6] [6]

TECHNICAL REPORTS SERIES No

International Atomic Energy Agency, Absorbed Dose Determination in External Beam Radiotherapy . TECHNICAL REPORTS SERIES No . 398 Rev. 1 . VIENNA (2024). https://www- pub.iaea.org/MTCD/Publications/PDF/p15048-DOC-010-398-Rev1_web.pdf

2024

[7] [7]

Kawrakow, E

I. Kawrakow, E. Mainegra -Hing, D. W. O. Rogers, F. Tessier, B. Walters, The EGSnrc Code System: Monte Carlo simulation of electron and photon transport . Technical Report PIRS -701, N ational Research Council Canada (2025). https://nrc-cnrc.github.io/EGSnrc/doc/pirs701-egsnrc.pdf

2025

[8] [8]

ICRU Report No

International Commission on Radiation Units and Measurements, Key data for ionizing – radiation dosimetry: Measurement standards and applications. ICRU Report No. 90. ICRU Vol. 14 No.1 (2016)

2016