Intercomparison of radiation protection devices in a high-energy stray neutron field. Part III: Instrument response

Joint EURADOS WG9-WG11 rem-counter intercomparison in a Mevion S250i proton therapy facility with Hyperscan pulsed synchrocyclotron

Objective Proton therapy is gaining popularity because of the improved dose delivery over conventional radiation therapy. The secondary dose to healthy tissues is dominated by secondary neutrons. Commercial rem-counters are valuable instruments for the on-line assessment of neutron ambient dose equivalent (H*(10)). In general, however, a priori knowledge of the type of facility and of the radiation field is required for the proper choice of any survey meter. The novel Mevion S250i Hyperscan synchrocyclotron mounts the accelerator directly on the gantry. It provides a scanned 227 MeV proton beam, delivered in pulses with a pulse width of 10 µs at 750 Hz frequency, which is afterwards degraded in energy by a range shifter modulator system. This environment is particularly challenging for commercial rem-counters; therefore, we tested the reliability of some of the most widespread rem-counters to understand their limits in the Mevion S250i stray neutron field. Approach This work, promoted by the European Radiation Dosimetry Group (EURADOS), describes a rem-counter intercomparison at the Maastro Proton Therapy centre in the Netherlands, which houses the novel Mevion S250i Hyperscan system. Several rem-counters were employed in the intercomparison (LUPIN, LINUS, WENDI-II, LB6411, NM2B-458, NM2B-495Pb), which included simulation of a patient treatment protocol employing a water tank phantom. The outcomes of the experiment were compared with models and data from the literature. Main results We found that only the LUPIN allowed for a correct assessment of H*(10) within a 20% uncertainty. All other rem-counters underestimated the reference H*(10) by factors from 2 to more than 10, depending on the detector model and on the neutron dose per pulse. In pulsed fields, the neutron dose per pulse is a fundamental parameter, while the average neutron dose rate is a secondary quantity. An average 150-200 µSv/GyRBE neutron H*(10) at various positions around the phantom and at distances between 186 cm and 300 cm from it was measured per unit therapeutic dose delivered to the target. Significance Our results are partially in line with results obtained at similar Mevion facilities employing passive energy modulation. Comparisons with facilities employing active energy modulation confirmed that the neutron H*(10) can increase up to more than a factor of 10 when passive energy modulation is employed. The challenging environment of the Mevion stray neutron field requires the use of specific rem-counters sensitive to high-energy neutrons (up to a few hundred MeV) and specifically designed to withstand pulsed neutron fields.

Quality assurance for the use of computational methods in dosimetry: activities of EURADOS Working Group 6 'Computational Dosimetry'

Working Group (WG) 6 'Computational Dosimetry' of the European Radiation Dosimetry Group promotes good practice in the application of computational methods for radiation dosimetry in radiation protection and the medical use of ionising radiation. Its cross-sectional activities within the association cover a large range of current topics in radiation dosimetry, including more fundamental studies of radiation effects in complex systems. In addition, WG 6 also performs scientific research and development as well as knowledge transfer activities, such as training courses. Monte Carlo techniques, including the use of anthropomorphic and other numerical phantoms based on voxelised geometrical models, play a strong part in the activities pursued in WG 6. However, other aspects and techniques, such as neutron spectra unfolding, have an important role as well. A number of intercomparison exercises have been carried out in the past to provide information on the accuracy with which computational methods are applied and whether best practice is being followed. Within the exercises that are still ongoing, the focus has changed towards assessing the uncertainty that can be achieved with these computational methods. Furthermore, the future strategy of WG 6 also includes an extension of the scope toward experimental benchmark activities and evaluation of cross-sections and algorithms, with the vision of establishing a gold standard for Monte Carlo methods used in medical and radiobiological applications.

Ambient neutron and photon dose equivalent H*(10) around a pencil beam scanning proton therapy facility

OBJECTIVES

To measure leakage ambient dose equivalent H*(10) from stray secondary neutron and photon radiation around proton therapy (PT) facility and evaluate adequacy of shielding design.

METHODS AND MATERIALS

H*(10) measurement were carried out at 149 locations around cyclotron vault (CV), beam transport system (BTS) and first treatment room (GTR3) of a multiroom PT facility using WENDI-II and SmartIon survey meter. Measurement were performed under extreme case scenarios wherein maximum secondary neutrons and photons were produced around CV, BTS and GTR3 by stopping 230MeV proton of 300nA on beam degrader, end of BTS and isocenter of GTR3. Weekly time average dose rate (TADR) were calculated from H*(10) value measured at selective hot spots by irradiating actual treatment plans of mix clinical sites.

