Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kernel function”

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

Helicity-dependent generalization of the JIMWLK evolution

We generalize the JIMWLK evolution equation to include the small-$x$ evolution of helicity distributions. To derive helicity JIMWLK, we use the method devised by A.H. Mueller in [Phys. Lett. B 523, 243 (2001)] for the usual unpolarized JIMWLK, and apply it to the helicity evolution equations derived in J. High Energy Phys. 01 (2016) 072; Phys. Rev. D 95, 014033 (2017); J. High Energy Phys. 10 (2017) 198; Phys. Rev. D 99, 054032 (2019)]. We obtain a functional evolution equation for the generalized JIMWLK weight functional. The functional now is a sum of a helicity-independent term and a helicity-dependent term. The kernel of the new evolution equation consists of the standard eikonal leading-order JIMWLK kernel plus new terms containing derivatives with respect to the sub-eikonal quark and gluon fields. The new kernel we derive can be used to generate the standard JIMWLK evolution and to construct small-$x$ evolution equations for flavor-singlet operators made of any number of light-cone Wilson lines along with one Wilson line containing sub-eikonal helicity-dependent local operator insertions (a "polarized Wilson line").

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Convergence of Weak-SINDy Surrogate Models

In this paper, we give an in-depth error analysis for surrogate models generated by a variant of the Sparse Identification of Nonlinear Dynamics (SINDy) method. We start with an overview of a variety of nonlinear system identification techniques, namely SINDy, weak-SINDy, and the occupation kernel method. Under the assumption that the dynamics are a finite linear combination of a set of basis functions, these methods establish a linear system to recover coefficients. We illuminate the structural similarities between these techniques and establish a projection property for the weak-SINDy technique. Following the overview, we analyze the error of surrogate models generated by a simplified version of weak-SINDy. In particular, under the assumption of boundedness of a composition operator given by the solution, we show that (i) the surrogate dynamics converges towards the true dynamics and (ii) the solution of the surrogate model is reasonably close to the true solution. Finally, as an application, we discuss the use of a combination of weak-SINDy surrogate modeling and proper orthogonal decomposition (POD) to build a surrogate model for partial differential equations (PDEs).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Experience of Migrating a Parallel Graph Coloring Program from CUDA to SYCL

We describe the experience of converting a CUDA implementation of a parallel graph coloring algorithm to SYCL. The goals are for our work to be useful to application and compiler developers by providing a detailed description of migration paths between CUDA and SYCL. We will describe how CUDA functions are mapped to SYCL functions. Evaluating the CUDA and SYCL implementations of the algorithm shows that the performance of SYCL and CUDA kernels are comparable over the test graph set on NVIDIA P100 and V100 GPUs. The SYCL program also allows for performance evaluation with the OpenCL and Level Zero interfaces and power profiling on an Intel GPU computing platform.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An Approach to Characterizing Spatial Aspects of Image System Blur

In imaging applications, such as X-ray radiography, the experimental system setup and data capture process introduce system blur and distortion. In images taken with thick scintillators, a spatially varying system blur is observed, requiring a nonstandard approach to blur and distortion estimation and removal. We demonstrate an approach to estimation of local point spread functions using checkerboard radiographs collected at the Naval Research Laboratory and extend that into a map defining the kernel over the camera field of view.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

A fast and accurate method for perturbative resummation of transverse momentum-dependent observables

Here, we propose a novel strategy for the perturbative resummation of transverse momentum-dependent (TMD) observables, using the q T spectra of gauge bosons (γ*, Higgs) in pp collisions in the regime of low (but perturbative) transverse momentum q T as a specific example. First we introduce a scheme to choose the factorization scale for virtuality in momentum space instead of in impact parameter space, allowing us to avoid integrating over (or cutting off) a Landau pole in the inverse Fourier transform of the latter to the former. The factorization scale for rapidity is still chosen as a function of impact parameter b, but in such a way designed to obtain a Gaussian form (in ln b) for the exponentiated rapidity evolution kernel, guaranteeing convergence of the b integral. We then apply this scheme to obtain the q T spectra for Drell-Yan and Higgs production at NNLL accuracy. In addition, using this scheme we are able to obtain a fast semi-analytic formula for the perturbative resummed cross sections in momentum space: analytic in its dependence on all physical variables at each order of logarithmic accuracy, up to a numerical expansion for the pure mathematical Bessel function in the inverse Fourier transform that needs to be performed just once for all observables and kinematics, to any desired accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Real Time Monitoring of the Dynamic Intracluster Diffusion of Single Gold Atoms into Silver Nanoclusters

