Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Field-induced magnetization of δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of magnetization, induced by an external magnetic field, for δ-phase plutonium. The fully relativistic electronic structure accounts for Zeeman splitting effects through a Hamiltonian that couples the magnetic field to both spin and orbital magnetic moments. The electronic-structure model is further improved by an extension to DFT in terms of an orbital-orbital coupling via the conventional orbital-polarization method, as has routinely been done for plutonium. The response to the applied magnetic field is shown to be weak in the DFT model, with induced magnetic moments of the order of in magnetic fields up to 30 T. These results are in accord with recent assessments from x-ray magnetic circular dichroism measurements in magnetic field on δ-plutonium. Somewhat surprisingly, a model assuming δ-plutonium to be absent magnetic moments shows stronger response to an external static magnetic field than models allowing for antiferromagnetism or magnetic disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Investigating 3,4-bis(3-nitrofurazan-4-yl)furoxan detonation with a rapidly tuned density functional tight binding model

In this work, we describe a machine learning approach to rapidly tune density functional tight binding models for the description of detonation chemistry in organic molecular materials. Resulting models enable simulations on the several 10s of ps scales characteristic to these processes, with “quantum-accuracy.” We use this approach to investigate early shock chemistry in 3,4-bis(3-nitrofurazan-4-yl)furoxan, a hydrogen-free energetic material known to form onion-like nanocarbon particulates following detonation. We find that the ensuing chemistry is significantly characterized by the formation of large C x N y O z species, which are likely precursors to the experimentally observed carbon condensates. Beyond utility as a means of investigating detonation chemistry, the present approach can be used to generate quantum-based reference data for the development of full machine-learned interatomic potentials capable of simulation on even greater time and length scales, i.e., for applications where characteristic time scales exceed the reach of methods including Kohn–Sham density functional theory, which are commonly used for reference data generation.

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

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. In this paper, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

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

Modal density function and number of propagating modes in ducts

Often raised questions in duct sound propagation studies involve the total number of propagating modes, the number of propagating radial modes for a particular spinning lobe number, and the number of modes possible between two given values of cutoff ratio or eigenvalue. These questions can be answered approximately by using the modal distribution function which is the integral of the modal density function for ducts in a manner similar to that previously published for architectural acoustics. The modal density functions are derived for rectangular and circular ducts with a uniform steady flow. Results from this continuous theory are compared to the actual (discrete) modal distributions.

Rice, E. J.↗
DOE OSTI2021

Global study of separable pairing interaction in covariant density functional theory

A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on the separable pairing interaction. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. For a given part of the nuclear chart the scaling factors for spherical nuclei are smaller than those for deformed ones; this feature exists also in nonrelativistic density functional theories. Furthermore, its origin is traced back to particle-vibration coupling in odd-A nuclei which is missing in all existing global studies of pairing. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1984

The non-Gaussian joint probability density function of slope and elevation for a nonlinear gravity wave field

On the basis of the mapping method developed by Huang et al. (1983), an analytic expression for the non-Gaussian joint probability density function of slope and elevation for nonlinear gravity waves is derived. Various conditional and marginal density functions are also obtained through the joint density function. The analytic results are compared with a series of carefully controlled laboratory observations, and good agreement is noted. Furthermore, the laboratory wind wave field observations indicate that the capillary or capillary-gravity waves may not be the dominant components in determining the total roughness of the wave field. Thus, the analytic results, though derived specifically for the gravity waves, may have more general applications.

Huang, N. E.↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
DOE OSTI2021

Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch–Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Finally, conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

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

Magnetization dynamics with time-dependent spin-density functional theory: Significance of exchange-correlation torques

In spin-density functional theory (SDFT) for noncollinear magnetic materials, the Kohn-Sham system features exchange-correlation (xc) scalar potentials and magnetic fields. The significance of the xc magnetic fields is not very well explored; in particular, they can give rise to local torques on the magnetization, which are absent in standard local and semilocal approximations. Here, exact benchmark solutions for a five-site extended Hubbard lattice at half filling and in the presence of spin-orbit coupling are compared with SDFT results obtained using orbital-dependent exchange-only approximations. The magnetization dynamics following short-pulse excitations is found to be reasonably well described in the exchange-only approximation for weak to moderate interactions. For stronger interactions and near transitions between magnetically ordered and frustrated phases, exchange and correlation torques tend to compensate each other and must both be accounted for.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Thermal hybrid exchange-correlation density functional for improving the description of warm dense matter

Finite-temperature density functional theory (FT-DFT) has become a standard tool for first principles calculations of the properties of warm dense matter (WDM) relevant to high-energy density physics (HEDP) applications. In this work we introduce theoretical grounds of thermal hybrid exchange-correlation (XC) functionals within the generalized Mermin-Kohn-Sham scheme for an improved description of WDM. Building on the previously developed KSDT (Karasiev-Sjostrom-Dufty-Trickey) [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)] local density approximation (LDA) and the KDT16 (Karasiev-Dufty-Trickey 2016) [Karasiev et al., Phys. Rev. Lett. 120, 076402 (2018)] generalized-gradient approximation (GGA) XC free-energy density functionals, we construct a new thermal hybrid XC functional, referred to here as KDT0. The KDT0 model at low temperature reduces to the popular ground-state PBE0 hybrid due to properties of the used KDT16 density functional approximation. Application to static calculations of electronic band gap and band structure at a wide range of temperatures for various systems of interest to HEDP show that KDT0 provides a significant improvement to the lower LDA and GGA rung XC functionals and to the ground-state PBE0 hybrid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations

In this research, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew–Burke–Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd–Scuseria–Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

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

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