Design and testing of a simulation framework for dosimetric motion studies integrating an anthropomorphic computational phantom into four-dimensional Monte Carlo

We have designed a simulation framework for motion studies in radiation therapy by integrating the anthropomorphic NCAT phantom into a 4D Monte Carlo dose calculation engine based on DPM. Representing an artifact-free environment, the system can be used to identify class solutions as a function of geometric and dosimetric parameters. A pilot dynamic conformal study for three lesions ( approximately 2.0 cm) in the right lung was performed (70 Gy prescription dose). Tumor motion changed as a function of tumor location, according to the anthropomorphic deformable motion model. Conformal plans were simulated with 0 to 2 cm margin for the aperture, with additional 0.5 cm for beam penumbra. The dosimetric effects of intensity modulated radiotherapy (IMRT) vs. conformal treatments were compared in a static case. Results show that the Monte Carlo simulation framework can model tumor tracking in deformable anatomy with high accuracy, providing absolute doses for IMRT and conformal radiation therapy. A target underdosage of up to 3.67 Gy (lower lung) was highlighted in the composite dose distribution mapped at exhale. Such effects depend on tumor location and treatment margin and are affected by lung deformation and ribcage motion. In summary, the complexity in the irradiation of moving targets has been reduced to a controlled simulation environment, where several treatment options can be accurately modeled and quantified The implemented tools will be utilized for extensive motion study in lung/liver irradiation.

Psychometric characteristics of the Spanish version of instruments to measure neck pain disability

BACKGROUND

The NDI, COM and NPQ are evaluation instruments for disability due to NP. There was no Spanish version of NDI or COM for which psychometric characteristics were known. The objectives of this study were to translate and culturally adapt the Spanish version of the Neck Disability Index Questionnaire (NDI), and the Core Outcome Measure (COM), to validate its use in Spanish speaking patients with non-specific neck pain (NP), and to compare their psychometric characteristics with those of the Spanish version of the Northwick Pain Questionnaire (NPQ).

METHODS

Translation/re-translation of the English versions of the NDI and the COM was done blindly and independently by a multidisciplinary team. The study was done in 9 primary care Centers and 12 specialty services from 9 regions in Spain, with 221 acute, subacute and chronic patients who visited their physician for NP: 54 in the pilot phase and 167 in the validation phase. Neck pain (VAS), referred pain (VAS), disability (NDI, COM and NPQ), catastrophizing (CSQ) and quality of life (SF-12) were measured on their first visit and 14 days later. Patients' self-assessment was used as the external criterion for pain and disability. In the pilot phase, patients' understanding of each item in the NDI and COM was assessed, and on day 1 test-retest reliability was estimated by giving a second NDI and COM in which the name of the questionnaires and the order of the items had been changed.

RESULTS

Comprehensibility of NDI and COM were good. Minutes needed to fill out the questionnaires [median, (P25, P75)]: NDI. 4 (2.2, 10.0), COM: 2.1 (1.0, 4.9). Reliability: [ICC, (95%CI)]: NDI: 0.88 (0.80, 0.93). COM: 0.85 (0.75,0.91). Sensitivity to change: Effect size for patients having worsened, not changed and improved between days 1 and 15, according to the external criterion for disability: NDI: -0.24, 0.15, 0.66; NPQ: -0.14, 0.06, 0.67; COM: 0.05, 0.19, 0.92.

VALIDITY

Results of NDI, NPQ and COM were consistent with the external criterion for disability, whereas only those from NDI were consistent with the one for pain. Correlations with VAS, CSQ and SF-12 were similar for NDI and NPQ (absolute values between 0.36 and 0.50 on day 1, between 0.38 and 0.70 on day 15), and slightly lower for COM (between 0.36 and 0.48 on day 1, and between 0.33 and 0.61 on day 15). Correlation between NDI and NPQ: r = 0.84 on day 1, r = 0.91 on day 15. Correlation between COM and NPQ: r = 0.63 on day 1, r = 0.71 on day 15.

CONCLUSION

