Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONFIGURATION INTERACTION”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Low-energy enhancement in the magnetic dipole γ-ray strength functions of heavy nuclei

A low-energy enhancement (LEE), which was observed experimentally in the gamma-ray strength function (γSF) describing the decay of compound nuclei, would have profound effects on r-process nucleosynthesis if it persists in heavy neutron-rich nuclei. The LEE was shown to be a feature of the magnetic dipole (M1) strength function in configuration-interaction shell-model calculations in medium-mass nuclei. However, its existence in heavy open-shell nuclei remains an open question. Here, using a combination of many-body methods, we identify a LEE in the M1 γSFs of heavy samarium nuclei. In particular, we use the static-path plus random-phase approximation (SPA+RPA), which includes static and small-amplitude quantal fluctuations beyond the mean field. Using the SPA+RPA strength as a prior, we apply the maximum-entropy method (MEM) to obtain finite-temperature M1 γSFs from exact imaginary-time response functions calculated with the shell model Monte Carlo (SMMC) method. We find that the slope of the LEE in samarium isotopes is roughly independent of the average initial energy over a wide range below the neutron separation energy. As the neutron number increases, strength transfers to a low-energy excitation, which we interpret as the scissors mode built on top of excited states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Structure of Ne 20 states in resonance O 16 + α elastic scattering

Background: The nuclear structure of the cluster bands in 20 Ne presents a challenge for different theoretical approaches. It is especially difficult to explain the broad 0 + , 2 + states at 9 MeV excitation energy. Simultaneously, it is important to obtain more reliable experimental data for these levels in order to quantitatively assess the theoretical framework. Purpose: To obtain new data on 20 Ne α- cluster structure. Method: Thick target inverse kinematics technique was used to study the 16 O+α resonance elastic scattering and the data were analyzed using an R matrix approach. The 20 Ne spectrum, the cluster and nucleon spectroscopic factors were calculated using cluster-nucleon configuration interaction model (CNCIM). Results: We determined the parameters of the broad resonances in 20 Ne: 0 + level at 8.77 ± 0.150 MeV with a width of 750 $^{+500}_{-220}$ keV; 2 + level at 8.75 ± 0.100 MeV with the width of 695 ± 120 keV; the width of 9.48 MeV level of 65 ± 20 keV and showed that 9.19 MeV, 2 + level (if exists) should have width ≤ 10 keV. The detailed comparison of the theoretical CNCIM predictions with the experimental data on cluster states was made. Conclusions: Our experimental results by the TTIK method generally confirm the adopted data on α-cluster levels in 20 Ne. The CNCIM gives a good description of the 20 Ne positive parity states up to an excitation energy of ~7 MeV, predicting reasonably well the excitation energy of the states and their cluster and single particle properties. At higher excitations, the qualitative disagreement with the experimentally observed structure is evident, especially for broad resonances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Search for excited states in O 25

Theoretical calculations suggest the presence of low-lying excited states in 25 O. Previous experimental searches by means of proton knockout on 26 F produced no evidence for such excitations. We search for excited states in 25 O using the 24 O(d,p) 25 O reaction. The theoretical analysis of excited states in unbound 25,27 O is based on the configuration interaction approach that accounts for couplings to the scattering continuum. We use invariant-mass spectroscopy to measure neutron-unbound states in 25 O. For the theoretical approach, we use the complex-energy Gamow Shell Model and Density Matrix Renormalization Group method with a finite-range two-body interaction optimized to the bound states and resonances of 23–26 O, assuming a core of 22 O. Here we predict energies, decay widths, and asymptotic normalization coefficients. Our calculations in a large spdf space predict several low-lying excited states in 25 O of positive and negative parity, and we obtain an experimental limit on the relative cross section of a possible J π = 1/2 + state with respect to the ground state of 25 O at σ 1/2+ /σ g.s. = 0.25$^{+1.0}_{-0.25}$. We also discuss how the observation of negative parity states in 25 O could guide the search for the low-lying negative parity states in 27 O. Previous experiments based on the proton knockout of 26 F suffered from the low cross sections for the population of excited states in 25 O because of low spectroscopic factors. In this respect, neutron transfer reactions carry more promise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Nuclear deformation in the laboratory frame

