Evaluation of Gafchromic EBT2 film for the measurement of anisotropy function for high-dose-rate 192Ir brachytherapy source with respect to thermoluminescent dosimetry

Aim

The aim of this work was to assess the suitability of the use of a Gafchromic EBT2 film for the measurement of anisotropy function for microSelectron HDR ¹⁹²Ir (classic) source with a comparative dosimetry method using a Gafchromic EBT2 film and thermoluminescence dosimeters (TLDs).

Background

Sealed linear radiation sources are commonly used for high dose rate (HDR) brachytherapy treatments. Due to self-absorption and oblique filtration of radiation in the source capsule material, an inherent anisotropy is present in the dose distribution around the source which can be described by a measurable two-dimensional anisotropy function, F(r, θ).

Materials and methods

Measurements were carried out in a specially designed and locally fabricated PMMA phantom with provisions to accommodate miniature LiF TLD rods and EBT2 film dosimeters at identical radial distances with respect to the ¹⁹²Ir source.

Results

The data of anisotropy function generated by the use of the Gafchromic EBT2 film method are in agreement with their TLD measured values within 4%. The produced data are also consistent with their experimental and Monte Carlo calculated results for this source available in the literature.

Conclusion

Gafchromic EBT2 film was found to be a feasible dosimeter in determining anisotropy in the dose distribution of ¹⁹²Ir source. It offers high resolution and is a viable alternative to TLD dosimetry at discrete points. The method described in this paper is useful for comparing the performances of detectors and can be applied for other brachytherapy sources as well.     

A Monte Carlo study on dose distribution validation of GZP6 60Co stepping source

Aim

Stepping source in brachytherapy systems is used to treat a target lesion longer than the effective treatment length of the source. Cancerous lesions in the cervix, esophagus and rectum are examples of such a target lesion.

Background

In this study, the stepping source of a GZP6 afterloading intracavitary brachytherapy unit was simulated using Monte Carlo (MC) simulation and the results were used for the validation of the GZP6 treatment planning system (TPS).

Materials and methods

The stepping source was simulated using MCNPX Monte Carlo code. Dose distributions in the longitudinal plane were obtained by using a matrix shift method for esophageal tumor lengths of 8 and 10 cm. A mesh tally has been employed for the absorbed dose calculation in a cylindrical water phantom. A total of 5 × 10⁸ photon histories were scored and the MC statistical error obtained was at the range of 0.008–3.5%, an average of 0.2%.

Results

The acquired MC and TPS isodose curves were compared and it was shown that the dose distributions in the longitudinal plane were relatively coincidental. In the transverse direction, a maximum dose difference of 7% and 5% was observed for tumor lengths of 8 and 10 cm, respectively.

Conclusion

Considering that the certified source activity is given with ±10% uncertainty, the obtained difference is reasonable. It can be concluded that the accuracy of the dose distributions produced by GZP6 TPS for the stepping source is acceptable for its clinical applications.     

Using matrix summation method for three dimensional dose calculation in brachytherapy

AimThe purpose of this study is to calculate radiation dose around a brachytherapy source in a water phantom for different seed locations or rotation the sources by the matrix summation method.BackgroundMonte Carlo based codes like MCNP are widely used for performing radiation transport calculations and dose evaluation in brachytherapy. But for complicated situations, like using more than one source, moving or rotating the source, the routine Monte Carlo method for dose calculation needs a long time running.Materials and methodsThe MCNPX code has been used to calculate radiation dose around a ¹⁹²Ir brachytherapy source and saved in a 3D matrix. Then, we used this matrix to evaluate the absorbed dose in any point due to some sources or a source which shifted or rotated in some places by the matrix summation method.ResultsThree dimensional (3D) dose results and isodose curves were presented for ¹⁹²Ir source in a water cube phantom shifted for 10 steps and rotated for 45 and 90° based on the matrix summation method. Also, we applied this method for some arrays of sources.ConclusionThe matrix summation method can be used for 3D dose calculations for any brachytherapy source which has moved or rotated. This simple method is very fast compared to routine Monte Carlo based methods. In addition, it can be applied for dose optimization study.     

A Monte Carlo evaluation for effects of probable dimensional uncertainties of low dose rate brachytherapy seeds on dose