Although most psychometric characteristics of NDI, NPQ and COM are similar, those from the latter one are worse and its use may lead to patients' evolution seeming more positive than it actually is. NDI seems to be the best instrument for measuring NP-related disability, since its results are the most consistent with patient's assessment of their own clinical status and evolution. It takes two more minutes to answer the NDI than to answer the COM, but it can be reliably filled out by the patient without assistance.

TRIAL REGISTRATION

Clinical Trials Register NCT00349544.

Effects of organ motion on IMRT treatments with segments of few monitor units

Interplay between organ (breathing) motion and leaf motion has been shown in the literature to have a small dosimetric impact for clinical conditions (over a 30 fraction treatment). However, previous studies did not consider the case of treatment beams made up of many few-monitor-unit (MU) segments, where the segment delivery time (1-2 s) is of the order of the breathing period (3-5 s). In this study we assess if breathing compromises the radiotherapy treatment with IMRT segments of low number of MUs. We assess (i) how delivered dose varies, from patient to patient, with the number of MU per segment, (ii) if this delivered dose is identical to the average dose calculated without motion over the path of the motion, and (iii) the impact of the daily variation of the delivered dose as a function of MU per segment. The organ motion was studied along two orthogonal directions, representing the left-right and cranial-caudal directions of organ movement for a patient setup in the supine position. Breathing motion was modeled as sin(x), sin4(x), and sin6(x), based on functions used in the literature to represent organ motion. Measurements were performed with an ionization chamber and films. For a systematic study of motion effects, a MATLAB simulation was written to model organ movement and dose delivery. In the case of a single beam made up of one single segment, the dose delivered to point in a moving target over 30 fractions can vary up to 20% and 10% for segments of 10 MU and 20 MU, respectively. This dose error occurs because the tumor spends most of the time near the edges of the radiation beam. In the case of a single beam made of multiple segments with low MU, we observed 2.4%, 3.3%, and 4.3% differences, respectively, for sin(x), sin4(x), and sin6(x) motion, between delivered dose and motion-averaged dose for points in the penumbra region of the beam and over 30 fractions. In approximately 5-10% of the cases, differences between the motion-averaged dose and the delivered 30-fraction dose could reach 6%, 8% and 10-12%, respectively for sin(x), sin4(x), and sin6(x) motion. To analyze a clinical IMRT beam, two patient plans were randomly selected. For one of the patients, the beams showed a likelihood of up to 25.6% that the delivered dose would deviate from the motion-averaged dose by more than 1%. For the second patient, there was a likelihood of up to 62.8% of delivering a dose that differs by more than 1% from the motion-averaged dose and a likelihood of up to approximately 30% for a 2% dose error. For the entire five-beam IMRT plan, statistical averaging over the beams reduces the overall dose error between the delivered dose and the motion-averaged dose. For both patients there was a likelihood of up to 7.0% and 33.9% that the dose error was greater than 1%, respectively. For one of the patients, there was a 12.6% likelihood of a 2% dose error. Daily intrafraction variation of the delivered dose of more than 10% is non-negligible and can potentially lead to biological effects. We observed [for sin(x), sin4(x), and sin6(x)] that below 10-15 MU leads to large daily variations of the order of 15-35%. Therefore, for small MU segments, non-negligible biological effects can be incurred. We conclude that for most clinical cases the effects may be small because of the use of many beams, it is desirable to avoid low-MU segments when treating moving targets. In addition, dose averaging may not work well for hypo-fractionation, where fewer fractions are used. For hypo-fractionation, PDF modeling of the tumor motion in IMRT optimization may not be adequate.

Amorphous silicon EPID calibration for dosimetric applications: comparison of a method based on Monte Carlo prediction of response with existing techniques

