Simulation of secondary particle production and absorbed dose to tissue in light ion beams

During radiation therapy with an ion beam, the production of secondary particles like neutrons, protons and heavier ions contribute to the dose delivered to tumour and healthy tissues outside the treated volume. Also, the secondary particles leaving the patient are of interest for radiation background around the ion-therapy facility. Calculations of secondary particle production and the dose absorbed by water, soft tissue and a multi-material phantom simulating the heterogeneous media of the patient body were performed for protons, helium, lithium and carbon ions in the energy range up to 400 MeV u(-1). The Monte Carlo code SHIELD-HIT for transport of protons and light ions in tissue-like media was used in these studies. The neutron ambient dose-equivalent, H*(10), was determined for neutrons leaving the water phantom irradiated with different light ion beams. The comparison of calculated secondary particle production in the water and PMMA phantoms irradiated with helium and carbon ions shows satisfactory agreement with experimental data.     

An accelerator-based epithermal neutron beam design for BNCT and dosimetric evaluation using a voxel head phantom

The beam shaping assembly design has been investigated in order to improve the epithermal neutron beam for accelerator-based boron neutron capture therapy in intensity and quality, and dosimetric evaluation for the beams has been performed using both mathematical and voxel head phantoms with MCNP runs. The neutron source was assumed to be produced from a conventional 2.5 MeV proton accelerator with a thick (7)Li target. The results indicate that it is possible to enhance epithermal neutron flux remarkably as well as to embody a good spectrum shaping to epithermal neutrons only with the proper combination of moderator and reflector. It is also found that a larger number of thermal neutrons can reach deeply into the brain and, therefore, can reduce considerably the treatment time for brain tumours. Consequently, the epithermal neutron beams designed in this study can treat more effectively deep-seated brain tumours.     

Comparison doses of secondary neutron with the heavy ions in a 75-Mev/n heavy ion beam

The angular distributions for neutrons of energy >6 MeV that are induced by 75 MeV/n 12C6+ and 16O8+ ions were measured with the activation method of Al threshold detectors at the radiobiological terminal of HIRFL. The data were obtained by a high-purity Ge(HpGe) detector. The results show that the neutron angular distributions produced by heavy ion beams are strongly peaked in the forward direction and decreased exponentially with angles in experimental area. The experimental conditions for these measurements were similar to those for biological experiments, so the results should be representative of neutrons produced by heavy ions during the biological experiments and tumour therapy. Comparing with the neutron doses produced by the heavy ion beam, the heavy ion dose is the main factor in biological effects and tumour therapy response, so the contribution of neutron dose can be neglected.     

Dosimetry of clinical neutron and proton beams: an overview of recommendations

Neutron therapy beams are obtained by accelerating protons or deuterons on Beryllium. These neutron therapy beams present comparable dosimetric characteristics as those for photon beams obtained with linear accelerators; for instance, the penetration of a p(65)+Be neutron beam is comparable with the penetration of an 8 MV photon beam. In order to be competitive with conventional photon beam therapy, the dosimetric characteristics of the neutron beam should therefore not deviate too much from the photon beam characteristics. This paper presents a brief summary of the neutron beams used in radiotherapy. The dosimetry of the clinical neutron beams is described. Finally, recent and future developments in the field of physics for neutron therapy is mentioned. In the last two decades, a considerable number of centres have established radiotherapy treatment facilities using proton beams with energies between 50 and 250 MeV. Clinical applications require a relatively uniform dose to be delivered to the volume to be treated, and for this purpose the proton beam has to be spread out, both laterally and in depth. The technique is called 'beam modulation' and creates a region of high dose uniformity referred to as the 'spread-out Bragg peak'. Meanwhile, reference dosimetry in these beams had to catch up with photon and electron beams for which a much longer tradition of dosimetry exists. Proton beam dosimetry can be performed using different types of dosemeters, such as calorimeters, Faraday cups, track detectors and ionisation chambers. National standard dosimetry laboratories will, however, not provide a standard for the dosimetry of proton beams. To achieve uniformity on an international level, the use of an ionisation chamber should be considered. This paper reviews and summarises the basic principles and recommendations for the absorbed dose determination in a proton beam, utilising ionisation chambers calibrated in terms of absorbed dose to water. These recommendations are based on the recent IAEA TRS398 Code of Practice: 'Absorbed Dose Determination in External Beam Radiotherapy: An International Code of Practice for Dosimetry based on Standards of Absorbed Dose to Water'.     

Analytical shielding calculations for a proton therapy facility

