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 1,621 records · Page 90

Structure formation in the Deser-Woodard nonlocal gravity model: a reappraisal

In this work, we extend previous analyses of the structure formation in the f (□{sup −1} R ) model of nonlocal gravity proposed by Deser and Woodard (DW), which reproduces the background expansion of ΛCDM with no need of a cosmological constant nor of any dimensional constant beside Newton's one. A previous analysis based on redshift-space distortions (RSD) data concluded that the model was ruled out. In this work we revisit the issue and find that, when recast in a localized model, the DW model is not ruled out and actually gives a better fit to RSD data than ΛCDM. In fact, the DW model presents a suppressed growth of matter perturbations with respect to ΛCDM and a slightly lower value of σ{sub 8}, as favored by observations. We also produce analytical approximations of the two modified gravity functions, i.e. the anisotropic stress η and the relative change of Newton's constant Y , and of f σ{sub 8}( z ) as a function of redshift. Finally, we also show how much the fit depends on initial conditions when these are generalized with respect to a standard matter-dominated era.

79 ASTRONOMY AND ASTROPHYSICS↗

Vorticity and Vertical Motions Diagnosed from Satellite Deep-Layer Temperatures

Spatial fields of satellite-measured deep-layer temperatures are examined in the context of quasigeostrophic theory. It is found that midtropospheric geostrophic vorticity and quasigeostrophic vertical motions can be diagnosed from microwave temperature measurements of only two deep layers. The lower- ( 1000-400 hPa) and upper- (400-50 hPa) layer temperatures are estimated from limb-corrected TIROS-N Microwave Sounding Units (MSU) channel 2 and 3 data, spatial fields of which can be used to estimate the midtropospheric thermal wind and geostrophic vorticity fields. Together with Trenberth's simplification of the quasigeostrophic omega equation, these two quantities can be then used to estimate the geostrophic vorticity advection by the thermal wind, which is related to the quasigeostrophic vertical velocity in the midtroposphere. Critical to the technique is the observation that geostrophic vorticity fields calculated from the channel 3 temperature features are very similar to those calculated from traditional, 'bottom-up' integrated height fields from radiosonde data. This suggests a lack of cyclone-scale height features near the top of the channel 3 weighting function, making the channel 3 cyclone-scale 'thickness' features approximately the same as height features near the bottom of the weighting function. Thus, the MSU data provide observational validation of the LID (level of insignificant dynamics) assumption of Hirshberg and Fritsch.

Spencer, Roy W.↗

Predicting impurity spectral functions using machine learning

The Anderson Impurity Model (AIM) is a canonical model of quantum many-body physics. Here we investigate whether machine learning models, both neural networks (NN) and kernel ridge regression (KRR), can accurately predict the AIM spectral function in all of its regimes, from empty orbital, to mixed valence, to Kondo. To tackle this question, we construct two large spectral databases containing approximately 410 000 and 600 000 spectral functions of the single-channel impurity problem. We show that the NN models can accurately predict the AIM spectral function in all of its regimes, with pointwise mean absolute errors down to 0.003 in normalized units. We find that the trained NN models outperform models based on KRR and enjoy a speedup on the order of 10 5 over traditional AIM solvers. Finally, the required size of the training set of our model can be significantly reduced using farthest point sampling in the AIM parameter space, which is important for generalizing our method to more complicated multichannel impurity problems of relevance to predicting the properties of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Validating first-principles phonon lifetimes via inelastic neutron scattering

Phonon lifetimes are a key component of quasiparticle theories of transport; yet first-principles lifetimes are rarely directly compared with inelastic neutron scattering (INS) results. Existing comparisons show discrepancies even at temperatures where perturbation theory is expected to be reliable. In this paper, we demonstrate that the reciprocal space voxel (q voxel), which is the finite region in reciprocal space required in INS data analysis, must be explicitly accounted for within theory in order to draw a meaningful comparison. Here, we demonstrate accurate predictions of peak widths of the scattering function when accounting for the q voxel in CaF 2 and ThO 2 . Passing this test implies high fidelity of the phonon interactions and the approximations used to compute the Green's function, serving as a critical benchmark of theory and indicating that other material properties should be accurately predicted, which we demonstrate for thermal conductivity.

