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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Density and pair-density scaling for deriving the Euler equation in density-functional and pair-density-functional theory

A link between density and pair density functional theories is presented. Density and pair density scaling are used to derive the Euler equation in both theories. Density scaling provides a constructive way of obtaining approximations for the Pauli potential. The Pauli potential (energy) of the density functional theory is expressed as the difference of the scaled and original exchange-correlation potentials (energies).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Dual energy CT for attenuation correction with PET/CT

Purpose: The authors evaluate the energy dependent noise and bias properties of monoenergetic images synthesized from dual-energy CT (DECT) acquisitions. These monoenergetic images can be used to estimate attenuation coefficients at energies suitable for positron emission tomography (PET) and single-photon emission computed tomography (SPECT) imaging. This is becoming more relevant with the increased use of quantitative imaging by PET/CT and SPECT/CT scanners. There are, however, potential variations in the noise and bias of synthesized monoenergetic images as a function of energy. Methods: The authors used analytic approximations and simulations to estimate the noise and bias of synthesized monoenergetic images of water-filled cylinders with different shapes and the NURBS-based cardiac-torso (NCAT) phantom from 40 to 520 keV, the range of SPECT and PET energies. The dual-kVp spectra were based on the GE Lightspeed VCT scanner at 80 and 140 kVp with added filtration of 0.5 mm Cu. The authors evaluated strategies of noise suppression with sinogram smoothing and dose minimization with reduction of tube currents at the two kVp settings. The authors compared the impact of DECT-based attenuation correction with single-kVp CT-based attenuation correction on PET quantitation for the NCAT phantom for soft tissue and high-Z materials of bone and iodine contrast enhancement. Results: Both analytic calculations and simulations displayed the expected minimum noise value for a synthesized monoenergetic image at an energy between the mean energies of the two spectra. In addition the authors found that the normalized coefficient of variation in the synthesized attenuation map increased with energy but reached a plateau near 160 keV, and then remained constant with increasing energy up to 511 keV and beyond. The bias was minimal, as the linear attenuation coefficients of the synthesized monoenergetic images were within 2.4% of the known true values across the entire energy range. Compared with no sinogram smoothing, sinogram smoothing can dramatically reduce noise in the DECT-derived attenuation map. Through appropriate selection of tube currents for high and low kVp scans, DECT can deliver roughly the same amount of radiation dose as that of a single kVp CT scan, but could be used for PET attenuation correction with reduced bias in contrast agent regions by a factor of ∼2.6 and slightly reduced RMSE for the total image. Conclusions: When DECT is used for attenuation correction at higher energies, there is a noise amplification that is dependent on the energy of the synthesized monoenergetic image of linear attenuation coefficients. Sinogram smoothing reduces the noise amplification in DECT-derived attenuation maps without increasing bias. With an appropriate selection of CT techniques, a DECT scan with the same radiation dose as a single CT scan can result in a PET image with improved quantitative accuracy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Structural, electronic and thermal properties of Mo{sub 3}Ir superconductors

Self consistent first principle calculations on superconducting material Mo{sub 3}Ir of A15 phase have been performed to understand their fundamental characteristics of the structural, electronic and thermal properties. The bulk modulus (B{sub 0}), Debye temperature (θ{sub D}), density of states (N (E{sub F})) and electronic specific heat coefficient (γ) have been computed in terms of the electronic structure results obtained by using the tight-binding linear muffin tin orbital (TB-LMTO) method based on the density functional theory (DFT) within the local density approximation (LDA). Structural, electronic and thermal properties calculated here are found to corroborate well with the experimental and theoretical results of literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Modelling charge transfer reactions with the frozen density embedding formalism

The frozen density embedding (FDE) subsystem formulation of density-functional theory is a useful tool for studying charge transfer reactions. In this work charge-localized, diabatic states are generated directly with FDE and used to calculate electronic couplings of hole transfer reactions in two {pi}-stacked nucleobase dimers of B-DNA: 5{sup '}-GG-3{sup '} and 5{sup '}-GT-3{sup '}. The calculations rely on two assumptions: the two-state model, and a small differential overlap between donor and acceptor subsystem densities. The resulting electronic couplings agree well with benchmark values for those exchange-correlation functionals that contain a high percentage of exact exchange. Instead, when semilocal GGA functionals are used the electronic couplings are grossly overestimated.

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

