Publications

@article{morcos_impact_2021,

title = {On the impact of absorbed dose specification, tissue heterogeneities, and applicator heterogeneities on Monte Carlo-based dosimetry of Ir-192, Se-75, and Yb-169 in conventional and intensity-modulated brachytherapy for the treatment of cervical cancer},

author = {Marc Morcos and Akila N. Viswanathan and Shirin A. Enger},

doi = {10.1002/mp.14802},

issn = {2473-4209},

year  = {2021},

date = {2021-05-01},

journal = {Medical Physics},

volume = {48},

number = {5},

pages = {2604--2613},

abstract = {PURPOSE: The purpose of this study was to evaluate the impact of dose reporting schemes and tissue/applicator heterogeneities for 192 Ir-, 75 Se-, and 169 Yb-based MRI-guided conventional and intensity-modulated brachytherapy. METHODS AND MATERIALS: Treatment plans using a variety of dose reporting and tissue/applicator segmentation schemes were generated for a cohort (n = 10) of cervical cancer patients treated with 192 Ir-based Venezia brachytherapy. Dose calculations were performed using RapidBrachyMCTPS, a Geant4-based research Monte Carlo treatment planning system. Ultimately, five dose calculation scenarios were evaluated: (a) dose to water in water (Dw,w ); (b) Dw,w taking the applicator material into consideration (Dw,wApp ); (c) dose to water in medium (Dw,m ); (d and e) dose to medium in medium with mass densities assigned either nominally per structure (Dm,m (Nom) ) or voxel-by-voxel (Dm,m ). 

RESULTS: Ignoring the plastic Venezia applicator (Dw,wApp ) overestimates Dm,m by up to 1% (average) with high energy source (192 Ir and 75 Se) and up to 2% with 169 Yb. Scoring dose to water (Dw,wApp or Dw,m ) generally overestimates dose and this effect increases with decreasing photon energy. Reporting dose other than Dm,m (or Dm,m Nom ) for 169 Yb-based conventional and intensity-modulated brachytherapy leads to a simultaneous overestimation (up to 4%) of CTVHR D90 and underestimation (up to 2%) of bladder D2cc due to a significant dip in the mass-energy absorption ratios at the depths of nearby targets and OARs. Using a nominal mass-density assignment per structure, rather than a CT-derived voxel-by-voxel assignment for MRI-guided brachytherapy, amounts to a dose error up to 1% for all radionuclides considered. 

CONCLUSIONS: The effects of the considered dose reporting schemes trend correspondingly between conventional and intensity-modulated brachytherapy. In the absence of CT-derived mass densities, MRI-only-based dosimetry can adequately approximate Dm,m by assigning nominal mass densities to structures. Tissue and applicator heterogeneities do not significantly impact dosimetry for 192 Ir and 75 Se, but do for 169 Yb; dose reporting must be explicitly defined since Dw,m and Dw,w may overstate the dosimetric benefits.},

keywords = {Brachytherapy, Computer-Assisted, dynamic shield brachytherapy, Female, Humans, IMBT, Intensity modulated brachytherapy, Iridium Radioisotopes, Monte Carlo Method, MR-guided brachytherapy, Radiotherapy Dosage, Radiotherapy Planning, rotating shield brachytherapy, RSBT, Uterine Cervical Neoplasms},

pubstate = {published},

tppubtype = {article}

}

@article{morcos_novel_2021,

title = {A novel minimally invasive dynamic-shield, intensity-modulated brachytherapy system for the treatment of cervical cancer},

author = {Marc Morcos and Majd Antaki and Akila N. Viswanathan and Shirin A. Enger},

doi = {10.1002/mp.14459},

issn = {2473-4209},

year  = {2021},

date = {2021-01-01},

journal = {Medical Physics},

volume = {48},

number = {1},

pages = {71--79},

abstract = {PURPOSE: To present a novel, MRI-compatible dynamicshield intensity modulated brachytherapy (IMBT) applicator and delivery system using 192 Ir, 75 Se, and 169 Yb radioisotopes for the treatment of locally advanced cervical cancer. Needle-free IMBT is a promising technique for improving target coverage and organs at risk (OAR) sparing. 

METHODS AND MATERIALS: The IMBT delivery system dynamically controls the rotation of a novel tungsten shield placed inside an MRI-compatible, 6-mm wide intrauterine tandem. Using 36 cervical cancer cases, conventional intracavitary brachytherapy (IC-BT) and intracavitary/interstitial brachytherapy (IC/IS-BT) (10Ci 192 Ir) plans were compared to IMBT (10Ci 192 Ir; 11.5Ci 75 Se; 44Ci 169 Yb). All plans were generated using the Geant4-based Monte Carlo dose calculation engine, RapidBrachyMC. Treatment plans were optimized then normalized to the same high-risk clinical target volume (HR-CTV) D90 and the D2cc for bladder, rectum, and sigmoid in the research brachytherapy planning system, RapidBrachyMCTPS. Plans were renormalized until either of the three OAR reached dose limits to calculate the maximum achievable HR-CTV D90 and D98 . RESULTS: Compared to IC-BT, IMBT with either of the three radionuclides significantly improves the HR-CTV D90 and D98 by up to 5.2% ± 0.3% (P textless 0.001) and 6.7% ± 0.5% (P textless 0.001), respectively, with the largest dosimetric enhancement when using 169 Yb followed by 75 Se and then 192 Ir. Similarly, D2cc for all OAR improved with IMBT by up to 7.7% ± 0.6% (P textless 0.001). For IC/IS-BT cases, needle-free IMBT achieved clinically acceptable plans with 169 Yb-based IMBT further improving HR-CTV D98 by 1.5% ± 0.2% (P = 0.034) and decreasing sigmoid D2cc by 1.9% ± 0.4% (P = 0.048). Delivery times for IMBT are increased by a factor of 1.7, 3.3, and 2.3 for 192 Ir, 75 Se, and 169 Yb, respectively, relative to conventional 192 Ir BT. 

CONCLUSIONS: Dynamic shield IMBT provides a promising alternative to conventional IC- and IC/IS-BT techniques with significant dosimetric enhancements and even greater improvements with intermediate energy radionuclides. The ability to deliver a highly conformal, OAR-sparing dose without IS needles provides a simplified method for improving the therapeutic ratio less invasively and in a less resource intensive manner.},

keywords = {Brachytherapy, Computer-Assisted, Female, Humans, Image-guided cervix brachytherapy, Intensity modulated brachytherapy, Monte Carlo based dosimetry, Monte Carlo Method, MRI-guided GYN brachytherapy, Organs at Risk, Radiotherapy Dosage, Radiotherapy Planning, Uterine Cervical Neoplasms},

pubstate = {published},

tppubtype = {article}

}

