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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Preferential hole defect formation in monolayer W Se 2 by electron-beam irradiation

Monolayer transition-metal dichalcogenides (TMDCs) have been extensively studied due to their wide range of physical properties and applications. It has been demonstrated that the electron beam in a transmission electron microscope (TEM) or scanning TEM (STEM) generates chalcogen vacancies that agglomerate into dispersed linelike or holelike defects. Here we employ a STEM electron beam and demonstrate that, in WSe 2 , beam-induced chalcogen vacancies initially form certain dispersed multivacancy structures, as seen in TMDCs in prior work. However, with suitable control of the STEM focused electron beam, these multivacancies gradually evolve into a dense network of ten-, 12-, 14-, and 16-member ring round holes, whereas the same process leads predominantly to chalcogen-vacancy line defects in other trigonal-prismatic TMDCs. Density functional theory calculations find that the underlying atomic-scale processes lead preferentially to defect structures that lower the total energy so that we are able to track the formation of the observed multivacancy complexes, which then lead to the formation of dense large round holes in WSe 2 . The same processes in WS 2 , however, are energetically unfavorable, while linear multivacancy defects are preferred, as observed. The demonstrated control of the formation of unique high-density round holes in WSe 2 has potential for applications such as atomic and molecular sieving.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Understanding Inlet Concentration Effects on the Electrocatalytic Conversion of CO 2 to Formic Acid in Gas-Fed Electrolyzers

The electrochemical CO 2 reduction reaction (CO2RR) to produce value-added products remains a developing technology for utilizing waste CO 2 streams. Most device-level CO2RR studies use pure CO 2 gas feeds; however, the effect of dilute CO 2 on the electrolyzer performance is an important consideration for large-scale electrolyzer operation, single-pass conversion, and real-world CO 2 source utilization. This work investigates the effect that the CO 2 concentration has on the performance of formic acid (HCOOH) producing tin oxide (SnO 2 ) and bismuth oxide (Bi 2 O 3 ) catalysts in an electrolyzer device setting. Surprisingly, SnO2 demonstrated an approximately 20% increase in HCOOH selectivity (Faradaic efficiency) when the CO 2 concentration decreased from 100 to 20%. In contrast, Bi 2 O 3 consistently demonstrated high selectivity toward HCOOH across the same CO 2 concentration range. The effects of the CO 2 concentration on selectivity were further investigated with half-cell experiments and in situ Raman spectroscopy, which revealed dynamic changes in the cathodic overpotential and chemical state of the catalyst that depended on the CO 2 concentration. Density functional theory calculations showed how changes in the surface oxidation state of Sn, varying from fully oxidized SnO 2 to metallic Sn(0), affect the thermodynamic barriers of the three main observed products: HCOOH, CO, and H 2 . Our results indicate that dilute CO 2 concentrations required larger cathodic overpotentials to sustain a fixed current density, which, in turn, pushed the Sn-based catalyst toward a more reduced surface that was favorable to HCOOH formation. On the other hand, the Bi-based catalyst remained in a metallic state at CO2RR-relevant potentials and demonstrated a consistent product selectivity regardless of CO 2 concentration. These findings highlight how varying the CO 2 inlet gas concentrations affects the chemical state of catalysts and the resulting performance metrics.

42 ENGINEERING↗
DOE OSTI2022

Meta-GGA exchange-correlation free energy density functional to increase the accuracy of warm dense matter simulations

We discuss strategies for thermalization of the ground-state meta-generalized gradient approximation (meta-GGA) exchange-correlation (XC) functionals. A simple but accurate scheme is implemented via universal additive thermal correction to XC using a perturbative-like self-consistent approach. The additive correction with explicit temperature dependence is applied to the ground-state deorbitalized, strongly constrained and appropriately normed (SCAN-L) meta-GGA XC leading to thermal XC functional denoted here as T-SCAN-L. Thermal T-SCAN-L meta-GGA functional shows significant improvement in density functional theory calculation accuracy for warm dense matter by a factor of 3 to 10, achieving unprecedented accuracy of total pressure between a few tenths and 1% when compared to traditional XC functionals, as demonstrated by the comparison to pathintegral Monte Carlo simulations for helium equation of state. Furthermore, the T-SCAN-L calculations of dc conductivity of warm dense aluminum also give better agreement with experiments over other XC functionals such as PBE and SCAN-L.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Single-crystalline orthorhombic GdAlGe as a rare-earth magnetic Dirac nodal-line metal

