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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Evaluating gyro-viscosity in the Kelvin-Helmholtz instability by kinetic simulations

In the present paper, the finite-Larmor-radius (gyro-viscous) term [K. V. Roberts and J. B. Taylor, Phys. Rev. Lett. 8, 197–198 (1962)] is evaluated by using a full kinetic Vlasov simulation result of the Kelvin-Helmholtz instability (KHI). The velocity field and the pressure tensor are calculated from the high-resolution data of the velocity distribution functions obtained by the Vlasov simulation, which are used to approximate the Finite-Larmor-Radius (FLR) term according to Roberts and Taylor [Phys. Rev. Lett. 8, 197–198 (1962)]. The direct comparison between the pressure tensor and the FLR term shows an agreement. It is also shown that the anisotropic pressure gradient enhanced the linear growth of the KHI when the inner product between the vorticity of the primary velocity shear layer and the magnetic field is negative, which is consistent with the previous FLR-magnetohydrodynamic simulation result. This result suggests that it is not sufficient for reproducing the kinetic simulation result by fluid simulations to include the FLR term (or the pressure tensor) only in the equation of motion for fluid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Communication: Effect of accidental mode degeneracy on Raman intensity in 2D materials: Hybrid functional study of bilayer phosphorene

Bulk black phosphorus has two optical phonon modes labeled as A{sub g}{sup 2} and B{sub 2u}, respectively, that are nearly degenerate in frequency. However, density functional theory calculations using local or semi-local functionals cannot reproduce this degeneracy. Here, we propose a hybrid functional approach aided by van der Waals (vdW) force fields, which can accurately describe the lattice dynamic and electronic properties of both bulk and few-layer black phosphorus (phosphorene). Using this approach we show that in bilayer phosphorene, the two Raman modes derived from the B{sub 2u} and A{sub g}{sup 2} modes could exhibit strong resonance as a result of the accidental degeneracy so that both modes could be observed in Raman experiment. Without the mode degeneracy, however, the Raman intensity of the B{sub 2u}-derived mode would be too weak to be observed. We further show that the accidental degeneracy is correlated to the applied strain, which enables Raman spectroscopy to be a powerful tool for characterizing built-in strains in 2D materials, e.g., due to the interaction with substrates, which has emerged as an important issue in vdW epitaxy.

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

Structure formation in the Deser-Woodard nonlocal gravity model: a reappraisal

In this work, we extend previous analyses of the structure formation in the f (□{sup −1} R ) model of nonlocal gravity proposed by Deser and Woodard (DW), which reproduces the background expansion of ΛCDM with no need of a cosmological constant nor of any dimensional constant beside Newton's one. A previous analysis based on redshift-space distortions (RSD) data concluded that the model was ruled out. In this work we revisit the issue and find that, when recast in a localized model, the DW model is not ruled out and actually gives a better fit to RSD data than ΛCDM. In fact, the DW model presents a suppressed growth of matter perturbations with respect to ΛCDM and a slightly lower value of σ{sub 8}, as favored by observations. We also produce analytical approximations of the two modified gravity functions, i.e. the anisotropic stress η and the relative change of Newton's constant Y , and of f σ{sub 8}( z ) as a function of redshift. Finally, we also show how much the fit depends on initial conditions when these are generalized with respect to a standard matter-dominated era.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE CLUSTERING OF GALAXIES IN THE SDSS-III BARYON OSCILLATION SPECTROSCOPIC SURVEY: LUMINOSITY AND COLOR DEPENDENCE AND REDSHIFT EVOLUTION

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over {approx}3300 deg{sup 2} in the redshift range 0.43 < z < 0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r{sub 0} = 7.93 {+-} 0.06 h {sup -1} Mpc and an index of {gamma} = 1.85 {+-} 0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k + e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

RF Gap Transformation in PARMILA

The PARMILA (Phase And Radial Motion in Ion Linear Accelerators) code is used to design ion linacs and simulate the beam dynamics therein. The original code was developed in the 1960’s and later the Los Alamos Accelerator Code Group developed a PC version that has been used worldwide. Recently, in an effort to modernize the code and use it as the basis for a high-performance beam dynamics simulator with the LANSCE linac, the details of the code have come under additional scrutiny. PARMILA is a “z code” that uses various transformations to approximately represent the effect of the beam line element on the beam. In the case of the RF gap, it treats the particle dynamics inside the gap as a combination of drift + kick at the mid-plane + drift. Recently, we uncovered an undocumented keyword “off-exact” in PARMILA that enables the code to use the full modified-Bessel functions instead of low-order approximations in the RF gap transformation. This change improved the accuracy of the transformation for off-axis particles. However, a discrepancy between the analytical formula and the code was found for the transformation of the transverse coordinates. The purpose of this tech note is to document the derivation of the gap transformation and the related DTL cell length calculation used in the PARMILA code.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Equations for the distributions of functionals of a random-walk trajectory in an inhomogeneous medium