The aim of this study is to determine effects of size deviations of brachytherapy seeds on two dimensional dose distributions around the seed. Although many uncertainties are well known, the uncertainties which stem from geometric features of radiation sources are weakly considered and predicted. Neither TG-43 report which is not completely in common consensus, nor individual scientific MC and experimental studies include sufficient data for geometric uncertainties. Sizes of seed and its components can vary in a manufacturing deviation. This causes geometrical uncertainties, too. In this study, three seeds which have different geometrical properties were modeled using EGSnrc-Code Packages. Seeds were designed with all their details using the geometry package. 5% deviations of seed sizes were assumed. Modified seeds were derived from original seed by changing sizes by 5%. Normalizations of doses which were calculated from three kinds of brachytherapy seed and their derivations were found to be about 3%–20%. It was shown that manufacturing differences of brachytherapy seed cause considerable changes in dose distribution.     

Calibration of 192Ir high dose rate brachytherapy source using different calibration procedures

Aim

To calibrate Ir-192 high dose rate (HDR) brachytherapy source using different calibration methods and to determine the accuracy and suitability of each method for routine calibrations.

Background

The source calibration is an essential part of the quality assurance programme for dosimetry of brachytherapy sources. The clinical use of brachytherapy source requires an independent measurement of the air kerma strength according to the recommendations of medical physics societies.

Materials and methods

The Ir-192 HDR brachytherapy source from Gammamed plus machine (Varian Medical Systems, Palo Alto, CA) was calibrated using three different procedures, one using the well-type ionization chamber, second by the in-air calibration method and third using solid water phantoms. The reference air kerma rate (RAKR) of the source was determined using Deutsche Gesellschaft fur Medizinische Physik (DGMP) recommendations.

Results

The RAKR determined using different calibration methods are in good agreement with the manufacturer stated value. The mean percentage variations of 0.21, −0.94, −0.62 and 0.58 in RAKR values with respect to the manufacturer quoted values were observed with the well-type chamber, in-air calibration, cylindrical phantom and slab phantom measurements, respectively.

Conclusion

Measurements with a well-type chamber are relatively simple to perform. For in-air measurements, the indigenously designed calibration jig provides an accurate positioning of the source and chamber with minimum scatter contribution. The slab phantom system has an advantage that no additional phantom and chamber are required other than those used for external beam therapy dosimetry. All the methods of calibration discussed in this study are effective to be used for routine calibration purposes.     

Tissue composition effect on dose distribution in neutron brachytherapy/neutron capture therapy

Aim

The aim of this study is to assess the effect of the compositions of various soft tissues and tissue-equivalent materials on dose distribution in neutron brachytherapy/neutron capture therapy.

Background

Neutron brachytherapy and neutron capture therapy are two common radiotherapy modalities.

Materials and methods

Dose distributions were calculated around a low dose rate ²⁵²Cf source located in a spherical phantom with radius of 20.0 cm using the MCNPX code for seven soft tissues and three tissue-equivalent materials. Relative total dose rate, relative neutron dose rate, total dose rate, and neutron dose rate were calculated for each material. These values were determined at various radial distances ranging from 0.3 to 15.0 cm from the source.

Results

Among the soft tissues and tissue-equivalent materials studied, adipose tissue and plexiglass demonstrated the greatest differences for total dose rate compared to 9-component soft tissue. The difference in dose rate with respect to 9-component soft tissue varied with compositions of the materials and the radial distance from the source. Furthermore, the total dose rate in water was different from that in 9-component soft tissue.

Conclusion

Taking the same composition for various soft tissues and tissue-equivalent media can lead to error in treatment planning in neutron brachytherapy/neutron capture therapy. Since the International Commission on Radiation Units and Measurements (ICRU) recommends that the total dosimetric uncertainty in dose delivery in radiotherapy should be within ±5%, the compositions of various soft tissues and tissue-equivalent materials should be considered in dose calculation and treatment planning in neutron brachytherapy/neutron capture therapy.     

Monte Carlo evaluation of the dose calculation algorithm of TomoTherapy for clinical cases in dynamic jaws mode

