A Method for Estimating the Lung Clinical Target Volume DVH from IMRT with and without Respiratory Gating

Motion of lung tumors from respiration has been reported in the literature to be as large as of 1-2 cm. This motion requires an additional margin between the Clinical Target Volume (CTV) and the Planning Target Volume (PTV). While such a margin is necessary, it may not be sufficient to ensure proper delivery of Intensity Modulated Radiotherapy (IMRT) to the CTV during the simultaneous movement of the DMLC. Gated treatment has been proposed to improve normal tissues sparing as well as to ensure accurate dose coverage of the tumor volume. The following questions have not been addressed in the literature: a) what is the dose error to a target volume without gated IMRT treatment\ulcorner b) what is an acceptable gating window for such treatment. In this study, we address these questions by proposing a novel technique for calculating the 3D dose error that would result if a lung IMRT plan were delivered without gating. The method is also generalized for gated treatment with an arbitrary triggering window. IMRT plans for three patients with lung tumor were studied. The treatment plans were generated with HELIOS for delivery with 6 MV on a CL2100 Varian linear accelerator with a 26 pair MLC. A CTV to PTV margin of 1 cm was used. An IMRT planning system searches for an optimized fluence map ${\Phi}$ (x,y) for each port, which is then converted into a dynamic MLC file (DMLC). The DMLC file contains information about MLC subfield shapes and the fractional Monitor Units (MUs) to be delivered for each subfield. With a lung tumor, a CTV that executes a quasi periodic motion z(t) does not receive ${\Phi}$ (x,y), but rather an Effective Incident Fluence EIF(x,y). We numerically evaluate the EIF(x,y) from a given DMLC file by a coordinate transformation to the Target's Eye View (TEV). In the TEV coordinate system, the CTV itself is stationary, and the MLC is seen to execute a motion -z(t) that is superimposed on the DMLC motion. The resulting EIF(x,y)is inputted back into the dose calculation engine to estimate the 3D dose to a moving CTV. In this study, we model respiratory motion as a sinusoidal function with an amplitude of 10 mm in the superior-inferior direction, a period of 5 seconds, and an initial phase of zero.