Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHELL MODELS”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Structure of high-lying excited states in 47 Ca with ℓ ≥ 3 populated in fast-beam 𝛾-ray-tagged one-neutron pickup reactions

Located just one neutron removed from doubly magic, stable 48 Ca, the adjacent isotope 47 Ca has been studied extensively with neutron-removing transfer reactions and interpreted within a shell-model picture. Neutron-adding transfer reactions onto the stable 46 Ca, however, have not been performed in nearly 60 years and never with the aid of high-resolution 𝛾-ray spectroscopy. Here, we report on the 12 C⁡( 46 , Ca 47 Ca+𝛾)⁢𝑋 one-neutron pickup reaction from 12 C onto 46 Ca at more than 50 MeV/nucleon, a regime of orbital angular momentum mismatch that selectively populates 47 Ca final states based on ℓ ≥ 3 neutron configurations. Newly observed excited states and 𝛾-ray transitions are discussed in comparison to 𝑓⁡𝑝 shell-model calculations in a variety of configuration spaces. Furthermore, it is shown that the nature of the observed spectrum of final states can be well interpreted, albeit with a significant fragmentation of the 𝜈⁢0⁢𝑓 5/2 strength observed in the measurement that is not borne out in the structure calculations.

39 ≤ A ≤ 58↗
DOE OSTI2026

Halo Nuclei from Ab Initio Nuclear Theory

A realistic description of halo nuclei, characterized by low-lying breakup thresholds, requires a proper treatment of continuum effects. We have developed an ab initio approach, the No-Core Shell Model with Continuum (NCSMC), capable of describing both bound and unbound states in light nuclei in a unified way. With chiral two- and three-nucleon interactions as the only input, we can predict the structure and dynamics of halo and other light nuclei and, by comparing to available experimental data, test the quality of chiral nuclear forces. We review NCSMC calculations of weakly bound states and resonances of the exotic halo nuclei 6He, 8B, 11Be, and 15C. For the latter, we discuss its production in the capture reaction 14C(n,𝛾 )15C. We highlight the challenges of a description of 6He as a Borromean n-n-4He system. Finally, we present our calculations of excited states in 10Be exhibiting a one-neutron halo structure and a large scale No-Core Shell Model investigation of 11Li as a precursor of a full n-n-9Li NCSMC study.

Navrátil, Petr↗
DOE OSTI2021

Wigner SU(4) symmetry, clustering, and the spectrum of $$^{12}$$C

Abstract We present lattice calculations of the low-lying spectrum of $$^{12}$$ 12 C using a simple nucleon–nucleon interaction that is independent of spin and isospin and therefore invariant under Wigner’s SU(4) symmetry. We find strong signals for all excited states up to $$\sim 15$$ ∼ 15 MeV above the ground state, and explore the structure of each state using a large variety of $$\alpha $$ α cluster and harmonic oscillator trial states, projected onto given irreducible representations of the cubic group. We are able to verify earlier findings for the $$\alpha $$ α clustering in the Hoyle state and the second $$2^+$$ 2 + state of $$^{12}$$ 12 C. The success of these calculations to describe the full low-lying energy spectrum using spin-independent interactions suggest that either the spin-orbit interactions are somewhat weak in the $$^{12}$$ 12 C system, or the effects of $$\alpha $$ α clustering are diminishing their influence. This is in agreement with previous findings from ab initio shell model calculations.

Shen, Shihang (ORCID:0000000280516466)↗
DOE OSTI2020

Toroidal States of the 12 C Nucleus

Among the states of 12 C, there is an important subset of $K$ = 0 and $K$ = I planar intrinsic states in which the intrinsic motion of the nucleons are confined in the planar region defined by the three-alpha cluster or by their generated toroid. The intrinsic nuclear densities of these states are toroidal in nature. In this work we study these 12 C toroidal states from the generator-coordinate viewpoints in both the alpha cluster model and the toroidal shell model. Numerical solutions in the toroidal mean field approximation are examined to pave the way for future extensions and refinements.

12C Hoyle state↗
DOE OSTI2025

In-beam 𝛾-ray spectroscopy toward the proton dripline: The curious case of 32 Ar

