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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Microsecond Molecular Dynamics Simulations of Proteins Using a Quasi-Equilibrium Solvation Shell Model

In this work, we describe the development and implementation of a quasi-equilibrium hydration shell model of biomolecular solvation with adaptive boundaries. Applying the model to microsecond-long molecular dynamics simulations of several protein systems of varying complexity, we find that the model simulation results are of comparable quality to those obtained from simulations of fully solvated systems, but at a reduced computational cost. We discuss the dominant sources of error in the model and outline directions for future improvements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

New insights into backbending in the symmetry-adapted shell-model framework

Here we provide insights into the backbending phenomenon within the symmetry-adapted framework which naturally describes the intrinsic deformation of atomic nuclei. For 20 Ne, the canonical example of backbending in light nuclei, the ab initio symmetry-adapted no-core shell model shows that while the energy spectrum replicates the backbending from experimental energies under the rotor-model assumption, there is no change in the intrinsic deformation or intrinsic spin of the yrast band around the backbend. For the traditional example of 48 Cr, computed in the valence shell with empirical interactions, we confirm a high-spin nucleus that is effectively near-spherical, in agreement with previous models. However, we find that this spherical distribution results, on average, from an almost equal mixing of deformed prolate shapes with deformed oblate shapes. Microscopic calculations confirm the importance of spin alignment and configuration mixing, but surprisingly unveil no anomalous increase in moment of inertia. Furthermore, this finding opens the path toward further understanding the rotational behavior and moment of inertia of medium-mass nuclei.

20 ≤ A ≤ 38↗
DOE OSTI2021

Shell model analysis of the B ( E 2,2 + → 0 + ) values in the A = 70, T = 1 triplet 70 Kr, 70 Br, and 70 Se

The B(E2, 2 + → 0 + ) transition strengths of the T=1 isobaric triplet 70 Kr, 70 Br, 70 Se, recently measured at the RIKEN Radioactive Isotope Beam Factory (RIBF), are discussed in terms of state-of-the-art large scale shell model calculations using the JUN45 and JUN45 + LNPS plus Coulomb interactions. In this Letter we argue that, depending on the effective charges used, the calculations are either in line with the experimental data within statistical uncertainties, or the anomaly happens in 70 Br, rather than 70 Kr. In the latter case, we suggest that it can be due to the presence of a hitherto undetected 1 + T = 0 state below the yrast 2 + T = 1 state. Furthermore, our results do not support a shape change of 70 Kr with respect to the other members of the isobaric multiplet.

59 ≤ A ≤ 89↗
DOE OSTI2022

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↗
DOE OSTI2020

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↗
DOE OSTI2020

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↗
DOE OSTI2020

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↗
NASA NTRS1985

Onion-shell model for cosmic ray electrons and radio synchrotron emission in supernova remnants

The spectrum of cosmic ray electrons, accelerated in the shock front of a supernova remnant (SNR), is calculated in the test-particle approximation using an onion-shell model. Particle diffusion within the evolving remnant is explicity taken into account. The particle spectrum becomes steeper with increasing radius as well as SNR age. Simple models of the magnetic field distribution allow a prediction of the intensity and spectrum of radio synchrotron emission and their radial variation. The agreement with existing observations is satisfactory in several SNR's but fails in other cases. Radiative cooling may be an important effect, especially in SNR's exploding in a dense interstellar medium.

Beck, R.↗
DOE OSTI2022

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↗
DOE OSTI2024

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↗
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 OSTI2024

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↗
DOE OSTI2020

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↗
NASA NTRS1977

The moon's permanent magnetic field - A cratered-shell model

The possibility is considered that the properties of the larger-scale lunar magnetic field may depend in a rather straightforward manner on the geometrical properties of the distribution of larger craters. The relevant effects of cratering are analyzed in the framework of a model of a spherically symmetric shell magnetized to some specific intensity by a concentric dipolar magnetic field. The remanent dipole moment of the model is determined by calculating the dipole moment of the material removed to form craters and them computing the remanent dipole moment for three mutually orthogonal orientations of the magnetizing dipole. The resulting absolute values for the remanent dipole moment are found to be comparable to the upper limit of 10 to the -19th power G-cu cm estimated from Apollo subsatellite measurements. It is concluded that for a crust with an empirically reasonable mean magnetization, the craters can account for a dipole moment as great as that indicated by the estimated upper limit.

Weiss, H.↗
NASA NTRS1988

Hierarchic plate and shell models based on p-extension

This paper is concerned with formulations of finite element models for beams, arches, plates and shells based on the principle of virtual work. The focus is on computer implementation of hierarchic sequences of finite element models suitable for numerical solution of a large variety of practical problems which may concurrently contain thin and thick plates and shells, stiffeners, and regions where three dimensional representation is required. The approximate solutions corresponding to the hierarchic sequence of models converge to the exact solution of the fully three dimensional model. The stopping criterion is based on (1) estimation of the relative error in energy norm; (2) equilibrium tests; and (3) observation of the convergence of quantities of interest.

Szabo, Barna A.↗
NASA NTRS1989

Hierarchic plate and shell models based on p-extension

Formulation of hierarchic sequences of finite element models for beams, arches, plates and shells based on the principle of virtual work is described. The exact solutions corresponding to models in the hierarchic sequence converge to the exact solution of the fully three-dimensional problem of linear elasticity. The stopping criterion is that the functionals of interest must be substantially independent of the choice of the model. This process is closely related to p-extensions. Aspects of implementation are discussed in connection with axisymmetric shells and an example is presented. An application of superconvergent extraction methods for the computation of stress resultants is demonstrated.

Szabo, B. A.↗
NASA NTRS1988

Hierarchic plate and shell models based on p-extension

Formulations of finite element models for beams, arches, plates and shells based on the principle of virtual work was studied. The focus is on computer implementation of hierarchic sequences of finite element models suitable for numerical solution of a large variety of practical problems which may concurrently contain thin and thick plates and shells, stiffeners, and regions where three dimensional representation is required. The approximate solutions corresponding to the hierarchic sequence of models converge to the exact solution of the fully three dimensional model. The stopping criterion is based on: (1) estimation of the relative error in energy norm; (2) equilibrium tests, and (3) observation of the convergence of quantities of interest.

Szabo, Barna A.↗
NASA NTRS1997

Variability in shell models of GRBs

Many cosmological models of gamma-ray bursts (GRBs) assume that a single relativistic shell carries kinetic energy away from the source and later converts it into gamma rays, perhaps by interactions with the interstellar medium or by internal shocks within the shell. Although such models are able to reproduce general trends in GRB time histories, it is difficult to reproduce the high degree of variability often seen in GRBs. The authors investigate methods of achieving this variability using a simplified external shock model. Since the model emphasizes geometric and statistical considerations, rather than the detailed physics of the shell, it is applicable to any theory that relies on relativistic shells. They find that the variability in GRBs gives strong clues to the efficiency with which the shell converts its kinetic energy into gamma rays.

Sumner, M. C.↗