The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

Pixel value to electron density conversion for a Megavoltage Cone Beam CT System

Introduction: The motivation for developing megavoltage (and kilovoltage) cone beam CT (MV CBCT) capabilities in the radiotherapy treatment room was primarily based on the need to improve patient set-up accuracy. There has recently been an interest in using the cone beam CT data for treatment planning. Accurate treatment planning, however, requires knowledge of the electron density of the tissues receiving radiation in order to calculate dose distributions. This is obtained from CT, utilising a conversion between CT number and electron density of various tissues. The use of MV CBCT has particular advantages compared to treatment planning with kilovoltage CT in the presence of high atomic number materials and requires the conversion of pixel values from the image sets to electron density. Therefore, a study was undertaken to characterise the pixel value to electron density relationship for the Siemens MV CBCT system, MVision, and determine the effect, if any, of differing the number of monitor units used for acquisition. If a significant difference with number of monitor units was seen then pixel value to ED conversions may be required for each of the clinical settings. The calibration of the MV CT images for electron density offers the possibility for a daily recalculation of the dose distribution and the introduction of new adaptive radiotherapy treatment strategies. Methods: A Gammex Electron Density CT Phantom was imaged with the MVCB CT system. The pixel value for each of the sixteen inserts, which ranged from 0.292 to 1.707 relative electron density to the background solid water, was determined by taking the mean value from within a region of interest centred on the insert, over 5 slices within the centre of the phantom. These results were averaged and plotted against the relative electron densities of each insert with a linear least squares fit was preformed. This procedure was performed for images acquired with 5, 8, 15 and 60 monitor units. Results: The linear relationship between MVCT pixel value and ED was demonstrated for all monitor unit settings and over a range of electron densities. The number of monitor units utilised was found to have no significant impact on this relationship. Discussion: It was found that the number of MU utilised does not significantly alter the pixel value obtained for different ED materials. However, to ensure the most accurate and reproducible MV to ED calibration, one MU setting should be chosen and used routinely. To ensure accuracy for the clinical situation this MU setting should correspond to that which is used clinically. If more than one MU setting is used clinically then an average of the CT values acquired with different numbers of MU could be utilized without loss in accuracy. Conclusions: No significant differences have been shown between the pixel value to ED conversion for the Siemens MV CT cone beam unit with change in monitor units. Thus as single conversion curve could be utilised for MV CT treatment planning. To fully utilise MV CT imaging for radiotherapy treatment planning further work will be undertaken to ensure all corrections have been made and dose calculations verified. These dose calculations may be either for treatment planning purposes or for reconstructing the delivered dose distribution from transit dosimetry measurements made using electronic portal imaging devices. This will potentially allow the cumulative dose distribution to be determined through the patient’s multi-fraction treatment and adaptive treatment strategies developed to optimize the tumour response.

Prostate radiotherapy treatment plan quality : effects of hip prostheses

Introduction This study aimed to examine the geometric and dosimetric results when radiotherapy treatment plans are designed for prostate cancer patients with hip prostheses. Methods Ten EBRT treatment plans for localised prostate cancer, in the presence of hip prostheses, were analysed and compared with a reference set of 196 treatment plans for localised prostate cancer in patients without prostheses. Crowe et al.’s TADA code [1] was used to extract treatment plan parameters and evaluate doses to target volumes and critical structures against recommended goals [2] and constraints [3, 4]. Results The need to avoid transmitting the radiation beam through the hip prostheses limited the range of gantry angles available for use in both the rotational (VMAT) and the non-rotational (3DCRT and IMRT) radiotherapy treatments. This geometric limitation (exemplified in the VMAT data shown in Fig. 1) reduced the overall quality of the treatment plans for patients with prostheses compared to the reference plans. All plans with prostheses failed the PTV dose homogeneity requirement [2], whereas only 4 % of the plans without prostheses failed this test. Several treatment plans for patients with hip prostheses also failed the QUANTEC requirements that less than 50 % of the rectum receive 50 Gy and less than 35 % of the rectum receive 60 Gy to keep the grade 3 toxicity rate below 10 % [3], or the Hansen and Roach requirement that less than 25 % of the bladder receive 75 Gy [4]. Discussion and conclusions The results of this study exemplify the difficulty of designing prostate radiotherapy treatment plans, where beams provide adequate doses to targeted tissues while avoiding nearby organs at risk, when the presence of hip prostheses limits the available treatment geometries. This work provides qualitative evidence of the compromised dose distributions that can result, in such cases.

Gel dosimetry—it’s been the next big thing in radiotherapy for the past 25 years

