Abstract
Prompt gamma ray in proton therapy is the product of a nuclear reaction between a proton and a target. The characteristic energies and intensities of prompt gamma lines can be used to determine the types of elements and their amounts in the target. In several previous experiments, it was demonstrated that no matter how complex the reaction cross-section is, once the energy of the incident proton and the irradiated element are determined, there is a definite linear relationship between the element concentration and the number of gamma-ray photons. However, this linear relationship is difficult to apply to medical imaging, and the nonlinear behavior of hydrogen has not been investigated so far. In this study, this linear relationship is extended to mixed elemental materials including a nonlinear case such as hydrogen, and a universal mathematical form, which is referred to as the prompt gamma spectroscopy retrieval algorithm (PGSRA), is developed. The basic assumption of the PGSRA is that the PGS of the sample material has a relationship with the molar gamma lines of the elements. For carbon and oxygen, this relationship is linear, while for hydrogen, this relationship is nonlinear. As the 2.23 MeV gamma line originates from neutron absorption radiation, the behavior of hydrogen is carefully investigated. The linear and nonlinear relationships are verified using Monte Carlo simulations with different combinations of carbon, oxygen, and hydrogen, such as PMMA, pentanediol, and ethanediol. The PGSRA developed in this work could be the first bridge between PGS and medical imaging.
1 Introduction
Prompt gamma (PG) is being widely investigated in proton and carbon ion therapy because it is a promising solution for online range verification [1, 2]. Different kinds of PG cameras were developed for 1D range measurement [3–6] and 3D image measurement [7–9]. Some of these cameras have already been tested in clinical trials [10, 11]. New PG timing systems, such as fast detectors based on CeBr3 [12] and LaBr3 [13], have been developed in recent years. Zhang et al. developed a new multi-slit PG imaging system based on BGO + SiPM [14, 15].
Monte Carlo code, such as FLUKA [16–18] and Geant4 [19, 20], are powerful tools for particle research. Some new ideas or new designs can be checked by conducting simulations on MC platforms; for example, the time-of-flight (TOF) method was adopted to study PG imaging using Geant4 [21, 22], and a new pixelated PG detector configuration was investigated using FLUKA [23]. Different configurations, such as multi-slit and knife-edge configurations, were compared using Geant4 [24]. Changran et al. investigated the Doppler shift effect for range verification using Geant4 [25, 26]. Recently, a new method for spatiotemporal emission reconstruction PG timing (SER-PGT) was developed and verified by FLUKA [27].
Protons deposit energy in matter through interactions when protons get into tissues during proton therapy. Proton–nuclear interactions involve both elastic and inelastic processes, including nuclear capture and nuclear scattering. For many scattering processes, the tissue nuclei remain intact and are left in an excited energy state. The decay of these excited nuclei typically produces a gamma ray with an energy ranging from 0 MeV to 11 MeV within a few nanoseconds of the proton–nucleus interaction, known as PG emission. The energy spectra of this PG emission depend on the specific nuclear energy states of the excited elemental nuclei, resulting in each element producing a unique spectrum, known as PG spectroscopy. Thus, tissues composed of different elements and elemental concentrations will produce different emission spectra during irradiation [28]. At present, prompt gamma spectroscopy (PGS) can be used to measure the beam range with high resolution [29–31]. In addition, Joost M. Verburg et al. investigated gamma lines and developed a novel optimization algorithm to determine the absolute range with a mean statistical precision of 1.1 mm [5, 11].
Range assessment is not the only capability of PGS. Similar to prompt gamma-ray neutron activation analysis (PGNAA) [32–35], the characteristic energies and intensities of gamma lines can be used to determine the types of elements and their amounts [36, 37]. The major difference between PGS and PGNAA is the kind of particle that initiates the nuclear reaction in the material. Verburg et al. assessed the complex cross-section and total photon yield for different incident proton energies on different target materials [35]. Polf et al. characterized how PG emission from tissue changes as a function of carbon and oxygen concentrations [38, 39]. Paulo Martins et al. demonstrated a feasible technique for proton and ion beam spectroscopy (PIBS) that was capable of determining the elemental concentrations of irradiated tissues during particle therapy [40]. In the experiments of Polf et al. and Paulo Martins et al., a linear relationship was demonstrated between the concentration of irradiated oxygen in tissue-equivalent samples and the total emission of PG radiation during irradiation by fixed energy protons.
Their experimental work demonstrated that no matter how complex the reaction cross-section is, once the energy of the incident proton and the irradiated element are determined, there is a definite relationship between the element content and the number of gamma-ray photons. However, this linear relationship is difficult to apply to medical imaging. These works did not extend the mathematical form of the linear relationship to more than one element. Our goal is to develop an algorithm that could be used to determine the content of multiple elements in the target, not just the concentration, from the PGS data obtained during irradiation. Additionally, a discussion on how to calculate the hydrogen amount as intensity from the neutron (slow neutron) absorption reaction (Eq. 1), which is not the direct production of proton/carbon irradiation, was not presented in the aforementioned works. Thus, the intensity is not a linear function of the hydrogen amount [41]. Jeyasingam et al. also measured PG-ray emissions from elements found in tissue; however, they also did not investigate the gamma lines from the neutron absorption reaction [42].
In the present study, the nonlinear behavior of the intensity from different samples is first investigated. In addition, a retrieval algorithm is developed to measure the amounts of all elements in homogeneous materials, including PMMA, pentanediol, and ethanediol. Their chemical formulas and basic molecular structures are presented in Figure 1. As a feasibility study, all the irradiation and gamma-ray transport processes are simulated by using the FLUKA v4 which is currently distributed by the FLUKA CERN collaboration [43].
FIGURE 1
2 Materials and methods
2.1 Basic assumptions
From an ideal simulation shown in Figure 2, we used a proton beam to irradiate a PMMA molecule, which yields gamma-ray spectroscopy (left in blue). Similarly, we used a proton beam to irradiate one carbon nucleus and one oxygen nucleus, which yields gamma-ray spectroscopy (right in dark) and (right in red), respectively. Finally, we used a neutron to bombard a hydrogen nucleus to obtain gamma-ray spectroscopy (right in yellow). Here, is spectroscopy from and reactions. There is a relationship between and , , . In this section, we aim to formalize a mathematical form to describe this relationship.
FIGURE 2
Normally, the unit of is photon number per incident proton per molecule or per mol at different gamma lines. In this study, gamma-ray spectroscopy is the quantitative of the energy spectra of gamma ray from the interaction between the proton and the specific matter or element. Generally, the x-axis of a spectroscopy is the energy channel, while the y-axis is the gamma-ray count.
As experiments from Polf et al. and Paulo Martins et al. proved that there is a definite relationship between the element content and the number of gamma-ray photons, we formalize a mathematical model from a phase-only perspective which could be solved by the numerical method in later sections. Thus, we assume a simple relationship, as in Eq. 2.
The spectroscopy of a target molecule is a sum of the spectroscopy from elemental carbon, oxygen, and hydrogen with corresponding weight ratios: , , and . The dominant factor is the number of element atoms in the molecule. From a microscopic point of view, the more atoms of an element in a molecule, the higher the probability is of emitting its characteristic gamma rays. Weight ratios and represent the numbers of carbon and oxygen, respectively, in the chemical formula; here, for example, and . However, for hydrogen, the major gamma line at 2.23 MeV strongly relies on neutron yields based on the absorption reaction (Eq. 1). Thus, the spectroscopy of hydrogen depends not only on the number of hydrogen atoms but also on the number of neutrons that cause reactions, which in turn come from reactions in other proton–nuclear interactions such as and . Thus, the weight ratio of hydrogen spectroscopy should be related to , , and . Considering the different sources of neutrons in absorption reaction, we have Eqs 3–5:
The first two terms of Eq. 3 are the same as Eq. 2, as neutrons for the absorption reaction come from the interaction of and . The third term of Eq. 2 hydrogen spectroscopy should be split into two parts: spectroscopy originates from reactions where the neutrons come from reactions, and the spectroscopy originates from reactions where the neutrons come from reactions. Correspondingly, their weight ratios Eq. 4 and Eq. 5 should be a nonlinear function of , , and .
In Eq. 4, the dominant parameter is the number of hydrogen atoms in the molecule. The ratio represents the factor of the neutron, that is, the neutron absorption reaction (Eq. 1), where the neutron is emitted from the reaction. The ratio is a hydrogen enhancement factor because when there are more hydrogens in a molecule with the same number of carbon or oxygen, the neutron absorption probability will increase. In Eq. 5, the dominant parameter is the number of hydrogen atoms in the molecule, and the ratio represents the factor of the neutron emitted from the reaction. is a hydrogen enhancement factor same as in Eq. 4.
In Section 3, we will illustrate that for hydrogen, there is only one peak at 2.23 MeV in its PGS, and it does not contribute to the other characteristic gamma lines in the carbon and oxygen spectra. Thus, for an unknown material we have the following changes: , , in Eq. 3. And Eq. 3 can be mathematically separated into two equations (Eq. 6):
In this article, we use to represent PG spectroscopy in a given material or an element per molecule or per mol, and represents the count at a given energy of 2.23 MeV in the material per molecule or per mol. The lower index is the material, and the upper index is the energy. and are 2.23 MeV gamma counts originating from carbon and oxygen per molecule or per mol. Eq. 1 to Eq. 6 are formalized based on a single molecule or 1 mol. In the next section, we will derive these equations to a large volume target.
2.2 PGS retrieval algorithm
In a real experiment, one cannot shoot a single molecule target with a single proton. If the target phantom has mol molecules, then the elements’ molar amount is
The total PG count from a phantom iswhere is total spectroscopy and is molar spectroscopy. The upper index FullEnergy means the gamma lines at every energy in consideration. For 2.23 MeV, the total gamma count of hydrogen is
In fact, is one peak point count at 2.23 MeV in . After both sides of Eq. 6 are multiplied by , we obtainHere, and are the PGS per mol carbon and oxygen, respectively. and are the PG counts per mol hydrogen with neutrons originating from carbon and oxygen, respectively. These molar PGS values can be measured or calculated from simple substances, such as hydrocarbons and water.
After both sides of Eq. 4 and Eq. 5 are multiplied by , we obtain volume weight ratios Eq. 8 for and Eq. 9 for in Eq. 7b:
As shown in Figure 3, the dark curve is the MC simulated PGS from 150 MeV proton beam, which irradiate a PMMA cylinder, and the red points are the peaks in the PGS spectrum. For example, the count at 2 MeV results from the contribution of 4 reactions; the first three reactions come from oxygen and the last one comes from carbon [44]. The count at 4.44 MeV also results from a contribution of both carbon and oxygen. The count at 5.25 MeV only comes from oxygen.
FIGURE 3
There are only two unknown variables and in Eq. 7a. Thus, for a binary linear equation, we need at least two peak points to solve the equation. From Figure 3, we always obtain more than 3 points from a PGS spectrum, and Eq. 7a can be solved by using the least squares method. These points are sorted into three different groups to compare the influence of the selection of different energy peaks on the results:
Low: low-energy peaks from 1.5 to 4 MeV.
High: high-energy peaks from 4 to 7 MeV.
Medium-high: peaks from medium energy 2 MeV to high energy 7 MeV.
To obtain a numerical solution, we discretize Eq. 7a:where matrix is the measured molar PGS point of carbon and oxygen, and is the measured PGS point of the sample material. Index is the selected gamma energy point of PGS from the target phantom. Then, the unknown element molar number matrix can be determined by the least square method.
The last step is to solve Eq. 7b and bring coefficients and into Eq. 7b; we have Eq. 12:
Now, only is an unknown variable, and the analytical solution of iswhere is the peak point count at 2.23 MeV for target material , which is determined by neutron distribution (energy, time, and lateral distribution) inside the phantom. and are the PG counts per mol hydrogen with neutrons originating from carbon and oxygen. If , and are properly measured or calculated, we can get hydrogen molar amount .
Finally, the target elements’ molar amounts , and are all retrieved. We refer to this method as the prompt gamma spectroscopy retrieval algorithm (PGSRA). In addition, PG can be used to measure the beam range in a target by using the TOF method [6] or a multi-slit camera [23, 45, 46]. Thus, the target volume and density can be measured.where is the lateral surface area and is the proton range in the target. Density can be determined by Eq. 15:
This mathematical algorithm is not limited to hydrogen, oxygen, and carbon. It can be extended to sulfur, nitrogen, or any other element. As a feasibility study, the focus of this study is hydrogen, oxygen, and carbon. We attempt to verify this algorithm on three molecules: PMMA (), pentanediol (), and ethanediol (). The purpose is to retrieve the molar numbers of the elements and calculate the density of these samples.
2.3 Cloud computing
To improve the statistics, the run primary number is for every MC simulation. The typical time consumption for a primary particle number of is almost 2 h on a standalone computer (Intel Core (TM) i7-8550U CPU @ 1.8 GHz and 16 GB memory). It takes 100 h, which is approximately 1 week, to complete one simulation. Thus, we have to use a parallel method to reduce the time consumption [47]. It is logical to break down a large primary number into several smaller ones with valid random numbers. The Python script is Fluka_p_run.py. Every simulation takes several hours; therefore, we develop a special Python script (Status.py) that terminates cloud computing once the simulation is finished. Another function of the script is to monitor statuses such as the number of processes that are still running. This script is cost-effective and time-saving. Another useful script is KillAll.py, which is used to terminate all the jobs on the cloud if a problem is discovered. In this simulation, we use Amazon Web Service’s HPC cluster instead of the traditional physical cluster.
2.4 Neutron distribution
As Eq. 13 described, neutron distribution is the key to the nonlinear solution of PGSRA. First, we must figure out where the neutrons are generated and where the capture reactions take place. We do not consider the hydrogen atoms that are far away from the beam path in lateral direction as they can reduce the accuracy of the algorithm, and nowadays pencil beams are adopted in most proton therapy centers. In practice, we are not likely to put a detector in the target to measure the locations of the reaction. We can roughly determine the location of the neutron by its time of flight from inside the target to detector outside. Thus, second, we calculate the time-of-flight for the neutron in different lateral layers in the target phantom. Based on data from the EXFOR library [48], slow neutrons (energy below ∼10 keV) have a higher cross-section (∼1 mb) for absorption reactions (Eq. 1). Thus, third, we need to determine the neutron distribution in energy. By the way, we also calculate the TOF of gamma photons in different layers in the same setup as neutrons.
A water phantom with 7 layers ranging from a radii of 2 cm to 14 cm was irradiated by 150 MeV protons, as shown in Figure 4. The USRYIELD card was adopted to score the neutron yield at different energies and time of flight for each radius. In this section, water is chosen as the material for calculation because the human body is mostly composed of water, which can be used as a rough estimation of the neutron distribution. All the calculation results in this section will finally give the recommended radius value of the phantom, which will be used in all the following simulations.
FIGURE 4
For the R = 2 cm layer, we set two slow neutron energy groups: 0–1 eV and 0–1 keV and one fast neutron group: 1 keV–150 MeV. From Figure 5, we can see that high-energy neutrons are generated immediately after the protons reach the target, and the slow neutron (0–1 keV group) yield increases rapidly and peaks at approximately 50 ns. The amount of 0–1 eV neutrons is much smaller than that from 0 to 1 keV neutrons before 200 ns.
FIGURE 5
The neutron (0–100 keV group) time-of-flight (TOF) distribution from R = 2 cm to R = 14 cm is shown in Figure 6. The neutron yield increases gradually and peaks at approximately 5 ns when R = 2 cm. As neutrons travel to the outside layers, the peak time increases and the neutron yield decreases. Thus, we choose R = 2 cm as the target radius for the absorption reaction region.
FIGURE 6
Correspondingly, we score the TOF of gamma photons (0–10 MeV group) in each layer from R = 2 cm to R = 14 cm. As shown in Figure 7, most photons are generated in the center layer R = 2 cm from 0.1 ns to 1.4 ns. The outside layers will start generating photons slightly later, within 0.6 ns, as neutron and other secondary particles require a short time to travel in the lateral direction in the water phantom. Thus, we will set an R = 2 cm slender cylindrical target in the feasibility study in the following simulation work.
FIGURE 7
2.5 MC simulations
To verify the retrieval algorithm in Section 2.2, a group of simulations are performed by using the FLUKA 4.0 code [16–18]. Figure 8 shows the basic simulation setup: one 150 MeV primary proton beam bombards a slender cylindrical target. The primary beam is Gaussian with FWHM = 5 mm and zero momentum spread. The radius of the target is 2 cm. The target length is larger than the beam range, as shown in columns 5 and 6 of Table 1. As the 150 MeV proton range in H0C is 17.98 cm, its length is set to 20 cm to cover a short neutron flight path. The lengths of H1C ∼ H20C and H0O ∼ H20O are set to 18 cm. There is a virtual gamma detector (FLUKA scorer) around the target surface, and all the photons from the target to the scorer are recorded by the USRBDX card in the range of 0–10 MeV with 1,000 energy bins. USRBDX is a virtual detector for a boundary crossing fluence or current estimator [16].
FIGURE 8
TABLE 1
| Molecular formula | Number of carbon atoms | Number of hydrogen atoms | Density (g/cm3) | Range (cm) | Length (cm) | Radius (cm) | USRBDX area (cm2) | Range volume (cm3) | Target volume (cm3) | Range mass (g) | Target mass (g) | Molecular molar mass (g/mol) | Molecular molar amount in range | Hydrogen molar amount in target |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| H0C | 1 | 0 | 1 | 17.98 | 20.00 | 2 | 251.33 | 225.94 | 251.33 | 225.94 | 251.33 | 12 | 18.83 | 0.00 |
| H1C | 1 | 1 | 1 | 16.39 | 18.00 | 2 | 226.19 | 205.96 | 226.19 | 205.96 | 226.19 | 13 | 15.84 | 17.40 |
| H2C | 1 | 2 | 1 | 15.25 | 18.00 | 2 | 226.19 | 191.64 | 226.19 | 191.64 | 226.19 | 14 | 13.69 | 32.31 |
| H3C | 1 | 3 | 1 | 14.27 | 18.00 | 2 | 226.19 | 179.32 | 226.19 | 179.32 | 226.19 | 15 | 11.95 | 45.24 |
| H4C | 1 | 4 | 1 | 13.47 | 18.00 | 2 | 226.19 | 169.27 | 226.19 | 169.27 | 226.19 | 16 | 10.58 | 56.55 |
| H5C | 1 | 5 | 1 | 12.78 | 18.00 | 2 | 226.19 | 160.60 | 226.19 | 160.60 | 226.19 | 17 | 9.45 | 66.53 |
| H6C | 1 | 6 | 1 | 12.38 | 18.00 | 2 | 226.19 | 155.57 | 226.19 | 155.57 | 226.19 | 18 | 8.64 | 75.40 |
| H8C | 1 | 8 | 1 | 11.84 | 18.00 | 2 | 226.19 | 148.79 | 226.19 | 148.79 | 226.19 | 20 | 7.44 | 90.48 |
| H10C | 1 | 10 | 1 | 11.22 | 18.00 | 2 | 226.19 | 140.99 | 226.19 | 140.99 | 226.19 | 22 | 6.41 | 102.82 |
| H12C | 1 | 12 | 1 | 10.81 | 18.00 | 2 | 226.19 | 135.84 | 226.19 | 135.84 | 226.19 | 24 | 5.66 | 113.10 |
| H15C | 1 | 15 | 1 | 10.36 | 18.00 | 2 | 226.19 | 130.19 | 226.19 | 130.19 | 226.19 | 27 | 4.82 | 125.66 |
| H20C | 1 | 20 | 1 | 9.77 | 18.00 | 2 | 226.19 | 122.77 | 226.19 | 122.77 | 226.19 | 32 | 3.84 | 141.37 |
Parameters for the MC simulation of combinations of carbon and hydrogen.
This setup means that we ignore the energy difference of proton travel through the target in the beam direction. In fact, the proton beam enters and loses energy as it travels. Nuclei near the front face of the target experience higher energy protons than those toward the back of the target. The reaction probabilities change as a function of energy, so different reactions dominate at different depths in the target. We divide the 20-cm water target into five layers and obtain the gamma spectra at different depths by using 5 USRBDXs, each of which is 4 cm. As shown in Figure 9, the characteristic gamma counts increase as protons travel along the target and suddenly decrease in layer 5. However, for this feasibility study, we only adopt one USRBDX to score spectra for the full length to simplify the calculation.
FIGURE 9
2.5.1 FLUKA setup for elemental molar PGS
The first group of simulations is completed to obtain molar spectroscopy results for carbon, oxygen, and hydrogen. Thus, two series of sample materials are simulated with 150 MeV proton irradiation. One series is hydrocarbons: , and the other is . Certainly, most of these compounds do not exist in our natural world except and . However, in our feasibility MC simulations, FLUKA can set up different combinations of atoms by using the COMPOUND card, which defines a mixture of different atoms. To simplify the calculation, the densities of all sample materials are set to the same as water .
As the focus of this study is not on range measurement, we measure the range by scoring proton fluence along the beam axis when 90% of the primary protons have stopped in the target. The USRBIN card, with a Cartesian grid of voxels, is adopted for such a measurement. The USRBIN card can score the distribution of proton fluence or other quantities in a regular spatial structure.
Not all the material in the target can contribute to the PGS spectra scored by USRBDX. After all the protons stop in the target, there will be no gamma radiation directly induced by protons. However, after reaching the longitudinal depth of the proton range, there are still neutrons from previous reactions along the beam path that can produce 2.23 MeV gamma rays. Thus, to calculate the molar PGS values for carbon and oxygen, we should use the range mass; for hydrogen, we should use the target mass. All the discussed parameters are shown in Tables 1, 2.
TABLE 2
| Molecular formula | Number of oxygen atoms | Number of hydrogen atoms | Density (g/cm3) | Range (cm) | Length (cm) | Radius (cm) | USRBDX area (cm2) | Range volume (cm3) | Target volume (cm3) | Range mass (g) | Target mass (g) | Molecular molar mass (g/mol) | Molecular molar amount in range | Hydrogen molar amount in target |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| H0O | 1 | 0 | 1 | 18.00 | 18.00 | 2 | 226.19 | 226.19 | 226.19 | 226.19 | 226.19 | 16 | 14.14 | 0.00 |
| H1O | 1 | 1 | 1 | 17.22 | 18.00 | 2 | 226.19 | 216.39 | 226.19 | 216.39 | 226.19 | 17 | 12.73 | 13.31 |
| H2O | 1 | 2 | 1 | 16.00 | 18.00 | 2 | 226.19 | 201.06 | 226.19 | 201.06 | 226.19 | 18 | 11.17 | 25.13 |
| H3O | 1 | 3 | 1 | 15.20 | 18.00 | 2 | 226.19 | 191.01 | 226.19 | 191.01 | 226.19 | 19 | 10.05 | 35.71 |
| H4O | 1 | 4 | 1 | 14.37 | 18.00 | 2 | 226.19 | 180.58 | 226.19 | 180.58 | 226.19 | 20 | 9.03 | 45.24 |
| H5O | 1 | 5 | 1 | 13.79 | 18.00 | 2 | 226.19 | 173.29 | 226.19 | 173.29 | 226.19 | 21 | 8.25 | 53.86 |
| H6O | 1 | 6 | 1 | 13.39 | 18.00 | 2 | 226.19 | 168.26 | 226.19 | 168.26 | 226.19 | 22 | 7.65 | 61.69 |
| H8O | 1 | 8 | 1 | 12.61 | 18.00 | 2 | 226.19 | 158.46 | 226.19 | 158.46 | 226.19 | 24 | 6.60 | 75.40 |
| H10O | 1 | 10 | 1 | 12.00 | 18.00 | 2 | 226.19 | 150.80 | 226.19 | 150.80 | 226.19 | 26 | 5.80 | 87.00 |
| H12O | 1 | 12 | 1 | 11.39 | 18.00 | 2 | 226.19 | 143.13 | 226.19 | 143.13 | 226.19 | 28 | 5.11 | 96.94 |
| H15O | 1 | 15 | 1 | 10.99 | 18.00 | 2 | 226.19 | 138.10 | 226.19 | 138.10 | 226.19 | 31 | 4.45 | 109.45 |
| H20O | 1 | 20 | 1 | 10.40 | 18.00 | 2 | 226.19 | 130.69 | 226.19 | 130.69 | 226.19 | 36 | 3.63 | 125.66 |
Parameters for the MC simulation of combinations of oxygen and hydrogen.
2.5.2 FLUKA setup for PGS
The second group of simulations is completed to obtain full energy PGS spectra of PMMA, pentanediol, and ethanediol. The configuration is similar to that in Section 2.5.1. The only difference is that the COMPOUND card simultaneously contains carbon, oxygen, and hydrogen with different numbers of atoms. In addition, their densities are set to , , and , respectively, in the MC simulations. All the discussed parameters are presented in Table 3.
TABLE 3
| Sample material | Number of carbon atoms | Number of hydrogen atoms | Number of oxygen atoms | Density (g/cm3) | Range (cm) | Length (cm) | Radius (cm) | USRBDX area (cm2) | Range volume (cm3) | Target volume (cm3) | Range mass (g) | Target mass (g) | Molecular molar mass (g/mol) | Molecular molar amount in range | Hydrogen molar amount in target |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 5 | 8 | 2 | 1.18 | 14.00 | 18.00 | 2 | 226.19 | 175.93 | 226.19 | 207.60 | 226.19 | 100 | 2.08 | 21.35 | |
| 5 | 12 | 2 | 0.9939 | 15.78 | 18.00 | 2 | 226.19 | 198.30 | 226.19 | 197.09 | 224.81 | 104 | 1.90 | 25.94 | |
| 2 | 6 | 2 | 1.1088 | 14.63 | 18.00 | 2 | 226.19 | 183.85 | 226.19 | 203.85 | 250.80 | 62 | 3.29 | 24.27 |
Parameters for the MC simulation of the sample materials.
3 Results
3.1 Element molar PGS
The PGS spectra from combinations of carbon and hydrogen bombarded by 150 MeV protons are shown in Figure 10. Characteristic gamma lines from carbon decrease with an increasing number of hydrogen atoms, while 2.23 MeV gamma lines increase as the number of hydrogen atoms increases. This is consistent with our intuition that the carbon concentration is high in low hydrogen combinations and the opposite for hydrogen concentrations. In addition, 2.23 MeV lines from neutron absorption radiation will be investigated carefully later.
FIGURE 10
For every combination, the beam range, which decreases as the number of hydrogen atoms increases, is measured. As shown in Table 1, the range volume (column 9) and range mass (column 11) can be determined. Then, the PGS spectra are divided by the molecular molar amount in range, which is the molar PGS, as presented in Figure 11. The difference is that the characteristic gamma lines from carbon come to the same level in molar PGS except at an energy of approximately 2.23 MeV.
FIGURE 11
We extract peak gamma lines (2.00 MeV, 2.23 MeV, 2.88 MeV, 4.44 MeV, and 6.12 MeV) from Figures 10, 11. As shown on the left of Figure 12, the 2.00 MeV, 2.88 MeV, 4.44 MeV, and 6.12 MeV counts have a linear relationship with carbon concentration, which has been confirmed by experiments [40]. As shown on the right of Figure 12, the molar count curves of these selected gamma lines are at the same level except for 2.23 MeV. Thus, the basic assumption of a linear carbon term in Eq. 2 is reasonable.
FIGURE 12
The molar count standard errors are at 2 MeV, at 2.23 MeV, at 2.88 MeV, at 4.44 MeV, and at 6.12 MeV as shown in Figure 13. Thus, the count deviation for 2.23 MeV is at least an order of magnitude higher than that of the other gamma lines. We will analyze the nonlinear behavior in the next section.
FIGURE 13
The PGS spectra from combinations of oxygen and hydrogen bombarded by 150 MeV protons are analyzed by the same pattern as carbon in the supplementary material. While their behavior is similar to those from the combinations of carbon and hydrogen, their characteristic gamma lines are more abundant.
Thus, matrix in Eq. 10 can be determined by the element molar PGS.
3.2 Molar PG count for hydrogen
We extract the 2.23 MeV gamma count from the PGS spectrum described in Section 3.1. Additionally, we obtain a nonlinear PG count per mol hydrogen with its number of atoms, as the dotted lines show in Figure 14. The count increases when the number of hydrogen atoms is 1–12 and then decreases. The reason is that the more hydrogens there are in a molecule, the higher the neutron absorption probability. However, the presence of too many hydrogens will stop neutrons from going further to the next molecule as energy is lost in the collision of neutrons and hydrogen nuclei. Thus, hydrogen appears to be an enhancer in low hydrogen combinations. The ratio term in Eqs 4, 5 is reasonable when is less than ∼12. is the upper limit of human organ tissue from the ICRP 110 report as list in its ANNEX B section [49].
FIGURE 14
Here, we must choose a molar PG count values to assign and in Eq. 13. In principle, the PG counts corresponding to HC and HO should be the best option. However, most of the combinations with a large number of hydrogen atoms are not natural; in this study, we choose the 2.23 MeV molar count from and . That is,
3.3 Element amount measurement by PGSRA
Figure 15 shows the full energy PGS spectra of PMMA, pentanediol, and ethanediol detected by the USRBDX card during the virtual MC simulations described in Section 2.5.2. To solve Eq. 10, we select the peak points to assign matrix . Most gamma lines are listed in Table 4. The MC simulated energy in column 3 is in good agreement with the experimental results [50].
FIGURE 15
TABLE 4
| Group | Energy experiment (MeV) | Energy MC (MeV) | Relative error (%) | Transition | Reaction | Mean life (s) |
|---|---|---|---|---|---|---|
| Deprecated | 0.718 | 0.72 | 0.28 | 10B∗0.718 → g.s. | 12C(p,x)10B∗ | 1.0 × 10−9 |
| 16O(p,x)10B∗ | 1.0 × 10−9 | |||||
| 12C(p,x)10C(ϵ)10B∗ | 27.8 | |||||
| 1.022 | 1.02 | 0.20 | 10B∗1.740→10B∗0.718 | 12C(p,x)10B∗ | 7.5 × 10−15 | |
| 16O(p,x)10B∗ | 7.5 × 10−15 | |||||
| Low | 1.635 | 1.63 | 0.31 | 14N∗3.948→14N∗2.313 | 16O(p,x)14N∗ | 6.9 × 10−15 |
| 2.000 | 2.00 | 0.00 | 11C∗2.000 → g.s. | 12C(p,x)11C∗ | 1.0 × 10−14 | |
| 2.124 | 2.12 | 0.19 | 11B∗2.125 → g.s. | 12C(p,x)11B∗ | 5.5 × 10−15 | |
| 2.313 | 3.31 | 0.13 | 14N∗2.313 → g.s. | 16O(p,x)14N∗ | 9.8 × 10−14 | |
| 2.742 | 2.82 | 2.84 | 16O∗8.872 → 16O∗6.130 | 16O(p,p’)16O∗ | 1.8 × 10−13 | |
| 3.684 | 3.69 | 0.16 | 13C∗3.685 → g.s. | 16O(p,x)13C∗ | 1.6 × 10−15 | |
| 3.853 | 3.88 | 0.70 | 13C∗3.854 → g.s. | 16O(p,x)13C∗ | 1.2 × 10−11 | |
| High | 4.438 | 4.43 | 0.18 | 12C∗4.439 → g.s. | 12C(p,p’)12C∗ | 6.1 × 10−14 |
| 16O(p,x)12C∗ | 6.1 × 10−14 | |||||
| 5.180 | 5.18 | 0.00 | 15O∗5.181 → g.s. | 16O(p,x)15O∗ | < 4.9 × 10−14 | |
| 6.129 | 6.11 | 0.31 | 16O∗6.130 → g.s. | 16O(p,p’)16O∗ | 2.7 × 10−11 | |
| 6.476 | 6.43 | 0.71 | 11C∗6.478 → g.s. | 12C(p,x)11C∗ | < 8.7 × 10−15 | |
| 6.741 | 6.76 | 0.28 | 11B∗6.743 → g.s. | 12C(p,x)11B∗ | 4.3 × 10−20 | |
| 6.790 | 6.79 | 0.00 | 11B∗6.792 → g.s. | 12C(p,x)11B∗ | 5.6 × 10−19 |
Gamma lines from proton reactions with oxygen and carbon. Experimental energy data, transition, reaction, and mean life from accelerated particle interactions [50].
These points are sorted into different groups:
All: All peaks from 1.5 to 7 MeV; energies lower than 1.5 MeV are not chosen in this algorithm, as some of these energies are not generated from pure (p,γ) channels. These gamma lines are not as distinct as those from Compton scattering.
Low: Low-energy peaks from 1.5 to 4 MeV.
High: High-energy peaks from 4 to 7 MeV.
Medium-high: Peaks of medium energy 2 MeV to high energy 7 MeV.
Every group will include at least 3 points, and thus, the numerical solution of the element molar amount matrix can be determined by the least square method for each group.
Table 5 lists the numerical solutions of and for each group. The molar amount of carbon is approximately 10, and the ratio of is close to 2.5. The last column is the error of with respect to the atomic ratio in a PMMA molecule. As a result, the solution of the medium-high group is the best solution, with smallest relative error of 0.78% as shown in Table 5 with bold number.
TABLE 5
| group | Mc [mol] | Mo [mol] | Mc/Mo | Mc/Mo relative error (%) |
|---|---|---|---|---|
| All | 10.1522 | 4.2419 | 2.3933 | 4.27 |
| Low | 10.1350 | 4.2509 | 2.3842 | 4.63 |
| High | 10.2351 | 4.2068 | 2.4330 | 2.68 |
| Medium-high | 10.3469 | 4.1066 | 2.5196 | 0.78 |
Numeric solution of PMMA.
For pentanediol, the best solution is the high-energy group, as shown in Table 6 with bold number.
TABLE 6
| group | Mc [mol] | Mo [mol] | Mc/Mo | Mc/Mo relative error (%) |
|---|---|---|---|---|
| All | 9.2649 | 4.0280 | 2.3001 | 7.99 |
| Low | 9.1582 | 4.0814 | 2.2439 | 10.24 |
| High | 9.7719 | 3.8373 | 2.5466 | 1.86 |
| Medium-high | 9.5617 | 3.9140 | 2.4430 | 2.28 |
Numeric solution of pentanediol.
For ethanediol, the best solution is the high-energy group, as shown in Table 7 with bold number.
TABLE 7
| group | Mc [mol] | Mo [mol] | Mc/Mo | Mc/Mo relative error (%) |
|---|---|---|---|---|
| All | 6.3268 | 6.6354 | 0.9535 | 4.65 |
| Low | 6.1875 | 6.6743 | 0.9271 | 7.29 |
| High | 6.8350 | 6.7923 | 1.0063 | 0.63 |
| Medium-high | 6.6864 | 6.7517 | 0.9903 | 0.97 |
Numeric solution of ethanediol.
Medium- and high-energy points are preferred for the numerical solution as these characteristic gamma lines are more distinct. Now, and are determined, and can be extracted from the PGS spectra of the sample materials. and have been determined in Section 3.2. Then, we can calculate the hydrogen molar amount by using Eq. 13. The sample mass and density can be calculated by using Eq. 15.
The results are summarized in Table 8 and compared with the MC setup values from Table 3. The relative errors of the hydrogen molar amount are approximately 10% for PMMA and pentanediol. However, their densities are consistent with the setup values, and the relative errors are smaller than 2%. On the contrary, for ethanediol, the hydrogen amount relative error is only 0.87%, and the density relative error is 3.31%. This is because although there is a high atom number of hydrogen, it has much lower weight than carbon and oxygen in materials. Thus, density error is dominated by carbon and oxygen.
TABLE 8
| Sample material | Variable | PGSPA result | MC setup value | Relative error (%) |
|---|---|---|---|---|
| PMMA = 0.0349 | [mol] | 4.1066 | 4.1519 | 1.09 |
| [mol] | 10.3469 | 10.3798 | 0.32 | |
| [mol] | 20.6003 | 21.3528 | 3.52 | |
| Sample Mass [g] | 210.4682 | 207.5964 | 1.38 | |
| Density [g/cm3] | 1.17 | 1.18 | 1.38 | |
| Pentanediol = 0.0365 | [mol] | 3.8373 | 3.7901 | 1.24 |
| [mol] | 9.7719 | 9.4754 | 3.13 | |
| [mol] | 20.4229 | 22.7409 | 10.19 | |
| Sample Mass [g] | 199.0820 | 197.0877 | 1.01 | |
| Density [g/cm3] | 1.0040 | 0.99 | 1.01 | |
| Ethanediol = 0.0377 | [mol] | 6.7923 | 6.5758 | 3.29 |
| [mol] | 6.8350 | 6.5758 | 3.94 | |
| [mol] | 19.8981 | 19.7273 | 0.87 | |
| Sample Mass [g] | 210.5943 | 203.8484 | 3.31 | |
| Density [g/cm3] | 1.1455 | 1.11 | 3.31 |
Numeric solutions of and density.
Finally, we can recalculate the PMMA PGS spectra by using our retrieval algorithm, as shown in Figure 16. In the 0–4 MeV range, all peaks are consistent with the MC curve. In the 4–7 MeV range, only the peak at 6.1 MeV shows a small deviation from the MC curve. The retrieved PGSs of pentanediol and ethanediol are in supplementary material.
FIGURE 16
4 Discussion and conclusion
Based on 3 basic assumptions: 1) the sample material has a relationship with its element molar gamma lines in PG spectroscopy, 2) for carbon and oxygen, this relationship is linear, and 3) for hydrogen, this relationship is nonlinear and is based on a neutron absorption reaction; an element retrieval algorithm is developed in the present study. The challenge is how to carefully deal with the nonlinear hydrogen 2.23 MeV gamma count. A series of simulations confirm that these assumptions are reasonable. The mathematical formulas of the PGSRA theoretically explain the linear relationship of the gamma count and element concentration detected by Polf et al. [38, 39] and Paulo et al. [40].
The element amount has been successfully calculated by the PGSRA for the following three sample materials: PMMA, pentanediol, and ethanediol. The relative errors of carbon and oxygen are smaller than 4% for all the sample materials. The relative error of hydrogen, based on its nonlinear behavior, is approximately 10%. In addition, the relative errors of the sample density are smaller than 3.5% for all the sample materials. This new retrieval algorithm is valid for predicting the element amount and density in an unknown material. It may be the first bridge between PGS and medical imaging.
Statements
Data availability statement
The original contributions presented in the study are included in the article. PGS spectra from combinations of oxygen and hydrogen, retrieval pentanediol PGS and ethanediol PGS are in the Supplementary Material, and further inquiries can be directed to the corresponding author.
Author contributions
J-L Wang: theoretical model, FLUKA simulation, and manuscript. X-GW: physics, prompt gamma analysis, and review. Z-FL: medical physics, proton therapy, and review. S-QX: data analysis. Z-HZ: detector. D-QH: neutron adsorption. AR: cloud computing. YZ: detector and physics. CL: nuclear reaction. T-XL: detector and physics. MZ: gamma detector. D-JD: data process. X-DW: data process. H-SR: data process.
Funding
This work was supported by the National Natural Science Foundation of China (U1932209, 11975315, U1867210, and 11905134), the Leading Innovation Project under Grant (LC192209000701 and LC202309000201), and the Continuous Basic Scientific Research Project (BJ20002501 and WDJC-2019-13).
Acknowledgments
We sincerely thank Michael Lu from Mevion Medical Systems for many discussions about imaging devices in proton therapy. The authors also kindly acknowledge the great support from AWS engineers Yin-Xiang Li and Xiao-Chen Ye.
Conflict of interest
A Rusanov was employed by the company Elekta Limited.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2022.961162/full#supplementary-material
References
1.
MacKayRI. Image guidance for proton therapy. Clin Oncol (2018) 30(5):293–8. 10.1016/j.clon.2018.02.004
2.
KnopfACLomaxA. In vivo proton range verification: A review. Phys Med Biol (2013) 58(15):R131–60. 10.1088/0031-9155/58/15/R131
3.
MinCHKimCHYounMYKimJW. Prompt gamma measurements for locating the dose falloff region in the proton therapy. Appl Phys Lett (2006) 89(18):183517. 10.1063/1.2378561
4.
SmeetsJRoellinghoffFPrieelsDStichelbautFBenilovABuscaPet alPrompt gamma imaging with a slit camera for real-time range control in proton therapy. Phys Med Biol (2012) 57(11):3371–405. 10.1088/0031-9155/57/11/3371
5.
VerburgJMSecoJ. Proton range verification through prompt gamma-ray spectroscopy. Phys Med Biol (2014) 59(23):7089–106. 10.1088/0031-9155/59/23/7089
6.
GolnikCHueso-GonzálezFMüllerADendoovenPEnghardtWFiedlerFet alRange assessment in particle therapy based on prompt γ-ray timing measurements. Phys Med Biol (2014) 59(18):5399–422. 10.1088/0031-9155/59/18/5399
7.
DraegerEMackinDPetersonSChenHAverySBeddarSet al3D prompt gamma imaging for proton beam range verification. Phys Med Biol (2018) 63(3):035019. 10.1088/1361-6560/aaa203
8.
RoellinghoffFRichardMHChevallierMConstanzoJDauvergneDFreudNet alDesign of a Compton camera for 3D prompt-γ imaging during ion beam therapy. Nucl Instr Methods Phys Res Section A: Acc Spectrometers Detectors Associated Equipment (2011) 648(Suppl. 1):20–3. 10.1016/j.nima.2011.01.069
9.
PanainoCMVMackayRIKirkbyKJTaylorMJ. A new method to reconstruct in 3D the emission position of the prompt gamma rays following proton beam irradiation. Sci Rep (2019) 9(1):18820–13. 10.1038/s41598-019-55349-7
10.
RichterCPauschGBarczykSPriegnitzMKeitzIThieleJet alFirst clinical application of a prompt gamma based in vivo proton range verification system. Radiother Oncol (2016) 118(2):232–7. 10.1016/j.radonc.2016.01.004
11.
Hueso-GonzálezFRabeMRuggieriTABortfeldTVerburgJM. A full-scale clinical prototype for proton range verification using prompt gamma-ray spectroscopy. Phys Med Biol (2018) 63(18):185019. 10.1088/1361-6560/aad513
12.
WernerTBertholdJHueso-GonzálezFKoeglerTPetzoldtJRoemerKet alProcessing of prompt gamma-ray timing data for proton range measurements at a clinical beam delivery. Phys Med Biol (2019) 64(10):105023. 10.1088/1361-6560/ab176d
13.
MarcatiliSCollotJCurtoniSDauvergneDHostachyJYKoumeirCet alUltra-fast prompt gamma detection in single proton counting regime for range monitoring in particle therapy. Phys Med Biol (2020) 65(24):245033–17. 10.1088/1361-6560/ab7a6c
14.
ZhangHFanPXuTWeiQMaTLiuYet alScintillation detector design study for prompt gamma photon detection in proton therapy monitoring. In: 2017 IEEE Nucl Sci Symp Med Imaging Conf NSS/MIC 2017 - Conf Proc; 21-28 October 2017; Atlanta, GA, USA (2018). p. 7–10. 10.1109/NSSMIC.2017.8532768
15.
ZhangHFanPWangSXiaYXuTWeiQet alDesign and performance evaluation of a BGO + SiPM detector for high-energy prompt gamma imaging in proton therapy monitoring. IEEE Trans Radiat Plasma Med Sci (2020) 4(2):184–93. 10.1109/trpms.2020.2972594
16.
FerrariASalaPRFassòARanftJ. Fluka: A multi-particle transport code. MC-INFN (2005) 773. 10.2172/877507
17.
BöhlenTTCeruttiFChinMPWFassoAFerrariAOrtegaPet alThe FLUKA Code: Developments and challenges for high energy and medical applications. Nucl Data Sheets (2014) 120:211–4. 10.1016/j.nds.2014.07.049
18.
KozłowskaWSBöhlenTTCuccagnaCFerrariAFracchiollaFMagroGet alFLUKA particle therapy tool for Monte Carlo independent calculation of scanned proton and carbon ion beam therapy. Phys Med Biol (2019) 64(7):075012. 10.1088/1361-6560/ab02cb
19.
AgostinelliSAllisonJAmakoKApostolakisJAraujoHArcePet alGEANT4 - a simulation toolkit. Nucl Instr Methods Phys Res Section A: Acc Spectrometers Detectors Associated Equipment (2003) 506(3):250–303. 10.1016/S0168-9002(03)01368-8
20.
CirroneGAPCuttoneGDi RosaFMazzagliaSERomanoFAttiliAet alHadrontherapy: An open source, Geant4-based application for proton-ion therapy studies. In: IEEE Nuclear Science Symposium Conference Record; 24 October 2009 - 01 November 2009; Orlando, FL, USA (2009). p. 4186–9. 10.1109/NSSMIC.2009.5402279
21.
BiegunAKSeravalliELopesPCRinaldiIPintoMOxleyDCet alTime-of-flight neutron rejection to improve prompt gamma imaging for proton range verification: A simulation study. Phys Med Biol (2012) 57(20):6429–44. 10.1088/0031-9155/57/20/6429
22.
JacquetMMarcatiliSGallin-MartelMLBoulyJLBoursierYDauvergneDet alA time-of-flight-based reconstruction for real-time prompt-gamma imaging in proton therapy. Phys Med Biol (2021) 66(13):135003–32. 10.1088/1361-6560/ac03ca
23.
WangJLCruzLALuM. Pixelated prompt gamma imaging detector for online measurement of proton beam: Monte Carlo feasibility study by FLUKA. Radiat Detect Technol Methods (2018) 2(1):4. 10.1007/s41605-017-0032-0
24.
ParkJHKimSHKuYLeeHSKimCHShinDHet alComparison of knife-edge and multi-slit camera for proton beam range verification by Monte Carlo simulation. Nucl Eng Technol (2019) 51(2):533–8. 10.1016/j.net.2018.10.002
25.
HanYTangXGengCShuDGongCZhangXet alInvestigation of in vivo beam range verification in carbon ion therapy using the Doppler shift effect of prompt gamma: A Monte Carlo simulation study. Radiat Phys Chem Oxf Engl 1993 (2019) 162:72–81. 10.1016/j.radphyschem.2019.04.036
26.
GengCHanYTangXShuDGongCAltieriS. Evaluation of using the Doppler shift effect of prompt gamma for measuring the carbon ion range in vivo for heterogeneous phantoms. Nucl Instr Methods Phys Res Section A: Acc Spectrometers Detectors Associated Equipment (2020) 959:163439. 10.1016/j.nima.2020.163439
27.
PennazioFFerreroVD’OnghiaGGarbolinoSFiorinaEMarti VillarrealOAet alProton therapy monitoring: Spatiotemporal emission reconstruction with prompt gamma timing and implementation with PET detectors. Phys Med Biol (2022) 67(6):065005. 10.1088/1361-6560/ac5765
28.
PolfJCPetersonSCiangaruGGillinMBeddarS. Prompt gamma-ray emission from biological tissues during proton irradiation: A preliminary study. Phys Med Biol (2009) 54(3):731–43. 10.1088/0031-9155/54/3/017
29.
TestaMBajardMChevallierMDauvergneDFreudNHenriquetPet alReal-time monitoring of the Bragg-peak position in ion therapy by means of single photon detection. Radiat Environ Biophys (2010) 49(3):337–43. 10.1007/s00411-010-0276-2
30.
PauschGPetzoldtJBerthelMEnghardtWFiedlerFGolnikCet alScintillator-based high-throughput fast timing spectroscopy for real-time range verification in particle therapy. IEEE Trans Nucl Sci (2016) 63(2):664–72. 10.1109/TNS.2016.2527822
31.
Dal BelloRMagalhaes MartinsPBronsSHermannGKihmTSeimetzMet alPrompt gamma spectroscopy for absolute range verification of 12C ions at synchrotron-based facilities. Phys Med Biol (2020) 65(9):095010. 10.1088/1361-6560/ab7973
32.
IdiriZMazrouHAmokraneABedekS. Characterization of an Am-Be PGNAA set-up developed for in situ liquid analysis: Application to domestic waste water and industrial liquid effluents analysis. Nucl Instr Methods Phys Res Section B: Beam Interactions Mater Atoms (2010) 268(2):213–8. 10.1016/j.nimb.2009.10.185
33.
KhelifiRAmokraneABodeP. Detection limits of pollutants in water for PGNAA using Am-Be source. Nucl Instr Methods Phys Res Section B: Beam Interactions Mater Atoms (2007) 262(2):329–32. 10.1016/j.nimb.2007.06.003
34.
NaqviAAGarwanMANagadiMMMaslehuddinMAl-AmoudiOSBKhateeb-ur-Rehman. Non-destructive analysis of chlorine in fly ash cement concrete. Nucl Instr Methods Phys Res Section A: Acc Spectrometers Detectors Associated Equipment (2009) 607(2):446–50. 10.1016/j.nima.2009.05.192
35.
NaqviAAAl-MatouqFAKhiariFZKhateeb-Ur-RehmanGondalMAIsabAA. Optimization of a prompt gamma setup for analysis of environmental samples. J Radioanal Nucl Chem (2013) 296:215–21. 10.1007/s10967-012-2045-y
36.
JiaWChengCHeiDLingYWangHChenD. Method for correcting thermal neutron self-shielding effect for aqueous bulk sample analysis by PGNAA technique. J Radioanal Nucl Chem (2015) 304(3):1133–7. 10.1007/s10967-015-3962-3
37.
HeiDJiangZJiaWChengCWangHLiJet alThe background influence of cadmium detection in saline water using PGNAA technique. J Radioanal Nucl Chem (2016) 310(1):27–31. 10.1007/s10967-016-4767-8
38.
VerburgJMShihHASecoJ. Simulation of prompt gamma-ray emission during proton radiotherapy. Phys Med Biol (2012) 57(17):5459–72. 10.1088/0031-9155/57/17/5459
39.
PolfJCPanthiRMackinDSMcCleskeyMSaastamoinenARoederBTet alMeasurement of characteristic prompt gamma rays emitted from oxygen and carbon in tissue-equivalent samples during proton beam irradiation. Phys Med Biol (2013) 58(17):5821–31. 10.1088/0031-9155/58/17/5821
40.
Magalhaes MartinsPDal BelloRAckermannBBronsSHermannGKihmTet alPIBS: Proton and ion beam spectroscopy for in vivo measurements of oxygen, carbon, and calcium concentrations in the human body. Sci Rep (2020) 10(1):7007–14. 10.1038/s41598-020-63215-0
41.
JiaWChengCShanQHeiDLingYZhangYet alStudy on the elements detection and its correction in aqueous solution. Nucl Instr Methods Phys Res Section B: Beam Interactions Mater Atoms (2015) 342:240–3. 10.1016/j.nimb.2014.10.010
42.
JeyasugiththanJCameroJNSymonsJJonesPBufflerAGeduldDet alMeasuring prompt gamma-ray emissions from elements found in tissue during passive-beam proton therapy. Biomed Phys Eng Express (2021) 7(2):025013. 10.1088/2057-1976/abe33d
43.
AhdidaCBozzatoDCalzolariDCeruttiFCharitonidisNCimminoAet alNew capabilities of the FLUKA multi-purpose code. Front Phys (2022) 9:1–14. 10.3389/fphy.2021.788253
44.
CleggABFoleyKJSalmonGLSegelRE. Gamma radiation from the medium energy proton bombardment of lithium, beryllium, boron, carbon and nitrogen. Proc Phys Soc (1961) 78(5):681–94. 10.1088/0370-1328/78/5/305
45.
LeeHRKimSHParkJHJungWGLimHKimCH. Prototype system for proton beam range measurement based on gamma electron vertex imaging. Nucl Instr Methods Phys Res Section A: Acc Spectrometers Detectors Associated Equipment (2017) 857:82–97. 10.1016/j.nima.2017.03.022
46.
ParkJHKimSHKuYKimCHLeeHRJeongJHet alMulti-slit prompt-gamma camera for locating of distal dose falloff in proton therapy. Nucl Eng Technol (2019) 51(5):1406–16. 10.1016/j.net.2019.03.008
47.
WangJLCruzLAWuQBWangQWeiYWangHK. Radiation shielding design of a compact single-room proton therapy based on synchrotron. Nucl Sci Tech (2020) 31(1):1. 10.1007/s41365-019-0712-1
48.
OtukaNDupontESemkovaVPritychenkoBBlokhinAIAikawaMet alToward more complete and accurate experimental nuclear reaction data library (EXFOR) - international collaboration between nuclear reaction data centres (NRDC). New York: Nuclear Data Sheets (2014). p. 8–11.
49.
SowbyFD. Annals of the ICRP 110. Ann ICRP (1981) 6(1):1. 10.1016/0146-6453(81)90127-5
50.
MurphyRJ. NUCLEAR DEEXCITATION GAMMA-RAY LINES FROM ACCELERATED PARTICLE INTERACTIONS benzion kozlovsky. The Astrophysical Journal Supplement Series (2002). 2, p. 523–41. Published online.
Summary
Keywords
proton therapy, prompt gamma spectroscopy retrieval algorithm, Monte Carlo, FLUKA, neutron absorption
Citation
Wang J-L, Wu X-G, Li Z-F, Xie S-Q, Hei D-Q, Zhao Z-H, Rusanov A, Zheng Y, Li C-B, Li T-X, Zheng M, Wang X-D, Ding D-J and Ruan H-S (2022) Prompt gamma spectroscopy retrieval algorithm for element and density measurements accelerated by cloud computing. Front. Phys. 10:961162. doi: 10.3389/fphy.2022.961162
Received
04 June 2022
Accepted
12 October 2022
Published
10 November 2022
Volume
10 - 2022
Edited by
Angelo Schiavi, Sapienza University of Rome, Italy
Reviewed by
Elisa Fiorina, National Institute of Nuclear Physics of Turin, Italy
Dominique Thers, IMT Atlantique Bretagne-Pays de la Loire, France
Updates
Copyright
© 2022 Wang, Wu, Li, Xie, Hei, Zhao, Rusanov, Zheng, Li, Li, Zheng, Wang, Ding and Ruan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jin-Long Wang, 1561103570@qq.com
This article was submitted to Medical Physics and Imaging, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.