Comparison of SRIM, MCNPX and GEANT simulations with experimental data for thick Al absorbers

Proton computerized tomography deals with relatively thick targets like the human head or trunk. In this case precise analytical calculation of the proton final energy is a rather complicated task, thus the Monte Carlo simulation stands out as a solution. We used the GEANT4.8.2 code to calculate the proton final energy spectra after passing a thick Al absorber and compared it with the same conditions of the experimental data. The ICRU49, Ziegler85 and Ziegler2000 models from the low energy extension pack were used. The results were also compared with the SRIM2008 and MCNPX2.4 simulations, and with solutions of the Boltzmann transport equation in the Fokker-Planck approximation. (C) 2009 Elsevier Ltd. All rights reserved.

Cálculo das funções de anisotropia e radial de dose para braquiterapia com semente de 192Ir em diferentes materiais de interesse em radiologia com o auxílio do programa baseado no método de Monte Carlo MCNPX

The brachytherapy braquiterapia is an x-ray modality radiotherapy in which stamped or halfstamped radioactive sources in format of seeds are used, wires or to one short distance, in contact or implanted to the fabric to be treated. All the treatment modalities require a previous planning. The formalism recommended for calculation of dose was considered by the AAPM for the first TG-43 Report. In it distributions of dose of isolated seeds measured and calculated by Monte Carlo method in water instead of with models half-empiricists. In this work we in accordance with present some preliminary results of the calculation of functions of radial anisotropy and of dose in the distance for seed of 192Ir, wide used in brachytherapy treatments of high tax of dose (HDR), with the aid of the program based on the Monte Carlo method MCNPX v2.50 (Mount Carlo N ParticleeXtended). The materials chosen in the simulation beyond water, had been MS20 and estriado muscle

Viabilidade do Código MCNPX para o Estudo da Técnica de Fluorescência de Raios-X por Reflexão Total

The X-ray fluorescence analysis (XRF) is an important technique for the qualitative and quantitative determination of chemical components in a sample. It is based on measurement of the intensity of the emitted characteristic radiation by the elements of the sample, after being properly excited. One of the modalities of this technique is the total reflection x-ray fluorescence (TXRF). In TXRF, the angle of refraction of the incident beam tends to zero and the refracted beam is tangent to the sample-support interface. Thus, there is a minimum angle of incidence that there is no refracted beam and all the incident radiation undergoes total reflection. As it is implemented in very small samples, in a film format, self-absorption effects should not very relevant. In this study, we evaluated the feasibility of using code MCNPX (Monte Carlo N - Particle eXtended), to simulate a measure implemented by the TXRF technique. In this way, it was verified the quality of response of a system by TXRF spectroscopy using synchrotron radiation as excitation beam for a simple setup, by retrieving the characteristic energies and the concentrations of the elements in the sample. The steps of data processing, after obtaining the excitation spectra, were the same as in a real experiment and included the obtaining of the sensitivity curve for the simulated system. The agreement between the theoretical and simulated values of Ka characteristic energies for different elements was lower than 1 % .The obtained concentration of the elements of the sample had high relatively errors ( between 6 and 60 % ) due mainly to lack of knowing about some realistic physical parameters of the sample , such as density . In this way, this result does not preclude the use of MCNPX code for this type of application

Comparison of proton energy loss in thick absorbers in terms of a reduced calibration curve

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Comparison of GEANT4 simulations with experimental data for thick al absorbers

Proton beams in medical applications deal with relatively thick targets like the human head or trunk. Therefore, relatively small differences in the total proton stopping power given, for example, by the different models provided by GEANT4 can lead to significant disagreements in the final proton energy spectra when integrated along lengthy proton trajectories. This work presents proton energy spectra obtained by GEANT4.8.2 simulations using ICRU49, Ziegler1985 and Ziegler2000 models for 19.68MeV protons passing through a number of Al absorbers with various thicknesses. The spectra were compared with the experimental data, with TRIM/SRIM2008 and MCNPX2.4.0 simulations, and with the Payne analytical solution for the transport equation in the Fokker-Plank approximation. It is shown that the MCNPX simulations reasonably reproduce well all experimental spectra. For the relatively thin targets all the methods give practically identical results but this is not the same for the thick absorbers. It should be noted that all the spectra were measured at the proton energies significantly above 2MeV, i.e., in the so-called Bethe-Bloch region. Therefore the observed disagreements in GEANT4 results, simulated with different models, are somewhat unexpected. Further studies are necessary for better understanding and definitive conclusions. © 2009 American Institute of Physics.

