SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Toward shell model interactions with credible uncertainties

Background: The nuclear shell model is a powerful framework for predicting nuclear structure observables, but relies on interaction matrix elements fit to experimental data as its inputs. Extending the shell model's applicability, particularly toward dripline nuclei, requires efficient fitting methods and credible uncertainty quantification. Traditional approaches face computational challenges and may underestimate uncertainties. Purpose: We develop and test a framework combining eigenvector continuation and Markov chain Monte Carlo to efficiently fit shell model interaction matrix elements and quantify their uncertainties. Methods: Eigenvector continuation is used to emulate shell model calculations, reducing computational costs. The emulator enables Markov chain Monte Carlo sampling to optimize interaction matrix elements and rigorously assess parametric uncertainties. Here, the framework is benchmarked using the USDB interaction in the 𝑠⁢𝑑 shell. Results: The emulator reproduces the USDB interaction with negligible error, validating its use in shell model fitting applications. However, we find that to obtain credible predictive intervals, the model defect of the shell model itself, rather than experimental or emulator error, must be taken into account in order to obtain credible uncertainties. Conclusions: The proposed framework provides an efficient and rigorous approach for fitting shell model interactions and quantifying uncertainties. Further, the normality assumption used in the past appears sufficient to describe the distribution of interaction matrix elements. However, it is crucial to account for model correlations to avoid underestimating uncertainties.

Nuclear forces↗

DOE OSTI2025

Response functions and giant monopole resonances for light to medium-mass nuclei from the ab initio symmetry-adapted no-core–shell model

Using the ab initio symmetry-adapted no-core–shell model, we compute sum rules and response functions for light to medium-mass nuclei, starting from interactions that are derived in the chiral effective field theory. Specifically, we investigate electromagnetic transitions of monopole, dipole and quadrupole nature for 4 He, and explore dominant features of giant monopole resonances in symmetric nuclei such as the closed-shell 4 He and 16 O light nuclei, the intermediate-mass open-shell 20 Ne and the medium-mass closed-shell 40 Ca. Furthermore, for the NNLO opt chiral potential, we determine parameter-free monopole sum rules, which can provide information on the incompressibility of symmetric nuclear matter. Here, we report 213(10) MeV as an estimate for the compression modulus for infinite nuclear matter, which overlaps with the lower range of values often used in current astrophysical applications.

ab initio nuclear structure↗

DOE OSTI2025

Advancing quantum simulations of the nuclear shell model with Gray-code–based resource-efficient protocols

Background: Some of the computational limitations in solving the nuclear many-body problem could be overcome by utilizing quantum computers. The nuclear shell-model calculations providing deeper insights into the properties of atomic nuclei are one such case with high demand for resources, as the size of the Hilbert space grows exponentially with the number of particles involved. Quantum algorithms are being developed to overcome these challenges and advance such calculations. Purpose: To develop quantum circuits for the nuclear shell-model, leveraging the capabilities of noisy intermediate-scale quantum (NISQ) devices. Here, we aim to minimize resource requirements (specifically in terms of qubits and gates) and strive to reduce the impact of noise by employing relevant mitigation techniques. Methods: We achieve noise resilience by designing an optimized Ansatz for the variational quantum eigensolver (VQE) based on Givens rotations and incorporating qubit-ADAPT-VQE in combination with variational quantum deflation (VQD) to compute ground and excited states, incorporating the zero-noise extrapolation mitigation technique. Furthermore, the qubit requirements are significantly reduced by mapping the basis states to qubits using Gray-code encoding and generalizing transformations of fermionic operators to efficiently represent many-body states. Results: By employing the resource-efficient protocols, we achieve the ground and excited state energy levels of 38 Ar and 6 Li with better accuracy. These energy levels are presented for noiseless simulations, noisy conditions, and after applying noise mitigation techniques. Results are compared for Jordan-Wigner and Gray-code encoding using VQE, qubit-ADAPT-VQE, and VQD. Conclusions: Our work highlights the potential of resource-efficient protocols to leverage the full potential of NISQ devices in scaling the nuclear shell model calculations, offering a pathway toward more complex quantum simulations in nuclear physics. This approach establishes a framework for studying other nuclear systems with improved quantum resource efficiency, marking a significant advancement in applying quantum computing to realistic nuclear physics applications.

Physics - Nuclear physics and radiation physics↗

DOE OSTI2025

