Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy 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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Semiconductor Nanostructures by Scientific Design (Final Technical Report)

We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

NNL.Fe.qSNAP-ZBL.2024.1: A Fe Spectral Neighbor Analysis Potential for Radiation Damage Simulations

The NNL.Fe.qSNAP-ZBL.2024.1 machine-learned potential (MLP) has been generated to support the development of an elemental body-centered cubic (BCC) Fe athermal recombination corrected neutron damage model and simulations of primary recoil atom (PRA) cascades in BCC Fe. This MLP is a quadratic spectral neighbor analysis potential (qSNAP) hybridized with the universal Ziegler-Beirsack-Littmark (ZBL) potential at short-range and is named according to Naval Nuclear Laboratory MLP naming conventions (NNL.material-system.MLP-type.year.version). Training set calculations for Fe are presented along with the subsequent MLP fitting procedure. A key criterion of the fitting procedure is that ZBL describes the short-range interaction with minimal impact on the MLP. The MLP is compared to density functional theory (DFT) predicted properties relevant to radiation damage simulation, including threshold displacement energies, for validation. The NNL.Fe.qSNAP-ZBL.2024.1 potential is considered suitable for molecular dynamics (MD) simulations of radiation defects up to 800 K and PRA cascades in BCC Fe up to around 10 keV. The potential can additionally be used on a limited basis for recoils of 10–20 keV, within which range the emergence of structures outside the training set in cascade simulations may cause system instabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Tunneling magnetoresistance in MgO tunnel junctions with Fe-based leads in empirically corrected density functional theory

The minority-spin Fe/MgO interface states are at the Fermi level in density functional theory (DFT), but experimental evidence and GW calculations place them slightly higher in energy. This small shift can strongly influence tunneling magnetoresistance (TMR) in junctions with a thin MgO barrier and its dependence on the concentration of Co in the electrodes. Here, in this study, an empirical potential correction to DFT is introduced to shift the interface states up to match the tunnel spectroscopy data. With this shift, TMR in Fe/MgO/Fe junctions exceeds 800% and 3000% at 3 and 4 monolayers (ML) of MgO, respectively. We further consider the effect of alloying of the Fe electrodes with up to 30% Co or 10% V, treating them in the coherent potential approximation (CPA). Alloying with Co broadens the interface states and brings a large incoherent minority-spin spectral weight to the Fermi level. Alloying with V brings the minority-spin resonant states close to the Fermi level. However, in both cases the minority-spin spectral weight at the Fermi level resides primarily at the periphery of the Brillouin zone, which is favorable for spin filtering. Using convolutions of k II -resolved barrier densities of states calculated in CPA, it is found that TMR is strongly reduced by alloying with Co or V but still remains above 500% at 4 ML of MgO up to 30% of Co or 5% V. At 5 ML, the TMR increases above 1000% in all systems considered. However, while TMR declines sharply with increasing bias up to 0.2 eV in the tunnel junctions with pure Fe leads, it remains almost constant up to 0.5 eV if leads are alloyed with Co.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Charge state transition levels of Ni in $β$-Ga 2 O 3 crystals from experiment and theory: An attractive candidate for compensation doping

Nickel-doped β-Ga 2 O 3 crystals were investigated by optical absorption and photoconductivity, revealing Ni-related deep levels. Here, the photoconductivity spectra were fitted using the phenomenological Kopylov and Pikhtin model to identify the energy of the zero-phonon transition (thermal ionization), Franck–Condon shift, and effective phonon energy. The resulting values are compared with the predicted ones by first-principle calculations based on the density functional theory (DFT). An acceptor level (0/-) of 1.9 eV and a donor level (+/0) of 1.1 eV above the valence band minimum are consistently determined for Ni Ga , which preferentially incorporates on the octahedrally coordinated Ga site. Temperature-dependent resistivity measurements yield a thermal activation energy of ~2.0 eV that agrees well with the determined Ni acceptor level. Conclusively, Ni is an eminently suitable candidate for compensation doping for producing semi-insulating β-Ga 2 O 3 substrates due to the position of the acceptor level (below and close to the mid-bandgap).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Using DAPPER to extract the photon strength function of 58Fe using the inverse Oslo and shape methods

