Convolution/superposition using the Monte Carlo method

A phase II study of weekly carboplatin and concurrent radiotherapy in older adults with locally advanced non-small cell lung cancer (LOGIK1902)

BACKGROUND

Concurrent chemoradiotherapy is the standard therapy for locally advanced non-small cell lung cancer (NSCLC). However, there is little evidence supporting its use in older adults. Low-dose daily carboplatin combined with thoracic radiotherapy is considered a standard regimen for this population. To establish a simple and feasible carboplatin administration method, we conducted a study of weekly carboplatin and concurrent radiotherapy for older adults with locally advanced NSCLC.

METHODS

This prospective, single-arm, multicenter, phase II clinical trial included patients aged ≥75 years with unresectable stage III NSCLC and Eastern Cooperative Oncology Group performance status 0-1. Patients received chemoradiotherapy (60 Gy/30 fractions plus concurrent weekly carboplatin at an area under curve of 2 mg mL-1 min-1). The primary endpoint was the overall response rate (ORR). Key secondary endpoints included progression-free survival (PFS), overall survival (OS), and safety.

RESULTS

From July 2020 to June 2022, 37 patients were enrolled from 15 institutions, and 36 patients were evaluable for efficacy and safety. The ORR was 63.9% (95% confidence interval [CI] = 47.6-77.5). Median PFS was 14.6 months (95% CI = 9.1-18.1). Median OS was 25.5 months (95% CI = 17.4-not reached). Grade 4 leucopenia, neutropenia, and thrombocytopenia were observed in one patient (2.8%) each.

CONCLUSION

Weekly carboplatin and concurrent radiation therapy was safe in older adults with locally advanced NSCLC, and promising activity was observed.

EGSnrc application for IMRT planning

The aim of this study was to describe a detailed instruction of intensity modulated radiotherapy (IMRT) planning simulation using BEAMnrc-DOSXYZnrc code system (EGSnrc package) and present a new graphical user interface based on MATLAB code (The MathWorks) to combine more than one. 3ddose file which were obtained from the IMRT plan. This study was performed in four phases: the commissioning of Varian Clinac iX6 MV, the simulation of IMRT planning in EGSnrc, the creation of in-house VDOSE GUI, and the analysis of the isodose contour and dose volume histogram (DVH) curve from several beam angles. The plan paramaters in sequence and control point files were extracted from the planning data in Tan Tock Seng Hospital Singapore (multileaf collimator (MLC) leaf positions - bank A and bank B, gantry angles, coordinate of isocenters, and MU indexes). VDOSE GUI which was created in this study can display the distribution dose curve in each slice and beam angle. Dose distributions from various MLC settings and beam angles yield different dose distributions even though they used the same number of simulated particles. This was due to the differences in the MLC leaf openings in every field. The value of the relative dose error between the two dose ditributions for "body" was 51.23 %. The Monte Carlo (MC) data was normalized with the maximum dose but the analytical anisotropic algorithm (AAA) data was normalized by the dose in the isocenter. In this study, we have presented a Monte Carlo simulation framework for IMRT dose calculation using DOSXYZnrc source 21. Further studies are needed in conducting IMRT simulations using EGSnrc to minimize the different dose error and dose volume histogram deviation.

A Simple Method for 2-D In Vivo Dosimetry by Portal Imaging

PURPOSE

To improve patient safety and treatment quality, verification of dose delivery in radiotherapy is desirable. We present a simple, easy-to-implement, open-source method for in vivo planar dosimetry of conformal radiotherapy by electronic portal imaging device (EPID).

METHODS

Correlation ratios, which relate dose in the mid-depth of slab phantoms to transit EPID signal, were determined for multiple phantom thicknesses and field sizes. Off-axis dose is corrected for by means of model-based convolution. We tested efficacy of dose reconstruction through measurements with off-reference values of attenuator thickness, field size, and monitor units. We quantified the dose calculation error in the presence of thickness changes to simulate anatomical or setup variations. An example of dose calculation on patient data is provided.

RESULTS

With varying phantom thickness, field size, and monitor units, dose reconstruction was almost always within 3% of planned dose. In the presence of thickness changes from planning CT, the dose discrepancy is exaggerated by up to approximately 1.5% for 1 cm changes upstream of the isocenter plane and 4% for 1 cm changes downstream.

CONCLUSION

