Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Nonempirical Interactions for the Nuclear Shell Model: An Update

The nuclear shell model has perhaps been the most signifcant conceptual and computational paradigm for the understanding of the structure of atomic nuclei. While the shell model has been used predominantly in a phenomenological context, there have been efforts stretching back more than half a century to derive shell model parameters based on a realistic interaction between nucleons. More recently, several ab initio many-body methods—in particular, many-body perturbation theory, the no-core shell model, the in-medium similarity renormalization group, and coupled-cluster theory—have developed the capability to provide effective shell model Hamiltonians. We offer an update on the status of these methods and investigate the connections between them and their potential strengths and weaknesses, with a particular focus on the in-medium similarity renormalization group approach. Three-body forces are demonstrated to be important for understanding the modifications needed in phenomenological treatments. We then review some applications of these methods to comparisons with recent experimental measurements, and conclude with some remaining challenges in ab initio shell model theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Proton-neutron entanglement in the nuclear shell model

Abstract We compute the proton-neutron entanglement entropy in the interacting nuclear shell model for a variety of nuclides and interactions. Some results make intuitive sense, for example, that the shell structure, as governed by single-particle and monopole energies, strongly affects the energetically available space and thus the entanglement entropy. We also find a surprising result: that the entanglement entropy at low excitation energy tends to decrease for nuclides when N ≠ Z . While we provide evidence this arises from the physical nuclear force by contrasting with random two-body interactions which shows no such decrease, the exact mechanism is unclear. Nonetheless, the low entanglement suggests that in models of neutron-rich nuclides, the coupling between protons and neutrons may be less computationally demanding than one might otherwise expect.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2020

The convergence of quadrupole rotational invariants from the nuclear shell model

Nuclei exhibit both single-particle and collective degrees of freedom, with the latter often subdivided into vibrational and rotational motions. Experimentally identifying the relative roles of these collective modes is extremely challenging, particularly in the face of possible shape coexistence. Model-independent, invariant quantities describing the deformation of a nucleus in the intrinsic frame have long been known but their determination potentially requires a large quantity of experimental data to achieve convergence. Through comparison with the nuclear shell model, the question of convergence is addressed. Shell-model calculations performed in the sd- and pf-shell model spaces are used to determine electric-quadrupole matrix elements for a multitude of low-lying states using the first 40 states of the relevant spins. Relative contributions to the rotationally invariant quantities from multiple states can therefore be determined. It is found that, on average, the inclusion of four intermediate states results in the leading-order invariant, $\langle\hat{Q^2}\rangle$, converging to within 10% of its true value and the triaxiality term, cos (3δ), converging to its true value, though some variance remains. Higher-order quantities relating to the softness of the nuclear shape are found to converge more slowly. The convergence of quadrupole rotationally invariant sum rules was quantified in the sd- and pf-shell model spaces and indicates the challenge inherent in a full determination of nuclear shape. The present study is limited to relatively small valence spaces. Finally, larger spaces, such as the rare-earth region, potentially offer faster convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Present Status of Nuclear Shell-Model Calculations of 0 νββ Decay Matrix Elements

Neutrinoless double beta (0νββ) decay searches are currently among the major foci of experimental physics. The observation of such a decay will have important implications in our understanding of the intrinsic nature of neutrinos and shed light on the limitations of the Standard Model. The rate of this process depends on both the unknown neutrino effective mass and the nuclear matrix element (M 0ν ) associated with the given 0νββ transition. The latter can only be provided by theoretical calculations, hence the need of accurate theoretical predictions of M 0ν for the success of the experimental programs. This need drives the theoretical nuclear physics community to provide the most reliable calculations of M 0ν . Among the various computational models adopted to solve the many-body nuclear problem, the shell model is widely considered as the basic framework of the microscopic description of the nucleus. Here, we review the most recent and advanced shell-model calculations of M 0ν considering the light-neutrino-exchange channel for nuclei of experimental interest. We report the sensitivity of the theoretical calculations with respect to variations in the model spaces and the shell-model nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Solving nuclear structure problems with the adaptive variational quantum algorithm

Here we use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the sd and pf shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2023

CoSMo - Shell Model Code

CoSMo is a comprehensive shell model code suite designed for nuclear shell model and configuration interaction calculations. Its highly structured and templated nature allows for flexibility and ease in applications, including those to open quantum systems with non-Hermitian Hamiltonians, clustering, and time-dependent dynamics.

Volya, Alexander↗
DOE OSTI2021

Entanglement rearrangement in self-consistent nuclear structure calculations

Background: Entanglement plays a central role in a diverse array of increasingly important research areas, including quantum computation, simulation, measurement, sensing, and communication. Extensive suites of investigations have been performed to better understand entanglement in atomic and molecular quantum many-body systems, while the exploration of entanglement in the structure of nuclei and their reactions is presently in its infancy. Purpose: The goal of this work is to begin investigating the entanglement properties of nuclei from first-principles nuclear many-body calculations. We attempt to identify common features and emergent structures of entanglement that could ultimately lead to new and natural many-body schemes. With an eye toward quantum accelerators in future hybrid-supercomputers, criteria for partitioning nuclear many-body calculations into quantum and classical components may provide advantages in future large-scale computations. Along the way we look for explanations of the relative success of phenomenological models such as the nuclear shell model, and for better ways to match to low-energy nuclear effective field theories and lattice QCD calculations to nuclear many-body techniques that are based upon entanglement. Method: We explore the entanglement between single-particle states in 4 He and 6 He. The patterns of entanglement emerging from different single-particle bases are compared, and possible links with the convergence of observables are explored, in particular, ground-state energies. The nuclear wave functions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases exhibiting different degrees of complexity are determined, in particular, harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases. Specifically, single-orbital entanglement entropy, two-orbital mutual information, and negativity are studied. Results: The entanglement structures in 4 He and 6 He are found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters we have worked with. In particular a core-valence structure clearly emerges from the full no-core calculation of 6 He. Here, the two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. Conclusions: Measures of one- and two-nucleon entanglement are found to be useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Time-dependent relaxation of observables in complex quantum systems

