Investigation of radiological properties of some shielding materials on charged and uncharged radiation interaction for neutron generator

NASA Astrophysics Data System (ADS)

Büyükyıldız, Mehmet

2017-04-01

Radiation interaction parameters such as total stopping power, projected range (longitudinal and lateral) straggling, mass attenuation coefficient, effective atomic number (Zeff) and electron density (Neff) of some shielding materials were investigated for photon and heavy charged particle interactions. The ranges, stragglings and mass attenuation coefficients were calculated for the high-density polyethylene(HDPE), borated polyethylene (BPE), brick (common silica), concrete (regular), wood, water, stainless steel (304), aluminum (alloy 6061-O), lead and bismuth using SRIM Monte Carlo software and WinXCom program. In addition, effective atomic numbers (Zeff) and electron densities (Neff) of HDPE, BPE, brick (common silica), concrete (regular), wood, water, stainless steel (304) and aluminum (alloy 6061-O) were calculated in the energy region 10 keV-100 MeV using mass stopping powers and mass attenuation coefficients. Two different methods namely direct and interpolation procedures were used to calculate Zeff for comparison and significant differences were determined between the methods. Variations of the ranges, longitudinal and lateral stragglings of water, concrete and stainless steel (304) were compared with each other in the continuous kinetic energy region and discussed with respect to their Zeffs. Moreover, energy absorption buildup factors (EABF) and exposure buildup factors (EBF) of the materials were determined for gamma rays as well and were compared with each other for different photon energies and different mfps in the photon energy region 0.015-15 MeV.

Evidence for age-related performance degradation of (241)Am foil sources commonly used in UK schools.

PubMed

Whitcher, R; Page, R D; Cole, P R

2014-06-01

The characteristics of alpha radiation have for decades been demonstrated in UK schools using small sealed (241)Am sources. There is a small but steady number of schools who report a considerable reduction in the alpha count rate detected by an end-window GM detector compared with when the source was new. This cannot be explained by incorrect apparatus or set-up, foil surface contamination, or degradation of the GM detector. The University of Liverpool and CLEAPSS collaborated to research the cause of this performance degradation. The aim was to find what was causing the performance degradation and the ramifications for both the useful and safe service life of the sources. The research shows that these foil sources have greater energy straggling with a corresponding reduction in spectral peak energy. A likely cause for this increase in straggling is a significant diffusion of the metals over time. There was no evidence to suggest the foils have become unsafe, but precautionary checks should be made on old sources.

A Student Experiment to Demonstrate the Energy Loss and Straggling of Electrons in Matter.

ERIC Educational Resources Information Center

de Bruin, M.; Huijgen, F. W. J.

1990-01-01

Described is an introductory experiment that allows students to directly observe and measure the linear energy transfer in matter. Illustrated are the experimental setup including the radioactive source, electronic equipment, and the detector; measurement and calculations; and the results. (CW)

High resolution Cerenkov and range detectors for balloon-borne cosmic-ray experiment

NASA Technical Reports Server (NTRS)

Ahlen, S. P.; Cartwright, B. G.; Tarle, G.

1975-01-01

A combination of an active Cerenkov detector and passive range detectors is proposed for the high resolution measurement of isotopic composition in the neighborhood of iron in the galactic cosmic rays. A large area (4,300 sq cm) Cerenkov counter and passive range detectors were tested. Tests with heavy ions (2.1 GeV/amu C-12, 289 MeV/amu Ar-40, and 594 MeV/amu Ne-20) revealed the spatial uniformity of response of the Cerenkov counter to be better than 1% peak-to-peak. Light collection efficiency is independent of projectile energy and incidence angle to within at least 0.5%. Passive Lexan track recorders to measure range in the presence of the nuclear interaction background which results from stopping particles through 0.9 interaction lengths of matter were also tested. It was found that nuclear interactions produce an effective range straggling distribution only approximately 75% wider than that expected from range straggling alone. The combination of these tested techniques makes possible high mass resolution in the neighborhood of iron.

The rotation of ripple pattern and the shape of the collision cascade in ion sputtered thin metal films

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mishra, P.; Ghose, D.

The sputter ripple formation in polycrystalline metal thin films of Al, Co, Cu, and Ag has been studied by 16.7 keV Ar{sup +} and O{sub 2}{sup +} ion bombardment as a function of angle of ion incidence. The experimental results show the existence of a critical angle of ion incidence ({theta}{sub c}) beyond which the ripples of wave vectors perpendicular to the projected ion beam direction appear. Monte Carlo simulation (SRIM) is carried out to calculate the depth, longitudinal and lateral straggling widths of energy deposition as these values are crucial in determining the critical angle {theta}{sub c}. It ismore » found that the radial energy distribution of the damage cascade has the maximum slightly away from the ion path in contradiction to the Gaussian distribution and the distribution is better characterized by an exponential function. The lower values of lateral straggling widths as those extracted from the measured critical angles using the Bradley and Harper theory indicate a highly anisotropic deposited-energy distribution.« less

Low-Energy Proton Testing Methodology

NASA Technical Reports Server (NTRS)

Pellish, Jonathan A.; Marshall, Paul W.; Heidel, David F.; Schwank, James R.; Shaneyfelt, Marty R.; Xapsos, M.A.; Ladbury, Raymond L.; LaBel, Kenneth A.; Berg, Melanie; Kim, Hak S.;

Active Radiation Detectors for Use in Space Beyond Low Earth Orbit: Spatial and Energy Resolution Requirements and Methods for Heavy Ion Charge Classification

NASA Astrophysics Data System (ADS)

McBeth, Rafe A.

Space radiation exposure to astronauts will need to be carefully monitored on future missions beyond low earth orbit. NASA has proposed an updated radiation risk framework that takes into account a significant amount of radiobiological and heavy ion track structure information. These models require active radiation detection systems to measure the energy and ion charge Z. However, current radiation detection systems cannot meet these demands. The aim of this study was to investigate several topics that will help next generation detection systems meet the NASA objectives. Specifically, this work investigates the required spatial resolution to avoid coincident events in a detector, the effects of energy straggling and conversion of dose from silicon to water, and methods for ion identification (Z) using machine learning. The main results of this dissertation are as follows: 1. Spatial resolution on the order of 0.1 cm is required for active space radiation detectors to have high confidence in identifying individual particles, i.e., to eliminate coincident events. 2. Energy resolution of a detector system will be limited by energy straggling effects and the conversion of dose in silicon to dose in biological tissue (water). 3. Machine learning methods show strong promise for identification of ion charge (Z) with simple detector designs.

Energy straggling eliminated as a limitation to charge resolution of transmission detectors. [used for particle identification

NASA Technical Reports Server (NTRS)

Tarle, G.; Ahlen, S. P.; Price, P. B.

1981-01-01

It is pointed out that detectors of the energy loss of penetrating charged particles are widely used for particle identification. These measurements are hampered, however, by fluctuations in the amount of energy deposited within the detector. It is shown that this limitation can be overcome with a new nuclear track detector, CR-39(DOP), and that the charge resolution of this detector exceeds that of any other, including semiconductor diodes.

Comparison of Cf-252 thin-film sources prepared by evaporation or self-transfer

DOE PAGES

Algutifan, Noor J.; Sherman, Steven R.; Alexander, Charles W.

2014-11-29

Californium-252 (Z = 98) is valued as a potent neutron source due to its spontaneous fission decay path. Thin film sources containing Cf-252 were prepared by two techniques: evaporation and self-transfer. The sources were analyzed by alpha and gamma spectroscopy. Results indicate that self-transfer sources exhibit less alpha energy straggling and energy loss than evaporative sources. Fission fragments may also self-transfer, and sources made by self-transfer may need some decay time to reach radioactive equilibrium.

Study of Gaussian Doped Double Gate JunctionLess (GD-DG-JL) transistor including source drain depletion length: Model for sub-threshold behavior

NASA Astrophysics Data System (ADS)

Kumari, Vandana; Kumar, Ayush; Saxena, Manoj; Gupta, Mridula

2018-01-01

The sub-threshold model formulation of Gaussian Doped Double Gate JunctionLess (GD-DG-JL) FET including source/drain depletion length is reported in the present work under the assumption that the ungated regions are fully depleted. To provide deeper insight into the device performance, the impact of gaussian straggle, channel length, oxide and channel thickness and high-k gate dielectric has been studied using extensive TCAD device simulation.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zieb, Kristofer James Ekhart; Hughes, Henry Grady III; Xu, X. George

The release of version 6.2 of the MCNP6 radiation transport code is imminent. To complement the newest release, a summary of the heavy charged particle physics models used in the 1 MeV to 1 GeV energy regime is presented. Several changes have been introduced into the charged particle physics models since the merger of the MCNP5 and MCNPX codes into MCNP6. Here, this article discusses the default models used in MCNP6 for continuous energy loss, energy straggling, and angular scattering of heavy charged particles. Explanations of the physics models’ theories are included as well.

Review of heavy charged particle transport in MCNP6.2

NASA Astrophysics Data System (ADS)

Zieb, K.; Hughes, H. G.; James, M. R.; Xu, X. G.

2018-04-01

The release of version 6.2 of the MCNP6 radiation transport code is imminent. To complement the newest release, a summary of the heavy charged particle physics models used in the 1 MeV to 1 GeV energy regime is presented. Several changes have been introduced into the charged particle physics models since the merger of the MCNP5 and MCNPX codes into MCNP6. This paper discusses the default models used in MCNP6 for continuous energy loss, energy straggling, and angular scattering of heavy charged particles. Explanations of the physics models' theories are included as well.

Changes in the frequency distribution of energy deposited in short pathlengths as a function of energy degradation of the primary beam.

NASA Technical Reports Server (NTRS)

Baily, N. A.; Steigerwalt, J. E.; Hilbert, J. W.

1972-01-01

The frequency distributions of event size in the deposition of energy over small pathlengths have been measured after penetration of 44.3 MeV protons through various thicknesses of tissue-equivalent material. Results show that particle energy straggling of an initially monoenergetic proton beam after passage through an absorber causes the frequency distributions of energy deposited in short pathlengths of low atomic number materials to remain broad. In all cases investigated, the ratio of the most probable to the average energy losses has been significantly less than unity.

Programming an Experiment Control System

NASA Technical Reports Server (NTRS)

Lange, Stuart

2004-01-01

As NASA develops plans for more and more ambitious missions into space, it is the job of NASA's researchers to develop the technologies that will make those planed missions feasible. One such technology is energy conversion. Energy is all around us; it is in the light that we see in the chemical bonds that hold compounds together, and in mass itself.Energy is the fundamental building block of our universe, yet it has always been straggle for humans to convert this energy into useable forms, like electricity. For space-based applications, NASA requires efficient energy conversion method that require little or no fuel.

Review of Heavy Charged Particle Transport in MCNP6.2

DOE PAGES

Zieb, Kristofer James Ekhart; Hughes, Henry Grady III; Xu, X. George; ...

2018-01-05

The release of version 6.2 of the MCNP6 radiation transport code is imminent. To complement the newest release, a summary of the heavy charged particle physics models used in the 1 MeV to 1 GeV energy regime is presented. Several changes have been introduced into the charged particle physics models since the merger of the MCNP5 and MCNPX codes into MCNP6. Here, this article discusses the default models used in MCNP6 for continuous energy loss, energy straggling, and angular scattering of heavy charged particles. Explanations of the physics models’ theories are included as well.

Monte Carlo Simulation of Spacecraft Particle Detectors to Assess the True Human Risk

NASA Technical Reports Server (NTRS)

O'Neill, Patrick M.

2002-01-01

Particle detectors (DOSTEL, CPDS, and TEPC) measure the energy deposition spectrum inside earth orbiting - manned spacecraft (shuttle, space station). These instruments attempt to emulate the deposition of energy in human tissue to evaluate the health risk. However, the measurements are often difficult to relate to tissue equivalent because nuclear fragmentation (internuclear cascade/evaporation), energy-loss straggling, heavy ions, spacecraft shielding and detector geometry/orientation, and coincidence thresholds significantly affect the measured spectrum. 'A le have developed a high fidelity Monte Carlo model addressing each of these effects that significantly improves interpretation of these instruments and the resulting assessment of radiation risk to humans.

Analytical and Numerical Solutions of Generalized Fokker-Planck Equations - Final Report

DOE Office of Scientific and Technical Information (OSTI.GOV)

Prinja, Anil K.

The overall goal of this project was to develop advanced theoretical and numerical techniques to quantitatively describe the spreading of a collimated beam of charged particles in space, in angle, and in energy, as a result of small deflection, small energy transfer Coulomb collisions with the target nuclei and electrons. Such beams arise in several applications of great interest in nuclear engineering, and include electron and ion radiotherapy, ion beam modification of materials, accelerator transmutation of waste, and accelerator production of tritium, to name some important candidates. These applications present unique and difficult modeling challenges, but from the outset aremore » amenable to the language of ''transport theory'', which is very familiar to nuclear engineers and considerably less-so to physicists and material scientists. Thus, our approach has been to adopt a fundamental description based on transport equations, but the forward peakedness associated with charged particle interactions precludes a direct application of solution methods developed for neutral particle transport. Unique problem formulations and solution techniques are necessary to describe the transport and interaction of charged particles. In particular, we have developed the Generalized Fokker-Planck (GFP) approach to describe the angular and radial spreading of a collimated beam and a renormalized transport model to describe the energy-loss straggling of an initially monoenergetic distribution. Both analytic and numerical solutions have been investigated and in particular novel finite element numerical methods have been developed. In the first phase of the project, asymptotic methods were used to develop closed form solutions to the GFP equation for different orders of expansion, and was described in a previous progress report. In this final report we present a detailed description of (i) a novel energy straggling model based on a Fokker-Planck approximation but which is adapted for a multigroup transport setting, and (ii) two unique families of discontinuous finite element schemes, one linear and the other nonlinear.« less

Advances in Monte-Carlo code TRIPOLI-4®'s treatment of the electromagnetic cascade

NASA Astrophysics Data System (ADS)

Mancusi, Davide; Bonin, Alice; Hugot, François-Xavier; Malouch, Fadhel

2018-01-01

TRIPOLI-4® is a Monte-Carlo particle-transport code developed at CEA-Saclay (France) that is employed in the domains of nuclear-reactor physics, criticality-safety, shielding/radiation protection and nuclear instrumentation. The goal of this paper is to report on current developments, validation and verification made in TRIPOLI-4 in the electron/positron/photon sector. The new capabilities and improvements concern refinements to the electron transport algorithm, the introduction of a charge-deposition score, the new thick-target bremsstrahlung option, the upgrade of the bremsstrahlung model and the improvement of electron angular straggling at low energy. The importance of each of the developments above is illustrated by comparisons with calculations performed with other codes and with experimental data.

Outbreak patterns of the novel avian influenza (H7N9)

NASA Astrophysics Data System (ADS)

Pan, Ya-Nan; Lou, Jing-Jing; Han, Xiao-Pu

2014-05-01

The attack of novel avian influenza (H7N9) in East China caused a serious health crisis and public panic. In this paper, we empirically analyze the onset patterns of human cases of the novel avian influenza and observe several spatial and temporal properties that are similar to other infectious diseases. More specifically, using the empirical analysis and modeling studies, we find that the spatio-temporal network that connects the cities with human cases along the order of outbreak timing emerges two-regime-power-law edge-length distribution, indicating the picture that several islands with higher and heterogeneous risk straggle in East China. The proposed method is applicable to the analysis of the spreading situation in the early stage of disease outbreak using quite limited dataset.

Invited article: Characterization of background sources in space-based time-of-flight mass spectrometers.

PubMed

Gilbert, J A; Gershman, D J; Gloeckler, G; Lundgren, R A; Zurbuchen, T H; Orlando, T M; McLain, J; von Steiger, R

2014-09-01

For instruments that use time-of-flight techniques to measure space plasma, there are common sources of background signals that evidence themselves in the data. The background from these sources may increase the complexity of data analysis and reduce the signal-to-noise response of the instrument, thereby diminishing the science value or usefulness of the data. This paper reviews several sources of background commonly found in time-of-flight mass spectrometers and illustrates their effect in actual data using examples from ACE-SWICS and MESSENGER-FIPS. Sources include penetrating particles and radiation, UV photons, energy straggling and angular scattering, electron stimulated desorption of ions, ion-induced electron emission, accidental coincidence events, and noise signatures from instrument electronics. Data signatures of these sources are shown, as well as mitigation strategies and design considerations for future instruments.

Comparisons between GRNTRN simulations and beam measurements of proton lateral broadening distributions

NASA Astrophysics Data System (ADS)

Mertens, Christopher; Moyers, Michael; Walker, Steven; Tweed, John

Recent developments in NASA's High Charge and Energy Transport (HZETRN) code have included lateral broadening of primary ion beams due to small-angle multiple Coulomb scattering, and coupling of the ion-nuclear scattering interactions with energy loss and straggling. The new version of HZETRN based on Green function methods, GRNTRN, is suitable for modeling transport with both space environment and laboratory boundary conditions. Multiple scattering processes are a necessary extension to GRNTRN in order to accurately model ion beam experiments, to simulate the physical and biological-effective radiation dose, and to develop new methods and strategies for light ion radiation therapy. In this paper we compare GRNTRN simulations of proton lateral scattering distributions with beam measurements taken at Loma Linda Medical University. The simulated and measured lateral proton distributions will be compared for a 250 MeV proton beam on aluminum, polyethylene, polystyrene, bone, iron, and lead target materials.

Description of a Generalized Analytical Model for the Micro-dosimeter Response

NASA Technical Reports Server (NTRS)

Badavi, Francis F.; Stewart-Sloan, Charlotte R.; Xapsos, Michael A.; Shinn, Judy L.; Wilson, John W.; Hunter, Abigail

2007-01-01

An analytical prediction capability for space radiation in Low Earth Orbit (LEO), correlated with the Space Transportation System (STS) Shuttle Tissue Equivalent Proportional Counter (TEPC) measurements, is presented. The model takes into consideration the energy loss straggling and chord length distribution of the TEPC detector, and is capable of predicting energy deposition fluctuations in a micro-volume by incoming ions through both direct and indirect ionic events. The charged particle transport calculations correlated with STS 56, 51, 110 and 114 flights are accomplished by utilizing the most recent version (2005) of the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport WZETRN), which has been extensively validated with laboratory beam measurements and available space flight data. The agreement between the TEPC model prediction (response function) and the TEPC measured differential and integral spectra in lineal energy (y) domain is promising.

