Impact of photon cross section uncertainties on Monte Carlo-determined depth-dose distributions.

This work studies the impact of systematic uncertainties associated to interaction cross sections on depth dose curves determined by Monte Carlo simulations. The corresponding sensitivity factors are quantified by changing cross sections by a given amount and determining the variation in the dose. The influence of total and partial photon cross sections is addressed. Partial cross sections for Compton and Rayleigh scattering, photo-electric effect, and pair production have been accounted for. The PENELOPE code was used in all simulations. It was found that photon cross section sensitivity factors depend on depth. In addition, they are positive and negative for depths below and above an equilibrium depth, respectively. At this depth, sensitivity factors are null. The equilibrium depths found in this work agree very well with the mean free path of the corresponding incident photon energy. Using the sensitivity factors reported here, it is possible to estimate the impact of photon cross section uncertainties on the uncertainty of Monte Carlo-determined depth dose curves.

The Influence of DNA Configuration on the Direct Strand Break Yield.

PURPOSE: To study the influence of DNA configuration on the direct damage yield. No indirect effect has been accounted for. METHODS: The GEANT4-DNA code was used to simulate the interactions of protons and alpha particles with geometrical models of the A-, B-, and Z-DNA configurations. The direct total, single, and double strand break yields and site-hit probabilities were determined. Certain features of the energy deposition process were also studied. RESULTS: A slight increase of the site-hit probability as a function of the incident particle linear energy transfer was found for each DNA configuration. Each DNA form presents a well-defined site-hit probability, independently of the particle linear energy transfer. Approximately 70% of the inelastic collisions and ~60% of the absorbed dose are due to secondary electrons. These fractions are slightly higher for protons than for alpha particles at the same incident energy. CONCLUSIONS: The total direct strand break yield for a given DNA form depends weakly on DNA conformation topology. This yield is practically determined by the target volume of the DNA configuration. However, the double strand break yield increases with the packing ratio of the DNA double helix; thus, it depends on the DNA conformation.

Experimental and Monte Carlo-simulated spectra of standard mammography-quality beams.

UNLABELLED: A spectrometric study of standard mammography-quality beams by using experimental and Monte Carlo simulation methods was carried out in this work. The qualities of these beams are described according to the International Electrotechical Commission 61267 standard and the Technical Report Series 457 International Atomic Energy Agency report. Specifically, the non-attenuated RQR-M beam series was studied. METHODS: A Si-PIN diode-based spectrometer and the PENELOPE Monte Carlo code (v. 2008F1) were used for experiments and simulations, respectively. In addition, an ionization chamber was used to determine the half-value layers (HVLs) of each beam quality. The measurements were done in the mammography dosimeter calibration setup of our laboratory, and the Monte Carlo simulations reproduced such conditions. RESULTS: The relative differences between the HVLs calculated from experimental and simulated spectra were lower than 2.4% for all the beam qualities studied. These differences are 1.2% and 3.1% when comparing the HVLs calculated from the experimental and simulated spectra to those determined by using the ionization chamber, respectively. A semi-empirical relation was found to obtain the nominal tube potential from the effective tube potential. CONCLUSION: According to our results, the mammography beams used in this work have energy spectra similar to clinical beams.

Estimation of the RBE of mammography-quality beams using a combination of a Monte Carlo code with a B-DNA geometrical model.

The PENELOPE code is used to determine direct strand break yields corresponding to photons from a (60)Co source and 28 and 30 kV x-ray beams impacting on a B-DNA geometrical model, which accounts for five organizational levels of the human genetic material. Direct single, double and total strand break probabilities are determined in a liquid water homogeneous medium with 1.06 g cm(-3) density. The spectra produced by the x-ray beams at various depths in the phantom have been used to study the dependence of the damage yield on the depth. The relative biological effectiveness (RBE) is also estimated using the (60)Co radiation qualities as the reference. According to this work, the damage probabilities and thus the RBE are, within the uncertainties, similar for both x-ray energies and are independent of the depth into the phantom. Furthermore, the total strand break yield is invariant with respect to the energy of the incident photons. The RBE for low-energy x-ray beams determined here (1.3 ± 0.1) is lower than that reported by Kellerer, taking into account that he used a 200 kV radiation as the reference quality. However, our RBE values are consistent with those determined by Kühne et al (2005 Radiat. Res. 164 669-76), which used the same biological endpoint and reference quality as our study. Also, our RBE values are similar to those determined by Verhaegen and Reniers (2004 Radiat. Res. 162 592-9).

The invariance of the total direct DNA strand break yield.