RESULTS

The maximum total H*(10) were within 2 µSv/hr around CV, 5 µSv/hr around outer wall of BTS which increases up to 62 µSv/hr at the end of inside BTS corridor. Maximum H*(10) of 20.8 µSv/hr in treatment control console (P125), 23.4 µSv/hr behind the common wall between GTR3 and GTR2 (P132) and 25.7 µSv/hr above isocenter (P99) were observed around GTR3. Reduction of beam current from 6 to 3 nA and 1 nA at nozzle exit lead to decrease in total H*(10) at P125 from 20.8 to 11.35 and 4.62 µSv/hr. In comparison to extreme case scenario, H*(10) value at P125, P132 and P99 from clinically relevant irradiation parameters were reduce by a factor ranging from 8.6 for high range cube to 46.4 for brain clinical plan. The maximum weekly TADR per fraction was highest for large volume, sacral chordoma patient at 8.5 µSv/hr compare to 0.3 µSv/hr for brain patient. The calculated weekly TADR for 30 mix clinical cases and 15 fractions of 1 L cube resulted total weekly TADR of 83-84 µSv/hr at P125, P132 and P99. The maximum annual dose level at these hot spots were estimated at 4.37 mSv/Yr.

CONCLUSION

We have carried out an extensive measurement of H*(10) under different conditions. The shielding thickness of our PT facility is adequate to limit the dose to occupational worker and general public within the permissible stipulated limit. The data reported here can bridge the knowledge gap in ambient dose around PT facility and can also be used as a reference for any new and existing proton facility for intercomparison and validation.

ADVANCES IN KNOWLEDGE

First extensive investigation of neutron and photon H*(10) around PT facility and can bridge the knowledge gap on ambient dose.

The work programme of EURADOS on internal and external dosimetry

Since the early 1980s, the European Radiation Dosimetry Group (EURADOS) has been maintaining a network of institutions interested in the dosimetry of ionising radiation. As of 2017, this network includes more than 70 institutions (research centres, dosimetry services, university institutes, etc.), and the EURADOS database lists more than 500 scientists who contribute to the EURADOS mission, which is to promote research and technical development in dosimetry and its implementation into practice, and to contribute to harmonisation of dosimetry in Europe and its conformance with international practices. The EURADOS working programme is organised into eight working groups dealing with environmental, computational, internal, and retrospective dosimetry; dosimetry in medical imaging; dosimetry in radiotherapy; dosimetry in high-energy radiation fields; and harmonisation of individual monitoring. Results are published as freely available EURADOS reports and in the peer-reviewed scientific literature. Moreover, EURADOS organises winter schools and training courses on various aspects relevant for radiation dosimetry, and formulates the strategic research needs in dosimetry important for Europe. This paper gives an overview on the most important EURADOS activities. More details can be found at www.eurados.org .

RECOMBINATION CHAMBERS-DO THE OLD IDEAS REMAIN USEFUL?

Recombination methods were invented almost contemporaneously with the introduction of the ICRU/ICRP concept of a radiation quality factor expressed in terms of linear energy transfer. Since then, recombination chambers and related methodologies have been extensively developed, along with new concepts of dosimetry and radiation protection. This article briefly presents the main features of recombination chambers and their possible applications in practical systems for measuring dosimetric quantities.

Radiation Protection Studies for Medical Particle Accelerators using Fluka Monte Carlo Code

Radiation protection (RP) in the use of medical cyclotrons involves many aspects both in the routine use and for the decommissioning of a site. Guidelines for site planning and installation, as well as for RP assessment, are given in international documents; however, the latter typically offer analytic methods of calculation of shielding and materials activation, in approximate or idealised geometry set-ups. The availability of Monte Carlo (MC) codes with accurate up-to-date libraries for transport and interaction of neutrons and charged particles at energies below 250 MeV, together with the continuously increasing power of modern computers, makes the systematic use of simulations with realistic geometries possible, yielding equipment and site-specific evaluation of the source terms, shielding requirements and all quantities relevant to RP at the same time. In this work, the well-known FLUKA MC code was used to simulate different aspects of RP in the use of biomedical accelerators, particularly for the production of medical radioisotopes. In the context of the Young Professionals Award, held at the IRPA 14 conference, only a part of the complete work is presented. In particular, the simulation of the GE PETtrace cyclotron (16.5 MeV) installed at S. Orsola-Malpighi University Hospital evaluated the effective dose distribution around the equipment; the effective number of neutrons produced per incident proton and their spectral distribution; the activation of the structure of the cyclotron and the vault walls; the activation of the ambient air, in particular the production of 41Ar. The simulations were validated, in terms of physical and transport parameters to be used at the energy range of interest, through an extensive measurement campaign of the neutron environmental dose equivalent using a rem-counter and TLD dosemeters. The validated model was then used in the design and the licensing request of a new Positron Emission Tomography facility.

