Search NASASearch

SEARCH · Search NASA

Results for “density functional calculations”

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.

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory

Carbon in GaN as a nonradiative recombination center

Trap-assisted nonradiative recombination has been shown to limit the efficiency of optoelectronic devices. While substitutional carbon (⁠C N ) has been suggested to be a nonradiative recombination center in GaN devices, a complete recombination cycle including the two charge-state transition levels has not been previously described. In this work, we investigate the trap-assisted recombination process due to C N in GaN, including multiphonon emission, radiative recombination, trap-assisted Auger–Meitner (TAAM) recombination as well as thermal emission of holes. Our study shows the key role of TAAM processes at the high carrier densities relevant for devices. We also reveal the carrier-density regimes where thermal emission and radiative recombination are expected to play an observable role. Our results highlight that carbon concentrations exceeding ~10 17 cm −3 can have a noticeable impact on device efficiency, not just in GaN active layers but also in InGaN and AlGaN. Furthermore, our comprehensive formalism not only offers detailed results for carbon but also provides a general framework for assessing the multiple processes that participate in trap-assisted recombination in semiconductors.

36 MATERIALS SCIENCE

Electronic structure of the kagome compound CaTi 3⁢ Bi 4 using high-field torque magnetometry and density functional theory

Here, we report systematic torque magnetometry measurements to investigate the electronic properties of the newly discovered kagome compound CaTi 3 ⁢Bi 4 . Electrical transport, magnetic susceptibility, and thermal measurements reveal no evidence of a magnetic ground state in this material. Torque data obtained in magnetic fields up to 41.5 T exhibit clear de Haas–van Alphen (dHvA) oscillations, with nine distinct frequencies ranging from 13 to 6164 T. Angular-dependent dHvA measurements show that, with the exception of the lowest frequency (13 T), all observed frequencies nearly follow a 1/cos ⁡𝜃 dependence, where 𝜃 is the angle between the crystallographic 𝑐 axis and the magnetic field direction. This behavior is characteristic of quasi-two-dimensional Fermi-surface sheets with nearly circular cross sections. To further elucidate the electronic structure, we performed density functional theory (DFT) calculations of the band structure and Fermi surface. The calculated bands reveal the presence of multiple Dirac points (DP), flat bands (FB), and van Hove singularities (VHS) near the Fermi level. The resulting Fermi surface consists of several quasi-two-dimensional cylindrical sheets, consistent with the experimentally observed 1/cos⁡ 𝜃 dependence. Notably, the theoretical dHvA frequencies, derived from extremal Fermi-surface cross-sectional areas, agree well with the experimental values and reproduce their angular dependence. Remarkably, all experimentally observed frequencies are captured by the DFT predictions. Pressure-dependent calculations up to 10 GPa show that the electronic features—DP, FB, and VHS—evolve systematically with pressure. In particular, the VHS shifts closer to the Fermi level, demonstrating that pressure acts as an effective tuning parameter in this material. These combined experimental and theoretical results provide a comprehensive understanding of the electronic structure of CaTi 3 ⁢Bi 4 and demonstrate how pressure can be used to tune its key electronic features, offering valuable guidance for exploring related kagome materials.

