Research Article  Open Access
Energy Modulated Photon Radiotherapy: A Monte Carlo Feasibility Study
Abstract
A novel treatment modality termed energy modulated photon radiotherapy (EMXRT) was investigated. The first step of EMXRT was to determine beam energy for each gantry angle/anatomy configuration from a pool of photon energy beams (2 to 10 MV) with a newly developed energy selector. An inverse planning system using gradient search algorithm was then employed to optimize photon beam intensity of various beam energies based on presimulated Monte Carlo pencil beam dose distributions in patient anatomy. Finally, 3D dose distributions in six patients of different tumor sites were simulated with Monte Carlo method and compared between EMXRT plans and clinical IMRT plans. Compared to current IMRT technique, the proposed EMXRT method could offer a better paradigm for the radiotherapy of lung cancers and pediatric brain tumors in terms of normal tissue sparing and integral dose. For prostate, head and neck, spine, and thyroid lesions, the EMXRT plans were generally comparable to the IMRT plans. Our feasibility study indicated that lower energy (<6 MV) photon beams could be considered in modern radiotherapy treatment planning to achieve a more personalized care for individual patient with dosimetric gains.
1. Introduction
High energy megavoltage photon beams (10 MV and higher) are suitable to treat deepseated tumors due to their penetration power and skin sparing. However, the large exit dose outside the target volume and the neutron contamination can be of great concerns in clinical practice. In addition, high energy photons have a pronounced dose buildup or builddown effect at the heterogeneous medium interface due to the increased range of secondary electrons, causing increase of beam penumbra, overdose of normal tissues, and underdose of target edge volumes [1–3]. On the other hand, low energy (6 MV and lower) photon beams have narrow penumbra which are suitable to deliver a tight dose distribution around the target. Low exit dose also benefits the critical structures adjacent to the target.
In recent years, there has been increasing interest in exploring the potential benefits of lower energy photon beams, such as Cobalt60 gamma ray (energy peaks at 1.17 and 1.33 MeV) and lower energy (<6 MV) Xrays. In 2007, Keller et al. showed that intermediate energy Xrays (1.2 MV) combined with small fields can reduce the radiological penumbra in intracranial stereotactic radiosurgery (SRS), which could be substantially beneficial for improved dose distribution homogeneity and better sparing of critical structures [4]. Fox et al. have compared Cobalt60 gamma ray with 6/18 MV photons and demonstrated that nearly identical intensitymodulated radiotherapy (IMRT) plans can be achieved between Co60 and 6 MV photons [5]. Dhanesar et al. have indicated that Co60 Tomotherapy is capable of providing the stateoftheart conformal dose delivery in both phantom and clinical planning studies [6]. If the buildup effect is of concern when the target volume is near air cavities, Behrens suggested that 4 MV should be preferred over both 6 MV and 8 MV [7]. More recently, Dong et al. investigated a 2 MV FFF beam for extracranial robotic IMRT. Their results demonstrated that the dual energy plan (2 and 6 MV) had the best dosimetry in terms of equivalent target coverage and improved organsatrisk (OARs) sparing, followed by 2 MV only and 6 MV only plans [8]. In addition, Zhang et al. studied the effects of 3 MV photon beams for lung cancer treatment [9]. By comparing 6 MV, 3 MV, and dual energy plans for 31 lung cancer patients, they concluded that 3 MV photon beams have statistically significant dosimetric benefits in terms of improved tumor coverage and reduced doses to the adjacent critical structures [9].
The potential of mixing photon beam energy has also been explored in some early investigations. Malhotra et al. showed that, for the prostate cases, the mixed energy photons of 6, 10, 15, or 18 MV would be better in terms of integral dose than using either low or high energy photons alone [10]. In a study by StHilaire et al., beam energy was added as an optimization parameter in an automatic aperturebased inverse planning system. Their work demonstrated that energy optimization could produce plans of better quality with less peripheral dose and fewer monitor units (MUs) for prostate and lung tumors [11]. Park et al. also showed the advantage of energy and intensitymodulated radiotherapy in terms of less integral dose while maintaining the plan quality [12]. In their other study, a cylindrical energy modulator with adjustable thickness of mercury along the beam axis was used to replace the flattening filter and to modulate the photon beam energy. They further commissioned a virtual machine in a treatment planning system (TPS) based on Monte Carlo simulated photon beams going through the energy modulator of various thickness of mercury [13]. While this novel energy modulator provided a way to generate mean energy adjustable photon beams, the incorporation of this modulator with toxic mercury into a modern Linac treatment head would pose a series of engineering, manufacturing, and safety challenges.
Based on all the previous observations, we propose a new treatment modality termed energy modulated photon radiotherapy (EMXRT) to systematically explore the benefits of lower energy photon beams for cancer radiotherapy. In EMXRT, a broad range of lower energy photons (i.e., 2, 3, 4, 5, 6, 10 MV) were considered as variables during the inverse planning process, with the consideration of tumor size, effective path length inside the patient for a certain gantry angle, and skin dose as well as beam exit dose. Higher energy photons (>10 MV) were not chosen in this study to avoid the overlapping with previous investigations [10, 12]. All the energy beams were modeled using an EGS4/BEAM Monte Carlo code based on the specifications of an actual linear accelerator. An inhouse inverse planning system based on gradient search algorithm has been developed and used to optimize photon beam intensity based on presimulated Monte Carlo pencil beam dose distributions in patient anatomy. The pros and cons of EMXRT have been discussed in six clinical cases with respect to the conventional IMRT treatments.
2. Material and Methods
2.1. Patient Selection and Workflow
Six cancer patients treated with IMRT were retrospectively selected to represent a variety of tumors (i.e., prostate, brain, lung, spine, thyroid, and head and neck) and locations. As shown in Figure 1, a patient CT phantom was first created from planning CT images, with the target and other critical structures delineated in a commercial treatment planning system (TPS). For each patient, the original IMRT treatment was planned with a Varian Eclipse TPS v10.0.28 on a Varian Trilogy Linac (Varian Oncology Systems, Palo Alto, CA). Accordingly, an EMXRT treatment was planned with our inhouse inverse planning system based on the presimulated Monte Carlo pencil beam dose distributions of a variety of lower energy photon beams. For the sake of comparison, we have also replanned the IMRT using the same inhouse inverse planning system with the same beam energy set as in the original Eclipse plan. Once the inverse planning was done, Monte Carlo method was used to simulate the dose distributions in patient anatomy for both IMRT and EMXRT plans. Specifically, a benchmarked EGS4 user code, MCSIM, was used to simulate the 3D dose distributions in the obtained patient CT phantoms with all the beam configurations, leaf leakage, and transmission being considered [14, 15]. Finally, the two plans were evaluated and compared in terms of target coverage and normal tissue sparing.
2.2. Monte Carlo Simulation and Modeling
It is widely accepted that the Monte Carlo method is one of the most accurate dose calculation methods available for radiotherapy. In this study, an EGS4/BEAM code has been used to simulate the particles emanating from a Varian Linac treatment head. The machine configurations were the same as those of the clinical 6 MV photon beams; only the energy of incident electron beam was set to be the nominal energy of 2, 3, 4, and 5 MeV, respectively [16, 17]. The configuration of 10 MV was different as the flattening filter for 10 MV was used. The phase space data was first scored above the photon jaws located at 28 cm downstream from the target and then EGS4/BEAMDP was used to derive the multiple source model for all the energy beams. The obtained multiple source models consisted of detailed numerical description of the energy spectrum, spatial distribution, fluence distribution, source location, shape, and size of each source for a particular treatment head [18, 19]. The multiple source models have been shown equivalent to the phase space data in representing the photon beams from the Linac treatment head and replicating the dose distributions in water, yet eliminating the inconvenience of large data transfer and latent variance related to the phase space [18–20].
EGS4/MCSIM has been used for 3D dose calculations in both water phantoms and patient CT anatomy in this study. To demonstrate the beam characteristics of low energy photons, a series of percent depth dose curves (PDDs) and dose profiles along axis at various depths has been simulated with MCSIM in a rectilinear water phantom. In addition, the PDDs and dose profiles of 6 and 10 MV photon beams in water were measured and compared with Monte Carlo simulations. To simulate 3D dose distributions in patient anatomy, an inhouse DICOM tool named DICOMan was used to construct a realistic CT phantom based on the planning CT images and contoured structures of individual patient exported from Eclipse TPS via DICOM RT protocol [21]. To enhance the efficiency of Monte Carlo simulations, a series of variance reduction techniques has been implemented in MCSIM such as photon interaction forcing, Russian roulette, particle splitting, and electron track repeating. In all Monte Carlo simulations, the energy cutoffs for electrons (ECUT) and photons (PCUT) and the energy thresholds for δray production (AE) and bremsstrahlung production (AP) were set as ECUT = AE = 700 keV and PCUT = AP = 10 keV, respectively. A better than 1% statistical uncertainty (1σ) was achieved for beam PDDs, profiles, and the target voxels in dose simulations. The benchmark results of EGS4/MCSIM have been reported previously [14, 15].
2.3. Energy Modulated Photon Radiotherapy Treatment Planning
2.3.1. Energy Selector
For EMXRT treatment planning, an energy selector was developed to provide a convenient way to determine an appropriate energy for photon beams at different gantry angles in the studied energy range. As shown in Figure 2, for any photon beam impinging on the patient contour from outside and transiting the tumor volume inside, there were four special points to be considered by the energy selector, that is, the entrance point to the contour (point 1, corresponding to skin dose), the entrance point to the target (point 2, target entrance dose) and exit point from the target (point 3, target exit dose), and finally the exit point from the contour (point 4, exit dose). The skin dose is one of the limiting factors during the treatment planning, especially when lower energies were considered. In addition, the effective path length (EP) and tumor size (TS) can be determined based on the distances between points 1, 2, and 3. The target dose homogeneity can be represented by the target entrance and exit doses. The exit dose at point 4 can be a general representative of the doses to the critical structures and healthy tissues downstream.
Ideally, a uniform dose distribution in the target with minimum doses delivered to the OARs would be expected. However, when treatments are delivered to the patients, the radiation beam may intersect different mediums such as bone, muscle, fat, lung, and air cavities. Because of the inhomogeneity of human body, we used the effective path length (EP) in this study, which was obtained by normalizing the physical path length to a waterequivalent path length from patient skin to the isocenter. Tumor size (TS) was calculated as the waterequivalent length from the target entrance point to the target exit point. A wide range of tumor size from 2 to 16 cm was investigated in this study to accommodate most of the clinical cases.
A spherical tumor model merged in water, whose diameter was equal to tumor size (TS) and whose center was effective path length (EP) away from the entrance, was used to imitate various clinical situations where impinging beams intersected four points along the beam direction with the contour and the tumor. Correlation coefficients (CC), as defined in (1), were calculated between the ideal doses and the realistic doses of different energies for these four points. Particularly, for the ideal situation, and were set to be 10% while = = 100% for the ideal homogeneous dose distribution in the tumor. For the realistic situation, were dose values at four points derived from PDD curves for a given beam energy, TS, and EP in water. Weighting factor was assigned to each point to show different importance when considering photon energy for a radiation treatment. According to our experience, we assigned 0.5 to skin dose but 1.0 to the other three points in this study. CC was calculated as below:
The energy corresponding to the maximum CC was chosen as the optimal energy for this EP/TS scenario. Following the same fashion, an energy lookup table has been generated for various EP/TS configurations, which can be readily used in energy selection for different patient cases. It should be pointed out that some common knowledge was also applied in the implementation of energy selector. For example, for some deepseated tumors with very large EP (>18 cm), regardless of tumor size, high energy (10 MV) photon beam would be chosen. Moreover, EP should always be larger than the radius of tumor (1/2 of TS).
2.3.2. Inverse Planning
With multiple source models as beam input, the dose deposition coefficients (DDCs) were scored to yield the dose distribution from each beamlet (1 cm × 1 cm) to each structure voxel in the patient for the specific beam energy at specific gantry angle. Therefore, for the dose to the voxel , is given by the summation of the DDCs from all beamlets weighted by their own weights as shown inwhere is the total number of beamlets for all beams, represents the DDC from the th beamlet to the th voxel, and is the weigh for the th beamlet.
Then, the weighting factors for all beamlets were optimized via an inhouse inverse planning code so that the desired dose distribution could be achieved. Further details of the optimization method used in this work can be found elsewhere [22]. The system utilizes a steepest descent search algorithm, with a quadratic objective function augmented by dose volume constraints of both target and OARs, as shown below:where , , and are the prescribed dose to the target, the th dose volume constraints for the target, and the th dose volume constraints for the th critical structure. and are the total number of dose points inside the target and the th critical structure, respectively. , , and are the weights assigned to the th dose volume constraints for target, the th critical structure, and the th dose volume constraint for the th critical structure, respectively. is a flag, defined as 1 when the constraint is violated and 0 when it is not. is the number of critical structures considered in the optimization and is the number of dose volume constraints for the th critical structure. The factor will be increased as the iterations proceed. Thus, the penalty assigned to the constraint violations becomes increasingly severe as represented by . In each iteration, the weights of beamlets were adjusted along the negative of the gradient to minimize the objective function. The optimization will stop once the value of the objective function keeps almost unchanged for several rounds or the predefined time limit is violated.
To simplify the optimization, the complete irradiated area outlines (CIAO) extracted from the IMRT plan in Eclipse, which were the general beam shape opened for fitting the target, were used as the initial guess such that the intensity of the beamlets inside the corresponding CIAO was set to be one while the outside was set to be zero.
2.4. Plan Evaluation
Dose volumetric analysis was performed to compare the new modality EMXRT with the conventional IMRT. To start with, the dose volume histograms (DVHs) of both IMRT and EMXRT plans were normalized such that 100% isodose lines would encompass 95% of target volume. For PTVs, the mean dose () and the dose levels covering 2% and 98% of the PTV volume ( and ) were calculated and compared. In addition, the homogeneity index (HI), defined as HI = ()/, was used to evaluate the homogeneity of the PTV coverage. Some specifically selected dosimetric parameters of OARs were generated from DVHs and analyzed. The integral dose was also calculated for surrounding normal tissue by integrating the doses over all the voxels within the patient volume while excluding the PTV. Finally, the isodose distributions in the axial, coronal, and sagittal views were compared between the two modalities.
3. Results
3.1. Benchmark Results
Figure 3 compared the percent depth doses of various beam energies in water between the Monte Carlo simulations (line with symbols) and the measurements (symbols) for the field of 10 cm × 10 cm at 100 cm sourcetosurface distance (SSD). Note that all the Monte Carlo simulations were performed with MCSIM with multiple source models as beam input, and the depth dose measurements were carried out in an IBA blue phantom water scanning system and scanned with two CC13 ionization chambers (IBA Dosimetry America, Bartlett, TN) with one as field detector and the other as reference detector. The agreement between the measurements and the Monte Carlo results for 6 and 10 MV have been validated in our previous work regarding multiple source modeling for accurate Monte Carlo dose calculations [20]. Part of the percentage depth dose data and the surface dose for all the beams were listed in Tables 1 and 2, respectively. As beam energy decreased, the TPR_{20,10} decreased from 0.627 for 10 MV to 0.446 for 2 MV. Although lower energy photons carried less penetrating power as indicated by the smaller TPR_{20,10}, it could be advantageous in some situations where critical structures were to be spared in the downstream.

 
a, b, c, and d are percentage doses at depths of 0, 10, 20, and 30 cm, respectively. e is the tissue phantom ratio calculated as . 
Beam profiles for all the energies were also generated from Monte Carlo simulations at 10 cm depth for the 10 cm × 10 cm (Figure 4(a)) and 40 cm × 40 cm (Figure 4(b)) field sizes. For 10 cm × 10 cm open field, the overall flatness of the beam profiles was considered clinically acceptable, while, for the 40 cm × 40 cm open field, horn shaped profiles were observed for the lower energies, such as 2 MV, 3 MV, and 4 MV. This uneven flatness was primarily caused by the nonoptimized design of flattening filter for these energies used in Monte Carlo beam simulation as well as the inherent fluctuation in the particle fluence scored in the phase space for beam modeling. However, in this work, the uneven fluence will be handled by assigning the appropriate weighting factors for the beamlet dose distributions during the inverse planning process.
(a)
(b)
3.2. Energy Selector
In this work, the correlation coefficients (CC) were calculated for all the available beams for each scenario with EP ranging from 2 to 21 cm and TS from 2 to 16 cm. The energy with the maximum CC was chosen as the estimated energy for the specific EP/TS configuration. An energy lookup table shown in Table 3 was determined by the proposed energy selector. Overall, higher energy photons were preferred for largesized tumors and/or deepseated ones, consistent with conventional clinical practice. 2 MV was not indicated in this energy table due to its very high skin dose and steep PDD curve. 5 MV was preferred in a large number of scenarios when the skin dose, the exit dose, and the tumor dose homogeneity were considered. For tumors larger than 11 cm in diameter or located as deep as 18 cm, the suitable photon energy was 10 MV as determined by the energy selector, which is consistent with our expectation.

