Introduction
Proton therapy shows the ability of reducing normal tissue doses while maintaining adequate dose coverage to tumour for many different sites. Reference Boehling, Grosshans and Bluett1–Reference Yock, Schneider, Friedmann and Adams4 Two beam delivery techniques were generally used in proton therapy, which include passively scattered proton therapy (PSPT) and scanning beam proton therapy (SBPT). 5 The PSPT utilises a single or double scattering foil/foils to spread the beam and achieve lateral beam coverage, while SBPT such as uniform scanning and pencil beam scanning (PBS) proton therapy involves directly moving a pencil shape-charged particle beam (here in this case proton) while depositing the dose throughout the target volume. The MEVION S250i with HYPERSCAN and Adaptive Aperture™ (Mevion Medical Systems Inc., Littleton, MA, USA) is a single room compact PBS system, recently installed and commissioned only in few institutions. This compact design utilising a gantry-mounted superconducting synchrocyclotron eliminates the need for complex beam transport system and provides lower cost option for implementing proton therapy. Comparing with the PSPT nozzle, the PBS nozzle does not need any scattering foils in the beamline. Patient-specific apertures were also eliminated and a new Adaptive Aperture (AA) system was introduced to collimate the peripheral proton spots. Patient compensator was also replaced by the range shifter plates to modify the beam energy and depths. A TOPAS Monte Carlo simulation of such system shown detailed system schematic in Figure 1. Reference Chiang, Bunker, Jin, Ahmad and Chen6
TOPAS Monte Carlo simulation schematic of Mevion S250i HYPERSCAN nozzle.
The production of radioactive materials plays an important role in radiation safety for proton therapy facilities. Reference Moritz7 A quantitative assessment of the production of long-live radioactive materials is also required for state licensing and radiation safety assessment. In our case, radioactive materials are mainly produced as a result of the operation of the synchrocyclotron accelerator. Cyclotron components, beam transport line and other room structure elements become radioactive due to the direct irradiations with proton beam as well as indirect irradiations with secondary particles. The main objective of the present study has been to provide a quantitatively assessment of the production of radioactive materials during a year-long clinical operation of a Mevion S250i Hyperscan proton therapy system. To the best of our knowledge, this report is the first of its kind to address residual radioactivity from a clinical-operated particle therapy system.
Materials and Methods
Radioactive material production
The possible source terms for the Mevion S250i Hyperscan proton system may include the following: (1) activation of cyclotron due to internal beam loss, (2) activation of collimating components after beam exit, (3) activation of carbon used as absorber material after beam exit, (4) activation of aluminium used in the frames and supports near the beamline, (5) activation of accelerator and vault by thermal neutron capture secondary to proton beam irradiation and (6) activation of air and water in the vault by high-energy protons. The activation of cyclotron will be consistent with those produced in high-energy accelerator when high-energy proton interactions on steel and brass components. The expected isotopes are principally including Be-7, C-11, F-18, Na-22, Na-24, Sc-46, Sc-48, V-48, Cr-51, Mn-52, Mn-54, Co-57, Co-60 and Zn-65. After the proton beam entering nozzle, the activation of collimating components might also happen when the proton is stopping by the AA, which is made of brass. The possible isotopes would be similar to those activations at cyclotron as listed above. The proton beam will also interact with carbon-based absorber material in range shifter plates and produce radioactive isotope. The products of proton interactions on carbon could be H-3, Be-7, C10 and C-11. Tritium, with half-life of 12·3 years, will not reach quantities in the μCi range.
On the other hand, the principal materials, such as copper and manganese, near an accelerator with significant thermal neutron activation cross-section will interact with secondary neutron too. The principal produced isotopes include Cu-64 and Mn-56. Eu-152 can also be produced by the neutron capture reaction in trace amounts of stable Europium of concrete. Reference Carroll8 The aluminium used in frame and supports near the beamline might also become activated by spallation reaction with secondary neutrons above 20 MeV. Activation products of spallation reaction are limited to F-18, Na-24 and Na-22. For last possible source, the radioactive isotopes with short half-life such as O-14, O-15, N-13 and C-11 will also be produced when high-energy beam passes through air or water.
Among all possible radioactive nuclides, Na-22 from spallation reaction and neutron captured produced Eu-152 with long half-life of 2·6 and 13·54 years were assumed to be the main contributor in this assessment; other possible radionuclide can also be found but Eu-152 and Na-22 are the most important. Reference Tesse, Boogert and Gnacadja9 Therefore, in this technical note, the accumulated activity of Na-22 and Eu-152 will be calculated based on the measured dose level at isocentre.
Year-long weekly room survey
To qualitatively assess the production of radioactive materials during the normal condition of clinical operation, we have conducted a year-long weekly room survey starting from 28 January 2019. The residual radioactive level was monitored inside the treatment room on every Sunday morning after roughly 40 hours of normal daily operation of patient treatments which ended around 5 PM Friday afternoon. Three locations of measurement that included room corner, room isocentre and nozzle surface are shown in Figure 2. A hand-held gamma spectrometer, identiFINDER-U (Thermo Scientific, Waltham, MA, USA) was used for the room survey for the production of accumulated radioactive materials during normal clinical operation.
Scheme of measurement points at treatment room of MEVION S250i proton therapy system.
Relationship between radionuclide gamma emission and exposure rate
The relationship between exposure (X) and energy fluence (Ψ) can be obtained by the following equation:
$$X = {{\Psi}}{\left( {{{{\mu _{en}}} \over \rho }} \right)_{air}}{\left( {{e \over {\overline W}}} \right)_{air}}$$
where
${\left( {{{{\mu _{en}}} \over \rho }} \right)_{air}}$
and
${\left( {{{\overline W} \over e}} \right)_{air}}$
are the mass energy absorption coefficients and the average energy required per unit charge of ionisation produced in air, respectively.
The energy fluence rate (Ψ) from a radioactive source can be determined by the following equation:
$${\rm{{\vskip -2.5pt\dot}\hskip-5pt\rm\Psi }} = {{AyE} \over {4\pi {r^2}}}$$
where ‘A’ represents the activity of the source in becquerels (Bq), ‘y’ and ‘E’ represent the yield and energy of photons, respectively and ‘r’ is the distance in cm.
Using the above two equations, we then obtained the relationships between the exposure rate (
$${\vskip -1.7pt\dot}\hskip-5pt{X}$$
) and radionuclide gamma emission by the following equation:
$${\vskip -1.7pt\dot}\hskip-5pt{X} = {{AyE} \over {4\pi {r^2}}}{\left( {{{{\mu _{en}}} \over \rho }} \right)_{air}}{\left( {{e \over {\overline W}}} \right)_{air}}$$
Using an average of 33·8 eV per ion pair, combining all the numerical terms, and considering gamma ray in different energies and yields, we obtained the following:
$${\vskip -1.7pt\dot}\hskip-5.5pt{X} = 5.263 \times {10^{ - 6}}A\sum {{{y_i}{E_i}{{\left( {{{{\mu _{en}}} \over \rho }} \right)}_i}} \over {{r^2}}}$$
where A is in Bq, yi and Ei represent the yield and energy of photon at each energy i (MeV). The decay information of Na-22 and Eu-152 including energies, yields and mass attenuation coefficients is shown in Tables 1 and 2, respectively. 10,11
The decay information of Na-22 including energies, yields and mass attenuation coefficients
The decay information of Eu-152 including energies, yields and mass attenuation coefficients
Results and Discussion
The year-long measurements for the locations of room isocentre and the corner are shown in Figure 3. The room activation levels were shown to gradually increase in the first 9 months and slowly plateaued thereafter. A few outlier data point readings were highlighted as they were taken during insufficient machine cold down time due to weekend patient treatment and maintenance work and were removed from further analysis. From these data, we could safely assume an average activation level of 35 urem/hour at room isocentre with normal treatment operation.
Year-long room survey measurements at two different locations, the outlier data points were highlight and removed from the results.
The production of radioactive materials is mainly due to the neutron capture, which produces Eu-152, and spallation reaction, which generates Na-22. The typical spectrum and possible radionuclides measured from our weekly room survey at nozzle surface are shown in Figure 4.
Energy spectrum measured at nozzle surface.
If we assume that the measured dose level was contributed mainly by Na-22 from secondary neutron-activated structural aluminium at 3 m away from the isocentre in the beamline, then the accumulated activity of Na-22 was calculated to be 274·9 μCi. On the other hand, if we assume that the measured dose level was contributed mainly by Eu-152 from neutron capture of concrete at about 2 m away from isocentre, then the accumulated activity of Eu-152 was calculated to be 75·1 μCi. The details of our calculation are shown in the Appendix.
In this study, we conservatively assumed that the contribution of measured dose level was only from one source, namely Na-22 or Eu-152. The accumulated activity might be lower for individual radionuclide if contributions from all other possible radionuclides are considered. However, comparing two main contribution radionuclides, the accumulated activity should be less than 300 μCi for all possible radionuclide for the assessment purpose of the production of radioactive materials.
Conclusions
At a worst-case scenario, the accumulated activity on aluminium structures was found to be less than 300 μCi after 1 year of patient treatment. We do not anticipate that the accumulated radiation level would rise in the future under normal condition of operation with patient treatments. We have thus reported for the first time the residual radioactivity from a clinical operation of a compact proton therapy system. We do not anticipate the accumulated radiation level to rise in future, and safe clinical operation will thus continue.
Acknowledgement
None.
Appendix
The relationship between radionuclide gamma emission and exposure rate is already shown in equation 4.
The decay radiation data including gamma energies and yields of decays can be found from National Nuclear Data website established by Brookhaven National Laboratory (https://www.nndc.bnl.gov/nudat2/indx_dec.jsp). The mass energy absorption coefficients for the photon in air can be obtained by interpolation from the table of National Institution of Standard and Technology (https://physics.nist.gov/PhysRefData/XrayMassCoef/ComTab/air.html).
The energies and yields of four possible decay photons and their mass energy absorption coefficients from parent nucleus Na-22 were listed in Table 1. Therefore, if we assume that the measured dose level was contributed mainly by Na-22 at 3 m away, the accumulated activity of Na-22 can be calculated to be 274·9 μCi using equation 5 shown below.
$$A\left( {Ci} \right) = {{{\vskip -1.7pt\dot}\hskip-5.5pt{X}\left( {R/hr} \right)} \over {5.263 \times {{10}^{ - 6}}\sum \left( {{{{y_i}{E_i}{{\left( {{{{\mu _{en}}} \over \rho }} \right)}_i}} \over {{r^2}}}} \right)}} \times {1 \over {3.7 \times {{10}^{10}}}}\left( {{{Ci} \over {Bq}}} \right)$$
The energies (above 1 MeV) and yields of four possible decay photons and their mass energy absorption coefficients from parent nucleus Eu-152 were listed in Table 2. Again, if we assume that the measured dose level was contributed mainly by Eu-152 at 2 m away from isocentre in the beamline, the accumulated activity of Eu-152 can be calculated to be 75·1 μCi using equation 5 shown above.
