Design and characterisation of a minibeam collimator utilising Monte Carlo simulation and a clinical linear accelerator

Objective.Spatially fractionated radiotherapy is showing promise as a treatment modality. Initial focus was on beams of photons at low energy produced from a synchrotron but more recently research has expanded to include applications in proton therapy. Interest in photon beams remains and this is the focus of this paperApproach.This study presents a 3D printed tungsten minibeam collimator intended to produce peak-to-valley dose ratios (PVDR) of between seven and ten with a 1 MV, bremsstrahlung generated, photon beam. The design of the collimator is motivated by a Monte Carlo study estimating the PVDR for different collimator designs at different energies. This collimator was characterised on a clinical linear accelerator (Elekta VersaHD) as well as an orthovoltage unit.Main results.The performance of the fabricated collimator was measured on Elekta VersaHD running in unflattened mode with a 6 MV beam. On the Elekta VersaHD units the PVDR was measured to be between approximately 1.5 and 2.0 at 3 cm deep. For measurements with the orthovoltage unit PVDRs of greater than 10 were observed at a depth of 4 cm.Significance.The results confirmed that the predictions from simulation could be reproduced on linear accelerators currently in clinical usage, producing PVDRs between 2-2.5. Using the model to predict PVDRs using 1 MV photon beams, the threshold considered to produce enhanced normal tissue dose tolerance (>7) was surpassed. This suggests the possibility of using such techniques with versions of existing Linac technology which have been modified to operate at low energy and high beam currents.

Monte-Carlo calculation of output correction factors for ionization chambers, solid-state detectors, and EBT3 film in small fields of high-energy photons

High-energy accelerators are often used in oncological practice, but the information on the small-field dosimetry for the photon beams with nominal energy above 10 MV is limited. The goal of the present work was to determine the values of the output correction factor ( k Q clin , Q ref f clin , f ref $k_{{Q}_{{\rm{clin}}},{Q}_{{\rm{ref}}}}^{{f}_{{\rm{clin}}},{f}_{{\rm{ref}}}}$ ) for solid-state detectors (Diode E, PTW 60017; microDiamond, PTW 60019), EBT3 film, and ionization chambers (Semiflex, PTW 31010; Semiflex 3D, PTW 31021; PinPoint, PTW 31015; PinPoint 3D, PTW 31016) in the small fields formed by 10, 15, 18, and 20 MV photon beams. The output correction factors were calculated by Monte-Carlo method using EGSnrc toolkit for six field sizes (from 0.5 × 0.5 cm 2 $0.5 \times 0.5\ {\rm{cm}}^2$ to 10 × 10 cm 2 $10 \times 10\ {\rm{cm}}^2$ ) for isocentric and constant source-to-surface distance (SSD) techniques. The decrease in the field size led to an increase in k Q clin , Q ref f clin , f ref $k_{{Q}_{{\rm{clin}}},{Q}_{{\rm{ref}}}}^{{f}_{{\rm{clin}}},{f}_{{\rm{ref}}}}$ for ionization chambers, while for solid-state detectors and radiochromic film, k Q clin , Q ref f clin , f ref $k_{{Q}_{{\rm{clin}}},{Q}_{{\rm{ref}}}}^{{f}_{{\rm{clin}}},{f}_{{\rm{ref}}}}$ were less than unity at the smallest field size. A larger sensitive volume of ionization chamber corresponded to a stronger deviation of output correction factor from unity: 1.847 (125 mm³ PTW 31010) versus up to 1.183 (16 mm³ PTW 31016) at the smallest field of 10 MV beam. The calculated output correction factors were used to correct the output factors for PTW 60017, PTW 60019, and EBT3. The deviation of the corrected output factor from the results of Monte-Carlo simulation did not exceed 3% in the fields from 1.0 × 1.0 cm 2 $1.0 \times 1.0\ {\rm{cm}}^2$ to 4.0 × 4.0 cm 2 $4.0 \times 4.0\ {\rm{cm}}^2$ for 10 and 18 MV beams. Thus, Diode E, microDiamond, and EBT3 film can be recommended for small-field dosimetry of high-energy photons.