Exploring isospin symmetry breaking in exotic nuclei: High-precision mass measurement of 23 Si and shell-model calculations of 𝑇 = 5/2 nuclei

Here, we present a high-precision mass measurement of the proton-rich nucleus 23 Si, performed with the LEBIT Penning trap at the Facility for Rare Isotope Beams (FRIB) utilizing the time-of-flight ion cyclotron resonance (TOF-ICR) technique. We determined a mass excess of 23362.9(5.8) keV, which agrees with a recent storage-ring measurement from the experimental Cooler-Storage Ring (CSRe) in Lanzhou but has a factor of 20 improved precision 23 Si is hence the nucleus with the most precisely known mass among all nuclei with an isospin projection of 𝑇 𝑧 = −5/2. We performed shell-model calculations with the USDC and USDCm Hamiltonians to study binding energy differences and Thomas-Ehrmann shifts in mirror systems with an isospin up to 𝑇 = 5/2. Our experimental result and other recently reported masses of neutron-deficient sd-shell nuclei agree well with the theoretical predictions, demonstrating that isospin symmetry breaking in sd-shell nuclei—even at high isospin values—is well described by modern shell-model calculations.

20 ≤ A ≤ 38↗

DOE OSTI2025

A critical verification of beam and shell models of wind turbine blades

Ever-increasing wind turbine size has challenged predictive capabilities on several fronts. Here, to address part of the blade structural modeling uncertainty, a systematic model fidelity comparison study was conducted on commonly used finite elements. pyNuMAD was utilized to create beam, shell, and solid models of a 100 m long blade undergoing large static deflections. The solid model avoided the use of layered-solid elements by resolving core and facesheet layers. An unprecedented model with 73.7 million elements revealed insights that have never been possible from prior experimental and numerical studies. As compared to the solid element model, the tip deflection from the shell and beam model was found to be about 2% and 4.3% too low, respectively. The twist from the beam model was found to be about 5.6% too high, while the twist from shell model was 24% too low, though improvement was demonstrated with mesh refinement. The beam model adhesive stresses were more accurate than the shell model. Out-of-plane stresses were of great significance near geometric and material discontinuities, and neither the shell nor beam model captured these effects well. Failure predictions from beam, shell, or layered-solid models are unlikely to be reliable at trailing edges, adhesives, ply-drops, spar-cap boundaries.

17 WIND ENERGY↗

DOE OSTI2025

Exploring the onset of collectivity approaching 𝑁=40 with masses of manganese isotopes

Isotopes in the region of the nuclear chart below 68 Ni have been the subject of intense experimental and theoretical effort due to the potential onset of a new “island of inversion” when crossing the harmonic oscillator subshell closure at 𝑁=40. Here, we have measured the masses of 64−68 Mn using TITAN's multiple-reflection time-of-flight mass spectrometer, resulting in the first precision mass measurements of 67 Mn and 68 Mn. These results are compared to ab initio calculations and modern shell model calculations and show an increase in collectivity approaching 𝑁=40.

binding energy & masses↗

DOE OSTI2026

Collapse of Nuclear Collectivity along the 𝑁 = 𝑍 Line

The lifetime of the 𝐽 𝜋 = 2 + state in the self-conjugate 88 Ru nucleus has been determined in an experiment performed using rare-isotope beams provided by the new Facility for Rare Isotope Beams. This is the heaviest 𝑁 = 𝑍 nucleus for which such a measurement has been achieved. 88 Ru was populated by both one-neutron knockout and charge-exchange reactions, and the lifetime of 14.3$^{+2.5}_{−3.4}$ ps was determined using the triple-foil plunger technique. The extracted electromagnetic transition strength shows that the quadrupole collectivity has dropped significantly compared with the highly deformed 𝑁 = 𝑍 region around 𝐴 ∼ 80. These results are compared with state-of-the-art large-scale shell-model and discrete nonorthogonal shell-model calculations. The theoretical calculations indicate a moderate triaxial deformation and suggest that low-lying states in this nucleus are no longer dominated by strong many-particle many-hole excitations, unlike the lighter, highly deformed 𝑁 = 𝑍 nuclei nearby.

59 ≤ A ≤ 89↗

DOE OSTI2026

How threshold effects in spectroscopic factors influence heavy-ion knockout reactions