Alloying metal materials with heterometal atoms is an efficient way to diversify the function of materials, but in-depth understanding of the dynamic heterometallic diffusion inside the alloying materials is rather limited, especially at the atomic level. Here, we report the real-time monitoring of the dynamic diffusion process of a single gold (Au) atom into an atomically precise silver nanocluster (Ag NC), Ag 25 (MHA) 18 (MHA = 6-mercaptohexanoic acid), by using in situ UV–vis absorption spectroscopy in combination with mass and tandem mass spectrometry. We found that the Au heteroatom first replaces the Ag atom at the surface Ag 2 (MHA) 3 motifs of Ag 25 (MHA) 18 . After that, the Au atom diffuses into the surface layer of the icosahedral Ag 13 kernel and finally occupies the center of the alloy NCs to form the thermodynamically stable [email protected] 24 (MHA) 18 product. Density functional theory (DFT) calculations reveal that the key thermodynamic driving force is the preference of the Au heteroatom for the central site of alloy NCs. Furthermore, the real-time monitoring approach developed in this study could also be extended to other metal alloy systems to reveal the reaction dynamics of intracluster diffusion of heteroatoms, as well as the formation mechanisms of metal alloy nanomaterials.

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

Functional characterization of three GH10 xylanases

Decreased availability of fossil fuels has driven an increased demand for sustainable energy sources such as biofuels. One area of bioenergy research involves development of biocatalysts able to efficiently degrade biomass “waste”‐products, such as substrates derived from corn stover, rather than easily‐fermentable biomass used for food, such as starch from corn kernels. Xylanases are of interest for their ability to hydrolyze hemicellulosic sugars found in plant cell walls, allowing utilization of waste‐derived biomass. In this study, three xylanases of glycoside hydrolase family 10 (GH10) are functionally characterized and optimized for temperature and pH, as well as kinetically analyzed for xylose monomer and oligomer production. These enzymes, GH10‐4, GH10‐9, and GH10‐15, demonstrate activity upon substrates xylan and arabinoxylan. For GH10‐4, optima are 50°C at pH 8 on xylan and 45°C at pH 8 on arabinoxylan. For GH10‐9, optima are 70°C at pH 5 on xylan and 70°C at pH 5 on arabinoxylan. For GH10‐15, optima are 50°C at pH 8 on xylan and 45°C at pH 5 on arabinoxylan. Kinetic analysis reveals GH10‐4 produces xylose (X1), xylobiose (X2), xylotriose (X3), xylotetraose (X4), and xylopentaose (X5) from xylan, while GH10‐9 and GH10‐15 produce X1, X2, and X3. From arabinoxylan, GH10‐4 produces X1, X2, X3, X4, while GH10‐9 and GH10‐15 produce X1, X2, X3. Conclusive kinetic data may suggest applicability of these enzymes for industrial use within biocatalysts and contribute further insight to the functionality of enzymes within GH10. Support or Funding Information The authors are grateful for financial support by the Integrated Biological Sciences Summer Research Program (IBS‐SRP, NSF) at the University of Wisconsin‐Madison and the Great Lakes Bioenergy Research Center (GLBRC) of the U.S. Department of Energy. Presentation travel support is provided by the Ronald E. McNair Scholars Program of Winthrop University. The authors also thank Lai Bergeman, Craig Bingman, Bob Smith, Emily Beebe, Kirk Vander Meulen, Amber Smith, David Wassarman, and John Greenler for professional support.

Park, Jesslyn↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Training quantum neural networks using the quantum information bottleneck method

Abstract We provide in this paper a concrete method for training a quantum neural network to maximize the relevant information about a property that is transmitted through the network. This is significant because it gives an operationally well founded quantity to optimize when training autoencoders for problems where the inputs and outputs are fully quantum. We provide a rigorous algorithm for computing the value of the quantum information bottleneck quantity within error ε that requires O ( log 2 ⁡ ( 1 / ϵ ) + 1 / δ 2 ) queries to a purification of the input density operator if its spectrum is supported on { 0 } ⋃ [ δ , 1 − δ ] for δ > 0 and the kernels of the relevant density matrices are disjoint. We further provide algorithms for estimating the derivatives of the QIB function, showing that quantum neural networks can be trained efficiently using the QIB quantity given that the number of gradient steps required is polynomial.

Çatlı, Ahmet Burak (ORCID:0000000152294141)↗
DOE OSTI2022

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. In this work, we introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Solving Inverse Stochastic Problems from Discrete Particle Observations Using the Fokker--Planck Equation and Physics-Informed Neural Networks

The Fokker--Planck (FP) equation governing the evolution of the probability density function (PDF) is applicable to many disciplines, but it requires specification of the coefficients for each case, which can be functions of space-time and not just constants and hence require the development of a data-driven modeling approach. When the data available is directly on the PDF, there exist methods for inverse problems that can be employed to infer the coefficients and thus determine the FP equation and subsequently obtain its solution. Herein, we address a more realistic scenario, where only sparse data are given on the particles' positions at a few time instants, which are not sufficient to accurately construct directly the PDF even at those times from existing methods, e.g., kernel estimation algorithms. To this end, we develop a general framework based on physics-informed neural networks (PINNs) that introduces a new loss function using the Kullback--Leibler divergence to connect the stochastic samples with the FP equation to simultaneously learn the equation and infer the multidimensional PDF at all times. In particular, we consider two types of inverse problems, type I, where the FP equation is known but the initial PDF is unknown, and type II, in which, in addition to the unknown initial PDF, the drift and diffusion terms are also unknown. In both cases, we investigate problems with either Brownian or Lévy noise or a combination of both. Here, we demonstrate the new PINN framework in detail in the one-dimensional (1D) case, but we also provide results for up to five dimensions demonstrating that we can infer both the FP equation and dynamics simultaneously at all times with high accuracy using only very few discrete observations of the particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

