Quasi-independent monitor unit calculation for intensity modulated sequential tomotherapy

Transcription

1 JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS, VOLUME 3, NUMBER 2, SPRING 2002 Quasi-independent monitor unit calculation for intensity modulated sequential tomotherapy Jen-San Tsai, 1,2 Mark J. Engler, 1,2 and James Liu 1 1 Department of Radiation Oncology, Tufts University School of Medicine, Boston, 136 Harrison Avenue, Boston, Massachusetts New England Medical Center Hospitals (a Lifespan Partner), Boston, Massachusetts Received 31 July 2001; accepted for publication 14 December 2002 The number of linac monitor units MU from intensity modulated sequential tomotherapy IMST is substantially larger than the MU delivered in conventional radiation therapy, and the relation between MU and dose is obscure due to complicated variation of the beam intensities. The purpose of this work was to develop a practical method of verifying the MU and dose from IMST so that the MU of each arced beam could be double-checked for accuracy. MU calculations for 41 arced beams from 14 IMST patients were performed using the variables of vane open fraction time, field size, target depth, output factor, TMR, and derived intensity distribution. Discrepancy between planned and checked MU was quantified as 100 (MU cal MU plan )/MU plan percent. All 41 discrepancies were clustered between 5% to 4%, illustrated in a Gaussian-shaped histogram centered at % standard deviation indicating the present MU calculations are in agreement with the planned expectations. To confirm the correctness of the present calculated MUs of the IMST plans, eight of the calculated IMST plans are performed dose verifications using their hybrid plans, which are created by transporting patient s IMST plan beams onto a spherical polystyrene Phantom for dose distribution within the Phantom. The dose was measured with a 0.07 cc ionization chamber inserted in the spherical Phantom during the hybrid plan irradiation. Average discrepancy between planned and measured doses was found to be % with single standard deviation uncertainty. The spread of the discrepancies of present calculated MUs relative to their planned ones are attributed to uncertainties of effective field size, effective planned dose corresponding to each arc, and inaccuracy of quantification of scattered dose from adjacent arced beams. Overall, the present calculation of MUs is consistent with what derived from treatment plans. Since the MUs are verified by actual dose measurements, therefore the present MU calculation technique is considered adequate for double-checking planned IMST MUs American College of Medical Physics. DOI: / PACS number s : j, a Key words: independent calculation, monitor unit, sequential tomotherapy, intensity modulated radiation therapy I. INTRODUCTION Intensity modulated sequential tomotherapy IMST, NOMOS Corp., Sewickley, PA is the technique using both field shaping and intensity variation to configure a dose distribution 3D conformal to the target to minimize the dose toxicity to the normal tissues adjacent to the target. Unfortunately, at the present time, there is no technology that enables the radiation source itself to achieve a well-controlled beam intensity pattern in both time sequence and specified space. Therefore, current beam intensity variation mainly relies on an external apparatus to manipulate the opening of the beam either partially or fully in time space. Since the radiation beam maintains its Õ2002Õ3 2 Õ135Õ19Õ$ Am. Coll. Med. Phys. 135

