Efficient photon beam dose calculations using DOSXYZnrc with BEAMnrc

An automation of Monte Carlo workflow for dosimetry study of an Elekta LINAC delivery system in radiotherapy

Purpose

This study aims to automate the Monte Carlo (MC) workflow utilized for radiotherapy dosimetry, focusing on an Elekta LINAC delivery system. It addresses the challenge of integrating MC simulations into routine clinical practice, making this accurate yet complex method more accessible and efficient for radiotherapy dosimetry.

Methods and Materials

We developed a user-friendly software featuring a graphical user interface (GUI) that integrates EGSnrc for MC simulations. The software streamlines the process from retrieving Digital Imaging and Communications in Medicine (DICOM) data to executing dose calculations and comparing dose distributions. To validate our proposed tool, we compared its computed doses for IMRT and VMAT plans from the Pinnacle TPS for an Elekta Versa HD linear accelerator against MC simulation results. This comparison utilized our in-house software and GUI as the tool, covering various treatment sites and prescriptions.

Results

The automated MC workflow demonstrated high accuracy in dose calculations and streamlined integration with clinical workflows. The comparison between the MC-simulated and TPS-calculated doses revealed excellent agreement, highlighting the reliability of MC for independent dose verification in complex treatment scenarios.

Conclusions

The automated MC workflow developed represents a substantial improvement in the practicality and efficiency of MC simulations in radiotherapy. This advancement not only simplifies the dosimetry process but also ensures high accuracy, establishing it as a valuable tool for routine patient-specific quality assurance and the development of specialized treatment procedures.

How efficient are metal-polymer and dual-metals-polymer non-lead radiation shields?

INTRODUCTION

Lead shields are often used to attenuate ionising radiations. However, to make lighter, recyclable and more efficient shields compared to lead, combinations of new metallic compounds together with polymer, for example, flexible polyvinyl chloride (PVC) have been developed recently. In this study, the capabilities of non-lead radiation shields made of one or two metallic compounds and polymer were evaluated.

METHODS

Monte Carlo (MC)-based BEAMnrc code was used to build a functional model based on a Philips X-ray machine in the range of radiographic energies. The MC model was then verified by IPEM Report 78 as a standardised global reference. The MC model was then used to evaluate the efficiency of non-lead-based garments made of metallic compound and polymer (MCP) including BaSO4 -PVC, Bi2 O3 -PVC, Sn-PVC and W-PVC, as well as dual-metallic compounds and polymer (DMCP) including Bi2 O3 -BaSO4 -PVC, Bi2 O3 -Sn-PVC, W-Sn-PVC and W-BaSO4 -PVC. The absorbed doses were determined at the surface of a water phantom and compared directly with the doses obtained for 0.5 mm pure lead (Pb).

RESULTS

Bi2 O3 -BaSO4 -PVC and W-BaSO4 -PVC were found to be efficient shields for most of the energies. In addition to the above radiation shields, Bi2 O3 -Sn-PVC was also found to be effective for the spectrum of 60 keV. Bi2 O3 -BaSO4 -PVC as a non-lead dual metals-PVC shield was shown to be more efficient than pure lead in diagnostic X-ray range.

CONCLUSION

Combination of two metals-PVC, a low atomic number (Z) metal together with a high atomic number metal, and also single-metal-PVC shields were shown to be efficient enough to apply as radiation protection shields instead of lead-based garments.

Viability of the virtual cone technique using a fixed small multi-leaf collimator field for stereotactic radiosurgery of trigeminal neuralgia

Dosimetric uncertainties in very small (≤1.5 × 1.5 cm2 ) photon fields are remarkably higher, which undermines the validity of the virtual cone (VC) technique with a diminutive and variable MLC fields. We evaluate the accuracy and reproducibility of the VC method with a very small, fixed MLC field setting, called a fixed virtual cone (fVC), for small target radiosurgery such as trigeminal neuralgia (TGN). The fVC is characterized by 0.5 cm x 0.5 cm high-definition (HD) MLC field of 10MV FFF beam defined at 100 cm SAD, while backup jaws are positioned at 1.5 cm x 1.5 cm. A spherical dose distribution equivalent to 5 mm (diameter) physical cone was generated using 10-14 non-coplanar, partial arcs. Dosimetric accuracy was validated using SRS diode (PTW 60018), SRS MapCHECK (SNC) measurements. As a quality assurance measure, 10 treatment plans (SRS) for TGN, consisting of various arc ranges at different collimator angles were analyzed using 6 MV FFF and 10 MV FFF beams, including a field-by-field study (n = 130 fields). Dose outputs were compared between the Eclipse TPS and measurements (SRS MapCHECK). Moreover, dosimetric changes in the field defining fVC, prompted by a minute (± 0.5-1.0 mm) leaf shift, was examined among TPS, diode measurements, and Monte Carlo (MC) simulations. The beam model for fVC was validated (≤3% difference) using SRS MapCHECK based absolute dose measurements. The equivalent diameters of the 50% isodose distribution were found comparable to that of a 5 mm cone. Additionally, the comparison of field output factors, dose per MU between the TPS and SRS diode measurements using the fVC field, including ± 1 mm leaf shift, yielded average discrepancies within 5.5% and 3.5% for 6 MV FFF and 10 MV FFF beams, respectively. Overall, the fVC method is a credible alternative to the physical cone (5 mm) that can be applied in routine radiosurgical treatment of TGN.

Monte-Carlo techniques for radiotherapy applications I: introduction and overview of the different Monte-Carlo codes

Introduction:

The dose calculation plays a crucial role in many aspects of contemporary clinical radiotherapy treatment planning process. It therefore goes without saying that the accuracy of the dose calculation is of very high importance. The gold standard for absorbed dose calculation is the Monte-Carlo algorithm.

Methods:

This first of two papers gives an overview of the main openly available and supported codes that have been widely used for radiotherapy simulations.

Results:

The paper aims to provide an overview of Monte-Carlo in the field of radiotherapy and point the reader in the right direction of work that could help them get started or develop their existing understanding and use of Monte-Carlo algorithms in their practice.

Conclusions:

It also serves as a useful companion to a curated collection of papers on Monte-Carlo that have been published in this journal.

Assessing the Effect of Directional Bremsstrahlung Splitting on the Output Spectra and Parameters Using BEAMnrc Monte Carlo Simulation Package

Introduction:

EGSnrc software package is one of the computational packages for Monte Carlo simulation in radiation therapy and has several subset codes. Directional bremsstrahlung splitting (DBS) is a technique that applies braking radiations in interactions in this software. This study aimed to evaluate the effect of this technique on the simulation time, uncertainty, particle number of phase-space data, and photon beam spectrum resulting from a medical linear accelerator (LINAC).

Materials and methods:

The gantry of the accelerator, including the materials and geometries of different parts, was simulated using the BEAMnrc code (a subset code in the EGSnrc package). The phase-space data were recorded in different parts of the LINAC. The DBS values (1, 10, 100, and 1000) were changed, and their effects were evaluated on the simulation parameters and output spectra.

Results:

Increasing the DBS value from 1 to 1000 resulted in an increase in the simulation time from 1.778 to 11.310 hours, and increasing the number of particles in the phase-space plane (5 590 732-180 328 382). When the DBS had been picked up from 1 to 100, the simulation uncertainty decreased by about 1.29%. In addition, the DBS increment value from 100 to 1000 leads to an increase in uncertainty and simulation time of about 0.71% and 315%, respectively.

