A procedure for benchmarking laboratory exposures with 1 A GeV iron ions.

A new version of the HZETRN code capable of simulating HZE ions with either laboratory or space boundary conditions is 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 asymptotic/Neumann expansions with non-perturbative corrections. The code contains energy loss with straggling, nuclear attenuation, nuclear fragmentation with energy dispersion and downshifts, and off-axis dispersion with multiple scattering under preparation. The present benchmark is for a broad directed beam for 1 A GeV iron ion beams with 2 A MeV width and four targets of polyethylene, polymethyl metachrylate, aluminum, and lead of varying thickness from 5 to 30 g/cm2. The benchmark quantities will be dose, track averaged LET, dose averaged LET, fraction of iron ion remaining, and fragment energy spectra after 23 g/cm2 of polymethyl metachrylate.     

Fragmentation studies of relativistic iron ions using plastic nuclear track detectors.

We measured fluence and fragmentation of high-energy (1 or 5 A GeV) 56Fe ions accelerated at the Alternating Gradient Synchrotron or at the NASA Space Radiation Laboratory (Brookhaven National Laboratory, NY, USA) using solid-state CR-39 nuclear track detectors. Different targets (polyethylene, PMMA, C, Al, Pb) were used to produce a large spectrum of charged fragments. CR-39 plastics were exposed both in front and behind the shielding block (thickness ranging from 5 to 30 g/cm2) at a normal incidence and low fluence. The radiation dose deposited by surviving Fe ions and charged fragments was measured behind the shield using an ionization chamber. The distribution of the measured track size was exploited to distinguish the primary 56Fe ions tracks from the lighter fragments. Measurements of projectile's fluence in front of the shield were used to determine the dose per incident particle behind the block. Simultaneous measurements of primary 56Fe ion tracks in front and behind the shield were used to evaluate the fraction of surviving iron projectiles and the total charge-changing fragmentation cross-section. These physical measurements will be used to characterize the beam used in parallel biological experiments.     

Statistical validation of HZETRN as a function of vertical cutoff rigidity using ISS measurements

Measurements taken in Low Earth Orbit (LEO) onboard the International Space Station (ISS) and transit vehicles have been extensively used to validate radiation transport models. Primarily, such comparisons were done by integrating measured data over mission or trajectory segments so that individual comparisons to model results could be made. This approach has yielded considerable information but is limited in its ability to rigorously quantify and differentiate specific model errors or uncertainties. Further, as exploration moves beyond LEO and measured data become sparse, the uncertainty estimates derived from these validation cases will no longer be applicable. Recent improvements in the underlying numerical methods used in HZETRN have resulted in significant decreases in code run time. Therefore, the large number of comparisons required to express error as a function of a physical quantity, like cutoff rigidity, are now possible. Validation can be looked at in detail over any portion of a flight trajectory (e.g. minute by minute) such that a statistically significant number of comparisons can be made. This more rigorous approach to code validation will allow the errors caused by uncertainties in the geometry models, environmental models, and nuclear physics models to be differentiated and quantified. It will also give much better guidance for future model development. More importantly, it will allow a quantitative means of extrapolating uncertainties in LEO to free space. In this work, measured data taken onboard the ISS during solar maximum are compared to results obtained with the particle transport code HZETRN. Comparisons are made at a large number (∼77,000) of discrete time intervals, allowing error estimates to be given as a function of cutoff rigidity. It is shown that HZETRN systematically underestimates exposure quantities at high cutoff rigidity. The errors are likely associated with increased angular variation in the geomagnetic field near the equator, the lack of pion production in HZETRN, and errors in high energy nuclear physics models, and will be the focus of future work.     

Nuclear fragmentation database for GCR transport code development

