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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

The simplest response theory methods for computing vertical excitation spectra in condensed-phase are configuration interaction with single excitations (CIS) and linear-response time-dependent density functional theory (TDDFT) within the Tamm–Dancoff approximation. In applications to X-ray absorption spectroscopy (XAS), methods like CIS and TDDFT that codify only single excitations into the wave function are prone to catastrophic errors in main-edge and post-edge features whose shapes act as a crucial fingerprint in structural analyses of liquids. Here, we show that these errors manifest primarily due to a lack of orbital relaxation in conventional linear-response theories and that core-ionized (n - 1-electron) references, like those of electron-affinity TDDFT, can eliminate the errors in the spectral profile, even in the highest-energy parts of the post-edge. Crucially, we find that single excitations atop core-ionized references are sufficient to elucidate liquid-phase XAS spectra with semi-quantitative accuracy, opening the door for methods like electron-affinity CIS/TDDFT to be used as efficient alternatives to higher-order wave function approaches.

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

The response of soil carbon storage and microbially mediated carbon turnover to simulated climatic disturbance in a northern peatland forest. Revisiting the concept of soil organic matter recalcitrance

The goal of this project was to investigate changes in the structure of dissolved and solid phase organic matter, the production of CO 2 and CH 4 , and the composition of decomposer microbial communities in response to the climatic forcing of environmental processes that determine the balance between carbon gas production versus storage and sequestration in peatlands. Cutting-edge analytical chemistry and next generation sequencing of microbial genes were been applied to habitats at the Marcell Experimental Forest (MEF), where the US DOE’s Oak Ridge National Laboratory and the USDA Forest Service are constructing a large-scale ecosystem study entitled, “Spruce and Peatland Responses Under Climatic and Environmental Change”(SPRUCE). Our study represented a comprehensive characterization of the sources, transformation, and decomposition of organic matter in the S1 bog at MEF. Multiple lines of evidence point to distinct, vertical zones of organic matter transformation: 1) the acrotelm consisting of living mosses, root material, and newly formed litter (0-30 cm), 2) the mesotelm, a mid-depth transition zone (30-75 cm) characterized by labile organic C compounds and intense decomposition, and 3) the underlying catotelm (below 75cm) characterized by refractory organic compounds as well as relatively low decomposition rates. These zones are in part defined by physical changes in hydraulic conductivity and water table depth. O-alkyl-C, which represents the carbohydrate fraction in the peat, was shown to be an excellent proxy for soil decomposition rates. The carbon cycle in deep peat was shown to be fueled by modern carbon sources further indicating that hydrology and surface vegetation play a role in belowground carbon cycling. We provide the first metagenomic study of an ombrotrophic peat bog, with novel insights into microbial specialization and functions in this unique terrestrial ecosystem. Vertical structuring of microbial communities closely paralleled the chemical evolution of peat, with large shifts in microbial populations occurring in the biogeochemical hotspot, the mesotelm, where the highest rates of decomposition were detected. Stable isotope geochemistry and potential rates of methane production paralleled vertical changes in methanogen community composition to indicate a predominance of acetoclastic methanogenesis mediated by the Methanosarcinales in the mesotelm, while hydrogen-utilizing methanogens dominated in the deeper catotelm. Evidence pointed to the availability of phosphorus as well as nitrogen limiting the microbially-mediated turnover of organic carbon at MEF. Prior to initiation of the experimental treatments, our study provided key baseline data for the SPRUCE site on the vertical stratification of peat decomposition, key enzymatic pathways, and microbial taxa containing these pathways. The sensitivity of soil carbon turnover to climate change is strongly linked to recalcitrant carbon stocks and the temperature sensitivity of decomposition is thought to increase with increasing molecular complexity of carbon substrates. This project delivered results on how climate change perturbations impact the microbially-mediated turnover of recalcitrant organic matter in peatland forest soils, both under controlled conditions in the laboratory and at the ecosystem-scale in the field. This project revisited the concept of “recalcitrance” in the regulation of soil carbon turnover using a combination of natural abundance radiocarbon and optical spectroscopic measurements on bulk DOM, and high resolution molecular characterization of DOM. The project elucidated how organic matter reactivity and decomposition will respond to climate change in a both a qualitative (organic matter lability) and quantitiative (increased rates) manner. An Aromaticity Index was developed to represent a more direct and accurate parameter for modeling of DOM reactivity in peatlands. The abundance and community composition of soil microorganisms that mediate C cycling were interrogated with depth in the peat, with season, and in manipulated climate enclosures at unprecedented resolution. Therefore this project delivered strategic new insights on the functioning of peatland ecosystems that collectively store approximately one-third of the world's soil carbon. Furthermore, results from the detailed characterization of DOM lability and microbial community structure/ function will be employed to further develop biogeochemical models to include microbial respiration pathways as well as to track carbon flow with a term that incorporates relative reactivity based on aromaticity index. As it stands now, detailed soil organic matter structure and microbial parameters are not included in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Outside the Continental United States International Travel and Contagion Impact Quick Look Tool