Conclusion:

Although using the DBS technique reduces the simulation time or uncertainty, increasing the DBS from a specific value, equal to 100 in our study, increases simulation uncertainties and times. Therefore, we propose considering a specific DBS value as we obtained for the Monte Carlo simulation of photon beams produced by linear accelerators.

Variations in size-specific effective dose with patient stature and beam width for kV cone beam CT imaging in radiotherapy

The facilities now available on linear accelerators for external beam radiotherapy enable radiation fields to be conformed to the shapes of tumours with a high level of precision. However, in order for the treatment delivered to take advantage of this, the patient must be positioned on the couch with the same degree of accuracy. Kilovoltage cone beam computed tomography systems are now incorporated into radiotherapy linear accelerators to allow imaging to be performed at the time of treatment, and image-guided radiation therapy is now standard in most radiotherapy departments throughout the world. However, because doses from imaging are much lower than therapy doses, less effort has been put into optimising radiological protection of imaging protocols. Standard imaging protocols supplied by the equipment vendor are often used with little adaptation to the stature of individual patients, and exposure factors and field sizes are frequently larger than necessary. In this study, the impact of using standard protocols for imaging anatomical phantoms of varying size from a library of 193 adult phantoms has been evaluated. Monte Carlo simulations were used to calculate doses for organs and tissues for each phantom, and results combined in terms of size-specific effective dose (SED). Values of SED from pelvic scans ranged from 11 mSv to 22 mSv for male phantoms and 8 mSv to 18 mSv for female phantoms, and for chest scans from 3.8 mSv to 7.6 mSv for male phantoms and 4.6 mSv to 9.5 mSv for female phantoms. Analysis of the results showed that if the same exposure parameters and field sizes are used, a person who is 5 cm shorter will receive a size SED that is 3%-10% greater, while a person who is 10 kg lighter will receive a dose that is 10%-14% greater compared with the average size.

Development of a GPU-accelerated Monte Carlo dose calculation module for nuclear medicine, ARCHER-NM: demonstration for a PET/CT imaging procedure

Objective. This paper describes the development and validation of a GPU-accelerated Monte Carlo (MC) dose computing module dedicated to organ dose calculations of individual patients undergoing nuclear medicine (NM) internal radiation exposures involving PET/CT examination. Approach. This new module extends the more-than-10-years-long ARCHER project that developed a GPU-accelerated MC dose engine by adding dedicated NM source-definition features. To validate the code, we compared dose distributions from the point ion source, including 18F, 11C, 15O, and 68Ga, calculated for a water phantom against a well-tested MC code, GATE. To demonstrate the clinical utility and advantage of ARCHER-NM, one set of 18F-FDG PET/CT data for an adult male NM patient is calculated using the new code. Radiosensitive organs in the CT dataset are segmented using a CNN-based tool called DeepViewer. The PET image intensity maps are converted to radioactivity distributions to allow for MC radiation transport dose calculations at the voxel level. The dose rate maps and corresponding statistical uncertainties were calculated at the acquisition time of PET image. Main results. The water-phantom results show excellent agreement, suggesting that the radiation physics module in the new NM code is adequate. The dose rate results of the 18F-FDG PET imaging patient show that ARCHER-NM’s results agree very well with those of the GATE within −2.45% to 2.58% (for a total of 28 organs considered in this study). Most impressively, ARCHER-NM obtains such results in 22 s while it takes GATE about 180 min for the same number of 5 × 108 simulated decay events. Significance. This is the first study presenting GPU-accelerated patient-specific MC internal radiation dose rate calculations for clinically realistic 18F-FDG PET/CT imaging case involving autosegmentation of whole-body PET/CT images. This study suggests that the proposed computing tools—ARCHER-NM— are accurate and fast enough for routine internal dosimetry in NM clinics.

Engaging medical physics students in active and authentic learning through the use of monte-carlo simulation and inverse treatment planning

PURPOSE

Engagement and participation of students with the learning process has been recognised as a growing problem across the higher education sector. The aim of this study was to investigate the value and impact of introducing Problem-based Learning (PBL) activities into a radiotherapy physics unit of a postgraduate medical physics course.

METHODS

Computer-based problem solving activities on 1) monte-carlo modelling of a linear accelerator and 2) inverse radiotherapy treatment planning were designed and implemented into a one-semester unit on radiotherapy physics. The value and impact of the activities on the student learning were evaluated through student surveys, a focus group, and peer observation of the sessions by members of the learning design team. Student attendance and grade profile data are also reported.

RESULTS

Overall the results indicated that students had a positive experience with the new problem solving activities that were implemented. Survey responses from a number of students indicated a desire for increased theoretical and technical support prior to and during activities. Another underlying theme that emerged from survey and focus group response was the perceived lack of reward in terms of marks for their efforts working on the learning activities. This may have influenced students' choices around attendance and participation. No significant changes were noted in the overall grades achieved in the unit.

CONCLUSIONS

Students appreciated the more hands-on approach to learning in the form of more authentic activities that they could directly relate to clinical radiotherapy. Further work is required to update and integrate assessment into the new learning delivery model more directly.

The relation between XR-QA2 and RT-QA2 GafchromicTMfilm optical density and absorbed dose in water produced by radionuclides

Purpose. In this study, Monte Carlo (MC) simulations were done to relate the dose-response of the film to that in water. The effect of backscattering materials (PMMA, lead, polystyrene, and air) was investigated on its influence on film density for radionuclides including Am-241, Tc-99m, I-131, Cs-137.Methods. A BEAMnrc MC simulation was designed to score a phase-space file (PSF) below the container of the radionuclide under consideration to use as an input file for the subsequent DOSXYZnrc MC simulation. The geometry of the container holding the radionuclide was built using the component modules available in BEAMnrc. BEAMDP was used to investigate the container effect on the radionuclide spectrum as well as the fluence. The DOSXYZnrc simulation produced the absorbed dose in XR-QA2 and RT-QA2 Gafchromic^TMfilms. The DOSXYZnrc simulations were repeated for the Gafchromic^TMfilm now replaced with water to get the absorbed dose in water. From these results, conversion factors for the dose in water to the film dose for the different radionuclides, Am-241, Tc-99m, I-131, and Cs-137 were obtained. The actual film dose was calculated using the specific gamma exposure constant (Γ) at a distance of 50 cm for a point source approximation. From the BEAMnrc simulations, the particle fluence was extracted from PSFs to correct for the fluence at 0.1 cm below the sources from the fluence 50 cm away since the inverse square law will not apply to finite-size sources. The absorbed dose profiles in the film were compared to the absorbed dose profiles from the MC simulations.Results. A fitting function based on the neutron depletion model fits the optical density versus absorbed film dose data well and can be used as a calibration tool to obtain the film dose from its optical density. Lead as a backscatter material results in a higher optical density change but a lower absorbed dose. The XR-QA2 Gafchromic^TMfilm is more sensitive than the RT-QA2 Gafchromic^TMfilm, showing a more responsive optical density (OD) change in the energy range of radionuclides used in this study. Conversion factors were determined to convert the dose in water to the dose in Gafchromic^TMfilm. The Am-241 and I-131 simulated absorbed dose in the film to dose in water does not fluctuate as much as the simulated absorbed dose in film and water when using Tc-99m and Cs-137. Validation was shown for the comparison of the film and MC simulation absorbed dose profiles.Conclusions. MC BEAMnrc simulations are useful to simulate radionuclides and their containers. BEAMDP extracted energy spectra showed that the radionuclide containers produced a Compton effect on the energy spectra and added filtration on the lower spectral photon components. Extracted fluence ratios from PSFs were used to calculate the absorbed dose value at 0.1 cm distance from the source. By using the fit function, the dose in the film can be determined for known optical density values. The effect of the backscatter materials showed that the XR-QA2 Gafchromic^TMfilm results in higher optical density values than the RT-QA2 Gafchromic^TMfilm. The absorbed dose in both the films is comparable but not for a radionuclide such as Am-241 with an activity of 74MBq. The lead backscatter material showed to be the most prominent in optical density enhancement, and the air equivalent material was the least prominent. The XR-QA2 Gafchromic^TMfilm is the most sensitive and will be the best option if working with low energies. The absorbed dose in the XR-QA2 Gafchromic^TMfilm also showed a good comparison to the absorbed dose in water for the Am-241 radionuclide with an activity of 74MBq. The absorbed dose in the films compares well to the MC simulated doses.