@article{antaki_fast_2020,

title = {Fast mixed integer optimization (FMIO) for high dose rate brachytherapy},

author = {Majd Antaki and Christopher L Deufel and Shirin A. Enger},

doi = {10.1088/1361-6560/aba317},

issn = {1361-6560},

year  = {2020},

date = {2020-12-01},

journal = {Physics in Medicine and Biology},

volume = {65},

number = {21},

pages = {215005},

abstract = {The purpose of this work was to develop an efficient quadratic mixed integer programming algorithm for high dose rate (HDR) brachytherapy treatment planning problems and integrate the algorithm into an open-source Monte Carlo based treatment planning software, RapidBrachyMCTPS. The mixed-integer algorithm yields a globally optimum solution to the dose volume histogram (DVH) based problem and, unlike other methods, is not susceptible to local minimum trapping. A hybrid linear-quadratic penalty model coupled to a mixed integer programming model was used to optimize treatment plans for 10 prostate cancer patients. Dose distributions for each dwell position were calculated with RapidBrachyMCTPS with type A uncertainties less than 0.2% in voxels within the planning target volume (PTV). The optimization process was divided into two parts. First, the data was preprocessed, in which the problem size was reduced by eliminating voxels that had negligible impact on the solution (e.g. far from the dwell position). Second, the best combination of dwell times to obtain a plan with the highest score was found. The dwell positions and dose volume constraints were used as input to a commercial mixed integer optimizer (Gurobi Optimization, Inc.). A penalty-based criterion was adopted for the scoring. The voxel-reduction technique successfully reduced the problem size by an average of 91%, without loss of quality. The preprocessing of the optimization process required on average 4 s and solving for the global maximum required on average 33 s. The total optimization time averaged 37 s, which is a substantial improvement over the ∼15 min optimization time reported in published literature. The plan quality was evaluated by evaluating dose volume metrics, including PTV D90, rectum and bladder D1cc and urethra D0.1cc. In conclusion, fast mixed integer optimization is an order of magnitude faster than current mixed-integer approaches for solving HDR brachytherapy treatment planning problems with DVH based metrics.},

keywords = {Algorithms, Brachytherapy, Computer-Assisted, Humans, Linear Models, Male, Monte Carlo Method, Prostatic Neoplasms, Radiation Dosage, Radiotherapy Dosage, Radiotherapy Planning, Software, Time Factors},

pubstate = {published},

tppubtype = {article}

}

@article{mao_rapidbrachydl_2020,

title = {RapidBrachyDL: Rapid Radiation Dose Calculations in Brachytherapy Via Deep Learning},

author = {Ximeng Mao and Joelle Pineau and Roy Keyes and Shirin A. Enger},

doi = {10.1016/j.ijrobp.2020.04.045},

issn = {1879-355X},

year  = {2020},

date = {2020-11-01},

journal = {International Journal of Radiation Oncology, Biology, Physics},

volume = {108},

number = {3},

pages = {802--812},

abstract = {PURPOSE: Detailed and accurate absorbed dose calculations from radiation interactions with the human body can be obtained with the Monte Carlo (MC) method. However, the MC method can be slow for use in the time-sensitive clinical workflow. The aim of this study was to provide a solution to the accuracy-time trade-off for 192Ir-based high-dose-rate brachytherapy by using deep learning. 

METHODS AND MATERIALS: RapidBrachyDL, a 3-dimensional deep convolutional neural network (CNN) model, is proposed to predict dose distributions calculated with the MC method given a patient's computed tomography images, contours of clinical target volume (CTV) and organs at risk, and treatment plan. Sixty-one patients with prostate cancer and 10 patients with cervical cancer were included in this study, with data from 47 patients with prostate cancer being used to train the model. 

RESULTS: Compared with ground truth MC simulations, the predicted dose distributions by RapidBrachyDL showed a consistent shape in the dose-volume histograms (DVHs); comparable DVH dosimetric indices including 0.73% difference for prostate CTV D90, 1.1% for rectum D2cc, 1.45% for urethra D0.1cc, and 1.05% for bladder D2cc; and substantially smaller prediction time, acceleration by a factor of 300. RapidBrachyDL also demonstrated good generalization to cervical data with 1.73%, 2.46%, 1.68%, and 1.74% difference for CTV D90, rectum D2cc, sigmoid D2cc, and bladder D2cc, respectively, which was unseen during the training. 

CONCLUSION: Deep CNN-based dose estimation is a promising method for patient-specific brachytherapy dosimetry. Desired radiation quantities can be obtained with accuracies arbitrarily close to those of the source MC algorithm, but with much faster computation times. The idea behind deep CNN-based dose estimation can be safely extended to other radiation sources and tumor sites by following a similar training process.},

keywords = {Brachytherapy, Colon, Computer, Computer-Assisted, Deep Learning, Female, Humans, Iridium Radioisotopes, Male, Monte Carlo Method, Neural Networks, Organs at Risk, Prostate, Prostatic Neoplasms, Radiotherapy Dosage, Radiotherapy Planning, Rectum, Retrospective Studies, Sigmoid, Urinary Bladder, Uterine Cervical Neoplasms},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_monte_2020,

title = {Monte Carlo dosimetric characterization of a new high dose rate 169 Yb brachytherapy source and independent verification using a multipoint plastic scintillator detector},

author = {Gabriel Famulari and Haydee M. Linares Rosales and Justine Dupere and David C. Medich and Luc Beaulieu and Shirin A. Enger},

doi = {10.1002/mp.14336},

issn = {2473-4209},

year  = {2020},

date = {2020-09-01},

journal = {Medical Physics},

volume = {47},

number = {9},

pages = {4563--4573},

abstract = {PURPOSE: A prototype 169 Yb source was developed in combination with a dynamic rotating platinum shield system (AIM-Brachy) to deliver intensity modulated brachytherapy (IMBT). The purpose of this study was to evaluate the dosimetric characteristics of the bare/shielded 169 Yb source using Monte Carlo (MC) simulations and perform an independent dose verification using a dosimetry platform based on a multipoint plastic scintillator detector (mPSD). 

METHODS: The TG-43U1 dosimetric parameters were calculated for the source model using RapidBrachyMCTPS. Real-time dose rate measurements were performed in a water tank for both the bare/shielded source using a custom remote afterloader. For each dwell position, the dose rate was independently measured by the three scintillators (BCF-10, BCF-12, and BCF-60). For the bare source, dose rate was measured at distances up to 3 cm away from the source over a range of 7 cm along the catheter. For the shielded source, measurements were performed with the mPSD placed at 1 cm from the source at four different azimuthal angles ( 0 ∘ , 9 0 ∘ , 18 0 ∘ , and 27 0 ∘ ). 

RESULTS: The dosimetric parameters were tabulated for the source model. For the bare source, differences between measured and calculated along-away dose rates were generally below 5-10%. Along the transverse axis, deviations were, on average (range), 3.3% (0.6-6.2%) for BCF-10, 1.7% (0.9-2.9%) for BCF-12, and 2.2% (0.3-4.4%) for BCF-60. The maximum dose rate reduction due to shielding at a radial distance of 1 cm was 88.8 ± 1.2%, compared to 83.5 ± 0.5% as calculated by MC. 

CONCLUSIONS: The dose distribution for the bare/shielded 169 Yb source was independently verified using mPSD with good agreement in regions close to the source. The 169 Yb source coupled with the partial-shielding system is an effective technique to deliver IMBT.},

keywords = {Brachytherapy, Geant4, IMBT, Monte Carlo Method, mPSD, Plastics, Radiometry, Radiotherapy Dosage, shielded applicator, TG-186, TG-43, Yb-169},

pubstate = {published},

tppubtype = {article}

}

