Biomedical subjects

Sung-Joon Ye

Publications and source records attributed to Sung-Joon Ye.

10 recordsLinked to original sources

PubMed2006

Doppler broadening effect on low-energy photon dose calculations using MCNP5 and PENELOPE.

Recent releases of the MCNP5 and PENELOPE Monte Carlo codes include the transport algorithm and momentum profiles that are necessary for accounting for Doppler broadening in Compton scattering processes. Such improvements might be particularly important in low-energy photon dose calculations. MCPLIB04 and PENDBASE (PENELOPE photon dataset) are based on the EPDL97 library with Compton momentum profiles, while MCPLIB03 and MCPLIB02 are based on the 1970's old library, with MCPLIB03 including the Compton momentum profiles. To isolate the dosimetric effects of Doppler broadening by the transport algorithm and Compton momentum profiles, we varied the choice of the above photon databases, in the same simulation geometry, using either version of MCNP5 or MCNP4 (no Doppler algorithm). We computed dose rate constants and dose distributions for r = 0.2-10 cm from a point source in a 50-cm-diameter sphere of water. Nine discrete energies for primary photon sources were chosen in the range of 10-150 keV. The results from both versions of MCNP with MCPLIB04 agreed with those of PENELOPE within statistical uncertainties (+/-1%) over the entire ranges of energies and radial distances investigated. MCNP5 with either MCPLIB03 or MCPLIB02 yielded almost identical data within statistical uncertainties (+/-1%) over the entire ranges of energies and radial distances investigated. This implies that in spite of the spectral broadening of scattered photons due to the orbital electron motion, the dosimetric effect of Doppler broadening for Compton interactions in water appears to be insignificant in the energy range investigated. The spectral dose analysis with and without the Doppler broadening supported this conclusion.

Algorithms↗

PubMed2006

Comprehensive evaluation of a commercial macro Monte Carlo electron dose calculation implementation using a standard verification data set.

A commercial electron dose calculation software implementation based on the macro Monte Carlo algorithm has recently been introduced. We have evaluated the performance of the system using a standard verification data set comprised of two-dimensional (2D) dose distributions in the transverse plane of a 15 X 15 cm2 field. The standard data set was comprised of measurements performed for combinations of 9-MeV and 20-MeV beam energies and five phantom geometries. The phantom geometries included bone and air heterogeneities, and irregular surface contours. The standard verification data included a subset of the data needed to commission the dose calculation. Additional required data were obtained from a dosimetrically equivalent machine. In addition, we performed 2D dose measurements in a water phantom for the standard field sizes, a 4 cm X 4 cm field, a 3 cm diameter circle, and a 5 cm X 13 cm triangle for the 6-, 9-, 12-, 15-, and 18-MeV energies of a Clinac 21EX. Output factors were also measured. Synthetic CT images and structure contours duplicating the measurement configurations were generated and transferred to the treatment planning system. Calculations for the standard verification data set were performed over the range of each of the algorithm parameters: statistical precision, grid-spacing, and smoothing. Dose difference and distance-to-agreement were computed for the calculation points. We found that the best results were obtained for the highest statistical precision, for the smallest grid spacing, and for smoothed dose distributions. Calculations for the 21EX data were performed using parameters that the evaluation of the standard verification data suggested would produce clinically acceptable results. The dose difference and distance-to-agreement were similar to that observed for the standard verification data set except for the portion of the triangle field narrower than 3 cm for the 6- and 9-MeV electron beams. The output agreed with measurements to within 2%, with the exception of the 3-cm diameter circle and the triangle for 6 MeV, which were within 5%. We conclude that clinically acceptable results may be obtained using a grid spacing that is no larger than approximately one-tenth of the distal falloff distance of the electron depth dose curve (depth from 80% to 20% of the maximum dose) and small relative to the size of heterogeneities. For judicious choices of parameters, dose calculations agree with measurements to better than 3% dose difference and 3-mm distance-to-agreement for fields with dimensions no less than about 3 cm.

Journal Article↗

PubMed2005

Using a photon phase-space source for convolution/superposition dose calculations in radiation therapy.

