ASSESSING THE CONTRIBUTION OF CROSS-SECTIONS TO THE UNCERTAINTY OF MONTE CARLO CALCULATIONS IN MICRO- AND NANODOSIMETRY.

Within EURADOS Working Group 6 'Computational Dosimetry', the micro and nanodosimetry task group 6.2 has recently conducted a Monte Carlo (MC) exercise open to participants around the world. The aim of this exercise is to quantify the contribution to the uncertainty of micro and nanodosimetric simulation results arising from the use of different electron-impact cross-sections, and hence physical models, employed by different MC codes (GEANT4-DNA, PENELOPE, MCNP6, FLUKA, NASIC and PHITS). Comparison of the participants' simulation results for both micro and nanodosimetric quantities using different MC codes was the first step of the exercise. The deviation between results is due to different cross-sections but also different tracking methods and particle transport cut-off energies. The second step of the exercise will involve using identical cross-section datasets to account only for the other variations in the first step, thus enabling the determination of the uncertainty contribution due to different cross-sections. This paper presents a comparison of the MC simulation results obtained in the first part of the exercise. For the microdosimetric simulations, particularly in the configuration where the electron source is contained within the micrometric target, the choice of MC code has a small influence on the results. For the nanodosimetric results, on the other hand, the mean ionisation cluster size distribution (ICSD) was sensitive to the physical models used in the MC codes. The ICSD was therefore chosen to study the influence of different cross-section data on the uncertainty of simulation results.

AN ALGORITHM TO DETERMINE THE NANODOSIMETRIC IMPACT OF GOLD NANOPARTICLES ON CELL MODELS.

High-Z nanomaterials, e.g. gold nanoparticles (GNPs), are being investigated worldwide for potential application in radiation imaging and therapy. Photon irradiation of cells containing GNP was shown to produce enhanced DNA damage which is believed to be related to the increased secondary electron (SE) yield and ionization density. In this work, an algorithm was developed for simulating the physical radiation damage inside the nucleus of a spherical cell model for the case of uniformly distributed GNPs within the cytoplasm. Previously calculated energy spectra of SE emerging from a single NP irradiated with different photon sources are used as input to obtain the SE energy spectrum at the surface of the cell nucleus. In a second step, the SE transport inside the cell nucleus is simulated with a track structure Monte Carlo code to obtain the spatial distribution of ionizations. The preliminary results presented here show that the developed algorithm allows for a fast calculation of the SE spectra at the cell nucleus surface, thus enabling a more realistic assessment of the ionization density inside the cell nucleus than that obtained by the simulation of a single GNP. Furthermore, the algorithm can be easily adapted to investigate both the effect of GNP clustering and the impact of GNP-GNP interactions on SE spectra.

Proton-impact ionisation cross sections for nanodosimetric track structure simulations.

Monte Carlo simulations of the particle track structure require accurate ion- and electron-impact cross-section data of the medium. These data are scarce and often inconsistent when measured by different groups. In this work, literature data on ionisation cross sections (CSs) of nitrogen and propane for protons with energies 0.1-10 MeV are reviewed and implemented in the code PTra. Methane data were used to obtain proton-impact CSs of propane due to their absence in the literature. PTra is benchmarked by comparing simulated particle-track parameters to experimental results, measured with an ion-counting nanodosemeter.

Comparison of measured and Monte Carlo simulated track structure parameters in nanometric volumes.

The track structure of ionising particles in biological matter can only be assessed by simulations, since neither the type of interaction and its products nor the interaction positions in biological matter can be detected with nanometer resolution. Hence, there is a need to benchmark the deployed computer codes using suitable experimental data. For this purpose, the frequency distributions of ionisation clusters produced in the sensitive volume of the PTB ion counting nanodosemeter by monoenergetic protons and alpha particles (with energies between 0.1 and 20 MeV) were measured. C3H8 and N2 were alternately used as the working gas. The measured data were compared with the results of simulations obtained with the PTB Monte Carlo code PTra. Measured and simulated characteristics of the particle track structure are in good agreement for protons over the entire energy range investigated. For alpha particles with energies above the Bragg peak a good agreement can also be seen, whereas for energies below the Bragg peak differences of as much as 25 % occur.

Comparison of nanodosimetric parameters of track structure calculated by the Monte Carlo codes Geant4-DNA and PTra.

The concept of nanodosimetry is based on the assumption that initial damage to cells is related to the number of ionizations (the ionization cluster size) directly produced by single particles within, or in the close vicinity of, short segments of DNA. The ionization cluster-size distribution and other nanodosimetric quantities, however, are not directly measurable in biological targets and our current knowledge is mostly based on numerical simulations of particle tracks in water, calculating track structure parameters for nanometric target volumes. The assessment of nanodosimetric quantities derived from particle-track calculations using different Monte Carlo codes plays, therefore, an important role for a more accurate evaluation of the initial damage to cells and, as a consequence, of the biological effectiveness of ionizing radiation. The aim of this work is to assess the differences in the calculated nanodosimetric quantities obtained with Geant4-DNA as compared to those of the ad hoc particle-track Monte Carlo code 'PTra' developed at Physikalisch-Technische Bundesanstalt (PTB), Germany. The comparison of the two codes was made for incident electrons of energy in the range between 50 eV and 10 keV, for protons of energy between 300 keV and 10 MeV, and for alpha particles of energy between 1 and 10 MeV as these were the energy ranges available in both codes at the time this investigation was carried out. Good agreement was found for nanodosimetric characteristics of track structure calculated in the high-energy range of each particle type. For lower energies, significant differences were observed, most notably in the estimates of the biological effectiveness. The largest relative differences obtained were over 50%; however, generally the order of magnitude was between 10% and 20%.

