Monte Carlo Dosimetry of the 60Co BEBIG High Dose Rate for Brachytherapy

Introduction

The use of high-dose-rate brachytherapy is currently a widespread practice worldwide. The most common isotope source is ¹⁹²Ir, but ⁶⁰Co is also becoming available for HDR. One of main advantages of ⁶⁰Co compared to ¹⁹²Ir is the economic and practical benefit because of its longer half-live, which is 5.27 years. Recently, Eckert & Ziegler BEBIG, Germany, introduced a new afterloading brachytherapy machine (MultiSource^®); it has the option to use either the ⁶⁰Co or ¹⁹²Ir HDR source. The source for the Monte Carlo calculations is the new ⁶⁰Co source (model Co0.A86), which is referred to as the new BEBIG ⁶⁰Co HDR source and is a modified version of the ⁶⁰Co source (model GK60M21), which is also from BEBIG.

Objective and Methods

The purpose of this work is to obtain the dosimetry parameters in accordance with the AAPM TG-43U1 formalism with Monte Carlo calculations regarding the BEBIG ⁶⁰Co high-dose-rate brachytherapy to investigate the required treatment-planning parameters. The geometric design and material details of the source was provided by the manufacturer and was used to define the Monte Carlo geometry. To validate the source geometry, a few dosimetry parameters had to be calculated according to the AAPM TG-43U1 formalism. The dosimetry studies included the calculation of the air kerma strength S _k, collision kerma in water along the transverse axis with an unbounded phantom, dose rate constant and radial dose function. The Monte Carlo code system that was used was EGSnrc with a new cavity code, which is a part of EGS++ that allows calculating the radial dose function around the source. The spectrum to simulate ⁶⁰Co was composed of two photon energies, 1.17 and 1.33 MeV. Only the gamma part of the spectrum was used; the contribution of the electrons to the dose is negligible because of the full absorption by the stainless-steel wall around the metallic ⁶⁰Co. The XCOM photon cross-section library was used in subsequent simulations, and the photoelectric effect, pair production, Rayleigh scattering and bound Compton scattering were included in the simulation. Variance reduction techniques were used to speed up the calculation and to considerably reduce the computer time. The cut-off energy was 10 keV for electrons and photons. To obtain the dose rate distributions of the source in an unbounded liquid water phantom, the source was immersed at the center of a cube phantom of 100 cm³. The liquid water density was 0.998 g/cm³, and photon histories of up to 10¹⁰ were used to obtain the results with a standard deviation of less than 0.5% (k = 1). The obtained dose rate constant for the BEBIG ⁶⁰Co source was 1.108±0.001 cGyh^-1U^-1, which is consistent with the values in the literature. The radial dose functions were compared with the values of the consensus data set in the literature, and they are consistent with the published data for this energy range.     

Plasmon Spectroscopy for Subnanometric Copper Particles: Dielectric Function and Core–Shell Sizing

In the last years, there has been a growing interest in the study of transition metal nanoparticles (Nps) due to their potential applications in several fields of science and technology. In particular, their optical properties are governed by the characteristics of the dielectric function of the metal, its size and environment. This work analyses the separated contribution of free and bound electrons on the optical properties of copper Nps. Usually, the contribution of free electrons to the dielectric function is corrected for particle size through the modification of the damping constant, which is changed as usual introducing a term inversely proportional to the particle’s radius to account for the extra collisions with the boundary when the size approaches the electronic mean free path limit (about 10 nm). For bound electron contribution, the interband transitions from the d-band to the conduction band are considered together with the fact that the electronic density of states in the conduction band must be made size-dependent to account for the larger spacing between electronic energy levels as the particle decreases in size below 2 nm. Taking into account these specific modifications of free and bound electron contributions to the dielectric function, it was possible to fit the bulk complex dielectric function, and consequently, determine optical parameters and band energy values such as the coefficient for bound electron contribution Q _bulk?=?2?×?10²⁴, gap energy E _g?=?1.95 eV, Fermi energy E _F?=?2.15 eV, and damping constant for bound electrons γ _b?=?1.15?×?10¹⁴ Hz. With both size-dependent contributions to the dielectric function, extinction spectra of copper Nps in the subnanometer radius range can be calculated using Mie’s theory and its behaviour with size can be analysed. These studies are applied to fit experimental extinction spectra of very small spherical core–shell Cu–Cu₂O Nps generated by ultrafast laser ablation of a solid target in water. Theoretical calculations for subnanometric core radius are in excellent agreement with experimental results obtained from core–shell colloidal Nps. From the fitting, it is possible determining core radius and shell thickness of the Nps, showing that optical extinction spectroscopy is a good complementary technique to standard high-resolution electron microscopy for sizing spherical nanometric-subnanometric Nps.     

