Development of an x-ray-opaque-marker system for quantitative phantom positioning in patient-specific quality assurance.

PURPOSE: We developed an x-ray-opaque-marker (XOM) system with inserted fiducial markers for patient-specific quality assurance (QA) in CyberKnife (Accuray) and a general-purpose linear accelerator (linac). The XOM system can be easily inserted or removed from the existing patient-specific QA phantom. Our study aimed to assess the utility of the XOM system by evaluating the recognition accuracy of the phantom position error and estimating the dose perturbation around a marker. METHODS: The recognition accuracy of the phantom position error was evaluated by comparing the known error values of the phantom position with the values measured by matching the images with target locating system (TLS; Accuray) and on-board imager (OBI; Varian). The dose perturbation was evaluated for 6 and 10 MV single-photon beams through experimental measurements and Monte Carlo simulations. RESULTS: The root mean squares (RMSs) of the residual position errors for the recognition accuracy evaluation in translations were 0.07 mm with TLS and 0.30 mm with OBI, and those in rotations were 0.13° with TLS and 0.15° with OBI. The dose perturbation was observed within 1.5 mm for 6 MV and 2.0 mm for 10 MV from the marker. CONCLUSIONS: Sufficient recognition accuracy of the phantom position error was achieved using our system. It is unnecessary to consider the dose perturbation in actual patient-specific QA. We concluded that the XOM system can be utilized to ensure quantitative and accurate phantom positioning in patient-specific QA with CyberKnife and a general-purpose linac.

Measurements of temporal response of luminescence of water at lower energy than Cerenkov-light threshold during carbon-ion irradiation.

Although the luminescence of water at lower energy than the Cerenkov-light threshold during carbon-ion irradiation was found and imaging was possible, the temporal response has not been measured, and so the difference from Cerenkov-light remains unclear. To clarify this point, we measured the temporal response of the luminescence of water at lower energy than the Cerenkov-light threshold and compared it with that of Cerenkov-light. We used silicon photomultiplier (Si-PM) modules to measure the temporal response at the Bragg peak area of a water phantom during irradiation of the carbon ion where the Cerenkov-light was not included. We also measured the temporal response at the shallow depth of the water phantom where the Cerenkov-light was included. In both areas, we measured the temporal waveforms of the light produced by the irradiation of the carbon ions in which the ripples of spills were clearly observed. We found no difference in the waveforms between the Bragg peak and the shallow depths of water. Our results do not contradict the hypothesis that the luminescence of water and Cerenkov-light are produced by the same mechanism.

OPTIMIZATION OF AN ADDITIONAL COLLIMATOR IN A BEAM DELIVERY SYSTEM FOR REDUCTION OF THE SECONDARY NEUTRON EXPOSURE IN PASSIVE CARBON-ION THERAPY.

The aim of this work is to optimize an additional collimator in a beam delivery system to reduce neutron exposure to patients in passive carbon-ion therapy. All studies were performed by Monte Carlo simulation assuming the beam delivery system at Heavy-Ion Medical Accelerator in Chiba. We calculated the neutron ambient dose equivalent at patient positions with an additional collimator, and optimized the position, aperture size and material of the collimator to reduce the neutron ambient dose equivalent. The collimator located 125 and 470 cm upstream from the isocenter could reduce the dose equivalent near the isocenter by 35%, while the collimator located 813 cm upstream from the isocenter was ineffective. As for the material of the collimator, iron and nickel could conduct reduction slightly better than aluminum and polymethyl methacrylate. The additional collimator is an effective method for the reduction of the neutron ambient dose equivalent near the isocenter.

Addition of luminescence process in Monte Carlo simulation to precisely estimate the light emitted from water during proton and carbon-ion irradiation.

Although luminescence of water lower in energy than the Cerenkov-light threshold during proton and carbon-ion irradiation has been found, the phenomenon has not yet been implemented for Monte Carlo simulations. The results provided by the simulations lead to misunderstandings of the physical phenomenon in optical imaging of water during proton and carbon-ion irradiation. To solve the problems, as well as to clarify the light production of the luminescence of water, we modified a Monte Carlo simulation code to include the light production from the luminescence of water and compared them with the experimental results of luminescence imaging of water. We used GEANT4 for the simulation of emitted light from water during proton and carbon-ion irradiation. We used the light production from the luminescence of water using the scintillation process in GEANT4 while those of Cerenkov light from the secondary electrons and prompt gamma photons in water were also included in the simulation. The modified simulation results showed similar depth profiles to those of the measured data for both proton and carbon-ion. When the light production of 0.1 photons/MeV was used for the luminescence of water in the simulation, the simulated depth profiles showed the best match to those of the measured results for both the proton and carbon-ion compared with those used for smaller and larger numbers of photons/MeV. We could successively obtain the simulated depth profiles that were basically the same as the experimental data by using GEANT4 when we assumed the light production by the luminescence of water. Our results confirmed that the inclusion of the luminescence of water in Monte Carlo simulation is indispensable to calculate the precise light distribution in water during irradiation of proton and carbon-ion.

