Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized functions”

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

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2024

Definition of fragmentation functions and the violation of sum rules

We point out a problem with the formulation and derivations of sum rules for quark fragmentation functions that impacts their validity in QCD, but which potentially points toward an improved understanding of final states in inclusive hard processes. Fragmentation functions give the distribution of final-state hadrons arising from a parton exiting a hard scattering, and the sum rules for momentum, electric charge, etc. express conservation of these quantities. The problem arises from a mismatch between the quark quantum numbers of the initial quark and the fact that all observed final-state hadrons are confined bound states with color zero. We point that, in a confining theory like QCD, the Wilson line in the operator definition of a fragmentation function entails that the final state in a fragmentation function includes a bound state in the external field generated by the Wilson line. We justify this with the aid of general features of string hadronization. The anomalous bound states are restricted to fractional momentum z = 0 . They tend to invalidate sum rules like the one for charge conservation when applied to the fragmentation functions inferred from experimental data but not the momentum sum rule. We propose to exploit our ideas in future studies as a way to relate the fragmentation functions extracted from inclusive cross sections to more detailed nonperturbative descriptions of final state hadronization. We also describe scenarios wherein the traditional sum rules might remain approximately valid with a reasonably high degree of accuracy. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

INSURE: An Information Theory iNspired diSentanglement and pURification modEl for Domain Generalization

Domain Generalization (DG) aims to learn a generalizable model on the unseen target domain by only training on the multiple observed source domains. Although a variety of DG methods have focused on extracting domain-invariant features, the domain-specific class-relevant features have attracted attention and been argued to benefit generalization to the unseen target domain. To take into account the class-relevant domain-specific information, in this paper we propose an Information theory iNspired diSentanglement and pURification modEl (INSURE) to explicitly disentangle the latent features to obtain sufficient and compact (necessary) class-relevant feature for generalization to the unseen domain. Specifically, we first propose an information theory inspired loss function to ensure the disentangled class-relevant features contain sufficient class label information and the other disentangled auxiliary feature has sufficient domain information. Additionally, we further propose a paired purification loss function to let the auxiliary feature discard all the class-relevant information and thus the class-relevant feature will contain sufficient and compact (necessary) class-relevant information. Moreover, instead of using multiple encoders, we propose to use a learnable binary mask as our disentangler to make the disentanglement more efficient and make the disentangled features complementary to each other. We conduct extensive experiments on five widely used DG benchmark datasets including PACS, VLCS, OfficeHome, TerraIncognita, and DomainNet. The proposed INSURE achieves state-of-the-art performance. We also empirically show that domain-specific class-relevant features are beneficial for domain generalization. The code is available at https://github.com/yuxi120407/INSURE .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-dimensional methods for quantum homodyne tomography

We provide optimized recursion relations for homodyne tomography. We improve previous methods by mitigating the divergences intrinsic in the calculation of the pattern functions used previously, and detail how to implement the data analysis through Monte Carlo simulations. Our refinements are necessary for the reconstruction of excited quantum states which populate a high-dimensional subspace of the electromagnetic field Hilbert space. •Stabilization of numerics for quantum homodyne tomography reconstructions.•High dimensional quantum homodyne tomography reconstruction.•Monte Carlo simulations of quantum homodyne tomography.•Quantum Homodyne tomography simulations.

Mosco, Nicola↗
DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

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

A resolution independent neural operator

The Deep operator network (DeepONet) is a powerful yet simple neural operator architecture that utilizes two deep neural networks to learn mappings between infinite-dimensional function spaces. This architecture is highly flexible, allowing the evaluation of the solution field at any location within the desired domain. However, it imposes a strict constraint on the input space, requiring all input functions to be discretized at the same locations; this limits its practical applications. Here, in this work, we introduce a general framework for operator learning from input–output data with arbitrary number and locations of sensors. This begins by introducing a resolution-independent DeepONet (RI-DeepONet), enabling it to handle input functions that are arbitrarily, but sufficiently finely, discretized. To this end, we propose two dictionary learning algorithms to adaptively learn a set of appropriate continuous basis functions, parameterized as implicit neural representations (INRs), from correlated signals defined on arbitrary point cloud data. These basis functions are then used to project arbitrary input function data as a point cloud onto an embedding space (i.e., a vector space of finite dimensions) with dimensionality equal to the dictionary size, which can be directly used by DeepONet without any architectural changes. In particular, we utilize sinusoidal representation networks (SIRENs) as trainable INR basis functions. The introduced dictionary learning algorithms are then used in a similar way to learn an appropriate dictionary of basis functions for the output function data, which defines a new neural operator architecture referred to as the R esolution I ndependent N eural O perator (RINO). In the RINO, the operator learning task simplifies to learning a mapping from the coefficients of input basis functions to the coefficients of output basis functions. We demonstrate the robustness and applicability of RINO in handling arbitrarily (but sufficiently richly) sampled input and output functions during both training and inference through several numerical examples.

