Macro Monte Carlo for dose calculation of proton beams

Although the Monte Carlo (MC) method allows accurate dose calculation for proton radiotherapy, its usage is limited due to long computing time. In order to gain efficiency, a new macro MC (MMC) technique for proton dose calculations has been developed. The basic principle of the MMC transport is a local to global MC approach. The local simulations using GEANT4 consist of mono-energetic proton pencil beams impinging perpendicularly on slabs of different thicknesses and different materials (water, air, lung, adipose, muscle, spongiosa, cortical bone). During the local simulation multiple scattering, ionization as well as elastic and inelastic interactions have been taken into account and the physical characteristics such as lateral displacement, direction distributions and energy loss have been scored for primary and secondary particles. The scored data from appropriate slabs is then used for the stepwise transport of the protons in the MMC simulation while calculating the energy loss along the path between entrance and exit position. Additionally, based on local simulations the radiation transport of neutrons and the generated ions are included into the MMC simulations for the dose calculations. In order to validate the MMC transport, calculated dose distributions using the MMC transport and GEANT4 have been compared for different mono-energetic proton pencil beams impinging on different phantoms including homogeneous and inhomogeneous situations as well as on a patient CT scan. The agreement of calculated integral depth dose curves is better than 1% or 1 mm for all pencil beams and phantoms considered. For the dose profiles the agreement is within 1% or 1 mm in all phantoms for all energies and depths. The comparison of the dose distribution calculated using either GEANT4 or MMC in the patient also shows an agreement of within 1% or 1 mm. The efficiency of MMC is up to 200 times higher than for GEANT4. The very good level of agreement in the dose comparisons demonstrate that the newly developed MMC transport results in very accurate and efficient dose calculations for proton beams.

Intracavity Generation of Radially Polarized CO2 Laser Beams Based on a Simple Binary Dielectric Diffraction Grating

We present experimental results on the intracavity generation of radially polarized light by incorporation of a polarization-selective mirror in a CO2 -laser resonator. The selectivity is achieved with a simple binary dielectric diffraction grating etched in the backsurface of the mirror substrate. Very high polarization selectivity was achieved, and good agreement of simulation and experimental results is shown. The overall radial polarization purity of the generated laser beam was found to be higher than 90% .

Irradiation of intracerebral 9L gliosarcoma by a single array of microplanar x-ray beams from a synchrotron: balance between curing and sparing

The purpose of this work was the understanding of microbeam radiation therapy at the ESRF in order to find the best compromise between curing of tumors and sparing of normal tissues, to obtain a better understanding of survival curves and to report its efficiency. This method uses synchrotron-generated x-ray microbeams. Rats were implanted with 9L gliosarcomas and the tumors were diagnosed by MRI. They were irradiated 14 days after implantation by arrays of 25 microm wide microbeams in unidirectional mode, with a skin entrance dose of 625 Gy. The effect of using 200 or 100 microm center-to-center spacing between the microbeams was compared. The median survival time (post-implantation) was 40 and 67 days at 200 and 100 microm spacing, respectively. However, 72% of rats irradiated at 100 microm spacing showed abnormal clinical signs and weight patterns, whereas only 12% of rats were affected at 200 microm spacing. In parallel, histological lesions of the normal brain were found in the 100 microm series only. Although the increase in lifespan was equal to 273% and 102% for the 100 and 200 microm series, respectively, the 200 microm spacing protocol provides a better sparing of healthy tissue and may prove useful in combination with other radiation modalities or additional drugs.

VMC++ versus BEAMnrc: a comparison of simulated linear accelerator heads for photon beams

