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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Crystal Structure Prediction of Binary Alloys via Deep Potential

Predicting crystal structure has been a challenging problem in physics and materials science for a long time. A reliable energy calculation engine combined with an efficient global search algorithm, such as particle swarm optimization algorithm or genetic algorithm, is needed to conduct crystal structure prediction. In recent years, machine learning-based interatomic potential energy surface models have been proposed, potentially allowing us to perform crystal structure prediction for systems with the accuracy of density functional theory (DFT) and the speed of empirical force fields. In this paper, we employ a previously developed Deep Potential model to predict the intermetallic compound of the aluminum–magnesium system, and find six meta-stable phases with negative or nearly zero formation energy. In particular, Mg 12 Al 8 shows excellent ductility and Mg 5 Al 27 has a high Young's modulus. Based on our benchmark results, we propose a relatively robust structure screening criterion that selects potentially stable structures from the Deep Potential-based convex hull and performs DFT refinement. By using this criterion, the computational cost needed to construct the convex hull with ab initio accuracy can be dramatically reduced.

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

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗
DOE OSTI2025

Examining Ni Coarsening in Solid Oxide Electrolysis Cells by Characterizing NiH on Ni (111) Using a Combined Theoretical Approach

Ni coarsening in the fuel electrode of solid oxide cells (SOCs) is an important degradation mechanism. In this talk, density-functional theory and kinetic Monte Carlo methods are used to explore the hypothesis that the surface diffusion of NiH on Ni may promote Ni coarsening in the SOC operated in electrolysis cell mode. Using both methods and defining the diffusivity as the product of the surface coverage and single-molecule diffusivity, the diffusivity of NiH on Ni (111) is found to be sufficiently large under a significant overpotential to support the above hypothesis. Also, the time between the formation and dissociation of NiH on Ni (111) is predicted to be short at low coverages of H on Ni (111). Thus, significant progress is made toward developing a model of Ni coarsening considering both molecular and dissociated forms of NiH on Ni (111).

density functional theory (DFT)↗
DOE OSTI2025

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anomalous stability of non-van der Waals bonded B4C nanosheets through surface reconstruction

Boron carbide (B4C) has been well studied both theoretically and experimentally in its bulk form due to its exceptional hardness and use as a high-temperature thermoelectric. However, the properties of its two-dimensional nanosheets are not well established. In this paper, using van der Waals-corrected density-functional theory simulations, we show that bulk B4C can be cleaved along different directions to form B4C nanosheets with low formation energies. We find that there is minimal dependence of formation energies on cleavage planes and surface terminations, even though the bulk is not van der Waals layered. This anomalous stability of B4C nanosheets is found to be a result of surface reconstructions that are unique to B-rich systems. While the density of states of the bulk B4C indicate that it is a semiconductor, the B4C nanosheets are found to be predominantly metallic. We attribute this metallic behavior to a redistribution of charges on the surface bonds of the films. The Seebeck coefficients of the B4C films remain comparable to those of the bulk and are nearly constant as a function of temperature. Our results provide guidance for experimental synthesis efforts and future application of B4C nanosheets in nanoelectronic and thermoelectric applications.

Physics↗
DOE OSTI2025

Static and dynamic properties of atomic nuclei with high-resolution potentials

Here, we compute ground-state and dynamical properties of 4 He and 16 O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schrodinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the 4 He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of 16 O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Pathways for Electron Transfer at MgO–Water Interfaces from Ab Initio Molecular Dynamics

The nature of electron transfer across metal oxide–water interfaces depends significantly on the band gap of the oxide and its band edge energies relative to the potentials of relevant aqueous redox couples. Here we focus on the water interface with MgO, a prototypical wide band gap oxide whose conduction band edge is close in energy to that of water. We investigate the behavior of an excess electron at and out of equilibrium near the interface using ab initio molecular dynamics based on hybrid density functional theory. Our simulations show that under equilibrium conditions the excess electron (donated by an Al impurity in MgO) localizes to a midgap defect state comparable in energy and shape to a hydrated electron in bulk water. To characterize the electron transfer from the conduction band of MgO to interfacial product states, we dope near-equilibrium configurations of the pristine MgO–water system with Al and run short trajectories of these instantaneously out-of-equilibrium systems. Here we observe two distinct products associated with the excess electron: a surface-localized electron (e surf –) and an aqueous hydrogen radical (H • ). The H • pathway exhibits a much higher activation barrier despite being more exoergic, making e surf – the kinetic product. Our characterization of the pathways on the basis of Marcus theory is consistent with the poor observed utility of MgO for water radiolysis. Moreover, we anticipate that the computational framework employed here will be broadly applicable to assessing electron transfer mechanisms at aqueous, photocatalytic interfaces.

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

Binary and Ternary Colloidal Cu–Sn–Te Nanocrystals for Thermoelectric Thin Films