PENELOPE simulations and experiment for 6 MV clinac iX accelerator for standard and small static fields

The goal of this work was to produce accurate data for use as a 'gold standard' and a valid tool for measurements in reference dosimetry for standard/small static field sizes from 0.5 × 0.5 to 10 × 10 cm². It is based on the accuracy of the phase space files (PSFs) as a key quantity. Because the IAEA general public database provides few PSFs for the Varian iX, we simulated the head through Monte Carlo (MC) simulations and calculated validated PSFs for 12 square field sizes including seven for small static fields. The resulting dosimetric calculations allowed us to reach a good level of agreement in comparison to our relative and absolute dose measurements performed on a Varian iX in water phantom. Measured and MC calculated output factors were investigated for different detectors. Based on the TRS 483 formalism and MC (PENELOPE/penEasy), we calculated output correction factors for the unshielded Diode-E (T60017) and the PinPoint-3D (T31016) micro-chamber according to manufacturers' blueprints. Our MC results were in agreement with the recommended data; they compete with recent measurements and MC simulations and in particular the TRS 483 MC data obtained from similar simulations. Moreover, our MC results provide supplemental data in comparison to TRS 483 data in particular for the PinPoint-3D (T31016). We suggest our MC output correction factors as new datasets for future TRS compilations. The work was substantial, used different robust MC strategies depending on the scoring regions, and led in most cases to uncertainties of less than 1%.

Reference dosimetry of modulated and dynamic photon beams

In the late 1980s, a new technique was proposed that would revolutionize radiotherapy. Now referred to as intensity-modulated radiotherapy, it is at the core of state-of-the-art photon beam delivery techniques, such as helical tomotherapy and volumetric modulated arc therapy. Despite over two decades of clinical application, there are still no established guidelines on the calibration of dynamic modulated photon beams. In 2008, the IAEA-AAPM work group on nonstandard photon beam dosimetry published a formalism to support the development of a new generation of protocols applicable to nonstandard beam reference dosimetry (Alfonso et al 2008 Med. Phys. 35 5179-86). The recent IAEA Code of Practice TRS-483 was published as a result of this initiative and addresses exclusively small static beams. But the plan-class specific reference calibration route proposed by Alfonso et al (2008 Med. Phys. 35 5179-86) is a change of paradigm that is yet to be implemented in radiotherapy clinics. The main goals of this paper are to provide a literature review on the dosimetry of nonstandard photon beams, including dynamic deliveries, and to discuss anticipated benefits and challenges in a future implementation of the IAEA-AAPM formalism on dynamic photon beams.

Reconstruction of the electron source intensity distribution of a clinical linear accelerator using in-air measurements and a genetic algorithm

Background and purpose

The electron source intensity distribution of a clinical linear accelerator has a great influence on the calculation of output factors for small radiation fields where source occlusion by the collimating devices takes place. The purpose of this study was to present a new method for the electron source reconstruction problem.

Materials and methods

The measurements were performed in-air using diode and 6 MV 1 × 1 cm² photon field in flattening filter-free mode. In Monte Carlo simulation, an electron target area was divided into a number of square subsources. Then, the in-air doses in 2D silicon chip array were calculated individually from each subsource. A genetic algorithm search was applied in order to determine the optimal weight factors for all subsources that provide the best agreement between simulated and measured doses.

Results

It was found that the reconstructed electron source intensity from a clinical linear accelerator has the two-dimensional elliptical double Gaussian distribution. The source intensity distribution consisted of two intensity components along the in-plane (x) and cross-plane (y) directions characterized by full width half-maximum (FWHM): FWHM_x1 = 0.27 cm, FWHM_x2 = 0.08 cm, FWHM_y1 = 0.24 cm, FWHM_y2 = 0.06 cm, where broader components are 81% and 53% of the total intensity along × and y axis respectively.

