You are viewing a javascript disabled version of the site. Please enable Javascript for this site to function properly.
Go to headerGo to navigationGo to searchGo to contentsGo to footer
In content section. Select this link to jump to navigation

Evaluating photoneutron dose equivalents for lung cancer using PMMA phantoms undergoing 15 MV IMRT

Abstract

BACKGROUND:

According to statistics of the Ministry of Health and Welfare in 2017, the second leading cause of death in Taiwan was lung cancer.

OBJECTIVE:

Routine treatment planning does not consider photoneutron dose equivalent (PNDE) of patient induced secondary radiation resulting from primary exposure of lung cancer. However, such treatment is potentially important for improving estimates of health risks.

METHODS:

This study used 10, 30, 50, 70, and 90 kg of polymethylmethacrylate (PMMA) phantoms as patient to measure PNDE varying anatomical area during lung cancer of intensity modulated radiotherapy (IMRT) treatment. Paired thermoluminescent dosimeters (TLD-600 and 700) were calibrated using university reactor neutrons. TLDs were inserted into phantom which was closely corresponded of the represented tissues or organs.

RESULTS:

Neutron doses (ND) of organ or tissue (NDT) were determined in these phantoms using paired TLDs approach. The risks of incurring fatal secondary malignancies, maximum statistical and total errors were estimated. Evaluated PNDE ranged from 0.80 ± 0.12 to 0.56 ± 0.08 mSv/Gy for these phantoms.

CONCLUSION:

The estimated NDT decreased with increasing distance that is from the central axis. Evaluated PNDE and ND𝑠𝑘𝑖𝑛 for these phantoms were discussed. This investigation also identified secondary risks associated with PNDE relating to radiation protection.

1.Introduction

Medical electron linear accelerators (linac) operating above 10 MV produce high-energy bremsstrahlung photon, then induce neutrons via (x, yn) reactions [1, 2, 3]. These reactions occur at the accelerator target (tungsten, iron), flattening filters, shielding materials, collimator jaws (inside the linac head), and in the patient body [1, 2, 3, 4]. The linac treatment room contains a mixed radiation field produced by the induced neutrons and photons. The medical linac, Siemens Primus (Siemens AG, Munich, Germany) at the Radiation Oncology Department of Changhua Christian Hospital (CCH) which could provide energies that accelerated 10 and 15 MV voltages. The linac is designed of intensity modulated radiotherapy (IMRT). The fast neutron flux (Φf) for 15 MV of linac was 1.75 × 104 (n/cm2s)/Gy-X, as well as flux of thermal neutron (Φth) was 1.26 × 105 (n/cm2s)/Gy-X at SSD of 100 cm, respectively [5]. Patient thus was exposed to significant undesirable radiations comprising photons and neutrons leakage (both direct and scattered components) [1]. Knowledge of the undesirable photoneutron dose equivalent (PNDE) in patients of different body weights is necessary. These data were necessary to determine safety and improved treatment plan of lung cancer in patients [1]. According to ICRP 60, the radiation weighting factor (WR) for neutrons ranging from 5 to 20 MeV is quite high, and one would expect significant PNDE to the critical organs, like gonads and lung [5, 6, 7]. A series of measurements and calculations were carried out to determine the PNDE and neutron dose of organ or tissue (NDT) in the patient undergoing 15 MV linac treatments for lung cancer. The authors are aware of limited studies assessing the secondary ND that results from differences in patient weight from step-and-shot IMRT [1]. The experiment described here is to use polymethylmethacrylate (PMMA) phantoms because of different body-weight with the same anthropometric-shaped skeletons as patient surrogates to assess the PNDE of various anatomical areas. Paired TLDs dosimeters (TLD-600 and 700) (Harshaw Chemical Company, Solon, OH, USA) were arranged and inserted. Simple equation could be evaluated PNDE during phantoms’ lung cancer. High ND𝑙𝑢𝑛𝑔 values were recorded up to 1.02 ± 0.13 mSv/Gy of Rando. This study is recommended to prevent extra ND to patients and could be satisfied the ALARA [7].

2.Materials and methods

2.1Polymethylmethacrylate (PMMA) phantom measurement