Power saturation of ESR signal in ammonium tartrate exposed to 60Co gamma-ray photons, electrons and protons

In this paper we present an investigation of the electron spin resonance (ESR) line shape of ammonium tartrate (AT) dosimeters exposed to radiation with different linear energy transfer (LET). We exposed our dosimeters to gamma-ray photons ((60)Co), 7 MeV and 14 MeV initial energy electrons, and 19.3 MeV initial energy protons. The differences in the power saturation behavior of ESR spectra of AT irradiated with photons, electrons and protons could be correlated to the effective LET of the radiation beams. We analyzed the behavior of peak-to-peak amplitude as a function of microwave power, and we developed a fitting procedure that permits us to obtain the dependence of the homogeneity parameter of the line shape on the LET of the radiation using the Castner saturation theory. This simple procedure allows us to distinguish the LET of the radiation beam.     

Effect of Cosmological Neutrinos on Discrimination Between the Two Enantiomers of a Chiral Molecule

In the framework of an extraterrestrial origin of biological homochirality, universal mechanisms are of particular interest. In this sense we consider the weak parity-violating neutrino-electron interaction through weak charged currents W ^± between the relic flux of cosmological neutrinos and the electrons of a chiral molecule. We use the known theoretical result of the split in energy of the two helicity sates of an electron in the cosmic neutrino bath, due to weak charged currents. In the case that electrons of a chiral molecule are submitted to a helicoidal potential due to the nuclear conformation, these electrons have opposite helicities for the two enantiomers of the molecule and consequently the mentioned neutrino-electron interaction would produce a splitting in energy between the two enantiomers. An estimation of this energy for the case of a single electron yields a small value of the order of 10⁻²⁶ eV. This value results amplified by the contribution of all the molecular electrons having helicity and other possible mechanisms.     

Qualitative and quantitative analyses of the decomposition products that arise from the exposure of thymine to monochromatic ultrasoft X rays and 60Co gamma rays in the solid state

HPLC analyses of condensed thymine irradiated with monochromatic synchrotron ultrasoft X rays in the energy region around nitrogen and oxygen K-shell edges were performed. Cobalt-60 gamma rays were used as a reference radiation. The radiation chemical dose response of each separated thymine decomposition product was also determined. Uracil (U), 5-(hydroxymethyl)uracil (HMU), 5,6-dihydrothymine (DHT), 5-formyluracil (foU) and four main unknown products were found in the HPLC chromatogram of the sample irradiated with ultrasoft X rays in vacuo. Similar spectra of the products were also found in the gamma-ray experiment; however, some unknown products that appeared after elution of the thymine peak were significantly larger than those in the ultrasoft X- ray experiment. This result indicates the difference in radiation quality. The G value of DHT produced by gamma radiation was 10 times larger than those produced by the ultrasoft X- ray photons with energies of 395 and 407 eV corresponding to below and on the nitrogen K-shell edge, respectively. This result suggests that the differences in the photon energy and/ or in the energy spectra of the secondary electron between ultrasoft X rays and gamma rays are causing differences in the process of the radiation chemistry. Moreover, the yields of all the thymine decomposition products induced by 538 eV photons (oxygen K-shell edge) were significantly smaller than those induced by photons around the nitrogen K-shell edge. The K-shell excitation of oxygen in thymine may efficiently promote the production of small thymine fragments susceptible to desorption from the sample.     

Method of simplified kinetic equation for simulating the electron kinetics in dense gases in strong electric fields