Simulation study of a coincidence detection system for non-invasive determination of arterial blood time-activity curve measurements

Background

Arterial sampling in PET studies for the purposes of kinetic modeling remains an invasive, time-intensive, and expensive procedure. Alternatives to derive the blood time-activity curve (BTAC) non-invasively are either reliant on large vessels in the field of view or are laborious to implement and analyze as well as being prone to many processing errors. An alternative method is proposed in this work by the simulation of a non-invasive coincidence detection unit.

Results

We utilized GATE simulations of a human forearm phantom with a blood flow model, as well as a model for dynamic radioactive bolus activity concentration based on clinical measurements. A fixed configuration of 14 and, also separately, 8 detectors were employed around the phantom, and simulations were performed to investigate signal detection parameters. Bismuth germanate (BGO) crystals proved to show the highest count rate capability and sensitivity to a simulated BTAC with a maximum coincidence rate of 575 cps. Repeatable location of the blood vessels in the forearm allowed a half-ring design with only 8 detectors. Using this configuration, maximum coincident rates of 250 cps and 42 cps were achieved with simulation of activity concentration determined from ¹⁵O and ¹⁸F arterial blood sampling. NECR simulated in a water phantom at 3 different vertical positions inside the 8-detector system (Y = − 1 cm, Y = − 2 cm, and Y = −3 cm) was 8360 cps, 13,041 cps, and 20,476 cps at an activity of 3.5 MBq. Addition of extra axial detection rings to the half-ring configuration provided increases in system sensitivity by a factor of approximately 10.

Conclusions

Initial simulations demonstrated that the configuration of a single half-ring 8 detector of monolithic BGO crystals could describe the simulated BTAC in a clinically relevant forearm phantom with good signal properties, and an increased number of axial detection rings can provide increased sensitivity of the system. The system would find use in the derivation of the BTAC for use in the application of kinetic models without physical arterial sampling or reliance on image-based techniques.

Calculated beam quality correction factors for ionization chambers in MV photon beams

The beam quality correction factor, [Formula: see text], which corrects for the difference in the ionization chamber response between the reference and clinical beam quality, is an integral part of radiation therapy dosimetry. The uncertainty of [Formula: see text] is one of the most significant sources of uncertainty in the dose determination. To improve the accuracy of available [Formula: see text] data, four partners calculated [Formula: see text] factors for 10 ionization chamber models in linear accelerator beams with accelerator voltages ranging from 6 MV to 25 MV, including flattening-filter-free (FFF) beams. The software used in the calculations were EGSnrc and PENELOPE, and the ICRU report 90 cross section data for water and graphite were included in the simulations. Volume averaging correction factors were calculated to correct for the dose averaging in the chamber cavities. A comparison calculation between partners showed a good agreement, as did comparison with literature. The [Formula: see text] values from TRS-398 were higher than our values for each chamber where data was available. The [Formula: see text] values for the FFF beams did not follow the same [Formula: see text], [Formula: see text] relation as beams with flattening filter (values for 10 MV FFF beams were below fits made to other data on average by 0.3%), although our FFF sources were only for Varian linacs.

New capabilities of the Monte Carlo dose engine ARCHER-RT: Clinical validation of the Varian TrueBeam machine for VMAT external beam radiotherapy

PURPOSE

The Monte Carlo radiation transport method is considered the most accurate approach for absorbed dose calculations in external beam radiation therapy. In this study, an efficient and accurate source model of the Varian TrueBeam 6X STx Linac is developed and integrated with a fast Monte Carlo photon-electron transport absorbed dose engine, ARCHER-RT, which is capable of being executed on CPUs, NVIDIA GPUs, and AMD GPUs. This capability of fast yet accurate radiation dose calculation is essential for clinical utility of this new technology. This paper describes the software and algorithmic developments made to the ARCHER-RT absorbed dose engine.

METHODS

AMD's Heterogeneous-Compute Interface for Portability (HIP) was implemented in ARCHER-RT to allow for device independent execution on NVIDIA and AMD GPUs. Architecture-specific atomic-add algorithms have been identified and both more accurate single-precision and double-precision computational absorbed dose calculation methods have been added to ARCHER-RT and validated through a test case to evaluate the accuracy and performance of the algorithms. The validity of the source model and the radiation transport physics were benchmarked against Monte Carlo simulations performed with EGSnrc. Secondary dose-check physics plans, and a clinical prostate treatment plan were calculated to demonstrate the applicability of the platform for clinical use. Absorbed dose difference maps and gamma analyses were conducted to establish the accuracy and consistency between the two Monte Carlo models. Timing studies were conducted on a CPU, an NVIDIA GPU, and an AMD GPU to evaluate the computational speed of ARCHER-RT.

RESULTS

Percent depth doses were computed for different field sizes ranging from 1.5 cm²  × 1.5 cm² to 22 cm²  × 40cm² and the two codes agreed for all points outside high gradient regions within 3%. Axial profiles computed for a 10 cm²  × 10 cm² field for multiple depths agreed for all points outside high gradient regions within 2%. The test case investigating the impact of native single-precision compared to double-precision showed differences in voxels as large as 71.47% and the implementation of KAS single-precision reduced the difference to less than 0.01%. The 3%/3mm gamma pass rates for an MPPG5a multileaf collimator (MLC) test case and a clinical VMAT prostate plan were 94.2% and 98.4% respectively. Timing studies demonstrated the calculation of a VMAT plan was completed in 50.3, 187.9, and 216.8 s on an NVIDIA GPU, AMD GPU, and Intel CPU, respectively.

CONCLUSION

ARCHER-RT is capable of patient-specific VMAT external beam photon absorbed dose calculations and its potential has been demonstrated by benchmarking against a well validated EGSnrc model of a Varian TrueBeam. Additionally, the implementation of AMD's HIP has shown the flexibility of the ARCHER-RT platform for device independent calculations. This work demonstrates the significant addition of functionality added to ARCHER-RT framework which has marked utility for both research and clinical applications and demonstrates further that Monte Carlo-based absorbed dose engines like ARCHER-RT have the potential for widespread clinical implementation.