High-resolution in-beam 𝛾-ray spectroscopy was used to study excited states of the neutron-deficient nucleus 32 Ar populated in fast-beam induced four- and six-nucleon removal reactions from 36,38 Ca. One new 𝛾-ray transition and indications for an additional two were found, allowing for a glimpse at the level scheme beyond the 2$^{+}_{1}$ state. The nature of the new 1900⁢(4)-keV transition is discussed in the context of the known energy spectrum of the mirror nucleus 32 Si and shell-model calculations using the FSU and SDPF-M cross-shell effective interactions. Its resulting parent state at 3767⁢(5)⁢keV, more than 1.3 MeV above the proton separation energy, is tentatively assigned to have mixed 𝑠⁢𝑑-shell and 2p-2h character. It might either be the mirror of the 𝐽 𝜋 = 2$^{+}_{2}$ state of 32 Si at 4230.8⁢(8)⁢keV, but with a decay branch favoring a transition to the 2$^{+}_{2}$ over the ground state, or the mirror of the 4983.9⁢(11)-keV state with quantum numbers 0 + . Furthermore, the resulting mirror-energy differences of −473⁢(5) and −1218⁢(5)⁢keV are both sizable when compared to systematics; in the latter case the result would, in fact, be among the largest reported to date in the nuclear chart or suggest the potential existence of an additional, hitherto unidentified, low-lying 0 + state of 32 Si.

20 ≤ A ≤ 38↗
DOE OSTI2022

Isomers in Tl 203 and core excitations built on a five-nucleon-hole structure

Isomers with three- and five-nucleon-hole configurations have been established in 203 Tl These include newly identified levels with a three-nucleon-hole structure: I π = (15/2 – ) with T 1/2 = 7.9(5) ns and I π = (35/2 – ) with T 1/2 = 4.0(5) ns. In addition, five-nucleon-hole states have also been established: I π = (39/2 – ) with T 1/2 = 1.9(2) ns and I π = (49/2 + ) with T 1/2 = 3.4(4) ns. The previously determined long-lived decay, T 1/2 = 6.6(3) μs from this work, is associated with isomerism of the I π = (29/2 + ) state. Levels above this long-lived isomer have been identified through a delayed-prompt γ–γ coincidence measurement. Five-nucleon-hole states with excitation energies E x ≈ 7 MeV have been established as well as possible octupole excitations of the 208Pb core built on these levels. Furthermore, the level scheme of 203 Tl is extended up to E x ≈ 11 MeV with the inclusion of 25 new transitions. Empirical and shell-model calculations have been performed to aid in the description of the observed states which are found to be of intrinsic character.

190 ≤ A ≤ 219↗
DOE OSTI2022

β − decay of exotic P and S isotopes with neutron number near 28

β - decay of very neutron-rich isotopes of P and S, studied at the National Superconducting Cyclotron Laboratory (NSCL) using the Beta Counting Station (BCS) consisting of a Double Sided Strip Detector surrounded by clovers detectors for observing delayed γ transitions is reported here. β-decay half-lifes and delayed neutron emission probabilities were extracted for 42,43,44 P and 44,46 S by analyzing spatial and temporal correlations between implants and decay events in the Si detector with further coincidence with γ transitions. Detection of delayed γ rays allowed for the identification of negative parity 1p1h states in 42 S for the first time, also constraining the parent ( 42 P) spin/parity to 2 - or 3 - . For the most exotic isotope studied, 46 S, no strong γ transition was observed unlike lighter even-even S isotopes, thus implying the shift of Gamow Teller (GT) strength distribution to higher energies. In conclusion, comparison of experimental observations to detailed shell-model calculations using the SDPFSDG-MU interaction allowed us to infer the importance and role of First Forbidden (FF) β transitions as the neutron number approaches and then exceeds N = 28.

39 ≤ A ≤ 58↗
DOE OSTI2023

Nuclear-level effective theory of 𝜇→𝑒 conversion: Formalism and applications