Conclusions

The obtained results demonstrated an elliptical double Gaussian intensity distribution of the incident electron source. We anticipate that the proposed method has universal applications independent of the type of linear accelerator, modality or energy.

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.

Hybrid Monte Carlo source model: Advantages and deficiencies

Monte Carlo (MC) simulation is the gold standard for dose calculation. An accurate mathematical source model can be used for the radiation beams. Source models can consist of sub-sources or fewer sources with data that need to be measured. This can speed up treatment plan verification without the need for a full simulation of the radiation treatment machine.

Aims: This study aimed to construct a novel hybrid source model for 6 MV photon beams for an Elekta Synergy accelerator and to commission it against measured beam data and treatments plans.

Methods and Material: The model comprised of a circular photon and planar electron contamination source. The modified Schiff formula provided off-axis variable bremsstrahlung spectra. Collimation and scatter were modelled with error functions. An exponential function modelled the transmitted fluence through the collimators. The source model was commissioned by comparing simulated and measured MC data. Dose data included profiles, depth dose and film measurements in a Rando phantom. Field sizes ranged from 1 × 1 cm2 to 40 × 40 cm2.

Results: Regular, wedged and asymmetrical fields could be modelled within 1.5% or 1.5 mm. More than 95% of all points lie within 3% or 3 mm for the multi-leaf collimators contours data. A gamma criterion of 3% or 3 mm was met for a complex treatment case.

Conclusions: The two sub-source model replicated clinical 6 MV Elekta Synergy photons beams and could calculate the dose accurately for conformal treatments in complex geometries such as a head-and-neck case.

Experimental Validation of Monte Carlo Simulations Based on a Virtual Source Model for TomoTherapy in a RANDO Phantom

A virtual source model for Monte Carlo simulations of helical TomoTherapy has been developed previously by the authors. The purpose of this work is to perform experiments in an anthropomorphic (RANDO) phantom with the same order of complexity as in clinical treatments to validate the virtual source model to be used for quality assurance secondary check on TomoTherapy patient planning dose. Helical TomoTherapy involves complex delivery pattern with irregular beam apertures and couch movement during irradiation. Monte Carlo simulation, as the most accurate dose algorithm, is desirable in radiation dosimetry. Current Monte Carlo simulations for helical TomoTherapy adopt the full Monte Carlo model, which includes detailed modeling of individual machine component, and thus, large phase space files are required at different scoring planes. As an alternative approach, we developed a virtual source model without using the large phase space files for the patient dose calculations previously. In this work, we apply the simulation system to recompute the patient doses, which were generated by the treatment planning system in an anthropomorphic phantom to mimic the real patient treatments. We performed thermoluminescence dosimeter point dose and film measurements to compare with Monte Carlo results. Thermoluminescence dosimeter measurements show that the relative difference in both Monte Carlo and treatment planning system is within 3%, with the largest difference less than 5% for both the test plans. The film measurements demonstrated 85.7% and 98.4% passing rate using the 3 mm/3% acceptance criterion for the head and neck and lung cases, respectively. Over 95% passing rate is achieved if 4 mm/4% criterion is applied. For the dose–volume histograms, very good agreement is obtained between the Monte Carlo and treatment planning system method for both cases. The experimental results demonstrate that the virtual source model Monte Carlo system can be a viable option for the accurate dose calculation of helical TomoTherapy.

Implementation of a double Gaussian source model for the BEAMnrc Monte Carlo code and its influence on small fields dose distributions