The University of Pennsylvania is building a proton therapy facility in collaboration with Walter Reed Army Medical Center. The proposed facility has four gantry rooms, a fixed beam room and a research room, and will use a cyclotron as the source of protons. In this study, neutron shielding considerations for the proposed proton therapy facility were investigated using analytical techniques and Monte Carlo simulated neutron spectra. Neutron spectra calculations were done using the GEANT4 (v6.2) simulation code for various materials: water, carbon, iron, nickel and tantalum to estimate the neutron production at proton beam energies of 100, 175 and 250 MeV. Dose equivalent calculations were performed using analytical methods at various critical points within the facility, by folding the GEANT4 produced neutron spectra with dose equivalent rate data from the National Council on Radiation Protection and Measurements (NCRP) Report #144.     

In phantom Dosimetric response of tooth enamel to neutrons

Electron Paramagnetic Resonance dosimetry based on tooth enamel has one important application in dose reconstruction of nuclear plant workers, where the contribution of neutrons to individual dose is often important. Evaluation of tooth enamel response to neutrons is thus an important goal. A few experimental data at thermal and fast neutron energies are available. A first evaluation of the tooth enamel relative response to 60Co in monoenergetic neutron flux of 2.8 and of 14 MeV, published elsewhere, has provided results apparently non-consistent with the results obtained at lower and higher energies. A comparison of those results in the 2.8 and 14 MeV beams with those available in the literature for other beams is reported and possible reasons for incongruities are discussed. Dose conversion factors of enamel to the water and air are also calculated and reported.     

Investigation of radiation fields outside the Sub-critical Assembly in Dubna

The radiation fields outside the planned experimental Sub-critical Assembly in Dubna (SAD) have been studied in order to provide a basis for the design of the concrete shielding that cover the reactor core. The effective doses around the reactor, induced by leakage of neutrons and photons through the shielding, have been determined for a shielding thickness varying from 100 to 200 cm. It was shown that the neutron flux and the effective dose is higher above the shielding than at the side of it, owing to the higher fraction of high-energy spallation neutrons emitted in the direction of the incident beam protons. At the top, the effective dose was found to be -150 microSv s(-1) for a concrete thickness of 100 cm, while -2.5 microSv s(-1) for a concrete thickness of 200 cm. It was also shown that the high-energy neutrons (> 10 MeV), which are created in the proton-induced spallation interactions in the target, contribute for the major part of the effective doses outside the reactor.     

In-phantom spectra and dose distributions from a high-energy neutron therapy beam

In radiotherapy with external beams, healthy tissues surrounding the target volumes are inevitably irradiated. In the case of neutron therapy, the estimation of dose to the organs surrounding the target volume is particularly challenging, because of the varying contributions from primary and secondary neutrons and photons of different energies. The neutron doses to tissues surrounding the target volume at the Louvain-la-Neuve (LLN) facility were investigated in this work. At LLN, primary neutrons have a broad spectrum with a mean energy of about 30 MeV. The transport of a 10×10 cm² beam through a water phantom was simulated by means of the Monte Carlo code MCNPX. Distributions of energy-differential values of neutron fluence, kerma and kerma equivalent were estimated at different locations in a water phantom. The evolution of neutron dose and dose equivalent inside the phantom was deduced. Measurements of absorbed dose and of dose equivalent were then carried out in a water phantom using an ionization chamber and superheated drop detectors (SDDs). On the beam axis, the calculations agreed well with the ionization chamber data, but disagreed significantly from the SDD data due to the detector's under-response to neutrons above 20 MeV. Off the beam axis, the calculated absorbed doses were significantly lower than the ionization chamber readings, since gamma fields were not accounted for. The calculated data are doses from neutron-induced charge particles, and these agreed with the values measured by the photon-insensitive SDDs. When exposed to the degraded spectra off the beam axis, the SDD offered reliable estimates of the neutron dose equivalent.     

Calculations of neutron shielding data for 10-100 MeV proton accelerators

The characteristics of neutron sources and their attenuation in concrete were investigated in detail for protons with energies ranging from 10 to 100 MeV striking on target materials of C, N, Al, Fe, Cu and W. A two-step approach was adopted: thick-target double-differential neutron yields were first calculated from the (p, xn) cross sections recommended in the ICRU Report 63; further, transport simulations of those neutrons in concrete were performed by using the FLUKA Monte Carlo code. The purpose of this study is to provide reasonably accurate parameters for shielding design for 10-100 MeV proton accelerators. Source terms and the corresponding attenuation lengths in concrete for several target materials are given as a function of proton energies and neutron emission angles.     

