Engineering Papers⌕ Search

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.

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

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Constraints on magnetism and correlations in RuO 2 from lattice dynamics and Mössbauer spectroscopy

We provide experimental evidence for the absence of a magnetic moment in bulk RuO 2 , a candidate altermagnetic material, by using a combination of Mössbauer spectroscopy, nuclear forward scattering, inelastic X-ray and neutron scattering, and density functional theory calculations. Using complementary Mössbauer and nuclear forward scattering, we determine the Ru magnetic hyperfine splitting to be negligible. Inelastic X-ray and neutron scattering-derived lattice dynamics of RuO 2 are compared to density functional theory calculations of varying flavors. Comparisons among theory with experiments indicate that electronic correlations, rather than magnetic order, are key in describing the lattice dynamics.

Mössbauer spectroscopy↗
DOE OSTI2022

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗
DOE OSTI2025

High-Pressure Phase Transition of Metastable Wurtzite-Like CuInSe 2 Nanocrystals

Ternary I–III–VI 2 semiconductors, such as CuInSe 2 , exhibit diverse polymorphs with unique structural characteristics and optoelectronic properties. This study investigates the pressure-induced phase transitions of metastable wurtzite-like CuInSe 2 nanocrystals. Using a combination of synchrotron X-ray diffraction, pair distribution function analysis, and density functional theory calculations, we reveal a transition from cation-ordered wurtzite-like (Pmc2 1 ) to cation-disordered NaCl-like (Fm3̅m) structures at 7.7 GPa. The cation-disordered NaCl-like phase persists upon decompression. Bulk modulus calculations highlight size-dependent deviations from bulk material behavior. These findings deepen our understanding of phase stability in colloidal I–III–VI 2 semiconductor nanocrystals, with implications for tailoring functional materials under extreme conditions.

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

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

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

Electron leakage through heterogeneous LiF on lithium–metal battery anodes

The solid–electrolyte interphase (SEI) that forms on lithium ion battery (LIB) anodes prevents degradation-causing transfer of electrons to the electrolyte. Grain boundaries (GBs) between different SEI components, like LiF, have been suggested to accelerate Li+ transport. However, using the non-equilibrium Green's function technique with density functional theory (NEGF-DFT), we find that GBs enhance electron tunneling in thin LiF films by 1–2 orders of magnitude, depending on the bias. Extrapolating to thicker films using the Wentzel–Kramers–Brillouin (WKB) method emphasizes that safer batteries require passivation of GBs in the SEI.

25 ENERGY STORAGE↗
DOE OSTI2022

Electron passivation in CaF 2 on calcium metal anodes

Current electrolytes in calcium-ion batteries suffer from a lack of stability and degradation caused by reduction from the anode. The solid-electrolyte interphase (SEI) that forms on the anodes during operation stems the flow of electrons from the anode to the electrolyte. CaF 2 is a common inorganic compound found in the SEI, and is derived from electrolyte salts such as Ca(PF 6 ) 2 . CaF 2 can exist in crystalline, polycrystalline, and amorphous phases in the SEI, and as recent work has shown, different phases of the same compound can have vastly different electronic conductivities. Using the non-equilibrium Green's function technique with density functional theory (NEGF-DFT), we find that amorphous phase systems enhance electron tunneling in thin CaF 2 films by 1–2 orders of magnitude when compared to crystalline and polycrystalline CaF 2 systems. Transport through several amorphous structures was considered showing that, despite their random structures, their conductance properties are similar. Finally, through analysis of the decay constant β and the low-bias conductance of each system, we show that crystalline and polycrystalline CaF 2 offer greater protection of the electrolyte than amorphous CaF 2 .

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

ExchCXX v0.1.0

ExchCXX is a modern C++ library for the evaluation of exchange-correlation (XC) functionals required for density functional theory (DFT) calculations. In addition to providing high-level wrappers for the ubiquitous Libxc library for the CPU evaluation of XC functionals, ExchCXX also provides locally developed implementations of a small subset of XC functionals which may be evaluated either on the host (CPU) or device (GPU, FPGA, etc).

Williams-Young, David↗
DOE OSTI2021

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4a Raw Data

Calculated formation energy of (0001 ̅) c-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (0001 ̅) c-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗
DOE OSTI2021

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4b Raw Data

Calculated formation energy of (1010 ̅) m-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (1010 ̅) m-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗
DOE OSTI2020

Mesh-Sequenced Realizations for Evaluation of Subgrid-Scale Models for Turbulent Combustion

This paper develops a new approach for analysis of subgrid closures for turbulent combustion as modeled using direct quadrature (finite-rate chemistry) techniques. The approach, termed multiresolution analysis through mesh-sequenced realizations (MRA-MSR), conducts simultaneous, constrained large-eddy simulations on a set of hierarchically coarsened meshes. Furthermore, the availability of underlying fine-mesh (subgrid) data corresponding to coarse-mesh locations allows a clearer assessment of the effects of unresolved fluctuations on apparent reactivity. A key to MRA-MSR is the correlation of eddy structures at coarser mesh levels, which is facilitated by the transfer of filtered fine-mesh velocity information. A seven-mesh MRA-MSR hierarchy using three resolution levels is applied to one of the Sydney bluff-body stabilized methane–hydrogen flames. Analysis of the simultaneously evolved data at different resolution levels reveals several interesting trends. First, at high Damköhler numbers, there is clear evidence of attenuation of apparent reactivity due to the effects of unresolved fluctuations. Secondly, single-point, single-time filtered density functions of a normalized subgrid Damköhler number show a characteristic beta probability density function (PDF) form and display evidence of scale similarity. Interrogation of the MRA-MSR database also shows that the recently-proposed least-squares minimization (LSM) turbulence–chemistry interaction model can account for the observed diminishment in reactivity at high Damköhler numbers but cannot reduce scatter significantly. A new form of the LSM model, which makes use of the normalized subgrid Damköhler number beta PDF distribution, performs slightly better than the original model, illustrating the potential of MRA-MSR both in assessing existing closure concepts and in developing new ones.

