Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory 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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermodynamics and kinetics of H adsorption and intercalation for graphene on 6 H -SiC(0001) from first-principles calculations

Previous experimental observations for H intercalation under graphene on SiC surfaces motivate the clarification of configuration stabilities and kinetic processes related to intercalation. From first-principles density-functional-theory calculations, we analyze H adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that in the low-H coverage regime (described by a single H atom within a sufficiently large supercell), intercalation into the gallery between TLG and BLG or into the gallery underneath BLG is more favorable thermodynamically than adsorption on top of TLG. However, intercalation into the gallery between TLG and BLG is most favorable. We obtain energy barriers of about 1.3 and 2.3 eV for a H atom diffusing on and under TLG, respectively. From an additional analysis of the energy landscape in the vicinity of a step on the TLG, we assess how readily one guest H atom on the TLG terrace can directly penetrate the TLG into the gallery between TLG and BLG versus crossing a TLG step to access the gallery. Finally, we also perform density functional theory calculations for higher H coverages revealing a shift in favorability to intercalation of H underneath BLG and characterizing the variation with H coverage in interlayer spacings.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tuning the reactivity of carbon surfaces with oxygen-containing functional groups

Oxygen-containing carbons are promising supports and metal-free catalysts for many reactions. However, distinguishing the role of various oxygen functional groups and quantifying and tuning each functionality is still difficult. Here we investigate the role of Brønsted acidic oxygen-containing functional groups by synthesizing a diverse library of materials. By combining acid-catalyzed elimination probe chemistry, comprehensive surface characterizations, 15N isotopically labeled acetonitrile adsorption coupled with magic-angle spinning nuclear magnetic resonance, machine learning, and density-functional theory calculations, we demonstrate that phenolic is the main acid site in gas-phase chemistries and unexpectedly carboxylic groups are much less acidic than phenolic groups in the graphitized mesoporous carbon due to electron density delocalization induced by the aromatic rings of graphitic carbon. The methodology can identify acidic sites in oxygenated carbon materials in solid acid catalyst-driven chemistry.

Zhou, Jiahua↗
DOE OSTI2025

Phonon Olympics: Phonon property and lattice thermal conductivity benchmarking from open-source packages

Three widely used open-source packages for determining phonon properties and lattice thermal conductivities (ALAMODE, phono3py, and ShengBTE) are benchmarked by teams of expert users and the package developers. The phonons for Ge, RbBr, monolayer MoSe 2 , and AlN are modeled at zero temperature, and they scatter through three-phonon and phonon-isotope processes, with thermal conductivities obtained from the linearized Peierls–Boltzmann transport equation with input from density functional theory calculations. Over a wide range of temperatures, the thermal conductivities calculated by the teams fall within at most ±15% of their mean values for each of the four materials. The phonon frequencies, obtained from the harmonic force constants, do not show large differences between the calculations, indicating that the modal heat capacities and group velocities are not responsible for the thermal conductivity variations. It is the lifetimes associated with three-phonon scattering, obtained from the cubic force constants, that drive the variations. The many decisions required to calculate the cubic force constants (e.g., supercell size, atomic displacement, neighbor cutoff, and application of symmetries) make identification of the precise origin of the thermal conductivity variations challenging. The calculated thermal conductivities do not generally show agreement with experimental measurements, which is attributed to the limitations of the density functional theory calculations. Guidance for the development of best practices is provided, which will help to standardize protocols needed for building thermal conductivity databases. The results provide a baseline for future benchmarking of other packages and more advanced calculations.