Gel dosimetry and plastic chemical dosimeters such as PresageTM are capable of very accurately mapping dose distributions in three dimensions. Combined with their near tissue equivalence one would expect that after several decades of development they would be the dosimeter of choice for dosimetry, however they have not achieve widespread clinical use. This presentation will include a brief description and history of developments in gels and 3D plastics for dosimetry, the limitations and advantages, and their role in the future.

Calibration of treatment time intensity modulated EPID images for anatomical imaging

Intensity Modulated Radiotherapy (IMRT) is a well established technique for delivering highly conformal radiation dose distributions. The complexity of the delivery techniques and high dose gradients around the target volume make verification of the patient treatment crucial to the success of the treatment. Conventional treatment protocols involve imaging the patient prior to treatment, comparing the patient set-up to the planned set-up and then making any necessary shifts in the patient position to ensure target volume coverage. This paper presents a method for calibrating electronic portal imaging device (EPID) images acquired during IMRT delivery so that they can be used for verifying the patient set-up.

Interplay of MLC, gantry and respiratory motion during DCAT delivery

This study investigated the possible interplay effects arising from the treatment of moving targets using the dynamic conformal arc therapy (DCAT) technique. Dose from a modulated test beam was measured, with and without phantom motion and with and without a 30o arc rotation, using a diode array placed on a sinusoidally moving platform. Measurements were repeated at five different collimator angles (0, 22.5, 45, 67.5 and 90o), at two different dose rates (300 and 600 MU/min). Results showed that the effect of respiratory motion on the measured dose distribution increased slightly when the beams were delivered as arcs, rather than with a static gantry angle, and that this effect increased substantially as the collimator angle was increased from 0o (MLC motion perpendicular to respiratory motion) to 90o (MLC motion parallel to respiratory motion). The dose oscillations arising from interplay between phantom and MLC motion were found to increase in magnitude when the dose rate was increased. These results led to the development of simple recommendations for minimizing the negative effects of motion interplay on DCAT dose distributions

The level and distribution of chromosomal aberration of tomato seeds at different penetration depths of carbon ions

The relationship between the penetration depth and the level and distribution of chromosomal aberration of the root tip cells were investigated by exposure of the superposed tomato seeds to 80 MeV/u carbon ions. The results showed that on the entrance of the beam the chromosomal aberration level was low. Damage such as breaks and gaps were dominant. At the Bragg peak, the chromosomal aberration level was high. The yields of dicentrics, rings and disintegrated small chromosomes increased but the yields of breaks and gaps decreased. These results are consistent with the distribution of the physical depth dose pro. le of carbon ions. It is effective to deposit the Bragg peak on the seeds to induce hereditary aberration in the mutation breeding with heavy ions.

Physical property measurement of therapeutic carbon ion beam in the shallow-seated tumor therapy terminal at HIRFL

For the first time the physical properties of therapeutic carbon-ion beam supplied by, the shallow-seated tumor therapy terminal at the Heavy Ion Research Facility in Lanzhou (HIRFL) are measured. For a 80.55MeV/u C-12 ion beam delivered to the therapy terminal, the homogeneity of irradiation fields is 73.48%, when the beam intensity varied in the range of 0.001-0.1nA (i.e. 1 X 10(6) - 1 X 10(8) particles per second). The stability of the beam intensity within a few minutes is estimated to be 80.87%. The depth-dose distribution of the beam at the isocenter of the therapy facility is measured, and the position of the high-dose Bragg peak is found to be located at the water-equivalent depth of 13.866mm. Based on the relationship between beam energy and Bragg peak position, the corresponding beam energy at the isocenter of the therapy terminal is evaluated to be 71.71MeV/u for the original 80.55MeV/u C-12 ion beam, which consisted basically with calculation. The readout of the previously-used air-free ionization chamber regarding absorbed dose is calibrated as well in this experiment. The results indicate that the performance of the therapy facility should be optimized further to meet the requirements of clinical trial.

Germ cell loss induced by C-12(6+) ion irradiation in young female mice

The ovaries of Kun-Ming strain mice (3 weeks) were irradiated with different doses of C-12(6+) ion in the Bragg peak or the plateau region. At 10th day after irradiation, ovarian and uterine weights were measured: normal and atretic (identified with the oocyte to be degenerating or absent) primordial, primary and preantral follicles were identified in the largest cross-section of each ovary. Percentage (%) of normal follicles of each developmental stage of oogenesis was calculated. The data showed that compared to controls, there was a dose-related decrease in percentage of normal follicles in each developmental stage. And the weights of ovary and uterus were also reduced with doses of irradiation. Moreover, these effects were much more significant in the Bragg peak region and the region close to the Bragg peak than in the beam's entrance (the plateau region). Radiosensitivity varied in different follicle maturation stages. Primordial follicles, which are thought to be extremely sensitive to ionizing irradiation, were reduced by 86.6%, while primary and preantral follicles reduced only by 72.5% and 61.8% respectively, by exposure with 6 Gy of C-12(6+) ion in the Bragg peak region and the region close to the Bragg peak. The data suggested that due to their optimal depth-dose distribution in the Bragg peak region, heavy ions are ones of the best particles for radiotherapy of tumors located next of vital organs or/and surrounded by normal tissues, especially radiosensitive tissues such as gonads.