Lung treatment plans of IMRT have been developed, delivered to anthropomorphic Rando phantom. This was made of Radiology Support Devices, Long Beach, Canada. Treatment planes of these phantoms were reviewed by professional medical doctor (Mu-Tai Liu, MD) and senior radiotherapist (Sing-Sheng Huang) [7, 8, 9, 10, 11]. To calculate the PNDE associated with IMRT treatment, this study conducted paired TLDs measurement in Rando and PMMA phantoms. The Rando phantom was 175 cm and 73.5 kg. This phantom was comprised 35 axial segments, those contained simulated human skeleton as well as lung soft tissue. Rando is suitable for IMRT radiation measurements [4]. The phantom was sliced at 2.5 cm intervals and holes were drilled in each slice for TLD-600 insertions. Paired TLDs were loaded into each slice of phantoms as well as can use to directly assess PNDE and secondary radiation to phantoms.

This study developed five phantoms using anthropometric skeletons made of PMMA and epoxy-resin to simulate average body-weights of 90, 70, 50, 30 and 10 kg [12]. The PMMA phantom was designed as published [9]. PMMA and epoxy-resin were chosen as the neutron and photon transport properties. Density were formulated as follows: skeleton-cortical-bone, 1.486; skin, 1.105; and lung-tissue 0.296 g/cm3 of PMMA, respectively [9, 12]. The PMMA phantoms were made of the human design; contained 31 slices. Those were represented the head, neck, torso and abdomen. These were without two arms and legs [11, 12]. Each slice was drilled with an 11 mm diameter hole to enable insertion of the TLDs into organs. The dimensions and physical properties of theses phantoms is listed in Table 1.

Table 1

Rando, polymethylmethacrylate (PMMA) phantoms’ dimension and physical properties

PhantomRandoPolymethylmethacrylate (PMMA)
Weight (kg)1701030507090
Height (cm)294.550788493112
Weight (kg)234.56.7519.031.544.157
Slices353131313131
cm/slices2.51.62.32.73.03.6

2.2Paired TLDs

For subtracting the signal (SN) contributed by photons from total SNs of TLD-600, and 700, all neutron measurements were performed using paired TLDs approach composed with a batch of 88 high sensitivity 3.2 × 3.2 × 0.89 mm3 dosimeters [10]. TLD-600 comprises 95.6% 6Li, and 4.4% 7Li. TLD-700 is 99.7% 7Li. 6Li and 7Li of these TLDs could provide different responses to neutron reactions of the 6Li (n, α)3H nuclear reaction [9]. Before irradiation, paired TLDs were annealed by heating to 400C (1 hr) and 200C (2 hr) in furnace 47900. Then TLDs were cooled to normal temperature before use [10]. Following IMRT exposure, a fully automated Harshaw 3500 reader made of Bicron NE Solon, OH was used to analyzed these TLDs. The readout SN involved a two steps, (A) heating to 300C at speed of 10C per second, holding for a further 1 second, (B) heating to 50C (1 second) [10].

2.3Calibrated paired TLDs

The ND calibrations were performed at the boron neutron capture therapy (BNCT) beams of National Tsing Hua University (NTHU). Using BNCT beams are 3.10 × 107 n/cm2s (Φf), as well as 1.31 × 108 n/cm2s (Φth) measured at the exit of the beam for reactor power of 1 MW [9]. Additionally, the Φth kerma factor of TLD-600s was calibrated as 1.40 × 10-10 mGy/(n/cm2s) by Monte Carlo simulation [9]. For paired TLDs’ photon dose calibration, linearity calibration was carried by Siemens Primus of CCH’ 6 MV therapeutic unit [6]. This machine typically provides a dose of 400 MU (monitor unit)/min, corresponding to 3.3 Gy/min. For 10 cm × 10 cm field size, 1MU equals approximately 1 cGy at SSD 100 cm. Paired TLDs were arranged in 10-cm-diameter methacrylate boxes. An ionization chamber (IC) like pencil-shaped, Wellhofer IC-69, CT Probes Models 350407 made by Nuclear Associates Victoreen, Inc., containing 0.6 cm3 was used to measure X-ray radiation. The IC was connected by way of 0.9 m low-noise cable to a SN digitizing preamplifier. Due to subtracting photon response, ND in these phantoms is calculated as

(1)
SN600n+γ=SN600n+SN600γ
(2)
SN700n+γ=SN700n+SN700γ

