安井 啓祐

Background This study evaluates dose perturbations caused by nonradioactive seeds in clinical cases by employing treatment planning system-based Monte Carlo (TPS-MC) simulation. Methodology We investigated dose perturbation using a water-equivalent phantom and 20 clinical cases of prostate cancer (10 cases with seeds and 10 cases without seeds) treated at Fujita Health University Hospital, Japan. First, dose calculations for a simple geometry were performed using the RayStation MC algorithm for a water-equivalent phantom with and without a seed. TPS-independent Monte Carlo (full-MC) simulations and film measurements were conducted to verify the accuracy of TPS-MC simulation. Subsequently, dose calculations using TPS-MC were performed on CT images of clinical cases of prostate cancer with and without seeds, and the dose distributions were compared. Results In clinical cases, dose calculations using MC simulations revealed hotspots around the seeds. However, the size of the hotspot was not correlated with the number of seeds. The maximum difference in dose perturbation between TPS-MC simulations and film measurements was 3.9%, whereas that between TPS-MC simulations and full-MC simulations was 3.7%. The dose error of TPS-MC was negligible for multiple beams or rotational irradiation. Conclusions Hotspots were observed in dose calculations using TPS-MC performed on CT images of clinical cases with seeds. The dose calculation accuracy around the seeds using TPS-MC simulations was comparable to that of film measurements and full-MC simulations, with differences within 3.9%. Although the clinical impact of hotspots occurring around the seeds is minimal, utilizing MC simulations on TPSs can be beneficial to verify their presence.

Deep learning-based CT image reconstruction (DLR) is a state-of-the-art method for obtaining CT images. This study aimed to evaluate the usefulness of DLR in radiotherapy. Data were acquired using a large-bore CT system and an electron density phantom for radiotherapy. We compared the CT values, image noise, and CT value-to-electron density conversion table of DLR and hybrid iterative reconstruction (H-IR) for various doses. Further, we evaluated three DLR reconstruction strength patterns (Mild, Standard, and Strong). The variations of CT values of DLR and H-IR were large at low doses, and the difference in average CT values was insignificant with less than 10 HU at doses of 100 mAs and above. DLR showed less change in CT values and smaller image noise relative to H-IR. The noise-reduction effect was particularly large in the low-dose region. The difference in image noise between DLR Mild and Standard/Strong was large, suggesting the usefulness of reconstruction intensities higher than Mild. DLR showed stable CT values and low image noise for various materials, even at low doses; particularly for Standard or Strong, the reduction in image noise was significant. These findings indicate the usefulness of DLR in treatment planning using large-bore CT systems.

In single-isocenter multiple-target stereotactic radiotherapy (SIMT-SRT), it is difficult to evaluate both the geometrical accuracy and absorbed dose measurement when irradiating off-isocenter targets. This study aimed to develop a simple quality assurance (QA) method to evaluate off-isocenter irradiation position accuracy in SIMT-SRT and compare its feasibility with that of a commercial device. First, we created two types of inserts and metallic balls with a diameter of 5 mm to be inserted into a commercially available phantom (SIMT phantom). Second, we developed a dedicated analysis software using Python for the Winston-Lutz test (WLT). Third, an image processing software, including the filtered back-projection algorithm, was developed to analyze the images obtained using an electronic portal imaging device (EPID). Fourth, the feasibility of our method was evaluated by comparing it with the results of WLT using two commercially available phantoms: WL-QA and MultiMet-WL cubes. Notably, 92% of the results in one-dimensional deviations were within 0.26 mm (EPID pixel width). The correlation coefficients were 0.52, 0.92, and 0.96 in the left-right, superior-inferior, and anterior-posterior directions, respectively. In the WLT, a maximum two-dimensional deviation of 0.70 mm was detected in our method, while the deviation in the other method was within 0.5 mm. The advantage of our method is that it can evaluate the geometrical accuracy at any gantry angle during dynamic rotation irradiation using a filtered back-projection algorithm, even if the target is located off the isocenter. Our method can perform WLT at arbitrary positions and is suitable for the QA of dynamic rotation irradiation using an EPID.

