SEARCH · Engineering Papers

Results for “Density functional”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Self-interaction-corrected Kohn–Sham effective potentials using the density-consistent effective potential method

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1992

Effect of fiber orientation on the fracture toughness of brittle matrix composites

The effective fracture toughness of brittle matrix materials can be increased through the addition of short, poorly bonded fibers which bridge the growing crack. The orientation distribution of the fibers is likely to be biased, and not in an ideal random or aligned state. A micromechanical model is formulated for the postcracking behavior using the force-displacement relation for an arbitrary fiber bridging a crack, the fiber orientation density function, and the fiber location density function. This model is then used to determine an effective traction law for the bridging fibers, as well as the steady state bridging toughness increment. In most cases, the results may be placed in the form of a product of the aligned fiber results times a modifying integrated orientation factor. The frictional shear stress on fiber pull-out is allowed to vary during pull-out, modeling the effects of matrix breakdown, fiber surface smoothing or wear debris accumulation. Results are presented for a variety of representative planar and three-dimensional fiber orientation states.

Jain, L. K.↗

NASA NTRS1992

Robust control algorithms for Mars aerobraking

Four atmospheric guidance concepts have been adapted to control an interplanetary vehicle aerobraking in the Martian atmosphere. The first two offer improvements to the Analytic Predictor Corrector (APC) to increase its robustness to density variations. The second two are variations of a new Liapunov tracking exit phase algorithm, developed to guide the vehicle along a reference trajectory. These four new controllers are tested using a six degree of freedom computer simulation to evaluate their robustness. MARSGRAM is used to develop realistic atmospheres for the study. When square wave density pulses perturb the atmosphere all four controllers are successful. The algorithms are tested against atmospheres where the inbound and outbound density functions are different. Square wave density pulses are again used, but only for the outbound leg of the trajectory. Additionally, sine waves are used to perturb the density function. The new algorithms are found to be more robust than any previously tested and a Liapunov controller is selected as the most robust control algorithm overall examined.

Shipley, Buford W., Jr.↗

NASA NTRS1995

A LEO concept for millimeter wave satellite communication

A conceptual 60 satellite LEO constellation for millimeter wave communication is discussed. It could be launched in segments, with the first 30 satellites providing high elevation angles for all time in the Northern latitudes between Miami and Thule. The second set of 30 satellites would complete the worldwide coverage with emphasis on high ground elevation angles in the densely populated temperate zones. Full earth searches for all time are used to generate probability density functions for elevation angle. The density functions are used to derive optimum frequencies for random elevation systems. The 55 degree average elevation angle and 14 degree standard deviation are seen to be acceptable for 0.997 rain availability in Washington, DC for the 40 to 47 GHz region. The 40 to 47 GHz region is nearly optimum, if 0.99 rain availability is acceptable.

Jackson, A. H.↗

NASA NTRS2001

Validation of TRMM Satellite Rainfall Products Over Oklahoma for a Three Year Period (1998-2000)

This study evaluates the level 2 and level 3 rainfall products from the Tropical Rainfall Measuring Mission gridded distribution of 118 high resolution tipping buckets from the Oklahoma Mesonet. The Tropical Rainfall Measuring Mission, a joint satellite mission between NASA and the National Space Development Agency (NASDA) of Japan, was designed to estimate global precipitation between 40 S and 40 N latitude. The TRMM satellite, consisting of two main precipitation sensors, a passive microwave (TMI) and precipitation radar (PR) sensors, was launched in 1997. Although the great advantage of TRMM has been its ability to sample precipitation globally over the tropical oceans in places where ground sensors do not exist, regional points over land still offer the best opportunity for validating these rain estimates. For this study, the gauge data was gridded at I x I degree resolution between 34 and 36 N and -95 and -100 W. The location of Oklahoma is somewhat unique, in that, it is located near the turning point for the satellite. This study investigates effects of temporal sampling on the satellite measurements and the resulting rainfall bias observed between space and ground based sensors. The first part of the study analyzes the monthly rainfall statistics of the 3AI2 (TMI), the 3A25 (PR) and the 31142 (combined sensors) with the rain gauges. It uses regression techniques and a computation of the probability density function (PDF) and cumulative probability density function (CDF) for each sensor to evaluate the effects of temporal sampling on the satellite products and the sufficiency of using two overpasses a day to estimate rainfall on a monthly scale. In further probing this issue, the study also looked at the direct coincidence of ground and satellite measurements by comparing instantaneous comparisons obtained from the level 2 orbital track data, e.g. 2A12, 2A25. Lastly, the effects of temporal sampling were further studied by sub-sampling the gauge data only at overpass times. This methodology provides a quantitative way of separately inferring the fraction of bias due to temporal sampling and the rainfall algorithm. In addition to showing the salient effects of temporal sampling, the study also revisits the historical problem of correlating point gauge estimates with area estimates, though, the problem is ameliorated by averaging and gridding the data to the area resolution scale of the satellite.