A method for calculating the electron kinetics in dense gases in strong electric fields is developed. The method differs from the forward-backward approximation proposed by Ra?zer and Shne?der for “high-energy” electrons in that it introduces the effective cosines of the scattering angles with respect to the electric force, μ⁺(?, E) and μ^?(?, E), which are different from +1 and ?1, as in the forward-backward approximation. The method was implemented numerically for atmospheric-pressure helium and molecular nitrogen for fields in the ranges 10–200 and 50–800 kV/cm, respectively. The cosines μ⁺(?, E) and μ^?(?, E), the frequency of “fatal” collisions making high-energy electrons to pass from an acceleration regime to a deceleration one, and the rate at which the electrons leave the low-energy reservoir with energies of ≤15 eV for nitrogen and of ≤20 eV for helium are calculated.     

Double layer in a cylindrical hollow-cathode discharge

A dc cylindrical coaxial glow discharge with an inner grid anode has been studied. The region between the two electrodes is seen dark, while a brightly glowing region forms inside the grid anode up to the center. The current-voltage characteristic of a dc cylindrical glow discharge in nitrogen is similar to that of a normal glow discharge, while the normal glow discharge voltage decreases with increasing pressure. The minimum plasma potentials are observed in the hollow cathode region due to the accumulation of electrons at the back of the grid anode. At the center, some of the passed electrons are converged, so their potential is decreased. These electrons have a sufficient time to be redistributed to form one group with a Maxwellian electron energy distribution function. The electron temperature measured by electric probes varies from 1.6 to 3.6 eV, while the plasma density varies from 3.9 × 10¹⁶ to 7 × 10¹³ m⁻³, depending on the discharge current and probe position. The plasma density increases as the electrons move radially from the grid toward the central region, while their temperature decreases.     

Dielectric Dispersion in Ga<Subscript>2</Subscript>S<Subscript>3</Subscript> Thin Films

In this work, the structural, compositional, optical, and dielectric properties of Ga₂S₃ thin films are investigated by means of X-ray diffraction, scanning electron microscopy, energy dispersion X-ray analysis, and ultraviolet—visible light spectrophotometry. The Ga₂S₃ thin films which exhibited amorphous nature in its as grown form are observed to be generally composed of 40.7 % Ga and 59.3 % S atomic content. The direct allowed transitions optical energy bandgap is found to be 2.96 eV. On the other hand, the modeling of the dielectric spectra in the frequency range of 270–1,000 THz, using the modified Drude-Lorentz model for electron-plasmon interactions revealed the electrons scattering time as 1.8 (fs), the electron bounded plasma frequency as ~0.76–0.94 (GHz) and the reduced resonant frequency as 2.20–4.60 ×10¹⁵ (Hz) in the range of 270–753 THz. The corresponding drift mobility of electrons to the terahertz oscillating incident electric field is found to be 7.91 (cm ²/Vs). The values are promising as they nominate the Ga₂S₃ thin films as effective candidates in thin-film transistor and gas sensing technologies.     

Radiation therapy at compact Compton sources

The principle of the compact Compton source is presented briefly. In collision with an ultrarelativistic electron bunch a laser pulse is back-scattered as hard X-rays. The radiation cone has an opening of a few mrad, and the energy bandwidth is a few percent. The electrons that have an energy of the order of a few tens of MeV either circulate in storage ring, or are injected to a linac at a frequency of 10–100 MHz. At the interaction point the electron bunch collides with the laser pulse that has been amplified in a Fabry-Perot resonator. There are several machines in design or construction phase, and projected fluxes are 10¹² to 10¹⁴ photons/s. The flux available at 80 keV from the ThomX machine is compared with that used in the Stereotactic Synchrotron Radiation Therapy clinical trials. It is concluded that ThomX has the potential of serving as the radiation source in future radiation therapy programs, and that ThomX can be integrated in hospital environment.     

Laser-generated plasmas at INFN-LNS

