PubMed Health⌕ Search

Biomedical subjects

S B Jiang

Publications and source records attributed to S B Jiang.

At least 19 recordsLinked to original sources

Effect of statistical uncertainties on Monte Carlo treatment planning.

This paper reviews the effect of statistical uncertainties on radiotherapy treatment planning using Monte Carlo simulations. We discuss issues related to the statistical analysis of Monte Carlo dose calculations for realistic clinical beams using various variance reduction or time saving techniques. We discuss the effect of statistical uncertainties on dose prescription and monitor unit calculation for conventional treatment and intensity-modulated radiotherapy (IMRT) based on Monte Carlo simulations. We show the effect of statistical uncertainties on beamlet dose calculation and plan optimization for IMRT and other advanced treatment techniques such as modulated electron radiotherapy (MERT). We provide practical guidelines for the clinical implementation of Monte Carlo treatment planning and show realistic examples of Monte Carlo based IMRT and MERT plans.

Computer Simulation↗

A quality assurance phantom for IMRT dose verification.

This paper investigates a quality assurance (QA) phantom specially designed to verify the accuracy of dose distributions and monitor units (MU) calculated by clinical treatment planning optimization systems and by the Monte Carlo method for intensity-modulated radiotherapy (IMRT). The QA phantom is a PMMA cylinder of 30 cm diameter and 40 cm length with various bone and lung inserts. A procedure (and formalism) has been developed to measure the absolute dose to water in the PMMA phantom. Another cylindrical phantom of the same dimensions, but made of water, was used to confirm the results obtained with the PMMA phantom. The PMMA phantom was irradiated by 4, 6 and 15 MV photon beams and the dose was measured using an ionization chamber and compared to the results calculated by a commercial inverse planning system (CORVUS, NOMOS, Sewickley, PA) and by the Monte Carlo method. The results show that the dose distributions calculated by both CORVUS and Monte Carlo agreed to within 2% of dose maximum with measured results in the uniform PMMA phantom for both open and intensity-modulated fields. Similar agreement was obtained between Monte Carlo calculations and measured results with the bone and lung heterogeneity inside the PMMA phantom while the CORVUS results were 4% different. The QA phantom has been integrated as a routine QA procedure for the patient's IMRT dose verification at Stanford since 1999.

Equipment Design↗

A Monte Carlo dose calculation tool for radiotherapy treatment planning.

A Monte Carlo user code, MCDOSE, has been developed for radiotherapy treatment planning (RTP) dose calculations. MCDOSE is designed as a dose calculation module suitable for adaptation to host RTP systems. MCDOSE can be used for both conventional photon/electron beam calculation and intensity modulated radiotherapy (IMRT) treatment planning. MCDOSE uses a multiple-source model to reconstruct the treatment beam phase space. Based on Monte Carlo simulated or measured beam data acquired during commissioning, source-model parameters are adjusted through an automated procedure. Beam modifiers such as jaws, physical and dynamic wedges, compensators, blocks, electron cut-outs and bolus are simulated by MCDOSE together with a 3D rectilinear patient geometry model built from CT data. Dose distributions calculated using MCDOSE agreed well with those calculated by the EGS4/DOSXYZ code using different beam set-ups and beam modifiers. Heterogeneity correction factors for layered-lung or layered-bone phantoms as calculated by both codes were consistent with measured data to within 1%. The effect of energy cut-offs for particle transport was investigated. Variance reduction techniques were implemented in MCDOSE to achieve a speedup factor of 10-30 compared to DOSXYZ.

Computer Simulation↗

The MLC tongue-and-groove effect on IMRT dose distributions.