Simulations to design an online motion compensation system for scanned particle beams

Respiration-induced target motion is a major problem in intensity-modulated radiation therapy. Beam segments are delivered serially to form the total dose distribution. In the presence of motion, the spatial relation between dose deposition from different segments will be lost. Usually, this results in over-and underdosage. Besides such interplay effects between target motion and dynamic beam delivery as known from photon therapy, changes in internal density have an impact on delivered dose for intensity-modulated charged particle therapy. In this study, we have analysed interplay effects between raster scanned carbon ion beams and target motion. Furthermore, the potential of an online motion strategy was assessed in several simulations. An extended version of the clinical treatment planning software was used to calculate dose distributions to moving targets with and without motion compensation. For motion compensation, each individual ion pencil beam tracked the planned target position in the lateral aswell as longitudinal direction. Target translations and rotations, including changes in internal density, were simulated. Target motion simulating breathing resulted in severe degradation of delivered dose distributions. For example, for motion amplitudes of +/- 15 mm, only 47% of the target volume received 80% of the planned dose. Unpredictability of resulting dose distributions was demonstrated by varying motion parameters. On the other hand, motion compensation allowed for dose distributions for moving targets comparable to those for static targets. Even limited compensation precision (standard deviation similar to 2 mm), introduced to simulate possible limitations of real-time target tracking, resulted in less than 3% loss in dose homogeneity.

Inactive and mutagenic effects induced by carbon beams of different LET values in a red yeast strain

To evaluate biological action of microorganism exposed to charged particles during the long distance space exploration. Induction of inactivation and mutation in a red yeast strain Rhodotorula glutinis AY 91015 by carbon beams of different LET values (14.9-120 0 keV mu m(-1)) was investigated It was found that survival curves were exponential, and mutation curves were linear for all LET values The dependence of inactivation cross section on LET approached saturation near 120 0 keV mu m(-1) The imitation cross section saturated when LET was higher than 582 keV mu m(-1) Meanwhile, the highest RBEI for inactivation located at 120 0 key mu m(-1) and the highest RBEm for mutation was at 58.2 key mu m(-1) The experiments imply that the most efficient mutagenic part of the depth dose profile of carbon ion is at the plateau region with intermediate LET value in which energy deposited is high enough to Induce mutagenic lesions but too low to induce over kill effect in the yeast cells (C) 2010 Elsevier B V All rights reserved

Optical-CT 3D Dosimetry Using Fresnel Lenses with Minimal Refractive-Index Matching Fluid.

Telecentric optical computed tomography (optical-CT) is a state-of-the-art method for visualizing and quantifying 3-dimensional dose distributions in radiochromic dosimeters. In this work a prototype telecentric system (DFOS-Duke Fresnel Optical-CT Scanner) is evaluated which incorporates two substantial design changes: the use of Fresnel lenses (reducing lens costs from $10-30K t0 $1-3K) and the use of a 'solid tank' (which reduces noise, and the volume of refractively matched fluid from 1 ltr to 10 cc). The efficacy of DFOS was evaluated by direct comparison against commissioned scanners in our lab. Measured dose distributions from all systems were compared against the predicted dose distributions from a commissioned treatment planning system (TPS). Three treatment plans were investigated including a simple four-field box treatment, a multiple small field delivery, and a complex IMRT treatment. Dosimeters were imaged within 2 h post irradiation, using consistent scanning techniques (360 projections acquired at 1 degree intervals, reconstruction at 2mm). DFOS efficacy was evaluated through inspection of dose line-profiles, and 2D and 3D dose and gamma maps. DFOS/TPS gamma pass rates with 3%/3mm dose difference/distance-to-agreement criteria ranged from 89.3% to 92.2%, compared to from 95.6% to 99.0% obtained with the commissioned system. The 3D gamma pass rate between the commissioned system and DFOS was 98.2%. The typical noise rates in DFOS reconstructions were up to 3%, compared to under 2% for the commissioned system. In conclusion, while the introduction of a solid tank proved advantageous with regards to cost and convenience, further work is required to improve the image quality and dose reconstruction accuracy of the new DFOS optical-CT system.

Intensity-modulated radiotherapy reduces lung irradiation in patients with carcinoma of the oesophagus.

