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

Spaceflight-Induced Cardiovascular Changes and Recovery During NASA's Functional Task Test

Microgravity-induced physiological changes could impair a crewmember s performance upon return to a gravity environment. The Functional Task Test (FTT) is designed to correlate these physiological changes to performance in mission-critical tasks. The Recovery from Fall/Stand Test (RFST) simulates one such task, measuring the ability to recover from a prone position and the cardiovascular response to orthostasis. The purpose of this study was to evaluate spaceflight-induced cardiovascular changes during the FTT. METHODS: Five astronauts participated in the FTT before 10-15 day missions, on landing day (R+0), and one (R+1), six (R+6) and thirty (R+30) days after landing. The RFST consisted of a 2-minute prone rest followed by a 3-minute stand during which heart rate (HR, Holter) and continuous blood pressure (BP, Finometer) were measured. Spectral heart rate variability (HRV) was calculated during the RFST to approximate autonomic function. Statistical analysis was performed with two-factor repeated measures ANOVA. RESULTS: During RFST, HR was higher on R+0 than preflight (p<0.004). This increase in HR persisted on R+1 and R+6 during the stand portion of RFST (p<0.026). BP was well-regulated on all test days. Parasympathetic activity was diminished on R+0 (p=0.035). Sympathovagal balance tended to be affected by spaceflight (main effect, p=0.072), appearing to be slightly elevated during postflight RFST except on R+30. Additionally, analysis of HR during the functional tasks yielded a higher HR on R+0 than preflight during 8 of 11 tasks analyzed, where all tasks had HR return to preflight values by R+30 (p<0.05). CONCLUSION: Spaceflight causes an increase in HR, decrease in parasympathetic activity, and increase in sympathovagal balance, which we confirmed during RFST. These spaceflight-induced changes seen in the RFST, along with the increased postflight HR in most functional tasks, can be used to assess functional performance after short-duration spaceflight.

Arzeno, N. M.↗

Correlations far from equilibrium in charged strongly coupled fluids subjected to a strong magnetic field

Within a holographic model, we calculate the time evolution of 2-point and 1-point correlation functions (of selected operators) within a charged strongly coupled system of many particles. That system is thermalizing from an anisotropic initial charged state far from equilibrium towards equilibrium while subjected to a constant external magnetic field. One main result is that thermalization times for 2-point functions are significantly (approximately three times) larger than those of 1-point functions. Magnetic field and charge amplify this difference, generally increasing thermalization times. However, there is also a competition of scales between charge density, magnetic field, and initial anisotropy, which leads to an array of qualitative changes on the 2- and 1-point functions. There appears to be a strong effect of the medium on 2-point functions at early times, but approximately none at later times. At strong magnetic fields, an apparently universal thermalization time emerges, at which all 2-point functions appear to thermalize regardless of any other scale in the system. Hence, this time scale is referred to as saturation time scale. As extremality is approached in the purely charged case, 2- and 1-point functions appear to equilibrate at infinitely late time. We also compute 2-point functions of charged operators. Our results can be taken to model thermalization in heavy ion collisions, or thermalization in selected condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multivariable quantum signal processing (M-QSP): prophecies of the two-headed oracle

Recent work shows that quantum signal processing (QSP) and its multi-qubit lifted version, quantum singular value transformation (QSVT), unify and improve the presentation of most quantum algorithms. QSP/QSVT characterize the ability, by alternating ansätze, to obliviously transform the singular values of subsystems of unitary matrices by polynomial functions; these algorithms are numerically stable and analytically well-understood. That said, QSP/QSVT require consistent access to a single oracle, saying nothing about computing joint properties of two or more oracles; these can be far cheaper to determine given an ability to pit oracles against one another coherently. This work introduces a corresponding theory of QSP over multiple variables: M-QSP. Surprisingly, despite the non-existence of the fundamental theorem of algebra for multivariable polynomials, there exist necessary and sufficient conditions under which a desired stable multivariable polynomial transformation is possible. Moreover, the classical subroutines used by QSP protocols survive in the multivariable setting for non-obvious reasons, and remain numerically stable and efficient. Up to a well-defined conjecture, we give proof that the family of achievable multivariable transforms is as loosely constrained as could be expected. The unique ability of M-QSP to obliviously approximate joint functions of multiple variables coherently leads to novel speedups incommensurate with those of other quantum algorithms, and provides a bridge from quantum algorithms to algebraic geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reconstructing high-dimensional Hilbert-valued functions via compressed sensing