Where SN600n+γ is denoted as total TLD-600 SN. SN600n is net neutron SN. SN600γ is γ-ray contribution to TLD-600 SN. SN700n+γ is denoted as total TLD-700 SN. SN700n and SN700γ are the respective contribution to TLD-700 SN of neutron and γ-ray.

(3)
SN600n=SN600n+γ×CF600n

CF600n is conversion factor (CF) of ND of TLD-600 calibrated by BNCT.

2.4Computed tomography (CT) scans and treatment planning for lung cancer

Before each treatment, numerous measurements were conducted using inhomogeneous PMMA phantoms. This study obtained archival Computed tomography (CT) scan images of the phantoms at CCH. A 3 cm × 3 cm × 3 cm lung tumor, in the center of left lung at a depth of 5 cm, is shown as Fig. 1a–c. Organ contours’ CT images were transferred into Siemens Primus for IMRT. CT protocol scans of these phantoms for three-dimensional visualizations of the lung region were made by way of a multi-slice CT (MSCT) (high speed, General Electric). The settings of mA and kV of the CT scanner were varied according to the phantoms’ weights. The lung block’ size and place was obtained by way of the treatment planning system. After loading these TLDs into phantoms, three repeat scans were conducted over the specified area. All technique parameters as well as displayed quantities were recorded. The IMRT approaches were delivered to the phantoms with 200 cGy prescribed to the planned target volume (PTV) to facilitate detection of the small doses outside the treatment field.

Figure 1.

Use polymethylmethacrylate (PMMA) phantoms as patient surrogates (a) outer appearances (b) from 9th to 20th scout view of 70 kg (c) 15th CT slice of 70 kg phantom.

Use polymethylmethacrylate (PMMA) phantoms as patient surrogates (a) outer appearances (b) from 9th to 20th scout view of 70 kg (c) 15th CT slice of 70 kg phantom.

2.5Paired TLDs inserted into Rando and PMMA phantoms

Calculation points are located in the lung and tumor center in the superior portion of the torso of these phantoms of radiotherapy 15 MV linac during lung cancer treatment [1]. Forty-four groups, each containing 2 TLD-600 and 700s, were inserted into phantoms at the positions of slices 3–35 (Rando) and 3–31 (PMMA), respectively. For all measuring locations of NDs, two TLD-600s were irradiated together as well as average signal determined. Four group TLDs were attached to 15th’s surface of the Rando’ slices in order to measure the ND𝑠𝑘𝑖𝑛.

2.6Evaluation of the PNDE

PNDE were determined for 16 critical organs and remainders recommended by ICRP 103 as published [12]. In addition, prostate of male and uterus of female selected as at particular risk of fatal secondary malignancy recommended by ICRP 103 [11, 13]. Table 2 lists the positions of 44 groups imparted within these phantoms and weighting factors (WT) [5]. The doses listed for the gonads and testes, as well as for the uterus and ovaries, were administered to males and females, respectively.

Table 2

44 ND measurement points, and the weighting factors WT recommended by ICRP 103

OrganMeasurement pointsICRP 103Practical evaluation
BrainBrain0.010.01
SalivarySalivary0.010.01
ThyroidThyroid0.040.04
Red bone marrow0.12
Clavicle vertebrae0.04
Thoracic vertebra0.04
Thighbone femur0.04
LungLung0.120.12
SkinSkin0.010.01
EsophagusEsophagus0.040.04
BreastBreast0.120.12
Bone surface0.01
Clavicle vertebrae0.025
Thoracic vertebra0.025
Thighbone0.025
Rib0.025
LiverLiver0.040.04
Colon0.12
Ascending0.03
Transverse0.03
Descending0.03
Sigmoid flexure0.03
StomachStomach0.120.12
BladderBladder0.040.04
GonadsTestes0.080.08
Remainder0.12
Lens0.017
Herat0.017
Pancreas0.017
Spleen0.017
Kidney0.017
Small intestine0.017
Prostate (male)/uterus (female)0.017
Total1.001.000

2.7Evaluating NDT

The PNDE was obtained for the position of chest PA. Meanwhile, PNDE denotes the total weighted equivalent doses of 16 organs and the “remainder” recommended by ICRP 103 as listed in Table 2, expressed as follows:

(4)
𝑃𝑁𝐷𝐸=TWTHT