Our novel electronic portal imaging device in vivo dosimetry allows clinically accurate 2-dimensional reconstruction of dose inside a phantom/patient at isocenter depth. Due to its simplicity, commissioning can be performed in a few hours per energy and may be modified to the user's needs. It may provide useful dose delivery information to detect harmful errors, guide adaptive radiotherapy, and assure quality of treatment.

New Developments in Intensity Modulated Radiation Therapy

As intensity modulated radiation therapy (IMRT) becomes routine clinical practice, its advantages and limitations are better understood. With these new understandings, some new developments have emerged in an effort to alleviate the limitations of the current IMRT practice. This article describes a few of these efforts made at the University of Maryland, including: i) improving IMRT efficiency with direct aperture optimization; ii) broadening the scope of optimization to include the mode of delivery and beam angles; and iii) new planning methods for intensity modulated arc therapy (IMAT).

Optimizing the MLC model parameters for IMRT in the RayStation treatment planning system

Unlike other commercial treatment planning systems (TPS) which model the rounded leaf end differently (such as the MLC dosimetric leaf gap (DLG) or rounded leaf‐tip radius), the RayStation TPS (RaySearch Laboratories, Stockholm, Sweden) models transmission through the rounded leaf end of the MLC with a step function, in which the radiation transmission through the leaf end is the square root of the average MLC transmission factor. We report on the optimization of MLC model parameters for the RayStation planning system. This (TPS) models the rounded leaf end of the MLC with the following parameters: leaf‐tip offset, leaf‐tip width, average transmission factor, and tongue and groove. We optimized the MLC model parameters for IMRT in the RayStation v. 4.0 planning system and for a Varian C‐series linac with a 120‐leaf Millennium MLC, and validated the model using measured data. The leaf‐tip offset is the geometric offset due to the rounded leaf‐end design and resulting divergence of the light/radiation field. The offset value is a function of the leaf‐tip position, and tabulated data are available from the vendor. The leaf‐tip width was iteratively evaluated by comparing computed and measured transverse dose profiles of MLC defined fields at dmax in water. In‐water profile comparisons were also used to verify the MLC leaf position (leaf‐tip offset). The average transmission factor and leaf tongue‐and‐groove width were derived iteratively by maximizing the agreement between measurements and RayStation TPS calculations for five clinical IMRT QA plans. Plan verifications were performed by comparing MapCHECK2 measurements and Monte Carlo calculations. The MLC model was validated using five test IMRT cases from the AAPM Task Group 119 report. Absolute gamma analyses (3 mm/3% and 2 mm/2%) were applied. In addition, computed output factors for MLC‐defined small fields (2×2,3×3,4×4,6×6 cm2) of both 6 MV and 18 MV photons were compared to those independently measured by the Imaging and Radiation Oncology Core (IROC), Houston, TX. 6 MV and 18 MV models were both determined to have the same MLC parameters: leaf‐tip offset=0.3 cm,2.5% transmission, and leaf tongue‐and‐groove width=0.05 cm. IMRT QA analysis for five test cases in TG‐119 resulted in a 100% passing rate with 3 mm/3% gamma analysis for 6 MV, and >97.5% for 18 MV. The passing rate was >94.6% for 6 MV and >90.9% for 18 MV when the 2 mm/2% gamma analysis criteria was applied. These results compared favorably with those published in AAPM Task Group 119. The reported MLC model parameters serve as a reference for other users.

PACS number(s): 87.55.D, 87.56.nk

Beam controlled arc therapy--a delivery concept for stationary targets

Volumetric modulated arc therapy (VMAT) presupposes that it is beneficial to deliver radiation from all beam angles as the gantry rotates, requiring the multi-leaf collimator to maintain continuity in shape from one angle to another. In turn, radiation from undesirable beam angles could compromise the dose distribution. In this work, we challenge the notion that the radiation beam must be held on as the gantry rotates around the patient. We propose a new approach for delivering intensity-modulated arc therapy, beam-controlled arc therapy (BCAT), during which the radiation beam is controlled on or off and the dose rate is modulated while the gantry rotates around the patient. We employ linear-programming-based dose optimization to each aperture weight, resulting in some zero weight apertures. During delivery, the radiation beam is held off at control points with zero weights as the MLC shape transits to the next non-zero weight shape. This was tested on ten head and neck cases. Plan quality and delivery efficiency were compared with VMAT. Improvements of up to 17% (p-value 0.001) and 57% (p-value 0.018) in organ-at-risk sparing and target dose uniformity, respectively, were achieved. Compared to the fixed number of apertures used in single-arc and double-arc VMAT, the BCAT used 109 and 175 apertures on average, respectively. The difference in total MUs for VMAT and BCAT plans was less than 4%. Plan quality improvement was confirmed after delivery with γ analysis resulting in over 99% agreement, or 4 in 1099 points that failed.

