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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

QCD factorization with multihadron fragmentation functions

Important aspects of quantum chromodynamics (QCD) factorization theorems are the properties of the objects involved that can be identified as universal. One example is that the definitions of parton densities and fragmentation functions for different types of hadrons differ only in the identity of the nonperturbative states that form the matrix elements, but are otherwise the same. This leads to independence of perturbative calculations on nonperturbative details of external states. It also lends support to interpretations of correlation functions as encapsulations of intrinsic nonperturbative properties. These characteristics have usually been presumed to still hold true in fragmentation functions even when the observed nonperturbative state is a small-mass cluster of n hadrons rather than simply a single isolated hadron. However, the multidifferential aspect of cross sections that rely on these latter types of fragmentation functions complicates the treatment of kinematical approximations in factorization derivations. That has led to recent claims that the operator definitions for fragmentation functions need to be modified from the single hadron case with nonuniversal prefactors. With such concerns as our motivation, we retrace the steps for factorizing the unpolarized semi-inclusive e + e − annihilation cross section and confirm that they do apply without modification to the case of a small-mass multihadron observed in the final state. In particular, we verify that the standard operator definition from single hadron fragmentation, with its usual prefactor, remains equally valid for the small-mass n -hadron case with the same hard parts and evolution kernels, whereas the more recently proposed definitions with nonuniversal prefactors do not. Our results reaffirm the reliability of most past phenomenological applications of dihadron fragmentation functions. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Simple Models of Lateral Heterostructures

General analytical expressions for densities of states of a lateral heterostructure, formed by the contact of two square semi-infinite lattices with single-band and two-band spectra, were obtained in the tight-binding approximation by the Green’s function method. The semi-elliptical density of states was used for numerical estimates, and the model of two interacting dimers was proposed to estimate the charge transfer. Application of this approach to description of lateral epitaxial and graphene-like heterostructures is discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Approximating and incorporating model uncertainty in an inversion for seismic source functions: Preliminary results

We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.

42 ENGINEERING↗
DOE OSTI2020

A Multi-Function AAA Algorithm Applied to Frequency-Dependent Line Modeling: Preprint

We present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages respect to the popular Vector Fitting (VF) algorithm using examples of approximations of admittance matrices discussed in a very recent paper and by the creators of VF. We show that not only MIAAA outperforms VF in terms of achieving better accuracy using a smaller number of poles, but also has no numerical issues to achieve convergence. In contrast to VF, MIAA is not sensitive to the location of input sample points and it does not require good estimates for the location of the desired approximation poles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface

In this work, we explore the suitability of several density functionals with the generalized gradient approximation (GGA) and beyond for describing the dissociative chemisorption of methane on the reconstructed Pt(110)-(2×1) surface. The bulk and surface structures of the metal, methane adsorption energy, and dissociation barrier are used to assess the functionals. A van der Waals corrected GGA functional (optPBE-vdW) and a metaGGA functional with van der Waals correction (MS PBEl-rVV10) are selected for ab initio molecular dynamics calculations of the sticking probability. Our results suggest that the use of these two functionals may lead to a better agreement with existing experimental results, thus serving as a good starting point for future development of reliable machine-learned potential energy surfaces for the dissociation of methane on the Pt(110)-(2×1) surface.

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

Status of free-energy representations for the homogeneous electron gas

Renewed interest in the homogeneous electron gas (HEG) has been stimulated by recent accurate simulations of it over a wide range of densities and temperatures. Those data, combined with known theoretical limits, have led to analytical representations of the free energy. Such a representation is, at least in principle, the complete HEG equation of state. The initial objective here is to establish that the two best representations of the simulation data and constraints are effectively the same in both functional form and accuracy of representation. The second objective is to disclose and delineate a significant difficulty. Despite their expected accuracy for the free energy, the underlying functional form is not adequate for derived thermodynamic properties. As an example, the specific heats obtained from the representations exhibit anomalies suggesting the need first for additional simulation data in critical regimes, then for refined fitting functions. Here, the existing representations are, however, almost certainly adequate for applications based on the free energy alone (e.g., density functional theory for warm dense matter). The third objective is to show that, despite their inability to provide a complete thermodynamic description of the HEG, the best analytical representations do provide a fully adequate exchange-correlation local density approximation for free energy density functional calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-principles investigation of uranium mononitride (UN): Effect of magnetic ordering, spin-orbit interactions and exchange correlation functional

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys. Various properties of UN have been previously studied using different density functional theory (DFT) methodologies; however, there are still inconsistencies when it comes to the dynamical stability and defect properties of UN. We address these inconsistencies by studying the UN phonons and defect properties using DFT calculations employing two generalized gradient approximation (GGA) exchange-correlation functionals: PBE and AM05, with and without an added on-site Coulomb repulsion term (+U). Furthermore, we investigate the importance of spin-orbit coupling (SOC) when calculating the properties of UN. We use the different methodologies to determine the preference of UN to have antiferromagnetic (AFM) ordering, as seen in experiments, or ferromagnetic (FM) ordering of the uranium spins. We compare the crystallographic properties, density of states, the DFT X-ray photoelectron spectra and phonon dispersions calculated using the different methodologies. We demonstrate that GGA + U reproduces the AFM ordering in UN, but the crystal structure is dynamically unstable. We also show that magnetic ordering is important in finding the lowest energy defective structure, and that SOC has a distinct influence on the energy of the different uranium interstitial defects. Lastly, we discuss the point defect formation energies under U-rich and N-rich conditions, and the stoichiometric formation energies calculated with the different methodologies, providing insight into the observed tendency for forming hypostoichiometric UN.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron scattering on A = 3 nuclei from quantum Monte Carlo based approaches

