Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Beam Simulations”

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.

At least 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Radiative Losses in the Chromosphere During A C-Class Flare

Context. Solar flares release an enormous amount of energy (~10(exp 32) erg) into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. Aims. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. Methods. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semiempirical model atmospheres, which were inferred from a multiline inversion of observed Stokes profiles obtained with the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. The radiative losses were computed taking into account the effect of partial redistribution and non-local thermodynamic equilibrium. To estimate the integrated radiative losses in the chromosphere, the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. We also compared our time series of radiative losses with those from the RADYN flare simulations. Results. We obtained a high spatial-resolution map of integrated radiative losses around the flare peak time. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H and K and Mgii h and k lines are strong and comparable to the Ca IR triplet (~32kW m(exp -2)). Since our flare is a relatively weak event, the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution (~7kW m(exp -2)) in the selected limits of the chromosphere. The temporal evolution of total integrated radiative losses exhibits sharply rising losses (0.4kW m(exp -2) (s(exp -1)) and a relatively slow decay (0.23kW m(exp -2) s(exp -1)). The maximum value of total radiative losses is reached around the flare peak time and can go up to 175kWm􀀀2 for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of 5 × 10(exp 10) erg s(exp -1) cm(exp -2).

Chromosphere↗
DOE OSTI2022

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2020

Feasibility of “Weak GPS” Real-Time Positioning/Timing at Lunar Distance

There are multiple Global Navigation Satellite Systems (GNSS’s), comprising over 100 navigation satellites in the Earth’s medium and high orbits. Most of these satellites have antennas that point Nadir to earth, and transmit navigation signals so vehicles on Earth’s surface and in its vicinity can perform trilateration and estimate its 3-dimenional (3D) positioning. The sidelobes of these antennas can occasionally point to the Moon. It is postulated that a lunar vehicle carrying a large enough receiving antenna can occasionally detects and receives four or more sidelobes of these weak GNSS signals, thus enabling the vehicle to perform 3D positioning using an onboard GNSS receiver. We propagate the orbits of the GNSS satellites from United States’ Global Positioning Satellite (GPS) constellation, the Europe’s Galileo constellation, and the Russia’s GLONASS constellation, a total of 81 satellites. We simulate the visibility of these satellites by a lunar vehicle in a Near Rectilinear Halo Orbit (NRHO), based on the assumption that the lunar vehicle is “in-view” of a GNSS satellite as long as it falls within the 40- degree beam-width of the satellite. We also simulate the 3D positioning performance as a function of satellites’ ephemeris errors and pseudo-range errors. The preliminary results show that the lunar vehicle can “see” 5 – 13 satellites, and achieve a 3D positioning error (one-sigma) of 200 – 300 meters based on reasonable ephemeris and pseudo- range error assumptions. We also consider the case of using relative positioning to mitigate the GNSS satellites’ ephemeris biases. That is, by assuming a reference receiver with accurately known positioning that is close to the lunar vehicle, and computing the relative position of the lunar vehicle relative to the reference.

Cheung, Kar-Ming↗
DOE OSTI2025

Enhanced Monte Carlo Simulations for Electron Energy Loss Mitigation in Real-Space Nanoimaging of Thick Biological Samples and Microchips

High-resolution imaging using Transmission Electron Microscopy (TEM) is essential for applications such as grain boundary analysis, microchip defect characterization, and biological imaging. However, TEM images are often compromised by electron energy spread and other factors. In TEM mode, where the objective and projector lenses are positioned downstream of the sample, electron–sample interactions cause energy loss, which adversely impacts image quality and resolution. This study introduces a simulation tool to estimate the electron energy loss spectrum (EELS) as a function of sample thickness, covering electron beam energies from 300 keV to 3 MeV. Leveraging recent advances in MeV-TEM/STEM technology, which includes a state-of-the-art electron source with 2-picometer emittance, an energy spread of 3 × 10 -5 , and optimized beam characteristics, we aim to minimize energy spread. By integrating EELS capabilities into the BNL Monte Carlo (MC) simulation code for thicker samples, we evaluate electron beam parameters to mitigate energy spread resulting from electron–sample interactions. Based on our simulations, we propose an experimental procedure for quantitively distinguishing between elastic and inelastic scattering. The findings will guide the selection of optimal beam settings, thereby enhancing resolution for nanoimaging of thick biological samples and microchips.

36 MATERIALS SCIENCE↗
NASA NTRS2019

NASA’s Galactic Cosmic Ray Simulator at Brookhaven National Laboratory: Enabling Human Exploration Missions to the Moon and Mars

With exciting new Agency plans for a sustainable return to the moon, astronauts will once again leave earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic rays (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPE’s with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever-penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include epithelial carcinogenesis and leukemia, central nervous system effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including cataracts, circulatory and heart disease, as well as, potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy ion exposures and the validation of countermeasures in a relevant space environment. NASA has developed the “GCR Simulator” at Brookhaven National Laboratory to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (~65-75%) and alpha particles (~10-20%) with heavier ions (Z≤3) contributing the remainder. The “GCR Simulator” exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the five heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires 75-90 minutes. To more closely simulate the low dose rates found in space, sequential field exposures can be divided into daily fractions over 2-4 weeks, with individual fractions as low as 0.1-0.2 mGy. In the large beam configuration (60 x 60 cm(exp 2)), 54 special housing cages can accommodate 2-3 mice each for a 70-75 min duration or ~15 individually housed rats. Emerging research results from our 2018 runs utilizing mixed heavy ion fields and protracted space exposures are forthcoming and deepen our understanding of the numerous health risks faced by our astronauts. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NASA Space Radiation Laboratory to enable future exploration missions.

Lisa C Simonsen↗
DOE OSTI2023

Atomistic simulations for investigation of substrate effects on lipid in-source fragmentation in secondary ion mass spectrometry

In beam-based ionization methods, the substrate plays an important role on the desorption mechanism of molecules from surfaces. Both the specific orientation that a molecule adopts at a surface and the strength of the molecule-surface interaction can greatly influence desorption processes, which in turn will affect the ion yield and the degree of in-source fragmentation of a molecule. In the beam-based method of secondary ion mass spectrometry (SIMS), in-source fragmentation can be significant and molecule specific due to the hard ionization method of using a primary ion beam for molecule desorption. To investigate the role of the substrate on orientation and in-source fragmentation, we have used atomistic simulations—molecular dynamics in combination with density functional theory calculations—to explore the desorption of a sphingolipid (palmitoylsphingomyelin) from a model surface (gold). We then compare SIMS data from this model system to our modeling findings. Using this approach, we found that the combined adsorption and binding energy of certain bonds associated with the headgroup fragments (C3H8N+, C5H12N+, C5H14NO+, and C5H15PNO4+) was a good predictor for fragment intensities (as indicated by relative ion yields). This is the first example where atomistic simulations have been applied in beam-based ionization of lipids, and it presents a new approach to study biointerfacial lipid ordering effects on SIMS imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Symplectic Particle Tracking in a Thick Nonlinear McMillan Lens for the Fermilab Integrable Optics Test Accelerator (IOTA)

The McMillan system is a novel method to increase the tune spread of a beam without decreasing its dynamic aperture due to the systems integrability. While the ideal system is based on an infinitely thin kick, the physical design requires a thick electron lens, including a solenoid. Particle transport through the lens is difficult to simulate due to the nature of the force on the circulating beam. This paper demonstrates accurate simulation of a thick McMillan lens in a solenoid using symplectic integrators derived from Yoshidas method.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Modeling of emittance growth due to Coulomb collisions in plasma-based accelerators

Coulomb collisions with background plasma are one source of emittance degradation in plasma accelerators. This work shows that the emittance growth due to Coulomb collisions can be correctly captured in particle-in-cell simulations, with a proper Monte Carlo binary collision module. The theory of the emittance growth due to Coulomb collisions is extended from a monoenergetic matched beam to a mismatched beam with energy spread and is compared with simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Facile Tensile Testing Platform for In Situ Transmission Electron Microscopy of Nanomaterials

In situ tensile testing using transmission electron microscopy (TEM) is a powerful technique to probe structure-property relationships of materials at the atomic scale. In this work, a facile tensile testing platform for in situ characterization of materials inside a transmission electron microscope is demonstrated. The platform consists of: 1) a commercially available, flexible, electron-transparent substrate (e.g., TEM grid) integrated with a conventional tensile testing holder, and 2) a finite element simulation providing quantification of specimen-applied strain. The flexible substrate (carbon support film of the TEM grid) mitigates strain concentrations usually found in free-standing films and enables in situ straining experiments to be performed on materials that cannot undergo localized thinning or focused ion beam lift-out. The finite element simulation enables direct correlation of holder displacement with sample strain, providing upper and lower bounds of expected strain across the substrate. The tensile testing platform is validated for three disparate material systems: sputtered gold-palladium, few-layer transferred tungsten disulfide, and electrodeposited lithium, by measuring lattice strain from experimentally recorded electron diffraction data. The results show good agreement between experiment and simulation, providing confidence in the ability to transfer strain from holder to sample and relate TEM crystal structural observations with material mechanical properties.