Development and validation of a 1.5 T MR-Linac full accelerator head and cryostat model for Monte Carlo dose simulations

PURPOSE

To develop, implement, and validate a full 1.5 T/7 MV magnetic resonance (MR)-Linac accelerator head and cryostat model in EGSnrc for high precision dose calculations accounting for magnetic field effects that are independent from the vendor treatment planning system.

METHODS

Primary electron beam parameters for the implemented model were adapted to be in accordance with measured dose profiles of the Elekta Unity (Elekta AB, Stockholm, Sweden). Parameters to be investigated were the mean electron energy as well as the Gaussian radial intensity and energy distributions. Energy tuning was done comparing depth dose profiles simulated with monoenergetic beams of varying energies to measurements. The optimum radial intensity distribution was found by varying the radial full width at half maximum (FWHM) and comparing simulated and measured lateral profiles. The influence of the energy distribution was investigated by comparing simulated lateral and depth dose profiles with varying energy spreads to measured data. Comparison of simulations and measurements was performed by calculating average and maximum local dose deviations. The model was validated recalculating a clinical intensity-modulated radiation therapy plan for the MR-Linac and comparing the resulting dose distribution with simulations from the commercial treatment planning system Monaco using the gamma criterion.

RESULTS

Comparison of simulated and measured data showed that the optimum initial electron beam for MR-Linac simulations was monoenergetic with an electron energy of (7.4 ± 0.2) MeV. The optimum Gaussian radial intensity distribution has a FWHM of (2.2 ± 0.3) mm. The average relative deviations were smaller than 1% for all simulated profiles with optimum electron parameters, whereas the largest maximum deviation of 2.07% was found for the 22 × 22 cm 2 cross-plane profile. Profiles were insensitive to energy spread variations. The IMRT plan recalculated with the final MR-Linac model with optimized initial electron beam parameters showed a gamma pass rate of 99.83 % using a gamma criterion of 3%/3 mm.

CONCLUSIONS

The EGSnrc MR-Linac model developed in this study showed good accordance with measurements and was successfully used to recalculate a first full clinical IMRT treatment plan. Thus, it shows the general possibility for future secondary dose calculations of full IMRT plans with EGSnrc, which needs further detailed investigations before clinical use.

Optimization of Phase Space files from clinical linear accelerators

This work proposes a methodology to produce an optimized phase-space (PhSp) for the Elekta Synergy linac by tuning the energy and direction of particles inside the 6-MV Elekta Precise PhSp, provided by the International Atomic Energy Agency (IAEA), for Monte Carlo (MC) simulations. First, the energies of the particles emerging from the original PhSp were increased by different factors, producing new PhSps. Percentage depth dose (PDD) profiles were simulated and compared to measured data from a Synergy linac for 6-MV photon beam. This process was repeated until a minimum difference was reached. Particles' directions were then manipulated following identified correlations to lateral profiles, resulting in two distinct perturbation factors based on inline and crossline profiles. Both factors were merged into one single optimal factor. For energy optimization, an increase of 0.32 MeV applied to all particles inside the original PhSp, but to 0.511 MeV annihilation photons, provided the best results. The direction optimization factor was the combination of the individual factors for inline (0.605%) and crossline (0.051%). The agreement between measured and simulated profiles, when using the optimized PhSp, improved considerably in comparison to simulations performed with the original IAEA PhSp. For all fields and depths analyzed, the discrepancies for PDD, inline and crossline profiles dropped from 11.2%, 15.7% and 27.5% to under 1.4%, 4.7% and 13.2%, respectively. The optimized PhSp should not replace the full linac modelling, however it offers an alternative for MC dose calculations when neither geometric details nor validated IAEA PhSp are available to the user.

Monte Carlo Study of Unflattened Photon Beams Shaped by Multileaf Collimator

INTRODUCTION

This study investigates basic dosimetric properties of unflattened 6 MV photon beam shaped by multileaf collimator and compares them with those of flattened beams.

MATERIALS AND METHODS

Monte Carlo simulation model using BEAM code was developed for a 6MV photon beam based on Varian Clinic 600 unique performance linac operated with and without a flattening filter in beam line. Dosimetric features including lateral profiles, central axis depth dose, photon and electron spectra were calculated for flattened and unflattened cases, separately.

RESULTS

An increase in absolute depth dose with a factor of more than 2.4 was observed for unflattened beam which was dependent on depth. PDDs values were found to be lower for unflattened beam for all field sizes. Significant decrease in calculated mlc leakage was observed when the flattening filter was removed from the beam line. The total scatter factor, S_CP was found to show less variation with field sizes for unflattened beam indicating a decrease in head scatter. The beam profiles for unflattened case are found to have lower relative dose value in comparison with flattened beam near the field edge, and it falls off faster with distance.

CONCLUSION

Our study showed that increase in the dose rate and lower peripheral dose could be considered as realistic advantages for unflattened 6MV photon beams.

Investigating the effect of dental implant materials with different densities on radiotherapy dose distribution using Monte-Carlo simulation and pencil beam convolution algorithm

OBJECTIVES

The aim of this study was to investigate the effect of dental implant materials with different physical densities on dose distribution for head and neck cancer radiotherapy planning.

METHODS

Titanium (Ti), Titanium alloy (Ti-6Al-4V), Zirconia (Y-TZP), Zirconium oxide (ZrO₂), Alumina (Al₂O₃) and polyetheretherketone (PEEK) dental implant materials were used for determination of implant material effect on dose distribution. Dental implant effect was investigated by using pencil beam convolution (PBC) algorithm of Eclipse treatment planning systems (TPS) and Monte Carlo (MC) simulation technique. 6 MV photon beam of the Varian 2300 C/D linear accelerator was simulated by EGSnrc-based BEAMnrc MC code system.

RESULTS

Reasonable consistency was determined for percentage depth dose (PDD) curves between MC simulation and water phantom measurements at 6.4 MeV initial electron energy. The consistency between modelled linear accelerator PDD curve calculations and water-phantom PDD measurements were compatible within 1 % range. The dose increase in front of the dental implant calculated by MC simulation is in the range of 0.4-20.2%. We found by MC and PBC calculations that the differences in dose increase in front of the dental implant materials is in the range of 0.1-17.2% and is dependent on the physical density of the dental implant.

CONCLUSIONS

Dose increase for Zirconia was noted to be maximum while PEEK implant dose increase was minimum among the whole dental implant materials studied. This study revealed that the Eclipse TPS PBC algorithm could not accurately estimate the backscatter radiation from dental implant materials.

Commissioning of a dedicated commercial Co-60 total body irradiation unit

