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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Parameter-free prediction of phase transition in PbTiO 3 through combination of quantum mechanics and statistical mechanics

Thermodynamics of ferroelectric materials and their ferroelectric to paraelectric (FE-PE) transitions is commonly described by the phenomenological Landau theory and more recently by effective Hamiltonian and various potentials, all with model parameters fitted to experimental or theoretical data. Here, in this work, we show that the zentropy theory, which considers the total entropy of a system as a weighted sum of entropies of configurations that the system may experience and the statistical entropy among the configurations, can predict the FE-PE transition without fitting parameters. For PbTiO 3 , the configurations are identified as the FE configurations with 90- or 180° domain walls in addition to the ground state FE configuration without domain wall. With the domain wall energies predicted from first-principles based on density functional theory in the literature as the only inputs, the FE-PE transition for PbTiO 3 is predicted showing remarkable agreement with experiments, unveiling the microscopic fundamentals of the transition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct Deoxygenation of Phenol over Fe-Based Bimetallic Surfaces Using On-the-Fly Surrogate Models

We present an accelerated nudged elastic band (NEB) study of phenol direct deoxygenation (DDO) on Fe-based bimetallic surfaces using a recently developed Gaussian process regression (GPR) calculator. Our test calculations demonstrate that the GPR calculator achieves up to 3 times speedup compared to conventional density functional theory calculations while maintaining high accuracy, with energy barrier errors below 0.015 eV. Using GPR-NEB, we systematically examine the DDO mechanism on pure Fe(110) and surfaces modified with Co and Ni in both top and subsurface layers. Our results show that subsurface Co and Ni substitutions preserve favorable thermodynamics and kinetics for both C–O bond cleavage and C–H bond formation, comparable to those on the pure Fe(110) surface. In contrast, top-layer substitutions generally increase the C–O bond cleavage barrier, render the step endothermic, and result in significantly higher reverse reaction rates, making DDO unfavorable on these surfaces. This work demonstrates the effectiveness of GRR-accelerated transition state searches for complex surface reactions and provides insights into rational design of bimetallic catalysts for selective deoxygenation.

Aromatic compounds↗
DOE OSTI2023

Binding of Li + to Negatively Charged and Neutral Ligands in Polymer Electrolytes

Conceptually, single-ion polymer electrolytes (SIPE) with the anion bound to the polymer could solve major issues in Li-ion batteries, but their conductivity is too low. Experimentally, weakly interacting anionic groups have the best conductivity. To provide a theoretical basis for this result, density functional theory calculations of the optimized geometries and energies are performed for charged ligands used in SIPE. Comparison is made to neutral ligands found in dual-ion conductors, which demonstrate higher conductivity. Further, the free energy differences between adding and subtracting a ligand are small enough for the neutral ligands to have the conductivity seen experimentally. However, charged ligands have large barriers, implying that lithium transport will coincide with the slow polymer diffusion, as observed in experiments. Overall, SIPE will require additional solvent to achieve a sufficiently high conductivity. Additionally, the binding of mono- and bidentate geometries varies, providing a simple and clear reason that polarizable force fields are required for detailed interactions.

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

High-Temperature Rotating Disk Electrode Study of Platinum Bimetallic Catalysts in Phosphoric Acid

Understanding the H 3 PO 4 effect on the catalyst’s activity under a relevant condition is important for high-temperature polymer electrolyte membrane fuel cell (HT-PEMFC) catalyst research. Here, we report a high-temperature rotating disk electrode (HT-RDE) study of oxygen reduction reaction (ORR) in H 3 PO 4 . With the regular electrochemical protocol, we found that H 3 PO 4 reduction could occur during cyclic voltammetry study and form a reductive species—phosphorus acid (H 3 PO 3 ). Further, to obtain reliable ORR measurement, we optimized the protocol to avoid the H 3 PO 3 generation. The ORR activity of carbon-supported PtM (M = Fe, Co, Ni, Ru, Pd, and Ir) bimetallic alloy catalysts measured with this HT-RDE method showed higher ORR activity than Pt. To understand the alloying effect, we combine experiments in diluted solutions to distinguish the alloying effect on Pt–O binding and Pt–H 3 PO 4 binding. The results indicate that H 3 PO 4 mainly reduces available sites for ORR, with little effect on neighboring site’s Pt–O binding via Pt–H 3 PO 4 interaction, which is also supported by the density functional theory calculation of the Pt–O binding energy with/without H 2 PO 4 . Further study in a phosphoric acid-doped quaternary ammonium-biphosphate ion pair coordinated polyphenylene (PA-QAPOH) membrane electrode assembly (MEA) shows that the active alloy catalyst has better performance in both the HT-RDE and MEA. Also, the MEA gives higher ORR activity than the HT-RDE because of the higher pressure and less phosphoric acid content of the MEA. Yet, the gap between the HT-RDE and MEA is significantly smaller than that between the room temperature (RT)-RDE and MEA, suggesting the importance of temperature and H 3 PO 4 concentration in understanding ORR in HT-PEMFCs.

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

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗
DOE OSTI2025