For a given linac design, the dosimetric characteristics of a photon beam are determined uniquely by the energy and radial distributions of the electron beam striking the x-ray target. However, in the usual commissioning of a beam from measured data, a large number of variables can be independently tuned, making it difficult to derive a unique and self-consistent beam model. For example, the measured dosimetric penumbra in water may be attributed in various proportions to the lateral secondary electron range, the focal spot size and the transmission through the tips of a non-divergent collimator; the head-scatter component in the tails of the transverse profiles may not be easy to resolve from phantom scatter and head leakage; and the head-scatter tails corresponding to a certain extra-focal source model may not agree self-consistently with in-air output factors measured on the central axis. To reduce the number of adjustable variables in beam modelling, we replace the focal and extra-focal sources with a single phase-space plane scored just above the highest adjustable collimator in a EGS/BEAM simulation of the linac. The phase-space plane is then used as photon source in a stochastic convolution/superposition dose engine. A photon sampled from the uncollimated phase-space plane is first propagated through an arbitrary collimator arrangement and then interacted in the simulation phantom. Energy deposition kernel rays are then randomly issued from the interaction points and dose is deposited along these rays. The electrons in the phase-space file are used to account for electron contamination. 6 MV and 18 MV photon beams from an Elekta SL linac are used as representative examples. Except for small corrections for monitor backscatter and collimator forward scatter for large field sizes (<0.5% with <20 x 20 cm2 field size), we found that the use of a single phase-space photon source provides accurate and self-consistent results for both relative and absolute dose calculations.

Algorithms↗

PubMed2005

Determination of field size-dependent wedge factors from a few selected measurements.

Some modern treatment-planning systems (TPSs) provide for input of wedge factor (WF) tables covering the entire range of square and elongated fields available on the LINAC. Depending on the field size increment chosen and the number of available wedge orientations, one may have to take more than 100 measurements per wedge and photon energy to commission the TPS. To expedite TPS commissioning while maintaining high accuracy, we demonstrate a simple method that requires only a few measurements per wedge, from which the remaining wedge factors can be found through linear interpolation based on field area. For the externally mounted wedges of two common LINACs, we have shown that WFs are proportional to field area and are nearly independent of field elongation and wedge orientation. Wedge factors computed from five to seven measurements comprised of square fields and a single, large rectangular field agreed with direct measurements throughout the entire range of achievable field dimensions within 0.6% at 6 MV and within 1% at 15 MV. Making the same set of measurements and using the equivalent square method to find WFs at other field sizes leads to errors up to 2%. Measuring the WF for a 10 x 10 cm2 field and applying the same value to all field sizes can lead to errors of up to 10% at both 6 MV and 15 MV.

Algorithms↗

PubMed2005

Monte Carlo techniques for scattering foil design and dosimetry in total skin electron irradiations.

Total skin electron irradiation (TSEI) with single fields requires large electron beams having good dose uniformity, dmax at the skin surface, and low bremsstrahlung contamination. To satisfy these requirements, energy degraders and scattering foils have to be specially designed for the given accelerator and treatment room. We used Monte Carlo (MC) techniques based on EGS4 user codes (BEAM, DOSXYZ, and DOSRZ) as a guide in the beam modifier design of our TSEI system. The dosimetric characteristics at the treatment distance of 382 cm source-to-surface distance (SSD) were verified experimentally using a linear array of 47 ion chambers, a parallel plate chamber, and radiochromic film. By matching MC simulations to standard beam measurements at 100 cm SSD, the parameters of the electron beam incident on the vacuum window were determined. Best match was achieved assuming that electrons were monoenergetic at 6.72 MeV, parallel, and distributed in a circular pattern having a Gaussian radial distribution with full width at half maximum = 0.13 cm. These parameters were then used to simulate our TSEI unit with various scattering foils. Two of the foils were fabricated and experimentally evaluated by measuring off-axis dose uniformity and depth doses. A scattering foil, consisting of a 12 x 12 cm2 aluminum plate of 0.6 cm thickness and placed at isocenter perpendicular to the beam direction, was considered optimal. It produced a beam that was flat within +/-3% up to 60 cm off-axis distance, dropped by not more than 8% at a distance of 90 cm, and had an x-ray contamination of <3%. For stationary beams, MC-computed dmax, Rp, and R50 agreed with measurements within 0.5 mm. The MC-predicted surface dose of the rotating phantom was 41% of the dose rate at dmax of the stationary phantom, whereas our calculations based on a semiempirical formula in the literature yielded a drop to 42%. The MC simulations provided the guideline of beam modifier design for TSEI and estimated the dosimetric performance for stationary and rotational irradiations.

Calibration↗

PubMed2004

A dynamic supraclavicular field-matching technique for head-and-neck cancer patients treated with IMRT.

