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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Excitation energies with linear response density matrix functional theory along the dissociation coordinate of an electron-pair bond in N-electron systems

Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ω{sub α} and oscillator strengths f{sub α} for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ω{sub α}(R) curves along the bond dissociation coordinate R for the molecules LiH, Li{sub 2}, and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Dynamically consistent method for mixed quantum-classical simulations: A semiclassical approach

We introduce a new semiclassical (SC) framework, the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), that can be tuned to reproduce existing quantum-limit and classical-limit SC approximations to quantum real-time correlation functions. Applying a modified Filinov transformation to a quantum-limit SC formulation leads to the association of a Filinov parameter with each degree of freedom in the system; varying this parameter from zero to infinity controls the extent of quantization of the corresponding mode. The resulting MQC-IVR expression provides a consistent dynamic framework for mixed quantum-classical simulations and we demonstrate its numerical accuracy in the calculation of real-time correlation functions for a model 1D system and a model 2D system over the full range of quantum- to classical-limit behaviors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

First-principles study of electronic and elastic properties of LuAl{sub 3}

A systematic theoretical study of electronic structure of rare earth intermetallic LuAl{sub 3} has been carried out using full potential linearized augmented plane wave (FP-LAPW) method based on density functional theory (DFT) within the generalized gradient approximation(GGA) for exchange and correlation potential. The ground state properties such as lattice constant (a{sub o}), bulk modulus (B) and pressure derivative of bulk modulus (B′) were evaluated. LuAl{sub 3} has the cubic AuCu{sub 3} type crystal structure. The electronic properties of this compound have been analyzed quantatively from band structure and DOS. It is clear from band structure that this compound is metallic in nature. The calculated elastic constants infer that this compound is mechanically stable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Thermodynamic and mechanical properties of TiC from ab initio calculation

The temperature-dependent thermodynamic and mechanical properties of TiC are systematically investigated by means of a combination of density-functional theory, quasi-harmonic approximation, and thermal electronic excitation. It is found that the quasi-harmonic Debye model should be pertinent to reflect thermodynamic properties of TiC, and the elastic properties of TiC decease almost linearly with the increase of temperature. Calculations also reveal that TiC possesses a pronounced directional pseudogap across the Fermi level, mainly due to the strong hybridization of Ti 3d and C 2p states. Moreover, the strong covalent bonding of TiC would be enhanced (reduced) with the decrease (increase) of temperature, while the change of volume (temperature) should have negligible effect on density of states at the Fermi level. The calculated results agree well with experimental observations in the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Influence of anharmonic phonon decay on self-heating in Si nanowire transistors

Anharmonic phonon-phonon scattering is incorporated into an electro-thermal quantum transport approach based on the nonequilibrium Green's function formalism. Electron-phonon and phonon-phonon interactions are taken into account through scattering self-energies solved in the self-consistent Born approximation. While studying self-heating effects in ultra-scaled Si nanowire transistors, it is found that the phonon decay process softens the artificial accumulation of high energy phonons caused by electron relaxations close to the drain region. This leads to an increase of the device current in the ON-state and a reduction of the effective lattice temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Spin-orbit coupling and the static polarizability of single-wall carbon nanotubes

We calculate the static longitudinal polarizability of single-wall carbon tubes in the long wavelength limit taking into account spin-orbit effects. We use a four-orbital orthogonal tight-binding formalism to describe the electronic states and the random phase approximation to calculate the dielectric function. We study the role of both the Rashba as well as the intrinsic spin-orbit interactions on the longitudinal dielectric response, i.e., when the probing electric field is parallel to the nanotube axis. The spin-orbit interaction modifies the nanotube electronic band dispersions, which may especially result in a small gap opening in otherwise metallic tubes. The bandgap size and state features, the result of competition between Rashba and intrinsic spin-orbit interactions, result in drastic changes in the longitudinal static polarizability of the system. We discuss results for different nanotube types and the dependence on nanotube radius and spin-orbit couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Function space requirements for the single-electron functions within the multiparticle Schroedinger equation

Our previously described method to approximate the many-electron wavefunction in the multiparticle Schroedinger equation reduces this problem to operations on many single-electron functions. In this work, we analyze these operations to determine which function spaces are appropriate for various intermediate functions in order to bound the output. This knowledge then allows us to choose the function spaces in which to control the error of a numerical method for single-electron functions. We find that an efficient choice is to maintain the single-electron functions in L{sup 2} Intersection L{sup 4}, the product of these functions in L{sup 1} Intersection L{sup 2}, the Poisson kernel applied to the product in L{sup 4}, a function times the Poisson kernel applied to the product in L{sup 2}, and the nuclear potential times a function in L{sup 4/3}. Due to the integral operator formulation, we do not require differentiability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Electronic structure and thermoelectric properties of half-Heusler compounds with eight electron valence count—KScX (X = C and Ge)