PurposeFor the TomoTherapy^® system, longitudinal conformation can be improved by selecting a smaller field width but at the expense of longer treatment time. Recently, the TomoEdge^® feature has been released with the possibility to move dynamically the jaws at the edges of the target volume, improving longitudinal penumbra and enabling faster treatments. Such delivery scheme requires additional modeling of treatment delivery. Using a previously validated Monte Carlo model (TomoPen), we evaluated the accuracy of the implementation of TomoEdge in the new dose engine of TomoTherapy for 15 clinical cases.MethodsTomoPen is based on PENELOPE. Particle tracking in the treatment head is performed almost instantaneously by 1) reading a particle from a phase-space file corresponding to the largest field and 2) correcting the weight of the particle depending on the actual jaw and MLC configurations using Monte Carlo pre-generated data. 15 clinical plans (5 head-and-neck, 5 lung and 5 prostate tumors) planned with TomoEdge and with the last release of the treatment planning system (VoLO^®) were re-computed with TomoPen. The resulting dose-volume histograms were compared.ResultsGood agreement was achieved overall, with deviations for the target volumes typically within 2% (D₉₅), excepted for small lung tumors (17 cm³) where a maximum deviation of 4.4% was observed for D₉₅. The results were consistent with previously reported values for static field widths.ConclusionsFor the clinical cases considered in the present study, the introduction of TomoEdge did not impact significantly the accuracy of the computed dose distributions.     

Monte Carlo dose calculation of GZP6 60Co stepping source based on a matrix shift technique

Background

As a routine method for stepping source simulation, a Monte Carlo program is run according to the number of steps and then the summation of dose from each run is taken to obtain total dose distribution. This method is time consuming.

Aim

As an alternative method, a matrix shift based technique was applied to simulate a stepping source for brachytherapy.

Materials and methods

The stepping source of GZP6 brachytherapy unit was simulated. In a matrix shift method, it is assumed that a radiation source is stationary and instead the data matrix is shifted based on the number of steps. In this study, by running MCNPX program for one point and calculation of the dose matrix using the matrix shift method, the isodose curves for the esophageal cancer tumor lengths of 4 and 6 cm were obtained and compared with the isodose curves obtained by running MCNPX programs in each step position separately (15 and 23 steps for esophageal cancer tumor lengths of 4 and 6 cm, respectively).

Results

The difference between the two dose matrixes for the stepping and matrix shift methods based on the average dose differences are 3.85 × 10⁻⁴ Gy and 5.19 × 10⁻⁴ Gy for treatment length of 4 cm and 6 cm, respectively. Dose differences are insignificant and these two methods are equally valid.

Conclusions

The matrix shift method presented in this study can be used for calculation of dose distribution for a brachytherapy stepping source as a quicker tool compared to other routine Monte Carlo based methods.     

Air kerma strength characterization of a GZP6 Cobalt-60 brachytherapy source

Background

Task group number 40 (TG-40) of the American Association of Physicists in Medicine (AAPM) has recommended calibration of any brachytherapy source before its clinical use. GZP6 afterloading brachytherapy unit is a ⁶⁰Co high dose rate (HDR) system recently being used in some of the Iranian radiotherapy centers.

Aim

In this study air kerma strength (AKS) of ⁶⁰Co source number three of this unit was estimated by Monte Carlo simulation and in air measurements.

Materials and methods

Simulation was performed by employing the MCNP-4C Monte Carlo code. Self-absorption of the source core and its capsule were taken into account when calculating air kerma strength. In-air measurements were performed according to the multiple distance method; where a specially designed jig and a 0.6 cm³ Farmer type ionization chamber were used for the measurements. Monte Carlo simulation, in air measurement and GZP6 treatment planning results were compared for primary air kerma strength (as for November 8th 2005).

Results

Monte Carlo calculated and in air measured air kerma strength were respectively equal to 17240.01 μGym² h⁻¹ and 16991.83 μGym² h⁻¹. The value provided by the GZP6 treatment planning system (TPS) was “15355 μGym² h⁻¹”.

Conclusion

The calculated and measured AKS values are in good agreement. Calculated-TPS and measured-TPS AKS values are also in agreement within the uncertainties related to our calculation, measurements and those certified by the GZP6 manufacturer. Considering the uncertainties, the TPS value for AKS is validated by our calculations and measurements, however, it is incorporated with a large uncertainty.     

The structure of polymer chains in confinement. A Monte Carlo study