@article{famulari_novel_2020,

title = {A novel 169 Yb-based dynamic-shield intensity modulated brachytherapy delivery system for prostate cancer},

author = {Gabriel Famulari and Marie Duclos and Shirin A. Enger},

doi = {10.1002/mp.13959},

issn = {2473-4209},

year  = {2020},

date = {2020-03-01},

journal = {Medical Physics},

volume = {47},

number = {3},

pages = {859--868},

abstract = {PURPOSE: Intensity modulated brachytherapy (IMBT) is a novel high dose rate brachytherapy (HDR BT) technique which incorporates static or dynamic shielding to increase tumor coverage and/or spare healthy tissues. The purpose of this study is to present a novel delivery system (AIM-Brachy) design that can enable dynamic-shield IMBT for prostate cancer. 

METHODS: The AIM-Brachy system dynamically controls the rotation of platinum shields, placed within interstitial catheters, which partially collimate the radiation emitted from an 169 Yb source. Conventional HDR BT (10 Ci 192 Ir) and IMBT (18 Ci 169 Yb) plans were generated for 12 patients using an in-house column generation-based optimizer, coupled to a Geant4-based dose calculation engine, RapidBrachyMC. Treatment plans were normalized to match the same PTV D90 coverage as the clinical plan. Intershield attenuation effects were taken into account. A sensitivity analysis was performed to evaluate the dosimetric impact of systematic longitudinal source positioning errors ( ± 1 mm, ± 2 mm, and ± 3 mm) and rotational errors ( ± 5 ∘ , ± 10 ∘ , and ± 15 ∘ ) on clinically relevant parameters (PTV D90 and urethra D10 ). 

RESULTS: The platinum shield reduced the dose rate on the shielded side at 1 cm to 18.1% of the dose rate on the unshielded side. For equal PTV D90 coverage, the urethral D10 was reduced by 13.3% ± 4.7%, without change to other plan quality indices (PTV V100 , V150, V200 , bladder V75 , rectum V75 , HI, COIN). Delivery times for HDR BT and IMBT were 9.2 ± 1.6 min and 18.6 ± 4.0 min, respectively. In general, the PTV D90 was more sensitive to source positioning errors than rotational errors, while the urethral D10 was more sensitive to rotational errors than source positioning errors. For a typical range of positioning errors ( ± 1 mm, ± 5 ∘ ), the overall tolerance was textless2%. 

CONCLUSIONS: The AIM-Brachy system was proposed to deliver dynamic-shield IMBT for prostate cancer with the potential to create a low dose tunnel within the urethra. The urethra-sparing properties are desirable to minimize the occurrence and severity of urethral strictures or, alternatively, to provide a method for dose escalation.},

keywords = {Brachytherapy, Cohort Studies, Computer-Assisted, Humans, IMBT, Intensity-Modulated, Male, Monte Carlo, Monte Carlo Method, prostate cancer, Prostatic Neoplasms, Radioisotopes, Radiotherapy, Radiotherapy Planning, Uncertainty, Yb-169, Ytterbium},

pubstate = {published},

tppubtype = {article}

}