We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electron scattering on ${\mathbf{A=3}}$ nuclei from quantum Monte Carlo based approaches

We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Predicting the thermodynamic stability of double-perovskite halides from density functional theory

Recently, a series of double-perovskite halide compounds such as Cs 2 AgBiCl 6 and Cs 2 AgBiBr 6 have attracted intensive interest as promising alternatives to the solar absorber material CH 3 NH 3 PbI 3 because they are Pb-free and may exhibit enhanced stability. The thermodynamic stability of a number of double-perovskite halides has been predicted based on density functional theory (DFT) calculations of compound formation energies. In this paper, we found that the stability prediction can be dependent on the approximations used for the exchange-correlation functionals, e.g., the DFT calculations using the widely used Perdew, Burke, Ernzerhof (PBE) functional predict that Cs 2 AgBiBr 6 is thermodynamically unstable against phase-separation into the competing phases such as AgBr, Cs 2 AgBr 3 , Cs 3 Bi 2 Br 9 , etc., obviously inconsistent with the good stability observed experimentally. The incorrect prediction by the PBE calculation results from its failure to predict the correct ground-state structures of AgBr, AgCl, and CsCl. By contrast, the DFT calculations based on local density approximation, optB86b-vdW, and optB88-vdW functionals predict the ground-state structures of these binary halides correctly. Furthermore, the optB88-vdW functional is found to give the most accurate description of the lattice constants of the double-perovskite halides and their competing phases. Given these two aspects, we suggest that the optB88-vdW functional should be used for predicting thermodynamic stability in the future high-throughput computational material design or the construction of the Materials Genome database for new double-perovskite halides. As a result, using different exchange-correlation functionals has little influence on the dispersion of the conduction and the valence bands near the electronic bandgap; however, the calculated bandgap can be affected indirectly by the optimized lattice constant, which varies for different functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Stochastic particle production in a de Sitter background

Here, we explore non-adiabatic particle production in a de Sitter universe for a scalar spectator field, by allowing the effective mass m 2 (t) of this field and the cosmic time interval between non-adiabatic events to vary stochastically. Two main scenarios are considered depending on the (non-stochastic) mass M of the spectator field: the conformal case with M 2 =2H 2 , and the case of a massless field. We make use of the transfer matrix formalism to parametrize the evolution of the system in terms of the "occupation number", and two phases associated with the transfer matrix; these are used to construct the evolution of the spectator field. Assuming short-time interactions approximated by Dirac-delta functions, we numerically track the change of these parameters and the field in all regimes: sub- and super-horizon with weak and strong scattering. In all cases a log-normally distributed field amplitude is observed, and the logarithm of the field amplitude approximately satisfies the properties of a Wiener process outside the horizon. We derive a Fokker-Planck equation for the evolution of the transfer matrix parameters, which allows us to calculate analytically non-trivial distributions and moments in the weak-scattering limit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Theoretical Investigation of Structural, Electronic, Magnetic, and Mechanical Properties of Quaternary Heusler Alloys CoVTiX (X = As, Si)

The structural, electronic, magnetic, and mechanical properties of CoVTiAs and CoVTiSi quaternary Heusler alloys are studied using the full-potential linearized augmented plane-wave (FP-LAPW) method in the framework of the density functional theory (DFT). The generalized gradient approximation (GGA) is chosen for the exchange-correlation energy, whereas the modified Becke-Johnson (mBJ) formalism is applied for the electronic properties. Both CoVTiAs and CoVTiSi are stable in the type 2 (FM) structure. The mechanical properties of CoVTiX (X = As, Si) are predicted from the calculated elastic constants. Results identify that these compounds are mechanically stable. Our results with the mBJ approximation predict that CoVTiAs and CoVTiSi alloys are half-metallic ferromagnets with band gaps of 0.450 and 0.424 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Effect of energy dependence of the density of states on pressure-dependent rate constants

