2020
Carroll, Liam; Croteau, Etienne; Kertzscher, Gustavo; Sarrhini, Otman; Turgeon, Vincent; Lecomte, Roger; Enger, Shirin A.
In: Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB), vol. 76, pp. 92–99, 2020, ISSN: 1724-191X.
Abstract | Links | BibTeX | Tags: Algorithms, Arterial input function, Arteries, Dynamic PET, Electrons, Humans, Imaging, Non-invasive detector development, Phantoms, Positron-Emission Tomography, Scintillation
@article{carroll_cross-validation_2020,
title = {Cross-validation of a non-invasive positron detector to measure the arterial input function for pharmacokinetic modelling in dynamic positron emission tomography},
author = {Liam Carroll and Etienne Croteau and Gustavo Kertzscher and Otman Sarrhini and Vincent Turgeon and Roger Lecomte and Shirin A. Enger},
doi = {10.1016/j.ejmp.2020.06.009},
issn = {1724-191X},
year = {2020},
date = {2020-08-01},
journal = {Physica medica: PM: an international journal devoted to the applications of physics to medicine and biology: official journal of the Italian Association of Biomedical Physics (AIFB)},
volume = {76},
pages = {92--99},
abstract = {Kinetic modeling of positron emission tomography (PET) data can assess index rate of uptake, metabolism and predict disease progression more accurately than conventional static PET. However, it requires knowledge of the time-course of the arterial blood radioactivity concentration, called the arterial input function (AIF). The gold standard to acquire the AIF is by invasive means. The purpose of this study was to validate a previously developed dual readout scintillating fiber-based non-invasive positron detector, hereinafter called non-invasive detector (NID), developed to determine the AIF for dynamic PET measured from the human radial artery. The NID consisted of a 3 m long plastic scintillating fiber with each end coupled to a 5 m long transmission fiber followed by a silicon photomultiplier. The scintillating fiber was enclosed inside the grooves of a plastic cylindrical shell. Two sets of experiments were performed to test the NID against a previously validated microfluidic positron detector. A closed-loop microfluidic system combined with a wrist phantom was used. During the first experiment, the three PET radioisotopes 18F, 11C and 68Ga were tested. After optimizing the detector, a second series of tests were performed using only 18F and 11C. The maximum pulse amplitude to electronic noise ratio was 52 obtained with 11C. Linear regressions showed a linear relation between the two detectors. These preliminary results show that the NID can accurately detect positrons from a patient's wrist and has the potential to non-invasively measure the AIF during a dynamic PET scan. The accuracy of these measurements needs to be determined.},
keywords = {Algorithms, Arterial input function, Arteries, Dynamic PET, Electrons, Humans, Imaging, Non-invasive detector development, Phantoms, Positron-Emission Tomography, Scintillation},
pubstate = {published},
tppubtype = {article}
}
2018
Famulari, Gabriel; Renaud, Marc-André; Poole, Christopher M.; Evans, Michael D. C.; Seuntjens, Jan; Enger, Shirin A.
RapidBrachyMCTPS: a Monte Carlo-based treatment planning system for brachytherapy applications Journal Article
In: Physics in Medicine and Biology, vol. 63, no. 17, pp. 175007, 2018, ISSN: 1361-6560.
Abstract | Links | BibTeX | Tags: Brachytherapy, Computer-Assisted, Humans, Imaging, Monte Carlo Method, Phantoms, Radiotherapy Dosage, Radiotherapy Planning, Software
@article{famulari_rapidbrachymctps_2018,
title = {RapidBrachyMCTPS: a Monte Carlo-based treatment planning system for brachytherapy applications},
author = {Gabriel Famulari and Marc-André Renaud and Christopher M. Poole and Michael D. C. Evans and Jan Seuntjens and Shirin A. Enger},
doi = {10.1088/1361-6560/aad97a},
issn = {1361-6560},
year = {2018},
date = {2018-08-01},
journal = {Physics in Medicine and Biology},
volume = {63},
number = {17},
pages = {175007},
abstract = {Despite being considered the gold standard for brachytherapy dosimetry, Monte Carlo (MC) has yet to be implemented into a software for brachytherapy treatment planning. The purpose of this work is to present RapidBrachyMCTPS, a novel treatment planning system (TPS) for brachytherapy applications equipped with a graphical user interface (GUI), optimization tools and a Geant4-based MC dose calculation engine, RapidBrachyMC. Brachytherapy sources and applicators were implemented in RapidBrachyMC and made available to the user via a source and applicator library in the GUI. To benchmark RapidBrachyMC, TG-43 parameters were calculated for the microSelectron v2 (192Ir) and SelectSeed (125I) source models and were compared against previously validated MC brachytherapy codes. The performance of RapidBrachyMC was evaluated for a prostate high dose rate case. To assess the accuracy of RapidBrachyMC in a heterogeneous setup, dose distributions with a cylindrical shielded/unshielded applicator were validated against film measurements in a Solid WaterTM phantom. TG-43 parameters calculated using RapidBrachyMC generally agreed within 1%-2% of the results obtained in previously published work. For the prostate case, clinical dosimetric indices showed general agreement with Oncentra TPS within 1%. Simulation times were on the order of minutes on a single core to achieve uncertainties below 2% in voxels within the prostate. The calculation time was decreased further using the multithreading features of Geant4. In the comparison between MC-calculated and film-measured dose distributions, at least 95% of points passed the 3%/3 mm gamma index criteria in all but one case. RapidBrachyMCTPS can be used as a post-implant dosimetry toolkit, as well as for MC-based brachytherapy treatment planning. This software is especially well suited for the development of new source and applicator models.},
keywords = {Brachytherapy, Computer-Assisted, Humans, Imaging, Monte Carlo Method, Phantoms, Radiotherapy Dosage, Radiotherapy Planning, Software},
pubstate = {published},
tppubtype = {article}
}
Famulari, Gabriel; Pater, Piotr; Enger, Shirin A.
Microdosimetric Evaluation of Current and Alternative Brachytherapy Sources-A Geant4-DNA Simulation Study Journal Article
In: International Journal of Radiation Oncology, Biology, Physics, vol. 100, no. 1, pp. 270–277, 2018, ISSN: 1879-355X.
Abstract | Links | BibTeX | Tags: Brachytherapy, Gadolinium, Imaging, Iodine Radioisotopes, Iridium Radioisotopes, Linear Energy Transfer, Monte Carlo Method, Phantoms, Radioisotopes, Radiometry, Radiotherapy Dosage, Relative Biological Effectiveness, Selenium Radioisotopes, Ytterbium
@article{famulari_microdosimetric_2018,
title = {Microdosimetric Evaluation of Current and Alternative Brachytherapy Sources-A Geant4-DNA Simulation Study},
author = {Gabriel Famulari and Piotr Pater and Shirin A. Enger},
doi = {10.1016/j.ijrobp.2017.09.040},
issn = {1879-355X},
year = {2018},
date = {2018-01-01},
journal = {International Journal of Radiation Oncology, Biology, Physics},
volume = {100},
number = {1},
pages = {270--277},
abstract = {PURPOSE: Radioisotopes such as 75Se, 169Yb, and 153Gd have photon energy spectra and half-lives that make them excellent candidates as alternatives to 192Ir for high-dose-rate brachytherapy. The aim of the present study was to evaluate the relative biological effectiveness (RBE) of current (192Ir, 125I, 103Pd) and alternative (75Se, 169Yb, 153Gd) brachytherapy radionuclides using Monte Carlo simulations of lineal energy distributions.
METHODS AND MATERIALS: Brachytherapy sources (microSelectron v2 [192Ir, 75Se, 169Yb, 153Gd], SelectSeed [125I], and TheraSeed [103Pd]) were placed in the center of a spherical water phantom with a radius of 40 cm using the Geant4 Monte Carlo simulation toolkit. The kinetic energy of all primary, scattered, and fluorescence photons interacting in a scoring volume were tallied at various depths from the source. Electron tracks were generated by sampling the photon interaction spectrum and tracking all the interactions down to 10 eV using the event-by-event capabilities of the Geant4-DNA models. The dose mean lineal energy (y¯D) values were obtained through random sampling of transfer points and overlaying spherical scoring volumes within the associated volume of the tracks. The scoring volume diameter was determined by fitting the y¯D ratio for 125I to its observed RBE.
RESULTS: y¯D increased with the increasing distance from the source for 192Ir, 75Se, and 169Yb, remained constant for 153Gd and 125I, and decreased for 103Pd. The diameter at which the y¯D ratio coincided with the RBE of 1.15 to 1.20 for 125I was ∼25 to 40 nm. The RBE (reference 1 MeV photons) at high doses and dose rates for 192Ir, 75Se, 169Yb, 153Gd, 125I, and 103Pd was 1.028 to 1.034, 1.05 to 1.07, 1.12 to 1.15, 1.16 to 1.21, 1.15 to 1.20, and 1.17 to 1.22, respectively.
CONCLUSIONS: The radiation quality of the radionuclides under investigation was greater than that of high-energy photons. The present study has provided a set of values to modify the prescription doses for brachytherapy to account for the variation in radiation quality among radionuclides.},
keywords = {Brachytherapy, Gadolinium, Imaging, Iodine Radioisotopes, Iridium Radioisotopes, Linear Energy Transfer, Monte Carlo Method, Phantoms, Radioisotopes, Radiometry, Radiotherapy Dosage, Relative Biological Effectiveness, Selenium Radioisotopes, Ytterbium},
pubstate = {published},
tppubtype = {article}
}
METHODS AND MATERIALS: Brachytherapy sources (microSelectron v2 [192Ir, 75Se, 169Yb, 153Gd], SelectSeed [125I], and TheraSeed [103Pd]) were placed in the center of a spherical water phantom with a radius of 40 cm using the Geant4 Monte Carlo simulation toolkit. The kinetic energy of all primary, scattered, and fluorescence photons interacting in a scoring volume were tallied at various depths from the source. Electron tracks were generated by sampling the photon interaction spectrum and tracking all the interactions down to 10 eV using the event-by-event capabilities of the Geant4-DNA models. The dose mean lineal energy (y¯D) values were obtained through random sampling of transfer points and overlaying spherical scoring volumes within the associated volume of the tracks. The scoring volume diameter was determined by fitting the y¯D ratio for 125I to its observed RBE.
RESULTS: y¯D increased with the increasing distance from the source for 192Ir, 75Se, and 169Yb, remained constant for 153Gd and 125I, and decreased for 103Pd. The diameter at which the y¯D ratio coincided with the RBE of 1.15 to 1.20 for 125I was ∼25 to 40 nm. The RBE (reference 1 MeV photons) at high doses and dose rates for 192Ir, 75Se, 169Yb, 153Gd, 125I, and 103Pd was 1.028 to 1.034, 1.05 to 1.07, 1.12 to 1.15, 1.16 to 1.21, 1.15 to 1.20, and 1.17 to 1.22, respectively.
CONCLUSIONS: The radiation quality of the radionuclides under investigation was greater than that of high-energy photons. The present study has provided a set of values to modify the prescription doses for brachytherapy to account for the variation in radiation quality among radionuclides.
2017
Famulari, Gabriel; Pater, Piotr; Enger, Shirin A.
Microdosimetry calculations for monoenergetic electrons using Geant4-DNA combined with a weighted track sampling algorithm Journal Article
In: Physics in Medicine and Biology, vol. 62, no. 13, pp. 5495–5508, 2017, ISSN: 1361-6560.
Abstract | Links | BibTeX | Tags: Algorithms, DNA, DNA Damage, Electrons, Imaging, Monte Carlo Method, Phantoms, Photons, Radiometry
@article{famulari_microdosimetry_2017,
title = {Microdosimetry calculations for monoenergetic electrons using Geant4-DNA combined with a weighted track sampling algorithm},
author = {Gabriel Famulari and Piotr Pater and Shirin A. Enger},
doi = {10.1088/1361-6560/aa71f6},
issn = {1361-6560},
year = {2017},
date = {2017-07-01},
journal = {Physics in Medicine and Biology},
volume = {62},
number = {13},
pages = {5495--5508},
abstract = {The aim of this study was to calculate microdosimetric distributions for low energy electrons simulated using the Monte Carlo track structure code Geant4-DNA. Tracks for monoenergetic electrons with kinetic energies ranging from 100 eV to 1 MeV were simulated in an infinite spherical water phantom using the Geant4-DNA extension included in Geant4 toolkit version 10.2 (patch 02). The microdosimetric distributions were obtained through random sampling of transfer points and overlaying scoring volumes within the associated volume of the tracks. Relative frequency distributions of energy deposition f(textgreaterE)/f(textgreater0) and dose mean lineal energy ([Formula: see text]) values were calculated in nanometer-sized spherical and cylindrical targets. The effects of scoring volume and scoring techniques were examined. The results were compared with published data generated using MOCA8B and KURBUC. Geant4-DNA produces a lower frequency of higher energy deposits than MOCA8B. The [Formula: see text] values calculated with Geant4-DNA are smaller than those calculated using MOCA8B and KURBUC. The differences are mainly due to the lower ionization and excitation cross sections of Geant4-DNA for low energy electrons. To a lesser extent, discrepancies can also be attributed to the implementation in this study of a new and fast scoring technique that differs from that used in previous studies. For the same mean chord length ([Formula: see text]), the [Formula: see text] calculated in cylindrical volumes are larger than those calculated in spherical volumes. The discrepancies due to cross sections and scoring geometries increase with decreasing scoring site dimensions. A new set of [Formula: see text] values has been presented for monoenergetic electrons using a fast track sampling algorithm and the most recent physics models implemented in Geant4-DNA. This dataset can be combined with primary electron spectra to predict the radiation quality of photon and electron beams.},
keywords = {Algorithms, DNA, DNA Damage, Electrons, Imaging, Monte Carlo Method, Phantoms, Photons, Radiometry},
pubstate = {published},
tppubtype = {article}
}
2011
Xu, Chen; Verhaegen, Frank; Laurendeau, Denis; Enger, Shirin A.; Beaulieu, Luc
An algorithm for efficient metal artifact reductions in permanent seed Journal Article
In: Medical Physics, vol. 38, no. 1, pp. 47–56, 2011, ISSN: 0094-2405.
Abstract | Links | BibTeX | Tags: Algorithms, Artifacts, Brachytherapy, Humans, Imaging, Metals, Monte Carlo Method, Phantoms, Tomography, X-Ray Computed
@article{xu_algorithm_2011,
title = {An algorithm for efficient metal artifact reductions in permanent seed},
author = {Chen Xu and Frank Verhaegen and Denis Laurendeau and Shirin A. Enger and Luc Beaulieu},
doi = {10.1118/1.3519988},
issn = {0094-2405},
year = {2011},
date = {2011-01-01},
journal = {Medical Physics},
volume = {38},
number = {1},
pages = {47--56},
abstract = {PURPOSE: In permanent seed implants, 60 to more than 100 small metal capsules are inserted in the prostate, creating artifacts in x-ray computed tomography (CT) imaging. The goal of this work is to develop an automatic method for metal artifact reduction (MAR) from small objects such as brachytherapy seeds for clinical applications.
METHODS: The approach for MAR is based on the interpolation of missing projections by directly using raw helical CT data (sinogram). First, an initial image is reconstructed from the raw CT data. Then, the metal objects segmented from the reconstructed image are reprojected back into the sinogram space to produce a metal-only sinogram. The Steger method is used to determine precisely the position and edges of the seed traces in the raw CT data. By combining the use of Steger detection and reprojections, the missing projections are detected and replaced by interpolation of non-missing neighboring projections.
RESULTS: In both phantom experiments and patient studies, the missing projections have been detected successfully and the artifacts caused by metallic objects have been substantially reduced. The performance of the algorithm has been quantified by comparing the uniformity between the uncorrected and the corrected phantom images. The results of the artifact reduction algorithm are indistinguishable from the true background value.
CONCLUSIONS: An efficient algorithm for MAR in seed brachytherapy was developed. The test results obtained using raw helical CT data for both phantom and clinical cases have demonstrated that the proposed MAR method is capable of accurately detecting and correcting artifacts caused by a large number of very small metal objects (seeds) in sinogram space. This should enable a more accurate use of advanced brachytherapy dose calculations, such as Monte Carlo simulations.},
keywords = {Algorithms, Artifacts, Brachytherapy, Humans, Imaging, Metals, Monte Carlo Method, Phantoms, Tomography, X-Ray Computed},
pubstate = {published},
tppubtype = {article}
}
METHODS: The approach for MAR is based on the interpolation of missing projections by directly using raw helical CT data (sinogram). First, an initial image is reconstructed from the raw CT data. Then, the metal objects segmented from the reconstructed image are reprojected back into the sinogram space to produce a metal-only sinogram. The Steger method is used to determine precisely the position and edges of the seed traces in the raw CT data. By combining the use of Steger detection and reprojections, the missing projections are detected and replaced by interpolation of non-missing neighboring projections.
RESULTS: In both phantom experiments and patient studies, the missing projections have been detected successfully and the artifacts caused by metallic objects have been substantially reduced. The performance of the algorithm has been quantified by comparing the uniformity between the uncorrected and the corrected phantom images. The results of the artifact reduction algorithm are indistinguishable from the true background value.
CONCLUSIONS: An efficient algorithm for MAR in seed brachytherapy was developed. The test results obtained using raw helical CT data for both phantom and clinical cases have demonstrated that the proposed MAR method is capable of accurately detecting and correcting artifacts caused by a large number of very small metal objects (seeds) in sinogram space. This should enable a more accurate use of advanced brachytherapy dose calculations, such as Monte Carlo simulations.
2006
Enger, Shirin A.; af Rosenschöld, Per Munck; Rezaei, Arash; Lundqvist, Hans
Monte Carlo calculations of thermal neutron capture in gadolinium: a comparison of GEANT4 and MCNP with measurements Journal Article
In: Medical Physics, vol. 33, no. 2, pp. 337–341, 2006, ISSN: 0094-2405.
Abstract | Links | BibTeX | Tags: Computer-Assisted, Fast Neutrons, Gadolinium, Humans, Imaging, Monte Carlo Method, Neutron Capture Therapy, Phantoms, Radiologic, Radiometry, Radiotherapy Planning, Reproducibility of Results, Technology
@article{enger_monte_2006,
title = {Monte Carlo calculations of thermal neutron capture in gadolinium: a comparison of GEANT4 and MCNP with measurements},
author = {Shirin A. Enger and Per Munck af Rosenschöld and Arash Rezaei and Hans Lundqvist},
doi = {10.1118/1.2150787},
issn = {0094-2405},
year = {2006},
date = {2006-02-01},
journal = {Medical Physics},
volume = {33},
number = {2},
pages = {337--341},
abstract = {GEANT4 is a Monte Carlo code originally implemented for high-energy physics applications and is well known for particle transport at high energies. The capacity of GEANT4 to simulate neutron transport in the thermal energy region is not equally well known. The aim of this article is to compare MCNP, a code commonly used in low energy neutron transport calculations and GEANT4 with experimental results and select the suitable code for gadolinium neutron capture applications. To account for the thermal neutron scattering from chemically bound atoms [S(alpha,beta)] in biological materials a comparison of thermal neutron fluence in tissue-like poly(methylmethacrylate) phantom is made with MCNP4B, GEANT4 6.0 patch1, and measurements from the neutron capture therapy (NCT) facility at the Studsvik, Sweden. The fluence measurements agreed with MCNP calculated results considering S(alpha,beta). The location of the thermal neutron peak calculated with MCNP without S(alpha,beta) and GEANT4 is shifted by about 0.5 cm towards a shallower depth and is 25%-30% lower in amplitude. Dose distribution from the gadolinium neutron capture reaction is then simulated by MCNP and compared with measured data. The simulations made by MCNP agree well with experimental results. As long as thermal neutron scattering from chemically bound atoms are not included in GEANT4 it is not suitable for NCT applications.},
keywords = {Computer-Assisted, Fast Neutrons, Gadolinium, Humans, Imaging, Monte Carlo Method, Neutron Capture Therapy, Phantoms, Radiologic, Radiometry, Radiotherapy Planning, Reproducibility of Results, Technology},
pubstate = {published},
tppubtype = {article}
}