On the quantification of the dosimetric accuracy of collapsed cone convolution superposition (CCCS) algorithm for small lung volumes using IMRT

Specialized techniques that make use of small field dosimetry are common practice in today's clinics. These new techniques represent a big challenge to the treatment planning systems due to the lack of lateral electronic equilibrium. Because of this, the necessity of planning systems to overcome such difficulties and provide an accurate representation of the true value is of significant importance. Pinnacle3 is one such planning system. During the IMRT optimization process, Pinnacle3 treatment planning system allows the user to specify a minimum segment size which results in multiple beams composed of several subsets of different widths. In this study, the accuracy of the engine dose calculation, collapsed cone convolution superposition algorithm (CCCS) used by Pinnacle3, was quantified by Monte Carlo simulations, ionization chamber, and Kodak extended dose range film (EDR2) measurements for 11 SBRT lung patients. Lesions were < 3.0 cm in maximal diameter and <27.0cm3 in volume. The Monte Carlo EGSnrc\BEAMnrc and EGS4\MCSIM were used in the comparison. The minimum segment size allowable during optimization had a direct impact on the number of monitor units calculated for each beam. Plans with the smallest minimum segment size (0.1 cm2 to 2.0 cm2) had the largest number of MUs. Although PTV coverage remained unaffected, the segment size did have an effect on the dose to the organs at risk. Pinnacle3-calculated PTV mean doses were in agreement with Monte Carlo-calculated mean doses to within 5.6% for all plans. On average, the mean dose difference between Monte Carlo and Pinnacle3 for all 88 plans was 1.38%. The largest discrepancy in maximum dose was 5.8%, and was noted for one of the plans using a minimum segment size of 1.0 cm2. For minimum dose to the PTV, a maximum discrepancy between Monte Carlo and Pinnacle3 was noted of 12.5% for a plan using a 6.0 cm2 minimum segment size. Agreement between point dose measurements and Pinnacle3-calculated doses were on average within 0.7% in both phantoms. The profiles show a good agreement between Pinnacle3, Monte Carlo, and EDR2 film. The gamma index and the isodose lines support the result.

GENERALIZED GEOMETRIC APPROACHES FOR LEAF SEQUENCING PROBLEMS IN RADIATION THERAPY

The 3-D static leaf sequencing (SLS) problem arises in radiation therapy for cancer treatments, aiming to deliver a prescribed radiation dose to a target tumor accurately and efficiently. The treatment time and machine delivery error are two crucial factors to the solution (i.e., a treatment plan) for the SLS problem. In this paper, we prove that the 3-D SLS problem is NP-hard, and present the first ever algorithm for the 3-D SLS problem that can determine a tradeoff between the treatment time and machine delivery error (also called the "tongue-and-groove" error in medical literature). Our new 3-D SLS algorithm with error control gives the users (e.g., physicians) the option of specifying a machine delivery error bound, and subject to the given error bound, the algorithm computes a treatment plan with the minimum treatment time. We formulate the SLS problem with error control as computing a k-weight shortest path in a directed graph and build the graph by computing g-matchings and minimum cost flows. Further, we extend our 3-D SLS algorithm to all the popular radiotherapy machine models with different constraints. In our extensions, we model the SLS problems for some of the radiotherapy systems as computing a minimum g-path cover of a directed acyclic graph. We implemented our new 3-D SLS algorithm suite and conducted an extensive comparison study with commercial planning systems and well-known algorithms in medical literature. Some of our experimental results based on real medical data are presented.

Correlation of dosimetric parameters obtained with the analytical anisotropic algorithm and toxicity of chest chemoradiation in lung carcinoma