The F E integral for the normalized Boltzmann-weighted number of unimolecular states above the threshold energy is the key quantity for computing the collision efficiency in the pressure-dependent unimolecular rate theory developed by Troe, who calls this the energy dependence factor of the density of states. By using the Whitten–Rabinovitch approximation and assuming that the Whitten–Rabinovitch a(E) function is independent of energy, F E can be approximated by an analytical formula; this approximate formula is widely used because of its convenience and computational efficiency. Here we test its validity by comparing the rate constants computed by using the approximate F E to the ones determined by using the numerically integrated F E . For small-sized molecules and for reactions with high threshold energies E 0 , the differences are negligible at all temperatures, but in other cases, the approximate formula tends to underestimate F E and thus overestimates the collision efficiency, and this leads to smaller pressure falloff. When a(E) at high energies differs appreciably from a(E 0 ), we find that the underestimation of pressure-dependent rate constants by using the approximate formula can be greater than a factor of 5 at high temperatures. The physical insight we draw from this study is that, for reactions with threshold energies below about 30 kcal mol -1 , the rate of collisional energy transfer can be appreciably slowed down by the increase in the density of states at higher energies, and this increases the falloff effect by which finite-pressure rate constants become lower than the high-pressure limit, especially at higher temperatures.

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

Find the Gap: Black Hole Population Analysis with an Astrophysically Motivated Mass Function

We introduce a novel black hole mass function that realistically models the physics of pair-instability supernovae with a minimal number of parameters. Applying this to all events in the LIGO-Virgo Gravitational-Wave Transient Catalog 2 (GWTC-2), we detect a peak at M BHMG = 46$^{+17}_{-6}$ M ⊙ . Repeating the analysis without the black holes from the event GW190521, we find this feature at M BHMG =54± 6 M ⊙ . These results establish the edge of the anticipated "black hole mass gap" at a value compatible with the expectation from standard stellar structure theory. The mass gap manifests itself as a discontinuity in the mass function and is populated by a distinct, less-abundant population of higher-mass black holes. We find that the primary black hole population scales with power-law index -1.95 ± 0.51 (-1.97 ± 0.44) with (without) GW190521, consistent with models of star formation. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A new computational approach for the solutions of generalized pantograph-delay differential equations

In this study, a new computational approach is presented to solve the generalized pantograph-delay differential equations (PDDEs). The solutions obtained by our scheme represent by a linear combination of a special kind of basis functions, and can be deduced in a straightforward manner. Firstly, using the least squares approximation method and the Lagrange-multiplier method, the given PDDE is converted to a linear system of algebraic equations, and those unknown coefficients of the solution of the problem are determined by solving this linear system. Secondly, a PDDE related to the error function of the approximate solution is constructed based on the residual error function technique, and error estimation is presented for the suggested method. The convergence of the approximate solution is proved. Several numerical examples are given to demonstrate the accuracy and efficiency. Comparisons are made between the proposed method and other existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Post-Newtonian kinetic theory

A kinetic theory for relativistic gases in the presence of gravitational fields is developed in the second post-Newtonian approximation. The corresponding Boltzmann equation is determined from the evolution of the one-particle distribution function with respect to the proper time along the world line of the particle. From the knowledge of the equilibrium Maxwell–Jüttner distribution function in the second post-Newtonian approximation the components of the particle four-flow and energy–momentum tensor are obtained. The Eulerian hydrodynamic equations for the mass density, mass-energy density and momentum density in the second post-Newtonian approximation are determined from the Boltzmann equation. It is shown that the combination of the hydrodynamic equations of mass and mass-energy densities leads to the hydrodynamic equation for the internal energy density in the first post-Newtonian approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multioutput Convolutional Neural Network for Improved Parameter Extraction in Time-Resolved Electrostatic Force Microscopy Data

Time-resolved scanning probe microscopy methods, like time-resolved electrostatic force microscopy (trEFM), enable imaging of dynamic processes ranging from ion motion in batteries to electronic dynamics in microstructured thin film semiconductors for solar cells. Reconstructing the underlying physical dynamics from these techniques can be challenging due to the interplay of cantilever physics with the actual transient kinetics of interest in the resulting signal. Previously, quantitative trEFM used empirical calibration of the cantilever or feed-forward neural networks trained on simulated data to extract the physical dynamics of interest. Both these approaches are limited by interpreting the underlying signal as a single exponential function, which serves as an approximation but does not adequately reflect many realistic systems. Here, we present a multi-branched, multi-output convolutional neural network (CNN) that uses the trEFM signal in addition to the physical cantilever parameters as input. The trained CNN accurately extracts parameters describing both single-exponential and bi-exponential underlying functions, and more accurately reconstructs real experimental data in the presence of noise. This article demonstrates an application of physics-informed machine learning to complex signal processing tasks, enabling more efficient and accurate analysis of trEFM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