Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear shell model”

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 415 records · Page 23

Extensions of a classical mechanics “piston-model” for understanding the impact of asymmetry on ICF implosions: The cases of mode 2, mode 2/1 coupling, time-dependent asymmetry, and the relationship to coast-time

As long suspected, low mode asymmetry in inertially confined fusion (ICF) implosions has been implicated as a performance limiting factor [Casey et al., “Evidence of three-dimensional asymmetries seeded by high-density carbon-ablator nonuniformity in experiments at the national ignition facility,” Phys. Rev. Lett. 126, 025002 (2021)]. Recently a non-linear, but solvable, theory [Hurricane et al., “An analytic asymmetric-piston model for the impact of mode-1 shell asymmetry on ICF implosions,” Phys. Plasmas 27, 062704 (2020)] based upon the simple picture of a pair of asymmetric pistons has generated new insights and provided some practical formulas for estimating the degradation of an implosion due to mode-1 asymmetry and demonstrated a previously unrecognized connection between measured hot-spot drift velocity, nuclear down-scatter ratio asymmetry, and the concept of residual kinetic energy (RKE). Asymmetry of the implosion “shell,” as opposed to asymmetry of the hot-spot, was key to the classical mechanics model because the majority of the kinetic energy in an implosion is carried by the shell. Herein, the two-piston model is extended to a six-piston model in order to capture mode-2 asymmetry and coupling between mode-1 and mode-2. A key result of this new six-piston model is that the weighted harmonic mean of shell areal density is the fundamental quantity that determines the RKE and performance degradations for a three-dimensional implosion. Agreement is found between the scalings coming from the theory and ICF implosion data from the National Ignition Facility and to large ensembles of detailed simulations. The connection between the piston model's dependence upon the radius of peak velocity and coast-time is also highlighted in this paper. Finally, by extending the two-piston model to include time-dependent “swing,” it is shown in the Appendix that the shell asymmetry at the time of stagnation dominates the solution for RKE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure of 30 Mg explored via in-beam γ-ray spectroscopy

In the island of inversion, ground states of neutron-rich sd-shell nuclei exhibit strong admixtures of intruder configurations from the fp shell. The nucleus 30 Mg, located at the boundary of the island of inversion, serves as a cornerstone to track the structural evolution as one approaches this region. Spin-parity assignments for excited states in 30 Mg, especially negative-parity levels, have yet to be established. In the present work, the nuclear structure of 30Mg was investigated by in-beam gamma-ray spectroscopy mainly focusing on firm spin-parity determinations. High-intensity rare-isotope beams of 31 Mg, 32 Mg, 34 Si, and 35 P bombarded a Be target to induce nucleon removal reactions populating states in 30 Mg. Gamma rays were detected by the state-of-the-art gamma-ray tracking array GRETINA. For the direct one-neutron removal reaction, final-state exclusive cross sections and parallel momentum distributions were deduced. Multi-nucleon removal reactions from different projectiles were exploited to gain complementary information. With the aid of the parallel momentum distributions, an updated level scheme with revised spin-parity assignments was constructed. Spectroscopic factors associated with each state were also deduced. Results were confronted with large-scale shell-model calculations using two different effective interactions, showing excellent agreement with the present level scheme. Furthermore, a marked difference in the spectroscopic factors indicates that the full delineation of the transition into the island of inversion remains a challenge for theoretical models.

20 ≤ A ≤ 38↗

A high-entropy-wind r-process study based on nuclear-structure quantities from the new finite-range droplet model FRDM(2012)

Attempts to explain the source of r-process elements in our solar system (S.S.) by particular astrophysical sites still face entwined uncertainties, stemming from the extrapolation of nuclear properties far from stability, inconsistent sources of different properties (e.g., nuclear masses and β-decay properties), and the poor understanding of astrophysical conditions, which are hard to disentangle. In this paper we present results from the investigation of r-process in the high-entropy wind (HEW) of core-collapse supernovae (here chosen as one of the possible scenarios for this nucleosynthesis process), using new nuclear-data input calculated in a consistent approach, for masses and β-decay properties from the new finite-range droplet model FRDM(2012). The accuracy of the new mass model is 0.56 MeV with respect to AME2003, to which it was adjusted. We compare the new HEW r-process abundance pattern to the latest S.S. r-process residuals and to our earlier calculations with the nuclear-structure quantities based on FRDM(1992). Substantial overall and specific local improvements in the calculated pattern of the r-process between A ≅ 110 and {sup 209}Bi, as well as remaining deficiencies, are discussed in terms of the underlying spherical and deformed shell structure far from stability.