A patient-specific aperture system with an energy absorber for spot scanning proton beams: Verification for clinical application.

PURPOSE: In the authors' proton therapy system, the patient-specific aperture can be attached to the nozzle of spot scanning beams to shape an irradiation field and reduce lateral fall-off. The authors herein verified this system for clinical application. METHODS: The authors prepared four types of patient-specific aperture systems equipped with an energy absorber to irradiate shallow regions less than 4 g/cm(2). The aperture was made of 3-cm-thick brass and the maximum water equivalent penetration to be used with this system was estimated to be 15 g/cm(2). The authors measured in-air lateral profiles at the isocenter plane and integral depth doses with the energy absorber. All input data were obtained by the Monte Carlo calculation, and its parameters were tuned to reproduce measurements. The fluence of single spots in water was modeled as a triple Gaussian function and the dose distribution was calculated using a fluence dose model. The authors compared in-air and in-water lateral profiles and depth doses between calculations and measurements for various apertures of square, half, and U-shaped fields. The absolute doses and dose distributions with the aperture were then validated by patient-specific quality assurance. Measured data were obtained by various chambers and a 2D ion chamber detector array. RESULTS: The patient-specific aperture reduced the penumbra from 30% to 70%, for example, from 34.0 to 23.6 mm and 18.8 to 5.6 mm. The calculated field width for square-shaped apertures agreed with measurements within 1 mm. Regarding patient-specific aperture plans, calculated and measured doses agreed within -0.06% ± 0.63% (mean ± SD) and 97.1% points passed the 2%-dose/2 mm-distance criteria of the γ-index on average. CONCLUSIONS: The patient-specific aperture system improved dose distributions, particularly in shallow-region plans.

High resolution Cerenkov light imaging of induced positron distribution in proton therapy.

PURPOSE: In proton therapy, imaging of the positron distribution produced by fragmentation during or soon after proton irradiation is a useful method to monitor the proton range. Although positron emission tomography (PET) is typically used for this imaging, its spatial resolution is limited. Cerenkov light imaging is a new molecular imaging technology that detects the visible photons that are produced from high-speed electrons using a high sensitivity optical camera. Because its inherent spatial resolution is much higher than PET, the authors can measure more precise information of the proton-induced positron distribution with Cerenkov light imaging technology. For this purpose, they conducted Cerenkov light imaging of induced positron distribution in proton therapy. METHODS: First, the authors evaluated the spatial resolution of our Cerenkov light imaging system with a (22)Na point source for the actual imaging setup. Then the transparent acrylic phantoms (100 × 100 × 100 mm(3)) were irradiated with two different proton energies using a spot scanning proton therapy system. Cerenkov light imaging of each phantom was conducted using a high sensitivity electron multiplied charge coupled device (EM-CCD) camera. RESULTS: The Cerenkov light's spatial resolution for the setup was 0.76 ± 0.6 mm FWHM. They obtained high resolution Cerenkov light images of the positron distributions in the phantoms for two different proton energies and made fused images of the reference images and the Cerenkov light images. The depths of the positron distribution in the phantoms from the Cerenkov light images were almost identical to the simulation results. The decay curves derived from the region-of-interests (ROIs) set on the Cerenkov light images revealed that Cerenkov light images can be used for estimating the half-life of the radionuclide components of positrons. CONCLUSIONS: High resolution Cerenkov light imaging of proton-induced positron distribution was possible. The authors conclude that Cerenkov light imaging of proton-induced positron is promising for proton therapy.

Assessment of spatial uncertainty in computed tomography-based Gamma Knife stereotactic radiosurgery process with automated positioning system.

BACKGROUND: In this study, we assessed the geometric accuracy of an automated positioning system in Gamma Knife (GK) surgery. Specifically, we looked at the total spatial uncertainty over the entire treatment range of GK stereotactic radiosurgery (SRS) procedures in both the GK model C and the Perfexion (PFX). METHODS: An originally-developed phantom and a radiochromic film were used for obtaining actual dose distributions. The phantom, with inserted films on different axial planes (z = 60, 75, 100, 125, 140 mm), sagittal planes (x = 60, 75, 100, 125, 140 mm), and coronal planes (y = 60, 75, 100, 125, 140 mm), was placed on a Leksell skull frame. Computed tomography (CT) was then performed with a stereotactic localizer box attached to the frame, and dose planning was made using the Leksell GammaPlan treatment planning system. The phantom finally received beam delivery using a single shot of a 4-mm collimator helmet. The discrepancy between the planned shot position and the irradiated center position was evaluated by a dedicated film analysis software. RESULTS: The total uncertainty of CT-based GK SRS was less than 1 mm for almost all measured points over the stereotactic space in both the model C and the PFX. In addition, the geometric accuracy of the automated positioning system was estimated to be less than 0.1 mm and equal to 0.5 mm in the central and peripheral areas, respectively. CONCLUSIONS: We confirmed that the total spatial uncertainties of both the GK model C and the PFX are acceptable for clinical use.