Over the next decade new 𝜇→𝑒 conversion searches at Fermilab (Mu2e) and J-PARC (COMET, DeeMe) are expected to advance limits on charged lepton flavor violation (CLFV) by more than four orders of magnitude. Here, by considering the consequence of 𝑃 and 𝐶⁢𝑃 on elastic 𝜇→𝑒 conversion and the structure of possible charge and current densities, we show that rates are governed by six nuclear responses and a single scale, 𝑞/𝑚 𝑁 , where 𝑞 ≈ 𝑚 𝜇 is the momentum transferred from the leptons to the nucleus. To relate this result to microscopic formulations of CLFV, we construct in nonrelativistic effective theory (NRET) the CLFV nucleon-level interaction, pointing out the relevance of the dimensionless scales 𝑦=($^{𝑞⁢𝑏}_2$) 2 ⁢ > |$\overrightarrow{𝑣}_N$⁢|>|$\overrightarrow{𝑣}_μ$⁢|>|$\overrightarrow{𝑣}_T$| , where 𝑏 is the nuclear size, $\overrightarrow{𝑣}_N$ and $\overrightarrow{𝑣}_μ$ are the nucleon and muon intrinsic velocities, and $\overrightarrow{𝑣}_T$ is the target recoil velocity. We discuss previous work, noting the lack of a systematic treatment of the various small parameters. Because the parameter 𝑦 is not small, a proper calculation of 𝜇→𝑒 conversion requires a full multipole expansion of the nuclear response functions, an apparently daunting task with Coulomb-distorted electron partial waves. We demonstrate that the multipole expansion can be carried out to high precision by introducing a simplifying local momentum 𝑞 eff for the electron. Previous work has been limited to simple charge or spin interactions, thereby treating the nucleus effectively as a point particle. We show that such formulations are not compatible with the general form of the 𝜇→𝑒 conversion rate, failing to generate three of the six allowed nuclear response functions. The inclusion of the nucleon velocity $\overrightarrow{𝑣}_N$ yields an NRET with 16 operators and a rate of the general form. Consequently, in the current discovery era for CLFV, it provides the most sensible starting point for experimental analysis, defining what can and cannot be determined about CLFV from the highly exclusive process of 𝜇→𝑒 conversion. Finally, we expand the NRET operator basis to account for the effects of $\overrightarrow{𝑣}_μ$, associated with the muon's lower component, generating corrections to the CLFV coefficients of the point-nucleus response functions. Using advanced shell-model methods, we compute 𝜇→𝑒 conversion rates for a series of experimental targets, deriving bounds on the coefficients of the CLFV operators. These calculations are the first to include a general basis of CLFV operators, full evaluation of the associated nuclear response functions, and an accurate treatment of electron and muon Coulomb effects. We discuss target selection as an experimental “knob” that can be turned to probe the microscopic origins of CLFV. We describe two types of coherence that enhance certain CLFV operators and selection rules that blind elastic 𝜇→𝑒 conversion to others. We discuss the matching of the NRET onto higher level effective field theories, such as those constructed at the light quark level, noting opportunities to build on existing work in direct detection of dark matter. We discuss the relation of 𝜇→𝑒 conversion to 𝜇→𝑒+𝛾 and 𝜇→3⁢𝑒, showing how MEG II and Mu3e results will complement those of Mu2e and COMET. Finally we describe a accompanying script—in Mathematica and Python versions—that can be used to compute 𝜇→𝑒 conversion rates in various nuclear targets for the full set of NRET operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Onset of band structure in 70 Ga

Excited states in the odd-odd nucleus 70 Ga were studied using the 62 Ni( 14 C, αpn) fusion-evaporation reaction at the John D. Fox Superconducting Accelerator Laboratory at Florida State University with a 50 MeV beam. Depopulating γ rays were measured in coincidence with a Compton-suppressed Ge array consisting of three Clover detectors and seven single-crystal detectors. Analysis of γ–γ coincidences led to the addition of 16 new transitions to the 70 Ga level scheme, including transitions associated with the onset of a new positive-parity band structure likely based on the πg 9/2 ⊗ νg 9/2 configuration. Spins and parities were assigned using directional correlation of oriented nuclei (DCO) ratios and linear polarization measurements. Shell-model calculations with the JUN45 effective interaction reproduce the excitation energies of the positive-parity states reasonably well but predict negative-parity states that are typically 400–500 keV lower than observed. Furthermore, total Routhian surface calculations for the lowest positive-parity configuration with signature α = 0 indicate significant triaxiality (γ ≈ −20°) and the development of competing non-collective excitations at J = 8, corresponding to the experimentally observed onset of the positive-parity band. Similar features are found in the calculated surfaces for the lowest negative-parity states with α = 0.

59 ≤ A ≤ 89↗
DOE OSTI2021

Nucleon-pair approximation for nuclei from spherical to deformed regions

Here we model low-lying states of atomic nuclei in the nucleon-pair approximation of the shell model, using three approaches to select collective nucleon pairs: the generalized seniority scheme, the conjugate gradient method, and the Hartree-Fock approach. We find the collective pairs obtained from the generalized seniority scheme provide a good description for nearly spherical nuclei, and those from the conjugate gradient method or the Hartree-Fock approach work well for transitional and deformed nuclei. Our NPA calculations using collective pairs with angular momenta 0, 2, and 4 (denoted by S D G pairs) reproduce the nuclear shape evolution in the N = 26 isotones, Ca 46 , Ti 48 , Cr 50 , and Fe 52 , and yield good agreement with full configuration-interaction calculations of low-lying states in medium-heavy transitional and deformed nuclei: Ti 44 – 48 , Cr 48 , Cr 50 , Fe 52 , Zn 60 – 64 , Ge 64 , 66 , Mo 84 , and Xe 108 – 112 . Finally, using the S D G I -pair approximation we describe low-lying states of Ba 112 , 114 , cases difficult to reach by conventional configuration-interaction methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Persistence of collectivity in the low-lying states of 30,31 Na inside the 𝑁 = 20 island of inversion

