Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional approximations”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Layer-dependent topological phase in a two-dimensional quasicrystal and approximant

The electronic and topological properties of materials are derived from the interplay between crystalline symmetry and dimensionality. Simultaneously introducing “forbidden” symmetries via quasiperiodic ordering with low dimensionality into a material system promises the emergence of new physical phenomena. Here, we isolate a two-dimensional (2D) chalcogenide quasicrystal and approximant, and investigate their electronic and topological properties. The 2D layers of the materials with a composition close to Ta 1.6 Te, derived from a layered transition metal dichalcogenide, are isolated with standard exfoliation techniques, and investigated with electron diffraction and atomic resolution scanning transmission electron microscopy. Density functional theory calculations and symmetry analysis of the large unit cell crystalline approximant of the quasicrystal, Ta 21 Te 13 , reveal the presence of symmetry-protected nodal crossings in the quasicrystalline and approximant phases, whose presence is tunable by layer number. Our study provides a platform for the exploration of physics in quasicrystalline, low-dimensional materials and the interconnected nature of topology, dimensionality, and symmetry in electronic systems.

36 MATERIALS SCIENCE↗
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 OSTI2016

Density functional theory investigation of the LiIn 1- x Ga x Se 2 solid solution

Here, the electronic structure and optical properties of the LiIn 1-x Ga x Se 2 (x=0, 0.25, 0.5, 0.75, 1) solid solution were studied by density functional theory (DFT) with pure functionals. The exchange-correlation is treated within the local density approximation (LDA) and generalized-gradient approximation (GGA). The electronic structures for each respective compound are discussed in detail. Calculations reveal that gallium incorporation can be used to tune the optical-electrical properties of the solid solution and correlates with the lattice parameter. The band gap trend of the LiIn 1-x Ga x Se 2 system follows a nonlinear behavior between the LiInSe 2 and LiGaSe 2 ternary boundaries. The bowing parameter is estimated to be on the order of 0.1- 0.3 eV at the point. Low-temperature optical absorption revealed a 30% change in the temperature dependence of the band gap for the intermediate compound LiIn 0.6 Ga 0.4 Se 2 compared to ternary boundaries and suggests the heat capacity to be another control element through strain.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multiplicity dependence of light-flavor hadron production in p p collisions at s = 7 TeV

Comprehensive results on the production of unidentified charged particles, $π^±, K^±, K^0_S, K^*(892)^0, p,$ $\bar{p}$, $\phi$ (1020), Λ, $\bar{Λ}$, $\Xi^-$, $\bar{\Xi}^+$, $Ω^-$, and $\bar{Ω}^+$ hadrons in proton-proton ($pp$) collisions at $\sqrt{s}$ = 7 TeV at midrapidity ($|y|$ < 0.5) as a function of charged-particle multiplicity density are presented. In order to avoid autocorrelation biases, the actual transverse momentum ($p$ T ) spectra of the particles under study and the event activity are measured in different rapidity windows. In the highest multiplicity class, the charged-particle density reaches about 3.5 times the value measured in inelastic collisions. While the yield of protons normalized to pions remains approximately constant as a function of multiplicity, the corresponding ratios of strange hadrons to pions show a significant enhancement that increases with increasing strangeness content. Furthermore, all identified particle-to-pion ratios are shown to depend solely on charged-particle multiplicity density, regardless of system type and collision energy. The evolution of the spectral shapes with multiplicity and hadron mass shows patterns that are similar to those observed in $p$-Pb and Pb-Pb collisions at Large Hadron Collider energies. The obtained $p$ T distributions and yields are compared to expectations from QCD-based $pp$ event generators as well as to predictions from thermal and hydrodynamic models. These comparisons indicate that traces of a collective, equilibrated system are already present in high-multiplicity $pp$ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Effect of Impurities on the Formation Energy of Point Defects in the γ-TiAl Alloy