A critical need for NASA is the ability to accurately model the transport of heavy ions in the Galactic Cosmic Rays (GCR) through matter, including spacecraft walls, equipment racks, etc. Nuclear interactions are of great importance in the GCR transport problem, as they can cause fragmentation of the incoming ion into lighter ions. Since the radiation dose delivered by a particle is proportional to the square of (charge/velocity), fragmentation reduces the dose delivered by incident ions. The other mechanism by which dose can be reduced is ionization energy loss, which can lead to some particles stopping in the shielding. This is the conventional notion of shielding, but it is not applicable to human spaceflight since the particles in the GCR tend to be too energetic to be stopped in the relatively thin shielding that is possible within payload mass constraints. Our group has measured a large number of fragmentation cross sections, intended to be used as input to, or for validation of, NASA’s radiation transport models. A database containing over 200 charge-changing cross sections and over 2000 fragment production cross sections has been compiled. In this report, we examine in detail the contrast between fragment measurements at large acceptance and small acceptance. We use output from the PHITS Monte Carlo code to test our assumptions using as an example ⁴⁰Ar data (and simulated data) at a beam energy of 650 MeV/nucleon. We also present preliminary analysis in which isotopic resolution was attained for beryllium fragments produced by beams of ¹⁰B and ¹¹B. Future work on the experimental data set will focus on extracting and interpreting production cross sections for light fragments.     

基于核主成分分析的多输出模型确认方法

目前对于不确定性环境下多个相关的复杂计算模型进行确认的方法存在计算困难及稳定性较差的问题。针对这类复杂计算模型,提出了一种新的基于核主成分分析(KPCA)的多输出模型确认方法。该方法将核主成分分析与面积法的思想相结合,构造了一个新的易于计算且稳定性高的模型确认指标。所提方法通过核主成分分析将相关的输出变量转化为不相关的核主成分,再对每一核主成分进行模型与实验的对比,从而避免了传统多输出模型确认方法中需要求解多个输出的联合累积分布函数的困难。由于核主成分分析(PCA)方法能够有效提取分析对象的非线性成分,因此基于核主成分分析的多输出模型确认方法较基于主成分分析的模型确认方法更为稳定,这表现在相同的实验样本数据下核主成分分析的方法具有更低的出错率。另外核主成分分析通过核主成分提取,可以实现多输出模型的降维,从而降低多输出模型确认的复杂度。所提方法既可以用于一般的多输出模型的确认,也可以用于多确认点的输出模型的确认。最后通过数值算例和工程算例证明了该方法的正确性与有效性。     

Nuclear fragmentation of high-energy heavy-ion beams in water.

As a part of the physical-technical program of the heavy-ion therapy project at GSI we have investigated the nuclear fragmentation of high-energy ion beams delivered by the heavy-ion synchrotron SIS, using water as a tissue-equivalent target. For a direct comparison of fragmentation properties, beams of 10B, 12C, 14N, and 16O were produced simultaneously as secondary beams from a primary 18O beam and separated in flight by magnetic beam analysis. The Z-distributions of beam fragments produced in the water target were measured via energy loss in a large ionisation chamber and a scintillator telescope. From these data we obtained both total and partial charge-changing cross sections. In addition we have performed Bragg measurements using two parallel-plate ionization chambers and a water target of variable length. The detailed shape of the measured Bragg curves and the measured cross sections are in good agreement with model calculations based on semi-empirical formulae.     

Computational methods for the HZETRN code.

Asymptotic expansion has been used to simplify the transport of high charge and energy ions for broad beam applications in the laboratory and space. The solution of the lowest order asymptotic term is then related to a Green's function for energy loss and straggling coupled to nuclear attenuation providing the lowest order term in a rapidly converging Neumann series for which higher order collisions terms are related to the fragmentation events including energy dispersion and downshift. The first and second Neumann corrections were evaluated numerically as a standard for further analytic approximation. The first Neumann correction is accurately evaluated over the saddle point whose width is determined by the energy dispersion and located at the downshifted ion collision energy. Introduction of the first Neumann correction leads to significant simplification of the second correction term allowing application of the mean value theorem and a second saddle point approximation. The regular dependence of the second correction spectral dependence lends hope to simple approximation to higher corrections. At sufficiently high energy nuclear cross-section variations are small allowing non-perturbative methods to all orders and renormalization of the second corrections allow accurate evaluation of the full Neumann series.     

