Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 55 records · Page 3

Machine learning approach to pattern recognition in nuclear dynamics from the ab initio symmetry-adapted no-core shell model

A novel machine learning approach is used to provide further insight into atomic nuclei and to detect orderly patterns amid a vast data of large-scale calculations. The method utilizes a neural network that is trained on ab initio results from the symmetry-adapted no-core shell model (SA-NCSM) for light nuclei. We show that the SA-NCSM, which expands ab initio applications up to medium-mass nuclei by using dominant symmetries of nuclear dynamics, can reach heavier nuclei when coupled with the machine learning approach. In particular, we find that a neural network trained on probability amplitudes for s- and p-shell nuclear wave functions not only predicts dominant configurations for heavier nuclei but in addition, when tested for the 20 Ne ground state, accurately reproduces the probability distribution. The non-negligible configurations predicted by the network provide an important input to the SA-NCSM for reducing ultralarge model spaces to manageable sizes that can be, in turn, utilized in SA-NCSM calculations to obtain accurate observables. The neural network is capable of describing nuclear deformation and is used to track the shape evolution along the 20-42 Mg isotopic chain, suggesting a shape coexistence that is more pronounced toward the very neutron-rich isotopes. We provide first descriptions of the structure and deformation of 24 Si and 40 Mg of interest to x-ray burst nucleosynthesis, and even of the extremely heavy nuclei such as 166,168 Er and 236 U, that build on first-principles considerations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Toward scalable many-body calculations for nuclear open quantum systems using the Gamow Shell Model

Drip-line nuclei have very different properties from those of the valley of stability, as they are weakly bound and resonant. Therefore, the models devised for stable nuclei can no longer be applied therein. Hence, a new theoretical tool, the Gamow Shell Model (GSM), has been developed to study the many-body states occurring at the limits of the nuclear chart. GSM is a configuration interaction model based on the use of the so-called Berggren basis, which contains bound, resonant and scattering states, so that inter-nucleon correlations are fully taken into account and the asymptotes of extended many-body wave functions are precisely handled. However, large complex symmetric matrices must be diagonalized in this framework, therefore the use of very powerful parallel machines is needed therein. So, in order to fully take advantage of their power, a 2D partitioning scheme using hybrid MPI/OpenMP parallelization has been developed in our GSM code. The specificities of the 2D partitioning scheme in the GSM framework will be described and illustrated with numerical examples. It will then be shown that the introduction of this scheme in the GSM code greatly enhances its capabilities.

2D partitioning↗

Uncertainty quantification of an empirical shell-model interaction using principal component analysis

Recent investigations have emphasized the importance of uncertainty quantification (UQ) in nuclear theory. Here, we carry out UQ for configuration-interaction shell-model calculations in the 1$\textit{s}$–0$\textit{d}$ valence space, investigating the sensitivity of observables to perturbations in the 66 parameters (matrix elements) of a high-quality empirical interaction. The large parameter space makes computing the corresponding Hessian numerically costly, so we compare a cost-effective approximation, using the Feynman-Hellmann theorem, to the full Hessian and find it works well. Diagonalizing the Hessian yields the principal components of the interaction: linear combinations of parameters ordered by sensitivity. This approximately decoupled distribution of parameters facilitates theoretical uncertainty propagation onto structure observables: electromagnetic transitions, Gamow-Teller decays, and dark matter-nucleus scattering matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutron-rich calcium isotopes within realistic Gamow shell model calculations with continuum coupling

Based on the realistic nuclear force of the high-precision CD-Bonn potential, we have performed comprehensive calculations for neutron-rich calcium isotopes using the Gamow shell model (GSM) which includes resonance and continuum. The realistic GSM calculations produce well binding energies, one- and two-neutron separation energies, predicting that 57 Ca is the heaviest bound odd isotope and 70 Ca is the dripline nucleus. Resonant states are predicted, which provides useful information for future experiments on particle emissions in neutron-rich calcium isotopes. Shell evolutions in the calcium chain around neutron numbers $\textit{N}$ = 32, 34, and 40 are understood by calculating effective single-particle energies, the excitation energies of the first 2 + states and two-neutron separation energies. The calculations support shell closures at 52 Ca ($\textit{N}$ = 32) and 54 Ca ($\textit{N}$ = 34) but show a weakening of shell closure at 60 Ca ($\textit{N}$ = 40). Finally, the possible shell closure at 70 Ca ($\textit{N}$ = 50) is predicted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Lifetime measurements of excited states in neutron-rich 53 Ti: Benchmarking effective shell-model interactions

