pith. sign in

arxiv: 2606.05807 · v1 · pith:MANBQGU7new · submitted 2026-06-04 · 🌌 astro-ph.IM · astro-ph.HE

Simulation of proton-induced activation for low-Earth orbit high energy astrophysics missions

Pith reviewed 2026-06-27 23:52 UTC · model grok-4.3

classification 🌌 astro-ph.IM astro-ph.HE
keywords proton-induced activationlow-Earth orbitinstrumental backgroundBateman equationsGeant4 simulationradioisotope productionspace astrophysics missionsbackground modeling
0
0 comments X

The pith

A three-step algorithm separates isotope production, decay evolution and detector response to efficiently simulate proton activation backgrounds in LEO astrophysics instruments.

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

The paper introduces an efficient three-step method for modeling the radioactive backgrounds created when trapped protons in low-Earth orbit activate spacecraft materials. Direct Monte Carlo tracking of rare production events followed by many decays is too slow, so the new approach first uses Geant4 to catalog every radioisotope made by monochromatic proton beams, then solves the Bateman equations numerically for each decay chain under arbitrary irradiation histories, and finally simulates the detector signals from those decays. Because the production step is precomputed once per energy and the decay step is analytic, the method reconstructs full time-dependent backgrounds orders of magnitude faster than brute-force simulation. Validation against direct runs shows agreement across many decades in activity and time for two mission concepts, HERMES and eXTP. The resulting framework lets designers quickly test orbit choices, shielding, and duty cycles without prohibitive computing costs.

Core claim

The central claim is that a decoupled three-step procedure—monochromatic proton irradiation libraries generated with Geant4, numerical integration of linearized Bateman equations for decay chains, and separate simulation of detector response to each isotope—accurately reproduces the activation-induced background for continuous, time-varying LEO proton spectra while greatly lowering computational cost compared with direct Monte Carlo.

What carries the argument

Three-step algorithm that decouples Geant4-based radioisotope production from monochromatic protons, numerical Bateman-equation solutions for decay evolution, and detector-response simulation.

If this is right

  • The method enables rapid testing of different orbits, radiation models, and operational duty cycles for background budgeting.
  • It applies across detector technologies, as shown for the HERMES and eXTP instruments.
  • Precomputed isotope libraries can be reused for many irradiation histories without repeating the expensive production step.
  • The framework supports optimization of shielding and material choices before flight.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same libraries could be reused for non-proton sources if analogous production tables are generated.
  • Integration with real-time space-weather proton spectra could yield predictive background forecasts during SAA passages.
  • The approach may extend to other activation channels such as secondary neutrons if production data are added.
  • Material databases could be ranked by predicted activation background using the same machinery.

Load-bearing premise

Linear superposition of monochromatic proton results plus numerical Bateman solutions will reproduce the activation and decay behavior of real continuous, time-varying LEO proton spectra without large systematic errors introduced by the decoupling.

What would settle it

A head-to-head comparison, for a realistic LEO orbit and full time history, in which the activity or count-rate time series predicted by the three-step method deviates from a direct Monte Carlo run by more than the combined statistical uncertainty over the full dynamic range.

read the original abstract

Proton-induced activation represents a major source of instrumental background for high-energy astrophysics missions in low-Earth orbit, where trapped protons, particularly during transits within the South Atlantic Anomaly region, irradiate spacecraft materials and generate radioactive isotopes. Direct Monte Carlo simulations of activation and of the ensuing decays are computationally inefficient, due to the low probability of nuclide production and the large number of decay events required for sufficient statistical accuracy. In this paper we provide a new implementation of an efficient three-step algorithm that decouples isotope production, radioactive-decay evolution, and background synthesis, enabling rapid reconstruction of activation-induced background for arbitrary irradiation histories. The method combines Geant4-based identification of all radioisotopes produced by monochromatic proton irradiations, numerical solutions of the Bateman equations for linearized decay chains, and simulation of the detector response to each isotope decay emissions. The approach greatly reduces the computational cost while maintaining accuracy, as demonstrated through validation against direct simulations, which show excellent agreement over many orders of magnitude in activity and time. This method is applied to two representative case studies: HERMES and eXTP/LAD and WFM, covering different detector technologies and orbital configurations. The presented framework enables fast exploration of design and operational scenarios (e.g., orbit selection, radiation models, or duty cycles) and is well suited for background budgeting and optimization of future high-energy space missions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The paper presents a three-step decoupled algorithm for simulating proton-induced activation backgrounds in LEO high-energy astrophysics missions. It uses Geant4 to compute radioisotope production yields from monochromatic proton irradiations, solves the Bateman equations numerically for decay-chain evolution under arbitrary irradiation histories, and simulates detector response to the resulting decays. The method is claimed to reduce computational cost substantially while preserving accuracy, with validation against direct Monte Carlo simulations showing excellent agreement over many orders of magnitude in activity and time; it is demonstrated on HERMES and eXTP (LAD/WFM) case studies.

