Using radial distribution functions to calculate cellular cross-absorbed dose forßemitters: comparison with reference methods and application for18F-FDG cell labeling.

To further improve the understanding ofin vitrobiological effects of incorporated radionuclides, it is essential to accurately determine cellular absorbed doses. In the case ofßemitters, the cross-dose is a major contribution, and can involve up to millions of cells. Realistic and efficient computational models are needed for that purpose. Conventionally, distances between each cell are calculated and the related dose contributions are cumulated to get the total cross-dose (standard method). In this work, we developed a novel approach for the calculation of the cross-absorbed dose, based on the use of the radial distribution function (rdf)) that describes the spatial properties of the cellular model considered. The dynamic molecular tool LAMMPS was used to create 3D cellular models and computerdfsfor various conditions of cell density, volume size, and configuration type (lattice and randomized geometry). The novel method is suitable for any radionuclide of nuclear medicine. Here, the model was applied for the labeling of cells with18F-FDG used for PET imaging, and first validated by comparison with other reference methods. MeanScrossvalues calculated with the novel approach versus the standard method agreed very well (relative differences less that 0.1%). Implementation of therdf-based approach with LAMMPS allowed to achieved results considerably faster than with the standard method, the computing time decreasing from hours to seconds for 106cells. Therdf-based approach was also faster and easier to accommodate more complex cellular models than the standard and other published methods. Finally, a comparative study of the meanScrossfor different types of configuration was carried out, as a function of the cell density and the volume size, allowing to better understand the impact of the configuration on the cross-absorbed dose.

Uncertainty induced by chest wall thickness assessment methods on lung activity estimation for plutonium and americium: a large population-based study.

In vivo lung counting aims at assessing the retained activity in the lungs. The calibration factor relating the measured counts to the worker's specific retained lung activity can be obtained by several means and strongly depends on the chest wall thickness. Here we compare, for 374 male nuclear workers, the activity assessed with a reference protocol, where the material equivalent chest wall thickness is known from ultrasound measurements, with two other protocols. The counting system is an array of four germanium detectors.It is found that non site-specific equations for the assessment of the chest wall thickness induce large biases in the assessment of activity. For plutonium isotopes or (241)Am the proportion of workers for whom the retained activity is within ± 10% of the reference one is smaller than 10%.The use of site-specific equations raises this proportion to 20% and 58% for plutonium and (241)Am, respectively.Finally, for the studied population, when site-specific equations are used for the chest wall thickness, the standard uncertainties for the lung activity are 42% and 12.5%, for plutonium and (241)Am, respectively. Due to the relatively large size of the studied population, these values are a relatively robust estimate of the uncertainties due to the assessment of the chest wall thickness for the current practice at this site.

Microdosimetry of alpha particles for simple and 3D voxelised geometries using MCNPX and Geant4 Monte Carlo codes.

Microdosimetry using Monte Carlo simulation is a suitable technique to describe the stochastic nature of energy deposition by alpha particle at cellular level. Because of its short range, the energy imparted by this particle to the targets is highly non-uniform. Thus, to achieve accurate dosimetric results, the modelling of the geometry should be as realistic as possible. The objectives of the present study were to validate the use of the MCNPX and Geant4 Monte Carlo codes for microdosimetric studies using simple and three-dimensional voxelised geometry and to study their limit of validity in this last case. To that aim, the specific energy (z) deposited in the cell nucleus, the single-hit density of specific energy f(1)(z) and the mean-specific energy were calculated. Results show a good agreement when compared with the literature using simple geometry. The maximum percentage difference found is <6 %. For voxelised phantom, the study of the voxel size highlighted that the shape of the curve f(1)(z) obtained with MCNPX for <1 µm voxel size presents a significant difference with the shape of non-voxelised geometry. When using Geant4, little differences are observed whatever the voxel size is. Below 1 µm, the use of Geant4 is required. However, the calculation time is 10 times higher with Geant4 than MCNPX code in the same conditions.

EURADOS intercomparison on measurements and Monte Carlo modelling for the assessment of americium in a USTUR leg phantom.

