Fluctuations in measured radioactive decay rates inside a modified Faraday cage: Correlations with space weather.

For several years, reports have been published about fluctuations in measured radioactive decay time-series and in some instances linked to astrophysical as well as classical environmental influences. Anomalous behaviors of radioactive decay measurement and measurement of capacitance inside and outside a modified Faraday cage were documented by our group in previous work. In the present report, we present an in-depth analysis of our measurement with regard to possible correlations with space weather, i.e. the geomagnetic activity (GMA) and cosmic-ray activity (CRA). Our analysis revealed that the decay and capacitance time-series are statistically significantly correlated with GMA and CRA when specific conditions are met. The conditions are explained in detail and an outlook is given on how to further investigate this important finding. Our discovery is relevant for all researchers investigating radioactive decay measurements since they point out that the space weather condition during the measurement is relevant for partially explaining the observed variability.

Study of the photoneutron generation caused by a LinAc Beryllium window with a 6 MeV treatment beam.

In most conventional radiation therapy treatments, special attention is payed for neutron contamination when working with energy beams above 8 MeV and generally it is only considered for shielding requirements, not for dose study in patients or employees. The present work is focused on studying the unwanted generated photoneutrons in a Medical Linear Accelerator (LinAc) Varian TrueBeam using a 6 MeV radiation treatment beam. To that, Monte Carlo (MC) simulation code MCNP6.1.1 was used. This version of the code allows the use of unstructured mesh geometries as a novelty, offering more reliable results and higher speed computation. The particularity of the studied LinAc is the presence of a beryllium filter at the treatment head. Since Beryllium causes photonuclear reactions (γ, n) at energies much lower than other LinAc composing materials, this work aims to analyze if this type of units, when using low energy treatment beams (6 MeV), produce neutron pollution and to ensure that this unwanted radiation can be considered negligible.

Brachytherapy organ dose estimation using Monte Carlo simulations of realistic patient models.

Radiation Therapy Planning Systems (RTPS) currently used in hospitals contain algorithms based on deterministic simplifications that do not properly consider electrons lateral transport in the areas where there are changes of density, and as a result, erroneous dose predictions could be produced. According to this, the present work proposes the use of Monte Carlo method in brachytherapy planning systems, which could affect positively on the radiotherapy treatment planning, since it provides results that are more accurate and takes into account the in homogeneities density variations. This paper presents a Monte Carlo (MC) simulation of a brachytherapy prostate treatment with I-125 seeds, using the latest version of MCNP, v.6.1.1. To that, a 3D model of the anatomy of a real anonymized patient is created from the segmentation of Computed Tomography (CT) images. Treatment over this 3D model is simulated and the dose given to each organ is obtained. These doses are compared with those calculated by deterministic planning system used in hospital demonstrating the effectiveness of MC method in the planning of brachytherapy treatments, because of not only the results precision but also regarding the affordable computing times.

Dose calculation in computerized tomography.

The purpose of this work is to develop an automatic methodology to obtain the dose received by a patient, (classified by organs), after being subjected to ionizing radiation because of CT images acquisition. The methodology starts from CT images, these images are automatically segmented and voxelized taking into account the CT numbers in order to obtain a 3D model used in Monte Carlo Simulations to calculate the dose inside the patient.Monte Carlo codes used in this work are MCNP.6.1.1, (whose results will be taken as a reference values) and MC-GPU, which appears to be a good candidate to be implemented in the methodology since his GPU parallelization offers a high speed calculation. Results show good agreement between simulated values obtained by MC-GPU and MCNP6.1.1.

Experimental validation of neutron activation simulation of a varian medical linear accelerator.

This work presents a Monte Carlo simulation using the last version of MCNP, v. 6.1.1, of a Varian CLinAc emitting a 15MeV photon beam. The main objective of the work is to estimate the photoneutron production and activated products inside the medical linear accelerator head. To that, the Varian LinAc head was modelled in detail using the manufacturer information, and the model was generated with a CAD software and exported as a mesh to be included in the particle transport simulation. The model includes the transport of photoneutrons generated by primary photons and the (n, γ) reactions which can result in activation products. The validation of this study was done using experimental measures. Activation products have been identified by in situ gamma spectroscopy placed at the jaws exit of the LinAc shortly after termination of a high energy photon beam irradiation. Comparison between experimental and simulation results shows good agreement.

Air radon concentration decrease in a waste water treatment plant.

