Search NASA⌕ Search

Engineering topics

Surman, R.

Publications and source records attributed to Surman, R..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Uncertainty quantification of mass models using ensemble Bayesian model averaging

Developments in the description of the masses of atomic nuclei have led to various nuclear mass models that provide predictions for masses across the whole chart of nuclides. These mass models play an important role in understanding the synthesis of heavy elements in the rapid neutron capture ( r ) process. However, it is still a challenging task to estimate the size of uncertainty associated with the predictions of each mass model. In this work, a method called ensemble Bayesian model averaging (EBMA) is introduced to quantify the uncertainty of one-neutron separation energies (S 1 n ) which are directly relevant in the calculations of r -process observables. Here, this Bayesian method provides a natural way to perform model averaging, selection, and uncertainty quantification, by combining the mass models as a mixture of normal distributions whose parameters are optimized against the experimental data, employing the Markov chain Monte Carlo method using the no-u-turn sampler. The EBMA model optimized with all the experimental S 1 n from the AME2003 nuclides are shown to provide reliable uncertainty estimates when tested with the new data in the AME2020.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measuring the β-decay Properties of Neutron-rich Exotic Pm, Sm, Eu, and Gd Isotopes to Constrain the Nucleosynthesis Yields in the Rare-earth Region

Abstract The β -delayed neutron-emission probabilities of 28 exotic neutron-rich isotopes of Pm, Sm, Eu, and Gd were measured for the first time at RIKEN Nishina Center using the Advanced Implantation Detector Array (AIDA) and the BRIKEN neutron detector array. The existing β -decay half-life ( T 1/2 ) database was significantly increased toward more neutron-rich isotopes, and uncertainties for previously measured values were decreased. The new data not only constrain the theoretical predictions of half-lives and β -delayed neutron-emission probabilities, but also allow for probing the mechanisms of formation of the high-mass wing of the rare-earth peak located at A ≈ 160 in the r -process abundance distribution through astrophysical reaction network calculations. An uncertainty quantification of the calculated abundance patterns with the new data shows a reduction of the uncertainty in the rare-earth peak region. The newly introduced variance-based sensitivity analysis method offers valuable insight into the influence of important nuclear physics inputs on the calculated abundance patterns. The analysis has identified the half-lives of 168 Sm and of several gadolinium isotopes as some of the key variables among the current experimental data to understand the remaining abundance uncertainty at A = 167–172.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Searching for the origin of the rare-earth peak with precision mass measurements across Ce–Eu isotopic chains

A nuclear mass survey of rare-earth isotopes has been conducted with the Canadian Penning Trap mass spectrometer using the most neutron-rich nuclei thus far extracted from the CARIBU facility. Here, we present a collection of 12 nuclear masses determined with a precision of ≤ 10 keV/c 2 for Z = 58-63 nuclei near N = 100. Independently, a detailed study exploring the role of nuclear masses in the formation of the r-process rare-earth abundance peak has been performed. Employing a Markov chain Monte Carlo (MCMC) technique, mass predictions of lanthanide isotopes have been made which uniquely reproduce the observed solar abundances near A = 164 under three distinct astrophysical outflow conditions. We demonstrate that the mass surface trends thus far mapped out by our measurements are most consistent with MCMC mass predictions given an r process that forms the rare-earth peak during an extended (n, γ) ⇌ (γ, n) equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Kilonovae Across the Nuclear Physics Landscape: The Impact of Nuclear Physics Uncertainties on r-process-powered Emission

Merging neutron stars produce “kilonovae”—electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their luminosities, with peak brightness varying by more than an order of magnitude. We also find variability in the shape of the kilonova light curves and their color, which in some cases runs counter to the expectation that increasing levels of lanthanide and/or actinide enrichment will be correlated with longer, dimmer, redder emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Following nuclei through nucleosynthesis: A novel tracing technique