For all six cases, the EP and TS of all beam angles were measured in Eclipse and listed in Table 4. In IMRT plans, 10 MV photon beams were used for the prostate patient and 6 MV photon beams for all the other patients. In EMXRT plans, beam energies were determined according to the energy lookup table in Table 3. Compared with IMRT plans, more beams of lower energy and multienergy photon beams were employed in EMXRT plans. Specifically, in the prostate EMXRT plan, 10 MV was selected for the left and right posterior oblique beams while 5 MV was preferred for the anterior posterior and two anterior oblique beams. For the brain case, even the effective path lengths for all the beam angles were relatively small, 5 MV instead of lower energy photon beams were determined by the energy selector as the large volume of tumor would require higher energy photons with stronger penetration power. For the other four cases, as tumor sizes largely fell into the range of 4 to 8 cm, 5 MV photon beams were dominantly chosen by the energy selector based on the consideration of tumor dose coverage and sparing of critical structures adjacent to the tumor.
 
is the volume of PTV. is the diameter of equivalent sphere of PTV volume. is the prescription dose for PTV. is effective path length. is tumor size. 
3.3. Plan Evaluation
The dose distributions along the axial, coronal, and sagittal planes for the lung and the prostate cases were shown in Figures 5 and 6 with 95%, 80%, 60%, 40%, and 20% isodose distributions represented by the red, yellow, green, light blue, and dark blue lines, respectively. As shown in Figure 5, the 95% red isodose line in EMXRT conformed to the tumor volume (red color organ) better than that in IMRT. In addition, the 40% and 20% isodose lines covered less volume in EMXRT. On the other hand, the IMRT plan delivered a better dose distribution in the prostate case than the EMXRT plan in terms of target dose conformity and normal tissue sparing, as shown in Figure 6.
(a)
(b)
(a)
(b)
Tables 5 and 6 compared the dosimetric parameters of PTV and OARs in IMRT plans with those in EMXRT plans for all the six cases. As mentioned early, all plans were normalized such that the 100% isodose lines encompassed 95% of PTV. Overall, the PTV coverage in both EMXRT and IMRT plans showed no big differences, with up to 1% difference in terms of and up to 1.7% difference on of PTV in the spine case. Compared to IMRT plans, the EMXRT plans showed the same homogeneity of PTV dose distributions indicated by HI values for the prostate, the head and neck, and thyroid cases. For the lung and brain cases, EMXRT plans yielded better PTV dose homogeneity as compared to IMRT, up by 5.9% and 7.1%, respectively. However, the PTV dose homogeneity was 10.5% worse in EMXRT plan of spine lesion when compared with IMRT.
 
