A diamond target for megavoltage cone-beam CT.

PURPOSE: To determine the properties of a megavoltage cone-beam CT system using the unflattened beam from a sintered diamond target at 4 and 6 MV. METHODS: A sintered diamond target was used in place of a graphite target as part of an imaging beam line (an unflattened beam from a graphite target) installed on a linear accelerator. The diamond target, with a greater density than the graphite target, permitted imaging at the lower beam energy (4 MV) required with the graphite target and the higher beam energy (6 MV) conventionally used with the tungsten/stainless steel target and stainless steel flattening filter. Images of phantoms and patients were acquired using the different beam lines and compared. The beam spectra and dose distributions were determined using Monte Carlo simulation. RESULTS: The diamond target allowed use of the same beam energy as for treatment, simplifying commissioning and quality assurance. Images acquired with the diamond target at 4 MV were similar to those obtained with the graphite target at 4 MV. The slight reduction in low energy photons due to the higher-Z sintering material in the diamond target had minimal effect on image quality. Images acquired at 6 MV with the diamond target showed a small decrease in contrast-to-noise ratio, resulting from a decrease in the fraction of photons in the beam in the energy range to which the detector is most sensitive. CONCLUSIONS: The diamond target provides images of a similar quality to the graphite target. Diamond allows use of the higher beam energy conventionally used for treatment, provides a higher dose rate for the same beam current, and potentially simplifies installation and maintenance of the beam line.

Simulation of large x-ray fields using independently measured source and geometry details.

PURPOSE: Obtain an accurate simulation of the dose from the 6 and 18 MV x-ray beams from a Siemens Oncor linear accelerator by comparing simulation to measurement. Constrain the simulation by independently determining parameters of the treatment head and incident beam, in particular, the energy and spot size. METHODS: Measurements were done with the treatment head in three different configurations: (1) The clinical configuration, (2) the flattening filter removed, and (3) the target and flattening filter removed. Parameters of the incident beam and treatment head were measured directly. Incident beam energy and spectral width were determined from the percent-depth ionization of the raw beam (as described previously), spot size was determined using a spot camera, and the densities of the flattening filters were determined by weighing them. Simulations were done with EGSnrc/BEAMnrc code. An asymmetric simulation was used, including offsets of the spot, primary collimator, and flattening filter from the collimator rotation axis. RESULTS: Agreement between measurement and simulation was obtained to the least restrictive of 1% or 1 mm at 6 MV, both with and without the flattening filter in place, except for the buildup region. At 18 MV, the agreement was 1.5%/1.5 mm with the flattening filter in place and 1%/1 mm with it removed, except for in the buildup region. In the buildup region, the discrepancy was 2%/2 mm at 18 MV and 1.5%/1.5 mm at 6 MV with the flattening filter either removed or in place. The methodology for measuring the source and geometry parameters for the treatment head simulation is described. Except to determine the density of the flattening filter, no physical modification of the treatment head is necessary to obtain those parameters. In particular, the flattening filter does not need to be removed as was done in this work. CONCLUSIONS: Good agreement between measured and simulated dose distributions was obtained, even in the buildup region. The simulation was tightly constrained by independent measurements of parameters of the incident beam and treatment head. The method of obtaining the input parameters is described, and can be carried out on a clinical linear accelerator.

Determination of electron energy, spectral width, and beam divergence at the exit window for clinical megavoltage x-ray beams.

Monte Carlo simulations of x-ray beams typically take parameters of the electron beam in the accelerating waveguide to be free parameters. In this paper, a methodology is proposed and implemented to determine the energy, spectral width, and beam divergence of the electron source. All treatment head components were removed from the beam path, leaving only the exit window. With the x-ray target and flattener out of the beam, uncertainties in physical characteristics and relative position of the target and flattening filter, and in spot size, did not contribute to uncertainty in the energy. Beam current was lowered to reduce recombination effects. The measured dose distributions were compared with Monte Carlo simulation of the electron beam through the treatment head to extract the electron source characteristics. For the nominal 6 and 18 MV x-ray beams, the energies were 6.51 +/- 0.15 and 13.9 +/- 0.2 MeV, respectively, with the uncertainties resulting from uncertainties in the detector position in the measurement and in the stopping power in the simulations. Gaussian spectral distributions were used, with full widths at half maximum ranging from 20 +/- 4% at 6 MV to 13 +/- 4% at 18 MV required to match the fall-off portion of the percent-depth ionization curve. Profiles at the depth of maximum dose from simulations that used the manufacturer-specified exit window geometry and no beam divergence were 2-3 cm narrower than measured profiles. Two simulation configurations yielding the measured profile width were the manufacturer-specified exit window thickness with electron source divergences of 3.3 degrees at 6 MV and 1.8 degrees at 18 MV and an exit window 40% thicker than the manufacturer's specification with no beam divergence. With the x-ray target in place (and no flattener), comparison of measured to simulated profiles sets upper limits on the electron source divergences of 0.2 degrees at 6 MV and 0.1 degrees at 18 MV. A method of determining source characteristics without mechanical modification of the treatment head, and therefore feasible in clinics, is presented. The energies and spectral widths determined using this method agree with those determined with only the exit window in the beam path.