@article{morcos_monte_2020,

title = {Monte Carlo dosimetry study of novel rotating MRI-compatible shielded tandems for intensity modulated cervix brachytherapy},

author = {Marc Morcos and Shirin A. Enger},

doi = {10.1016/j.ejmp.2020.02.014},

issn = {1724-191X},

year  = {2020},

date = {2020-03-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 = {71},

pages = {178--184},

abstract = {PURPOSE: Intensity modulated brachytherapy (IMBT) with rotating metal shields enables dose modulation that can better conform to the tumor while reducing OAR doses. In this work, we investigate novel rotating shields, compatible with MRI-compatible tandems used for cervix brachytherapy. Three unique shields were evaluated using the traditional 192Ir source. Additionally, 75Se and 169Yb isotopes were investigated as alternative sources. 

METHODS: Three different IMBT shields were modeled and simulated in RapidBrachyMCTPS. Each tungsten shield was designed to fit inside a 6 mm-wide MRI-compatible tandem. The active core of the source was replaced with 192Ir, 75Se and 169Yb. Transmission factors (TFs) were calculated and defined as the dose ratio at 1 cm on opposite sides of the shielded tandem on the transverse plane. Polar and azimuthal anisotropy plots were extracted from simulations. Dose homogeneities V200%V100% were calculated for all radionuclide-shield combinations. 

RESULTS: TFs are favorable for IMBT and ranged between 12.9% and 32.2% for 192Ir, 4.0%-16.1% for 75Se and 1.2-6.4% for 169Yb for all shield designs. Average beam-widths in the polar and azimuthal directions were reduced to the range of 42°-112° and 27°-107°, respectively, for all shield-radionuclide combinations. Dose homogeneities for all the radionuclide-shield combinations were within 12% of the non-IMBT tandem. 

CONCLUSIONS: This study has quantitatively assessed the influence of various rotating cervical cancer-specific IMBT tandem shields on dosimetry. The dynamic single-channel shields and narrow beam-widths in the polar and azimuthal direction give rise to highly anisotropic distributions. Intermediate-to-high energy radionuclides, 75Se and 169Yb substantially improve the modulation capacity of IMBT and pave the way for treating large and complex cervical cancer without interstitial needle implantation.},

keywords = {Anisotropy, Brachytherapy, Female, Humans, Image-guided cervix brachytherapy, Intensity modulated brachytherapy, Intensity-Modulated, Iridium Radioisotopes, Magnetic Resonance Imaging, Monte Carlo based dosimetry, Monte Carlo Method, MRI-guided GYN brachytherapy, Radiometry, Radiotherapy, Selenium Radioisotopes, Uterine Cervical Neoplasms, Ytterbium},

pubstate = {published},

tppubtype = {article}

}

@article{shoemaker_dosimetric_2019,

title = {Dosimetric Considerations for Ytterbium-169, Selenium-75, and Iridium-192 Radioisotopes in High-Dose-Rate Endorectal Brachytherapy},

author = {Tristan Shoemaker and Té Vuong and Harry Glickman and Samar Kaifi and Gabriel Famulari and Shirin A. Enger},

doi = {10.1016/j.ijrobp.2019.07.003},

issn = {1879-355X},

year  = {2019},

date = {2019-11-01},

journal = {International Journal of Radiation Oncology, Biology, Physics},

volume = {105},

number = {4},

pages = {875--883},

abstract = {PURPOSE: To investigate differences between prescribed and postimplant calculated dose in 192Ir high-dose-rate endorectal brachytherapy (HDR-EBT) by evaluating dose to clinical target volume (CTV) and organs at risk (OARs) calculated with a Monte Carlo-based dose calculation software, RapidBrachyMC. In addition, dose coverage, conformity, and homogeneity were compared among the radionuclides 192Ir, 75Se, and 169Yb for use in HDR-EBT. 

METHODS AND MATERIALS: Postimplant dosimetry was evaluated using 23 computed tomography (CT) images from patients treated with HDR-EBT using the 192Ir microSelectron v2 (Elekta AB, Stockholm, Sweden) source and the Intracavitary Mold Applicator Set (Elekta AB, Stockholm, Sweden), which is a flexible applicator capable of fitting a tungsten rod for OAR shielding. Four tissue segmentation schemes were evaluated: (1) TG-43 formalism, (2) materials and nominal densities assigned to contours of foreign objects, (3) materials and nominal densities assigned to contoured organs in addition to foreign objects, and (4) materials specified as in (3) but with voxel mass densities derived from CT Hounsfield units. Clinical plans optimized for 192Ir were used, with the results for 75Se and 169Yb normalized to the D90 of the 192Ir clinical plan. RESULTS: In comparison to segmentation scheme 4, TG-43-based dosimetry overestimates CTV D90 by 6% (P = .00003), rectum D50 by 24% (P = .00003), and pelvic bone D50 by 5% (P = .00003) for 192Ir. For 169Yb, CTV D90 is overestimated by 17% (P = .00003) and rectum D50 by 39% (P = .00003), and pelvic bone D50 is significantly underestimated by 27% (P = .007). Postimplant dosimetry calculations also showed that a 169Yb source would give 20% (P = .00003) lower rectum V60 and 17% (P = .00008) lower rectum D50. 

CONCLUSIONS: Ignoring high-Z materials in dose calculation contributes to inaccuracies that may lead to suboptimal dose optimization and disagreement between prescribed and calculated dose. This is especially important for low-energy radionuclides. Our results also show that with future magnetic resonance imaging-based treatment planning, loss of CT density data will only affect calculated dose in nonbone OARs by 2% or less and bone OARs by 13% or less across all sources if material composition and nominal mass densities are correctly assigned.},

keywords = {Brachytherapy, Femur, Humans, Iridium Radioisotopes, Monte Carlo Method, Organs at Risk, Pelvic Bones, Radioisotopes, Radiotherapy Dosage, Rectal Neoplasms, Rectum, Selenium Radioisotopes, Tomography, Urinary Bladder, X-Ray Computed, Ytterbium},

pubstate = {published},

tppubtype = {article}

}

@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}

}