2 136 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 1. IMST beam ares, couch indices, and the coordinate system, T: target. constant monitor units MU delivery rate, this implies that some radiation beams, which are presumably to irradiate some specific anatomic site and the neighboring tissues, are sometimes blocked off or partially attenuated. Thus to deliver a prescribed dose to the patient, the IMST linac runs a substantially larger MU and simultaneously far more unpredictably than the MU delivered in conventional radiation therapy. In fact, the relationship between MU and dose for IMST is determined by a set of multiple intertwined procedures ranging from imaging, computer networking, treatment planning, optimization, and treatment delivery 1 9 and is entirely obscured to the radiation physicist, radiation oncologist, and therapist delivering the treatment. In compliance with the QA of dose delivery, especially in IMST dose escalation treatments, AAPM report #38, 10 accompanied by the policy of Radiological Physics Center, invokes the importance and necessity of double checking of MUs prior to treatment in radiotherapy. In the past, Tsai et al., 11 Low et al., 12 Teh et al., 13 and Verellen et al., 14 have performed integrally IMST plan verification, indirectly fulfilling some MUs second check, using a humanoid phantom prior to IMST treatment. This integrated plan verification procedure requires about four hours per patient. To reduce the time required for IMST plan verification, some QA procedure may not have to proceed on every individual patient because of known slowly varying quality condition of the radiotherapy mode without pose any significant risk in a short period of time. However, to maintain the high vigilant QA, MUs second check, being directly related to dose delivery and varying from patient to patient, is perhaps the non-compromising alternative QA needed. This article addresses this issue. Reports of MU calculations for verifications of multi-leaf collimator MLC intensitymodulated radiation therapy IMRT as a dosimetry QA have been done by Boyer et al. 15 and Kung et al. 16 Recently, Ayyangar et al. 17 have done an independent dose calculation for MIMiC collimated IMRT NOMOS Corp., Sewickley, PA by using the plan optimized intensity beamlets and subsequently repeated dose calculation on the same anatomy by simple dosimetric formulas to check the plan dose distribution. This article, in the same wavelength, independently develops a MU calculation based on some modeled parameters is to dedicate for IMST using the MIMiC collimator. In the present article, the MU calculation is based on the input data composed of measured values and some modeled estimations, for example the average dose, which is inevitably embedded in some intrinsic uncertainty. This implies that we are dealing with the macroscopic, instead of voxel-wised microscopic quantities. Since the purpose of the present work is for vigilant checking of the MUs derived from the IMST plan, our task is to scrutinize the MUs of IMST plan so that they are not far off from the present calculated ones. In other words, a high accuracy of the calculated MU may not be essential albeit we are further exploring for better accuracy.

3 137 Tsai, Engler, and Liu: Quasi-independent monitor unit II. MATERIALS AND METHODS A. Glossary of terms To help readers understand the article, some terminology, acronyms, and symbols involved with the article are specified, including their units as follows: Coordinates system The coordinate system (x, y,z) of NOMOS IMST and couch index are sketched in Fig. 1. The couchindex beam arc j, with index coordinate at z j, being synonymous to z coordinate, is used to describethe couch movement of the patient relative to the isocenter. D Dose in Gray Gy. d Depth in cm. DVH Dose volume histogram. Gx Gantry rotational axis. IMRT Intensity modulated radiation therapy. IMST Intensity modulated sequential tomotherapy, a special category of IMRT. IVSC Inverse square correction factor relative to the SAD source to axial distance. MIMiC Multi-vane intensity modulated collimator. 40 vanes, in two rows, operated in binary mode, are used for radiation intensity controlling. m 1/2 Average field width spanned by the MIMiC at the Gx in the cross-plane direction. MU The monitor unit when beam is on. Normally this MU is calibrated so that 1MU 1 cgy at the specific depth. O/P Output factor, is the linear accelerator linac dose ratio at the maximum dose depth, cm 2 field size to any field size. PCF Practical Calibration factor. r TMR A point in 3D space with implicit coordinates (x,y,z) see Fig. 1 Tissue maximum ratio, a dose ratio for a given field sizes at any depth to the maximum dose depth. (t i ) j Fraction of time that remains opened in the rotation of gantry angle indices i 1,2,3,...,31 30,40,50,...,,...,330 for a given vane j, j 1,2,3,...,40 of the MIMiC. This is equivalentto the relative beam intensity delivered by vane j at gantry angle i. W 1cm,W 2cm The radiation field widths of the MIMiC at 1 cm, 2 cm mode respectively projected along the gantryrotational axis at the isocenter height. These are the widths that are measured at the couch height of isocenter, and the couch should move along the gantry rotational axis for the adjacent arcs irradiation. Gantry angle. Couch angle. The average field size in cm 2. %DD Percentage of depth dose. B. Mathematical equations for arc MU calculation In the gantry rotating IMST, dose delivery to the patient is performed using a multi-vane switching MIMiC collimator, and is treated arc by arc. 4,5 Each arc has its own MU derived from the inverse plan optimization Corvus version 3.0, NOMOS Corp., Sewickley, PA. In the present article, we temporarily adopt a hypothesis, to be justified subsequently, that the dose distributed in the transverse plane coplanar to the beam arc central axial plane is overwhelmingly contributed by the disc beam arc as implicitly illustrated in Fig. 1. Theoretically, each beam arc contributes some dose to an anatomic location r. However, due to the drastical drop-off of the beam profile beyond a half radiation field width i.e., either 0.5 W 1cm, or 0.5W 2cm along the gantry rotational