Astrophysical nucleosynthesis is a family of diverse processes by which atomic nuclei undergo nuclear reactions and decay to form new nuclei. The complex nature of nucleosynthesis, which can involve as many as tens of thousands of interactions between thousands of nuclei, makes it difficult to study any one of these interactions in isolation using standard approaches. In this work, we present a new technique, nucleosynthesis tracing, that we use to quantify the relative fraction of nuclear abundances that pass through individual nuclear reaction, decay, and fission processes at any point during nucleosynthesis. Here, we apply this technique to study fission and β – decay as they occur in the rapid neutron capture (r) process of nucleosynthesis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing β-decay, α-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Following Fission Products in Explosive Astrophysical Environments

The astrophysical process by which the heaviest elements are formed in the universe is known as the rapid neutron capture process, or r process, of nucleosynthesis. The r process is characterized by the neutron capture and β - decay of short-lived, neutron-rich atomic nuclei; in suitably extreme environments, nuclear fission can also play a major role in determining the ensuing nucleosynthesis. In this work, we present the application of our recently developed nucleosynthesis tracing framework to precisely quantify the impact that neutron-induced and β - -delayed fission processes have in r-process environments that produce fissioning nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Propagation of statistical uncertainties of Skyrme mass models to simulations of r -process nucleosynthesis

Uncertainties in nuclear models have a major impact on simulations that aim at understanding the origin of heavy elements in the universe through the rapid neutron capture process ( r process) of nucleosynthesis. Within the framework of the nuclear density functional theory, we use results of Bayesian statistical analysis to propagate uncertainties in the parameters of energy density functionals to the predicted r -process abundance pattern, by way not only of the nuclear masses but also through the influence of the masses on β -decay and neutron capture rates. We point out the importance of the nonequilibrium end stage of the r process in determining the width of the resulting abundance pattern uncertainty bands. We additionally make the first identifications of specific parameters of Skyrme-like energy density functionals which show tentative correlations with particular aspects of the r -process abundance pattern. While previous studies have explored the reduction in the abundance pattern uncertainties due to anticipated new measurements of neutron-rich nuclei, here we point out that an even larger reduction will occur when these new measurements are used to reduce the uncertainty of model predictions of masses, which are then propagated through to the abundance pattern. We make a quantitative prediction for how large this reduction will be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Exploring the mass surface near the rare-earth abundance peak via precision mass measurements at JYFLTRAP

The JYFLTRAP double Penning trap at the Ion Guide Isotope Separator On-Line facility has been used to measure the atomic masses of 13 neutron-rich rare-earth isotopes. Eight of the nuclides, 161 Pm, 163 Sm, 164,165 Eu, 167 Gd, and 165,167,168 Tb, were measured for the first time. The systematics of the mass surface has been studied via one- and two-neutron separation energies as well as neutron pairing-gap and shell-gap energies. The proton-neutron pairing strength has also been investigated. Furthermore, the impact of the new mass values on the astrophysical rapid neutron capture process has been studied. The calculated abundance distribution results in a better agreement with the solar abundance pattern near the top of the rare-earth abundance peak at around A≈165.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Strong one-neutron emission from two-neutron unbound states in β decays of the r -process nuclei Ga 86 , 87