Crystal engineering is a method for discovering new quantum materials and phases, which may be achieved using external pressure or strain. Chemical pressure is unique in that it generates internal pressure perpetually to the lattice. As an example, GdAlSi from the rare-earth (𝑅) 𝑅⁢Al⁢𝑋 (𝑋=Si or Ge) family of Weyl semimetals is considered. Replacing Si with the larger isovalent element Ge creates sufficiently large chemical pressure to induce a structural transition from the tetragonal structure of GdAlSi, compatible with a Weyl semimetallic state, to an orthorhombic phase in GdAlGe, resulting in an inversion-symmetry-protected nodal-line metal. We find that GdAlGe hosts an antiferromagnetic ground state with two successive orderings, at 𝑇 N ⁢1=35K and 𝑇 N ⁢2=30K. In-plane isothermal magnetization shows a magnetic field induced metamagnetic transition at 6.2 T for 2 K. Furthermore, electron-hole compensation gives rise to a large magnetoresistance of ∼100% at 2 K and 14 T. Angle-resolved photoemission spectroscopy measurements and density functional theory calculations reveal a Dirac-like linear band dispersion over an exceptionally large energy range of ∼1.5eV with a high Fermi velocity of ∼10 6 m/s, a rare feature not observed in any magnetic topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Computational Investigation of a CO 2 Conversion Strategy via Diels–Alder Reaction in a Carbon Capture Solvent

Molecular-level insights into reactive separations are crucial for the design of new conversion pathways of carbon dioxide (CO 2 ). This work explores a postulated pathway that directs CO 2 to undergo inverse-electron-demand Diels–Alder reactions to produce heterocycles using the CO 2 chemically fixed on water-lean solvent molecules. Density functional theory calculations are applied to evaluate the lowest unoccupied molecular orbital (LUMO) energies of three types of reactants (1,3-butadiene, 1,3-cyclohexadiene, and 1,2,4,5-tetrazine) with various functional substituents. These calculations also provide a data set (5.8k data) for developing a machine learning model to efficiently predict LUMO energies. A computational screening of LUMO energies for an additional 47k diene and tetrazine candidates is performed, and a list of candidates with lowered LUMO energies by electron-withdrawing substituents is provided. These candidates are further examined by their reaction energy barriers computed from the interatomic potential or density functional theory. Two major energy barriers are identified, one for the proton transfer within the water-lean solvent and the other for the CO 2 transfer from the solvent molecule to the reactant candidate (diene or tetrazine). The functional substituents have a more significant impact on the second barrier but a very slight one on the first barrier. This exploratory work demonstrates a new possibility for guiding experimental efforts toward the chemical conversion of fixated CO 2 to value-added compounds.

Chemical reactions↗
DOE OSTI2022

High-pressure structural systematics of dysprosium metal compressed in a neon pressure medium to 182 GPa

Here we present an experimental and theoretical study of dysprosium metal compressed in the soft pressure transmitting medium Ne up to 182 and 300 GPa, respectively. Angle-dispersive x-ray powder diffraction data from each of the high-pressure polymorphs shows anisotropic compression behavior indicating changes to the electron density distribution throughout its polymorphic landscape. We compare the monoclinic (mC4) and orthorhombic (oF16) structures for the collapsed structure for Dy above 82 GPa and verify that the oF16 structure offers a better fit to our data than the previously reported mC4 structure. Further, we have found that the oF16 structure undergoes similar anisotropic compression of its lattice parameters, with a turning point above 160 GPa; suggesting a potential phase transition at pressures much higher than achieved in this study. Density functional theory calculations show the likely candidate for this new high-pressure phase is the isosymmetric oF8 structure, which is predicted to be lower in energy than the oF16 structure above 275 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Molecular Fe-N 4 Moieties Coupled with Atomic Co-N 4 Sites Toward Improved Oxygen Reduction Performance