A two-decade-old puzzle in heavy-ion one-nucleon knockout reactions is the strong correlation between the reduction factor R S = σ exp /σ t h and the Fermi surface asymmetry ΔS. Theoretical cross sections typically rely on spectroscopic factors (SFs) from shell model (SM) calculations, which neglect continuum coupling effects. Here, we employ the Gamow shell model (GSM), which explicitly incorporates continuum coupling, to compute SFs for p-shell nuclei and predict corresponding theoretical cross sections. Systematic calculations demonstrate that using GSM-derived SFs substantially reduces discrepancies between theoretical and experimental results. This improvement is particularly significant for deeply bound nucleon knockout in nuclei near the dripline, where traditional SM-based calculations fall short. As a result, using GSM SFs, the ratio R s exhibits no pronounced dependence on ΔS. Furthermore, both the ratio of GSM SFs to SM SFs and their corresponding reaction cross sections ratios exhibit a strong ΔS dependence. We have also compared GSM SFs and cross sections with those from the no-core shell model calculations, giving a similar pronounced sensitivity to ΔS. Detailed analysis attributes these correlations to threshold effects for SFs in weakly bound systems. Overall, incorporating continuum coupling via GSM enhances the reliability of SF predictions for exotic, weakly bound nuclei and provides key insights toward resolving the enduring puzzle in heavy-ion knockout reactions from a nuclear structure perspective.

Gamow shell model↗

DOE OSTI2025

Triaxiality and shape dynamics in 70 Ge

The electromagnetic properties of low-lying states in 70 Ge were investigated via multi-step Coulomb excitation of a 70 Ge beam impinging on a 208 Pb target at the ATLAS facility of the Argonne National Laboratory. A total of 27 transitional elements and six diagonal matrix elements coupling 11 low-lying states, were extracted from the measured cross sections. These were used to calculate reduced transition probabilities, spectroscopic quadrupole moments, and rotational invariant shape parameters, providing enhanced precision and expanding on previous studies. The experimental data were compared within several theoretical frameworks, including the generalized triaxial rotor model, configuration interaction shell-model calculations, and computations within the combined frameworks of relativistic density functional theory and the five-dimensional collective Hamiltonian. The results demonstrate a good agreement with the experimental data and, in conjunction with calculations using a two-state mixing model, support significant triaxiality and strong mixing between the 0$^{+}_{1}$ and 0$^{+}_{2}$ states. This results in the magnitudes of their respective quadrupole deformations [β rms (0$^{+}_{1}$) = 0.228 (3), β rms (0$^{+}_{2}$) = 0.273 (1)] being more similar than previously observed. Furthermore, the implications of these results for understanding the complex shape coexistence phenomena, the role of triaxiality, and shape evolution along the Ge isotopic chain are discussed.

59 ≤ A ≤ 89↗

DOE OSTI2025

Measurements of Gamow-Teller transitions from 59 Co via the 59 Co ⁢(𝑡, 3 He +𝛾) charge-exchange reaction and its application to the stellar electron-capture rates

Electron-capture reactions on iron-group nuclei play a crucial role in the late stages of massive star evolution. Since stellar evolution simulations depend on accurate electron-capture rates—which are highly sensitive to the detailed Gamow-Teller (GT) strength distributions—reliable theoretical models are essential. However, experimental data on GT strength distributions are scarce. High-resolution measurements are therefore vital for benchmarking and improving these theoretical calculations. To provide high-resolution data on Gamow-Teller strength distributions of iron-group nuclei and to compare these results with theoretical calculations within this mass region. Differential cross sections for the 59 Co ⁢(𝑡, 3 He)⁢ 59 Fe charge-exchange reaction at 115 MeV/u were measured using the S800 spectrometer. Furthermore, to resolve individual levels that are not distinguishable in the S800 particle singles data, coincident 𝛾 rays from the 59 Fe residual nucleus were detected by using the Gamma-Ray Energy Tracking In-beam Nuclear Array 𝛾-ray tracking array. Here, the Gamow-Teller transition strength distribution from the ground state of 59 Co to 59 Fe was extracted up to an excitation energy of 10 MeV. Additionally, transition strengths for several low-lying states were determined from coincident 𝛾-ray measurements. Electron-capture rates calculated using the present data indicate that these low-lying states contribute significantly to the overall rates in relevant stellar environments. The experimental results show reasonable agreement with theoretical predictions based on both shell-model and projected shell-model calculations. High-resolution data on Gamow-Teller strength distributions—particularly for individual low-lying states—are essential for accurately determining electron-capture rates in iron-group nuclei. Coincident 𝛾-ray measurements provide a powerful tool for obtaining such detailed information. While the present work demonstrates that shell-model calculations successfully reproduce the experimental results, such comparisons are scarce and more experimental data are desirable.