A collaboration of the EURADOS working group on 'Internal Dosimetry' and the United States Transuranium and Uranium Registries (USTUR) has taken place to carry out an intercomparison on measurements and Monte Carlo modelling determining americium deposited in the bone of a USTUR leg phantom. Preliminary results and conclusions of this intercomparison exercise are presented here.

EURADOS coordinated action on research, quality assurance and training of internal dose assessments.

EURADOS working group on 'Internal Dosimetry (WG7)' represents a frame to develop activities in the field of internal exposures as coordinated actions on quality assurance (QA), research and training. The main tasks to carry out are the update of the IDEAS Guidelines as a reference document for the internal dosimetry community, the implementation and QA of new ICRP biokinetic models, the assessment of uncertainties related to internal dosimetry models and their application, the development of physiology-based models for biokinetics of radionuclides, stable isotope studies, biokinetic modelling of diethylene triamine pentaacetic acid decorporation therapy and Monte-Carlo applications to in vivo assessment of intakes. The working group is entirely supported by EURADOS; links are established with institutions such as IAEA, US Transuranium and Uranium Registries (USA) and CEA (France) for joint collaboration actions.

Female workers and in vivo lung monitoring: a simple model for morphological dependence of counting efficiency curves.

This paper addresses the question of the morphological dependence of counting efficiency curves for in vivo lung monitoring of workers, with a particular focus on the case of female workers for whom different chest girth and cup size are considered. A library of 24 female torsos, with chest girth varying from 85 to 120 and cup size from A to F, was constructed using mesh and NURBS formats. The anatomical realism and usefulness of these models for simulating in vivo counting measurements are illustrated and simulations are reported for a typical 4-germanium (Ge) counting system. A simple analytic formula describing the relation between efficiency curves obtained for each female phantom is given. This formula uses the mass attenuation coefficient for adipose tissue and two parameters which are dependant on lung volume and breast weight. The model is tested against Monte Carlo simulated data, experimental data obtained with the Livermore phantom and published data. The model correctly describes the efficiency curve and, since the parameters depend on the counting geometry, it is shown how to estimate them from experimental measurements.

Application of the ICRP/ICRU reference computational phantoms to internal dosimetry: calculation of specific absorbed fractions of energy for photons and electrons.

The emission of radiation from a contaminated body region is connected with the dose received by radiosensitive tissue through the specific absorbed fractions (SAFs) of emitted energy, which is therefore an essential quantity for internal dose assessment. A set of SAFs were calculated using the new adult reference computational phantoms, released by the International Commission on Radiological Protection (ICRP) together with the International Commission on Radiation Units and Measurements (ICRU). Part of these results has been recently published in ICRP Publication 110 (2009 Adult reference computational phantoms (Oxford: Elsevier)). In this paper, we mainly discuss the results and also present them in numeric form. The emission of monoenergetic photons and electrons with energies ranging from 10 keV to 10 MeV was simulated for three source organs: lungs, thyroid and liver. SAFs were calculated for four target regions in the body: lungs, colon wall, breasts and stomach wall. For quality assurance purposes, the simulations were performed simultaneously at the Helmholtz Zentrum München (HMGU, Germany) and at the Institute for Radiological Protection and Nuclear Safety (IRSN, France), using the Monte Carlo transport codes EGSnrc and MCNPX, respectively. The comparison of results shows overall agreement for photons and high-energy electrons with differences lower than 8%. Nevertheless, significant differences were found for electrons at lower energy for distant source/target organ pairs. Finally, the results for photons were compared to the SAF values derived using mathematical phantoms. Significant variations that can amount to 200% were found. The main reason for these differences is the change of geometry in the more realistic voxel body models. For electrons, no SAFs have been computed with the mathematical phantoms; instead, approximate formulae have been used by both the Medical Internal Radiation Dose committee (MIRD) and the ICRP due to the limitations imposed by the computing power available at this time. These approximations are mainly based on the assumption that electrons are absorbed locally in the source organ itself. When electron SAFs are calculated explicitly, discrepancies with this simplifying assumption are notable, especially at high energies and for neighboring organs where the differences can reach the same order of magnitude as for photon SAFs.

Analytical and Monte Carlo assessment of activity and local dose after a wound contamination by activation products.