HT (equivalent dose), is multiplied by WR, and obtained by

(5)
HT=RWRNDT

Where NDT is the averaged ND of the organs. PNDE is in J/kg unit and thus is termed Gray (Gy). Complete prescribed photon doses (200 cGy) of 15MV were delivered thus enabling detection of small ND outside the treatment field owing to radiation safety concerns. All measurements were conducted in triplicate for minimize fluctuations of TLDs. For small organs like the thyroid and gonad, ND𝑡ℎ𝑦𝑟𝑜𝑖𝑑 and ND𝑔𝑜𝑛𝑎𝑑 obtained from the TLD-600 inserted in these organ centers served as the organ dose. For larger organs, like lung and colon, many TLDs were placed into each slices. The mean ND values were regarded as ND𝑙𝑢𝑛𝑔 as well as ND𝑐𝑜𝑙𝑜𝑛 (cf. Table 2).

3.Results

3.1Calibration of TLD-600

The response linearity of photon dose for paired TLDs was tested ranging from μSv to mSv. Results demonstrated a good linearity at CCH. The dose-response of TLD-700 is y (mSv) = 2.05 + 0.144 ×SN700γ (nC), R2 is equal to 0.99893 for TLD-600 as well as 0.9994 for TLD-700, respectively. From measurements to calculate paired TLDs response to BNCT neutron, CF600n from Eq. (3), is equal to 1.167 × 10-4 mSv/nC, herein. This is in agreement with Hsu [9].

3.2Photoneutron dose equivalent (PNDE)

Figure 2.

Estimate of photoneutron dose equivalent (PNDE) and skin dose per Gy (mSv/Gy) undergoing lung cancer of 15 MV IMRT treatments.

Estimate of photoneutron dose equivalent (PNDE) and skin dose per Gy (mSv/Gy) undergoing lung cancer of 15 MV IMRT treatments.

Specifically, no standardized whole-body phantoms exist for evaluating the PNDE of different weights. Thus, this investigation developed five phantoms to simulate different body-weights, from 10 to 90 kg. As illustrated in Fig. 2, PNDE at the different weights decrease abruptly with body-weight. Meanwhile, the maximum ND𝑠𝑘𝑖𝑛 measurements at the 15th slice were 1.98 ± 0.30, 2.03 ± 0.26, 1.85 ± 0.24, 1.44 ± 0.23, 1.18 ± 0.27, and 1.01 ± 0.23 mSv/Gy for the Rando and PMMA phantoms, respectively. Fitting an exponential function (Eq. (6)) to this trend yielded the PNDE and ND𝑠𝑘𝑖𝑛’s relationship. Furthermore, these values reduced exponentially with increasing body-weight.

(6)
𝑃𝑁𝐷𝐸=0.85-0.00519M+2.15×10-5M2,R2=0.97493
(7)
ND𝑠𝑘𝑖𝑛=2.23-0.0166M+3.04×10-5M2,R2=0.98265

Where PNDE and ND𝑠𝑘𝑖𝑛 were administered in mSv/Gy, while M denotes PMMA phantoms in kg with R2 is coefficient of determination.

4.Discussions

4.1Estimates of the secondary NDT of other organs

Most neutrons those originate in the gentry head, then move from the tumor center to the patients lens and gonad during IMRT. Due to the ND of the 44 measurement (Table 2), estimates were evaluated that reflected the ND to anywhere in bodies of different body-weights for lung cancer treatment. Figure 3a–g present the estimates of the total NDT on the distance away from the five phantoms’ lung edge. Estimated NDT decreased with increasing distance from central axis. This phenomenon is also the same as the results published by Chao et al. and Lin et al. [2, 6]. Comparisons with the measured values of these PMMA phantoms revealed that the small differences identified can be explained using seven different angles of multiple field incident these phantoms [2, 6].