The purpose of this study was to evaluate the effect of quality assurance (QA)-related setup errors in passive proton therapy for prostate cancer with and without a hydrogel spacer. We used 20 typical computed tomography (CT) images of prostate cancer: 10 patients with and 10 patients without spacers. The following 12 model errors were assumed: output error ± 2%, range error ± 1 mm, setup error ± 1 mm for three directions, and multileaf collimator (MLC) position error ± 1 mm. We created verification plans with model errors and compared the prostate-rectal (PR) distance and dose indices with and without the spacer. The mean PR distance at the isocenter was 1.1 ± 1.3 mm without the spacer and 12.9 ± 2.9 mm with the spacer (P < 0.001). The mean rectum V53.5 GyE, V50 GyE, and V34.5 GyE in the original plan were 2.3%, 4.1%, and 12.1% without the spacer and 0.1%, 0.4%, and 3.3% with the spacer (P = 0.0011, < 0.001, and < 0.001). The effects of the range and lateral setup errors were small; however, the effects of the vertical/long setup and MLC error were significant in the cases without the spacer. The means of the maximum absolute change from original plans across all scenarios in the rectum V53.5 GyE, V50 GyE, and V34.5 GyE were 1.3%, 1.5%, and 2.3% without the spacer, and 0.2%, 0.4%, and 1.3% with the spacer (P < 0.001, < 0.001, and = 0.0019). This study indicated that spacer injections were also effective in reducing the change in the rectal dose due to setup errors.

A radiophotoluminescent glass dosimeter (RGD) is used for a postal audit of a photon beam because of its various excellent characteristics. However, it has not been used for scanning proton beams because its response characteristics have not been verified. In this study, the response of RGD to scanning protons was investigated to develop a dosimetry protocol using the linear energy transfer (LET)-based correction factor. The responses of RGD to four maximum-range-energy-pattern proton beams were verified by comparing it with ionization chamber (IC) dosimetry. The LET at each measurement depth was calculated via Monte Carlo (MC) simulation. The LET correction factor ( k LET RGD ) was the ratio between the uncorrected RGD dose ( D raw RGD ) and the IC dose at each measurement depth. k LET RGD can be represented as a function of LET using the following equation: k LET RGD LET = - 0.035 LET + 1.090 . D raw RGD showed a linear under-response with increasing LET, and the maximum dose difference between the IC dose and D raw RGD was 15.2% at an LET of 6.07 keV/μm. The LET-based correction dose ( D LET RGD ) conformed within 3.6% of the IC dose. The mean dose difference (±SD) of D raw RGD and D LET RGD was -2.5 ± 6.9% and 0.0 ± 1.6%, respectively. To achieve accurate dose verification for scanning proton beams using RGD, we derived a linear regression equation based on LET. The results show that with appropriate LET correction, RGD can be used for dose verification of scanning proton beams.

<title>Abstract</title> To understand the current state of flattening filter-free (FFF) beam implementation in C-arm linear accelerators (LINAC) in Japan, the quality assurance (QA)/quality control (QC) 2018–2019 Committee of the Japan Society of Medical Physics (JSMP) conducted a 37-question survey, designed to investigate facility information and specifications regarding FFF beam adoption and usage. The survey comprised six sections: facility information, devices, clinical usage, standard calibration protocols, modeling for treatment planning (TPS) systems and commissioning and QA/QC. A web-based questionnaire was developed. Responses were collected between 18 June and 18 September 2019. Of the 846 institutions implementing external radiotherapy, 323 replied. Of these institutions, 92 had adopted FFF beams and 66 had treated patients using them. FFF beams were used in stereotactic radiation therapy (SRT) for almost all disease sites, especially for the lungs using 6 MV and liver using 10 MV in 51 and 32 institutions, respectively. The number of institutions using FFF beams for treatment increased yearly, from eight before 2015 to 60 in 2018. Farmer-type ionization chambers were used as the standard calibration protocol in 66 (72%) institutions. In 73 (80%) institutions, the beam-quality conversion factor for FFF beams was calculated from TPR20,10, via the same protocol used for beams with flattening filter (WFF). Commissioning, periodic QA and patient-specific QA for FFF beams also followed the procedures used for WFF beams. FFF beams were primarily used in high-volume centers for SRT. In most institutions, measurement and QA was conducted via the procedures used for WFF beams.