BEAMnrc, a code for simulating medical linear accelerators based on EGSnrc, has been bench-marked and used extensively in the scientific literature and is therefore often considered to be the gold standard for Monte Carlo simulations for radiotherapy applications. However, its long computation times make it too slow for the clinical routine and often even for research purposes without a large investment in computing resources. VMC++ is a much faster code thanks to the intensive use of variance reduction techniques and a much faster implementation of the condensed history technique for charged particle transport. A research version of this code is also capable of simulating the full head of linear accelerators operated in photon mode (excluding multileaf collimators, hard and dynamic wedges). In this work, a validation of the full head simulation at 6 and 18 MV is performed, simulating with VMC++ and BEAMnrc the addition of one head component at a time and comparing the resulting phase space files. For the comparison, photon and electron fluence, photon energy fluence, mean energy, and photon spectra are considered. The largest absolute differences are found in the energy fluences. For all the simulations of the different head components, a very good agreement (differences in energy fluences between VMC++ and BEAMnrc <1%) is obtained. Only a particular case at 6 MV shows a somewhat larger energy fluence difference of 1.4%. Dosimetrically, these phase space differences imply an agreement between both codes at the <1% level, making VMC++ head module suitable for full head simulations with considerable gain in efficiency and without loss of accuracy.

Testing the GlaaS algorithm for dose measurements on low- and high-energy photon beams using an amorphous silicon portal imager

The GLAaS algorithm for pretreatment intensity modulation radiation therapy absolute dose verification based on the use of amorphous silicon detectors, as described in Nicolini et al. [G. Nicolini, A. Fogliata, E. Vanetti, A. Clivio, and L. Cozzi, Med. Phys. 33, 2839-2851 (2006)], was tested under a variety of experimental conditions to investigate its robustness, the possibility of using it in different clinics and its performance. GLAaS was therefore tested on a low-energy Varian Clinac (6 MV) equipped with an amorphous silicon Portal Vision PV-aS500 with electronic readout IAS2 and on a high-energy Clinac (6 and 15 MV) equipped with a PV-aS1000 and IAS3 electronics. Tests were performed for three calibration conditions: A: adding buildup on the top of the cassette such that SDD-SSD = d(max) and comparing measurements with corresponding doses computed at d(max), B: without adding any buildup on the top of the cassette and considering only the intrinsic water-equivalent thickness of the electronic portal imaging devices device (0.8 cm), and C: without adding any buildup on the top of the cassette but comparing measurements against doses computed at d(max). This procedure is similar to that usually applied when in vivo dosimetry is performed with solid state diodes without sufficient buildup material. Quantitatively, the gamma index (gamma), as described by Low et al. [D. A. Low, W. B. Harms, S. Mutic, and J. A. Purdy, Med. Phys. 25, 656-660 (1998)], was assessed. The gamma index was computed for a distance to agreement (DTA) of 3 mm. The dose difference deltaD was considered as 2%, 3%, and 4%. As a measure of the quality of results, the fraction of field area with gamma larger than 1 (%FA) was scored. Results over a set of 50 test samples (including fields from head and neck, breast, prostate, anal canal, and brain cases) and from the long-term routine usage, demonstrated the robustness and stability of GLAaS. In general, the mean values of %FA remain below 3% for deltaD equal or larger than 3%, while they are slightly larger for deltaD = 2% with %FA in the range from 3% to 8%. Since its introduction in routine practice, 1453 fields have been verified with GLAaS at the authors' institute (6 MV beam). Using a DTA of 3 mm and a deltaD of 4% the authors obtained %FA = 0.9 +/- 1.1 for the entire data set while, stratifying according to the dose calculation algorithm, they observed: %FA = 0.7 +/- 0.9 for fields computed with the analytical anisotropic algorithm and %FA = 2.4 +/- 1.3 for pencil-beam based fields with a statistically significant difference between the two groups. If data are stratified according to field splitting, they observed %FA = 0.8 +/- 1.0 for split fields and 1.0 +/- 1.2 for nonsplit fields without any significant difference.

Utilizing Gaussian-Schell model beams to mitigate atmospheric turbulence in free space optical communications / by Kyle R. Drexler.