Assessment of the neutron dose field around a biomedical cyclotron: FLUKA simulation and experimental measurements

In the planning of a new cyclotron facility, an accurate knowledge of the radiation field around the accelerator is fundamental for the design of shielding, the protection of workers, the general public and the environment. Monte Carlo simulations can be very useful in this process, and their use is constantly increasing. However, few data have been published so far as regards the proper validation of Monte Carlo simulation against experimental measurements, particularly in the energy range of biomedical cyclotrons. In this work a detailed model of an existing installation of a GE PETtrace 16.5MeV cyclotron was developed using FLUKA. An extensive measurement campaign of the neutron ambient dose equivalent H^∗(10) in marked positions around the cyclotron was conducted using a neutron rem-counter probe and CR39 neutron detectors. Data from a previous measurement campaign performed by our group using TLDs were also re-evaluated. The FLUKA model was then validated by comparing the results of high-statistics simulations with experimental data. In 10 out of 12 measurement locations, FLUKA simulations were in agreement within uncertainties with all the three different sets of experimental data; in the remaining 2 positions, the agreement was with 2/3 of the measurements. Our work allows to quantitatively validate our FLUKA simulation setup and confirms that Monte Carlo technique can produce accurate results in the energy range of biomedical cyclotrons.

Overview of the Radiation Dosimetry Experiment (RaD-X) flight mission

The NASA Radiation Dosimetry Experiment (RaD-X) stratospheric balloon flight mission addresses the need to reduce the uncertainty in predicting human exposure to cosmic radiation in the aircraft environment. Measurements were taken that characterize the dosimetric properties of cosmic ray primaries, the ultimate source of aviation radiation exposure, and the cosmic ray secondary radiations that are produced and transported to aviation altitudes. In addition, radiation detectors were flown to assess their potential application to long-term, continuous monitoring of the aircraft radiation environment. RaD-X was successfully launched from Fort Sumner, New Mexico (34.5°N, 104.2°W), on 25 September 2015. Over 18 h of science data were obtained from a total of four different type dosimeters at altitudes above 20 km. The RaD-X flight mission was supported by laboratory radiation exposure testing of the balloon flight dosimeters and also by coordinated radiation measurements taken on ER-2 and commercial aircraft. This paper provides the science background and motivation for the RaD-X flight mission, a brief description of the balloon flight profile and the supporting aircraft flights, and a summary of the articles included in the RaD-X special collection and their contributions to the science goals of the RaD-X mission.

SECONDARY NEUTRON DOSES IN A PROTON THERAPY CENTRE

The formation of secondary high-energy neutrons in proton therapy can be a concern for radiation protection of staff. In this joint intercomparative study (CERN, SCK•CEN and IBA/IRISIB/ULB), secondary neutron doses were assessed with different detectors in several positions in the Proton Therapy Centre, Essen (Germany). The ambient dose equivalent H(*)(10) was assessed with Berthold LB 6411, WENDI-2, tissue-equivalent proportional counter (TEPC) and Bonner spheres (BS). The personal dose equivalent Hp(10) was measured with two types of active detectors and with bubble detectors. Using spectral and basic angular information, the reference Hp(10) was estimated. Results concerning staff exposure show H(*)(10) doses between 0.5 and 1 nSv/monitoring unit in a technical room. The LB 6411 showed an underestimation of H(*)(10), while WENDI-2 and TEPC showed good agreement with the BS data. A large overestimation for Hp(10) was observed for the active personal dosemeters, while the bubble detectors showed only a slight overestimation.

An analytical model of leakage neutron equivalent dose for passively-scattered proton radiotherapy and validation with measurements

Exposure to stray neutrons increases the risk of second cancer development after proton therapy. Previously reported analytical models of this exposure were difficult to configure and had not been investigated below 100 MeV proton energy. The purposes of this study were to test an analytical model of neutron equivalent dose per therapeutic absorbed dose (H/D) at 75 MeV and to improve the model by reducing the number of configuration parameters and making it continuous in proton energy from 100 to 250 MeV. To develop the analytical model, we used previously published H/D values in water from Monte Carlo simulations of a general-purpose beamline for proton energies from 100 to 250 MeV. We also configured and tested the model on in-air neutron equivalent doses measured for a 75 MeV ocular beamline. Predicted H/D values from the analytical model and Monte Carlo agreed well from 100 to 250 MeV (10% average difference). Predicted H/D values from the analytical model also agreed well with measurements at 75 MeV (15% average difference). The results indicate that analytical models can give fast, reliable calculations of neutron exposure after proton therapy. This ability is absent in treatment planning systems but vital to second cancer risk estimation.