Here we develop a formalism for calculating the distribution of the axial quadrupole operator in the laboratory frame within the rotationally invariant framework of the configuration-interaction shell model. The calculation is carried out using a finite-temperature auxiliary-field quantum Monte Carlo method. We apply this formalism to isotope chains of even-mass samarium and neodymium nuclei and show that the quadrupole distribution provides a model-independent signature of nuclear deformation. Two technical advances are described that greatly facilitate the calculations. The first is to exploit the rotational invariance of the underlying Hamiltonian to reduce the statistical fluctuations in the Monte Carlo calculations. The second is to determine quadruple invariants from the distribution of the axial quadrupole operator in the laboratory frame. This allows us to extract effective values of the intrinsic quadrupole shape parameters without invoking an intrinsic frame or a mean-field approximation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Few-nucleon and many-nucleon systems with semilocal coordinate-space regularized chiral nucleon-nucleon forces

We employ a variety of ab initio methods, including Faddeev-Yakubovsky equations, no-core configuration interaction approach, coupled-cluster theory, and in-medium similarity renormalization group, to perform a comprehensive analysis of the nucleon-deuteron elastic and breakup reactions and selected properties of light and medium-mass nuclei up to 48 C a using the recently constructed semilocal coordinate-space regularized chiral nucleon-nucleon potentials. We compare the results with those based on selected phenomenological and chiral EFT two-nucleon potentials, discuss the convergence pattern of the chiral expansion, and estimate the achievable theoretical accuracy at various chiral orders using an approach to quantify truncation errors of the chiral expansion without relying on cutoff variation. We also address the robustness of this method and explore alternative ways to estimate the theoretical uncertainty from the truncation of the chiral expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Statistical theory of deformation distributions in nuclear spectra

The dependence of the nuclear level density on intrinsic deformation is an important input to dynamical nuclear processes such as fission. The auxiliary-field Monte Carlo (AFMC) method is a powerful method for computing state densities. However, the statistical distribution of intrinsic shapes is not readily accessible due to the formulation of the AFMC method in a spherical configuration-interaction shell-model approach. Instead, the theory of deformation has largely relied on a mean-field approximation which breaks rotational symmetry. We show here how the distributions of the intrinsic quadrupole deformation parameters can be calculated within the AFMC method and present results for a chain of even-mass samarium nuclei ( 148 Sm, 150 Sm, 152 Sm, 154 Sm) which includes spherical, transitional, and strongly deformed isotopes. The method relies on a Landau-like expansion of the Helmholtz free energy in invariant polynomials of the quadrupole tensor. Here we find that an expansion to fourth order provides an excellent description of the AFMC results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Triaxiality in selenium-76

Background: Selenium and germanium nuclei are associated with both triaxiality and shape coexistence. The relative influence of these deformation effects on the low-lying nuclear structure remains the subject of much discussion with additional attention drawn to 76 Se and 76 Ge due to the potential for the observation of neutrinoless double- β decay. Purpose: Experimental observables related to the deformation of 76 Se are lacking in precision. The purpose of the present paper is to provide electric quadrupole matrix elements with improved precision in order to determine the deformation of low-lying states in a model-independent manner. Methods: Sub-barrier Coulomb excitation was employed at the reaccelerated beam facility of the National Superconducting Cyclotron Laboratory using the JANUS setup. Using this method nineteen E2 matrix elements were extracted. Results: Extracted matrix elements agree within uncertainties with those in the literature but with improved precision. Through both a comparison with geometric models and a model-independent evaluation of E2 matrix elements using rotational invariants, the ground state of 76 Se is best described as having a significant triaxial component, but is not maximally triaxially deformed. Conclusions: Selenium-76 exhibits a significant degree of triaxiality in its ground state. Finally, a detailed comparison with configuration-interaction calculations indicates that this can be well reproduced theoretically.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Finite-temperature many-body perturbation theory in the canonical ensemble