We have investigated the tongue-and-groove effect on the IMRT dose distributions for a Varian MLC. We have compared the dose distributions calculated using the intensity maps with and without the tongue-and-groove effect. Our results showed that, for one intensity-modulated treatment field, the maximum tongue-and-groove effect could be up to 10% of the maximum dose in the dose distributions. For an IMRT treatment with multiple gantry angles (> or = 5), the difference between the dose distributions with and without the tongue-and-groove effect was hardly visible, less than 1.6% for the two typical clinical cases studied. After considering the patient setup errors, the dose distributions were smoothed with reduced and insignificant differences between plans with and without the tongue-and-groove effect. Therefore, for a multiple-field IMRT plan (> or = 5), the tongue-and-groove effect on the IMRT dose distributions will be generally clinically insignificant due to the smearing effect of individual fields. The tongue-and-groove effect on an IMRT plan with small number of fields (< 5) will vary depending on the number of fields in a plan (coplanar or non-coplanar), the MLC leaf sequences and the patient setup uncertainty, and may be significant (> 5% of maximum dose) in some cases, especially when the patient setup uncertainty is small (< or = 2 mm).

Bone Neoplasms↗

Monte Carlo simulation of electron beams from an accelerator head using PENELOPE.

The Monte Carlo code PENELOPE has been used to simulate electron beams from a Siemens Mevatron KDS linac with nominal energies of 6, 12 and 18 MeV. Owing to its accuracy, which stems from that of the underlying physical interaction models, PENELOPE is suitable for simulating problems of interest to the medical physics community. It includes a geometry package that allows the definition of complex quadric geometries, such as those of irradiation instruments, in a straightforward manner. Dose distributions in water simulated with PENELOPE agree well with experimental measurements using a silicon detector and a monitoring ionization chamber. Insertion of a lead slab in the incident beam at the surface of the water phantom produces sharp variations in the dose distributions, which are correctly reproduced by the simulation code. Results from PENELOPE are also compared with those of equivalent simulations with the EGS4-based user codes BEAM and DOSXYZ. Angular and energy distributions of electrons and photons in the phase-space plane (at the downstream end of the applicator) obtained from both simulation codes are similar, although significant differences do appear in some cases. These differences, however, are shown to have a negligible effect on the calculated dose distributions. Various practical aspects of the simulations, such as the calculation of statistical uncertainties and the effect of the 'latent' variance in the phase-space file, are discussed in detail.

Dose-Response Relationship, Radiation↗

Derivation of electron and photon energy spectra from electron beam central axis depth dose curves.

A method for deriving the electron and photon energy spectra from electron beam central axis percentage depth dose (PDD) curves has been investigated. The PDD curves of 6, 12 and 20 MeV electron beams obtained from the Monte Carlo full phase space simulations of the Varian linear accelerator treatment head have been used to test the method. We have employed a 'random creep' algorithm to determine the energy spectra of electrons and photons in a clinical electron beam. The fitted electron and photon energy spectra have been compared with the corresponding spectra obtained from the Monte Carlo full phase space simulations. Our fitted energy spectra are in good agreement with the Monte Carlo simulated spectra in terms of peak location, peak width, amplitude and smoothness of the spectrum. In addition, the derived depth dose curves of head-generated photons agree well in both shape and amplitude with those calculated using the full phase space data. The central axis depth dose curves and dose profiles at various depths have been compared using an automated electron beam commissioning procedure. The comparison has demonstrated that our method is capable of deriving the energy spectra for the Varian accelerator electron beams investigated. We have implemented this method in the electron beam commissioning procedure for Monte Carlo electron beam dose calculations.

Algorithms↗

Monte Carlo based treatment planning for modulated electron beam radiation therapy.