β-delayed one-neutron and two-neutron branching ratios (P$_{1n}$ and P$_{2n}$) have been measured in the decay of A = 84 to 87 Ga isotopes at the Radioactive-Isotope Beam Factory (RIBF) at the RIKEN Nishina Center using a high-efficiency array of 3 He neutron counters (BRIKEN). Two-neutron emission was observed in the decay of 84, 85, 87 Ga for the first time and the branching ratios were measured to be P$_{2n}$ = 1.6 (2)%, 1.3 (2)%, and 10.2 ( 28 ) stat ( 5 ) sys %, respectively. One-neutron branching ratio of 87 Ga (P$_{1n}$= 81 (9) stat (8) sys %) and half-life of 29(4) ms were measured for the first time. The branching ratios of 86 Ga were also measured to be P$_{1n}$ = 74 (2) stat (8) sys % and 16.2 (9) stat ( 6 ) sys % with better precision than a previous study. The observation that P$_{1n}$ > P$_{2n}$ for both 86,87 Ga was unexpected and is interpreted as a signature of dominating one-neutron emission from the two-neutron unbound excited states in 86, 87 Ge. In order to interpret the experimental results, shell-model and Hauser-Feshbach statistical model calculations of delayed particle and γ-ray emission probabilities were performed. This model framework reproduces the experimental results. The shell model alone predicts P$_{2n}$ significantly larger than P$_{1n}$ for the 87 Ga decay, and it is necessary to invoke a statistical description to successfully explain the observation that P$_{1n}$ > P$_{2n}$. Our new results demonstrate the relevance and importance of a statistical description of neutron emission for the prediction of the decay properties of multineutron emitters and that it must be included in the r -process modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Using excitation-energy dependent fission yields to identify key fissioning nuclei in r -process nucleosynthesis

The possibility that nucleosynthesis in neutron star mergers may reach fissioning nuclei introduces significant uncertainties in predicting the relative abundances of $r$-process material from such events. We evaluate the impact of using sets of fission yields given by the 2016 GEF code for spontaneous (sf), neutron-induced $((n, f))$, and $β$-delayed (βdf) fission processes which take into account the approximate initial excitation energy of the fissioning compound nucleus. We further explore energy-dependent fission dynamics in the $r$ process by considering the sensitivity of our results to the treatment of the energy sharing and de-excitation of the fission fragments using the FREYA code. We show that the asymmetric-to-symmetric yield trends predicted by GEF 2016 can reproduce the high-mass edge of the second $r$-process peak seen in solar data and examine the sensitivity of this result to the mass model and astrophysical conditions applied. We consider the effect of fission yields and barrier heights on the nuclear heating rates used to predict kilonova light curves. We find that fission barriers influence the contribution of 254 Cf spontaneous fission to the heating at ~100 d, such that a light curve observation consistent with such late-time heating would both confirm that actinides were produced in the event and imply the fission barriers are relatively high along the 254 Cf β-feeding path. We lastly determine the key nuclei responsible for setting the r-process abundance pattern by averaging over thirty trajectories from a 1.2–1.4 M ⊙ neutron star merger simulation. Finally, we show it is largely the odd-N nuclei undergoing $(Z, N)(n, f)$ and $(Z, N)β$df that control the relative abundances near the second peak. We find the 'hot spots' for β-delayed and neutron-induced fission given all mass models considered and show most of these nuclei lie between the predicted $N$ = 184 shell closure and the location of currently available experimental decay data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

β-delayed Fission in r-process Nucleosynthesis

In this work, we present β-delayed neutron emission and β-delayed fission (βdf) calculations for heavy, neutron-rich nuclei using the coupled Quasi-Particle Random Phase Approximation plus Hauser-Feshbach (QRPA+HF) approach. From the initial population of a compound nucleus after β-decay, we follow the statistical decay, taking into account competition between neutrons, γ-rays, and fission. We find a region of the chart of nuclides where the probability of βdf is ~100%, which likely prevents the production of superheavy elements in nature. For a subset of nuclei near the neutron dripline, neutron multiplicity and the probability of fission are both large, leading to the intriguing possibility of multi-chance βdf, a decay mode for extremely neutron-rich heavy nuclei. In this decay mode, β-decay can be followed by multiple neutron emission, leading to subsequent daughter generations that each have a probability to fission. We explore the impact of βdf in rapid neutron-capture process (r-process) nucleosynthesis in the tidal ejecta of a neutron star–neutron star merger and show that it is a key fission channel that shapes the final abundances near the second r-process peak.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Californium-254 and Kilonova Light Curves

Neutron star mergers o er unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the r-process. We have performed dynamical nucleosynthesis calculations and identi ed a single isotope, 254 Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the e ciency of ssion thermalization compared to other nuclear channels. We estimate the ssion fragment yield of this nucleus and outline the astrophysical conditions under which 254 Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Precision Mass Measurements of Neutron-Rich Neodymium and Samarium Isotopes and Their Role in Understanding Rare-Earth Peak Formation