4 138 Tsai, Engler, and Liu: Quasi-independent monitor unit axis from the beam central line the location of the couch index coordinate z j of arc j, itis reasonable to say that the dose distributed on the transversal plane z j is primarily contributed by the beam arc j. The dose contribution from its adjacent arcs is very minor mainly the scattering and leakage. In other words, to evaluate the MU j value for beam arc j, we have to identify the z coordinate in the IMST plan of anatomic CT image so that it matches the beam arc j index coordinate z j. Thus, at any given anatomic location r in the transverse plane z j, the dose (r) j and the MU j from this arc j have the relationship depicted by the general dosimetry equation. Dose r j D r j MU j Output factor field size relative to cm 2 j t fraction of time vane opened j TMR depth, field size j Beam profile r j or IVSC O scattered dose from adjacent arcs O total leakage i MU j i O/P i j t i j TMR d i, i j Beam profile i,r j IVSC O scattered dose from adj. arcs leakage MU F I 20 / F I j O/P j t j TMR d, j where Beam profile,r j IVSC O dose scattered from adj. arcs O j leakage j, 1 i the running index for gantry angle, in 10 step. That is, for example, 1,2,3,...,i,...,31 30,40,50,...,,...,330. The average value of the parameter, etc. D(r) j Total dose at the anatomic location r corresponding to arc j index coordinate z j in the z direction. The factor ( F I 20 )/( F I ) j is accounted for the offset of MU j because Corvus IMST system the first 10 interval of initial i and final 10 interval of final F of gantry angles, all vanes of MIMiC are still in close status due to the angular acceleration and deceleration of the gantry rotation. TMR is assumed that TMR(d i, i ) j TMR( d, ) j, which is valid within the uncertainty 1% based on our calculation. The zero orders O dose scattering dose from adjacent beams and O ( j leakage j ) can be explored by inserting a film at coordinate z j to the humanoid phantom and performing the IMST plan irradiation without the beam arc j. For IMST using the NOMOS MIMiC collimator, each term in Eq. 1 cannot be quantified uniquely because of variable field size i, beam intensities ( t i ), and depth d i. Therefore, the present MU calculation is done implicitly by some modeled parameters expressed in the Eq. 1. For simplicity, the subscript j for arc j will be dropped by default. As mentioned above, the purpose of the present MU calculation is used as a secondary verification of MU derived from the IMST plan for radiotherapy QA. We did not analyze in detail the high order dose components of O dose scattering and O leakage for every arc. It is also understood that the dose is accumulated only when the vanes open, which define the dose area. C. Estimation of average vane-open time The average vane-open time t is quantified using the percentage sum of each vane i s open time fraction t i and averaged over the total vanes number that have involved in the beam intensity

5 139 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 2. The fraction of each vane s open time. The darker indicates longer open time, while the white vanes indicate all closed. The numerical on the left side shows the numerical fraction of time opened in percentage. The sketch shows the tumor depth variation as gantry rotates during IMST treatment. modulation and field shaping through the entire arc. This can be illustrated in Fig. 2. The summation is performed from gantry initial angle I to gantry stop angle F : F t k t k / k m k, 2 I where t i s are numerical percentage printout from the IMST plan and are displayed side by side besides the corresponding vane pattern, and m k 1 if t k 0 k 1,2,3,..., if t k 0 This extracted average time t is, in fact, the beam on fraction.