We describe the commissioning of the first dedicated commercial total body irradiation (TBI) unit in clinical operation. The Best Theratronics GammaBeam 500 is a Co-60 teletherapy unit with extended field size and imaging capabilities. Radiation safety, mechanical and imaging systems, and radiation output are characterized. Beam data collection, calibration, and external dosimetric validation are described. All radiation safety and mechanical tests satisfied relevant requirements and measured dose distributions meet recommendations of American Association of Physicists in Medicine (AAPM) Report #17. At a typical treatment distance, the dose rate in free space per unit source activity using the thick flattening filter is 1.53 × 10-3 cGy*min-1 *Ci-1 . With a 14,000 Ci source, the resulting dose rate at the midplane of a typical patient is approximately 17 and 30 cGy/min using the thick and thin flattening filters, respectively, using the maximum source to couch distance. The maximum useful field size, defined by the 90% isodose line, at this location is 225 × 78 cm with field flatness within 5% over the central 178 × 73 cm. Measured output agreed with external validation within 0.5%. End-to-end testing was performed in a modified Rando phantom. In-house MATLAB software was developed to calculate patient-specific dose distributions using DOSXYZnrc, and fabricate custom 3D-printed forms for creating patient-specific lung blocks. End-to-end OSLD and diode measurements both with and without lung blocks agreed with Monte Carlo calculated doses to within 5% and in-phantom film measurements validated dose distribution uniformity. Custom lung block transmission measurements agree well with design criteria and provide clinically favorable dose distributions within the lungs. Block placement is easily facilitated using the flat panel imaging system with an exposure time of 0.01 min. In conclusion, a novel dedicated TBI unit has been commissioned and clinically implemented. Its mechanical, dosimetric, and imaging capabilities are suitable to provide state-of-the-art TBI for patients as large as 220 cm in height and 78 cm in width.

Study of unflattened photon beams shaped by multileaf collimator using BEAMnrc code

Purpose

In our study basic dosimetric properties of a flattening filter free 6 MV photon beam shaped by multileaf collimators (MLC) is examined using the Monte Carlo (MC) method.

Methods and Materials

BEAMnrc code was used to make a MC simulation model for 6 MV photon beam based on Varian Clinic 600 unique performance linac, operated with and without a flattening filter in beam line. Dosimetric features including central axis depth dose, beam profiles, photon and electron spectra were calculated and compared for flattened and unflattened cases.

Results

Dosimetric field size and penumbra were found to be smaller for unflattened beam, and the decrease in field size was less for MLC shaped in comparison with jaw-shaped unflattened beam. Increase in dose rate of >2·4 times was observed for unflattened beam indicating a shorter beam delivery time for treatment. MLC leakage was found to decrease significantly when the flattening filter was removed from the beam line. The total scatter factor showed slower deviation with field sizes for unflattened beam indicating a reduced head scatter.

Conclusions

Our study demonstrated that improved accelerator characteristics can be achieved by removing flattening filter from beam line.

Monte Carlo systems used for treatment planning and dose verification

General-purpose radiation transport Monte Carlo codes have been used for estimation of the absorbed dose distribution in external photon and electron beam radiotherapy patients since several decades. Results obtained with these codes are usually more accurate than those provided by treatment planning systems based on non-stochastic methods. Traditionally, absorbed dose computations based on general-purpose Monte Carlo codes have been used only for research, owing to the difficulties associated with setting up a simulation and the long computation time required. To take advantage of radiation transport Monte Carlo codes applied to routine clinical practice, researchers and private companies have developed treatment planning and dose verification systems that are partly or fully based on fast Monte Carlo algorithms. This review presents a comprehensive list of the currently existing Monte Carlo systems that can be used to calculate or verify an external photon and electron beam radiotherapy treatment plan. Particular attention is given to those systems that are distributed, either freely or commercially, and that do not require programming tasks from the end user. These systems are compared in terms of features and the simulation time required to compute a set of benchmark calculations.

Increasing efficiency of BEAMnrc-simulated Co-60 beams using directional source biasing

PURPOSE

This study describes the implementation of a directional source biasing (DSB) scheme for efficiently simulating Cobalt-60 treatment heads using the BEAMnrc Monte Carlo code. Previous simulation of Co-60 beams with BEAMnrc was impractical because of the time required to track photons not directed into the treatment field and to simulate secondary charged particles.

METHODS

In DSB, efficiency is increased by splitting each photon emitted by the Co-60 source a user-defined number of times. Only those split primary photons directed into a user-defined splitting field (encompassing the treatment field) are sampled, yielding many low-weight photons directed into the field. Efficiency can be further increased by taking advantage of radial symmetry at the top of the treatment head to reduce the number of split primary photons tracked in this portion. There is also an option to generate contaminant electrons in DSB.

RESULTS

The DSB scheme in BEAMnrc increases the photon fluence calculation efficiency in a 10 × 10 cm(2) Co-60 beam by a factor of 1800 with a concurrent increase in contaminant electron fluence calculation efficiency by a factor of 1200. Implementation of DSB in beampp, a C++ code for accelerator simulations based on EGSnrc and the C++ class library, egspp, increases photon fluence efficiency by a factor of 2800 and contaminant electron fluence efficiency by a factor of 1600. Optimum splitting numbers are in the range of 20,000-40,000. For dose calculations in a water phantom (0.5 × 0.5 × 0.5 cm(3) voxels) this translates into a factor of ∼400 increase in dose calculation efficiency (all doses > 0.5 × Dmax). An example calculation of the ratio of dose to water to dose to chamber (the basis of the beam quality correction factor) to within 0.2% in a realistic chamber using a full simulation of a Co-60 treatment head as a source indicates the practicality of Co-60 simulations with DSB.

CONCLUSIONS

The efficiency improvement resulting from DSB makes Monte Carlo commissioning of Co-60 beams and calculation of beam quality correction factors feasible.

Approaching Oxygen-Guided Intensity-Modulated Radiation Therapy

The outcome of cancer radiation treatment is strongly correlated with tumor oxygenation. The aim of this study is to use oxygen tension distributions in tumors obtained using Electron Paramagnetic Resonance (EPR) imaging to devise better tumor radiation treatment. The proposed radiation plan is delivered in two steps. In the first step, a uniform 50% tumor control dose (TCD50) is delivered to the whole tumor. For the second step an additional dose boost is delivered to radioresistant, hypoxic tumor regions. FSa fibrosarcomas grown in the gastrocnemius of the legs of C3H mice were used. Oxygen tension images were obtained using a 250 MHz pulse imager and injectable partially deuterated trityl OX63 (OX71) spin probe. Radiation was delivered with a novel animal intensity modulated radiation therapy (IMRT) XRAD225Cx microCT/radiation therapy delivery system. In a simplified scheme for boost dose delivery, the boost area is approximated by a sphere, whose radius and position are determined using an EPR O2 image. The sphere that irradiates the largest fraction of hypoxic voxels in the tumor was chosen using an algorithm based on Receiver Operator Characteristic (ROC) analysis. We used the fraction of irradiated hypoxic volume as the true positive determinant and the fraction of irradiated normoxic volume as the false positive determinant in the terms of that analysis. The most efficient treatment is the one that demonstrates the shortest distance from the ROC curve to the upper left corner of the ROC plot. The boost dose corresponds to the difference between TCD90 and TCD50 values. For the control experiment an identical radiation dose to the normoxic tumor area is delivered.

Comparative evaluation of modern dosimetry techniques near low- and high-density heterogeneities

