A method for converting dose-to-medium to dose-to-tissue in Monte Carlo studies of gold nanoparticle-enhanced radiotherapy

The Effect of Dose Enhancement in Tumor With Silver Nanoparticles on Surrounding Healthy Tissues: A Monte Carlo Study

Objectives: Cancer-related death rates account for approximately one-third of all deaths, and this rate is increasing remarkably every year. In this study, we examined the dose enhancement factor (DEF) in the tumor and surrounding tissues by adding different concentrations of silver nanoparticles (AgNPs) to the brain tumor using the Monte Carlo (MC) technique. Methods: This study used MCNP6.2 simulation software. A Planning Target Volume (PTV) of 1 × 1 × 1 cm3 was placed in the center of a cubic cranial model with dimensions of 5 × 5 × 5 cm3. Five different simulations were initially generated using the simple method. These simulations included pure PTV and PTV consisting of 4 different silver concentrations (5, 10, 20, and 30 mg/g). Additionally, a model was created using the nanolattice method, considering the size, position, and distribution of the AgNPs. Irradiation was performed using a source with a 6 MV linac photon spectrum. Measurements were performed using the *f8 tally, and DEF values were calculated. Results: In the simulation study using the simple method, the DEF value of PTV increased linearly with concentration, whereas the DEF values were lower than the simulation results with the nanolattice model (1.9 vs 1.4 for 30 mg/g NP concentration). Performing the simple method, we observed no remarkable dose increase in lateral OARs surrounding PTV. While a remarkable dose decrease was observed in distal OARs, a dose increase in the proximal OAR was observed, which was consistent with that of PTV. However, according to the results obtained by performing the nanolattice method, the dose increase was observed in both the proximal OAR and the distal OAR and was similar to that of PTV. Conclusion: While enhancing the dose in the tumor by adding NPs into the tumor, it is essential to consider whether it also increases the OAR dose. In addition, simulation studies on NPs showed that the dose increase varied significantly with particle size, position, and distribution. Hence, these factors should be considered carefully.

Multiscale Monte Carlo simulations of gold nanoparticle dose-enhanced radiotherapy I: Cellular dose enhancement in microscopic models

BACKGROUND

The introduction of Gold NanoParticles (GNPs) in radiotherapy treatments necessitates considerations such as GNP size, location, and quantity, as well as patient geometry and beam quality. Physics considerations span length scales across many orders of magnitude (nanometer-to-centimeter), presenting challenges that often limit the scope of dosimetric studies to either micro- or macroscopic scales.

PURPOSE

To investigate GNP dose-enhanced radiation Therapy (GNPT) through Monte Carlo (MC) simulations that bridge micro-to-macroscopic scales. The work is presented in two parts, with Part I (this work) investigating accurate and efficient MC modeling at the single cell level to calculate nucleus and cytoplasm Dose Enhancement Factors (n,cDEFs), considering a broad parameter space including GNP concentration, GNP intracellular distribution, cell size, and incident photon energy. Part II then evaluates cell dose enhancement factors across macroscopic (tumor) length scales.

METHODS

Different methods of modeling gold within cells are compared, from a contiguous volume of either pure gold or gold-tissue mixture to discrete GNPs in a hexagonal close-packed lattice. MC simulations with EGSnrc are performed to calculate n,cDEF for a cell with radiusr cell = 7.35 $r_{\rm cell}=7.35$  µm and nucleusr nuc = 5 $r_{\rm nuc} = 5$  µm considering 10 to 370 keV incident photons, gold concentrations from 4 to 24 mgAu /gtissue , and three different GNP configurations within the cell: GNPs distributed around the surface of the nucleus (perinuclear) or GNPs packed into one (or four) endosome(s). Select simulations are extended to cells with different cell (and nucleus) sizes: 5 µm (2, 3, and 4 µm), 7.35 µm (4 and 6 µm), and 10 µm (7, 8, and 9 µm).

RESULTS

n,cDEFs are sensitive to the method of modeling gold in the cell, with differences of up to 17% observed; the hexagonal lattice of GNPs is chosen (as the most realistic model) for all subsequent simulations. Across cell/nucleus radii, source energies, and gold concentrations, both nDEF and cDEF are highest for GNPs in the perinuclear configuration, compared with GNPs in one (or four) endosome(s). Across all simulations of the (rcell , rnuc ) = (7.35, 5) µm cell, nDEFs and cDEFs range from unity to 6.83 and 3.87, respectively. Including different cell sizes, nDEFs and cDEFs as high as 21.5 and 5.5, respectively, are observed. Both nDEF and cDEF are maximized at photon energies above the K- or L-edges of gold by 10 to 20 keV.