Shtefiienko, Kyryl [West Texas A & M University, C

Density Functional Theory Analysis of Candidate Tritium Permeation Barriers for Vanadium Alloys as a Fusion Structural Material

Vanadium alloys are being considered as a candidate structural material for the breeder blanket component of future fusion reactors due to their good thermal conductivity, low activation, and chemical compatibility with FLiBe. However, vanadium has a high hydrogen solubility, resulting in hydrogen diffusion and trapping in the vanadium structural alloy and lower than expected tritium breeding ratios. To mitigate the permeation of hydrogen into the structural alloy, a tritium permeation barrier will be coated onto the vanadium alloy. Multiple materials have been proposed as candidate tritium permeation barriers, but the extent to which each of these materials impacts hydrogen diffusion and trapping is still unknown. In this project, we perform atomic-scale calculations to narrow down the list of potential tritium permeation barrier candidates and to understand the mechanisms through which hydrogen diffusion is affected in these materials.

36 MATERIALS SCIENCE

Three-dimensional dispersion in the type-II Dirac semimetals PtTe 2 and PdTe 2 revealed through circular dichroism in angle-resolved photoemission spectroscopy

PtTe 2 and PdTe 2 are among the first transition metal dichalcogenides that were predicted to host type-II Dirac fermions, exotic particles prohibited in free space. These materials are layered and air stable, which makes them top candidates for technological applications that take advantage of their anisotropic magnetotransport properties. Here, in this work, we provide a detailed characterization of the electronic structure of PtTe 2 and PdTe 2 using angle-resolved photoemission spectroscopy (ARPES) and density functional theory calculations, offering an alternative interpretation for one of the Dirac-like dispersions in these materials. Through the use of circularly polarized light, we report a different behavior of such dispersion in PdTe 2 compared to PtTe 2 , that we relate to a symmetry analysis of the dipole matrix element. Such analysis reveals a link between the observed circular dichroism and the different momentum-dependent terms in the dispersion of these two compounds, despite their close similarity in crystal structure. Additionally, our data show a clear difference in the circular dichroic signal for the type-II Dirac cones characteristic of these materials, compared to their topologically protected surface states. Our paper provides a useful reference for the ARPES characterization of other transition metal dichalcogenides with topological properties and illustrates the use of circular dichroism as a guide to identify the topological character and attributes of two otherwise equivalent band dispersions.

angle-resolved photoemission spectroscopy

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems

Derivation of low-energy Hamiltonians for heavy-fermion materials

Here, by utilizing a multiorbital periodic Anderson model with parameters obtained from ab initio band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both nonmagnetic 4⁢𝑓 0 and magnetic 4⁢𝑓 2 virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn 3 . The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of 4⁢𝑓 materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Thermally Stable Co@C3N4 Single-Atom Catalysts for CO Oxidation: Atomic-Level Insights into Structure and Activity

Single Co atoms supported on C3N4 (Co@C3N4) have demonstrated high activity and selectivity in photocatalysis. However, the investigation of structure–function relationships and reaction mechanisms under photocatalytic conditions is very challenging due to the complex conditions of light absorption, charge transfer, and catalysis. In this study, we employed thermal CO oxidation as a prototypical probe reaction to benchmark the intrinsic catalytic performance and track the active-site evolution of Co@C3N4. Single Co atoms were identified and shown to be the catalytically active sites for CO oxidation based on control experiments and isotope-labeling experiments. The Co sites remained atomically dispersed before, during, and after the reaction with temperatures up to 400 °C, as established by in situ X-ray absorption fine structure (XAFS) combined with density functional theory (DFT), FDMNES simulations, and dynamic-time-warping (DTW)-assisted X-ray absorption near edge structure (XANES) matching. Together with theoretical calculations, the integrated analysis reveals a stable coordination environment under reaction conditions, which correlates with sustained activity, establishing Co@C3N4 single-atom catalysts as thermally stable CO oxidation catalysts. Beyond these findings, the current study provides a workflow for unambiguously assigning active sites in Co@C3N4 for thermal CO oxidation. This workflow will aid the understanding of their behavior in photocatalysis in the future, where light-driven dynamics obscure direct structure–function links. Notably, this study provides fundamental insights for the rational design of robust single-atom catalysts and a foundation for the broader application of Co@C3N4 catalysts in oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics

Propane Activation on Pt Electrodes at Room Temperature: Quantification of Adsorbate Identity and Coverage

Abstract C−H bond activation is the first step in manufacturing chemical products from readily available light alkane feedstock and typically proceeds via carbon‐intensive thermal processes. The ongoing emphasis on decarbonization via electrification motivates low‐temperature electrochemical alternatives that could lead to sustainable chemicals production. Platinum (Pt) electrocatalysts have shown activity towards reacting alkanes; however, little is known about propane electrocatalytic activation and conditions suitable for enabling selective oxidation to valuable products. Herein, we utilize a combination of electrochemical mass spectrometry (ECMS) and density functional theory (DFT) calculations to elucidate the potential dependence of propane activation on Pt electrocatalysts. Results show a strong dependence of adsorption on the applied potential in room‐temperature aqueous acidic electrolyte, with a maximum coverage of propane‐derived adsorbates at 0.30 V vs RHE. Using charge deconvolution and deuterated experiments, the mechanism of adsorption was elucidated, and C 3 H 2 * was determined as the average dehydrogenated propane‐derived adsorbate species. DFT calculations further corroborate these results, showing that the formation of deeply dehydrogenated species is energetically accessible at room temperature. The combined theoretical and experimental findings yield insights for selective activation of paraffinic C−H bonds at room temperature, aqueous conditions—a critical step towards decarbonized chemical manufacturing.

Chemistry

Microscopic Scattering Approach to In-Gap States

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of beta Bi2Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov-de Gennes formalism[1]. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic dIdV maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate C4v representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr d-Bi p hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data for magnetic in-gap states on superconductors. Further advances examining layered 2D material surfaces, such as NbSe2, will be described[2]. For this system the superconducting properties are obtained from a full anisotropic Eliashberg calculation of the superconducting order parameter along with the charge density wave gap. Additional features associated with proposals to measure the dynamics of these individual YSR states will be presented. [1] arXiv:2507.08740 [2] arXiv:2507.11856

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND