Search NASA⌕ Search

SEARCH · Search NASA

Results for “fermi operator expansion”

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.

9 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective

Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

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

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

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

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multipolar Fermi Surface Deformations in Sr 2 ⁢RuO 4 Probed by Resistivity and Sound Attenuation: A Window into Electron Viscosity and the Collision Operator

Recent developments in electron hydrodynamics have demonstrated the importance of considering the full structure of the electron-electron scattering operator, which encodes a sequence of lifetimes, one for each component of the Fermi surface deformation in a multipolar expansion. In this context, the dipolar lifetime is measured by resistivity, whereas the quadrupolar component probes the viscosity and can be measured in the bulk via sound attenuation. We introduce a framework to extract the collision operator of an arbitrary metal by combining resistivity and sound attenuation measurements with a realistic calculation of the scattering operator that includes multiband and umklapp effects. The collision operator allows for the prediction of a plethora of properties, including the nonlocal conductivity, and can be used to predict hydrodynamic behavior for bulk metals. As a first application, we apply this framework to Sr 2 ⁢RuO 4 in a temperature range where electron-electron scattering is dominant. Furthermore, we find quantitative agreement between our model and the temperature dependence of both the resistivity and the sound attenuation, we find the quadrupolar (𝐵 1⁢𝑔 ) relaxation rate to be 30% higher than the dipolar one due to the presence of hot spots on the 𝛾 band, and we predict a strongly anisotropic viscosity arising from the 𝛼 and 𝛽 bands.

Boltzmann theory↗
NASA NTRS2001

Polarizabilities and Other Properties of the td Muons Molecular Ion

Wavefunctions of Hylleraas type were used earlier to calculate energy levels of muonic systems. Recently, we found in the case of the molecular ions H2+, D2+ and HD+ that it was necessary to include high powers of the internuclear distance in the Hylleraas functions to localize the nuclear motion when treating the ions as three-body systems without invoking the Born-Oppenheimer approximation. We try the same approach in a muonic system, td(mu-). Improved convergence is obtained for J = 0 and 1 states for shorter expansions when we use this type of generalized Hylleraas function, but as the expansion length increases the high powers are no longer useful. We obtain good energy values for the two lowest J = 0 and J = 1 states and compare them with the best earlier calculations. Expectation values are obtained for various operators, the Fermi contact parameters, and the permanent quadrupole moment. The cusp conditions are also calculated. The polarizability of the ground state is then calculated using second-order perturbation theory with intermediate J = 1 pseudostates. It should be possible to measure the polarizability by observing Rydberg states of atoms with td(mu-) acting as the nucleus.

Bhatia, A. K.↗
NASA NTRS2009

Concluding Remarks: The Current Status and Future Prospects for GRB Astronomy

We are in a remarkable period of discovery in GRB astronomy. The current satellites including Swift, Fermi. AGILE and INTEGRAL are detecting and observing bursts of all varieties. Increasing capabilities for follow-up observations on the ground and in space are leading to rapid and deep coverage across the electromagnetic spectrum, The future will see continued operation of the current experiments and with future missions like SVOM plus possible rni_Ssions like JANUS and EXIST. An exciting expansion of capabilities is occurring in areas of gravitational waves and neutrinos that could open new windows on the GRB phenomenon. Increased IR capabilities on the ground and with missions like JWST will enable further exploration of high redshift bursts. The future is bright.

Gehrels, Neil↗
DOE OSTI2025

Energy-enhanced expansion of the standard model effective field theory

We formalize energy-scaling arguments in the standard model effective field theory (SMEFT) to estimate the effects of operators up to dimension ten. Our approach relies on weakly coupled UV completions with no presumed large hierarchies between the Wilson coefficients. We introduce a classification based on the number of external legs and an energy-counting parameter. We establish a dual expansion in 𝑣/Λ and 𝐸/Λ. Extending to four-, five-, and six-particle vertices, our framework highlights energy-enhanced operators that dominate high-energy processes at the High Luminosity-Large Hadron Collider. This organization streamlines experimental analyses to only include operators with energetic impact in their analyses and enhances the discoverability of new physics within the SMEFT framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Parton physics from a heavy-quark operator product expansion: Lattice QCD calculation of the fourth moment of the pion distribution amplitude

The pion light-cone distribution amplitude (LCDA) is an essential nonperturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of 𝑚 𝜋 = 550 MeV. This calculation finds ⟨𝜉 2 ⟩ = 0.202⁢(8)⁢(9) and ⟨𝜉 4 ⟩ = 0.039⁢(28)⁢(11) where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the $\overline{\textrm{MS}}$ scheme at a renormalization scale of 𝜇 = 2 GeV.

Detmold, William [Massachusetts Inst. of Technolog↗