Sources of inner radiation zone energetic helium ions: Cross-field transport versus in-situ nuclear reactions

Radial transport theory for inner radiation zone MeV He ions has been extended by combining radial diffusive transport, losses due to Coulomb friction and charge exchange reaction with local generation of ³He and ⁴He ions due to nuclear reactions taking place on the inner edge of the inner radiation zone. From interactions between high energy trapped protons and upper atmospheric constituents we have included a nuclear reaction yield helium flux source that was numerically derived from a nuclear reaction model originally developed at the Institute of Nuclear Researches of Moscow, Russia and implemented in the computer system at the University of Campinas, Brazil. Magnetospheric transport computations have been made covering the L-shell range L=1.0 to 1.6 and the resulting MeV He ion flux distributions show a strong influence of the local nuclear source mechanism on the inner zone energetic He ion content.     

试验模态分析最佳测试点选取方法的优化

作为有限元模型修正的重要依据,试验模态分析的结果受多方面因素的影响,其中很重要的一方面就是测试点的选取.基于目前试验模态分析中最常用的两种测试点选择方法——有效独立法（EI,Effective Independent）和模态动能法（MKE,Modal Kinetic Energy）,提出了一种新的测试点优化选择方法——有效独立-模态动能变异系数法（EI-CVMKE,Effective Independent-Coefficient Variance of Modal Kinetic Energy）.并提出了一种新的评估测试点选择方法优劣的准则——模态动能变异系数准则（CVMKE）.运用模态置信准则（MAC,Modal Assurance Criterion）、CVMKE和奇异值分解比（SVD ration,Singular Value Decomposition ratio）评价并比较了EI,MKE和EI-CVMKE的测试点优化结果.以印刷电路板作为算例,并通过实物模态试验验证,结果表明EI-CVMKE选取的测试点在考虑了截断模态线性独立的同时又具有较高的模态动能,有较强的抗噪声性能,在3种方法中最佳.      

Fragmentation of 1 GeV/nucleon iron ions in thick targets relevant for space exploration.

We have measured charged nuclear fragments produced by 1 GeV/nucleon 56Fe ions interacting with aluminium, polyethylene and lead. These materials are relevant for assessment of radiation risk for manned space flight. The data will be presented in a form suitable for comparison with models of nuclear fragmentation and transport, including linear energy transfer (LET) spectrum, fluence for iron and fragments, event-tack- and event-dose-averaged LET, total dose and iron contribution to dose.     

On the origin of trapped heavy ions at L = 1.4-1.6.

Aboard the NASA satellite Long Duration Exposure Facility (LDEF) heavy ions of nuclear charge Z = 8-26 were detected with energies between 15 and 50 MeV/nuc which are far below the cutoff energy required of fully stripped ions to reach the LDEF orbit. The arrival directions and the falling energy spectra of these particles are consistent with a trapped component incident in the South Atlantic Anomaly at L = 1.4-1.6. The trapped oxygen, neon and argon ions probably originate from the anomalous cosmic rays, whereas the origin of the other particles like magnesium, silicon and iron is not yet solved but may be associated with the October 89 solar energetic particle events.     

Validation of the new trapped environment AE9/AP9/SPM at low Earth orbit

The completion of the international space station (ISS) in 2011 has provided the space research community an ideal proving ground for future long duration human activities in space. Ionizing radiation measurements in ISS form the ideal tool for the validation of radiation environmental models, nuclear transport codes and nuclear reaction cross sections. Indeed, prior measurements on the space transportation system (STS; shuttle) provided vital information impacting both the environmental models and the nuclear transport code developments by indicating the need for an improved dynamic model of the low Earth orbit (LEO) trapped environment. Additional studies using thermo-luminescent detector (TLD), tissue equivalent proportional counter (TEPC) area monitors, and computer aided design (CAD) model of earlier ISS configurations, confirmed STS observations that, as input, computational dosimetry requires an environmental model with dynamic and directional (anisotropic) behavior, as well as an accurate six degree of freedom (DOF) definition of the vehicle attitude and orientation along the orbit of ISS.     