The 2-year survival rate after conventional radiotherapy for carcinoma of the oesophagus is around 10–20% [8]. Concomitant chemoradiation schedules have produced survival figures of 25–30% at 5 years, and this is now considered standard treatment [1]. Conformal radiotherapy techniques offer the potential to deliver higher doses of radiation to oesophageal tumours [5], and this may improve local tumour control. However, concerns regarding late normal tissue damage to the lung parenchyma and spinal cord remain a concern. Intensitymodulated radiotherapy (IMRT) allows complex dose distributions to be produced, and can reduce the dose to radiosensitive organs close to the tumour [2]. The present study was designed to investigate the impact of beam intensity modulation on treatment planning for carcinoma of the oesophagus, by comparing a standard three-dimensional conformal radiotherapy (3DCRT) technique to an IMRT technique using the same number and orientation of treatment fields.

Optimisation of radiotherapy for carcinoma of the parotid gland: a comparison of conventional, three-dimensional conformal, and intensity-modulated techniques

Background and purpose: To compare external beam radiotherapy techniques for parotid gland tumours using conventional radiotherapy (RT), three-dimensional conformal radiotherapy (3DCRT), and intensity-modulated radiotherapy (IMRT). To optimise the IMRT techniques, and to produce an IMRT class solution.Materials and methods: The planning target volume (PTV), contra-lateral parotid gland, oral cavity, brain-stem, brain and cochlea were outlined on CT planning scans of six patients with parotid gland tumours. Optimised conventional RT and 3DCRT plans were created and compared with inverse-planned IMRT dose distributions using dose-volume histograms. The aim was to reduce the radiation dose to organs at risk and improve the PTV dose distribution. A beam-direction optimisation algorithm was used to improve the dose distribution of the IMRT plans, and a class solution for parotid gland IMRT was investigated.Results: 3DCRT plans produced an equivalent PTV irradiation and reduced the dose to the cochlea, oral cavity, brain, and other normal tissues compared with conventional RT. IMRT further reduced the radiation dose to the cochlea and oral cavity compared with 3DCRT. For nine- and seven-field IMRT techniques, there was an increase in low-dose radiation to non-target tissue and the contra-lateral parotid gland. IMRT plans produced using three to five optimised intensity-modulated beam directions maintained the advantages of the more complex IMRT plans, and reduced the contra-lateral parotid gland dose to acceptable levels. Three- and four-field non-coplanar beam arrangements increased the volume of brain irradiated, and increased PTV dose inhomogeneity. A four-field class solution consisting of paired ipsilateral coplanar anterior and posterior oblique beams (15, 45, 145 and 170o from the anterior plane) was developed which maintained the benefits without the complexity of individual patient optimisation.Conclusions: For patients with parotid gland tumours, reduction in the radiation dose to critical normal tissues was demonstrated with 3DCRT compared with conventional RT. IMRT produced a further reduction in the dose to the cochlea and oral cavity. With nine and seven fields, the dose to the contra-lateral parotid gland was increased, but this was avoided by optimisation of the beam directions. The benefits of IMRT were maintained with three or four fields when the beam angles were optimised, but were also achieved using a four-field class solution. Clinical trials are required to confirm the clinical benefits of these improved dose distributions.

Out-of-field cell survival following exposure to intensity-modulated radiation fields

PURPOSE:
To determine the in-field and out-of-field cell survival of cells irradiated with either primary field or scattered radiation in the presence and absence of intercellular communication.
METHODS AND MATERIALS:
Cell survival was determined by clonogenic assay in human prostate cancer (DU145) and primary fibroblast (AGO1552) cells following exposure to different field configurations delivered using a 6-MV photon beam produced with a Varian linear accelerator.
RESULTS:
Nonuniform dose distributions were delivered using a multileaf collimator (MLC) in which half of the cell population was shielded. Clonogenic survival in the shielded region was significantly lower than that predicted from the linear quadratic model. In both cell lines, the out-of-field responses appeared to saturate at 40%-50% survival at a scattered dose of 0.70 Gy in DU-145 cells and 0.24 Gy in AGO1522 cells. There was an approximately eightfold difference in the initial slopes of the out-of-field response compared with the a-component of the uniform field response. In contrast, cells in the exposed part of the field showed increased survival. These observations were abrogated by direct physical inhibition of cellular communication and by the addition of the inducible nitric oxide synthase inhibitor aminoguanidine known to inhibit intercellular bystander effects. Additional studies showed the proportion of cells irradiated and dose delivered to the shielded and exposed regions of the field to impact on response.
CONCLUSIONS:
These data demonstrate out-of-field effects as important determinants of cell survival following exposure to modulated irradiation fields with cellular communication between differentially irradiated cell populations playing an important role. Validation of these observations in additional cell models may facilitate the refinement of existing radiobiological models and the observations considered important determinants of cell survival.