A coarse-grained model of polymer star chains confined in two parallel impenetrable surfaces, which were attractive for polymer beads was studied. The flexible homopolymer chains were built of united atoms whose positions in space were restricted to vertices of a simple cubic lattice. The chains were modeled in good solvent conditions and, thus, there were no long-range specific interactions between polymer beads—only the excluded volume was present. The influence of the polymer density and the distances between the confining surfaces on the properties of star-branched polymers was studied. It is shown that the chains adsorbed on one surface could change their position so that they swap between both surfaces with frequency depending on the size of the slit and on the density of the system only. The increase of the polymer density diminished the frequency of jumps and caused that chains became only partially adsorbed. The analysis of structural elements of chains showed that the increase of the density of the system leads to increase of the number of bridges connecting the two adsorbing surfaces, thus, the frequency of jumps between them decreases.     

A depth dose study between AAA and AXB algorithm against Monte Carlo simulation using AIP CT of a 4D dataset from a moving phantom

Aim

To identifying depth dose differences between the two versions of the algorithms using AIP CT of a 4D dataset.

Monte Carlo docking simulations of cyclomaltoheptaose and dimethyl cyclomaltoheptaose with paclitaxel

The molecular basis for the remarkable enhancement of the solubility of paclitaxel by O-dimethylcyclomaltoheptaose (DM-beta-CD) over cyclomaltoheptaose (beta-cyclodextrin, beta-CD) was investigated with Monte Carlo docking-minimization simulation. As possible guests of inclusion complexation for the host cyclic oligosaccharides, two functional moieties of the suggested solution structure of paclitaxel were used where one is the C-3'N benzoyl moiety (B-ring) and the other is a hydrophobic (HP) cluster site among the C-3' phenyl (C-ring), C-2 benzoate (A-ring), and C-4 acetoxy moieties. The energetic preference of inclusion complexation of DM-beta-CD over beta-CD was analyzed on the basis of more efficient partitioning process of DM-beta-CD into the hydrophobic cluster site of the paclitaxel.     

Mixtures of Associating and Non-associating Chains on Activated Surfaces: A Monte Carlo Approach

Abstract

The behavior of mixtures of associating and non-associating chains confined in pores with activated surfaces is studied by means of molecular simulation. The fluid molecules are modeled as a chain of four tangent Lennard-Jones spheres. Some of the chains have an additional associating square-well site placed in an end sphere. The activated surfaces of the slit pore are modeled via an integrated Lennard-Jones (10-4-3) potential with specific association sites protruding from the surface. We present Gibbs ensemble Monte Carlo simulation results for the partitioning of mixtures of chains in the bulk and confined phases for this particular model. The chain-wall association governs the adsorption behavior of the system. The preferential adsorption of associating chains is seen to strongly depend on temperature and pore width. Selectivities obtained are in the range of those seen in experiments of alkane-alkanol mixtures.     

A Monte Carlo method for Bayesian analysis of linkage between single markers and quantitative trait loci. II. A simulation study

A Bayesian approach to the statistical mapping of Quantitative Trait Loci (QTLs) using single markers was implemented via Markov Chain Monte Carlo (MCMC) algorithms for parameter estimation and hypothesis testing. Parameters were estimated by marginal posterior means computed with a Gibbs sampler with data augmentation. Variables sampled included the augmented data (marker-QTL genotypes, polygenic effects), the event of linkage or nonlinkage, and the parameters (allele frequencies, QTL substitution effect, recombination rate, polygenic and residual variances). The analysis was evaluated empirically via application to simulated granddaughter designs consisting of 2000 sons, 20 related sires and their ancestors. Results obtained in this study and preliminary work on multiple linked markers and multiple QTLs support the usefulness of the Bayesian method for the statistical mapping of QTLs.     

A P-NMR study of structural and functional aspects of phosphate and phosphonate distribution in Tetrahymena

³¹P-NMR has been used to study the chemical nature of cytoplasmic components of live Tetrahymena in a non-invasive manner. The technique has further been used to characterize the physical behaviour of lipids extracted from this organism. In particular, we have shown the presence of large quantities of pyrophosphate and of tripolyphosphate in acid extracts of the organism. These are not detectable in the live cell due to the motionally rigid nature of the storage granules. We have characterized the distribution of phosphonic acids in the organism and followed the phase behaviour of the extracted cell lipids. Aqueous dispersions of extracted lipid show both bilayer and non-bilayer behaviour in the range of the growth temperature. The phosphonolipid in Tetrahymena appears to play a role similar to that of phosphatidylethanolamine in regulating the phase behaviour of the membrane. The high degree of unsaturation in the fatty acids of Tetrahymena is most likely responsible for the polymorphic phase behaviour observed near the growth temperature.     