The Canadian Penning Trap mass spectrometer at the Californium Rare Isotope Breeder Upgrade (CARIBU) facility was used to measure the masses of eight neutron-rich isotopes of Nd and Sm. These measurements are the first to push into the region of nuclear masses relevant to the formation of the rare-earth abundance peak at A similar to 165 by the rapid neutron-capture process. We compare our results with theoretical predictions obtained from "reverse engineering" the mass surface that best reproduces the observed solar abundances in this region through a Markov chain Monte Carlo technique. Our measured masses are consistent with the reverse-engineering predictions for a neutron star merger wind scenario.

Orford, R.↗
DOE OSTI2017

Conceptual design of a hybrid neutron-gamma detector for study of β-delayed neutrons at the RIB facility of RIKEN

Here, the conceptual design of the BRIKEN neutron detector at the radioactive ion beam factory (RIBF) of the RIKEN Nishina Center is reported. The BRIKEN setup is a complex system aimed at detecting heavy-ion implants, β particles, γ rays and β-delayed neutrons. The whole setup includes the Advanced Implantation Detection Array (AIDA), two HPGe Clover detectors and up to 166 3 He-filled counters embedded in a high-density polyethylene moderator. The design is quite complex due to the large number and different types of 3 He-tubes involved and the additional constraints introduced by the ancillary detectors for charged particles and γ rays. This article reports on a novel methodology developed for the conceptual design and optimisation of the 3 He-counter array, aiming for the best possible performance in terms of neutron detection. The algorithm is based on a geometric representation of two selected detector parameters of merit, namely, the average neutron detection efficiency and the efficiency flatness as a function of a reduced number of geometric variables. The response of the neutron detector is obtained from a systematic Monte Carlo simulation implemented in GEANT4. The robustness of the algorithm allowed us to design a versatile detection system, which operated in hybrid mode includes the full neutron counter and two clover detectors for high-precision gamma spectroscopy. In addition, the system can be reconfigured into a compact mode by removing the clover detectors and re-arranging the 3 He tubes in order to maximize the neutron detection performance. Both operation modes shows a rather flat and high average efficiency. In summary, we have designed a system which shows an average efficiency for hybrid mode ( 3 He tubes + clovers) of 68.6% and 64% for neutron energies up to 1 and 5 MeV, respectively. For compact mode (only 3 He tubes), the average efficiency is 75.7% and 71% for neutron energies up to 1 and 5 MeV, respectively. The performance of the BRIKEN detection system has been also quantified by means of Monte Carlo simulations with different neutron energy distributions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Reverse engineering nuclear properties from rare earth abundances in the r process

The bulk of the rare earth elements are believed to be synthesized in the rapid neutron capture process or r process of nucleosynthesis. The solar r-process residuals show a small peak in the rare earths around $A\sim 160$, which is proposed to be formed dynamically during the end phase of the r process by a pileup of material. This abundance feature is of particular importance as it is sensitive to both the nuclear physics inputs and the astrophysical conditions of the main r process. Here, we explore the formation of the rare earth peak from the perspective of an inverse problem, using Monte Carlo studies of nuclear masses to investigate the unknown nuclear properties required to best match rare earth abundance sector of the solar isotopic residuals. When nuclear masses are changed, we recalculate the relevant β-decay properties and neutron capture rates in the rare earth region. The feedback provided by this observational constraint allows for the reverse engineering of nuclear properties far from stability where no experimental information exists. We investigate a range of astrophysical conditions with this method and show how these lead to different predictions in the nuclear properties influential to the formation of the rare earth peak. Finally, we conclude that targeted experimental campaigns in this region will help to resolve the type of conditions responsible for the production of the rare earth nuclei, and will provide new insights into the longstanding problem of the astrophysical site(s) of the r process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