Kilovoltage beam Monte Carlo dose calculations in submillimeter voxels for small animal radiotherapy

Absorbed dose calculation for a realistic CT-derived mouse phantom irradiated with a standard Cs-137 cell irradiator using a Monte Carlo method

Computed tomography (CT) derived Monte Carlo (MC) phantoms allow dose determination within small animal models that is not feasible with in-vivo dosimetry. The aim of this study was to develop a CT-derived MC phantom generated from a mouse with a xenograft tumour that could then be used to calculate both the dose heterogeneity in the tumour volume and out of field scattered dose for pre-clinical small animal irradiation experiments. A BEAMnrc Monte-Carlo model has been built of our irradiation system that comprises a lead collimator with a 1 cm diameter aperture fitted to a Cs-137 gamma irradiator. The MC model of the irradiation system was validated by comparing the calculated dose results with dosimetric film measurement in a polymethyl methacrylate (PMMA) phantom using a 1D gamma-index analysis. Dose distributions in the MC mouse phantom were calculated and visualized on the CT-image data. Dose volume histograms (DVHs) were generated for the tumour and organs at risk (OARs). The effect of the xenographic tumour volume on the scattered out of field dose was also investigated. The defined gamma index analysis criteria were met, indicating that our MC simulation is a valid model for MC mouse phantom dose calculations. MC dose calculations showed a maximum out of field dose to the mouse of 7% of Dmax. Absorbed dose to the tumour varies in the range 60%-100% of Dmax. DVH analysis demonstrated that tumour received an inhomogeneous dose of 12 Gy-20 Gy (for 20 Gy prescribed dose) while out of field doses to all OARs were minimized (1.29 Gy-1.38 Gy). Variation of the xenographic tumour volume exhibited no significant effect on the out of field scattered dose to OARs. The CT derived MC mouse model presented here is a useful tool for tumour dose verifications as well as investigating the doses to normal tissue (in out of field) for preclinical radiobiological research.

Monte Carlo simulations of EBT3 film dose deposition for percentage depth dose (PDD) curve evaluation

Purpose

To use Monte Carlo (MC) calculations to evaluate the effects of Gafchromic EBT3 film orientation on percentage depth dose (PDD) curves.

Methods

Dose deposition in films placed in a water phantom, and oriented either parallel or perpendicular with respect to beam axis, were simulated with MC and compared to PDDs scored in a homogenous water phantom. The effects of introducing 0.01–1.00 mm air gaps on each side of the film as well as a small 1°‐3° tilt for film placed in parallel orientation were studied. PDDs scored based on two published EBT3 film compositions were compared. Three photon beam energies of 120 kVp, 220 kVp, and 6 MV and three field sizes between 1 × 1 and 5 × 5 cm² were considered. Experimental PDDs for a 6‐MV 3 × 3 cm² beam were acquired.

Results

PDD curves for films in perpendicular orientation more closely agreed to water PDDs than films placed in parallel orientation. The maximum difference between film and water PDD for films in parallel orientation was −12.9% for the 220 kVp beam. For the perpendicular film orientation, the maximum difference decreased to 5.7% for the 120 kVp beam. The inclusion of an air gap had the largest effect on the 6‐MV 1 × 1 cm² beam, for which the dose in the buildup region was underestimated by 21.2% compared to the simulation with no air gap. A 2° film tilt decreased the difference between the parallel film and homogeneous water phantom PDDs from −5.0% to −0.5% for the 6 MV 3 × 3 cm² beam. The “newer” EBT3 film composition resulted in larger PDD discrepancies than the previous composition. Experimental film data qualitatively agreed with MC simulations.

Conclusions

PDD measurements with films should either be performed with film in perpendicular orientation to the beam axis or in parallel orientation with a ~ 2º tilt and no air gaps.

A Multimodality Image Guided Precision Radiation Research Platform: Integrating X-ray, Bioluminescence, and Fluorescence Tomography With Radiation Therapy

PURPOSE

Small animal irradiation is crucial to the investigation of radiobiological mechanisms. The paradigm of clinical radiation therapy is trending toward high-precision, stereotactic treatment. However, translating this scheme to small animal irradiation is challenging owing to the lack of high-quality image guidance. To overcome this obstacle, we developed a multimodality image guided precision radiation platform.

METHODS AND MATERIALS

The platform consists of 4 modules: x-ray computed tomography (CT), bioluminescence tomography (BLT), fluorescence molecular tomography (FMT), and radiation therapy. CT provides animal anatomy and material density for radiation dose calculation, as well as body contour for BLT and FMT reconstruction. BLT and FMT provide tumor localization to guide radiation beams and molecular activity to evaluate treatment outcome. Furthermore, we developed a Monte Carlo-based treatment planning system (TPS) for 3-dimensional dose calculation, calibrated it using radiochromic films sandwiched in a water-equivalent phantom, and validated it using in vivo dosimeters surgically implanted into euthanized mice (n = 4). Finally, we performed image guided irradiation on mice bearing orthotopic breast and prostate tumors and confirmed radiation delivery using γH2AX histology.

RESULTS

The Monte Carlo-based TPS was successfully calibrated by benchmarking simulation dose against film measurement. For in vivo dosimetry measured in the euthanized mice, the average difference between the TPS calculated dose and measured dose was 3.86% ± 1.12%. Following the TPS-generated treatment plan, we successfully delivered 20 Gy dose to an animal bearing an orthotopic prostate tumor using 4 BLT-guided radiation beams and 5 Gy dose to an animal bearing an orthotopic breast tumor using a single FMT-guided radiation beam. γH2AX histology presented significantly more DNA damage in irradiated tumors and thus validated the dose delivery accuracy.

CONCLUSIONS

Combined with Monte Carlo TPS, this multimodality CT/BLT/FMT image guided small animal radiation platform can specifically localize tumors, accurately calculate dose distribution, precisely guide radiation delivery, and molecularly evaluate treatment response. It provides an advanced toolset for radiobiology and translational cancer research.

Evaluation of a Novel Liquid Fiducial Marker, BioXmark®, for Small Animal Image-Guided Radiotherapy Applications

BioXmark^® (Nanovi A/S, Denmark) is a novel fiducial marker based on a liquid, iodine-based and non-metallic formulation. BioXmark^® has been clinically validated and reverse translated to preclinical models to improve cone-beam CT (CBCT) target delineation in small animal image-guided radiotherapy (SAIGRT). However, in phantom image analysis and in vivo evaluation of radiobiological response after the injection of BioXmark^® are yet to be reported. In phantom measurements were performed to compare CBCT imaging artefacts with solid fiducials and determine optimum imaging parameters for BioXmark^®. In vivo stability of BioXmark^® was assessed over a 5-month period, and the impact of BioXmark^® on in vivo tumour response from single-fraction and fractionated X-ray exposures was investigated in a subcutaneous syngeneic tumour model. BioXmark^® was stable, well tolerated and detectable on CBCT at volumes ≤10 µL. Our data showed imaging artefacts reduced by up to 84% and 89% compared to polymer and gold fiducial markers, respectively. BioXmark^® was shown to have no significant impact on tumour growth in control animals, but changes were observed in irradiated animals injected with BioXmark^® due to alterations in dose calculations induced by the sharp contrast enhancement. BioXmark^® is superior to solid fiducials with reduced imaging artefacts on CBCT. With minimal impact on the tumour growth delay, BioXmark^® can be implemented in SAIGRT to improve target delineation and reduce set-up errors.

Impact of backscatter material thickness on the depth dose of orthovoltage irradiators for radiobiology research

The orthovoltage x-ray energy frequently used in radiation research is prone to dosimetry errors due to insufficient backscatter conditions. In many radiobiology studies, especially for cell irradiations, precise dose calculation algorithms such as Convolution-Superposition or Monte Carlo are impractical and as such, less accurate hand calculation methods are used for dose estimation. These dose estimation methods typically assume full backscatter conditions. The purpose of this study is to demonstrate the magnitude of the dose error that results from insufficient backscatter, and to provide lookup tables to account this issue. The beam spectra of several widely used commercial systems (XRAD-225, XRAD-320, SARRP) were used in Monte Carlo (MC) simulations on a series of phantom setups to investigate the impact of varying backscatter conditions on dosimetry. The depth dose curves for different field sizes, water phantom thicknesses and beam qualities were generated. In addition, depth dependent backscatter factors for different field sizes and different beam qualities were calculated. It is demonstrated that as much as a 50% dose difference exists for different backscatter conditions at the beam qualities studied. The choice of cell dish size as well as other changes in the experiment setup can have more than 10% impact on the dose. The impact of backscatter is reduced with a decrease in field size. Further, the thickness needed to provide full backscatter can be approximated as being equal to the field size. It is imperative to ensure full backscatter conditions during system and dosimeter calibration, or to use the look-up table provided in this study.