Evaluation of the new electron-transport algorithm in MCNP6.1 for the simulation of dose point kernel in water

NASA Astrophysics Data System (ADS)

Antoni, Rodolphe; Bourgois, Laurent

2017-12-01

In this work, the calculation of specific dose distribution in water is evaluated in MCNP6.1 with the regular condensed history algorithm the "detailed electron energy-loss straggling logic" and the new electrons transport algorithm proposed the "single event algorithm". Dose Point Kernel (DPK) is calculated with monoenergetic electrons of 50, 100, 500, 1000 and 3000 keV for different scoring cells dimensions. A comparison between MCNP6 results and well-validated codes for electron-dosimetry, i.e., EGSnrc or Penelope, is performed. When the detailed electron energy-loss straggling logic is used with default setting (down to the cut-off energy 1 keV), we infer that the depth of the dose peak increases with decreasing thickness of the scoring cell, largely due to combined step-size and boundary crossing artifacts. This finding is less prominent for 500 keV, 1 MeV and 3 MeV dose profile. With an appropriate number of sub-steps (ESTEP value in MCNP6), the dose-peak shift is almost complete absent to 50 keV and 100 keV electrons. However, the dose-peak is more prominent compared to EGSnrc and the absorbed dose tends to be underestimated at greater depths, meaning that boundaries crossing artifact are still occurring while step-size artifacts are greatly reduced. When the single-event mode is used for the whole transport, we observe the good agreement of reference and calculated profile for 50 and 100 keV electrons. Remaining artifacts are fully vanished, showing a possible transport treatment for energies less than a hundred of keV and accordance with reference for whatever scoring cell dimension, even if the single event method initially intended to support electron transport at energies below 1 keV. Conversely, results for 500 keV, 1 MeV and 3 MeV undergo a dramatic discrepancy with reference curves. These poor results and so the current unreliability of the method is for a part due to inappropriate elastic cross section treatment from the ENDF/B-VI.8 library in those energy ranges. Accordingly, special care has to be taken in setting choice for calculating electron dose distribution with MCNP6, in particular with regards to dosimetry or nuclear medicine applications.

Efficiency improvement in proton dose calculations with an equivalent restricted stopping power formalism

NASA Astrophysics Data System (ADS)

Maneval, Daniel; Bouchard, Hugo; Ozell, Benoît; Després, Philippe

2018-01-01

The equivalent restricted stopping power formalism is introduced for proton mean energy loss calculations under the continuous slowing down approximation. The objective is the acceleration of Monte Carlo dose calculations by allowing larger steps while preserving accuracy. The fractional energy loss per step length ɛ was obtained with a secant method and a Gauss-Kronrod quadrature estimation of the integral equation relating the mean energy loss to the step length. The midpoint rule of the Newton-Cotes formulae was then used to solve this equation, allowing the creation of a lookup table linking ɛ to the equivalent restricted stopping power L eq, used here as a key physical quantity. The mean energy loss for any step length was simply defined as the product of the step length with L eq. Proton inelastic collisions with electrons were added to GPUMCD, a GPU-based Monte Carlo dose calculation code. The proton continuous slowing-down was modelled with the L eq formalism. GPUMCD was compared to Geant4 in a validation study where ionization processes alone were activated and a voxelized geometry was used. The energy straggling was first switched off to validate the L eq formalism alone. Dose differences between Geant4 and GPUMCD were smaller than 0.31% for the L eq formalism. The mean error and the standard deviation were below 0.035% and 0.038% respectively. 99.4 to 100% of GPUMCD dose points were consistent with a 0.3% dose tolerance. GPUMCD 80% falloff positions (R80 ) matched Geant’s R80 within 1 μm. With the energy straggling, dose differences were below 2.7% in the Bragg peak falloff and smaller than 0.83% elsewhere. The R80 positions matched within 100 μm. The overall computation times to transport one million protons with GPUMCD were 31-173 ms. Under similar conditions, Geant4 computation times were 1.4-20 h. The L eq formalism led to an intrinsic efficiency gain factor ranging between 30-630, increasing with the prescribed accuracy of simulations. The L eq formalism allows larger steps leading to a O(constant) algorithmic time complexity. It significantly accelerates Monte Carlo proton transport while preserving accuracy. It therefore constitutes a promising variance reduction technique for computing proton dose distributions in a clinical context.

Method of measuring field funneling and range straggling in semiconductor charge-collecting junctions

NASA Technical Reports Server (NTRS)

Zoutendyk, John A. (Inventor); Malone, Carl J. (Inventor)

1987-01-01

Electric-field funneling length is measured while irradiating a semiconductor charge-collecting junction with electron-hole-pair generating charged particles at a first junction bias voltage. The bias voltage is then reduced to a second level in order to reduce the depth of the depletion region such that the total charge can no longer be collected by drift and measured in the energy band previously displayed in the multichannel analyzer. This is representative of the maximum electric field funnelling length which may be calculated by measuring the difference at the second bias voltage level of the depletion width and the ion penetration range. The bias voltage is further lowered to a third level at which the particles are collected over a spread of energy levels while at least some of the particles are still collected at the selected energy level. From this the different depths of penetration of the particles are determined while additional effects due to diffusion are minimized.

Method of measuring field funneling and range straggling in semiconductor charge-collecting junctions

NASA Technical Reports Server (NTRS)

Zoutendyk, J. A. (Inventor)

1985-01-01

Electric-field funneling length is measured while irradiating a semiconductor charge-collecting junction with electron-hole-pair generating charged particles at a first junction bias voltage. The bias voltage is then reduced to a second level in order to reduce the depth of the depletion region such that the total charge can no longer be collected by drift and measured in the energy band previously displayed in the multichannel analyzer. This is representative of the maximum electric field funneling length which may be calculated by measuring the difference at the second bias voltage level of the depletion width and the ion penetration range. The bias voltage is further lowered to a third level at which the particles are collected over a spread of energy levels while at least some of the particles are still collected at the selected energy level. From this the different depths of penetration of the particles are determined while additional effects due to diffusion are minimized.

A -100 kV Power Supply for Ion Acceleration in Space-based Mass Spectrometers

NASA Astrophysics Data System (ADS)

Gilbert, J. A.; Zurbuchen, T.; Battel, S.

2017-12-01

High voltage power supplies are used in many space-based time-of-flight (TOF) mass spectrometer designs to accelerate incoming ions and increase the probability of their measurement and proper identification. Ions are accelerated in proportion to their charge state, so singly charged ions such as pickup ions are accelerated less than their multiple-charge state solar wind counterparts. This lack of acceleration results in pickup ion measurements with lower resolution and without determinations of absolute energy. Acceleration reduces the effects of angular scattering and energy straggling when ions pass through thin membranes such as carbon foils, and it brings ion energies above the detection threshold of traditional solid state detectors. We have developed a power supply capable of operating at -100 kV for ion acceleration while also delivering up to 10 W of power for the operation of a floating TOF system. We also show results of benchtop calibration and ion beam tests to demonstrate the functionality and success of this approach.

Extended wave-packet model to calculate energy-loss moments of protons in matter

NASA Astrophysics Data System (ADS)

Archubi, C. D.; Arista, N. R.

2017-12-01

In this work we introduce modifications to the wave-packet method proposed by Kaneko to calculate the energy-loss moments of a projectile traversing a target which is represented in terms of Gaussian functions for the momentum distributions of electrons in the atomic shells. These modifications are introduced using the Levine and Louie technique to take into account the energy gaps corresponding to the different atomic levels of the target. We use the extended wave-packet model to evaluate the stopping power, the energy straggling, the inverse mean free path, and the ionization cross sections for protons in several targets, obtaining good agreements for all these quantities on an extensive energy range that covers low-, intermediate-, and high-energy regions. The extended wave-packet model proposed here provides a method to calculate in a very straightforward way all the significant terms of the inelastic interaction of light ions with any element of the periodic table.

Recent Progress in the Development of a Multi-Layer Green's Function Code for Ion Beam Transport

NASA Technical Reports Server (NTRS)

Tweed, John; Walker, Steven A.; Wilson, John W.; Tripathi, Ram K.

2008-01-01

To meet the challenge of future deep space programs, an accurate and efficient engineering code for analyzing the shielding requirements against high-energy galactic heavy radiation is needed. To address this need, a new Green's function code capable of simulating high charge and energy ions with either laboratory or space boundary conditions is currently under development. The computational model consists of combinations of physical perturbation expansions based on the scales of atomic interaction, multiple scattering, and nuclear reactive processes with use of the Neumann-asymptotic expansions with non-perturbative corrections. The code contains energy loss due to straggling, nuclear attenuation, nuclear fragmentation with energy dispersion and downshifts. Previous reports show that the new code accurately models the transport of ion beams through a single slab of material. Current research efforts are focused on enabling the code to handle multiple layers of material and the present paper reports on progress made towards that end.

A study of trends and techniques for space base electronics

NASA Technical Reports Server (NTRS)

Trotter, J. D.; Wade, T. E.; Gassaway, J. D.; Mahmood, Q.

1978-01-01

A sputtering system was developed to deposit aluminum and aluminum alloys by the dc sputtering technique. This system is designed for a high level of cleanliness and for monitoring the deposition parameters during film preparation. This system is now ready for studying the deposition and annealing parameters upon double-level metal preparation. A technique recently applied for semiconductor analysis, the finite element method, was studied for use in the computer modeling of two dimensional MOS transistor structures. It was concluded that the method has not been sufficiently well developed for confident use at this time. An algorithm was developed for confident use at this time. An algorithm was developed for implementing a computer study which is based upon the finite difference method. The program which was developed was modified and used to calculate redistribution data for boron and phosphorous which had been predeposited by ion implantation with range and straggle conditions. Data were generated for 111 oriented SOS films with redistribution in N2, dry O2 and steam ambients.

Freestanding ultrathin single-crystalline SiC substrate by MeV H ion-slicing

NASA Astrophysics Data System (ADS)