Benchmark data are presented for the zeroth- through third-order many-body perturbation corrections to the electronic Helmholtz energy, internal energy, and entropy in the canonical ensemble in a wide range of temperature. They are determined as numerical λ-derivatives of the respective quantities computed by thermal full configuration interaction with a perturbation-scaled Hamiltonian, $\hat{H} = \hat{H}_0 + λ\hat{V}$. Sum-over-states analytical formulas for up to the third-order corrections to these properties are also derived as analytical λ-derivatives. These formulas, which are verified by exact numerical agreement with the benchmark data, are given in terms of the Hirschfelder–Certain degenerate perturbation energies and should be valid for both degenerate and nondegenerate reference states at any temperature down to zero. Further, the results in the canonical ensemble are compared with the same in the grand canonical ensemble.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Emergent Sp(3,$\mathbb{R}$) Dynamical Symmetry in the Nuclear Many-Body System from an Ab Initio Description

Ab initio nuclear theory provides not only a microscopic framework for quantitative description of the nuclear many-body system, but also a foundation for deeper understanding of emergent collective correlations. In this work, a symplectic Sp(3, $\mathbb{R}$) $\supset$ U(3) dynamical symmetry is identified in ab initio predictions, from a no-core configuration interaction approach, and found to provide a qualitative understanding of the spectrum of 7 Be. Low-lying states form an Elliott SU(3) spectrum, while an Sp(3, $\mathbb{R}$) excitation gives rise to an excited rotational band with strong quadrupole connections to the ground state band.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

The QCAD framework for quantum device modeling

We present the Quantum Computer Aided Design (QCAD) simulator that targets modeling quantum devices, particularly Si double quantum dots (DQDs) developed for quantum computing. The simulator core includes Poisson, Schrodinger, and Configuration Interaction solvers which can be run individually or combined self-consistently. The simulator is built upon Sandia-developed Trilinos and Albany components, and is interfaced with the Dakota optimization tool. It is being developed for seamless integration, high flexibility and throughput, and is intended to be open source. The QCAD tool has been used to simulate a large number of fabricated silicon DQDs and has provided fast feedback for design comparison and optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

NeuroCoreX: An Open-Source FPGA-Based Spiking Neural Network Emulator with On-Chip Learning

Spiking Neural Networks (SNNs) are computational models inspired by the event-driven communication and connectivity patterns of biological neural circuits. They enable high energy efficiency and natural support for diverse architectures ranging from layered networks to small-world and graphstructured topologies. In this work, we introduce NeuroCoreX, an open-source, FPGA-based spiking neural network emulator that provides real-time, on-chip learning and flexible network organization. NeuroCoreX supports both feedforward sensory inputs streamed directly from sensors or PCs via UART and recurrent on-chip connectivity, enabling simultaneous processing and learning from external stimuli and internal network dynamics-capabilities rarely available in existing FPGA SNN platforms. The system implements a Leaky Integrate-and-Fire (LIF) neuron model with current-based synapses and supports pair-based STDP learning on both feedforward and recurrent synapses. A lightweight Python interface enables interactive configuration, live monitoring, weight read-back, and experiment control. Importantly, NeuroCoreX is tightly integrated with the SuperNeuroMAT simulator, allowing SNN models to be transferred seamlessly from software to hardware for hardware-in-the-loop development. By combining real-time plasticity, flexible connectivity, and an open-source VHDL implementation, NeuroCoreX provides an extensible and accessible platform for neuromorphic research, algorithm-hardware co-design, and energy-efficient edge intelligence.

Gautam, Ashish [ORNL]↗
DOE OSTI2021

A Visual Designer of Layer‐wise Relevance Propagation Models