@article{decunha_new_2018,

title = {A new delivery system to resolve dosimetric issues in intravascular brachytherapy},

author = {Joseph M. DeCunha and Shirin A. Enger},

doi = {10.1016/j.brachy.2018.01.012},

issn = {1873-1449},

year  = {2018},

date = {2018-06-01},

journal = {Brachytherapy},

volume = {17},

number = {3},

pages = {634--643},

abstract = {PURPOSE: Renewed interest is being expressed in intravascular brachytherapy (IVBT). A number of unresolved issues exist in the discipline. Providing a homogeneous and adequate dose to the target remains difficult in IVBT. The guidewire that delivers the device to the target, arterial plaques, and stent struts are all known to reduce the dose delivered to target. The viability and efficacy of a proposed IVBT delivery system designed to resolve the issue of guidewire attenuation is evaluated and compared to that of a popular and commercially available IVBT device. 

METHODS AND MATERIALS: Monte Carlo simulations are conducted to determine distributions of absorbed dose around an existing and proposed IVBT delivery system. 

RESULTS: For the Novoste Beta-Cath 3.5F (TeamBest®), dose in water varies by 10% as a function of angle in the plane perpendicular to the delivery catheter due to off-centering of seeds in the catheter. Dose is reduced by 52% behind a stainless steel guidewire and 64% behind a guidewire, arterial plaque, and stent strut for the Novoste Beta-Cath 3.5F. Dose is not perturbed by the presence of a guidewire for the proposed device and is reduced by 46% by an arterial plaque and stent strut. 

CONCLUSIONS: Dose attenuation by guidewire is likely the single greatest source of dose attenuation in IVBT in terms of absolute dose reduction and is greater than previously reported for the Novoste Beta-Cath 3.5F. The Novoste Beta-Cath 3.5F delivers an inhomogeneous dose to target. A delivery system is proposed, which resolves the issue of guidewire attenuation in IVBT and should reduce treatment times.},

keywords = {Brachytherapy, Catheterization, Catheters, Computer Simulation, Coronary Vessels, Humans, Intravascular, Monte Carlo Method, Physics, Radiation Dosage, Radiometry, Restenosis, Stents, Strontium Radioisotopes},

pubstate = {published},

tppubtype = {article}

}

@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}

}

@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}

}

@article{poole_determination_2015,

title = {Determination of subcellular compartment sizes for estimating dose variations in radiotherapy},

author = {Christopher M. Poole and Anders Ahnesjö and Shirin A. Enger},

doi = {10.1093/rpd/ncv305},

issn = {1742-3406},

year  = {2015},

date = {2015-09-01},

journal = {Radiation Protection Dosimetry},

volume = {166},

number = {1-4},

pages = {361--364},

abstract = {The variation in specific energy absorbed to different cell compartments caused by variations in size and chemical composition is poorly investigated in radiotherapy. The aim of this study was to develop an algorithm to derive cell and cell nuclei size distributions from 2D histology samples, and build 3D cellular geometries to provide Monte Carlo (MC)-based dose calculation engines with a morphologically relevant input geometry. Stained and unstained regions of the histology samples are segmented using a Gaussian mixture model, and individual cell nuclei are identified via thresholding. Delaunay triangulation is applied to determine the distribution of distances between the centroids of nearest neighbour cells. A pouring simulation is used to build a 3D virtual tissue sample, with cell radii randomised according to the cell size distribution determined from the histology samples. A slice with the same thickness as the histology sample is cut through the 3D data and characterised in the same way as the measured histology. The comparison between this virtual slice and the measured histology is used to adjust the initial cell size distribution into the pouring simulation. This iterative approach of a pouring simulation with adjustments guided by comparison is continued until an input cell size distribution is found that yields a distribution in the sliced geometry that agrees with the measured histology samples. The thus obtained morphologically realistic 3D cellular geometry can be used as input to MC-based dose calculation programs for studies of dose response due to variations in morphology and size of tumour/healthy tissue cells/nuclei, and extracellular material.},

keywords = {Algorithms, Breast Neoplasms, Cell Nucleus, Computer Simulation, Computer-Assisted, ErbB-2, Female, Humans, Image Processing, Imaging, Immunoenzyme Techniques, Male, Monte Carlo Method, Prostatic Neoplasms, Radiotherapy Dosage, Radiotherapy Planning, Receptor, Signal Processing, Subcellular Fractions, Three-Dimensional},

pubstate = {published},

tppubtype = {article}

}

@article{adams_interstitial_2014,

title = {Interstitial rotating shield brachytherapy for prostate cancer},

author = {Quentin E. Adams and Jinghzu Xu and Elizabeth K. Breitbach and Xing Li and Shirin A. Enger and William R. Rockey and Yusung Kim and Xiaodong Wu and Ryan T. Flynn},

doi = {10.1118/1.4870441},

issn = {2473-4209},

year  = {2014},

date = {2014-05-01},

journal = {Medical Physics},

volume = {41},

number = {5},

pages = {051703},

abstract = {PURPOSE: To present a novel needle, catheter, and radiation source system for interstitial rotating shield brachytherapy (I-RSBT) of the prostate. I-RSBT is a promising technique for reducing urethra, rectum, and bladder dose relative to conventional interstitial high-dose-rate brachytherapy (HDR-BT). 

METHODS: A wire-mounted 62 GBq(153)Gd source is proposed with an encapsulated diameter of 0.59 mm, active diameter of 0.44 mm, and active length of 10 mm. A concept model I-RSBT needle/catheter pair was constructed using concentric 50 and 75 μm thick nickel-titanium alloy (nitinol) tubes. The needle is 16-gauge (1.651 mm) in outer diameter and the catheter contains a 535 μm thick platinum shield. I-RSBT and conventional HDR-BT treatment plans for a prostate cancer patient were generated based on Monte Carlo dose calculations. In order to minimize urethral dose, urethral dose gradient volumes within 0-5 mm of the urethra surface were allowed to receive doses less than the prescribed dose of 100%. 

RESULTS: The platinum shield reduced the dose rate on the shielded side of the source at 1 cm off-axis to 6.4% of the dose rate on the unshielded side. For the case considered, for the same minimum dose to the hottest 98% of the clinical target volume (D(98%)), I-RSBT reduced urethral D(0.1cc) below that of conventional HDR-BT by 29%, 33%, 38%, and 44% for urethral dose gradient volumes within 0, 1, 3, and 5 mm of the urethra surface, respectively. Percentages are expressed relative to the prescription dose of 100%. For the case considered, for the same urethral dose gradient volumes, rectum D(1cc) was reduced by 7%, 6%, 6%, and 6%, respectively, and bladder D(1cc) was reduced by 4%, 5%, 5%, and 6%, respectively. Treatment time to deliver 20 Gy with I-RSBT was 154 min with ten 62 GBq (153)Gd sources. 

CONCLUSIONS: For the case considered, the proposed(153)Gd-based I-RSBT system has the potential to lower the urethral dose relative to HDR-BT by 29%-44% if the clinician allows a urethral dose gradient volume of 0-5 mm around the urethra to receive a dose below the prescription. A multisource approach is necessary in order to deliver the proposed (153)Gd-based I-RSBT technique in reasonable treatment times.},

keywords = {Brachytherapy, Catheters, Computer-Assisted, Equipment Design, Gadolinium, Humans, Iridium Radioisotopes, Male, Monte Carlo Method, Needles, Nickel, Platinum Compounds, Prostatic Neoplasms, Radiation Protection, Radioisotopes, Radiotherapy Dosage, Radiotherapy Planning, Rectum, Time Factors, Titanium, Urethra, Urinary Bladder},

pubstate = {published},

tppubtype = {article}

}