Relative difference = 100% × ( − )/. , , and are the mean, maximum, and minimum doses of PTV, respectively. and are the dose levels covering 98% and 2% of the PTV volume, respectively. HI = ( − )/. 
 
Relative difference = 100% × ( − )/. The “+” and “−” indicate whether the dose index of the EMXRT plan is larger or smaller than that of IMRT plan. is the mean dose, and is the maximum dose. 
For the prostate case, rectum and bladder received equivalent doses in EMXRT compared with IMRT, in terms of mean dose and , with up to 3.1% and 1% relative differences, respectively (Table 6). The mean dose to the femur head increased by 29.2% and the was up by 42.2% in EMXRT plan, although both indexes were still within the tolerances.
In the head and neck case, all critical structures showed comparable DVHs as indicated in Figure 7(b). Some OARs such as spinal cord, brainstem, and mandible were slightly more spared in EMXRT while, for left and right parotids, the mean doses were 11.2% and 7.1% higher in EMXRT than in IMRT.
(a)
(b)
(c)
(d)
(e)
(f)
In terms of lung cancer radiotherapy, our study indicated that EMXRT was a much better modality than IMRT, not only on the target coverage (Table 5), but also on the OAR sparing (Table 6). The mean dose and maximum dose of the spinal cord were 18.2% and 11.3% lower in EMXRT than in IMRT. The right and left lungs were much more spared with the mean doses decreased from 1.7 Gy and 3.5 Gy to 1.1 Gy and 2.6 Gy, a 35.3% and 25.7% reduction, respectively.
The spinal cord dose in the spine case was a big challenge due to its proximity to the target. In EMXRT plan, with almost the same target coverage, the mean dose to the spinal cord was reduced by 4.1% and the maximum dose got a reduction of 8.2%. Liver in this case also received 7.7% less dose on average as compared to IMRT. However, the target dose coverage was 10.5% less homogeneous in EMXRT plan than in IMRT plan, due to the 5 MV used in the majority of beams for a lesion wrapping around the bony structure.
The brain case in this study was a pediatric case, so any dose reduction on the OARs can be of great significance. Replacing the 6 MV photons in IMRT with the 5 MV photons in EMXRT, it was found that the OARs got better spared yet the PTV dose coverage was 7.1% more homogeneous in EMXRT. Specifically, the mean dose to the brainstem and the maximum doses to the eyes and chiasm were 6.3%, 12.9%, and 7.7% lower in EMXRT than in IMRT. In combination with better target dose coverage, our results suggested that EMXRT with 5 MV photon beams was a better treatment modality for the pediatric brain tumors than the conventional IMRT with 6 MV photon beams.
For the case of thyroid cancer, except for the 21.1% higher mean dose to the pharynx, all the other OARs such as the spinal cord, the esophagus, and the trachea were more spared in EMXRT than in IMRT by 1.7% to 9%. The PTV coverage for the two plans was basically the same.
The integral dose defined as the product of absorbed dose with the volume of normal tissues was calculated and compared between the two modalities. Figure 8 showed the relative difference of integral dose for the six cases. It was obvious that EMXRT outperformed IMRT in the lung, the brain, and the thyroid cases with the reduction of integral dose by 8.8%, 2.3%, and 4.0%, respectively. As to the spine case, both modalities delivered the same amount of integral dose to the total healthy tissues. For the prostate and head and neck cases, IMRT delivered 4.5% and 2.2% less integral doses to the normal tissues than EMXRT did.
4. Discussion
In this work, the impact of using mixed lower photon beam for IMRT on six various sites is investigated. Compared to the clinical IMRT energy settings, the proposed mixed lower photon energies plans showed improved target coverage and homogeneity, better sparing of the OARs, and reduced integral dose for the lung cancer and pediatric brain tumor patient. For other lesions such as prostate, spine, thyroid, and head and neck, the lower energy beams were less beneficial than the clinical energy settings.
Beam energies were determined easily by our energy selector according to tumor size and effective path length along beam direction. The lower energy photon beams (<6 MV) help to reduce the doses deposited to the adjacent critical structures in lung cancer EMXRT plan due to the rapid dose falloff with smaller penumbra at the field edge. The brain case investigated in this work was a pediatric case. Cranial radiation therapy has been associated with the highest risk of longterm cognitive morbidity particularly in younger children. There has been an established doseresponse relationship correlating the high OAR dose in cranial radiotherapy with the poor performance on intellectual measures in the long term [23]. The improvements on normal tissue sparing and integral dose reduction are meaningful. Further study will be carried out to explore the role of lower energy photon beams in pediatric cancer radiotherapy.
To verify the primary conclusion of the potential dosimetric benefits of EMXRT for the lung and pediatric brain patients, one more case for each kind of cancer was analyzed and compared. Following the same procedure mentioned above, the energy settings for the EMXRT plans were 4, 5, 5, 3, 3, 4, and 5 MV and 5, 5, 4, 5, 4, 4, and 4 MV for the additional lung and brain cases, respectively, while the corresponding IMRT plans used 6 MV for all beams (Table 7).
 