Abstract We consider time-dependent relaxation of observables in quantum systems of chaotic and regular type. Using statistical arguments and exact numerical solutions we show that the spread of the initial wave function in the Hilbert space and the main characteristics of evolution of observables have certain generic features. The study compares examples of regular dynamics, a completely chaotic case of the Gaussian orthogonal ensemble, a bosonic system with random interactions, and a fully realistic case of the time evolution of various initial non-stationary states in the nuclear shell model. In the case of the Gaussian orthogonal ensemble we show that the survival probability obtained analytically also fully defines the relaxation timescale of observables. This is not the case in general. Using the realistic nuclear shell model and the quadrupole moment as an observable we demonstrate that the relaxation time is significantly longer than defined by the survival probability of the initial state. The full analysis does not show the presence of an analog of the Lyapunov exponent characteristic for examples of classical chaos.

Volya, Alexander (ORCID:0000000217656466)↗
DOE OSTI2019

Is the Structure of Si 42 Understood?

A more detailed test of the implementation of nuclear forces that drive shell evolution in the pivotal nucleus 42 Si—going beyond earlier comparisons of excited-state energies—is important. The two leading shell-model effective interactions, SDPF-MU and SDPF-U-Si, both of which reproduce the low-lying 42 Si(2+1) energy, but whose predictions for other observables differ significantly, are interrogated by the population of states in neutron-rich 42 Si with a one-proton removal reaction from 43 P projectiles at 81 MeV/nucleon. The measured cross sections to the individual 42 Si final states are compared to calculations that combine eikonal reaction dynamics with these shell-model nuclear structure overlaps. The differences in the two shell-model descriptions are examined and linked to predicted low-lying excited 0 + states and shape coexistence. Based on the present data, which are in better agreement with the SDPF-MU calculations, the state observed at 2150(13) keV in 42 Si is proposed to be the second 0$ ^+ _2 $ level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Exact isovector pairing in a shell-model framework: Role of proton-neutron correlations in isobaric analog states

Here, we utilize a nuclear shell model Hamiltonian with only two adjustable parameters to generate, for the first time, exact solutions for pairing correlations for light to medium-mass nuclei, including the challenging proton-neutron pairs, while also identifying the primary physics involved. In addition to single-particle energy and Coulomb potential terms, the shell model Hamiltonian consists of an isovector T = 1 pairing interaction and an average proton-neutron isoscalar T = 0 interaction, where the T = 0 term describes the average interaction between nonpaired protons and neutrons. This Hamiltonian is exactly solvable, where, utilizing three to seven single-particle energy levels, we reproduce experimental data for 0 + state energies for isotopes with mass A = 10 through A = 62 exceptionally well including isotopes from He to Ge. Additionally, we isolate effects due to like-particle and proton-neutron pairing, provide estimates for the total and proton-neutron pairing gaps, and reproduce N (neutron) = Z (proton) irregularity. These results provide a further understanding for the key role of proton-neutron pairing correlations in nuclei, which is especially important for waiting-point nuclei on the $\mathcal{rp}$ path of nucleosynthesis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Search for Weak Side Branches in the Electromagnetic Decay Paths of the 6526-keV 10 + Isomer in 54 Fe

High-spin nuclear isomers in $N≈Z$ nuclei between doubly magic 40 Ca and 56 Ni provide an excellent testing ground for the nuclear shell model and questions related to isospin symmetry breaking in the vicinity of the proton drip line. The purpose of the present study is to investigate the possibility of weak electromagnetic decay branches along the decay paths of the 6526-keV 10 + isomer in 54 Fe. The isomer was strongly populated by means of the fusion-evaporation reaction 24 Mg( 36 Ar,$α$2$p$) 54$m$ Fe. The Gammasphere array was used to detect $γ$-ray cascades emitted from the isomeric state. By means of $γγγ$ coincidences, weak non-yrast decay branches can be discriminated, with the isomer’s half-life confirmed at $T$ 1/2 = 363 (4) ns. The yrast ${6}_{1}^{+}$ → ${2}_{1}^{+}$ $E4$ cross-over transition was interrogated. The observations are compared with shell-model calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Many-body factorization and position–momentum equivalence of nuclear short-range correlations

While mean-field approximations, such as the nuclear shell model, provide a good description of many bulk nuclear properties, they fail to capture the important effects of nucleon–nucleon correlations such as the short-distance and high-momentum components of the nuclear many-body wave function1. Here, we study these components using the effective pair-based generalized contact formalism and ab initio quantum Monte Carlo calculations of nuclei from deuteron to 40 Ca. We observe a universal factorization of the many-body nuclear wave function at short distance into a strongly interacting pair and a weakly interacting residual system. The residual system distribution is consistent with that of an uncorrelated system, showing that short-distance correlation effects are predominantly embedded in two-body correlations. Spin- and isospin-dependent ‘nuclear contact terms’ are extracted in both coordinate and momentum space for different realistic nuclear potentials. The contact coefficient ratio between two different nuclei shows very little dependence on the nuclear interaction model. These findings thus allow extending the application of mean-field approximations to short-range correlated pair formation by showing that the relative abundance of short-range pairs in the nucleus is a long-range (that is, mean field) quantity that is insensitive to the short-distance nature of the nuclear force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