For EPID dosimetry, the calibration should ensure that all pixels have a similar response to a given irradiation. A calibration method (MC), using an analytical fit of a Monte Carlo simulated flood field EPID image to correct for the flood field image pixel intensity shape, was proposed. It was compared with the standard flood field calibration (FF), with the use of a water slab placed in the beam to flatten the flood field (WS) and with a multiple field calibration where the EPID was irradiated with a fixed 10x10 field for 16 different positions (MF). The EPID was used in its normal configuration (clinical setup) and with an additional 3 mm copper slab (modified setup). Beam asymmetry measured with a diode array was taken into account in MC and WS methods. For both setups, the MC method provided pixel sensitivity values within 3% of those obtained with the MF and WS methods (mean difference<1%, standard deviation<2%). The difference of pixel sensitivity between MC and FF methods was up to 12.2% (clinical setup) and 11.8% (modified setup). MC calibration provided images of open fields (5x5 to 20x20 cm2) and IMRT fields to within 3% of that obtained with WS and MF calibrations while differences with images calibrated with the FF method for fields larger than 10x10 cm2 were up to 8%. MC, WS and MF methods all provided a major improvement on the FF method. Advantages and drawbacks of each method were reviewed.

Electron dosimetry of angular fields

Shaping electron fields through the use of lead cut-outs may result in there being acute angles in part of the field. Using both experimental techniques and EGSnrc Monte Carlo simulations an investigation was carried out to determine the dosimetric consequences of this. Measurements were made to investigate how the field dose was related to the angle between adjacent sides in the cut-outs. The study involved two electron energies (9 MeV and 12 MeV) and source-skin distances (SSDs) in the range 1000-1100 mm. For angles less than about 120 degrees the dose received in the angular region decreased significantly, the effect being more pronounced at 12 MeV than at 9 MeV, and at longer SSDs. The planar shapes of the Monte Carlo dose distributions agreed with those experimentally determined to within +/-1.5 mm at 9 MeV and +/-1.0 mm at 12 MeV, demonstrating the validity of using such calculations for this purpose. Graphs are presented which may help in the prospective assessment of the dose reductions likely to be incurred.

An investigation of dose calculation accuracy in intensity-modulated radiotherapy of sites in the head & neck

Knowledge of the accuracy of dose calculations in intensity-modulated radiotherapy of the head and neck is essential for clinical confidence in these highly conformal treatments. High dose gradients are frequently placed very close to critical structures, such as the spinal cord, and good coverage of complex shaped nodal target volumes is important for long term-local control. A phantom study is presented comparing the performance of standard clinical pencil-beam and collapsed-cone dose algorithms to Monte Carlo calculation and three-dimensional gel dosimetry measurement. Calculations and measurements are individually normalized to the median dose in the primary planning target volume, making this a purely relative study. The phantom simulates tissue, air and bone for a typical neck section and is treated using an inverse-planned 5-field IMRT treatment, similar in character to clinically used class solutions. Results indicate that the pencil-beam algorithm fails to correctly model the relative dose distribution surrounding the air cavity, leading to an cverestimate of the target coverage. The collapsed-cone and Monte Carlo results are very similar, indicating that the clinical collapsed-cone algorithm is perfectly sufficient for routine clinical use. The gel measurement shows generally good agreement with the collapsed-cone and Monte Carlo calculated dose, particularly in the spinal cord dose and nodal target coverage, thus giving greater confidence in the use of this class solution.

A quantitative study of IMRT delivery effects in commercial planning systems for the case of oesophagus and prostate tumours