Jia, Qi; Huang, Kai; You, Tiangui; Yi, Ailun; Lin, Jiajie; Zhang, Shibin; Zhou, Min; Zhang, Bin; Zhang, Bo; Yu, Wenjie; Ou, Xin; Wang, Xi

2018-05-01

SiC is a widely used wide-bandgap semiconductor, and the freestanding ultrathin single-crystalline SiC substrate provides the material platform for advanced devices. Here, we demonstrate the fabrication of a freestanding ultrathin single-crystalline SiC substrate with a thickness of 22 μm by ion slicing using 1.6 MeV H ion implantation. The ion-slicing process performed in the MeV energy range was compared to the conventional case using low-energy H ion implantation in the keV energy range. The blistering behavior of the implanted SiC surface layer depends on both the implantation temperature and the annealing temperature. Due to the different straggling parameter for two implant energies, the distribution of implantation-induced damage is significantly different. The impact of implantation temperature on the high-energy and low-energy slicing was opposite, and the ion-slicing SiC in the MeV range initiates at a much higher temperature.

An analytical model for calculating microdosimetric distributions from heavy ions in nanometer site targets.

PubMed

Czopyk, L; Olko, P

2006-01-01

The analytical model of Xapsos used for calculating microdosimetric spectra is based on the observation that straggling of energy loss can be approximated by a log-normal distribution of energy deposition. The model was applied to calculate microdosimetric spectra in spherical targets of nanometer dimensions from heavy ions at energies between 0.3 and 500 MeV amu(-1). We recalculated the originally assumed 1/E(2) initial delta electrons spectrum by applying the Continuous Slowing Down Approximation for secondary electrons. We also modified the energy deposition from electrons of energy below 100 keV, taking into account the effective path length of the scattered electrons. Results of our model calculations agree favourably with results of Monte Carlo track structure simulations using MOCA-14 for light ions (Z = 1-8) of energy ranging from E = 0.3 to 10.0 MeV amu(-1) as well as with results of Nikjoo for a wall-less proportional counter (Z = 18).

An Improved Green's Function for Ion Beam Transport

NASA Technical Reports Server (NTRS)

Tweed, J.; Wilson, J. W.; Tripathi, R. K.

2003-01-01

Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a Green's function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple solutions can be found for the high charge and energy (HZE) by ignoring nuclear energy downshifts and dispersion. Such solutions were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we have found global solutions including these effects providing a broader class of HZE ion solutions.

Calcium, phosphorus and protein levels as factors in the distribution of the pheasant

USGS Publications Warehouse

Dale, F.H.; DeWitt, J.B.

1958-01-01

Summary of work on pheasant nutrition conducted since 1949 at the Patuxent Research Refuge. Pheasant chicks fed experimental diets failed to develop normally on protein levels of 15 and 18%. With 22% protein they grew at a reduced rate as compared to those on 28%. Protein level of the reproductive diet was shown to be important; low production of eggs and young resulted from levels below 25%. Calcium was found to be even more critical than protein level for reproduction; birds on a winter diet that furnished 145 mg./kg. per day had poor reproductive success the following spring. About 600 mg./kg. of Ca per day was necessary in the reproduction diet. Birds on an intermediate level of Ca (about 0.5% of diet) showed evidence of cumulative deficiency. It was concluded that pheasants receiving levels of Ca no higher than 0.5% in nature might display 'straggling failure' such as has been observed in several midwestern areas.

Simulation of neutron production using MCNPX+MCUNED.

PubMed

Erhard, M; Sauvan, P; Nolte, R

2014-10-01

In standard MCNPX, the production of neutrons by ions cannot be modelled efficiently. The MCUNED patch applied to MCNPX 2.7.0 allows to model the production of neutrons by light ions down to energies of a few kiloelectron volts. This is crucial for the simulation of neutron reference fields. The influence of target properties, such as the diffusion of reactive isotopes into the target backing or the effect of energy and angular straggling, can be studied efficiently. In this work, MCNPX/MCUNED calculations are compared with results obtained with the TARGET code for simulating neutron production. Furthermore, MCUNED incorporates more effective variance reduction techniques and a coincidence counting tally. This allows the simulation of a TCAP experiment being developed at PTB. In this experiment, 14.7-MeV neutrons will be produced by the reaction T(d,n)(4)He. The neutron fluence is determined by counting alpha particles, independently of the reaction cross section. © The Author 2013. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.

Computer simulation program for medium-energy ion scattering and Rutherford backscattering spectrometry

NASA Astrophysics Data System (ADS)

Nishimura, Tomoaki

2016-03-01

A computer simulation program for ion scattering and its graphical user interface (MEISwin) has been developed. Using this program, researchers have analyzed medium-energy ion scattering and Rutherford backscattering spectrometry at Ritsumeikan University since 1998, and at Rutgers University since 2007. The main features of the program are as follows: (1) stopping power can be chosen from five datasets spanning several decades (from 1977 to 2011), (2) straggling can be chosen from two datasets, (3) spectral shape can be selected as Gaussian or exponentially modified Gaussian, (4) scattering cross sections can be selected as Coulomb or screened, (5) simulations adopt the resonant elastic scattering cross section of 16O(4He, 4He)16O, (6) pileup simulation for RBS spectra is supported, (7) natural and specific isotope abundances are supported, and (8) the charge fraction can be chosen from three patterns (fixed, energy-dependent, and ion fraction with charge-exchange parameters for medium-energy ion scattering). This study demonstrates and discusses the simulations and their results.

