No Arabic abstract
Purpose: This is a theoretical simulation study for proof of concept of radiochromic film dosimetry to measure physical and biological doses without plan-based quenching correction for patient-specific quality assurance of carbon-ion radiotherapy. Methods: We took a layer-stacking carbon-ion beam comprised of range-shifted beamlets. The dosimeter response was simulated according to an experimental quenching model. The beam model followed a treatment planning system. The beam was decomposed into finely arranged beamlets with weights estimated by deconvolution of longitudinal dosimeter responses. The distributions of physical and biological doses were reconstructed from the estimated weights, and were compared with the plan. We also evaluated the sensitivity to measurement errors and to erratic delivery with an undelivered beamlet. Results: The reconstructed physical and biological doses accurately reproduced the simulated delivery with errors approximately corresponding to the measurement errors. The erratic beam delivery was easily detectable by comparison of biological dose distribution to the plan. Conclusions: We have developed a method to measure physical and biological doses by longitudinal dosimetry of quenched response without using plan data. The method only involves a general optimization algorithm, a radiobiology model, and experimental beamlet data, and requires no extra corrections. Theoretically, this approach is applicable to various dosimeters and to proton and ion beams of any delivery method, regardless of quenching or biological effectiveness.
In carbon-ion radiotherapy, single-beam delivery each day in alternate directions has been commonly practiced for operational efficiency, taking advantage of the Bragg peak and the relative biological effectiveness (RBE) for uniform dose conformation to a tumor. The treatment plans are usually evaluated with total RBE-weighted dose, which is however deficient in relevance to the biological effect in the linear-quadratic model due to its quadratic-dose term, or the dose-fractionation effect. In this study, we reformulate the extrapolated response dose (ERD), or synonymously BED, which normalizes the dose-fractionation and cell-repopulation effects as well as the RBE of treating radiation, based on inactivation of a single model cell system and a typical treating radiation in carbon-ion RT. The ERD distribution virtually represents the biological effect of the treatment regardless of radiation modality or fractionation scheme. We applied the ERD formulation to simplistic model treatments and to a preclinical survey for hypofractionation based on an actual prostate-cancer treatment of carbon-ion radiotherapy. The proposed formulation was demonstrated to be practical and to offer theoretical implications. In the prostate-cancer case, the ERD distribution was very similar to the RBE-weighted-dose distribution of the actual treatment in 12 fractions. With hypofractionation, while the RBE-weighted-dose distribution varied significantly, the ERD distribution was nearly invariant, implying that the carbon-ion radiotherapy would be insensitive to fractionation. However, treatment evaluation with simplistic biological dose is intrinsically limited and must be complemented in practice somehow by clinical experiences and biology experiments.
Purpose: This study aims to optimize and characterize the response of a mPSD for in vivo dosimetry in HDR brachytherapy. Methods: An exhaustive analysis was carried out in order to obtain an optimized mPSD design that maximize the scintillation light collection produced by the interaction of ionizing photons. Several mPSD prototypes were built and tested in order to determine the appropriate order of scintillators relative to the photodetector, as well as their length as a function of the scintillation light emitted. Scintillators BCF-60, BCF-12 and BCF-10 constituted the mPSD sensitive volume.Each scintillator contribution to the total spectrum was determined by irradiations in the low energy range.For the best mPSD design, a numerical optimization was done in order to select the optical components that better match the light emission profile. The optimized dosimetric system was used for HDR brachytherapy dose determination. The system performance was quantified in term of signal to noise ratio and signal to background ratio. Results: It was determined that BCF-60 should be placed at the distal position, BCF-12 in the center and BCF-10 at proximal position with respect to the photodetector.This configuration allowed for optimized light transmission through the collecting fiber, avoiding inter-scintillator excitation and self-absorption effects.The optimized luminescence system allowed for signal deconvolution using a multispectral approach, extracting the dose to each element while taking into account Cerenkov stem effect.Differences between the mPSD measurements and TG-43 remain below 5%. In all measurement conditions, the system was able to properly differentiate the produced scintillation signal from the background one. Conclusions: A mPSD was constructed and optimized for HDR brachytherapy dosimetry, enabling real time dose determination, up to 6.5cm from the 192Ir source.
Depth distributions of positron-emitting nuclei in PMMA phantoms are calculated within a Monte Carlo model for Heavy-Ion Therapy (MCHIT) based on the GEANT4 toolkit (version 8.0). The calculated total production rates of $^{11}$C, $^{10}$C and $^{15}$O nuclei are compared with experimental data and with corresponding results of the FLUKA and POSGEN codes. The distributions of e$^+$ annihilation points are obtained by simulating radioactive decay of unstable nuclei and transporting positrons in surrounding medium. A finite spatial resolution of the Positron Emission Tomography (PET) is taken into account in a simplified way. Depth distributions of $beta^+$-activity as seen by a PET scanner are calculated and compared to available data for PMMA phantoms. The calculated $beta^+$-activity profiles are in good agreement with PET data for proton and $^{12}$C beams at energies suitable for particle therapy. The MCHIT capability to predict the $beta^+$-activity and dose distributions in tissue-like materials of different chemical composition is demonstrated.
This article is intended to supplement our 2015 paper in Medical Physics titled Noise properties of CT images reconstructed by use of constrained total-variation, data-discrepancy minimization, in which ordered subsets methods were employed to perform total-variation constrained data-discrepancy minimization for image reconstruction in X-ray computed tomography. Here we provide details regarding implementation of the ordered subsets algorithms and suggestions for selection of algorithm parameters. Detailed pseudo-code is included for every algorithm implemented in the original manuscript.
MESA (Mainz Energy recovery Superconducting Accelerator) is an energy recovery linac (ERL) which is under construction at Johannes Gutenberg University in Mainz. It will be operated in external beam (EB) mode with 150 $mu$A electron beam at 155 MeV and energy recovery (ER) mode with 1 mA (first stage) and later 10 mA (second stage) electron beam at 105 MeV. An important factor which may limit performance of the machine is a beam breakup (BBU) instability which may occur due to excitation of higher-order modes (HOMs) in superconducting RF cavities. This effect occurs only when the injected beam current exceeds a threshold value. The aim of the present work is to develop a software for reliable determination of the threshold current in MESA, find main factors which may change its value and finally make a decision concerning capability of MESA operation at 10 mA and need for additional measures for suppressing BBU instability.