Based on the random-trap model and using the mean-field approximation, we derive an equation that allows the distribution of a functional of the trajectory of a particle making random walks over inhomogeneous-lattice site to be calculated. The derived equation is a generalization of the Feynman-Kac equation to an inhomogeneous medium. We also derive a backward equation in which not the final position of the particle but its position at the initial time is used as an independent variable. As an example of applying the derived equations, we consider the one-dimensional problem of calculating the first-passage time distribution. We show that the average first-passage times for homogeneous and inhomogeneous media with identical diffusion coefficients coincide, but the variance of the distribution for an inhomogeneous medium can be many times larger than that for a homogeneous one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Predicting impurity spectral functions using machine learning

The Anderson Impurity Model (AIM) is a canonical model of quantum many-body physics. Here we investigate whether machine learning models, both neural networks (NN) and kernel ridge regression (KRR), can accurately predict the AIM spectral function in all of its regimes, from empty orbital, to mixed valence, to Kondo. To tackle this question, we construct two large spectral databases containing approximately 410 000 and 600 000 spectral functions of the single-channel impurity problem. We show that the NN models can accurately predict the AIM spectral function in all of its regimes, with pointwise mean absolute errors down to 0.003 in normalized units. We find that the trained NN models outperform models based on KRR and enjoy a speedup on the order of 10 5 over traditional AIM solvers. Finally, the required size of the training set of our model can be significantly reduced using farthest point sampling in the AIM parameter space, which is important for generalizing our method to more complicated multichannel impurity problems of relevance to predicting the properties of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Validating first-principles phonon lifetimes via inelastic neutron scattering

Phonon lifetimes are a key component of quasiparticle theories of transport; yet first-principles lifetimes are rarely directly compared with inelastic neutron scattering (INS) results. Existing comparisons show discrepancies even at temperatures where perturbation theory is expected to be reliable. In this paper, we demonstrate that the reciprocal space voxel (q voxel), which is the finite region in reciprocal space required in INS data analysis, must be explicitly accounted for within theory in order to draw a meaningful comparison. Here, we demonstrate accurate predictions of peak widths of the scattering function when accounting for the q voxel in CaF 2 and ThO 2 . Passing this test implies high fidelity of the phonon interactions and the approximations used to compute the Green's function, serving as a critical benchmark of theory and indicating that other material properties should be accurately predicted, which we demonstrate for thermal conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Connection formulas for thermal density functional theory

We show that the adiabatic connection formula of ground-state density functional theory relates the correlation energy to a coupling-constant integral over a purely potential contribution, and is widely used to understand and improve approximations. The corresponding formula for thermal density functional theory is cast as an integral over temperatures instead, ranging upward from the system's physical temperature. We also show how to relate different correlation components to each other, either in terms of temperature or coupling-constant integrations. Lastly, we illustrate our results on the uniform electron gas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

On the electronic, structural, and thermodynamic properties of Au supported on α-Fe{sub 2}O{sub 3} surfaces and their interaction with CO

Extensive first principles calculations are carried out to investigate Au monomers and dimers supported on α-Fe{sub 2}O{sub 3}(0001) surfaces in terms of structure optimizations, electronic structure analyses, and ab initio thermodynamics calculations of surface phase diagrams. All computations rely on density functional theory in the generalized gradient approximation (Perdew-Burke-Ernzerhof (PBE)) and account for on-site Coulomb interactions via inclusion of a Hubbard correction (PBE+U). The relative stability of Au monomers/dimers on the stoichiometric termination of α-Fe{sub 2}O{sub 3}(0001) decorated with various vacancies (multiple oxygen vacancies, iron vacancy, and mixed iron-oxygen vacancies) has been computed as a function of the oxygen chemical potential. The charge rearrangement induced by Au at the oxide contact is analyzed in detail and discussed. On one hand, ab initio thermodynamics predicts that under O-rich conditions, structures obtained by replacing a surface Fe atom with a Au atom are thermodynamically stable over a wide range of temperatures. On the other hand, the complex of a CO molecule on a Au atom substituting surface Fe atoms is thermodynamically stable only in a much more narrow range of values of the O chemical potential under O-rich conditions. In the case of a Au dimer, under O-rich conditions, supported Au atoms at an O-Fe di-vacancy are more stable. However, upon CO adsorption, the complex of a CO molecule and 2 Au atoms located at a single Fe vacancy is more favorable.

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