This study focuses on understanding the impact of intensity-modulated radiotherapy (IMRT) delivery effects when applied to plans generated by commercial treatment-planning systems such as Pinnacle (ADAC Laboratories Inc.) and CadPlan/Helios (Varian Medical Systems). These commercial planning systems have had several version upgrades (with improvements in the optimization algorithm), but the IMRT delivery effects have not been incorporated into the optimization process. IMRT delivery effects include head-scatter fluence from IMRT fields, transmission through leaves and the effect of the rounded shape of the leaf ends. They are usually accounted for after optimization when leaf sequencing the "optimal" fluence profiles, to derive the delivered fluence profile. The study was divided into two main parts: (a) analysing the dose distribution within the planning-target volume (PTV), produced by each of the commercial treatment-planning systems, after the delivered fluence had been renormalized to deliver the correct dose to the PTV; and (b) studying the impact of the IMRT delivery technique on the surrounding critical organs such as the spinal cord, lungs, rectum, bladder etc. The study was performed for tumours of (i) the oesophagus and (ii) the prostate and pelvic nodes. An oesophagus case was planned with the Pinnacle planning system for IMRT delivery, via multiple-static fields (MSF) and compensators, using the Elekta SL25 with a multileaf collimator (MLC) component. A prostate and pelvic nodes IMRT plan was performed with the Cadplan/Helios system for a dynamic delivery (DMLC) using the Varian 120-leaf Millennium MLC. In these commercial planning systems, since IMRT delivery effects are not included into the optimization process, fluence renormalization is required such that the median delivered PTV dose equals the initial prescribed PTV dose. In preparing the optimum fluence profile for delivery, the PTV dose has been "smeared" by the IMRT delivery techniques. In the case of the oesophagus, the critical organ, spinal cord, received a greater dose than initially planned, due to the delivery effects. The increase in the spinal cord dose is of the order of 2-3 Gy. In the case of the prostate and pelvic nodes, the IMRT delivery effects led to an increase of approximately 2 Gy in the dose delivered to the secondary PTV, the pelvic nodes. In addition to this, the small bowel, rectum and bladder received an increased dose of the order of 2-3 Gy to 50% of their total volume. IMRT delivery techniques strongly influence the delivered dose distributions for the oesophagus and prostate/pelvic nodes tumour sites and these effects are not yet accounted for in the Pinnacle and the CadPlan/Helios planning systems. Currently, they must be taken into account during the optimization stage by altering the dose limits accepted during optimization so that the final (sequenced) dose is within the constraints.

Assessing the effect of electron density in photon dose calculations

Photon dose calculation algorithms (such as the pencil beam and collapsed cone, CC) model the attenuation of a primary photon beam in media other than water, by using pathlength scaling based on the relative mass density of the media to water. In this study, we assess if differences in the electron density between the water and media, with different atomic composition, can influence the accuracy of conventional photon dose calculations algorithms. A comparison is performed between an electron-density scaling method and the standard mass-density scaling method for (i) tissues present in the human body (such as bone, muscle, etc.), and for (ii) water-equivalent plastics, used in radiotherapy dosimetry and quality assurance. We demonstrate that the important material property that should be taken into account by photon dose algorithms is the electron density, and not the mass density. The mass-density scaling method is shown to overestimate, relative to electrondensity predictions, the primary photon fluence for tissues in the human body and water-equivalent plastics, where 6%-7% and 10% differences were observed respectively for bone and air. However, in the case of patients, differences are expected to be smaller due to the large complexity of a treatment plan and of the patient anatomy and atomic composition and of the smaller thickness of bone/air that incident photon beams of a treatment plan may have to traverse. Differences have also been observed for conventional dose algorithms, such as CC, where an overestimate of the lung dose occurs, when irradiating lung tumors. The incorrect lung dose can be attributed to the.incorrect modeling of the photon beam attenuation through the rib cage (thickness of 2-3 cm in bone upstream of the lung tumor) and through the lung and the oversimplified modeling of electron transport in convolution algorithms. In the present study, the overestimation of the primary photon fluence, using the mass-density scaling method, was shown to be a consequence of the differences in the hydrogen content between the various media studied and water. On the other hand, the electron-density scaling method was shown to predict primary photon fluence in media other than water to within 1%-2% for all the materials studied and for energies up to 5 MeV. For energies above 5 MeV, the accuracy of the electron-density scaling method was shown to depend on the photon energy, where for materials with a high content of calcium (such as bone, cortical bone) or for primary photon energies above 10 MeV, the pair-production process could no longer be neglected. The electron-density scaling method was extended to account for pair-production attenuation of the primary photons. Therefore the scaling of the dose distributions in media other than water became dependent on the photon energy. The extended electron-scaling method was shown to estimate the photon range to within 1% for all materials studied and for energies from 100 keV to 20 MeV, allowing it to be used to scale dose distributions to media other than water and generated by clinical radiotherapy photon beams with accelerator energies from 4 to 20 MV.