Photoelectron spectra of early 3 d -transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3 d -transition metal dioxide anions, Sc O 2 - , Ti O 2 - , V O 2 - , Cr O 2 - , and Mn O 2 - , are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point ( G 0 W 0 @PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV–0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew–Burke–Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G 0 W 0 @PBE0, the G -only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3 d -transition metal dioxide anions. Adding eigenvalue self-consistency also in W , on the other hand, worsens the agreement with experiment overall. Overall, our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

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

Applicability of Density Functional Theory in Reproducing Accurate Vibrational Spectra of Surface Bound Species

The structural equilibrium parameters, the adsorption energies, and the vibrational frequencies of the nitrogen molecule and the hydrogen atom adsorbed on the (111) surface of rhodium have been investigated using different generalized-gradient approximation (GGA), nonlocal correlation, meta-GGA, and hybrid functionals, namely, Perdew, Burke, and Ernzerhof (PBE), Revised-RPBE, vdW-DF, Tao, Perdew, Staroverov, and Scuseria functional (TPSS), and Heyd, Scuseria, and Ernzerhof (HSE06) functional in the plane wave formalism. Among the five tested functionals, nonlocal vdW-DF and meta-GGA TPSS functionals are most successful in describing energetics of dinitrogen physisorption to the Rh(111) surface, while the PBE functional provides the correct chemisorption energy for the hydrogen atom. It was also found that TPSS functional produces the best vibrational spectra of the nitrogen molecule and the hydrogen atom on rhodium within the harmonic formalism with the error of 22.62 and 21.1% for the NAN stretching and RhAH stretching frequency. Thus, TPSS functional was proposed as a method of choice for obtaining vibrational spectra of low weight adsorbates on metallic surfaces within the harmonic approximation. At the anharmonic level, by decoupling the RhAH and NAN stretching modes from the bulk phonons and by solving one- and two-dimensional Schr€odinger equation associated with the RhAH, RhAN, and NAN potential energy we calculated the anharmonic correction for NAN and RhAH stretching modes as 231 cm21 and 277 cm21 at PBE level. Anharmonic vibrational frequencies calculated with the use of the hybrid HSE06 function are in best agreement with available experiments.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

A practical model for fluid flow in discrete-fracture porous media by using the numerical manifold method

Here in this study, a numerical manifold method (NMM) model is developed to analyze flow in porous media with discrete fractures in a non-conforming mesh. This new model is based on a two-cover-mesh system with a uniform triangular mathematical mesh and boundary/fracture-divided physical covers, where local independent cover functions are defined. The overlapping parts of the physical covers are elements where the global approximation is defined by the weighted average of the physical cover functions. The mesh is generated by a tree-cutting algorithm. A new model that does not introduce additional degrees of freedom (DOF) for fractures was developed for fluid flow in fractures. The fracture surfaces that belong to different physical covers are used to represent fracture flow in the direction of the fractures. In the direction normal to the fractures, the fracture surfaces are regarded as Dirichlet boundaries to exchange fluxes with the rock matrix. Furthermore, fractures that intersect with Dirichlet or Neumann boundaries are considered. Simulation examples are designed to verify the efficiency of the tree-cutting algorithm, the calculation's independency from the mesh orientation, and accuracy when modeling porous media that contain fractures with multiple intersections and different orientations. The simulation results show good agreement with available analytical solutions. Finally, the model is applied to cases that involve nine intersecting fractures and a complex network of 100 fractures, both of which achieve reasonable results. The new model is very practical for modeling flow in fractured porous media, even for a geometrically complex fracture network with large hydraulic conductivity contrasts between fractures and the matrix.

58 GEOSCIENCES↗
DOE OSTI2021

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Adaptive Sparse Grid Discretization (ASGarD)

This code solves the Vlasov-Poisson equations in 6-dimensions without approximation. It uses sparse grid discretization to solve the high-dimensional PDEs. The code functions as a discontinuous Galerkin finite element solver. Advantages/benefits include no approximations, supports high-dimensionality, leadership-scale capability, and implicit time advance.

Hahn, Steven [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2013

Statistics of general functions of a Gaussian field-application to non-Gaussianity from preheating

We provide a general formula for calculating correlators of arbitrary function of a Gaussian field. This work extends the standard leading-order approximation based on the δN formalism to the case where truncation of the δN at some low order does not yield the correct answer. As an application of this formula, we investigate 2, 3 and 4-point functions of the primordial curvature perturbation generated in the massless preheating model by approximating the mapping between the curvature perturbation and the Gaussian field as a sum of the many spiky normal distribution functions as suggested by lattice calculations. We also discuss observational consequences of this case and show that trispectrum would be a key observable to search signature of preheating in the CMB map. It is found the forms of the curvature correlation functions for any δN, at the leading order in the correlator of the Gaussian field, coincide with the standard local type ones. Within this approximation, it is also found that the standard formula for the non-linearity parameters given by the product of the derivatives of the e-folding number still holds after we replace the bare e-folding number appearing in the original δN expansion with the one smoothed in the field space with a Gaussian window function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Combined hybrid functional and DFT+U calculations for metal chalcogenides

In the density-functional studies of materials with localized electronic states, the local/semilocal exchange-correlation functionals are often either combined with a Hubbard parameter U as in the LDA+U method or mixed with a fraction of exactly computed (Fock) exchange energy yielding a hybrid functional. Although some inaccuracies of the semilocal density approximations are thus fixed to a certain extent, the improvements are not sufficient to make the predictions agree with the experimental data. Here, we put forward the perspective that the hybrid functional scheme and the LDA+U method should be treated as complementary, and propose to combine the range-separated Heyd-Scuseria-Ernzerhof (HSE) hybrid functional with the Hubbard U. We thus present a variety of HSE+U calculations for a set of II-VI semiconductors, consisting of zinc and cadmium monochalcogenides, along with comparison to the experimental data. Our findings imply that an optimal value U{sup *} of the Hubbard parameter could be determined, which ensures that the HSE+U{sup *} calculation reproduces the experimental band gap. It is shown that an improved description not only of the electronic structure but also of the crystal structure and energetics is obtained by adding the U{sup *} term to the HSE functional, proving the utility of HSE+U{sup *} approach in modeling semiconductors with localized electronic states.

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

A COMPARATIVE STUDY OF OPTICAL/ULTRAVIOLET VARIABILITY OF NARROW-LINE SEYFERT 1 AND BROAD-LINE SEYFERT 1 ACTIVE GALACTIC NUCLEI

The ensemble optical/ultraviolet (UV) variability of narrow-line Seyfert 1 (NLS1)-type active galactic nuclei (AGNs) is investigated, based on a sample selected from the Sloan Digital Sky Survey (SDSS) Stripe 82 region with multi-epoch photometric scanning data. As a comparison, a control sample of broad-line Seyfert 1 (BLS1)-type AGNs is also incorporated. To quantify properly the intrinsic variation amplitudes and their uncertainties, a novel method of parametric maximum likelihood is introduced that has, as we argued, certain virtues over previously used methods. The majority of NLS1-type AGNs exhibit significant variability on timescales from about 10 days to a few years with, however, smaller amplitudes on average compared to BLS1-type AGNs. About 20 NLS1-type AGNs that show relatively large variations are presented and may deserve future monitoring observations, for instance, reverberation mapping. The averaged structure functions of variability, constructed using the same maximum likelihood method, show remarkable similarity in shape for the two types of AGNs on timescales longer than about 10 days, which can be approximated by a power law or an exponential function. This, along with other similar properties, such as the wavelength-dependent variability, is indicative of a common dominant mechanism responsible for the long-term optical/UV variability of both NLS1- and BLS1-type AGNs. Toward the short timescales, however, there is tentative evidence that the structure function of NLS1-type AGNs continues to decline, whereas that of BLS1-type AGNs flattens with some residual variability on timescales of days. If this can be confirmed, it may suggest that an alternative mechanism, such as X-ray reprocessing, starts to dominate in BLS1-type AGNs, but not in NLS1-type AGNs, on such timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

First-principles modeling of plasmons in aluminum under ambient and extreme conditions

The theoretical understanding of plasmon behavior is crucial for an accurate interpretation of inelastic scattering diagnostics in many experiments. In this work, we highlight the utility of linear response time-dependent density functional theory (LR-TDDFT) as a first-principles framework for consistently modeling plasmon properties. We provide a comprehensive analysis of plasmons in aluminum from ambient to warm dense matter conditions and assess typical properties such as the dynamical structure factor, the plasmon dispersion, and the plasmon lifetime. We compare our results with scattering measurements and with other TDDFT results as well as models such as the random phase approximation, the Mermin approach, and the dielectric function obtained using static local field corrections of the uniform electron gas parametrized from path-integral Monte Carlo simulations. We conclude that results for the plasmon dispersion and lifetime are inconsistent between experiment and theories and that the common practice of extracting and studying plasmon dispersion relations is an insufficient procedure to capture the complicated physics contained in the dynamic structure factor in its full breadth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Regarding the use and misuse of retinal protonated Schiff base photochemistry as a test case for time-dependent density-functional theory

The excited-state relaxation of retinal protonated Schiff bases (PSBs) is an important test case for biological applications of time-dependent (TD) density-functional theory (DFT). While well-known shortcomings of approximate TD-DFT might seem discouraging for application to PSB relaxation, progress continues to be made in the development of new functionals and of criteria allowing problematic excitations to be identified within the framework of TD-DFT itself. Furthermore, experimental and theoretical ab initio advances have recently lead to a revised understanding of retinal PSB photochemistry, calling for a reappraisal of the performance of TD-DFT in describing this prototypical photoactive system. Here, we re-investigate the performance of functionals in (TD-)DFT calculations in light of these new benchmark results, which we extend to larger PSB models. We focus on the ability of the functionals to describe primarily the early skeletal relaxation of the chromophore and investigate how far along the out-of-plane pathways these functionals are able to describe the subsequent rotation around formal single and double bonds. Conventional global hybrid and range-separated hybrid functionals are investigated as the presence of Hartree-Fock exchange reduces problems with charge-transfer excitations as determined by the Peach-Benfield-Helgaker-Tozer Λ criterion and by comparison with multi-reference perturbation theory results. While we confirm that most functionals cannot render the complex photobehavior of the retinal PSB, do we also observe that LC-BLYP gives the best description of the initial part of the photoreaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Multi-Function AAA Algorithm Applied to Frequency Dependent Line Modeling

Modeling of power distribution system components that are valid for a wide range of frequencies are crucial for highly accurate modeling of electromagnetic transient (EMT) events. This has recently become of interest due to the improvements needed for the resilient operation of distribution systems. Vector fitting (VF) is a very popular and commonly used algorithm for wide band representations of power system components in EMT simulations. In this research, we present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages with respect to VF using examples of approximations of admittance matrices discussed in the literature. We show that MIAAA not only outperforms VF in terms of achieving better accuracy using lesser number of poles, but also has no numerical issues achieving convergence. In contrast to VF, MIAAA 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. The novelty of this research work is the use of recent mathematical results to solve existing challenges in distribution system modeling and to develop rational approximations for power system models that intend to be optimal in terms of accuracy and performance.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2013

Lipid composition can determine membrane protein orientation during reconstitution into proteoliposomes

Artificial membrane systems allow researchers to study the structure and function of membrane proteins in a matrix that approximates their natural environment and to integrate these proteins in ex vivo devices such as electronic biosensors, thin-film protein arrays, or biofuel cells. Given that most membrane proteins have vectorial functions, both functional studies and applications require effective control over protein orientation within a lipid bilayer. In this work, we explored the role of the bilayer surface charge in determining transmembrane protein orientation and functionality during formation of proteoliposomes. We reconstituted a model vectorial ion pump, proteorhodopsin, in liposomes of opposite charges and varying charge densities and determined the resultant protein orientation. Antibody-binding assay and proteolysis of proteoliposomes showed physical evidence of preferential orientation, and functional assays verified the vectorial nature of ion transport in this system. Our results indicate that the manipulation of lipid composition can indeed control orientation of an asymmetrically charged membrane protein, proteorhodopsin, in liposomes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Developing new and understanding old approximations in TDDFT

When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange–correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn–Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange–correlation potential from a coupling-constant path integration.

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

Self-consistent RPA based on a many-body vacuum

Self-Consistent RPA is extended in a way so that it is compatible with a variational ansatz for the ground-state wave function as a fermionic many-body vacuum. Employing the usual equation-of-motion technique, we arrive at extended RPA equations of the Self-Consistent RPA structure. In principle the Pauli principle is, therefore, fully respected. However, the correlation functions entering the RPA matrix can only be obtained from a systematic expansion in powers of some combinations of RPA amplitudes. We demonstrate for a model case that this expansion may converge rapidly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