36 MATERIALS SCIENCE↗

On the Gibbs phenomenon 4: Recovering exponential accuracy in a sub-interval from a Gegenbauer partial sum of a piecewise analytic function

We continue our investigation of overcoming Gibbs phenomenon, i.e., to obtain exponential accuracy at all points (including at the discontinuities themselves), from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function. We show that if we are given the first N Gegenbauer expansion coefficients, based on the Gegenbauer polynomials C(sub k)(sup mu)(x) with the weight function (1 - x(exp 2))(exp mu - 1/2) for any constant mu is greater than or equal to 0, of an L(sub 1) function f(x), we can construct an exponentially convergent approximation to the point values of f(x) in any subinterval in which the function is analytic. The proof covers the cases of Chebyshev or Legendre partial sums, which are most common in applications.

Gottlieb, David↗

The Prediction of Jet Noise Ground Effects Using an Acoustic Analogy and a Tailored Green's Function

An assessment of an acoustic analogy for the mixing noise component of jet noise in the presence of an infinite surface is presented. The reflection of jet noise by the ground changes the distribution of acoustic energy and is characterized by constructive and destructive interference patterns. The equivalent sources are modeled based on the two-point cross- correlation of the turbulent velocity fluctuations and a steady Reynolds-Averaged Navier-Stokes (RANS) solution. Propagation effects, due to reflection by the surface and refaction by the jet shear layer, are taken into account by calculating the vector Green's function of the linearized Euler equations (LEE). The vector Green's function of the LEE is written in relation to Lilley's equation; that is, approximated with matched asymptotic solutions and the Green's function of the convective Helmholtz equation. The Green's function of the convective Helmholtz equation for an infinite flat plane with impedance is the Weyl-van der Pol equation. Predictions are compared with an unheated Mach 0.95 jet produced by a nozzle with an exit diameter of 0.3302 meters. Microphones are placed at various heights and distances from the nozzle exit in the peak jet noise direction above an acoustically hard and an asphalt surface. The predictions are shown to accurately capture jet noise ground effects that are characterized by constructive and destructive interference patterns in the mid- and far-field and capture overall trends in the near-field.

Miller, Steven A. E.↗

Prediction of a Dynamically Stable New Half-Metallic Phase for the BaN and BaC Compounds

The structural, magnetic, elastic, mechanical, and thermodynamic properties of BaC and BaN compounds in different phases were studied using first-principle calculations based on spin-polarized density functional theory within the generalized gradient approximation (GGA-PBEsol) and the modified Becke–Johnson approach (mBJ-GGA-PBEsol) for the exchange-correlation energy and potential. The following phases—rock-salt (NaCl), CsCl, zinc blende (ZB), NiAs- and WZ-type hexagonal, tetragonal (P4/nmm), and orthorombic (Pnma) phases of BaC and BaN compounds—were considered. We obtained that Pnma phase has the lowest energy configuration as a function of the volume for both the BaN and BaC compounds. The ferromagnetic phase is energetically favored with respect to the non-magnetic phase in the BaN and BaC compounds, except for the CsCl phase in the BaC compound. Considering the phonon dynamics of BaN and BaC compounds in the Pnma, NaCl, ZB, and WZ phases, we observed that the BaN and BaC compounds in the Pnma, NaCl, and ZB phases are dynamically stable. The calculated elastic properties for the Pnma, NaCl, and ZB phases show that they are elastically stable. The Pnma phase for the BaN and BaC compounds, which is a new phase was found to be dynamically and elastically stable. The BaN and BaC compounds exhibit half-metallic behavior in the Pnma, NaCl, and ZB phases. The half-metallic and magnetic character found in the BaN and BaC compounds are attributed to the presence of spin-polarized 2p orbitals of the nitrogen and carbon atoms, respectively. We found that BaN and BaC compounds are half-metallic ferromagnets with magnetic moment of 1 μ{sub B} and 2 μ{sub B} per formula unit, respectively. Using the GGA-PBEsol (mBJ-GGA-PBEsol) approach, our calculated half-metallic gaps for BaN and BaC compounds are 0.22 eV (0.54 eV) and 0.32 eV (0.44 eV) in the Pnma phase, 0.23 eV (1.32 eV) and 0.35 eV (1.00 eV) in the NaCl phase, and 0.38 eV (1.54 eV) and 0.50 eV (1.57 eV) in the ZB phase, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Testing Orions Fairing Separation System

