Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function approximation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Half-Metallic Ferromagnetism Character in Cr-Doped CaS Diluted Magnetic Insulator and Semiconductor: an Ab Initio Study

The overall aim of this study is to investigate theoretically the structural, electronic, and magnetic properties of calcium sulfide (CaS) doped with chromium (Cr) impurity, in order to conduct a new search dilute magnetic semiconductors (DMS) suitable for different applications in electronics and spintronics. For measuring, the physical property of this compound is implemented using the first principles approach employed in WIEN2K code. The structural characteristics are optimized using the Generalized Gradient Approximation established by Perdew-Burk-Ernzerhof (PBE-GGA). We calculate and minimize the total energy of the three ternary compounds (Ca{sub 0.75}Cr{sub 0.25}S, Ca{sub 0.50}Cr{sub 0.50}S, and Ca{sub 0.25}Cr{sub 0.75}S) in the paramagnetic (PM), ferromagnetic (FM), and antiferromagnetic (AFM) phase. We find all compounds stable in (FM) structure, whereas the modified Becke and Johnson local density approximation (mBJ-LDA) functional has been employed to evaluate the electronic and magnetic properties. Based on our findings, indicate that this system revealed a half-metallic ferromagnetic behavior with half-metallic gap (HM) and 100% spin-polarized at the fermi level for all chromium (Cr) concentrations. This advantageous set of properties is due to the half-metallic behavior, where the majority spin and minority spin exhibit metallic and semiconducting behaviors respectively. The chromium atom is the most important source of the total magnetic moment in these compounds (4 μβ) by comparison with magnetic moments produced by Ca and S atoms, which have minor contribution. Finally, our prediction results require an experimental confirmation in the future.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Decision-Making Machine Learning Approach in Hermite Spectral Approximations of Partial Differential Equations

The accuracy and effectiveness of Hermite spectral methods for the numerical discretization of partial differential equations on unbounded domains are strongly affected by the amplitude of the Gaussian weight function employed to describe the approximation space. This is particularly true if the problem is under-resolved, i.e., there are no enough degrees of freedom. The issue becomes even more crucial when the equation under study is time-dependent, forcing in this way the choice of Hermite functions where the corresponding weight depends on time. In order to adapt dynamically the approximation space, it is here proposed an automatic decision-making process that relies on machine learning techniques, such as deep neural networks and support vector machines. The algorithm is numerically tested with success on a simple 1D problem, but the main goal is its exportability in the context of more serious applications. Here we also show at the end an application in the framework of plasma physics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Empirical Study of the Effect of Thermal Loading on the Heating Efficiency of Variable-Speed Air Source Heat Pumps

Heating buildings with air source heat pumps (ASHPs) has the potential to save energy compared to utilizing conventional heat sources. Accurate understanding of the efficiency of ASHPs is important to maximize the energy savings. While it is well understood that, in general, ASHP efficiency decreases with decreasing outdoor temperature, it is not well understood how the ASHP efficiency changes with different levels of thermal loading, even though it is an important consideration for sizing and controlling ASHPs. The goal of this study was to create an empirical model of the ASHP efficiency as a function of two independent variables–outside temperature and level of thermal loading. Four ductless mini-split ASHPs were evaluated in a cold chamber where the temperature (representing the outdoor temperature) was varied over a wide range. For each temperature, the ASHP performance data were collected at several levels of thermal loading. The data for all four ASHPs were combined and approximated with an analytical function that can be used as a general model for the ASHP steady-state efficiency as a function of the outside temperature and level of thermal loading. To the knowledge of the authors, no such empirical model that is solely based on third-party test data has been published before. While limitations exist, the model can be used to help guide future selection and operation of ASHPs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Existence of La-site antisite defects in LaMO 3 (M = Mn, Fe, and Co) predicted with many-body diffusion quantum Monte Carlo