Research on high-efficiency and cost-efficient catalysts for oxygen reduction reaction (ORR) is still a vital but challenging issue for commercializing metal–air batteries. Herein, a single-molecule/atom hybrid catalyst is developed to boost the ORR, in which iron phthalocyanine molecules containing molecular Fe-N 4 moieties couple with atomic Co-N 4 sites on the surface of polyhedral carbon. Density functional theory calculations reveal that face-to-face laminated construction of Fe-N 4 and Co-N 4 in the hybrid catalyst can effectively modulate the electronic structure of active iron atoms and reduce the energy barrier of the rate-determining step for ORR. As a result, this hybrid catalyst demonstrates excellent ORR performance, featuring a half-wave potential of 0.904 V, a peak power density of 238.3 mW cm -2 for zinc–air battery, and outstanding electrocatalytic stability. Here, this work offers a distinctive and robust molecular/atomic engineering approach to creating efficient electrocatalysts, advancing the fields of metal–air batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spin qubit properties of the boron-vacancy/carbon defect in the two-dimensional hexagonal boron nitride

Spin qubit defects in two-dimensional materials have a number of advantages over those in three-dimensional hosts including simpler technologies for defect creation and control, as well as qubit accessibility. In this work, we select the V B C B defect in the hexagonal boron nitride (hBN) as a possible optically controllable spin qubit and explain its triplet ground state and neutrality. In this defect a boron vacancy is combined with a carbon dopant substituting the closest boron atom to the vacancy. Our density-functional-theory calculations confirmed that the system has dynamically stable spin triplet and singlet ground states. As revealed from our linear response GW calculations, the spin-sensitive electronic states are localized around the three undercoordinated N atoms and make local peaks in the density of electronic states within the bandgap. Using the triplet and singlet ground state energies, as well as the energies of the optically excited states, obtained from solution to the Bethe–Salpeter equation, we construct the spin-polarization cycle, which is found to be favorable for the spin qubit initialization. The calculated zero-field splitting parameters ensure that the splitting energy between the spin projections in the triplet ground state is comparable to that of the known spin qubits. We thus propose the V B C B defect in hBN as a promising spin qubit.

2D BN↗
DOE OSTI2023

Revisiting the structures and phase transitions of Ba 2 NaNb 5 O 15

The room-temperature and low-temperature structure(s) of Ba 2 NaNb 5 O 15 (BNN) have been debated since the structure was proposed in the 1960s. This work revisits the structures and phase transitions of BNN, combining high-resolution X-ray and neutron powder diffraction with density functional theory calculations. Temperature-dependent high-resolution X-ray powder diffraction patterns are collected from 4 to 918 K, and sequential batch Rietveld refinement using a symmetry mode approach to describe the structure is used to extract the main structural changes as a function of temperature. The data show that the average structure of BNN is best described by the Ama2 space group, and no other structural phase transitions were observed below the ferroelastic transition. The symmetry mode analysis, combining results from diffraction and density functional theory, shows significant octahedral tilting and corrugations of both the A1 and A2 sites along the c direction. A strong correlation between the spontaneous strain and the octahedral tilting was observed, and a potential connection with emerging microstructure at low temperatures is proposed, all enabled by the symmetry mode approach used in this work.

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

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2026

Hollow-structured Ni-N-C catalysts for highly selective CO 2 electroreduction

Atomically dispersed single-atom catalysts have emerged as promising non-precious catalyst alternatives to expensive Ag and Au catalysts for electrochemical CO 2 reduction reaction (CO 2 RR). In particular, nickel-nitrogen-carbon (Ni-N-C) catalysts have demonstrated a high faradaic efficiency (FE) toward CO formation at low overpotentials. Nonetheless, the exact nature of Ni active sites under CO 2 RR remains elusive and conventional Ni-N-C catalysts are limited by microporosity and low density of Ni single atoms, hindering performance in CO 2 electrolyzers. Here, we report the synthesis of hollow-structured Ni-N-C ( hs -Ni-N-C) catalysts via a post-synthesis modification (PSM) strategy using partial ligand exchange of 2-methylimidazole with 3-amino-1,2,4-triazole. This approach enables the formation of a hollow structure, resulting in more than a twofold increase in Ni atom density compared to regular Ni-N-C (r-Ni-N-C). In a zero-gap CO 2 electrolyzer, the optimized hs -Ni-N-C allows for achieving an FE CO of 97% at a current density of > 100 mA cm⁻ 2 , while maintaining high CO selectivity with stable performance over 100 h at 2.5 V. hs-Ni-N-C shows a more than sevenfold increase in the CO partial current density relative to r-Ni-N-C resulting from the combined effects of a higher density of Ni single-atom sites, improved kinetics, and lower transport resistance under the operating conditions, as indicated by electrochemical impedance spectra and distribution of relaxation times analysis. Operando high energy-resolution X-ray absorption spectroscopy (XAS) reveals that atop-bonded CO on Ni single sites induces dynamic transformations of the Ni–N coordination environment, leading to a symmetric coordination structure of hs -Ni-N-C. Under CO 2 RR, the catalysts undergo a more pronounced structural change and form a minor fraction of Ni nanoparticles. Density functional theory calculations are consistent with the XAS results and provide molecular insights showing that the interplay between protonation and CO adsorption leads to adsorbate-induced restructuring of the Ni single atom. This work demonstrates the synergistic role of hollow structure and high-density Ni atoms in governing CO 2 RR selectivity and provides mechanistic insights into the structural dynamics of single-atom catalysts under operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An improved oxygen reduction reaction activity and CO 2 -tolerance of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ achieved by a surface modification with barium cobaltite coatings