79 ASTRONOMY AND ASTROPHYSICS↗

Solving The Long-Standing Problem Of Nuclear Reactions At The Highest Microscopic Level: Annual Continuation And Progress Report

The aim of this project is to develop a comprehensive framework that will lead to a fundamental description of both structural properties and reactions of light nuclei in terms of constituent protons and neutrons interacting through nucleon-nucleon (NN) and three-­nucleon (NNN) forces. This project will provide the research community with the theoretical and computational tools that will enable: (1) an accurate prediction for fusion reactions that power stars and Earth-­based fusion facilities; (2) an improved description of the spectroscopy of exotic nuclei, including light Borromean systems; and (3) a fundamental understanding of the three-­nucleon force in nuclear reactions and nuclei at the drip line. To achieve this goal, we build upon a promising technique emerged recently as a candidate to reach a fundamental description of low-­energy binary reactions between light ions, that is the ab initio no-­core shell model combined with the resonating-­group method (NCSM/RGM). This approach has demonstrated the capability to describe binary reactions below the three-­body breakup threshold based, up to now, on similarity-­renormalization-­group (SRG) 5 evolved NN only potentials. To advance the understanding of nuclear reactions at low energies and light exotic nuclei, this project aims at extending the NCSM/RGM approach to include the full range of NNN interactions as well as the treatment of three-­cluster bound and continuum states. Three-­nucleon interactions are unavoidable components of a fundamental nuclear Hamiltonian obtained in a low-­energy effective theory. In addition, three-­nucleon force terms are induced by the SRG procedure and have to be taken into account for such a transformation to be unitary in many-­body calculations. At the same time, the introduction of three-­body cluster states is key to achieve a microscopic description of Borromean systems as well as three-­body breakup reactions. This project will both enhance the fundamentality and enlarge the scope of our microscopic description of nuclear properties. A successful completion of this project will result in improved accuracy of the 3 He( 3 He,2p) 4 He and 3 He(α,γ) 7 Be reaction rates and consequently, in enhancement of the predictive capability of the standard solar model. In addition, we will study also the mirror reactions 3 H( 3 H,2n) 4 He and 3 H(α,γ) 7 Li (a key reaction for the production of 7 Li in the standard big-­bang nucleosynthesis), and the spectroscopy of the 6 He, 6 Be, and 11 Li nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Beta-Decay Half-Lives beyond 54 Ca: A Systematic Survey of Decay Properties Approaching the Neutron Dripline

In an experiment performed at the Facility for Rare Isotope Beams (FRIB) using the FRIB Decay Station initiator, 15 new half-lives of isotopes near 54 Ca were measured. A new method of extracting lifetimes from experimental data, taking into account the unknown 𝛽-delayed neutron emission branches of very neutron-rich nuclei, was developed to enable systematic uncertainty analysis. The experiment observed a dramatic change in the half-life systematics for the isotopes with neutron number 𝑁 = 34. Beyond 𝑁 = 34, the decline of nuclear lifetime is much slower, leading to longer than anticipated lifetimes for near-dripline nuclei. State-of-the-art shell-model calculations can explain the experimental results, revealing the imprint of shell effects and the need for modification of single-particle neutron states. In conclusion, the results from a newly developed quasiparticle random phase approximation model with potential for making global predictions were also tested against the experimental results and good agreement was found.

Physics - Nuclear physics and radiation physics↗

Nuclear ground-state masses and deformations: FRDM(2012)

We tabulate the atomic mass excesses and binding energies, ground-state shell-plus-pairing corrections, ground-state microscopic corrections, and nuclear ground-state deformations of 9318 nuclei ranging from {sup 16}O to A=339. The calculations are based on the finite-range droplet macroscopic and the folded-Yukawa single-particle microscopic nuclear-structure models, which are completely specified. Relative to our FRDM(1992) mass table in Möller et al. (1995), the results are obtained in the same model, but with considerably improved treatment of deformation and fewer of the approximations that were necessary earlier, due to limitations in computer power. The more accurate execution of the model and the more extensive and more accurate experimental mass data base now available allow us to determine one additional macroscopic-model parameter, the density-symmetry coefficient L, which was not varied in the previous calculation, but set to zero. Because we now realize that the FRDM is inaccurate for some highly deformed shapes occurring in fission, because some effects are derived in terms of perturbations around a sphere, we only adjust its macroscopic parameters to ground-state masses. The values of ten constants are determined directly from an optimization to fit ground-state masses of 2149 nuclei ranging from {sup 16}O to {sub 106}{sup 265}Sg and {sub 108}{sup 264}Hs. The error of the mass model is 0.5595 MeV for the entire region of nuclei included in the adjustment, but is only 0.3549 MeV for the region N≥65. We also provide masses in the FRLDM, which in the more accurate treatments now has an error of 0.6618 MeV, with 0.5181 MeV for nuclei with N≥65, both somewhat larger than in the FRDM. But in contrast to the FRDM, it is suitable for studies of fission and has been extensively so applied elsewhere, with FRLDM(2002) constants. The FRLDM(2012) fits 31 fission-barrier heights from {sup 70}Se to {sup 252}Cf with a root-mean-square deviation of 1.052 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Monte Carlo calculations of weak transitions in A = 6 – 10 nuclei