@article{enger_gadolinium-153_2013,

title = {Gadolinium-153 as a brachytherapy isotope},

author = {Shirin A. Enger and Darrell R. Fisher and Ryan T. Flynn},

doi = {10.1088/0031-9155/58/4/957},

issn = {1361-6560},

year  = {2013},

date = {2013-02-01},

journal = {Physics in Medicine and Biology},

volume = {58},

number = {4},

pages = {957--964},

abstract = {The purpose of this work was to present the fundamental dosimetric characteristics of a hypothetical (153)Gd brachytherapy source using the AAPM TG-43U1 dose-calculation formalism. Gadolinium-153 is an intermediate-energy isotope that emits 40-100 keV photons with a half-life of 242 days. The rationale for considering (153)Gd as a brachytherapy source is for its potential of patient specific shielding and to enable reduced personnel shielding requirements relative to (192)Ir, and as an isotope for interstitial rotating shield brachytherapy (I-RSBT). A hypothetical (153)Gd brachytherapy source with an active core of 0.84 mm diameter, 10 mm length and specific activity of 5.55 TBq of (153)Gd per gram of Gd was simulated with Geant4. The encapsulation material was stainless steel with a thickness of 0.08 mm. The radial dose function, anisotropy function and photon spectrum in water were calculated for the (153)Gd source. The simulated (153)Gd source had an activity of 242 GBq and a dose rate in water 1 cm off axis of 13.12 Gy h(-1), indicating that it would be suitable as a low-dose-rate or pulsed-dose-rate brachytherapy source. The beta particles emitted have low enough energies to be absorbed in the source encapsulation. Gadolinium-153 has an increasing radial dose function due to multiple scatter of low-energy photons. Scattered photon dose takes over with distance from the source and contributes to the majority of the absorbed dose. The anisotropy function of the (153)Gd source decreases at low polar angles, as a result of the long active core. The source is less anisotropic at polar angles away from the longitudinal axes. The anisotropy function increases with increasing distance. The (153)Gd source considered would be suitable as an intermediate-energy low-dose-rate or pulsed-dose-rate brachytherapy source. The source could provide a means for I-RSBT delivery and enable brachytherapy treatments with patient specific shielding and reduced personnel shielding requirements relative to (192)Ir.},

keywords = {Anisotropy, Brachytherapy, Equipment Design, Gadolinium, Humans, Iridium Radioisotopes, Male, Monte Carlo Method, Photons, Prostatic Neoplasms, Radiation, Radiation Protection, Radioisotopes, Radiotherapy Dosage, Scattering},

pubstate = {published},

tppubtype = {article}

}