We present and analyze a novel sparse polynomial technique for approximating high-dimensional Hilbert-valued functions, with application to parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our theoretical framework treats the function approximation problem as a joint sparse recovery problem, where the set of jointly sparse vectors is possibly infinite. To achieve the simultaneous reconstruction of Hilbert-valued functions in both parametric domain and Hilbert space, we propose a novel mixed-norm based `1 regularization method that exploits both energy and sparsity. Our approach requires extensions of concepts such as the restricted isometry and null space properties, allowing us to prove recovery guarantees for sparse Hilbert-valued function reconstruction. We complement the enclosed theory with an algorithm for Hilbert-valued recovery, based on standard forward-backward algorithm, meanwhile establishing its strong convergence in the considered infinite-dimensional setting. Finally, we demonstrate the minimal sample complexity requirements of our approach, relative to other popular methods, with numerical experiments approximating the solutions of high-dimensional parameterized elliptic PDEs.

Dexter, Nick↗

Electronic structure and thermoelectric properties of CoAsSb with post-DFT approaches

Here, we study the electronic structure and thermoelectric properties of recently synthesized CoAsSb. The calculated bandgap becomes more accurate for increasingly-complex electronic structure methods: generalized gradient approximation, hybrid functionals, selfconsistent linearized quasiparticle GW method (LQSGW), and LQSGW with simplified vertex corrections. By equating the bandgaps of each method from a rigid shift of the bands, we evaluate the contributions made to thermoelectric properties beyond the bandgap. In doing so, we evaluate the efficacy of a commonpractice: a rigid shift applied to less-costly electronic structure methods in order to gain some accuracy of the more-costly methods. We find that while the shift made the Seebeck coefficients much closer to one another than from the original electronic structures, there remain differences in the Goldsmid-Sharp (thermoelectric) bandgap between the methods and from the intended electronic bandgap. Additionally, some lasting differences in temperature dependence remain between the methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hydrodynamic instability and coalescence of binary neutron stars

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configuration, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of confirgurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, and spins. The influence of the stiffness of the equation of state is also explored by varying the effective polytopic index. Typically, we find that the radial infall velocity just prior to contact is about 10% of the tangential orbital velocity. Once the stability limit is reached, the final evolution only takes another orbit. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibiltiy of mass transfer form one neutron star to the other. We show that stable mass transfer is unlikely except when the mass of one of the components is very small (M is less than or approximately 0.4 solar mass) and the viscosity is high enough to maintain corotation. Otherwise, either the two stars come into contact or the dynamical instability sets in before a Roche limit can be reached.

Lai, Dong↗

Skin and Proximity Effects in Electrodes and Furnace Shells

A review of two-dimensional (2D) analytical models of skin and proximity effects in large industrial furnaces with three electrodes arranged in an equilateral triangle is given. The models cover three different cases: one electrode only, three electrodes where two are approximated by line currents, and induced shell currents where all electrodes are approximated by line currents. The first two models show how the skin and proximity effects depend on electrode material properties and size, and the distance between the electrodes. The third model shows how the strength of the induced shell currents will depend on electrode position and furnace size. These models are compared to numerical studies including distributed electrodes and shell currents. The analytical models are accurate when induced shell currents can be disregarded. However, strong shell currents may have a significant impact on the current distribution within the electrodes. This electrode-shell proximity effect competes with the electrode-electrode proximity effect. Finally, the 2D models have been compared with three-dimensional (3D) case studies of large industrial furnaces. In 3D, the shell currents are significantly smaller than what are predicted by the 2D models, but they are sufficiently strong to cause a significant correction of the electrode current density.

36 MATERIALS SCIENCE↗

Overview of the High Performance Antiproton (HiPAT) Experiment

The annihilation of matter with antimatter represents the highest energy density of any known reaction, producing 10(exp 8) MJ/g, approximately 10 orders of magnitude more energy per unit mass than chemical based combustion. To take the first step towards using this energy for propulsion applications the NASA MSFC Propulsion Research Center (PRC) has initiated a research activity examining the storage of low energy antiprotons. Storage was identified as a key enabling technology since it builds the experience base necessary to understand the handling of antiprotons for virtually all utilization and high-density storage concepts. To address this need, a device referred to as the High Performance Antiproton Trap (HiPAT) is under development at the NASA MSFC PRC. The HiPAT is an electromagnetic system (Penning-Malmberg design) consisting of a 4 Tesla superconductor, a high voltage confinement electrode system (operation up to 20 KV), and an ultra high vacuum test section (operating in the 10(exp -12) torr range). The system was designed to be portable with an ultimate goal of maintaining 10(exp 12) charged particles with a half-life of 18 days. Currently, this system is being experimentally evaluated using normal matter ions which are cheap to produce and relatively easy to handle. These normal ions provide a good indication of overall trap behavior, with the exception of assessing annihilation losses. The ions are produced external to HiPAT using two hydrogen ion sources, with adjustable beam energy and current. Ion are transported in a beam line and controlled through the use of electrostatic optics. These optics serve to both focus and gate the incoming ions, providing microsecond-timed pulses that are dynamically captured by cycling the HiPAT electric containment field like a 'trap door'. The layout of this system more closely simulates the operations expected at an actual antiproton production facility where 'packets' of antiprotons with pulse widths measured in 100's of nanoseconds could be provided. Initial dynamic capture experiments have been performed with both trap and ton source system functioning at approximately 10% of maximum levels, minimizing the potential for extraneous effects. Dynamic trapping techniques demonstrated the successful capture of millions of hydrogen ions with good agreement with the predicted loading based on the timing sequence, trap electric field, and ion beam current. These techniques will be expanded to examine multiple shot capture or 'stacking' to increase the total number of ions stored within HiPAT.

Martin, James J.↗

Recent Developments in DFTB+, a Software Package for Efficient Atomistic Quantum Mechanical Simulations

DFTB+ is a flexible, open-source software package developed by its community, designed for fast and efficient atomistic quantum mechanical simulations. It employs various methods that approximate density functional theory (DFT), such as density functional-based tight binding (DFTB) and the extended tight binding (xTB) approach allowing simulations of large systems over extended time scales with reasonable accuracy, while being significantly faster than traditional ab initio methods. In recent years, several new extensions of the DFTB method have been developed and implemented in the DFTB+ program package in order to improve the accuracy and generality of the available simulation results. In this paper, we review those enhancements, show several use case examples and discuss the strengths and limitations of its features.

36 MATERIALS SCIENCE↗

Probabilistic partition of unity networks for high–dimensional regression problems

We explore the probabilistic partition of unity network (PPOU-Net) model in the context of high-dimensional regression problems and propose a general framework focusing on adaptive dimensionality reduction. With the proposed framework, the target function is approximated by a mixture of experts model on a low-dimensional manifold, where each cluster is associated with a fixed-degree polynomial. We present a training strategy that leverages the expectation maximization (EM) algorithm. During the training, we alternate between (i) applying gradient descent to update the DNN coefficients; and (ii) using closed-form formulae derived from the EM algorithm to update the mixture of experts model parameters. Under the probabilistic formulation, step (ii) admits the form of embarrassingly paralleliazable weighted least-squares solves. The PPOU-Nets consistently outperform the baseline fully-connected neural networks of comparable sizes in numerical experiments of various data dimensions. Here, we also explore the proposed model in applications of quantum computing, where the PPOU-Nets act as surrogate models for cost landscapes associated with variational quantum circuits.

97 MATHEMATICS AND COMPUTING↗

Improved Architectures and Training Algorithms for Deep Operator Networks

Operator learning techniques have recently emerged as a powerful tool for learning maps between infinite-dimensional Banach spaces. Trained under appropriate constraints, they can also be effective in learning the solution operator of partial differential equations (PDEs) in an entirely self-supervised manner. In this work we analyze the training dynamics of deep operator networks (DeepONets) through the lens of Neural Tangent Kernel theory, and reveal a bias that favors the approximation of functions with larger magnitudes. To correct this bias we propose to adaptively re-weight the importance of each training example, and demonstrate how this procedure can effectively balance the magnitude of back-propagated gradients during training via gradient descent. We also propose a novel network architecture that is more resilient to vanishing gradient pathologies. Taken together, our developments provide new insights into the training of DeepONets and consistently improve their predictive accuracy by a factor of 10-50x, demonstrated in the challenging setting of learning PDE solution operators in the absence of paired input-output observations.

97 MATHEMATICS AND COMPUTING↗

Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method

Fermi orbital descriptors (FODs) play a key role in Fermi-Löwdin orbital self-interaction correction (FLOSIC) calculations used to remove self-interaction from approximate density functionals. Optimal FODs are obtained by minimizing the self-interaction-corrected total energy, and, in this process, identifying initial sets of FODs becomes crucial. Here we propose, implement and test a novel method for automatically initializing FODs, quick-FOD, based on the minimization of an empirical energy expression that involves a Coulomb-like FOD-electron density attraction, an FOD-FOD short-range repulsion, and an exchange-like FOD-density repulsion. Quick-FOD successfully reproduces FOD arrangements in qualitatively good agreement with Lewis theory and with full-fledged FLOSIC calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gradient-based optimization for regression in the functional tensor-train format

Predictive analysis of complex computational models, such as uncertainty quantification (UQ), must often rely on using an existing database of simulation runs. In this paper we consider the task of performing low-multilinear-rank regression on such a database. Specifically we develop and analyze an efficient gradient computation that enables gradient-based optimization procedures, including stochastic gradient descent and quasi-Newton methods, for learning the parameters of a functional tensor-train (FT). We compare our algorithms with 22 other nonparametric and parametric regression methods on 10 real-world data sets and show that for many physical systems, exploiting low-rank structure facilitates efficient construction of surrogate models. Here, we use a number of synthetic functions to build insight into behavior of our algorithms, including the rank adaptation and group-sparsity regularization procedures that we developed to reduce overfitting. Finally we conclude the paper by building a surrogate of a physical model of a propulsion plant on a naval vessel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

We constructed and tested a new set of norm-conserving pseudopotentials, as well as companion Gaussian basis sets, for the actinide series (Ac–Lr) using the Goedecker, Teter and Hutter (GTH) formalism within the generalized gradient approximation (GGA) functional for exchange-correlation of Perdew, Burke and Ernzerhof (PBE). This fills a critical void in the literature for studying the chemistry of 5f–block elements in condensed phase. The accuracy and reliability of the newly parameterized An–GTH pseudopotentials and basis sets, a variety of tests on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium ( [Xe]4f 14 ) and large core ( [Xe]4f 14 5d 10 ) options that successfully reproduced structure and energetics, particularly redox processes. The medium core size set, in particular, reproduced all electron calculations over multiple oxidation states from 0 to VII, whereas the large core set was suitable only for the early series elements and low oxidation states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pareto Optimization of Oligomer Polarizability and Dipole Moment Using a Genetic Algorithm

High-performance electronic components are highly sought after in order to produce increasingly smaller and cheaper electronic devices. Drawing inspiration from inorganic dielectric materials, in which both polarizability and polarization contribute, organic materials can also maximize both. For a large set of small molecules drawn from PubChem, a Pareto-like front appears between the polarizability and dipole moment, indicating the presence of an apparent trade-off between these two properties. We tested this balance in π-conjugated materials by searching for novel conjugated hexamers with simultaneously large polar- izabilities and dipole moments with potential use for dielectric materials. Using a genetic algorithm (GA) screening technique in conjunction with an approximate density functional tight-binding method for property calculations, we were able to efficiently search chemical space for optimal hexamers. Given the scope of chemical space, using the GA technique saves considerable time and resources by speeding up molecular searches compared to a systematic search. Here, we also explored the underlying structure–function relationships, including sequence and monomer properties, that characterize large polarizability and dipole moment regimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Libration of hydroxyl groups in layered aluminum (oxy)hydroxides and other material analogs: insights from inelastic neutron scattering and theory

We analyzed the hydroxyl librational signatures of five structurally related aluminum (oxy)hydroxides, using inelastic neutron scattering (INS) and plane-wave lattice dynamics simulations. A clear trend across these aluminum-containing phases illustrates the relationship between hydrogen bonding, local atomic structure, and the spectral location and profile of the librational bands. The INS spectra have been compared to previous optical spectroscopy and computational studies, highlighting the complementary nature of the INS technique. Taking into account other structurally or chemically related material analogs, we have identified a correlation between a blueshift (to higher energy) of the upper librational band edge and the geometry of the hydrogen bond interactions, mirroring (with opposite correlation) the well-known redshift in the intramolecular O–H stretching energy with increasing hydrogen bond strength. For hydroxyl groups that do not participate in hydrogen bonding effectively, the bending librations occur at lower energies and hybridize with metal–oxygen lattice modes. Standard density functional theory approximations, including dispersion corrections, struggle to correctly predict vibrational frequencies of motions dominated by H but perform well for metal–oxygen modes, allowing us to make detailed mode assignments in several cases, including a demonstration of how layer-to-layer disorder in boehmite hydrogen bond orientations is reflected in the sharp but minor low energy peaks (at ∼70–80 meV) of the INS spectrum.

Wang, Hsiu-Wen [Oak Ridge National Laboratory (ORN↗

A near X-point charge exchange neutral spectroscopy (CENS) system for DIII-D

A 16 channel spectroscopy system has been installed on DIII-D to provide information about the energy distribution of the atomic neutrals using the Doppler shift and broadening of passive Balmer-α emission. Here, the lines of sight are dominantly in the toroidal direction, with tangencies near the lowest point on closed magnetic flux surfaces moving from the lower divertor upward into the confined plasma. This allows the system to provide details of the neutrals as they undergo various atomic physics processes while traveling upward from the lower divertor. The spectrally resolved measurements provide several advantages that complement typical optical filter based measurements of hydrogenic spectral lines. These include more direct measurements of the neutrals that fuel deeper in the plasma and capturing neutral velocity distribution information via the Doppler broadened and shifted line. In addition, the spectral separation of higher energy “thermal” neutrals and bright cold emission from the scrape-off layer allows for a more straightforward calculation of the underlying neutral densities based on the emission because of the dependence of the thermal emission on confined plasma properties that are approximately flux functions and well measured.

Atomic and molecular physics↗