Fisher, Brad↗

NASA NTRS1977

The correlation function for density perturbations in an expanding universe. I - Linear theory

The evolution of the two-point correlation function for adiabatic density perturbations in the early universe is studied. Analytical solutions are obtained for the evolution of linearized spherically symmetric adiabatic density perturbations and the two-point correlation function for these perturbations in the radiation-dominated portion of the early universe. The results are then extended to the regime after decoupling. It is found that: (1) adiabatic spherically symmetric perturbations comparable in scale with the maximum Jeans length would survive the radiation-dominated regime; (2) irregular fluctuations are smoothed out up to the scale of the maximum Jeans length in the radiation era, but regular fluctuations might survive on smaller scales; (3) in general, the only surviving structures for irregularly shaped adiabatic density perturbations of arbitrary but finite scale in the radiation regime are the size of or larger than the maximum Jeans length in that regime; (4) infinite plane waves with a wavelength smaller than the maximum Jeans length but larger than the critical dissipative damping scale could survive the radiation regime; and (5) black holes would also survive the radiation regime and might accrete sufficient mass after decoupling to nucleate the formation of galaxies.

Mcclelland, J.↗

DOE OSTI2021

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗

DOE OSTI2023

Aerosol Particle Deposition on a Spent Nuclear Fuel Assembly Spacer Grid

The flow and particle deposition patterns on surfaces in an idealized spacer grid for a 17x17 pressurized water reactor (PWR) assembly in a spent fuel canister are modeled using computational fluid dynamics (CFD) with laminar flow. The effects of gravitational settling, non-Stokesian flow, and particle slip are first rigorously analyzed. From the analysis, non-Stokesian effects and slip may be neglected for the particle sizes and conditions expected in a canister. For particles that do not settle out, a swirling flow pattern at the corners of a spacer grid channel directs particles to the leeward side of the flow vanes where much of the deposition occurs. Particle deposition increases with increasing particle diameter. Deposition also increases with decreasing flow velocity as this provides more time for particles to settle and deposit on the leeward side of the flow vanes. The fraction of particles that are transmitted through a spacer grid is determined as a function of inlet gas velocity and particle diameter by running the CFD calculation for each set of conditions and for each particle diameter. Curve fits of the transmission curve as a function of particle diameter for a specified spacer grid and flow velocity are applied to a lognormal particle mass density function for the inlet particles. The resulting mass density function and aerosol mass fraction that passes through the spacer grid can be determined analytically without resorting to numerical iteration. A sample calculation of the analytical solution is demonstrated for a lognormal particle mass density function.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure—An experimental and theoretical study

The 57Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α→ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth’s inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe2O3. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1983

A non-Gaussian statistical model for surface elevation of nonlinear random wave fields

Probability density function of the surface elevation of a nonlinear random wave field is obtained. The wave model is based on the Stokes expansion carried to the third order for both deep water waves and waves in finite depth. The amplitude and phase of the first-order component of the Stokes wave are assumed to be Rayleigh and uniformly distributed and slowly varying, respectively. The probability density function for the deep water case was found to depend on two parameters: the root-mean-square surface elevation and the significant slope. For water of finite depth, an additional parameter, the nondimensional depth, is also required. An important difference between the present result and the Gram-Charlier representation is that the present probability density functions are always nonnegative. It is also found that the 'constant' term in the Stokes expansion, usually neglected in deterministic studies, plays an important role in determining the details of the density function. The results compare well with laboratory and field experiment data.

Huang, N. E.↗

DOE OSTI2020

Incorporating Electronic Information into Machine Learning Potential Energy Surfaces via Approaching the Ground-State Electronic Energy as a Function of Atom-Based Electronic Populations

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). Additionally, the new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on LinHn clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Global functional connectivity density alterations in patients with bipolar disorder with auditory verbal hallucinations and modest short‐term effects of transcranial direct current stimulation augmentation treatment—Baseline and follow‐up study

Abstract Objectives To investigate the neuroimaging characteristics of auditory verbal hallucinations (AVHs) in patients with bipolar disorder (BP) experiencing depressive episodes with and without AVHs, and alterations in those characteristics after transcranial direct current stimulation (tDCS). Methods For a baseline pilot study, we recruited 80 patients with BP and depressive status (40 with and 40 without AVHs), and 40 healthy controls (HCs). Their global functional connectivity density (gFCD) was screened by functional magnetic resonance imaging (fMRI). Voxel‐wise one‐way analysis of covariance (ANCOVA) was conducted to detect intergroup differences in gFCD. In a follow‐up study, the effects of 5 weeks of tDCS augmentation treatment on clinical symptoms and gFCD were assessed in the 40 BP patients with AVHs. Results Compared to HCs, BP patients with and without AVHs exhibited increased gFCD in the central parietal lobe, insular lobe, and middle cingulate cortex, with decreased gFCD in the posterior parietal cortex, lateral prefrontal cortex, and occipital lobe (all bilateral). Only patients with AVHs showed increased gFCD in the Broca and Wernicke regions, and decreased gFCD in the hippocampus (all bilateral). After 5 weeks of tDCS, AVHs were slightly alleviated and gFCD abnormalities in the hippocampus were mildly attenuated. Conclusions Patients with BP and AVHs showed disturbances in the brain's communication capacity mainly in the left frontoparietal network, control network, and memory circuitry. Five weeks of tDCS alleviated AVHs slightly, without improving depressive symptoms, and attenuated hippocampal gFCD alterations in these patients.

