A model based on a specific phantom, called QuAArC, has been designed for the evaluation of planning and verification systems of complex radiotherapy treatments, such as volumetric modulated arc therapy (VMAT). This model uses the high accuracy provided by the Monte Carlo (MC) simulation of log files and allows the experimental feedback from the high spatial resolution of films hosted in QuAArC. This cylindrical phantom was specifically designed to host films rolled at different radial distances able to take into account the entrance fluence and the 3D dose distribution. Ionization chamber measurements are also included in the feedback process for absolute dose considerations. In this way, automated MC simulation of treatment log files is implemented to calculate the actual delivery geometries, while the monitor units are experimentally adjusted to reconstruct the dose-volume histogram (DVH) on the patient CT. Prostate and head and neck clinical cases, previously planned with Monaco and Pinnacle treatment planning systems and verified with two different commercial systems (Delta4 and COMPASS), were selected in order to test operational feasibility of the proposed model. The proper operation of the feedback procedure was proved through the achieved high agreement between reconstructed dose distributions and the film measurements (global gamma passing rates > 90% for the 2%/2 mm criteria). The necessary discretization level of the log file for dose calculation and the potential mismatching between calculated control points and detection grid in the verification process were discussed. Besides the effect of dose calculation accuracy of the analytic algorithm implemented in treatment planning systems for a dynamic technique, it was discussed the importance of the detection density level and its location in VMAT specific phantom to obtain a more reliable DVH in the patient CT. The proposed model also showed enough robustness and efficiency to be considered as a pre-treatment VMAT verification system.
PURPOSE:The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice.METHODS:The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called "biophysical" map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process.RESULTS:Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements.CONCLUSIONS:A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.
The manufacturer has introduced the new EBT2 film model so as to improve its predecessor, the EBT radiochromic film model. According to the manufacturer, some of its main advantages include a higher tolerance to light exposure and it can correct non-uniformity of the active layer thickness using a marker dye. However, the equivalence in uniformity between both models was questioned by some authors, and the asymmetrical configuration of layers of the EBT2 film model produces a new dependence on the film side being scanned (front and back orientation). In this study, the EBT2 radiochromic film model was compared with the EBT model and the new marker dye feature was assessed. We also compared this correction method with a pre-irradiated pixel value correction method. An Epson Expression 10000XL scanner in transmission mode was used to scan the films and the red channel response was analyzed. We confirmed the lower-measured signal dependence on the visible light exposure of the EBT2 film model. Differences in pixel values remained below 0.5% for a minimum of 15 days. In regard to the uniformity, similar results for EBT2 and EBT film models were obtained; in both cases inhomogeneity was found to be less than 1%, in relative pixel value from the mean. However, we found that the signal-to-noise ratio was reduced for low doses by 37% for old EBT2 batch and by 21% for new EBT2 batch compared to signal-to-noise ratio for EBT. The EBT2 film model's pixel value difference for the front and back orientation reached up to 1.0% in the red channel. Our results did not show a clear advantage between to use a pre-irradiated pixel value correction and to use the manufacturer's correction.
Monte Carlo simulation is an accurate way of assessing radiotherapy dose distribution in nonhomogeneous volumes, but it requires long processing times. A new distribution model simulates radiotherapy treatments and runs on a PC network, which reduces the processing time and makes for a powerful treatment-verification tool.
Total skin electron therapy (TSET) is a complex technique which requires non-standard measurements and dosimetric procedures. This paper investigates an essential first step towards TSET Monte Carlo (MC) verification. The non-standard 6 MeV 40 x 40 cm2 electron beam at a source to surface distance (SSD) of 100 cm as well as its horizontal projection behind a polymethylmethacrylate (PMMA) screen to SSD = 380 cm were evaluated. The EGS4 OMEGA-BEAM code package running on a Linux home made 47 PCs cluster was used for the MC simulations. Percentage depth-dose curves and profiles were calculated and measured experimentally for the 40 x 40 cm2 field at both SSD = 100 cm and patient surface SSD = 380 cm. The output factor (OF) between the reference 40 x 40 cm2 open field and its horizontal projection as TSET beam at SSD = 380 cm was also measured for comparison with MC results. The accuracy of the simulated beam was validated by the good agreement to within 2% between measured relative dose distributions, including the beam characteristic parameters (R50, R80, R100, Rp, E0) and the MC calculated results. The energy spectrum, fluence and angular distribution at different stages of the beam (at SSD = 100 cm, at SSD = 364.2 cm, behind the PMMA beam spoiler screen and at treatment surface SSD = 380 cm) were derived from MC simulations. Results showed a final decrease in mean energy of almost 56% from the exit window to the treatment surface. A broader angular distribution (FWHM of the angular distribution increased from 13 degrees at SSD = 100 cm to more than 30 degrees at the treatment surface) was fully attributable to the PMMA beam spoiler screen. OF calculations and measurements agreed to less than 1%. The effect of changing the electron energy cut-off from 0.7 MeV to 0.521 MeV and air density fluctuations in the bunker which could affect the MC results were shown to have a negligible impact on the beam fluence distributions. Results proved the applicability of using MC as a treatment verification tool for complex radiotherapy techniques.
BACKGROUND AND PURPOSE:In treatments where it is necessary to conform the field shape yielding a very small effective beam area, dosimetry and conventional treatment planning may be inaccurate. The Monte Carlo (MC) method can be an alternative to verify dose calculations. A conjunctival mucosa-associated lymphoid tissues lymphoma is presented, to show the importance of an independent assessment in critical situations.MATERIALS AND METHODS:In this work, the MC technique has been employed using the program BEAM (based on EGS4 code). Electron beam simulation has been performed and the results have been compared with those obtained with films. The patient dose distribution has been obtained by two methods: the full Monte Carlo (FMC) simulation and a conventional planning system (PLATO).RESULTS:Concerning dosimetry, some differences have been observed in the comparison of profiles obtained with film and those obtained with the MC method. Moreover, significant differences were found in the patient isodose distribution between both calculation methods.CONCLUSIONS:The results highlight that, in treatments where small beams are needed, conventional dosimetry and planning systems have some limitations. Therefore, an independent and more accurate assessment, such as MC, would be desirable.
Small field dosimetry is complicated due to the lack of electronic equilibrium and to the high steep dose gradients. This works compares PDD curves, profiles and output factors measured with conventional detectors (film, diode, TLD and ionisation chamber) and calculated with Monte Carlo. The 6 MV nominal energy from a Philips SL-18 linac has been simulated by using the OMEGA code. MC calculation reveals itself as a convenient method to validate OF and profiles in special conditions, such as small fields.
The wide range of possibilities available in Radiotherapy with conformai fields cannot be covered experimentally. For this reason, dosimetrical and planning procedures are based on approximate algorithms or systematic measurements. Dose distribution calculations based on Monte Carlo (MC) simulations can be used to check results. In this work, two examples of conformai field treatments are shown: A prostate carcinoma and an ocular lymphoma. The dose distributions obtained with a conventional Planning System and with MC have been compared. Some significant differences have been found.
The precision of Radiosurgery Treatment planning systems is limited by the approximations of their algorithms and by their dosimetrical input data. This fact is especially important in small fields. However, the Monte Carlo methods is an accurate alternative as it considers every aspect of particle transport. In this work an acoustic neurinoma is studied by comparing the dose distribution of both a planning system and Monte Carlo. Relative shifts have been measured and furthermore, Dose-Volume Histograms have been calculated for target and adjacent organs at risk.
CPU time is one of the most important inconveniences in the general use of Monte Carlo calculations applied to Radiotherapy. Although some variance reduction techniques have been developed, the time lessening achieved is not enough for the practical use of Full Monte Carlo (FMC). Moreover, the expected evolution of processors is not as fast as desirable. Our choice has then been the optimisation of CPU time by means of the distributed processing technique by which, the model developed divides the simulation of the treatment in several parts, having each a given physical sense. The EGS4 code implemented under OMEGA-BEAM has been used. This code has been compiled on a HP X-Class SPP 2000 machine with 24 processors PA-RISC 8000 type. The distribution process is carried out by means of customised script which run under graphical interface. Most treatments are solved in less than one hour, while complex multi-isocentric are treatments (10 arcs with three different isocentres on a 256x256 voxel geometry applied to 72 CT slices with 4096 density levels) are calculated between one and two hours.
A distribution model to simulate physical dosimetry measurements with Monte Carlo (MC) techniques has been developed. This approach is indicated to solve the simulations where there are continuous changes of measurement conditions (and hence of the input parameters) such as a TPR curve or the estimation of the resolution limit of an optical densitometer in the case of small field profiles. As a comparison, a high resolution scan for narrow beams with no iterative process is presented.The model has been installed on a network PCs without any resident software. The only requirement for these PCs has been a small and temporal Linux partition in the hard disks and to be connecting by the net with our server PC.
The choice of the most appropriate strategy in a Radiotherapy treatment is mainly based on the use of a planning system. With the introduction of new techniques (conformal and/or small fields, asymmetrical and non coplanar beams, true 3D calculations, IMRT) the trustworthiness of the algorithms is being questioned. An alternative verification procedure is every time more necessary to warranty a treatment delivery. The reliability of Monte Carlo is generally accepted. However, its clinical use has not been operative due to the high CPU times needed. During the last few years our objective has been focussed to reduce this time by means of new process distribution techniques. Tnis drop has made it feasible, not only the physical dosimetry under special conditions, but also a numerous variety of clinical cases: photon and electron conformal fields, Radiosurgery and IMRT. The carried out procedure is presented. Furthermore, experimental dosimetry data as well as conventional TPS calculations are compared with Monte Carlo simulations.