PURPOSE: The conventional single-isocenter and half-beam (SIHB) technique for matching supraclavicular fields with head-and-neck (HN) intensity-modulated radiotherapy (IMRT) fields is subject to substantial dose inhomogeneities from imperfect accelerator jaw/MLC calibration. It also limits the isocenter location and restricts the useful field size for IMRT. We propose a dynamic field-matching technique to overcome these limitations. METHODS AND MATERIALS: The proposed dynamic field-matching technique makes use of wedge junctions for the abutment of supraclavicular and HN IMRT fields. The supraclavicular field was shaped with a multileaf collimator (MLC), which was orientated such that the leaves traveled along the superoinferior direction. The leaves that defined the superior field border moved continuously during treatment from 1.5 cm below to 1.5 cm above the conventional match line to generate a 3-cm-wide wedge-shaped junction. The HN IMRT fields were optimized by taking into account the dose contribution from the supraclavicular field to the junction area, which generates a complementary wedge to produce a smooth junction in the abutment region. This technique was evaluated on a polystyrene phantom and 10 HN cancer patients. Treatment plans were generated for the phantom and the 10 patients. Dose profiles across the abutment region were measured in the phantom on films. For patient plans, dose profiles that passed through the center of the neck lymph nodes were calculated using the proposed technique and the SIHB technique, and dose uniformity in the abutment region was compared. Field mismatches of +/- 1 mm and +/- 2 mm because of imperfect jaw/MLC calibration were simulated, and the resulting dose inhomogeneities were studied for the two techniques with film measurements and patient plans. Three-dimensional volumetric doses were analyzed, and equivalent uniform doses (EUD) were computed. The effect of field mismatches on EUD was compared for the two match techniques. RESULTS: For a perfect jaw/MLC calibration, dose profiles for the 10 patients in the 3-cm match zone had an average inhomogeneity range of -1.6% to +1.6% using the dynamic-matching technique and -3.7% to +3.8% according to the SIHB technique. Measurements showed that dose inhomogeneities that resulted from 1-mm and 2-mm jaw/MLC calibration errors were reduced from as large as 27% and 45% with the SIHB technique to less than 2% and 5.7% with the dynamic technique, respectively. For -1-mm, -2-mm, +1-mm, and +2-mm jaw/MLC calibration errors, respectively, treatment plans for the 10 patients yielded average dose inhomogeneities of -5.9%, -3.0%, +2.7%, and +5.8% with the dynamic technique as compared to -22.8%, -11.1%, +9.8%, and +22.1% with the SIHB technique. Calculation based on a dose-volume histogram (DVH) showed that the SIHB technique resulted in larger changes in EUD of the PTV in the junction area than did the dynamic technique. CONCLUSION: Compared with the conventional SIHB technique, the dynamic field-matching technique provides superior dose homogeneity in the abutment region between the supraclavicular and HN IMRT fields. The dynamic feathering mechanism substantially reduces dose inhomogeneities that result from imperfect jaw/MLC calibration. In addition, isocenter location in the dynamic field-matching technique can be chosen for reproducible patient setup and for adequate IMRT field size rather than being dictated by the match position. It also allows angling of the supraclavicular field to reduce the volume of healthy lung irradiated, which is impractical with the SIHB technique. In principle, this technique should be applicable to any treatment site that requires the abutment of static and intensity-modulated fields.

Calibration↗

PubMed2004

Dose errors due to inhomogeneities in balloon catheter brachytherapy for breast cancer.

PURPOSE: To evaluate dose errors in balloon catheter brachytherapy of breast cancer due to inhomogeneities, such as iodine-containing radiographic contrast medium in the balloon, the lack of scattering medium, and the low density of lung that are not considered by commercial treatment planning systems (TPS). METHODS AND MATERIALS: By accounting for these inhomogeneities in breast/lung phantoms, Monte Carlo simulations were performed to calculate doses in the breast and lung. Doses were also calculated by a commercial TPS. The Monte Carlo doses and the TPS doses were compared along the transverse and longitudinal axes of the source. RESULTS: The Monte Carlo doses were lower by 4-10% on the prescription line than the TPS doses, depending on the concentration (5-25% by volume) of the contrast medium, and on the direction from the source. The lack of scattering medium around the breast contributes to the differences more than the attenuation by the contrast medium. Attenuation contributed approximately 1.0-4.8% at the concentrations investigated in this study. CONCLUSIONS: Current treatment planning systems, which assume a source in a large homogeneous water-equivalent medium, significantly overestimate doses in breast brachytherapy.

Brachytherapy↗

PubMed2004

Benchmark of PENELOPE code for low-energy photon transport: dose comparisons with MCNP4 and EGS4.