Comparison of electron scattering algorithms in Geant4.

Electron scattering algorithms in Geant4 versions 9.4 and 9.5 were benchmarked by comparing scattered distributions against previously measured values at 13 and 20 MeV, for low, intermediate, and high atomic number materials. Several scattering models were used: Versions 93 and 95 of the Urban model, with different step size limits near boundaries; Goudsmit-Saunderson multiple scattering; and single scattering. The Urban93 and Urban95 models with a large step size limit (as in the Option 0 physics list) were found to give results most closely matching the experimental results. Scattered distributions using the Urban models were all narrower than measured by up to 6%, consistent with previous published simulations using EGSnrc. This is suggestive of a systematic difference between simulations and measurement. The magnitudes of the differences were similar to previously published results using Geant4, although there were differences in detail. In particular, the current results were typically 2% narrower than values. Results with the more restrictive step size limit in Option 3 were even more narrow, and close to those with single scattering. The Goudsmit-Saunderson multiple scattering model produced distributions up to 15% different from measured in Geant4 version 9.5 and up to 45% different in Geant4 version 9.4.

Evaluation of motion management strategies based on required margins.

Strategies for delivering radiation to a moving lesion each require a margin to compensate for uncertainties in treatment. These motion margins have been determined here by separating the total uncertainty into components. Probability density functions for the individual sources of uncertainty were calculated for ten motion traces obtained from the literature. Motion margins required to compensate for the center of mass motion of the clinical treatment volume were found by convolving the individual sources of uncertainty. For measurements of position at a frequency of 33 Hz, system latency was the dominant source of positional uncertainty. Averaged over the ten motion traces, the motion margin for tracking with a latency of 200 ms was 4.6 mm. Gating with a duty cycle of 33% required a mean motion margin of 3.2-3.4 mm, and tracking with a latency of 100 ms required a motion margin of 3.1 mm. Feasible reductions in the effects of the sources of uncertainty, for example by using a simple prediction algorithm to anticipate the lesion position at the end of the latency period, resulted in a mean motion margin of 1.7 mm for tracking with a latency of 100 ms, 2.4 mm for tracking with a latency of 200 ms, and 2.1-2.2 mm for the gating strategies with duty cycles of 33%. A crossover tracking latency of 150 ms was found, below which tracking strategies could take advantage of narrower motion margins than gating strategies. The methods described here provide a means to guide selection of a motion management strategy for a given patient.

SU-E-T-499: Validation of the Varian Generic Phase Space Files for Monte Carlo Calculations of Dose Distributions for the TrueBeam Linac Head.

PURPOSE: To validate the generic phase space files for Varian TrueBeam linac head simulations. METHODS: The generic phase space files include the simulation results of 6MV, 10MV, 6MV FFF, and 10MV FFF (flattening-filter free) operating modes of TrueBeam for patient-independent linac head components. Using the generic phase space files as the radiation sources, the BEAMnrc Monte Carlo codes are used to simulate the patient-dependent parts of the TrueBeam linac and the resulting phase space files are generated at a plane just before entering a water phantom for 4 different field sizes (5×5, 10×10, 20×20, and 40×40 cm2 ). Dose distributions are calculated by DOSXYZnrc in the water phantom of size 50×50×40 cm3 . The percentage-depth-dose (PDD) curves and lateral dose profiles at three different depths (dmax, 10cm, 20cm) are obtained. Comprehensive comparisons have been made for a total of 64 dose profiles (including PDDs) between the Monte Carlo calculations and the measured data. The gamma index analysis is performed for all the comparisons. RESULTS: The matching of the calculated dose distributions to the measured ones is analyzed by the gamma index method with a criterion of 2% dose tolerance and 2 mm distance-to-agreement. Of the 64 comparisons, the minimum gamma index passing rate is at least 92%, after taking into account the statistical nature of the Monte Carlo calculated dose values. Despite the existence of latent variance of phase space files, the phantom dose calculation uncertainty can be less than 1% for field sizes as small as 5×5 cm2 . The computing time saved by using phase space files could be a factor of 5-10. CONCLUSIONS: The Varian generic phase space files are accurate and efficient radiation sources for Monte Carlo calculations of radiation dose distributions for TrueBeam linac head. This work was supported in part by Varian Medical Systems and the NIH (1R01 CA104205 and 1R21 CA153587).

Thermal conductivity of spin-polarized liquid 3He.

We present the first measurements of the thermal conductivity of spin-polarized normal liquid 3He. Using the rapid melting technique to produce nuclear polarizations up to 0.7, and a vibrating wire both as a heater and a thermometer, we show that, unlike the viscosity, the conductivity increases much less than predicted for s-wave scattering. We suggest that this might be due to a small probability for head-on collisions between quasiparticles.