Background: The purpose of this study was to investigate the dosimetric characteristics of a new type of two-dimensional diode detector array used for quality assurance of stereotactic radiosurgery (SRS). Materials and Methods: The devices used in this study were the SRS MapCHECK detector and the StereoPHAN. The detector has 1013 diode detectors over an area of 77 x 77 mm2. The reproducibility, dose linearity, dose rate dependencies, output factors (OPFs) and angular dependencies were investigated as dosimetric characteristics. The OPFs were measured and compared between AP and PA direction ranging from 0.5 x 0.5 to 7 x 7 cm2. The angular dependencies were measured using 19 gantry angles. Results: The dose reproducibility and linearities showed sufficient performance of 6 MV and 10 MV. At 40 MU/min, there was a 1.3% difference from the ionization chamber measurements. For the flattening filter-free beam, there was no dose rate dependency from the 400 MU/minute to 2400 MU/minute, and the variation was within 0.5%. For small irradiation fields of 1 cm or less, the measured value of the SMC differed in AP and PA directions by up to 4.5%. The maximum gantry angle dependency of the detector was 5.3%. A maximum difference of -3.1% occurred between the measurements and TPS calculations. Conclusion: Results indicate that the new 2D diode detector is stable and useful for QA and end-to-end testing of SRS due to its excellent dose characteristics, high resolution and ease of handling when combined with the StereoPHAN.

The purpose of this study was to propose a verification method and results of intensity-modulated proton therapy (IMPT), using a commercially available heterogeneous phantom. We used a simple simulated head and neck and prostate phantom. An ionization chamber and radiochromic film were used for measurements of absolute dose and relative dose distribution. The measured doses were compared with calculated doses using a treatment planning system. We defined the uncertainty of the measurement point of the ionization chamber due to the effective point of the chamber and mechanical setup error as 2 mm and estimated the dose variation base on a 2 mm error. We prepared a HU-relative stopping power conversion table and fluence correction factor that were specific to the heterogeneous phantom. The fluence correction factor was determined as a function of depth and was obtained from the ratio of the doses in water and in the phantom at the same effective depths. In the simulated prostate plan, composite doses of measurements and calculations agreed within ±1.3% and the maximum local dose differences of each field were 10.0%. Composite doses in the simulated head and neck plan agreed within 4.0% and the maximum local dose difference for each field was 12.0%. The dose difference for each field came within 2% when taking the measurement uncertainty into consideration. In the composite plan, the maximum dose uncertainty was estimated as 4.0% in the simulated prostate plan and 5.8% in the simulated head and neck plan. Film measurements showed good agreement, with more than 92.5% of points passing a gamma value (3%/3 mm). From these results, the heterogeneous phantom should be useful for verification of IMPT by using a phantom-specific HU-relative stopping power conversion, fluence correction factor, and dose error estimation due to the effective point of the chamber.