Abstract Layer‐wise Relevance Propagation (LRP) is an emerging and widely‐used method for interpreting the prediction results of convolutional neural networks (CNN). LRP developers often select and employ different relevance backpropagation rules and parameters, to compute relevance scores on input images. However, there exists no obvious solution to define a “best” LRP model. A satisfied model is highly reliant on pertinent images and designers' goals. We develop a visual model designer, named as VisLRPDesigner, to overcome the challenges in the design and use of LRP models. Various LRP rules are unified into an integrated framework with an intuitive workflow of parameter setup. VisLRPDesigner thus allows users to interactively configure and compare LRP models. It also facilitates relevance‐based visual analysis with two important functions: relevance‐based pixel flipping and neuron ablation. Several use cases illustrate the benefits of VisLRPDesigner. The usability and limitation of the visual designer is evaluated by LRP users.

Huang, Xinyi↗
DOE OSTI2025

Learning local and semi-local density functionals from exact exchange-correlation potentials and energies

Finding accurate exchange-correlation (XC) functionals remains the defining challenge in density functional theory (DFT). Despite 40 years of active development, attaining general purpose chemical accuracy is still elusive with existing functionals. We present a data-driven pathway to learn the XC functional by using the exact density, XC energy, and XC potential. While the exact densities are obtained from accurate configuration interaction (CI), the exact XC energies and XC potentials are obtained via inverse DFT calculations on the CI densities. We demonstrate how simple neural network (NN)–based local density approximation (LDA) and generalized gradient approximation (GGA), trained on just five atoms and two molecules, provide remarkable improvement in total energies and densities. Particularly, the NN-based GGA functional attains similar accuracy as the higher rung SCAN meta-GGA on various thermochemistry datasets. These results underscore the promise of using the XC potential in modeling XC functionals and can pave the way for systematic learning of increasingly accurate XC functionals.

Science & Technology - Other Topics↗
DOE OSTI2022

Multiparticle-hole excitations in nuclei near N = Z = 20: $^{41}$K

This experimental study of high-spin structure near the N = Z = 20 region focuses on 41 K and also reports three newly observed γ transitions in 41 Ca from the same reaction. High-spin states were populated using the 26 Mg( 18 O, p2nγ) 41 K and 26 Mg( 18 O, 3nγ) 41 Ca reactions at a beam energy of 50 MeV at the Florida State University John D. Fox Superconducting Linear Accelerator Laboratory, employing the FSU high-purity germanium detector array. The level scheme of 41 K was extended to 12325 keV, possibly up to J π = 25/2 − or 27/2 + , by means of 25 new γ-ray transitions, and that of 41 Ca to 9916 keV. Linear polarization and angular-distribution measurements were used to provide spin and parity information for several states in the 41 K decay scheme. The results are compared with spsdpf cross-shell shell-model calculations using the FSU interaction. Configurations involving zero or one nucleon promoted from the sd to the fp shell reproduce the energies of many known states reasonably well, while multi-particle excitations reveal a more complex interplay of single-particle motion, collective pairing, and deformation, posing an interesting challenge for future theoretical work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

NuLattice: Ab initio computations of atomic nuclei on lattices

Here, we introduce NuLattice, a Python software package for ab initio computations of atomic nuclei on lattices. The computational tools consist of Hartree Fock, the coupled cluster method, the in-medium similarity renormalization group, and full configuration interaction. At present, the employed interactions are from pion-less effective field theory at leading order and consist of two-body and three-body contacts. We present results for light nuclei 2 H, 3,4 He, 8 Be, 12 C, and 16 O. NuLattice algorithms exploit the sparsity and locality of lattice interactions, and as a result computations can be run on laptops.

Rothman, Maxwell [Univ. of Tennessee, Knoxville, T↗
DOE OSTI2020

CSBCOO SPMM: Compressed Sparse Block Coordinate format sparse matrix dense vector multiply benchmark (CSBCOO SPMM) v1.0

The purpose of this microbenchmark is to provide a means to explore different programming methodologies using a simple, but not trivial, mathematical kernel. The kernel is based on the matrix vector product in the MFDn nuclear configuration interaction code's LOBPCG eigensolver. This particular mathematical operation is a key element in the solution of systems of equations and block eigensolvers.

Cook, Brandon↗
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↗