The expanding clinical use of low-energy photon emitting 125I and 103Pd seeds in recent years has led to renewed interest in their dosimetric properties. Numerous papers pointed out that higher accuracy could be obtained in Monte Carlo simulations by utilizing newer libraries for the low-energy photon cross-sections, such as XCOM and EPDL97. The recently developed PENELOPE 2001 Monte Carlo code is user friendly and incorporates photon cross-section data from the EPDL97. The code has been verified for clinical dosimetry of high-energy electron and photon beams, but has not yet been tested at low energies. In the present work, we have benchmarked the PENELOPE code for 10-150 keV photons. We computed radial dose distributions from 0 to 10 cm in water at photon energies of 10-150 keV using both PENELOPE and MCNP4C with either DLC-146 or DLC-200 cross-section libraries, assuming a point source located at the centre of a 30 cm diameter and 20 cm length cylinder. Throughout the energy range of simulated photons (except for 10 keV), PENELOPE agreed within statistical uncertainties (at worst +/- 5%) with MCNP/DLC-146 in the entire region of 1-10 cm and with published EGS4 data up to 5 cm. The dose at 1 cm (or dose rate constant) of PENELOPE agreed with MCNP/DLC-146 and EGS4 data within approximately +/- 2% in the range of 20-150 keV, while MCNP/DLC-200 produced values up to 9% lower in the range of 20-100 keV than PENELOPE or the other codes. However, the differences among the four datasets became negligible above 100 keV.

Benchmarking↗

PubMed2004

Attenuation of intracavitary applicators in 192Ir-HDR brachytherapy.

Unlike the penetrating monoenergetic 662 keV gamma rays emitted by 137Cs LDR sources, the spectrum of 192Ir used in HDR brachytherapy contains low-energy components. Since these are selectively absorbed by the high-atomic number materials of which intracavitary applicators are made, the traditional neglect of applicator attenuation can lead to appreciable dose errors. We investigated the attenuation effects of a uterine applicator, and of a set of commonly used vaginal cylinders. The uterine applicator consists of a stainless steel source guide tube with a wall thickness of 0.5 mm and a density of 8.02 g/cm3, whereas the vaginal cylinders consist of the same stainless steel tube plus concentric polysulfone cylinders with a radius of 1 or 2 cm and a density of 1.40 g/cm3. Monte Carlo simulations were performed to compute dose distributions for a bare 192Ir-HDR source, and for the same source located within the applicators. Relative measurements of applicator attenuation using ion-chambers (0.125 cm3) confirmed the Monte Carlo results within 0.5%. We found that the neglect of the applicator attenuation overestimates the dose along the transverse plane by up to 3.5%. At oblique angles, the longer photon path within applicators worsens the error. We defined attenuation-corrected radial dose and anisotropy functions, and applied them to a treatment having multiple dwell positions inside a vaginal cylinder.

Algorithms↗

PubMed2003

Dosimetric characteristics of a linear array of gamma or beta-emitting seeds in intravascular irradiation: Monte Carlo studies for the AAPM TG-43/60 formalism.

In principle, the AAPM TG-43/60 formalism for intravascular brachytherapy (IVBT) dosimetry of catheter-based sources is fully valid with a single seed of cylindrical symmetry and in the region comparable to or larger than the mean-free path of emitting radiation. However, for the geometry of a linear array of seeds within the few millimeter range of interest in IVBT, the suitability of the AAPM TG-43/60 formalism has not been fully addressed yet. We have meticulously investigated the dosimetric characteristics of catheter-based gamma (192Ir) and beta (90Sr/Y) sources using Monte Carlo methods before applying the AAPM TG-43/60 formalism. The dosimetric perturbation due to radiation interactions with neighboring seeds is at most 2% over the entire region of interest for the 192Ir source, while it increases to about 5% for the 90Sr/Y source. As the transaxial distance (y) increases beyond 3 mm, the sum of the dose contributions from neighboring seeds exceeds the dose contribution from the center seed for both sources. However, it continues to increase with the increasing y for 192Ir but is saturated beyond y = 5 mm for 9Sr/Y. Even within a few millimeters from the seeds, the dose from the low-energy betas of 192Ir is still less than 1% of the total dose. The radial dose and anisotropy functions are reformulated in reduced cylindrical coordinate with the reference point at y = 2 mm. The dose rate constant of 192Ir and the dose rate of 90Sr/Y at the reference point showed a fairly good agreement (within +/- 2%) with earlier studies and the NIST-traceable value, respectively. We conclude that the dosimetric perturbation caused by close proximity of neighboring seeds is nearly negligible so that the AAPM TG-43/60 formalism can be applied to a linear array of seeds.

Alpha Particles↗