6 140 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 3. The tumor depth is varied as the gantry rotates. This figure illustrates how the average depth of the tumor is extracted in the present work of MU calculation. The radial lines are drawn diverging from the tumor center. All the radial lines are in equal angular 10 intervals. Each radial quantifies a depth from the tumor center to the skin. The average depth is obtained by summing up all the radial line depths and divided by the total number of radial line depths. D. Average field size, ŠA, estimation For a given specific vane pattern within a 10 interval of gantry rotation, and definition of m k ( ) in Eq. 3, the average number of vane pairs open across the MIMiC port of the average field size during gantry rotation from gantry angle I to F, as illustrated in Fig. 2 would be m 1/2 40 k 1 t k 40 k 1 m k k 1 t k Thus the averaged field size, for the 2 cm mode MIMiC would be, W 2cm m 1/2 cm 2, where W 2cm 3.4 cm is the index width of 2 cm mode MIMiC along the gantry rotational axis, Gx. The real time fraction, t, comes from the fact that any observed vane open pattern size is probabilistic in current IMST, and a function of the vane open time. Similarly, for 1 cm mode MIMiC with index width W 1cm 1.7 cm the averaged field size would be W 1cm m 1/2 cm 2. E. Average depth, Šd The average depth d for a given anatomic arc is illustrated in the sketch of Fig. 3. An equal angular divergent lines emitting from the common point of the tumor geometric center is fabricated in a transparent sheet. The transparent sheet prints the external anatomic contour and the tumor obtained by overlaying on the transparent sheet on the CT image of some specific couch index coordinate z i along the superior-inferior (S/I) direction. From the transparent sheet, the external anatomic contour is delineated. From the external anatomic contour, the depth, d( ) ata given angle from the tumor geometric center to the anatomic surface is evaluated by a rule. The average depth, d, is thus obtained by the formula, d d scale factor of the image. Total divergent lines assessed 5 The uncertainty of d includes measurement of d( ), noncoincidence of Gx and the target center, and adjacent image slices. The overall uncertainty of d is expected up to 5%, taking consideration of contour difference of its adjacent CT image slices.

7 141 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 4. In the present modeled D, the average dose around the tumor within the field size, we assumed D da D ( m 1/2 /2) 2. F. Average dose, ŠD, estimation To quantify the average tumor dose, D j, on the axial CT image plane directly corresponding to the arc j, we use the image contour coplanar to the beam arc j at z coordinate index z j and display its isodose distribution derived from the IMST plan. The procedure is sketched in Fig. 4. The average dimension of the isodose distribution m 1/2 is equal to the average vane pattern size in the transverse MIMiC direction of the arc, as shown in Fig. 2. This isodose distribution is assumed entirely radiated by the shown intensity vane patterns of the arc j. Next, in this axial isodose distribution plane of arc j, we go to find out the maximum isodose level, D h, by assigning an isodose level such that a point or a tiny spot, instead of a spread out isodose distribution, appears. Afterwards, we go to find the isodose level D l such that the isodose contour area A is equal to ( m 1/2 /2) 2. This is based on an assumption that an average width of field size spanned from the MIMiC projected and perpendicular to Gx is m 1/2. With this averaged field width m 1/2, the arc for a gantry 300 rotation will fabricate an isodose area equivalent to a circle of radius m 1/2. Our quantification of the average dose, D, is the isodose between D h and D l at the arc j,

8 142 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 5.D l, D, andd h isodose distributions. FIG. 6. Measured output factors of the MIMiC for IMST beam dose calculation.

9 143 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 7. A typical beam profile of 2 cm mode MIMiC. The small penumbra indicates that the radiation extending to the adjacent index positions is negligible. D A D pixel area A D Total A l D h D l /2 D h D l /2. However, of the D, a small fraction comes from scattered dose due to adjacent arcs. So, the real dose effectively contributed from the arc j beam is D D. This fractional dose D will be estimated from film measurements. The pursued average dose, D, is simply the dose of the mean 6 FIG. 8. The setup of film and Rando Phantom for the study of the radiation dose contributing to the adjacent axial plane.

10 144 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 9. a The setup of the spherical phantom and mini-ionization chamber for absolute dose verification between the IMST plans and that is measured. b The dose distribution of the hybrid plan using a Spherical Phantom. isodose levels of D h and D l. This is illustrated in Fig. 5, in which the isodose line of D h is shown as near a spot circle, and isodose line D l is a contour whose area is qualitatively equal to ( m 1/2 /2) 2. G. Output factor The output factors of the IMST beam are shown in Fig. 6 for 1 cm and 2 cm mode of MIMiC respectively. These output factors for specific field sizes are measured using films, Si-diode, and TLD for mutual consistency. If the equivalent IMST field size is not equal to any measured points shown in Fig. 5, then an interpolated output factor is made by the dashed curve. H. Tissue maximum ration TMR The TMR of IMST is deduced from the percent depth dose (%DD) curves not shown in this article, which are stored in the dosimetry base of the Corvus IMST computer. The relationship between TMR and %DD is expressed as