The photon strength function of 58 Fe has been extracted using both the Oslo and Shape methods from particle-γ coincidence data measured using the Detector Array for Photons, Protons, and Exotic Residues, which probes nuclei utilizing (d,p) reactions in inverse kinematics. Four particle-γ coincidence matrices, each constructed with different treatments of the γ-ray energies, are explored in order to observe the impact on the resulting nuclear level density and photon strength. The final photon strength function reported is found to agree well with previous Oslo measurements of other iron isotopes. Systematic uncertainties are included, using different model parameters and their reported errors to perform the Oslo method normalization. The model-independent Shape method is explored and the functional form of the photon strength function obtained is in agreement with the Oslo method results. A low-energy enhancement is not reported for 58 Fe in this work given possible subtraction issues originating from strongly populated states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Photodriven electron-transfer dynamics in a series of heteroleptic Cu(I)–anthraquinone dyads

Solar fuels catalysis is a promising route to efficiently harvesting, storing, and utilizing abundant solar energy. To achieve this promise, however, molecular systems must be designed with sustainable components that can balance numerous photophysical and chemical processes. To that end, here we report on the structural and photophysical characterization of a series of Cu(I)–anthraquinone-based electron donor–acceptor dyads. The dyads utilized a heteroleptic Cu(I) bis-diimine architecture with a copper(I) bis-phenanthroline chromophore donor and anthraquinone electron acceptor. We characterized the structures of the complexes using x-ray crystallography and density functional theory calculations and the photophysical properties via resonance Raman and optical transient absorption spectroscopy. The calculations and resonance Raman spectroscopy revealed that excitation of the Cu(I) metal-to-ligand charge-transfer (MLCT) transition transfers the electron to a delocalized ligand orbital. The optical transient absorption spectroscopy demonstrated that each dyad formed the oxidized copper–reduced anthraquinone charge-separated state. Unlike most Cu(I) bis-phenanthroline complexes where increasingly bulky substituents on the phenanthroline ligands lead to longer MLCT excited-state lifetimes, here, we observe a decrease in the long-lived charge-separated state lifetime with increasing steric bulk. The charge-separated state lifetimes were best explained in the context of electron-transfer theory rather than with the energy gap law, which is typical for MLCT excited states, despite the complete conjugation between the phenanthroline and anthraquinone moieties.

14 SOLAR ENERGY↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core– electron binding energy (CEBE) analysis. We extend the development of realspace pseudopotential methods to calculating 1s, 2s, and 2p 3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn–Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

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

Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo

Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative—particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods.

Chemistry↗
DOE OSTI2023

Defect graph neural networks for materials discovery in high-temperature clean-energy applications

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. In conclusion, it will help to tackle future materials discovery problems in clean energy and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

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

Photophysical properties of zero-dimensional perovskites studied by PBE0 and GW+BSE methods

Due to the exceptional photophysical properties and great potential applications in solid state lighting, zero-dimensional halide perovskites with broadband emitting have recently attracted much attention. In this work, to reveal the origin of the broadband absorption and emission spectra of 0D (C 4 N 2 H 14 Br) 4 SnBr x I 6-x (x = 3), the ground state and exciton dynamics properties have been studied by hybrid PBE0 functional and the GW+BSE method. Different self-trapped excitons (STEs) are obtained, and their emission energies by PBE0 are consistent well with experimental data. By considering the vibrational broadening of emission, our calculated full width at half maximum of emission accords well with the experiment. The absorption of free exciton and STEs is also calculated by considering the quasiparticle and excitonic effect. Our results show that these light-induced trapped states can be observed by experiments, as confirmed in other similar systems by the femtosecond transient absorption spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Energetics of Base–Acid Pairs for the Design of High-Temperature Fuel Cell Polymer Electrolytes

