| Home | << 1 2 3 4 5 6 >> |
|
Assam, I., Vijande, J., Ballester, F., Perez-Calatayud, J., Poppe, B., & Siebert, F. A. (2022). Evaluation of dosimetric effects of metallic artifact reduction and tissue assignment on Monte Carlo dose calculations for I-125 prostate implants. Med. Phys., 49, 6195–6208.
Abstract: Purpose Monte Carlo (MC) simulation studies, aimed at evaluating the magnitude of tissue heterogeneity in I-125 prostate permanent seed implant brachytherapy (BT), customarily use clinical post-implant CT images to generate a virtual representation of a realistic patient model (virtual patient model). Metallic artifact reduction (MAR) techniques and tissue assignment schemes (TAS) are implemented on the post-implant CT images to mollify metallic artifacts due to BT seeds and to assign tissue types to the voxels corresponding to the bright seed spots and streaking artifacts, respectively. The objective of this study is to assess the combined influence of MAR and TAS on MC absorbed dose calculations in post-implant CT-based phantoms. The virtual patient models used for I-125 prostate implant MC absorbed dose calculations in this study are derived from the CT images of an external radiotherapy prostate patient without BT seeds and prostatic calcifications, thus averting the need to implement MAR and TAS. Methods The geometry of the IsoSeed I25.S17plus source is validated by comparing the MC calculated results of the TG-43 parameters for the line source approximation with the TG-43U1S2 consensus data. Four MC absorbed dose calculations are performed in two virtual patient models using the egs_brachy MC code: (1) TG-43-based D-w,w-TG(43), (2) D-w,D-w-MBDC that accounts for interseed scattering and attenuation (ISA), (3) D-m,D-m that examines ISA and tissue heterogeneity by scoring absorbed dose in tissue, and (4) D-w,D-m that unlike D-m,D-m scores absorbed dose in water. The MC absorbed doses (1) and (2) are simulated in a TG-43 patient phantom derived by assigning the densities of every voxel to 1.00 g cm(-3) (water), whereas MC absorbed doses (3) and (4) are scored in the TG-186 patient phantom generated by mapping the mass density of each voxel to tissue according to a CT calibration curve. The MC absorbed doses calculated in this study are compared with VariSeed v8.0 calculated absorbed doses. To evaluate the dosimetric effect of MAR and TAS, the MC absorbed doses of this work (independent of MAR and TAS) are compared to the MC absorbed doses of different I-125 source models from previous studies that were calculated with different MC codes using post-implant CT-based phantoms generated by implementing MAR and TAS on post-implant CT images. Results The very good agreement of TG-43 parameters of this study and the published consensus data within 3% validates the geometry of the IsoSeed I25.S17plus source. For the clinical studies, the TG-43-based calculations show a D-90 overestimation of more than 4% compared to the more realistic MC methods due to ISA and tissue composition. The results of this work generally show few discrepancies with the post-implant CT-based dosimetry studies with respect to the D-90 absorbed dose metric parameter. These discrepancies are mainly Type B uncertainties due to the different I-125 source models and MC codes. Conclusions The implementation of MAR and TAS on post-implant CT images have no dosimetric effect on the I-125 prostate MC absorbed dose calculation in post-implant CT-based phantoms.
|
|
Ma, Y. Z., Vijande, J., Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., et al. (2017). A generic TG-186 shielded applicator for commissioning model-based dose calculation algorithms for high-dose-rate Ir-192 brachytherapy. Med. Phys., 44(11), 5961–5976.
Abstract: PurposeA joint working group was created by the American Association of Physicists in Medicine (AAPM), the European Society for Radiotherapy and Oncology (ESTRO), and the Australasian Brachytherapy Group (ABG) with the charge, among others, to develop a set of well-defined test case plans and perform calculations and comparisons with model-based dose calculation algorithms (MBDCAs). Its main goal is to facilitate a smooth transition from the AAPM Task Group No. 43 (TG-43) dose calculation formalism, widely being used in clinical practice for brachytherapy, to the one proposed by Task Group No. 186 (TG-186) for MBDCAs. To do so, in this work a hypothetical, generic high-dose rate (HDR) Ir-192 shielded applicator has been designed and benchmarked. MethodsA generic HDR Ir-192 shielded applicator was designed based on three commercially available gynecological applicators as well as a virtual cubic water phantom that can be imported into any DICOM-RT compatible treatment planning system (TPS). The absorbed dose distribution around the applicator with the TG-186 Ir-192 source located at one dwell position at its center was computed using two commercial TPSs incorporating MBDCAs (Oncentra((R)) Brachy with Advanced Collapsed-cone Engine, ACE, and BrachyVision ACUROS) and state-of-the-art Monte Carlo (MC) codes, including ALGEBRA, BrachyDose, egs_brachy, Geant4, MCNP6, and Penelope2008. TPS-based volumetric dose distributions for the previously reported source centered in water and source displaced test cases, and the new source centered in applicator test case, were analyzed here using the MCNP6 dose distribution as a reference. Volumetric dose comparisons of TPS results against results for the other MC codes were also performed. Distributions of local and global dose difference ratios are reported. ResultsThe local dose differences among MC codes are comparable to the statistical uncertainties of the reference datasets for the source centered in water and source displaced test cases and for the clinically relevant part of the unshielded volume in the source centered in applicator case. Larger local differences appear in the shielded volume or at large distances. Considering clinically relevant regions, global dose differences are smaller than the local ones. The most disadvantageous case for the MBDCAs is the one including the shielded applicator. In this case, ACUROS agrees with MC within [-4.2%, +4.2%] for the majority of voxels (95%) while presenting dose differences within [-0.12%, +0.12%] of the dose at a clinically relevant reference point. For ACE, 95% of the total volume presents differences with respect to MC in the range [-1.7%, +0.4%] of the dose at the reference point. ConclusionsThe combination of the generic source and generic shielded applicator, together with the previously developed test cases and reference datasets (available in the Brachytherapy Source Registry), lay a solid foundation in supporting uniform commissioning procedures and direct comparisons among treatment planning systems for HDR Ir-192 brachytherapy.
|
|
Zhang, X., Xiao, Y. T., & Gimeno, B. (2020). Multipactor Suppression by a Resonant Static Magnetic Field on a Dielectric Surface. IEEE Trans. Electron Devices, 67(12), 5723–5728.
Abstract: In this article, we study the suppression of the multipactor phenomenon on a dielectric surface by a resonant static magnetic field. A homemade Monte Carlo algorithm is developed for multipactor simulations on a dielectric surface driven by two orthogonal radio frequency (RF) electric field components. When the static magnetic field is perpendicular to the tangential and normal RF electric fields, it is shown that if the normal electric field lags the tangential electric field by pi/2, the superposition of the normal and tangential electric fields will trigger a gyro-acceleration of the electron cloud and restrain the multipactor discharge effectively. By contrast, when the normal electric field is in advance of the tangential electric field by pi/2, the difference between the normal and tangential electric fields drives gyro-motion of the electron cloud. Consequently, two enhanced discharge zones are inevitable. The suppression effects of the resonant static magnetic field that is parallel to the tangential RF electric field or to the normal RF electric field are also presented.
|
|
Aguiar, P., Rafecas, M., Ortuño, J. E., Kontaxakis, G., Santos, A., Pavia, J., et al. (2010). Geometrical and Monte Carlo projectors in 3D PET reconstruction. Med. Phys., 37(11), 5691–5702.
Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.
|
|
Valdes-Cortez, C., Niatsetski, Y., Perez-Calatayud, J., Ballester, F., & Vijande, J. (2022). A Monte Carlo study of the relative biological effectiveness in surface brachytherapy. Med. Phys., 49, 5576–5588.
Abstract: Purpose This work aims to simulate clustered DNA damage from ionizing radiation and estimate the relative biological effectiveness (RBE) for radionuclide (rBT)- and electronic (eBT)-based surface brachytherapy through a hybrid Monte Carlo (MC) approach, using realistic models of the sources and applicators. Methods Damage from ionizing radiation has been studied using the Monte Carlo Damage Simulation algorithm using as input the primary electron fluence simulated using a state-of-the-art MC code, PENELOPE-2018. Two Ir-192 rBT applicators, Valencia and Leipzig, one Co-60 source with a Freiburg Flap applicator (reference source), and two eBT systems, Esteya and INTRABEAM, have been included in this study implementing full realizations of their geometries as disclosed by the manufacturer. The role played by filtration and tube kilovoltage has also been addressed. Results For rBT, an RBE value of about 1.01 has been found for the applicators and phantoms considered. In the case of eBT, RBE values for the Esteya system show an almost constant RBE value of about 1.06 for all depths and materials. For INTRABEAM, variations in the range of 1.12-1.06 are reported depending on phantom composition and depth. Modifications in the Esteya system, filtration, and tube kilovoltage give rise to variations in the same range. Conclusions Current clinical practice does not incorporate biological effects in surface brachytherapy. Therefore, the same absorbed dose is administered to the patients independently on the particularities of the rBT or eBT system considered. The almost constant RBE values reported for rBT support that assumption regardless of the details of the patient geometry, the presence of a flattening filter in the applicator design, or even significant modifications in the photon energy spectra above 300 keV. That is not the case for eBT, where a clear dependence on the eBT system and the characteristics of the patient geometry are reported. A complete study specific for each eBT system, including detailed applicator characteristics (size, shape, filtering, among others) and common anatomical locations, should be performed before adopting an existing RBE value.
|
|
Peppa, V., Thomson, R. M., Enger, S. A., Fonseca, G. P., Lee, C. N., Lucero, J. N. E., et al. (2023). A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy. Med. Phys., 50(7), 4675–4687.
Abstract: PurposeTo provide the first clinical test case for commissioning of Ir-192 brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow. Acquisition and Validation MethodsA computational patient phantom model was generated from a clinical multi-catheter Ir-192 HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic Ir-192 HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution. Data Format and Usage NotesThe dataset is available online at ,. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations. Potential ApplicationsThe dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.
|
|
Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., Mourtada, F., Fonseca, G. P., et al. (2015). A generic high-dose rate Ir-192 brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism. Med. Phys., 42(6), 3048–3062.
Abstract: Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) Ir-192 source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR Ir-192 source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic Ir-192 source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra (R) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision AcuRos (TM)]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and pENELopE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR Ir-192 source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 +/- 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR Ir-192 source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.
|
|
Granero, D., Candela-Juan, C., Vijande, J., Ballester, F., Perez-Calatayud, J., Jacob, D., et al. (2016). Technical Note: Dosimetry of Leipzig and Valencia applicators without the plastic cap. Med. Phys., 43(5), 2087–4pp.
Abstract: Purpose: High dose rate (HDR) brachytherapy for treatment of small skin lesions using the Leipzig and Valencia applicators is a widely used technique. These applicators are equipped with an attachable plastic cap to be placed during fraction delivery to ensure electronic equilibrium and to prevent secondary electrons from reaching the skin surface. The purpose of this study is to report on the dosimetric impact of the cap being absent during HDR fraction delivery, which has not been explored previously in the literature. Methods: GEANT4 Monte Carlo simulations (version 10.0) have been performed for the Leipzig and Valencia applicators with and without the plastic cap. In order to validate the Monte Carlo simulations, experimental measurements using radiochromic films have been done. Results: Dose absorbed within 1 mm of the skin surface increases by a factor of 1500% for the Leipzig applicators and of 180% for the Valencia applicators. Deeper than 1 mm, the overdosage flattens up to a 10% increase. Conclusions: Differences of treating with or without the plastic cap are significant. Users must check always that the plastic cap is in place before any treatment in order to avoid overdosage of the skin. Prior to skin HDR fraction delivery, the timeout checklist should include verification of the cap placement. (C) 2016 American Association of Physicists in Medicine.
|
|
Ortega, P. G., Torres-Espallardo, I., Cerutti, F., Ferrari, A., Gillam, J. E., Lacasta, C., et al. (2015). Noise evaluation of Compton camera imaging for proton therapy. Phys. Med. Biol., 60(5), 1845–1863.
Abstract: Compton Cameras emerged as an alternative for real-time dose monitoring techniques for Particle Therapy (PT), based on the detection of prompt-gammas. As a consequence of the Compton scattering process, the gamma origin point can be restricted onto the surface of a cone (Compton cone). Through image reconstruction techniques, the distribution of the gamma emitters can be estimated, using cone-surfaces backprojections of the Compton cones through the image space, along with more sophisticated statistical methods to improve the image quality. To calculate the Compton cone required for image reconstruction, either two interactions, the last being photoelectric absorption, or three scatter interactions are needed. Because of the high energy of the photons in PT the first option might not be adequate, as the photon is not absorbed in general. However, the second option is less efficient. That is the reason to resort to spectral reconstructions, where the incoming. energy is considered as a variable in the reconstruction inverse problem. Jointly with prompt gamma, secondary neutrons and scattered photons, not strongly correlated with the dose map, can also reach the imaging detector and produce false events. These events deteriorate the image quality. Also, high intensity beams can produce particle accumulation in the camera, which lead to an increase of random coincidences, meaning events which gather measurements from different incoming particles. The noise scenario is expected to be different if double or triple events are used, and consequently, the reconstructed images can be affected differently by spurious data. The aim of the present work is to study the effect of false events in the reconstructed image, evaluating their impact in the determination of the beam particle ranges. A simulation study that includes misidentified events (neutrons and random coincidences) in the final image of a Compton Telescope for PT monitoring is presented. The complete chain of detection, from the beam particle entering a phantom to the event classification, is simulated using FLUKA. The range determination is later estimated from the reconstructed image obtained from a two and three-event algorithm based on Maximum Likelihood Expectation Maximization. The neutron background and random coincidences due to a therapeutic-like time structure are analyzed for mono-energetic proton beams. The time structure of the beam is included in the simulations, which will affect the rate of particles entering the detector.
Keywords: proton therapy; Compton camera; Monte Carlo methods; FLUKA; prompt gamma; range verification; MLEM
|
|
n_TOF Collaboration(Mendoza, E. et al), Giubrone, G., & Tain, J. L. (2011). Improved Neutron Capture Cross Section Measurements with the n_TOF Total Absorption Calorimeter. J. Korean Phys. Soc., 59(2), 1813–1816.
Abstract: The n_TOF collaboration operates a Total Absorption Calorimeter (TAC) [1] for measuring neutron capture cross-sections of low-mass and/or radioactive samples. The results obtained with the TAC have led to a substantial improvement of the capture cross sections of (237)Np and (240)Pu [2]. The experience acquired during the first measurements has allowed us to optimize the performance of the TAC and to improve the capture signal to background ratio, thus opening the way to more complex and demanding measurements on rare radioactive materials. The new design has been reached by a series of detailed Monte Carlo simulations of complete experiments and dedicated test measurements. The new capture setup will be presented and the main achievements highlighted.
|