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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗
DOE OSTI2025

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗
DOE OSTI2017

Self-interaction effects on charge-transfer collisions

In this article, we investigate the role of the self-interaction error in the simulation of collisions using time-dependent density functional theory (TDDFT) and Ehrenfest dynamics. In addition, we compare many different approximations of the exchange and correlation potential, using as a test system the collision of H + + CH 4 at 30 eV. We find that semilocal approximations, like the Perdew-Burke- Ernzerhof (PBE), and even hybrid functionals, such as the Becke, 3-parameter, Lee-Yang-Parr (B3LYP), produce qualitatively incorrect predictions for the scattering of the proton. This discrepancy appears because the self-interaction error allows the electrons to jump too easily to the proton, leading to radically different forces with respect to the non-self-interacting case. Lastly, from our results, we conclude that using a functional that is self-interaction free is essential to properly describing charge-transfer collisions between ions and molecules in TDDFT.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Investigation of Transport Parameters of Graphene-Based Nanostructures

The paper presents results of computer simulation of the main transport parameters of nanostructures obtained through the row-by-row removal of carbon atoms from graphene ribbon. Research into the electrical parameters is carried out within the density functional theory using the non-equilibrium Green functions in the local-density approximation. Virtual NanoLab based on Atomistix ToolKit is used to construct structures and analyze simulation results. Current-voltage characteristics, differential conductivity and transmittance spectra of nanostructures are calculated at different values of bias voltage. It is found that there is a large region of negative differential resistance in current-voltage characteristics of nanostructures caused by resonant tunneling of quasi-particles. Differential (dI/dV) characteristic also has similar changes. The obtained results can be useful for building novel electronic devices in the field of nanoelectronics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

A phase-based interior penalty discontinuous Galerkin method for the Helmholtz equation with spatially varying wavenumber

This paper is concerned with an interior penalty discontinuous Galerkin (IPDG) method based on a flexible type of non-polynomial local approximation space for the Helmholtz equation with varying wavenumber. The local approximation space consists of multiple polynomial-modulated phase functions which can be chosen according to the phase information of the solution. We obtain some approximation properties for this space and a priori L 2 error estimates for the h-convergence of the IPDG method using duality argument. We also provide ample numerical examples to show that, building phase information into the local spaces often gives more accurate results comparing to using the standard polynomial spaces.

42 ENGINEERING↗
DOE OSTI2019

Hydrogen Effect on Electron-Phonon Interactions in L1{sub 0} FePd

The effect of H position in L1{sub 0}FePd magnetic equiatomic layered structure on electron-phonon interactions is studied via ab initio pseudopotentials density functional approach in the general gradient approximation. The structures investigated are L1{sub 0}FePd (parent structure), and the specific hydrides 2(FePd)1H (H@Fe layer) and 2(FePd)2H (H@Pd layer). At ground state (0 K), electronic band structure results demonstrate the emergence of Fermi surfaces (FS) in topologies of small isolated multipockets in the hydrides. Using the harmonic approximation I find the phonon energy dispersion of parent and H@Fe simple and classical but for H@Pd complex. I linked this latter structure to its topology of FS, known as Kohn anomaly. Using McMillan-Eliashberg model, I get superconducting transition temperature (T{sub c}) values ≈ 0.00 K, 3.48 K, and 19.75 K for the parent, H@Fe and H@Pd respectively. Therefore, these hypothetical hydrides prove that H position in the structure has a direct influence on T{sub c} values and consequently the suppression of magnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Approximating Photo- z PDFs for Large Surveys

Modern galaxy surveys produce redshift probability density functions (PDFs) in addition to traditional photometric redshift (photo-z) point estimates. However, the storage of photo-z PDFs may present a challenge with increasingly large catalogs, as we face a trade-off between the accuracy of subsequent science measurements and the limitation of finite storage resources. This work presents qp, a Python package for manipulating parameterizations of one-dimensional PDFs, as suitable for photo-z PDF compression. We use qp to investigate the performance of three simple PDF storage formats (quantiles, samples, and step functions) as a function of the number of stored parameters on two realistic mock data sets, representative of upcoming surveys with different data qualities. We propose some best practices for choosing a photo-z PDF approximation scheme and demonstrate the approach on a science case using performance metrics on both ensembles of individual photo-z PDFs and an estimator of the overall redshift distribution function. We show that both the properties of the set of PDFs we wish to approximate and the fidelity metric(s) chosen affect the optimal parameterization. Additionally, we find that quantiles and samples outperform step functions, and we encourage further consideration of these formats for PDF approximation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion

There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1D, 2D, and 3D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, t d0 . In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, t d0 . The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the first two terms for high-accuracy approximations (with less than 10-7 relative error) for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2D/3D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial in √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.

42 ENGINEERING↗
DOE OSTI2019

Improved approximations of non-Gaussian probability, force, and energy of a single polymer chain

The mechanical behavior of polymers such as their probability density function (PDF), strain energy, and entropic force, has long been described by the non-Gaussian statistical model. Non-Gaussian models are often approximated by the Kuhn-Grün (KG) distribution function, which is derived from the first-order approximation of the complex Rayleigh's exact Fourier integral distribution. Furthermore, the KG function is widely accepted in polymer physics, where the non-Gaussian theory is often used to describe energy of the chains with various flexibility ratios. However, the KG function is shown to be relevant only for long chains and becomes extremely inaccurate for chains with fewer than 40 segments. In comparison to the KG model, other approximations of the non-Gaussian statistical model are often less accurate, and those with higher accuracy are usually too complex to be implemented in large-scale simulations. Here a new accurate approximation of the non-Gaussian PDF, entropic force, and strain energy of a single chain subsequently is developed to describe the mechanics of a polymer chain. With a similar level of complexity, the presented approximations of non-Gaussian PDF, strain energy, and entropic force are at least 10 times more accurate than KG approximations and thus are an excellent alternative option to be used in micromechanical constitutive models.