© 2018 The Authors. Journal of Applied Clinical Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine. The purpose of this study was to provide periodic quality assurance (QA) methods for respiratory-gated proton beam with a range modulation wheel (RMW) and to clarify the characteristics and long-term stability of the respiratory-gated proton beam. A two-dimensional detector array and a solid water phantom were used to measure absolute dose, spread-out Bragg peak (SOBP) width and proton range for monthly QA. SOBP width and proton range were measured using an oblique incidence beam to the lateral side of a solid water phantom and compared between with and without a gating proton beam. To measure the delay time of beam-on/off for annual QA, we collected the beam-on/off signals and the dose monitor-detected pulse. We analyzed the results of monthly QA over a 15-month period and investigated the delay time by machine signal analysis. The dose deviations at proximal, SOBP center and distal points were −0.083 ± 0.25%, 0.026 ± 0.20%, and −0.083 ± 0.35%, respectively. The maximum dose deviation between with and without respiratory gating was −0.95% at the distal point and other deviations were within ±0.5%. Proximal and SOBP center doses showed the same trend over a 15-month period. Delay times of beam-on/off for 200 MeV/SOBP 16 cm were 140.5 ± 0.8 ms and 22.3 ± 13.0 ms, respectively. Delay times for 160 MeV/SOBP 10 cm were 167.5 ± 15.1 ms and 19.1 ± 9.8 ms. Our beam delivery system with the RMW showed sufficient stability for respiratory-gated proton therapy and the system did not show dependency on the energy and the respiratory wave form. The delay times of beam-on/off were within expectations. The proposed QA methods will be useful for managing the quality of respiratory-gated proton beams and other beam delivery systems.

© 2018 American Association of Physicists in Medicine. In this study, we evaluate dosimetric advantages of using patient-specific aperture system with intensity-modulated proton therapy (IMPT) for head and neck tumors at the shallow depth. We used four types of patient-specific aperture system (PSAS) to irradiate shallow regions less than 4 g/cm2 with a sharp lateral penumbra. Ten head and neck IMPT plans with or without aperture were optimized separately with the same 95% prescription dose and same dose constraint for organs at risk (OARs). The plans were compared using dose volume histograms (DVHs), dose distributions, and some dose indexes such as volume receiving 50% of the prescribed dose (V50), mean or maximum dose (Dmean and Dmax) to the OARs. All examples verified in this study had decreased V50 and OAR doses. Average, maximum, and minimum relative reductions of V50 were 15.4%, 38.9%, and 1.0%, respectively. Dmax and Dmean of OARs were decreased by 0.3% to 25.7% and by 1.0% to 46.3%, respectively. The plans with the aperture over more than half of the field showed decreased V50 or OAR dose by more than 10%. The dosimetric advantage of patient-specific apertures with IMPT was clarified in many cases. The PSAS has some dosimetric advantages for clinical use, and in some cases, it enables to fulfill dose constraints.

To measure the absorbed dose to water D w in proton beams using a radiophotoluminescent glass dosimeter (RGD), a method with the correction for the change of the mass stopping power ratio (SPR) and the linear energy transfer (LET) dependence of radiophotoluminescent efficiency [Formula: see text] is proposed. The calibration coefficient in terms of D w for RGDs (GD-302M, Asahi Techno Glass) was obtained using a 60Co γ-ray. The SPR of water to the RGD was calculated by Monte Carlo simulation, and [Formula: see text] was investigated experimentally using a 70 MeV proton beam. For clinical usage, the residual range R res was used as a quality index to determine the correction factor for the beam quality [Formula: see text] and the LET quenching effect of the RGD [Formula: see text]. The proposed method was evaluated by measuring D w at different depths in a 200 MeV proton beam. For both non-modulated and modulated proton beams, [Formula: see text] decreases rapidly where R res is less than 4 cm. The difference in [Formula: see text] between a non-modulated and a modulated proton beam is less than 0.5% for the R res range from 0 cm to 22 cm. [Formula: see text] decreases rapidly at a LET range from 1 to 2 keV µm-1. In the evaluation experiments, D w using RGDs, [Formula: see text] showed good agreement with that obtained using an ionization chamber and the relative difference was within 3% where R res was larger than 1 cm. The uncertainty budget for [Formula: see text] in a proton beam was estimated to investigate the potential of RGD postal dosimetry in proton therapy. These results demonstrate the feasibility of RGD dosimetry in a therapeutic proton beam and the general versatility of the proposed method. In conclusion, the proposed methodology for RGDs in proton dosimetry is applicable where R res  >  1 cm and the RGD is feasible as a postal audit dosimeter for proton therapy.