Electronic band structure and structural properties of two representative half-Heusler (HH) compounds with 8 electron valence count (VC), KScC and KScGe, have been studied using first principles methods within density functional theory and generalized gradient approximation. These systems differ from the well studied class of HH compounds like ZrNiSn and ZrCoSb which have VC = 18 because of the absence of d electrons of the transition metal atoms Ni and Co. Electronic transport properties such as Seebeck coefficient (S), electrical conductivity (σ), electronic thermal conductivity (κ{sub e}) (the latter two scaled by electronic relaxation time), and the power factor (S{sup 2}σ) have been calculated using semi-classical Boltzmann transport theory within constant relaxation time approximation. Both the compounds are direct band gap semiconductors with band extrema at the X point. Their electronic structures show a mixture of heavy and light bands near the valance band maximum and highly anisotropic conduction and valence bands near the band extrema, desirable features of good thermoelectric. Optimal p- or n-type doping concentrations have been estimated based on thermopower and maximum power factors. The optimum room temperature values of S are ∼1.5 times larger than that of the best room temperature thermoelectric Bi{sub 2}Te{sub 3}. We also discuss the impact of the band structure on deviations from Weidemann-Franz law as one tunes the chemical potential across the band gap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Formation of the dopant-oxygen vacancy complexes and its influence on the photoluminescence emissions in Gd-doped HfO{sub 2}

Rare earth doping is widely used to improve the desired properties of high-k dielectric oxides. However, whether rare earth doping can suppress the formation of oxygen vacancies is still debated. By using the first-principles calculations with the generalized gradient approximation and more advanced hybrid functional, we have investigated the structural and electronic properties of the dopant-oxygen vacancy complexes in Gd-doped HfO{sub 2}. Our calculations indicate that the Gd dopants interacting with oxygen vacancies can substantially shift up the V{sub O} energy states towards the conduction band edge. This together with other effects, such as capturing the localized electrons at the oxygen vacancy by Gd dopants and suppressing the randomicity of oxygen vacancy formation, improves the reliability of the devices made from Gd-doped HfO{sub 2}. Based on our calculated results, we have presented an explanation for the experimentally observed decrease of the V{sub O}-related photoluminescence intensities upon Gd doping in HfO{sub 2}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

First-Principles Calculations to Investigate the Structural, Electronic, and Half-Metallic Properties of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} (x = 0, 0.25, 0.5, 0.75, and 1) Quaternary Heusler Alloys

The structural, electronic, and magnetic properties of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} (x = 0, 0.25, 0.5, 0.75, and 1) quaternary Heusler alloys have been investigated using first-principles density functional theory within generalized gradient approximation (GGA) and local spin-density approximation with Hubbard-U corrections (LSDA+U). The calculations have been performed using the full-potential linearized augmented Plane wave plus local orbitals method. The crystal structures of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} were found to be cubic in (F-43m) symmetry for x = 0, 1 and cubic in (P-43m) symmetry for x = 0.25, 0.75, while tetragonal in (P-4m2) symmetry for x = 0.5. The results show that all Ti{sub 2}RhSn1-xSix, Ti{sub 2}RhGe{sub 0.75}Si{sub 0.25}, and Ti{sub 2}RhGe alloys are half-metallic ferrimagnetic. These new materials are good candidates for potential applications in spintronic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: Using density functional theory to accelerate wavefunction methods

The development and implementation of increasingly accurate methods for electronic structure calculations mean that, for many atomistic simulation problems, treating light nuclei as classical particles is now one of the most serious approximations. Even though recent developments have significantly reduced the overhead for modeling the quantum nature of the nuclei, the cost is still prohibitive when combined with advanced electronic structure methods. Here we present how multiple time step integrators can be combined with ring-polymer contraction techniques (effectively, multiple time stepping in imaginary time) to reduce virtually to zero the overhead of modelling nuclear quantum effects, while describing inter-atomic forces at high levels of electronic structure theory. This is demonstrated for a combination of MP2 and semi-local DFT applied to the Zundel cation. The approach can be seamlessly combined with other methods to reduce the computational cost of path integral calculations, such as high-order factorizations of the Boltzmann operator or generalized Langevin equation thermostats.

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

Magneto-Electronic, Thermodynamic, and Thermoelectric Properties of 5f-Electron System BaBkO{sub 3}