59 ≤ A ≤ 89↗

DOE OSTI2025

Electromagnetic moments of 215,217 Bi: Probing shell evolution beyond N = 126

The nuclear properties of bismuth isotopes (Z = 83) , with just one valence proton above the closed spherical shell at (Z = 82) , are expected to be governed by a single unpaired proton. However, already in semimagic 209 Bi (Z = 83, N = 126) , , the magnetic moment (μ) strongly deviates from the single-particle Schmidt value. A near linear decrease in μ with the increase of N after the N = 126 magic number was observed up to N = 130 . In order to test whether this trend is kept at N > 130 and to reveal the underlying mechanisms, an investigation of 215,217 Bi (N = 132, 134) has been undertaken. The magnetic dipole and electric quadrupole moments of the I π = 9/2 - nuclear ground states in these isotopes have been measured for the first time using the in-source resonance-ionization spectroscopy technique at ISOLDE (CERN). It has been shown that the linearly decreasing trend of μ( 209,211,213 Bi g ) is broken in 215,217 Bi with a nearly constant value of μ observed. Experimental data have been compared to calculations in the framework of the configuration-interaction shell model with the monopole-based universal V MU +LS interaction. The peculiarities in the behavior of μ(Bi, 9/2 - ) with increasing neutron number are explained as being due to the shell evolution, change of the neutron orbitals occupancies and strong configuration mixing beyond N = 130 . Also, the difference in the μ trends for bismuth (Z = 83) and astatine (Z = 85) isotopes with N > 126 are reproduced by the shell-model calculations. It is shown that monopole interaction plays noticeable role in the description of the peculiarities of the μ behaviour. Additionally, the extension of the application of the V MU interaction to the μ isotopic trends for heavy nuclei is important for further study of the capabilities of this promising version of the shell-model calculations.

Dipole magnetic moments↗

DOE OSTI2026

Isospin symmetry holds beyond the proton drip line: First observation of excited states in 55Cu

Isospin symmetry, combined with shell-model description, has long been successful in reproducing nuclear structure and predicting mirror energy differences. However, these descriptions have been primarily limited to mirror nuclei with well-bound ground states. As one approaches the proton drip line, coupling to the continuum becomes important. Despite this, little is known about medium-mass nuclei at the drip line. Here we report the measurement of the excitation energies of several unbound states in a nucleus beyond the proton drip line, 55Cu, performed via γ-ray spectroscopy at the Radioactive Isotope Beam Factory (RIKEN). Several low-lying yrast states are observed to decay predominantly by γ-ray emission. In contrast, the state–predicted to be strongly populated in the one-proton knockout reaction–is notably absent, suggesting that it decays primarily via prompt proton emission. We show that even for such an unbound system, shell-model calculations based on a harmonic oscillator potential reproduce the observed mirror energy differences remarkably well in all the observed states. Based on this result, we propose that a general pattern might emerge where, in medium-mass systems at the proton drip line, the stronger Coulomb barrier together with the centrifugal barrier effectively confines the wave functions, so that the shell model is able to describe low-lying states, enabling accurate predictions of mirror energy differences between isobaric analogue states.