11 145 Tsai, Engler, and Liu: Quasi-independent monitor unit FiG. 9 Continued. FIG. 10. The film dose distribution accumulated from its adjacent arcs.

12 146 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 11. Discrepancies of all the IMST plans. FIG. 12. The summed extrapatient MU discrepancies between the calculated MU and that from the IMST planned. The dashed curve is the Gaussian-fitted curve and the vertical dashed line is the center of the Gaussian distribution.

13 147 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 13. The average dose D quantification is sometimes complicated by 2 the high dose gradient of the image which is near the tumor edge, and 2 multiple targets contributed the dose to their high dose areas. TMR d, %DD d, S p d /S p m 100 d / 100 d max 2, where (S p d /S p m ) is the phantom scattering ratio for the field size projected to different depths of d and d max. For simplicity, we have set (S p d /S p m ) 1. The squares of spherical parentheses account for the term IVSC. 7 I. Dose contribution from adjacent arcs In current IMST, the separations between each couch index at coordinates (z j ) and others are in equal distance, either W 1cm or W 2cm depending on the MIMiC mode used. According to the IMST beam profiles in z-direction, Fig. 7, the radiation intensity at any z-coordinate beyond the beam center by W 1cm or W 2cm is negligibly small albeit the scattering dose from the adjacent arcs is not totally neglected as expressed in Eq. 1. In order to study the dose contribution from the adjacent arcs to the MU assessing arc, a film Kodak, X-Omat V is inserted to Rando Phantom as shown in Fig. 8. The Rando Phantom is set up in such a manner as described in Tsai et al. 1 The

14 148 Tsai, Engler, and Liu: Quasi-independent monitor unit TABLE I. Absolute dose verified from the hybrid plans of the spherical polystyrene phantom. Trial # Verified dose Gy Hybrid planned dose Gy Ratio verified/hybrid film is located at the beam arc s couch index z coordinate. The phantom is then irradiated by IMST plan for arcs j 1, j 1. The dose at tumor target T g contributed from arcs j 1, j 1 is quantified by film scan. J. IMST plan dose verification At the moment, we do the IMST arc MU calculations and then compare them with the planned obtained MU. Logically, this is based on the assumption that the MUs derived from the IMST plans are accurate, and correct. To confirm that the MUs from the IMST plans are true, several IMST plans have their MUs verified using a hybrid plan on a spherical phantom. The general IMST hybrid plan technique and its components have been described in Tsai et al. 10 In the present work, a spherical polystyrene phantom used for dose calibration in Leksell Gamma Knife stereotactic radiosurgery facility in conjunction with a Capintec mini ionization chamber Model PR- 05P, 0.07 cc are employed see sketch, Fig. 9 a. The spherical phantom and mini-ionization FIG. 14. IMST beam arcs and its couch indices.

15 149 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 15. IMRT treatment sheet. chamber are used to verify the absolute dose achieved from the hybrid plan Fig. 9b. The measured dose from this hybrid plan is compared with the dose expected from the IMST hybrid plan. This will test the reliability and accuracy of the association of the plan MUs with the absolute dose. III. RESULTS AND DISCUSSION A. Scattered dose from adjacent arcs This has wide uncertainty depending on the number of arcs for each plan, tumor size, and the index width mode of MIMiC. Figure 10 shows one of our studies of the films embedded within the Rando phantom, the scattered dose focused to the area, where the D is estimated between 6 8 cgy/fx, accounting for D 6 8/270 cgy or 2 3% of the D. In fact, the scattered dose from the adjacent arcs is slightly spread widely over the anatomic space. This implies the quantification of the scattered dose from adjacent arcs is slightly not well defined, and thus incurring a wider range in discrepancy up to 2 4% depending on the target size, target number and dose prescription option to targets. Nevertheless, this discrepancy due to dose scattering from adjacent arcs is smaller than the propagated or accumulated uncertainties from variable fields size i, depth d i, time open fraction t i, quantified dose D(r) j, and O/P due to vane switching. 5 B. Calculated MU relative to the planned Fourteen plans were used to test the current model MU calculations. All the cases and their MU discrepancies between the present calculation work and the IMST plans were shown in Fig. 11. A typical example of the present MU calculation is shown in the Appendix. In each case, we graphically display its value in percentage / 2% as experimental error. For example, in the first arc of the GBM patient, the obtained discrepancy 4% will be displayed in the horizontal axis of histogram between 2 to 6%. All such discrepancies are summed up to a combined histogram for the extrapatient MU discrepancy, as shown in Fig. 12. Overall, the average of the total extrapatient MU discrepancy yields %, where 3.5% is the standard deviation. Overall, the global MU discrepancy between the present calculations and the planned is in agreement. The wide spread of discrepancy is attributed to the following reasons: 1 Average dose, D, uncertainty. This D is quantified based on the model that the dose volume histogram DVH between the region of D l and D h, shown in Fig. 4, is triangle in shape. As a result, the D is a mean value of D l, and D h. Any DVH deviated from the aforementioned shape will alter the D, and automatically change the calculated MU value. In several cases, the beam index position lies either near the tumor edge, where the dose distribution shows a high