McGaughey, Alan J. H. [Carnegie Mellon Univ., Pitt↗
DOE OSTI2021

Finite-size correction for slab supercell calculations of materials with spontaneous polarization

Abstract The repeated slab approach has become a de facto standard to accurately describe surface properties of materials by density functional theory calculations with periodic boundary conditions. For materials exhibiting spontaneous polarization, we show that the conventional scheme of passivation with pseudo hydrogen is unable to realize a charge-neutral surface. The presence of a net surface charge induces via Gauss’s law a macroscopic electric field through the slab and results in poor size convergence with respect to the thickness of the slab. We propose a modified passivation method that accounts for the effect of spontaneous polarization, describes the correct bulk limits and boosts convergence with respect to slab thickness. The robustness, reliability, and superior convergence of energetics and electronic structure achieved by the proposed method are demonstrated using the example of polar ZnO surfaces.

Chemistry↗
DOE OSTI2023

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Confinement of quasi-atomic structures in Ti 2 N and Ti 3 N 2 MXene electrides

Metal carbides, nitrides, or carbonitrides of early transition metals, better known as MXenes, possess notable structural, electrical, and magnetic properties. Analyzing electronic structures by calculating structural stability, band structure, density of states, Bader charge transfer, and work functions utilizing first principle calculations, we revealed that titanium nitride MXenes, namely Ti 2 N and Ti 3 N 2 , have excess anionic electrons in their lattice voids, making them MXene electrides. Bulk Ti 3 N 2 has competing antiferromagnetic (AFM) and ferromagnetic (FM) configurations with slightly more stable AFM configuration, while the Ti 2 N MXene is nonmagnetic. Although Ti 3 N 2 favors AFM configuration with hexagonal crystal systems having 6/ mmm point group symmetry, Ti 3 N 2 does not support altermagnetism. The monolayer of the Ti 3 N 2 MXene is a ferromagnetic electride. These unique properties of having non-nuclear interstitial anionic electrons in the electronic structure of titanium nitride MXene have not yet been reported in the literature. Density functional theory calculations show TiN is neither an electride, MXene, or magnetic.

Anionic electrons↗
DOE OSTI2026

Ultrafast One-Dimensional Peierls-Distortion Dynamics in 1⁢T′−Re⁢S 2 Revealed by 4D Electron Microscopy

Rhenium disulfide (ReS 2 ), a prototypical 2D semiconductor with an anisotropic 1⁢T′ structure due to the pronounced Peierls distortion, has demonstrated great potential for polarization-dependent optoelectronics and photonics. Here, we report an ultrafast phase transition occurring within 1 ps, accompanied by a one-dimensional Peierls-distortion relaxation in 1⁢T′−ReS 2 revealed with combined 4D electron microscopy and time-dependent density functional theory calculations. Upon femtosecond laser pulse excitation, the Re-Re dimerization morphology rapidly transforms from a diamond cluster to zigzag chains and persists for several nanoseconds. Here, this ultrafast Peierls-distortion relaxation is further verified by transient changes in optical anisotropy via polarization-dependent transient absorption spectroscopy. Time-dependent density functional theory calculations attribute this novel phase transition to strong correlation between Peierls distortion and impulsive photoexcited carrier doping, predicting a transient band-gap collapse. This Letter opens up an exciting avenue for ultrafast control of Peierls-distortion-induced anisotropy and the metal-insulator transition in 1⁢T′−ReS 2 .

1T’-ReS2↗
DOE OSTI2021

A Polar Magnetic and Insulating Double Corundum Oxide: Mn2MnSbO6 with Ordered Mn(II) and Mn(III) Ions

A new magnetic insulator Mn2MnSbO6 with a polar crystal structure and an ordered Mn2+ and Mn3+ arrangement was synthesized under a high pressure of 7.5 GPa and 1300 degrees C. The crystal structure of Mn2MnSbO6, investigated by synchrotron powder X-ray diffraction, was found to be isomorphous with that of Ni3TeO6-type, space group R3. The non-centrosymmetric structure was confirmed by the second-harmonic generation measurements. The X-ray absorption near-edge spectroscopy measurement confirmed the nominal oxidation states of Mn22+Mn3+SbO6. Magnetic measurements indicate that Mn2MnSbO6 orders antiferromagnetically below 44 K and undergoes a field-induced spin-flop transition at 5 K. First-principles calculations indicate an antiferromagnetic ground state with up/up/up/down/down/down (uuuddd) spin configuration of the six crystallographically unique Mn ions in the c-axis doubled magnetic structure. The density functional theory calculations also substantiate the experimentally observed charge ordering of the Mn2+/Mn3+ ions and the insulating behavior due to a bandgap of 0.52 eV. To the best of our knowledge, this is the first double corundum oxide containing Jahn-Teller active Mn3+ ions.

Feng, Hai L.↗
DOE OSTI2022

Strongly correlated itinerant magnetism near superconductivity in NiTa 4 Se 8

Metallic ferromagnets with strongly interacting electrons often exhibit remarkable electronic phases such as ferromagnetic superconductivity, complex spin textures, and nontrivial topology. Here, in this report, we discuss the synthesis of a layered magnetic metal NiTa 4 Se 8 (or Ni 1/4 TaSe 2 ) with a Curie temperature of 58 Kelvin. Magnetization data and density functional theory calculations indicate that the nickel atoms host uniaxial ferromagnetic order of about 0.7µ B per atom, while an even smaller moment is generated in the itinerant tantalum conduction electrons. Strong correlations are evident in flat bands near the Fermi level, a high heat capacity coefficient, and a high Kadowaki-Woods ratio. Density functional theory calculations suggest that electron and hole Fermi surfaces in the ferromagnetic phase are associated with opposite spin polarization. When the system is diluted of magnetic ions, the samples become superconducting below about 2 Kelvin. We discuss possible mechanisms for superconductivity in this family.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

Single-crystal growth and superconductivity in RbNi 2 Se 2

Herein we report the synthesis and characterization of RbNi 2 Se 2 , an analog of the iron chalcogenide superconductor Rb x Fe 2 Se 2 , via transport, angle-resolved photoemission spectroscopy, and density functional theory calculations. A superconducting transition at T c =1.20 K is identified. In the normal state, RbNi 2 Se 2 shows paramagnetic and Fermi-liquid behaviors. A large Sommerfeld coefficient yields an effective electron mass of m* ≈ 6m e . In the superconducting state, zero-field electronic specific-heat data C es can be described by a two-gap BCS model, indicating that RbNi 2 Se 2 is a possible multigap superconductor. Our density functional theory calculations and angle-resolved photoemission spectroscopy measurements demonstrate that RbNi 2 Se 2 exhibits relatively weak correlations and multiband characteristics, consistent with the multigap superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Probing Active-Site Relocation in Cu/SSZ-13 SCR Catalysts during Hydrothermal Aging by in situ EPR Spectroscopy, Kinetic Studies, and DFT Calculations

Cu/SSZ-13 selective catalytic reduction catalysts display activity loss after mild hydrothermal aging without support degradation or loss of isolated Cu ion content. By applying electron paramagnetic resonance spectroscopy to compare hydrated, dehydrated and NH3-saturated catalysts before and after hydrothermal aging, as well as density functional theory calculations, it is found that hydrothermal aging induces Cu relocation, leading to stronger Cu-support interactions. This is particularly evidenced by ZCuOH conversion to Z2Cu, which leads to the generation of a unique structure with two Cu(II) ions positioned in a double 6-membered ring prism. The Cu(II) ions that do not undergo such dramatic changes also appear to stay closer to framework windows as evidenced by the overall increased anisotropy during EPR measurements. In situ EPR measurements under SCR conditions demonstrate that, in comparison to the fresh catalyst, higher percentages of Cu species stay as Cu(II) in the aged catalyst. This suggests that the reduction half-cycle of the SCR redox, i.e., Cu(II) ? Cu(I), becomes slower for the hydrothermally aged catalyst.

42 ENGINEERING↗
DOE OSTI2024

Design and Synthesis of PtPdNiCoMn High‐Entropy Alloy Electrocatalyst for Enhanced Alkaline Hydrogen Evolution Reaction: A Theoretically Supported Predictive Design Approach

Electrocatalytic hydrogen generation requires a multifunctional electrocatalyst with abundant active sites to drive multielectron transfer reactions. High entropy alloys (HEA) are five or more-elements with high configurational entropy are considered unique materials for next-generation electrocatalysts. Here, in this work, based on new screening guidelines for catalyst selections that combine density-functional theory calculated Gibbs formation-enthalpy with bond length and electronegativity variance, a novel HEA electrocatalyst consisting of five elements, namely, Pt, Pd, Ni, Co, and Mn has been designed. By simple room temperature electrodeposition, the designed catalyst is prepared and its hydrogen evolution reaction (HER) is explored and validated through experimental and theoretical approaches. The HEA demonstrated a superior HER activity with an overpotential of 22.6 mV at -10 mA cm -2 which outperforms Pt/C commercial catalyst. No evident degradation of the material is detected even after 100 hours of continuous operation under high current density. Moreover, the HEA has shown exceptional performance in harsh electrolyte conditions such as in simulated seawater and actual seawater. Remarkably, the density-functional theory calculated Gibbs formation-enthalpy is small (≈0 eV) compared to Pt/C placing the new HEA near the apex of Trasatti's model of Volcano plot, which is also suggestive of superior HER activity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Angular Momentum in Rotating Superfluid Droplets

The angular momentum of rotating superfluid droplets originates from quantized vortices and capillary waves, the interplay between which remains to be uncovered. Here, the rotation of isolated submicrometer superfluid 4 He droplets is studied by ultrafast x-ray diffraction using a free electron laser. The diffraction patterns provide simultaneous access to the morphology of the droplets and the vortex arrays they host. In capsule-shaped droplets, vortices form a distorted triangular lattice, whereas they arrange along elliptical contours in ellipsoidal droplets. The combined action of vortices and capillary waves results in droplet shapes close to those of classical droplets rotating with the same angular velocity. The findings are corroborated by density functional theory calculations describing the velocity fields and shape deformations of a rotating superfluid cylinder.

O'Connell, Sean M.O.↗
DOE OSTI2020

A Comparative Study of Redox Mediators for Improved Performance of Li-Oxygen Batteries.

Redox meditators (RMs) are soluble catalysts located in an electrolyte that can improve the energy efficiency (reduced overpotential) and cyclability of Li-oxygen (Li-O-2) batteries. In this work, 20 RMs within a Li-O-2 system with dimethyl sulfoxide and tetraethylene glycol dimethyl ether electrolytes are studied and their electrochemical features such as redox potential, the separation of cathodic and anodic peaks, and their current intensities are measured using cyclic voltammetry (CV) experiments. Six RMs are selected as "primary" choices based on their electrochemical performance, and stability tests are then performed to examine their electrochemical responses after consecutive cycles. Moreover, galvanostatic cycling tests are performed within a Li-O-2 battery system assembled with selected six RMs for real case consistency investigations. It is found that results from CV to galvanostatic cycling tests are consistent for halides and organometallic RMs, where the former exhibit much higher stability. However, the organic RMs show high reversibility in CV but low in battery cycling results. Density functional theory calculations are carried out to gain more understanding of the stability and redox potentials of the RMs. This study provides comparative information to select the most reliable RMs for Li-O-2 batteries along with new fundamental understanding of their electrochemical activity and stability.

cyclic voltammetry↗
DOE OSTI2020

Computationally aided, entropy-driven synthesis of highly efficient and durable multi-elemental alloy catalysts

Multi-elemental alloy nanoparticles (MEA-NPs) hold great promise for catalyst discovery in a virtually unlimited compositional space. However, rational and controllable synthesize of these intrinsically complex structures remains a challenge. Here, we report the computationally aided, entropy-driven design and synthesis of highly efficient and durable catalyst MEA-NPs. The computational strategy includes prescreening of millions of compositions, prediction of alloy formation by density functional theory calculations, and examination of structural stability by a hybrid Monte Carlo and molecular dynamics method. Selected compositions can be efficiently and rapidly synthesized at high temperature (e.g., 1500 K, 0.5 s) with excellent thermal stability. We applied these MEA-NPs for catalytic NH 3 decomposition and observed outstanding performance due to the synergistic effect of multi-elemental mixing, their small size, and the alloy phase. We anticipate that the computationally aided rational design and rapid synthesis of MEA-NPs are broadly applicable for various catalytic reactions and will accelerate material discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomically Precise PdSe 2 Pentagonal Nanoribbons

We report atomically precise pentagonal PdSe 2 nanoribbons (PNRs) fabricated on a pristine PdSe 2 substrate with a hybrid method of top-down and bottom-up processes. The PNRs form a uniform array of dimer structure with a width of 2.4 nm and length of more than 200 nm. In situ four-probe scanning tunneling microscopy (STM) reveals metallic behavior of PNRs with ballistic transport for at least 20 nm in length. Density functional theory calculations produce a semiconducting density of states of isolated PNRs and find that the band gap narrows and disappears quickly once considering coupling between PNR stacking layers or interaction with the PdSe 2 substrate. The coupling of PNRs is further corroborated by Raman spectroscopy and field-effect transistor measurements. The facile method of fabricating atomically precise PNRs offers an air-stable functional material for dimensional control.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nanoconfinement of Molecular Magnesium Borohydride Captured in a Bipyridine-Functionalized Metal-Organic Framework

The lower limit of metal hydride nanoconfinement is demonstrated through the coordination of a molecular hydride species to binding sites inside the pores of a metal-organic framework (MOF). Magnesium borohydride, which has a high hydrogen capacity, is incorporated into the pores of UiO-67bpy (Zr6O4(OH)4(bpydc)6 with bpydc2- = 2,2'-bipyridine-5,5'- dicarboxylate) by solvent impregnation. The MOF retained its long-range order, and transmission electron microscopy and elemental mapping confirmed the retention of the crystal morphology and revealed a homogeneous distribution of the hydride within the MOF host. Notably, the B-, N-, and Mg-edge XAS data confirm the coordination of Mg(II) to the N atoms of the chelating bipyridine groups. In situ 11B MAS NMR studies helped elucidate the reaction mechanism and revealed that complete hydrogen release from Mg(BH4)2 occurs as low as 200 °C. Sieverts and thermogravimetric measurements indicate an increase in the rate of hydrogen release, with the onset of hydrogen desorption as low as 120 °C, which is approximately 150 °C lower than that of the bulk material. Furthermore, density functional theory calculations support the improved dehydrogenation properties and confirm the drastically lower activation energy for B-H bond dissociation.

Nanoconfinement, Metal Hydrides, Metal-Organic Fra↗