42 ENGINEERING↗
DOE OSTI2018

Electronic, Magnetic and Optical Properties of TM (Ti, V, Cr, Mn and Co)-Doped Rare-Earth Nitride GdN: First-Principles Theory

The half-metallic ferromagnetic behavior of rare-earth nitride Gd{sub 0.95} TM{sub 0.05}N (TM = Ti, V, Cr, Mn and Co), based on diluted magnetic semiconductors (DMSs), is investigated using the Korringa–Kohn–Rostoker (KKR) method combined with the coherent potential approximation (CPA) within a framework of density functional theory (DFT). The energy difference between the ferromagnetic and disorder local moment states has been evaluated. The exchange interactions obtained from first-principles calculations resulted in ferromagnetic states with Curie temperatures within the ambient conditions. Moreover, the optical absorption spectra obtained by ab initio calculations confirm the ferromagnetic stability based on the charge state of magnetic impurities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

While time-dependent density functional theory has emerged as a method of choice for computing electronic spectra and response of molecules and materials, its reliability is hindered by the adiabatic approximation that is commonly made. In this work, we address one problematic aspect that arises from this approximation: an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. Here, we demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

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

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Insight into Structural, Electronic, Magnetic, Mechanical, and Thermodynamic Properties of Actinide Perovskite BaPuO{sub 3}

In this paper, we have performed a systematic investigation on structural, magnetic, electronic, mechanical, and thermodynamic properties of BaPuO{sub 3} perovskite within density functional theory (DFT) using generalized gradient approximations (GGA), onsite coulomb repulsion (GGA + U), and modified Becke-Johnson (mBJ). The calculated structural parameters were found in good agreement with the experimental results. A large value of magnetic moment equal to integer value of 4 μ{sub B} was obtained for the compound. The spin-polarized electronic band structure and density of states present 100% of spin polarization at Fermi level, resulting in half-metallic nature for the compound with spin-up states as metallic and spin- down states as a semiconducting. The elastic and mechanical properties have also been predicted. Moreover, we have calculated thermodynamic properties like Debye temperature (𝜃{sub D}), specific heat (C{sub V}), entropy (S), etc. using quasi-harmonic Debye model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

First Principle Mechanical and Thermodynamic Properties of Some TbX (X = S, Se) Compounds

The mechanical and thermodynamic results of some terbium chalcogenides, TbX (X = S, Se) compounds have been examined using full-potential linear augmented plane wave (FP-LAPW) method based on density functional theory. The local spin density approximation with Hubbard potential (LSDA+U) has been used for proper treatment of 4f-electrons of Tb. The magnetic moment of TbS and TbSe was calculated to be 5.44μ{sub B} and 5.57μ{sub B}, respectively. The ground state properties such as lattice parameter, bulk modulus, and pressure derivative of bulk modulus were obtained and found in good agreement with already available theory and experiment results. The mechanical properties such as Young modulus (E), shear modulus (G), Poisson ratio (υ), Cauchy pressure and B/G ratio were also predicted. Further, in thermodynamic properties, Debye temperatures (θ{sub D}), specific heats at constant volume(C{sub V}), thermal expansion co-efficient (α), and entropies (S) and their behavior with high temperature(0–1200 K) and pressure (0–50 GPa) were also predicted.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Correlated particle transport enables biological free energy transduction

Studies of biological transport frequently neglect the explicit statistical correlations among particle site occupancies (i.e., they use a mean-field approximation). Neglecting correlations sometimes captures biological function, even for out-of-equilibrium and interacting systems. Here, we show that neglecting correlations fails to describe free energy transduction, mistakenly predicting an abundance of slippage and energy dissipation, even for networks that are near reversible and lack interactions among particle sites. Interestingly, linear charge transport chains are well described without including correlations, even for networks that are driven and include site-site interactions typical of biological electron transfer chains. We examine three specific bioenergetic networks: a linear electron transfer chain (as found in bacterial nanowires), a near-reversible electron bifurcation network (as in complex III of respiration and other recently discovered structures), and a redox-coupled proton pump (as in complex IV of respiration).

bacterial nanowires↗
DOE OSTI2025

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Low-frequency Raman spectrum of 2D layered perovskites: Local atomistic motion or superlattice modes?

We report the low-frequency Raman spectrum (ω = 10 cm -1 –150 cm -1 ) of a wide variety of alkylammonium iodide based 2D lead halide perovskites (2D LHPs) as a function of A-site cation (MA = methylammonium and FA = formamidinium), octahedral layer thickness (n = 2–4), organic spacer chain length (butyl-, pentyl-, hexyl-), and sample temperature (T = 77 K–293 K). Using density functional theory calculations under the harmonic approximation for n = 2 BA:MAPbI, we assign several longitudinal/transverse optical phonon modes between 30 cm -1 and 100 cm -1 , the eigen displacements of which are analogous to that observed previously for octahedral twists/distortions in bulk MAPbI. Additionally, we propose an alternative assignment for low-frequency modes below this band (<30 cm -1 ) as zone-folded longitudinal acoustic phonons corresponding to the periodicity of the entire layered structure. Furthermore, we compare measured spectra to predictions of the Rytov elastic continuum model for zone-folded dispersion in layered structures. Our results are consistent across the various 2D LHPs studied herein, with energetic shifts of optical phonons corresponding to microscopic structural differences between materials and energetic shifts of acoustic phonons according to changes in the periodicity and elastic properties of the perovskite/organic subphases. This study highlights the importance of both the local atomic order and the superlattice structure on the vibrational properties of layered 2D materials.

2D materials↗
DOE OSTI2022

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