Measured ND𝑙𝑖𝑣𝑒𝑟 differed among different points of the liver inside the Rando and PMMA phantoms. This study measured ND𝑙𝑖𝑣𝑒𝑟 at the 22nd slice of the phantoms. The averaged ND𝑙𝑖𝑣𝑒𝑟 was 1.05 ± 0.20, and ranged from 1.24 ± 0.21 to 0.81 ± 0.19 mSv/Gy of PMMA phantoms. The ND𝑡𝑢𝑚𝑜𝑟 of Rando was 2.70 mSv/Gy, 0.135% of 200 cGy prescribed to PTV at the tumor center. Generally, in medical radiation exposure, patients are most worried about the gonads. Gonads are much more sensitive to radiation [7]. Averaged ND𝑔𝑜𝑛𝑎𝑑 obtained in this study was 0.25 ± 0.04 mSv/Gy, and ranged from 0.19 ± 0.03 to 0.32 ± 0.05 mSv/Gy (cf Fig. 2). The ND𝑙𝑢𝑛𝑔 for lung exposed to 15 MV stated large values, up to 1.02 ± 0.13 mSv/Gy of lung cancer herein. 1.23 ± 0.16 (10 kg), 1.12 ± 0.24 (30 kg), 1.08 ± 0.18 (50 kg), 0.98 ± 0.21 (70 kg), 0.96 ± 0.20 (90 kg) mSv/Gy for these PMMA phantoms. Thyroid, red bone marrow, and remainder were comparatively low NDT during IMRT. The averaged PNDEs were 0.67 ± 0.10 and 0.71 ± 0.12 mSv/Gy for male and female of Rando phantoms, respectively. The averaged ND𝑡𝑢𝑚𝑜𝑟 were 2.70 ± 0.40, 2.23 ± 0.38, 1.98 ± 0.36, 1.87 ± 0.32, 1.12 ± 0.20, and 0.92 ± 0.17 mSv/Gy for Rando and 10 to 90 kg PMMA phantoms, respectively. ND𝑙𝑢𝑛𝑔 at the lung center decreased with increasing distance to the tumor center (cf Fig. 3). Therefore, at the end of a 200 cGy lung cancer treatment, Rando’s PNDE was 0.67 ± 0.10 mSv. PNDE from the 10 kg phantom was approximately 1.43 ± 0.21 times more than that of the 90 kg phantom.

4.2Comparison with other experimental data

This is the first to use different body-weight of PMMA to evaluate PNDE herein, and compared the results presented here under lung cancer treatment with organ absorbed ND investigated with bubble detector in water used for 15 MV x-ray linac-based radiotherapy of Siemens Primus, CCH under a 10 × 10 cm2 [1, 6]. The effective doses of fast and thermal neutrons on the water phantom surface were 1.42 mSv/Gy, and thus were approximately 1.21 times higher than that of 70 kg phantom. These above differences can be attributed to different treatment planning, and phantom compositions of PMMA, epoxy-resin and water used for the conversion of absorbed dose.

Table 3

PNDE (mSv/Gy) and estimated risk (%) of secondary cancer coming from a total 72 Gy for Rando and PMMA phantoms

Phantom (kg)PNDE (mSv/Gy)Risk (%)
Rando0.67 ± 0.102.85
100.80 ± 0.128.35
300.74 ± 0.117.73
500.63 ± 0.092.68
700.61 ± 0.092.59
900.56 ± 0.082.38

Figure 3.

Secondary ND at distal edge of lung on the central axis varying distances during IMRT treatment with 15 MV (a) 90 kg, (b) 70 kg irradiation, (c) 70 kg, (d) 50 kg, (e) 30 kg, (f) 10 kg, (g) Rando phantoms.

Secondary ND at distal edge of lung on the central axis varying distances during IMRT treatment with 15 MV (a) 90 kg, (b) 70 kg irradiation, (c) 70 kg, (d) 50 kg, (e) 30 kg, (f) 10 kg, (g) Rando phantoms.

4.3Risk of secondary malignancies

The risk coming from neutron-induced malignancies is evaluated based on the differences between in the effective dose of the PMMA and Rando phantoms calculated by the BEIR VII [14]. BEIR VII stated that radiotherapy patients’ risk should be approximately 5.9% per Sv for lung cancer patients given a large population spanning all ages and 14.5% per Sv for 10 year old male [4, 14]. The prescription used for phantoms was 72 Gy of the PTV during IMRT. Comparing calculated PNDE of five PMMA phantoms, Table 3 demonstrates that the risk also increases with decreasing body-weight, as well as with decreasing PNDE. This study is in good agreement with Kry et al. stated the risk following prostate radiotherapy was 2–5% for all linac [15].

5.Conclusion