The projector augmented-wave method within density functional theory is used to study the effect of 4d elements and Group IIIA and IVA elements on the energetics of point defect formation and elastic moduli of the γ-TiAl alloy. To calculate the point defect formation energy, the grand canonical formalism was used. It is shown that the formation energy of aluminum vacancies decreases by approximately 1.3 eV with increasing its content in the series of alloys Ti{sub 3}Al–TiAl–TiAl{sub 3}, whereas the formation energy of titanium vacancies changes insignificantly. In general, the formation of titanium vacancies in these alloys is preferred to the formation of aluminum vacancies. It has been found that Nb, Mo, Tc, Ru, Rh, and Pd impurities on the aluminum sublattice contribute to an increase in the formation energy of an aluminum vacancies, and Mo and Tc also lead to an increase in the formation energy of titanium vacancies. All 4d impurities, if they substitute for titanium, reduce the formation energy of aluminum vacancies, and Nb and Mo are favorable for increasing the formation energy of titanium vacancies. The influence of impurities on the chemical bond in the γ-TiAl alloy and its elastic moduli and characteristics based on them is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Multicomponent mutual diffusion in the warm, dense matter regime

In this work, we present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential. Onsager coefficients then connect the macroscopic diffusion to microscopic particle motions, evinced in trajectories characterized by positions and velocities, through various autocorrelation functions (ACFs). These trajectories are generated by molecular dynamics (MD) simulations either through the Born-Oppenheimer approximation, which treats the ions classically and the electrons quantum-mechanically by an orbital-free density-functional theory, or through a classical MD approach with Yukawa pair-potentials, whose effective ionizations and electron screening length derive from quantal considerations. We employ the reference-mean form of the ACFs and determine the center-of-mass coefficients through a simple reference-frame-dependent similarity transformation. The Onsager terms in turn determine the mutual diffusion coefficients. We examine a representative sample of ternary mixtures as a function of density and temperature from those with only light elements (D-Li-C, D-Li-Al) to those with highly asymmetric mass components (D-Li-Cu, D-Li-Ag, H-C-Ag). We also follow trends in the diffusion as a function of number concentration and evaluated the efficacy of various approximations such as the Darken approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Sharp edges to neutral hydrogen disks in galaxies and the extragalactic radiation field

It is shown that the very sharp truncation of the neutral hydrogen distribution seen in NGC 3198 (and probably M33) is well modeled as the result of ionization of the atomic gas by the extragalactic radiation field. Below a critical column density of about a few times 10 exp 19/sq cm the gas is dominantly ionized and undetectable in the 21-cm line. It is inferred from the photoionization models that the total disk gas distribution in NGC 3198 is actually fairly axisymmetric. The critical column density for ionization is not a strong function of galaxy mass or mass distribution; thus, all galaxies should show a cutoff at approximately the same column density. Specific models of 3198 suggest that the extragalactic ionizing photon flux is 5000-10,000 photons/sq cm s.

Maloney, Philip↗
DOE OSTI2021

Electronic properties of single-layer CoO 2 /Au(111)

Here we report direct measurements via angle-resolved photoemission spectroscopy (ARPES) of the electronic dispersion of single-layer (SL) CoO 2 . The Fermi contour consists of a large hole pocket centered at the $\overline{\Gamma}$ point. To interpret the ARPES results, we use density functional theory (DFT) in combination with the multi-orbital Gutzwiller Approximation (DFT+GA), basing our calculations on crystalline structure parameters derived from x-ray photoelectron diffraction and low-energy electron diffraction. Our calculations are in good agreement with the measured dispersion. We conclude that the material is a moderately correlated metal. We also discuss substrate effects, and the influence of hydroxylation on the CoO 2 SL electronic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Path-Integral Monte Carlo Simulation of the Warm Dense Homogeneous Electron Gas

We perform calculations of the 3D finite-temperature homogeneous electron gas in the warm-dense regime $(r_s ≡ (3/4πn)^{1/3}a^{-1}_0 = 1.0 – 40.0$ and $Θ ≡ T/T_F = 0.0625 – 8.0)$ using restricted path-integral Monte Carlo simulations. Precise energies, pair correlation functions, and structure factors are obtained. For all densities, we find a significant discrepancy between the ground state parametrized local density approximation and our results around $T_F$. These results can be used as a benchmark for developing finite-temperature density functionals, as well as input for orbital-free density function theory formulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Lastly, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

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