Reduced calibration curve for proton computed tomography

The pCT deals with relatively thick targets like the human head or trunk. Thus, the fidelity of pCT as a tool for proton therapy planning depends on the accuracy of physical formulas used for proton interaction with thick absorbers. Although the actual overall accuracy of the proton stopping power in the Bethe-Bloch domain is about 1%, the analytical calculations and the Monte Carlo simulations with codes like TRIM/SRIM, MCNPX and GEANT4 do not agreed with each other. A tentative to validate the codes against experimental data for thick absorbers bring some difficulties: only a few data is available and the existing data sets have been acquired at different initial proton energies, and for different absorber materials. In this work we compare the results of our Monte Carlo simulations with existing experimental data in terms of reduced calibration curve, i.e. the range - energy dependence normalized on the range scale by the full projected CSDA range for given initial proton energy in a given material, taken from the NIST PSTAR database, and on the final proton energy scale - by the given initial energy of protons. This approach is almost energy and material independent. The results of our analysis are important for pCT development because the contradictions observed at arbitrary low initial proton energies could be easily scaled now to typical pCT energies. © 2010 American Institute of Physics.

Monte Carlo and analytical calculation of lateral deflection of proton beams in homogeneous targets

Proton radiation therapy is a precise form of radiation therapy, but the avoidance of damage to critical normal tissues and the prevention of geographical tumor misses require accurate knowledge of the dose delivered to the patient and the verification of his position demand a precise imaging technique. In proton therapy facilities, the X-ray Computed Tomography (xCT) is the preferred technique for the planning treatment of patients. This situation has been changing nowadays with the development of proton accelerators for health care and the increase in the number of treated patients. In fact, protons could be more efficient than xCT for this task. One essential difficulty in pCT image reconstruction systems came from the scattering of the protons inside the target due to the numerous small-angle deflections by nuclear Coulomb fields. The purpose of this study is the comparison of an analytical formulation for the determination of beam lateral deflection, based on Molière's theory and Rutherford scattering with Monte Carlo calculations by SRIM 2008 and MCNPX codes. © 2010 American Institute of Physics.

Comparison of the GEANT4 releases 8.2 and 9.2 in terms of a pCT reduced calibration curve

The GEANT4 simulations are essential for the development of medical tomography with proton beams pCT. In the case of thin absorbers the latest releases of GEANT4 generate very similar final spectra which agree well with the results of other popular Monte Carlo codes like TRIM/SRIM, or MCNPX. For thick absorbers, however, the disagreements became evident. In a part, these disagreements are due to the known contradictions in the NIST PSTAR and SRIM reference data. Therefore, it is interesting to compare the GEANT4 results with each other, with experiment, and with diverse code results in a reduced form, which is free from this kind of doubts. In this work such comparison is done within the Reduced Calibration Curve concept elaborated for the proton beam tomography. © 2010 IEEE.

Monte Carlo simulation applied in total reflection X-ray fluorescence: Preliminary results

The X-ray Fluorescence (XRF) analysis is a technique for the qualitative and quantitative determination of chemical constituents in a sample. This method is based on detection of the characteristic radiation intensities emitted by the elements of the sample, when properly excited. A variant of this technique is the Total Reflection X-ray Fluorescence (TXRF) that utilizes electromagnetic radiation as excitation source. In total reflection of X-ray, the angle of refraction of the incident beam tends to zero and the refracted beam is tangent to the sample support interface. Thus, there is a minimum angle of incidence at which no refracted beam exists and all incident radiation undergoes total reflection. In this study, we evaluated the influence of the energy variation of the beam of incident x-rays, using the MCNPX code (Monte Carlo NParticle) based on Monte Carlo method. © 2013 AIP Publishing LLC.

Estudo teórico da relação entre número de CT de raios-X e Stopping Power de prótons

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo da interação de prótons com alvos não homogêneos, aplicados a tomografia com feixes de prótons