The interaction energy of base–acid plays a key role in acid retention of phosphoric acid (PA)-doped polymer electrolytes under fuel cell operating conditions. Here, we investigate the energetics of proton-accepting and hydroxide-donating organic bases using density functional theory calculations. Because of their weak basicity, proton-accepting organic bases such as benzimidazole have relatively low interaction energy with the acid in the absence of water (15.3–28.0 kcal mol –1 ). Energetics of the proton-accepting base–PA complex increases by adding water, indicating that the interactions in the base–acid complex strengthen in the presence of water. On the other hand, hydroxide-donating organic bases, such as tetramethylammonium hydroxide, have high interaction energy with PA (~110 kcal mol –1 ), which remains high in the presence of water. The chemical shifts of 31P NMR support the energetics of the base–acid complexes. Overall, this study further discusses the benefit of incorporating hydroxide-donating organic bases into the polymeric structure over proton-accepting bases as a way to increase acid retention.

25 ENERGY STORAGE↗
DOE OSTI2025

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

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

Intrinsic Layer-Dependent Surface Energy and Exfoliation Energy of van der Waals Materials

Stacking and twisting 2D van der Waals (vdW) layers have become versatile platforms to tune the electron correlation. These platforms rely on exfoliating vdW materials down to a single vdW layer and a few vdW layers. We calculate the intrinsic layer-dependent surface and exfoliation energies of typical vdW materials such as graphite, h-BN, black P, MX 2 (M = Mo or W; X = S, Se, or Te), MX (M = Ga or In; X = S, Se, or Te), Bi 2 Te 3 , and MnBi 2 Te 4 using density functional theory. For exchange-correlation functionals with explicit vdW interaction, a single vdW layer always has the smallest surface energy, giving a surface energy reduction when compared to that of thicker vdW layers. Furthermore, the magnitude of this surface energy reduction quickly decreases with an increase in the number of atomic layers inside the single vdW layer for different vdW materials. Such atomic-layer dependence in surface energy reduction helps explain the different effectiveness of exfoliation for different vdW materials down to a single vdW layer.

Anode materials↗
DOE OSTI2025

Ba 1−x Sr x FeO 3−δ as an improved oxygen storage material for chemical looping air separation: a computational and experimental study

Chemical looping air separation (CLAS) is a promising technology to generate oxygen-rich gas streams to enable efficient carbon dioxide capture during fossil fuel combustion or gasification. CLAS relies on the capture and release of oxygen from the atmosphere using the redox properties of an oxygen-selective solid oxide carrier. This study investigates the redox characteristics of Ba 1−x Sr x FeO 3−δ (0.0 ≤ x ≤ 0.417, 0.0 ≤ δ ≤ 0.5) using a combination of density functional theory (DFT) calculations and experimental verification using X-ray diffraction, thermogravimetric analysis, and oxygen-temperature-programmed desorption. The DFT computed energies of the Ba 1−x Sr x FeO 3−δ perovskites reveal a composition-dependent transition from hexagonal to cubic phases as the Sr-concentration or oxygen vacancy concentration increases. Oxygen vacancy formation energies of the cubic perovskites are found to be lower than those of their hexagonal counterparts. A low oxygen diffusion barrier of ∼1 eV combined with the thermodynamic preference of Ba 1−x Sr x FeO 3−δ compositions that form in a cubic phase suggests them as promising candidates for oxygen storage applications. The experimental results corroborate this finding by identifying Ba 0.75 Sr 0.25 FeO 3−δ in the cubic phase as an optimal composition offering low-temperature oxygen storage capacities comparable to that of the state-of-the-art Sr 0.75 Ca 0.25 FeO 3−δ perovskite oxygen storage material at 325 °C and 350 °C.

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

Comparative Critical Mass Calculations for NNL and ENDF/B-VIII.0 Zirconium Hydride Thermal Neutron Scattering Laws

Zirconium hydride (ZrH x ) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high- enriched uranium (HEU) and ZrH x in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of 235 U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH x .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Comparative critical mass calculations for NNL and ENDF/B-VIII.0 - Zirconium hydride thermal neutron scattering laws

Zirconium hydride (ZrH{sub x}) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high-enriched uranium (HEU) and ZrH{sub x} in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of {sup 235}U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH{sub x}. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