2D materials↗
DOE OSTI2024

Performance of the Fermilab Linac Injector

The Fermilab Linac injection line consists of a 35 keV magnetron-type H- ion source, two-solenoid Low Energy Beam Transport (LEBT), 201 MHz 4-rod 750 keV Radio Frequency Quadrupole (RFQ), and a Medium Energy Transport (MEBT) containing 4 quadrupoles and a bunching cavity. The injector delivers 25 mA, 48 s pulses to drift-tube Linac at a repetition rate of 15Hz. The transmission efficiency has been lower than expected since commissioning. Recent beam current measurements suggest that the beam is primarily lost upstream of the RFQ exit. Numerical simulations indicate that ions passing through the non-linear field region of the solenoids could produce a beam with an increased emittance resulting in up to 50 % of the LEBT beam current failing to meet the RFQ acceptance. An aperture restriction was installed upstream of the first solenoid to remove these ions. This report describes the results of measurements and simulations as well as the LEBT tuning.

Jones, Daniel C.↗
DOE OSTI2024

Performance of the Fermilab Linac Injector

The Fermilab Linac injection line consists of a 35 keV magnetron-type H- ion source, two-solenoid Low Energy Beam Transport (LEBT), 201 MHz 4-rod 750 keV Radio Frequency Quadrupole (RFQ), and a Medium Energy Transport (MEBT) containing 4 quadrupoles and a bunching cavity. The injector delivers 25 mA, 48 s pulses to drift-tube Linac at a repetition rate of 15Hz. The transmission efficiency has been lower than expected since commissioning. Recent beam current measurements suggest that the beam is primarily lost upstream of the RFQ exit. Numerical simulations indicate that ions passing through the non-linear field region of the solenoids could produce a beam with an increased emittance resulting in up to 50 % of the LEBT beam current failing to meet the RFQ acceptance. An aperture restriction was installed upstream of the first solenoid to remove these ions. This report describes the results of measurements and simulations as well as the LEBT tuning.