PURPOSE: The invariance of the total direct strand break yield when DNA is irradiated by different types of particles and energies has been reported by previous works. This study is intended to explain the physical causes of this behavior. METHODS: The GEANT4-DNA extension of the GEANT4 general purpose Monte Carlo simulation toolkit has been used to determine direct strand break yields induced by protons and alpha particles impacting on a B-DNA geometrical model, including five organization levels of the human genetic material. The linear energy transfer (LET) of such particles ranges from 4.8 keV/microm (10 MeV protons) to about 235 keV/microm (2 MeV alpha particles), at 5.225 pm depth (near the center of the region of interest). Direct total, single and double strand break probabilities have been determined in a liquid water homogeneous medium with a 1.06 g/cm3 density. The energetic spectra of single strand breaks (SSB), the number of energy deposition events, and the SSB/event ratio were determined. RESULTS: The target-hit probability was found to be independent of both the type and the energy of the incident particle, even if this latter is a secondary electron. This probability is determined by the geometrical properties of the system. The total strand break yield and the number of energy deposition events required to reach a certain absorbed dose were found nearly independent of the type and energy of the incident ion (proton or alpha). In contrast, the double strand break (DSB) yield was found strongly dependent on the LET of the incident radiation. CONCLUSIONS: The SSB generation process is homogeneous and independent of the LET of the particles involved, at least within the proton and alpha particle energy range here studied. The target-hit probability is only determined by the ratio between the total volume occupied by targets and that of the ROI where the radiation deposits its energy. The maximum separation distance between two adjacent SSBs to produce a DSB is the parameter that breaks the homogeneity of the target-hit process, making the DSB production process strongly heterogeneous.

Evaluation of radiographic image quality parameters obtained with the REX simulator.

The PTW REX phantom was used to study the radiographic image quality parameters in X-ray devices in the X-ray Diagnostics Department, as well as the system of film processing at the University Hospital of Rio de Janeiro State. X-ray devices were evaluated by performing tests on 11 screen-film combinations from X-ray devices in 3 rooms. The results showed that six film-screen combinations exhibited poor performances. For determination of air kerma output in the X-ray field, two devices presented significant variation >2 %. The grid attenuation factor in three devices had been approved, while two films were within the limits of sensitometric specifications. The modulation transfer function, which evaluates the level of image degradation, revealed that five film-screen combinations exhibited bad performance. The tests with the REX phantom revealed that the X-ray equipment and the system of processing at the University Hospital presented discrepancies in relation to the expected values, contributing to loss of quality of the radiographs.

Experimental derivation of the fluence non-uniformity correction for air kerma near brachytherapy linear sources.

In brachytherapy, one of the elements to take into account for measurements free in air is the non-uniformity of the photon fluence due to the beam divergence that causes a steep dose gradient near the source. The correction factors for this phenomenon have been usually evaluated by two available theories by Kondo and Randolph [Radiat. Res. 13, 37-60 (1960)] and Bielajew [Phys. Med. Biol. 35, 517-538 (1990)], both conceived for point sources. This work presents the experimental validation of the Monte Carlo calculations made by Rodriguez and deAlmeida [Phys. Med. Biol. 49, 1705-1709 (2004)] for the non-uniformity correction specifically for a Cs-137 linear source measured using a Farmer type ionization chamber. The experimental values agree very well with the Monte Carlo calculations and differ from the results predicted by both theoretical models widely used. This result confirms that for linear sources there are some important differences at short distances from the source and emphasizes that those theories should not be used for linear sources. The data provided in this study confirm the limitations of the mentioned theories when linear sources are used. Considering the difficulties and uncertainties associated with the experimental measurements, it is recommended to use the Monte Carlo data to assess the non-uniformity factors for linear sources in situations that require this knowledge.

The measurement of photoneutron dose in the vicinity of clinical linear accelerators.

In Brazil, the replacement of rather old cobalt and cesium teletherapy machines with high-energy (E > 10 MV) medical linear accelerators (linacs) started in the year 2000, as part of an effort by the Ministry of Health to update radiotherapy installations. Since then, the contamination of undesirable neutrons in the therapeutic beam generated by these high-energy photons has become an issue of concern when considering patient and occupational doses. The walls of the treatment room are shielded to attenuate the primary and secondary X-ray fluence, and this shielding is generally considered adequate also to attenuate neutrons. However, these neutrons are scattered through the treatment room maze and might result in a radiological problem at the door entrance, an area of high occupancy by the workers of a radiotherapy facility. This paper presents and discusses the results of ambient dose equivalent measurements of neutron using bubble detectors. The measurements were made at different points inside the treatment rooms, including the isocentre and the maze. Several radiation oncology centres, which are users of Varian Clinac or Siemens machines, have agreed to allow measurements to be taken at their facilities. The measured values were compared with the results obtained through the semi-empirical Kersey method of neutron dose equivalent calculation at maze entrances, with reported values provided by the manufacturers as well as values published in the literature. It was found that the measured values were below the dose limits adopted by the Brazilian Regulatory Agency (CNEN), requiring no additional shielding in any of the points measured.