Dosimetric comparison of absolute and relative dose distributions between tissue maximum ratio and convolution algorithms for acoustic neurinoma plans in Gamma Knife radiosurgery.

BACKGROUND: The treatment planning for Gamma Knife (GK) stereotactic radiosurgery (SRS) that performs dose calculations based on tissue maximum ratio (TMR) algorithm has disadvantages in predicting dose in tissue heterogeneity. The latest version of the planning software is equipped with a convolution dose algorithm as an optional extra and the new algorithm is able to compensate for head inhomogeneity. However, the effect of this improved calculation method requires detailed validation in clinical cases. In this study, we compared absolute and relative dose distributions of treatment plans for acoustic neurinoma between TMR and the convolution calculation. METHODS: Twenty-nine clinically used plans created by TMR algorithm were recalculated by convolution method. Differences between TMR and convolution were evaluated in terms of absolute dose (beam-on time), dosimetric parameters including target coverage, selectivity, conformity index, gradient index, radical homogeneity index and the dose-volume relationship. RESULTS: The discrepancy in estimated absolute dose to the target ranged from 1 to 7 % between TMR and convolution. In addition, dosimetric parameters of the two methods achieved statistical significance. However, it was difficult to see the change of relative dose distribution by visual assessment on a monitor. CONCLUSIONS: Convolution, heterogeneity correction calculation, and the algorithm are necessary to reduce the dosimetric uncertainty of each case in GK SRS.

[Assessment of overall spatial accuracy in image guided stereotactic body radiotherapy using a spine registration method].

Stereotactic body radiotherapy (SBRT) for lung and liver tumors is always performed under image guidance, a technique used to confirm the accuracy of setup positioning by fusing planning digitally reconstructed radiographs with X-ray, fluoroscopic, or computed tomography (CT) images, using bony structures, tumor shadows, or metallic markers as landmarks. The Japanese SBRT guidelines state that bony spinal structures should be used as the main landmarks for patient setup. In this study, we used the Novalis system as a linear accelerator for SBRT of lung and liver tumors. The current study compared the differences between spine registration and target registration and calculated total spatial accuracy including setup uncertainty derived from our image registration results and the geometric uncertainty of the Novalis system. We were able to evaluate clearly whether overall spatial accuracy is achieved within a setup margin (SM) for planning target volume (PTV) in treatment planning. After being granted approval by the Hospital and University Ethics Committee, we retrospectively analyzed eleven patients with lung tumor and seven patients with liver tumor. The results showed the total spatial accuracy to be within a tolerable range for SM of treatment planning. We therefore regard our method to be suitable for image fusion involving 2-dimensional X-ray images during the treatment planning stage of SBRT for lung and liver tumors.

Microdosimetric calculation of penumbra for biological dose in wobbled carbon-ion beams with Monte Carlo Method.

In carbon-ion radiotherapy, it is important to evaluate the biological dose because the relative biological effectiveness values vary greatly in a patient's body. The microdosimetric kinetic model (MKM) is a method of estimating the biological effect of radiation by use of microdosimetry. The lateral biological dose distributions were estimated with a modified MKM, in which we considered the overkilling effect in the high linear-energy-transfer region. In this study, we used the Monte Carlo calculation of the Geant4 code to simulate a horizontal port at the Heavy Ion Medical Accelerator in Chiba of the National Institute of Radiological Sciences. The lateral biological dose distributions calculated by Geant4 were almost flat as the lateral absorbed dose in the flattened area. However, in the penumbra region, the lateral biological dose distributions were sharper than the lateral absorbed dose distributions. Furthermore, the differences between the lateral absorbed dose and biological dose distributions were dependent on the depth for each multi-leaf collimator opening size. We expect that the lateral biological dose distribution presented here will enable high-precision calculations for a treatment-planning system.

Analytical design method of optimum ridge filters for wobbled and collimated proton beams.

A novel design method of ridge filters (RFs) has been developed for general proton beam lines which use a single-radius beam wobbling method. It can be applied to beam lines that transport both protons and carbon ions which are about three times longer than regular beam lines dedicated to protons. We designed an RF with an SOBP (spread-out Bragg peak) width of 60 mm in water for the 160-MeV proton beam of the HIMAC (Heavy Ion Medical Accelerator in Chiba) biology beam line using an existing model of the RF. Yet we observed a slope in the SOBP region when we used the RF. To elucidate the source of the slope, we have developed a new calculation model taking into account the geometry of the RF and a beam-limiting device. The source for the slope was found to be the large scattering effect of protons in the RF and beam restriction by a ring collimator (aperture diameter: 160 mm) placed just before the RF. When both fluence reduction by the scattering effect of protons in the RF and the beam-collimation effect are taken into account, proper RFs can be designed universally for a given beam line arrangement using the single-radius beam-wobbling method from the start without any trial-and-error process. This will serve to reduce the commissioning time of newly designed beam delivery systems.