ABSTRACT This paper describes a tool that will allow public health analysts to estimate infectious disease risk at the country level as a function of different international transportation modes. The prototype focuses on a cholera epidemic originating within Latin America or the Caribbean, but it can be expanded to consider other pathogens as well. This effort leverages previous work in collaboration with the Centers for Disease Control and Prevention to develop the International Travel to Community Impact (IT-CI) model, which analyzes and assesses potential international disease outbreaks then estimates the associated impacts to U.S. communities and the nation as a whole and orient it for use Outside the Continental United States (OCONUS). For brevity, we refer to this refined model as OIT-CI. First, we developed an operationalized meta-population spatial cholera model for Latin America and the Caribbean at the secondary administrative-level boundary. Secondly, we developed a robust function of human airline critical to approximating mixing patterns in the meta- population model. In the prototype version currently presented here, OIT-CI models a cholera epidemic originating in a Latin American or Caribbean country and spreading via airline transportation routes. Disease spread is modeled at the country level using a patch model with a connectivity function based on demographic, geospatial, and human transportation data. We have also identified data to estimate the water and health-related infrastructure capabilities of each country to include this potential impact on disease transmission.

Biosurveillance↗
DOE OSTI2017

Development of the self-consistent approximation and its application to the problem of magnetoelastic resonance in an inhomogeneous medium

Our previously proposed approximation involving both the first and second terms of the expansion of the vertex function is generalized to the system of two interacting wavefields of different physical nature. A system of self-consistent equations for the matrix Green’s function and matrix vertex function is derived. On the basis of this matrix generalization of the new self-consistent approximation, a theory of magnetoelastic resonance is developed for a ferromagnetic model, where the magnetoelastic coupling parameter ε(x) is inhomogeneous. Equations for magnetoelastic resonance are analyzed for one-dimensional inhomogeneities of the coupling parameter. The diagonal and off-diagonal elements of the matrix Green’s function of the system of coupled spin and elastic waves are calculated with the change in the ratio between the average value ε and rms fluctuation Δε of the coupling parameter between waves from the homogeneous case (ε ≠ 0, Δε = 0) to the extremely randomized case (ε = 0, Δε ≠ 0) at various correlation wavenumbers of inhomogeneities k{sub c}. For the limiting case of infinite correlation radius (k{sub c} = 0), in addition to approximate expressions, exact analytical expressions corresponding to the summation of all diagrams of elements of the matrix Green’s function are obtained. The results calculated for an arbitrary k{sub c} value in the new self-consistent approximation are compared to the results obtained in the standard self-consistent approximation, where only the first term of the expansion of the vertex function is taken into account. It is shown that the new approximation corrects disadvantages of the Green’s functions calculated in the standard approximation such as the dome shape of resonances and bends on the sides of resonance peaks. The appearance of a fine structure of the spectrum in the form of a narrow resonance on the Green’s function of spin waves and a narrow antiresonance on the Green’s function of elastic waves, which was previously predicted in the standard self-consistent approximation, is confirmed. With an increase in the parameter k{sub c}, the Green’s functions calculated in the standard and new approximations approach each other and almost coincide with each other at k{sub c}/k ≥ 0.5. At the same time, the results of this work indicate that the new self-consistent approximation has a certain advantage for studying the problems of stochastic radiophysics in media with long-wavelength inhomogeneities (small k{sub c} values), because it describes both the shape and width of peaks much better than the standard approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Efficient sampling of free energy landscapes with functions in Sobolev spaces