Online virtual isocenter based radiation field targeting for high performance small animal microirradiation

Advances in precision microirradiators for small animal radiation oncology studies have provided the framework for novel translational radiobiological studies. Such systems target radiation fields at the scale required for small animal investigations, typically through a combination of on-board computed tomography image guidance and fixed, interchangeable collimators. Robust targeting accuracy of these radiation fields remains challenging, particularly at the millimetre scale field sizes achievable by the majority of microirradiators. Consistent and reproducible targeting accuracy is further hindered as collimators are removed and inserted during a typical experimental workflow. This investigation quantified this targeting uncertainty and developed an online method based on a virtual treatment isocenter to actively ensure high performance targeting accuracy for all radiation field sizes. The results indicated that the two-dimensional field placement uncertainty was as high as 1.16 mm at isocenter, with simulations suggesting this error could be reduced to 0.20 mm using the online correction method. End-to-end targeting analysis of a ball bearing target on radiochromic film sections showed an improved targeting accuracy with the three-dimensional vector targeting error across six different collimators reduced from [Formula: see text] mm (mean  ±  SD) to [Formula: see text] mm for an isotropic imaging voxel size of 0.1 mm.

Adaptation, Commissioning, and Evaluation of a 3D Treatment Planning System for High-Resolution Small-Animal Irradiation

Although spatially precise systems are now available for small-animal irradiations, there are currently limited software tools available for treatment planning for such irradiations. We report on the adaptation, commissioning, and evaluation of a 3-dimensional treatment planning system for use with a small-animal irradiation system. The 225-kV X-ray beam of the X-RAD 225Cx microirradiator (Precision X-Ray) was commissioned using both ion-chamber and radiochromic film for 10 different collimators ranging in field size from 1 mm in diameter to 40 × 40 mm(2) A clinical 3-dimensional treatment planning system (Metropolis) developed at our institution was adapted to small-animal irradiation by making it compatible with the dimensions of mice and rats, modeling the microirradiator beam orientations and collimators, and incorporating the measured beam data for dose calculation. Dose calculations in Metropolis were verified by comparison with measurements in phantoms. Treatment plans for irradiation of a tumor-bearing mouse were generated with both the Metropolis and the vendor-supplied software. The calculated beam-on times and the plan evaluation tools were compared. The dose rate at the central axis ranges from 74 to 365 cGy/min depending on the collimator size. Doses calculated with Metropolis agreed with phantom measurements within 3% for all collimators. The beam-on times calculated by Metropolis and the vendor-supplied software agreed within 1% at the isocenter. The modified 3-dimensional treatment planning system provides better visualization of the relationship between the X-ray beams and the small-animal anatomy as well as more complete dosimetric information on target tissues and organs at risk. It thereby enhances the potential of image-guided microirradiator systems for evaluation of dose-response relationships and for preclinical experimentation generally.

A novel technique for image-guided local heart irradiation in the rat