Configuration and validation of an analytical model predicting secondary neutron radiation in proton therapy using Monte Carlo simulations and experimental measurements

PURPOSE

This study focuses on the configuration and validation of an analytical model predicting leakage neutron doses in proton therapy.

METHODS

Using Monte Carlo (MC) calculations, a facility-specific analytical model was built to reproduce out-of-field neutron doses while separately accounting for the contribution of intra-nuclear cascade, evaporation, epithermal and thermal neutrons. This model was first trained to reproduce in-water neutron absorbed doses and in-air neutron ambient dose equivalents, H*(10), calculated using MCNPX. Its capacity in predicting out-of-field doses at any position not involved in the training phase was also checked. The model was next expanded to enable a full 3D mapping of H*(10) inside the treatment room, tested in a clinically relevant configuration and finally consolidated with experimental measurements.

RESULTS

Following the literature approach, the work first proved that it is possible to build a facility-specific analytical model that efficiently reproduces in-water neutron doses and in-air H*(10) values with a maximum difference less than 25%. In addition, the analytical model succeeded in predicting out-of-field neutron doses in the lateral and vertical direction. Testing the analytical model in clinical configurations proved the need to separate the contribution of internal and external neutrons. The impact of modulation width on stray neutrons was found to be easily adjustable while beam collimation remains a challenging issue. Finally, the model performance agreed with experimental measurements with satisfactory results considering measurement and simulation uncertainties.

CONCLUSION

Analytical models represent a promising solution that substitutes for time-consuming MC calculations when assessing doses to healthy organs.

Comparing Geant4 hadronic models for the WENDI-II rem meter response function

The WENDI-II rem meter is one of the most popular neutron dosemeters used to assess a useful quantity of radiation protection, namely the ambient dose equivalent. This is due to its high sensitivity and its energy response that approximately follows the conversion function between neutron fluence and ambient dose equivalent in the range of thermal to 5 GeV. The simulation of the WENDI-II response function with the Geant4 toolkit is then perfectly suited to compare low- and high-energy hadronic models provided by this Monte Carlo code. The results showed that the thermal treatment of hydrogen in polyethylene for neutron <4 eV has a great influence over the whole detector range. Above 19 MeV, both Bertini Cascade and Binary Cascade models show a good correlation with the results found in the literature, while low-energy parameterised models are not suitable for this application.

Neutron spectrometry with Bonner Spheres for area monitoring in particle accelerators

Selecting the instruments to determine the operational quantities in the neutron fields produced by particle accelerators involves a combination of aspects, which is peculiar to these environments: the energy distribution of the neutron field, the continuous or pulsed time structure of the beam, the presence of other radiations to which the neutron instruments could have significant response and the large variability in the dose rate, which can be observed when moving from areas near the beam line to free-access areas. The use of spectrometric techniques in support of traditional instruments is highly recommended to improve the accuracy of dosimetric evaluations. The multi-sphere or Bonner Sphere Spectrometer (BSS) is certainly the most used device, due to characteristics such as the wide energy range, large variety of active and passive detectors suited for different workplaces, good photon discrimination and the simple signal management. Disadvantages are the poor energy resolution, weight and need to sequentially irradiate the spheres, leading to usually long measurement sessions. Moreover, complex unfolding analyses are needed to obtain the neutron spectra. This work is an overview of the BSS for area monitoring in particle accelerators.

Workplace monitoring of mixed neutron-photon radiation fields and its contribution to external dosimetry

Workplace monitoring is a common procedure for determining measures for routine radiation protection in a particular working environment. For mixed radiation fields consisting of neutrons and photons, it is of increased importance because it contributes to the improved accuracy of individual monitoring. An example is the determination of field-specific correction factors, which can be applied to the readings of personal dosemeters. This paper explains the general problems associated with individual dosimetry of neutron radiation, and describes the various options for workplace monitoring. These options cover a range from the elaborate field characterisation using transport calculations or spectrometers to the simpler approach using area monitors. Examples are given for workplaces in nuclear industry, at particle accelerators and at flight altitudes.