Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hartree–Fock”

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 91 records · Page 5

Efficient Monte Carlo event generation for neutrino-nucleus exclusive cross sections

Modern neutrino-nucleus cross section computations need to incorporate sophisticated nuclear models to achieve greater predictive precision. However, the computational complexity of these advanced models often limits their practicality for experimental analyses. To address this challenge, we introduce a new Monte Carlo method utilizing normalizing flows to generate surrogate cross sections that closely approximate those of the original model while significantly reducing computational overhead. As a case study, we built a Monte Carlo event generator for the neutrino-nucleus cross section model developed by the Ghent group. This model employs a Hartree-Fock procedure to establish a quantum mechanical framework in which both the bound and scattering nucleon states are solutions to the mean-field nuclear potential. The surrogate cross sections generated by our method demonstrate excellent accuracy with a relative effective sample size of more than 98.4%, providing a computationally efficient alternative to traditional Monte Carlo sampling methods for differential cross sections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of 𝐷-line energies in sodiumlike Ir

We report measurements of the 𝐷 1 and 𝐷 2 transition energies in Na-like Ir (Ir 66+ ). The 𝐷 1 3⁢𝑠−3⁢𝑝 1/2 transition energy, measured in the extreme ultraviolet (EUV) region at 169.977 ± 0.005 eV with a relative uncertainty of 29 ppm, represents one of the most precise 𝐷-line energy measurements to date for high atomic numbers 𝑍 ≥ 30. The 𝐷 2 3⁢𝑠−3⁢𝑝 3/2 transition energy, measured in the x-ray region, is 621.11 ± 0.06 eV. Therefore, there is a fine-structure separation of 451.13 ± 0.06 eV for the Na-like Ir 3⁢𝑝 levels. We present theoretical calculations of these energies using quantum electrodynamic-extended relativistic many-body perturbation theory and multiple-configuration Dirac-Hartree-Fock methods, and discuss the discrepancies between experimental and theoretical values. The prospect of determining the absolute nuclear charge radii of heavy elements, including rare isotopes, from these measurements is explored.

Atomic spectra↗

Revisiting the assignment of atomic charges in metal oxides based on core-level x-ray photoelectron spectra: The case of Ti in SrTiO3(001)

We demonstrate that assigning formal charges to transition metal (TM) cations based on core-level (CL) x-ray photoemission binding energies in oxides leads to physically inconsistent pictures of electronic structure. O 2p–TM 3d hybridization is well known to result in significant covalency in TM–O bonds, thereby reducing TM cation charges from their fully ionic values. However, the ionic bonding model remains the working paradigm for assigning TM CL features, and the resulting cation charges are often taken to be representative of the material under study. Here, we show that a more physically meaningful way to assign charges is to extract information about charge distributions utilizing Dirac–Hartree–Fock theory to calculate CL spectra from first principles and then use the resulting wave functions to determine charges based on orbital occupancies. TM cation charges can also be determined using density functional theory and Bader population analysis. We illustrate these two methods using the Ti 2p spectrum for SrTiO3(001) and show that the agreement between them is excellent. Significantly, the resulting Ti charge is considerably lower than the formal charge. The high degree of similarity between the Ti 2p spectrum for SrTiO3 and those for the rutile and anatase polymorphs of TiO2 suggests that the charge densities surrounding Ti in the latter materials are similar to that in SrTiO3. Taking a broader perspective, oxides containing other first-row transition metals also exhibit covalent character, leading to TM cation charges lower than the analogous fully ionic values in these materials as well.

Chambers, Scott A. (ORCID:000000025415043X)↗

Valley-spin polarization at zero magnetic field induced by strong hole-hole interactions in monolayer WSe 2

Monolayer transition metal dichalcogenides have emerged as prominent candidates to explore the complex interplay between spin and valley degrees of freedom. Their strong spin-orbit interaction and broken inversion symmetry lead to the spin-valley locking effect, in which carriers occupying theKandK′ valleys of the reciprocal space must have opposite spins. This effect is particularly strong for holes due to a larger spin-orbit gap in the valence band. By reducing the dimensionality of a monolayer of WSe 2 to 1D via electrostatic confinement, we demonstrate that spin-valley locking and strong hole-hole interactions lead to a ferromagnetic state where hole transport is spin-valley polarized, even without an applied magnetic field. A massive Dirac fermion model in the Hartree-Fock approximation reveals that many-body hole-exchange interactions lead to this polarized ground-state. This observation opens the possibility of implementing a robust and stable valley-polarized system, essential in valleytronic applications.

Science & Technology - Other Topics↗

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↗

Lifetimes of low-lying levels in 158 Gd