A Monte Carlo based treatment planning system for modulated electron radiation therapy (MERT) is presented. This new variation of intensity modulated radiation therapy (IMRT) utilizes an electron multileaf collimator (eMLC) to deliver non-uniform intensity maps at several electron energies. In this way, conformal dose distributions are delivered to irregular targets located a few centimetres below the surface while sparing deeper-lying normal anatomy. Planning for MERT begins with Monte Carlo generation of electron beamlets. Electrons are transported with proper in-air scattering and the dose is tallied in the phantom for each beamlet. An optimized beamlet plan may be calculated using inverse-planning methods. Step-and-shoot leaf sequences are generated for the intensity maps and dose distributions recalculated using Monte Carlo simulations. Here, scatter and leakage from the leaves are properly accounted for by transporting electrons through the eMLC geometry. The weights for the segments of the plan are re-optimized with the leaf positions fixed and bremsstrahlung leakage and electron scatter doses included. This optimization gives the final optimized plan. It is shown that a significant portion of the calculation time is spent transporting particles in the leaves. However, this is necessary since optimizing segment weights based on a model in which leaf transport is ignored results in an improperly optimized plan with overdosing of target and critical structures. A method of rapidly calculating the bremsstrahlung contribution is presented and shown to be an efficient solution to this problem. A homogeneous model target and a 2D breast plan are presented. The potential use of this tool in clinical planning is discussed.

Electrons↗

Modeling the extrafocal radiation and monitor chamber backscatter for photon beam dose calculation.

A simple analytical approach has been developed to model extrafocal radiation and monitor chamber backscatter for clinical photon beams. Model parameters for both the extrafocal source and monitor chamber backscatter are determined simultaneously using conventional measured data, i.e., in-air output factors for square and rectangular fields defined by the photon jaws. The model has been applied to 6 MV and 15 MV photon beams produced by a Varian Clinac 2300C/D accelerator. Contributions to the in-air output factor from the extrafocal radiation and monitor chamber backscatter, as predicted by the model, are in good agreement with the measurements. The model can be used to calculate the in-air output factors analytically, with an accuracy of 0.2% for symmetric or asymmetric rectangular fields defined by jaws when the calculation point is at the isocenter and 0.5% when the calculation point is at an extended SSD. For MLC-defined fields, with the jaws at the recommended positions, calculated in-air output factors agree with the measured data to within 0.3% at the isocenter and 0.7% at off-axis positions. The model has been incorporated into a Monte Carlo dose algorithm to calculate the absolute dose distributions in patients or phantoms. For three MLC-defined irregular fields (triangle shape, C-shape, and L-shape), the calculations agree with the measurements to about 1% even for points at off-axis positions. The model will be particularly useful for IMRT dose calculations because it accurately predicts beam output and penumbra dose.

Biophysical Phenomena↗

Er:YAG laser for the treatment of actinic keratoses.

BACKGROUND: There is no single optimal treatment for multiple facial actinic keratoses. The existing therapies such as topical 5-fluorouracil, chemical peels, cryotherapy, dermabrasion, and CO2 laser resurfacing can produce prolonged recovery time or are often operator dependent. OBJECTIVE: The purpose of this study was to investigate another therapeutic modality which provides a shorter recovery time with uniform results. We performed a prospective pilot study investigating the use of the Er:YAG laser for the treatment of multiple facial actinic keratoses. METHODS: Five patients with multiple facial actinic keratoses were treated with two to three passes of Er:YAG laser. Anesthesia was achieved in all cases by topical application and local infiltration when indicated. All patients were treated with 2.0 J, 5 mm spot size, and a fluence of 10 J/cm2. Clinical and histologic evaluations were performed both pre- and postoperatively. RESULTS: All patients showed a decrease in the total number of clinical actinic keratoses on the face ranging from 86 to 96%. In addition to the reversal of actinic damage in the epidermis, histologic evidence revealed increased fibroplasia and decreased superficial solar elastosis 3 months after the laser resurfacing. Reepithelialization occurred in 5-8 days, and erythema lasted for about 3-6 weeks after the procedure. There was no evidence of scarring or pigmentary changes in any of the patients during the follow-up period. CONCLUSION: Er:YAG laser skin resurfacing is a safe and effective treatment for multiple facial actinic keratoses. Histologic data suggest a new zone of collagen deposition occurs in the superficial papillary dermis. Under our current parameters, Er:YAG laser skin resurfacing has a relatively short recovery period and a low risk of scarring. Unlike the CO2 laser, Er:YAG laser skin resurfacing can be performed with topical anesthesia alone.

Aged↗

Validation of a Monte Carlo dose calculation tool for radiotherapy treatment planning.