CONCLUSIONS

Considering 5000 unique simulation scenarios, this work comprehensively investigates many physics trends on DEFs at the cellular level, including demonstrating that cellular DEFs are sensitive to gold modeling approach, intracellular GNP configuration, cell/nucleus size, gold concentration, and incident source energy. These data should prove especially useful in research as well as treatment planning, allowing one to optimize or estimate DEF using not only GNP uptake, but also account for average tumor cell size, incident photon energy, and intracellular configuration of GNPs. Part II will expand the investigation, taking the Part I cell model and applying it in cm-scale phantoms.

Accelerated brachytherapy with the Xoft electronic source used in association with iodine, gold, bismuth, gadolinium, and hafnium nano-radioenhancers

PURPOSE

The current study was designed to calculate the dose enhancement factor (DEF) of iodine (I), gold (Au), bismuth (Bi), gadolinium (Gd), and hafnium (Hf) nanoparticles (NP)s by Monte Carlo (MC) modeling of an electronic brachytherapy source in resection cavities of breast tumors.

METHODS AND MATERIALS

The GEANT4 MC code was used for simulation of a phantom containing a water-filled balloon and a Xoft source (50 kVp) to irradiate the margins of a resected breast tumor. NPs with a diameter of 20 nm and concentrations from 1 to 5% w/w were simulated in a tumor margin with 5 mm thickness as well as a hypothetical breast model consisting of spherical island-like residual tumor-remnants. The DEFs for all NPs were calculated in both models.

RESULTS

In the margin-loaded model, for the concentration of 1% w/w heavy atom, DEFs of 2.5, 2.3, 2.1, 2, and 1.7 were calculated for Bi, Au, I, Hf, and Gd NPs (descending order), which increased, almost linearly with concentration for all NPs. Moreover, normal tissue dose behind the NP-loaded margin declined significantly depending on NP type and concentration. When modeling residual tumor islands, DEF values were very close to the margin-loaded values except for Bi and I, where DEFs of 2.55 and 1.7 were seen, respectively.

CONCLUSIONS

Considerable dose enhancements were obtained for the heavy atom NPs studied in the partial breast brachytherapy with a Xoft electronic source. In addition, normal tissue doses were lowered in the points beyond the NP-loaded margin. The findings revealed promising outcomes and the probability of improved tumor control for NP-aided brachytherapy with the Xoft electronic source.

Consistency checks of results from a Monte Carlo code intercomparison for emitted electron spectra and energy deposition around a single gold nanoparticle irradiated by X-rays

Organized by the European Radiation Dosimetry Group (EURADOS), a Monte Carlo code intercomparison exercise was conducted where participants simulated the emitted electron spectra and energy deposition around a single gold nanoparticle (GNP) irradiated by X-rays. In the exercise, the participants scored energy imparted in concentric spherical shells around a spherical volume filled with gold or water as well as the spectral distribution of electrons leaving the GNP. Initially, only the ratio of energy deposition with and without GNP was to be reported. During the evaluation of the exercise, however, the data for energy deposition in the presence and absence of the GNP were also requested. A GNP size of 50 nm and 100 nm diameter was considered as well as two different X-ray spectra (50 kVp and 100kVp). This introduced a redundancy that can be used to cross-validate the internal consistency of the simulation results. In this work, evaluation of the reported results is presented in terms of integral quantities that can be benchmarked against values obtained from physical properties of the radiation spectra and materials involved. The impact of different interaction cross-section datasets and their implementation in the different Monte Carlo codes is also discussed.

On the Equivalence of the Biological Effect Induced by Irradiation of Clusters of Heavy Atom Nanoparticles and Homogeneous Heavy Atom-Water Mixtures

