Biomedical subjects

C M Ma

Publications and source records attributed to C M Ma.

At least 55 records · Page 3Linked to original sources

PubMed2000

Calculation of x-ray transmission through a multileaf collimator.

A ray tracing based method has been developed to calculate the x-ray transmission through a multileaf collimator (MLC) for beam delivery verification and dose calculation in intensity modulated radiotherapy (IMRT). The path length of a ray line in the MLC is accurately calculated using the exact geometry of the MLC leaves. The fluence distribution of an IMRT field is calculated first using a point source. The fluence distribution for a realistic beam model is obtained, as an approximation, by convolving the point source fluence distribution with the distribution of source strength. Full ray tracing calculations are performed using analytic and Monte Carlo simulated beam models to verify the accuracy of the convolution method. The calculation is in better agreement with measurements using either film or a beam imaging system (BIS) than previous calculations for MLC transmission using a simplified model. This ray tracing calculation can be applied to the problem of verifying dynamic MLC leaf sequences as part of a patient-specific quality assurance process for IMRT.

Algorithms↗

PubMed2000

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↗

PubMed2000

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↗

PubMed2000

Biotransformation of a C-glycosylflavone, abrusin 2''-O-beta-D-apioside, by human intestinal bacteria.

After anaerobic incubation of abrusin 2''-O-beta-D-apioside (1) with a human fecal suspension, five metabolites were isolated and identified as abrusin (2), 1-(2',6'-dihydroxy-3',4'-dimethoxyphenyl)-3-(4''-hydroxyphenyl)propan-1- one (5), 5,6-dimethoxybenzene-1,3-diol (6), 3-(4'-hydroxyphenyl)propionic acid (7) and 3-phenylpropionic acid (8). However, methyl ether derivatives of abrusin (4'-O-methylabrusin and 4'-O-, 5-O-dimethylabrusin) resisted degradation under the same conditions.

Biotransformation↗

PubMed2000

Four new saponins from the root bark of Aralia elata.

Four new saponins, 3-O-[beta-D-glucopyranosyl(1-->3)-alpha-L-arabinopyranosyl]-16a lpha-hydroxyoleanolic acid 28-O-beta-D-glucopyranosyl ester (called aralia-saponin I), 3-O-[beta-D-glucopyranosyl(1-->3)-alpha-L-arabinopyranosyl]-16a lpha-hydroxyhederagenin 28-O-beta-D-glucopyranosyl ester (aralia-saponin II), 3-O-[beta-D-glucopyranosyl(1-->3)-beta-D-glucopyranosyl(1-->3)-alpha-L-+ ++arabinopyranosyl]-16alpha-hydroxyoleanolic acid 28-O-beta-D-glucopyranosyl ester (aralia-saponin III), 3-O-[beta-D-glucopyranosyl(1-->3)-beta-D-gucopyranosyl(1-->3)-beta -D-glucucopyranosyl]-16alpha-hydroxyoleanolic acid 28-O-beta-D-glucopyranosyl ester (aralia-saponin IV), were isolated from the root bark of Aralia elata (Miq.) Seem., together with nineteen known compounds including glycosides of (20S)-protopanaxadiol and (20S)-protopanaxatriol. Their structures were determined on the basis of chemical and spectroscopy methods.

Carbohydrate Conformation↗

PubMed1999

Synchronizing dynamic multileaf collimators for producing two-dimensional intensity-modulated fields with minimum beam delivery time.

PURPOSE: Leaf motion synchronization of dynamic multileaf collimators (DMLC) for intensity-modulated radiotherapy (IMRT) is important in improving dose distribution and reducing "tongue-and-groove" effects for a prescribed intensity profile. Leaf synchronization could also be used in transforming a one-dimensional leaf-setting algorithm into a two-dimensional leaf-setting algorithm. In this work, we aim to develop a generalized leaf synchronization method for delivering IMRT with the minimized beam delivery time and the optimized subfield variations for a leaf-setting sequence. METHODS AND MATERIALS: With the leaf synchronization procedure, all active MLC leaf pairs start and finish off a leaf sequence simultaneously. In this work, the MLC leaf pairs were synchronized under the condition that the resulting leaf sequence produces the desired intensity profile with the minimum beam delivery time. The parameter of the leaf synchronization function was determined through the least-square minimization of the area variations of all subfields within a leaf sequence. The leaf synchronization and optimization procedure were applied and analyzed for clinical relevant intensity profiles for treating the head-and-neck cancer patients using IMRT. RESULTS: The total monitor units and the optimized beam delivery time of generating a two-dimensional intensity profile was proven through this work to be the global minimum of all leaf-setting sequences including the unsynchronized leaf-setting sequences. The optimized parameter for subfield variations of the synchronized leaf trajectories was found to be dependent on individual intensity profiles. For all our studied cases, the unsynchronized leaf trajectories always have significantly larger subfield variations than the synchronized leaf trajectories. CONCLUSION: It is important and also feasible to synchronize and optimize dynamic MLC leaf motions while still keeping the total beam delivery time minimum for delivering arbitrary two-dimensional intensity-modulated fields.

Algorithms↗

PubMed1999

Anti-HIV-1 protease triterpenoid saponins from the seeds of Aesculus chinensis.

Eight bioactive triterpenoid saponins (1-8) were isolated from the seeds of Aesculus chinensis, four of which are novel compounds. The major saponins were identified as escin Ia (1), Ib (2), isoescin Ia (3) and Ib (4), while the new compounds were identified as 22alpha-tigloyl-28-acetylprotoaescigenin-3beta-O-¿beta -D-glucopyranos yl (1-2) ¿beta-D-glucopyranosyl (1-4)-beta-D-glucopyranosiduronic acid (escin IVc, 5), 22alpha-angeloyl-28-acetylprotoaescigenin-3beta-O-¿bet a-D-glucopyrano syl (1-2) ¿beta-D-glucopyranosyl (1-4)-beta-D-glucopyranosiduronic acid (escin IVd, 6), 28-tigloylprotoaescigenin-3beta-O-¿beta-D-glucopyranosyl (1-2) ¿beta-D-glucopyranosyl (1-4)-beta-D-glucopyranosiduronic acid (escin IVe, 7), and 28-angeloylprotoaescigenin-3beta-O-¿beta-D-glucopyranosyl (1-2) ¿beta-D-glucopyranosyl (1-4)-beta-D-glucopyranosiduronic acid (escin IVf, 8). The structures were determined by chemical and spectroscopic methods. All the above compounds were evaluated for their inhibitory activity against HIV-1 protease.

Carbohydrate Sequence↗

PubMed1999

Inhibitory effect of triterpenes from Crataegus pinatifida on HIV-I protease.

The methanol extracts of the leaves of Crataegus pinnatifida showed potent inhibitory activities against HIV-1 protease at a concentration of 100 micrograms/ml. The subsequent fractionation and isolation of the extract gave two active compounds. Their structures were identified as uvaol (1) and ursolic acid (2) by spectral data. These active compounds inhibit HIV-1 protease with IC50 values of 5.5 and 8.0 microM, respectively.

Anti-HIV Agents↗

PubMed1999

Mass-energy absorption coefficient and backscatter factor ratios for kilovoltage x-ray beams.

For low-energy (up to 150 kV) x-rays, the ratio of mass-energy absorption coefficients for water to air, (mu(en)/rho)w.air, and the backscatter factor B are used in the conversion of air kerma, measured free-in-air, to water kerma on the surface of a water phantom. For clinical radiotherapy, similar conversion factors are needed for the determination of the absorbed dose to biological tissues on (or near) the surface of a human body. We have computed the mu(en)/rho ratios and B factor ratios for different biological tissues including muscle, soft tissue, lung, skin and bone relative to water. The mu(en)/rho ratios were obtained by integrating the respective mass-energy absorption coefficients over the in-air primary photon spectra. We have also calculated the mu(en)/rho ratios at different depths in a water phantom in order to convert the measured in-phantom water kerma to the absorbed dose to various biological tissues. The EGS4/DOSIMETER Monte Carlo code system has been used for the simulation of the energy fluence at different depths in a water phantom irradiated by a kilovoltage x-ray beam of variable beam quality (HVL: 0.1 mm Al-5 mm Cu), field size and source-surface distance (SSD). The same code was also used in the calculation of the B factor ratios, soft tissue to water and bone to water. The results show that the B factor for bone differs from the B factor for water by up to 20% for a 100 kV beam (HVL: 2.65 mm Al) with a 100 cm2 field. On the other hand, the difference in the B factor between water and soft tissue is insignificant (well within 1% generally). This means that the B factors for water may be directly used to convert the 'in-air' water kerma to surface kerma for human soft tissues.

Humans↗

PubMed1999

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↗

PubMed1999

Stopping-power ratios for clinical electron beams from a scatter-foil linear accelerator.

Restricted mass collision stopping-power ratios for electron beams from a scatter-foil medical linear accelerator (Varian Clinac 2100C) were calculated for various combinations of beams, phantoms and detector materials using the Monte Carlo method. The beams were of nominal energy 6, 12 or 20 MeV, with square dimensions 1 x 1 cm2 to 10 x 10 cm2. They were incident at nominal SSDs of 100 or 120 cm and inclined at 90 degrees or 30 degrees to the surface of homogeneous water phantoms or water phantoms interspersed with layered lung or bone-like materials. The broad beam water-to-air stopping-power ratios were within 1.3% of the AAPM TG21 protocol values and consistent with the results of Ding et al to within 0.2%. On the central axis the stopping-power ratio variations for narrow beams compared with normally incident broad beams were 0.1% or less for water-to-LiF-100, graphite, ferrous sulfate dosimeter solution, polystyrene and PMMA, 0.5% for water-to-silicon and 1% for water-to-air and water-to-photographic-film materials. The transverse variations of the stopping-power ratios were up to 4% for water-to-silicon, 7% for water-to-photographic-film materials and 10% for water-to-air in the penumbral regions (where the dose was 10% of the global dose maximum) at shallow depths compared with the values at the same depths on the central axis. In the inhomogeneous phantoms studied, the stopping-power ratio correction factors varied more significantly for air, followed by photographic materials and silicon, at various depths on the central axis in the heterogeneous regions. For the simple layered phantoms studied, the estimation of the stopping-power ratio correction factors based on the relative electron-density derived effective depth approach yielded results that were within 0.5% of the Monte Carlo derived values for all the detector materials studied.

Electrons↗

PubMed1999

Theoretical considerations of monitor unit calculations for intensity modulated beam treatment planning.

A treatment planning system to compute intensity modulated radiotherapy (IMRT) treatments using inverse planning was investigated. The system was designed to optimize the intensity patterns required to treat a specified target volume with specified normal structure constraints. A beam model that uses the convolution of pencil beams was used to compute the dose distributions. A multileaf collimator leaf-setting sequence intended to produce the intensity pattern was computed along with the monitor units required to deliver each of a number of fixed-gantry modulated fields. Computer calculations are commonly verified using an independent manual procedure. It is difficult to calculate treatment delivery monitor units for this variant of IMRT using manual methods. Since manual calculations are not feasible, it is important both to understand and to verify the calculation of treatment monitor units by the planning system algorithm. A formal analysis was made of the dose calculation model and the monitor unit calculation embedded in the algorithm. Experimental verification of the dose delivered by plans computed with the methodology demonstrated an agreement of better than 4% between the dose model and measurements.

Algorithms↗

PubMed1999

Clinical implementation of a Monte Carlo treatment planning system.

The purpose of this study was to implement the Monte Carlo method for clinical radiotherapy dose calculations. We used the EGS4/BEAM code to obtain the phase-space data for 6-20 MeV electron beams and 4, 6, and 15 MV photon beams for Varian Clinac 1800, 2100C, and 2300CD accelerators. A multiple-source model was used to reconstruct the phase-space data for both electron and photon beams, which retained the accuracy of the Monte Carlo beam data. The multiple-source model reduced the phase-space data storage requirement by a factor of 1000 and the accelerator simulation time by a factor of 10 or more. Agreement within 2% was achieved between the Monte Carlo calculations and measurements of the dose distributions in homogeneous and heterogeneous phantoms for various field sizes, source-surface distances, and beam modulations. The Monte Carlo calculated electron output factors were within 2% of the measured values for various treatment fields while the heterogeneity correction factors for various lung and bone phantoms were within 1% for photon beams and within 2% for electron beams. The EGS4/DOSXYZ Monte Carlo code was used for phantom and patient dose calculations. The results were compared to the dose distributions produced by a conventional treatment planning system and an intensity-modulated radiotherapy inverse-planning system. Significant differences (>5% in dose and >5 mm shift in isodose lines) were found between Monte Carlo calculations and the analytical calculations implemented in the commercial systems. Treatment sites showing the largest dose differences were for head and neck, lung, and breast cases.

Algorithms↗

PubMed1999

Treatment head design for multileaf collimated high-energy electrons.

This paper describes how a conventional treatment head can be modified for use of multileaf collimated electron beams. Automatic and dynamic beam delivery are possible for both electrons and photons by using the computer controlled multileaf collimator (MLC) for both photon and electron beams. Thereby, the electron beams can be mixed more freely into the treatment to take advantage of the specific depth modulation characteristics of electrons. The investigation was based on Monte Carlo calculations using the software package BEAM. The physical parameters used in this optimization were the beam penumbra and the virtual/effective point source position. These parameters are essential for shaping beams, beam matching and for dosimetry calculations. The optimization was carried out by modifying a number of parameters: replacing the air atmosphere in the treatment head with helium, adding a helium bag below the MLC, changing the position of the scattering foils, modifying the monitor chamber, and adjusting the position of the MLC. The beam characteristics for some of these designs were found to fulfil our criteria for clinically useful beams down to at least 9 MeV.

Computer Simulation↗

PubMed1998

Monte Carlo calculations of electron beam output factors for a medical linear accelerator.

The purpose of this study was to investigate the application of the Monte Carlo technique to the calculation and analysis of output factors for electron beams used in radiotherapy. The code EGS4/BEAM was used to obtain phase-space files for 6, 12 and 20 MeV clinical electron beams from a scattering-foil linac (Varian Clinac 2100C) for a clinically representative range of applicator and square or rectangular insert combinations. The source-to-surface distance used was 100 cm. The field sizes ranged from 1 x 1 cm2 to 20 x 20 cm2. These phase-space files were analysed to study the intrinsic beam characteristics and used as source input for relative dose and output factor computations in homogeneous water phantoms using the code EGS4/DOSXYZ. The calculated relative central-axis depth-dose and transverse dose profiles at various depths of clinical interest agreed with the corresponding measured dose profiles to within 2% of the maximum dose. Calculated output factors for the fields studied agreed with measured output factors to about 2%. This demonstrated that for the Varian Clinac 2100C linear accelerator, electron beam dose calculations in homogeneous water phantoms can be performed accurately at the 2% level using Monte Carlo simulations.

Computer Simulation↗

PubMed1998

An optimized leaf-setting algorithm for beam intensity modulation using dynamic multileaf collimators.

A leaf-setting algorithm is developed for generating arbitrary beam intensity profiles in discrete levels using dynamic multileaf collimators (DMLCs). The algorithm starts with the algebraic expression for the area under the beam profile. It is shown that the coefficients in this expression can be transformed into the specifications for the leaf-setting sequence. It is proven that the algorithm optimizes beam delivery time and total monitor units for the DMLC leaf setting for intensity modulated radiotherapy (IMRT). The algorithm is demonstrated to be applicable to both the 'step-and-shoot' and 'dynamic' type of beam delivery. The graphical interpretation and numerical implementation scheme of the algorithm is illustrated using a simplified example.

Algorithms↗

PubMed1998

Dosimetric evaluation of a widely used kilovoltage x-ray unit for endocavitary radiotherapy.

In this paper we present the dosimetric data of a Therapax DTX300 kilovoltage x-ray unit for endocavitary rectal irradiation. The unit if operated at tube voltage of 40-60 kVp (30 mA) with an added filtration of 0.2-0.4 mm Al generates acceptable beam qualities comparable to those of the original Papillon technique. Relative dosimetric measurements were performed at the cone end (37.2 cm SSD) of a 3 cm diameter rectal cone using various detectors to ensure the accuracy. A Monte Carlo method was used to calculate correction factors for the diode used in the percentage depth-dose (PDD) measurement, and to study the effect of the detector size on the beam profile. The PDD data were determined using the diode measurement corrected for its energy and angular response. It was found that the PTW N23342 and Markus parallel-plate chamber can be used directly to measure the PDD for this beam quality with 2% uncertainty. Measurement and Monte Carlo results have shown that the detector size has a significant effect on the penumbral profile. Film and diode detectors have a better spatial resolution compared to ionization chambers, but they may give an incorrect profile tail due to either nonlinear response at low energy or angular dependence. This can be corrected using the ionization-chamber measurement, based on the Monte Carlo analysis. The isodose distributions for this x-ray unit are presented.

Brachytherapy↗

PubMed1998

Study of dosimetry consistency for kilovoltage x-ray beams.

In this paper, the consistency of kilovoltage (tube potentials between 40 and 300 kV) x-ray beam dosimetry using the "in-air" method and the in-phantom measurement has been studied. The procedures for the measurement of the central-axis depth-dose curve, which serve as a link between the dose at the reference depth to the dose elsewhere in a phantom, were examined. The uncertainties on the measured dose distributions were analyzed with the emphasis on the surface dose measurement. The Monte Carlo method was used to calculate the perturbation correction factors for a photon diode and a NACP plane-parallel ionization chamber at different depths in a water phantom irradiated by 100-300 kV (2.43 mm Al-3.67 mm Cu half-value layer) x-ray beams. The depth-dose curves measured with these two detectors, after correcting for the perturbation effect (up to 15% corrections), agreed with each other to within 1.5%. Comparisons of the doses at the phantom surface and at 2 cm depth in water for photon beams of 100-300 kV tube potential obtained using the "backscatter" method and those using the "in-phantom" measurement have shown that the "in-air" method can be equally applied to this energy range if the depth-dose curve can be measured accurately. To this end, measured depth ionization curves require depth-dependent correction factors.

Humans↗