The properties of LaMO 3 (M: 3d transition metal) perovskite crystals are significantly dependent on point defects, whether introduced accidentally or intentionally. The most studied defects in La-based perovskites are the oxygen vacancies and doping impurities on the La and M sites. Here, we identify that intrinsic antisite defects, the replacement of La by the transition metal, M, can be formed under M-rich and O-poor growth conditions, based on results of an accurate many-body ab initio approach. Our fixed-node diffusion Monte Carlo (FNDMC) calculations of LaMO 3 (M = Mn, Fe, and Co) find that such antisite defects can have low formation energies and are magnetized. Complementary density functional theory (DFT)-based calculations show that Mn antisite defects in LaMnO 3 may cause the p-type electronic conductivity. These features could affect spintronics, redox catalysis, and other broad applications. Our bulk validation studies establish that FNDMC reproduces the antiferromagnetic state of LaMnO 3 , whereas DFT with PBE (Perdew–Burke–Ernzerhof), SCAN (strongly constrained and appropriately normed), and the LDA+U (local density approximation with Coulomb U) functionals all favor ferromagnetic states, at variance with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Weak Form Scientific Machine Learning: Test Function Construction for System Identification

Weak form Scientific Machine Learning (WSciML) is a recently developed framework for data-driven modeling and scientific discovery. It leverages the weak form of equation error residuals to provide enhanced noise robustness in system identification via convolving model equations with test functions, reformulating the problem to avoid direct differentiation of data. The performance, however, relies on wisely choosing a set of compactly supported test functions. In this work, we mathematically motivate a novel data-driven method for constructing Single-scale-Local reference functions for creating the set of test functions. Our approach numerically approximates the integration error introduced by the quadrature and identifies the support size for which the error is minimal, without requiring access to the model parameter values. Through numerical experiments across various models, noise levels, and temporal resolutions, we demonstrate that the selected supports consistently align with regions of minimal parameter estimation error. We also compare the proposed method against the strategy for constructing Multi-scale-Global (and orthogonal) test functions introduced in our prior work, demonstrating the improved computational efficiency.

FOS: Computer and information sciences↗
DOE OSTI2026

Augmenting subspace optimization methods with linear bandits

In this work, we consider the framework of methods for unconstrained minimization that are, in each iteration, restricted to a model that is only a valid approximation to the objective function on some affine subspace containing an incumbent point. These methods are of practical interest in computational settings where derivative information is either expensive or impossible to obtain. Recent attention has been paid in the literature to employing randomized matrix sketching for generating the affine subspaces within this framework. We consider a relatively straightforward, deterministic augmentation of such a generic subspace optimization method. In particular, we consider a sequential optimization framework where actions consist of one-dimensional linear subspaces and rewards consist of (approximations to) the magnitudes of directional derivatives computed in the direction of the action subspace. Reward maximization in this context is consistent with maximizing lower bounds on descent guaranteed by first-order Taylor models. This sequential optimization problem can be analysed through the lens of dynamic regret. We modify an existing linear upper confidence bound (UCB) bandit method and prove sublinear dynamic regret in the subspace optimization setting. We demonstrate the efficacy of employing this linear UCB method in a setting where forward-mode algorithmic differentiation can provide directional derivatives in arbitrary directions and in a derivative-free setting. For the derivative-free setting, we propose SS-POUNDers, an extension of the derivative-free optimization method POUNDers that employs the linear UCB mechanism to identify promising subspaces. Our numerical experiments suggest a preference, in either computational setting, for employing a linear UCB mechanism within a subspace optimization method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

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

Production of light-flavor hadrons in pp collisions at $\sqrt{s}~=~7\text { and }\sqrt{s} = 13 \, \text { TeV} $