A multiscale local effect model (LEM)-based framework was implemented to study the cell damage caused by the irradiation of clusters of gold nanoparticles (GNPs) under clinically relevant conditions. The results were compared with those obtained by a homogeneous mixture of water and gold (MixNP) irradiated under similar conditions. To that end, Monte Carlo simulations were performed for the irradiation of GNP clusters of different sizes and MixNPs with a 6 MV Linac spectrum to calculate the dose enhancement factor in water. The capabilities of our framework for the prediction of cell damage trends are examined and discussed. We found that the difference of the main parameter driving the cell damage between a cluster of GNPs and the MixNP was less than 1.6% for all cluster sizes. Our results demonstrate for the first time a simple route to intuit the radiobiological effects of clusters of nanoparticles through the consideration of an equivalent homogenous gold/water mixture. Furthermore, the negligible difference on cell damage between a cluster of GNPs and MixNP simplifies the modelling for the complex geometries of nanoparticle aggregations and saves computational resources.

Intercomparison of Monte Carlo calculated dose enhancement ratios for gold nanoparticles irradiated by X-rays: Assessing the uncertainty and correct methodology for extended beams

Results of a Monte Carlo code intercomparison exercise for simulations of the dose enhancement from a gold nanoparticle (GNP) irradiated by X-rays have been recently reported. To highlight potential differences between codes, the dose enhancement ratios (DERs) were shown for the narrow-beam geometry used in the simulations, which leads to values significantly higher than unity over distances in the order of several tens of micrometers from the GNP surface. As it has come to our attention that the figures in our paper have given rise to misinterpretation as showing 'the' DERs of GNPs under diagnostic X-ray irradiation, this article presents estimates of the DERs that would have been obtained with realistic radiation field extensions and presence of secondary particle equilibrium (SPE). These DER values are much smaller than those for a narrow-beam irradiation shown in our paper, and significant dose enhancement is only found within a few hundred nanometers around the GNP. The approach used to obtain these estimates required the development of a methodology to identify and, where possible, correct results from simulations whose implementation deviated from the initial exercise definition. Based on this methodology, literature on Monte Carlo simulated DERs has been critically assessed.

The potential role of borophene as a radiosensitizer in boron neutron capture therapy (BNCT) and particle therapy (PT)

The potential role of borophene as a radiosensitizer in PT and BNCT was investigated. Our study focused on two aspects: (1) the synthesis and characterization of borophene nanomaterials; and (2) biocompatibility and dose enhancement. To overcome the limitation of vapor-based technology, we successfully deployed the liquid-phase exfoliation (LPE) method to produce borophene targeting for biomedical applications. Bringing together spatial distribution and dose deposition, the in vitro microdosimetry study was carried out in the presence of borophene. A quantitative study of the dose enhancement ratio (DER) was performed with Monte-Carlo simulation. The synthesized borophene showed good biocompatibility with less than 10% cell death at a concentration of up to 0.2 mg ml-1. The uptake of borophene within individual cells penetrated through cell membranes but outside the nucleus. For proton PT, no significant change in the DER is found. For carbon PT, the DER increases by about 5% as the concentration of 10B reaches 1 mg g-1. For BNCT, a DER of more than 2 can be obtained for a concentration as low as 100 μg g-1. This study lays a foundation for utilizing novel borophene-based nanomaterials as radiosensitizers as well as imaging probes in cancer treatment.

Evaluation of dose enhancement with gold nanoparticles in kilovoltage radiotherapy using the new EGS geometry library in Monte Carlo simulation

<abstract><sec> <title>Purpose</title> <p>This study compared the dose enhancement predicted in kilovoltage gold nanoparticle-enhanced radiotherapy using the newly developed EGS lattice and the typical gold-water mixture method in Monte Carlo simulation. This new method considered the gold nanoparticle-added volume consisting of solid nanoparticles instead of a gold-water mixture. In addition, this particle method is more realistic in simulation.</p> </sec><sec> <title>Methods</title> <p>A heterogeneous phantom containing bone and water was irradiated by the 105 and 220 kVp x-ray beams. Gold nanoparticles were added to the tumour volume with concentration varying from 3–40 mg/mL in the phantom. The dose enhancement ratio (DER), defined as the ratio of dose at the tumour with and without adding gold nanoparticles, was calculated by the gold-water mixture and particle method using Monte Carlo simulation for comparison.</p> </sec><sec> <title>Results</title> <p>It is found that the DER was 1.44–4.71 (105 kVp) and 1.27–2.43 (220 kVp) for the gold nanoparticle concentration range of 3–40 mg/mL, when they were calculated by the gold-water mixture method. The DER was slightly larger and equal to 1.47–4.84 (105 kVp) and 1.29–2.5 (220 kVp) for the same concentration range, when the particle method was used. Moreover, the DER predicted by both methods increased with an increase of nanoparticle concentration, and a decrease of x-ray beam energy.</p> </sec><sec> <title>Conclusion</title> <p>The deviation of DER determined by the particle and gold-water mixture method was insignificant when considering the uncertainty in the calculation of DER (2%) in the nanoparticle concentration range of 3–40 mg/mL. It is therefore concluded that the gold-water mixture method could predict the dose enhancement as accurate as the newly developed particle method.</p> </sec></abstract>