Hot plasmas can be generated by fast and intense laser pulses ablating solids placed in vacuum. A Nd:Yag laser operating at the fundamental and second harmonics with 9-ns pulses (maximum energy of 900 mJ) focused on metallic surfaces produces high ablation yields of the order of μg/pulse and dense plasma that expands adiabatically at supersonic velocity along the normal to the target surface. The plasma emits neutral and charged particles. Charge states up to 10+ have been measured in heavy elements ablated with intensities of the order of 10¹⁰ W/cm². The ion temperature of the plasma is evaluated from the ion energy distributions measured with an ion energy analyzer. The electron temperature is measured through Faraday cups placed at the end of long drift tubes by using time-of-flight technique. The neutral temperature is measured with a special mass quadrupole spectrometer placed along the normal to the target surface. The plasma temperature increases with the laser pulse intensity. The ion temperature reaches values of the order of 400 eV, the electron temperature is of the order of 1 keV for hot electrons and 0.1 eV for thermal electrons, and the neutral temperature is of the order of 200 eV. The experimental apparatus, the diagnostic techniques, and the procedures for the plasma temperature characterization will be presented and discussed in detail. Published in Russian in Fizika Plazmy, 2006, Vol. 32, No. 6, pp. 558–564. The text was submitted by the authors in English.     

High-energy electron irradiation of proteins and nucleic acids: collisional stopping power and average energy loss

Inactivation of proteins due to the direct action of ionizing radiation and the electron energy loss spectra of organic materials indicate that an average of 60–66 eV of energy is lost from high energy electrons in each inelastic collision with target molecules. The average energy loss per inelastic collision with high energy electrons in solid, carbon-based materials, proteins and nucleic acids is calculated from mass collisional stopping powers and empirical total inelastic cross-sections. Bragg’s Additivity Law is used for the calculation of the mean excitation energy of molecules. For simple organic compounds, the calculated average energy loss is close to that obtained by direct observation of the energy loss suffered by electrons as they pass through thin films of organic material. The density effect correction for the rate of energy loss, important in the more complex case of proteins irradiated with 10 MeV electrons, is determined using the comparable mass collisional stopping power of water and proteins. In this manner, a value is obtained for the average energy per inelastic collision of high energy electrons with proteins, which is similar to the average energy per inactivating event of proteins. Analogous calculations for nucleic acids are also presented.     

Evaluation of energy deposition and secondary particle production in proton therapy of brain using a slab head phantom

AimEvaluation of energy deposition of protons in human brain and calculation of the secondary neutrons and photons produced by protons in proton therapy.BackgroundRadiation therapy is one of the main methods of treating localized cancer tumors. The use of high energy proton beam in radiotherapy was proposed almost 60 years ago. In recent years, there has been a revival of interest in this subject in the context of radiation therapy. High energy protons suffer little angular deflection and have a well-defined penetration range, with a sharp increase in the energy loss at the end of their trajectories, namely the Bragg peak.Materials and methodsA slab head phantom was used for the purpose of simulating proton therapy in brain tissue. In this study simulation was carried out using the Monte Carlo MCNPX code.ResultsBy using mono energetic proton pencil beams, energy depositions in tissues, especially inside the brain, as well as estimating the neutron and photon production as a result of proton interactions in the body, together with their energy spectra, were calculated or obtained. The amount of energy escaped from the head by secondary neutrons and photons was determined.ConclusionsIt was found that for high energy proton beams the amount of escaped energy by neutrons is almost 10 times larger than that by photons. We estimated that at 110 MeV beam energy, the overall proton energy “leaked” from the head by secondary photons and neutrons to be around 1%.     

Electron transfer between horse heart and Candida krusei cytochromes c in the free and bound states

Electron transfer between horse heart and Candida krusei cytochromes c in the free and phosvitin-bound states was examined by difference spectrum and stopped-flow methods. The difference spectra in the wavelength range of 540–560 nm demonstrated that electrons are exchangeable between the cytochromes c of the two species. The equilibrium constants of the electron transfer reaction for the free and phosvitin-bound forms, estimated from these difference spectra, were close to unity at 20°C in 20 mM Tris-HCl buffer (pH 7.4). The electron transfer rate for free cytochrome c was (2–3) · 10⁴ M^?1 · s^?1 under the same conditions. The transfer rate for the bound form increased with increase in the binding ratio at ratios below half the maximum, and was almost constant at higher ratios up to the maximum. The maximum electron exchange rate was about 2 · 10⁶ M^?1 · s^?1, which is 60–70 times that for the free form at a given concentration of cytochrome c. The activation energy of the reaction for the bound cytochrome c was equal to that for the free form, being about 10 kcal/mol. The dependence of the exchange rate on temperature, cytochrome c concentration and solvent viscosity suggests that enhancement of the electron transfer rate between cytochromes c on binding to phosvitin is due to increase in the collision frequency between cytochromes c concentrated on the phosvitin molecule.     