@article{enger_layered_2012,

title = {Layered mass geometry: a novel technique to overlay seeds and applicators onto patient geometry in Geant4 brachytherapy simulations},

author = {Shirin A. Enger and Guillaume Landry and Michel D'Amours and Frank Verhaegen and Luc Beaulieu and Makoto Asai and Joseph Perl},

doi = {10.1088/0031-9155/57/19/6269},

issn = {1361-6560},

year  = {2012},

date = {2012-10-01},

journal = {Physics in Medicine and Biology},

volume = {57},

number = {19},

pages = {6269--6277},

abstract = {A problem faced by all Monte Carlo (MC) particle transport codes is how to handle overlapping geometries. The Geant4 MC toolkit allows the user to create parallel geometries within a single application. In Geant4 the standard mass-containing geometry is defined in a simulation volume called the World Volume. Separate parallel geometries can be defined in parallel worlds, that is, alternate three dimensional simulation volumes that share the same coordinate system with the World Volume for geometrical event biasing, scoring of radiation interactions, and/or the creation of hits in detailed readout structures. Until recently, only one of those worlds could contain mass so these parallel worlds provided no solution to simplify a complex geometric overlay issue in brachytherapy, namely the overlap of radiation sources and applicators with a CT based patient geometry. The standard method to handle seed and applicator overlay in MC requires removing CT voxels whose boundaries would intersect sources, placing the sources into the resulting void and then backfilling the remaining space of the void with a relevant material. The backfilling process may degrade the accuracy of patient representation, and the geometrical complexity of the technique precludes using fast and memory-efficient coding techniques that have been developed for regular voxel geometries. The patient must be represented by the less memory and CPU-efficient Geant4 voxel placement technique, G4PVPlacement, rather than the more efficient G4NestedParameterization (G4NestedParam). We introduce for the first time a Geant4 feature developed to solve this issue: Layered Mass Geometry (LMG) whereby both the standard (CT based patient geometry) and the parallel world (seeds and applicators) may now have mass. For any area where mass is present in the parallel world, the parallel mass is used. Elsewhere, the mass of the standard world is used. With LMG the user no longer needs to remove patient CT voxels that would include for example seeds. The patient representation can be a regular voxel grid, conducive to G4NestedParam, and the patient CT derived materials remain exact, avoiding the inaccuracy of the backfilling technique. Post-implant dosimetry for one patient with (125)I permanent seed implant was performed using Geant4 version 9.5.p01 using three different geometrical techniques. The first technique was the standard described above (G4PVPlacement). The second technique placed patient voxels as before, but placed seeds with LMG (G4PVPlacement+LMG). The third technique placed patient voxels through G4NestedParam and seeds through LMG (G4NestedParam+LMG). All the scenarios were calculated with 3 different image compression factors to manipulate the number of voxels. Additionally, the dosimetric impact of the backfilling technique was investigated for the case of calcifications in close proximity of sources. LMG eliminated the need for backfilling and simplified geometry description. Of the two LMG techniques, G4PVPlacement+LMG had no benefit to calculation time or memory use, actually increasing calculation time, but G4NestedParam+LMG reduced both calculation time and memory. The benefits of G4NestedParam+LMG over standard G4PVPlacement increased with increasing voxel numbers. For the case of calcifications in close proximity to sources, LMG not only increased efficiency but also yielded more accurate dose calculation than G4PVPlacement. G4NestedParam in combination with LMG present a new, efficient approach to simulate radiation sources that overlap patient geometry. Cases with brachytherapy applicators would constitute a direct extension of the method.},

keywords = {Brachytherapy, Computer-Assisted, Humans, Monte Carlo Method, Radiotherapy Dosage, Radiotherapy Planning},

pubstate = {published},

tppubtype = {article}

}

@article{enger_dose_2012,

title = {Dose to tissue medium or water cavities as surrogate for the dose to cell nuclei at brachytherapy photon energies},

author = {Shirin A. Enger and Anders Ahnesjö and Frank Verhaegen and Luc Beaulieu},

doi = {10.1088/0031-9155/57/14/4489},

issn = {1361-6560},

year  = {2012},

date = {2012-07-01},

journal = {Physics in Medicine and Biology},

volume = {57},

number = {14},

pages = {4489--4500},

abstract = {It has been suggested that modern dose calculation algorithms should be able to report absorbed dose both as dose to the local medium, D(m,m,) and as dose to a water cavity embedded in the medium, D(w,m), using conversion factors from cavity theory. Assuming that the cell nucleus with its DNA content is the most important target for biological response, the aim of this study is to investigate, by means of Monte Carlo (MC) simulations, the relationship of the dose to a cell nucleus in a medium, D(n,m,) to D(m,m) and D(w,m), for different combinations of cell nucleus compositions and tissue media for different photon energies used in brachytherapy. As D(n,m) is very impractical to calculate directly for routine treatment planning, while D(m,m) and D(w,m) are much easier to obtain, the questions arise which one of these quantities is the best surrogate for D(n,m) and which cavity theory assumptions should one use for its estimate. The Geant4.9.4 MC code was used to calculate D(m,m,) D(w,m) and D(n,m) for photon energies from 20 (representing the lower energy end of brachytherapy for ¹⁰³Pd or ¹²⁵I) to 300 keV (close to the mean energy of (¹⁹²Ir) and for the tissue media adipose, breast, prostate and muscle. To simulate the cell and its nucleus, concentric spherical cavities were placed inside a cubic phantom (10 × 10 × 10 mm³). The diameter of the simulated nuclei was set to 14 µm. For each tissue medium, three different setups were simulated; (a) D(n,m) was calculated with nuclei embedded in tissues (MC-D(n,m)). Four different published elemental compositions of cell nuclei were used. (b) D(w,m) was calculated with MC (MC-D(w,m)) and compared with large cavity theory calculated D(w,m) (LCT-D(w,m)), and small cavity theory calculated D(w,m) (SCT-D(w,m)). (c) D(m,m) was calculated with MC (MC-D(m,m)). MC-D(w,m) is a good substitute for MC-D(n,m) for all photon energies and for all simulated nucleus compositions and tissue types. SCT-D(w,m) can be used for most energies in brachytherapy, while LCT-D(w,m) should only be considered for source spectra well below 50 keV, since contributions to the absorbed dose inside the nucleus to a large degree stem from electrons released in the surrounding medium. MC-D(m,m) is not an appropriate substitute for MC-D(n,m) for the lowest photon energies for adipose and breast tissues. The ratio of MC-D(m,m) to MC-D(n,m) for adipose and breast tissue deviates from unity by 34% and 15% respectively for the lowest photon energy (20 keV), whereas the ratio is close to unity for higher energies. For prostate and muscle tissue MC-D(m,m) is a good substitute for MC-D(n,m). However, for all photon energies and tissue types the nucleus composition with the highest hydrogen content behaves differently than other compositions. Elemental compositions of the tissue and nuclei affect considerably the absorbed dose to the cell nuclei for brachytherapy sources, in particular those at the low-energy end of the spectrum. Thus, there is a need for more accurate data for the elemental compositions of tumours and healthy cells. For the nucleus compositions and tissue types investigated, MC-D(w,m) is a good substitute to MC-D(n,m) for all simulated photon energies. Whether other studied surrogates are good approximations to MC-D(n,m) depends on the target size, target composition, composition of the surrounding tissue and photon energy.},

keywords = {Brachytherapy, Cell Line, Cell Nucleus, Humans, Monte Carlo Method, Photons, Radiation Dosage, Radiotherapy Dosage, Water},

pubstate = {published},

tppubtype = {article}

}

@article{enger_exploring_2012,

title = {Exploring (57)Co as a new isotope for brachytherapy applications},

author = {Shirin A. Enger and Hans Lundqvist and Michel D'Amours and Luc Beaulieu},

doi = {10.1118/1.3700171},

issn = {0094-2405},

year  = {2012},

date = {2012-05-01},

journal = {Medical Physics},

volume = {39},

number = {5},

pages = {2342--2345},

abstract = {PURPOSE: The characteristics of the radionuclide (57)Co make it interesting for use as a brachytherapy source. (57)Co combines a possible high specific activity with the emission of relatively low-energy photons and a half-life (272 days) suitable for regular source exchanges in an afterloader. (57)Co decays by electron capture to the stable (57)Fe with emission of 136 and 122 keV photons. 

METHODS: A hypothetical (57)Co source based on the Flexisource brachytherapy encapsulation with the active core set as a pure cobalt cylinder (length 3.5 mm and diameter 0.6 mm) covered with a cylindrical stainless-steel capsule (length 5 mm and thickness 0.125 mm) was simulated using Geant4 Monte Carlo (MC) code version 9.4. The radial dose function, g(r), and anisotropy function F(r,θ), for the line source approximation were calculated following the TG-43U1 formalism. The results were compared to well-known (192)Ir and (125)I radionuclides, representing the higher and the lower energy end of brachytherapy, respectively. 

RESULTS: The mean energy of photons in water, after passing through the core and the encapsulation material was 123 keV. This hypothetical (57)Co source has an increasing g(r) due to multiple scatter of low-energy photons, which results in a more uniform dose distribution than (192)Ir. 

CONCLUSIONS: (57)Co has many advantages compared to (192)Ir due to its low-energy gamma emissions without any electron contamination. (57)Co has an increasing g(r) that results in a more uniform dose distribution than (192)Ir due to its multiple scattered photons. The anisotropy of the (57)Co source is comparable to that of (192)Ir. Furthermore, (57)Co has lower shielding requirements than (192)Ir.},

keywords = {Anisotropy, Brachytherapy, Cobalt Radioisotopes, Monte Carlo Method, Radiation, Radiometry, Scattering},

pubstate = {published},

tppubtype = {article}

}

@article{enger_modeling_2011,

title = {Modeling a hypothetical 170Tm source for brachytherapy applications},

author = {Shirin A. Enger and Michel D'Amours and Luc Beaulieu},

doi = {10.1118/1.3626482},

issn = {0094-2405},

year  = {2011},

date = {2011-10-01},

journal = {Medical Physics},

volume = {38},

number = {10},

pages = {5307--5310},

abstract = {PURPOSE: To perform absorbed dose calculations based on Monte Carlo simulations for a hypothetical (170)Tm source and to investigate the influence of encapsulating material on the energy spectrum of the emitted electrons and photons. 

METHODS: GEANT4 Monte Carlo code version 9.2 patch 2 was used to simulate the decay process of (170)Tm and to calculate the absorbed dose distribution using the GEANT4 Penelope physics models. A hypothetical (170)Tm source based on the Flexisource brachytherapy design with the active core set as a pure thulium cylinder (length 3.5 mm and diameter 0.6 mm) and different cylindrical source encapsulations (length 5 mm and thickness 0.125 mm) constructed of titanium, stainless-steel, gold, or platinum were simulated. The radial dose function for the line source approximation was calculated following the TG-43U1 formalism for the stainless-steel encapsulation. 

RESULTS: For the titanium and stainless-steel encapsulation, 94% of the total bremsstrahlung is produced inside the core, 4.8 and 5.5% in titanium and stainless-steel capsules, respectively, and less than 1% in water. For the gold capsule, 85% is produced inside the core, 14.2% inside the gold capsule, and a negligible amount (textless1%) in water. Platinum encapsulation resulted in bremsstrahlung effects similar to those with the gold encapsulation. The range of the beta particles decreases by 1.1 mm with the stainless-steel encapsulation compared to the bare source but the tissue will still receive dose from the beta particles several millimeters from the source capsule. The gold and platinum capsules not only absorb most of the electrons but also attenuate low energy photons. The mean energy of the photons escaping the core and the stainless-steel capsule is 113 keV while for the gold and platinum the mean energy is 160 keV and 165 keV, respectively. 

CONCLUSIONS: A (170)Tm source is primarily a bremsstrahlung source, with the majority of bremsstrahlung photons being generated in the source core and experiencing little attenuation in the source encapsulation. Electrons are efficiently absorbed by the gold and platinum encapsulations. However, for the stainless-steel capsule (or other lower Z encapsulations) electrons will escape. The dose from these electrons is dominant over the photon dose in the first few millimeter but is not taken into account by current standard treatment planning systems. The total energy spectrum of photons emerging from the source depends on the encapsulation composition and results in mean photon energies well above 100 keV. This is higher than the main gamma-ray energy peak at 84 keV. Based on our results, the use of (170)Tm as a brachytherapy source presents notable challenges.},

keywords = {Algorithms, Brachytherapy, Computer Simulation, Computer-Assisted, Electrons, Equipment Design, Gold, Humans, Models, Monte Carlo Method, Photons, Platinum, Radioisotopes, Radiotherapy Planning, Stainless Steel, Theoretical, Thulium, Titanium},

pubstate = {published},

tppubtype = {article}

}

@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}

}

@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}

}