Near the 𝛽-stability line, nuclei with magic numbers are expected to exhibit spherical ground states. However, mass measurements of neutron-rich nuclei near the 𝑁 = 20 island of inversion have indicated an excess in the binding energies. These results, together with subsequent spectroscopic studies and theoretical works, point to the structural evolution from normal configurations to intruder-deformed configurations characterized by particle-hole excitations across the shell gap. Despite evidence of intruder-dominant ground states for the 30,31 Na isotopes in the vicinity of 𝑁 = 20, a question remains as to whether the deformation persists into the ground-state bands beyond the first excited states. Here, in this work, we report on a heavy-ion inelastic scattering measurement performed with GRETINA, the TRIPLEX device, and the S800 spectrograph to study the low-lying excitations in neutron-rich 30,31 Na . The observed gamma rays for the (4 + ) → (3 + ), (3 + ) → 2$^+_{g.s.}$, and (4 + ) → 2$^+_{g.s.}$ decays in 30 Na and for the (7/2 + ) → (5/2 + ) and (5/2 + )→3/2$^+_{g.s.}$ decays in 31 Na were used to extract the reduced transition probabilities 𝐵⁡(𝐸⁢2↑). We report the first measurement of the 𝐵⁡(𝐸⁢2↑) between the 3/2$^+_{g.s.}$ and the (7/2 + ) states in 31 Na . The results from this study are compared to shell-model calculations, confirming the persistence of large collectivity in the low-lying excited states consistent with the formation of the well-deformed ground-state bands in 30,31 Na .

collective models↗
DOE OSTI2020

Shell Element Material Model Verification Problems for DYNA3D: Part II

A suite of supplementary shell element material model verification tests were developed for the explicit finite element program DYNA3D, in continuation of the work performed in LLNL-TR-792469 (Shell Element Material Model Verification Problems for DYNA3D). The testing procedure developed in the preceding report is extended, and used to verify all remaining untested shell model features and inputs. A collection of 78 feature-specific verification tests are proposed, some of which are applicable to multiple material models, though many are specialized to a particular model. The collective suite of tests cover all 27 currently available material models in DYNA3D. In the course of developing the proposed test suite, 24 separate shell model-related issues were identified and have subsequently been resolved.

42 ENGINEERING↗
DOE OSTI2025

Puzzling 𝐵⁡(𝐸⁢2;0 + → 2 + ) strength in the proton dripline nucleus 36 Ca

Recent measurements of the 𝐸⁢2 transition rate from the ground state to the first 2 + excited state of the proton dripline nucleus 36 Ca show an unusual pattern when compared to its isotopic neighbor 38 Ca: despite having a higher 𝐸 𝑥 ⁡(2$^{+}_{1}$) excitation energy, the 𝐵⁡(𝐸⁢2;0$^{+}_{1}$ → 2$^{+}_{1}$) rate in 36 Ca is larger. The question that naturally arises is to what extent this observation can be attributed to the unbound character of the 2$^{+}_{1}$ state. To understand the influence of the continuum space on the low-energy properties of 36 Ca, we carried out Gamow shell model calculations that can account for the continuum coupling effects associated with the occupation of unbound 𝑓⁡𝑝 shells. We found that in the threshold 2 + state, 36 Ca is spatially diffused, which impacts the observed 𝐵⁡(𝐸⁢2) trend.

20 ≤ A ≤ 38↗
DOE OSTI2022

α-cluster structure of 18 Ne

In this work, we investigate α-clustering in 18 Ne and compare it with what is known about clustering in the mirror nucleus 18 O. The excitation function for resonant elastic scattering of 14 O + α was measured in inverse kinematics using the TexAT active-target detector system, covering excitation energies from 8 to 17 MeV. The data were analyzed using a multi-channel R-matrix approach. Detailed spectroscopic information—including excitation energies, spin and parity assignments, and partial α and total widths—was extracted from the R-matrix analysis and compared with theoretical models and previously available data. Correspondence between levels in 18 O and 18 Ne is established. An extensive shell-model analysis of the 18 O– 18 Ne mirror pair shows very good agreement between theory and experiment, especially in the distribution of clustering strength. The comparison reveals that certain states, particularly at higher excitation energies, exhibit significantly stronger α-cluster character than predicted. This suggests that these states have a collective structure aligned with the corresponding α-reaction channel.

