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 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 OSTI2013

Shell model nuclear matrix elements for competing mechanisms contributing to double beta decay

Recent progress in the shell model approach to the nuclear matrix elements for the double beta decay process are presented. This includes nuclear matrix elements for competing mechanisms to neutrionless double beta decay, a comparison between closure and non-closure approximation for {sup 48}Ca, and an updated shell model analysis of nuclear matrix elements for the double beta decay of {sup 136}Xe.

73 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 OSTI2017

Description of resonant states in the shell model

A technique for describing scattering states within the nuclear shell model is proposed. This technique is applied to scattering of nucleons by particles based on ab initio No-Core Shell Model calculations of {sup 5}He and {sup 5}Li nuclei with JISP16 NN interaction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Level density of the sd-nuclei—Statistical shell-model predictions

Accurate knowledge of the nuclear level density is important both from a theoretical viewpoint as a powerful instrument for studying nuclear structure and for numerous applications. For example, astrophysical reactions responsible for the nucleosynthesis in the universe can be understood only if we know the nuclear level density. We use the configuration-interaction nuclear shell model to predict nuclear level density for all nuclei in the sd-shell, both total and for individual spins (only with positive parity). To avoid the diagonalization in large model spaces we use the moments method based on statistical properties of nuclear many-body systems. In the cases where the diagonalization is possible, the results of the moments method practically coincide with those from the shell-model calculations. Using the computed level densities, we fit the parameters of the Constant Temperature phenomenological model, which can be used by practitioners in their studies of nuclear reactions at excitation energies appropriate for the sd-shell nuclei.

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 OSTI2015

Relativistic nuclear theory—nucleons and mesons: origin, current status, and trends

The nuclear shell model (NSM) is a fundamental model of nuclear theory. At the initial stage, the NSM was developed on the basis of the Schrödinger equation, in particular because it was not clear which Dirac matrices should be associated with various components of the (relativistic) shell model potential. In the early 1970s, the relativistic version of the NSM with meson fields as its main components was developed on the basis of the Dirac equation. The relativistic nuclear shell model (RNSM) includes meson fields with different space–time transformation properties (scalar, four-vector, etc.) that determine the behavior of the (corresponding) meson fields under Lorentz transformations. This fact directly indicates that nuclear theory should be relativistic and based on the Dirac equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