The activity and local dose following a right index finger wound contamination by activation products are assessed. Measurements with a high purity germanium detector and a four positions measurement protocol enabled a better localization of the contaminant source. From the source location and detector calibration, the remaining wound activity and local absorbed dose were deduced. An analytical model, based on a two dimensional simplification of the problem, is presented. It is shown to provide a fast and quite accurate activity assessment when the contaminants are described as a point source. The contaminants' location and activity were then more accurately assessed using Monte Carlo calculations based on the OEDIPE software and a voxelized phantom of the index finger. Describing the contaminant mixture as a point source resulted in an agreement of experimental and computed data around 6% for most of the radionuclides. The total activity, due to 11 radionuclides, was estimated to be (9.5 +/- 0.4) kBq at measurement day. Since the point source is found to be less than 1 mm under the skin, the equivalent skin dose is calculated and found to be around 680 mSv in the first year after the contamination, and this value decreases to 250 mSv in the second year. The relevance of equivalent skin dose as an estimate of the sanitary risk is discussed, and it is concluded that for this case it gives the upper end estimate of the risk.

Internal dosimetry: towards harmonisation and coordination of research.

The CONRAD Project is a Coordinated Network for Radiation Dosimetry funded by the European Commission 6th Framework Programme. The activities developed within CONRAD Work Package 5 ('Coordination of Research on Internal Dosimetry') have contributed to improve the harmonisation and reliability in the assessment of internal doses. The tasks carried out included a study of uncertainties and the refinement of the IDEAS Guidelines associated with the evaluation of doses after intakes of radionuclides. The implementation and quality assurance of new biokinetic models for dose assessment and the first attempt to develop a generic dosimetric model for DTPA therapy are important WP5 achievements. Applications of voxel phantoms and Monte Carlo simulations for the assessment of intakes from in vivo measurements were also considered. A Nuclear Emergency Monitoring Network (EUREMON) has been established for the interpretation of monitoring data after accidental or deliberate releases of radionuclides. Finally, WP5 group has worked on the update of the existing IDEAS bibliographic, internal contamination and case evaluation databases. A summary of CONRAD WP5 objectives and results is presented here.

Analysis of a computational problem involving complex voxel geometries.

This communication briefly summarises the results obtained from the 'International comparison on MC modeling for in vivo measurement of Americium in a knee phantom' organised within the EU Coordination Action CONRAD (Coordinated Network for Radiation Dosimetry) as a joint initiative of EURADOS working groups 6 (computational dosimetry) and 7 (internal dosimetry). Monte Carlo simulations using the knee voxel phantom proved to be a viable approach to provide the calibration factor needed for in vivo measurements.

OEDIPE: a new graphical user interface for fast construction of numerical phantoms and MCNP calculations.

Although great efforts have been made to improve the physical phantoms used to calibrate in vivo measurement systems, these phantoms represent a single average counting geometry and usually contain a uniform distribution of the radionuclide over the tissue substitute. As a matter of fact, significant corrections must be made to phantom-based calibration factors in order to obtain absolute calibration efficiencies applicable to a given individual. The importance of these corrections is particularly crucial when considering in vivo measurements of low energy photons emitted by radionuclides deposited in the lung such as actinides. Thus, it was desirable to develop a method for calibrating in vivo measurement systems that is more sensitive to these types of variability. Previous works have demonstrated the possibility of such a calibration using the Monte Carlo technique. Our research programme extended such investigations to the reconstruction of numerical anthropomorphic phantoms based on personal physiological data obtained by computed tomography. New procedures based on a new graphical user interface (GUI) for development of computational phantoms for Monte Carlo calculations and data analysis are being developed to take advantage of recent progress in image-processing codes. This paper presents the principal features of this new GUI. Results of calculations and comparison with experimental data are also presented and discussed in this work.

Monte Carlo modelling for in vivo measurements of Americium in a knee voxel phantom: general criteria for an international comparison.

The general criteria and the scientific approach adopted for an 'International comparison on Monte Carlo modelling for in vivo measurement of Americium in a knee phantom' that is being organised within the EU Coordination Action CONRAD (Coordinated Network for Radiation Dosimetry) are described her. Detection system and a knee voxel phantom based on a computerised axial tomography of the Spitz anthropometric knee phantom with a homogeneous distribution of 241Am in bone have been considered for the simulation of three specific situations: (a) a single Low Energy Germanium detector for a point 241Am source in air; (b) the calculation of photon fluence spectra in air around the voxel phantom; and (c) the calculation of the energy distribution of pulses and peak detection efficiency in the real detection system geometry.