Trends and techniques for space base electronics. [mathematical models, ion implantation, and semiconductors

NASA Technical Reports Server (NTRS)

Gassaway, J. D.; Mahmood, Q.; Trotter, J. D.

1978-01-01

A system was developed for depositing aluminum and aluminum alloys by the D.C. sputtering technique. This system which was designed for a high level of cleanliness and ion monitoring the deposition parameters during film preparation is ready for studying the deposition and annealing parameters upon double level metal preparation. The finite element method was studied for use in the computer modeling of two dimensional MOS transistor structures. An algorithm was developed for implementing a computer study which is based upon the finite difference method. The program was modified and used to calculate redistribution data for boron and phosphorous which had been predeposited by ion implantation with range and straggle conditions typical of those used at MSFC. Data were generated for 111 oriented SOS films with redistribution in N2, dry O2 and steam ambients. Data are given showing both two dimensional effects and the evolution of the junction depth, sheet resistance and integrated dose with redistribution time.

Hybrid Methods for Muon Accelerator Simulations with Ionization Cooling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kunz, Josiah; Snopok, Pavel; Berz, Martin

Muon ionization cooling involves passing particles through solid or liquid absorbers. Careful simulations are required to design muon cooling channels. New features have been developed for inclusion in the transfer map code COSY Infinity to follow the distribution of charged particles through matter. To study the passage of muons through material, the transfer map approach alone is not sufficient. The interplay of beam optics and atomic processes must be studied by a hybrid transfer map--Monte-Carlo approach in which transfer map methods describe the deterministic behavior of the particles, and Monte-Carlo methods are used to provide corrections accounting for the stochasticmore » nature of scattering and straggling of particles. The advantage of the new approach is that the vast majority of the dynamics are represented by fast application of the high-order transfer map of an entire element and accumulated stochastic effects. The gains in speed are expected to simplify the optimization of cooling channels which is usually computationally demanding. Progress on the development of the required algorithms and their application to modeling muon ionization cooling channels is reported.« less

A fast and complete GEANT4 and ROOT Object-Oriented Toolkit: GROOT

NASA Astrophysics Data System (ADS)

Lattuada, D.; Balabanski, D. L.; Chesnevskaya, S.; Costa, M.; Crucillà, V.; Guardo, G. L.; La Cognata, M.; Matei, C.; Pizzone, R. G.; Romano, S.; Spitaleri, C.; Tumino, A.; Xu, Y.

2018-01-01

Present and future gamma-beam facilities represent a great opportunity to validate and evaluate the cross-sections of many photonuclear reactions at near-threshold energies. Monte Carlo (MC) simulations are very important to evaluate the reaction rates and to maximize the detection efficiency but, unfortunately, they can be very cputime-consuming and in some cases very hard to reproduce, especially when exploring near-threshold cross-section. We developed a software that makes use of the validated tracking GEANT4 libraries and the n-body event generator of ROOT in order to provide a fast, realiable and complete MC tool to be used for nuclear physics experiments. This tool is indeed intended to be used for photonuclear reactions at γ-beam facilities with ELISSA (ELI Silicon Strip Array), a new detector array under development at the Extreme Light Infrastructure - Nuclear Physics (ELI-NP). We discuss the results of MC simulations performed to evaluate the effects of the electromagnetic induced background, of the straggling due to the target thickness and of the resolution of the silicon detectors.

Optical properties of P ion implanted ZnO

NASA Astrophysics Data System (ADS)

Pong, Bao-Jen; Chou, Bo-Wei; Pan, Ching-Jen; Tsao, Fu-Chun; Chi, Gou-Chung

2006-02-01

Red and green emissions are observed from P ion implanted ZnO. Red emission at ~680 nm (1.82 eV) is associated with the donor-acceptor pair (DAP) transition, where the corresponding donor and acceptor are interstitial zinc (Zn i) and interstitial oxygen (O i), respectively. Green emission at ~ 516 nm (2.40 eV) is associated with the transition between the conduction band and antisite oxygen (O Zn). Green emission at ~516nm (2.403 eV) was observed for ZnO annealed at 800 oC under ambient oxygen, whereas, it was not visible when it was annealed in ambient nitrogen. Hence, the green emission is most likely not related to oxygen vacancies on ZnO sample, which might be related to the cleanliness of ZnO surface, a detailed study is in progress. The observed micro-strain is larger for N ion implanted ZnO than that for P ion implanted ZnO. It is attributed to the larger straggle of N ion implanted ZnO than that of P ion implanted ZnO. Similar phenomenon is also observed in Be and Mg ion implanted GaN.

An improved Green's function for ion beam transport

NASA Technical Reports Server (NTRS)

Tweed, J.; Wilson, J. W.; Tripathi, R. K.

2004-01-01

Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a Green's function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple solutions can be found for high charge and energy (HZE) ions by ignoring nuclear energy downshifts and dispersion. Such solutions were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we have found global solutions including these effects providing a broader class of HZE ion solutions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Silicon patterning using ion blistering and e-beam lithography

NASA Astrophysics Data System (ADS)

Giguere, A.; Terreault, B.; Beerens, J.; Aimez, V.; Beauvais, J.

2004-03-01

We explore the limits of silicon patterning using ion blistering in conjunction with e-beam lithography. In a first approach, we implanted 3.5E16 H/cm**2 at 5 keV through variable width (0.1-10 micron) e-beam written PMMA masks. The resist was then removed and the samples were rapid-thermal-annealed (RTA) up to 650 °C. In the wider trenches, round blisters with 800-900 nm diameter and 15 nm height and a few exfoliations are observed, which are similar to those observed on an unmasked surface. In submicron trenches (500-1000 nm), there is a transition in morphology created by the proximity to the border; the blisters are smaller and they are densely aligned along the trench direction ("pearl-string" pattern). No effect is observed in the lowest dimension trenches. The results are discussed in terms of stress/strain fields, defect configuration, and mask shadowing and charging effects. Ultimate pattern resolution will be limited by lateral straggling of the ions in and by the mechanics of lateral crack propagation.

Use of SRIM and Garfield with Geant4 for the characterization of a hybrid 10B/3He neutron detector

NASA Astrophysics Data System (ADS)

van der Ende, B. M.; Rand, E. T.; Erlandson, A.; Li, L.

2018-06-01

This paper describes a method for more complete neutron detector characterization using Geant4's Monte Carlo methods for characterizing overall detector response rate and Garfield interfaced with SRIM for the simulation of the detector's raw pulses, as applied to a hybrid 10B/3He detector. The Geant4 models characterizing the detector's interaction with a 252Cf point source and parallel beams of mono-energetic neutrons (assuming ISO 8529 reference energy values) compare and agree well with calibrated 252Cf measurements to within 6.4%. Validated Geant4 model outputs serve as input to Garfield+SRIM calculations to provide meaningful pulse height spectra. Modifications to Garfield for this work were necessary to account for simultaneous tracking of electrons resulting from proton and triton reaction products from a single 3He neutron capture event, and it was further necessary to interface Garfield with the energy loss, range, and straggling calculations provided by SRIM. Individual raw pulses generated by Garfield+SRIM are also observed to agree well with experimentally measured raw pulses from the detector.

Proton Lateral Broadening Distribution Comparisons Between GRNTRN, MCNPX, and Laboratory Beam Measurements

NASA Technical Reports Server (NTRS)

Mertens, Christopher J.; Moyers, Michael F.; Walker, Steven A.; Tweed, John

2010-01-01

Recent developments in NASA s deterministic High charge (Z) and Energy TRaNsport (HZETRN) code have included lateral broadening of primary ion beams due to small-angle multiple Coulomb scattering, and coupling of the ion-nuclear scattering interactions with energy loss and straggling. This new version of HZETRN is based on Green function methods, called GRNTRN, and is suitable for modeling transport with both space environment and laboratory boundary conditions. Multiple scattering processes are a necessary extension to GRNTRN in order to accurately model ion beam experiments, to simulate the physical and biological-effective radiation dose, and to develop new methods and strategies for light ion radiation therapy. In this paper we compare GRNTRN simulations of proton lateral broadening distributions with beam measurements taken at Loma Linda University Proton Therapy Facility. The simulated and measured lateral broadening distributions are compared for a 250 MeV proton beam on aluminum, polyethylene, polystyrene, bone substitute, iron, and lead target materials. The GRNTRN results are also compared to simulations from the Monte Carlo MCNPX code for the same projectile-target combinations described above.

Statistical exchange-coupling errors and the practicality of scalable silicon donor qubits

NASA Astrophysics Data System (ADS)

Song, Yang; Das Sarma, S.

2016-12-01

Recent experimental efforts have led to considerable interest in donor-based localized electron spins in Si as viable qubits for a scalable silicon quantum computer. With the use of isotopically purified 28Si and the realization of extremely long spin coherence time in single-donor electrons, the recent experimental focus is on two-coupled donors with the eventual goal of a scaled-up quantum circuit. Motivated by this development, we simulate the statistical distribution of the exchange coupling J between a pair of donors under realistic donor placement straggles, and quantify the errors relative to the intended J value. With J values in a broad range of donor-pair separation ( 5 <|R |<60 nm), we work out various cases systematically, for a target donor separation R0 along the [001], [110] and [111] Si crystallographic directions, with |R0|=10 ,20 or 30 nm and standard deviation σR=1 ,2 ,5 or 10 nm. Our extensive theoretical results demonstrate the great challenge for a prescribed J gate even with just a donor pair, a first step for any scalable Si-donor-based quantum computer.

An accurate model for the computation of the dose of protons in water.

PubMed

Embriaco, A; Bellinzona, V E; Fontana, A; Rotondi, A

2017-06-01

The accurate and fast calculation of the dose in proton radiation therapy is an essential ingredient for successful treatments. We propose a novel approach with a minimal number of parameters. The approach is based on the exact calculation of the electromagnetic part of the interaction, namely the Molière theory of the multiple Coulomb scattering for the transversal 1D projection and the Bethe-Bloch formula for the longitudinal stopping power profile, including a gaussian energy straggling. To this e.m. contribution the nuclear proton-nucleus interaction is added with a simple two-parameter model. Then, the non gaussian lateral profile is used to calculate the radial dose distribution with a method that assumes the cylindrical symmetry of the distribution. The results, obtained with a fast C++ based computational code called MONET (MOdel of ioN dosE for Therapy), are in very good agreement with the FLUKA MC code, within a few percent in the worst case. This study provides a new tool for fast dose calculation or verification, possibly for clinical use. Copyright © 2017 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.

Tailoring the structural and magnetic properties of masked CoPt thin films using ion implantation

NASA Astrophysics Data System (ADS)

Kumar, Durgesh; Gupta, Surbhi; Jin, Tianli; Nongjai, R.; Asokan, K.; Piramanayagam, S. N.

2018-05-01

The effects of ion implantations through a mask on the structural and magnetic properties of Co80Pt20 films were investigated. The mask was patterned using the self-assembly of diblock copolymers. For implantation, high (40 keV for 14N+ and 100 keV for 40Ar+) and low (7.5 keV for 14N+ and 4.5 keV for 40Ar+) energy 14N+ and 40Ar+ ions were used to modify the structural and magnetic properties of these films. X-ray diffraction and TRIM simulations were performed for understanding the structural changes due to ion implantations. These results revealed the intermixing of Co atoms in lower layers and lattice expansion in Co80Pt20 magnetic and Ru layers. A lateral straggling of Co caused an increase in the exchange coupling in the masked region. Depletion of Co atoms in Co80Pt20 layer caused a decrease in the anisotropy constant, which were further confirmed by the alternating gradient force magnetometer and magnetic force microscopy results. The magnetic force microscopy images showed an increase in domain width and domain wall width confirming the above-mentioned effects.

Comparison of various stopping gases for 3He-based position sensitive neutron detectors

NASA Astrophysics Data System (ADS)

Doumas, A.; Smith, G. C.

2012-05-01

A range of solid state, scintillator and gas based detectors are being developed for use at the next generation of high flux neutron facilities. Since gas detectors are expected to continue to play a key role in future specific thermal neutron experiments, a comparison of the performance characteristics of prospective stopping gases is beneficial. Gas detectors typically utilize the reaction 3He(n,p)t to detect thermal neutrons; the 3He gas is used in a mixture containing a particular stopping gas in order to maintain relatively short ranges for the proton and triton pair emitted from the n-3He reaction. Common stopping gases include hydrocarbons (e.g. propane), carbon tetrafluoride, and noble gases such as argon and xenon. For this study, we utilized the Monte Carlo simulation code "Stopping and Range of Ions in Matter" to analyze the expected behavior of argon, xenon, carbon dioxide, difluoroethane and octafluoropropane as stopping gases for thermal neutron detectors. We also compare these findings to our previously analyzed performance of propane, butane and carbon tetrafluoride. A discussion of these gases includes their behavior in terms of proton and triton range, ionization distribution and straggle.

The MONET code for the evaluation of the dose in hadrontherapy

NASA Astrophysics Data System (ADS)

Embriaco, A.

2018-01-01

The MONET is a code for the computation of the 3D dose distribution for protons in water. For the lateral profile, MONET is based on the Molière theory of multiple Coulomb scattering. To take into account also the nuclear interactions, we add to this theory a Cauchy-Lorentz function, where the two parameters are obtained by a fit to a FLUKA simulation. We have implemented the Papoulis algorithm for the passage from the projected to a 2D lateral distribution. For the longitudinal profile, we have implemented a new calculation of the energy loss that is in good agreement with simulations. The inclusion of the straggling is based on the convolution of energy loss with a Gaussian function. In order to complete the longitudinal profile, also the nuclear contributions are included using a linear parametrization. The total dose profile is calculated in a 3D mesh by evaluating at each depth the 2D lateral distributions and by scaling them at the value of the energy deposition. We have compared MONET with FLUKA in two cases: a single Gaussian beam and a lateral scan. In both cases, we have obtained a good agreement for different energies of protons in water.

Studies on heavy charged particle interaction, water equivalence and Monte Carlo simulation in some gel dosimeters, water, human tissues and water phantoms

NASA Astrophysics Data System (ADS)

Kurudirek, Murat

2015-09-01

Some gel dosimeters, water, human tissues and water phantoms were investigated with respect to their radiological properties in the energy region 10 keV-10 MeV. The effective atomic numbers (Zeff) and electron densities (Ne) for some heavy charged particles such as protons, He ions, B ions and C ions have been calculated for the first time for Fricke, MAGIC, MAGAT, PAGAT, PRESAGE, water, adipose tissue, muscle skeletal (ICRP), muscle striated (ICRU), plastic water, WT1 and RW3 using mass stopping powers from SRIM Monte Carlo software. The ranges and straggling were also calculated for the given materials. Two different set of mass stopping powers were used to calculate Zeff for comparison. The water equivalence of the given materials was also determined based on the results obtained. The Monte Carlo simulation of the charged particle transport was also done using SRIM code. The heavy ion distribution along with its parameters were shown for the given materials for different heavy ions. Also the energy loss and damage events in water when irradiated with 100 keV heavy ions were studied in detail.

Validation of a multi-layer Green's function code for ion beam transport

NASA Astrophysics Data System (ADS)

Walker, Steven; Tweed, John; Tripathi, Ram; Badavi, Francis F.; Miller, Jack; Zeitlin, Cary; Heilbronn, Lawrence

To meet the challenge of future deep space programs, an accurate and efficient engineering code for analyzing the shielding requirements against high-energy galactic heavy radiations is needed. In consequence, a new version of the HZETRN code capable of simulating high charge and energy (HZE) ions with either laboratory or space boundary conditions is currently under development. The new code, GRNTRN, is based on a Green's function approach to the solution of Boltzmann's transport equation and like its predecessor is deterministic in nature. The computational model consists of the lowest order asymptotic approximation followed by a Neumann series expansion with non-perturbative corrections. The physical description includes energy loss with straggling, nuclear attenuation, nuclear fragmentation with energy dispersion and down shift. Code validation in the laboratory environment is addressed by showing that GRNTRN accurately predicts energy loss spectra as measured by solid-state detectors in ion beam experiments with multi-layer targets. In order to validate the code with space boundary conditions, measured particle fluences are propagated through several thicknesses of shielding using both GRNTRN and the current version of HZETRN. The excellent agreement obtained indicates that GRNTRN accurately models the propagation of HZE ions in the space environment as well as in laboratory settings and also provides verification of the HZETRN propagator.

A fast - Monte Carlo toolkit on GPU for treatment plan dose recalculation in proton therapy

NASA Astrophysics Data System (ADS)

Senzacqua, M.; Schiavi, A.; Patera, V.; Pioli, S.; Battistoni, G.; Ciocca, M.; Mairani, A.; Magro, G.; Molinelli, S.

2017-10-01

In the context of the particle therapy a crucial role is played by Treatment Planning Systems (TPSs), tools aimed to compute and optimize the tratment plan. Nowadays one of the major issues related to the TPS in particle therapy is the large CPU time needed. We developed a software toolkit (FRED) for reducing dose recalculation time by exploiting Graphics Processing Units (GPU) hardware. Thanks to their high parallelization capability, GPUs significantly reduce the computation time, up to factor 100 respect to a standard CPU running software. The transport of proton beams in the patient is accurately described through Monte Carlo methods. Physical processes reproduced are: Multiple Coulomb Scattering, energy straggling and nuclear interactions of protons with the main nuclei composing the biological tissues. FRED toolkit does not rely on the water equivalent translation of tissues, but exploits the Computed Tomography anatomical information by reconstructing and simulating the atomic composition of each crossed tissue. FRED can be used as an efficient tool for dose recalculation, on the day of the treatment. In fact it can provide in about one minute on standard hardware the dose map obtained combining the treatment plan, earlier computed by the TPS, and the current patient anatomic arrangement.

Measurement of the 234U(n, f ) cross-section with quasi-monoenergetic beams in the keV and MeV range using a Micromegas detector assembly

NASA Astrophysics Data System (ADS)

Stamatopoulos, A.; Kanellakopoulos, A.; Kalamara, A.; Diakaki, M.; Tsinganis, A.; Kokkoris, M.; Michalopoulou, V.; Axiotis, M.; Lagoyiannis, A.; Vlastou, R.

2018-01-01

The 234U neutron-induced fission cross-section has been measured at incident neutron energies of 452, 550, 651 keV and 7.5, 8.7, 10 MeV using the 7Li ( p, n) and the 2H( d, n) reactions, respectively, relative to the 235U( n, f ) and 238U( n, f ) reference reactions. The measurement was performed at the neutron beam facility of the National Center for Scientific Research "Demokritos", using a set-up based on Micromegas detectors. The active mass of the actinide samples and the corresponding impurities were determined via α-spectroscopy using a surface barrier silicon detector. The neutron spectra intercepted by the actinide samples have been thoroughly studied by coupling the NeuSDesc and MCNP5 codes, taking into account the energy and angular straggling of the primary ion beams in the neutron source targets in addition to contributions from competing reactions ( e.g. deuteron break-up) and neutron scattering in the surrounding materials. Auxiliary Monte Carlo simulations were performed making combined use of the FLUKA and GEF codes, focusing particularly on the determination of the fission fragment detection efficiency. The developed methodology and the final results are presented.

Secondary ion mass spectroscopy study of Au trapping and migration in the Au-irradiated YBa2Cu3O7 - delta film

NASA Astrophysics Data System (ADS)

Li, Yupu; Kilner, J. A.; Liu, J. R.; Chu, W. K.; Wagner, G. A.; Somekh, R. E.

1996-05-01

The range data and migration of Au in YBa2Cu3O7-δ film were studied with implanted 197Au (1.5 MeV 5×1015 Au+/cm2) as a tracer. The film was a c-axis oriented film, ˜750 nm thick, deposited by high-pressure planar dc sputtering on <100> LaAlO3. Analysis by secondary ion mass spectroscopy shows that the as-implanted Au concentration distribution is essentially Gaussian-like and the depth (R̂p) of maximum Au concentration (˜1.2 wt %) is 201 nm. The projected range (R¯p) and (R̂p) are found to be in very good agreement with the simulated data by TRIM-95, whereas the measured ``straggle'' (ΔRp*) is about 20% larger than that by TRIM-95 simulation. It has also been found that the implanted 197Au starts to migrate within the film at a temperature between 650 and 700 °C, which is much higher than that for the implanted 2H (˜175 °C) and the implanted 18O (between 250 and 300 °C) in c-oriented YBa2Cu3O7-δ films.

Thick target bremsstrahlung spectra for 1.00-, 1.25-, and 1.40-Mev electrons

USGS Publications Warehouse

Miller, W.; Motz, J.W.; Cialella, C.

1954-01-01

The spectrum of radiation produced by 1.0-, 1.25-, and 1.40-Mev electrons incident on a thick tungsten target was measured at 0A????and 90A????with the incident beam by a method involving the magnetic analysis of Compton electrons. The effects of electron scattering and energy loss in the target preclude any simple interpretation of this data to yield a differential bremsstrahlung cross section. However, an estimate of the spectra to be expected at 0A????and 90A????was obtained by combining the Sauter expression for the bremsstrahlung cross section with the available information on electron scatter and energy loss in the target and backscatter from the target. The reliability of the estimate is limited because the Sauter formula was calculated by using the Born approximation, the electron scattering calculations are applicable to an infinite medium only, and the backscatter was estimated empirically from Bothe's experimental data which were obtained with lower energy electrons (370 kev). Furthermore electron energy straggling was neglected. Nevertheless, the predicted spectral shapes at 0A????and 90A????and the relative intensities at these two angles are in qualitative agreement with the measured values. The absolute magnitudes of the measured intensities at both angles are about a factor of two greater than the predicted values. ?? 1954 The American Physical Society.

Influence of surface topography on RBS measurements: case studies of (Cu/Fe/Pd) multilayers and FePdCu alloys nanopatterned by self-assembly

NASA Astrophysics Data System (ADS)

Krupinski, M.; Perzanowski, M.; Zabila, Y.; Zarzycki, A.; Marszałek, M.

2017-03-01

In this paper the influence of surface topography on Rutherford backscattering spectrometry (RBS) is discussed. (Cu/Fe/Pd) multilayers with total thickness of about 10 nm were deposited by physical vapor deposition on self-organized array of SiO2 nanoparticles with the size of 50 nm and 100 nm. As a reference, the multilayered systems were also prepared on flat substrates under the same conditions. After the deposition, morphology of the systems was studied by scanning electron microscopy (SEM), while chemical analysis was performed using Rutherford backscattering spectrometry. It was found that the RBS spectra and determined compositions for flat and patterned multilayers differ. The difference is discussed by taking into account the effect of additional inelastic scattering and energy straggling occurring due to developed topography of patterned systems. Then, the multilayers were annealed in 600 °C in order to obtain FePdCu alloy. The phenomenon of solid-state dewetting resulted in the formation of isolated alloy islands on the top of SiO2 nanoparticles. The SEM and RBS analysis were repeated showing correlation between the size distribution of obtained alloy islands and broadening of peaks appearing in RBS spectra. Invited talk at 8th International Workshop on Advanced Materials Science and Nanotechnology (IWAMSN2016), 8-12 November 2016, Ha Long City, Vietnam.

Alpha-spectrometry and fractal analysis of surface micro-images for characterisation of porous materials used in manufacture of targets for laser plasma experiments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Aushev, A A; Barinov, S P; Vasin, M G

2015-06-30

We present the results of employing the alpha-spectrometry method to determine the characteristics of porous materials used in targets for laser plasma experiments. It is shown that the energy spectrum of alpha-particles, after their passage through porous samples, allows one to determine the distribution of their path length in the foam skeleton. We describe the procedure of deriving such a distribution, excluding both the distribution broadening due to statistical nature of the alpha-particle interaction with an atomic structure (straggling) and hardware effects. The fractal analysis of micro-images is applied to the same porous surface samples that have been studied bymore » alpha-spectrometry. The fractal dimension and size distribution of the number of the foam skeleton grains are obtained. Using the data obtained, a distribution of the total foam skeleton thickness along a chosen direction is constructed. It roughly coincides with the path length distribution of alpha-particles within a range of larger path lengths. It is concluded that the combined use of the alpha-spectrometry method and fractal analysis of images will make it possible to determine the size distribution of foam skeleton grains (or pores). The results can be used as initial data in theoretical studies on propagation of the laser and X-ray radiation in specific porous samples. (laser plasma)« less

Transport Measurements on Si Nanostructures with Counted Sb Donors

NASA Astrophysics Data System (ADS)

Singh, Meenakshi; Bielejec, Edward; Garratt, Elias; Ten Eyck, Gregory; Bishop, Nathaniel; Wendt, Joel; Luhman, Dwight; Carroll, Malcolm; Lilly, Michael

2014-03-01

Donor based spin qubits are a promising platform for quantum computing. Single qubits using timed implant of donors have been demonstrated.1 Extending this to multiple qubits requires precise control over the placement and number of donors. Such control can be achieved by using a combination of low-energy heavy-ion implants (to reduce depth straggle), electron-beam lithography (to define position), focused ion beam (to localize implants to one lithographic site) and counting the number of implants with a single ion detector.2 We report transport measurements on MOS quantum dots implanted with 5, 10 and 20 Sb donors using the approach described above. A donor charge transition is identified by a charge offset in the transport characteristics. Correlation between the number of donors and the charge offsets is studied. These results are necessary first steps towards fabricating donor nanostructures for two qubit interactions. This work was performed, in part, at the Center for Integrated Nanotechnologies, a U.S. DOE Office of Basic Energy Sciences user facility. The work was supported by Sandia National Laboratories Directed Research and Development Program. Sandia National Laboratories is a multi-program laboratory operated by Sandia Corporation, a Lockheed-Martin Company, for the U. S. Department of Energy under Contract No. DE-AC04-94AL85000. 1J. J. Pla et al., Nature 496, 334 (2013) 2J. A. Seamons et al., APL 93, 043124 (2008).

Accelerator mass spectrometry of the heaviest long-lived radionuclides with a 3-MV tandem accelerator

NASA Astrophysics Data System (ADS)

Vockenhuber, Christof; Golser, Robin; Kutschera, Walter; Priller, Alfred; Steier, Peter; Winkler, Stephan; Liechtenstein, Vitaly

2002-12-01

A 3-MV pelletron tandem accelerator is the heart of the Vienna environmental research accelerator (VERA). The original design of the beam transport components allows the transport of ions of all elements, from the lightest to the heaviest. For light ions the suppression of neighboring masses was sufficient to measure isotopic ratios of {(14}) C/{(12}) C and {(26}) Al/{(27}) Al as low as 10{(-15}) and {(10}) Be/{(9}) Be down to 10{(-13}) . To suppress neighboring masses for the heaviest radionuclides in the energy range of 10-20 MeV, the resolution of VERA was increased both by improving the ion optics of existing elements at the injection side and by installing a new high-resolution electrostatic separator at the high-energy side. Interfering ions which pass all beam filters are identified with a Bragg-type ionization detector and a high-resolution time-of-flight system. Two ultra-thin diamond-like carbon (DLC) foils are used in the start and stop detector, which substantially reduces losses due to beam straggling. This improved set up enables us to measure even the heaviest long-lived radionuclides, where stable isobaric interferences are absent (e.g. {(236}) U and {(244}) Pu), down to environmental levels. Moreover, the advantage of a `small' and well manageable machine like VERA lies in its higher stability and reliability which allows to measure these heavy radionuclides more accurately, and also a large number of samples.

Xe- and U-tracks in apatite and muscovite near the etching threshold

NASA Astrophysics Data System (ADS)

Wauschkuhn, Bastian; Jonckheere, Raymond; Ratschbacher, Lothar

2015-01-01

Ion irradiation of a wedge-shaped Durango apatite backed by a mica detector allows investigating ion track ranges and etching properties at different points along the tracks. Transmission profiles obtained by irradiation with 2 × 106 cm-2 11.1 MeV/amu 132Xe and 2 × 106 cm-2 11.1 MeV/amu 238U parallel to the apatite c-axis correspond to ranges calculated with SRIM (Xe: 76.3 μm; U: 81.1 μm). However, the measured profiles show much greater etchable track-length variations than the calculated longitudinal straggles. The probable cause is that the length deficit exhibits significant variation from track to track. The measured length deficit in muscovite is in agreement with most existing data. In contrast, the length deficit in apatite appears to be close to zero, which is in conflict with all earlier estimates. This probably results from the etching properties of the apatite basal face, which permit surface-assisted sub-threshold etching of track sections in the nuclear stopping regime. These sections are not accessible from the opposite direction, i.e. by etching towards the endpoint of the tracks or in the direction of the ion beam. This conclusion is supported by the fact that linear dislocations are revealed in apatite basal faces and by the observation of imperfect etch pits that are separated from the etched ion track channel by a section that appears unetched under the microscope.

The center for production of single-photon emitters at the electrostatic-deflector line of the Tandem accelerator of LABEC (Florence)

NASA Astrophysics Data System (ADS)

Lagomarsino, Stefano; Sciortino, Silvio; Gelli, Nicla; Flatae, Assegid M.; Gorelli, Federico; Santoro, Mario; Chiari, Massimo; Czelusniac, Caroline; Massi, Mirko; Taccetti, Francesco; Agio, Mario; Giuntini, Lorenzo

2018-05-01

The line for the pulsed beam of the 3 MeV Tandetron accelerator at LABEC (Florence) has been upgraded for ion implantation experiments aiming at the fabrication of single-photon emitters in a solid-state matrix. A system based on Al attenuators has been calibrated in order to extend the energy range of the implanted ions from MeV down to the tens of keV. A new motorized XY stage has been installed in the implantation chamber for achieving ultra-fine control on the position of each implanted ion, allowing to reach the scale imposed by lateral straggling. A set-up for the activation of the implanted ions has been developed, based on an annealing furnace operating under controlled high-vacuum conditions. The first experiments have been performed with silicon ions implanted in diamond and the luminescent signal of the silicon-vacancy (SiV) center, peaked at 738 nm, has been observed for a wide range of implantation fluences (108 ÷ 1015 cm-2) and implantation depths (from a few nm to 2.4 μm). Studies on the efficiency of the annealing process have been performed and the activation yield has been measured to range from 1% to 3%. The implantation and annealing facility has thus been tuned for the production of SiV centers in diamond, but is in principle suitable for other ion species and solid-state matrices.

Hypolipidemic Activity of Chloroform Extract of Mimosa pudica Leaves

PubMed Central

Rajendran, Rekha; Krishnakumar, Ekambaram

2010-01-01

Mimosa pudica Lin., known as chue Mue, is a stout straggling prostrate shrubby plant, with spinous stipules and globose pinkish flower heads, and grows as weed in almost all parts of the country. It is traditionally used for its various properties and hence in the present study, chloroform extract of Mimosa pudica leaves has been screened for its hypolipidemic activity. Hypolipidemic activity is screened by inducing hyperlipidemia with the help of atherogenic diet in wistar albino rats and serum levels of various biochemical parameters such as total cholesterol, triglycerides, LDL, VLDL and HDL cholesterol were determined. Atherogenic index shows the measure of the athero-genic potential of the drugs. Chloroform extract showed significant (p < 0.05) hypolipidemic effect by lowering the serum levels of biochemical parameters such as significant reduction in the level of serum cholesterol, triglyceride, LDL, VLDL and increase in HDL level which was similar to the standard drug Atorvastatin. Chloroform extract exhibited significant atherogenic index and percentage protection against hyperlipidemia. These biochemical observations were in turn confirmed by histopathological examinations of aorta, liver and kidney sections and are comparable with the standard hypolipidemic drug Atorvastatin. Preliminary phytochemical analysis revealed the presence of phytoconstituents such as steroids, flavonoids, glycosides, alkaloids, phenolic compounds which is further confirmed by the thin layer chromatography, High Performance Thin Layer Chromatography (HPTLC). The overall experimental results suggests that the biologically active phytoconstituents such as flavonoids, glycosides alkaloids present in the chloroform extract of Mimosa pudica, may be responsible for the significant hypolipidemic activity and the results justify the use of Mimosa pudica as a significant hypolipidemic agent. PMID:23408779

Beam tracking with micromegas & wire chambers in secondary electron detection configuration

NASA Astrophysics Data System (ADS)

Voštinar, M.; Fernández, B.; Pancin, J.; Alvarez, M. A. G.; Chaminade, T.; Damoy, S.; Doré, D.; Drouart, A.; Druillole, F.; Frémont, G.; Kebbiri, M.; Materna, T.; Monmarthe, E.; Panebianco, S.; Papaevangelou, T.; Riallot, M.; Savajols, H.; Spitaels, C.

2013-12-01

The focal plane of S3 (Super Separator Spectrometer), a new experimental area of SPIRAL2 at GANIL, will be used for identification of nuclei, and requires the reconstruction of their trajectories and velocities by the Time Of Flight (TOF) method. Classical tracking detectors used in-beam would generate a lot of angular and energy straggling due to their thickness. One solution is the use of a SED (Secondary Electron Detection), which consists of a thin emissive foil in beam coupled to a low pressure gaseous detector out of the beam, for the detection of secondary electrons ejected from the foil. Moreover, this type of detector can be used for classical beam tracking at low energies, or for example at NFS (GANIL) for the FALSTAFF experiment for the reconstruction of fission fragments trajectories. Several low pressure gaseous detectors such as wire chambers and Micromegas have been constructed and tested since 2008. High counting rate capabilities and good time resolution obtained in previous tests motivated the construction of a new real-size 2D prototype wire chamber and a 2D bulk Micromegas at low pressure. For the first time, spatial resolution of the Micromegas at low pressure (below 20 mbar) in the SED configuration was measured. Different tests have been performed in order to characterize time and spatial properties of both prototypes, giving spatial resolution in the horizontal (X) direction of 0.90(0.02) mm FWHM for the real size prototype and 0.72(0.08) mm FWHM for Micromegas, and a time resolution of ~ 110(25) ps for the real size prototype.

SHIELD-HIT12A - a Monte Carlo particle transport program for ion therapy research

NASA Astrophysics Data System (ADS)

Bassler, N.; Hansen, D. C.; Lühr, A.; Thomsen, B.; Petersen, J. B.; Sobolevsky, N.

2014-03-01

Purpose: The Monte Carlo (MC) code SHIELD-HIT simulates the transport of ions through matter. Since SHIELD-HIT08 we added numerous features that improves speed, usability and underlying physics and thereby the user experience. The "-A" fork of SHIELD-HIT also aims to attach SHIELD-HIT to a heavy ion dose optimization algorithm to provide MC-optimized treatment plans that include radiobiology. Methods: SHIELD-HIT12A is written in FORTRAN and carefully retains platform independence. A powerful scoring engine is implemented scoring relevant quantities such as dose and track-average LET. It supports native formats compatible with the heavy ion treatment planning system TRiP. Stopping power files follow ICRU standard and are generated using the libdEdx library, which allows the user to choose from a multitude of stopping power tables. Results: SHIELD-HIT12A runs on Linux and Windows platforms. We experienced that new users quickly learn to use SHIELD-HIT12A and setup new geometries. Contrary to previous versions of SHIELD-HIT, the 12A distribution comes along with easy-to-use example files and an English manual. A new implementation of Vavilov straggling resulted in a massive reduction of computation time. Scheduled for later release are CT import and photon-electron transport. Conclusions: SHIELD-HIT12A is an interesting alternative ion transport engine. Apart from being a flexible particle therapy research tool, it can also serve as a back end for a MC ion treatment planning system. More information about SHIELD-HIT12A and a demo version can be found on http://www.shieldhit.org.

Chasing Ghosts in Space Radiobiology Research: The Lost Focus on Non-Targeted Effects

NASA Astrophysics Data System (ADS)

Cucinotta, Francis; Saganti, Premkumar; Cacao, Eliedonna

2016-07-01

The doses and dose-rates of astronaut exposures to galactic cosmic rays (GCR) are accurately known, and lead to particle hits per cell nucleus from high charge and energy (HZE) particles of much less than one hit per cell per week. A large number of experiments have shown that additivity of biological effects is a valid assumption for space radiation exposures, while experiments at higher doses and dose-rates than occur in space continue to be a focus of the majority of space radiobiology research. Furthermore HZE particle exposures with mono-energetic particles manifest themselves as a mixed-radiation field due to the contributions of delta-rays and the random impact parameter of a particles track core to DNA and non-DNA targets in cells and tissues. The mixed-field manifestation of mono-energetic HZE particle exposures is well known from theoretical studies of microdosimetry and track structure. Additional mixed-field effects occur for single species experiments due to nuclear fragmentation in particle accelerator beam-lines and biological samples along with energy straggling. In contrast to these well known aspects of space radiobiology there are many open questions on the contribution of non-targeted effects to low dose and dose-rate exposures. Non-targeted effects (NTEs) include bystander effects and genomic instability, and have been shown to be the most important outstanding question for reducing uncertainties in space radiation cancer risk assessment. The dose-rate and radiation quality dependence of NTE's has not been established, while there is an over-arching need to develop 21st century experimental models of human cancer risk. We review possible mechanisms of NTE's and how new experiments to address these issues could be designed.

The Plasma Proton Environment within Saturn's inner magnetosphere as Observed by the Cassini Plasma Spectrometer (CAPS) during Saturn Orbit Insertion

NASA Astrophysics Data System (ADS)

Sittler, E. C., Jr.; Elrod, M. K.; Johnson, R. E.; Cooper, J. F.; Tseng, W. L.; Smith, H. T.; Chornay, D. J.; Shappirio, M.; Simpson, D. G.

2017-12-01

In analyzing the Cassini data between Saturn's A-ring outer edge and Mimas' L shell numerous inconsistencies have been reported in estimates of total ionic charge and electron density. The primary focus of our work is to understand these inconsistencies. We present our recent discovery of plasma protons during Saturn Orbit Insertion (SOI) outbound pass of the magnetospheric region between the F and G rings. We also searched for H2+ ions but no such events were found. The discovery of protons was made possible by a recent analysis of the CAPS Ion Mass Spectrometer's (IMS's) time-of-flight (TOF) composition data in a mode of reduced post-acceleration voltage at 6 kV instead of the usual 14.6 kV. All previous work for this region had not considered the TOF data. The new proton analysis was enabled by minimum scattering of 6 kV protons in the instrument's ultrathin carbon foils (CF), in comparison to larger scattering for the heavier ions such as for O+ and O2+. We use a SIMION model of the CAPS IMS including the effects of energy straggling and scattering by the instrument's CFs in an attempt to understand the TOF composition data for the heavier ions. This analysis within the uncertainties of the instrument allows us to estimate the relative abundances of the heavier ions and thus run our 2D velocity ion moments code to get ion densities, temperatures and velocities during the SOI outbound pass through the F-ring and G-ring gap. Comparisons with other data sets will be made.

Neutron-induced fission cross-section measurement of 234U with quasi-monoenergetic beams in the keV and MeV range using micromegas detectors

NASA Astrophysics Data System (ADS)

Tsinganis, A.; Kokkoris, M.; Vlastou, R.; Kalamara, A.; Stamatopoulos, A.; Kanellakopoulos, A.; Lagoyannis, A.; Axiotis, M.

2017-09-01

Accurate data on neutron-induced fission cross-sections of actinides are essential for the design of advanced nuclear reactors based either on fast neutron spectra or alternative fuel cycles, as well as for the reduction of safety margins of existing and future conventional facilities. The fission cross-section of 234U was measured at incident neutron energies of 560 and 660 keV and 7.5 MeV with a setup based on `microbulk' Micromegas detectors and the same samples previously used for the measurement performed at the CERN n_TOF facility (Karadimos et al., 2014). The 235U fission cross-section was used as reference. The (quasi-)monoenergetic neutron beams were produced via the 7Li(p,n) and the 2H(d,n) reactions at the neutron beam facility of the Institute of Nuclear and Particle Physics at the `Demokritos' National Centre for Scientific Research. A detailed study of the neutron spectra produced in the targets and intercepted by the samples was performed coupling the NeuSDesc and MCNPX codes, taking into account the energy spread, energy loss and angular straggling of the beam ions in the target assemblies, as well as contributions from competing reactions and neutron scattering in the experimental setup. Auxiliary Monte-Carlo simulations were performed with the FLUKA code to study the behaviour of the detectors, focusing particularly on the reproduction of the pulse height spectra of α-particles and fission fragments (using distributions produced with the GEF code) for the evaluation of the detector efficiency. An overview of the developed methodology and preliminary results are presented.

Allelic interaction of F1 pollen sterility loci and abnormal chromosome behaviour caused pollen sterility in intersubspecific autotetraploid rice hybrids.

PubMed

He, J H; Shahid, M Q; Li, Y J; Guo, H B; Cheng, X A; Liu, X D; Lu, Y G

2011-08-01

The intersubspecific hybrids of autotetraploid rice has many features that increase rice yield, but lower seed set is a major hindrance in its utilization. Pollen sterility is one of the most important factors which cause intersubspecific hybrid sterility. The hybrids with greater variation in seed set were used to study how the F(1) pollen sterile loci (S-a, S-b, and S-c) interact with each other and how abnormal chromosome behaviour and allelic interaction of F(1) sterility loci affect pollen fertility and seed set of intersubspecific autotetraploid rice hybrids. The results showed that interaction between pollen sterility loci have significant effects on the pollen fertility of autotetraploid hybrids, and pollen fertility further decreased with an increase in the allelic interaction of F(1) pollen sterility loci. Abnormal ultra-structure and microtubule distribution patterns during pollen mother cell (PMC) meiosis were found in the hybrids with low pollen fertility in interphase and leptotene, suggesting that the effect-time of pollen sterility loci interaction was very early. There were highly significant differences in the number of quadrivalents and bivalents, and in chromosome configuration among all the hybrids, and quadrivalents decreased with an increase in the seed set of autotetraploid hybrids. Many different kinds of chromosomal abnormalities, such as chromosome straggling, chromosome lagging, asynchrony of chromosome disjunction, and tri-fission were found during the various developmental stages of PMC meiosis. All these abnormalities were significantly higher in sterile hybrids than in fertile hybrids, suggesting that pollen sterility gene interactions tend to increase the chromosomal abnormalities which cause the partial abortion of male gametes and leads to the decline in the seed set of the autotetraploid rice hybrids. © 2011 The Author(s).

Dosimetric comparisons of carbon ion treatment plans for 1D and 2D ripple filters with variable thicknesses

NASA Astrophysics Data System (ADS)

Printz Ringbæk, Toke; Weber, Uli; Santiago, Alina; Simeonov, Yuri; Fritz, Peter; Krämer, Michael; Wittig, Andrea; Bassler, Niels; Engenhart-Cabillic, Rita; Zink, Klemens

2016-06-01

A ripple filter (RiFi)—also called mini-ridge filter—is a passive energy modulator used in particle beam treatments that broadens the Bragg peak (BP) as a function of its maximum thickness. The number of different energies requested from the accelerator can thus be reduced, which significantly reduces the treatment time. A new second generation RiFi with 2D groove shapes was developed using rapid prototyping, which optimizes the beam-modulating material and enables RiFi thicknesses of up to 6 mm. Carbon ion treatment plans were calculated using the standard 1D 3 mm thick RiFi and the new 4 and 6 mm 2D RiFis for spherical planning target volumes (PTVs) in water, eight stage I non-small cell lung cancer cases, four skull base chordoma cases and three prostate cancer cases. TRiP98 was used for treatment planning with facility-specific base data calculated with the Monte Carlo code SHIELD-HIT12A. Dose-volume-histograms, spatial dose distributions and dosimetric indexes were used for plan evaluation. Plan homogeneity and conformity of thinner RiFis were slightly superior to thicker RiFis but satisfactory results were obtained for all RiFis investigated. For the 6 mm RiFi, fine structures in the dose distribution caused by the larger energy steps were observed at the PTV edges, in particular for superficial and/or very small PTVs but performances for all RiFis increased with penetration depth due to straggling and scattering effects. Plans with the new RiFi design yielded for the studied cases comparable dosimetric results to the standard RiFi while the 4 and 6 mm RiFis lowered the irradiation time by 25-30% and 45-49%, respectively.

Review on Monte-Carlo Tools for Simulating Relativistic Runaway Electron Avalanches and the Propagation of TerretrialTerrestrial-Gamma Ray Flashes in the Atmosphere

NASA Astrophysics Data System (ADS)

Sarria, D.

2016-12-01

The field of High Energy Atmospheric Physics (HEAP) includes the study of energetic events related to thunderstorms, such as Terrestrial Gamma-ray Flashes (TGF), associated electron-positron beams (TEB), gamma-ray glows and Thunderstorm Ground Enhancements (TGE). Understanding these phenomena requires accurate models for the interaction of particles with atmospheric air and electro-magnetic fields in the <100 MeV energy range. This study is the next step of the work presented in [C. Rutjes et al., 2016] that compared the performances of various codes in the absence of electro-magnetic fields. In the first part, we quantify simple but informative test cases of electrons in various electric field profiles. We will compare the avalanche length (of the Relativistic Runaway Electron Avalanche (RREA) process), the photon/electron spectra and spatial scattering. In particular, we test the effect of the low-energy threshold, that was found to be very important [Skeltved et al., 2014]. Note that even without a field, it was found to be important because of the straggling effect [C. Rutjes et al., 2016]. For this first part, we will be comparing GEANT4 (different flavours), FLUKA and the custom made code GRRR. In the second part, we test the propagation of these high energy particles in the atmosphere, from production altitude (around 10 km to 18 km) to satellite altitude (600 km). We use a simple and clearly fixed set-up for the atmospheric density, the geomagnetic field, the initial conditions, and the detection conditions of the particles. For this second part, we will be comparing GEANT4 (different flavours), FLUKA/CORSIKA and the custom made code MC-PEPTITA. References : C. Rutjes et al., 2016. Evaluation of Monte Carlo tools for high energy atmospheric physics. Geosci. Model Dev. Under review. Skeltved, A. B. et al., 2014. Modelling the relativistic runaway electron avalanche and the feedback mechanism with geant4. JGRA, doi :10.1002/2014JA020504.

Influence of bill and foot morphology on the ectoparasites of barn owls.

PubMed

Bush, Sarah E; Villa, Scott M; Boves, Than J; Brewer, Dallas; Belthoff, James R

2012-04-01

Preening is the principle behavioral defense used by birds to combat ectoparasites. Most birds have a small overhang at the tip of their bills that is used to shear through the tough cuticle of ectoparasitic arthropods, making preening much more efficient. Birds may also scratch with their feet to defend against ectoparasites. This is particularly important for removing ectoparasites on the head, which birds cannot preen. Scratching may be enhanced by the comb-like serrations that are found on the claws of birds in many avian families. We examined the prevalence and intensity of ectoparasites of barn owls (Tyto alba pratincola) in southern Idaho in relation to bill hook length and morphological characteristics of the pectinate claw. The barn owls in our study were infested with 3 species of lice (Phthiraptera: Ischnocera): Colpocephalum turbinatum , Kurodaia subpachygaster, and Strigiphilus aitkeni . Bill hook length was associated with the prevalence of these lice. Owls with longer hooks were more likely to be infested with lice. Conventional wisdom suggests that the bill morphology of raptors has been shaped by selection for efficient foraging; our data suggest that hook morphology may also play a role in ectoparasite defense. The number of teeth on the pectinate claw was also associated with the prevalence of lice. Owls that had claws with more teeth were less likely to be infested with lice, which suggests that larger pectinate claws may offer relatively more protection against ectoparasitic lice. Experiments that manipulate the bill hook and pectinate claw are needed to confirm whether these host characters are involved in ectoparasite defense. Finally, we recovered mammalian ectoparasites from 4 barn owls. We recovered species of mammalian lice (Phthiraptera:Anoplura) and fleas (Siphonaptera) that are commonly found on microtine rodents. The owls probably acquired these parasites from recently eaten prey. This represents 1 of the few documented cases of parasites "straggling" from prey to predator.

SU-C-BRC-06: OpenCL-Based Cross-Platform Monte Carlo Simulation Package for Carbon Ion Therapy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Qin, N; Tian, Z; Pompos, A

2016-06-15

Purpose: Monte Carlo (MC) simulation is considered to be the most accurate method for calculation of absorbed dose and fundamental physical quantities related to biological effects in carbon ion therapy. Its long computation time impedes clinical and research applications. We have developed an MC package, goCMC, on parallel processing platforms, aiming at achieving accurate and efficient simulations for carbon therapy. Methods: goCMC was developed under OpenCL framework. It supported transport simulation in voxelized geometry with kinetic energy up to 450 MeV/u. Class II condensed history algorithm was employed for charged particle transport with stopping power computed via Bethe-Bloch equation. Secondarymore » electrons were not transported with their energy locally deposited. Energy straggling and multiple scattering were modeled. Production of secondary charged particles from nuclear interactions was implemented based on cross section and yield data from Geant4. They were transported via the condensed history scheme. goCMC supported scoring various quantities of interest e.g. physical dose, particle fluence, spectrum, linear energy transfer, and positron emitting nuclei. Results: goCMC has been benchmarked against Geant4 with different phantoms and beam energies. For 100 MeV/u, 250 MeV/u and 400 MeV/u beams impinging to a water phantom, range difference was 0.03 mm, 0.20 mm and 0.53 mm, and mean dose difference was 0.47%, 0.72% and 0.79%, respectively. goCMC can run on various computing devices. Depending on the beam energy and voxel size, it took 20∼100 seconds to simulate 10{sup 7} carbons on an AMD Radeon GPU card. The corresponding CPU time for Geant4 with the same setup was 60∼100 hours. Conclusion: We have developed an OpenCL-based cross-platform carbon MC simulation package, goCMC. Its accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon therapy.« less

Comparisons of Integrated Radiation Transport Models with Microdosimetry Data in Spaceflight

NASA Technical Reports Server (NTRS)

Cucinotta, Francis A.; Nikjoo, H.; Kim, M. Y.; Hu, X.; Dicello, J. F.; Pisacane, V. L.

2006-01-01

Astronauts are exposed to galactic cosmic rays (GCR), trapped protons, and possible solar particle events (SPE) during spaceflight. For such complicated mixtures of radiation types and kinetic energies, tissue equivalent proportional counters (TEPC's) represent a simple time-dependent approach for radiation monitoring. Of interest in radiation protection is the average quality factor of a radiation field defined as a function of linear energy transfer, LET, Q(sub ave)(LET). However TEPC's measure the average quality factors as a function of lineal energy (y), Q(sub ave)(y) defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector. Using integrated space radiation models that includes the transport code HZETRN/BRYNTRN, the quantum nuclear interaction model, QMSFRG, and results from Monte-Carlo track simulations of TEPC's response to ions, we consider comparisons of model calculations to TEPC results from NASA missions in low Earth orbit and make predictions for lunar and Mars missions. Good agreement between the model and measured spectra from past NASA missions is found. A finding of this work is that TEPC's values for trapped or solar protons of Q(sub ave)(y) range from 1.9-2.5, overestimating Q(sub ave)(LET), which ranges from 1.4-1.6 with both quantities increasing with shielding depth due to nuclear secondaries Comparisons for the complete GCR spectra show that Q(sub ave)(LET) for GCR is approximately 3.5-4.5, while TEPC's measure 2.9-3.4 for Q(sub ave)(y) with the GCR values decreasing with depth as heavy ions are absorbed in shielding material. Our results support the use of TEPC's for space radiation environmental monitoring when computational analysis is used for proper data interpretation.

Modulation power of porous materials and usage as ripple filter in particle therapy.

PubMed

Printz Ringbæk, Toke; Simeonov, Yuri; Witt, Matthias; Engenhart-Cabillic, Rita; Kraft, Gerhard; Zink, Klemens; Weber, Uli

2017-04-07

Porous materials with microscopic structures like foam, sponges, lung tissues and lung substitute materials have particular characteristics, which differ from those of solid materials. Ion beams passing through porous materials show much stronger energy straggling than expected for non-porous solid materials of the same thickness. This effect depends on the microscopic fine structure, the density and the thickness of the porous material. The beam-modulating effect from a porous plate enlarges the Bragg peak, yielding similar benefits in irradiation time reduction as a ripple filter. A porous plate can additionally function as a range shifter, which since a higher energy can be selected for the same penetration depth in the body reduces the scattering at the beam line and therefore improves the lateral fall-off. Bragg curve measurements of ion beams passing through different porous materials have been performed in order to determine the beam modulation effect of each. A mathematical model describing the correlation between the mean material density, the porous pore structure size and the strength of the modulation has been developed and a new material parameter called 'modulation power' is defined as the square of the Gaussian sigma divided by the mean water-equivalent thickness of the porous absorber. Monte Carlo simulations have been performed in order to validate the model and to investigate the Bragg peak enlargement, the scattering effects of porosity and the lateral beam width at the end of the beam range. The porosity is found to only influence the lateral scattering in a negligible way. As an example of a practical application, it is found that a 20 mm and 50 mm plate of Gammex LN300 performs similar to a 3 mm and 6 mm ripple filter, respectively, and at the same time can improve the sharpness of the lateral beam due to its multifunctionality as a ripple filter and a range shifter.

A fast GPU-based Monte Carlo simulation of proton transport with detailed modeling of nonelastic interactions.

PubMed

Wan Chan Tseung, H; Ma, J; Beltran, C

2015-06-01

Very fast Monte Carlo (MC) simulations of proton transport have been implemented recently on graphics processing units (GPUs). However, these MCs usually use simplified models for nonelastic proton-nucleus interactions. Our primary goal is to build a GPU-based proton transport MC with detailed modeling of elastic and nonelastic proton-nucleus collisions. Using the cuda framework, the authors implemented GPU kernels for the following tasks: (1) simulation of beam spots from our possible scanning nozzle configurations, (2) proton propagation through CT geometry, taking into account nuclear elastic scattering, multiple scattering, and energy loss straggling, (3) modeling of the intranuclear cascade stage of nonelastic interactions when they occur, (4) simulation of nuclear evaporation, and (5) statistical error estimates on the dose. To validate our MC, the authors performed (1) secondary particle yield calculations in proton collisions with therapeutically relevant nuclei, (2) dose calculations in homogeneous phantoms, (3) recalculations of complex head and neck treatment plans from a commercially available treatment planning system, and compared with (GEANT)4.9.6p2/TOPAS. Yields, energy, and angular distributions of secondaries from nonelastic collisions on various nuclei are in good agreement with the (GEANT)4.9.6p2 Bertini and Binary cascade models. The 3D-gamma pass rate at 2%-2 mm for treatment plan simulations is typically 98%. The net computational time on a NVIDIA GTX680 card, including all CPU-GPU data transfers, is ∼ 20 s for 1 × 10(7) proton histories. Our GPU-based MC is the first of its kind to include a detailed nuclear model to handle nonelastic interactions of protons with any nucleus. Dosimetric calculations are in very good agreement with (GEANT)4.9.6p2/TOPAS. Our MC is being integrated into a framework to perform fast routine clinical QA of pencil-beam based treatment plans, and is being used as the dose calculation engine in a clinically applicable MC-based IMPT treatment planning system. The detailed nuclear modeling will allow us to perform very fast linear energy transfer and neutron dose estimates on the GPU.

Interpretation of TEPC Measurements in Space Flights for Radiation Monitoring

NASA Technical Reports Server (NTRS)

Kim, Myung-Hee Y.; Nikjoo, Hooshang; Dicello, John F.; Pisacane, Vincent; Cucinotta, Francis A.

2007-01-01

For the proper interpretation of radiation data measured in space, the results of integrated radiation transport models were compared with the tissue equivalent proportional counter (TEPC) measurements. TEPC is a simple, time-dependent approach to radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Qave(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Qave(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. The deviation of y from LET is caused by energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. The response distribution functions of the wall-less and walled TEPCs were calculated from Monte-Carlo track simulations. Using an integrated space radiation model (which includes the transport codes HZETRN and BRYNTRN, and the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions from Monte-Carlo track simulations, we compared model calculations with the walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. Good agreement was found for Qave(y) between the model and measured spectra from past NASA missions. The Qave(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Qave(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Qave(y). The GCR values are decreasing with the shield thickness. Our analysis of the measurements of TEPCs can be used for a proper interpretation of observed data of monitoring the space radiation environment.

Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy

PubMed Central

Qin, Nan; Pinto, Marco; Tian, Zhen; Dedes, Georgios; Pompos, Arnold; Jiang, Steve B.; Parodi, Katia; Jia, Xun

2017-01-01

Monte Carlo (MC) simulation is considered as the most accurate method for calculation of absorbed dose and fundamental physics quantities related to biological effects in carbon ion therapy. To improve its computational efficiency, we have developed a GPU-oriented fast MC package named goCMC, for carbon therapy. goCMC simulates particle transport in voxelized geometry with kinetic energy up to 450 MeV/u. Class II condensed history simulation scheme with a continuous slowing down approximation was employed. Energy straggling and multiple scattering were modeled. δ-electrons were terminated with their energy locally deposited. Four types of nuclear interactions were implemented in goCMC, i.e., carbon-hydrogen, carbon-carbon, carbon-oxygen and carbon-calcium inelastic collisions. Total cross section data from Geant4 were used. Secondary particles produced in these interactions were sampled according to particle yield with energy and directional distribution data derived from Geant4 simulation results. Secondary charged particles were transported following the condensed history scheme, whereas secondary neutral particles were ignored. goCMC was developed under OpenCL framework and is executable on different platforms, e.g. GPU and multi-core CPU. We have validated goCMC with Geant4 in cases with different beam energy and phantoms including four homogeneous phantoms, one heterogeneous half-slab phantom, and one patient case. For each case 3 × 107 carbon ions were simulated, such that in the region with dose greater than 10% of maximum dose, the mean relative statistical uncertainty was less than 1%. Good agreements for dose distributions and range estimations between goCMC and Geant4 were observed. 3D gamma passing rates with 1%/1 mm criterion were over 90% within 10%) isodose line except in two extreme cases, and those with 2%/1 mm criterion were all over 96%. Efficiency and code portability were tested with different GPUs and CPUs. Depending on the beam energy and voxel size, the computation time to simulate 107 carbons was 9.9–125 sec, 2.5–50 sec and 60–612 sec on an AMD Radeon GPU card, an NVidia GeForce GTX 1080 GPU card and an Intel Xeon E5-2640 CPU, respectively. The combined accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon ion therapy. PMID:28140352

Gaseous time projection chambers for rare event detection: results from the T-REX project. I. Double beta decay

NASA Astrophysics Data System (ADS)

Irastorza, I. G.; Aznar, F.; Castel, J.; Cebrián, S.; Dafni, T.; Galán, J.; Garcia, J. A.; Garza, J. G.; Gómez, H.; Herrera, D. C.; Iguaz, F. J.; Luzon, G.; Mirallas, H.; Ruiz, E.; Seguí, L.; Tomás, A.

2016-01-01

As part of the T-REX project, a number of R&D and prototyping activities have been carried out during the last years to explore the applicability of gaseous Time Projection Chambers (TPCs) with Micromesh Gas Structures (Micromegas) in rare event searches like double beta decay, axion research and low-mass WIMP searches. In both this and its companion paper, we compile the main results of the project and give an outlook of application prospects for this detection technique. While in the companion paper we focus on axions and WIMPs, in this paper we focus on the results regarding the measurement of the double beta decay (DBD) of 136Xe in a high pressure Xe (HPXe) TPC. Micromegas of the microbulk type have been extensively studied in high pressure Xe and Xe mixtures. Particularly relevant are the results obtained in Xe + trimethylamine (TMA) mixtures, showing very promising results in terms of gain, stability of operation, and energy resolution at high pressures up to 10 bar. The addition of TMA at levels of ~ 1% reduces electron diffusion by up to a factor of 10 with respect to pure Xe, improving the quality of the topological pattern, with a positive impact on the discrimination capability. Operation with a medium size prototype of 30 cm diameter and 38 cm of drift (holding about 1 kg of Xe at 10 bar in the fiducial volume, enough to contain high energy electron tracks in the detector volume) has allowed to test the detection concept in realistic experimental conditions. Microbulk Micromegas are able to image the DBD ionization signature with high quality while, at the same time, measuring its energy deposition with a resolution of at least a ~ 3% FWHM @ Qββ. This value was experimentally demonstrated for high-energy extended tracks at 10 bar, and is probably improvable down to the ~ 1% FWHM levels as extrapolated from low energy events. In addition, first results on the topological signature information (one straggling track ending in two blobs) show promising background discrimination capabilities out of reach of other experimental implementations. Moreover, microbulk Micromegas have very low levels of intrinsic radioactivity, and offer cost-effective scaling-up options. All these results demonstrate that Micromegas-read HPXe TPC remains a very competitive technique for the next generation DBD experiments.

Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy

NASA Astrophysics Data System (ADS)

Qin, Nan; Pinto, Marco; Tian, Zhen; Dedes, Georgios; Pompos, Arnold; Jiang, Steve B.; Parodi, Katia; Jia, Xun

2017-05-01

Monte Carlo (MC) simulation is considered as the most accurate method for calculation of absorbed dose and fundamental physics quantities related to biological effects in carbon ion therapy. To improve its computational efficiency, we have developed a GPU-oriented fast MC package named goCMC, for carbon therapy. goCMC simulates particle transport in voxelized geometry with kinetic energy up to 450 MeV u-1. Class II condensed history simulation scheme with a continuous slowing down approximation was employed. Energy straggling and multiple scattering were modeled. δ-electrons were terminated with their energy locally deposited. Four types of nuclear interactions were implemented in goCMC, i.e. carbon-hydrogen, carbon-carbon, carbon-oxygen and carbon-calcium inelastic collisions. Total cross section data from Geant4 were used. Secondary particles produced in these interactions were sampled according to particle yield with energy and directional distribution data derived from Geant4 simulation results. Secondary charged particles were transported following the condensed history scheme, whereas secondary neutral particles were ignored. goCMC was developed under OpenCL framework and is executable on different platforms, e.g. GPU and multi-core CPU. We have validated goCMC with Geant4 in cases with different beam energy and phantoms including four homogeneous phantoms, one heterogeneous half-slab phantom, and one patient case. For each case 3× {{10}7} carbon ions were simulated, such that in the region with dose greater than 10% of maximum dose, the mean relative statistical uncertainty was less than 1%. Good agreements for dose distributions and range estimations between goCMC and Geant4 were observed. 3D gamma passing rates with 1%/1 mm criterion were over 90% within 10% isodose line except in two extreme cases, and those with 2%/1 mm criterion were all over 96%. Efficiency and code portability were tested with different GPUs and CPUs. Depending on the beam energy and voxel size, the computation time to simulate {{10}7} carbons was 9.9-125 s, 2.5-50 s and 60-612 s on an AMD Radeon GPU card, an NVidia GeForce GTX 1080 GPU card and an Intel Xeon E5-2640 CPU, respectively. The combined accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon ion therapy.

Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy.

PubMed

Qin, Nan; Pinto, Marco; Tian, Zhen; Dedes, Georgios; Pompos, Arnold; Jiang, Steve B; Parodi, Katia; Jia, Xun

2017-05-07

Monte Carlo (MC) simulation is considered as the most accurate method for calculation of absorbed dose and fundamental physics quantities related to biological effects in carbon ion therapy. To improve its computational efficiency, we have developed a GPU-oriented fast MC package named goCMC, for carbon therapy. goCMC simulates particle transport in voxelized geometry with kinetic energy up to 450 MeV u -1 . Class II condensed history simulation scheme with a continuous slowing down approximation was employed. Energy straggling and multiple scattering were modeled. δ-electrons were terminated with their energy locally deposited. Four types of nuclear interactions were implemented in goCMC, i.e. carbon-hydrogen, carbon-carbon, carbon-oxygen and carbon-calcium inelastic collisions. Total cross section data from Geant4 were used. Secondary particles produced in these interactions were sampled according to particle yield with energy and directional distribution data derived from Geant4 simulation results. Secondary charged particles were transported following the condensed history scheme, whereas secondary neutral particles were ignored. goCMC was developed under OpenCL framework and is executable on different platforms, e.g. GPU and multi-core CPU. We have validated goCMC with Geant4 in cases with different beam energy and phantoms including four homogeneous phantoms, one heterogeneous half-slab phantom, and one patient case. For each case [Formula: see text] carbon ions were simulated, such that in the region with dose greater than 10% of maximum dose, the mean relative statistical uncertainty was less than 1%. Good agreements for dose distributions and range estimations between goCMC and Geant4 were observed. 3D gamma passing rates with 1%/1 mm criterion were over 90% within 10% isodose line except in two extreme cases, and those with 2%/1 mm criterion were all over 96%. Efficiency and code portability were tested with different GPUs and CPUs. Depending on the beam energy and voxel size, the computation time to simulate [Formula: see text] carbons was 9.9-125 s, 2.5-50 s and 60-612 s on an AMD Radeon GPU card, an NVidia GeForce GTX 1080 GPU card and an Intel Xeon E5-2640 CPU, respectively. The combined accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon ion therapy.

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

NASA Technical Reports Server (NTRS)

Cucinotta, Francis A.; Plante, Ianik; Ponomarev, Artem L.; Kim, Myung-Hee Y.

2013-01-01

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

NASA Technical Reports Server (NTRS)

Badavi, Francis F.; Xapsos, Mike

2008-01-01

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Sparse-view proton computed tomography using modulated proton beams.

PubMed

Lee, Jiseoc; Kim, Changhwan; Min, Byungjun; Kwak, Jungwon; Park, Seyjoon; Lee, Se Byeong; Park, Sungyong; Cho, Seungryong

2015-02-01

Proton imaging that uses a modulated proton beam and an intensity detector allows a relatively fast image acquisition compared to the imaging approach based on a trajectory tracking detector. In addition, it requires a relatively simple implementation in a conventional proton therapy equipment. The model of geometric straight ray assumed in conventional computed tomography (CT) image reconstruction is however challenged by multiple-Coulomb scattering and energy straggling in the proton imaging. Radiation dose to the patient is another important issue that has to be taken care of for practical applications. In this work, the authors have investigated iterative image reconstructions after a deconvolution of the sparsely view-sampled data to address these issues in proton CT. Proton projection images were acquired using the modulated proton beams and the EBT2 film as an intensity detector. Four electron-density cylinders representing normal soft tissues and bone were used as imaged object and scanned at 40 views that are equally separated over 360°. Digitized film images were converted to water-equivalent thickness by use of an empirically derived conversion curve. For improving the image quality, a deconvolution-based image deblurring with an empirically acquired point spread function was employed. They have implemented iterative image reconstruction algorithms such as adaptive steepest descent-projection onto convex sets (ASD-POCS), superiorization method-projection onto convex sets (SM-POCS), superiorization method-expectation maximization (SM-EM), and expectation maximization-total variation minimization (EM-TV). Performance of the four image reconstruction algorithms was analyzed and compared quantitatively via contrast-to-noise ratio (CNR) and root-mean-square-error (RMSE). Objects of higher electron density have been reconstructed more accurately than those of lower density objects. The bone, for example, has been reconstructed within 1% error. EM-based algorithms produced an increased image noise and RMSE as the iteration reaches about 20, while the POCS-based algorithms showed a monotonic convergence with iterations. The ASD-POCS algorithm outperformed the others in terms of CNR, RMSE, and the accuracy of the reconstructed relative stopping power in the region of lung and soft tissues. The four iterative algorithms, i.e., ASD-POCS, SM-POCS, SM-EM, and EM-TV, have been developed and applied for proton CT image reconstruction. Although it still seems that the images need to be improved for practical applications to the treatment planning, proton CT imaging by use of the modulated beams in sparse-view sampling has demonstrated its feasibility.

Gaseous time projection chambers for rare event detection: results from the T-REX project. I. Double beta decay

DOE Office of Scientific and Technical Information (OSTI.GOV)

Irastorza, I.G.; Aznar, F.; Castel, J., E-mail: igor.irastorza@cern.ch, E-mail: faznar@unizar.es, E-mail: jfcastel@unizar.es

2016-01-01

As part of the T-REX project, a number of R and D and prototyping activities have been carried out during the last years to explore the applicability of gaseous Time Projection Chambers (TPCs) with Micromesh Gas Structures (Micromegas) in rare event searches like double beta decay, axion research and low-mass WIMP searches. In both this and its companion paper, we compile the main results of the project and give an outlook of application prospects for this detection technique. While in the companion paper we focus on axions and WIMPs, in this paper we focus on the results regarding the measurementmore » of the double beta decay (DBD) of {sup 136}Xe in a high pressure Xe (HPXe) TPC. Micromegas of the microbulk type have been extensively studied in high pressure Xe and Xe mixtures. Particularly relevant are the results obtained in Xe + trimethylamine (TMA) mixtures, showing very promising results in terms of gain, stability of operation, and energy resolution at high pressures up to 10 bar. The addition of TMA at levels of ∼ 1% reduces electron diffusion by up to a factor of 10 with respect to pure Xe, improving the quality of the topological pattern, with a positive impact on the discrimination capability. Operation with a medium size prototype of 30 cm diameter and 38 cm of drift (holding about 1 kg of Xe at 10 bar in the fiducial volume, enough to contain high energy electron tracks in the detector volume) has allowed to test the detection concept in realistic experimental conditions. Microbulk Micromegas are able to image the DBD ionization signature with high quality while, at the same time, measuring its energy deposition with a resolution of at least a ∼ 3% FWHM @ Q{sub ββ}. This value was experimentally demonstrated for high-energy extended tracks at 10 bar, and is probably improvable down to the ∼ 1% FWHM levels as extrapolated from low energy events. In addition, first results on the topological signature information (one straggling track ending in two blobs) show promising background discrimination capabilities out of reach of other experimental implementations. Moreover, microbulk Micromegas have very low levels of intrinsic radioactivity, and offer cost-effective scaling-up options. All these results demonstrate that Micromegas-read HPXe TPC remains a very competitive technique for the next generation DBD experiments.« less

Effect of Ion Flux (Dose Rate) in Source-Drain Extension Ion Implantation for 10-nm Node FinFET and Beyond on 300/450mm Platforms

NASA Astrophysics Data System (ADS)

Shen, Ming-Yi

The improvement of wafer equipment productivity has been a continuous effort of the semiconductor industry. Higher productivity implies lower product price, which economically drives more demand from the market. This is desired by the semiconductor manufacturing industry. By raising the ion beam current of the ion implanter for 300/450mm platforms, it is possible to increase the throughput of the ion implanter. The resulting dose rate can be comparable to the performance of conventional ion implanters or higher, depending on beam current and beam size. Thus, effects caused by higher dose rate must be investigated further. One of the major applications of ion implantation (I/I) is source-drain extension (SDE) I/I for the silicon FinFET device. This study investigated the dose rate effects on the material properties and device performance of the 10-nm node silicon FinFET. In order to gain better understanding of the dose rate effects, the dose rate study is based on Synopsys Technology CAD (TCAD) process and device simulations that are calibrated and validated using available structural silicon fin samples. We have successfully shown that the kinetic monte carlo (KMC) I/I simulation can precisely model both the silicon amorphization and the arsenic distribution in the fin by comparing the KMC simulation results with TEM images. The results of the KMC I/I simulation show that at high dose rate more activated arsenic dopants were in the source-drain extension (SDE) region. This finding matches with the increased silicon amorphization caused by the high dose-rate I/I, given that the arsenic atoms could be more easily activated by the solid phase epitaxial regrowth process. This increased silicon amorphization led to not only higher arsenic activation near the spacer edge, but also less arsenic atoms straggling into the channel. Hence, it is possible to improve the throughput of the ion implanter when the dopants are implanted at high dose rate if the same doping level with a lower wafer dose can be achieved. In addition, the leakage current might also be reduced due to less undesired dopants in the channel. However, the twin defects from the problematic Si{111} recrystallization is well-known to cause excessive leakage current to the FinFET. This drawback can offset the benefits of the high dose rate I/I mentioned above. This work produced the first attempt at simulating the electrical impact of twin defects on advanced-node (10 nm) FinFET device performance. It was found that the high dose-rate I/I causes more twin defects in the silicon fin, and the physical locations of these defects were close to the channel. The defects undesirably induced trap-assisted band-to-band tunneling near the drain, which increased the leakage current. This issue could be mitigated by using asymmetrical gate overlap/underlap design or thicker spacer for SDE I/I so that the twin defects are not located in the depletion region near the drain.

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

NASA Technical Reports Server (NTRS)

Kim, Myung-Hee Y.; Nikjoo, Hooshang; Dicello, John F.; Pisacane, Vincent; Cucinotta, Francis A.

2007-01-01

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

About the Nature of a Coherence of Light Waves

NASA Astrophysics Data System (ADS)

Demyaneko, P. O.; Zinkovskiy, Y. F.; Savenko, Y. V.

The confrontation of corpuscular and wave hypotheses was not partly stacked in frameworks of the uniform theory. Fundamental works of Fresnel and Fraunhofer on a diffraction of light have erected a wave theory in a rank of dominant true. The wave theory did not so effectively explain developments of corpuscular properties of a light. Its feeble place was also necessity for concept "coherence", introduced for explanation of a light interference. The wave coherence is properly understood as waves ability to interfere. The problem of a light coherence continues to be interesting for investigators [L], but clear understanding of its nature is not yet appeared. Because, it is unconvincing to consider an attribution to the act of separate atom radiation of electromagnetic waves "zug" of a few meters length as explanation of the coherence nature, when it has become possible to generate light pulses by duration ˜ 10 -15 s. Let's note there is the spatial separation of a primary luminous flux on two secondary ones in a basis of all ways of deriving of coherent luminous fluxes. And these fluxes are able to interfere when are brought together. Their coherence was explained that at partitioning each "zug" was bisected, and at joining again met mutual coherent "its halves". There was not disputed the question, how happens " separation of each waves "zug" in halves". M. Plank postulated (1900) heated bodies radiate electromagnetic waves not continuously, but by separate portions he called "quantums" of energy. Its have a quantity is proportional to frequency of electromagnetic waves. A. Einstein has entered (1905) a hypothesis of light quantums -- light is indeed generated as quantums, and in further it exists as a flux of quantums and interacts with matter also, as a flux of separate quantums. The term "photon" was entered by G.N. Lewis (1929) properly for emphasising of light quantums and that one underlined corpuscularity of a light. At investigation of an atom structure there was set (E. Rutherford, N. Bohr, 1911) quantums are generated in atoms at transitions of excited electrons from higher energy levels onto lowest levels. At that, there are radiated quantums-waves of electromagnetic energy into environmental space. In different light sources "working body" has the "own" structure of energy levels defining spectral characteristics of these sources. So, the development of representations about the nature of a light returns to a corpuscular hypothesis. It has become clear, that the light organically combines in itself both property of waves and the properties of particles. It depends on requirements of experiment which one from developments will be prevalent. Inseparable unity of corpuscular and wave properties is proper for all microparticles (a hypothesis De Brogle, 1924) and has received a title of "wave-corpuscle dualism". Let's make a common view about "sizes" of a photon. As was mentioned, the light pulses can have duration ˜ 10 -15 s. Spatial length of such pulses in direction of motion ˜ 10 -6 m, that comparable with a light wavelength λ . It is possible to suspect that it will be a size of a photon in direction of its propagation. An estimate of "cross" of the sizes of a photon we shall obtain by analyzing of light diffraction on a narrow slot. The angular size of central diffraction peak at decreasing of width of a slot b is increased, and it reaches 180 at b = λ . Then the light intensity promptly impinges behind it. From this it is possible to assume, that the cross sectional dimensions of photon also is comparable with λ . It is necessary to clear understand, that photon, as the wave formation, does not have sharp borders. It is possible to speak only about the sizes of area containing a dominant share of photon energy. So, photon is a spatial localized electromagnetic perturbation, that allows to allot it with properties of a particle. Essential properties of a photon are indivisibility and existence only in a motion. So, the light is a photons flux: both light wave and light electromagnetic field consist of final number of photons. At that, it is important to remember that in any light source along with spontaneous mechanism it also operates a mechanism of induced radiation, generating identical (coherent) quantums. Due to it, there is radiated a partially coherent flux, consisting of large or small groups of quantums ("quantum packets"), from any light source. In limits of a separate packet its component quantums are coherent, because all of them are originated by one quantum which has appeared spontaneously, which induced occurrence of other quantums of this packet, passing by other excited atoms. The representation about quantum packets gives clear physical explanation to concept of "light waves zugs". Quantum packet is that "zug of waves". "Quality" of a light source (in sense of its coherence) is determined by sizes of quantum packets -- the larger they the more qualitative source, radiating them. There are understandable a better coherence of a gas light sources: the atoms in gas are arranged on large distances and do not hinder for spontaneously generated light quantum to overcome without absorption or dispersion that large distance, challenging on it an induced radiation of other excited atoms. The low coherence of glow-discharge tubes is stipulated by that the radiation in them goes only from surface layer of atoms and the requirements for development of the mechanism of induced radiation are unfavorable. It is also obvious the high coherence of a laser radiation due to a positive optical back coupling. The coherent quantums of one quantum packet exist a long time inside the resonator; they are reproducing there during all this time. Due to this the lasers are capable to generate multi-km quantum packets ("zugs"). By the way, it could not to explain "by emissive opportunities" of one atom. It is understandable a division of quantum packets on semi transparent mirrors: the part of quantums of each packet simply transits through a semi transparent mirrors, and remaining ones are simply reflected from it. The model of quantum packets gives clear explanation of coherence parameters of light flux. A length of coherence is a spatial extent of a quantum packet in direction of its propagation. A coherence time is a time of flight of quantum packet by a fixed spectator. A coherence radius (size) is a spatial extent of a quantum packet in direction, perpendicular to direction of its propagation. A volume of coherence is simply a volume of quantum packet. Separately it is necessary to tell about the fact of increasing of coherence radius of a light flux, propagating in space. Iterated, including by us, assertion: "at induced transitions there are generated the same quantums as ones induced them". It is not necessary to understand it too literally. What perfect was a light source, the spectral line of its radiation always has final width. That means, there is a certain frequency dispersion of quantums, generated by source, or modules of their wave vectors. Apparently, it is necessary to expect as well certain dispersion of particular straggling of wave vectors directions inside separate quantum packets. Beginning with experimentally obtained radius of sunlight coherence on surface of the Earth, it was determined a value of angular divergence of quantum packets. With the help of the obtained thus value, there were calculated values of coherence radiuses of light, coming on the Earth from more remote stars. Obtained calculated values are well compared with experimentally obtained values of light coherence radiuses for these stars. Starting from proposed concept of quantum packets, we have given explanation to such development of wave properties of the light as interference, in particular, its variety, when superimposed coherent fluxes interfere. It is not less important from a point of view of the coordination of their explanations with our representations about luminous flux structure, there is an analysis also such developments of wave properties of light, as its interference on thin films, "Newton's ringes", etc. For explanation of this variety of interference there is no need for concept coherence, as in such interference is watched always and for a light from any sources. There is a special interest to phenomenas bound with diffraction of light, from a point of view of quantum packets model. The prime task here is to give a corresponding explanation to the content both senses of Huygens' and Huygens-Fresnel principles. These problems will be considered in following our works. [L] Mandel L., Wolf E. Optical Coherence and Quantum Optics / Cambrige, 1995