Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculations”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗
DOE OSTI2024

BEAST DB: Grand-Canonical Database of Electrocatalyst Properties

We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze the electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain a granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage.

database↗
NASA NTRS1975

Short-term time variability of Cygnus X-1

Autocorrelation functions, power-density spectra, and cross-correlation functions were calculated for 71 observations of Cyg X-1 by the Uhuru satellite. The average autocorrelation function and power-density spectrum were found to be independent of energy and consistent with a shot-noise source with a characteristic time constant of 0.45 s. The data in different energy bands were found to be 100 percent correlated.

Weisskopf, M. C.↗
DOE OSTI2025

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗
DOE OSTI2024

Origin of light-induced metastability in ZrTe 5

Here we study the nonequilibrium electronic structure of a model Dirac semimetal ZrTe 5 by using time-and-angle resolved photoemission spectroscopy and density functional theory–based electron and phonon calculations. By measuring the electronic dispersion near the Γ point at time delays up to 10 picoseconds, we discovered that the band spectral weight does not recover during the measured temporal window, revealing the existence of a light-induced metastable state in the electronic structure of this material. Our calculations find that the photoexcited A 1⁢g phonon mode leads to a band renormalization that both supports our experimental observations at the zone center and predicts changes to the band structure outside of our experimental window, ultimately showing the evolution from a direct to an indirect gap semimetal; such band renormalization dramatically reduces the electron-hole recombination rate giving rise to the metastability in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Machine-Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

To cope with legacy greenhouse gas emissions and to achieve net-zero emissions by 2050, the U.S. Department of Energy (DOE) is funding efforts to develop direct air capture (DAC), a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been well studied as DAC sorbent materials due to their tunable structural and compositional properties. Thermodynamic simulations using force fields are often used to provide predictions of a material’s performance in many separations. However, these force fields often make assumptions about bonds and the physics of the adsorption process. A new class of force fields called machine-learned force fields (MLFFs) use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). In this work, models were developed to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using MLFFs. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
NASA NTRS1986

The statistics of peaks of Gaussian random fields

A set of new mathematical results on the theory of Gaussian random fields is presented, and the application of such calculations in cosmology to treat questions of structure formation from small-amplitude initial density fluctuations is addressed. The point process equation is discussed, giving the general formula for the average number density of peaks. The problem of the proper conditional probability constraints appropriate to maxima are examined using a one-dimensional illustration. The average density of maxima of a general three-dimensional Gaussian field is calculated as a function of heights of the maxima, and the average density of 'upcrossing' points on density contour surfaces is computed. The number density of peaks subject to the constraint that the large-scale density field be fixed is determined and used to discuss the segregation of high peaks from the underlying mass distribution. The machinery to calculate n-point peak-peak correlation functions is determined, as are the shapes of the profiles about maxima.

Bardeen, J. M.↗
NASA NTRS2006

Intergalactic Photon Spectra from the Far-IR to the UV Lyman Limit for 0 < z < 6 and the Optical Depth of the Universe to High-Energy Gamma Rays

We calculate the intergalactic photon density as a function of both energy and redshift for 0<z<6 for photon energies from.003 eV to the Lyman limit cutoff at 13.6 eV in a (Omega)CDM universe with (Omega)(Lambda)=0.7 and (Omega)m=0.3. The basic features of our backward-evolution model for galaxies were developed in earlier papers by Malkan & Stecker. With a few improvements, we find that this evolutionary model gives predictions of new deep number counts from Spitzer, as well as a calculation of the spectral energy distribution of the diffuse infrared background, which are in good agreement with the data. We then use our calculated intergalactic photon densities to extend previous work on the absorption of high-energy Gamma-rays in intergalactic space owing to interactions with low-energy photons and the 2.7 K cosmic microwave background radiation. We calculate the optical depth of the universe, Tau , for Gamma-rays having energies from 4 GeV to 100 TeV emitted by sources at redshifts from 0 to 5. We also give an analytic fit with numerical coefficients for approximating (E(Gamma), z). As an example of the application of our results, we calculate the absorbed spectrum of the blazar PKS 2155-304 at z=0.117 and compare it with the spectrum observed by the HESS air Cerenkov Gamma-ray telescope array.

Stecker, F. W.↗
DOE OSTI2026

Enhanced stability and ordering of 2D MnBi2Se2Te2