Molecular simulations of biological and physical phenomena generally involve sampling complicated, rough energy landscapes characterized by multiple local minima. In this work, we introduce a new family of methods for advanced sampling that draw inspiration from functional representations used in machine learning and approximation theory. As shown here, such representations are particularly well suited for learning free energies using artificial neural networks. As a system evolves through phase space, the proposed methods gradually build a model for the free energy as a function of one or more collective variables, from both the frequency of visits to distinct states and generalized force estimates corresponding to such states. Implementation of the methods is relatively simple and, more importantly, for the representative examples considered in this work, they provide computational efficiency gains of up to several orders of magnitude over other widely used simulation techniques.

Approximation theory↗
DOE OSTI2019

Effects of Co Doping and Point Defect on the Ferromagnetism of ZnO

The effects of different doping proportions of Co doping and oxygen vacancies (V{sub O}) on the ferromagnetism of zinc oxide (ZnO) are controversial. To solve this problem, we investigated the effects of Co doping and point defects (V{sub O}, V{sub Zn}, H{sub i}, and Zn{sub i}) on the ferromagnetism of ZnO through first-principle calculations by using generalized gradient approximation + U (GGA + U) under density functional theory. Results show that the closer the distance between Co and vacancies (V{sub O} or V{sub Zn}), the lower the formation energy and the more stable the system will be. Co-doped ZnO with V{sub O} or V{sub Zn} exhibits long-range ordered ferromagnetism and Curie temperature above the room temperature. Especially, Co-doped ZnO with V{sub Zn} shows a half-metallic behavior, resulting in 100% conduction hole polarizability, which is highly beneficial for dilute magnetic semiconductors (DMSs) as a hole injection source. The ferromagnetism of Zn{sub 14}CoO{sub 16} is caused by the double exchange effects among the electrons of the O–2p, Co–3d, and Zn–3d orbits mediated by the hole carriers after complexes were formed by the Co doping and Zn vacancy. Similarly, the origin of ferromagnetism in Zn{sub 15}CoO{sub 15} is derived from the double exchange effects among the electrons of O–2p, Co–3d, and Zn–4s states mediated by the electron carriers after complexes were formed by the Co doping and O vacancy. In addition, the result shows that Co-doped ZnO with interstitial H (H{sub i}) or Zn (Zn{sub i}) is unfavorable for ferromagnetism and should be avoided experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Meta-GGA exchange-correlation free energy density functional to increase the accuracy of warm dense matter simulations

We discuss strategies for thermalization of the ground-state meta-generalized gradient approximation (meta-GGA) exchange-correlation (XC) functionals. A simple but accurate scheme is implemented via universal additive thermal correction to XC using a perturbative-like self-consistent approach. The additive correction with explicit temperature dependence is applied to the ground-state deorbitalized, strongly constrained and appropriately normed (SCAN-L) meta-GGA XC leading to thermal XC functional denoted here as T-SCAN-L. Thermal T-SCAN-L meta-GGA functional shows significant improvement in density functional theory calculation accuracy for warm dense matter by a factor of 3 to 10, achieving unprecedented accuracy of total pressure between a few tenths and 1% when compared to traditional XC functionals, as demonstrated by the comparison to pathintegral Monte Carlo simulations for helium equation of state. Furthermore, the T-SCAN-L calculations of dc conductivity of warm dense aluminum also give better agreement with experiments over other XC functionals such as PBE and SCAN-L.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Clusters in nuclear matter and the equation of state for astrophysical applications

Thermodynamic properties and the composition of dense matter are affected strongly by correlations leading to the formation of nuclei at densities below nuclear saturation density. These features of the equation of state (EoS) are essential in astrophysical applications, e.g. the simulation of core-collapse supernovae. Theoretical approaches to consider correlation effects on different levels of approximation are discussed. A generalized relativistic density functional approach is presented that exhibits the correct limits at low and high densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

An Empirical Fitting Method for Type Ia Supernova Light Curves: A Case Study of SN 2011fe