The purpose of this paper is to describe an outline of a proton therapy system in Nagoya Proton Therapy Center (NPTC). The NPTC has a synchrotron with a linac injector and three treatment rooms: two rooms are equipped with a gantry and the other one is equipped with a fixed horizontal beamline. One gantry treatment room has a pencil beam scanning treatment delivery nozzle. The other two treatment rooms have a passive scattering treatment delivery nozzle. In the scanning treatment delivery nozzle, an energy absorber and an aperture system to treat head and neck cancer have been equipped. In the passive treatment delivery nozzle, a multi-leaf collimator is equipped. We employ respiratory gating to treat lung and liver cancers for passive irradiation. The proton therapy system passed all acceptance tests. The first patient was treated on February 25, 2013, using passive scattering fixed beams. Respiratory gating is commonly used to treat lung and liver cancers in the passive scattering system. The MLCs are our first choice to limit the irradiation field. The use of the aperture for scanning irradiation reduced the lateral fall off by half or less. The energy absorber and aperture system in scanning delivery is beneficial to treat head and neck cancer.

© 2016 Elsevier Inc. Purpose To determine the relative biological effectiveness (RBE), oxygen enhancement ratio (OER), and contribution of the indirect effect of spot scanning proton beams, passive scattering proton beams, or both in cultured cells in comparison with clinically used photons. Methods and Materials The RBE of passive scattering proton beams at the center of the spread-out Bragg peak (SOBP) was determined from dose-survival curves in 4 cell lines using 6-MV X rays as controls. Survival of 2 cell lines after spot scanning and passive scattering proton irradiation was then compared. Biological effects at the distal end region of the SOBP were also investigated. The OER of passive scattering proton beams and 6 MX X rays were investigated in 2 cell lines. The RBE and OER values were estimated at a 10% cell survival level. The maximum degree of protection of radiation effects by dimethyl sulfoxide was determined to estimate the contribution of the indirect effect against DNA damage. All experiments comparing protons and X rays were made under the same biological conditions. Results The RBE values of passive scattering proton beams in the 4 cell lines examined were 1.01 to 1.22 (average, 1.14) and were almost identical to those of spot scanning beams. Biological effects increased at the distal end of the SOBP. In the 2 cell lines examined, the OER was 2.74 (95% confidence interval, 2.56-2.80) and 3.08 (2.84-3.11), respectively, for X rays, and 2.39 (2.38-2.43) and 2.72 (2.69-2.75), respectively, for protons (P<.05 for both cells between X rays and protons). The maximum degree of protection was significantly higher for X rays than for proton beams (P<.05). Conclusions The RBE values of spot scanning and passive scattering proton beams were almost identical. The OER was lower for protons than for X rays. The lower contribution of the indirect effect may partly account for the lower OER of protons.

Purpose: The main purpose in this study was to present the results of beam modeling and how the authors systematically investigated the influence of double and triple Gaussian proton kernel models on the accuracy of dose calculations for spot scanning technique. Methods: The accuracy of calculations was important for treatment planning software (TPS) because the energy, spot position, and absolute dose had to be determined by TPS for the spot scanning technique. The dose distribution was calculated by convolving in-air fluence with the dose kernel. The dose kernel was the in-water 3D dose distribution of an infinitesimal pencil beam and consisted of an integral depth dose (IDD) and a lateral distribution. Accurate modeling of the low-dose region was important for spot scanning technique because the dose distribution was formed by cumulating hundreds or thousands of delivered beams. The authors employed a double Gaussian function as the in-air fluence model of an individual beam. Double and triple Gaussian kernel models were also prepared for comparison. The parameters of the kernel lateral model were derived by fitting a simulated in-water lateral dose profile induced by an infinitesimal proton beam, whose emittance was zero, at various depths using Monte Carlo (MC) simulation. The fitted parameters were interpolated as a function of depth in water and stored as a separate look-up table. These stored parameters for each energy and depth in water were acquired from the look-up table when incorporating them into the TPS. The modeling process for the in-air fluence and IDD was based on the method proposed in the literature. These were derived using MC simulation and measured data. The authors compared the measured and calculated absolute doses at the center of the spread-out Bragg peak (SOBP) under various volumetric irradiation conditions to systematically investigate the influence of the two types of kernel models on the dose calculations. Results: The authors investigated the difference between double and triple Gaussian kernel models. The authors found that the difference between the two studied kernel models appeared at mid-depths and the accuracy of predicting the double Gaussian model deteriorated at the low-dose bump that appeared at mid-depths. When the authors employed the double Gaussian kernel model, the accuracy of calculations for the absolute dose at the center of the SOBP varied with irradiation conditions and the maximum difference was 3.4%. In contrast, the results obtained from calculations with the triple Gaussian kernel model indicated good agreement with the measurements within ?1.1%, regardless of the irradiation conditions. Conclusions: The difference between the results obtained with the two types of studied kernel models was distinct in the high energy region. The accuracy of calculations with the double Gaussian kernel model varied with the field size and SOBP width because the accuracy of prediction with the double Gaussian model was insufficient at the low-dose bump. The evaluation was only qualitative under limited volumetric irradiation conditions. Further accumulation of measured data would be needed to quantitatively comprehend what influence the double and triple Gaussian kernel models had on the accuracy of dose calculations.