Dosimetric impact of statistical uncertainty on Monte Carlo dose calculation algorithm in volumetric modulated arc therapy using Monaco TPS for three different clinical cases

AimTo study the dosimetric impact of statistical uncertainty (SU) per plan on Monte Carlo (MC) calculation in Monaco? treatment planning system (TPS) during volumetric modulated arc therapy (VMAT) for three different clinical cases.BackgroundDuring MC calculation SU is an important factor to decide dose calculation accuracy and calculation time. It is necessary to evaluate optimal acceptance of SU for quality plan with reduced calculation time.Materials and methodsThree different clinical cases as the lung, larynx, and prostate treated using VMAT technique were chosen. Plans were generated with Monaco? V5.11 TPS with 2% statistical uncertainty. By keeping all other parameters constant, plans were recalculated by varying SU, 0.5%, 1%, 2%, 3%, 4%, and 5%. For plan evaluation, conformity index (CI), homogeneity index (HI), dose coverage to PTV, organ at risk (OAR) dose, normal tissue receiving dose ≥5 Gy and ≥10 Gy, integral dose (NTID), calculation time, gamma pass rate, calculation reproducibility and energy dependency were analyzed.ResultsCI and HI improve as SU increases from 0.5% to 5%. No significant dose difference was observed in dose coverage to PTV, OAR doses, normal tissue receiving dose ≥5 Gy and ≥10 Gy and NTID. Increase of SU showed decrease in calculation time, gamma pass rate and increase in PTV max dose. No dose difference was seen in calculation reproducibility and dependent on energy.ConclusionFor VMAT plans, SU can be accepted from 1% to 3% per plan with reduced calculation time without compromising plan quality and deliverability by accepting variations in point dose within the target.     

Analyzing indirect secondary electron contrast of unstained bacteriophage T4 based on SEM images and Monte Carlo simulations

The indirect secondary electron contrast (ISEC) condition of the scanning electron microscopy (SEM) produces high contrast detection with minimal damage of unstained biological samples mounted under a thin carbon film. The high contrast image is created by a secondary electron signal produced under the carbon film by a low acceleration voltage. Here, we show that ISEC condition is clearly able to detect unstained bacteriophage T4 under a thin carbon film (10-15 nm) by using high-resolution field emission (FE) SEM. The results show that FE-SEM provides higher resolution than thermionic emission SEM. Furthermore, we investigated the scattered electron area within the carbon film under ISEC conditions using Monte Carlo simulation. The simulations indicated that the image resolution difference is related to the scattering width in the carbon film and the electron beam spot size. Using ISEC conditions on unstained virus samples would produce low electronic damage, because the electron beam does not directly irradiate the sample. In addition to the routine analysis, this method can be utilized for structural analysis of various biological samples like viruses, bacteria, and protein complexes.     

Effects of macro elements and nitrogen source on adventitious root growth and ginsenoside production in ginseng (Panax ginseng C. A. Meyer)

We investigated the effects on ginseng adventitious root growth and ginsenoside production when macro-element concentrations and nitrogen source were manipulated in the culture media. Biomass growth was greatest in the medium supplemented with 0.5-strength NH₄PO₃, whereas ginsenoside accumulation was highest (9.90 mg g^-1 DW) in the absence of NH₄PO₃. At levels of 1.0-strength KNO₃, root growth was maximum, but a 2.0 strength of KNO₃ led to the greatest ginsenoside content (9.85 mg g^-l). High concentrations of MgSO₄ were most favorable for both root growth and ginsenoside accumulation (up to 8.89 mg g^-1 DW). Root growth and ginsenoside content also increased in proportion to the concentration of CaCI₂ in the medium, with the greatest accumulation of ginsenoside (8.91 mg g^-1 DW) occurring at a 2.0 strength. The NH₄/NO₃ ^-- ratio also influenced adventitious root growth and ginsenoside production; both parameters were greater when the NO₃ ^- concentration was higher than that of NH₄ ⁺. Maximum root growth was achieved at an NH₄ ⁺/NO₃ ^- ratio of 7.19/18.50, while ginsenoside production was greatest (83.37 mg L^-1) when NO₃ ^- was used as the sole N source.