Jones, Daniel C.↗
DOE OSTI2024

Performance of the Fermilab Linac Injector

The Fermilab Linac injection line consists of a 35 keV magnetron-type H\textsuperscript{-} ion source, two-solenoid Low Energy Beam Transport (LEBT), 201 MHz 4-rod 750 keV Radio Frequency Quadrupole (RFQ), and a Medium Energy Transport (MEBT) containing 4 quadrupoles and a bunching cavity. The injector delivers 25 mA, 48 µs pulses to the drift-tube Linac at a repetition rate of 15Hz. The transmission efficiency of the injector has been lower than expected since commissioning. Recent beam current measurements suggest that the beam is primarily lost upstream of the RFQ exit. Numerical simulations indicate that ions passing through the non-linear field region of the solenoids could produce a beam with an increased emittance resulting in up to 50\% of the LEBT beam current failing to meet the RFQ acceptance. An aperture restriction was installed upstream of the first solenoid to remove these ions. This report describes the results of measurements and simulations as well as the LEBT tuning.

43 PARTICLE ACCELERATORS↗
NASA NTRS1992

Lidar multiple scattering: Dependence on atmospheric parameters

Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).

Wiegner, Matthias↗
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 19

The topics include: 1) The abundances of Iron-60 in Pyroxene Chondrules from Unequilibrated Ordinary Chondrites; 2) LL-Ordinary Chondrite Impact on the Moon: Results from the 3.9 Ga Impact Melt at the Landing Site of Appolo 17; 3) Evaluation of Chemical Methods for Projectile Identification in Terrestrial and Lunar Impactites; 4) Impact Cratering Experiments in Microgravity Environment; 5) New Achondrites with High-Calcium Pyroxene and Its implication for Igneous Differentiation of Asteroids; 6) Climate History of the Polar Regions of Mars Deduced form Geologic Mapping Results; 7) The crater Production Function for Mars: A-2 Cumulative Power-Law Slope for Pristine Craters Greater than 5 km in Diameter Based on Crater Distribution for Northern Plains Materials; 8) High Resolution Al-26 Chronology: Resolved Time Interval Between Rim and Interior of a Highly Fractionated Compact Type a CAI from Efremovka; 9) Assessing Aqueous Alteration on Mars Using Global Distributions of K and Th; 10) FeNi Metal Grains in LaPaz Mare Basalt Meteorites and Appolo 12 Basalts; 11) Unique Properties of Lunar Soil for In Situ Resource Utilization on the Moon; 12) U-Pb Systematics of Phosphates in Nakhlites; 13) Measurements of Sound Speed in Granular Materials Simulated Regolith; 14) The Effects of Oxygen, Sulphur and Silicon on the Dihedral Angles Between Fe-rich Liquid Metal and Olivine, Ringwoodite and Silicate Perovskite: Implications for Planetary Core Formation; 15) Seismic Shaking Removal of Craters 0.2-0.5 km in Diameter on Asteroid 433 Eros; 16) Focused Ion Beam Microscoopy of ALH84001 Carbonate Disks; 17) Simulating Micro-Gravity in the Laboratory; 18) Mars Atmospheric Sample Return Instrument Development; 19) Combined Remote LIBS and Raman Spectroscopy Measurements; 20) Unusual Radar Backscatter Properties Along the Northern Rim of Imbrium Basin; 21) The Mars Express/NASAS Project at JPL; 22) The Geology of the Viking 2 Lander Site Revisited; 23) An Impact Genesis for Loki Patera? 24) Mars Polar Cap Edges Tracked over 3 Full Mars Years; 25) Elemental Abundance in Presolar SiC: Comparing Grains Separated by Acid Residue and Gently Separation Procedures; 26) First Results from the Descent Imager/Spectral Radiometer (DISR) Experiment on the Huygens Entry Probe of Titan; 27) Minor Element Behavior of Pallasite Olivine: Understanding Pallasite Thermal History and Chronology; 28) Canonical Anorthite in a Grosnaja Forsterite-bearing CAI; 29) Experimental Evidence for Condensation of 'Astrophysical' Carbonate; 30) Distribution and Classification of Multiple Coronae on Venus; 31) Recognition of Rayed Craters on Mars in THEMIS Thermal Infrared Imagery: Implications for Martian Meteorite Source Regions; 32) Geochemical Modeling of Evaporites on Mars: Insight from Meridiani Planum; 33) Hadean Crustal Processes Revealed from Oxygen Isotopes and U-Th-Pb Depth Profiling of Pre-4.0 Ga Detrital Zircons from Western Australia; 34) On Modeling the Seepage of Water into the Martian Subsurface; 35) Martial Gullies and Groundwater: A Series of Unfortunate Exceptions; 36) Olivine and Carbonate Globules in ALH84001: A Terrestrial Analog, and Implications for Water on Mars; 37) A Reevaluation of Mass Movements Within the Valles Marineris Region of Mars Using MOLA and MOC Data; 38) Evidence of Hydrated 109P/Swift-Tuttle Meteoroids from Meteor Spectroscopy; 39) Cr-54 Anomalies in the Solar System: Their Extent and Origin; 40) Reevaluation of the Mn-53-Cr-53 Systematic in the Basaltic Achondrites; 41) Effective Liquid Metal-Silicate Mixing Upon Shock by Power-Law Droplet Size Scaling in Richtmyer-Meshkov Like Perturbations; 42) Post-Impact Deformation of Impact Craters: Towards a Better Understanding Through the Study of Mjolnir Crater; 43) Cutting Silica Aerogel for Particle Extraction; 44) Liquid Hydrocarbons on Titan's Surface? How Cassini ISS Observations Fit into the Story (So Far); and 45) Mesoscale Simulations of Polar Circulations: Late Spring to Late Summe

