Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 505 records · Page 28

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of over-magnetization of elemental transition metal solids from the SCAN density functional

Current investigations have revealed that the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation exchange-correlation functional significantly over-magnetizes elemental Fe, Co, and Ni solids. For the paradigmatic case, bcc Fe, the error relative to experiment is ≳ 20 % . A comparative analysis of magnetization results from SCAN and its deorbitalized counterpart, SCAN-L, leads to identification of the source of the discrepancy. It is not from the difference between Kohn-Sham (SCAN-L) and generalized Kohn-Sham (SCAN) procedures. The key is the iso-orbital indicator α (the ratio of the local Pauli and Thomas-Fermi kinetic energy densities). Its deorbitalized counterpart α L has more dispersion in both spin channels with respect to magnetization in an approximate region between 0.6 and 1.2 bohrs around an Fe nucleus. The overall effect is that the SCAN switching function evaluated with α L reduces the energetic disadvantage of the down channel with respect to up compared to the original α , which in turn reduces the magnetization. This identifies the origin of the SCAN magnetization error as insensitivity of the SCAN switching function to α values in the approximate range 0.5 ≲ α ≲ 0.8 and oversensitivity for α ≳ 0.8 .

36 MATERIALS SCIENCE↗

Shock-induced metallization of polystyrene along the principal Hugoniot investigated by advanced thermal density functionals

To date, none of the ab-initio molecular dynamics simulations of polystyrene, often used as an ablator material in inertial confinement fusion targets, with the standard ground-state exchange-correlation (XC) functional in density-functional theory can satisfactorily agree with experiments in terms of reflectivity measurements. Here, we use recently developed thermal strongly constrained and appropriately normed Laplacian dependent meta-generalized gradient approximation XC density functional (T-SCAN-L) and thermal hybrid XC density functional (KDT0) to show that the inclusion of thermal and inhomogeneity effects is crucial for accurate prediction of structural evolution and corresponding insulator-metal transition (IMT) during shock compression. Optical reflectivity calculated as an indicator of IMT is in perfect accord with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Time-dependent current-density-functional theory of charge, energy and spin transport and dynamics in nanoscale systems (Final Report)

The focus of this research program is on developing density functional methods for the calculation of non-equilibrium phenomena of the transport type (finite wavelength, low frequency) and of the dynamical type (long wavelength, high frequency). One such method is the time-dependent current density functional theory, in which the evolution of the physical system is mimicked in a fictitious non-interacting system subjected to an effective one-particle potential, which we approximate as a function of the equilibrium density and non-equilibrium current density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimal Power Flow of Radial Networks and Its Variations: A Sequential Convex Optimization Approach

This paper proposes a sequential convex optimization method to solve broader classes of optimal power flow (OPF) problems over radial networks. The non-convex branch power flow equation is decomposed as a second-order cone inequality and a non-convex constraint involving the difference of two convex functions. Provided with an initial solution offered by an inexact second-order cone programming relaxation model, this approach solves a sequence of convexified penalization problems, where concave terms are approximated by linear functions and updated in each iteration. It could recover a feasible power flow solution, which usually appears to be very close, if not equal, to the global optimal one. Two variations of the OPF problem, in which non-cost related objectives are optimized subject to power flow constraints and the convex relaxation is generally inexact, are elaborated in detail. One is the maximum loadability problem, which is formulated as a special OPF problem that seeks the maximal distance to the boundary of power flow insolvability. The proposed method is shown to outperform commercial nonlinear solvers in terms of robustness and efficiency. The other is the bi-objective OPF problem. A non-parametric scalarization model is suggested, and is further reformulated as an extended OPF problem by convexifying the objective function. It provides a single trade-off solution without any subjective preference. The proposed computation framework also helps retrieve the Pareto front of the bi-objective OPF via the e-constraint method or the normal boundary intersection method. This paper also discusses extensions for OPF problems over meshed networks based on the semidefinite programming relaxation method.

bi-objective optimization↗

Vibrational and Thermal Properties of Oxyanionic Crystals

The vibrational and thermal properties of dolomite and alkali chlorates and perchlorates were studied in the gradient approximation of density functional theory using the method of a linear combination of atomic orbitals (LCAO). Long-wave vibration frequencies, IR and Raman spectra, and mode Gruneisen parameters were calculated. Equation-of-state parameters, thermodynamic potentials, entropy, heat capacity, and thermal expansion coefficient were also determined. The thermal expansion coefficient of dolomite was established to be much lower than for chlorates and perchlorates. The temperature dependence of the heat capacity at T > 200 K was shown to be generally governed by intramolecular vibrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Conditional probability density functional theory

Here we present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this paper we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

36 MATERIALS SCIENCE↗

Discrete Hahn polynomials for numerical solution of two-dimensional variable-order fractional Rayleigh–Stokes problem

The main aim of this paper is to find the numerical solutions of 2D Rayleigh–Stokes problem with the variable-order fractional derivatives in the Riemann–Liouville sense. The presented method is based on collocation procedure in combination with the new operational matrix of the variable-order fractional derivatives, in the Caputo sense, for the discrete Hahn polynomials. The main advantage of the proposed method is obtaining a global approximation for spatial and temporal discretizations, and it reduced the problem to an algebraic system, which is easier to solve. Also, the profit of approximating a continuous function by Hahn polynomials is that for computing the coefficients of the expansion, we only have to compute a summation and the calculation of coefficients is exact. The error bound for the approximate solution is estimated. Finally, we evaluate results of the presented method with other numerical methods.