Level lifetimes of the yrast (5/2 - ) to 13/2 - states in the neutron-rich nucleus 53 Ti, produced in a multinucleon-transfer reaction, have been measured for the first time. The recoil distance Doppler-shift method was employed and lifetimes of the excited states were extracted by a lineshape analysis aided by GEANT4-based Monte-Carlo simulations. The experiment was performed at the Grand Accelerateur National d'Ions Lourds facility in Caen, France, by using the Advanced Gamma Tracking Array for the gamma-ray detection coupled to the large-acceptance variable mode spectrometer for an event-by-event particle identification and the Cologne plunger for deep-inelastic reactions. Reduced transition probabilities, deduced from the lifetimes, give new information on the nuclear structure of 53 Ti, and are used to benchmark different shell-model calculations using established interactions in the fp shell.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab initio no-core Gamow shell-model calculations of multineutron systems

The existence of multineutron systems has always been a debatable question. Indeed, both internucleon correlations and a large continuum coupling occur in these states. We then employ the ab initio no-core Gamow shell model to calculate the resonant energies and widths of the trineutron and tetraneutron systems with realistic interactions. In this work, our results indicate that trineutron and tetraneutron are both unbound and bear broad widths. The calculated energy and width of tetraneutron are also comparable with recent experimental data. Moreover, our calculations suggest that the energy of trineutron is lower than that of tetraneutron, while its resonance width is also narrower. This strongly suggests that trineutron is more likely to be experimentally observed than tetraneutron. We thus suggest experimentalists to search for trineutron at low energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutrinoless Double-Beta Decay Investigations of 82Se Using Three Shell Model Hamiltonians

Neutrinoless double-beta decay is considered one of the most promising processes that would help clarify some of the symmetry-breaking problems in our understanding of the observable universe. Recent studies of neutrinoless double-beta decay matrix elements have employed statistical approaches based on modified shell model effective Hamiltonians for 48Ca (Phys. Rev. C 106, 054302 (2022)) and 136Xe (Phys. Rev. C 107, 045501 (2023)). The analyses rely on inducing perturbations in the starting effective Hamiltonians to observe the behavior of a wide range of observables, besides the 0νββ) NME, that are compared with experimental data. Following a Bayesian model averaging approach, the range of probable values for the neutrinoless double-beta decay matrix elements is presented. In this paper, we present a similar study for 82Se, which is described in the same model space as 76Ge that is under experimental observation. Due to its faster calculation time compared to 76Ge, 82Se can be used as an appropriate substitute in our complex statistical study. Using the calculations performed for the statistical analysis of the neutrinoless double-beta decay matrix elements, we also search for the correlations between the observables that we can compare to experimental data.

Science & Technology - Other Topics↗

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↗

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗

Gamow-shell-model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we develop the complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 ≤ $\textit{A}$ ≤ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid 4 He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the $\textit{sp}$ space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: We report our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei 10 Li – 10 N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Be 7 and Li 7 nuclei within the no-core shell model with continuum

The production of $^7Be$ and $^7Li$ nuclei plays an important role in primordial nucleosynthesis, nuclear astrophysics, and fusion energy generation. The $^3He(α,γ) ^7Be$ and $^3H(α,γ) 7Li$ radiative-capture processes are important to determine the $^7Li$ abundance in the early universe and to predict the correct fraction of $pp$ -chain branches resulting in $^7Be$ versus $^8B$ neutrinos. The $^6Li(p,γ) ^7Be$ has been investigated recently hinting at a possible cross section enhancement near the threshold. The $^6Li (n,^3H) ^4He$ process can be utilized for tritium breeding in machines dedicated to fusion energy generation through the deuteron-tritium reaction, and is a neutron cross section standard used in the measurement and evaluation of fission cross sections. In this work we study the properties of $^7Be$ and $^7Li$ within the no-core shell model with continuum (NCSMC) method, using chiral nucleon-nucleon interactions as the only input, and analyze all the binary mass partitions involved in the formation of these systems. The NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states and thus is well suited to investigate systems with many resonances and pronounced clustering like $^7Be$ and $^7Li$. Our calculations reproduce all the experimentally known states of the two systems and provide predictions for several new resonances of both parities. Some of these new possible resonances are built on the ground states of $^6Li$ and $^6He$, and thus represent a robust prediction. We do not find any resonance in the $p + ^6Li$ mass partition near the threshold. On the other hand, in the $p + ^6He$ mass partition of $^7Li$ we observe an $S$ -wave resonance near the threshold producing a very pronounced peak in the calculated $S$ factor of the $^6He (p,γ)$ $^7Li$ radiative-capture reaction. While we do not find a resonance near the threshold in the $p + ^6Li$ channel, in the case of $^6He + p$ reaction a resonant S -wave state is predicted at a very low energy above the reaction threshold, which could be relevant for astrophysics and its implications should be investigated. We note though that this state lies above the three-body breakup threshold not included in our method and may be influenced by three-body continuum correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure of Cl 38 and the quest for a comprehensive shell model interaction