The purpose of this study is to compare performance of several dosimetric methods in heterogeneous phantoms irradiated by 6 and 18 MV beams. Monte Carlo (MC) calculations were used, along with two versions of Acuros XB, anisotropic analytical algorithm (AAA), EBT2 film, and MOSkin dosimeters. Percent depth doses (PDD) were calculated and measured in three heterogeneous phantoms. The first two phantoms were a 30×30×30 cm3 solid‐water slab that had an air‐gap of 20×2.5×2.35 cm3. The third phantom consisted of 30×30×5 cm3 solid water slabs, two 30×30×5 cm3 slabs of lung, and one 30×30×1 cm3 solid water slab. Acuros XB, AAA, and MC calculations were within 1% in the regions with particle equilibrium. At media interfaces and buildup regions, differences between Acuros XB and MC were in the range of +4.4% to −12.8%. MOSkin and EBT2 measurements agreed to MC calculations within ∼2.5%, except for the first centimeter of buildup where differences of 4.5% were observed. AAA did not predict the backscatter dose from the high‐density heterogeneity. For the third, multilayer lung phantom, 6 MV beam PDDs calculated by all TPS algorithms were within 2% of MC. 18 MV PDDs calculated by two versions of Acuros XB and AAA differed from MC by up to 2.8%, 3.2%, and 6.8%, respectively. MOSkin and EBT2 each differed from MC by up to 2.9% and 2.5% for the 6 MV, and by −3.1% and ∼2% for the 18 MV beams. All dosimetric techniques, except AAA, agreed within 3% in the regions with particle equilibrium. Differences between the dosimetric techniques were larger for the 18 MV than the 6 MV beam. MOSkin and EBT2 measurements were in a better agreement with MC than Acuros XB calculations at the interfaces, and they were in a better agreement to each other than to MC. The latter is due to their thinner detection layers compared to MC voxel sizes.

PACS numbers: 87.55.K‐, 87.55.kd, 87.55.km, 87.53.Bn, 87.55.k

Dose reduction of scattered photons from concrete walls lined with lead: Implications for improvement in design of megavoltage radiation therapy facility mazes

PURPOSE

This study explores the possibility of using lead to cover part of the radiation therapy facility maze walls in order to absorb low energy photons and reduce the total dose at the maze entrance of radiation therapy rooms.

METHODS

Experiments and Monte Carlo simulations were utilized to establish the possibility of using high-Z materials to cover the concrete walls of the maze in order to reduce the dose of the scatteredphotons at the maze entrance. The dose of the backscatteredphotons from a concrete wall was measured for various scattering angles. The dose was also calculated by the FLUKA and EGSnrc Monte Carlo codes. The FLUKA code was also used to simulate an existing radiotherapy room to study the effect of multiple scattering when adding lead to cover the concrete walls of the maze. Monoenergetic photons were used to represent the main components of the x ray spectrum up to 10 MV.

RESULTS

It was observed that when the concrete wall was covered with just 2 mm of lead, the measured dose rate at all backscattering angles was reduced by 20% for photons of energy comparable to Co-60 emissions and 70% for Cs-137 emissions. The simulations with FLUKA and EGS showed that the reduction in the dose was potentially even higher when lead was added. One explanation for the reduction is the increased absorption of backscatteredphotons due to the photoelectric interaction in lead. The results also showed that adding 2 mm lead to the concrete walls and floor of the maze reduced the dose at the maze entrance by up to 90%.

CONCLUSIONS

This novel proposal of covering part or the entire maze walls with a few millimeters of lead would have a direct implication for the design of radiation therapy facilities and would assist in upgrading the design of some mazes, especially those in facilities with limited space where the maze length cannot be extended to sufficiently reduce the dose.

ARCHERRT - a GPU-based and photon-electron coupled Monte Carlo dose computing engine for radiation therapy: software development and application to helical tomotherapy

PURPOSE

Using the graphical processing units (GPU) hardware technology, an extremely fast Monte Carlo (MC) code ARCHERRT is developed for radiation dose calculations in radiation therapy. This paper describes the detailed software development and testing for three clinical TomoTherapy® cases: the prostate, lung, and head & neck.

METHODS

To obtain clinically relevant dose distributions, phase space files (PSFs) created from optimized radiation therapy treatment plan fluence maps were used as the input to ARCHERRT. Patient-specific phantoms were constructed from patient CT images. Batch simulations were employed to facilitate the time-consuming task of loading large PSFs, and to improve the estimation of statistical uncertainty. Furthermore, two different Woodcock tracking algorithms were implemented and their relative performance was compared. The dose curves of an Elekta accelerator PSF incident on a homogeneous water phantom were benchmarked against DOSXYZnrc. For each of the treatment cases, dose volume histograms and isodose maps were produced from ARCHERRT and the general-purpose code, GEANT4. The gamma index analysis was performed to evaluate the similarity of voxel doses obtained from these two codes. The hardware accelerators used in this study are one NVIDIA K20 GPU, one NVIDIA K40 GPU, and six NVIDIA M2090 GPUs. In addition, to make a fairer comparison of the CPU and GPU performance, a multithreaded CPU code was developed using OpenMP and tested on an Intel E5-2620 CPU.

RESULTS

For the water phantom, the depth dose curve and dose profiles from ARCHERRT agree well with DOSXYZnrc. For clinical cases, results from ARCHERRT are compared with those from GEANT4 and good agreement is observed. Gamma index test is performed for voxels whose dose is greater than 10% of maximum dose. For 2%/2mm criteria, the passing rates for the prostate, lung case, and head & neck cases are 99.7%, 98.5%, and 97.2%, respectively. Due to specific architecture of GPU, modified Woodcock tracking algorithm performed inferior to the original one. ARCHERRT achieves a fast speed for PSF-based dose calculations. With a single M2090 card, the simulations cost about 60, 50, 80 s for three cases, respectively, with the 1% statistical error in the PTV. Using the latest K40 card, the simulations are 1.7-1.8 times faster. More impressively, six M2090 cards could finish the simulations in 8.9-13.4 s. For comparison, the same simulations on Intel E5-2620 (12 hyperthreading) cost about 500-800 s.

CONCLUSIONS

ARCHERRT was developed successfully to perform fast and accurate MC dose calculation for radiotherapy using PSFs and patient CT phantoms.

A virtual source model for Monte Carlo simulation of helical tomotherapy

The purpose of this study was to present a Monte Carlo (MC) simulation method based on a virtual source, jaw, and MLC model to calculate dose in patient for helical tomotherapy without the need of calculating phase‐space files (PSFs). Current studies on the tomotherapy MC simulation adopt a full MC model, which includes extensive modeling of radiation source, primary and secondary jaws, and multileaf collimator (MLC). In the full MC model, PSFs need to be created at different scoring planes to facilitate the patient dose calculations. In the present work, the virtual source model (VSM) we established was based on the gold standard beam data of a tomotherapy unit, which can be exported from the treatment planning station (TPS). The TPS‐generated sinograms were extracted from the archived patient XML (eXtensible Markup Language) files. The fluence map for the MC sampling was created by incorporating the percentage leaf open time (LOT) with leaf filter, jaw penumbra, and leaf latency contained from sinogram files. The VSM was validated for various geometry setups and clinical situations involving heterogeneous media and delivery quality assurance (DQA) cases. An agreement of <1% was obtained between the measured and simulated results for percent depth doses (PDDs) and open beam profiles for all three jaw settings in the VSM commissioning. The accuracy of the VSM leaf filter model was verified in comparing the measured and simulated results for a Picket Fence pattern. An agreement of <2% was achieved between the presented VSM and a published full MC model for heterogeneous phantoms. For complex clinical head and neck (HN) cases, the VSM‐based MC simulation of DQA plans agreed with the film measurement with 98% of planar dose pixels passing on the 2%/2 mm gamma criteria. For patient treatment plans, results showed comparable dose‐volume histograms (DVHs) for planning target volumes (PTVs) and organs at risk (OARs). Deviations observed in this study were consistent with literature. The VSM‐based MC simulation approach can be feasibly built from the gold standard beam model of a tomotherapy unit. The accuracy of the VSM was validated against measurements in homogeneous media, as well as published full MC model in heterogeneous media.

