VMC++ validation for photon beams in the energy range of 20-1000 keV

In high energy teletherapy, VMC++ is known to be a very accurate and efficient Monte Carlo (MC) code. In principle, the MC method is also a powerful dose calculation tool in other areas in radiation oncology, e.g., brachytherapy or orthovoltage radiotherapy. However, VMC++ is not validated for the low-energy range of such applications. This work aims in the validation of the VMC++ MC code for photon beams in the energy range between 20 and 1000 keV.

Comparison of monte carlo collimator transport methods for photon treatment planning in radiotherapy

The aim of this work was a Monte Carlo (MC) based investigation of the impact of different radiation transport methods in collimators of a linear accelerator on photon beam characteristics, dose distributions, and efficiency. Thereby it is investigated if it is possible to use different simplifications in the radiation transport for some clinical situations in order to save calculation time.

The Swiss IMRT dosimetry intercomparison using a thorax phantom

In 2008, a national intensity modulated radiation therapy (IMRT) dosimetry intercomparison was carried out for all 23 radiation oncology institutions in Switzerland. It was the aim to check the treatment chain focused on the planning, dose calculation, and irradiation process.

Susceptibility of Aethina tumida (Coleoptera: Nitidulidae) Larvae and Pupae to Entomopathogenic Nematodes

In this study, we evaluated the potential use of entomopathogenic nematodes as a control for the beetle Aethina tumida Murray (Coleoptera: Nitidulidae). In particular, we conducted 1) four screening bioassays to determine nematode (seven species, 10 total strains tested) and application level effects on A. tumida larvae and pupae, 2) a generational persistence bioassay to determine whether single inoculations with nematodes would control multiple generations of A. tumida larvae in treated soil, and 3) a field bioassay to determine whether the nematodes would remain efficacious in the field. In the screening bioassays, nematode efficacy varied significantly by tested nematode and the infective juvenile (IJ) level at which they were applied. Although nematode virulence was moderate in screening bioassays 1-3 (0 - 68% A. tumida mortality), A. tumida mortality approached higher levels in screening bioassay 4 (nearly 100% after 39 d) that suggest suitable applicability of some of the test nematodes as field controls for A. tumida. In the generational persistence bioassay, Steinernema Hobrave Cabanillas, Poinar & Raulston 7-12 strain and Heterorhabditis indica Poinar, Karunaka & David provided adequate A. tumida control for 19 wk after a single soil inoculation (76-94% mortality in A. tumida pupae). In the field bioassay, the same two nematode species also showed high virulence toward pupating A. tumida (88-100%) mortality. Our data suggest that nematode use may be an integral component of an integrated pest management scheme aimed at reducing A. tumida populations in bee colonies to tolerable levels.

Winter Losses of Honeybee Colonies (Hymenoptera: Apidae): The Role of Infestations With Aethina tumida (Coleoptera: Nitidulidae) and Varroa destructor (Parasitiformes: Varroidae)

Multiple infections of managed honeybee, Apis mellifera, colonies are inevitable due to the ubiquitous ectoparasitic mite Varroa destructor and might be an underlying cause of winter losses. Here we investigated the role of adult small hive beetles, Aethina tumida, alone and in combination with V. destructor for winter losses and for infections with the microsporidian endoparasite Nosema ceranae. We found no significant influence of A. tumida and V destructor alone or in combination on the numbers of N. ceranae spores. Likewise, A. tumida alone had no significant effects on winter losses, which is most likely due to the observed high winter mortality of the adult beetles. Therefore, our data suggest that A. tumida is unlikely to contribute to losses of overwintering honeybee colonies. However, high losses occurred in all groups highly infested with V. destructor, supporting the central role of the mite for colony losses.

Monte Carlo implementation, validation, and characterization of a 120 leaf MLC

Recently, the new high definition multileaf collimator (HD120 MLC) was commercialized by Varian Medical Systems providing high resolution in the center section of the treatment field. The aim of this work is to investigate the characteristics of the HD120 MLC using Monte Carlo (MC) methods.

A virtual source model for kilo-voltage cone beam CT: source characteristics and model validation