Automatic application of ICRP biokinetic models in voxel phantoms for in vivo counting and internal dose assessment.

As part of the improvement of calibration techniques of in vivo counting, the Laboratory of Internal Dose Assessment of the Institute of Radiological Protection and Nuclear Safety has developed a computer tool, 'OEDIPE', to model internal contamination, to simulate in vivo counting and to calculate internal dose. The first version of this software could model sources located in a single organ. As the distribution of the contamination evolves from the time of intake according to the biokinetics of the radionuclide, a new facility has been added to the software first to allow complex heterogeneous source modelling and then to automatically integrate the distribution of the contamination in the different tissues estimated by biokinetic calculation at any time since the intake. These new developments give the opportunity to study the influence of the biokinetics on the in vivo counting, leading to a better assessment of the calibration factors and the corresponding uncertainties.

Implementation of bioassay methods to improve assessment of incorporated radionuclides.

The present work which was carried out in the framework of an EU project (IDEA: Internal Dosimetry-Enhancements in Application; Contract Number: FIKR CT2001 00164) shall provide commonly acceptable guidelines for optimum performance of ICP-MS measurements with focus on urinary measurements of uranium, thorium and actinides. From the results of this work it is recommended that, whenever feasible, 24 h urine sampling should be conducted to avoid large uncertainties in the quantitation of daily urinary excretion values. For storage, urine samples should be acidified and kept frozen before analysis. Measurement of total uranium in urine by ICP-MS at physiological levels (<10 ng.l(-1)) requires no sample preparation besides UV photolysis and/or dilution. For the measurement of thorium in urine by ICP-MS, it can be concluded, that salt removal from the urine samples is not recommended. For the measurement of actinides in urine it is shown that ICP-MS is well-suited and a good alternative to alpha-spectrometry for isotopes with T1/2>5x10(4) years. In general, ICP-MS measurements are an easy, fast and cost-saving methodology. New improved measuring techniques (HR-SF-ICP-MS) with detection limits in urine of 150 pg.l(-1) (1.9 microBq.l(-1)) for 238U, 30 pg.l(-1) (2.4 microBq.l(-1)) for 235U and 100 pg.l(-1) (0.4 microBq.l(-1)) for (232)Th, respectively, meet all necessary requirements. This method should therefore become the routine technique for incorporation monitoring of workers and of members of the general public, in particular for uranium contamination.

Application of voxel phantoms in whole-body counting for the validation of calibration phantoms and the assessment of uncertainties.

This article is dedicated to the application of voxel phantoms in whole-body counting calibration. The first study was performed to validate this approach using IGOR, a physical phantom dedicated to fission and activation product (FAP) measurement, and a graphical user interface, developed at the IRSN internal dose assessment laboratory, called OEDIPE (French acronym for the tool for personalised internal dose assessment) associated with the Monte Carlo code MCNP. The method was validated by comparing the results of real measurements and simulations using voxel phantoms obtained from CT scan images of IGOR. To take this application further, two studies were carried out and are presented in this article. First, a comparison was made between the IGOR voxel based phantom (IGOVOX) and a voxel human body (Zubal Phantom) to confirm whether IGOR could be considered as a realistic representation of a human. Second, the errors made when considering sources homogeneously distributed in the body were assessed against real contamination by taking into account the biokinetic behaviour of the radioactive material for two modes of exposure: the ingestion of 137Cs in soluble form and the inhalation of insoluble 60Co several days after acute incorporation.

Validation of a personalized dosimetric evaluation tool (Oedipe) for targeted radiotherapy based on the Monte Carlo MCNPX code.