The purpose of this study was to analyze and revisit toxicity related to chest chemoradiotherapy and to correlate these side effects with dosimetric parameters obtained using analytical anisotropic algorithm (AAA) in locally unresectable advanced lung cancer. We retrospectively analyzed data from 47 lung cancer patients between 2005 and 2008. All received conformal 3D radiotherapy using high-energy linear accelerator plus concomitant chemotherapy. All treatment planning data were transferred into Eclipse 8.05 (Varian Medical Systems, Palo Alto, CA) and dosimetric calculations were performed using AAA. Thirty-three patients (70.2%) developed acute pneumopathy after radiotherapy (grades 1 and 2). One patient (2.1%) presented with grade 3 pneumopathy. Thirty-one (66%) presented with grades 1-2 lung fibrosis, and 1 patient presented with grade 3 lung fibrosis. Thirty-four patients (72.3%) developed grade 1-2 acute oesophagic toxicity. Four patients (8.5%) presented with grades 3 and 4 dysphagia, necessitating prolonged parenteral nutrition. Median prescribed dose was 64 Gy (range 50-74) with conventional fractionation (2 Gy per fraction). Dose-volume constraints were respected with a median V20 of 23.5% (maximum 34%) and a median V30 of 17% (maximum 25%). The median dose delivered to healthy contralateral lung was 13.1 Gy (maximum 18.1 Gy). At univariate analysis, larger planning target volume and V20 were significantly associated with the probability of grade ≥2 radiation-induced pneumopathy (p = 0.022 and p = 0.017, respectively). No relation between oesophagic toxicity and clinical/dosimetric parameters could be established. Using AAA, the present results confirm the predictive value of the V20 for lung toxicity as already demonstrated with the conventional pencil beam convolution approach.

Intensity-modulated arc therapy: principles, technologies and clinical implementation

Intensity-modulated arc therapy (IMAT) was proposed by Yu (1995 Phys. Med. Biol. 40 1435-49) as an alternative to tomotherapy. Over more than a decade, much progress has been made. The advantages and limitations of the IMAT technique have also been better understood. In recent years, single-arc forms of IMAT have emerged and become commercially adopted. The leading example is the volumetric-modulated arc therapy (VMAT), a single-arc form of IMAT that delivers apertures of varying weights with a single-arc rotation that uses dose-rate variation of the treatment machine. With commercial implementation of VMAT, wide clinical adoption has quickly taken root. However, there remains a lack of general understanding for the planning of such arc treatments, as well as what delivery limitations and compromises are made. Commercial promotion and competition add further confusion for the end users. It is therefore necessary to provide a summary of this technology and some guidelines on its clinical implementation. The purpose of this review is to provide a summary of the works from the radiotherapy community that led to wide clinical adoption, and point out the issues that still remain, providing some perspective on its further developments. Because there has been vast experience in IMRT using multiple intensity-modulated fields, comparisons between IMAT and IMRT are also made in the review within the areas of planning, delivery and quality assurance.

Comparing radiation treatments using intensity-modulated beams, multiple arcs, and single arcs

PURPOSE

A dosimetric comparison of multiple static-field intensity-modulated radiation therapy (IMRT), multiarc intensity-modulated arc therapy (IMAT), and single-arc arc-modulated radiation therapy (AMRT) was performed to evaluate their clinical advantages and shortcomings.

METHODS AND MATERIALS

Twelve cases were selected for this study, including three head-and-neck, three brain, three lung, and three prostate cases. An IMRT, IMAT, and AMRT plan was generated for each of the cases, with clinically relevant planning constraints. For a fair comparison, the same parameters were used for the IMRT, IMAT, and AMRT planning for each patient.

RESULTS

Multiarc IMAT provided the best plan quality, while single-arc AMRT achieved dose distributions comparable to those of IMRT, especially in the complicated head-and-neck and brain cases. Both AMRT and IMAT showed effective normal tissue sparing without compromising target coverage and delivered a lower total dose to the surrounding normal tissues in some cases.

CONCLUSIONS

IMAT provides the most uniform and conformal dose distributions, especially for the cases with large and complex targets, but with a delivery time similar to that of IMRT; whereas AMRT achieves results comparable to IMRT with significantly faster treatment delivery.

A two-stage sequential linear programming approach to IMRT dose optimization

The conventional IMRT planning process involves two stages in which the first stage consists of fast but approximate idealized pencil beam dose calculations and dose optimization and the second stage consists of discretization of the intensity maps followed by intensity map segmentation and a more accurate final dose calculation corresponding to physical beam apertures. Consequently, there can be differences between the presumed dose distribution corresponding to pencil beam calculations and optimization and a more accurately computed dose distribution corresponding to beam segments that takes into account collimator-specific effects. IMRT optimization is computationally expensive and has therefore led to the use of heuristic (e.g., simulated annealing and genetic algorithms) approaches that do not encompass a global view of the solution space. We modify the traditional two-stage IMRT optimization process by augmenting the second stage via an accurate Monte Carlo-based kernel-superposition dose calculations corresponding to beam apertures combined with an exact mathematical programming-based sequential optimization approach that uses linear programming (SLP). Our approach was tested on three challenging clinical test cases with multileaf collimator constraints corresponding to two vendors. We compared our approach to the conventional IMRT planning approach, a direct-aperture approach and a segment weight optimization approach. Our results in all three cases indicate that the SLP approach outperformed the other approaches, achieving superior critical structure sparing. Convergence of our approach is also demonstrated. Finally, our approach has also been integrated with a commercial treatment planning system and may be utilized clinically.