Source record↗
DOE OSTI2024

Optimize Electron Beam Energy toward In Situ Imaging of Thick Frozen Bio-Samples with Nanometer Resolution Using MeV-STEM

To optimize electron energy for in situ imaging of large biological samples up to 10 μm in thickness with nanoscale resolutions, we implemented an analytical model based on elastic and inelastic characteristic angles. This model has been benchmarked by Monte Carlo simulations and can be used to predict the transverse beam size broadening as a function of electron energy while the probe beam traverses through the sample. As a result, the optimal choice of the electron beam energy can be realized. In addition, the impact of the dose-limited resolution was analysed. While the sample thickness is less than 10 μm, there exists an optimal electron beam energy below 10 MeV regarding a specific sample thickness. However, for samples thicker than 10 μm, the optimal beam energy is 10 MeV or higher depending on the sample thickness, and the ultimate resolution could become worse with the increase in the sample thickness. Moreover, a MeV-STEM column based on a two-stage lens system can be applied to reduce the beam size from one micron at aperture to one nanometre at the sample with the energy tuning range from 3 to 10 MeV. In conjunction with the state-of-the-art ultralow emittance electron source that we recently implemented, the maximum size of an electron beam when it traverses through an up to 10 μm thick bio-sample can be kept less than 10 nm . This is a critical step toward the in situ imaging of large, thick biological samples with nanometer resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of the Role of the Role of Nuclear Physics in Heavy Element Nucleosynthesis, through the Study of Key Reactions, and the Improvement of Theoretical Reaction Rates (Final Report)