Technical Note: Taking EGSnrc to new lows: Development of egs++ lattice geometry and testing with microscopic geometries

PURPOSE

This work introduces a new lattice geometry library, egs_lattice, into the EGSnrc Monte Carlo code, which can be used for both modeling very large (previously unfeasible) quantities of geometries (e.g., cells or gold nanoparticles (GNPs)) and establishing recursive boundary conditions. The reliability of egs_lattice, as well as EGSnrc in general, is cross-validated and tested at short length scales and low energies.

METHODS

New Bravais, cubic, and hexagonal lattice geometries are defined in egs_lattice and their transport algorithms are described. Simulations of cells and GNP-containing cavities are implemented to compare to independent, published Geant4-DNA and PENELOPE results. Recursive boundary conditions, implemented through a cubic lattice, are used to perform electron Fano cavity tests. The Fano test is performed on three different-sized cells containing GNPs in the region around the nucleus for three source energies.

RESULTS

Lattices are successfully implemented in EGSnrc, and are used for validation. EGSnrc calculated the dose to cell cytoplasm and nucleus when irradiated by an internal electron source with a median difference of 0.6% compared to published Geant4-DNA results. EGSnrc calculated the ratio of dose to a microscopic cavity containing GNPs over dose to a cavity containing a homogeneous mixture of gold, and results generally agree (within 1%) with published PENELOPE results. The electron Fano cavity test is passed for all energies and cells considered, with sub-0.1% discrepancies between EGSnrc-calculated and expected values. Additionally, the recursive boundary conditions used for the Fano test provided a factor of over a million increase in efficiency in some cases.

CONCLUSIONS

The egs_lattice geometry library, currently available as a pull request on the EGSnrc GitHub "develop" branch, is now freely accessible as open-source code. Lattice geometry implementations cross-validated with independent simulations in other MC codes and verified with the electron Fano cavity test demonstrate not only the reliability of egs_lattice, but also, by extension, EGSnrc's ability to simulate transport in nanometer geometries and score in microscopic cavities.

On the importance of modeling gold nanoparticles distribution in dose-enhanced radiotherapy

Purpose: To investigate the effect of precise modeling for Monte Carlo simulations of gold nanoparticles (GNPs) dose-enhanced radiotherapy, two models characterized by their distribution of GNPs in a simulated macroscopic cubic tumor were introduced. The motivation was the widely documented tendency of GNPs to localize around the cell nucleus.

Methods: The introduced models composed of 2.7×107 ellipsoid cells, each of them containing a centrally located nucleus as the target for dose evaluation. In the first model, the spheres of GNP are homogeneously distributed in the whole tumor volume, and in the latter, GNPs are localized in the cytoplasms surrounded the nuclei.

Results: The results achieved through applying Monte Carlo radiation transports using the Mont Carlo N-Particle eXtended code (MCNPX) show an underestimation of nuclear dose enhancement caused by homogeneous model compared with that of heterogeneous distribution. By investigating various quantities, it was found that subcellular location of GNPs strongly governs the sensitivity of dose enhancement to the number and concentration of GNPs targeted in the tumor. Other obvious differences are revealed by studying the dose enhancement curves in depth of the tumor. While the heterogeneous model predicts an approximately constant dose enhancement in depth for primary photon energies of 50 keV and more, the homogeneous model estimates an energy-dependent increase of about 11 to 30%.

Conclusion: It can be concluded that defining a model in accordance with the experimental observations can effectively account for accurate prediction of macroscopic dose enhancement in the target of interest.

Monte Carlo dosimetry modeling of focused kV x-ray radiotherapy of eye diseases with potential nanoparticle dose enhancement

PURPOSE

Eye plaque brachytherapy is the most common approach for intraocular cancer treatment. It is, however, invasive and subject to large setup uncertainty due to the surgical operation. We propose a novel-focused kV x-ray technique with potential nanoparticle (NP) enhancement and evaluate its application in treating choroidal melanoma and iris melanoma by Monte Carlo (MC) dosimetry modeling.

