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 253 records · Page 14

Quantum Monte Carlo calculations in configuration space with three-nucleon forces

Neutron matter, through its connection to neutron stars as well as systems like cold atom gases, is one of the most interesting yet computationally accessible systems in nuclear physics. The Configuration-Interaction Monte Carlo (CIMC) method is a stochastic many-body technique allowing to tackle strongly coupled systems. In contrast to other Quantum Monte Carlo methods employed in nuclear physics, the CIMC method can be formulated directly in momentum space allowing for an efficient use of non-local interactions. In this work we extend CIMC method to include three-nucleon interactions through the normal-ordered two-body approximation. Here, we present results for the equation of state of neutron matter in line with other many-body calculations that employ low resolution chiral interactions, and provide predictions for the momentum distribution and the static structure factor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled-cluster theory for strong entanglement in nuclei

Atomic nuclei can exhibit shape coexistence and multireference physics that enters in their ground states, and accurately capturing the ensuing correlations and entanglement is challenging. Here, we address this problem by applying single-reference coupled-cluster theory based on spherical and deformed reference states and the tailored coupled-cluster method. The latter combines configuration interactions to capture static correlations with coupled-cluster theory for dynamic correlations. We compute the atomic nuclei 12 C, 28 Si, and 56 Ni and find that the tailored coupled-cluster method and the single-reference approach based on a deformed Hartree-Fock state yield the most accurate results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

𝛽 decay of 36 Mg and 36 Al: Identification of a 𝛽-decaying isomer in 36 Al

The level structure of 36 Al has been studied via 𝛽 decay of 36 Mg at the Facility for Rare Isotope Beams (FRIB) and the National Superconducting Cyclotron Laboratory (NSCL). A long-lived isomer in 36 Al was identified which decays by 𝛽 to an excited state of 36 Si. The ground state and the isomeric state of 36 Al were found to populate different energy levels of 36 Si. Furthermore, the results from the two data sets in the present work complement each other. Configuration interaction calculations performed with the FSU shell-model Hamiltonians provide reasonable descriptions to the experimental observations and offer insight into future improvements of the theoretical interpretation.

20 ≤ A ≤ 38↗

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↗

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↗

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↗

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↗

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

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↗

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↗

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↗

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↗

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↗

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↗

Mapping Rydberg States of H 2 with the Halfium R-Matrix Method

In this article, we use the Halfium R-matrix method to investigate the Rydberg states of the H 2 molecule up to n = 20, filling the gap above the low-lying bound states already calculated with configuration interaction packages. Moreover, we show that the use of Quantum Defect Theory scaling laws, allows for a comprehensive analysis of the regular patterns resulting from the coupling between Rydberg series and doubly excited states. The results should open the door for more efficient quasi-diabatization of the potential energy curves which is required for calculating cross sections and rate coefficients of the (e + H 2 + ) collisional processes, involved in the plasma modeling for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Multireference Approach to Electron and Electron–Nuclear Dynamics in Nanomaterials (Final Report)

Many important chemical and physical phenomena involve dynamics on large number of electronic states. Thus, there is a critical need to develop methods to simulate dynamics in dense manifolds of states. Towards this end, we have: a) developed the multiple cloning in dense manifolds of states (MCDMS) method, which is capable of accurately modeling the quantum mechanical coherence between populations on a large number of electronic states, b) implemented MCDMS into the free, open-source PySpawn software package, c) developed graphics processing unit-accelerated algorithms modeling electron dynamics in light fields via Floquet time-dependent configuration interaction (F-TDCI), and d) critically compared different orbital bases in order to achieve an accurate and efficient F-TDCI expansion. This grant ended in August 2020, when our group moved to from Michigan State University to Stony Brook University, where this project continues under grant number DE-SC0021643.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Systematic, Polynomial-Cost Approach to Exact Correlation Energies (Final Technical Report)

The full configuration interaction (FCI) wave function provides the exact solution to the Schrödinger equation in a given basis set. While FCI is intractable due to exponential computational costs, the many-body expansion (or the method of increments) can reduce scaling to a low-order polynomial with system size. This project entailed advances in the incremental FCI (iFCI) approach, designed to allow iFCI to reach larger system sizes and maintain its intrinsic high accuracy. This document describes advances in solvers for iFCI, strategies to treat multiple charge and spin states, and virtual state management methods to reduce memory requirements. Overall, this project allows iFCI to correlate (for the first time) 142 valence electrons in 444 orbitals in a realistic model of a transition metal complex.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