Shielding variation effects for 250 MeV protons on tissue targets

This paper provides results of computer simulation studies with the goal to analyse issues regarding radiation protection for personnel, patients and third persons involved in hadron therapy treatment. The treatment room and the patient are modelled by simple cylindrical geometries at incident proton energies of 250 MeV. Monte Carlo simulations of the energy and angular dependence of proton, neutron and photon radiation fields and resulting ambient dose equivalent distributions outside the shielding walls are performed. In order to investigate systematic uncertainties due to the shielding materials and inherent to the computer models, various concrete compositions, densities and water contents are modelled, and the influence of simulation parameters on the results obtained is determined. Generally, good agreement is found between results provided by MCNPX and FLUKA computer codes. Variations in neutron ambient dose attenuation from -50 to +/-30% are found due to varying concrete composition. Changes in the water content of the concrete in the order of 8% may cause variations up to 20%.     

Dosimetric measurements with a brain equivalent plastic walled ionization chamber in an epithermal neutron beam

The tissue substitute A-181 plastic, which has an elemental composition matching both the constituent hydrogen and nitrogen of brain tissue, was assessed for dosimetry in boron neutron capture therapy (BNCT). The sensitivity of an A-181 walled ionization chamber relative to photons for all neutrons in a clinical epithermal beam was calculated to vary between 0.79 +/- 0.04 in-air and 0.95 +/- 0.01 at depths of 4 cm and greater in-phantom. Differences in the total neutron doses measured with A-150 and A-181 plastic-walled chambers were attributed, within experimental error, to the dose produced by thermal neutron capture reactions from the different concentrations of nitrogen in the two tissue substitutes. The response of the A-181 chamber was converted to total neutron dose with an uncertainty increasing with depth in-phantom from 13 to 23% the magnitude of which is determined by the subtraction of a relatively large photon dose. The use of A-181 in place of A-150 plastic will no longer require partitioning the measured neutron dose by energy and should simplify dose reporting in BNCT.     

Differential absorbed dose distributions in lineal energy for neutrons and gamma rays at the mono-energetic neutron calibration facility

Absorbed dose distributions in lineal energy for neutrons and gamma rays of mono-energetic neutron sources from 140 keV to 15 MeV were measured in the Fast Neutron Laboratory at Tohoku University. By using both a tissue-equivalent plastic walled counter and a graphite-walled low-pressure proportional counter, absorbed dose distributions in lineal energy for neutrons were obtained separately from those for gamma rays. This method needs no knowledge of energy spectra and dose distributions for gamma rays. The gamma-ray contribution in this neutron calibration field >1 MeV neutron was <3%, while for <550 keV it was >40%. The measured neutron absolute absorbed doses per unit neutron fluence agreed with the LA150 evaluated kerma factors. By using this method, absorbed dose distributions in lineal energy for neutrons and gamma rays in an unknown neutron field can be obtained separately.     

A study of neutron radiation quality with a tissue-equivalent proportional counter for a low-energy accelerator-based in vivo neutron activation facility

The accelerator-based in vivo neutron activation facility at McMaster University has been used successfully for the measurement of several minor and trace elements in human hand bones due to their importance to health. Most of these in vivo measurements have been conducted at a proton beam energy (E(p)) of 2.00 MeV to optimise the activation of the selected element of interest with an effective dose of the same order as that received in chest X rays. However, measurement of other elements at the same facility requires beam energies other than 2.00 MeV. The range of energy of neutrons produced at these proton beam energies comes under the region where tissue-equivalent proportional counters (TEPCs) are known to experience difficulty in assessing the quality factor and dose equivalent. In this study, the response of TEPCs was investigated to determine the quality factor of neutron fields generated via the (7)Li(p, n)(7)Be reaction as a function of E(p) in the range 1.884-2.56 MeV at the position of hand irradiation in the facility. An interesting trend has been observed in the quality factor based on ICRP 60, Q(ICRP60), such that the maximum value was observed at E(p)=1.884 MeV (E(n)=33±16 keV) and then continued to decline with increasing E(p) until achieving a minimum value at E(p)=2.0 MeV despite a continuous increase in the mean neutron energy with E(p). This observation is contrary to what has been observed with direct fast neutrons where the quality factor was found to increase continuously with an increase in E(p) (i.e. increasing E(n)). The series of measurements conducted with thermal and fast neutron fields demonstrate that the (14)N(n, p)(14)C produced 580 keV protons in the detector play an important role in the response of the counter under 2.0 MeV proton energy (E(n) ≤ 250 keV). In contrast to the lower response of TEPCs to low-energy neutrons, the quality factor is overestimated in the range 1-2 depending on beam energy <2.0 MeV. This study provides an insight to understanding the response of TEPCs in low-energy neutron fields where the neutrons are moderated using a polyethylene moderator.     