PACS numbers: 87.53.‐j, 87.55.K‐

A hybrid phase-space and histogram source model for GPU-based Monte Carlo radiotherapy dose calculation

In GPU-based Monte Carlo simulations for radiotherapy dose calculation, source modelling from a phase-space source can be an efficiency bottleneck. Previously, this has been addressed using phase-space-let (PSL) sources, which provided significant efficiency enhancement. We propose that additional speed-up can be achieved through the use of a hybrid primary photon point source model combined with a secondary PSL source. A novel phase-space derived and histogram-based implementation of this model has been integrated into gDPM v3.0. Additionally, a simple method for approximately deriving target photon source characteristics from a phase-space that does not contain inheritable particle history variables (LATCH) has been demonstrated to succeed in selecting over 99% of the true target photons with only ~0.3% contamination (for a Varian 21EX 18 MV machine). The hybrid source model was tested using an array of open fields for various Varian 21EX and TrueBeam energies, and all cases achieved greater than 97% chi-test agreement (the mean was 99%) above the 2% isodose with 1% / 1 mm criteria. The root mean square deviations (RMSDs) were less than 1%, with a mean of 0.5%, and the source generation time was 4-5 times faster. A seven-field intensity modulated radiation therapy patient treatment achieved 95% chi-test agreement above the 10% isodose with 1% / 1 mm criteria, 99.8% for 2% / 2 mm, a RMSD of 0.8%, and source generation speed-up factor of 2.5.

Development and reproducibility evaluation of a Monte Carlo-based standard LINAC model for quality assurance of multi-institutional clinical trials

Technical developments in radiotherapy (RT) have created a need for systematic quality assurance (QA) to ensure that clinical institutions deliver prescribed radiation doses consistent with the requirements of clinical protocols. For QA, an ideal dose verification system should be independent of the treatment-planning system (TPS). This paper describes the development and reproducibility evaluation of a Monte Carlo (MC)-based standard LINAC model as a preliminary requirement for independent verification of dose distributions. The BEAMnrc MC code is used for characterization of the 6-, 10- and 15-MV photon beams for a wide range of field sizes. The modeling of the LINAC head components is based on the specifications provided by the manufacturer. MC dose distributions are tuned to match Varian Golden Beam Data (GBD). For reproducibility evaluation, calculated beam data is compared with beam data measured at individual institutions. For all energies and field sizes, the MC and GBD agreed to within 1.0% for percentage depth doses (PDDs), 1.5% for beam profiles and 1.2% for total scatter factors (Scps.). Reproducibility evaluation showed that the maximum average local differences were 1.3% and 2.5% for PDDs and beam profiles, respectively. MC and institutions' mean Scps agreed to within 2.0%. An MC-based standard LINAC model developed to independently verify dose distributions for QA of multi-institutional clinical trials and routine clinical practice has proven to be highly accurate and reproducible and can thus help ensure that prescribed doses delivered are consistent with the requirements of clinical protocols.

Monte Carlo study of radiation dose enhancement by gadolinium in megavoltage and high dose rate radiotherapy

MRI is often used in tumor localization for radiotherapy treatment planning, with gadolinium (Gd)-containing materials often introduced as a contrast agent. Motexafin gadolinium is a novel radiosensitizer currently being studied in clinical trials. The nanoparticle technologies can target tumors with high concentration of high-Z materials. This Monte Carlo study is the first detailed quantitative investigation of high-Z material Gd-induced dose enhancement in megavoltage external beam photon therapy. BEAMnrc, a radiotherapy Monte Carlo simulation package, was used to calculate dose enhancement as a function of Gd concentration. Published phase space files for the TrueBeam flattening filter free (FFF) and conventional flattened 6MV photon beams were used. High dose rate (HDR) brachytherapy with Ir-192 source was also investigated as a reference. The energy spectra difference caused a dose enhancement difference between the two beams. Since the Ir-192 photons have lower energy yet, the photoelectric effect in the presence of Gd leads to even higher dose enhancement in HDR. At depth of 1.8 cm, the percent mean dose enhancement for the FFF beam was 0.38±0.12, 1.39±0.21, 2.51±0.34, 3.59±0.26, and 4.59±0.34 for Gd concentrations of 1, 5, 10, 15, and 20 mg/mL, respectively. The corresponding values for the flattened beam were 0.09±0.14, 0.50±0.28, 1.19±0.29, 1.68±0.39, and 2.34±0.24. For Ir-192 with direct contact, the enhanced were 0.50±0.14, 2.79±0.17, 5.49±0.12, 8.19±0.14, and 10.80±0.13. Gd-containing materials used in MRI as contrast agents can also potentially serve as radiosensitizers in radiotherapy. This study demonstrates that Gd can be used to enhance radiation dose in target volumes not only in HDR brachytherapy, but also in 6 MV FFF external beam radiotherapy, but higher than the currently used clinical concentration (>5 mg/mL) would be needed.

An automated voxelized dosimetry tool for radionuclide therapy based on serial quantitative SPECT/CT imaging

PURPOSE

To create an accurate map of the distribution of radiation dose deposition in healthy and target tissues during radionuclide therapy.

METHODS

Serial quantitative SPECT∕CT images were acquired at 4, 24, and 72 h for 28 (177)Lu-octreotate peptide receptor radionuclide therapy (PRRT) administrations in 17 patients with advanced neuroendocrine tumors. Deformable image registration was combined with an in-house programming algorithm to interpolate pharmacokinetic uptake and clearance at a voxel level. The resultant cumulated activity image series are comprised of values representing the total number of decays within each voxel's volume. For PRRT, cumulated activity was translated to absorbed dose based on Monte Carlo-determined voxel S-values at a combination of long and short ranges. These dosimetric image sets were compared for mean radiation absorbed dose to at-risk organs using a conventional MIRD protocol (OLINDA 1.1).

RESULTS

Absorbed dose values to solid organs (liver, kidneys, and spleen) were within 10% using both techniques. Dose estimates to marrow were greater using the voxelized protocol, attributed to the software incorporating crossfire effect from nearby tumor volumes.

CONCLUSIONS

The technique presented offers an efficient, automated tool for PRRT dosimetry based on serial post-therapy imaging. Following retrospective analysis, this method of high-resolution dosimetry may allow physicians to prescribe activity based on required dose to tumor volume or radiation limits to healthy tissue in individual patients.

Dosimetric impact of gold markers implanted closely to lung tumors: a Monte Carlo simulation