(222)Rn is a naturally occurring gas created from the decay of (226)Ra. The long-term health risk of breathing radon is lung cancer. One particular place where indoor radon concentrations can exceed national guidelines is in wastewater treatment plants (WWTPs) where treatment processes may contribute to ambient airborne concentrations. The aim of this paper was to study the radon concentration decrease after the application of corrective measures in a Spanish WWTP. According to first measures, air radon concentration exceeded International Commission Radiologica1 Protection (ICRP) normative (recommends intervention between 400 and 1000 Bq m(-3)). Therefore, the WWTP improved mechanical forced ventilation to lower occupational exposure. This measure allowed to increase the administrative controls, since the limitation of workers access to the plant changed from 2 h d(-1) (considering a maximum permissible dose of 20 mSv y(-1) averaged over 5 y) to 7 h d(-1).

Neutron distribution and induced activity inside a Linac treatment room.

Induced radioactivity and photoneutron contamination inside a radiation therapy bunker of a medical linear accelerator (Linac) is investigated in this work. The Linac studied is an Elekta Precise electron accelerator which maximum treatment photon energy is 15 MeV. This energy exceeds the photonuclear reaction threshold (around 7 MeV for high atomic number metals). The Monte Carlo code MCNP6 has been used for quantifying the neutron contamination inside the treatment room for different gantry rotation configuration. Walls activation processes have also been simulated. The approach described in this paper is useful to prevent the overexposure of patients and medical staff.

Neutron activation processes simulation in an Elekta medical linear accelerator head.

Monte Carlo estimation of the giant-dipole-resonance (GRN) photoneutrons inside the Elekta Precise LINAC head (emitting a 15 MV photon beam) were performed using the MCNP6 (general-purpose Monte Carlo N-Particle code, version 6). Each component of LINAC head geometry and materials were modelled in detail using the given manufacturer information. Primary photons generate photoneutrons and its transport across the treatment head was simulated, including the (n, γ) reactions which undergo activation products. The MCNP6 was used to develop a method for quantifying the activation of accelerator components. The approach described in this paper is useful in quantifying the origin and the amount of nuclear activation.

Comparison between transmission and scattering spectrum reconstruction methods based on EPID images.

Numerous improved physics-based methods for Linac photon spectra reconstruction have been published; some of them are based on transmission data analysis and others on scattering data. In this work, the two spectrum unfolding approaches are compared in order to experimentally validate its robustness and to determine which is the optimal methodology for application on a clinical quality assurance routine. Both studied methods are based on EPID images generated when the incident photon beam impinges onto plastic blocks. The distribution of transmitted/scatter radiation produced by this object centered at the beam field size was measured. Measurements were performed using a 6 MeV photon beam produced by the linear accelerator. The same radiation distribution conditions were also simulated with Monte Carlo code for a series of monoenergetic identical geometry photon beams for both cases. Two systems of linear equations were generated to combine the polyenergetic EPID measurements with the monoenergetic simulation results. Regularization techniques were applied to solve the systems for obtaining the incident photon spectrum. We present a comparison between the well-known photon Spectral Reconstruction based on Transmission Data (Trans-based) technology and the Spectral Reconstruction based on Scattering Data (Scatt-based), which we both developed using EPID images. It is shown that Trans-based reconstruction results display much better agreement with photon spectrum theoretical predictions.

Uncertainty and sensitivity analysis of the effect of the mean energy and FWHM of the initial electron fluence on the Bremsstrahlung photon spectra of linear accelerators.

A calculation of the correct dose in radiation therapy requires an accurate description of the radiation source because uncertainties in characterization of the linac photon spectrum are propagated through the dose calculations. Unfortunately, detailed knowledge of the initial electron beam parameters is not readily available, and many researchers adjust the initial electron fluence values by trial-and-error methods. The main goal of this work was to develop a methodology to characterize the fluence of initial electrons before they hit the tungsten target of an Elekta Precise medical linear accelerator. To this end, we used a Monte Carlo technique to analyze the influence of the characteristics of the initial electron beam on the distribution of absorbed dose from a 6 MV linac photon beam in a water phantom. The technique is based on calculations with Software for Uncertainty and Sensitivity Analysis (SUSA) and Monte Carlo simulations with the MCNP5 transport code. The free parameters used in the SUSA calculations were the mean energy and full-width-at-half-maximum (FWHM) of the initial electron distribution. A total of 93 combinations of these parameters gave initial electron fluence configurations. The electron spectra thus obtained were used in a simulation of the electron transport through the target of the linear accelerator, which produced different photon (Bremsstrahlung) spectra. The simulated photon spectra were compared with the 6-MV photon spectrum provided by the linac manufacturer (Elekta). This comparison revealed how the mean energy and FWHM of the initial electron fluence affect the spectrum of the generated photons. This study has made it possible to fine-tune the examined electron beam parameters to obtain the resulted absorbed doses with acceptable accuracy (error<1%).

Assessment in vitro of cytogenetic and genotoxic effects of propolis on human lymphocytes.