A new EGS4/PRESTA Monte Carlo user code, MCDOSE, has been developed as a routine dose calculation tool for radiotherapy treatment planning. It is suitable for both conventional and intensity modulated radiation therapy. Two important features of MCDOSE are the inclusion of beam modifiers in the patient simulation and the implementation of several variance reduction techniques. Before this tool can be used reliably for clinical dose calculation, it must be properly validated. The validation for beam modifiers has been performed by comparing the dose distributions calculated by MCDOSE and the well-benchmarked EGS4 user codes BEAM and DOSXYZ. Various beam modifiers were simulated. Good agreement in the dose distributions was observed. The differences in electron cutout factors between the results of MCDOSE and measurements were within 2%. The accuracy of MCDOSE with various variance reduction techniques was tested by comparing the dose distributions in different inhomogeneous phantoms with those calculated by DOSXYZ without variance reduction. The agreement was within 1.0%. Our results demonstrate that MCDOSE is accurate and efficient for routine dose calculation in radiotherapy treatment planning, with or without beam modifiers.

Bone and Bones↗

Photon beam characterization and modelling for Monte Carlo treatment planning.

Photon beams of 4, 6 and 15 MV from Varian Clinac 2100C and 2300C/D accelerators were simulated using the EGS4/BEAM code system. The accelerators were modelled as a combination of component modules (CMs) consisting of a target, primary collimator, exit window, flattening filter, monitor chamber, secondary collimator, ring collimator, photon jaws and protection window. A full phase space file was scored directly above the upper photon jaws and analysed using beam data processing software, BEAMDP, to derive the beam characteristics, such as planar fluence, angular distribution, energy spectrum and the fractional contributions of each individual CM. A multiple-source model has been further developed to reconstruct the original phase space. Separate sources were created with accurate source intensity, energy, fluence and angular distributions for the target, primary collimator and flattening filter. Good agreement (within 2%) between the Monte Carlo calculations with the source model and those with the original phase space was achieved in the dose distributions for field sizes of 4 cm x 4 cm to 40 cm x 40 cm at source surface distances (SSDs) of 80-120 cm. The dose distributions in lung and bone heterogeneous phantoms have also been found to be in good agreement (within 2%) for 4, 6 and 15 MV photon beams for various field sizes between the Monte Carlo calculations with the source model and those with the original phase space.

Models, Theoretical↗

Removing the effect of statistical uncertainty on dose-volume histograms from Monte Carlo dose calculations.

Dose-volume histograms (DVHs) of the dose distributions calculated by the Monte Carlo method contain statistical uncertainties. The Monte Carlo DVH can be considered as blurred from the noiseless DVH by the statistical uncertainty. The focus of the present work is on the removal of the statistical uncertainty effect on the Monte Carlo DVHs and the reconstruction of the noiseless DVHs. We first study the effect of statistical uncertainty. It is found that the steeper the DVH, the more significant the effect. For typical critical structure DVHs the effect is usually negligible. For the target DVHs the effect could be clinically significant, depending on the value of uncertainty and the slope of the DVH. We then propose an iterative reconstruction algorithm. Using the DVHs and statistical uncertainties from the Monte Carlo simulations, we are able to reconstruct the noiseless DVHs. A hypothetical example and a number of clinical cases have been used to test the proposed algorithm. For each clinical case, two Monte Carlo simulations (denoted A and B) were performed. Simulation A has very large statistical uncertainties (about 10% of dose in the target volume) while simulation B has very small uncertainties (about 1%). DVHs from simulation B were used to approximate the noiseless DVHs. Using the proposed algorithm, the effect of statistical uncertainty can be removed from the DVHs of simulation A. The reconstructed DVHs were in good agreement with the DVHs from simulation B. The proposed approach is expected to be useful in removing the blurring effect on a quickly calculated Monte Carlo DVH when performing the iterative forward treatment planning.

Algorithms↗

Energy- and intensity-modulated electron beams for radiotherapy.