This dataset includes density functional theory and quantum Monte Carlo calculations on MnBi(Se(1-x)Te(x))4 system for its phase stability and defect formation energies. The data supports publication "Increased Defect Resistance and Ordering in MnBi2(Se(1–x)Te(x))4 via Accurate Diffusion Monte Carlo" (https://doi.org/10.1021/acs.jpcc.5c04299).

36 MATERIALS SCIENCE↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Further steps toward the next generation of covariant energy density functionals

The present study aims at further development of covariant energy density functionals (CEDFs) towards more accurate description of binding energies across the nuclear chart. Infinite basis corrections to binding energies in the fermionic and bosonic sectors of the covariant density functional theory are taken into account in the fitting protocol within the covariant density functional theory. In addition, total electron binding energies are used in the conversion of atomic binding energies into nuclear ones. Their dependence on neutron excess is investigated across the nuclear chart within the atomic approach. Furthermore, these factors were disregarded in the previous generation of covariant energy density functionals, but their omission leads to substantial global calculation errors for physical quantities of interest. For example, these errors for binding energies are of the order of 0.8 MeV or higher for the three major classes of covariant energy density functionals.

Binding energy & masses↗
DOE OSTI2024

First-Principles Study on the Role of Cu and Cl-Based Dopants in NiO

Utilization of wide band gap oxide-based materials in thin-film solar energy technologies has increased in recent years. Among the numerous candidate oxide materials, NiO has shown many desirable optoelectronic properties that are applicable to thin-film PV technologies such as cadmium telluride PV. However, one critical factor requiring further investigation is the p-type doping behavior of NiO, specifically when the cadmium telluride solar cell undergoes conventional processes such as copper doping and chlorine-based activation treatment. The previous literature has shown a large degree of variability in hole concentrations in NiO when copper is used as the primary dopant. This study uses first-principles computational modeling based on density functional theory coupled with defect equilibria calculations to quantitatively explore the role of copper and chlorine-based extrinsic dopants in the p-type doping activity of NiO. The study reveals the importance of extrinsic dopants and their binding interactions with nickel vacancies to effectively p-dope NiO. It is suggested that both the formation of V Ni + Cu Ni and V Ni + Cl O defect pairs under a supersaturated state of NiO are potential mechanisms for increasing hole densities. On the other hand, the production of 2Cu Ni severely limits the effectiveness of p-doping in NiO, even in the presence of the aforementioned defect pairs. Furthermore, the study provides a guideline for experimentalists interested in using copper or chlorine species to understand how to controllably p-dope NiO during thin-film synthesis.

chlorine↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
NASA NTRS1998

Operational Data Reduction Procedure for Determining Density and Vertical Structure of the Martian Upper Atmosphere from Mars Global Surveyor Accelerometer Measurements

The success of aerobraking by the Mars Global Surveyor (MGS) spacecraft was partly due to the analysis of MGS accelerometer data. Accelerometer data was used to determine the effect of the atmosphere on each orbit, to characterize the nature of the atmosphere, and to predict the atmosphere for future orbits. To interpret the accelerometer data, a data reduction procedure was developed to produce density estimations utilizing inputs from the spacecraft, the Navigation Team, and pre-mission aerothermodynamic studies. This data reduction procedure was based on the calculation of aerodynamic forces from the accelerometer data by considering acceleration due to gravity gradient, solar pressure, angular motion of the MGS, instrument bias, thruster activity, and a vibration component due to the motion of the damaged solar array. Methods were developed to calculate all of the acceleration components including a 4 degree of freedom dynamics model used to gain a greater understanding of the damaged solar array. The total error inherent to the data reduction procedure was calculated as a function of altitude and density considering contributions from ephemeris errors, errors in force coefficient, and instrument errors due to bias and digitization. Comparing the results from this procedure to the data of other MGS Teams has demonstrated that this procedure can quickly and accurately describe the density and vertical structure of the Martian upper atmosphere.

Cancro, George J.↗
DOE OSTI2024

Ab-initio vibrational properties of WTe2 and Li0.5WTe2

We present ab-initio calculations of the Td and 1T’ phases of WTe2 and the Td’ phase of Li 0.5 WTe 2 using density functional theory (DFT). These include lattice relaxations, total energy calculations, phonon calculations, and calculations of elastic constants. We also provide a detailed view of the transformation from Td to Td’ using a non-primitive representation of the Td lattice that more closely resembles the Td’ structure. The details presented here are intended as a supplement to a forthcoming experimental-computational paper on lithiated WTe 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

The role of defect charge, crystal chemistry, and crystal structure on positron lifetimes of vacancies in oxides

Density functional theory based positron lifetime (PL) calculations for cation and oxygen monovacancies in a range of oxides—hematite, magnetite, hercynite, and alumina—have been conducted to compare the impact of defect chemistry and crystal structure on the predicted lifetimes. The role of defect charge state has also been examined. A comparison across the same type of crystalline structure but different composition shows that oxygen vacancies only induce a slight increase in the positron-electron overlap and thus barely modify the PL as compared to the bulk. A much more substantial increase of PL is observed for cation monovacancies, regardless of crystal structure or the elemental nature of the vacancy, which we ascribe to an enhanced localization of charge density around the vacant site. The structural and compositional richness of the oxide leads to longer defect PLs, with defected hercynite exhibiting the longest PLs. The charge state of cation monovacancies modifies only by a small percentage the positron localization, relegating to secondary importance the metal defect's oxidation state in modifying the lifetime of positrons within vacancy traps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

Theoretical temperature-dependent branching ratios and laser thresholds of the 3F4 to 3H6 levels of Tm(3+) in ten garnets

The calculated energy levels, the branching ratios, and the estimated thresholds for thulium operating on the 3F4 to 3H6 transitions are reported. Garnet materials with the general formula A3B2C3O12 are evaluated. Calculations are performed for the A side under the assumption of D2 symmetry. X-ray data available in the literature are used to evaluate the crystal-field components, A sub nm. Even-n components are employed to calculate the crystal-field splittings within the manifold. Thermal occupation factors are determined in a straightforward manner using a Boltzmann distribution for the respective manifolds. Odd-n components are applied to calculate the transition probabilities for electric field transitions. It is determined that the magnetic dipole contributions to the transition probability are comparable to the electric dipole contributions in some cases. Thresholds as a function of the density of thulium atoms are calculated.

Filer, Elizabeth D.↗