{\it Ab initio} calculations of the Gamow-Teller (GT) matrix elements in the $\beta$ decays of $^6$He and $^{10}$C and electron captures in $^7$Be are carried out using both variational and Green's function Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange phenomenology or chiral effective field theory. The agreement with experimental data is excellent for the electron captures in $^7$Be, while theory overestimates the $^6$He and $^{10}$C data by $\sim 2\%$ and $\sim 10\%$, respectively. We show that for these systems correlations in the nuclear wave functions are crucial to explain the data, while many-body currents increase by $\sim 2$--$3\%$ the one-body GT contributions. These findings suggest that the longstanding $g_A$-problem, {\it i.e.}, the systematic overprediction ($\sim 20 \%$ in $A\le 18$ nuclei) of GT matrix elements in shell-model calculations, may be resolved, at least partially, by correlation effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Possible determination of high-lying single-particle components with ( p , d ) reactions

A detailed feasibility study on deducing the high-lying single particle components (HLSPCs), which are important but used to be ignored, in the ground and low-lying excited states of even-even light nuclei is performed by analyses of (p, d) reactions with 12 C, 24 Mg, 28 Si, and 40 Ca targets at 51.93 MeV. Coupled reaction channel analyses have been made for (p, d) transitions to the j-forbidden excited states in 11 C (5/2 – ; 4.32 MeV), 23 Mg (7/2 + ; 2.05 MeV), 27 Si (7/2 + ; 2.16 MeV), and 39 Ca (9/2 – ; 3.64 MeV), including the major allowed transition components together with direct components of HLSPCs. Spectroscopic amplitudes of the HLSPCs are deduced by fitting the angular distributions of the ground and the j-forbidden excited states simultaneously. Lastly, the present analysis demonstrates for the first time that information about HLSPCs in atomic nuclei can be obtained from analysis of (p, d) reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Uncertainty quantification of collective nuclear observables from the chiral potential parametrization

We perform an uncertainty estimate of quadrupole moments and B(E2) transition rates that inform nuclear collectivity. In particular, we study the low-lying states of 6 Li and 12 C using the ab initio symmetry-adapted no-core–shell model. For a narrow standard deviation of approximately 1% on the low-energy constants which parametrize high-precision chiral potentials, we find output standard deviations in the collective observables ranging from approximately 3%–6%. The results mark the first step towards a rigorous uncertainty quantification of collectivity in nuclei that aims to account for all sources of uncertainty in ab initio descriptions of challenging collective and clustering observables.

ab initio↗

Evolving R Coronae Borealis stars with mesa

The R Coronae Borealis (RCB) stars are rare hydrogen-deficient, carbon-rich supergiants. They undergo extreme, irregular declines in brightness of many magnitudes due to the formation of thick clouds of carbon dust. It is thought that RCB stars result from the mergers of CO/He white dwarf (WD) binaries. We constructed post-merger spherically symmetric models computed with the mesa code, and then followed the evolution into the region of the Hertzsprung-Russell (H-R) diagram where the RCB stars are located. We also investigated nucleosynthesis in the dynamically accreting material of CO/He WD mergers which may provide a suitable environment for significant production of 18 O and the very low 16 O/ 18 O values observed. We have also discovered that the N abundance depends sensitively on the peak temperature in the He-burning shell. Our mesa modelling consists of engineering the star by adding He-WD material to an initial CO-WD model, and then following the post-merger evolution using a nuclear-reaction network to match the observed RCB abundances as it expands and cools to become an RCB star. These new models are more physical because they include rotation, mixing, mass-loss, and nucleosynthesis within mesa . We follow the later evolution beyond the RCB phase to determine the stars’ likely lifetimes. The relative numbers of known RCB and extreme helium stars correspond well to the lifetimes predicted from the mesa models. In addition, most of computed abundances agree very well with the observed range of abundances for the RCB class.