Dosimetric studies are necessary for all patients treated with targeted radiotherapy. In order to attain the precision required, we have developed Oedipe, a dosimetric tool based on the MCNPX Monte Carlo code. The anatomy of each patient is considered in the form of a voxel-based geometry created using computed tomography (CT) images or magnetic resonance imaging (MRI). Oedipe enables dosimetry studies to be carried out at the voxel scale. Validation of the results obtained by comparison with existing methods is complex because there are multiple sources of variation: calculation methods (different Monte Carlo codes, point kernel), patient representations (model or specific) and geometry definitions (mathematical or voxel-based). In this paper, we validate Oedipe by taking each of these parameters into account independently. Monte Carlo methodology requires long calculation times, particularly in the case of voxel-based geometries, and this is one of the limits of personalized dosimetric methods. However, our results show that the use of voxel-based geometry as opposed to a mathematically defined geometry decreases the calculation time two-fold, due to an optimization of the MCNPX2.5e code. It is therefore possible to envisage the use of Oedipe for personalized dosimetry in the clinical context of targeted radiotherapy.

New method of voxel phantom creation: application for whole-body counting calibration and perspectives in individual internal dose assessment.

The purpose of this work is to present an innovative approach for the creation and application of voxel phantoms associated with the Monte Carlo calculation (MCNP) for the calibration of whole-body counting systems dedicated to the measurement of fission and activation products. The new method is based on a graphical user interface called 'OEDIPE' that allows to simulate a whole measurement process using all measurement parameters, the final goal being to approach a numerical calibration of the facilities. The creation of voxel phantoms and validation of the method are presented in this paper using the IGOR phantom. Finally, the efficiency of the method is discussed, in particular, with the perspective of validating IGOR as a suitable human-equivalent phantom and for the assessment of uncertainties in dose estimation due to the inhomogeneous distribution of activity in the body, correlated to the bio-kinetic behaviour of the radionuclides.

Dosimetric comparison of Monte Carlo codes (EGS4, MCNP, MCNPX) considering external and internal exposures of the Zubal phantom to electron and photon sources.

This paper aims at comparing dosimetric assessments performed with three Monte Carlo codes: EGS4, MCNP4c2 and MCNPX2.5e, using a realistic voxel phantom, namely the Zubal phantom, in two configurations of exposure. The first one deals with an external irradiation corresponding to the example of a radiological accident. The results are obtained using the EGS4 and the MCNP4c2 codes and expressed in terms of the mean absorbed dose (in Gy per source particle) for brain, lungs, liver and spleen. The second one deals with an internal exposure corresponding to the treatment of a medullary thyroid cancer by 131I-labelled radiopharmaceutical. The results are obtained by EGS4 and MCNPX2.5e and compared in terms of S-values (expressed in mGy per kBq and per hour) for liver, kidney, whole body and thyroid. The results of these two studies are presented and differences between the codes are analysed and discussed.

A computational tool based on voxel geometry for dose reconstruction of a radiological accident due to external exposure.

In the case of overexposure to ionising radiation, estimation of the absorbed dose in the organism is an important indicator for evaluating the biological consequences of this exposure. The physical dosimetry approach is based either on real reconstruction of the accident, using physical phantoms, or on calculation techniques. Tools using Monte Carlo simulations associated with geometric models are very powerful since they offer the possibility to simulate faithfully the victim and the environment for dose calculations in various accidental situations. Their work presents a new computational tool, called SESAME, dedicated to dose reconstruction of radiological accidents based on anthropomorphic voxel phantoms built from real medical images of the victim in association with the MCNP Monte Carlo code. The utility was, as a first step, validated for neutrons by experimental means using a physical tissue-equivalent phantom.

Application of new imaging and calculation techniques to activity and dose assessment in the case of a 106Ru contaminated wound.

The aim of this paper is to describe the dosimetric evaluation of a point contamination that occurred in a laboratory during the examination of an irradiated sample. The incident led to point contamination of the operator's finger due to the presence of mainly 106Ru, with its progeny, 106Rh. The paper reports on the activity and dose assessment, performed using several methods. The measured activity was obtained using a conventional device based on a germanium detector and confirmed using software developed at IRSN, based on reconstruction of voxel phantom associated with the Monte Carlo N-Particle code (MCNP) for in vivo measurement. Two dose assessment calculations were performed using both analytical and Monte Carlo methods, applying the same approach as for activity assessment based on the personal computational phantom of the finger. The results are compared, followed by a discussion on the suitability of the tools described in this study.