Traditional fairing systems are designed to fully encapsulate and protect their payload from the harsh ascent environment including acoustic vibrations, aerodynamic forces and heating. The Orion fairing separation system performs this function and more by also sharing approximately half of the vehicle structural load during ascent. This load-share condition through launch and during jettison allows for a substantial increase in mass to orbit. A series of component-level development tests were completed to evaluate and characterize each component within Orion's unique fairing separation system. Two full-scale separation tests were performed to verify system-level functionality and provide verification data. This paper summarizes the fairing spring, Pyramidal Separation Mechanism and forward seal system component-level development tests, system-level separation tests, and lessons learned.

Martinez, Henry↗

A persistent adjoint method with dynamic time-scaling and an application to mass action kinetics

In this article, we consider an optimization problem where the objective function is evaluated at the fixed-point of a contraction mapping parameterized by a control variable, and optimization takes place over this control variable. Since the derivative of the fixed-point with respect to the parameter can usually not be evaluated exactly, an adjoint dynamical system can be used to estimate gradients. Using this estimation procedure, the optimization algorithm alternates between derivative estimation and an approximate gradient descent step. We analyze a variant of this approach involving dynamic time-scaling, where after each parameter update the adjoint system is iterated until a convergence threshold is passed. Here, we prove that, under certain conditions, the algorithm can find approximate stationary points of the objective function. We demonstrate the approach in the settings of an inverse problem in chemical kinetics, and learning in attractor networks.

97 MATHEMATICS AND COMPUTING↗

Optimized structure and electronic band gap of monolayer GeSe from quantum Monte Carlo methods

Here, we have used highly accurate quantum Monte Carlo methods to determine the chemical structure and electronic band gaps of monolayer GeSe. Two-dimensional (2D) monolayer GeSe has received a great deal of attention due to its unique thermoelectric, electronic, and optoelectronic properties with a wide range of potential applications. Density functional theory (DFT) methods have usually been applied to obtain optical and structural properties of bulk and 2D GeSe. For the monolayer, DFT typically yields a larger band-gap energy than for bulk GeSe but cannot conclusively determine if the monolayer has a direct or indirect gap. Moreover, the DFT-optimized lattice parameters and atomic coordinates for monolayer GeSe depend strongly on the choice of approximation for the exchange-correlation functional, which makes the ideal structure-and its electronic properties-unclear. In order to obtain accurate lattice parameters and atomic coordinates for the monolayer, we use a surrogate Hessian-based parallel line search within diffusion Monte Carlo to fully optimize the GeSe monolayer structure. The DMC-optimized structure is different from those obtained using DFT, as are calculated band gaps. The potential energy surface has a shallow minimum at the optimal structure. This, combined with the sensitivity of the electronic structure to strain, suggests that the optical properties of monolayer GeSe are highly tunable by strain.

36 MATERIALS SCIENCE↗

Partition function corrections