6 ≤ A ≤ 19↗
DOE OSTI2025

Quantum magic and multipartite entanglement in the structure of nuclei

Motivated by the Gottesman-Knill theorem, we present a detailed study of the quantum complexity of p -shell and s d -shell nuclei. Valence-space nuclear shell-model wave functions generated by the code are mapped to qubit registers using the Jordan-Wigner mapping (12 qubits for the p shell and 24 qubits for the s d shell), from which measures of the many-body entanglement ( n -tangles) and magic (nonstabilizerness) are determined. While exact evaluations of these measures are possible for nuclei with a modest number of active nucleons, Monte Carlo simulations are required for the more complex nuclei. The broadly applicable Pauli-string I ̂ Z ̂ exact (PSIZe) Markov chain Monte Carlo (MCMC) technique is introduced to accelerate the evaluation of measures of magic in deformed nuclei (with hierarchical wave functions), by factors of ≈ 8 for some nuclei. Significant multinucleon entanglement is found in the s d shell, dominated by proton-neutron configurations, along with significant measures of magic. This is evident not only for the deformed states, but also for nuclei on the path to instability via regions of shape coexistence and level inversion. These results indicate that quantum-computing resources will accelerate precision simulations of such nuclei and beyond. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Proton and neutron contributions to the quadrupole transition strengths in 39 Ca and 39 K studied by lifetime measurements of mirror transitions

The E2 transition matrix elements of isobaric multiplets are expected to follow a linear trend as a function of isospin projection. However, measurements of the 2 + → 0 + transitions in the A = 38 triplet of Ca, K, and Ar show a deviation from this trend with an enhanced transition strength in 38 Ca with respect to its mirror 38 Ar. We have studied analogue 11/2 – → 7/2 – E2 transitions in 39 Ca and its mirror partner 39 K to determine if this enhancement persists in neighboring Ca isotopes. Recoil-distance lifetime measurements of 39 Ca and 39 K were performed utilizing a 42 Sc secondary beam, the TRIPLEX plunger, the GRETINA array, and the S800 spectrograph. Our data provide a lifetime measurement of the (11/2 – ) state in 39 Ca as well as an improved lifetime result for the (9/2 – ) state, while the 39 K data are used to validate the present analysis. Furthermore, a comparison of the present data to shell-model calculations suggests an enhanced transition strength in 39 Ca, pointing to both proton and neutron contributions to core excitations across the Z = N = 20 shell gaps in close proximity to 40 Ca.

39 ≤ A ≤ 58↗
DOE OSTI2023

Projection algorithm for state preparation on quantum computers

We present an efficient method to prepare states of a many-body system on quantum hardware, first isolating individual quantum numbers and then using time evolution to isolate the energy. Our method in its simplest form requires only one additional auxiliary qubit. The total time evolved for an accurate solution is proportional to the ratio of the spectrum range of the trial state to the gap to the lowest excited state, a substantial improvement over other projection algorithms, and the accuracy increases exponentially with the time evolved. Isolating the quantum numbers is efficient because of the known eigenvalues and increases the gap thus shortening the propagation time required. The success rate of the algorithm, or the probability of producing the desired state, is a simple function of measurement times and phases and is dominated by the square overlap of the original state to the desired state. We present examples from the nuclear shell model and the Heisenberg model. In conclusion, we compare this algorithm to previous algorithms for short evolution times and discuss potential further improvements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Evolution of low-lying M1 modes in germanium isotopes

Magnetic dipole strength functions are determined for the series of germanium isotopes from N=Z=32 to N = 48 on the basis of a large number of transition strengths calculated within the shell model. The evolution of the strength with increasing neutron number in the 1g 9/2 orbital is analyzed. A bimodal structure comprising an enhancement toward low transition energy and a resonance in the region of the scissors mode is identified. The low-energy enhancement is strongest near closed shells, in particular at the almost completely filled 1g 9/2 orbital, while the scissorslike resonance is most pronounced in the middle of the open shell, which correlates with the magnitude of the also deduced electric quadrupole transition strengths. The results are consistent with previous findings for the shorter series of iron isotopes and prove the occurrence and correlation of the two low-lying magnetic dipole modes as a global structural feature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