We evaluated the genetic damage by ethanolic extract of propolis (EEP) induced to human lymphocytes which were exposed to increasing concentrations (0-2000µgml(-1)). The results indicated that EEP reduced significantly the mitotic index (MI) and proliferation index (PI) when high concentrations of EEP were used. Sister chromatid exchange (SCE) rates indicated that EEP could have genotoxic effects at high concentrations. Exposure of the cells to the amount of ethanol used as solvent did not alter either the MI and cell proliferation kinetics (CPK), or the rate of SCE. The results showed: (a) statistical increase in the percentage the cells with CAs and in the frequency of SCE at the highest concentrations, (b) a decrease in MI and in the CPK values was observed, (c) no effect was noticed in negative controls. In conclusion, it can be assumed that high concentrations of EEP have a cyto and genotoxic effect, in vitro, for human peripheral lymphocytes.

X-ray simulation with the Monte Carlo code PENELOPE. Application to Quality Control.

A realistic knowledge of the energy spectrum is very important in Quality Control (QC) of X-ray tubes in order to reduce dose to patients. However, due to the implicit difficulties to measure the X-ray spectrum accurately, it is not normally obtained in routine QC. Instead, some parameters are measured and/or calculated. PENELOPE and MCNP5 codes, based on the Monte Carlo method, can be used as complementary tools to verify parameters measured in QC. These codes allow estimating Bremsstrahlung and characteristic lines from the anode taking into account specific characteristics of equipment. They have been applied to simulate an X-ray spectrum. Results are compared with theoretical IPEM 78 spectrum. A sensitivity analysis has been developed to estimate the influence on simulated spectra of important parameters used in simulation codes. With this analysis it has been obtained that the FORCE factor is the most important parameter in PENELOPE simulations. FORCE factor, which is a variance reduction method, improves the simulation but produces hard increases of computer time. The value of FORCE should be optimized so that a good agreement of simulated and theoretical spectra is reached, but with a reduction of computer time. Quality parameters such as Half Value Layer (HVL) can be obtained with the PENELOPE model developed, but FORCE takes such a high value that computer time is hardly increased. On the other hand, depth dose assessment can be achieved with acceptable results for small values of FORCE.

MCNP5 Monte Carlo simulation of amorphous silicon EPID dosimetry from MLC radiation therapy treatment beams.

The present work is focused on a MCNP Monte Carlo (MC) simulation of a multi-leaf collimator (MLC) radiation therapy treatment unit including its corresponding Electronic Portal Imaging Device (EPID). We have developed a methodology to perform a spatial calibration of the EPID signal to obtain dose distribution using MC simulations. This calibration is based on several images acquisition and simulation considering different thicknesses of solid water slabs, using a 6 MeV photon beam and a square field size of 20 cm x 20 cm. The resulting relationship between the EPID response and the MC simulated dose is markedly linear. This signal to dose EPID calibration was used as a dosimetric tool to perform the validation of the MLC linear accelerator MCNP model. Simulation results and measurements agreed within 2% of dose difference. The methodology described in this paper potentially offers an optimal verification of dose received by patients under complex multi-field conformal or intensity-modulated radiation therapy (IMRT).

Concentration-Dependent Protection by Ethanol Extract of Propolis against γ-Ray-Induced Chromosome Damage in Human Blood Lymphocytes.

Radioprotection with natural products may be relevant to the mitigation of ionizing radiation-induced damage in mammalian systems; in this sense, propolis extracts have shown effects such as antioxidant, antitumoral, anti-inflammatory, and immunostimulant. We report for the first time a cytogenetic study to evaluate the radioprotective effect, in vitro, of propolis against radiation-induced chromosomal damage. Lymphocytes were cultured with increasing concentrations of ethanol extract of propolis (EEP), including 20, 40, 120, 250, 500, 750, 1000, and 2000 µg mL(-1) and then exposed to 2 Gy γ-rays. A significant and concentration-dependent decrease is observed in the frequency of chromosome aberrations in samples treated with EEP. The protection against the formation of dicentrics was concentration-dependent, with a maximum protection at 120 µg mL(-1) of EEP. The observed frequency of dicentrics is described as negative exponential function, indicating that the maximum protectible fraction of dicentrics is approximately 44%. Free radical scavenging and antioxidant activities are the mechanisms that these substances use to protect cells from ionizing radiation.

Parametric study of the X-ray primary spectra obtained with the MTSVD unfolding method.

The modified truncated singular value decomposition (MTSVD) unfolding method is applied to obtain primary spectra for X-ray tubes in radiodiagnostic. Three parameters - voltage, anode angle and filter thickness - of the tube are tested. Unfolded spectra are compared with theoretical extracted from IPEM-78 catalogue. A 2σ error criterion is applied to assess the minimum variations in tested parameters that permits distinguishing between close spectra.