EP is effective path length. TS is tumor size. 
As shown in Figure 9(a), the EMXRT plan showed better plan quality in the second lung case with steeper PTV dose curve and lower OARs dose curves. The PTV dose was more homogeneous in the EMXRT plan with slightly lower HI value (0.22) compared to the IMRT plan (0.23). The of the whole lung tissue, that is, the percent volume of lung receiving 5 Gy, was decreased from 18.8% to 16.9%, with 10.1% reduction in EMXRT compared to that in IMRT, which means more lung tissues were spared by the low dose irradiation. For the EMXRT plan, the mean dose of the spinal cord, esophagus, great vessel, and trachea was decreased by 22.2%, 10.7%, 5.4%, and 4.1%, respectively (Table 8). The integral dose for the second lung case was decreased from 3.02 Gy to 2.73 Gy, with 9.7% reduction for the EMXRT plan. Overall, our conclusion was in agreement with the previous works [8, 9]; that is, the EMXRT and the application of lower energy photons were able to improve the target dose distribution and reduce the OARs dose and the integral dose for the lung cancer treatment.
 
Relative difference = 100% × ( − )/. is the mean doses of each organ/tissue. HI = ( − )/. , , and are the dose levels covering 98%, 2%, and 50% of the PTV volume, respectively. and are the percent of whole lung volume receiving 5 Gy and 20 Gy dose, respectively. 
(a)
(b)
For the second brain case, the PTV dose distribution was identical for both IMRT and EMXRT with same PTV mean dose and HI value (Figure 9(b) and Table 8). The mean dose to the brainstem was decreased from 9.8 Gy to 9.2 Gy, with a 5.7% reduction in EMXRT compared to that in IMRT. The left and right optic nerves received slightly less doses in EMXRT (2.8% and 2.4% reduction, resp.). The integral dose was 4.6% lower in EMXRT plan (1.44 Gy) than in IMRT (1.51 Gy). Our results suggest that EMXRT has potential dosimetric benefits for the pediatric brain cancer patients.
Although there are currently no technologies available to modulate photon energy, fast energy switching is possible as seen in some of the modern linear accelerators. For example, clinically it only takes less than 10 seconds for the Varian Trilogy to switch energy between 6 MV and 10 MV, making it possible to change energy for different gantry angles during EMXRT treatment. Furthermore, theoretically speaking, the EMXRT modality should provide increased flexibility as photon energy can be treated as a new variable in the planning optimization. Of course, the additional flexibility may increase the difficulty in the optimization and in quality assurance and quality control. Nevertheless, this study has demonstrated the feasibility of EMXRT and clinical benefits of lower energy (<6 MV) photon beams in the lung cancer radiotherapy and pediatric patients, opening a door for a more personalized treatment planning for individual patient.
As mentioned above, there are a growing number of studies on the potential advantages and applications of low energy photon beams (<6 MV) [6, 8, 9]. These energies (<6 MV) were especially considered in this new modality treatment study [24]. With advanced intensity modulation technology, it has been shown that radiotherapy has become less restricted by the weak penetration power of lower energy [25, 26]. The feasibility and potential benefits of EMXRT in the energy range of 2–10 MV have been demonstrated in this study. Higher energy photons (> 10 MV) for EMXRT have not been investigated so far, due primarily to their inherent physical characteristics, such as high exit dose to the critical structures downstream, longer secondary electron range, and neutron contamination.
To study the effect of higher energy photons for EMXRT, we compared a HEEMXRT (higher energy EMXRT) plan with the LEEMXRT (lower energy EMXRT) and IMRT plans of prostate cancer. In the LEEMXRT plan, 10 MV was selected for the left and right posterior oblique beams while 5 MV was used for the anterior posterior and two anterior oblique beams. In the HEEMXRT plan, the two 10 MV beams were replaced with 15 MV photons. Identical beam configuration and optimization constraints were used in all the three plans. As shown in Figure 10, the target dose distributions were identical for all the plans. However, the rectum received more irradiation in HEEMXRT plan with mean dose increased by 4.7% compared to that in IMRT plan. The mean dose to the femur head increased by 41.9% in HEEMXRT plan (13.7 Gy) compared to IMRT plan (10.6 Gy). The increased OAR doses were mostly due to the increased secondary electron range and penetrating power of the 15 MV photons. Our results indicated that, with current energy modulation method, it is not necessary to involve the photon energies higher than 10 MV. While it is possible to further improve the HEEMXRT plan by adjusting the optimization constraints, we kept them identical in this study for the sake of comparison.
In addition, it has been demonstrated that lower energy megavoltage photon beams can improve the image quality. The MV fan beam CT (MVCT) with effective energy of 3.5 MV from a Helical Tomotherapy unit has been shown to provide sufficient contrast for softtissue delineation [27, 28]. Several studies have further shown that, with lowZ targets in linear accelerators producing photon beams as low as 1.9 MV, the image quality could be greatly enhanced compared to 6 MV photons [29, 30]. Hence, photon beams of lower energy (<6 MV) could be potentially useful for the target localization and radiation delivery in the imageguided IMRT treatments. Furthermore, with the actual treatment beam being the source for MV portal or MVCT image acquisition, the imaging isocenter is identical to the treatment isocenter, eliminating the need to perform a comprehensive calibration and registration procedure usually required if a kilovoltage imaging system is installed perpendicular to the treatment beam axis [24]. The expanded beam energy range will not only increase the flexibility for personalized radiotherapy, but also provide more options for personalized imaging protocol [31].
5. Conclusion
In this study, we have proposed a new treatment modality termed energy modulated photon radiotherapy (EMXRT), which adaptively set the beam energy from a wide range of photon beams (2 to 10 MV) based on patient anatomy and beam angles, and demonstrated its feasibility in six clinical cases by using a Monte Carlo based inverse planning platform. Compared to current IMRT technique, the proposed EMXRT offered much improved OAR sparing in the radiotherapy of lung cancers and pediatric cases. For the prostate, head and neck, spine, and thyroid cases, the EMXRT plans achieved comparable quality to the IMRT plans. Our pilot study suggested that photon beam energy could be added as a new variable in treatment planning for a more personalized radiation treatment and lower energy photons (<6 MV) have the potential benefits in treating lung cancer and pediatric patients. The proposed EMXRT method indicated that photon beam energy may be considered as a technical parameter in current treatment approaches to increase the flexibility of radiotherapy and provide more possibility to perform personalized radiotherapy.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 I. J. Das and K. R. Kase, “Higher energy: is it necessary, is it worth the cost for radiation oncology?” Medical Physics, vol. 19, no. 4, pp. 917–925, 1992. View at: Publisher Site  Google Scholar
 B. L. Werner, I. J. Das, F. M. Khan, and A. S. Meigooni, “Dose perturbations at interfaces in photon beams,” Medical Physics, vol. 14, no. 4, pp. 585–595, 1987. View at: Publisher Site  Google Scholar
 P. J. White, R. D. Zwicker, and D. T. Huang, “Comparison of dose homogeneity effects due to electron equilibrium loss in lung for 6 MV and 18 MV photons,” International Journal of Radiation Oncology Biology Physics, vol. 34, no. 5, pp. 1141–1146, 1996. View at: Publisher Site  Google Scholar
 B. M. Keller, D. J. Beachey, and J.P. Pignol, “Experimental measurement of radiological penumbra associated with intermediate energy xrays (1 MV) and small radiosurgery field sizes,” Medical Physics, vol. 34, no. 10, pp. 3996–4002, 2007. View at: Publisher Site  Google Scholar
 C. Fox, H. E. Romeijn, B. Lynch, C. Men, D. M. Aleman, and J. F. Dempsey, “Comparative analysis of ^{60}Co intensitymodulated radiation therapy,” Physics in Medicine and Biology, vol. 53, no. 12, pp. 3175–3188, 2008. View at: Publisher Site  Google Scholar
 S. Dhanesar, J. Darko, C. P. Joshi, A. Kerr, and L. John Schreiner, “Cobalt60 tomotherapy: clinical treatment planning and phantom dose delivery studies,” Medical Physics, vol. 40, no. 8, Article ID 081710, 2013. View at: Publisher Site  Google Scholar
 C. F. Behrens, “Dose buildup behind air cavities for Co60, 4, 6 and 8 MV. Measurements and Monte Carlo simulations,” Physics in Medicine and Biology, vol. 51, no. 22, pp. 5937–5950, 2006. View at: Publisher Site  Google Scholar
 P. Dong, V. Yu, D. Nguyen et al., “Feasibility of using intermediate xray energies for highly conformal extracranial radiotherapy,” Medical Physics, vol. 41, no. 4, Article ID 041709, 2014. View at: Publisher Site  Google Scholar
 Y. Zhang, Y. Feng, M. Ahmad et al., “Intermediate megavoltage photon beams for improved lung cancer treatments,” PLoS ONE, vol. 10, no. 12, Article ID e0145117, 2015. View at: Publisher Site  Google Scholar
 H. Malhotra, J. Avadhani, S. de Boer, and M. Podgorsak, “SUFFT81: mixed beam energy for IMRT treatment of prostate carcinoma,” Medical Physics, vol. 32, no. 6, p. 1968, 2005. View at: Publisher Site  Google Scholar
 J. StHilaire, C. Sévigny, F. Beaulieu, L. Gingras, D. Tremblay, and L. Beaulieu, “Optimization of photon beam energy in aperturebased inverse planning,” Journal of Applied Clinical Medical Physics, vol. 10, no. 4, p. 3012, 2009. View at: Google Scholar
 J. M. Park, C. H. Choi, S. W. Ha, and S.J. Ye, “The dosimetric effect of mixedenergy IMRT plans for prostate cancer,” Journal of Applied Clinical Medical Physics, vol. 12, no. 4, p. 3563, 2011. View at: Google Scholar
 J. M. Park, J.I. Kim, C. Heon Choi, E. K. Chie, I. H. Kim, and S.J. Ye, “Photon energymodulated radiotherapy: Monte Carlo simulation and treatment planning study,” Medical Physics, vol. 39, no. 3, pp. 1265–1277, 2012. View at: Publisher Site  Google Scholar
 J. S. Li, T. Pawlicki, J. Deng, S. B. Jiang, E. Mok, and C.M. Ma, “Validation of a Monte Carlo dose calculation tool for radiotherapy treatment planning,” Physics in Medicine and Biology, vol. 45, no. 10, pp. 2969–2985, 2000. View at: Publisher Site  Google Scholar
 C.M. Ma, J. S. Li, T. Pawlicki et al., “A Monte Carlo dose calculation tool for radiotherapy treatment planning,” Physics in Medicine and Biology, vol. 47, no. 10, pp. 1671–1689, 2002. View at: Publisher Site  Google Scholar
 W. R. Nelson, H. Hirayama, and D. W. O. Rogers, “EGS4 code system,” Stanford Linear Accelerator Center Report SLAC265, Stanford Linear Accelerator Center, Menlo Park,, Calif, USA, 1985. View at: Google Scholar
 D. W. O. Rogers, B. A. Faddegon, G. X. Ding, C.M. Ma, J. We, and T. R. Mackie, “BEAM: a Monte Carlo code to simulate radiotherapy treatment units,” Medical Physics, vol. 22, no. 5, pp. 503–524, 1995. View at: Publisher Site  Google Scholar
 C.M. Ma and D. W. O. Rogers, “Beam characterization: a multiplesource model,” NRC Report PIRS 509d, 1995. View at: Google Scholar
 C.M. Ma, “Characterization of computer simulated radiotherapy beams for MonteCarlo treatment planning,” Radiation Physics and Chemistry, vol. 53, no. 3, pp. 329–344, 1998. View at: Publisher Site  Google Scholar
 J. Deng, S. B. Jiang, A. Kapur, J. Li, T. Pawlicki, and C.M. Ma, “Photon beam characterization and modelling for Monte Carlo treatment planning,” Physics in Medicine and Biology, vol. 45, no. 2, pp. 411–427, 2000. View at: Publisher Site  Google Scholar
 DICOMan, “The Universal DICOM Wizard for Radiation Oncology and Radiology,” Version 3.4.4, 2011, http://radonc.uams.edu/research/dicoman/. View at: Google Scholar
 S. B. Jiang, Development of a compensatorbased intensitymodulated radiation therapy system [Ph.D. thesis], Medical College of Ohio at Toledo, Ann Arbor, Mich, USA, 1998.
 T. S. Reimers, S. Ehrenfels, E. L. Mortensen et al., “Cognitive deficits in longterm survivors of childhood brain tumors: identification of predictive factors,” Medical and Pediatric Oncology, vol. 40, no. 1, pp. 26–34, 2003. View at: Publisher Site  Google Scholar
 Y. Zhang, Y. M. Feng, and J. Deng, “Applications of low energy megavoltage Xray beams in cancer radiotherapy,” British Journal of Medicine and Medical Research, vol. 6, no. 7, pp. 661–674, 2015. View at: Publisher Site  Google Scholar
 I. Madani, B. Vanderstraeten, S. Bral et al., “Comparison of 6 MV and 18 MV photons for IMRT treatment of lung cancer,” Radiotherapy & Oncology, vol. 82, no. 1, pp. 63–69, 2007. View at: Publisher Site  Google Scholar
 X.Q. Lu, C. Burman, B. Mychalczak et al., “Feasibility study of using lowenergy intensity modulated radiation for the treatment of localized prostate carcinoma to high doses,” in Proceedings of the 22nd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 1, pp. 449–451, Chicago, Ill, USA, 2000. View at: Publisher Site  Google Scholar
 S. Yartsev, T. Kron, and J. Van Dyk, “Tomotherapy as a tool in imageguided radiation therapy (IGRT): current clinical experience and outcomes,” Biomedical Imaging and Intervention Journal, vol. 3, no. 1, article e17, 2007. View at: Publisher Site  Google Scholar
 S. L. Meeks, J. F. Harmon Jr., K. M. Langen, T. R. Willoughby, T. H. Wagner, and P. A. Kupelian, “Performance characterization of megavoltage computed tomography imaging on a helical tomotherapy unit,” Medical Physics, vol. 32, no. 8, pp. 2673–2681, 2005. View at: Publisher Site  Google Scholar
 D. Parsons and J. L. Robar, “Beam generation and planar imaging at energies below 2.40 MeV with carbon and aluminum linear accelerator targets,” Medical Physics, vol. 39, no. 7, pp. 4568–4578, 2012. View at: Publisher Site  Google Scholar
 J. L. Robar, T. Connell, W. Huang, and R. G. Kelly, “Megavoltage planar and conebeam imaging with lowZ targets: dependence of image quality improvement on beam energy and patient separation,” Medical Physics, vol. 36, no. 9, pp. 3955–3963, 2009. View at: Publisher Site  Google Scholar
 Y. Zhang, Y. Feng, Y. Zhang et al., “Is it the time for personalized imaging protocols in cancer radiation therapy?” International Journal of Radiation Oncology Biology Physics, vol. 91, no. 3, pp. 659–660, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Ying Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.