The shape of the radiation source of a linac has a direct impact on the delivered dose distributions, especially in the case of small radiation fields. Traditionally, a single Gaussian source model is used to describe the electron beam hitting the target, although different studies have shown that the shape of the electron source can be better described by a mixed distribution consisting of two Gaussian components. Therefore, this study presents the implementation of a double Gaussian source model into the BEAMnrc Monte Carlo code. The impact of the double Gaussian source model for a 6 MV beam is assessed through the comparison of different dosimetric parameters calculated using a single Gaussian source, previously commissioned, the new double Gaussian source model and measurements, performed with a diode detector in a water phantom. It was found that the new source can be easily implemented into the BEAMnrc code and that it improves the agreement between measurements and simulations for small radiation fields. The impact of the change in source shape becomes less important as the field size increases and for increasing distance of the collimators to the source, as expected. In particular, for radiation fields delivered using stereotactic collimators located at a distance of 59 cm from the source, it was found that the effect of the double Gaussian source on the calculated dose distributions is negligible, even for radiation fields smaller than 5 mm in diameter. Accurate determination of the shape of the radiation source allows us to improve the Monte Carlo modeling of the linac, especially for treatment modalities such as IMRT, were the radiation beams used could be very narrow, becoming more sensitive to the shape of the source.

PACS number(s): 87.53.Bn, 87.55.K, 87.56.B‐, 87.56.jf

Direct reconstruction of the source intensity distribution of a clinical linear accelerator using a maximum likelihood expectation maximization algorithm

Direct determination of the source intensity distribution of clinical linear accelerators is still a challenging problem for small field beam modeling. Current techniques most often involve special equipment and are difficult to implement in the clinic. In this work we present a maximum-likelihood expectation-maximization (MLEM) approach to the source reconstruction problem utilizing small fields and a simple experimental set-up. The MLEM algorithm iteratively ray-traces photons from the source plane to the exit plane and extracts corrections based on photon fluence profile measurements. The photon fluence profiles were determined by dose profile film measurements in air using a high density thin foil as build-up material and an appropriate point spread function (PSF). The effect of other beam parameters and scatter sources was minimized by using the smallest field size ([Formula: see text] cm(2)). The source occlusion effect was reproduced by estimating the position of the collimating jaws during this process. The method was first benchmarked against simulations for a range of typical accelerator source sizes. The sources were reconstructed with an accuracy better than 0.12 mm in the full width at half maximum (FWHM) to the respective electron sources incident on the target. The estimated jaw positions agreed within 0.2 mm with the expected values. The reconstruction technique was also tested against measurements on a Varian Novalis Tx linear accelerator and compared to a previously commissioned Monte Carlo model. The reconstructed FWHM of the source agreed within 0.03 mm and 0.11 mm to the commissioned electron source in the crossplane and inplane orientations respectively. The impact of the jaw positioning, experimental and PSF uncertainties on the reconstructed source distribution was evaluated with the former presenting the dominant effect.

Rounded leaf end effect of multileaf collimator on penumbra width and radiation field offset: an analytical and numerical study

Background

Penumbra characteristics play a significant role in dose delivery accuracy for radiation therapy. For treatment planning, penumbra width and radiation field offset strongly influence target dose conformity and organ at risk sparing.

Methods

In this study, we present an analytical and numerical approach for evaluation of the rounded leaf end effect on penumbra characteristics. Based on the rule of half-value layer, algorithms for leaf position calculation and radiation field offset correction were developed, which were advantageous particularly in dealing with large radius leaf end. Computer simulation was performed based on the Monte Carlo codes of EGSnrc/BEAMnrc, with groups of leaf end radii and source sizes. Data processing technique of curve fitting was employed for deriving penumbra width and radiation field offset.

Results

Results showed that penumbra width increased with source size. Penumbra width curves for large radius leaf end were U-shaped. This observation was probably related to the fact that radiation beams penetrated through the proximal and distal leaf sides. In contrast, source size had negligible impact on radiation field offset. Radiation field offsets were found to be constant both for analytical method and numerical simulation. However, the overall resulting values of radiation field offset obtained by analytical method were slightly smaller compared with Monte Carlo simulation.

Conclusions

The method we proposed could provide insight into the investigation of rounded leaf end effects on penumbra characteristics. Penumbra width and radiation field offset calibration should be carefully performed to commission multileaf collimator for intensity modulated radiotherapy.