METHODS

A polycapillary x-ray lens was used to focus 45 kVp x rays to achieve pinpoint accuracy of dose delivery to small tumors near critical structures. In addition to allowing for beam focusing, the use of kV x rays takes advantage of the strong photoelectric absorption of metallic NPs in that energy regime and hence strong radiosensitization. We constructed an MC simulation program that takes into account the x-ray optic modeling and used GEANT4 for dosimetric calculation. Extensive phantom measurements using a prototype-focused x-ray system were carried out. The MC simulation of simple geometry phantom irradiation was first compared to measurements to verify the x-ray optic lens modeling in conjunction with the Geant4 dosimetric calculation. To simulate tumor treatment, a geometric eye model and two tumor models were constructed. Dose distributions of the simulated treatments were then calculated. NP radiosensitization was also simulated for two concentrations of 2 nm gold NP (AuNP) uniformly distributed in the tumor.

RESULTS

The MC-simulated full width at half maximum (FWHM) and central-axis depth dose of the focused kV x-ray beam match those measured on EBT3 films within ~10% around the depth of focus of the beam. Dose distributions of the simulated ocular tumor treatments show that focused x-ray beams can concentrate the high-dose region in or close to the tumor plus margin. For the simulated posterior choroidal tumor treatment, with sufficient tumor coverage, the doses to the optic disc and fovea are substantially reduced with focused x-ray therapy compared to eye plaque treatment (3.8 vs 39.8 Gy and 11.1 vs 53.8 Gy, respectively). The eye plaque treatment was calculated using an Eye Physics plaque with I-125 seeds under TG43 assumption. For the energy spectrum used in this study, the average simulated dose enhancement ratios (DERs) are roughly 2.1 and 1.1 for 1.0% and 0.1% AuNP mass concentration in the tumor, respectively.

CONCLUSION

Compared to eye plaque brachytherapy, the proposed focused kV x-ray technique is noninvasive and shows great advantage in sparing healthy critical organs without sacrificing the tumor control. The NP radiation dose enhancement is considerable at our proposed kV range even with a low NP concentration in the tumor, providing better critical structure protection and more flexibility for treatment planning.

Gold Nanoparticles in Radiotherapy and Recent Progress in Nanobrachytherapy

Over the last few decades, gold nanoparticles (GNPs) have emerged as "radiosensitizers" in oncology. Radiosensitizers are additives that can enhance the effects of radiation on biological tissues treated with radiotherapy. The interaction of photons with GNPs leads to the emission of low-energy and short-range secondary electrons, which in turn increase the dose deposited in tissues. In this context, GNPs are the subject of intensive theoretical and experimental studies aiming at optimizing the parameters leading to greater dose enhancement and highest therapeutic effect. This review describes the main mechanisms occurring between photons and GNPs that lead to dose enhancement. The outcome of theoretical simulations of the interactions between GNPs and photons is presented. Finally, the findings of the most recent in vivo studies about interactions between GNPs and photon sources (e.g., external beams, brachytherapy sources, and molecules labeled with radioisotopes) are described. The advantages and challenges inherent to each of these approaches are discussed. Future directions, providing new guidelines for the successful translation of GNPs into clinical applications, are also highlighted.

Energy optimization in gold nanoparticle enhanced radiation therapy

Gold nanoparticles (GNPs) have been demonstrated as radiation dose enhancing agents. Kilovoltage external photon beams have been shown to yield the largest enhancement due to the high interaction probability with gold. While orthovoltage irradiations are feasible and promising, they suffer from a reduced tissue penetrating power. This study quantifies the effect of varying photon beam energies on various beam arrangements, body, tumor, and cellular GNP uptake geometries. Cell survival was modeled based on our previously developed GNP-local effect model with radial doses calculated using the TOPAS-nBio Monte Carlo code. Cell survival curves calculated for tumor sites with GNPs were used to calculate the relative biological effectiveness (RBE)-weighted dose. In order to evaluate the plan quality, the ratio of the mean dose between the tumor and normal tissue for 50-250 kVp beams with GNPs was compared to the standard of care using 6 MV photon beams without GNPs for breast and brain tumors. For breast using a single photon beam, kV  +  GNP was found to yield up to 2.73 times higher mean RBE-weighted dose to the tumor than two tangential megavoltage beams while delivering the same dose to healthy tissue. For irradiation of brain tumors using multiple photon beams, the GNP dose enhancement was found to be effective for energies above 50 keV. A small tumor at shallow depths was found to be the most effective treatment conditions for GNP enhanced radiation therapy. GNP uptake distributions in the cell (with or without nuclear uptake) and the beam arrangement were found to be important factors in determining the optimal photon beam energy.