Recent advances in copper chalcogenide-based nanocrystals (NCs), copper sulfide, and copper selenide derived nanostructures, have drawn considerable attention. However, reports of crystal phase and shape engineering of binary or ternary copper telluride NCs remain rare. Here, a colloidal hot-injection approach for producing binary copper/tin telluride, and ternary copper tin telluride NCs with controllable compositions, crystal structures, and morphologies is reported. Here, the crystal phase and growth behavior of these tellurides are systematically studied from both experimental and theoretical perspectives.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deciphering the multiple deformation mechanisms responsible for sustained work hardening in a FeCrCoNi medium entropy alloy

Two important and desirable properties of materials for most structural applications are high tensile strength and ductility, which typically require high work hardening to delay necking. Here, in this work, we designed and tensile tested a face-centered cubic (fcc) Fe-Cr-Co-Ni medium-entropy alloy in which multiple deformation mechanisms are triggered during tensile loading at different temperatures to induce sustained work hardening. Our strategy involved control of the relative stabilities of the fcc, hcp (hexagonal close-packed), and bcc (body-centered cubic) phases in this quaternary system via high-throughput thermodynamic calculations. This alloy not only exhibits extensive deformation-induced nanotwinning at room temperature, but also displays a two-step sequential phase transformation [γ (fcc) → ε (hcp) martensite → α’ (bcc) martensite] at 77 K, which contrasts with the single-step phase transformation [γ → ε martensite] observed in many other fcc high/medium entropy alloys with a low stacking fault energy. The sequence of phase transformation at 77K was supported by first-principles density functional theory calculations. This work provides new templates for the design of alloys capable of multiple deformation mechanisms for sustained work hardening.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Excitonic Effects in X-ray Absorption Spectra of Fluoride Salts and Their Surfaces

Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode–electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. Here we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.

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

First-Principles Examination of Multiple Criteria of Organic Solvent Oxidative Stability in Batteries

Oxidative instability of the liquid electrolyte at or near battery cathode oxide surfaces has significant detrimental effects on batteries. Organic solvent molecules are often the fuel and precursors of such degradation processes, releasing electrons and protons that react with cathode oxides and electrolyte anions. These reactions contribute to cathode–electrolyte interphase (CEI) film formation, transition-metal ion dissolution, and phase transformation of the surface regions of the cathode. Here we apply density functional theory calculations to examine four criteria of oxidative stability (oxidation potential, hydrogen removal energies, and initial reactivity on two types of oxide facets) using four different solvent/additive molecules (ethylene carbonate, fluoroethylene carbonate, 1,3-dioxolane, and dimethyl ether). The ranking of molecular stability differs with each criterion. Surprisingly, the all-oxygen-terminated basal planes of layered oxides exhibit lower reaction barriers than spinel surface facets with exposed transition-metal cations, especially for ether solvents; the calculations also suggest basal planes contribute to the dissolution of transition-metal ions. Further, the structure–degradation relation complexity underscores the challenge of understanding the function of the CEI but also offers a guide to future degradation-mitigation strategies including facet engineering. Our predictions and models help establish a framework for future studies relevant to high-voltage conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Flash Communication: Boron K-edge XAS and TDDFT Studies of Covalent Metal–Ligand Bonding in Ni(C 2 B 9 H 11 ) 2

Ligand K-edge X-ray absorption spectroscopy (XAS), a technique that can measure variations in covalent metal–ligand bonding, has rarely been used to assess covalency in complexes containing metal–boron bonds. Here we describe ligand K-edge XAS and TDDFT studies of the Ni dicarbollide complex Ni(C 2 B 9 H 11 ) 2 (1) and the Ni-free salt (HNMe 3 )(C 2 B 9 H 12 ) (L1). The XAS spectrum for 1 reveals a pre-edge feature indicative of covalent Ni–B bonding, which is corroborated by time-dependent density functional theory (TDDFT) calculations and comparative analysis to L1 and inner-shell electron energy loss spectroscopy (ISEELS) collected on the same Ni complex.

Boron↗
DOE OSTI2023

Synergistic effects of mixing and strain in high entropy spinel oxides for oxygen evolution reaction

Developing stable and efficient electrocatalysts is vital for boosting oxygen evolution reaction (OER) rates in sustainable hydrogen production. High-entropy oxides (HEOs) consist of five or more metal cations, providing opportunities to tune their catalytic properties toward high OER efficiency. This work combines theoretical and experimental studies to scrutinize the OER activity and stability for spinel-type HEOs. Density functional theory confirms that randomly mixed metal sites show thermodynamic stability, with intermediate adsorption energies displaying wider distributions due to mixing-induced equatorial strain in active metal-oxygen bonds. The rapid sol-flame method is employed to synthesize HEO, comprising five 3d-transition metal cations, which exhibits superior OER activity and durability under alkaline conditions, outperforming lower-entropy oxides, even with partial surface oxidations. The study highlights that the enhanced activity of HEO is primarily attributed to the mixing of multiple elements, leading to strain effects near the active site, as well as surface composition and coverage.

08 HYDROGEN↗