We are developing an innovative dynamic tumor tracking irradiation technique using gold markers implanted around a tumor as a surrogate signal, a real‐time marker detection system, and a gimbaled X‐ray head in the Vero4DRT. The gold markers implanted in a normal organ will produce uncertainty in the dose calculation during treatment planning because the photon mass attenuation coefficient of a gold marker is much larger than that of normal tissue. The purpose of this study was to simulate the dose variation near the gold markers in a lung irradiated by a photon beam using the Monte Carlo method. First, the single‐beam and the opposing‐beam geometries were simulated using both water and lung phantoms. Subsequently, the relative dose profiles were calculated using a stereotactic body radiotherapy (SBRT) treatment plan for a lung cancer patient having gold markers along the anteriorposterior (AP) and right‐left (RL) directions. For the single beam, the dose at the gold marker‐phantom interface laterally along the perpendicular to the beam axis increased by a factor of 1.35 in the water phantom and 1.58 in the lung phantom, respectively. Furthermore, the entrance dose at the interface along the beam axis increased by a factor of 1.63 in the water phantom and 1.91 in the lung phantom, while the exit dose increased by a factor of 1.00 in the water phantom and 1.12 in the lung phantom, respectively. On the other hand, both dose escalations and dose de‐escalations were canceled by each beam for opposing portal beams with the same beam weight. For SBRT patient data, the dose at the gold marker edge located in the tumor increased by a factor of 1.30 in both AP and RL directions. In clinical cases, dose escalations were observed at the small area where the distance between a gold marker and the lung tumor was ≤ 5 mm, and it would be clinically negligible in multibeam treatments, although further investigation may be required.

PACS number: 87.10.Rt

Red bone marrow dose calculations in radiotherapy of prostate cancer based on the updated VCH adult male phantom

Red bone marrow (RBM) is an important dose-limiting tissue that has high radiosensitivity but is difficult to identify on clinical medical images. In this study, we investigated dose distribution in RBM for prostate cancer radiotherapy. Four suborgans were identified in the skeleton of the visible Chinese human phantom: cortical bone (CB), trabecular bone (TB), RBM, and yellow bone marrow (YBM). Dose distributions in the phantom were evaluated by the Monte Carlo method. When the left os coxae was taken as the organ-at-risk (OAR), the difference in absorbed dose between RBM and each CB and TB was up to 20%, but was much less (≤3.1%) between RBM and YBM. When the left os coxae and entire bone were both taken as OARs, RBM dose also increased with increasing planning target volume size. The results indicate the validity of using dose to homogeneous bone marrow mixture for estimating dose to RBM when RBM is not available in computational phantoms. In addition, the human skeletal system developed in this study provides a model for considering RBM dose in radiotherapy planning.

Monte Carlo comparison of superficial dose between flattening filter free and flattened beams

This study investigates the superficial dose from FFF beams in comparison with the conventional flattened ones using a Monte Carlo (MC) method. Published phase-space files which incorporated real geometry of a TrueBeam accelerator were used for the dose calculation in phantom and clinical cases. The photon fluence on the central axis is 3 times that of a flattened beam for a 6 MV FFF beam and 5 times for a 10 MV beam. The mean energy across the field in air at the phantom surface is 0.92-0.95 MeV for the 6 MV FFF beam and 1.18-1.30 MeV for the corresponding flattened beam. At 10 MV, the values are 1.52-1.72 and 2.15-2.87 MeV for the FFF and flattened beams, respectively. The phantom dose at the depth of 1 mm in the 6 MV FFF beam is 6% ± 2.5% (of the maximum dose) higher compared to the flattened beam for a 25 × 25 cm(2) field and 14.6% ± 1.9% for the 2 × 2 cm(2) field. For the 10 MV beam, the corresponding differences are 3.4% ± 1.5% and 10.7% ± 0.6%. The skin dose difference at selected points on the patient's surface between the plans using FFF and flattened beams in the head-and-neck case was 6.5% ± 2.3% (1SD), and for the breast case it was 6.4% ± 2.3%. The Monte Carlo simulations showed that due to the lower mean energy in the FFF beam, the clinical superficial dose is higher without the flattening filter compared to the flattened beam.

Commissioning of 6 MV medical linac for dynamic MLC-based IMRT on Monte Carlo code GEANT4

Monte Carlo simulation is the most accurate tool for calculating dose distributions. In particular, the Electron Gamma shower computer code has been widely used for multi-purpose research in radiotherapy, but Monte Carlo GEANT4 (GEometry ANd Tracking) is rare for radiotherapy with photon beams and needs to be verified further under various irradiation conditions, particularly multi-leaf collimator-based intensity-modulated radiation therapy (MLC-based IMRT). In this study, GEANT4 was used for modeling of a 6 MV linac for dynamic MLC-based IMRT. To verify the modeling of our linac, we compared the calculated data with the measured depth-dose for a 10 × 10 cm(2) field and the measured dose profile for a 35 × 35 cm(2) field. Moreover, 120 MLCs were modeled on the GEANT4. Five tests of MLC modeling were performed: (I) MLC transmission, (II) MLC transmission profile including intra- and inter-leaf leakage, (III) tongue-and-groove leakage, (IV) a simple field with different field sizes by use of MLC and (V) a dynamic MLC-based IMRT field. For all tests, the calculations were compared with measurements of an ionization chamber and radiographic film. The calculations agreed with the measurements: MLC transmissions by calculations and measurements were 1.76 ± 0.01 and 1.87 ± 0.01 %, respectively. In gamma evaluation method (3 %/3 mm), the pass rates of the (IV) and (V) tests were 98.5 and 97.0 %, respectively. Furthermore, tongue-and-groove leakage could be calculated by GEANT4, and it agreed with the film measurements. The procedure of commissioning of dynamic MLC-based IMRT for GEANT4 is proposed in this study.

GammaPod-a new device dedicated for stereotactic radiotherapy of breast cancer

PURPOSE

This paper introduces a new external beam radiotherapy device named GammaPod that is dedicated for stereotactic radiotherapy of breast cancer.

METHODS

The design goal of the GammaPod as a dedicated system for treating breast cancer is the ability to deliver ablative doses with sharp gradients under stereotactic image guidance. Stereotactic localization of the breast is achieved by a vacuum-assisted breast immobilization cup with built-in stereotactic frame. Highly focused radiation is achieved at the isocenter due to the cross-firing from 36 radiation arcs generated by rotating 36 individual Cobalt-60 beams. The dedicated treatment planning system optimizes an optimal path of the focal spot using an optimization algorithm borrowed from computational geometry such that the target can be covered by 90%-95% of the prescription dose and the doses to surrounding tissues are minimized. The treatment plan is intended to be delivered with continuous motion of the treatment couch. In this paper the authors described in detail the gamma radiation unit, stereotactic localization of the breast, and the treatment planning system of the GammaPod system.

RESULTS

A prototype GammaPod system was installed at University of Maryland Medical Center and has gone through a thorough functional, geometric, and dosimetric testing. The mechanical and functional performances of the system all meet the functional specifications.

CONCLUSIONS

An image-guided breast stereotactic radiotherapy device, named GammaPod, has been developed to deliver highly focused and localized doses to a target in the breast under stereotactic image guidance. It is envisioned that the GammaPod technology has the potential to significantly shorten radiation treatments and even eliminate surgery by ablating the tumor and sterilizing the tumor bed simultaneously.