Significance. If the accuracy of the linear superposition and Bateman integration holds across continuous, time-varying LEO proton spectra (including SAA passages), the framework would enable rapid exploration of orbital, material, and operational scenarios for background budgeting—addressing a key practical bottleneck in mission design for instruments like scintillators or silicon detectors in high-energy astrophysics.

major comments (2)
  1. [Abstract] The central claim that the three-step decoupling 'maintains accuracy' for continuous time-varying proton spectra rests on unshown quantitative validation. The abstract asserts 'excellent agreement over many orders of magnitude' but provides no error budgets, R² values, maximum relative deviations, or references to specific figures/tables comparing decoupled vs. direct simulations across mission-duration timescales.
  2. [Method (implied in abstract description)] The linear-superposition step for non-monochromatic, time-dependent fluxes (including SAA transits) is load-bearing for the method's applicability. No details are given on energy binning, time-stepping of production rates, or handling of numerical stiffness in the Bateman solver that would confirm absence of accumulating systematic bias relative to fully coupled simulations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments and for recognizing the potential utility of the three-step algorithm for activation background simulations. We address each major comment below.

read point-by-point responses
  1. Referee: [Abstract] The central claim that the three-step decoupling 'maintains accuracy' for continuous time-varying proton spectra rests on unshown quantitative validation. The abstract asserts 'excellent agreement over many orders of magnitude' but provides no error budgets, R² values, maximum relative deviations, or references to specific figures/tables comparing decoupled vs. direct simulations across mission-duration timescales.

    Authors: We agree that the abstract would benefit from explicit references to the quantitative validation results. The manuscript contains direct comparisons between the decoupled method and full Monte Carlo simulations in the validation section. In the revised manuscript we will update the abstract to reference the relevant figures and tables and to note the error metrics employed in those comparisons. revision: yes

  2. Referee: [Method (implied in abstract description)] The linear-superposition step for non-monochromatic, time-dependent fluxes (including SAA transits) is load-bearing for the method's applicability. No details are given on energy binning, time-stepping of production rates, or handling of numerical stiffness in the Bateman solver that would confirm absence of accumulating systematic bias relative to fully coupled simulations.

    Authors: We acknowledge that additional implementation details would strengthen the description of the linear-superposition approach. We will add a dedicated subsection to the Methods section that specifies the energy binning used for proton spectra, the time-stepping scheme for production rates under time-varying fluxes (including SAA passages), and the numerical methods chosen for the Bateman solver to control stiffness. We will also include a brief discussion of tests performed to verify the absence of accumulating systematic bias. revision: yes

Circularity Check

0 steps flagged

No significant circularity; method uses external Geant4 and standard Bateman solvers with independent validation

full rationale

The paper describes a three-step decoupling of isotope production (via Geant4 monochromatic runs), decay evolution (Bateman equations), and detector response simulation. Validation compares outputs to direct Monte Carlo simulations over orders of magnitude in activity and time, providing an external benchmark. No parameters are fitted to the target LEO background data, no self-citations form the load-bearing justification, and the linear superposition assumption is stated explicitly rather than derived from the method itself. The derivation chain is self-contained against external physics lists and numerical solvers.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on the standard assumption that Geant4 accurately models proton-induced isotope production cross-sections and that the Bateman equations provide an exact solution for linear decay chains; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Geant4 physics lists correctly identify all relevant radioisotopes produced by proton interactions in spacecraft materials
    Invoked in the first step of the algorithm; accuracy of the entire pipeline depends on this.
  • standard math Decay chains are linear and can be solved independently via Bateman equations without significant branching or nonlinear effects
    Core of the second step; standard but assumes no unmodeled interactions.