Deep operator network (DeepONet)↗
DOE OSTI2024

Generalized Quantum Signal Processing

Quantum signal processing (QSP) and quantum singular value transformation (QSVT) currently stand as the most efficient techniques for implementing functions of block-encoded matrices, a central task that lies at the heart of most prominent quantum algorithms. However, current QSP approaches face several challenges, such as the restrictions imposed on the family of achievable polynomials and the difficulty of calculating the required phase angles for specific transformations. In this paper, we present a generalized quantum signal processing (GQSP) approach, employing general SU(2) rotations as our signal-processing operators, rather than relying solely on rotations in a single basis. Our approach lifts all practical restrictions on the family of achievable transformations, with the sole remaining condition being that | P | ≤ 1 , a restriction necessary due to the unitary nature of quantum computation. Furthermore, GQSP provides a straightforward recursive formula for determining the rotation angles needed to construct the polynomials in cases where P and Q are known. In cases where only P is known, we provide an efficient optimization algorithm capable of identifying in under a minute of GPU time, a corresponding Q for polynomials of degree on the order of 10 7 . We further illustrate GQSP simplifies QSP-based strategies for Hamiltonian simulation, offer an optimal solution to the ϵ -approximate fractional query problem that requires O ( ( 1 / δ ) + log ( 1 / ϵ ) ) queries to perform where O ( 1 / δ ) is a proved lower bound, and introduces novel approaches for implementing bosonic operators. Moreover, we propose a novel framework for the implementation of normal matrices, demonstrating its applicability through synthesis of diagonal matrices, as well as the development of a new algorithm for convolution through synthesis of circulant matrices using only O ( d log N + log 2 N ) 1 and 2-qubit gates for a filter of lengths d . Published by the American Physical Society 2024

Motlagh, Danial↗
DOE OSTI2022

How Well Does Kohn–Sham Regularizer Work for Weakly Correlated Systems?

Kohn–Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn–Sham density functional theory that works for strongly correlated systems. Here we test KSR for a weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, and Be 2+ and testing on 1D hydrogen chains, LiH, BeH 2 , and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 mH.

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

A General Binary Isotherm Model for Amines Interacting with CO 2 and H 2 O

CO 2 capture by primary or secondary amines has been of great research interests for a century because of its industrial importance. Interest has grown even more, because of the need to eliminate the CO 2 emissions that drive global warming. Experimental evidence shows that CO 2 sorption in primary or secondary amines is accompanied by co-absorption of H 2 O. A quantitative analysis of such CO 2 -H 2 O co-absorption behavior is important for practical process design and theoretical understanding. Even though there is almost an experimental consensus that water enhances CO 2 uptake capacity, an analytic model to explain this phenomenon is not well established. Instead, some empirical models such as the Toth model are used to describe the isotherm without accounting for the presence of water. Recently, we have demonstrated that the isotherm equation of CO 2 sorption into strong-base anion exchange materials with quaternary ammonium can be derived from that of strong-base aqueous alkaline solutions by correcting for the drastic change in water activity and by including an appropriate parameterization of the water activity terms. In this paper, we generalize this model from quaternary ammonium to primary, secondary and tertiary amines either in solutions or as functional groups in polymer resins. For primary, secondary and tertiary amines, the isotherm equation can be derived by extending that of a weak-base aqueous alkaline solution such as aqueous ammonia. The model has been validated using experimental data of aqueous ammonia in literature. Furthermore, this general model even includes quaternary ammonium as a special limit. Hence, this general model offers a platform that can treat the isotherms of solid amines, aqueous amines and aqueous alkaline solutions in a unified way.

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

Progress in modelling fast-ion D-alpha spectra and neutral particle analyzer fluxes using FIDASIM

FIDASIM is a code that models signals produced by charge-exchange reactions between neutrals and ions (both fast and thermal) in magnetically confined plasmas. With the ion distribution function as input, the code predicts the efflux to a neutral particle analyzer diagnostic and the photon radiance of Balmer-alpha light to a fast-ion D α diagnostic, in addition to many other related quantities. A new, parallelized version of the Monte Carlo code FIDASIM has been developed in Fortran90 that is substantially faster than the original interactive data language version. Modified algorithms include more accurate treatments of the time dependent collisional-radiative equations that describe neutral energy levels, of the cloud of ‘halo’ neutrals that surround the injected neutral beam, and of finite Larmor radius effects. Enhanced physics capabilities include modelling ‘passive’ signals from cold edge neutrals, the ability to treat general three-dimensional magnetic confinement configurations, and calculations of diagnostic-specific weight functions that enable tomographic reconstructions of the fast-ion distribution function. Neutral beam attenuation, beam emission, and fast-ion birth profiles are also modelled. Finally, the new algorithms have been successfully validated against experimental data and new features have been tested through benchmarks between two independently developed versions of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Lie algebraic theory of barren plateaus for deep parameterized quantum circuits