Experimental and theoretical study of the neutron dose produced by carbon ion therapy beams

High-energy (12)C ions offer favourable conditions for the treatment of deep-seated local tumours. Several facilities for the heavy ion therapy are planned or under construction, for example the new clinical ion-therapy unit HIT at the Radiological University Clinics in Heidelberg. In order to improve existing treatment planning models, it is essential to evaluate the secondary fragment production and to include these contributions to the therapy dose with higher accuracy. Secondary neutrons are most abundantly produced in the reactions between (12)C beams and tissues. The dose contribution to tissues by a neutron is fairly small compared with the projectile and the other charged fragments due to no ionisation and the small reaction cross-sections; however, it distributes in a considerably wider region beyond the bragg-peak because of the strong penetrability. Systematic data on energy spectra and doses of secondary neutrons produced by (12)C beams using water targets of different thicknesses for various detection angles have therefore been measured in this study at GSI Darmstadt.     

Measurements of the response functions of a large size NE213 organic liquid scintillator for neutrons up to 800 MeV

The response functions of 25.4 cm (length) x 25.4 cm (diameter) NE213 organic liquid scintillator have been measured for neutrons in the energy range from 20 to 800 MeV at the Heavy-Ion Medical Accelerator in Chiba (HIMAC) and at the Research Center for Nuclear Physics (RCNP) of Osaka University. At HIMAC, white (continuous) energy spectrum neutrons were produced by the 400 MeV per nucleon carbon ion bombardment on a thick graphite target, whose energy spectrum has already been measured by Kurosawa et al., [Nucl. Sci. Eng. 132, 30 (1999)] and the response functions of the time-of-flight-gated monoenergetic neutrons in a wide energy range from 20 to 800 MeV were simultaneously measured. At RCNP, the quasi-monoenergetic neutrons were produced via 7Li(p,n)7Be reaction by 250 MeV proton beam bombardment on a thin 7Li target, and the TOF-gated 245 MeV peak neutrons were measured. The absolute peak neutron yield was obtained by the measurement of 478 keV gamma rays from the 7Be nuclei produced in a Li target. The measured results show that the response functions for monoenergetic neutrons < 250 MeV have a recoil proton plateau and an edge around the maximum light output, which increases with increasing incident neutron energy, on the other hand > 250 MeV, the plateau and the edge become unclear because the proton range becomes longer than the detector size and the escaping protons increase. It can be found that the efficiency of the 24.5 cm (diameter) x 25.4 cm (length) NE213 for the 250 MeV neutrons is -10 times larger than the 12.7 cm (length) x 12.7 cm (diameter) NE213, which is widely used as a neutron spectrometer.     

The neutron component of two high-energy photon reference fields