16 150 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 16. The CT image, including its isodose distribution, closest to the couch index location (Z mm) of the beam arc #1. gradient, or the between two tumor targets see Fig. 13, which further complicates D quantification. In Fig. 13, the arcs delivering dose for upper target share portioned MUs from the arcs that deliver dose to lower target, and verse versa. Under such circumstances, higher discrepancy between our calculated MUs and the IMST planned is expected because of less accurate dose quantification from the arcs being responsible for the dose coverage of the lower target. 2 Average fields size dimension,, uncertainty. The, being probabilistic in nature as mentioned above and depicted in Fig. 2, has incorporated some extent of uncertainty. The uncertainty propagates subsequent uncertainty of the average dose assessment. 3 Index z coordinate uncertainty for D quantification. In all 41 extrapatient IMST arcs calculations of MUs, we evaluate the D from the isodose distribution of the plan on an axial image at the z-coordinate specified by the IMST treatment sheet. This is based on the assumption that the vane pattern opened in MIMiC, no matter 1 or 2 cm mode, includes both rows. This however, is not true in some cases near the tumor edge, where only one row, either superior or inferior of the MIMiC is involved for beam field intensity modulating and shaping. In order to minimize the z-coordinate deviation for one row MIMiC open, the axial image of its z coordinate is adjusted. For example, in 1 cm mode MIMiC, if the arc index z coordinate is z 1, and only the superior row is involved, then the D quantification is not at axial image of z 1, instead, it is of z mm. 4 The average dose, D j, in the CT image of z j, being overwhelmingly contributed by the beam arc j. This is true for single couch angle setup. However, for multiple couch angle setup, the same D j may be crossed over by beam arc j of couch angle, and beam arc j of couch angle setup. The portion contributions to D j j from beam arc j and j is more complicated than we expected. Therefore, the present quasi-calculation of MU using Eq. 1 is infeasible for multiple couch IMST setup. 5 As aforementioned, we have simplified the present calculation in Eq. 7 that the phantom scattering ratio S d p /S m p 1. This may incur an uncertainty up to about 2%. 6 In our MUs calculations using Eq. 1, no heterogeneous density correction is made. Nev-

17 151 Tsai, Engler, and Liu: Quasi-independent monitor unit FIG. 17. Vane patterns of beam arc #1 at different gantry angle. ertheless, NOMOS s Corvus planning software, albeit its potential, still remains unproven for the correction. Thus, the MU discrepancy due to density heterogeneity remains inconclusive. Detailed dose calculation algorithm of Corvus was reported by Bleier et al. 18 For a given patient, and given beam arc i.e., given couch index, MU discrepancy between the present calculations and the planned does raise comparative uncertainty. It doesn t mean any flaw of the present methodology. Instead, it mainly attributes to the modeled and probabilistic quantities for some parameters proceeded in the present work. C. Dose verification of hybrid plan The result of absolute dose measurements from spherical polystyrene phantom of eight hybrid plans is listed in Table I. The averaged discrepancy between measured doses and the planned is