Astronomy & Astrophysics↗

Total kinetic energy release in the fast neutron induced fission of actinide nuclei

The total kinetic energy release and fission mass distributions for the fast neutron (En = 3–100 MeV) induced fission of 232 Th, 233 U, 235 U, 237 Np, 239 Pu, 240 Pu, and 242 Pu have been measured using the LANSCE facility. The neutron energies were deduced from time-of- flight measurements. The fission fragments were detected using Si PIN diode detectors, giving us the fragment energies. The actinide targets were made by vapor deposition leading to high-quality targets, that were thin and uniform with reduced impurities. Corrections were made to the data for pulse height defect and the fragment energy loss in the target and its backing. The TKE distributions were Gaussian in shape and their mean value as a function of incoming neutron energy could be fitted with second order polynomials. In the case of 233 U and 235 U, our measurements agree with prior work. Our measurements for 232 Th are unique. Our data agree with Viola scaling. The constant position of the heavy mass peak is interpreted as being due to the influence of the N = 88 and Z = 50 shells. The GEF model predictions agree with the data in general as do the CGMF model predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

THE STRUCTURE FUNCTION OF THE FREE NEUTRON AT HIGH X-BJORKEN

Understanding the internal structure of nucleons is one of the primary goal of nuclear physicists. As protons and neutrons are only the bound state solution of the QCD lagrangian (at least inside atomic nuclei), studying protons and neutrons helps uncover nuclear struc ture. Due to its easy availability, many studies on protons have been done on a wide range of kinematics. However, free neutron targets are not readily achievable. So, any information on neutrons has to be extracted from neutron-rich nuclei, and some nuclear models have to be used to subtract the contributions from other nucleons to extract the information on neutrons. So, the Barely Off-shell Nucleon Structure (BONuS12) experiment at Jefferson Lab was conducted to overcome these challenges by using spectator tagging. The experiment effectively created a quasi-free neutron target by scattering electrons off a deuterium target and detecting low-momentum, backward-moving protons using a custom-built Radial Time Projection Chamber (RTPC). Selecting the low momentum and backward-moving spectators would enable us to minimize the model-dependent effects due to final state interactions and target fragmentation. The RTPC was a 40 cm-long cylindrical detector that works on the principle of gaseous ionization. It had three layers of Gas Electron Multipliers (GEMs) for charge amplification and a surrounding readout pad. The scattered electrons were measured using the CLAS12 detector, and data were collected using a 10.4 GeV electron beam dur ing Spring and Summer 2020. Using spectator tagging, we extracted the structure function ratio Fn 2 of the quasi-free neutron in the deep inelastic scattering at high x, upto x ~ 0.8. The result was extracted in the region with the invariant mass W > 1.8 GeV/c2, and Q2 in the range 1.3 to 11 GeV2. This dissertation presents the methodology, event selection criteria and refinements, estimation and subtraction of backgrounds, and complete analysis of extraction of Fn 2/Fp 2 in a model-independent way. Also, systematic uncertainties in our final analysis will be discussed in detail.