Neutrino-driven winds in core-collapse supernovae have been identified as a possible site for the production of elements heavier than iron. Traditionally, these neutrino-driven winds have been proposed as the site of the main r-process. Recent simulations fail to reproduce the conditions required for the main r-process. while they remain a promising site for producing the lightest elements beyond iron, e.g., Sr, Y, and Zr through the νp process. The efficiency of the νp process depends on the hydrodynamical conditions, the electron fraction (which is related to the neutrino properties), and the nuclear reactions on many short-lived nuclei with limited (if at all) experimental information. The reaction rates on these nuclei are based on theoretical predictions using the Hauser-Feshbach model. Recent sensitivity studies have highlighted the importance of neutron-induced reactions on these nuclei along the νp process path. This work aimed to experimentally constrain reaction rates that are known to play a key role in the neutrino-p process nucleosynthesis. A secondary subsequently-added objective was to start the implementation of techniques that improve the description of nuclear properties in the Hauser-Feshbach model by extending the microscopic nuclear level density description offered via the shell model to high excitation energies without using experiment-based renormalizations. The main objective of this work was the experimental constraint of the 56 Ni(n,p) 56 Co reaction rate via a measurement of the inverse reaction 56 Co(p,n) 56 Ni at the National Superconducting Cyclotron Laboratory (NSCL) and later the Facility for Rare Isotope Beams (FRIB). This reaction is considered the key one for determining the yields possible by the neutrino-p process. A technique for this type of measurement in inverse kinematics at low energies did not exist before this work. The work also had two secondary objectives. First, to contribute to efforts to measure the same reaction in direct kinematics using a radioactive target at Los Alamos National Laboratory (LANL), and second, to advance work to implement shell-model-deduced microscopic level densities in Hauser-Feshbach calculations. The project has resulted in the development of the first technique to perform (p,n) cross-section measurements in relevant-for-astrophysics low energies in inverse kinematics using a magnetic spectrometer or separator, and neutron detectors for neutron tagging. It has also resulted in the precise measurement of the cross-section of the 40 Ar(p,n) 40 K reaction in a proof-of-principle experiment realized by using a beam-line quadrupole of the ReA3 accelerator of NSCL/FRIB. As part of this project the technique was successfully adapted to make use of the superior acceptance of the Separator for Capture Reactions (SECAR) at FRIB. In this project, the required experimental setup simulations and beam optics were developed and tested with the measurement of the 58 Fe(p,n) 58 Cu reaction cross-section. Additionally, this project contributed with simulation work to the development of a technique to measure (n,p) reactions with radioactive targets at LANL, and which resulted in the measurement of the key 56 Ni(n,p) 56 Co reaction cross-section at neutron energies above ≈1 MeV. Last, the project initiated work in the development of shell model based level densities using the moments method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