Shell Model Methods, Statistical Nuclear Reactions, and Beta-delayed Neutron Emission
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, implying a slope parameter L = 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Finally, combining astrophysical data with PREX-II and chiral effective field theory constraints yields S 0 = $33.0^{+2.0}_{-1.8}$ MeV , L = $53^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin} = 0.17^{+0.04}_{-0.04}$ fm.
Understanding the radiogenic neutron production rate through the (α, n) reaction is crucial in many areas of physics, including dark matter searches, neutrino studies, and nuclear astrophysics. In addition to its relevance for fundamental research, the (α, n) reaction also plays a significant role in nuclear energy technologies, for example by contributing to neutron production in subcritical systems using UO2, as well as in applications such as medical physics. This review examines the current state of (α, n) yield calculations and neutron spectra, describes the computational tools used for their estimation, and discusses the available cross-section data. We discuss the uncertainties affecting (α, n) yield estimations and propose a strategy to enhance their accuracy. Furthermore, this paper discusses and emphasizes the need for new measurements of (α, n) cross-sections for a variety of relevant materials. Such measurements are essential for improving neutron flux predictions, which are crucial for reducing uncertainties in sensitivity estimates for next-generation physics experiments operating in the keV—MeV range.
The knowledge of how neutrinos interact in nuclei is critical to an understanding of finite nuclei, neutrino physics and also to astrophysical environments like neutron stars and supernovae. Within this project we studied how nuclei decay through the emission of a neutrino and a lepton, i.e. β-decay. We implemented realistic treatments of many-nucleon correlations and currents to enable high-precision studies of neutrino physics.
The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.
The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.
The Goos-Hänchen (GH) shift describes a phenomenon in which a specularly reflected beam is translated along the reflecting surface such that the incident and reflected rays no longer intersect at the surface. Using a neutron spin-echo technique and a specially designed magnetic multilayer mirror, we have measured the relative phase between the reflected up and down neutron spin states in total reflection. The relative GH shift calculated from this phase shows a strong resonant enhancement at a particular incident neutron wave vector, which is due to a waveguiding effect in one of the magnetic layers. Calculations based on the observed phase difference between the neutron states indicate a propagation distance along the waveguide layer of 0.65 mm for the spin-down state, which we identify with the magnitude of the giant GH shift. The existence of a physical GH shift is confirmed by the observation of neutron absorption in the waveguide layer. Here, we propose ways in which our experimental method may be exploited for neutron quantum-enhanced sensing of thin magnetic layers.
The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.
In September 1980, the initial delayed critical configuration of the Solution High-Energy Burst Assembly (SHEBA I) was assembled, and the US Department of Energy’s Oak Ridge National Laboratory (ORNL) staff performed near critical and subcritical measurements using the californium source-driven noise analysis (CSDNA) method at the Los Alamos National Laboratory (LANL) Critical Experiments Facility. An unreflected 56 cm outside diameter stainless-steel cylindrical tank was partially filled with uranyl fluoride solution (235U enrichment was 4.95 wt %) until delayed criticality was achieved. Then, measurements were performed for various fuel solution heights from delayed critical to 60% of the height required for delayed criticality. The stainless-steel tank had an inside diameter of 54.6 cm and a height of 105 cm and was partially (20–36.5 cm) filled with an aqueous solution of uranyl fluoride (with a density of 2.162 g/cm 3 ). The density of the uranium in solution was 1.042 g/cm 3 , and the solution had a H/U atomic ratio of 20.43. The tank had a 6.35 mm wall thickness and an axial reentrant tube with an inside diameter of 6.02 cm and wall thickness of 0.165 cm. The tank was essentially unreflected on the top and sides because it was in a thin metal shed. The reactivity of the near–delayed critical configuration was −10.4 cents, which corresponds to a k eff value of 0.99922. In addition to the CSDNA measurements, the prompt neutron decay constants were determined from break frequency noise analysis (BNFA) measurements. The subcritical neutron multiplication factors from CSDNA and BFNA compared extremely well. These data can be used as the basis of International Nuclear Criticality Safety Evaluation benchmark for the near-critical configuration, and the k eff values at various subcritical configurations can be used as subcritical benchmarks. Furthermore, the measured prompt neutron decay constants can be used as reactor physics benchmarks. This report documents the experimental information for the measurements performed so that, at a later date, researchers could perform the required uncertainty and calculational analyses and documentation to use these data for an International Criticality Safety Benchmark Evaluation Program (ICSBEP) or Nuclear Energy Agency benchmarks. The measured prompt neutron decay constants can be used as the basis of a benchmark for the International Reactor Physics Evaluation Program (IRPhEP). The data from these measurements are available from the ORNL Records Management Services Department, and the logbook is available from ICSBEP at Idaho National Laboratory.
Tetrakis(phenylethynyl)benzene (TPEB), a high char yield arylacetylene resin, was blended with a resorcinol-based PEEK™-like oligomeric phthalonitrile resin (Res). Here, we observed the influence of each polymer precursor on the structural changes and material dynamics of different mixtures of these resins during thermal curing to produce thermosets with high char retention. We provide insight to help understand fundamental processes that control and tune this complex thermally driven curing mechanism, which, to date, has not been well understood. We show that an equivalent blend of Res and TPEB mitigates the exothermic curing process and exhibits a low viscosity, which may facilitate processing of these thermosetting blends for composite manufacturing. To obtain a more comprehensive overview of the different dynamics of these polymers during curing, we correlated thermodynamic changes with rheology, materials characterization, and molecular dynamics derived from neutron scattering. Our efforts shed more light on the thermally driven chemical changes, physical transformation, and gelation dynamics that drive or accompany softening and curing of these two distinct crosslinking resin groups.
We extend the excluded-volume model of isospin symmetric two-flavor dense quarkyonic matter [Phys. Rev. C 101, 035201 (2020)] including strange particles and address its implications for neutron stars. The effective sizes of baryons are defined from the diverging hard-core potentials in the short interdistance regime. Around the hard-core density, the repulsive core between baryons at short distances leads to a saturation in the number density of baryons and generates perturbative quarks from the lower phase space, which leads to the shell-like distribution of baryons by the Pauli exclusion principle. The strange-quark Fermi sea always appears at high densities but the Λ hyperon shell only appears when the effective size of the Λ hyperon is smaller than the effective size of nucleons. We find that the pressure of strange quarkyonic matter can be large enough to support neutron stars with two times solar mass and can have a large sound speed, $c^2_s ≃ 0.7$. Finally, the fraction of the baryon number carried by perturbative quarks is about 30% at the inner core of most massive neutron stars.
The principal goal of this work was to better understand the fusion of neutron-rich nuclei a topic relevant to the fields of both nuclear physics and nuclear astrophysics. The work provides insight into the structure and reactions of neutron-rich nuclei namely the extent of their neutron density distribution and its polarizability. By comparing the fusion excitation functions for a chain of isotopes with a common target nucleus changes in the attractive nuclear potential are assessed. This change in the attractive potential is related to changes in the neutron density distribution with increasing number of neutrons or changes in fusion dynamics with increasing neutron number. The impact of the pairing of valence neutrons and protons on the fusion cross-section is also examined. Measurement of an isotopic chain is a powerful tool to address this topic. The experimental program made use of several different accelerator facilities. The core of the experimental work involved experiments at the ReA3 accelerator and the Facility for Rare Isotope Beams (FRIB) situated at Michigan State University, at GANIL, the French national nuclear physics laboratory, and the University of Notre Dame. At ReA3 the degree to which α-clusters associated with fusion of 28,30,32 Si + 28 Si result from the collision dynamics or reflect an initial α-cluster structure was explored. Alpha clusters observed in fusion reactions exceed the predictions of the standard statistical model. This experiment provides a measure of how clusterization is impacted by the increasing neutron-richess of the system. In the experiment at GANIL we investigated fusion in 19 O + 12 C and 20 O + 12 C. This experiment utilized the recently developed active-target detector MuSIC@Indiana. Using this proven, highly efficient, active-target detector we measured the fusion excitation function for these neutron-rich systems. Grounded by these experimental measurements, theoretical models were used to examine the role of unpaired valence neutrons on the fusion cross-section at energies just above the fusion barrier.
This presentation discusses project goals, including: develop new UWBG ternary nitride materials for radiation tolerant neutron detectors, understand stability and defect behavior, and demonstrate device functionality.
Templates of uncertainties expected in specific measurement types were recently developed. One aim of these templates is to help evaluators in identifying (1) missing or suspiciously low uncertainties and (2) missing correlations between uncertainties of the same and different experiments, when estimating covariances for experimental data employed in their evaluations. These templates also provide realistic estimates of standard deviations and correlations for a particular uncertainty source and measurement type that can be used by evaluators in situations where they are not supplied by the experimenters. This information allows for a more comprehensive uncertainty analysis across all measurements considered in an evaluation and, thus, more realistic evaluated covariances. Here, in this work, we extend a template that is applicable to uncertainties expected in neutron-induced fission, (n,f), cross-section measurements. It is applied to improving covariances of 239 Pu(n,f) cross-section measurements in the database underlying the Neutron Data Standards evaluations. This particular example was chosen since this evaluation is primarily based on experimental information. Also, some uncertainties of individual 239 Pu(n,f) cross-section experiments in this database were suspected to be underestimated. The evaluated uncertainties obtained after updating the covariances in the database by means of the template indeed do increase compared to their original values. Even more importantly, the evaluated mean values change noticeably. These modified cross sections impact application calculations significantly, as is demonstrated by employing them in simulations of the effective neutron multiplication factor for a few selected critical assemblies. However, this updated evaluated 239 Pu(n,f) cross section should not be interpreted as the final one that should replace values of the current Neutron Data Standards project. Evaluations for the Neutron Data Standards of the 239 Pu(n,f) cross section must be linked to many other observables included in the associated database, most notably to cross sections for 235 U(n,f), but also to those for 10 B(n, α ), 6 Li(n,t), 238 U(n,f), and 238 U(n, γ ), because of included measurements of the 239 Pu(n,f) cross section that appear as ratios to these reactions. Some of these other reactions are correlated to further observables in the database. Hence, updating uncertainties of data sets of any of these observables can potentially impact the 239 Pu(n,f) cross section. Uncertainties for all measurements of these linked physical observables have to be updated before a comprehensive evaluation of the 239 Pu(n,f) cross section and its corresponding uncertainties can be provided.
Abstract not provided.
The Facility for Antiproton and Ion Research (FAIR) is in its final construction stage next to the campus of the Gesellschaft für Schwerionenforschung Helmholtzzentrum for heavy-ion research in Darmstadt, Germany. Once it starts its operation, it will be the main nuclear physics research facility in many basic sciences and their applications in Europe for the coming decades. Owing to the ability of the new fragment separator, Super-FRagment Separator, to produce high-intensity radioactive ion beams in the energy range up to about 2 GeV/nucleon, these can be used in various nuclear reactions. This opens a unique opportunity for various nuclear structure studies across a range of fields and scales: from low-energy physics via the investigation of multi-neutron systems and halos to high-density nuclear matter and the equation of state, following heavy-ion collisions, fission and study of short-range correlations in nuclei and hypernuclei. Here, the newly developed reactions with relativistic radioactive beams (R3B) set up at FAIR would be the most suitable and versatile for such studies. An overview of highlighted physics cases foreseen at R3B is given, along with possible future opportunities, at FAIR.
This Engineering Calculations and Analysis Report (ECAR) is a summary of the detailed Monte Carlo N-Particle (MCNP) physics analysis performed for the Purdue JW 16297/Fracture of Irradiated Neutron Embrittled Structural Steel and Examination (FINESSE) experiment irradiated in the B-10 position in Advanced Test Reactor (ATR), Cycle 171A.
We report the symmetry energy and its density dependence are pivotal for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283 ± 0.071 fm for the neutron-skin thickness of 208 Pb, $R^{^{208}Pb}_{skin}$, implying a symmetry-energy slope parameter L of 106 ± 37 MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S 0 , L, and $R^{^{208}Pb}_{skin}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with chiral effective field theory (χEFT) and PREX-II constraints yields S 0 = 33.0$^{+2.0}_{-1.8}$ MeV, L = 53$^{+14}_{-15}$ MeV, and $R^{^{208}Pb}_{skin}$ = 0.17 $^{+0.04}_{-0.04}$ fm. We also examine the consistency of several individual χEFT calculations with astrophysical observations and terrestrial experiments. We find that there is only mild tension between χ EFT, astrophysical data, and PREX-II's $R^{^{208}Pb}_{skin}$ measurement (p value = 12.3%) and that there is excellent agreement between χEFT, astrophysical data, and other nuclear experiments.