@article{enger_gadolinium_2006,

title = {Gadolinium neutron capture brachytherapy (GdNCB), a new treatment method for intravascular brachytherapy},

author = {Shirin A. Enger and Arash Rezaei and Per Munck af Rosenschöld and Hans Lundqvist},

doi = {10.1118/1.2146050},

issn = {0094-2405},

year  = {2006},

date = {2006-01-01},

journal = {Medical Physics},

volume = {33},

number = {1},

pages = {46--51},

abstract = {Restenosis is a major problem after balloon angioplasty and stent implantation. The aim of this study is to introduce gadolinium neutron capture brachytherapy (GdNCB) as a suitable modality for treatment of stenosis. The utility of GdNCB in intravascular brachytherapy (IVBT) of stent stenosis is investigated by using the GEANT4 and MCNP4B Monte Carlo radiation transport codes. To study capture rate, Kerma, absorbed dose and absorbed dose rate around a Gd-containing stent activated with neutrons, a 30 mm long, 5 mm diameter gadolinium foil is chosen. The input data is a neutron spectrum used for clinical neutron capture therapy in Studsvik, Sweden. Thermal neutron capture in gadolinium yields a spectrum of high-energy gamma photons, which due to the build-up effect gives an almost flat dose delivery pattern to the first 4 mm around the stent. The absorbed dose rate is 1.33 Gy/min, 0.25 mm from the stent surface while the dose to normal tissue is in order of 0.22 Gy/min, i.e., a factor of 6 lower. To spare normal tissue further fractionation of the dose is also possible. The capture rate is relatively high at both ends of the foil. The dose distribution from gamma and charge particle radiation at the edges and inside the stent contributes to a nonuniform dose distribution. This will lead to higher doses to the surrounding tissue and may prevent stent edge and in-stent restenosis. The position of the stent can be verified and corrected by the treatment plan prior to activation. Activation of the stent by an external neutron field can be performed days after catherization when the target cells start to proliferate and can be expected to be more radiation sensitive. Another advantage of the nonradioactive gadolinium stent is the possibility to avoid radiation hazard to personnel.},

keywords = {Biological, Blood Vessel Prosthesis, Brachytherapy, Computer Simulation, Computer-Assisted, Gadolinium, Graft Occlusion, Humans, Models, Monte Carlo Method, Neutron Capture Therapy, Radiometry, Radiotherapy Dosage, Radiotherapy Planning, Relative Biological Effectiveness, Statistical, Stents, Vascular},

pubstate = {published},

tppubtype = {article}

}