Variational quantum computing schemes train a loss function by sending an initial state through a parametrized quantum circuit, and measuring the expectation value of some operator. Despite their promise, the trainability of these algorithms is hindered by barren plateaus (BPs) induced by the expressiveness of the circuit, the entanglement of the input data, the locality of the observable, or the presence of noise. Up to this point, these sources of BPs have been regarded as independent. In this work, we present a general Lie algebraic theory that provides an exact expression for the variance of the loss function of sufficiently deep parametrized quantum circuits, even in the presence of certain noise models. Our results allow us to understand under one framework all aforementioned sources of BPs. This theoretical leap resolves a standing conjecture about a connection between loss concentration and the dimension of the Lie algebra of the circuit’s generators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Learning Functions Varying along a Central Subspace

Many functions of interest are in a high-dimensional space but exhibit low-dimensional structures. This paper studies regression of an s-Hölder function in $R^D$ which varies along a central subspace of dimension $d$ while $d \ll D$. A direct approximation of $f$ in $R^D$ with an accuracy $\varepsilon$ requires the number of samples in the order of $\varepsilon^{-(2s+D)/s}$. In this paper, we analyze the generalized contour regression (GCR) algorithm for the estimation of the central subspace and use piecewise polynomials for function approximation. GCR is among the best estimators for the central subspace, but its sample complexity is an open question. In this paper, we partially answer this questions by proving that if a variance quantity is exactly known, GCR leads to a mean squared estimation error of $O(n^{-1})$ for the central subspace. The estimation error of this variance quantity is also given in this paper. The mean squared regression error of $f$ is proved to be in the order of $(n/\log n)^{-\frac{2s}{2s+d}}$, where the exponent depends on the dimension of the central subspace instead of the ambient space . This result demonstrates that GCR is effective in learning the low-dimensional central subspace. We also propose a modified GCR with improved efficiency. Here, the convergence rate is validated through several numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Neutral Pion Electroproduction and development of a Neutral Particle Spectrometer

Protons and neutrons, i.e. nucleons, are the basic building blocks of the matter in the visible universe. The strong force binds the nucleons to form nuclei. The electromagnetic force forms the atoms by binding the electrons with the nuclei. The electromagnetic interaction is well understood by Quantum Electrodynamics (QED), which shows the most precise predictability amongst all the theories in physics. In QED, charges interact with each other by exchanging photons. The Quantum Chromodynamics (QCD) describes the strong interaction. Its degrees of freedom are quarks and gluons, the fundamental constituents of the nucleons. The quarks interact with each other by exchanging gluons. However, unlike QED, the gluons interact amongst themselves. This feature of self bindings of the gluons confines the quarks and gluons in the nucleons/hadrons, never to be seen as free. In order to study some of the features of QCD, such as confinement or the structure of the nucleon, one usually needs to rely on experiments. Electromagnetic probes, governed by the well-understood QED, are excellent tools to probe the nucleon. In general, different scales, e.g. electron beam energies, probe different regions of the nucleon. At low energy, of the order of a few GeV, the electron probes the nucleon in the valence quark region. As its energy increases, the electron probes the sea quark and gluon regions. The study of the nucleon structure in all these regions is needed to fully understand QCD. Form factors and parton distribution functions measured from elastic scattering and deep inelastic scattering of leptons off nucleons have provided a partial view of the internal structure of the nucleon. In the mid-1990s, Generalized Parton Distributions (GPDs) were developed. These new objects are a generalization of the form factors and parton distribution functions, but contain richer information on the nucleon internal structure. GPDs are accessible experimentally by deep exclusive reactions. Deeply virtual Compton scattering (DVCS) and deeply virtual meson production (DVMP) are some examples. The first dedicated DVCS/DVMP experiment took place in 2004 in Hall A at The Thomas Jefferson National Accelerator Facility, i.e. Jefferson Lab, in Virginia, U.S.A. A new DVCS/DVMP experiment, after the beam energy upgrade of Jefferson Lab, was carried out in Hall A in a wider kinematic range. Its data were taken from 2014 to 2016. In Hall C at Jefferson Lab, the next DVCS/DVMP experiment will take place. The Hall C experiment will further exploit the kinematic range with higher precision. A Neutral Particle Spectrometer (NPS) is in development to measure DVCS/DVMP events under high background conditions. Jefferson Lab will provide the highest precision data in the valence quark region for various exclusive reactions. The Electron-Ion Collider (EIC) is a future experimental facility currently planned to start operations around 2030 in the U.S.A. Its high energy and high luminosity will probe the sea quark and gluon regions providing answers to the outstanding questions of QCD, in particular in the region where matter is dominated by gluons. First of all, this document describes the data analysis and results of the Hall A neutral pion electroproduction off the proton, from the data taken in 2014-2016. Later, some of the developments towards the construction of the electromagnetic calorimeter of the NPS for the upcoming DVCS/DVMP experiment in Hall C are presented. Finally, one of the candidate materials for the EIC calorimeter, a glass scintillator, will be briefly introduced. I have participated to all these projects, in collaboration with many colleagues. I present in this thesis my contributions to each of these projects. My contributions to the neutral pion data analysis were focused on background subtractions on the calorimeter, acceptance calculations, and the estimation of the systematic uncertainty associated to the event selection cuts. Some necessary information on calibrations of the detectors and data analysis methods are also described. In the NPS project, I performed background dose calculations and energy and position resolution studies of the calorimeter, all using Monte Carlo simulations, with realistic geometries of the experimental apparatus. Characterization of the crystals of the calorimeter was also done. Additionally, I measured the radiation hardness of some glass scintillator in its early stage of development. In order to have a future reference when the glass calorimeter prototype will be tested, I simulated the energy resolution of the prototype.