18 152 Tsai, Engler, and Liu: Quasi-independent monitor unit less than 1%, indicating, overall, in good agreement. The standard deviation, roughly about 4%, displays a slight spread of the measured data relative to the planned ones. This is due to the fact that some of the locations of dose measurements within the phantom are in high dose gradient region that makes the discrepancies between the measured ones and the planned expectations slightly high. Any discrepancy if exists, however, may include many factors. First, in many cases, the hybrid plan isodose distribution near the spherical phantom center, where the ionization chamber is located, is in a high dose gradient. Second, IMST planning software uses a Practical Calibration Factor PCF to scale the absolute dose readings. PCF was assigned by actual dose measurements within the humanoid phantom for several simulated targets and compare with the planned ones. PCF is slightly correlated to the target s location relative to the gantry rotational axis Gx. In our initial assignment of PCF during the commission, several targets at various locations were measured to offset the deviation and obtain an average PCF. However, the patient s tumor locations are still randomly distributed within the anatomic space, and it is not surprising that there exists some dose uncertainty between the real measured dose and the planned one due to slight deviation of PCF. Third, the ionization chamber, more or less, has stem effect for dose detection efficiency on the IMST beams. IV. CONCLUSIONS As IMST dose, dose distributions, and MU have been extensively verified with measurements in anthropomorphic and other phantoms using the Corvus hybrid plan tool, the quasi-independent monitor unit verification may be applied to simplify QA cost overhead in the implementation of IMST. The current MU calculation can be used for an IMST treatment second check and thus as a tool of QA for IMST dose delivery. Because the MU calculation cannot vigilantly pick up other hidden errors, such as patient setup during treatment, hybrid phantom verification will continue to be performed as supplementary for patients not belonging to classes of anatomy and/or delivery techniques already verified in phantom for prior patients. ACKNOWLEDGMENTS We greatly acknowledge the NOMOS Corporation for providing information that leads to get MIMiC time open fraction. Additionally, we thank for Mark J. Rivard, Ph.D. for English corrections, and Ali Mahmoudieh, M.D. for some calculations checking. APPENDIX: ONE EXAMPLE OF QUASI MU CALCULATION FOR BEAM ARC #1 OF IMST PLAN FOR PT. G.K. Step 1. Identifying the beam arcs and their corresponding couch indices in the Z direction, as shown in Fig. 14. Step 2. Print a copy of the treatment sheet of the IMST plan #1519, as shown in Fig. 15. The treatment sheet displays the alignment point relative couch movement of each beam arc couch indices in terms of the Z coordinate. The couch index of beam arc 1 is Z mm. This means that the transverse central beam of the MIMiC is in coincidence with transverse plane Z mm, i.e., the CT image of Z mm, for gantry 30 to 330 rotation. The MU for this beam arc is 331. Step 3. The CT image whose Z coordinate closes most to the couch index of beam arc # mm is the axial CT image #297 the image registration number from Corvus planning algorithm at Z 58.0 mm, see Fig. 16. Step 4. Calculation of the average depth d of the target based on the CT image Fig. 16 is 6.06 cm. The calculation steps are illustrated in Fig. 3 using Eq. 5. Step 5. Vane patterns of beam arc #1 is shown in Fig. 17, from which the average fraction time of vanes open is t , shown in Fig. 2 from Eq. 2.

