SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
There is no human epidemiology data available to estimate the heavy ion cancer risks experienced by astronauts in space. Studies of tumor induction in mice are a necessary step to estimate risks to astronauts. Previous experimental data can be better utilized to model dose response for heavy ion tumorigenesis and plan future low dose studies.
Causality is violated in the early stages of state-of-the-art heavy-ion hydrodynamic simulations. Such violations are present in up to 75% of the fluid cells in the initial time and only after 2–3 fm/c of evolution do we find that 50% of the fluid cells are definitely causal. Superluminal propagation reaches up to 15% the speed of light in some of the fluid cells. The inclusion of pre-equilibrium evolution significantly reduces the number of acausal cells. Here our findings suggests that relativistic causality may place constraints on the available parameter space of heavy-ion collision simulations when factored into more thorough statistical analyses.
Heavy-ion storage rings have relatively large momentum acceptance which allows for multiple ion species to circulate at the same time. This needs to be considered in radioactive decay measurements of highly charged ions, where atomic charge exchange reactions can significantly alter the intensities of parent and daughter ions. In this study, we investigate this effect using the decay curves of ion numbers in the recent 205 Tl 81+ bound-state beta decay experiment conducted using the Experimental Storage Ring at GSI Darmstadt. To understand the intricate dynamics of ion numbers, we present a set of differential equations that account for various atomic and nuclear reaction processes—bound state beta decay, atomic electron recombination and capture, and electron ionization. By incorporating appropriate boundary conditions, we develop a set of differential equations that accurately simulate the decay curves of various simultaneously stored ions in the storage ring: 205 Tl 81+ , 205 Pb 81+ , 205 Pb 82+ , 200 Hg 79+ , and 200 Hg 80+ . Through a quantitative comparison between simulations and experimental data, we provide insights into the detailed reaction mechanisms governing stored heavy ions within the storage ring. Our approach effectively models charge-changing processes, reduces the complexity of the experimental setup, and provides a simpler method for measuring the decay half-lives of highly charged ions in storage rings.
The IS3D particlization module simulates the emission of hadrons from heavy-ion collisions via Monte-Carlo sampling of the Cooper–Frye formula which converts fluid dynamical information into local phase-space distributions for hadrons. The code package includes multiple choices for the non-equilibrium correction to these distribution functions: the 14-moment approximation, first-order Chapman–Enskog expansion, and two types of modified equilibrium distributions. This makes it possible to explore to what extent heavy-ion experimental data are sensitive to different choices for δf n , presently the main source of theoretical uncertainty in the particlization stage. Here, we validate our particle sampler with a high degree of precision by generating several million hadron emission events from a longitudinally boost-invariant hypersurface and comparing the event-averaged particle spectra and space–time distributions to the Cooper–Frye formula.
Simulating the Space Radiation environment with Monte Carlo Codes, such as FLUKA, requires the ability to model the interactions of heavy ions as they penetrate spacecraft and crew member's bodies. Monte-Carlo-type transport codes use total interaction cross sections to determine probabilistically when a particular type of interaction has occurred. Then, at that point, a distinct event generator is employed to determine separately the results of that interaction. The space radiation environment contains a full spectrum of radiation types, including relativistic nuclei, which are the most important component for the evaluation of crew doses. Interactions between incident protons with target nuclei in the spacecraft materials and crew member's bodies are well understood. However, the situation is substantially less comfortable for incident heavier nuclei (heavy ions). We have been engaged in developing several related heavy ion interaction models based on a Quantum Molecular Dynamics-type approach for energies up through about 5 GeV per nucleon (GeV/A) as part of a NASA Consortium that includes a parallel program of cross section measurements to guide and verify this code development.
In the hydrodynamic model description of heavy ion collisions, the final-state anisotropic flows v n are linearly related to the strengths of the multipole shape of the distribution of nucleons in the transverse plane, ε n : v n ∝ε n . The ε n , for n=1, 2, 3, 4, are sensitive to the shapes of the colliding ions, characterized by the quadrupole β 2 , octupole β 3 , and hexadecapole β 4 deformations. This sensitivity is investigated analytically and also in a Monte Carlo Glauber model. One observes a robust linear relation, $\langleε^{2}_{n}\rangle$ = $a^{'}_{n} + b^{'}_{n}β^{2}_{n}$, for events in a fixed centrality. The $\langleε^{2}_{1}\rangle$ has a contribution from β 3 and β 4 , and $\langleε^{2}_{3}\rangle$ from β 4 . In ultracentral collisions, there are little cross contributions between β 2 and ε 3 and between β 3 and ε 2 , but clear cross contributions are present in noncentral collisions. Additionally, $\langleε^{2}_{n}\rangle$ are insensitive to nonaxial shape parameters such as the triaxiality. This is good news because the measurements of v 2 , v 3 , and v 4 can be used to constrain simultaneously the β 2 , β 3 , and β 4 values. This is best done by comparing two colliding ions with similar mass numbers and therefore nearly identical $a^{'}_{n}$, to obtain a simple equation that relates the β n of the two species. Finally, this opens up the possibility to map the shape of the atomic nuclei at a timescale (<10 –24 s) much shorter than probed by low-energy nuclear structure physics (<10 –21 s), which ultimately may provide information complementary to that obtained in the nuclear structure experiments.
The track structure model of heavy ion cross sections was developed by Katz and co-workers in the 1960s. In this model the action cross section is evaluated by mapping the dose-response of a detector to gamma rays (modeled from biological target theory) onto the radial dose distribution from delta rays about the path of the ion. This is taken to yield the radial distribution of probability for a "hit" (an interaction leading to an observable end-point). Radial integration of the probability yields the cross section. When different response from ions of different Z having the same stopping power is observed this model may be indicated. Since the 1960s there have been several developments in the computation of the radial dose distribution, in the measurement of these distributions, and in new radiobiological data against which to test the model. The earliest model, by Butts and Katz made use of simplified delta ray distribution functions, of simplified electron range-energy relations, and neglected angular distributions. Nevertheless it made possible the calculation of cross sections for the inactivation of enzymes and viruses, and allowed extension to tracks in nuclear emulsions and other detectors and to biological cells. It set the pattern for models of observable effects in the matter through which the ion passed. Here we outline subsequent calculations of radial dose which make use of improved knowledge of the electron emission spectrum, the electron range-energy relation, the angular distribution, and some considerations of molecular excitation, of particular interest both close to the path of the ion and the outer limits of electron penetration. These are applied to the modeling of action cross sections for the inactivation of several strains of E-coli and B. subtilis spores where extensive measurements in the "thin-down" region have been made with heavy ion beams. Such calculations serve to test the radial dose calculations at the outer limit of electron penetration. We lack data from which to test these calculations in regions close to the path of the ion aside from our earliest work on latent tracks in plastics, though it appears that the criterion then suggested for the threshold of track formation, of a minimal dose at a minimal distance (of about 20 angstroms, in plastics), remains valid.
Observation of electromagnetic waves in the low-altitude auroral zone at frequencies between the proton and helium gyrofrequencies suggests that Alfven-ion cyclotron waves modified by the presence of helium ions are being excited. Estimates of the growth rates for this mode indicate that the auroral electron beam can provide the free energy for the instability. The effect of the heavy ions is to decrease the group velocity of the waves, leading to larger convective growth. Theoretical wave spectra are computed in the local approximation, which assumes that the gradient scale lengths in density and magnetic field are constant over the ray paths. Narrow banded spectral peaks similar to observations may be produced when the thickness of the electron beam is small (200 m at 3000 km altitude). Narrow beams also limit growth of whistler mode waves, which compete for the free energy of the electron beam.
The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. Furthermore, we confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this framework, it is possible to separate the effects of early time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as 𝑝−𝑝, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable.
The Coulomb dissociation contributions to fragmentation cross sections in relativistic heavy ion collisions, where more than one nucleon is removed, are estimated using the Weizsacker-Williams method of virtual quanta. Photonuclear cross sections taken from experimental results were used to fold into target photon number spectra calculated with the Weizsacker-Williams method. Calculations for several projectile target combinations over a wide range of charge numbers, and a wide range of incident projectile energies, are reported. These results suggest that multiple nucleon knockout by the Coulomb field may be of negligible importance in galactic heavy ion studies for projectiles lighter than Fe-56.
Proton data from the GOES 6 and 7 satellites and heavy ion data from the IMP-8 satellite have been compared to the expected results of Nymmik's new model for solar particle event fluences. This model calculates the energy spectra of Ions for protons through nickel for solar particle events, based upon the observed proton integral fluence above 30 MeV. Based upon 27 observed proton events of solar cycle 22, and three large historical events, with integral fluences above 30 MeV of greater than 106 particles/cm2, a reasonable agreement with model predictions is seen for more than half of the events. However, several events show a marked departure from the model predictions, leading to the conclusion that there may exist more than a single class of event, or that it may be necessary to include additional parameters within the model, such as solar disk position of the source flare, or height of disturbance in the solar corona. Data for heavy Ions, (oxygen and iron), were limited to a total of six solar particle events, of which only two occurred in solar cycle 22. The agreement between data and the model predictions appeared to be quite good, however this agreement was sensitively dependent upon the value taken for the proton fluence above 30 MeV.
It has been proposed that strong electromagnetic fields produced in the early stages of heavy-ion collisions can lead to splitting of the rapidity-odd directed flow of positive and negative hadrons. For light hadrons, the interpretation of such measurements is complicated by the low magnitude of directed flow as well as by ambiguities arising from transported quarks. To overcome these complications, we propose measurements using only hadrons carrying produced quarks (u¯, d¯, s s¯ ). Here, we discuss how to identify the kinematics where such hadrons are produced via the coalescence mechanism and therefore their flow is the sum of the flow of their constituent quarks. With this sum rule verified for certain combinations of hadrons, the expected systematic violation of this rule with increasing electric charge can be measured, which could be a consequence of the electromagnetic fields produced in the collisions. Our approach can be tested with the high statistics data from Phase II of the Beam Energy Scan program at the Relativistic Heavy Ion Collider.
Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.
Polarized quarks and antiquarks in high-energy heavy-ion collisions can lead to the spin alignment of vector mesons formed by quark coalescence. Here, using the relativistic spin Boltzmann equation for vector mesons derived from Kadanoff-Baym equations with an effective quark-meson model for strong interaction and quark coalescence model for hadronizaton, we calculate the spin density matrix element $ρ_{00}$ for Φ mesons and show that anisotropies of local field correlations with respect to the spin quantization direction lead to Φ meson's spin alignment. We propose that the local correlation or fluctuation of Φ fields is the dominant mechanism for the observed Φ meson's spin alignment and its strength can be extracted from experimental data as functions of collision energies. The calculated transverse momentum dependence of $ρ_{00}$ agrees with STAR's data. We further predict the azimuthal angle dependence of $ρ_{00}$ which can be tested in future experiments.
BACKGROUND: There is no human epidemiology data available to estimate the heavy ion cancer risks experienced by astronauts in space. Studies of tumor induction in mice are a necessary step to estimate risks to astronauts. Previous experimental data can be better utilized to model dose response for heavy ion tumorigenesis and plan future low dose studies. DOSE RESPONSE MODELS: The Harderian Gland data of Alpen et al.[1-3] was re-analyzed [4] using non-linear least square regression. The data set measured the induction of Harderian gland tumors in mice by high-energy protons, helium, neon, iron, niobium and lanthanum with LET s ranging from 0.4 to 950 keV/micron. We were able to strengthen the individual ion models by combining data for all ions into a model that relates both radiation dose and LET for the ion to tumor prevalence. We compared models based on Targeted Effects (TE) to one motivated by Non-targeted Effects (NTE) that included a bystander term that increased tumor induction at low doses non-linearly. When comparing fitted models to the experimental data, we considered the adjusted R2, the Akaike Information Criteria (AIC), and the Bayesian Information Criteria (BIC) to test for Goodness of fit.In the adjusted R2test, the model with the highest R2values provides a better fit to the available data. In the AIC and BIC tests, the model with the smaller values of the summary value provides the better fit. The non-linear NTE models fit the combined data better than the TE models that are linear at low doses. We evaluated the differences in the relative biological effectiveness (RBE) and found the NTE model provides a higher RBE at low dose compared to the TE model. POWER ANALYSIS: The final NTE model estimates were used to simulate example data to consider the design of new experiments to detect NTE at low dose for validation. Power and sample sizes were calculated for a variety of radiation qualities including some not considered in the Harderian Gland data set and with different background tumor incidences. We considered different experimental designs with varying number of doses and varying low doses dependant on the LET of the radiation. The optimal design to detect a NTE for an individual ion had 4 doses equally spaced below a maximal dose where bending due to cell sterilization was < 2%. For example at 100 keV/micron we would irradiate at 0.03 Gy, 0.065 Gy, 0.13 Gy, and 0.26 Gy and require 850 mice including a control dose for a sensitivity to detect NTE with 80% power. Sample sizes could be improved by combining ions similar to the methods used with the Harderian Gland data.
Baryon numbers are theorized to be carried by valence quarks in the standard QCD picture of the baryon structure. Another theory proposed an alternative baryon number carrier, a non-perturbative Y-shaped configuration of the gluon field, called the baryon junction in the 1970s. However, neither of these theories has been verified experimentally. Recently, searching for the baryon junction by investigating the correlation of net-charge and net-baryon yields at midrapidity in heavy-ion collisions has been suggested. Here, this paper presents studies of such correlations in collisions of various heavy ions from oxygen to uranium with the UrQMD Monte Carlo model. The UrQMD model implements valence quark transport as the primary means of charge and baryon stopping at midrapidity. Detailed studies are also conducted for isobaric $_{40}^{96}{\rm{Zr}}$ + $_{40}^{96}{\rm{Zr}}$ and $_{44}^{96}{\rm{Ru}}$ + $_{44}^{96}{\rm{Ru}}$ collisions. We found a universal trend of charge stopping with respect to baryon stopping and discovered that the charge stopping is always greater than the baryon stopping. This study provides a model baseline in valence quark transport for what is expected in net-charge and net-baryon yields at the midrapidity of relativistic heavy-ion collisions.
In this work, we study the flavor dependence of jet quenching in heavy-ion collisions. A next-to-leading-order perturbative QCD framework is applied to account for both quark and gluon contributions to light and heavy flavor hadron productions. A linear Boltzmann transport model coupled to hydrodynamic simulation is utilized to study the nuclear modification of heavy and light flavor jets in the dynamically evolving quark-gluon plasma. Within our jet quenching framework, we obtain a nice description of the nuclear modification factors of charged hadrons, D mesons, B mesons and B–decayed D mesons over a wide range of transverse momentum (8-300 GeV). Our result shows that perturbative QCD is sufficient to explain the color, mass and energy dependence of parton energy loss and jet quenching in heavy-ion collisions.
Heavy-ion-induced degradation in the reverse leakage current of SiC Schottky power diodes exhibits a strong dependence on the ion angle of incidence. This effect is studied experimentally for several different bias voltages applied during heavy-ion exposure. In addition, TCAD simulations are used to give insight on the physical mechanisms involved.