Two types of double-strand breaks in electron and photon tracks and their relation to exchange-type chromosome aberrations

Yields of DNA double-strand breaks (dsb), i. e. the average number of dsb, N–, per relative molar mass, M _r , and dose, D, produced by electrons and photons in the energy range 50 eV – 1 MeV were calculated. The experimental data of dsb induction by ultrasoft x-rays and by photons agree well with the calculated yields of dsb as a function of photon energy. The dsb are classified into simple and complex ones. Energy transfers of less than about 200 eV producing at least two ionizations generate mainly simple dsb, while low-energy electrons with an initial energy between 200 and 500 eV induce preferentially complex dsb. Assuming that dsb is the main DNA lesion leading to exchange-type chromosome aberrations (etca), three different mechanisms have to be considered: 1) complex dsb on its own; 2) interaction between two dsb induced by the same primary particle; and 3) interaction between two dsb induced by different primary particles. Mechanisms 1) and 2) produce a linear term, whereas mechanism 3) leads to a quadratic term for the yield of etca. The sum of contributions 1) and 2) to the yield of dicentrics describes fairly well the non-trivial structure of the experimental data. The results suggest that interaction between complex dsb does not contribute significantly to the formation of dicentrics via mechanism 3). Received: 30 July 1995 / Accepted in revised form: 28 March 1996     

Self-consistent band theoretic models of DNA

Previously investigated models are further elaborated here, with particular application to DNA. Accepting the premise that DNA may be characterized as an intrinsic semiconductor with a band gap energy of approximately 2 eV, a self-consistent method reveals term energies of reasonable magnitude. The treatment also determines that the mobilities are of the order of 10 to 10² cm²/V sec, indicating delocalization of the electrons with respect to the base pairs.     

Radiation quality and rat motor performance

The effects of bremsstrahlung, electron, gamma, and neutron radiations were investigated on the motor performance of male Sprague-Dawley rats. Rats were irradiated at a midline tissue dose rate of 20 Gy/min +/- 1 with one of the following: 18.6-MeV electrons (N = 40) or 18.1-MVp bremsstrahlung (N = 57) from a linear accelerator, 60Co 1.25-MeV gamma-ray photons (N = 48), or reactor neutrons at 1.67 MeV tissue-kerma weighted-mean energy (N = 43). Radiation effects were determined by establishing median effective doses (ED50) for rats trained on an accelerod, a shock-avoidance motor performance test. ED50's were based on 10-min postexposure performance. The ED50's were 61 Gy for electrons, 81 Gy for bremsstrahlung, 89 Gy for gamma-ray photons, and 98 Gy for neutrons. In terms of relative biological effectiveness to produce early performance decrement (10 min from the start of irradiation), significant differences existed between the electrons and the other three fields and between the bremsstrahlung and neutron fields. These differences could not be explained by macroscopic dose distribution patterns in the irradiated animals. The data imply that different radiation qualities are not equally effective at disrupting performance, with high-energy electrons being the most effective and neutrons the least.     

Radiation quality of tritium beta-rays: microdosimetric distributions for nanometer-size targets in a medium containing tritiated water

To elucidate the characteristics of the action of tritium beta-rays, the following parameters are derived: electron slowing down spectra of primary electrons (beta-rays) and delta-rays in a medium containing tritiated water; and frequency distributions for the microdosimetric quantity j (number of effective primary events per track per target), fj, for nanometer-size targets exposed to tritiated water. Features of the radiation quality of tritium beta-rays are discussed by comparing the present results with those for 60Co gamma-rays and 7 MeV electrons. It is concluded that, although tritium beta-rays, 60Co gamma-rays, and 7 MeV electrons are classified as the same low l.e.t. radiation, the radiation quality of tritium beta-rays is considerably different from those of 60Co gamma-rays and 7 MeV electrons, and has specific features such as a high average l.e.t., a small total electron fluence per unit absorbed dose, and a different microdosimetric distribution, fj, for nanometer-size targets.     

Electron and photon spectra for three gadolinium-based cancer therapy approaches