Ko, Ho-San↗
DOE OSTI2020

Parton Energy Loss in the Generalized High-Twist Approach

In this work, we calculate the radiative parton energy loss in the deeply inelastic scattering (DIS) off a large nucleus within a generalized high-twist approach. The final gluon radiation spectrum is a convolution of the hard partonic part and the transverse momentum dependent (TMD) quark-gluon correlation function, without the twist expansion in the transverse momentum of the initial gluons used in original high twist approach. The TMD quark-gluon correlation function can be factorized approximately as the product of initial quark distribution and TMD gluon distribution which can be used to define the generalized or TMD jet transport coefficient. The radiation spectrum will recover the Gylassy-Levai-Vitev (GLV) result in the first order of the opacity expansion, under the static scattering center and soft gluon radiation approximation. We also investigate numerically the difference as a result of the soft gluon radiation approximation, under the static scattering center approximation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories

Significance The ground state of a quantum mechanical system is the lowest-energy eigenstate of the Hamiltonian. In isolation, it persists unchanged forever, with symmetries dictated by those of the Hamiltonian. But near-eigenstates of broken symmetry can persist for long times, even on the scale of human measurement. The appearance of broken symmetries of the electron density or spin density in a density functional calculation can reveal strong correlations among the electrons that are present in a symmetry-unbroken wavefunction. Symmetry breaking can arise when a wave-like fluctuation drops to zero frequency. The presented examples are the stretched hydrogen molecule, antiferromagnetism in solids, and the static charge-density wave in a low-density jellium, which is shown quantitatively to be a zero-frequency plasma wave.

42 ENGINEERING↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
DOE OSTI2023

Function-based Taxonomy Implementation

Function-based taxonomies are relational mapping tools used to clearly illustrate how an organization delivers Critical Functions by using various people, processes, technologies, information, and infrastructure (PPTII). Critical Functions are the actions or activities that make up the organization’s primary purpose. Enabling Functions are the combination of PPTII used by the organization to deliver their Critical Functions. A function-based taxonomy provides a framework to categorize information and related artifacts that document an organization’s unique implementation of PPTII for Critical Function delivery. The ideal state of a function-based taxonomy is to organize the existing knowledge the organization already has throughout their numerous systems, policies, people, procedures, and configurations. A well-organized taxonomy helps organizations to effectively leverage their knowledge by clearly identifying dependencies and connections. This document is the result of combined engineering and analytical experience and describes a repeatable method for producing function-based taxonomies.

42 ENGINEERING↗
DOE OSTI2025

First Global Extraction of Generalized Parton Distributions from Experiment and Lattice Data with Next-to-Leading-Order Accuracy

We report the first global extraction of generalized parton distributions, GUMP 1.0, by combining deeply virtual Compton scattering and 𝜌-meson production data from Jefferson Lab and the Hadron-Electron Ring Accelerator with global fits of parton distribution functions, charge form factors, and lattice quantum chromodynamics simulations. Using a conformal moment space parametrization, we achieve a unified description across low- and high-𝑥 regions at next to leading order accuracy in perturbative corrections. The results provide state-of-the-art generalized parton distributions consistent with almost all known facts, enabling three-dimensional nucleon imaging in impact parameter space and, at the same time, establishing a benchmark for future theoretical and experimental studies of the nucleon structure.

Form factors↗