We present a new empirical fitting method for the optical light curves of Type Ia supernovae (SNe Ia). We find that a variant broken-power-law function provides a good fit, with the simple assumption that the optical emission is approximately the blackbody emission of the expanding fireball. This function is mathematically analytic and is derived directly from the photospheric velocity evolution. When deriving the function, we assume that both the blackbody temperature and photospheric velocity are constant, but the final function is able to accommodate these changes during the fitting procedure. Applying it to the case study of SN 2011fe gives a surprisingly good fit that can describe the light curves from the first-light time to a few weeks after peak brightness, as well as over a large range of fluxes (∼5 mag, and even ∼7 mag in the g band). Since SNe Ia share similar light-curve shapes, this fitting method has the potential to fit most other SNe Ia and characterize their properties in large statistical samples such as those already gathered and in the near future as new facilities become available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.

Jin, Shi↗
DOE OSTI2019

Crystal-momentum-resolved contributions to high-order harmonic generation in solids

We analytically and numerically investigate the emission of high-order harmonic radiation from model solids by intense few-cycle midinfrared laser pulses. In single-active-electron approximation, we expand the active electron's wave function in a basis of adiabatic Houston states and describe the solid's electronic band structure in terms of an adjustable Kronig-Penney model potential. For high-order harmonic generation (HHG) from MgO crystals, we examine spectra from two-band and converged multiband numerical calculations. We discuss the characteristics of intra- and interband contributions to the HHG spectrum for computations including initial crystal momenta either from the Γ point at the center of the first Brioullin zone (BZ) only or from the entire first BZ. For sufficiently high intensities of the driving laser field, we find relevant contributions to HHG from the entire first BZ. Here, based on numerically calculated spectra, we scrutinize the cutoff harmonic orders as a function of the laser peak intensity and find good qualitative agreement with our analytical saddle-point-approximation predictions and published theoretical data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

The Dirac equation in electronic structure calculations: Accurate evaluation of DFT predictions for actinides

Brooks, Johansson, and Skriver, using the LMTO-ASA method and considerable insight, were able to explain many of the ground state properties of the actinides. In the many years since this work was done, electronic structure calculations of increasing sophistication have been applied to actinide elements and compounds, attempting to quantify the applicability of DFT to actinides and actinide compounds and to try to incorporate other methodologies (i.e. DMFT) into DFT calculations. Through these calculations, the limits of both available density functionals and ad hoc methodologies are starting to become clear. However, it has also become clear that approximations used to incorporate relativity are not adequate to provide rigorous tests of the underlying equations of DFT, not to mention ad hoc additions. In this talk, we describe the result of full-potential LMTO calculations for the elemental actinides, comparing results obtained with a full Dirac basis with those obtained from scalar-relativistic bases, with and without variational spin-orbit. This comparison shows that the scalar relativistic treatment of actinides does not have sufficient accuracy to provide a rigorous test of theory and that variational spin-orbit introduces uncontrolled errors in the results of electronic structure calculations on actinide elements.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Transport and optical properties of warm dense aluminum in the two-temperature regime: Ab initio calculation and semiempirical approximation

This work is devoted to the investigation of transport and optical properties of liquid aluminum in the two-temperature case. At first optical properties, static electrical, and thermal conductivities were obtained in the ab initio calculation which is based on the quantum molecular dynamics, density functional theory, and the Kubo-Greenwood formula. Then the semiempirical approximation was constructed based on the results of our simulation. This approximation yields the dependences σ{sub 1{sub D{sub C}}}∝1/T{sub i}{sup 0.25} and K∝T{sub e}/T{sub i}{sup 0.25} for the static electrical conductivity and thermal conductivity, respectively, for liquid aluminum at ρ = 2.70 g/cm{sup 3}, 3 kK ≤ T{sub i} ≤ T{sub e} ≤ 20 kK. Our results are well described by the Drude model with the effective relaxation time τ∝T{sub i}{sup −0.25}. We have considered a number of other models for the static electrical and thermal conductivities of aluminum, they are all reduced in the low-temperature limit to the Drude model with different expressions for the relaxation time τ. Our results are not consistent with the models in which τ∝T{sub i}{sup −1} and support the models which use the expressions with the slower decrease of the relaxation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Electronic and magnetic properties of small rhodium clusters