Pós-graduação em Biologia Geral e Aplicada - IBB

Estudo teórico da contribuição à dose de partículas secundárias geradas por prótons de energias entre 100 e 200 MeV num Phantom de substituto de tecido MS20

The goal of this project is the reproduction, through the simulation code based on the MCNPX (Monte Carlo N-Particle eXtended) v2.50 method, of the proton beam interaction with the material, since, in proton therapy, only the particle ionization and excitation are analyzed and the occurence of nuclear interactive inelastic process are not considered. This work will help the development of studies concerning the contribution to the total dosis of secondary particles generated by nuclear interaction in proton therapy. They are: alpha particles ( ), deuterium(2H), tritium (3H), neutron (n) and helium (3He). A MS20 tissue substitute phantom was used as the target and the energy of the proton beams was within an interest range of 100 to 200MeV. With the results obtained, it was possible to generate graphics which allows the analysis of the dosis deposition relation with and without nuclear interaction, the percentage of secondary particles deposited dosis, the radial dispersion of neutrons in the material, the secondary particles multiplicity, as well as the relation between the secondary particles spectrum with the próton generated spectrum

Estudo teórico das funções de anisotropia e radial de dose para braquiterapia com semente de Ir192 em diferentes materiais de interesse em radiologia

The brachytherapy braquiterapia is an x-ray modality radiotherapy in which stamped or halfstamped radioactive sources in format of seeds are used, wires or to one short distance, in contact or implanted to the fabric to be treated. All the treatment modalities require a previous planning.The formalism recommended for calculation of dose was considered by the AAPM for the first TG-43 Report. In it distributions of dose of isolated seeds measured and calculated by Monte Carlo method in water instead of with models half-empiricists.In this work we in accordance with present some preliminary results of the calculation of functions of radial anisotropy and of dose in the distance for seed of 192Ir, wide used in brachytherapy treatments of high tax of dose (HDR), with the aid of the program based on the Monte Carlo method MCNPX v2.50 (Mount Carlo N ParticleeXtended).The materials chosen in the simulation beyond water, had been MS20 and estriado muscle

Estudo da contribuição do espalhamento nuclear não elástico de prótons na dose depositada em alvos com composição homogênea

Due the differences between interaction physics process with matter for protons and photons, the proton beam tomography (pCT) has some vantages to comparison with conventional tomography. Also it is confirmed that usually pCT has better dose distribution and highest contrast resolution. The pCT allow not only view the internal structure of an object without destroying it, but also directly measure of volume density of electrons. Also it is confirmed that usually pCT has better dose distribution and highest contrast resolution. At the same time, there are many scientific and technical aspects to a detailed study: the capacity and limitations of the pCT methods are not well clarified. Through computations, based on Monte Carlo Method was carried out a detailed study of the contribution of non-elastic nuclear spreading, and together was compared with an analytical model for the deflection angle and the lateral deflection of protons in the target volume. The programs used were SRIM 2006 code and MCNPX v.2.50 code

Estudo teórico da contribuição a dose de partículas secundárias geradas por prótons de energias entre 100-200 MeV num phantom de músculo estriado

In the treatment plans in conventional Proton therapy are considered only the elastic interactions of protons with electrons and/or nuclei, it means, mainly ionization and coulomb excitation processes. As the energy needed to reach the deep tumors should be of several hundred of MeVs, certainly the nuclear inelastic channels are open. Only some previous studies of the contribution of these processes in the full dose have been made towards targets composed of water. In this study will be presented the results of the simulation of the processes of interaction of beams of protons in the range of 100-200 MeV of energy with a cylindrical phantom composed by striated muscle (ICRU), emphasizing in the contribution to total dose due to the deposition of energy by secondary particles alpha (α), deuterium (2H), tritium (3H), neutron (n) and hélio3 (3He), originated by nuclear inelastic processes. The simulations were performed by using the method of Monte Carlo, via the computer code MCNPX v2.50 (Monte Carlo N-Particle eXtended). The results will be shown demonstrated through the graphics of the deposited dose with or without nuclear interaction, the percentual of dose deposited by secondary particles, the radial dispersion of neutrons, as well as the multiplicity of secondary particles