Implementation of a generalized coherent elastic scattering formulation for thermal neutron scattering analysis - 14487

The scattering interactions of thermal neutrons in matter are classified into inelastic and elastic components. In crystalline materials, coherent elastic scattering represents a dominant form of interaction. Traditionally, the treatment of this interaction in standard computer codes (e.g., the LEAPR module of the NJOY code system) was performed within the cubic approximation and restricted the analyst to a few materials that were hardwired into the code. In general, these restrictions can hinder the production of ENDF style thermal neutron scattering libraries for new materials outside the current portfolio of treated materials. In this work, a general formulation based on the exact theory of coherent elastic scattering was developed. This formulation is capable of processing the crystalline structure of any material. The cubic approximation is relaxed by accounting for the partial forces from asymmetric crystal structures in the calculation. The atom sites approximation is relaxed by utilizing a site-wise Debye-Waller factor. Specifically, the Debye-Waller factor is no longer regarded as a uniform constant across the crystal but treated as a matrix tensor that uses the integral of the partial phonon density of states. Therefore, the calculation of the coherent elastic scattering cross section is executed by applying the matrix-form of the Debye-Waller factor on each reciprocal space point and summing them up. Furthermore, the developed algorithm supports compounds with any number of elements and of any ratio to each other. This general formulation was successfully implemented and contributed to the production of ENDF style thermal neutron scattering libraries for silicon dioxide (α-SiO{sub 2}) and silicon carbide (3C-SiC). The needed input for the analysis was generated using density functional theory (DFT) atomistic simulations and dynamical matrix calculations for both materials. (authors)

42 ENGINEERING↗
DOE OSTI2017

Competing magnetic instabilities in the weak itinerant antiferromagnetic TiAu

The binary intermetallic compound TiAu is distinguished by an extremely sharp and narrow density of states peak at the Fermi level that has been proposed, via scattering between mirrored van Hove singularities, to be the mechanism for antiferromagnetic ordering versus the more fundamental ferromagnetic Stoner instability. Here we study, using density functional theory methods, magnetic tendencies and effects of doping, the latter within the virtual crystal approximation (VCA). Ferromagnetic tendencies are quantified using the fixed spin moment approach, illustrating the strong Stoner instability that does not, however, provide the ground state. Use of VCA results allows the identification of the value of the Stoner exchange constant I = 0.74eV for Ti. Magnetic fluctuations not included with semilocal density functionals are quantified with the procedure provided by Ortenzi and coauthors, with alloy concentrations corresponding to the quantum critical points reduced by a factor of three to five. Furthermore, our results provide useful guidelines for experimental doping studies of TiAu.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Fully coupled simulation of cosmic reionization. II. Recombinations, clumping factors, and the photon budget for reionization