The purpose of this investigation was to study the source characteristics of a clinical kilo-voltage cone beam CT unit and to develop and validate a virtual source model that could be used for treatment planning purposes.

X-ray phase contrast imaging of metallized and non-metallized biological samples

Imaging of biological samples has been performed with a variety of techniques for example electromagnetic waves, electrons, neutrons, ultrasound and X-rays. Also conventional X-ray imaging represents the basis of medical diagnostic imaging, it remains of limited use in this application because it is based solely on the differential absorption of X-rays by tissues. Coherent and bright photon beams, such as those produced by third-generation synchrotron X-ray sources, provide further information on subtle X-ray phase changes at matter interfaces. This complements conventional X-ray absorption by edge enhancement phenomena. Thus, phase contrast imaging has the potential to improve the detection of structures on images by detecting those structures that are invisible with X-ray absorption imaging. Images of a weakly absorbing nylon fibre were recorded in in-line holography geometry using a high resolution low-noise CCD camera at the ESRF in Grenoble. The method was also applied to improve image contrast for images of biological tissues. This paper presents phase contrast microradiographs of vascular tree casts and images of a housefly. These reveal very fine structures, that remain invisible with conventional absorption contrast only.

Leaf transmission reduction using moving jaws for dynamic MLC IMRT

The aim of this work is to investigate to what extent it is possible to use the secondary collimator jaws to reduce the transmitted radiation through the multileaf collimator (MLC) during an intensity modulated radiation therapy (IMRT). A method is developed and introduced where the jaws follow the open window of the MLC dynamically (dJAW method). With the aid of three academic cases (Closed MLC, Sliding-gap, and Chair) and two clinical cases (prostate and head and neck) the feasibility of the dJAW method and the influence of this method on the applied dose distributions are investigated. For this purpose the treatment planning system Eclipse and the Research-Toolbox were used as well as measurements within a solid water phantom were performed. The transmitted radiation through the closed MLC leads to an inhomogeneous dose distribution. In this case, the measured dose within a plane perpendicular to the central axis differs up to 40% (referring to the maximum dose within this plane) for 6 and 15 MV. The calculated dose with Eclipse is clearly more homogeneous. For the Sliding-gap case this difference is still up to 9%. Among other things, these differences depend on the depth of the measurement within the solid water phantom and on the application method. In the Chair case, the dose in regions where no dose is desired is locally reduced by up to 50% using the dJAW method instead of the conventional method. The dose inside the chair-shaped region decreased up to 4% if the same number of monitor units (MU) as for the conventional method was applied. The undesired dose in the volume body minus the planning target volume in the clinical cases prostate and head and neck decreased up to 1.8% and 1.5%, while the number of the applied MU increased up to 3.1% and 2.8%, respectively. The new dJAW method has the potential to enhance the optimization of the conventional IMRT to a further step.

Dosimetric impact of geometric errors due to respiratory motion prediction on dynamic multileaf collimator-based four-dimensional radiation delivery