The low-lying structure of the well-deformed nucleus 158 Gd has been revisited to elucidate the nature of the low-lying states in 158 Gd. Earlier (p, t) studies identified numerous 0 + states below 4.3 MeV, prompting questions about whether these states correspond to collective vibrations or shape coexistence. New and previously reported (n, n' γ) measurements are combined, including γ-γ coincidences, excitation functions, and angular distributions, to extract lifetimes and transition probabilities for 44 excited states up to 2.7 MeV, including 32 previously unmeasured levels. Our results confirm or revise γ-ray placements and provide detailed transition strengths, revealing both weakly collective and strongly enhanced B(E2) and B(E1) transition probabilities. In particular, a tentative 0 + state at 2437.8 keV exhibits a strong interband B(E2) transition, which may be a candidate for a possible two-phonon (ββ) excitation. Systematic comparisons with neighboring Gd isotopes, Hartree–Fock–Bogoliubov, and interacting-boson model predictions suggest that the first excited 0 + state in 158 Gd is predicted to be a β-vibration, although it is weakly collective.

158Gd↗

Mojo: MLIR-based Performance-Portable HPC Science Kernels on GPUs for the Python Ecosystem

We explore the performance and portability of the novel Mojo language for scientific computing workloads on GPUs. As the first language based on the LLVM’s Multi-Level Intermediate Representation (MLIR) compiler infrastructure, Mojo aims to close performance and productivity gaps by combining Python’s interoperability and CUDA-like syntax for compile-time portable GPU programming. We target four scientific workloads: a seven-point stencil (memory-bound), BabelStream (memory-bound), miniBUDE (compute-bound), and Hartree–Fock (compute-bound with atomic operations); and compare their performance against vendor baselines on NVIDIA H100 and AMD MI300A GPUs. We show that Mojo’s performance is competitive with CUDA and HIP for memory-bound kernels, whereas gaps exist on AMD GPUs for atomic operations and for fast-math compute-bound kernels on both AMD and NVIDIA GPUs. Although the learning curve and programming requirements are still fairly low-level, Mojo can close significant gaps in the fragmented Python ecosystem in the convergence of scientific computing and AI.

Godoy, William [ORNL] (ORCID:0000000225905178)↗

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB+QRPA calculations with the Gogny interaction

The Zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98Zr. This region includes a benchmarking case and offers insights into nuclear deformation phenomena. To investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially-symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the Quasiparticle Random Phase Approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the K π = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the K π = 0 − and K π = 1 − components. Intrinsic transition densities and associ ated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e. we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Recursive algorithm for constructing antisymmetric fermionic states in first quantization mapping

We devise a deterministic quantum algorithm to produce antisymmetric states of single-particle orbitals in the first quantization mapping. Unlike sorting-based antisymmetrization algorithms, which require ordered input states and high Clifford-gate overhead, our approach initializes the state of each particle independently. For a system of $η$ particles and $N$ single-particle states, our algorithm prepares antisymmetrized states of non-trivial localized (e.g., Hartree-Fock) orbitals using $O(η^2\sqrt{N})$ $T$-gates, outperforming alternative algorithms when $η ≲ \sqrt{N}$. To achieve such scaling, we require $O(\sqrt{N})$ dirty ancilla qubits for intermediate calculations. Knowledge of the single-particle states to be antisymmetrized can be leveraged to further improve the efficiency of the circuit, and a measurement-based variant reduces gate cost by roughly a factor of two. We show example circuits for two- and three-particle systems and discuss the generalization to an arbitrary number of particles. For a specific three-particle example, we decompose the circuit into Clifford $+T$ gates and study the impact of noise on the prepared state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimate of the Total Energy of the Beh(sub 2) Molecule

The existence of gaseous BeH 2 , predicted from simple valence theory, has still not been established experimentally. We have estimated the total energy of this molecule, assumed linear symmetric, as the sum of four negative terms; the calculated SCF energy, the SCF error to the Hartree-Fock limit, the correlation energy to the exact nonrelativistic energy, and the relativistic energy.

Molecular energy↗

Polarizability of molecular hydrogen.

Perturbation theory summation method for calculating molecular hydrogen static dipole electric polarizability, using coupled Hartree- Fock approximation

Taylor, H. S.↗

Theoretical study of HSiO/+/ and HOSi/+/

Results are reported for quantum-mechanical studies of the ions HSiO(+) and HOSi(+). The equilibrium geometries, rotation constants, and relative stabilities of these ions are determined using the matrix Hartree-Fock model with a basis set of Slater exponential functions. Optimum bond lengths are calculated, and an empirical correction to the computed bond lengths is introduced. The isomer HOSi(+) is found to be more stable than HSiO(+). A frequency of 35.77 GHz is obtained for the J = 1-0 transition of HOSi(+) by assuming that the vibration-rotation interaction constants for this isomer are the same as those for HCP.

Wilson, S.↗

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗

Stochastic mean-field theory and applications to multinucleon transfer and kinetic energy dissipation processes in heavy-ion collisions

In this Review article, a brief description of the stochastic mean-field (SMF) theory for describing reaction dynamics in low-energy heavy-ion collisions at bombarding energies in the vicinity of the Coulomb barrier is presented. In these collisions, as a result of strong Pauli blocking, binary nucleon collisions do not have a significant effect on the dissipation and fluctuations. At low energies, the mean-field fluctuations, due to initial correlations, have a dominant effect on fluctuations of macroscopic variables. The SMF theory proposes the determination of an ensemble of single-particle density matrices by specifying random initial fluctuations according to a distribution law. Employing an ensemble of single-particle density matrices, not only the mean values but also the distribution functions of the one-body observables can be determined. If the di-nuclear structure is maintained in heavy-ion collisions, such as deep inelastic collisions and fast quasi-fission reactions, a much simpler description of the reaction mechanism can be derived in terms of several macroscopic variables such as mass and charge asymmetry, and relative linear and relative angular momentum. In this case, by geometric projection of the SMF equations, it is possible to derive the quantal Langevin equations for macroscopic variables. As an application of quantal transport description, an analysis of multinucleon transfers and kinetic energy dissipation and fluctuations is presented for selected quasi-fission reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