Pigliapoco, S [Dipartimento di Fisica e Astronomia↗

DOE OSTI2026

Impact of the continuum on the γ decay of the lowest 2 + states in 14 C

The half-life, T 1/2 =14.6(33) fs, of the 7012-keV 2$^{+}_{1}$ state in 14 C was measured in an experiment employing the 9 Be( 6 Li,pγ) fusion-evaporation reaction and the GODDESS setup. The γ decay from the 2$^{+}_{2}$ near-threshold resonance, located 142 keV above the 8176-keV neutron-emission threshold, was also investigated. An upper limit of 4.0 × 10 –5 was established for the γ-decay branching ratio with respect to the neutron-decay channel. The B(E2) transition probabilities for these 2$^{+}_{1}$ and 2$^{+}_{2}$ states were compared to predictions from the Shell Model Embedded in the Continuum (SMEC). Significant modifications of these B(E2) probabilities, with respect to the standard shell model, are expected due to the coupling to the continuum. For calculations using the YSOX interaction, agreement was found for large negative values of V 0 , the coupling constant to the continuum. The central value V 0 = −645 MeV fm 3 results in a retardation by a factor ∼ 3.5 for the 2$^{+}_{1}$ → 0$^{+}_{1}$ transition, and an enhanced probability by a factor ∼ 2.5 for the γ-decay out of the 2$^{+}_{2}$ state. The latter factor reflects the effect of collectivization of the 2$^{+}_{2}$ excitation when the coupling to the continuum is taken into account.

14C↗

DOE OSTI2025

Total absorption spectroscopy for the 𝛽 + decay strength distribution of 60 Ga

𝛽-decay properties play an important role in most astrophysical processes. In the absence of experimental data, astrophysical models rely on global theoretical calculations to provide the relevant properties. It is therefore important to provide strong experimental constraints when possible. In the case of 𝛽-decay, the most sensitive probe is the 𝛽-decay strength distribution. We report here on the first measurement of the latter quantity for the 𝛽 + decay of 60 Ga using the total absorption spectroscopy technique. The experimental results are compared to theoretical calculations often used in astrophysical models, namely the shell model and the quasiparticle random phase approximation (QRPA), as well as an extension of QRPA that includes higher-order nucleonic calculations. Both models are in reasonable agreement with the experimental results.

59 ≤ A ≤ 89↗

DOE OSTI2026

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Quenching of the octupole rotational band in 71 Ge

The main goal of this study was to enhance the knowledge of the properties of the known octupole rotational band in 71 Ge, including transition rates, and to investigate whether the band persists to higher spin. An experiment at Florida State University using the 62 Ni ( 14 C, 𝛼⁢𝑛) reaction at 50 MeV was used to produce 71 Ge at high spin. An array of ten Compton-suppressed Ge detectors, consisting of three Clover detectors and seven single-crystal detectors, was used to measure the 𝛾 decays in coincidence. An analysis of the resulting coincidence spectra resulted in the placement of 13 additional transitions in the 71 Ge level scheme, one of which (1092 keV) extends the octupole band to a (35/2 − ) state at 8206 keV. However, calculations of the kinematic moment of inertia show that the 1092-keV transition disrupts the rigid rotational pattern and is likely associated with a band crossing, potentially quenching the octupole deformation. Transition strengths inferred from lifetime measurements in the octupole band are not well reproduced by either shell-model calculations using the JUN45 interaction or a semimicroscopic cluster model, while those for the 𝜈⁢𝑔 9/2 band are in good agreement with the corresponding shell-model predictions. Comparisons between the octupole band in 71 Ge and negative-parity bands based on 3 − octupole states in some neighboring even-even nuclei show both similarities and differences. Systematic trends within the 𝜈⁢𝑔 9/2 bands among odd-𝐴 Ge isotopes point to increased collectivity near 𝑁 = 40.

59 ≤ A ≤ 89↗

DOE OSTI2025

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗

DOE OSTI2026

Radiative strength functions from the energy-localized Brink-Axel hypothesis

Radiative strength functions (RSFs) model the bulk electromagnetic response of highly excited nuclei and are critical inputs for statistical reaction codes. In this paper, we present a definition of the RSF that is consistent with Hauser-Feshbach reaction codes and that can be efficiently computed with the shell model using the Lanczos strength-function (LSF) method. Here, we introduce a variant of the shell-model LSF method that exploits the energy-localized Brink-Axel hypothesis, which makes it possible to compute both electric and magnetic RSFs across all energies relevant to capture reactions. We verify agreement with the conventional definition of RSFs with benchmark calculations of 24 Mg and then present novel results for 56 Fe. For 56 Fe we find that (i) the 𝑀⁢1 RSF shape evolves smoothly with excitation energy, consistent with the energy-localized Brink-Axel hypothesis; (ii) both 𝑀⁢1 and 𝐸⁢1 transitions contribute significantly to the radiative strength below the photoabsorption threshold; and (iii) within the sdpf model space, the strength below 3 MeV observed in Oslo-type experiments cannot be fully reproduced. These results pave the way for a coherent microscopic description of RSFs and further motivate the use of energy-dependent RSFs in modern reaction codes.

Physics - Nuclear physics and radiation physics↗