The higher-spin structure of 38 Cl ($N=21$) was investigated following the 26 Mg( 14 C, $pn$) reaction at 30 and 37 MeV beam energies. The outgoing protons were detected in an $E-ΔE$ Si telescope placed at 0° close to the target with a Ta beam stopper between the target and telescope. Multiple $γ$ rays were detected in time coincidence with the protons using an enhanced version of the FSU $γ$-detection array. The level scheme was extended up to 8420 keV with a likely spin of 10 $ℏ$. A new multishell interaction was developed guided by the experimental information. This FSU interaction was built by fitting to the energies of 270 experimental states from 13 C to 51 Ti. Calculations using the FSU interaction reproduce observed properties of 38 Cl rather well, including the spectroscopic factors. The interaction has been successfully used to interpret the 1p1h and 2p2h configurations in some nearby nuclei as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Resonances of A = 4T = 1 isospin triplet states within the ab initio no-core Gamow shell model

Here, the A = 4 nuclei, i.e., 4 H, 4 He, and 4 Li , establish an interesting isospin T = 1 isobaric system. 4 H and 4 Li are unbound broad resonances, whereas 4 He is deeply bound in its ground state but unbound in all its excited states. The present situation is that experiments so far have not given consistent data on the resonances. Few-body calculations have well studied the scatterings of the 4N systems. In the present work, we provide many-body calculations of the broad resonance structures, in an ab initio framework with modern realistic interactions. It occurs that, indeed, 4 H, 4 Li, and excited 4 He are broad resonances, which is in accordance with experimental observations. The calculations also show that the first 1 - excited state almost degenerates with the 2 - ground state in the pair of mirror isobars of 4 H and 4 Li, which may suggest that the experimental data on energy and width are the mixture of the ground state and the first excited state. The T = 1 isospin triplet formed with an excited state of 4 He and ground states of 4 H and 4 Li is studied, focusing on the effect of isospin symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multi-shock design tools in double shell models

The essence of the project was to implement a program in python that demonstrated certain models put forth by fusion physicists. Across the time spent at Los Alamos National Laboratory, much was spent programming. Though fundamentally programmatical in nature, the task contained considerable difficulty hid within mathematically deriving and manipulating the models while maintaining data discretization. The result of solving these complications gave invaluable experience in applied mathematics. The final product, a piece of lengthy code, ended up having the full range of desired capabilities, with plans for more. Given the laser power history from a pulse shot on NIF, Omega, or any other Inertial Confinement Fusion facility, the code maps the terrain of radiation temperature inside the hohlraum throughout the period of the shot. Additionally, shock wave convergence prediction was implemented as additional capability. As the name of the project suggests, this tool acts as a compass of direction for researchers designing laser power profiles. Further development of this tool will undoubtedly assist in the research conducted on single and double-shell indirect ICF, consequently expanding the capabilities of operating stewardship over the national nuclear stockpile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of the core shell model for strengthening polymer filament interfaces

While surface segregation has been intensively studied in thin films, relatively little has been investigated in filament geometries. In contrast to thin films, the filament surface forms a continuum, with interactions both at the free air and internal filament interface, which makes it difficult to achieve an equilibrium configuration. This asymmetry has become particularly apparent with the increasing use of additive manufacturing where entire structures are constructed via layered deposition of filamentous layers from a moving thermal nozzle. In addition, this process is highly non-equilibrated with poor thermal retention and is responsible for insufficient filament–filament interfacial fusion and internal pore formation along the printing direction. These features all contribute to poor mechanical and structural properties, which are exacerbated in semi-crystalline polymers, such as Polylactic Acid (PLA), where the melting point, T m , is much higher than the glass transition, T g . Even though system is liquid the molecules become ordered within the crystal structure and this ordering interferes with interdiffusion, which in turn prevents proper interfacial fusion. Capitalizing on the differential in surface tension, we show that addition of less than 1%, or trace amounts, of longer chain polymers, an oriented core shell structure, will be formed with a self-limiting thickness of three time the radius of gyration, R g , of the tracer chains. The lower T g of the tracers and their orientation makes them poised for interdiffusion when subsequent layers are deposited producing structures with significant improvement of the mechanical properties and internal structural integrity.

36 MATERIALS SCIENCE↗