Some recent neutron capture therapy research has focused on using compounds containing the element gadolinium, which produces internal conversion and Auger cascade electrons. The low-energy, short-range Auger electrons are absorbed locally and increase cell killing dramatically as the gadolinium compounds are introduced into the cell nucleus and bind to the DNA. Detailed electron and photon spectra are needed for biophysical modeling and Monte Carlo calculations of damage to DNA. This paper presents calculated electron and photon spectra for three cases: thermal neutron absorption by (157)Gd, the beta-particle decay of (159)Gd, and the K-shell photoelectric event in gadolinium. The Monte Carlo sampling of atomic and nuclear transitions for each of the three cases was used to calculate a large number of representative decays. The sampled decays were used to determine average emissions and energy deposited in small spheres of tissue. The kinetic energy nuclear recoil from gamma-ray and electron emissions was calculated and found to be more than 10 eV for 26% of all (157)Gd neutron capture reactions.     

Characteristics of Bremsstrahlung emissions of 177Lu, 188Re,and 90Y for SPECT/CT quantification in radionuclide therapy

PurposeBeta particles emitted by radioisotopes used in targeted radionuclide therapies (TRT) create Bremsstrahlung (BRS) which may affect SPECT quantification when imaging these isotopes. The purpose of the current study was to investigate the characteristics of Bremsstrahlung produced in tissue by three β-emitting radioisotopes used in TRT.MethodsMonte Carlo simulations of ¹⁷⁷Lu, ¹⁸⁸Re, and ⁹⁰Y sources placed in water filled cylinders were performed. BRS yields, mean energies and energy spectra for (a) all photons generated in the decays, (b) photons that were not absorbed and leave the cylinder, and (c) photons detected by the camera were analyzed. Next, the results of simulations were compared with those from experiments performed on a clinical SPECT camera using same acquisition conditions and phantom configurations as in simulations.ResultsSimulations reproduced relatively well the shapes of the measured spectra, except for ⁹⁰Y which showed an overestimation in the low energy range. Detailed analysis of the results allowed us to suggest best collimators and imaging conditions for each of the investigated isotopes. Finally, our simulations confirmed that the BRS contribution to the energy spectra in quantitative imaging of ¹⁷⁷Lu and ¹⁸⁸Re could be ignored.ConclusionsFor ¹⁷⁷Lu and ¹⁸⁸Re, BRS contributes only marginally to the total spectra recorded by the camera. Our analysis shows that MELP and HE collimators are the best for imaging these two isotopes. For ⁹⁰Y, HE collimator should be used.     

Mechanism of radiation-induced reactions in aqueous solution of coumarin-3-carboxylic acid: Effects of concentration,gas and additive on fluorescent product yield

Abstract

The radiation-induced reactions of a water-soluble coumarin derivative, coumarin-3-carboxyl acid (C3CA), have been investigated in aqueous solutions by pulse radiolysis with a 35 MeV electron beam, final product analysis following ⁶⁰Co γ-irradiations and deterministic model simulations. Pulse radiolysis revealed that C3CA reacted with both hydroxyl radicals (^?OH) and hydrated electrons (e^? _aq) with near diffusion-controlled rate constants of 6.8 × 10⁹ and 2.1 × 10¹⁰ M^?1 s^?1, respectively. The reactivity of C3CA towards O₂^{? ?} was not confirmed by pulse radiolysis. Production of the fluorescent molecule, 7-hydroxy-coumarin-3-carboxylic acid (7OH-C3CA), was confirmed by final product analysis with a fluorescence spectrometer coupled to a high performance liquid chromatography (HPLC) system. Production yields of 7OH-C3CA following ⁶⁰Co γ-irradiations depended on the irradiation conditions and ranged from 0.025 to 0.18 (100 eV) ^?1. Yield varied with saturating gas, additive and C3CA concentration, implying the presence of at least two pathways capable of providing 7OH-C3CA as a stable product following the scavenging reaction of C3CA with ^?OH, including a peroxidation/elimination sequence and a disproportionation pathway. A reaction mechanism for the two pathways was proposed and incorporated into a deterministic simulation, showing that the mechanism can explain experimentally measured 7OH-C3CA yields with a constant conversion factor of 4.7% from ^?OH scavenging to 7OH-C3CA production, unless t-BuOH was added.