The vibration-rotation interaction coupling in diatomic molecules is analyzed and first-order expressions for the energy levels are deduced. Corrections to the partition function to account for this coupling are derived. The effects of long-range intermolecular forces due to induced dipole-induced dipole interactions, dipole-induced dipole interactions, charge-dipole interactions, and charge-charge interactions are considered. London's quantized oscillator model of the induced dipole-induced dipole interaction is developed. The partition function including effects of such dispersion forces is derived and is shown to lead to the virial equation of state. For a model combining the hard sphere and long-range attractive type potentials, this is found to lead to the van der Waals equation of state. The equation of state near the critical point is presented in reduced form (thermodynamic variables expressed in units of critical point values). Finally, approximate corrections to the electronic partition functions are derived to account for perturbations of excited electronic states induced by neighboring neutral particles, heavy ions, and fast electron collisions.

Source record↗

Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics

The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effective propagation constants in dense random media under effective medium approximation

The effective medium approximation is applied to study the effective propagation constants in a dense random medium. The dyadic Green's function is introduced to establish the effective medium approximation formalism for electromagnetic waves. The multiple scattering equations and the Lippmann-Schwinger equations for the transition operator are obtained in the configuration average form. The dispersion equations of multiple scattering is derived by using a standard method in quantum mechanics. To obtain an expression for the effective propagation constants the matrix elements of the configuration average dyadic transiton operator are calculated in momentum representation. Numerical illustations are carried out to demonstrate the difference in the effective propagation constants between the use of this approximation and the well-known quasi-crystalline approximation. A comparison is made with measured loss tangent in dry snow.

Zhu, P. Y.↗

Nuclear properties for astrophysical and radioactive-ion-beam applications (II)

We tabulate the ground-state odd-proton and odd-neutron spins and parities, proton and neutron pairing gaps, one- and two-neutron separation energies, quantities related to β-delayed one- and two-neutron emission probabilities, average energy and average number of emitted neutrons, β-decay energy release and half-life with respect to Gamow–Teller decay with a phenomenological treatment of first-forbidden decays, one- and two-proton separation energies, and α-decay energy release and half-life for 9318 nuclei ranging from {sup 16}O to {sup 339}136 and extending from the proton drip line to the neutron drip line. This paper is a new and improved version of Atomic Data And Nuclear Data Tables [66 131 (1997)]. The starting point of our present work is the new study (FRDM(2012)) of nuclear ground-state masses and deformations based on the finite-range droplet model and folded-Yukawa single-particle potential published in a previous issue of Atomic Data And Nuclear Data Tables [109–110, 1 (2016)]. The β-delayed neutron-emission probabilities and Gamow–Teller β-decay rates are obtained from a quasi-particle random-phase approximation with single-particle levels and wave functions at the calculated nuclear ground-state shapes as input quantities. A development since 1997 is we now use a Hauser–Feshbach approach to account for (n, γ) competition and treat first-forbidden decay in a phenomenological approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A new modified BFGS method for unconstrained optimization problems

Using chain rule, we propose a modified secant equation to get a more accurate approximation of the second curvature of the objective function. Then, based on this modified secant equation we present a new BFGS method for solving unconstrained optimization problems. The proposed method makes use of both gradient and function values, and utilizes information from two most recent steps, while the usual secant relation uses only the latest step information. Under appropriate conditions, we show that the proposed method is globally convergent without convexity assumption on the objective function. Comparative numerical results show computational efficiency of the proposed method in the sense of the Dolan–Moré performance profiles.

97 MATHEMATICS AND COMPUTING↗

Two dimensional eye tracking: Sampling rate of forcing function

A study was conducted to determine the minimum update rate of a forcing function display required for the operator to approximate the tracking performance obtained on a continuous display. In this study, frequency analysis was used to determine whether there was an associated change in the transfer function characteristics of the operator. It was expected that as the forcing function display update rate was reduced, from 120 to 15 samples per second, the operator's response to the high frequency components of the forcing function would show a decrease in gain, an increase in phase lag, and a decrease in coherence.

Hornseth, J. P.↗