97 MATHEMATICS AND COMPUTING↗

Practical algorithms for multivariate rational approximation

We present two approaches for computing rational approximations to multivariate functions, motivated by their effectiveness as surrogate models for high-energy physics (HEP) applications. Our first approach builds on the Stieltjes process to efficiently and robustly compute the coefficients of the rational approximation. Our second approach is based on an optimization formulation that allows us to include structural constraints on the rational approximation (in particular, constraints demanding the absence of singularities), resulting in a semi-infinite optimization problem that we solve using an outer approximation approach. We present results for synthetic and real-life HEP data, and we compare the approximation quality of our approaches with that of traditional polynomial approximations.

97 MATHEMATICS AND COMPUTING↗

First-principles study of electronic states in LiBe{sub 2}

By coupling structure prediction methods with first-principles total energy calculations we find lattice constant and bulk modulus of hcp LiBe{sub 2}. Structural and electronic properties are investigated using Linear Combination of Atomic Orbital’s method. The generalized gradient approximation based density functional calculations are attempted to look for the ground state. Thereafter ground state electronic band structure, total density of states and anisotropies in directional Compton profiles are computed. Step like behavior seen broadly in density of states (DOS) is appreciated by band structure suggesting 2D features. Anisotropic behavior of electronic states around Fermi level is also reflected in anisotropies of directional Compton profiles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Survival of the most transferable at the top of Jacob’s ladder: Defining and testing the ω B97M(2) double hybrid density functional

A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resu|ting functional, ωB97M(2) is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 datq points clearly surpasses that of all of the tested density funstionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wavenumber shift due to nonlinear plasma and wave interaction

Wavenumber shift of the ion Bernstein wave has been observed in the particle-in-cell simulations when the input power of the injected wave is sufficiently large. It is demonstrated that the increase of the total kinetic energy of ions, including both the thermal energy related to the random thermal motion and the oscillation energy due to the coherent motion with the wave, gives rise to such change of the wavenumber. However, the velocity distribution function of the ions can approximately be fitted as a Maxwellian distribution function, and thus, the linear dispersion relation still holds, provided that the initial ion temperature is replaced by the effective temperature measured in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗

Thermodynamics of Anharmonic Systems: Uncoupled Mode Approximations for Molecules

The partition functions, heat capacities, entropies, and enthalpies of selected molecules were calculated using uncoupled mode (UM) approximations, where the full-dimensional potential energy surface for internal motions was modeled as a sum of independent one-dimensional potentials for each mode. The computational cost of such approaches scales the same with molecular size as standard harmonic oscillator vibrational analysis using harmonic frequencies (HO hf ). To compute thermodynamic properties, a computational protocol for obtaining the energy levels of each mode was established. The accuracy of the UM approximation depends strongly on how the one-dimensional potentials of each modes are defined. If the potentials are determined by the energy as a function of displacement along each normal mode (UM-N), the accuracies of the calculated thermodynamic properties are not significantly improved versus the HO hf model. Significant improvements can be achieved by constructing potentials for internal rotations and vibrations using the energy surfaces along the torsional coordinates and the remaining vibrational normal modes, respectively (UM-VT). For hydrogen peroxide and its isotopologs at 300 K, UM-VT captures more than 70% of the partition functions on average. By con trast, the HO hf model and UM-N can capture no more than 50%. For a selected test set of C2 to C8 linear and branched alkanes and species with different moieties, the enthalpies calculated using the HO hf model, UM-N, and UM-VT are all quite accurate comparing with reference values though the RMS errors of the HO model and UM-N are slightly higher than UM-VT. However, the accuracies in entropy calculations differ significantly between these three models. For the same test set, the RMS error of the standard entropies calculated by UM-VT is 2.18 cal mol -1 K -1 at 1000 K. By contrast, the RMS error obtained using the HO model and UM-N are 6.42 and 5.73 cal mol -1 K -1 , respectively. For a test set composed of nine alkanes ranging from C5 to C8, the heat capacities calculated with the UM-VT model agree with the experimental values to within a RMS error of 0.78 cal mol -1 K -1 , which is less than one-third of the RMS error of the HO hf (2.69 cal mol -1 K -1) and UM-N (2.41 cal mol -1 K -1) models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗

Construction and Investigation of the Third-Order Approximation to the Solution of the Heat-Conduction Equation for Thin Shallow Shells by Using Legendre Polynomials in the Case of Stationary Heat Exchange

By using the N th order approximation of temperature function and its first derivative by Legendre polynomials, we construct the solution of problem of heat conduction for a thin-walled shallow isotropic shell and deduce the system of resolving equations for the N th approximation. For the first and third approximations, we solve this problem for the case of concentrated heat sources. In the third approximation, we also construct the plots of the dependences of temperature on the distance to the heat source and on the curvature of the shell under the conditions of symmetric or asymmetric stationary heat exchange.

42 ENGINEERING↗