Pokhrel, Madhusudhan [Old Dominion Univ., Norfolk,↗

New 26 P( p, γ ) 27 S Thermonuclear Reaction Rate and Its Astrophysical Implications in the rp -process

Accurate nuclear reaction rates for 26 P(p, γ) 27 S are pivotal for a comprehensive understanding of the rp-process nucleosynthesis path in the region of proton-rich sulfur and phosphorus isotopes. However, large uncertainties still exist in the current rate of 26 P(p, γ) 27 S because of the lack of nuclear mass and energy level structure information for 27 S. We reevaluate this reaction rate using the experimentally constrained 27 S mass, together with the shell model predicted level structure. It is found that the 26 P(p, γ) 27 S reaction rate is dominated by a direct capture reaction mechanism despite the presence of three resonances at E = 1.104, 1.597, and 1.777 MeV above the proton threshold in 27 S. The new rate is overall smaller than the other previous rates from the Hauser–Feshbach statistical model by at least 1 order of magnitude in the temperature range of X-ray burst interest. In addition, we consistently update the photodisintegration rate using the new 27 S mass. The influence of new rates of forward and reverse reaction in the abundances of isotopes produced in the rp-process is explored by postprocessing nucleosynthesis calculations. The final abundance ratio of 27 S/ 26 P obtained using the new rates is only 10% of that from the old rate. The abundance flow calculations show that the reaction path 26 P(p, γ) 27 S(β + ,ν) 27 P is not as important as previously thought for producing 27 P. The adoption of the new reaction rates for 26 P(p, γ) 27 S only reduces the final production of aluminum by 7.1% and has no discernible impact on the yield of other elements.

79 ASTRONOMY AND ASTROPHYSICS↗

Determination of the 60 Zn level density from neutron evaporation spectra

Nuclear reactions of interest for astrophysics and applications often rely on statistical model calculations for nuclear reaction rates, particularly for nuclei far from β stability. However, statistical model parameters are often poorly constrained, where experimental constraints are particularly sparse for exotic nuclides. For example, our understanding of the breakout from the NiCu cycle in the astrophysical rp-process is currently limited by uncertainties in the statistical properties of the proton-rich nucleus 60 Zn. We have determined the nuclear level density of 60 Zn using neutron evaporation spectra from 58 Ni( 3 He,n) measured at the Edwards Accelerator Laboratory. We compare our results to a number of theoretical predictions, including phenomenological, microscopic, and shell-model-based approaches. Notably, we find the 60 Zn level density is somewhat lower than expected for excitation energies populated in the 59 Cu(p,γ) 60 Zn reaction under rp-process conditions. This includes a level density plateau from roughly 5 to 6 MeV excitation energy, which is counter to the usual expectation of exponential growth and all theoretical predictions that we explore. Here, a determination of the spin distribution at the relevant excitation energies in 60 Zn is needed to confirm that the Hauser-Feshbach formalism is appropriate for the 59 Cu(p,γ) 60 Zn reaction rate at x-ray burst temperatures

59 ≤ A ≤ 89↗

Form factors for dark matter capture by the Sun in effective theories

In the effective theory of isoscalar and isovector dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle, 8 isotope-dependent nuclear response functions can be generated in the dark matter scattering by nuclei. We compute the 8 nuclear response functions for the 16 most abundant elements in the Sun, i.e. H, {sup 3}He, {sup 4}He, {sup 12}C, {sup 14}N, {sup 16}O, {sup 20}Ne, {sup 23}Na, {sup 24}Mg, {sup 27}Al, {sup 28}Si, {sup 32}S, {sup 40}Ar, {sup 40}Ca, {sup 56}Fe, and {sup 59}Ni, through numerical shell model calculations. We use our response functions to compute the rate of dark matter capture by the Sun for all isoscalar and isovector dark matter-nucleon effective interactions, including several operators previously considered for dark matter direct detection only. We study in detail the dependence of the capture rate on specific dark matter-nucleon interaction operators, and on the different elements in the Sun. We find that a so far neglected momentum dependent dark matter coupling to the nuclear vector charge gives a larger contribution to the capture rate than the constant spin-dependent interaction commonly included in dark matter searches at neutrino telescopes. Our investigation lays the foundations for model independent analyses of dark matter induced neutrino signals from the Sun. The nuclear response functions obtained in this study are listed in analytic form in an appendix, ready to be used in other projects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

β -delayed neutron emissions from N > 50 gallium isotopes

β-delayed γ-neutron spectroscopy has been performed on the decay of A = 84 to 87 gallium isotopes at the RI-beam Factory at the RIKEN Nishina Center using a high-efficiency array of 3 He neutron counters (BRIKEN). β-2n-γ events were measured in the decays of all of the four isotopes for the first time, which is direct evidence for populating the excited states of two-neutron daughter nuclei. Detailed decay schemes with the γ branching ratios were obtained for these isotopes, and the neutron emission probabilities (P xn ) were updated from the previous study. Hauser-Feshbach statistical model calculations were performed to understand the experimental branching ratios. We found that the P 1n and P 2n values are sensitive to the nuclear level densities of 1n daughter nuclei and showed that the statistical model reproduced the P 2n /P 1n ratio better when experimental levels plus shell-model level densities fit by the Gilbert-Cameron formula were used as the level-density input. We also showed the neutron and γ branching ratios are sensitive to the ground-state spin of the parent nucleus. Our statistical model analysis suggested J ≤ 3 for the unknown ground-state spin of the odd-odd nucleus 86 Ga, from the I γ (4 + → 2 + )/I γ (2 + → 0 + ) ratio of 84 Ga and the P 2n /P 1n ratio. In conclusion, these results show the necessity of detailed understanding of the decay scheme, including data from neutron spectroscopy, in addition to γ measurements of the multineutron emitters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