The production of \(\pi ^{\pm }\) , \(\mathrm{K}^{\pm }\) , \(\mathrm{K}^{0}_{S}\) , \(\mathrm{K}^{*}(892)^{0}\) , \(\mathrm{p}\) , \(\phi (1020)\) , \(\Lambda \) , \(\Xi ^{-}\) , \(\Omega ^{-}\) , and their antiparticles was measured in inelastic proton–proton (pp) collisions at a center-of-mass energy of \(\sqrt{s}\) = 13 TeV at midrapidity ( \(|y|<0.5\) ) as a function of transverse momentum ( \(p_{\mathrm{T}}\) ) using the ALICE detector at the CERN LHC. Furthermore, the single-particle \(p_{\mathrm{T}}\) distributions of \(\mathrm{K}^{0}_{S}\) , \(\Lambda \) , and \(\overline{\Lambda }\) in inelastic pp collisions at \(\sqrt{s} = 7\) TeV are reported here for the first time. The \(p_{\mathrm{T}}\) distributions are studied at midrapidity within the transverse momentum range \(0\le p_{\mathrm{T}}\le 20\) GeV/ c , depending on the particle species. The \(p_{\mathrm{T}}\) spectra, integrated yields, and particle yield ratios are discussed as a function of collision energy and compared with measurements at lower \(\sqrt{s}\) and with results from various general-purpose QCD-inspired Monte Carlo models. A hardening of the spectra at high \(p_{\mathrm{T}}\) with increasing collision energy is observed, which is similar for all particle species under study. The transverse mass and \(x_{\mathrm{T}}\equiv 2p_{\mathrm{T}}/\sqrt{s}\) scaling properties of hadron production are also studied. As the collision energy increases from \(\sqrt{s}\) = 7–13 TeV, the yields of non- and single-strange hadrons normalized to the pion yields remain approximately constant as a function of \(\sqrt{s}\) , while ratios for multi-strange hadrons indicate enhancements. The \(p_\mathrm{{T}}\) -differential cross sections of \(\pi ^{\pm }\) , \(\mathrm {K}^{\pm }\) and \(\mathrm {p}\) ( \(\overline{\mathrm{p}}\) ) are compared with next-to-leading order perturbative QCD calculations, which are found to overestimate the cross sections for \(\pi ^{\pm }\) and \(\mathrm{p}\) ( \(\overline{\mathrm{p}}\) ) at high \(p_\mathrm{{T}}\) .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Measurements of dihadron correlations relative to the event plane in Au+Au collisions at $\sqrt{s_{NN}} =200$ GeV*

Dihadron azimuthal correlations containing a high transverse momentum ($ p_{T} $) trigger particle are sensitive to the properties of the nuclear medium created at RHIC through the strong interactions occurring between the traversing parton and the medium, i.e. jet-quenching. Previous measurements revealed a strong modification to dihadron azimuthal correlations in Au+Au collisions with respect to p+p and d+Au collisions. The modification increases with the collision centrality, suggesting a path-length or energy density dependence to the jet-quenching effect. This paper reports STAR measurements of dihadron azimuthal correlations in mid-central (20%-60%) Au+Au collisions at $ \sqrt{s_{\rm{NN}}} = 200 $ GeV as a function of the trigger particle's azimuthal angle relative to the event plane, $ \phi_{s} = | \phi_{t}- \psi_{{\rm{EP}}}| $. The azimuthal correlation is studied as a function of both the trigger and associated particle $ p_{T} $. The subtractions of the combinatorial background and anisotropic flow, assuming Zero Yield At Minimum (ZYAM), are described. The correlation results are first discussed with subtraction of the even harmonic (elliptic and quadrangular) flow backgrounds. The away-side correlation is strongly modified, and the modification varies with $ \phi_{s} $, with a double-peak structure for out-of-plane trigger particles. The near-side ridge (long range pseudo-rapidity $ \Delta\eta $ correlation) appears to drop with increasing $ \phi_{s} $ while the jet-like component remains approximately constant. The correlation functions are further studied with the subtraction of odd harmonic triangular flow background arising from fluctuations. It is found that the triangular flow, while responsible for the majority of the amplitudes, is not sufficient to explain the $ \phi_{s} $-dependence of the ridge or the away-side double-peak structure. The dropping ridge with $ \phi_{s} $ could be attributed to a $ \phi_{s} $-dependent elliptic anisotropy; however, the physics mechanism of the ridge remains an open question. Even with a $ \phi_{s} $-dependent elliptic flow, the away-side correlation structure is robust. These results, with extensive systematic studies of the dihadron correlations as a function of $ \phi_{s} $, trigger and associated particle $ p_{T} $, and the pseudo-rapidity range $ \Delta\eta $, should provide stringent inputs to help understand the underlying physics mechanisms of jet-medium interactions in high energy nuclear collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Linearization errors in discrete goal-oriented error estimation

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.

42 ENGINEERING↗
DOE OSTI2023

Probing reaction channels via reinforcement learning

Abstract Chemical reactions are dynamical processes involving the correlated reorganization of atomic configurations, driving the conversion of an initial reactant into a result product. By virtue of the metastability of both the reactants and products, chemical reactions are rare events, proceeding fleetingly. Reaction pathways can be modelled probabilistically by using the notion of reactive density in the phase space of the molecular system. Such density is related to a function known as the committor function, which describes the likelihood of a configuration evolving to one of the nearby metastable regions. In theory, the committor function can be obtained by solving the backward Kolmogorov equation (BKE), which is a partial differential equation (PDE) defined in the full dimensional phase space. However, using traditional methods to solve this problem is not practical for high dimensional systems. In this work, we propose a reinforcement learning based method to identify important configurations that connect reactant and product states along chemical reaction paths. By shooting multiple trajectories from these configurations, we can generate an ensemble of states that concentrate on the transition path ensemble. This configuration ensemble can be effectively employed in a neural network-based PDE solver to obtain an approximation solution of a restricted BKE, even when the dimension of the problem is very high. The resulting solution provides an approximation for the committor function that encodes mechanistic information for the reaction, paving a new way for understanding of complex chemical reactions and evaluation of reaction rates.

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

