An analytical model to calculate absorbed fractions for internal dosimetry with alpha, beta and gamma emitters

We developed a general model for the calculation of absorbed fractions in ellipsoidal volumes of soft tissue uniformly filled with alpha, beta and gamma emitting radionuclides. The approach exploited Monte Carlo simulations with the Geant4 code to determine absorbed fractions in ellipsoids characterized by a wide range of dimensions and ellipticities, for monoenergetic emissions of each radiation type. The so-obtained absorbed fractions were put in an analytical relationship with the ’generalized radius’, calculated as 3V /S, where V is the ellipsoid volume and S its surface. Radiation-specific parametric functions were obtained in order to calculate the absorbed fraction of a given radiation in a generic ellipsoidal volume. The dose from a generic radionuclide can be calculated through a process of summation and integration over the whole radionuclide emission spectrum, profitably implemented in an electronic spreadsheet. We compared the results of our analytical calculation approach with those obtained from the OLINDA/EXM computer software, finding a good agreement in a wide range of sphere radii, for the high-energy pure beta emitter 90 Y, the commonly employed beta-gamma emitter 131 I, and the pure alpha emitter 213 Po. The generality of our approach makes it useful and easy to implement in clinical dosimetry calculations as well as in radiation safety estimations when doses from internal radionuclide uptake are to be taken into account.