In radiotherapy treatment of thoracic, breast and chest wall tumors, the heart may be included (partially or fully) in the radiation field. As a result, patients may develop radiation-induced heart disease (RIHD) several years after exposure to radiation. There are few methods available to prevent or reverse RIHD and the biological mechanisms remain poorly understood. In order to further study the effects of radiation on the heart, we developed a model of local heart irradiation in rats using an image-guided small animal conformal radiation therapy device (SACRTD) developed at our institution. First, Monte Carlo based simulations were used to design an appropriate collimator. EBT-2 films were used to measure relative dosimetry, and the absolute dose rate at the isocenter was measured using the AAPM protocol TG-61. The hearts of adult male Sprague-Dawley rats were irradiated with a total dose of 21 Gy. For this purpose, rats were anesthetized with isoflurane and placed in a custom-made vertical rat holder. Each heart was irradiated with a 3-beam technique (one AP field and 2 lateral fields), with each beam delivering 7 Gy. For each field, the heart was visualized with a digital flat panel X-ray imager and placed at the isocenter of the 1.8 cm diameter beam. In biological analysis of radiation exposure, immunohistochemistry showed γH2Ax foci and nitrotyrosine throughout the irradiated hearts but not in the lungs. Long-term follow-up of animals revealed histopathological manifestations of RIHD, including myocardial degeneration and fibrosis. The results demonstrate that the rat heart irradiation technique using the SACRTD was successful and that surrounding untargeted tissues were spared, making this approach a powerful tool for in vivo radiobiological studies of RIHD. Functional and structural changes in the rat heart after local irradiation are ongoing.

Modality comparison for small animal radiotherapy: a simulation study

PURPOSE

Small animal radiation therapy has advanced significantly in recent years. Whereas in the past dose was delivered using a single beam and a lead shield for sparing of healthy tissue, conformal doses can be now delivered using more complex dedicated small animal radiotherapy systems with image guidance. The goal of this paper is to investigate dose distributions for three small animal radiation treatment modalities.

METHODS

This paper presents a comparison of dose distributions generated by the three approaches-a single-field irradiator with a 200 kV beam and no image guidance, a small animal image-guided conformal system based on a modified microCT scanner with a 120 kV beam developed at Stanford University, and a dedicated conformal system, SARRP, using a 220 kV beam developed at Johns Hopkins University. The authors present a comparison of treatment plans for the three modalities using two cases: a mouse with a subcutaneous tumor and a mouse with a spontaneous lung tumor. A 5 Gy target dose was calculated using the EGSnrc Monte Carlo codes.

RESULTS

All treatment modalities generated similar dose distributions for the subcutaneous tumor case, with the highest mean dose to the ipsilateral lung and bones in the single-field plan (0.4 and 0.4 Gy) compared to the microCT (0.1 and 0.2 Gy) and SARRP (0.1 and 0.3 Gy) plans. The lung case demonstrated that due to the nine-beam arrangements in the conformal plans, the mean doses to the ipsilateral lung, spinal cord, and bones were significantly lower in the microCT plan (2.0, 0.4, and 1.9 Gy) and the SARRP plan (1.5, 0.5, and 1.8 Gy) than in single-field irradiator plan (4.5, 3.8, and 3.3 Gy). Similarly, the mean doses to the contralateral lung and the heart were lowest in the microCT plan (1.5 and 2.0 Gy), followed by the SARRP plan (1.7 and 2.2 Gy), and they were highest in the single-field plan (2.5 and 2.4 Gy). For both cases, dose uniformity was greatest in the single-field irradiator plan followed by the SARRP plan due to the sensitivity of the lower energy microCT beam to target heterogeneities and image noise.

CONCLUSIONS

The two treatment planning examples demonstrate that modern small animal radiotherapy techniques employing image guidance, variable collimation, and multiple beam angles deliver superior dose distributions to small animal tumors as compared to conventional treatments using a single-field irradiator. For deep-seated mouse tumors, however, higher-energy conformal radiotherapy could result in higher doses to critical organs compared to lower-energy conformal radiotherapy. Treatment planning optimization for small animal radiotherapy should therefore be developed to take full advantage of the novel conformal systems.

Toward an organ based dose prescription method for the improved accuracy of murine dose in orthovoltage x-ray irradiators

PURPOSE

Accurate dosimetry is essential when irradiating mice to ensure that functional and molecular endpoints are well understood for the radiation dose delivered. Conventional methods of prescribing dose in mice involve the use of a single dose rate measurement and assume a uniform average dose throughout all organs of the entire mouse. Here, the authors report the individual average organ dose values for the irradiation of a 12, 23, and 33 g mouse on a 320 kVp x-ray irradiator and calculate the resulting error from using conventional dose prescription methods.

METHODS