Zhuo, Chuanjun↗

NASA NTRS1972

A comparative study of stability methods for flexible satellites.

This paper compares three approaches to the stability of hybrid dynamical systems, all three methods being based on the Liapunov direct method. The first method uses testing density functions, whereas the second involves defining certain integral coordinates. Both the method using testing density functions and the method of integral coordinates lead to closed-form stability criteria in terms of the system parameters. Criteria obtained using the method of integral coordinates are in general less restrictive than those derived by the method using testing density functions. On the other hand, the latter method is easier to apply and requires less work than the former. The third method is the standard modal analysis. The modal analysis generally yields more involved criteria, depending on the number of modes used to represent the elastic displacements. As an application, the attitude stability of an earth-pointing satellite with multi-elastic domains is investigated.

Meirovitch, L.↗

DOE OSTI2022

Understanding the fundamental driver of semiconductor radiation tolerance with experiment and theory

Space is the operating environment of a multitude of systems that our society is heavily reliant on today, however, maintaining operability there necessitates special consideration of the electronic systems' tolerance of space radiation. Electronic systems are critically dependent on the electronic properties of their semiconductor components, which are modified by space radiation with an adverse impact on the space system performance. What innate property allows some semiconductors to sustain little damage while others accumulate defects rapidly with dose is poorly understood, which limits the extent to which radiation tolerance can be implemented as a design criterion. Here, to gain insight into what properties are drivers of semiconductor radiation tolerance, the first step is to generate a dataset of the relative radiation tolerance of a broad sampling of semiconductors. To accomplish this, Rutherford backscatter channeling experiments are used to compare the displaced lattice atom buildup in InAs, InP, GaP, GaN, ZnO, MgO, and Si as a function of stepwise alpha particle dose. With this experimental information on radiation-induced incorporation of interstitial defects in hand, hybrid density functional theory electron densities (and their derived quantities) are calculated and their gradient and Laplacian are evaluated to obtain key fundamental information about the interactions in each material. It is shown that simple, undifferentiated values (which are typically used to describe bond strength) are insufficient to predict radiation tolerance. Instead, the curvature of the electron density at bond critical points provides a measure of radiation tolerance consistent with the experimental results obtained. This curvature and associated forces surrounding bond critical points have the potential to disfavor the localization of displaced lattice atoms at these points, favoring their diffusion toward perfect lattice positions. With this criterion to predict radiation tolerance, simple density functional theory simulations can be conducted on potential new materials to gain insight into how they may operate in demanding high radiation environments.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Orbital contributions in the element-resolved valence electronic structure of Bi 2 Se 3

In this work, we studied the bulk band structure of a topological insulator (TI) Bi 2 Se 3 and determined the contributions of the Bi and Se orbital states to the valence bands using standing-wave excited hard x-ray photoemission spectroscopy (SW-HAXPES). This SW technique can provide the element-resolved information and extract individual Bi and Se contributions to the Bi2Se3 valence band. Comparisons with density-functional theory calculations (local density approximation and GW) reveal that the Bi 6s, Bi 6p, and Se 4p states are dominant in the Bi 2 Se 3 HAXPES valence band. Furthermore, these findings pave a way for studying the element-resolved band structure and orbital contributions of this class of TIs.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Atomic-Site-Specific Surface Valence-Band Structure from X-Ray Standing-Wave Excited Photoemission

X-ray standing-wave (XSW) excited photoelectron emission was used to measure the site-specific valence band (VB) for ½ monolayer (ML) Pt grown on a SrTiO 3 (001) surface. The XSW induced modulations in the core level (CL), and VB photoemission from the surface and substrate atoms were monitored for three ℎ⁡kl substrate Bragg reflections. The XSW CL analysis shows the Pt to have a face-centered-cubic-like cube-on-cube epitaxy with the substrate. The XSW VB information compares well to a density functional theory calculated projected density of states from the surface and substrate atoms. Altogether, this Letter represents a novel method for determining the contribution to the density of states by valence electrons from specific atomic surface sites.

74 ATOMIC AND MOLECULAR PHYSICS↗