THREE-POINT PHASE CORRELATIONS: A NEW MEASURE OF NONLINEAR LARGE-SCALE STRUCTURE

We derive an analytical expression for a novel large-scale structure observable: the line correlation function. The line correlation function, which is constructed from the three-point correlation function of the phase of the density field, is a robust statistical measure allowing the extraction of information in the nonlinear and non-Gaussian regime. We show that, in perturbation theory, the line correlation is sensitive to the coupling kernel F{sub 2}, which governs the nonlinear gravitational evolution of the density field. We compare our analytical expression with results from numerical simulations and find a 1σ agreement for separations r ≳ 30 h{sup −1} Mpc. Fitting formulae for the power spectrum and the nonlinear coupling kernel at small scales allow us to extend our prediction into the strongly nonlinear regime, where we find a 1σ agreement with the simulations for r ≳ 2 h{sup −1} Mpc. We discuss the advantages of the line correlation relative to standard statistical measures like the bispectrum. Unlike the latter, the line correlation is independent of the bias, in the regime where the bias is local and linear. Furthermore, the variance of the line correlation is independent of the Gaussian variance on the modulus of the density field. This suggests that the line correlation can probe more precisely the nonlinear regime of gravity, with less contamination from the power spectrum variance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

While time-dependent density functional theory has emerged as a method of choice for computing electronic spectra and response of molecules and materials, its reliability is hindered by the adiabatic approximation that is commonly made. In this work, we address one problematic aspect that arises from this approximation: an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. Here, we demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2020

Transverse momentum dependent PDFs at N3LO

We compute the quark and gluon transverse momentum dependent parton distribution functions at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. Our calculation is based on an expansion of the differential Drell-Yan and gluon fusion Higgs production cross sections about their collinear limit. This method allows us to employ cutting edge multiloop techniques for the computation of cross sections to extract these universal building blocks of the collinear limit of QCD. The corresponding perturbative matching kernels for all channels are expressed in terms of simple harmonic polylogarithms up to weight five. As a byproduct, we confirm a previous computation of the soft function for transverse momentum factorization at N 3 LO. Our results are the last missing ingredient to extend the q T subtraction methods to N 3 LO and to obtain resummed q T spectra at N 3 LL' accuracy both for gluon as well as for quark initiated processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Building workflows for an interactive human-in-the-loop automated experiment (hAE) in STEM-EELS

Exploring the structural, chemical, and physical properties of matter on the nano- and atomic scales has become possible with the recent advances in aberration-corrected electron energy-loss spectroscopy (EELS) in scanning transmission electron microscopy (STEM). However, the current paradigm of STEM-EELS relies on the classical rectangular grid sampling, in which all surface regions are assumed to be of equal a priori interest. However, this is typically not the case for real-world scenarios, where phenomena of interest are concentrated in a small number of spatial locations, such as interfaces, structural and topological defects, and multi-phase inclusions. One of the foundational problems is the discovery of nanometer- or atomic-scale structures having specific signatures in EELS spectra. Herein, we systematically explore the hyperparameters controlling deep kernel learning (DKL) discovery workflows for STEM-EELS and identify the role of the local structural descriptors and acquisition functions in experiment progression. In agreement with the actual experiment, we observe that for certain parameter combinations the experiment path can be trapped in the local minima. We demonstrate the approaches for monitoring the automated experiment in the real and feature space of the system and knowledge acquisition of the DKL model. Based on these, we construct intervention strategies defining the human-in-the-loop automated experiment (hAE). This approach can be further extended to other techniques including 4D STEM and other forms of spectroscopic imaging. The hAE library is available on Github at https://github.com/utkarshp1161/hAE/tree/main/hAE.

Pratiush, Utkarsh [Univ. of Tennessee, Knoxville, ↗
DOE OSTI2022

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Implementation of a small-angle scattering model in MCNPX for very cold neutron reflector studies

Current neutron moderator media do not sufficiently moderate neutrons below the cold neutron regime into the very cold neutron (VCN) regime that is desirable for some physics applications. Nesvizhevsky et al [1] have demonstrated that nanodiamond powder efficiently reflect VCN via small angle scattering. He suggests that these effects could be exploited to boost the neutron output of a VCN moderator. Simulation studies of nanoparticle reflectors are being investigated as part of the development of a VCN source option for the SNS second target station. We are pursuing an expansion of the MCNPX code by implementation of an analytical small-angle scattering function [2], which is adaptable by scattering particle sizes, distributions, and packing fractions in order to supplement currently existing scattering kernels. The analytical model and preliminary studies using MCNPX will be discussed.

Grammer, Kyle↗