Effect of age-dependent bone electron density on the calculated dose distribution from kilovoltage and megavoltage photon and electron radiotherapy in paediatric MRI-only treatment planning.

OBJECTIVE: MRI-only treatment planning (TP) can be advantageous in paediatric radiotherapy. However, electron density extraction is necessary for dose calculation. Normally, after bone segmentation, a bulk density is assigned. However, the variation of bone bulk density in patients makes the creation of pseudo CTs challenging. This study aims to assess the effects of bone density variations in children on radiation attenuation and dose calculation for MRI-only TP. METHODS: Bone contents of <15-year-old children were calculated, and substituted in the Oak Ridge National Laboratory paediatric phantoms. The percentage depth dose and beam profile of 150 kVp and 6 MV photon and 6 MeV electron beams were then calculated using Xcom, MCNPX (Monte Carlo N-particle version X) and ORLN phantoms. RESULTS: Using 150 kVp X-rays, the difference in attenuation coefficient was almost 5% between an 11-year-old child and a newborn, and ~8% between an adult and a newborn. With megavoltage radiation, the differences were smaller but still important. For an 18 MV photon beam, the difference of radiation attenuation between an 11-year-old child and a newborn was 4% and ~7.4% between an adult and a newborn. For 6 MeV electrons, dose differences were observed up to the 2 cm depth. The percentage depth dose difference between 1 and 10-year-olds was 18.5%, and between 10 and 15-year-olds was 24%. CONCLUSION: The results suggest that for MRI-only TP of photon- or electron-beam radiotherapy, the bone densities of each age group should be defined separately for accurate dose calculation. Advances in knowledge: This study highlights the need for more age-specific determination of bone electron density for accurate dose calculations in paediatric MRI-only radiotherapy TP.

Revision of orthovoltage chest wall treatment using Monte Carlo simulations.

PURPOSE: Given the high local control rates observed in breast cancer patients undergoing chest wall irradiation by kilovoltage x-rays, we aimed to revisit this treatment modality by accurate calculation of dose distributions using Monte Carlo simulation. METHODS AND MATERIAL: The machine components were simulated using the MCNPX code. This model was used to assess the dose distribution of chest wall kilovoltage treatment in different chest wall thicknesses and larger contour or fat patients in standard and mid sternum treatment plans. Assessments were performed at 50 and 100 cm focus surface distance (FSD) and different irradiation angles. In order to evaluate different plans, indices like homogeneity index, conformity index, the average dose of heart, lung, left anterior descending artery (LAD) and percentage target coverage (PTC) were used. Finally, the results were compared with the indices provided by electron therapy which is a more routine treatment of chest wall. RESULT: These indices in a medium chest wall thickness in standard treatment plan at 50 cm FSD and 15 degrees tube angle was as follows: homogeneity index 2.57, conformity index 7.31, average target dose 27.43 Gy, average dose of heart, lung and LAD, 1.03, 2.08 and 1.60 Gy respectively and PTC 11.19%. Assessments revealed that dose homogeneity in planning target volume (PTV) and conformity between the high dose region and PTV was poor. To improve the treatment indices, the reference point was transferred from the chest wall skin surface to the center of PTV. The indices changed as follows: conformity index 7.31, average target dose 60.19 Gy, the average dose of heart, lung and LAD, 3.57, 6.38 and 5.05 Gy respectively and PTC 55.24%. Coverage index of electron therapy was 89% while it was 22.74% in the old orthovoltage method and also the average dose of the target was about 50 Gy but in the given method it was almost 30 Gy. CONCLUSION: The results of the treatment study show that the optimized standard and mid sternum treatment for different chest wall thicknesses is with 50 cm FSD and zero (vertical) tube angle, while in large contour patients, it is with 100 cm FSD and zero tube angle. Finally, chest wall kilovoltage and electron therapies were compared, which revealed that electron therapy produces a better dose distribution than kilovoltage therapy.

Fast and accurate Monte Carlo modeling of a kilovoltage X-ray therapy unit using a photon-source approximation for treatment planning in complex media.