We report a theoretical study of the electronic and magnetic properties of rhodium-atomic clusters. The lowest energy structures at the semi-empirical level of rhodium clusters are first obtained from a novel global-minimum search algorithm, known as PTMBHGA, where Gupta potential is used to describe the atomic interaction among the rhodium atoms. The structures are then re-optimized at the density functional theory (DFT) level with exchange-correlation energy approximated by Perdew-Burke-Ernzerhof generalized gradient approximation. For the purpose of calculating the magnetic moment of a given cluster, we calculate the optimized structure as a function of the spin multiplicity within the DFT framework. The resultant magnetic moments with the lowest energies so obtained allow us to work out the magnetic moment as a function of cluster size. Rhodium atomic clusters are found to display a unique variation in the magnetic moment as the cluster size varies. However, Rh{sub 4} and Rh{sub 6} are found to be nonmagnetic. Electronic structures of the magnetic ground-state structures are also investigated within the DFT framework. The results are compared against those based on different theoretical approaches available in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Structural and electronic properties of Sr{sub x}Ba{sub 1-x}SnO{sub 3} from first principles calculations

Neutron diffraction data for Sr{sub x}Ba{sub 1-x}SnO{sub 3} (x=0.0, 0.2, 0.4, 0.6, 0.8 and 1.0) solid solutions were used as inputs to obtain optimized geometries and electronic properties using the density functional theory (DFT) formalism considering both the local density and generalized gradient approximations, LDA and GGA, respectively. The crystal structures and SnO{sub 6} octahedra tilting angles found after total energy minimization agree well with experiment, specially for the GGA data. Elastic constants were also obtained and compared with theoretical and experimental results for cubic BaSnO{sub 3}. While the alloys with cubic unit cell have an indirect band gap, tetragonal and orthorhombic alloys exhibit direct band gaps (exception made to x=1.0). The Kohn-Sham minimum electronic band gap oscillates from 1.52 eV (cubic x=0.0, LDA) to 2.61 eV (orthorhombic x=1.0, LDA), and from 0.74 eV (cubic BaSnO{sub 3}, GGA) to 1.97 eV (orthorhombic SrSnO{sub 3}, GGA). Parabolic interpolation of bands has allowed us to estimate the effective masses for charge carriers, which are shown to be anisotropic and larger for holes. - Graphical Abstract: Highlights: Black-Right-Pointing-Pointer DFT calculations were performed on Sr{sub x}Ba{sub 1-x}SnO{sub 3} solid solutions. Black-Right-Pointing-Pointer Calculated crystal structures agree well with experiment. Black-Right-Pointing-Pointer Alloys have direct or indirect gaps depending on the Sr molar fraction. Black-Right-Pointing-Pointer The Kohn-Sham gap variation from x=0.0 to x=1.0 is close to the experimental value. Black-Right-Pointing-Pointer Carrier effective masses are very anisotropic, specially for holes.

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

Investigation of structural stability and elastic properties of Zrh and Zrh{sub 2}: A first principles study

The electronic, structural and mechanical properties of ZrH and ZrH{sub 2} are investigated by means of first principles calculation based on density functional theory as implemented in VASP code with generalized gradient approximation. The calculated ground state properties are in good agreement with previous experimental and other theoretical results. Among the six crystallographic structures considered for ZrH, ZB phase is found to be the most stable phase, whereas ZrH{sub 2} is energetically stable in tetragonal structure at ambient condition. A structural phase transition from ZB→NaCl at a pressure 10 GPa is predicted for ZrH.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

The structural and electronic properties of cubic AgMO{sub 3} (M=Nb, Ta) by first principles calculations

We report the electronic structure of the AgMO{sub 3}(M=Nb, Ta) within the frame work of density functional theory and calculations are performed within the generalized gradient approximation (GGA) by using ultrasoft pseudopotentials. The calculated equilibrium lattice parameters and volumes are extracted from fitting of Birch third order equation of state and which are reasonable agreement with the available experimental results. The density of states,band structure of Ag(Nb,Ta)O{sub 3} reveals that the valance bands mostly occupied with O-2p and O-2s states and whereas conduction band occupied with Nb (Ta) 4d(5d) states including less contribution from Ag 5s states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