Turbulence affects traditional free space optical communication by causing speckle to appear in the received beam profile. This occurs due to changes in the refractive index of the atmosphere that are caused by fluctuations in temperature and pressure, resulting in an inhomogeneous medium. The Gaussian-Schell model of partial coherence has been suggested as a means of mitigating these atmospheric inhomogeneities on the transmission side. This dissertation analyzed the Gaussian-Schell model of partial coherence by verifying the Gaussian-Schell model in the far-field, investigated the number of independent phase control screens necessary to approach the ideal Gaussian-Schell model, and showed experimentally that the Gaussian-Schell model of partial coherence is achievable in the far-field using a liquid crystal spatial light modulator. A method for optimizing the statistical properties of the Gaussian-Schell model was developed to maximize the coherence of the field while ensuring that it does not exhibit the same statistics as a fully coherent source. Finally a technique to estimate the minimum spatial resolution necessary in a spatial light modulator was developed to effectively propagate the Gaussian-Schell model through a range of atmospheric turbulence strengths. This work showed that regardless of turbulence strength or receiver aperture, transmitting the Gaussian-Schell model of partial coherence instead of a fully coherent source will yield a reduction in the intensity fluctuations of the received field. By measuring the variance of the intensity fluctuations and the received mean, it is shown through the scintillation index that using the Gaussian-Schell model of partial coherence is a simple and straight forward method to mitigate atmospheric turbulence instead of traditional adaptive optics in free space optical communications.

Accelerating Diffraction-free Beams in Photonic Lattices

We study nondiffracting accelerating paraxial optical beams in periodic potentials, in both the linear and the nonlinear domains. In particular, we show that only a unique class of z-dependent lattices can support a true accelerating diffractionless beam. Accelerating lattice solitons, autofocusing beams and accelerating bullets in optical lattices are systematically examined.

Characterization of Optically Stimulated Luminescent Detectors in Photon & Proton Beams for Use in Anthropomorphic Phantoms

This study investigated characteristics of optically stimulated luminescent detectors (OSLDs) in protons, allowing comparison to thermoluminescent detectors, and to be implemented into the Radiological Physics Center’s (RPC) remote audit quality assurance program for protons, and for remote anthropomorphic phantom irradiations. The OSLDs used were aluminum oxide (Al2O3:C) nanoDots from Landauer, Inc. (Glenwood, Ill.) measuring 10x10x2 mm3. A square, 20(L)x20(W)x0.5(H) cm3 piece of solid water was fabricated with pockets to allow OSLDs and TLDs to be irradiated simultaneously and perpendicular to the beam. Irradiations were performed at 5cm depth in photons, and in the center of a 10 cm SOBP in a 200MeV proton beam. Additionally, the Radiological Physics Center’s anthropomorphic pelvic phantom was used to test the angular dependence of OSLDs in photons and protons. A cylindrical insert in the phantom allows the dosimeters to be rotated to any angle with a fixed gantry angle. OSLDs were irradiated at 12 angles between 0 and 360 degrees. The OSLDs were read out with a MicroStar reader from Landauer, Inc. Dose response indicates that at angles where the dosimeter is near parallel with the radiation beam response is reduced slightly. Measurements in proton beams do not show significant angular dependence. Post-irradiation fading of OSLDs was studied in proton beams to determine if the fading was different than that of photons. The fading results showed no significant difference from results in photon beams. OSLDs and TLDs are comparable within 3% in photon beams and a correction factor can be posited for proton beams. With angular dependence characteristics defined, OSLDs can be implemented into multiple-field treatment plans in photons and protons and used in the RPC’s quality assurance program.

Collimated Field Emission Beams from Metal Double-Gate Nanotip Arrays Optically Excited via Surface Plasmon Resonance

The generation of collimated electron beams from metal double-gate nanotip arrays excited by near infrared laser pulses is studied. Using electromagnetic and particle tracking simulations, we showed that electron pulses with small rms transverse velocities are efficiently produced from nanotip arrays by laser-induced field emission with the laser wavelength tuned to surface plasmon polariton resonance of the stacked double-gate structure. The result indicates the possibility of realizing a metal nanotip array cathode that outperforms state-of-the-art photocathodes.