To accurately recompute dose distributions in chest-wall radiotherapy with 120 kVp kilovoltage X-rays, an MCNP4C Monte Carlo model is presented using a fast method that obviates the need to fully model the tube components. To validate the model, half-value layer (HVL), percentage depth doses (PDDs) and beam profiles were measured. Dose measurements were performed for a more complex situation using thermoluminescence dosimeters (TLDs) placed within a Rando phantom. The measured and computed first and second HVLs were 3.8, 10.3 mm Al and 3.8, 10.6 mm Al, respectively. The differences between measured and calculated PDDs and beam profiles in water were within 2 mm/2% for all data points. In the Rando phantom, differences for majority of data points were within 2%. The proposed model offered an approximately 9500-fold reduced run time compared to the conventional full simulation. The acceptable agreement, based on international criteria, between the simulations and the measurements validates the accuracy of the model for its use in treatment planning and radiobiological modeling studies of superficial therapies including chest-wall irradiation using kilovoltage beam.

Monte Carlo N Particle code - Dose distribution of clinical electron beams in inhomogeneous phantoms.

Electron dose distributions calculated using the currently available analytical methods can be associated with large uncertainties. The Monte Carlo method is the most accurate method for dose calculation in electron beams. Most of the clinical electron beam simulation studies have been performed using non- MCNP [Monte Carlo N Particle] codes. Given the differences between Monte Carlo codes, this work aims to evaluate the accuracy of MCNP4C-simulated electron dose distributions in a homogenous phantom and around inhomogeneities. Different types of phantoms ranging in complexity were used; namely, a homogeneous water phantom and phantoms made of polymethyl methacrylate slabs containing different-sized, low- and high-density inserts of heterogeneous materials. Electron beams with 8 and 15 MeV nominal energy generated by an Elekta Synergy linear accelerator were investigated. Measurements were performed for a 10 cm × 10 cm applicator at a source-to-surface distance of 100 cm. Individual parts of the beam-defining system were introduced into the simulation one at a time in order to show their effect on depth doses. In contrast to the first scattering foil, the secondary scattering foil, X and Y jaws and applicator provide up to 5% of the dose. A 2%/2 mm agreement between MCNP and measurements was found in the homogenous phantom, and in the presence of heterogeneities in the range of 1-3%, being generally within 2% of the measurements for both energies in a "complex" phantom. A full-component simulation is necessary in order to obtain a realistic model of the beam. The MCNP4C results agree well with the measured electron dose distributions.

EchoSeed Model 6733 Iodine-125 brachytherapy source: improved dosimetric characterization using the MCNP5 Monte Carlo code.

This study primarily aimed to obtain the dosimetric characteristics of the Model 6733 (125)I seed (EchoSeed) with improved precision and accuracy using a more up-to-date Monte-Carlo code and data (MCNP5) compared to previously published results, including an uncertainty analysis. Its secondary aim was to compare the results obtained using the MCNP5, MCNP4c2, and PTRAN codes for simulation of this low-energy photon-emitting source. The EchoSeed geometry and chemical compositions together with a published (125)I spectrum were used to perform dosimetric characterization of this source as per the updated AAPM TG-43 protocol. These simulations were performed in liquid water material in order to obtain the clinically applicable dosimetric parameters for this source model. Dose rate constants in liquid water, derived from MCNP4c2 and MCNP5 simulations, were found to be 0.993 cGyh(-1) U(-1) (±1.73%) and 0.965 cGyh(-1) U(-1) (±1.68%), respectively. Overall, the MCNP5 derived radial dose and 2D anisotropy functions results were generally closer to the measured data (within ±4%) than MCNP4c and the published data for PTRAN code (Version 7.43), while the opposite was seen for dose rate constant. The generally improved MCNP5 Monte Carlo simulation may be attributed to a more recent and accurate cross-section library. However, some of the data points in the results obtained from the above-mentioned Monte Carlo codes showed no statistically significant differences. Derived dosimetric characteristics in liquid water are provided for clinical applications of this source model.

