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.
Approximating and incorporating model uncertainty in an inversion for seismic source functions: Preliminary results
We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.
A Multi-Function AAA Algorithm Applied to Frequency-Dependent Line Modeling: Preprint
We present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages respect to the popular Vector Fitting (VF) algorithm using examples of approximations of admittance matrices discussed in a very recent paper and by the creators of VF. We show that not only MIAAA outperforms VF in terms of achieving better accuracy using a smaller number of poles, but also has no numerical issues to achieve convergence. In contrast to VF, MIAA is not sensitive to the location of input sample points and it does not require good estimates for the location of the desired approximation poles.
Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface
In this work, we explore the suitability of several density functionals with the generalized gradient approximation (GGA) and beyond for describing the dissociative chemisorption of methane on the reconstructed Pt(110)-(2×1) surface. The bulk and surface structures of the metal, methane adsorption energy, and dissociation barrier are used to assess the functionals. A van der Waals corrected GGA functional (optPBE-vdW) and a metaGGA functional with van der Waals correction (MS PBEl-rVV10) are selected for ab initio molecular dynamics calculations of the sticking probability. Our results suggest that the use of these two functionals may lead to a better agreement with existing experimental results, thus serving as a good starting point for future development of reliable machine-learned potential energy surfaces for the dissociation of methane on the Pt(110)-(2×1) surface.
Methane dissociation on Pt(111): Searching for a specific reaction parameter density functional
The theoretical description of methane dissociating on metal surfaces is a current frontier in the field of gas-surface dynamics. Dynamical models that aim at achieving a highly accurate description of this reaction rely on potential energy surfaces based on density functional theory calculations at the generalized gradient approximation. We focus here on the effect that the exchange-correlation functional has on the reactivity of methane on a metal surface, using CHD{sub 3} + Pt(111) as a test case. We present new ab initio molecular dynamics calculations performed with various density functionals, looking also at functionals that account for the van der Waals (vdW) interaction. While searching for a semi-empirical specific reaction parameter density functional for this system, we find that the use of a weighted average of the PBE and the RPBE exchange functionals together with a vdW-corrected correlation functional leads to an improved agreement with quantum state-resolved experimental data for the sticking probability, compared to previous PBE calculations. With this semi-empirical density functional, we have also investigated the surface temperature dependence of the methane dissociation reaction and the influence of the rotational alignment on the reactivity, and compared our results with experiments.
Status of free-energy representations for the homogeneous electron gas
Renewed interest in the homogeneous electron gas (HEG) has been stimulated by recent accurate simulations of it over a wide range of densities and temperatures. Those data, combined with known theoretical limits, have led to analytical representations of the free energy. Such a representation is, at least in principle, the complete HEG equation of state. The initial objective here is to establish that the two best representations of the simulation data and constraints are effectively the same in both functional form and accuracy of representation. The second objective is to disclose and delineate a significant difficulty. Despite their expected accuracy for the free energy, the underlying functional form is not adequate for derived thermodynamic properties. As an example, the specific heats obtained from the representations exhibit anomalies suggesting the need first for additional simulation data in critical regimes, then for refined fitting functions. Here, the existing representations are, however, almost certainly adequate for applications based on the free energy alone (e.g., density functional theory for warm dense matter). The third objective is to show that, despite their inability to provide a complete thermodynamic description of the HEG, the best analytical representations do provide a fully adequate exchange-correlation local density approximation for free energy density functional calculations.
Low-rank separated representation surrogates of high-dimensional stochastic functions: Application in Bayesian inference
This study introduces a non-intrusive approach in the context of low-rank separated representation to construct a surrogate of high-dimensional stochastic functions, e.g., PDEs/ODEs, in order to decrease the computational cost of Markov Chain Monte Carlo simulations in Bayesian inference. The surrogate model is constructed via a regularized alternative least-square regression with Tikhonov regularization using a roughening matrix computing the gradient of the solution, in conjunction with a perturbation-based error indicator to detect optimal model complexities. The model approximates a vector of a continuous solution at discrete values of a physical variable. The required number of random realizations to achieve a successful approximation linearly depends on the function dimensionality. The computational cost of the model construction is quadratic in the number of random inputs, which potentially tackles the curse of dimensionality in high-dimensional stochastic functions. Furthermore, this vector-valued separated representation-based model, in comparison to the available scalar-valued case, leads to a significant reduction in the cost of approximation by an order of magnitude equal to the vector size. The performance of the method is studied through its application to three numerical examples including a 41-dimensional elliptic PDE and a 21-dimensional cavity flow.
First-principles investigation of uranium mononitride (UN): Effect of magnetic ordering, spin-orbit interactions and exchange correlation functional
Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys. Various properties of UN have been previously studied using different density functional theory (DFT) methodologies; however, there are still inconsistencies when it comes to the dynamical stability and defect properties of UN. We address these inconsistencies by studying the UN phonons and defect properties using DFT calculations employing two generalized gradient approximation (GGA) exchange-correlation functionals: PBE and AM05, with and without an added on-site Coulomb repulsion term (+U). Furthermore, we investigate the importance of spin-orbit coupling (SOC) when calculating the properties of UN. We use the different methodologies to determine the preference of UN to have antiferromagnetic (AFM) ordering, as seen in experiments, or ferromagnetic (FM) ordering of the uranium spins. We compare the crystallographic properties, density of states, the DFT X-ray photoelectron spectra and phonon dispersions calculated using the different methodologies. We demonstrate that GGA + U reproduces the AFM ordering in UN, but the crystal structure is dynamically unstable. We also show that magnetic ordering is important in finding the lowest energy defective structure, and that SOC has a distinct influence on the energy of the different uranium interstitial defects. Lastly, we discuss the point defect formation energies under U-rich and N-rich conditions, and the stoichiometric formation energies calculated with the different methodologies, providing insight into the observed tendency for forming hypostoichiometric UN.
Introducing constricted variational density functional theory in its relaxed self-consistent formulation (RSCF-CV-DFT) as an alternative to adiabatic time dependent density functional theory for studies of charge transfer transitions
We have applied the relaxed and self-consistent extension of constricted variational density functional theory (RSCF-CV-DFT) for the calculation of the lowest charge transfer transitions in the molecular complex X-TCNE between X = benzene and TCNE = tetracyanoethylene. Use was made of functionals with a fixed fraction (α) of Hartree-Fock exchange ranging from α = 0 to α = 0.5 as well as functionals with a long range correction (LC) that introduces Hartree-Fock exchange for longer inter-electronic distances. A detailed comparison and analysis is given for each functional between the performance of RSCF-CV-DFT and adiabatic time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation. It is shown that in this particular case, all functionals afford the same reasonable agreement with experiment for RSCF-CV-DFT whereas only the LC-functionals afford a fair agreement with experiment using TDDFT. We have in addition calculated the CT transition energy for X-TCNE with X = toluene, o-xylene, and naphthalene employing the same functionals as for X = benzene. It is shown that the calculated charge transfer excitation energies are in as good agreement with experiment as those obtained from highly optimized LC-functionals using adiabatic TDDFT. We finally discuss the relation between the optimization of length separation parameters and orbital relaxation in the RSCF-CV-DFT scheme.
Ab-initio study of electronic structure and elastic properties of ZrC
The electronic and elastic properties of ZrC have been investigated using the linear combination of atomic orbitals method within the framework of density functional theory. Different exchange-correlation functionals are taken into account within generalized gradient approximation. We have computed energy bands, density of states, elastic constants, bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, lattice parameters and pressure derivative of the bulk modulus by calculating ground state energy of the rock salt structure type ZrC.
Brownian motion of a plasma crystal
The dynamics of a plasma crystal under the action of random external forces is considered. The pair correlation functions of the particle displacements are calculated in the harmonic approximation by using the Langevin equations. The case of a planar hexagonal lattice is analyzed in more detail. Analogues of the Van Hove singularities in the spectral densities are discovered.
Computationally predicted energies and properties of defects in GaN
Recent developments in theoretical techniques have significantly improved the predictive power of density-functional-based calculations. In this article, we discuss how such advancements have enabled improved understanding of native point defects in GaN. We review the methodologies for the calculation of point defects, and discuss how techniques for overcoming the band-gap problem of density functional theory affect native defect calculations. In particular, we examine to what extent calculations performed with semilocal functionals (such as the generalized gradient approximation), combined with correction schemes, can produce accurate results. The properties of vacancy, interstitial, and antisite defects in GaN are described, as well as their interaction with common impurities. We also connect the first-principles results to experimental observations, and discuss how native defects and their complexes impact the performance of nitride devices. Overall, we find that lower-cost functionals, such as the generalized gradient approximation, combined with band-edge correction schemes can produce results that are qualitatively correct. However, important physics may be missed in some important cases, particularly for optical transitions and when carrier localization occurs.
A wave based method to predict the absorption, reflection and transmission coefficient of two-dimensional rigid frame porous structures with periodic inclusions
This paper presents an extension to the Wave Based Method to predict the absorption, reflection and transmission coefficients of a porous material with an embedded periodic set of inclusions. The porous unit cell is described using the Multi-Level methodology and by embedding Bloch–Floquet periodicity conditions in the weighted residual scheme. The dynamic pressure field in the semi-infinite acoustic domains is approximated using a novel wave function set that fulfils the Helmholtz equation, the Bloch–Floquet periodicity conditions and the Sommerfeld radiation condition. The method is meshless and computationally efficient, which makes it well suited for optimisation studies.
Excitation energies of molecules within time-independent density functional theory
Recently proposed exchange energy functional for excited-states is tested for obtaining excitation energies of diatomic molecules. The functional is the ground-state counterpart of the local-density approximation, the modified local spin density (MLSD). The MLSD functional is tested for the N{sub 2} and CO diatomic molecules. The excitation energy obtained with the MLSD functional for the N{sub 2} molecule is in close vicinity to that obtained from the exact exchange orbital functional, Krieger, Li and Iafrate (KLI). For the CO molecule, the departure in excitation energy is observed and is due to the overcorrection of self-interaction.
Effect of mechanical strain on electronic properties of bulk MoS{sub 2}
Ab-initio density functional theory based calculations of electronic properties of bulk and monolayer Molybdenum di-Sulfide (MoS{sub 2}) have been performed using all electron Full Potential Linearised Augmentad Plane Wave (FPLAPW) method using Elk code. We have used Generalised Gradient Approximation (GGA) for exchange and correlation functionals and performed calculaitons of Lattice parameters, Density Of States (DOS) and Band Structure (BS). Band structure calculations revealed that bulk MoS{sub 2} has indirect band gap of 0.97 eV and mono-layer MoS{sub 2} has direct band gap which has increased to 1.71 eV. These are in better agreement with experimental values as compared with the other calculations using pseudo-potential code. The effect of mechanical strain on the electronic properties of bulk MoS{sub 2} has also been studied. For the different values of compressive strain (varying from 2% to 8% in steps of 2%) along the c-axis, the corresponding DOS and BS are obtained. We observed that the band gap decreases by about 15% for every 2% increase in strain along the c-axis.
Electron scattering on A = 3 nuclei from quantum Monte Carlo based approaches
We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
Electron scattering on ${\mathbf{A=3}}$ nuclei from quantum Monte Carlo based approaches
We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
Predicting the thermodynamic stability of double-perovskite halides from density functional theory
Recently, a series of double-perovskite halide compounds such as Cs 2 AgBiCl 6 and Cs 2 AgBiBr 6 have attracted intensive interest as promising alternatives to the solar absorber material CH 3 NH 3 PbI 3 because they are Pb-free and may exhibit enhanced stability. The thermodynamic stability of a number of double-perovskite halides has been predicted based on density functional theory (DFT) calculations of compound formation energies. In this paper, we found that the stability prediction can be dependent on the approximations used for the exchange-correlation functionals, e.g., the DFT calculations using the widely used Perdew, Burke, Ernzerhof (PBE) functional predict that Cs 2 AgBiBr 6 is thermodynamically unstable against phase-separation into the competing phases such as AgBr, Cs 2 AgBr 3 , Cs 3 Bi 2 Br 9 , etc., obviously inconsistent with the good stability observed experimentally. The incorrect prediction by the PBE calculation results from its failure to predict the correct ground-state structures of AgBr, AgCl, and CsCl. By contrast, the DFT calculations based on local density approximation, optB86b-vdW, and optB88-vdW functionals predict the ground-state structures of these binary halides correctly. Furthermore, the optB88-vdW functional is found to give the most accurate description of the lattice constants of the double-perovskite halides and their competing phases. Given these two aspects, we suggest that the optB88-vdW functional should be used for predicting thermodynamic stability in the future high-throughput computational material design or the construction of the Materials Genome database for new double-perovskite halides. As a result, using different exchange-correlation functionals has little influence on the dispersion of the conduction and the valence bands near the electronic bandgap; however, the calculated bandgap can be affected indirectly by the optimized lattice constant, which varies for different functionals.