Organ doses were simulated in the Geant4 application for tomographic emission toolkit using the MOBY mouse whole-body phantom. Dosimetry was performed for three beams utilizing filters A (1.65 mm Al), B (2.0 mm Al), and C (0.1 mm Cu + 2.5 mm Al), respectively. In addition, simulated x-ray spectra were validated with physical half-value layer measurements.

RESULTS

Average doses in soft-tissue organs were found to vary by as much as 23%-32% depending on the filter. Compared to filters A and B, filter C provided the hardest beam and had the lowest variation in soft-tissue average organ doses across all mouse sizes, with a difference of 23% for the median mouse size of 23 g.

CONCLUSIONS

This work suggests a new dose prescription method in small animal dosimetry: it presents a departure from the conventional approach of assigninga single dose value for irradiation of mice to a more comprehensive approach of characterizing individual organ doses to minimize the error and uncertainty. In human radiation therapy, clinical treatment planning establishes the target dose as well as the dose distribution, however, this has generally not been done in small animal research. These results suggest that organ dose errors will be minimized by calibrating the dose rates for all filters, and using different dose rates for different organs.

High throughput film dosimetry in homogeneous and heterogeneous media for a small animal irradiator

PURPOSE

We have established a high-throughput Gafchromic film dosimetry protocol for narrow kilovoltage beams in homogeneous and heterogeneous media for small-animal radiotherapy applications. The kV beam characterization is based on extensive Gafchromic film dosimetry data acquired in homogeneous and heterogeneous media. An empirical model is used for parameterization of depth and off-axis dependence of measured data.

METHODS

We have modified previously published methods of film dosimetry to suit the specific tasks of the study. Unlike film protocols used in previous studies, our protocol employs simultaneous multi-channel scanning and analysis of up to nine Gafchromic films per scan. A scanner and background correction were implemented to improve accuracy of the measurements. Measurements were taken in homogeneous and inhomogeneous phantoms at 220 kVp and a field size of 5 × 5 mm(2). The results were compared against Monte Carlo simulations.

RESULTS

Dose differences caused by variations in background signal were effectively removed by the corrections applied. Measurements in homogeneous phantoms were used to empirically characterize beam data in homogeneous and heterogeneous media. Film measurements in inhomogeneous phantoms and their empirical parameterization differed by about 2%-3%. The model differed from MC by about 1% (water, lung) to 7% (bone). Good agreement was found for measured and modelled off-axis ratios.

CONCLUSIONS

EBT2 films are a valuable tool for characterization of narrow kV beams, though care must be taken to eliminate disturbances caused by varying background signals. The usefulness of the empirical beam model in interpretation and parameterization of film data was demonstrated.

Dependences of mucosal dose on photon beams in head-and-neck intensity-modulated radiation therapy: a Monte Carlo study

Dependences of mucosal dose in the oral or nasal cavity on the beam energy, beam angle, multibeam configuration, and mucosal thickness were studied for small photon fields using Monte Carlo simulations (EGSnrc-based code), which were validated by measurements. Cylindrical mucosa phantoms (mucosal thickness = 1, 2, and 3 mm) with and without the bone and air inhomogeneities were irradiated by the 6- and 18-MV photon beams (field size = 1 × 1 cm(2)) with gantry angles equal to 0°, 90°, and 180°, and multibeam configurations using 2, 4, and 8 photon beams in different orientations around the phantom. Doses along the central beam axis in the mucosal tissue were calculated. The mucosal surface doses were found to decrease slightly (1% for the 6-MV photon beam and 3% for the 18-MV beam) with an increase of mucosal thickness from 1-3 mm, when the beam angle is 0°. The variation of mucosal surface dose with its thickness became insignificant when the beam angle was changed to 180°, but the dose at the bone-mucosa interface was found to increase (28% for the 6-MV photon beam and 20% for the 18-MV beam) with the mucosal thickness. For different multibeam configurations, the dependence of mucosal dose on its thickness became insignificant when the number of photon beams around the mucosal tissue was increased. The mucosal dose with bone was varied with the beam energy, beam angle, multibeam configuration and mucosal thickness for a small segmental photon field. These dosimetric variations are important to consider improving the treatment strategy, so the mucosal complications in head-and-neck intensity-modulated radiation therapy can be minimized.

The importance of tissue segmentation for dose calculations for kilovoltage radiation therapy

PURPOSE

The aim of this work was to evaluate the effect of tissue segmentation on the accuracy of Monte Carlo (MC) dose calculations for kilovoltage radiation therapy, which are commonly used in preclinical radiotherapy studies and are also being revisited as a clinical treatment modality. The feasibility of tissue segmentation routinely done on the basis of differences in tissue mass densities was studied and a new segmentation scheme based on differences in effective atomic numbers was developed.

METHODS

MC dose calculations in a cylindrical mouse phantom with small cylindrical inhomogeneities consisting of 34 ICRU-44 tissues were performed using the EGSnrc/BEAMnrc and DOSXYZnrc codes. The dose to tissue was calculated for five different kilovoltage beams currently used in small animal radiotherapy: a microCT 120 kV beam, two 225 kV beams filtered with either 4 mm of Al or 0.5 mm of Cu, a heavily filtered 320 kV beam, and a 192Ir beam. The mean doses to the 34 ICRU-44 tissues as a function of tissue mass density and effective atomic number and beam energy were studied. A treatment plan for an orthotopic lung tumor model was created, and the dose distribution was calculated for three tissue segmentation schemes using 4, 8, and 39 tissue bins to assess the significance of the simulation results for kilovoltage radiotherapy.

RESULTS

In our model, incorrect assignment of adipose tissue to muscle caused dose calculation differences of 27%, 13%, and 7% for the 120 kV beam and the 225 kV beams filtered with 4 mm Al and 0.5 mm Cu, respectively. For the heavily filtered 320 kV beam and a 192Ir source, potential dose calculation differences due to tissue mis-assignment were below 4%. There was no clear relationship between the dose to tissue and its mass density for x-ray beams generated by tube potentials equal or less than 225 kV. A second order polynomial fit approximated well the absorbed dose to tissue as a function of effective atomic number for these beams. In the mouse study, the 120 kV beam dose to bone was overestimated by 100% and underestimated by 10% for the 4 and 8-tissue segmentation schemes compared to the 39-tissue segmentation scheme, respectively. Dose to adipose tissue was overestimated by 30% and underestimated by 10%, respectively. In general, organ at risk (OAR) doses were overestimated in the 4-tissue and the 8-tissue segmentation schemes compared to the 39-tissue segmentation.

CONCLUSIONS

Tissue segmentation was shown to be a key parameter for dose calculations with kilovoltage beams used in small animal radiotherapy when an x-ray tube with a potential < or = 225 kV is used as a source. A new tissue segmentation scheme with 39 tissues based on effective number differences derived from mass density differences has been implemented.

Dosimetric characterization of an image-guided stereotactic small animal irradiator

Small animal irradiation provides an important tool used by preclinical studies to assess and optimize new treatment strategies such as stereotactic ablative radiotherapy. Characterization of radiation beams that are clinically and geometrically scaled for the small animal model is uniquely challenging for orthovoltage energies and minute field sizes. The irradiator employs a commercial x-ray device (XRAD 320, Precision x-ray, Inc.) with a custom collimation system to produce 1-10 mm diameter beams and a 50 mm reference beam. Absolute calibrations were performed using the AAPM TG-61 methodology. Beam's half-value layer (HVL) and timer error were measured with an ionization chamber. Percent depth dose (PDD), output factors (OFs) and off-axis ratios were measured using radiochromic film, a diode and a pinpoint ionization chamber at 19.76 and 24.76 cm source-to-surface distance (SSD). PDD measurements were also compared with Monte Carlo (MC) simulations. In-air and in-water absolute calibrations for the reference 50 mm diameter collimator at 19.76 cm SSD were measured as 20.96 and 20.79 Gy min(-1), respectively, agreeing within 0.8%. The HVL at 250 kVp and 15 mAs was measured to be 0.45 mm Cu. The reference field PDD MC simulation results agree with measured data within 3.5%. PDD data demonstrate typical increased penetration with increasing field size and SSD. For collimators larger than 5 mm in diameter, OFs measured using film, an ion chamber and a diode were within 3% agreement.