Approaching periodic systems in ensemble density functional theory via finite one-dimensional models

Ensemble density functional theory (EDFT) is a generalization of ground-state DFT, which is based on an exact formal theory of finite collections of a system's ground and excited states. EDFT in various forms has been shown to improve the accuracy of calculated energy level differences in isolated model systems, atoms, and molecules, but it is not yet clear how EDFT could be used to calculate band gaps for periodic systems. We extend the application of EDFT toward periodic systems by estimating the thermodynamic limit with increasingly large finite one-dimensional 'particle in a box' systems, which approach the uniform electron gas (UEG). Using ensemble-generalized Hartree and local spin density approximation exchange-correlation functionals, we find that corrections go to zero in the infinite limit, as expected for a metallic system. However, there is a correction to the effective mass, with results comparable to other calculations on 1D, 2D, and 3D UEGs, which indicates promise for non-trivial results from EDFT on periodic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Machine-learning Kohn–Sham potential from dynamics in time-dependent Kohn–Sham systems

Abstract The construction of a better exchange-correlation potential in time-dependent density functional theory (TDDFT) can improve the accuracy of TDDFT calculations and provide more accurate predictions of the properties of many-electron systems. Here, we propose a machine learning method to develop the energy functional and the Kohn–Sham potential of a time-dependent Kohn–Sham (TDKS) system is proposed. The method is based on the dynamics of the Kohn–Sham system and does not require any data on the exact Kohn–Sham potential for training the model. We demonstrate the results of our method with a 1D harmonic oscillator example and a 1D two-electron example. We show that the machine-learned Kohn–Sham potential matches the exact Kohn–Sham potential in the absence of memory effect. Our method can still capture the dynamics of the Kohn–Sham system in the presence of memory effects. The machine learning method developed in this article provides insight into making better approximations of the energy functional and the Kohn–Sham potential in the TDKS system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Deep Learning without Global Optimization by Random Fourier Neural Networks

Here we introduce a new training algorithm for deep neural networks that utilize random complex exponential activation functions. Our approach employs a Markov chain Monte Carlo sampling procedure to iteratively train network layers, avoiding global and gradient-based optimization while maintaining error control. It consistently attains the theoretical approximation rate for residual networks with complex exponential activation functions, determined by network complexity. Additionally, it enables efficient learning of multiscale and high-frequency features, producing interpretable parameter distributions. Despite using sinusoidal basis functions, we do not observe Gibbs phenomena in approximating discontinuous target functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Generation, dynamics, and correlations of the fission fragments' angular momenta

The generation of angular momentum in fissioning nuclei is not well understood. The predictions of different models disagree, particularly concerning the correlation between the fragments' angular momenta. In this article, a time-dependent collective Hamiltonian model is proposed to treat the generation of the angular momentum in the fission fragments due to the quantum uncertainty principle as well as the dynamics of the collective wave function during and after scission. The model is constructed in the framework of the frozen Hartree-Fock approximation using a Skyrme energy functional to extract deformations of the fission fragments as well as the interactions in a derived collective Hamiltonian. The fission reactions studied are $^{240}$Pu $\rightarrow$ $^{132}$Sn+$^{108}$Ru and $^{240}$Pu $\rightarrow$ $^{144}$Ba+$^{96}$Sr. The model can account for a large part of the angular momentum found in experimental data. The orientation of the angular momentum of each fragment is found to be mainly in the plane perpendicular to the fission axis, in agreement with the experiment. The magnitudes of the angular momenta in the two fragments are nearly uncorrelated, in agreement with the recent experimental data of Wilson et al., Nature (London) 590, 566 (2021). Some of the conclusions of the traditional collective vibration model are supported by the present model but some are not. Surprisingly, it is found that the angular momenta of the fragments are slightly correlated positively as in a wriggling mode. As a result, it is also found that the presence of an octupole deformation in a fragment can significantly increase the generated angular momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