The synchronization of dynamic multileaf collimator (DMLC) response with respiratory motion is critical to ensure the accuracy of DMLC-based four dimensional (4D) radiation delivery. In practice, however, a finite time delay (response time) between the acquisition of tumor position and multileaf collimator response necessitates predictive models of respiratory tumor motion to synchronize radiation delivery. Predicting a complex process such as respiratory motion introduces geometric errors, which have been reported in several publications. However, the dosimetric effect of such errors on 4D radiation delivery has not yet been investigated. Thus, our aim in this work was to quantify the dosimetric effects of geometric error due to prediction under several different conditions. Conformal and intensity modulated radiation therapy (IMRT) plans for a lung patient were generated for anterior-posterior/posterior-anterior (AP/PA) beam arrangements at 6 and 18 MV energies to provide planned dose distributions. Respiratory motion data was obtained from 60 diaphragm-motion fluoroscopy recordings from five patients. A linear adaptive filter was employed to predict the tumor position. The geometric error of prediction was defined as the absolute difference between predicted and actual positions at each diaphragm position. Distributions of geometric error of prediction were obtained for all of the respiratory motion data. Planned dose distributions were then convolved with distributions for the geometric error of prediction to obtain convolved dose distributions. The dosimetric effect of such geometric errors was determined as a function of several variables: response time (0-0.6 s), beam energy (6/18 MV), treatment delivery (3D/4D), treatment type (conformal/IMRT), beam direction (AP/PA), and breathing training type (free breathing/audio instruction/visual feedback). Dose difference and distance-to-agreement analysis was employed to quantify results. Based on our data, the dosimetric impact of prediction (a) increased with response time, (b) was larger for 3D radiation therapy as compared with 4D radiation therapy, (c) was relatively insensitive to change in beam energy and beam direction, (d) was greater for IMRT distributions as compared with conformal distributions, (e) was smaller than the dosimetric impact of latency, and (f) was greatest for respiration motion with audio instructions, followed by visual feedback and free breathing. Geometric errors of prediction that occur during 4D radiation delivery introduce dosimetric errors that are dependent on several factors, such as response time, treatment-delivery type, and beam energy. Even for relatively small response times of 0.6 s into the future, dosimetric errors due to prediction could approach delivery errors when respiratory motion is not accounted for at all. To reduce the dosimetric impact, better predictive models and/or shorter response times are required.

Photon-beam subsource sensitivity to the initial electron-beam parameters

One limitation to the widespread implementation of Monte Carlo (MC) patient dose-calculation algorithms for radiotherapy is the lack of a general and accurate source model of the accelerator radiation source. Our aim in this work is to investigate the sensitivity of the photon-beam subsource distributions in a MC source model (with target, primary collimator, and flattening filter photon subsources and an electron subsource) for 6- and 18-MV photon beams when the energy and radial distributions of initial electrons striking a linac target change. For this purpose, phase-space data (PSD) was calculated for various mean electron energies striking the target, various normally distributed electron energy spread, and various normally distributed electron radial intensity distributions. All PSD was analyzed in terms of energy, fluence, and energy fluence distributions, which were compared between the different parameter sets. The energy spread was found to have a negligible influence on the subsource distributions. The mean energy and radial intensity significantly changed the target subsource distribution shapes and intensities. For the primary collimator and flattening filter subsources, the distribution shapes of the fluence and energy fluence changed little for different mean electron energies striking the target, however, their relative intensity compared with the target subsource change, which can be accounted for by a scaling factor. This study indicates that adjustments to MC source models can likely be limited to adjusting the target subsource in conjunction with scaling the relative intensity and energy spectrum of the primary collimator, flattening filter, and electron subsources when the energy and radial distributions of the initial electron-beam change.

Monte Carlo source model for photon beam radiotherapy: photon source characteristics

A major barrier to widespread clinical implementation of Monte Carlo dose calculation is the difficulty in characterizing the radiation source within a generalized source model. This work aims to develop a generalized three-component source model (target, primary collimator, flattening filter) for 6- and 18-MV photon beams that match full phase-space data (PSD). Subsource by subsource comparison of dose distributions, using either source PSD or the source model as input, allows accurate source characterization and has the potential to ease the commissioning procedure, since it is possible to obtain information about which subsource needs to be tuned. This source model is unique in that, compared to previous source models, it retains additional correlations among PS variables, which improves accuracy at nonstandard source-to-surface distances (SSDs). In our study, three-dimensional (3D) dose calculations were performed for SSDs ranging from 50 to 200 cm and for field sizes from 1 x 1 to 30 x 30 cm2 as well as a 10 x 10 cm2 field 5 cm off axis in each direction. The 3D dose distributions, using either full PSD or the source model as input, were compared in terms of dose-difference and distance-to-agreement. With this model, over 99% of the voxels agreed within +/-1% or 1 mm for the target, within 2% or 2 mm for the primary collimator, and within +/-2.5% or 2 mm for the flattening filter in all cases studied. For the dose distributions, 99% of the dose voxels agreed within 1% or 1 mm when the combined source model-including a charged particle source and the full PSD as input-was used. The accurate and general characterization of each photon source and knowledge of the subsource dose distributions should facilitate source model commissioning procedures by allowing scaling the histogram distributions representing the subsources to be tuned.