Solid oxide fuel cells (SOFCs) cathode often suffers from the poisoning effect of the contaminants commonly encountered in air such as CO 2 . Here we report an effective approach to enhancing the activity and CO 2 tolerance of the state-of-the-art La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF) cathode enabled by a coating of BaCoO 3-δ (BCO), as verified by the electrochemical testings, Raman analyses, and density functional theory calculations. When surface modified with a thin-film BCO coating, LSCF displays a much enhanced ORR activity and an improved durability against CO 2 . For example, anode supported SOFCs with the LSCF cathode coated with BCO coatings show a remarkable peak power density (Pmax) of 0.41 Wcm -2 and a significantly reduced degradation rate in current density of ~0.08% h -1 at 0.8 V and 700 °C for a period of 300 hs when humidified H2 (with 3 vol%H2O) was used as fuel and air with 8 vol% CO 2 as oxidant. The demonstrated performance is improved when compared with those of the cells with a blank LSCF electrode (a Pmax of ~0.36 Wcm -2 and a degradation rate of ~0.15% h -1 ) under the same conditions. Furthermore, the adsorption energy calculations suggests that BCO coating makes CO 2 adsorption much weaker than LSCF (-0.54 eV versus -1.07 eV).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Effects of thin metal contacts on few-layer van der Waals ferrielectric CuInP 2 S 6

Out-of-plane polarized ferroelectric materials in a capacitive structure provide a key component for several technological applications. Furthermore, two-dimensional materials are expected to aid in the quest for both ultrathin and flexible electronics. Of the various two-dimensional ferroelectrics with out-of-plane polarization, CuInP 2 S 6 is special in that the Cu atoms are highly mobile and it has been shown to possess both low- and high-polarization states. Here, using density-functional-theory calculations, we explore the stabilization of the ferroelectric state for several prototypical metal contacts (Gr, Ni, Cu, Au, and Ag). In all cases, we find that the ferroelectric state can be stabilized at fewer layers than in the freestanding case. For all of the considered conventional metal contacts, we also find the existence of a quasi-ferroelectric state that stabilizes a polar phase for thicknesses greater than two layers of CIPS. In the cases of Au and Ag, interfacial alignment and strain can be used to stabilize ferroelectricity at the bilayer limit. Furthermore, we find that the strength of the interaction between the contact and CuInP 2 S 6 also leads to stabilization of the high-polarization state when ferroelectricity is stabilized. Lastly, energy-barrier calculations show that the system is still switchable in the presence of contact doping from the metal contacts.

2D materials↗
DOE OSTI2024

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
DOE OSTI2023

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Importance of bulk excitations and coherent electron-photon-phonon scattering in photoemission from PbTe(111): Ab initio theory with experimental comparisons

This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Protein C-GeM: A Coarse-Grained Electron Model for Fast and Accurate Protein Electrostatics Prediction

The electrostatic potential (ESP) is a powerful property for understanding and predicting electrostatic charge distributions that drive interactions between molecules. In this study, we compare various charge partitioning schemes including fitted charges, density-based quantum mechanical (QM) partitioning schemes, charge equilibration methods, and our recently introduced coarse-grained electron model, C-GeM, to describe the ESP for protein systems. Furthermore, when benchmarked against high quality density functional theory calculations of the ESP for tripeptides and the crambin protein, we find that the C-GeM model is of comparable accuracy to ab initio charge partitioning methods, but with orders of magnitude improvement in computational efficiency since it does not require either the electron density or the electrostatic potential as input.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