This feasibility study involves the evaluated PNDE of weight-dose relationship for each patient in the lung cancer treatment based on photoneutron contamination using the paired TLDs approach. The highest NDT was close to the tumor, and decreased with increasing distance from the tumor center. PNDE reduced with increasing phantom weight. Evaluated PNDE ranged from 0.80 ± 0.12 to 0.56 ± 0.08 mSv/Gy for these phantoms. The calculated PNDEs clearly show that the risk decreased with increasing body-weight. Results demonstrated that paired TLDs method exhibits high stability as well as sensitivity while evaluating PNDE.

High ND𝑙𝑢𝑛𝑔 values were up to 1.02 ± 0.13 mSv/Gy of Rando. To prevent unnecessarily extra ND to patients as well as satisfy the ALARA principle are strongly recommended. The in vivoNDT measurements reported in this study reveal a measurable difference in these PMMA phantoms that should be considered in cases involving extra secondary radiation. However, the risk also increases with decreasing body-weight and found the extra PNDE to be negligible.

Acknowledgments

The authors would like to thank senior radiotherapist Sing-Sheng Huang for many fruitful discussions, and radiotherapists of CCH and the staff of NTHU for their efficient support.

Conflict of interest

None to report.

References

[1] 

Hashemi SM, Hashemi-Malayeri B, Raisali G, et al. Measurement of photoneutron dose produced by wedge filters of a high energy linac using polycarbonate films. Journal of Radiation Research. (2008) ; 49: (3): 279-283.

[2] 

Chao JH, Liu WS, Chen CY. Estimation of 41Ar concentrations in the vicinity of a high-energy medical accelerator. Radiatation Measurments. (2007) ; 42: : 1538-1544.

[3] 

National Council on Radiation Protection and Measurements, ((1984) )Neutron contamination from medical accelerators, NCRP-79. Washington, DC

[4] 

Harrison RM. Introduction to dosimetry and risk estimation of second cancer induction following radiotherapy. Radiation Measurements. (2013) ; 57: : 1-8.

[5] 

Lin JP, Chua TC, Lin SY, et al. The measurement of photoneutrons in the vicinity of a Siemens Primus linear accelerator. Applied Radiation and Isotopes. (2001) ; 55: : 315-321.

[6] 

Lin JP, Liu WC, Lin CC, Investigation of photoneutron dose equivalent from high-energy photons in radiotherapy. Applied Radiation and Isotopes. (2007) ; 65: : 599-604.

[7] 

International Commission on Radiological Protection. ((1991) )Recommendation of the ICRP, ICRP Publication 60, Annals of the ICRP, No. 1–3; Pergamon Press, Oxford, UK

[8] 

Ministry of Health and Welfare, ((2019) )Statistics of General Health Welfare of 2019, available at https://dep.mohw.gov.tw/ DOS/mp-113.html

[9] 

Hsu FY, Chiu MC, Chang YL, et al. Estimation of photon and neutron dose distributions in the THOR BNCT treatment room using dual TLD method. Radiation Measurements. (2008) ; 43: (2-6): 1089-1094.

[10] 

International Commission on Radiation Units and Measurement, ((1992) )Phantoms and computational models in therapy, diagnosis and protection ICRU-48. Bethesda

[11] 

Matsumoto T. Monte Carlo simulation of depth-dose distribution in several organic models for boron neutron capture therapy. Nuclear Instruments and Methods in Physics Research, A. (2007) ; 508: : 552-557.

[12] 

Lin HC, Lai TJ, Tseng HC, et al. Radiation doses with various body weights of phantom in brain 128-slice MDCT Examination. J Radiation Research. (2019) ; 60: (4): 466-475.

[13] 

International Commission on Radiological Protection, ((2007) )Recommendations of the International Commission on Radiological Protection, ICRP-103, Pergamon Press, Oxford

[14] 

Health Risks from Exposure to Low Levels of Ionizing Radiation, ((2006) ) BEIR VII-Phase 2. Committee to Assess Health Risks from Exposure to Low Levels of Ionizing Radiation. National Research Council The National Academies Press, Washington DC.

[15] 

Kry SF, Salehpour M, Followill DS, et al. The calculated risk of fatal secondary malignancies from intensity-modulated radiation therapy. International Journal of Radiation Oncology Biology Physics. (2005) ; 62: (4): 1195-1203.