The 4.4 MeV photon reference field described in ISO 4037 is produced by the (12)C(p,p')(12)C (E(x) = 4.4389 MeV) reaction using a thick elemental carbon target and a proton beam with an energy of 5.7 MeV. The relative abundance of the isotope (13)C in elemental carbon is 1.10%. Therefore, the 4.4 MeV photon field is contaminated by neutrons produced by the (13)C(p,n) (13)N reaction (Q = -3.003 MeV). The ambient dose equivalent H*(10) produced by these neutrons is of the same order of magnitude as the ambient dose equivalent produced by the 4.4 MeV photons. For the calibration of dosemeters, especially those also sensitive to neutrons, the spectral fluence distribution of these neutrons has to be known in detail. On the other hand, a mixed photon/neutron field is very useful for the calibration of tissue-equivalent proportional counters (TEPC), if this field combines a high-linear energy transfer (LET) component produced by low-energy neutrons and a low-LET component resulting from photons with about the same ambient dose equivalent and energies up to 7 MeV. Such a mixed field was produced at the PTB accelerator facility using a thin CaF(2) + (nat)C target and a 5.7 MeV proton beam.     

Analysis of the neutron component at high altitude mountains using active and passive measurement devices

The European Council directive 96/29/Euratom requires dosimetric precautions if the effective dose exceeds 1 mSv/a. On an average, this value is exceeded by aircrew members. Roughly half of the radiation exposure at flight altitudes is caused by cosmic ray-induced neutrons. Active (6LiI(Eu)-scintillator) and passive (TLDs) Bonner sphere spectrometers were used to determine the neutron energy spectra atop Mt. Sonnblick (3105 m) and Mt. Kitzsteinhorn (3029 m). Further measurements in a mixed radiation field at CERN as well as in a proton beam of 62 MeV at Paul Scherrer Institute, Switzerland, confirmed that not only neutrons but also charged particles contribute to the readings of active detectors, whereas TLD-600 and TLD-700 in pair allow the determination of the thermal neutron flux. Unfolding of the detector data obtained atop both mountains shows two relative maxima around 1 MeV and 85 MeV, which have to be considered for the assessment of the biologically relevant dose equivalent. By convoluting the spectra with appropriate conversion functions the neutron dose equivalent rate was determined to be 150 +/- 15 nSv/h. The total dose equivalent rate determined by the HTR-method was 210 +/- 15 nSv/h. The results are in good agreement with LET-spectrometer and Sievert counter measurements carried out simultaneously.     

High-fidelity MCNP modeling of a D-T neutron generator for active interrogation of special nuclear material

Fast and robust methods for interrogation of special nuclear material (SNM) are of interest to many agencies and institutions in the United States. It is well known that passive interrogation methods are typically sufficient for plutonium identification because of a relatively high neutron production rate from ²⁴⁰Pu [1]. On the other hand, identification of shielded uranium requires active methods using neutron or photon sources [2]. Deuterium-deuterium (2.45 MeV) and deuterium-tritium (14.1 MeV) neutron-generator sources have been previously tested and proven to be relatively reliable instruments for active interrogation of nuclear materials [3] and [4]. In addition, the newest generators of this type are small enough for applications requiring portable interrogation systems.Active interrogation techniques using high-energy neutrons are being investigated as a method to detect hidden SNM in shielded containers [4] and [5]. Due to the thickness of some containers, penetrating radiation such as high-energy neutrons can provide a potential means of probing shielded SNM. In an effort to develop the capability to assess the signal seen from various forms of shielded nuclear materials, the University of Michigan Neutron Science Laboratory’s D-T neutron generator and its shielding were accurately modeled in MCNP. The generator, while operating at nominal power, produces approximately 1×10¹⁰ neutrons/s, a source intensity which requires a large amount of shielding to minimize the dose rates around the generator. For this reason, the existing shielding completely encompasses the generator and does not include beam ports. Therefore, several MCNP simulations were performed to estimate the yield of uncollided 14.1-MeV neutrons from the generator for active interrogation experiments. Beam port diameters of 5, 10, 15, 20, and 25 cm were modeled to assess the resulting neutron fluxes. The neutron flux outside the beam ports was estimated to be approximately 2×10⁴ n/cm² s.     

Arrangement of high-energy neutron irradiation field and shielding experiment using 4 m concrete at KENS

An irradiation field of high-energy neutrons produced in the forward direction from a thick tungsten target bombarded by 500 MeV protons was arranged at the KENS spallation neutron source facility. In this facility, shielding experiment was performed with an ordinary concrete shield of 4 m thickness assembled in the irradiation room, 2.5 m downstream from the target centre. Activation detectors of bismuth, aluminium, indium and gold were inserted into eight slots inside the shield and attenuations of neutron reaction rates were obtained by measurements of gamma-rays from the activation detectors. A MARS14 Monte Carlo simulation was also performed down to thermal energy, and comparisons between the calculations and measurements show agreements within a factor of 3. This neutron field is useful for studies of shielding, activation and radiation damage of materials for high-energy neutrons, and experimental data are useful to check the accuracies of the transmission and activation calculation codes.     

Microdosimetric characteristics of the clinical proton beams at JINR,Dubna

High-energy proton radiotherapy beams give rise to secondary heavy charged particles with elevated linear energy transfer (LET), which contribute to the dose in a patient. This contribution to the characteristics of radiotherapy proton beams was experimentally studied by means of a LET spectrometer based on a track detector. The spectrometer permits LET spectra to be established in the region above 10 keV.micron-1 in tissue. Sets of track detectors were exposed in the various depths of a phantom irradiated with protons of two energies, 150 and 205 MeV. It was observed that the contribution of particles with the values of LET mentioned increases with the depth, representing from about 2 (at the surface) up to few tens% close to Bragg peak region of the total dose. There, some of primary protons contribute also above 10 keV.micron-1. Using the 'biological weighted function' proposed, the clinical RBE was calculated, it could approach 1.3. This effect has to be taken into account during the clinical beam production and the radiotherapy.