Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local density 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Pair distribution function analysis of discrete nanomaterials in PDFgui

Pair distribution functions (PDFs) are a leading tool for atomic structure analysis of nanomaterials. However, the most widely used programs for refining atomic structure against PDF data are based on extended crystallographic models, which cannot be applied to discrete, whole nanoparticles. Furthermore, this work describes a straightforward approach to simulate and refine atomistic models of discrete clusters and nanoparticles employing widely used PDF modelling programs such as PDFgui that utilize extended crystallographic models. In this approach, the whole particle to be modelled is contained within an expanded, and otherwise empty, unit cell that is sufficiently large to avoid correlations between atoms in neighbouring unit cells over the r range analyzed. The PDF of the particle is simulated as a composite using two conventional `phases': one that calculates the atom–atom correlations and one that approximates the local number density. This approach is first validated for large nanoparticles that are well modelled by a conventional shape factor model, and then applied to simulate the PDF of discrete particles and low-dimensional materials (graphene and MXene) and to model the experimental PDF data for single-layer FeS nanosheets. A comparison of this approach with the DiffPy-CMI program, which calculates the PDF of discrete species, shows that the composite modelling approach is equally or more accurate. Example input files for implementing this approach within PDFgui and TOPAS, and recommendations for selecting model parameters for reliable application of this refinement strategy, are provided.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Local thermodynamic mapping for effective liquid density-functional theory

The structural-mapping approximation introduced by Lutsko and Baus (1990) in the generalized effective-liquid approximation is extended to include a local thermodynamic mapping based on a spatially dependent effective density for approximating the solid phase in terms of the uniform liquid. This latter approximation, called the local generalized effective-liquid approximation (LGELA) yields excellent predictions for the free energy of hard-sphere solids and for the conditions of coexistence of a hard-sphere fcc solid with a liquid. Moreover, the predicted free energy remains single valued for calculations with more loosely packed crystalline structures, such as the diamond lattice. The spatial dependence of the weighted density makes the LGELA useful in the study of inhomogeneous solids.

Kyrlidis, Agathagelos↗
DOE OSTI2022

Sensitivity of the electronic and magnetic structures of cuprate superconductors to density functional approximations

Abstract We discuss the crystal, electronic, and magnetic structures of La 2− x Sr x CuO 4 (LSCO) for x = 0.0 and x = 0.25 employing 13 density functional approximations, representing the local, semi-local, and hybrid exchange-correlation approximations within the Perdew–Schmidt hierarchy. The meta-generalized gradient approximation (meta-GGA) class of functionals is found to perform well in capturing the key properties of LSCO, a prototypical high-temperature cuprate superconductor. In contrast, the localspin-density approximation, GGA, and the hybrid density functional fail to capture the metal-insulator transition under doping.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Orbital dependent complications for close vs well-separated electrons in diradicals

We investigate two limits in open-shell diradical systems: O3, in which the interesting orbitals are in close proximity to one another, and (C21H13)2, where there is a significant spatial separation between the two orbitals. In accord with earlier calculations, we find that standard density-functional approximations do not predict the open-shell character for the former case but uniformly predict the open-shell character for the latter case. We trace the qualitatively incorrect behavior in O3 predicted by these standard density functional approximations to self-interaction error and use the Fermi–Löwdin-orbital-self-interaction-corrected formalism to determine accurate triplet, closed-shell singlet, and open-shell broken-spin-symmetry electronic configurations. Analysis of the resulting many-electron overlap matrices allows us to unambiguously show that the broken-spin-symmetry configurations do not participate in the representation of the Ms = 0 triplet states and allows us to reliably extract the singlet–triplet splitting in O3 by analyzing the energy as a function of Fermi-orbital-descriptor permutations. The results of these analyses predict the percentage of open-shell character in O3, which agrees well with conventional wavefunction-based methods. While these techniques are expected to be required in cases near the Coulson–Fischer point, we find that they will be less necessary in diradical systems with well-separated electrons, such as (C21H13)2. Results based on energies from self-interaction-corrected generalized gradient, local density, and Hartree–Fock approximations and experimental results are in generally good agreement for O3. These results help form the basis for deriving extended Heisenberg-like Hamiltonians that are needed for descriptions of molecular magnets when there are competing low-energy electronic configurations.

Chemistry↗
DOE OSTI2022

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Deep View of the Large Magellanic Cloud with Six Years of Fermi-LAT Observations

Context. The nearby Large Magellanic Cloud (LMC) provides a rare opportunity of a spatially resolved view of an external star-forming galaxy in gamma-rays. The LMC was detected at 0.1-100 GeV as an extended source with CGRO/EGRET and using early observations with the Fermi-LAT. The emission was found to correlate with massive star-forming regions and to be particularly bright towards 30 Doradus. Aims. Studies of the origin and transport of cosmic rays (CRs) in the Milky Way are frequently hampered by line-of-sight confusion and poor distance determination. The LMC offers a complementary way to address these questions by revealing whether and how the gamma-ray emission is connected to specific objects, populations of objects, and structures in the galaxy. Methods. We revisited the gamma-ray emission from the LMC using about 73 months of Fermi-LAT P7REP data in the 0.2-100 GeV range. We developed a complete spatial and spectral model of the LMC emission, for which we tested several approaches: a simple geometrical description, template-fitting, and a physically driven model for CR-induced interstellar emission. Results. In addition to identifying PSR J0540-6919 through its pulsations, we find two hard sources positionally coincident with plerion N 157B and supernova remnant N 132D, which were also detected at TeV energies with H.E.S.S. We detect an additional soft source that is currently unidentified. Extended emission dominates the total flux from the LMC. It consists of an extended component of about the size of the galaxy and additional emission from three to four regions with degree-scale sizes. If it is interpreted as CRs interacting with interstellar gas, the large-scale emission implies a large-scale population of approximately 1-100 GeV CRs with a density of approximately 30% of the local Galactic value. On top of that, the three to four small-scale emission regions would correspond to enhancements of the CR density by factors 2 to 6 or higher, possibly more energetic and younger populations of CRs compared to the large-scale population. An alternative explanation is that this is emission from an unresolved population of at least two dozen objects, such as pulsars and their nebulae or supernova remnants. This small-scale extended emission has a spatial distribution that does not clearly correlate with known components of the LMC, except for a possible relation to cavities and supergiant shells. Conclusions. The Fermi-LAT GeV observations allowed us to detect individual sources in the LMC. Three of the newly discovered sources are associated with rare and extreme objects. The 30 Doradus region is prominent in GeV gamma-rays because PSR J0540-6919 and N 157B are strong emitters. The extended emission from the galaxy has an unexpected spatial distribution, and observations at higher energies and in radio may help to clarify its origin.

gamma rays: galaxies – Magellanic Clouds – cos↗
DOE OSTI2022

Experimental distinction of the molecularly induced Balmer emission contribution and its application for inferring molecular divertor density with 2D filtered camera measurements during detachment in JET L-mode plasmas

A previously presented model for generating 2D estimates of the divertor plasma conditions at JET from deuterium Balmer line intensity ratios, obtained from tomographic reconstructions of divertor camera images, was amended to consider also the Balmer emission arising from molecular processes. Utilizing the AMJUEL and H2VIBR atomic and molecular databases of EIRENE enabled also inference of the molecular divertor density from the distinguished molecularly induced emission. Analysis of a JET L-mode density scan suggests the molecularly induced emission accounting for up to 60%–70% and 10%–20% of the Balmer Dα and Dγ intensities, respectively, at the onset of detachment, while electron-ion recombination becomes increasingly dominant with deepening detachment. Similar observations were made by post-processing EDGE2D-EIRENE simulations, which indicated significant roles of molecular D$^{+}_{2}$ ions and vibrational excitation of the D 2 molecules as precursors for the molecularly induced emission. The experimentally inferred molecular density at the outer strike point was found to increase monotonously with decreasing strike point temperature, reaching approximately 30%–50% of the local electron density at n mol, osp = 1–2 X 10 20 m -3 at T e,osp ≈ 0.7 eV. A further steep increase by a factor of 3–5 was observed with decrease of T e,osp to 0.5 eV. The observations are in qualitative and reasonable quantitative agreement with EDGE2D-EIRENE predictions of n mol, osp within the uncertainties of the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Linear density drift instabilities in very low beta plasmas A different approach

This paper reports a study of the linear Vlasov electromagnetic dispersion equation for density drift instabilities in very low beta plasmas. A uniform magnetic field and a weak uniform density gradient are assumed. This paper differs from most other studies of this topic in three important ways: (1) no low frequency or long wavelength approximations are made, (2) no gauge condition is imposed and (3) a modification of the local approximation is used which is argued to be less arbitrary than the usual local approximation. Numerical solution of the resulting dispersion equation yields stabilization of the universal density drift instability at values of beta somewhat higher than those obtained from the local approximation, and near maximum growth shows no evidence of coupling between the universal and Alfven modes.

Gary, S. P.↗
DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1977

Numerical simulation of spacecraft charging phenomena

A numerical simulation program is being constructed having the following features: (1) infinite circular cylindrical geometry with angle-dependence, (2) inclusion of incident particles, photoelectrons, secondary electrons, backscattered electrons, any gun emissions, and any internal current pathways including surface conductive layers, (3) quasistatic time-dependent iteration, in which sheath potential changes during particle transit times are ignored, (4) use of approximate, locally-dependent space charge density expressions in solving Poisson's equation for sheath potentials, with use of numerical orbit-following to determine surface currents, and (5) incident particle velocity distributions isotropic or beam-like, or some superposition of these. Rationales for each of these features are discussed.

Laframboise, J. G.↗
DOE OSTI2021

How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations?

Herein, we examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys. , 2018, 20 , 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew–Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi–Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys. , 2019, 151 , 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.

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

Characterization of core neutrals using passive measurements of the D α spectrum near the X -point in the DIII-D tokamak

New spectroscopic measurements of deuterium Balmer-α emission are used to infer the spatial distribution of neutral particles near the X-point in a diverted high confinement mode plasma. The charge exchange neutral spectroscopy (CENS) diagnostic—recently installed on the DIII-D tokamak—uses 15 lines of sight extending from the edge of the confined region (ρ ≈ 0.7) to the X-point in lower single-null plasmas. Each CENS line of sight is spectrally resolved around the deuterium Balmer-α wavelength (6561 Å) to measure the Doppler shift, line broadening, and intensity of Dα emission. Thermal line broadening is used to identify emission from neutrals undergoing charge-exchange (CX) with high temperature ions in the confined plasma. This spectral information allows a more accurate determination of the neutral density deeper inside the confined plasma compared to traditional filter-based diagnostics. We present two methods of analyzing CENS measurements: (1) fitting the spectrum from each view-chord separately, and (2) a tomographic method for inverting the neutral density over a 2D region of space using all CENS views collectively. The neutral density profile is found to decay exponentially in the radial direction with two scale-lengths. In the pedestal region the neutral density decays at a rate approximately equal to the local mean-free-path for CX collisions, $L_{n_\mathrm{D0}}\approx \lambda_\text{CX}$. Further inside the plasma the neutral density decays at a rate equal to the mean-free-path for ionization, $L_{n_\mathrm{D0}}\approx \lambda_\text{inz.}$. The separatrix $n_\mathrm{D0}$ value is found to be approximately $2\times10^{15}$ m−3. Based on results from the 2D inversion, the density of neutrals is found to increase along the separatrix approaching the X-point by a factor of $\approx 5\times$ over the region covered by the CENS diagnostic.

X-point neutrals↗
NASA NTRS1979

Sunspots and the physics of magnetic flux tubes. I - The general nature of the sunspot. II - Aerodynamic drag

Analysis of the dynamical stability of a large flux tube suggests that the field of a sunspot must divide into many separate tubes within the first 1000 km below the surface. Buoyancy of the Wilson depression at the visible surface and probably also a downdraft beneath the sunspot hold the separate tubes in a loose cluster. Convective generation of Alfven waves, which are emitted preferentially downward, cools the tubes. Aerodynamic drag on a slender flux tube stretched vertically across a convective cell is also studied. Since the drag is approximately proportional to the local kinetic energy density, the density stratification weights the drag in favor of the upper layers. Horizontal motions concentrated in the bottom of the convective cell may reverse this density effect. A downdraft of about two km/sec through the flux tubes beneath the sunspot is hypothesized.

Parker, E. N.↗
NASA NTRS1983

Space radiation interaction mechanisms in materials

Models of charged-particle impact under conditions typical of the space environment are reported, with a focus on impact excitation and nuclear reactions, especially for heavy ions. Impact excitation is studied by using a global model for electronic excitation based on formal relations through the classical dielectric function to derive an approximation related to the local plasma (electron density distribution) within the atoms and molecules; and corrections to the model resulting from the nonfluid nature of this plasma are discussed. Nuclear reactions are studied by reducing quantum-mechanical treatment of this general N-body problem to an equivalent two-body problem that is solvable, and by comparing the results with experimental data. The equations for heavy-charged-particle transport are derived and solution techniques demonstrated. Finally, these methods of analysis are applied to study the change in the electrical properties of a GaAs semiconductor for photovoltaic applications; proton damage to GaAs crystals is found to arise from stable replacement defects and to be nonannealable, in contrast to electron-induced damage.

Wilson, J. W.↗
NASA NTRS2015

Modeling the Swift BAT Trigger Algorithm with Machine Learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. (2014) is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of approximately greater than 97% (approximately less than 3% error), which is a significant improvement on a cut in GRB flux which has an accuracy of 89:6% (10:4% error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of eta(sub 0) approximately 0.48(+0.41/-0.23) Gpc(exp -3) yr(exp -1) with power-law indices of eta(sub 1) approximately 1.7(+0.6/-0.5) and eta(sub 2) approximately -5.9(+5.7/-0.1) for GRBs above and below a break point of z(sub 1) approximately 6.8(+2.8/-3.2). This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online.

gamma rays: general↗
NASA NTRS1985

Geocoronal structure - The effects of solar radiation pressure and the plasmasphere interaction

The theory of planetary exospheres is extended to incorporate solar radiation pressure in a rigorous manner, and an evaporative geocoronal prototype (classical, motionless exobase) is constructed using Liouville's theorem. Model calculations for density and kinetic temperature at points along the earth-sun axis (solar and antisolar directions) reveal an extensive satellite component, comprising approximately 2/3 of the total hydrogen density near 10 earth radii, and a temperature profile suggestive of an isotropic quasi-Maxwellian velocity distribution for the bound component. A geotail is also evident as an enhancement of the density at local midnight compared to local noon that increases outward (from approximately 25 percent at 10 earth radii to over 60 percent at 20 earth radii). Additional mechanisms acting upon the geocorona alter the basic evaporative case in notable ways. Solar ionization has been included in a simple fashion; the effect is to partially deplete the density without otherwise altering the structure. Interaction with a simple plasmasphere via the Boltzmann equation results in 'heating' the geocorona and enhancing the escape flux at the expense of the density of the bound component, an effect not appreciated in earlier studies; the geotail survives this interaction.

Bishop, J.↗
DOE OSTI2024

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dynamically Collected Local Density using Low-Cost Lidar and its Application to Traffic Models

This article demonstrates the use of traffic density observations collected dynamically in the vicinity of probe vehicles. Fixed position sensors cannot capture the longitudinal evolution of local traffic density in the corridor. In this research, dynamic traffic density observations were collected in a naturalistic driving setting that was free of any controlled experiment biases. Speed from global positioning system and space headway from a light detection and ranging module was collected on one arterial and one freeway segment, 2 and 4mi long, respectively. The combined data frequency was approximately 3Hz. Space headway was used to estimate the local density and consequently to identify the density of a specific location in a corridor. Besides, driver behavior was characterized using the relationship between instantaneous speed and local density under different regimes of the Wiedemann car-following model. Macroscopic traffic stream models were used to investigate the relationship between dynamically collected instantaneous speed and local density. Using the longitudinal evolution of density, precise local density across the corridor can be obtained along with the leader and follower trajectories. A method to identify driver behavior across density ranges was developed for different facility types using a microscopic relationship between instantaneous speed and local density. Overall driving behavior on the freeway segment can be represented by translating the instantaneous speed and local density relationship to macroscopic stream models.

Engineering↗