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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Optimal Bandwidth for Multitaper Spectrum Estimation

A systematic method for bandwidth parameter selection is desired for Thomson multitaper spectrum estimation. We give a method for determining the optimal bandwidth based on a mean squared error (MSE) criterion. When the true spectrum has a second-order Taylor series expansion, one can express quadratic local bias as a function of the curvature of the spectrum, which can be estimated by using a simple spline approximation. This is combined with a variance estimate, obtained by jackknifing over individual spectrum estimates, to produce an estimated MSE for the log spectrum estimate for each choice of time-bandwidth product. The bandwidth that minimizes the estimated MSE then gives the desired spectrum estimate. Additionally, the bandwidth obtained using our method is also optimal for cepstrum estimates. We give an example of a damped oscillatory (Lorentzian) process in which the approximate optimal bandwidth can be written as a function of the damping parameter. Furthermore, the true optimal bandwidth agrees well with that given by minimizing estimated the MSE in these examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

The adiabatic semiclassical perturbation theory for vibrationally inelastic scattering. I - Collinear calculations. II - Three-dimensional treatment

A semiclassical approximation to treat vibrationally inelastic scattering is developed. The vibrational basis set used is adiabatic with respect to a reference potential which is chosen to be as close as possible to the true potential and also gives easily obtainable solutions to the vibrational wave equation. The radial wave functions are obtained using the WKB approximation, and the coupled Schroedinger equations are solved by a first-order perturbation method to yield a phase shift matrix which is exponentiated to give the full scattering matrix. Results were obtained for all the cases computed by Secrest and Johnson and by Clark and Dickinson, and the agreement is better than 10% for half of the cross-sections and rarely off by more than a factor of 2.

Cross, R. J., Jr.↗
NASA NTRS1992

Scattering from arbitrarily shaped microstrip patch antennas

The scattering properties of arbitrarily shaped microstrip patch antennas are examined. The electric field integral equation for a current element on a grounded dielectric slab is developed for a rectangular geometry based on Galerkin's technique with subdomain rooftop basis functions. A shape function is introduced that allows a rectangular grid approximation to the arbitrarily shaped patch. The incident field on the patch is expressed as a function of incidence angle theta(i), phi(i). The resulting system of equations is then solved for the unknown current modes on the patch, and the electromagnetic scattering is calculated for a given angle. Comparisons are made with other calculated results as well as with measurements.

Shively, David G.↗
DOE OSTI2022

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

The simplest response theory methods for computing vertical excitation spectra in condensed-phase are configuration interaction with single excitations (CIS) and linear-response time-dependent density functional theory (TDDFT) within the Tamm–Dancoff approximation. In applications to X-ray absorption spectroscopy (XAS), methods like CIS and TDDFT that codify only single excitations into the wave function are prone to catastrophic errors in main-edge and post-edge features whose shapes act as a crucial fingerprint in structural analyses of liquids. Here, we show that these errors manifest primarily due to a lack of orbital relaxation in conventional linear-response theories and that core-ionized (n - 1-electron) references, like those of electron-affinity TDDFT, can eliminate the errors in the spectral profile, even in the highest-energy parts of the post-edge. Crucially, we find that single excitations atop core-ionized references are sufficient to elucidate liquid-phase XAS spectra with semi-quantitative accuracy, opening the door for methods like electron-affinity CIS/TDDFT to be used as efficient alternatives to higher-order wave function approaches.

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

Development of the self-consistent approximation and its application to the problem of magnetoelastic resonance in an inhomogeneous medium

Our previously proposed approximation involving both the first and second terms of the expansion of the vertex function is generalized to the system of two interacting wavefields of different physical nature. A system of self-consistent equations for the matrix Green’s function and matrix vertex function is derived. On the basis of this matrix generalization of the new self-consistent approximation, a theory of magnetoelastic resonance is developed for a ferromagnetic model, where the magnetoelastic coupling parameter ε(x) is inhomogeneous. Equations for magnetoelastic resonance are analyzed for one-dimensional inhomogeneities of the coupling parameter. The diagonal and off-diagonal elements of the matrix Green’s function of the system of coupled spin and elastic waves are calculated with the change in the ratio between the average value ε and rms fluctuation Δε of the coupling parameter between waves from the homogeneous case (ε ≠ 0, Δε = 0) to the extremely randomized case (ε = 0, Δε ≠ 0) at various correlation wavenumbers of inhomogeneities k{sub c}. For the limiting case of infinite correlation radius (k{sub c} = 0), in addition to approximate expressions, exact analytical expressions corresponding to the summation of all diagrams of elements of the matrix Green’s function are obtained. The results calculated for an arbitrary k{sub c} value in the new self-consistent approximation are compared to the results obtained in the standard self-consistent approximation, where only the first term of the expansion of the vertex function is taken into account. It is shown that the new approximation corrects disadvantages of the Green’s functions calculated in the standard approximation such as the dome shape of resonances and bends on the sides of resonance peaks. The appearance of a fine structure of the spectrum in the form of a narrow resonance on the Green’s function of spin waves and a narrow antiresonance on the Green’s function of elastic waves, which was previously predicted in the standard self-consistent approximation, is confirmed. With an increase in the parameter k{sub c}, the Green’s functions calculated in the standard and new approximations approach each other and almost coincide with each other at k{sub c}/k ≥ 0.5. At the same time, the results of this work indicate that the new self-consistent approximation has a certain advantage for studying the problems of stochastic radiophysics in media with long-wavelength inhomogeneities (small k{sub c} values), because it describes both the shape and width of peaks much better than the standard approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

The role of under-determined approximations in engineering and science application

There is currently a great deal of interest in using response surfaces in the optimization of aircraft performance. The objective function and/or constraint equations involved in these optimization problems may come from numerous disciplines such as structures, aerodynamics, environmental engineering, etc. In each of these disciplines, the mathematical complexity of the governing equations usually dictates that numerical results be obtained from large computer programs such as a finite element method program. Thus, when performing optimization studies, response surfaces are a convenient way of transferring information from the various disciplines to the optimization algorithm as opposed to bringing all the sundry computer programs together in a massive computer code. Response surfaces offer another advantage in the optimization of aircraft structures. A characteristic of these types of optimization problems is that evaluation of the objective function and response equations (referred to as a functional evaluation) can be very expensive in a computational sense. Because of the computational expense in obtaining functional evaluations, the present study was undertaken to investigate under-determinined approximations. An under-determined approximation is one in which there are fewer training pairs (pieces of information about a function) than there are undetermined parameters (coefficients or weights) associated with the approximation. Both polynomial approximations and neural net approximations were examined. Three main example problems were investigated: (1) a function of one design variable was considered; (2) a function of two design variables was considered; and (3) a 35 bar truss with 4 design variables was considered.

Carpenter, William C.↗
DOE OSTI2025

Efficient sampling of free energy landscapes with functions in Sobolev spaces

Molecular simulations of biological and physical phenomena generally involve sampling complicated, rough energy landscapes characterized by multiple local minima. In this work, we introduce a new family of methods for advanced sampling that draw inspiration from functional representations used in machine learning and approximation theory. As shown here, such representations are particularly well suited for learning free energies using artificial neural networks. As a system evolves through phase space, the proposed methods gradually build a model for the free energy as a function of one or more collective variables, from both the frequency of visits to distinct states and generalized force estimates corresponding to such states. Implementation of the methods is relatively simple and, more importantly, for the representative examples considered in this work, they provide computational efficiency gains of up to several orders of magnitude over other widely used simulation techniques.

Approximation theory↗
DOE OSTI2019

Effects of Co Doping and Point Defect on the Ferromagnetism of ZnO

The effects of different doping proportions of Co doping and oxygen vacancies (V{sub O}) on the ferromagnetism of zinc oxide (ZnO) are controversial. To solve this problem, we investigated the effects of Co doping and point defects (V{sub O}, V{sub Zn}, H{sub i}, and Zn{sub i}) on the ferromagnetism of ZnO through first-principle calculations by using generalized gradient approximation + U (GGA + U) under density functional theory. Results show that the closer the distance between Co and vacancies (V{sub O} or V{sub Zn}), the lower the formation energy and the more stable the system will be. Co-doped ZnO with V{sub O} or V{sub Zn} exhibits long-range ordered ferromagnetism and Curie temperature above the room temperature. Especially, Co-doped ZnO with V{sub Zn} shows a half-metallic behavior, resulting in 100% conduction hole polarizability, which is highly beneficial for dilute magnetic semiconductors (DMSs) as a hole injection source. The ferromagnetism of Zn{sub 14}CoO{sub 16} is caused by the double exchange effects among the electrons of the O–2p, Co–3d, and Zn–3d orbits mediated by the hole carriers after complexes were formed by the Co doping and Zn vacancy. Similarly, the origin of ferromagnetism in Zn{sub 15}CoO{sub 15} is derived from the double exchange effects among the electrons of O–2p, Co–3d, and Zn–4s states mediated by the electron carriers after complexes were formed by the Co doping and O vacancy. In addition, the result shows that Co-doped ZnO with interstitial H (H{sub i}) or Zn (Zn{sub i}) is unfavorable for ferromagnetism and should be avoided experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

Optimizing the Zeldovich approximation

