No-core Shell Model Research Articles

Uncertainty quantification of collective nuclear observables from the chiral potential parametrization

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

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

<i>Ab initio</i> leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Ab initio study of Z(N) = 6 magicity* *This work is partially supported by the National Natural Science Foundation of China (Grant Nos. 12175280, 12250610193, 12375143, 11975282, 11705240, 11435014), the Natural Science Foundation of Gansu Province, China (Grant Nos. 20JR10RA067, 23JRRA675), the Chinese Academy of Sciences "Light of West China" Program, the Key Research Program of the Chinese Academy of

The existence of magic numbers of protons and neutrons in nuclei is essential for understanding the nuclear structure and fundamental nuclear forces. Over decades, researchers have conducted theoretical and experimental studies on a new magic number, Z(N)=6, focusing on observables such as radii, binding energy, electromagnetic transition, and nucleon separation energies. We performed ab initio no-core shell model calculations for the occupation numbers of the lowest single particle states in the ground states of Z(N)=6 and Z(N)=8 isotopes (isotones). The results of our calculations do not support Z(N)=6 as a magic number over a range of atomic numbers. However, 14C and 14O exhibit the characteristics of double-magic nuclei.

Direct Observation of Competing M1 and M3 Transitions in ^{10}B.

Excited states in ^{10}B were populated with the ^{10}B(p,p^{'}γ)^{10}B^{*} reaction at 8.5MeV and their γ decay was investigated via coincidence γ-ray spectroscopy. The emitted γ rays were measured using large-volume LaBr_{3}:Ce and CeBr_{3} detectors placed in anti-Compton shields. This allowed the observation of weak γ-ray transitions, such as the M3 transition between the J^{π},T=0^{+},1 isobaric analog state (IAS) and the J^{π},T=3^{+},0 ground state and the E2 transition between the J^{π},T=2_{1}^{+},0 state and the IAS, i.e., performing measurements of branching ratios at the level of λ≥10^{-4}. For the first time in ^{10}B, the competing M1 and M3 transitions from the decay of the IAS have been observed in a γ spectroscopy experiment. The experimental results are compared with abinitio no-core shell model calculation using the newest version of the local position-space chiral N^{3}LO nucleon-nucleon interaction. The calculations reproduce correctly the ordering of the bound states in ^{10}B, and are in reasonable agreement with the observed branching ratios and reduced transition probabilities.

Precise neural network predictions of energies and radii from the no-core shell model

Precise neural network predictions of energies and radii from the no-core shell model

Trineutron resonances in the SS-HORSE extension of the no-core shell model

Trineutron resonances in the SS-HORSE extension of the no-core shell model

Intruder band mixing in an ab initio description of 12Be

The spectrum of ▪ exhibits exotic features, e.g., an intruder ground state and shape coexistence, normally associated with the breakdown of a shell closure. While previous phenomenological treatments indicated the ground state has substantial contributions from intruder configurations, it is only with advances in computational abilities and improved interactions that this intruder mixing is observed in ab initio no-core shell model (NCSM) predictions. In this work, we extract electromagnetic observables and symmetry decompositions from the NCSM wave functions to demonstrate that the low-lying positive parity spectrum can be explained in terms of mixing of rotational bands with very different intrinsic structure coexisting within the low-lying spectrum. These observed bands exhibit an approximate SU(3) symmetry and are qualitatively consistent with Elliott model predictions.

Machine Learning in the Problem of Extrapolating Variational Calculations in Nuclear Physics

A modified machine learning method is proposed, utilizing an ensemble of artificial neural networks for the extrapolation of energies obtained in variational calculations, specifically in the No-core Shell Model (NCSM), to the case of the infinite basis. A new neural network topology is employed, and criteria for selecting both the data used for training and the trained neural networks for statistical analysis of the results are formulated. The approach is tested by extrapolating the deutron ground state energy in calculations with the Nijmegen II NN interaction and provides statistically significant results. This technique is applied to obtain extrapolated ground state energies of 6Li and 6He nuclei based on the NCSM calculations with Daejeon16 NN interaction.

Direct ab initio calculation of the 4He nuclear electric dipole polarizability

The calculation of nuclear electromagnetic sum rules by directly solving for numerous eigenstates in a large basis is numerically challenging and has not been performed for A>2 nuclei. With the significant progress of high performance computing, we show that calculating sum rules using numerous discretized continuum states obtained with the ab initio no-core shell model in the harmonic oscillator basis is achievable numerically. Specifically, we calculate the 4He electric dipole (E1) polarizability, that is an inverse energy weighted sum rule, employing the Daejeon16 NN interaction. We demonstrate that the calculations are numerically tractable as the dimension of the basis increases and are convergent. Our results for the 4He electric dipole polarizability are consistent with the most recent experimental data and are compared with those of other theoretical studies employing different techniques and various interactions.

Separation energies of light Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varLambda $$\\end{document} hypernuclei and their theoretical uncertainties

Separation energies of light Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varLambda $$\\end{document} hypernuclei (A≤5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$A\\le 5$$\\end{document}) and their theoretical uncertainties are investigated. Few-body calculations are performed within the Faddeev-Yakubovsky scheme and the no-core shell model. Thereby, modern and up-to-date nucleon-nucleon, three-nucleon and hyperon-nucleon potentials derived within chiral effective field theory are employed. It is found that the numerical uncertainties of the few-body methods are well under control and an accuracy of around 1 keV for the hypertriton and of better than 20 keV for the separation energies of the Λ4He\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{\\,4}_{\\varLambda }\ extrm{He}$$\\end{document} and Λ5He\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{\\,5}_{\\varLambda }\ extrm{He}$$\\end{document} hypernuclei can be achieved. Variations caused by differences in the nucleon-nucleon interaction are in the order of 10 keV for Λ3H\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{\\,3}_{\\varLambda }\ extrm{H}$$\\end{document} and no more than 110 keV for A=4,5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$A=4,\\,5$$\\end{document}Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varLambda $$\\end{document} hypernuclei, when recent high-precision potentials up to fifth order in the chiral expansion are employed. The variations are smaller than the expected contributions from chiral hyperon-nucleon-nucleon forces which arise at the chiral order of state-of-the-art hyperon-nucleon potentials. Estimates for those three-body forces are deduced from a study of the truncation uncertainties in the chiral expansion.

New studies and predictions for nuclear clustering and dynamics in the ab initio symmetry-adapted framework

We discuss recent studies and predictions for nuclear clustering and dynamics within the framework of the ab initio symmetry-adapted no-core shell model, which has opened new domains of the nuclear chart. In this framework, we show the emergence from first principles of collectivity and clustering in light to mediummass nuclei, with implications for constructing ab initio optical potentials, for studying clustering in stable and unstable nuclei, for reproducing enhanced deformations without effective charges, and for the formation of clusters and its sensitivity to the underlying inter-nucleon force.

γ-ray spectroscopic study of self-conjugate 10B nucleus with inelastic proton scattering

A γ-ray spectroscopic study of 10B is reported. Excited states in 10B were populated in inelastic proton scattering and their γ-decays were measured with the ELIFANT array, a state-of-the-art spectrometer with large volume LaBr3:Ce and CeBr3 detectors placed in anti-Compton shields and having unprecedented efficiency for high-energy γ rays. A few weak transitions were observed in the experiment, namely M3 transition between Jπ, T = 0+1, 1 isobaric analog state to the Jπ, T = 3+1,0 ground state and E2 transition between the Jπ, T = 2+1,0 state and the isobaric analog state. The results are compared to a new set of ab initio no-core shell model calculation using the newest version of the local position-space chiral N3LO nucleon-nucleon interaction, which correctly reproduces the spectrum of the excited states in 10B, a problem of the theory in previous years.

Description of the Proton-Decaying 0_{2}^{+} Resonance of the α Particle.

The recent precise experimental determination of the monopole transition form factor from the ground state of ^{4}He to its 0_{2}^{+} resonance via electron scattering has reinvigorated discussions about the nature of this first excited state of the α particle. The 0_{2}^{+} state has been traditionally interpreted in the literature as the isoscalar monopole resonance (breathing mode) or, alternatively, as a particle-hole shell-model excitation. To better understand the nature of this state, which lies only ∼410 keV above the proton emission threshold, we employ the coupled-channel representation of the no-core Gamow shell model. By considering the [^{3}H+p], [^{3}He+n], and [^{2}H+^{2}H] reaction channels, we explain the excitation energy and monopole form factor of the 0_{2}^{+} state. We argue that the continuum coupling strongly impacts the nature of this state, which carries characteristics of the proton decay threshold.

Hyperon-nucleon interaction constrained by light hypernuclei

Ab initio structure calculations for p-shell hypernuclei have recently become accessible through extensions of nuclear many-body methods, such as the no-core shell model, in combination with hyperon-nucleon interactions from chiral effective field theory. However, the low-energy constants in these hyperon-nucleon interactions are poorly constraint due to the very limited amount of experimental scattering data available. We present a hyperon-nucleon interaction that is additionally constrained by experimental ground-state and spectroscopic data for selected p-shell hypernuclei and, thus, optimized for hypernuclear structure calculations. We show that the previous overestimation of the hyperon separation energies in the p-shell is remedied and discuss the significantly improved description of the ▪ isotopic chain. We further discuss the uncertainty quantification for hypernuclear observables on the many-body level, obtained through a novel machine-learning tool.

Ab initio investigations of A=8 nuclei: α–α scattering, deformation in 8He, radiative capture of protons on 7Be and 7Li and the X17 boson

We apply the No-Core Shell Model with Continuum (NCSMC) that is 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. The NCSMC can predict structure and dynamics of light nuclei and, by comparing to available experimental data, test the quality of chiral nuclear forces. We discuss applications of NCSMC to the α–α scattering and the structure of 8Be, the p+7Be and p+7Li radiative capture and the production of the hypothetical X17 boson claimed in ATOMKI experiments. The 7Be(p, γ)8B reaction plays a role in Solar nucleosynthesis and Solar neutrino physics and has been subject of numerous experimental investigations. We also highlight our investigation of the neutron rich exotic 8He that has been recently studied experimentally at TRIUMF with an unexpected deformation reported.

Sizes of the Neutron–Proton Halo in Nucleon-Stable States of the 6Li Nucleus

The matter, neutron, and proton radii of nucleon-stable 1+ and 0+ states of the 6Li nucleus are studied theoretically within the no-core shell model. The results have been comparatively analyzed with the radii of the 0+ state of the 6He nucleus. To increase the accuracy of calculations, we have developed an extrapolation procedure. A new definition of the quantitative measure has been proposed and justified to describe the properties of a halo formed by loosely bound neutrons and protons in the A-nucleon problem. The sizes of the halo in the indicated states of 6Li have been calculated for the first time. It has been demonstrated that the sizes of their halos are close to those of the two-neutron halo in 6He. Thus, additional reliable evidence of the existence of the neutron–proton halo in the discussed states of 6Li has been obtained.

NuHamil : A numerical code to generate nuclear two- and three-body matrix elements from chiral effective field theory

The applicability of nuclear ab initio calculations has rapidly extended over the past decades. However, starting research projects is still challenging due to the required numerical expertise in the generation of underlying nuclear interaction matrix elements and many-body calculations. To ease the first issue, in this paper we introduce the numerical code NuHamil to generate the nucleon-nucleon (NN) and three-nucleon (3N) matrix elements expressed in a spherical harmonic-oscillator basis, inputs of many-body calculations. The ground-state energies for the selected doubly closed shell nuclei are calculated with the no-core shell-model (NCSM) and in-medium similarity renormalization group (IMSRG). The code is written in modern Fortran, and OpenMP+MPI hybrid parallelization is available for the 3N matrix-element calculations.

Parallel multithreaded deduplication of data sequences in nuclear structure calculations

High performance computing (HPC) applications that work with redundant sequences of data can benefit from their deduplication. We study this problem on the symmetry-adapted no-core shell model (SA-NCSM), where redundant sequences of different kinds naturally emerge in the data of the basis of the Hilbert space physically relevant to a modeled nucleus. For a fast solution of this problem on multicore architectures, we propose and present three multithreaded algorithms, which employ either concurrent hash tables or parallel sorting methods. Furthermore, we present evaluation and comparison of these algorithms based on experiments performed with real-world SA-NCSM calculations. The results indicate that the fastest option is to use a concurrent hash table, provided that it supports sequences of data as a type of table keys. If such a hash table is not available, the algorithm based on parallel sorting is a viable alternative.

Study of S, Cl and Ar isotopes with $$N \ge Z$$ using microscopic effective sd-shell interactions

In the present work, newly developed microscopic effective $sd$-valence shell interactions such as chiral next-to-next-to-next-to-leading order (N3LO), $J$-matrix inverse scattering potential (JISP16), Daejeon16 (DJ16), and monopole-modified DJ16 (DJ16A) are employed to study the nuclear structural properties of sulphur, chlorine, and argon isotopes with $N \geq Z$. These interactions are derived using the \textit{ab initio} no-core shell-model and the OLS unitary transformation method. We calculate energy spectra and electromagnetic properties to test the predictive strength of the effective interactions for these heavier $sd$-shell nuclei. For a complete systematic study, we compare the microscopic results with the phenomenological USDB results and experimental data. By looking at the excitation energies of these nuclei, the DJ16A interaction is found to be {the} most suitable for these $sd$-shell nuclei among all microscopic interactions. The electric quadrupole transition strength and excitation energy of the first $2^+$ state data of even-even sulphur isotopes indicate the presence of the $N=20$ shell closure. Quadrupole moment predictions are also made using these interactions where experimental data are unknown. Magnetic moments are in excellent agreement with the experimental values. The root-mean-square deviations are also calculated to provide an idea of how accurate the interactions are.

Nuclear structure properties of Si and P isotopes with the microscopic effective interactions

In the present work, we have reported comprehensive study of the sd-shell nuclei Si and P with neutron number varying from N=9 to N=20 using the microscopic effective valence shell interactions, namely, N3LO, JISP16 and DJ16A. These effective sd-shell interactions are developed using the ab initio no-core shell model wave functions and the Okubo-Lee-Suzuki transformation method. For comparison, we have also performed shell model calculations with the empirical USDB interaction. Energy spectra and electromagnetic properties of these isotopic chains have been studied. Theoretically calculated shell model results are compared with the available experimental data, to check the predictive strength of these microscopic interactions. It is found that the binding energies of the ground states are better reproduced with the DJ16A interaction as compared to other microscopic interactions and a proton subshell closure at Z=14 is obtained in Si. Spin-tensor decomposition of two-body interaction is presented to understand the contributions from central, vector and tensor components into these interactions. Spectroscopic strengths of 23Al(d,n)24Si are examined for the newly performed experiment at NSCL.