© 2015 American Association of Physicists in Medicine. 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/cm2. 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/cm2. 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 mmdistance criteria of the γ-index on average. Conclusions: The patient-specific aperture system improved dose distributions, particularly in shallow-region plans.

To investigate the effects of scattered radiation when a thin titanium plate (thickness, 0.05 cm) used for skull fixation in cerebral nerve surgery is irradiated by a 4-MV photon beam. We investigated the dose distribution of radiation inside a phantom that simulates a human head fitted with a thin titanium plate used for post-surgery skull fixation and compared the distribution data measured using detectors, obtained by Monte Carlo (MC) simulations, and calculated using a radiation treatment planning system (TPS). Simulations were shown to accurately represent measured values. The effects of scattered radiation produced by high-Z materials such as titanium are not sufficiently considered currently in TPS dose calculations. Our comparisons show that the dose distribution is affected by scattered radiation around a thin high-Z material. The depth dose is measured and calculated along the central beam axis inside a water phantom with thin titanium plates at various depths. The maximum relative differences between simulation and TPS results on the entrance and exit sides of the plate were 23.1% and - 12.7%, respectively. However, the depth doses do not change in regions deeper than the plate in water. Although titanium is a high-Z material, if the titanium plate used for skull fixation in cerebral nerve surgery is thin, there is a slight change in the dose distribution in regions away from the plate. In addition, we investigated the effects of variation of photon energies, sizes of radiation field and thickness of the plate. When the target to be irradiated is far from the thin titanium plate, the dose differs little from what it would be in the absence of a plate, though the dose escalation existed in front of the metal plate.

Recently, dose delivery techniques are evolved rapidly in radiotherapy. Intensity-modulated radiotherapy (IMRT) enables a non-uniform complex dose distribution. Therefore, it is necessary to obtain data of two dimensional (2D) dose distribution. Radiographic films are more commonly used in this situation. The radiographic film features high resolution, and obtains 2D data easily. However, the standard method of the radiographic film dosimetry has not been established in radiotherapy. In this study, the radiographic film EDR2 (Kodak) and Gafchromic Film RT-QA (ISP. INC) were used and compared by measurements and Monte Carlo simulations. Both films were digitalized using a scanner (EPSON Offirio ES 100000G) and analyzed by the DD-system (R-TECH. INC). Density-dose conversion curves were plotted and percentage depth dose (PDD) curves were calculated by the DD-system. The simulation geometry was the same as the film measurement arrangement. In comparison of PDD curves for EDR2, the calculated curve was larger than the curve obtained by an ion chamber and the curve measured by the film dosimetry and PDD curves for RT-QA was smaller than the curve by the ion chamber. Thus, the calculated curve was overestimated and the curve by the film dosimetry was underestimated. On the other hand, calculated PDD for RT-QA was smaller than the PDD by the ion chamber and the curve measured by RT-QA agreed with the PDD by the ion chamber.