Determination of absorbed dose in water at the reference point d(r0, theta0) for an 192Ir HDR brachytherapy source using a Fricke system.

A ring-shaped Fricke device was developed to measure the absolute dose on the transverse bisector of a 192Ir high dose rate (HDR) source at 1 cm from its center in water, D(r0, theta0). It consists of a polymethylmethacrylate (PMMA) rod (axial axis) with a cylindrical cavity at its center to insert the 192Ir radioactive source. A ring cavity around the source with 1.5 mm thickness and 5 mm height is centered at 1 cm from the central axis of the source. This ring cavity is etched in a disk shaped base with 2.65 cm diameter and 0.90 cm thickness. The cavity has a wall around it 0.25 cm thick. This ring is filled with Fricke solution, sealed, and the whole assembly is immersed in water during irradiations. The device takes advantage of the cylindrical geometry to measure D(r0, theta0). Irradiations were performed with a Nucletron microselectron HDR unit loaded with an 192Ir Alpha Omega radioactive source. A Spectronic 1001 spectrophotometer was used to measure the optical absorbance using a 1 mL quartz cuvette with 1.00 cm light pathlength. The PENELOPE Monte Carlo code (MC) was utilized to simulate the Fricke device and the 192Ir Alpha Omega source in detail to calculate the perturbation introduced by the PMMA material. A NIST traceable calibrated well type ionization chamber was used to determine the air-kerma strength, and a published dose-rate constant was used to determine the dose rate at the reference point. The time to deliver 30.00 Gy to the reference point was calculated. This absorbed dose was then compared to the absorbed dose measured by the Fricke solution. Based on MC simulation, the PMMA of the Fricke device increases the D(r0, theta0) by 2.0%. Applying the corresponding correction factor, the D(r0, theta0) value assessed with the Fricke device agrees within 2.0% with the expected value with a total combined uncertainty of 3.43% (k=1). The Fricke device provides a promising method towards calibration of brachytherapy radiation sources in terms of D(r0, theta0) and audit HDR source calibrations.

The radiation field characteristics of a 137Cs source used for calibration of radiation protection instruments.

In this work, an experimental study of the radiation beam characteristics of the photon beams produced by commercial irradiator that uses a single collimated 137Cs source is performed. A set of lead attenuators are placed at the exit window of the irradiator to vary the air Kerma that is required to cover the instrument scales at a particular calibration distance. The results strongly enhance the findings of previous studies that a single source with appropriate attenuation and collimation is a practical solution for calibration laboratories with an acceptable overall accuracy.

The photon fluence non-uniformity correction for air kerma near Cs-137 brachytherapy sources.

The use of brachytherapy sources in radiation oncology requires their proper calibration to guarantee the correctness of the dose delivered to the treatment volume of a patient. One of the elements to take into account in the dose calculation formalism is the non-uniformity of the photon fluence due to the beam divergence that causes a steep dose gradient near the source. The correction factors for this phenomenon have been usually evaluated by the two theories available, both of which were conceived only for point sources. This work presents the Monte Carlo assessment of the non-uniformity correction factors for a Cs-137 linear source and a Farmer-type ionization chamber. The results have clearly demonstrated that for linear sources there are some important differences among the values obtained from different calculation models, especially at short distances from the source. The use of experimental values for each specific source geometry is recommended in order to assess the non-uniformity factors for linear sources in clinical situations that require special dose calculations or when the correctness of treatment planning software is verified during the acceptance tests.

The design of a new insert for calibration of LDR and HDR 192Ir sources in a well-type ionization chamber.

Some well-type ionization chambers, present a very small sweet spot that are sufficient for small HDR sources. However, if a longer HDR source or LDR wires are calibrated, the positional uncertainty increases and an approximated correction factor must be applied, resulting in an increased uncertainty. One of the ways to avoid this problem would be to flatten the well chamber response by increasing its sweet spot region. This work uses the Monte Carlo code PENELOPE to simulate the response of a well-type chamber HDR-1000, with its original insert, by using an HDR 192Ir source and proposes a new insert design that increases its flatness region from 1.0 cm to approximately 4.0 cm (+/- 2.0 cm about the peak response).