Microdosimetric assessment of the radiation quality of a therapeutic proton beam: comparison between numerical simulation and experimental measurements.

Using protons for the treatment of ocular melanoma (especially of posterior pole tumours), the radiation quality of the beam must be precisely assessed to preserve the vision and to minimise the damage to healthy tissue. The radiation quality of a therapeutic proton beam at the Centre Antoine Lacassagne in Nice (France) was measured using microdosimetric techniques, i.e. a miniaturised version of a tissue-equivalent proportional counter. Measurements were performed in a 1-µm site at different depths in a Lucite phantom. Experimental data showed a significant increase in the beam quality at the distal edge of the spread-out Bragg peak (SOBP). In this paper, the numerical simulation of the experimental setup is done with the FLUKA Monte Carlo radiation transport code. The calculated microdosimetric spectra are compared with the measured ones at different depths in tissue for a monoenergetic proton beam (E=62 MeV) and for a modulated SOBP. Numerically and experimentally predicted relative biological effectiveness values are in good agreement. The calculated frequency-averaged and dose-averaged lineal energy mean values are consistent with measured data.

Proton-induced frequency distributions of ionization cluster size in propane.

The frequency distribution of clustered ionizations produced by a proton beam was measured in a nanodosimetric volume of the size of a DNA segment by means of an ion-counting nanodosimeter in the energy range from 0.4 to 3.5 MeV. In order to meet the needs of the ion-counting nanodosimeter, the accelerator's primary beam was reduced in intensity by means of Rutherford scattering. The comparison between experimental results and Monte Carlo simulations show a good agreement in the energy dependence of the mean cluster size, while the experimental cluster size distributions show a higher amount of large ionization clusters compared with those obtained with the simulations.

Dependence of nanodosimetric spectra on the sensitive volume length and ion drift in an ion-counting nanodosemeter.

Nanodosimetric spectra, measured in a well-defined ionisation sensitive volume of an ion-counting gaseous nanodosemeter, may have a valuable predictive value of radiation damage to DNA. In such devices, the distributions of radiation-induced ions are measured after their drift in gas. The sensitive-volume size, corresponding to a DNA segment length, can be tuned by selecting an appropriate time window for ion counting; the method's accuracy depends on the velocity distribution of the drifting ions. The results of ion-drift measurements in an ion-counting nanodosemeter were used for the precise calculation of its sensitive volume length. Monte Carlo simulations of nanodosimetric spectra, performed with the obtained data, are in good agreement with experimental data. The method's limitations, arising from the spread of drift velocities, are discussed.

A simple method for measuring thoron spatial distributions.

A simple but effective method that allows the measurement of the 220Rn spatial distribution in working or living environments using a solid-state detector is presented in this paper. The method is based on measurements of the alpha particles emitted by 216Po (the first 220Rn progeny) directly deposited on the detector surface at different distances from a 220Rn exhalation source. The validity of the method is shown by comparing the results of an experiment, where the 220Rn activity concentration is measured under conditions of diffusion at constant temperature, with finite-element calculations.

Mixed high energy photon and electron radiation fields for calibrating radiation protection dosemeters.

According to ISO 4037-3, calibrations of radiation protection dosemeters with photon radiation of energies above 3 MeV are performed under conditions of charged particle equilibrium. No information is provided concerning how to determine the response of dosemeters to radiation fields in the more general case when these conditions are not fulfilled. This paper deals with the production of mixed high energy photon and electron fields characterised by a lack or an excess of charged particles relative to conditions of equilibrium and describes a new procedure for the dosimetry in such fields. Through variation of the charged particle fluence fraction with respect to a nearly constant photon fluence, Hp(10) and H'(10) values varied by up to a factor of 1.74. The above mentioned basic study was utilised in the recent IAEA intercomparison (Co-ordinated Research Project 1996-1998) and EURADOS 'trial performance test' (1996-1998) for individual monitoring of photon radiation in testing response characteristics of individual dosemeters in non-charged particle equilibrium conditions.

Response of radiation protection dosemeters in mixed high-energy photon and electron radiation fields.

The response of radiation protection dosemeters in terms of the phantom-related operational quantities Hp(10) and H'(10.0 degrees) was measured for personal and area monitoring systems in mixed high-energy electron and photon radiation fields with energies up to 7 MeV. Using mixed radiation fields composed of different fractions of charged particle and photon fluence, three conditions were produced at the point of measurement: charged particle equilibrium (CPE) (a), a lack (b) and an excess (c) of charged particles relative to the conditions of CPE. Personal and area dosemeters of different types were investigated under conditions (a)-(c). A large variability of the response of the different dosemeter types was observed. The results are presented and discussed.