We have recently learned that the Zeldovich approximation can be successfully used for a far wider range of gravitational instability scenarios than formerly proposed; we study here how to extend this range. In previous work (Coles, Melott and Shandarin 1993, hereafter CMS) we studied the accuracy of several analytic approximations to gravitational clustering in the mildly nonlinear regime. We found that what we called the 'truncated Zeldovich approximation' (TZA) was better than any other (except in one case the ordinary Zeldovich approximation) over a wide range from linear to mildly nonlinear (sigma approximately 3) regimes. TZA was specified by setting Fourier amplitudes equal to zero for all wavenumbers greater than k(sub nl), where k(sub nl) marks the transition to the nonlinear regime. Here, we study the cross correlation of generalized TZA with a group of n-body simulations for three shapes of window function: sharp k-truncation (as in CMS), a tophat in coordinate space, or a Gaussian. We also study the variation in the crosscorrelation as a function of initial truncation scale within each type. We find that k-truncation, which was so much better than other things tried in CMS, is the worst of these three window shapes. We find that a Gaussian window e(exp(-k(exp 2)/2k(exp 2, sub G))) applied to the initial Fourier amplitudes is the best choice. It produces a greatly improved crosscorrelation in those cases which most needed improvement, e.g. those with more small-scale power in the initial conditions. The optimum choice of kG for the Gaussian window is (a somewhat spectrum-dependent) 1 to 1.5 times k(sub nl). Although all three windows produce similar power spectra and density distribution functions after application of the Zeldovich approximation, the agreement of the phases of the Fourier components with the n-body simulation is better for the Gaussian window. We therefore ascribe the success of the best-choice Gaussian window to its superior treatment of phases in the nonlinear regime. We also report on the accuracy of particle positions and velocities produced by TZA.

Melott, Adrian L.↗
DOE OSTI2022

Meta-GGA exchange-correlation free energy density functional to increase the accuracy of warm dense matter simulations

We discuss strategies for thermalization of the ground-state meta-generalized gradient approximation (meta-GGA) exchange-correlation (XC) functionals. A simple but accurate scheme is implemented via universal additive thermal correction to XC using a perturbative-like self-consistent approach. The additive correction with explicit temperature dependence is applied to the ground-state deorbitalized, strongly constrained and appropriately normed (SCAN-L) meta-GGA XC leading to thermal XC functional denoted here as T-SCAN-L. Thermal T-SCAN-L meta-GGA functional shows significant improvement in density functional theory calculation accuracy for warm dense matter by a factor of 3 to 10, achieving unprecedented accuracy of total pressure between a few tenths and 1% when compared to traditional XC functionals, as demonstrated by the comparison to pathintegral Monte Carlo simulations for helium equation of state. Furthermore, the T-SCAN-L calculations of dc conductivity of warm dense aluminum also give better agreement with experiments over other XC functionals such as PBE and SCAN-L.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1973

Short-term frequency stability of the Rb-87 maser.

Measurements of the short-term stability of the Rb-87 maser as a function of the maser power output and the receiver cutoff frequency are reported. The experimental data are compared to theoretical results obtained from an approximate theory. In this theory the transfer function of the maser for thermal noise is derived, and the spectral density of the phase fluctuations is calculated. An analytical expression for the 'Allan variance' is also given. A comparison of the stability of the Rb-87 maser with existing frequency standards shows its superiority for averaging times less than 1 sec.

Tetu, M.↗
DOE OSTI2017

An Empirical Fitting Method for Type Ia Supernova Light Curves: A Case Study of SN 2011fe

We present a new empirical fitting method for the optical light curves of Type Ia supernovae (SNe Ia). We find that a variant broken-power-law function provides a good fit, with the simple assumption that the optical emission is approximately the blackbody emission of the expanding fireball. This function is mathematically analytic and is derived directly from the photospheric velocity evolution. When deriving the function, we assume that both the blackbody temperature and photospheric velocity are constant, but the final function is able to accommodate these changes during the fitting procedure. Applying it to the case study of SN 2011fe gives a surprisingly good fit that can describe the light curves from the first-light time to a few weeks after peak brightness, as well as over a large range of fluxes (∼5 mag, and even ∼7 mag in the g band). Since SNe Ia share similar light-curve shapes, this fitting method has the potential to fit most other SNe Ia and characterize their properties in large statistical samples such as those already gathered and in the near future as new facilities become available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.

Jin, Shi↗
DOE OSTI2019

Crystal-momentum-resolved contributions to high-order harmonic generation in solids

We analytically and numerically investigate the emission of high-order harmonic radiation from model solids by intense few-cycle midinfrared laser pulses. In single-active-electron approximation, we expand the active electron's wave function in a basis of adiabatic Houston states and describe the solid's electronic band structure in terms of an adjustable Kronig-Penney model potential. For high-order harmonic generation (HHG) from MgO crystals, we examine spectra from two-band and converged multiband numerical calculations. We discuss the characteristics of intra- and interband contributions to the HHG spectrum for computations including initial crystal momenta either from the Γ point at the center of the first Brioullin zone (BZ) only or from the entire first BZ. For sufficiently high intensities of the driving laser field, we find relevant contributions to HHG from the entire first BZ. Here, based on numerically calculated spectra, we scrutinize the cutoff harmonic orders as a function of the laser peak intensity and find good qualitative agreement with our analytical saddle-point-approximation predictions and published theoretical data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

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↗
NASA NTRS1994

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.↗
DOE OSTI2021

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↗