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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Modeling the physisorption of graphene on metals

Many processes of technological and fundamental importance occur on surfaces. Adsorption is one of these phenomena that has received the most attention. However, it presents a great challenge to conventional density functional theory. By starting with the Lifshitz-Zaremba-Kohn second-order perturbation theory, we develop a long-range van der Waals (vdW) correction for physisorption of graphene on metals. The model importantly includes quadrupole-surface interaction and screening effects. The results show that, when the vdW correction is combined with the Perdew-Burke-Enzerhof functional, it yields adsorption energies in good agreement with the random-phase approximation, significantly improving upon other vdW methods. We also find that, compared with the leading-order interaction, the higher-order quadrupole-surface correction accounts for about 25$\%$ of the total vdW correction, suggesting the importance of the higher-order term.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Including screening in van der Waals corrected density functional theory calculations: The case of atoms and small molecules physisorbed on graphene

The Density Functional Theory (DFT)/van der Waals-Quantum Harmonic Oscillator-Wannier function (vdW-QHO-WF) method, recently developed to include the vdW interactions in approximated DFT by combining the quantum harmonic oscillator model with the maximally localized Wannier function technique, is applied to the cases of atoms and small molecules (X=Ar, CO, H{sub 2}, H{sub 2}O) weakly interacting with benzene and with the ideal planar graphene surface. Comparison is also presented with the results obtained by other DFT vdW-corrected schemes, including PBE+D, vdW-DF, vdW-DF2, rVV10, and by the simpler Local Density Approximation (LDA) and semilocal generalized gradient approximation approaches. While for the X-benzene systems all the considered vdW-corrected schemes perform reasonably well, it turns out that an accurate description of the X-graphene interaction requires a proper treatment of many-body contributions and of short-range screening effects, as demonstrated by adopting an improved version of the DFT/vdW-QHO-WF method. We also comment on the widespread attitude of relying on LDA to get a rough description of weakly interacting systems.

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

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗
DOE OSTI2015

Structural, bonding, and electronic properties of the hexagonal ferroelectric and paraelectric phases of LuMnO{sub 3} compound: A density functional theory study

We have investigated the structural, bonding, and electronic properties of both ferroelectric (FE) and paraelectric (PE) phases of the hexagonal LuMnO{sub 3} compound using calculations based on density functional theory. The structural properties have been determined by employing the generalized gradient approximation with Perdew-Burke-Ernzerhof and Wu-Cohen parameterization. The bonding and electronic properties have been treated by recently developed modified Becke-Johnson exchange potential, which succeeded to open a band gap for both PE and FE phases, in agreement with experimental predictions. The Bader’s topological analysis of electronic density showed that the character of the Lu–O axial bonds changes when the crystal exhibits the PE → FE structural transition. This fact is in agreement with experimental findings. The covalent character of the Lu–O bond significantly increases due to orbital hybridization between the Lu 5d{sub z}{sup 2} and O 2p{sub z}-states. This bonding mechanism causes the ferroelectricity in the hexagonal LuMnO{sub 3} compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

System-dependent exchange–correlation functional with exact asymptotic potential and ε{sub HOMO} ≈ − I

Density scaling considerations are used to derive an exchange–correlation explicit density functional that is appropriate for the electron deficient side of the integer and which recovers the exact r → ∞ asymptotic behaviour of the exchange–correlation potential. The functional has an unconventional mathematical form with parameters that are system-dependent; the parameters for an N-electron system are determined in advance from generalised gradient approximation (GGA) calculations on the N- and (N − 1)-electron systems. Compared to GGA results, the functional yields similar exchange–correlation energies, but HOMO energies that are an order of magnitude closer to the negative of the vertical ionisation potential; for anions, the HOMO energies are negative, as required. Rydberg excitation energies are also notably improved and the exchange–correlation potential is visibly lowered towards the near-exact potential. Further development is required to improve valence excitations, static isotropic polarisabilities, and the shape of the potential in non-asymptotic regions. The functional is fundamentally different to conventional approximations.

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