In the present work, the structural, electronic, magnetic, mechanical, thermodynamic, and thermoelectric properties of cubic perovskite BaBkO{sub 3} were investigated within the framework of density functional theory using generalized gradient approximation, onsite Coulomb interaction, and modified Becke-Johnson methods. BaBkO{sub 3} is stable in the ferromagnetic phase with a magnetic moment of 7 μ{sub B}. The spin-polarized electronic band structure reveals 100% spin polarization of the material with metallic character in spin-up state and semiconducting in spin-down state. Pugh and Cauchy relation displays that the material is highly ductile. The quasi-harmonic Debye model was used to determine the thermodynamic properties of BaBkO{sub 3} under high pressure and temperature. BaBkO{sub 3} presents high value of Seebeck coefficient and power factor at room temperature with a value of 135 μVK{sup − 1} and 1.5 × 10{sup 12}μWcm{sup − 1} K{sup − 2} s{sup − 1}, respectively. The possession of half metallicity and high Seebeck coefficient makes BaBkO{sub 3} suitable for spintronic and thermoelectric device applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

ON THE GRAVITATIONAL STABILITY OF GRAVITO-TURBULENT ACCRETION DISKS

Low mass, self-gravitating accretion disks admit quasi-steady, “gravito-turbulent” states in which cooling balances turbulent viscous heating. However, numerical simulations show that gravito-turbulence cannot be sustained beyond dynamical timescales when the cooling rate or corresponding turbulent viscosity is too large. The result is disk fragmentation. We motivate and quantify an interpretation of disk fragmentation as the inability to maintain gravito-turbulence due to formal secondary instabilities driven by: (1) cooling, which reduces pressure support; and/or (2) viscosity, which reduces rotational support. We analyze the axisymmetric gravitational stability of viscous, non-adiabatic accretion disks with internal heating, external irradiation, and cooling in the shearing box approximation. We consider parameterized cooling functions in 2D and 3D disks, as well as radiative diffusion in 3D. We show that generally there is no critical cooling rate/viscosity below which the disk is formally stable, although interesting limits appear for unstable modes with lengthscales on the order of the disk thickness. We apply this new linear theory to protoplanetary disks subject to gravito-turbulence modeled as an effective viscosity, and cooling regulated by dust opacity. We find that viscosity renders the disk beyond ∼60 au dynamically unstable on radial lengthscales a few times the local disk thickness. This is coincident with the empirical condition for disk fragmentation based on a maximum sustainable stress. We suggest turbulent stresses can play an active role in realistic disk fragmentation by removing rotational stabilization against self-gravity, and that the observed transition in behavior from gravito-turbulent to fragmenting may reflect instability of the gravito-turbulent state itself.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

An improved multilevel Monte Carlo method for estimating probability distribution functions in stochastic oil reservoir simulations

In this paper, we develop an improved multilevel Monte Carlo (MLMC) method for estimating cumulative distribution functions (CDFs) of a quantity of interest, coming from numerical approximation of large-scale stochastic subsurface simulations. Compared with Monte Carlo (MC) methods, that require a significantly large number of high-fidelity model executions to achieve a prescribed accuracy when computing statistical expectations, MLMC methods were originally proposed to significantly reduce the computational cost with the use of multifidelity approximations. The improved performance of the MLMC methods depends strongly on the decay of the variance of the integrand as the level increases. However, the main challenge in estimating CDFs is that the integrand is a discontinuous indicator function whose variance decays slowly. To address this difficult task, we approximate the integrand using a smoothing function that accelerates the decay of the variance. In addition, we design a novel a posteriori optimization strategy to calibrate the smoothing function, so as to balance the computational gain and the approximation error. The combined proposed techniques are integrated into a very general and practical algorithm that can be applied to a wide range of subsurface problems for high-dimensional uncertainty quantification, such as a fine-grid oil reservoir model considered in this effort. The numerical results reveal that with the use of the calibrated smoothing function, the improved MLMC technique significantly reduces the computational complexity compared to the standard MC approach. Finally, we discuss several factors that affect the performance of the MLMC method and provide guidance for effective and efficient usage in practice.

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

Configuration splitting of the Gamow-Teller resonance in antimony isotopes: Is this a real or a virtual effect?

The Gamow-Teller strength functions for a number of antimony isotopes were calculated within a semimicroscopic approach based on the continuum version of the charge-exchange quasiparticle random-phase approximation and on phenomenologically taking into account the fragmentation effect. The structural effect of splitting of the main maximum of the Gamow-Teller resonance in these isotopes was confirmed. Experimental data on the excitation of this resonance in a direct and a resonance reaction for the {sup 118}Sn parent nucleus were analyzed with allowance for this effect.

74 ATOMIC AND MOLECULAR PHYSICS↗