Generation of prescriptions robust against geometric uncertainties in dose painting by numbers

BACKGROUND

In the context of dose painting by numbers delivered with intensity-modulated radiotherapy, the robustness of dose distributions against geometric uncertainties can be ensured by robust optimization. As robust optimization is seldom available in treatment planning systems (TPS), we propose an alternative method that reaches the same goal by modifying the heterogeneous dose prescription (based on (18)FDG-PET) and guarantees coverage in spite of systematic and random errors with known standard deviations Σ and σ, respectively.

MATERIAL AND METHODS

The objective was that 95% of all voxels in the GTVPET received at least 95% of the prescribed dose despite geometric errors. The prescription was modified by a geometric dilation of αΣ for systematic errors and a deconvolution by a Gaussian function of width σ for random errors. For a 90% confidence interval, α = 2.5. Planning was performed on a TomoTherapy system, such that 95% of the voxels received at least 95% of the modified prescription and less than 5% of the voxels received more than 105% of the modified prescription. The applicability of the method was illustrated for two head-and-neck tumors.

RESULTS

Systematic and random displacements larger than αΣ and σ degraded coverage. Down to 62.8% of the points received at least 95% of prescribed dose for the largest considered displacements (5 mm systematic translation and 3 mm standard deviation for random errors). When systematic and random displacements were smaller than αΣ and σ, no degradation of target coverage was observed.

CONCLUSIONS

The method led to treatment plans with target coverage robust against geometric uncertainties without the need to incorporate these in the optimizer of the TPS. The methodology was illustrated for head-and-neck cancer but can be potentially extended to all treatment sites.

Angle dependent focal spot size of a conical X-ray target

Misaligned phantoms may severely affect the focal spot calculations. A method is proposed to determine the geometry of the X-ray target and the position of the image radiograph around the X-ray target to get a relatively smaller focal spot size. Results reveal that the focal spot size is not always isotropic around the target but it decreases as the point of observation shifts radially away from the center line of the conical X-ray target. This research will help in producing high quality X-ray images in multi-directions by properly aligning the phantoms and the radiograph tallies.

A slit method to determine the focal spot size and shape of TomoTherapy system

PURPOSE

To obtain accurate x-ray source profile measurements using a slit-collimator, the slit-collimator should have a narrow width, large height, and be positioned near the source. However, these conditions may not always be met. In this paper, the authors provide a detailed analysis of the slit measurement geometry and the relationship between the slit parameters and the measured x-ray source profile. The slit model allows the use of a shorter and more easily available slit-collimator, while accurate source profile measurements can still be obtained.

METHODS

Measurements were performed with a variety of slit widths and/or slit to source distances. The relationship derived between the slit parameters and the measured profile was used to determine the true focal spot profile through a least square fit of the profile data. The model was verified by comparing the predicted profiles at a variety of slit-collimator parameters with the measured results on the TomoTherapy Hi-Art system.

RESULTS

Both the treatment beam and the imaging beam were measured. For treatment mode, it was found that a source consisting of one Gaussian with a 0.75 mm full-width-half-maximum (FWHM) and 72% peak amplitude and a second Gaussian with a 2.27 mm FWHM and 18% peak amplitude matched measurement profiles. The overall source profile has a FWHM of 0.93 mm, but with a higher amplitude in the tail region than a single Gaussian. For imaging mode, the source consists of one Gaussian with a 0.68 mm FWHM and 82% peak amplitude and a second Gaussian with a 1.83 mm FWHM and 18% peak amplitude. The overall source profile has a FWHM of 0.77 mm.

CONCLUSIONS

Our study of the focal spot measurement using slit-collimators showed that accurate source profile measurements can be achieved through fitting of measurement results at different slit widths and source-to-slit distances (SSD). Quantitative measurements of the TomoTherapy linac focal spot showed that the source distribution could be better described with a model consisting of two Gaussian components rather than a single Gaussian model as assumed in previous studies.