Spin excitations and flat electronic bands in a Cr-based kagome superconductor

In the quest for topology- and correlation-driven quantum states, kagome lattice materials have garnered significant interest for their band structures, featuring flat bands (FBs) from the quantum destructive interference of the electronic wavefunction. Tuning an FB to the chemical potential could induce electronic instabilities and emergent orders. Despite extensive studies, direct evidence of FBs tuned to the chemical potential and their role in emergent orders in bulk materials remains lacking. Using angle-resolved photoemission spectroscopy, resonant inelastic X-ray scattering, and density functional theory, we show that the low-energy structure of the Cr-based kagome metal superconductor CsCr 3 Sb 5 is dominated by FBs at the Fermi level. We also observe low-energy magnetic excitations evolving across the low-temperature transition, largely consistent with the FB shift. Our results suggest that the low-temperature order contains a magnetic origin and that the kagome FBs may play a role in the emergence of this order.

electronic properties and materials↗
DOE OSTI2022

Vibrationally resolved optical excitations of the nitrogen-vacancy center in diamond

A comprehensive description of the optical cycle of spin defects in solids requires the understanding of the electronic and atomistic structure of states with different spin multiplicity, including singlet states which are particularly challenging from a theoretical standpoint. We present a general framework, based on spin-flip time-dependent density function theory, to determine the excited state potential energy surfaces of the many-body singlet states of spin defects; we then predict the vibrationally resolved absorption spectrum between singlet shelving states of a prototypical defect, the nitrogen-vacancy center in diamond. Our results, which are in very good agreement with experiments, provide an interpretation of the measured spectra and reveal the key role of specific phonons in determining absorption processes, and the notable influence of non-adiabatic interactions. The insights gained from our calculations may be useful in defining strategies to improve infrared-absorption-based magnetometry and optical pumping schemes. The theoretical framework developed here is general and applicable to a variety of other spin defects and materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

First-Principles insights into group-V impurities and their impact on germanium detector performance

The outstanding properties of high-purity germanium (HPGe) detectors, such as excellent energy resolution, high energy sensitivity, and a low background-to-signal ratio, make them essential and ideal candidates for detecting particle signatures in nuclear processes such as neutrino-less double beta decay (0νββ). However, the presence of defects and impurities in HPGe crystals can lead to charge trapping, which affects carrier mobility and results in significant energy resolution degradation. In this work, we employ density functional theory with a hybrid functional to study the energetics of possible point defects in Ge. Our findings indicate that group-V impurities form more readily in Ge compared to vacancy and interstitial of Ge. Unlike N dopants, which yield deep trap states, P, As, and Sb create shallow traps close to the conduction band edge of Ge. Furthermore, we predict that group-V defects can condense into defect complexes with Ge vacancies. These vacancy-impurity complexes form deep traps in Ge, similar to Ge vacancies, suggesting that both vacancies and vacancy-impurity complexes contribute to charge trapping in these detectors, thereby diminishing their performance.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly fluorescent purine-containing conjugated copolymers with tailored optoelectronic properties

Conjugated copolymers containing electron donor and acceptor units in their main chain have emerged as promising materials for organic electronic devices due to their tunable optoelectronic properties. In this work, we describe the use of direct arylation polymerization to create a series of fully π-conjugated copolymers containing the highly tailorable purine scaffold as a key design element. To create efficient coupling sites, dihalopurines are flanked by alkylthiophenes to create a monomer that is readily copolymerized with a variety of conjugated comonomers, ranging from electron-donating 3,4-dihydro-2H-thieno[3,4-b][1,4]dioxepine to electron-accepting 4,7-bis(5-bromo-3-hexylthiophen-2-yl)benzo[c][1,2,5]thiadiazole. The comonomer choice and electronic nature of the purine scaffold allow the photophysical properties of the purine-containing copolymers to be widely varied, with optical bandgaps ranging from 1.96–2.46 eV, and photoluminescent quantum yields as high as Φ = 0.61. Frontier orbital energy levels determined for the various copolymers using density functional theory tight binding calculations track with experimental results, and the geometric structures of the alkylthiophene-flanked purine monomer and its copolymer are found to be nearly planar. The utility of direct arylation polymerization and intrinsic tailorability of the purine scaffold highlight the potential of these fully conjugated polymers to establish structure–property relationships based on connectivity pattern and comonomer type, which may broadly inform efforts to advance purine-containing conjugated copolymers for various applications.

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

AENET–LAMMPS and AENET–TINKER : Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials

Machine-learning potentials (MLPs) trained on data from quantum-mechanics based first-principles methods can approach the accuracy of the reference method at a fraction of the computational cost. To facilitate efficient MLP-based molecular dynamics and Monte Carlo simulations, an integration of the MLPs with sampling software is needed. Here, we develop two interfaces that link the atomic energy network (ænet) MLP package with the popular sampling packages TINKER and LAMMPS. The three packages, ænet, TINKER, and LAMMPS, are free and open-source software that enable, in combination, accurate simulations of large and complex systems with low computational cost that scales linearly with the number of atoms. Scaling tests show that the parallel efficiency of the ænet–TINKER interface is nearly optimal but is limited to shared-memory systems. The ænet–LAMMPS interface achieves excellent parallel efficiency on highly parallel distributed memory systems and benefits from the highly optimized neighbor list implemented in LAMMPS. We demonstrate the utility of the two MLP interfaces for two relevant example applications: the investigation of diffusion phenomena in liquid water and the equilibration of nanostructured amorphous battery materials.

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

Review of radiation-induced defects in GaAs

Radiation-resilient optoelectronic materials are highly desired in space and nuclear applications, and understanding the relevant defects and electronic processes in those materials is crucial for both current and next-generation applications. GaAs is a prototypical semiconductor that serves as a foundational system for describing and understanding optoelectronic devices. In this review, both experimental studies and molecular dynamics (MD) simulations of irradiated GaAs are reviewed, with particular emphasis on the deep-level transient spectroscopy, irradiation-induced amorphization zones, and MD predictions of damage structures. The MD results are also compared to predictions of the non-ionized energy-loss model. Recent theoretical studies, in particular, density functional theory based calculations on the simple intrinsic defects and defect clusters in GaAs, are also reviewed. These defects have an important role in dictating the evolution of GaAs in radiation damage environments, as they impact the coupled dynamics of charge carriers. Finally, possible gaps and challenges toward the general understanding of defect evolution in GaAs are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ZENN: A thermodynamics-inspired computational framework for heterogeneous data–driven modeling

Traditional entropy-based methods—such as cross-entropy loss in classification problems—have long been essential tools for representing the information uncertainty and physical disorder in data and for developing artificial intelligence algorithms. However, the rapid growth of data across various domains has introduced new challenges, particularly the integration of heterogeneous datasets with intrinsic disparities. To address this, we introduce a zentropy-enhanced neural network (ZENN), extending zentropy theory into the data science domain via intrinsic entropy, enabling more effective learning from heterogeneous data sources. ZENN simultaneously learns both energy and intrinsic entropy components, capturing the underlying structure of multisource data. To support this, we redesign the neural network architecture to better reflect the intrinsic properties and variability inherent in diverse datasets. We demonstrate the effectiveness of ZENN on classification tasks and energy landscape reconstructions, showing its superior generalization capabilities and robustness-particularly in predicting high-order derivatives. In image and text classification tasks, ZENN demonstrates superior generalization by introducing a learnable temperature variable that models latent multisource heterogeneity, allowing it to surpass state-of-the-art models on CIFAR-10/100, BBC News, and AG News. As a practical application in materials science, we employ ZENN to reconstruct the Helmholtz energy landscape of Fe3Pt using data generated from density functional theory and capture key material behaviors, including negative thermal expansion and the critical point in the temperature–pressure space. Overall, this work presents a zentropy-grounded framework for data-driven machine learning, positioning ZENN as a versatile and robust approach for scientific problems involving complex, heterogeneous datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Diabatic Hamiltonian matrix elements made simple

With a view to applying the generator coordinate method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. Here, the method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Understanding Cu incorporation in the Cu 2x Hg 2-x GeTe 4 structure using resonant x-ray diffraction

The ability to control carrier concentration based on the extent of Cu solubility in the Cu 2x Hg 2-x GeTe 4 alloy compound (where 0 ≤ x ≤ 1) makes Cu 2x Hg 2-x GeTe 4 an interesting case study in the field of thermoelectrics. While Cu clearly plays a role in this process, it is unknown exactly how Cu incorporates into the Cu 2x Hg 2-x GeTe 4 crystal structure and how this affects the carrier concentration. In this work, we use a combination of resonant energy x-ray diffraction (REXD) experiments and density functional theory (DFT) calculations to elucidate the nature of Cu incorporation into the Cu 2x Hg 2-x GeTe 4 structure. REXD across the Cu k edge facilitates the characterization of Cu incorporation in the Cu 2x Hg 2-x GeTe 4 alloy and enables direct quantification of antisite defects. Here, we find that Cu substitutes for Hg at a 2:1 ratio, wherein Cu annihilates a vacancy and swaps with a Hg atom. DFT calculations confirm this result and further indicate that the incorporation of Cu occurs preferentially on one of the z = 1/4 or z = 3/4 planes before filling the other plane. Furthermore, the amount of Cu Hg antisite defects quantified by REXD was found to be directly proportional to the experimentally measured hole concentration, indicating that the Cu Hg defects are the driving force for tuning carrier concentration in the Cu 2x Hg 2-x GeTe 4 alloy. The link uncovered here between crystal structure, or more specifically antisite defects, and carrier concentration can be extended to similar cation-disordered material systems and will aid the development of improved thermoelectric and other functional materials through defect engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