Effect of Functional and Electron Correlation on the Structure and Spectroscopy of the Al 2 O 3 (001)–H 2 O Interface

Oxide–water interfaces are ubiquitous, with many applications in industry and the environment, yet there is a great deal of controversy over their properties and microscopic structure. This controversy stems, in part, from the unique H-bond networks formed at different surface terminations and mineral compositions. Density functional theory simulations of these interfaces require an accurate description of both the oxide mineral and water in diverse H-bond environments. Thus, herein we simulate the Al 2 O 3 (001)–H 2 O interface using the PBE, PBE-TS, RPBE, SCAN, and HSE06-TS functionals to determine how calculated interfacial properties depend on the choice of functional. We find that the structure of the first few layers of water at the surface is determined by electron correlation in a way that cannot be approximated using semiemipirical van der Waals corrections. Of the functionals investigated, we find that SCAN yields the most accurate interfacial structure, dynamics, and sum frequency generation spectrum. Furthermore, SCAN leads to a reduction in the order of the 2D H-bond network of water at the alumina surface predicted by GGA functionals, leading to a significant decrease in the anisotropy of the diffusion coefficient at the surface. We emphasize the importance of using a functional which accurately describes electron correlation for more complex oxides, such as transition-metal oxides, where electron correlation may play an even greater role in determining the structure and dynamics of the oxide–water interface.

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

The Effect of Core-Hole Shape on Attosecond Valence Electron Dynamics

Rapid X-ray ionization of a core electron is known to induce the attosecond motion of valence electrons; however, the effect of core-hole shape on the triggered dynamics remains relatively unknown. In this work, the sub-four fs response of prototypical functionalized/heterocyclic/polycyclic molecules was simulated using real-time time-dependent density functional theory (RT-TDDFT), a sudden approximation core-hole, and phenomenologically Auger–Meitner (AM) decay. These molecules included fluorobenzene, chlorobenzene, bromobenzene, phenol, thiophenol, pyridine, phosphorus, and azulene. It is observed that the valence electron dynamics are essentially independent of the core-hole created, provided that it is ionized from an inner-shell orbital and not an inner-valence orbital. This has broad implications for free-electron laser studies of X-ray pumped attosecond processes since the flexibility in edge allows for a wide range of experimental modalities, core-holes with longer AM lifetimes, and molecular targets.

aromatic compounds↗
DOE OSTI2016

Structural phase transition and electronic properties in samarium chalcogenides

The electronic structure and high pressure properties of samarium monochalcogenides SmS, SmSe and SmTe have been reported by using tight binding linear muffin-tin-orbital (TB-LMTO) method within the local density approximation (LDA). The total energy as a function of volume is evaluated. It is found that these monochalcogenides are stable in NaCl-type structure under ambient pressure. We predict a structural phase transition from NaCl-type (B{sub 1}-phase) structure to CsCl-type (B{sub 2}-type) structure for these compounds. Phase transition pressures were found to be 1.7, 4.4 and 6.6 GPa, for SmS, SmSe and SmTe respectively. Apart from this, the lattice parameter (a{sub 0}), bulk modulus (B{sub 0}), band structure (BS) and density of states (DOS) are calculated. From energy band diagram we observed that these compounds exhibit metallic character. The calculated values of equilibrium lattice parameter and phase transition pressure are in general good agreement with available data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Electronic and Magnetic Investigations of Rare-Earth Tm-doped AlGaN Ternary Alloy

Electronic structure and magnetic interaction of substitutional thulium rare earth-doped wurtzite Al{sub 0.5}Ga{sub 0.5}N ternary alloy have been performed using density-functional theory within local spin-density approximation with Hubbard-U corrections (LSDA+U) approach. The LSDA+U method is applied to the rare earth (RE) 4f states. The calculation of formation energy shows that it is more energetically favorable for a substitutional Tm atom to replace the Al atom than the Ga atom. For AlGaN:Tm, the lattices parameters are expanded due to larger ionic radius of Tm than that of Al atoms. The energy band gap of AlGaN:Tm has direct character and its width becomes small compared with that of AlGaN. The magnetic coupling between Tm ions in the nearest neighbor sites is ferromagnetic. Magnetic interaction of rare earth ion with the host states at the valence and conduction band edges has been investigated and compared to those of GaN:Mn and has been found to be relatively small.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Transition matrices and orbitals from reduced density matrix theory