This work investigates the feasibility of optimizing energy- and intensity-modulated electron beams for radiation therapy. A multileaf collimator (MLC) specially designed for modulated electron radiotherapy (MERT) was investigated both experimentally and by Monte Carlo simulations. An inverse-planning system based on Monte Carlo dose calculations was developed to optimize electron beam energy and intensity to achieve dose conformity for target volumes near the surface. The results showed that an MLC with 5 mm leaf widths could produce complex field shapes for MERT. Electron intra- and inter-leaf leakage had negligible effects on the dose distributions delivered with the MLC, even at shallow depths. Focused leaf ends reduced the electron scattering contributions to the dose compared with straight leaf ends. As anticipated, moving the MLC position toward the patient surface reduced the penumbra significantly. There were significant differences in the beamlet distributions calculated by an analytic 3-D pencil beam algorithm and the Monte Carlo method. The Monte Carlo calculated beamlet distributions were essential to the accuracy of the MERT dose distribution in cases involving large air gaps, oblique incidence and heterogeneous treatment targets (at the tissue-bone and bone-lung interfaces). To demonstrate the potential of MERT for target dose coverage and normal tissue sparing for treatment of superficial targets, treatment plans for a hypothetical treatment were compared using photon beams and MERT.

Algorithms↗

Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system.

The purpose of this work was to use Monte Carlo simulations to verify the accuracy of the dose distributions from a commercial treatment planning optimization system (Corvus, Nomos Corp., Sewickley, PA) for intensity-modulated radiotherapy (IMRT). A Monte Carlo treatment planning system has been implemented clinically to improve and verify the accuracy of radiotherapy dose calculations. Further modifications to the system were made to compute the dose in a patient for multiple fixed-gantry IMRT fields. The dose distributions in the experimental phantoms and in the patients were calculated and used to verify the optimized treatment plans generated by the Corvus system. The Monte Carlo calculated IMRT dose distributions agreed with the measurements to within 2% of the maximum dose for all the beam energies and field sizes for both the homogeneous and heterogeneous phantoms. The dose distributions predicted by the Corvus system, which employs a finite-size pencil beam (FSPB) algorithm, agreed with the Monte Carlo simulations and measurements to within 4% in a cylindrical water phantom with various hypothetical target shapes. Discrepancies of more than 5% (relative to the prescribed target dose) in the target region and over 20% in the critical structures were found in some IMRT patient calculations. The FSPB algorithm as implemented in the Corvus system is adequate for homogeneous phantoms (such as prostate) but may result in significant under or over-estimation of the dose in some cases involving heterogeneities such as the air-tissue, lung-tissue and tissue-bone interfaces.

Computer Simulation↗

Monte Carlo and experimental investigations of multileaf collimated electron beams for modulated electron radiation therapy.

Modulated electron radiation therapy (MERT) has been proposed as a means of delivering conformal dose to shallow tumors while sparing distal structures and surrounding tissues. Conventional systems for electron beam collimation are labor and time intensive in their construction and are therefore inadequate for use in the sequential delivery of multiple complex fields required by MERT. This study investigates two proposed methods of electron beam collimation: the use of existing photon multileaf collimators (MLC) in a helium atmosphere to reduce in-air electron scatter, and a MLC specifically designed for electron beam collimation. Monte Carlo simulations of a Varian Clinac 2100C were performed using the EGS4/BEAM system and dose calculations performed with the MCDOSE code. Dose penumbras from fields collimated by photon MLCs both with air and with helium at 6, 12, and 20 MeV at a range of SSDs from 70 to 90 cm were examined. Significant improvements were observed for the helium based system. Simulations were also performed on an electron specific MLC located at the level of the last scraper of a 25x25 cm2 applicator. A number of leaf materials, thicknesses, end shapes, and widths were simulated to determine optimal construction parameters. The results demonstrated that tungsten leaves 15 mm thick and 5 mm wide with unfocused ends would provide sufficient collimation for MERT fields. A prototype electron MLC was constructed and comparisons between film measurements and simulation demonstrate the validity of the Monte Carlo model. Further simulations of dose penumbras demonstrate that such an electron MLC would provide improvements over the helium filled photon MLC at all energies, and improvements in the 90-10 penumbra of 12% to 45% at 20 MeV and 6 MeV, respectively. These improvements were also seen in isodose curves when a complex field shape was simulated. It is thus concluded that an MLC specific for electron beam collimation is required for MERT.