Independent evaluation of an in-house brachytherapy treatment planning system using simulation, measurement and calculation methods.

Accuracy of treatment planning systems may significantly influence the efficacy of brachytherapy. The purpose of this work is a detailed, varied and independent evaluation of an in-house brachytherapy treatment planning software called STPS. Operational accuracy of STPS was investigated. Geometric tests were performed to validate entry and reconstruction of positional information from scanned orthogonal films. MCNP4C Monte Carlo code and TLDs were used for simulation and experimental measurement, respectively. STPS data were also compared with those from a commercial planning system (Nucletron PLATO). Discrepancy values between MCNP and STPS data and also those of PLATO and STPS at Manchester system dose prescription points (AL and AR) of tandem and ovoid configurations were 2.5% ± 0.5% and 5.4% ± 0.4%, respectively. Similar results were achieved for other investigated configurations. Observed discrepancies between MCNP and STPS at the dose prescription point and at 1 cm from the tip of the vaginal applicator were 4.5% and 25.6% respectively, while the discrepancy between the STPS and PLATO data at those points was 2.3%. The software showed submillimeter accuracy in its geometrical reconstructions. In terms of calculation accuracy, similar to PLATO, as attenuation of the sources and applicator body is not considered, dose was overestimated at the tip of the applicator, but based on the available criteria, dose accuracy at most points were acceptable. Our results confirm STPS's geometrical and operational reliability, and show that its dose computation accuracy is comparable to an established commercial TPS using the same algorithm.

Impact of the vaginal applicator and dummy pellets on the dosimetry parameters of Cs-137 brachytherapy source.

In this study, dose rate distribution around a spherical 137Cs pellet source, from a low-dose-rate (LDR) Selectron remote afterloading system used in gynecological brachytherapy, has been determined using experimental and Monte Carlo simulation techniques. Monte Carlo simulations were performed using MCNP4C code, for a single pellet source in water medium and Plexiglas, and measurements were performed in Plexiglas phantom material using LiF TLD chips. Absolute dose rate distribution and the dosimetric parameters, such as dose rate constant, radial dose functions, and anisotropy functions, were obtained for a single pellet source. In order to investigate the effect of the applicator and surrounding pellets on dosimetric parameters of the source, the simulations were repeated for six different arrangements with a single active source and five non-active pellets inside central metallic tubing of a vaginal cylindrical applicator. In commercial treatment planning systems (TPS), the attenuation effects of the applicator and inactive spacers on total dose are neglected. The results indicate that this effect could lead to overestimation of the calculated F(r,θ), by up to 7% along the longitudinal axis of the applicator, especially beyond the applicator tip. According to the results obtained in this study, in a real situation in treatment of patients using cylindrical vaginal applicator and using several active pellets, there will be a large discrepancy between the result of superposition and Monte Carlo simulations.

Monte Carlo and Lambertian light guide models of the light output from scintillation crystals at megavoltage energies.

A new model of the light output from single-crystal scintillators in megavoltage energy x-ray beams has been developed, based on the concept of a Lambertian light guide model (LLG). This was evaluated in comparison with a Monte Carlo (MC) model of optical photon transport, previously developed and reported in the literature, which was used as a gold standard. The LLG model was developed to enable optimization of scintillator detector design. In both models the dose deposition and light propagation were decoupled, the scintillators were cuboids, split into a series of cells as a function of depth, with Lambertian side and entrance faces, and a specular exit face. The signal in a sensor placed 1 and 1000 mm beyond the exit face was calculated. Cesium iodide (CSI) crystals of 1.5 and 3 mm square cross section and 1, 5, and 10 mm depth were modeled. Both models were also used to determine detector signal and optical gain factor as a function of CsI scintillator thickness, from 2 to 10 mm. Results showed a variation in light output with position of dose deposition of a factor of up to approximately 5, for long, thin scintillators (such as 10 X 1.5 x 1.5 mm3). For short, fat scintillators (such as 1 X 3 X 3 mm3) the light output was more uniform with depth. MC and LLG generally agreed to within 5%. Results for a sensor distance of 1 mm showed an increase in light output the closer the light originates to the exit face, while a distance of 1000 mm showed a decrease in light output the closer the light originates to the exit face. For a sensor distance of 1 mm, the ratio of signal for a 10 mm scintillator to that for a 2 mm scintillator was 1.98, whereas for the 1000 mm distance the ratio was 3.00. The ratio of quantum efficiency (QE) between 10 and 2 mm thicknesses was 4.62. We conclude that these models may be used for detector optimization, with the light guide model suitable for parametric study.