Heterogeneous multiscale Monte Carlo simulations for gold nanoparticle radiosensitization

PURPOSE

To introduce the heterogeneous multiscale (HetMS) model for Monte Carlo simulations of gold nanoparticle dose-enhanced radiation therapy (GNPT), a model characterized by its varying levels of detail on different length scales within a single phantom; to apply the HetMS model in two different scenarios relevant for GNPT and to compare computed results with others published.

METHODS

The HetMS model is implemented using an extended version of the EGSnrc user-code egs_chamber; the extended code is tested and verified via comparisons with recently published data from independent GNP simulations. Two distinct scenarios for the HetMS model are then considered: (a) monoenergetic photon beams (20 keV to 1 MeV) incident on a cylinder (1 cm radius, 3 cm length); (b) isotropic point source (brachytherapy source spectra) at the center of a 2.5 cm radius sphere with gold nanoparticles (GNPs) diffusing outwards from the center. Dose enhancement factors (DEFs) are compared for different source energies, depths in phantom, gold concentrations, GNP sizes, and modeling assumptions, as well as with independently published values. Simulation efficiencies are investigated.

RESULTS

The HetMS MC simulations account for the competing effects of photon fluence perturbation (due to gold in the scatter media) coupled with enhanced local energy deposition (due to modeling discrete GNPs within subvolumes). DEFs are most sensitive to these effects for the lower source energies, varying with distance from the source; DEFs below unity (i.e., dose decreases, not enhancements) can occur at energies relevant for brachytherapy. For example, in the cylinder scenario, the 20 keV photon source has a DEF of 3.1 near the phantom's surface, decreasing to less than unity by 0.7 cm depth (for 20 mg/g). Compared to discrete modeling of GNPs throughout the gold-containing (treatment) volume, efficiencies are enhanced by up to a factor of 122 with the HetMS approach. For the spherical phantom, DEFs vary with time for diffusion, radionuclide, and radius; DEFs differ considerably from those computed using a widely applied analytic approach.

CONCLUSIONS

By combining geometric models of varying complexity on different length scales within a single simulation, the HetMS model can effectively account for both macroscopic and microscopic effects which must both be considered for accurate computation of energy deposition and DEFs for GNPT. Efficiency gains with the HetMS approach enable diverse calculations which would otherwise be prohibitively long. The HetMS model may be extended to diverse scenarios relevant for GNPT, providing further avenues for research and development.

Optimization of photon beam energies in gold nanoparticle enhanced arc radiation therapy using Monte Carlo methods

As a recent area of development in radiation therapy, gold nanoparticle (GNP) enhanced radiation therapy has shown potential to increase tumour dose while maintaining acceptable levels of healthy tissue toxicity. In this study, the effect of varying photon beam energy in GNP enhanced arc radiation therapy (GEART) is quantified through the introduction of a dose scoring metric, and GEART is compared to a conventional radiotherapy treatment. The PENELOPE Monte Carlo code was used to model several simple phantoms consisting of a spherical tumour containing GNPs (concentration: 15 mg Au g^-1 tumour, 0.8 mg Au g^-1 normal tissue) in a cylinder of tissue. Several monoenergetic photon beams, with energies ranging from 20 keV to 6 MeV, as well as 100, 200, and 300 kVp spectral beams, were used to irradiate the tumour in a 360° arc treatment. A dose metric was then used to compare tumour and tissue doses from GEART treatments to a similar treatment from a 6 MV spectrum. This was also performed on a simulated brain tumour using patient computed tomography data. GEART treatments showed potential over the 6 MV treatment for many of the simulated geometries, delivering up to 88% higher mean dose to the tumour for a constant tissue dose, with the effect greatest near a source energy of 50 keV. This effect is also seen with the inclusion of bone in a brain treatment, with a 14% increase in mean tumour dose over 6 MV, while still maintaining acceptable levels of dose to the bone and brain.