Air↗

Electron beam modeling and commissioning for Monte Carlo treatment planning.

A hybrid approach for commissioning electron beam Monte Carlo treatment planning systems has been studied. The approach is based on the assumption that accelerators of the same type have very similar electron beam characteristics and the major difference comes from the on-site tuning of the electron incident energy at the exit window. For one type of accelerator, a reference machine can be selected and simulated with the Monte Carlo method. A multiple source model can be built on the full Monte Carlo simulation of the reference beam. When commissioning electron beams from other accelerators of the same type, the energy spectra in the source model are tuned to match the measured dose distributions. A Varian Clinac 2100C accelerator was chosen as the reference machine and a four-source beam model was established based on the Monte Carlo simulations. This simplified beam model can be used to generate Monte Carlo dose distributions accurately (within 2%/2 mm compared to those calculated with full phase space data) for electron beams from the reference machine with various nominal energies, applicator sizes, and SSDs. Three electron beams were commissioned by adjusting the energy spectra in the source model. The dose distributions calculated with the adjusted source model were compared with the dose distributions calculated using the phase space data for these beams. The agreement is within 1% in most of cases and 2% in all situations. This preliminary study has shown the capability of the commissioning approach for handling large variation in the electron incident energy. The possibility of making the approach more versatile is also discussed.

Algorithms↗

Microcystic adnexal carcinoma: collaborative series review and update.

BACKGROUND: Microcystic adnexal carcinoma (MAC) is a malignant appendageal tumor first described in 1982. It can be clinically and histologically confused with other malignant and benign cutaneous neoplasms, leading to inadequate initial treatment. This neoplasm is locally aggressive and deeply infiltrating, characterized by high morbidity and frequent recurrence. Mohs micrographic surgery has been used to conserve tissue and improve the likelihood for cure. OBJECTIVE: We report our experience using Mohs micrographic surgery for the treatment of MAC and compare with earlier reports in the literature. In addition, we review the epidemiology, clinical and histologic characteristics, and optimal treatment of this rare neoplasm. We also describe a 15-year-old white male patient with MAC on the scalp occurring only 7 years after radiation exposure. METHODS: The medical records of 11 patients with MAC who were treated by Mohs micrographic surgery were reviewed at both departments, and follow-up data were obtained. RESULTS: In all patients treated with Mohs micrographic surgery, there were no recurrences after a mean follow-up of 5 years. CONCLUSION: Mohs technique enables the detection of clinically unrecognizable tumor spread and perineural invasion often encountered with MAC. Aggressive initial treatment by microscopically controlled excision appears to offer the greatest likelihood of cure for this neoplasm, while providing conservation of normal tissue. In addition, we describe the second youngest patient with MAC and readdress the issue of previous radiotherapy as an important predisposing factor.

Adolescent↗

Monte Carlo modelling of electron beams from medical accelerators.

Monte Carlo simulation of radiation transport is considered to be one of the most accurate methods of radiation therapy dose calculation. With the rapid development of computer technology, Monte Carlo based treatment planning for radiation therapy is becoming practical. A basic requirement for Monte Carlo treatment planning is a detailed knowledge of the radiation beams from medical accelerators. A practical approach to obtain the above is to perform Monte Carlo simulation of radiation transport in the medical accelerator. Additionally, Monte Carlo modelling of the treatment machine head can also improve our understanding of clinical beam characteristics, help accelerator design and improve the accuracy of clinical dosimetry by providing more realistic beam data. This paper summarizes work over the past two decades on Monte Carlo simulation of clinical electron beams from medical accelerators.

Electrons↗