eIMRT: a web platform for the verification and optimization of radiation treatment plans

The eIMRT platform is a remote distributed computing tool that provides users with Internet access to three different services: Monte Carlo optimization of treatment plans, CRT & IMRT treatment optimization, and a database of relevant radiation treatments/clinical cases. These services are accessible through a user-friendly and platform independent web page. Its flexible and scalable design focuses on providing the final users with services rather than a collection of software pieces. All input and output data (CT, contours, treatment plans and dose distributions) are handled using the DICOM format. The design, implementation, and support of the verification and optimization algorithms are hidden to the user. This allows a unified, robust handling of the software and hardware that enables these computation-intensive services. The eIMRT platform is currently hosted by the Galician Supercomputing Center (CESGA) and may be accessible upon request (there is a demo version at http://eimrt.cesga.es:8080/eIMRT2/demo; request access in http://eimrt.cesga.es/signup.html). This paper describes all aspects of the eIMRT algorithms in depth, its user interface, and its services. Due to the flexible design of the platform, it has numerous applications including the intercenter comparison of treatment planning, the quality assurance of radiation treatments, the design and implementation of new approaches to certain types of treatments, and the sharing of information on radiation treatment techniques. In addition, the web platform and software tools developed for treatment verification and optimization have a modular design that allows the user to extend them with new algorithms. This software is not a commercial product. It is the result of the collaborative effort of different public research institutions and is planned to be distributed as an open source project. In this way, it will be available to any user; new releases will be generated with the new implemented codes or upgrades.

Effects of tumor motion in GRID therapy

Clinical and biological evidence suggest that the success of GRID therapy in debulking large tumors depends on the high peak-to-valley contrast in the dose distribution. In this study, we show that the peaks and valleys can be significantly blurred out by respiration-induced tumor motion, possibly affecting the clinical outcome. Using a kernel-based Monte Carlo dose engine that incorporates phantom motion, we calculate the dose distributions for a GRID with hexagonally arranged holes. The holes have a diameter of 1.3 cm and a minimum center-to-center separation of 2.1 cm (projected at the isocenter). The phantom moves either in the u parallel direction, which is parallel to a line joining any two nearest neighbors, or in the perpendicular u perpendicular direction. The displacement-time waveform is modeled with a cosn function, with n assigned 1 for symmetric motion, or 6 to simulate a large inhale-exhale asymmetry. Dose calculations are performed on a water phantom for a 6 MV x-ray beam. Near dmax, the static valley dose is 0.12D0, where D0 is the peak static dose. For motion in the u parallel direction, the peak and valley doses vary periodically with the amplitude of motion a and the transverse dose profiles are maximally flat near a=0.8 cm and a=1.9 cm. For the cos waveform, the minimum peak dose (Dpmin) is 0.67D0 and the maximum valley dose (Dvmax) is 0.60D0. Less dose blurring is seen with the cos6 waveform, with Dpmin=0.77D0 and Dvmax=0.45D0. For motion in the u perpendicular direction, the maximum flattening of dose profiles occurs at a=1.5 cm. GRIDs with smaller hole separations produce similar blurring at proportionally smaller amplitudes. The reported clinical response data from GRID therapy seem to indicate that mobile tumors, such as those in the thorax and abdomen, respond worse to GRID treatments than stationary tumors, such as those in the head and neck. To establish a stronger correlation between clinical response and tumor motion, and possibly improve the clinical response rates, it is recommended that prospective GRID therapy trials be conducted with motion compensation strategies, such as respiratory gating.

Arc-modulated radiation therapy (AMRT): a single-arc form of intensity-modulated arc therapy

Arc-modulated radiation therapy (AMRT) is a novel rotational intensity-modulated radiation therapy (IMRT) technique developed for a clinical linear accelerator that aims to deliver highly conformal radiation treatment using just one arc of gantry rotation. Compared to fixed-gantry IMRT and the multiple-arc intensity-modulated arc therapy (IMAT) techniques, AMRT promises the same treatment quality with a single-arc delivery. In this paper, we present a treatment planning scheme for AMRT, which addresses the challenges in inverse planning, leaf sequencing and dose calculation. The feasibility and performance of this AMRT treatment planning scheme have been verified with multiple clinical cases of various sites on Varian linear accelerators.

Automatic determination of primary electron beam parameters in Monte Carlo simulation

In order to obtain realistic and reliable Monte Carlo simulations of medical linac photon beams, an accurate determination of the parameters that define the primary electron beam that hits the target is a fundamental step. In this work we propose a new methodology to commission photon beams in Monte Carlo simulations that ensures the reproducibility of a wide range of clinically useful fields. For such purpose accelerated Monte Carlo simulations of 2 x 2, 10 x 10, and 20 x 20 cm2 fields at SSD = 100 cm are carried out for several combinations of the primary electron beam mean energy and radial FWHM. Then, by performing a simultaneous comparison with the correspondent measurements for these same fields, the best combination is selected. This methodology has been employed to determine the characteristics of the primary electron beams that best reproduce a Siemens PRIMUS and a Varian 2100 CD machine in the Monte Carlo simulations. Excellent agreements were obtained between simulations and measurements for a wide range of field sizes. Because precalculated profiles are stored in databases, the whole commissioning process can be fully automated, avoiding manual fine-tunings. These databases can also be used to characterize any accelerators of the same model from different sites.

Direct aperture optimization of breast IMRT and the dosimetric impact of respiration motion

We have studied the application of direct aperture optimization (DAO) as an inverse planning tool for breast IMRT. Additionally, we have analysed the impact of respiratory motion on the quality of the delivered dose distribution. From this analysis, we have developed guidelines for balancing the desire for a high-quality optimized plan with the need to create a plan that will not degrade significantly in the presence of respiratory motion. For a DAO optimized breast IMRT plan, the tangential fields incorporate a flash field to cover the range of respiratory motion. The inverse planning algorithm then optimizes the shapes and weights of additional segments that are delivered in combination with the open fields. IMRT plans were generated using DAO with the relative weights of the open segments varied from 0% to 95%. To assess the impact of breathing motion, the dose distribution for the optimized IMRT plan was recalculated with the isocentre sampled from a predefined distribution in a Monte Carlo convolution/superposition dose engine with the breast simulated as a rigid object. The motion amplitudes applied in this study ranged from 0.5 to 2.0 cm. For a range of weighting levels assigned to the open field, comparisons were made between the static plans and the plans recalculated with motion. For the static plans, we found that uniform dose distributions could be generated with relative weights for the open segments equal to and below 80% and unacceptable levels of underdosage were observed with the weights larger than 80%. When simulated breathing motion was incorporated into the dose calculation, we observed a loss in dose uniformity as the weight of the open field was decreased to below 65%. More quantitatively, for each 1% decrease in the weight, the per cent volume of the target covered by at least 95% of the prescribed dose decreased by approximately 0.10% and 0.16% for motion amplitudes equal to 1.5 cm and 2.0 cm, respectively. When taking into account the motion effects, the most uniform and conformal dose distributions were achieved when the open segment weights were in the range of 65-80%. Within this range, high-quality IMRT plans were produced for each case. The study demonstrates that DAO with tangential fields provides a robust and efficient technique for breast IMRT planning and delivery when the open segment weight is selected between 65% and 80%.

Direct aperture optimization for IMRT using Monte Carlo generated beamlets

This work introduces an EGSnrc-based Monte Carlo (MC) beamlet does distribution matrix into a direct aperture optimization (DAO) algorithm for IMRT inverse planning. The technique is referred to as Monte Carlo-direct aperture optimization (MC-DAO). The goal is to assess if the combination of accurate Monte Carlo tissue inhomogeneity modeling and DAO inverse planning will improve the dose accuracy and treatment efficiency for treatment planning. Several authors have shown that the presence of small fields and/or inhomogeneous materials in IMRT treatment fields can cause dose calculation errors for algorithms that are unable to accurately model electronic disequilibrium. This issue may also affect the IMRT optimization process because the dose calculation algorithm may not properly model difficult geometries such as targets close to low-density regions (lung, air etc.). A clinical linear accelerator head is simulated using BEAMnrc (NRC, Canada). A novel in-house algorithm subdivides the resulting phase space into 2.5 X 5.0 mm2 beamlets. Each beamlet is projected onto a patient-specific phantom. The beamlet dose contribution to each voxel in a structure-of-interest is calculated using DOSXYZnrc. The multileaf collimator (MLC) leaf positions are linked to the location of the beamlet does distributions. The MLC shapes are optimized using direct aperture optimization (DAO). A final Monte Carlo calculation with MLC modeling is used to compute the final dose distribution. Monte Carlo simulation can generate accurate beamlet dose distributions for traditionally difficult-to-calculate geometries, particularly for small fields crossing regions of tissue inhomogeneity. The introduction of DAO results in an additional improvement by increasing the treatment delivery efficiency. For the examples presented in this paper the reduction in the total number of monitor units to deliver is approximately 33% compared to fluence-based optimization methods.

A difference-matrix metaheuristic for intensity map segmentation in step-and-shoot IMRT delivery

At an intermediate stage of radiation treatment planning for IMRT, most commercial treatment planning systems for IMRT generate intensity maps that describe the grid of beamlet intensities for each beam angle. Intensity map segmentation of the matrix of individual beamlet intensities into a set of MLC apertures and corresponding intensities is then required in order to produce an actual radiation delivery plan for clinical use. Mathematically, this is a very difficult combinatorial optimization problem, especially when mechanical limitations of the MLC lead to many constraints on aperture shape, and setup times for apertures make the number of apertures an important factor in overall treatment time. We have developed, implemented and tested on clinical cases a metaheuristic (that is, a method that provides a framework to guide the repeated application of another heuristic) that efficiently generates very high-quality (low aperture number) segmentations. Our computational results demonstrate that the number of beam apertures and monitor units in the treatment plans resulting from our approach is significantly smaller than the corresponding values for treatment plans generated by the heuristics embedded in a widely use commercial system. We also contrast the excellent results of our fast and robust metaheuristic with results from an 'exact' method, branch-and-cut, which attempts to construct optimal solutions, but, within clinically acceptable time limits, generally fails to produce good solutions, especially for intensity maps with more than five intensity levels. Finally, we show that in no instance is there a clinically significant change of quality associated with our more efficient plans.

Feasibility of delivering grid therapy using a multileaf collimator

The feasibility of using a multileaf collimator (MLC) for grid therapy is demonstrated in this study. Grids with the projected field openings of 10 mm x 10 mm and 5 mm x 5 mm were created using multiple MLC-shaped fields. The deposited doses were measured with films at different depths in a solid water phantom and compared to those of Cerrobend grid collimators of similar hole sizes and hole separations. At the depth of maximum dose (dmax), the valley-to-peak dose ratios of the MLC grids were found to be about 11% and 19% for the respective 10 mm x 10 mm and 5 mm X 5 mm grid openings, and those of the corresponding grid blocks were about 15% and 20%. To quantify the dose contributed by transmission in the blocked areas due to the limited leaf thickness, Monte Carlo simulations (based on convolution/superposition method) were performed to calculate the doses in the solid water phantom using an ideal MLC with no leakage and perfect divergence in both the leaf end and side. About 7% reduction in the valley-to-peak dose ratio was found for both grid sizes at dmax. The results clearly showed that MLCs can be used to provide grid treatments with at least as good dosimetric properties as those of the Cerrobend grid blocks, though the former would in general require a longer delivery time.

Effect of beamlet step-size on IMRT plan quality

We have studied the degree to which beamlet step-size impacts the quality of intensity modulated radiation therapy (IMRT) treatment plans. Treatment planning for IMRT begins with the application of a grid that divides each beam's-eye-view of the target into a number of smaller beamlets (pencil beams) of radiation. The total dose is computed as a weighted sum of the dose delivered by the individual beamlets. The width of each beamlet is set to match the width of the corresponding leaf of the multileaf collimator (MLC). The length of each beamlet (beamlet step-size) is parallel to the direction of leaf travel. The beamlet step-size represents the minimum stepping distance of the leaves of the MLC and is typically predetermined by the treatment planning system. This selection imposes an artificial constraint because the leaves of the MLC and the jaws can both move continuously. Removing the constraint can potentially improve the IMRT plan quality. In this study, the optimized results were achieved using an aperture-based inverse planning technique called direct aperture optimization (DAO). We have tested the relationship between pencil beam step-size and plan quality using the American College of Radiology's IMRT test case. For this case, a series of IMRT treatment plans were produced using beamlet step-sizes of 1, 2, 5, and 10 mm. Continuous improvements were seen with each reduction in beamlet step size. The maximum dose to the planning target volume (PTV) was reduced from 134.7% to 121.5% and the mean dose to the organ at risk (OAR) was reduced from 38.5% to 28.2% as the beamlet step-size was reduced from 10 to 1 mm. The smaller pencil beam sizes also led to steeper dose gradients at the junction between the target and the critical structure with gradients of 6.0, 7.6, 8.7, and 9.1 dose%/mm achieved for beamlet step sizes of 10, 5, 2, and 1 mm, respectively.

Using a photon phase-space source for convolution/superposition dose calculations in radiation therapy

For a given linac design, the dosimetric characteristics of a photon beam are determined uniquely by the energy and radial distributions of the electron beam striking the x-ray target. However, in the usual commissioning of a beam from measured data, a large number of variables can be independently tuned, making it difficult to derive a unique and self-consistent beam model. For example, the measured dosimetric penumbra in water may be attributed in various proportions to the lateral secondary electron range, the focal spot size and the transmission through the tips of a non-divergent collimator; the head-scatter component in the tails of the transverse profiles may not be easy to resolve from phantom scatter and head leakage; and the head-scatter tails corresponding to a certain extra-focal source model may not agree self-consistently with in-air output factors measured on the central axis. To reduce the number of adjustable variables in beam modelling, we replace the focal and extra-focal sources with a single phase-space plane scored just above the highest adjustable collimator in a EGS/BEAM simulation of the linac. The phase-space plane is then used as photon source in a stochastic convolution/superposition dose engine. A photon sampled from the uncollimated phase-space plane is first propagated through an arbitrary collimator arrangement and then interacted in the simulation phantom. Energy deposition kernel rays are then randomly issued from the interaction points and dose is deposited along these rays. The electrons in the phase-space file are used to account for electron contamination. 6 MV and 18 MV photon beams from an Elekta SL linac are used as representative examples. Except for small corrections for monitor backscatter and collimator forward scatter for large field sizes (<0.5% with <20 x 20 cm2 field size), we found that the use of a single phase-space photon source provides accurate and self-consistent results for both relative and absolute dose calculations.

A stochastic convolution/superposition method with isocenter sampling to evaluate intrafraction motion effects in IMRT

Current methods to calculate dose distributions with organ motion can be broadly classified as "dose convolution" and "fluence convolution" methods. In the former, a static dose distribution is convolved with the probability distribution function (PDF) that characterizes the motion. However, artifacts are produced near the surface and around inhomogeneities because the method assumes shift invariance. Fluence convolution avoids these artifacts by convolving the PDF with the incident fluence instead of the patient dose. In this paper we present an alternative method that improves the accuracy, generality as well as the speed of dose calculation with organ motion. The algorithm starts by sampling an isocenter point from a parametrically defined space curve corresponding to the patient-specific motion trajectory. Then a photon is sampled in the linac head and propagated through the three-dimensional (3-D) collimator structure corresponding to a particular MLC segment chosen randomly from the planned IMRT leaf sequence. The photon is then made to interact at a point in the CT-based simulation phantom. Randomly sampled monoenergetic kernel rays issued from this point are then made to deposit energy in the voxels. Our method explicitly accounts for MLC-specific effects (spectral hardening, tongue-and-groove, head scatter) as well as changes in SSD with isocentric displacement, assuming that the body moves rigidly with the isocenter. Since the positions are randomly sampled from a continuum, there is no motion discretization, and the computation takes no more time than a static calculation. To validate our method, we obtained ten separate film measurements of an IMRT plan delivered on a phantom moving sinusoidally, with each fraction starting with a random phase. For 2 cm motion amplitude, we found that a ten-fraction average of the film measurements gave an agreement with the calculated infinite fraction average to within 2 mm in the isodose curves. The results also corroborate the existing notion that the interfraction dose variability due to the interplay between the MLC motion and breathing motion averages out over typical multifraction treatments. Simulation with motion waveforms more representative of real breathing indicate that the motion can produce penumbral spreading asymmetric about the static dose distributions. Such calculations can help a clinician decide to use, for example, a larger margin in the superior direction than in the inferior direction. In the paper we demonstrate that a 15 min run on a single CPU can readily illustrate the effect of a patient-specific breathing waveform, and can guide the physician in making informed decisions about margin expansion and dose escalation.