19 153 Tsai, Engler, and Liu: Quasi-independent monitor unit Step 6. From the vane pattern of Fig. 17, we are able to calculate the m 6.0 cm using Eq. 4 as shown in Fig. 4. Step 7. The average dose D isodose line whose area (6.0) 2 cm 2 is cgy/fx, shown in Fig. 16. From this dose, we assume 4 5 % is attributed from the scattering doses from beam arcs #2, 3, and 4 and radiation leakage. Step 8. The O/P factor for field size of cm 2 would be The beam width, i.e., the field length projected with MIMiC field port along the gantry rotational axis is 3.4 cm for 2 cm MIMiC set-up, which was used for the IMST plan #1519. Step 9. The TMR ( d 6.0 cm, cm 2 ) Step 10. Beam profile central-axis-ratio 1, i.e., assuming the average of the dose distributed to the adjacent axial CT image #297 by the beam arc #1 has a beam profile central-axis-ratio of 1. Step 11. MU calculated /( ) MU using Eq. 1. Step 12. The MU of this beam arc #1, derived from the IMST plan #1519, is 331. Therefore, the discrepancy is 5%. 1 R. Mohan, G. S. Mageras, B. Baldwin, L. Brewster, G. J. Kutcher, S. Leibel, C. M. Burman, C. C. Ling, and Z. Fuks, Clinically relevant optimization of 3-Dradiation treatments, Med. Phys. 19, G. S. Mageras, K. C. Podmaniczky, R. Mohan, A model for computer-controlled delivery of 3D conformal treatments, Med. Phys. 19, S. A. Leibel, R. Heimann, G. J. Kutcher, M. J. Zelefssky, C. M. Burman, E. Melian, J. P. Orazem, R. Mohan, T. T. LoSasso, Y.-C. Lo, J. A. Wiseberg, D. S. Chapman, L. C. C. Ling, and Z. Fuks, Three-dimensional conformal radiation therapy in locally advanced carcinoma of the prostate: Preliminary results of a phase I dose escalation study, Int. J. Radiat. Oncol., Biol., Phys. 28, M. P. Carol, Peacock: A system for planning and rotational delivery of intensity-modulated fields, Int. J. Imaging Syst. Technol. 6, M. P. Carol, Integrated 3-D conformal multivane intensity modulation delivery system for radiotherapy, in Proceedings of the Eleventh International Conference on the Use of Computers in Radiation Therapy, Stockport, Cheshire, UK, edited by A. R. Hounsell, pp M. Oldham and S. Webb, Intensity-modulated radiotherapy by means of static tomotherapy, Med. Phys. 24, T. R. Mackie, T. Holmes, S. Swerdloff, P. Reckwerdt, J. O. Deasy, J. Yang, B. Paliwal, and T. Kinsella, Tomotherapy: A new concept for delivery of dynamic conformal radiotherapy, Med. Phys. 20, J. M. Michaliski, R. K. Sur, W. B. Harms, and J. A. Purdy, Three dimensional conformal radiation therapy in pediatric parameningeal rhabdomyosarcomas, Int. J. Radiat. Oncol., Biol., Phys. 33, J.-S. Tsai, D. E. Wazer, M. N. Ling, J. K. Wu, M. Fagundes, T. DiPetrillo, B. Kramer, M. Koistinen, and M. J. Engler, Dosimetric verification of the dynamic intensity-modulated radiation therapy of 92 patients, Int. J. Radiat. Oncol., Biol., Phys. 40, AAPM Report No. 38, The role of the Physicist in Radiation Oncology, Professional Information and Clinical Relation Committee TG# J.-S. Tsai, M. J. Engler, M. N. Ling, J. K. Wu, B. Kramer, T. DiPetrillo, and D. E. Wazer, A non-invasive immobilization system and related quality assurance for dynamic modulated radiation therapy of intracranial and head and neck disease, Int. J. Radiat. Oncol., Biol., Phys. 43, D. A. Low, J. F. Dempsey, R. Venkatesan, S. Mutic, J. Markman, E. M. Haacke, and J. A. Purdy, Evaluation of polymer gels and MRI as a 3-D dosimeter for intensity-modulated radiation, Med. Phys. 26 8, B. S. Teh, W. Y. Mai, B. M. Uhl, M. E. Augspurger, W. H. Grant, H. H. Lu, S. Y. Woo, L. S. Carpenter, J. K. Chiu, and E. B. Butler, Intensity-modulated radiation therapy IMRT for prostate cancer with the use of rectal balloon for prostate immobilization: acute toxicity and dose-volume analysis, Int. J. Radiat. Oncol., Biol., Phys. 49, D. Verellen and N. Linthout, et al. Initial experience with intensity-modulated conformal radiation therapy for treatment of the head and neck region, Int. J. Radiat. Oncol., Biol., Phys. 39, A. Boyer et al., Theoretical considerations of monitor unit calculations for intensity modulated beam treatment planning, Med. Phys. 26, J. H. Kung, G. T. Y. Chen, and F. K. Kuchnir, A monitor unit verification calculation in intensity modulated radiotherapy as a dosimetry quality assurance, Med. Phys. 27, K. M. Ayyangar et al., Independent dose calculations for the Peacock System, Med. Dosim. 26, A. R. Bleier, M. P. Carol, and B. H. Curran, et al., Dose calculation for intensity-modulated radiotherapy using finite size pencil beam derived from standard measured data, Proceedings of the XII International Conference on the Use of Computers in Radiation Therapy, Salt Lake City, 1997, pp