Ground-based measurements of galactic cosmic ray fragmentation in shielding.

The mean free path for nuclear interactions of galactic cosmic-rays is comparable to shielding and tissue thicknesses present in human interplanetary exploration, resulting in a significant fraction of nuclear reaction products at depth. In order to characterize the radiation field, the energy spectrum, the angular distribution, and the multiplicity of each type of secondary particles must also be known as a function of depth. Reactions can take place anywhere in a thick absorber; therefore, it is necessary to know these quantities as a function of particle energy for all particles produced. HZE transport methods are used to predict the radiation field; they are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The most cost-effective method of validation is a comparison with ground-based experimental measurements. A research program to provide such validation measurements using neon, iron and other accelerated heavy ion beams will be discussed and illustrated using results from ongoing experiments and their comparison with current transport codes. The extent to which physical measurements yield radiobiological predictions will be discussed.     

Results of dosimetric measurements in space missions.

Detector packages consisting of plastic nuclear track detectors, nuclear emulsions, and theromoluminescence detectors were exposed at different locations inside the space laboratory Spacelab and at the astronauts' body and in different sections of the MIR space station. Total dose, particle fluence rate and linear energy transfer (LET) spectra of heavy ions, number of nuclear disintegrations and fast neutron fluence rates were determined of each exposure. The dose equivalent received by the Payload specialists (PSs) were calculated from the measurements, they range from 190 microSv d-1 to 770 microSv d-1. Finally, a preliminary investigation of results from a particle telescope of two silicon detectors, first used in the last BIORACK mission on STS 76, is reported.     

Inactivation, mutation induction and repair in Bacillus subtilis spores irradiated with heavy ions.

Studies on the response of bacterial spores to accelerated heavy ions (HZE particles) help in understanding problems of space radiobiology and exobiology. Layers of spores of Bacillus subtilis strains, differing in repair capabilities, were irradiated with accelerated boron, carbon and neon ions of linear energy transfer (LET) values up to 14000 MeV cm2/g. Inactivation as measured by loss of colony forming ability and induction of mutations as measured by reversion to histidine prototrophy and resistance to 150 micrograms/ml sodium azide were tested, as well as the influence of repair processes on these effects. For inactivation, the cross-sectional values sigma plotted as a function of LET follow a saturation curve. The plateau, which is reached around a LET of 2000 MeV cm2/g, occurs at 2.5 x 10(-9) cm2, a value in good agreement with the dimensions of the spore protoplast. Lethal damage produced at LET values < 2000 MeV cm2/g is reparable. Recombination repair is more effective than excision repair. At higher LET values, lethal damage could not be reconstituted by the repair mechanisms studied. In addition, at these high LET values, the frequency of induced mutations was drastically decreased. The data support the assumption of at least two qualitatively different types of lesion, depending on the LET of the affecting heavy ion.     

Trapped particle energy spectrum in shuttle middeck

We describe the differential energy spectrum of trapped particles measured by a solid-state charged particle telescope in the mid-deck of the Space Shuttle during the period of solar maximum. The telescope was flown in two high altitude flights at 28.5° and 57° inclination. Assuming, as is normally done, that the variations of Shuttle orientation during the missions lead to average isotropic incident spectra, the observed spectrum disagrees significantly from AP8 model calculations. This indicates the need to take into consideration the variations of solid-angle direction relative to the magnetic field. The measurements show that there is a very significant flux of secondary light ions. The energy spectra of these ions does not agree with the production spectrum from radiation transport calculations based on omni-directional AP8 Max model as an input energy spectrum.

We also describe measurements of linear energy transfer spectra using a tissue equivalent proportional counter (TEPC) flown both in the mid-deck and the payload bay of the Space Shuttle. Comparisons are made between linear energy transfer spectral measurements AP8 model-based radiation transport predictions, and thermoluminescent dosimeter (TLD) measurements. The absorbed dose-rate measurements using TLD's are roughly 25% lower than the TEPC-measured dose rate measurements.