Commissioning and quality assurance of the Pinnacle(3) radiotherapy treatment planning system for external beam photons.

The commissioning of a Pinnacle(3) treatment planning system is described. Four Elekta linear accelerators were commissioned for external beam photons. Measured data were used to derive parameter values for the Pinnacle(3) beam model by (1). fitting a Monte Carlo model of the accelerator head to measured data and then extracting the parameters for the Pinnacle(3) beam model, and by (2). using the auto-modelling facility within Pinnacle(3). Both of these methods yielded dose distributions in accord with published recommendations. A separate small-field beam model, customized for an in-house compact blocking system, was also created, which satisfied appropriate acceptance criteria for stereotactically guided conformal brain treatments. Inhomogeneous, oblique, asymmetrical and irregular fields were also assessed, with calculated and measured doses agreeing to within +/-3%. Dose-volume histogram calculation was found to be accurate to within +/-5% dose or volume for a grid size of 4 mm x 4 mm x 4 mm, with better accuracy being achieved for finer grids. Isocentric doses were compared between Pinnacle(3)'s collapsed cone convolution algorithm and the Bentley-Milan algorithm within the Target-2 treatment planning system. Dose differences were generally less than 3% in the dose prescribed, with larger values for breast plans, where the Pinnacle(3) algorithm calculated scatter more accurately. Pelvic and thoracic plans were also verified using an anthropomorphic phantom, with local dose differences between calculated and delivered dose of up to 8%, but mainly less than 3%, and with no systematic difference. Ionization chamber verifications using START and RT-01 trial procedures demonstrated differences between calculated and measured doses of less than 2%. Following satisfactory performance in the commissioning process, Pinnacle(3) has now been introduced into routine clinical use.

Optimization of the scintillation detector in a combined 3D megavoltage CT scanner and portal imager.

A parametric study is described leading to the optimization of a custom-made scintillation detector with a relatively high quantum efficiency (QE) for megavoltage photons and light output toward a remote lens. This detector allows low-dose portal imaging and continuous cone-beam megavoltage CT acquisition. The EGS4 Monte Carlo code was used to simulate the x-ray and electron transport in the detector. A Monte Carlo model of optical photon transport in a detector element was devised and used as well as various irradiation experiments on scintillators. Different detector materials and configurations were compared in terms of the optical photon irradiance on the lens from on- and off-axis detector elements and the practical constraints regarding detector construction and weight. Effects of scintillator material, detector element size, crystal coating type, and reflectivity, combinations of different coatings on detector faces, scintillator doping level, and crystal transparency were studied. With scintillator thickness adjusted to give an 18% x-ray QE at 6 MV, the light output of CsI(Tl) was at least eight times higher than ZnWO4, BGO and NE118 plastic. Further, CsI(Tl) showed the smallest decrease in QE going from 6 to 24 MV. The off-axis reduction in emittance from the periphery of the detector was relatively small with a slight dependence on the type and reflectivity of the coating and the crystal thickness for a fixed detector element cross section. Light output was more strongly dependent on the reflectivity of lambertian coatings than specular ones. For a fixed detector element cross section, optimum coating type depended on crystal thickness. Typical CsI(Tl) crystals showed a relatively small variation in light output with changes in optical attenuation length. The optimum detector element was found to be CsI(Tl) coated on five faces with TiO2-loaded epoxy resin offering about a ten-fold improvement in light output per incident photon compared to typical metal/phosphor screens.