42 ENGINEERING↗
DOE OSTI2022

Local Coordination Environment of 3d and 4d Transition Metal Ions in LiCl-KCl Eutectic Mixture

In this study, the structure and coordination environment of two 3d transition elements (Ni and Cr) is investigated in a molten chloride salt system. Electronic absorption spectroscopy was employed to elucidate their coordination environment in 3LiCl-2KCl eutectic salt, as a function of temperature. Density functional theory (DFT) modeling was used to determine the coordination environment of the transition metal species in the eutectic composition as well as the optical spectra computationally. The Ni 2+ and Cr 3+ exist in a tetrahedral and octahedral coordination environment, respectively, in eutectic salt. The spectra thus obtained were compared with the experimental data; a reasonable qualitative agreement was obtained between experimental and computational Ni 2+ and Cr 3+ spectra, and the coordination of both elements in the eutectic composition were in excellent agreement with the experimentally determined results. Computational results were also obtained for two 4d elements, Mo 3+ and Nb 3+ , with both quantum molecular dynamics (QMD) and hybrid functional optical spectra indicating octahedral coordination.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Revisiting the structures and phase transitions of Ba 2 NaNb 5 O 15

The room-temperature and low-temperature structure(s) of Ba 2 NaNb 5 O 15 (BNN) have been debated since the structure was proposed in the 1960s. This work revisits the structures and phase transitions of BNN, combining high-resolution X-ray and neutron powder diffraction with density functional theory calculations. Temperature-dependent high-resolution X-ray powder diffraction patterns are collected from 4 to 918 K, and sequential batch Rietveld refinement using a symmetry mode approach to describe the structure is used to extract the main structural changes as a function of temperature. The data show that the average structure of BNN is best described by the Ama2 space group, and no other structural phase transitions were observed below the ferroelastic transition. The symmetry mode analysis, combining results from diffraction and density functional theory, shows significant octahedral tilting and corrugations of both the A1 and A2 sites along the c direction. A strong correlation between the spontaneous strain and the octahedral tilting was observed, and a potential connection with emerging microstructure at low temperatures is proposed, all enabled by the symmetry mode approach used in this work.

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

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2021

Information Content of the Parity-Violating Asymmetry in Pb 208

The parity-violating asymmetry A PV in 208 Pb, recently measured by the PREX-2 Collaboration, is studied using modern relativistic (covariant) and nonrelativistic energy density functionals. We first assess the theoretical uncertainty on A PV which is intrinsic to the adopted approach. To this end, we use quantified functionals that are able to accommodate our previous knowledge on nuclear observables such as binding energies, charge radii, and the dipole polarizability α D of 208 Pb. We then add the quantified value of A PV together with α D to our calibration dataset to optimize new functionals. Based on these results, we predict a neutron skin thickness in 208 Pb r skin = 0.19 ± 0.02 fm and the symmetry-energy slope L = 54 ± 8 MeV. These values are consistent with other estimates based on astrophysical data and are significantly lower than those recently reported using a particular set of relativistic energy density functionals. Here, we also make a prediction for the A PV value in 48 Ca that will be soon available from the CREX measurement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Unusual Exchange Couplings and Intermediate Temperature Weyl State in Co 3 Sn 2 S 2

Understanding magnetism and its possible correlations to topological properties has emerged to the forefront as a difficult topic in studying magnetic Weyl semimetals. Co 3 Sn 2 S 2 is a newly discovered magnetic Weyl semimetal with a kagome lattice of cobalt ions and has triggered intense interest for rich fantastic phenomena. In this study, we report the magnetic exchange couplings of Co 3 Sn 2 S 2 using inelastic neutron scattering and two density functional theory (DFT) based methods: constrained magnetism and multiple-scattering Green’s function methods. Co 3 Sn 2 S 2 exhibits highly anisotropic magnon dispersions and linewidths below T C , and paramagnetic excitations above T C . The spin-wave spectra in the ferromagnetic ground state is well described by the dominant third-neighbor “across-hexagon” J d model. Our density functional theory calculations reveal that both the symmetry-allowed 120° antiferromagnetic orders support Weyl points in the intermediate temperature region, with distinct numbers and the locations of Weyl points. Our study highlights the important role Co 3 Sn 2 S 2 can play in advancing our understanding of kagome physics and exploring the interplay between magnetism and band topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Correlating the three-dimensional atomic defects and electronic properties of two-dimensional transition metal dichalcogenides

The electronic, optical and chemical properties of two-dimensional transition metal dichalcogenides strongly depend on their three-dimensional atomic structure and crystal defects. Using Re-doped MoS 2 as a model system, here we present scanning atomic electron tomography as a method to determine three-dimensional atomic positions as well as positions of crystal defects such as dopants, vacancies and ripples with a precision down to 4 pm. We further measure the three-dimensional bond distortion and local strain tensor induced by single dopants. By directly offering these experimental three-dimensional atomic coordinates to density functional theory, we obtain more accurate electronic band structures than derived from conventional density functional theory calculations that relies on relaxed three-dimensional atomic coordinates. We anticipate that scanning atomic electron tomography not only will be generally applicable to determine the three-dimensional atomic coordinates of two-dimensional materials, but also will enable ab initio calculations to better predict the physical, chemical and electronic properties of these materials.

36 MATERIALS SCIENCE↗