Search NASA⌕ Search

SEARCH · Search NASA

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Atomic ionization: sd energy imbalance and Perdew–Zunger self-interaction correction energy penalty in 3d atoms

To accurately describe the energetics of transition metal systems, density functional approximations (DFAs) must provide a balanced description of s- and d- electrons. One measure of this is the sd transfer error, which has previously been defined as E ( 3 d n − 1 4 s 1 ) − E ( 3 d n − 2 4 s 2 ) . Theoretical concerns have been raised about this definition due to its evaluation of excited-state energies using ground-state DFAs. A more serious concern appears to be strong correlation in the 4s 2 configuration. Here, we define a ground-state measure of the sd energy imbalance, based on the errors of s- and d-electron second ionization energies of the 3d atoms, that effectively circumvents the aforementioned problems. We find an improved performance as we move from the local spin density approximation (LSDA) to the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) to the regularized and restored Strongly Constrained and Appropriately Normed (r 2 SCAN) meta-GGA for first-row transition metal atoms. However, we find large (∼2 eV) ground-state sd energy imbalances when applying a Perdew–Zunger 1981 self-interaction correction. This is attributed to an “energy penalty” associated with the noded 3d orbitals. A local scaling of the self-interaction correction to LSDA results in a balance of s- and d-errors.

Science & Technology - Other Topics↗
DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental and computational studies on high-entropy carbide MoNbTaVWC 5 under high pressures

High-entropy carbide, MoNbTaVWC 5 , was synthesized from oxide precursors of the constituent metals, mixed with graphite powder in a microwave-generated hydrogen plasma at 26.66 kPa and 2100 °C. Ambient x-ray diffraction analysis confirms the full conversion of oxide precursors into a single-phase, face-centered cubic structure with a lattice parameter a = 4.3309 Å. Nanoindentation measured a hardness of 24.5 ± 1.3 GPa and an elastic modulus of 386 ± 22 GPa. The synthesized sample, mixed with a copper pressure marker, was studied by the radial x-ray diffraction technique with beryllium gasketing in a diamond anvil cell up to 70 GPa. The experimentally measured pressure–volume curve and shear strength were compared with theoretical predictions using the special quasi-random structure technique and density functional theory. MoNbTaVWC 5 achieved a 12% volume compression at 70 GPa and exhibited a high shear strength of 6.6 GPa. The present study demonstrates that the high-entropy carbide MoNbTaVWC 5 exhibits exceptional incompressibility and high strength under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The electronic structure, crystal fields, and magnetic anisotropy in RECo 5 magnets

The current progress in describing rare-earth-based magnets' electronic structure and magnetic properties is discussed. We use several currently popular electronic structure methods to show the typical values of critical parameters that define the physics of RECo 5 (RE = rare earth atom) materials. The magnetic moments and magnetic anisotropy of 4f atoms are obtained using several approaches, including anisotropic 4f-charge density-constrained DFT and DFT+HI methods. We also suggest the introduction of "penalty" functional for obtaining correct variational total energy in the traditional Hund's rule-constrained DFT-based techniques. The applicability and future extensions are discussed. The proposed combination of methods is potentially suitable for high-throughput computational searches of new rare-earth-containing magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

EGRET Observations of the Diffuse Gamma-Ray Emission in Orion: Analysis Through Cycle 6

We present a study of the high-energy diffuse emission observed toward Orion by the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray Observatory. The total exposure by EGRET in this region has increased by more than a factor of two since a previous study. A simple model for the diffuse emission adequately fits the data; no significant point sources are detected in the region studied (1 = 195 deg to 220 deg and b = -25 deg to -10 deg) in either the composite dataset or in two separate groups of EGRET viewing periods considered. The gamma-ray emissivity in Orion is found to be (1.65 +/- 0.11) x 10(exp -26)/s.sr for E > 100 MeV, and the differential emissivity is well-described as a combination of contributions from cosmic-ray electrons and protons with approximately the local density. The molecular mass calibrating ratio is N(H2)/W(sub CO) = (1.35 +/- 0.15) x 10(exp 20)/sq cm.(K.km/s).

Digel, S. W.↗
NASA NTRS1982

Stellar contributions to the hard X-ray galactic ridge

The number density of serendipitous sources in galactic plane Einstein Observatory IPC fields are compared with predictions based on the intensity of the HEAO-1 A2 unresolved hrd X-ray galactic ridge emission. It is concluded that theoretically predicted X-ray source populations of luminosity 8 x 10 to the 32nd power to 3 x 10 to the 34th power ergs s have 2 KeV to 10 KeV local surface densities of less than approximately .0008 L(32) pc/2 and are unlikely to be the dominant contributors to the hard X-ray ridge. An estimate for Be/neutron star binary systems, such as X Persei, gives a 2 keV to 10 keV local surface density of approximately 26 x 10 to the -5 power L(32) pc/2. Stellar systems of low luminosity, are more likely contributors. Both RS CVn and cataclysmic variable systems contribute 43% + or - 18% of the ridge. A more sensitive measurement of the ridge's hard X-ray spectrum should reveal Fe-line emission. We speculate that dM stars are further major contributors.

Worrall, S. M.↗
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↗
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↗
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.↗
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.↗
NASA NTRS1982

An electrostatic parabolic density drift instability

Assuming a uniform magnetic field and using the local approximation, the Vlassov theory for an electrostatic instability driven by a parabolic density gradient is considered. It is found that the instability grows only in the case of a sufficiently large, positive second derivative of the density. Linear growth rate parametric dependences are given, in addition to weakly nonlinear calculations on wave-particle transport due to this instability.

Gary, S. P.↗
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.↗
NASA NTRS1980

Primordial nucleosynthesis and Dirac's large numbers hypothesis

Consideration is given to the analysis of Falik (1979) which attempted to show that the cosmological model proposed by Canuto and Hsieh (1978) in which the gravitational constant varies with time contradicts observations of primordial helium. It is shown that the analysis was based on the assumptions that (1) the energy density of radiation in local thermodynamic equilibrium is approximately equal to the fourth power of the equilibrium temperature, where the product of the equilibrium temperature with the scale factor of the Robertson-Walker metric is constant, and (2) the gravitational constant is approximately equal to the inverse of the time even at early cosmological epochs. These assumptions are demonstrated to be invalid in the scale covariant theory of gravitation used to develop the model, thus negating the conclusion that the Canuto and Hsieh model excludes the primordial synthesis of helium.

Canuto, V.↗