In this contribution, we report two different methodologies for characterizing the electronic structure reorganization occurring when a chromophore undergoes an electronic transition. For the first method, we start by setting the theoretical background necessary to the reinterpretation through simple tensor analysis of (i) the transition density matrix and (ii) the natural transition orbitals in the scope of reduced density matrix theory. This novel interpretation is made more clear thanks to a short compendium of the one-particle reduced density matrix theory in a Fock space. The formalism is further applied to two different classes of excited states calculation methods, both requiring a single-determinant reference, that express an excited state as a hole-particle mono-excited configurations expansion, to which particle-hole correlation is coupled (time-dependent Hartree-Fock/time-dependent density functional theory) or not (configuration interaction single/Tamm-Dancoff approximation). For the second methodology presented in this paper, we introduce a novel and complementary concept related to electronic transitions with the canonical transition density matrix and the canonical transition orbitals. Their expression actually reflects the electronic cloud polarisation in the orbital space with a decomposition based on the actual contribution of one-particle excitations from occupied canonical orbitals to virtual ones. This approach validates our novel interpretation of the transition density matrix elements in terms of the Euclidean norm of elementary transition vectors in a linear tensor space. A proper use of these new concepts leads to the conclusion that despite the different principles underlying their construction, they provide two equivalent excited states topological analyses. This connexion is evidenced through simple illustrations of (in)organic dyes electronic transitions analysis.

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

Calculation of absorption spectra involving multiple excited states: Approximate methods based on the mixed quantum classical Liouville equation

We investigate the calculation of absorption spectra based on the mixed quantum classical Liouville equation (MQCL) methods. It has been shown previously that, for a single excited state, the averaged classical dynamics approach to calculate the linear and nonlinear spectroscopy can be derived using the MQCL formalism. This work focuses on problems involving multiple coupled excited state surfaces, such as in molecular aggregates and in the cases of coupled electronic states. A new equation of motion to calculate the dipole-dipole correlation functions within the MQCL formalism is first presented. Two approximate methods are then proposed to solve the resulted equations of motion. The first approximation results in a mean field approach, where the nuclear dynamics is governed by averaged forces depending on the instantaneous electronic states. A modification to the mean field approach based on first order moment expansion is also proposed. Numerical examples including calculation of the absorption spectra of Frenkel exciton models of molecular aggregates, and the pyrazine molecule are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Electronic density-of-states of amorphous vanadium pentoxide films: Electrochemical data and density functional theory calculations

Thin films of V{sub 2}O{sub 5} were prepared by sputter deposition onto transparent and electrically conducting substrates and were found to be X-ray amorphous. Their electrochemical density of states was determined by chronopotentiometry and displayed a pronounced low-energy peak followed by an almost featureless contribution at higher energies. These results were compared with density functional theory calculations for amorphous V{sub 2}O{sub 5}. Significant similarities were found between measured data and computations; specifically, the experimental low-energy peak corresponds to a split-off part of the conduction band apparent in the computations. Furthermore, the calculations approximately reproduce the experimental band gap observed in optical measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

The interacting quasiparticle–phonon picture and odd–even nuclei. Overview and perspectives

The role of the nucleon correlations in the ground states of even–even nuclei on the properties of low-lying states in odd–even spherical and transitional nuclei is studied. We reason about this subject using the language of the quasiparticle–phonon model which we extend to take account of the existence of quasiparticle⊗phonon configurations in the wave functions of the ground states of the even–even cores. Of paramount importance to the structure of the low-lying states happens to be the quasiparticle–phonon interaction in the ground states which we evaluated using both the standard and the extended random phase approximations. Numerical calculations for nuclei in the barium and cadmium regions are performed using pairing and quadrupole–quadrupole interaction modes which have the dominant impact on the lowest-lying states’ structure. It is found that states with same angular momentum and parity become closer in energy as compared to the predictions of models disregarding the backward amplitudes, which turns out to be in accord with the experimental data. In addition we found that the interaction between the last quasiparticle and the ground-state phonon admixtures produces configurations which contribute significantly to the magnetic dipolemoment of odd-A nuclei. It also reveals a potential for reproducing their experimental values which proves impossible if this interaction is neglected.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion

We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A near-stationary subspace for ridge approximation

Response surfaces are common surrogates for expensive computer simulations in engineering analysis. However, the cost of fitting an accurate response surface increases exponentially as the number of model inputs increases, which leaves response surface construction intractable for high-dimensional, nonlinear models. Here we describe ridge approximation for fitting response surfaces in several variables. A ridge function is constant along several directions in its domain, so fitting occurs on the coordinates of a low-dimensional subspace of the input space. We review essential theory for ridge approximation – e.g., the best mean-squared approximation and an optimal low-dimensional subspace – and we prove that the gradient-based active subspace is near-stationary for the least-squares problem that defines an optimal subspace. Motivated by the theory, we propose a computational heuristic that uses an estimated active subspace as an initial guess for a ridge approximation fitting problem. We show a simple example where the heuristic fails, which reveals a type of function for which the proposed approach is inappropriate. We then propose a simple alternating heuristic for fitting a ridge function, and we demonstrate the effectiveness of the active subspace initial guess applied to an airfoil model of drag as a function of its 18 shape parameters.

42 ENGINEERING↗
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

Stationary electron velocity distribution function in crossed electric and magnetic fields with collisions

Analytical studies and numerical simulations show that the electron velocity distribution function in a Hall thruster discharge with crossed electric and magnetic fields is not Maxwellian. This is due to the fact that the mean free path between collisions is greater than both the Larmor radius and the characteristic dimensions of the discharge channel. However in numerical models of Hall thrusters, a hydrodynamic approach is often used to describe the electron dynamics, because discharge simulation in a fully kinetic approach requires large computing resources and is time consuming. A more accurate modeling of the electron flow in the hydrodynamic approximation requires taking into account the non-Maxwellian character of the distribution function and finding its moments, an approach that reflects the properties of electrons drifting in crossed electric and magnetic fields better than the commonly used Euler or Navier-Stokes approximations. In the present paper, an expression for the electron velocity distribution function in rarefied spatially homogeneous stationary plasma with crossed electric and magnetic fields and predominance of collisions with heavy particles is derived in the relaxation approximation. The main moments of the distribution function including longitudinal and transversal temperatures, the components of the viscous stress tensor, and of the heat flux vector are calculated. Distinctive features of the hydrodynamic description of electrons with a strongly non-equilibrium distribution function and the prospects for further development of the proposed approach for calculating the distribution function in spatially inhomogeneous plasma are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Half-metallic Ferromagnetism in Novel Rh{sub 2}-based Full-Heusler Alloys Rh{sub 2}FeZ (Z = Ga and In)

In this approach, we have employed the ab initio method of full-potential linearized plane waves plus local orbital based on the framework of density functional theory, within the spin generalized gradient approximation plus U (GGA+U), where U is the Hubbard correlation term for exchange-correlation potential. This method is used to predict the equilibrium structural parameters, the electronic structure, and the magnetic properties of Rh{sub 2}FeGa and Rh{sub 2}FeIn full-Heusler alloys. The structural properties result show that both Rh{sub 2}FeGa and Rh{sub 2}FeIn compounds are stable in ferromagnetic states. The spin-polarized electronic structure demonstrate that both full-Heusler compounds are half-metals with half-metallic band gap (E{sub HM}) of 0.326 eV for Rh{sub 2}FeGa alloy and 0.245 eV for Rh{sub 2}FeIn alloy, where the apparition of an indirect gap from Γ to X in minority spin is shown by these Heusler alloys. Furthermore, calculations of total magnetic moment give the value of 5 µ{sub B}, which agrees with Slater-Pauling rule; also, the atomic magnetic moment of both Rh and Fe transition metals have the same sign which depicts the ferromagnetic behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Refine your search to explore more results.