Construction of exchange-correlation functionals through interpolation between the non-interacting and the strong-correlation limit

Drawing on the adiabatic connection of density functional theory, exchange-correlation functionals of Kohn-Sham density functional theory are constructed which interpolate between the extreme limits of the electron-electron interaction strength. The first limit is the non-interacting one, where there is only exchange. The second limit is the strong correlated one, characterized as the minimum of the electron-electron repulsion energy. The exchange-correlation energy in the strong-correlation limit is approximated through a model for the exchange-correlation hole that is referred to as nonlocal-radius model [L. O. Wagner and P. Gori-Giorgi, Phys. Rev. A 90, 052512 (2014)]. Using the non-interacting and strong-correlated extremes, various interpolation schemes are presented that yield new approximations to the adiabatic connection and thus to the exchange-correlation energy. Some of them rely on empiricism while others do not. Several of the proposed approximations yield the exact exchange-correlation energy for one-electron systems where local and semi-local approximations often fail badly. Other proposed approximations generalize existing global hybrids by using a fraction of the exchange-correlation energy in the strong-correlation limit to replace an equal fraction of the semi-local approximation to the exchange-correlation energy in the strong-correlation limit. The performance of the proposed approximations is evaluated for molecular atomization energies, total atomic energies, and ionization potentials.

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

Best recovery of the Laplace operator of a function from incomplete spectral data

This paper is concerned with the problem of best recovery for a fractional power of the Laplacian of a smooth function on R{sup d} from an exact or approximate Fourier transform for it, which is known on some convex subset of R{sup d}. A series of optimal recovery methods is constructed. Information about the Fourier transform outside some ball centred at the origin proves redundant--it is not used by the optimal methods. These optimal methods differ in the way they 'process' key information. Bibliography: 12 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Prediction of a Dynamically Stable New Half-Metallic Phase for the BaN and BaC Compounds

The structural, magnetic, elastic, mechanical, and thermodynamic properties of BaC and BaN compounds in different phases were studied using first-principle calculations based on spin-polarized density functional theory within the generalized gradient approximation (GGA-PBEsol) and the modified Becke–Johnson approach (mBJ-GGA-PBEsol) for the exchange-correlation energy and potential. The following phases—rock-salt (NaCl), CsCl, zinc blende (ZB), NiAs- and WZ-type hexagonal, tetragonal (P4/nmm), and orthorombic (Pnma) phases of BaC and BaN compounds—were considered. We obtained that Pnma phase has the lowest energy configuration as a function of the volume for both the BaN and BaC compounds. The ferromagnetic phase is energetically favored with respect to the non-magnetic phase in the BaN and BaC compounds, except for the CsCl phase in the BaC compound. Considering the phonon dynamics of BaN and BaC compounds in the Pnma, NaCl, ZB, and WZ phases, we observed that the BaN and BaC compounds in the Pnma, NaCl, and ZB phases are dynamically stable. The calculated elastic properties for the Pnma, NaCl, and ZB phases show that they are elastically stable. The Pnma phase for the BaN and BaC compounds, which is a new phase was found to be dynamically and elastically stable. The BaN and BaC compounds exhibit half-metallic behavior in the Pnma, NaCl, and ZB phases. The half-metallic and magnetic character found in the BaN and BaC compounds are attributed to the presence of spin-polarized 2p orbitals of the nitrogen and carbon atoms, respectively. We found that BaN and BaC compounds are half-metallic ferromagnets with magnetic moment of 1 μ{sub B} and 2 μ{sub B} per formula unit, respectively. Using the GGA-PBEsol (mBJ-GGA-PBEsol) approach, our calculated half-metallic gaps for BaN and BaC compounds are 0.22 eV (0.54 eV) and 0.32 eV (0.44 eV) in the Pnma phase, 0.23 eV (1.32 eV) and 0.35 eV (1.00 eV) in the NaCl phase, and 0.38 eV (1.54 eV) and 0.50 eV (1.57 eV) in the ZB phase, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A persistent adjoint method with dynamic time-scaling and an application to mass action kinetics