pith-pipeline@v0.9.1-grok · 5771 in / 1450 out tokens · 15508 ms · 2026-06-27T23:52:41.598863+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 30 canonical work pages · 5 internal anchors

  1. [1]

    (eds.) In-Orbit Background for X-Ray Detectors, p

    Campana, R.: In: Bambi, C., Santangelo, A. (eds.) In-Orbit Background for X-Ray Detectors, p. 39. Springer Nature, Singapore (2022). https://doi.org/10. 1007/978-981-16-4544-0 28-1

  2. [2]

    Nuclear Instruments and Methods in Physics Research A506(3), 250–303 (2003) https://doi.org/10.1016/S0168-9002(03)01368-8

    Agostinelli, S., Allison, J., Amako, K., Apostolakis, J., Araujo, H., Arce, P., Asai, M., Axen, D., Banerjee, S., Barrand, G., Behner, F., Bellagamba, L., Boudreau, J., Broglia, L., Brunengo, A., Burkhardt, H., Chauvie, S., Chuma, J., Chytracek, R., Cooperman, G., Cosmo, G., Degtyarenko, P., Dell’Acqua, A., Depaola, G., Dietrich, D., Enami, R., Feliciello...

  3. [3]

    MGGPOD: a Monte Carlo Suite for Modeling Instrumental Line and Continuum Backgrounds in Gamma-Ray Astronomy

    Weidenspointner, G., Harris, M.J., Sturner, S., Teegarden, B.J., Ferguson, C.: MGGPOD: a Monte Carlo Suite for Modeling Instrumental Line and Continuum Backgrounds in Gamma-Ray Astronomy. The Astrophysical Jour- nal Supplement Series156(1), 69–91 (2005) https://doi.org/10.1086/425577 arXiv:astro-ph/0408399 [astro-ph]

  4. [4]

    New Astronomy Review50(7-8), 629–632 (2006) https://doi.org/10.1016/j.newar.2006.06.049 25

    Zoglauer, A., Andritschke, R., Schopper, F.: MEGAlib The Medium Energy Gamma-ray Astronomy Library. New Astronomy Review50(7-8), 629–632 (2006) https://doi.org/10.1016/j.newar.2006.06.049 25

  5. [5]

    Modeling of proton-induced radioactivation background in hard X-ray telescopes: Geant4-based simulation and its demonstration by Hitomi's measurement in a low Earth orbit

    Odaka, H., Asai, M., Hagino, K., Koi, T., Madejski, G., Mizuno, T., Ohno, M., Saito, S., Sato, T., Wright, D.H., Enoto, T., Fukazawa, Y., Hayashi, K., Kataoka, J., Katsuta, J., Kawaharada, M., Kobayashi, S.B., Kokubun, M., Laurent, P., Lebrun, F., Limousin, O., Maier, D., Makishima, K., Mimura, T., Miyake, K., Mori, K., Murakami, H., Nakamori, T., Nakano,...

  6. [6]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Galgoczi, G., Ripa, J., Dilillo, P., Ohno, M., Campana, R., Werner, N.: A software toolkit to simulate activation background for high energy detectors onboard satel- lites. In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K. (eds.) Space Telescopes and Instrumentation 2020: Ultraviolet to Gamma Ray. Society of Photo-Optical Instrumentation Engineers (SPIE) ...

  7. [7]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Fiore, F., Burderi, L., Lavagna, M., Bertacin, R., Evangelista, Y., Campana, R., Fuschino, F., Lunghi, P., Monge, A., Negri, B., Pirrotta, S., Puccetti, S., Sanna, A., Amarilli, F., Ambrosino, F., Amelino-Camelia, G., Anitra, A., Auric- chio, N., Barbera, M., Bechini, M., Bellutti, P., Bertuccio, G., Cao, J., Ceraudo, F., Chen, T., Cinelli, M., Citossi, M...

  8. [8]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Feroci, M., Ambrosi, G., Antonelli, M., Argan, A., Babinec, V., Barbera, M., Bas- tia, P., Bayer, J., Bellutti, P., Bertucci, B., Bertuccio, G., Bonfitto, F., Bonvicini, V., Bozzo, E., Baudin, D., Bouyjou, F., Brienza, D., Cadoux, F., Campana, R., 26 Candia, R., Cao, J., Cavazzuti, E., Ceraudo, F., Chen, T., Chen, W., Coimbra, A., Colucci, A., Cong, X., C...

  9. [9]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Hernanz, M., Feroci, M., Evangelista, Y., Meuris, A., Schanne, S., Zampa, G., Tenzer, C., Bayer, J., Nowosielski, W., Michalska, M., Kalemci, E., Sungur, M., Brandt, S., Kuvvetli, I., Alvarez Franco, D., Carmona, A., G´ alvez, J.-L., Patruno, A., In’t Zand, J., Zwart, F., Santangelo, A., Bozzo, E., Zhang, S.-N., Lu, F., Xu, Y., Campana, R., Del Monte, E.,...

  10. [10]

    Journal of Jet Propulsion 28(9), 588–592 (1958) https://doi.org/10.2514/8.7396

    Van Allen, J.A., Ludwig, G.H., Ray, E.C., McIlwain, C.E.: Observation of high intensity radiation by satellites 1958 alpha and gamma. Journal of Jet Propulsion 28(9), 588–592 (1958) https://doi.org/10.2514/8.7396

  11. [11]

    Technical Report NSSDC/WDC-A-R&S 76-06, Natl

    Sawyer, D.M., Vette, J.I.: AP-8 trapped proton model environment for solar max- imum and minimum. Technical Report NSSDC/WDC-A-R&S 76-06, Natl. Space Sci. Data Cent (1976) 27

  12. [12]

    Technical Report NSSDC/WDC-A-R&S 91-24, NASA Goddard Space Flight Center (1991)

    Vette, J.I.: The AE-8 trapped electron model environment. Technical Report NSSDC/WDC-A-R&S 91-24, NASA Goddard Space Flight Center (1991)

  13. [13]

    IEEE Transactions on Nuclear Science46(6), 1475–1480 (1999) https://doi.org/10.1109/23.819110

    Heynderickx, D., Kruglanski, M., Pierrard, V., Lemaire, J., Looper, M.D., Blake, J.B.: A low altitude trapped proton model for solar minimum conditions based on SAMPEX/PET data. IEEE Transactions on Nuclear Science46(6), 1475–1480 (1999) https://doi.org/10.1109/23.819110

  14. [14]

    Space Science Review179(1-4), 579–615 (2013) https://doi.org/10.1007/s11214-013-9964-y

    Ginet, G.P., O’Brien, T.P., Huston, S.L., Johnston, W.R., Guild, T.B., Friedel, R., Lindstrom, C.D., Roth, C.J., Whelan, P., Quinn, R.A., Madden, D., Morley, S., Su, Y.-J.: AE9, AP9 and SPM: New Models for Specifying the Trapped Energetic Particle and Space Plasma Environment. Space Science Review179(1-4), 579–615 (2013) https://doi.org/10.1007/s11214-013-9964-y

  15. [15]

    2022, in Handbook of X-ray and Gamma-ray Astrophysics, ed

    Fiore, F., Guzman, A., Campana, R., Evangelista, Y.: HERMES-Pathfinder. In: Bambi, C., Sangangelo, A. (eds.) Handbook of X-ray and Gamma-ray Astrophysics, p. 38 (2022). https://doi.org/10.1007/978-981-16-4544-0 35-1

  16. [16]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Evangelista, Y., Fiore, F., Campana, R., Ceraudo, F., Della Casa, G., Demenev, E., Dilillo, G., Fiorini, M., Grassi, M., Guzman, A., Hedderman, P., Marchesini, E.J., Morgante, G., Mele, F., Nogara, P., Nuti, A., Piazzolla, R., Pliego Caballero, S., Rashevskaya, I., Russo, F., Sottile, G., Labanti, C., Baroni, G., Bellutti, P., Bertuccio, G., Cao, J., Chen...

  17. [17]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Evangelista, Y., Fiore, F., Campana, R., Baroni, G., Ceraudo, F., Della Casa, G., Demenev, E., Dilillo, G., Fiorini, M., Ghirlanda, G., Grassi, M., Guzm´ an, A., Hedderman, P., Marchesini, E.J., Morgante, G., Mele, F., Nava, L., Nogara, P., Nuti, A., Pliego Caballero, S., Rashevskaya, I., Russo, F., Sottile, G., Lavagna, M., Colagrossi, A., Silvestrini, S...

  18. [18]

    In: 2019 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), pp

    Gandola, M., Grassi, M., Mele, F., Malcovati, P., Bertuccio, G.: LYRA: A Multi- Chip ASIC Designed for HERMES X and Gamma Ray Detector. In: 2019 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), pp. 1–3 (2019). https://doi.org/10.1109/NSS/MIC42101.2019.9059616

  19. [19]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Guzman, A., Pliego, S., Bayer, J., Evangelista, Y., La Rosa, G., Sottile, G., Curzel, S., Campana, R., Fiore, F., Fuschino, F., Colagrossi, A., Fiorito, M., Nog- ara, P., Piazzolla, R., Russo, F., Santangelo, A., Tenzer, C.: The Payload Data Handling Unit (PDHU) on-board the HERMES-TP and HERMES-SP CubeSat Missions. In: den Herder, J.-W.A., Nikzad, S., Na...

  20. [20]

    Experimental Astronomy 58(3), 14 (2024) https://doi.org/10.1007/s10686-024-09958-4 arXiv:2410.06056 [astro-ph.IM]

    Dilillo, G., Marchesini, E.J., Baroni, G., Della Casa, G., Campana, R., Evan- gelista, Y., Guzm´ an, A., Hedderman, P., Bellutti, P., Bertuccio, G., Ceraudo, F., Citossi, M., Cirrincione, D., Dedolli, I., Demenev, E., Feroci, M., Ficorella, F., Fiorini, M., Gandola, M., Grassi, M., Rosa, G.L., Lombardi, G., Malco- vati, P., Mele, F., Nogara, P., Nuti, A.,...

  21. [21]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    Campana, R., Evangelista, Y., Fiore, F., Guzm´ an, A., Baroni, G., Della Casa, G., Dilillo, G., Hedderman, P., Marchesini, E.J., Bertuccio, G., Ceraudo, F., Demenev, E., Fiorini, M., Grassi, M., Malcovati, P., Mele, F., Nogara, P., Nuti, A., Perri, M., Pirrotta, S., Pliego-Caballero, S., Puccetti, S., Sottile, G., Russo, F., Trevisan, S.: Design and devel...

  22. [22]

    arXiv e-prints, 2407–14034 (2024) https://doi.org/10.48550/arXiv.2407.14034 arXiv:2407.14034 [astro-ph.IM]

    Trenti, M., Ortiz del Castillo, M., Mearns, R., McRobbie, J., Therakam, C., Chapman, A., Woods, A., Morgan, J., Barraclough, S., Rodriguez Mallo, I., Baroni, G., Fiore, F., Evangelista, Y., Campana, R., Guzman, A., Hedder- man, P.: SpIRIT Mission: In-Orbit Results and Technology Demonstrations. arXiv e-prints, 2407–14034 (2024) https://doi.org/10.48550/ar...

  23. [23]

    In: den Herder, J.-W.A., Nikzad, 29 S., Nakazawa, K

    Campana, R., Fuschino, F., Evangelista, Y., Dilillo, G., Fiore, F.: The HERMES- TP/SP background and response simulations. In: den Herder, J.-W.A., Nikzad, 29 S., Nakazawa, K. (eds.) Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 11444, p. 11444...

  24. [24]

    In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K

    ˇR´ ıpa, J., Dilillo, G., Campana, R., Galg´ oczi, G.: A comparison of trapped particle models in low Earth orbit. In: den Herder, J.-W.A., Nikzad, S., Nakazawa, K. (eds.) Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 11444, p. 114443 (2021). h...

  25. [25]

    Della Casa, G., Zampa, N., Cirrincione, D., Monzani, S., Baruzzo, M., Cam- pana, R., Cauz, D., Citossi, M., Crupi, R., Dilillo, G., Pauletta, G., Fiore, F., Vacchi, A.: New detailed characterization of the residual luminescence emit- ted by the GAGG:Ce scintillator crystals for the HERMES Pathfinder mission. Nuclear Instruments and Methods in Physics Rese...

  26. [26]

    The enhanced X-ray Timing and Polarimetry mission - eXTP

    Zhang, S., Santangelo, A., Feroci, M., Xu, Y., Lu, F., Chen, Y., Feng, H., Zhang, S., Brandt, S., Hernanz, M., Baldini, L., Bozzo, E., Campana, R., De Rosa, A., Dong, Y., Evangelista, Y., Karas, V., Meidinger, N., Meuris, A., Nandra, K., Pan, T., Pareschi, G., Orleanski, P., Huang, Q., Schanne, S., Sironi, G., Spiga, D., Svoboda, J., Tagliaferri, G., Tenz...

  27. [27]

    The Large Observatory for X-ray Timing (LOFT)

    Feroci, M., Stella, L., van der Klis, M., Courvoisier, T.J.-L., Hernanz, M., Hudec, 30 R., Santangelo, A., Walton, D., Zdziarski, A., Barret, D., Belloni, T., Braga, J., Brandt, S., Budtz-Jørgensen, C., Campana, S., den Herder, J.-W., Huovelin, J., Israel, G.L., Pohl, M., Ray, P., Vacchi, A., Zane, S., Argan, A., Attin` a, P., Bertuccio, G., Bozzo, E., Ca...

  28. [28]

    Science China Physics, Mechanics, and Astronomy68(11), 119502 (2025) https://doi.org/10.1007/s11433-025-2786-6 arXiv:2506.08101 [astro-ph.HE]

    Zhang, S.-N., Santangelo, A., Xu, Y., Feng, H., Lu, F., Chen, Y., Ge, M., Nandra, K., Wu, X., Feroci, M., Hernanz, M., Liu, C., He, H., Wang, Y., Jiang, W., Cui, W., Yang, Y., Wang, J., Li, W., Li, H., Du, Y., Liu, X., Meng, B., Wen, X., Zhang, A., Ma, J., Li, M., Li, G., Qi, L., Sun, J., Luo, T., Liu, H., Liu, X., Zhang, F., Luo, L., Zhu, Y., Zhao, Z., S...

  29. [29]

    Background simulations for the Large Area Detector onboard LOFT

    Campana, R., Feroci, M., Del Monte, E., Mineo, T., Lund, N., Fraser, G.W.: Background simulations for the Large Area Detector onboard LOFT. Experimen- tal Astronomy36(3), 451–477 (2013) https://doi.org/10.1007/s10686-013-9341-6 arXiv:1305.3789 [astro-ph.IM]

  30. [30]

    Proceedings of the Cambridge Philosophical Society15(pt V), 423–427 (1910)

    Bateman, H.: The solution of a system of differential equations occurring in the theory of radioactive transformations. Proceedings of the Cambridge Philosophical Society15(pt V), 423–427 (1910)

  31. [31]

    Ameri- can Journal of Physics70(4), 444–445 (2002) https://doi.org/10.1119/1.1427084

    Pressyanov, D.S.: Short solution of the radioactive decay chain equations. Ameri- can Journal of Physics70(4), 444–445 (2002) https://doi.org/10.1119/1.1427084

  32. [32]

    American Journal of Physics37(10), 1019–1022 (1969) https://doi.org/10.1119/1.1975187

    Onega, R.J.: Radioactivity Calculations. American Journal of Physics37(10), 1019–1022 (1969) https://doi.org/10.1119/1.1975187

  33. [33]

    American Journal of Physics71(7), 684–686 (2003) https://doi.org/10.1119/1

    Moral, L., Pacheco, A.F.: Algebraic approach to the radioactive decay equations. American Journal of Physics71(7), 684–686 (2003) https://doi.org/10.1119/1. 1571834

  34. [34]

    Computer Physics Communications181(1), 21–23 (2010) https://doi.org/10.1016/j.cpc.2009.08.011

    Amaku, M., Pascholati, P.R., Vanin, V.R.: Decay chain differential equations: Solution through matrix algebra. Computer Physics Communications181(1), 21–23 (2010) https://doi.org/10.1016/j.cpc.2009.08.011

  35. [35]

    IEEE Transactions on Nuclear Science60(4), 2966–2983 (2013) https://doi.org/ 10.1109/TNS.2013.2270894 arXiv:1307.0996 [physics.comp-ph] 32

    Hauf, S., Kuster, M., Batiˇ c, M., Bell, Z.W., Hoffmann, D.H.H., Lang, P.M., Neff, S., Pia, M.G., Weidenspointner, G., Zoglauer, A.: Radioactive Decays in Geant4. IEEE Transactions on Nuclear Science60(4), 2966–2983 (2013) https://doi.org/ 10.1109/TNS.2013.2270894 arXiv:1307.0996 [physics.comp-ph] 32