Dosimetric variation due to the photon beam energy in the small-animal irradiation: a Monte Carlo study

PURPOSE

The impact of photon beam energy and tissue heterogeneities on dose distributions and dosimetric characteristics such as point dose, mean dose, and maximum dose was investigated in the context of small-animal irradiation using Monte Carlo simulations based on the EGSnrc code.

METHODS

Three Monte Carlo mouse phantoms, namely, heterogeneous, homogeneous, and bone homogeneous were generated based on the same mouse computed tomography image set. These phantoms were generated by overriding the tissue type of none of the voxels (heterogeneous), all voxels (homogeneous), and only the bone voxels (bone homogeneous) to that of soft tissue. Phase space files of the 100 and 225 kVp photon beams based on a small-animal irradiator (XRad225Cx, Precision X-Ray Inc., North Branford, CT) were generated using BEAMnrc. A 360 degrees photon arc was simulated and three-dimensional (3D) dose calculations were carried out using the DOSXYZnrc code through DOSCTP in the above three phantoms. For comparison, the 3D dose distributions, dose profiles, mean, maximum, and point doses at different locations such as the isocenter, lung, rib, and spine were determined in the three phantoms.

RESULTS

The dose gradient resulting from the 225 kVp arc was found to be steeper than for the 100 kVp arc. The mean dose was found to be 1.29 and 1.14 times higher for the heterogeneous phantom when compared to the mean dose in the homogeneous phantom using the 100 and 225 kVp photon arcs, respectively. The bone doses (rib and spine) in the heterogeneous mouse phantom were about five (100 kVp) and three (225 kVp) times higher when compared to the homogeneous phantom. However, the lung dose did not vary significantly between the heterogeneous, homogeneous, and bone homogeneous phantom for the 225 kVp compared to the 100 kVp photon beams.

CONCLUSIONS

A significant bone dose enhancement was found when the 100 and 225 kVp photon beams were used in small-animal irradiation. This dosimetric effect, due to the presence of the bone heterogeneity, was more significant than that due to the lung heterogeneity. Hence, for kV photon energies of the range used in small-animal irradiation, the increase of the mean and bone dose due to the photoelectric effect could be a dosimetric concern.

Investigation of the effects of treatment planning variables in small animal radiotherapy dose distributions

PURPOSE

Methods used for small animal radiation treatment have yet to achieve the same dose targeting as in clinical radiation therapy. Toward understanding how to better plan small animal radiation using a system recently developed for this purpose, the authors characterized dose distributions produced from conformal radiotherapy of small animals in a microCT scanner equipped with a variable-aperture collimator.

METHODS

Dose distributions delivered to a cylindrical solid water phantom were simulated using a Monte Carlo algorithm. Phase-space files for 120 kVp x-ray beams and collimator widths of 1-10 mm at isocenter were generated using BEAMnrc software, and dose distributions for evenly spaced beams numbered from 5 to 80 were generated in DOSXYZnrc for a variety of targets, including centered spherical targets in a range of sizes, spherical targets offset from centered by various distances, and various ellipsoidal targets. Dose distributions were analyzed using dose volume histograms. The dose delivered to a mouse bearing a spontaneous lung tumor was also simulated, and dose volume histograms were generated for the tumor, heart, left lung, right lung, and spinal cord.

RESULTS

Results indicated that for centered, symmetric targets, the number of beams required to achieve a smooth dose volume histogram decreased with increased target size. Dose distributions for noncentered, symmetric targets did not exhibit any significant loss of conformality with increasing offset from the phantom center, indicating sufficient beam penetration through the phantom for targeting superficial targets from all angles. Even with variable collimator widths, targeting of asymmetric targets was found to have less conformality than that of spherical targets. Irradiation of a mouse lung tumor with multiple beam widths was found to effectively deliver dose to the tumor volume while minimizing dose to other critical structures.

CONCLUSIONS

Overall, this method of generating and analyzing dose distributions provides a quantitative method for developing practical guidelines for small animal radiotherapy treatment planning. Future work should address methods to improve conformality in asymmetric targets.