Biological and physical methods for risk estimation in interventional radiology: a detrimental effect approach.

Interventional radiologists and staff members are frequently exposed to the effects of direct and scattered radiation, which undergo in deterministic effects (radiodermitis, aged skin, cataracts, telangiectasia in nasal region, vasocellular epitelioms, hands depilation) and/or stochastic ones (cancer incidence). A methodology has been proposed for estimating the radiation risk or detriment from a group of six exposed interventional radiologists of the Hospital Universitario La Fe (Valencia, Spain), which had developed general exposition symptoms attributable to deterministic effects of ionizing radiation. Equivalent doses have been periodically registered using termoluminiscence dosimeters (TLD's) and wrist dosimeters, H(p)(10) and H(p)(0.07), respectively, and estimated through the observation of translocations in lymphocytes of peripheral blood (biological methods), by extrapolating the yield of translocations to their respective dose-effect curves. The software RADRISK has been applied for estimating radiation risks in these occupational radiation exposures. The minimum and maximum average excess ratio for skin cancer has been, using wrist physical doses, of [1.03 × 10(-3), 5.06 × 10(-2)], concluding that there is not an increased risk of skin cancer incidence. The minimum and maximum average excess ratio for leukemia has been, using TLD physical doses, of [7.84 × 10(-2), 3.36 × 10(-1)], and using biological doses, of [1.40 × 10(-1), 1.51], which is considerably higher than incidence rates, showing an excess radio-induced risk of leukemia in the group under study. Finally, the maximum radiological detriment in the group, evaluated as the total number of radio-induced cancers using physical dosimetry, has been of 2.18 per 1000 person-year (skin and leukemia), and using biological dosimetry of 9.20 per 1000 PY (leukemia). As a conclusion, this study has provided an assessment of the non-deterministic effects (rate of radio-induced cancer incidence) attributable to the group under study due to their professional activity.

6 and 15 MeV photon spectra reconstruction using an unfolding depth dose gradient methodology.

An accurate knowledge of the spectral distribution emission is essential for precise dose calculations in radiotherapy treatment planning. Reconstruction of photon spectra emitted by medical accelerators from measured depth dose distributions in a water cube is an important tool for commissioning a Monte Carlo treatment planning system. However, the reconstruction problem is an inverse radiation transport function which is poorly conditioned and its solution may become unstable due to small perturbations in the input data. In this paper we present a more stable spectral reconstruction method which can be used to provide an independent confirmation of source models for a given machine without any prior knowledge of the spectral distribution. This technique involves measuring the depth dose curve in a water phantom and applying an unfolding method using Monte Carlo simulated depth dose gradient curves for consecutives mono-energetic beams. We illustrate this theory to calculate a 6 and a 15 MeV photon beam emitted from an Elekta Precise radiotherapy unit using the gradient of depth dose curves in a cube-shaped water tank.

Uncertainty analysis in the simulation of X-ray spectra in the diagnostic range using the MCNP5 code.

An accurate knowledge of the photonic spectra emitted by X-ray tubes in radiodiagnostics is essential to better estimate the imparted dose to patients and to improve the image quality obtained with these devices. In this work, several X-ray spectra have been simulated using the MCNP5 code to simulate X-ray production in a commercial device. To validate the Monte Carlo results, simulated spectra have been compared to those extracted from the IPEM 78 database. The uncertainty associated to some geometrical features of the tube and its effect on the simulated spectra has been analyzed using the Noether-Wilks formula. This analysis has been focused on the thickness of collimators, filters, shielding and barrel shutter. Furthermore, results show that the uncertainty due to geometrical parameters (0.98% in terms of Root Mean Squared) is higher than the statistical uncertainty associated to the MCNP5 calculations.

Analysis of image quality parameter of conventional and dental radiographic digital images.

The image quality obtained by a radiographic equipment is very useful to characterize the physical properties of the image radiographic chain, in a quality control of the radiographic equipment. In the radiographic technique it is necessary that the evaluation of the image can guarantee the constancy of its quality to carry out a suitable diagnosis. In this work we have designed some radiographic phantoms for different radiographic digital devices, as dental, conventional, equipments with computed radiography (phosphor plate) and direct radiography (sensor) technology. Additionally, we have developed a software to analyse the image obtained by the radiographic equipment with digital processing techniques, as edge detector, morphological operators, statistical test for the detected combinations‥ The design of these phantoms let the evaluation of a wide range of operating conditions of voltage, current and time of the digital equipments. Moreover, the image quality analysis by the automatic software, let study it with objective parameters.