We use a fully coupled cosmological simulation including dark matter dynamics, multispecies hydrodynamics, nonequilibrium chemical ionization, flux-limited diffusion radiation transport, and a parameterized model of star formation and feedback (thermal and radiative) to investigate the epoch of hydrogen reionization in detail. In this paper, the first of several application papers, we investigate the mechanics of reionization from stellar sources forming in high-z galaxies, the utility of various formulations for the gas clumping factor on accurately estimating the effective recombination time in the intergalactic medium (IGM), and the photon budget required to achieve reionization. We also test the accuracy of the static and time-dependent models of Madau et al. as predictors of reionization completion/maintenance. We simulate a WMAP7 ΛCDM cosmological model in a 20 comoving Mpc cube, resolved with 800{sup 3} uniform fluid cells and dark matter particles. By tuning our star formation recipe to approximately match the observed high-redshift star formation rate density and galaxy luminosity function, we have created a fully coupled radiation hydrodynamical realization of hydrogen reionization, which begins to ionize at z ≈ 10 and is completed at z ≈ 5.8 without further tuning. We find that roughly two ionizing photons per H atom are required to convert the neutral IGM to a highly ionized state. After reionization concludes, we find that the quantity n-dot{sub ion}×(1 Gyr)/n{sub H} is ∼9 at z = 5, in rough agreement with measurements of the ionizing emissivity by Becker and Bolton. The complicated events during reionization that lead to this number can be generally described as inside-out, but in reality, the narrative depends on the level of ionization of the gas one attributes as being ionized. We find that the formula for the ionizing photon production rate needed to maintain the IGM in an ionized state derived by Madau et al. should not be used to predict the epoch of reionization completion because it ignores history-dependent terms in the global ionization balance which are not ignorable. We find that the time-dependent model for the ionized volume fraction Q{sub HII} is more predictive, but overestimates the redshift of reionization completion z{sub reion} by Δz ≈ 1. We propose a revised formulation of the time-dependent model that agrees with our simulation to a few percent accuracy. Finally, we use our simulation data to measure the absorption of ionizing radiation due to circumgalactic gas resolved on our mesh and find f-bar{sub esc}(CGM)≈0.7.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Continuation of probability density functions using a generalized Lyapunov approach

Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the efficient solution of a generalized Lyapunov equation using an iterative method involving low-rank approximations. We apply and illustrate the capabilities of the method using a problem in physical oceanography, i.e. the occurrence of multiple steady states of the Atlantic Ocean circulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Correlation effects on spin-polarized electron-hole quantum bilayer

We present a numerical calculation for the intra- and interlayer pair-correlation functions, g{sub ll’}(r), of spin-polarized electron-hole quantum bilayers at zero temperature. The calculations of g{sub ll’}(r) are performed by including electron correlations within the dynamical version of the self-consistent mean-field approximation of Singwi, Tosi, Land and Sjölander (qSTLS). Our study reveals that the critical layer density decreases (increases) due to the inclusion of finite width (mass-asymmetry) effect during the phase-transition from charge-density wave to Wigner crystal ground-state by yielding the pronounced oscillatory behavior ing{sub ll}(r). The results are compared with recent findings of spin-polarized electron-hole quantum bilayers with mass-symmetry and zero width effects. To highlight the importance of dynamical character of correlations, we have also compared our results with the STLS results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Microscopic Approach

Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field.

Li, Jianzhong↗
NASA NTRS2002

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗
DOE OSTI2014

The adsorption of α-cyanoacrylic acid on anatase TiO{sub 2} (101) and (001) surfaces: A density functional theory study

The adsorption of α-cyanoacrylic acid (CAA) on anatase TiO{sub 2} (101) and (001) surfaces, including adsorption energies, structures, and electronic properties, have been studied by means of density functional theory calculations in connection with ultrasoft pseudopotential and generalized gradient approximation based upon slab models. The most stable structure of CAA on anatase TiO{sub 2} (101) surface is the dissociated bidentate configuration where the cyano N and carbonyl O bond with two adjacent surface Ti atoms along [010] direction and the dissociated H binds to the surface bridging O which connects the surface Ti bonded with carbonyl O. While for the adsorption of CAA on (001) surface, the most stable structure is the bidentate configuration through the dissociation of hydroxyl in carboxyl moiety. The O atoms of carboxyl bond with two neighbor surface Ti along [100] direction, and the H from dissociated hydroxyl interacts with surface bridging O, generating OH species. The adsorption energies are estimated to be 1.02 and 3.25 eV for (101) and (001) surfaces, respectively. The analysis of density of states not only suggests the bonds between CAA and TiO{sub 2} surfaces are formed but also indicates that CAA adsorptions on TiO{sub 2} (101) and (001) surfaces provide feasible mode for photo-induced electron injection through the interface between TiO{sub 2} and CAA. This is resulted from that, compared with the contribution of CAA orbitals in valence bands, the conduction bands which are mainly composed of Ti 3d orbitals have remarkable reduction of the component of CAA orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