In this article, we consider an optimization problem where the objective function is evaluated at the fixed-point of a contraction mapping parameterized by a control variable, and optimization takes place over this control variable. Since the derivative of the fixed-point with respect to the parameter can usually not be evaluated exactly, an adjoint dynamical system can be used to estimate gradients. Using this estimation procedure, the optimization algorithm alternates between derivative estimation and an approximate gradient descent step. We analyze a variant of this approach involving dynamic time-scaling, where after each parameter update the adjoint system is iterated until a convergence threshold is passed. Here, we prove that, under certain conditions, the algorithm can find approximate stationary points of the objective function. We demonstrate the approach in the settings of an inverse problem in chemical kinetics, and learning in attractor networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Optimized structure and electronic band gap of monolayer GeSe from quantum Monte Carlo methods

Here, we have used highly accurate quantum Monte Carlo methods to determine the chemical structure and electronic band gaps of monolayer GeSe. Two-dimensional (2D) monolayer GeSe has received a great deal of attention due to its unique thermoelectric, electronic, and optoelectronic properties with a wide range of potential applications. Density functional theory (DFT) methods have usually been applied to obtain optical and structural properties of bulk and 2D GeSe. For the monolayer, DFT typically yields a larger band-gap energy than for bulk GeSe but cannot conclusively determine if the monolayer has a direct or indirect gap. Moreover, the DFT-optimized lattice parameters and atomic coordinates for monolayer GeSe depend strongly on the choice of approximation for the exchange-correlation functional, which makes the ideal structure-and its electronic properties-unclear. In order to obtain accurate lattice parameters and atomic coordinates for the monolayer, we use a surrogate Hessian-based parallel line search within diffusion Monte Carlo to fully optimize the GeSe monolayer structure. The DMC-optimized structure is different from those obtained using DFT, as are calculated band gaps. The potential energy surface has a shallow minimum at the optimal structure. This, combined with the sensitivity of the electronic structure to strain, suggests that the optical properties of monolayer GeSe are highly tunable by strain.

36 MATERIALS SCIENCE↗
DOE OSTI2014

L-asparagine crystals with wide gap semiconductor features: Optical absorption measurements and density functional theory computations

Results of optical absorption measurements are presented together with calculated structural, electronic, and optical properties for the anhydrous monoclinic L-asparagine crystal. Density functional theory (DFT) within the generalized gradient approximation (GGA) including dispersion effects (TS, Grimme) was employed to perform the calculations. The optical absorption measurements revealed that the anhydrous monoclinic L-asparagine crystal is a wide band gap material with 4.95 eV main gap energy. DFT-GGA+TS simulations, on the other hand, produced structural parameters in very good agreement with X-ray data. The lattice parameter differences Δa, Δb, Δc between theory and experiment were as small as 0.020, 0.051, and 0.022 Å, respectively. The calculated band gap energy is smaller than the experimental data by about 15%, with a 4.23 eV indirect band gap corresponding to Z → Γ and Z → β transitions. Three other indirect band gaps of 4.30 eV, 4.32 eV, and 4.36 eV are assigned to α3 → Γ, α1 → Γ, and α2 → Γ transitions, respectively. Δ-sol computations, on the other hand, predict a main band gap of 5.00 eV, just 50 meV above the experimental value. Electronic wavefunctions mainly originating from O 2p–carboxyl, C 2p–side chain, and C 2p–carboxyl orbitals contribute most significantly to the highest valence and lowest conduction energy bands, respectively. By varying the lattice parameters from their converged equilibrium values, we show that the unit cell is less stiff along the b direction than for the a and c directions. Effective mass calculations suggest that hole transport behavior is more anisotropic than electron transport, but the mass values allow for some charge mobility except along a direction perpendicular to the molecular layers of L-asparagine which form the crystal, so anhydrous monoclinic L-asparagine crystals could behave as wide gap semiconductors. Finally, the calculations point to a high degree of optical anisotropy for the absorption and complex dielectric function, with more structured curves for incident light polarized along the 100 and 101 directions.

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

Communication: The description of strong correlation within self-consistent Green's function second-order perturbation theory

We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H{sub 32} finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.

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

Modeling direct interband tunneling. II. Lower-dimensional structures

We investigate the applicability of the two-band Hamiltonian and the widely used Kane analytical formula to interband tunneling along unconfined directions in nanostructures. Through comparisons with k·p and tight-binding calculations and quantum transport simulations, we find that the primary correction is the change in effective band gap. For both constant fields and realistic tunnel field-effect transistors, dimensionally consistent band gap scaling of the Kane formula allows analytical and numerical device simulations to approximate non-equilibrium Green's function current characteristics without arbitrary fitting. This allows efficient first-order calibration of semiclassical models for interband tunneling in nanodevices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