Search NASA⌕ Search

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.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of External Conditions on the Black-to-Yellow Phase Transition of CsPbI 3 Based on First-Principles Calculations: Pressure and Moisture

All-inorganic CsPbI 3 perovskite has emerged as a promising candidate for next-generation solar cells. However, owing to the poor structural stability of black phase CsPbI 3 perovskite at room temperature, it spontaneously transforms to the yellow, photoinactive, nonperovskite phase at room temperature, limiting further development of CsPbI 3 perovskite solar cells. Here, to better understand the mechanism driving such undesirable phase transformation, we examine the thermodynamics and kinetics associated with γ-to-δ phase transformation using density functional theory calculation. A solid-state nudged elastic band method was employed to find minimum energy pathways assuming a concerted solid-state phase transformation between these two phases. The gas molecules representing atmospheric (H 2 O, O 2 , and N 2 ) and inert (Ar) ambient environments, as well as applied external pressure, were considered to examine the influence of external conditions on the black-to-yellow phase transition. Our calculation reveals that, with increasing pressure, the thermodynamic driving force for converting the γ-phase to the δ-phase increases, whereas a larger activation energy must be overcome for the phase transition to occur. We also investigate the moisture-induced phase transformation with an H 2 O molecule and its dissociated species (H + /OH – ) and demonstrate that the reaction energy barriers can be significantly lowered in the presence of H 2 O or OH – . On the other hand, other nonpolar species (O 2 , N 2 , and Ar) have a negligible effect on the phase transformation kinetics, while they might reduce the thermodynamic driving force for the phase transformation and suppress the undesirable phase transformation. Our theoretical prediction could support recent observations that the γ-to-δ phase transformation occurs rapidly and catalytically in the presence of moisture, whereas in dry argon and dry oxygen atmospheres, the γ-CsPbI 3 remains stable.

14 SOLAR ENERGY↗
DOE OSTI2024

Thermally induced structural evolution and nanoscale interfacial dynamics in Bi-Sb-Te layered nanostructures

Layered chalcogenides, including Bi-Sb-Te ternary alloys and heterostructures, are renowned as thermoelectric and topological insulators and have recently been highlighted as plasmonic building blocks beyond noble metals. Here, we conduct joint in situ transmission electron microscopy and density functional theory calculations to investigate the temperature-dependent nanoscale dynamics and interfacial properties, identifying the role of native defects and edge configurations in the anisotropic sublimation of Bi 2 Te 3 -Sb 2 Te 3 heterostructures and Sb 2-x Bi x Te 3 alloys. We report structural dynamics, including edge evolution, layer-by-layer sublimation, and the formation and coalescence of thermally induced polygonal nanopores. These nanopores are initiated by preferential dissociation of tellurium, reducing thermal stability in heterostructures. Triangular and quasi-hexagonal configurations dominate nanopore structures in heterostructures. Our calculations reveal antisite defects (Te Sb and Te Bi ) as key players in defect-assisted sublimation. These findings enhance our understanding of nanoscale dynamics and assist in designing tunable low-dimensional chalcogenides.

Bi2Te3-Sb2Te3 heterostructure↗
DOE OSTI2023

Machine learning the relationship between Debye temperature and superconducting transition temperature

Recently a relationship between the Debye temperature $Θ_D$ and the superconducting transition temperature $T_c$ of conventional superconductors has been proposed [Esterlis et al., npj Quantum Mater. 3, 59 (2018)]. The relationship indicates that $T_c$ ≤ $AΘ_D$ for phonon-mediated BCS superconductors, with $A$ being a prefactor of order ~ $0.1$. In order to verify this bound, we train machine learning (ML) models with 10 330 samples in the Materials Project database to predict $Θ_D$. Here, by applying our ML models to 9860 known superconductors in the NIMS SuperCon database, we find that the conventional superconductors in the database indeed follow the proposed bound. We also perform first-principles phonon calculations for $\mathrm{H_3S}$ and $\mathrm{LaH_{10}}$ at 200 GPa. The calculation results indicate that these high-pressure hydrides essentially saturate the bound of $T_c$ versus $Θ_D$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Unusual Exchange Couplings and Intermediate Temperature Weyl State in Co 3 Sn 2 S 2

Understanding magnetism and its possible correlations to topological properties has emerged to the forefront as a difficult topic in studying magnetic Weyl semimetals. Co 3 Sn 2 S 2 is a newly discovered magnetic Weyl semimetal with a kagome lattice of cobalt ions and has triggered intense interest for rich fantastic phenomena. In this study, we report the magnetic exchange couplings of Co 3 Sn 2 S 2 using inelastic neutron scattering and two density functional theory (DFT) based methods: constrained magnetism and multiple-scattering Green’s function methods. Co 3 Sn 2 S 2 exhibits highly anisotropic magnon dispersions and linewidths below T C , and paramagnetic excitations above T C . The spin-wave spectra in the ferromagnetic ground state is well described by the dominant third-neighbor “across-hexagon” J d model. Our density functional theory calculations reveal that both the symmetry-allowed 120° antiferromagnetic orders support Weyl points in the intermediate temperature region, with distinct numbers and the locations of Weyl points. Our study highlights the important role Co 3 Sn 2 S 2 can play in advancing our understanding of kagome physics and exploring the interplay between magnetism and band topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Observation and Characterization of Vibrationally Active Surface Species Accessed with Nonthermal Nitrogen Plasmas

Polycrystalline Ni, Pd, Cu, Ag, and Au foils exposed to nonthermal plasma (NTP)-activated N 2 are found to exhibit a vibrational feature near 2200 cm –1 in polarization-modulation infrared reflection-absorption spectroscopy (PM-IRAS) observations that are not present in the same materials exposed to N 2 under nonplasma conditions. The feature is similar to that reported elsewhere and is typically assigned to chemisorbed N 2 . We employ a combination of temperature-dependent experiments, sequential dosing, X-ray photoelectron spectroscopy, isotopic labeling, and density functional theory calculations to characterize the feature. Results are most consistent with a triatomic species, likely NCO, with the C and O likely originating from ppm-level impurities in the ultrahigh-purity (UHP) Ar and/or N 2 gas cylinders. Here, this work highlights the potential for nonthermal plasmas to access adsorbates inaccessible thermally as well as the potential contributions of ppm-level impurities to corrupt the interpretation of plasma catalytic chemistry.

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

Factors influencing quantum evaporation of helium from polar semiconductors from first principles

While there is much indirect evidence for the existence of dark matter (DM), to date it has evaded detection. Current efforts focus on DM masses over ~GeV—to push the sensitivity of DM searches to lower masses, new DM targets and detection schemes are needed. In this work, we focus on the latter—a novel detection scheme recently proposed to detect 10–100 meV phonons in polar target materials. Previous work showed that well-motivated models of DM can interact with polar semiconductors to produce an athermal population of phonons. This new sensing scheme proposes that these phonons then facilitate quantum evaporation of 3 He from a van der Waals film deposited on the target material. However, a fundamental understanding of the underlying process is still unclear, with several uncertainties related to the precise rate of evaporation and how it can be controlled. In this work, we use ab initio density functional theory calculations to compare the adsorption energies of helium atoms on a polar target material, sodium iodide, to understand the underlying evaporation physics. We explore the role of surface termination, monolayer coverage, and elemental species on the rate of He evaporation from the target material. Using this, we discuss the optimal target features for He-evaporation experiments and their range of tunability through chemical and physical modifications such as applied field and surface termination.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The role of surfaces in flexoelectricity

We develop the relationship between the strain derivative of the mean-inner potential and surface contributions to flexoelectricity, identifying the true surface-specific component of the flexoelectric response of finite samples. Density functional theory calculations on a range of experimentally observed, low energy SrTiO 3 , MgO, and Si surfaces demonstrate that the mean-inner potential and its contributions to flexoelectricity are sensitive to small differences in surface structure, chemistry, and adsorbates. We also introduce a method to estimate mean-inner potential contributions to flexoelectricity using electron scattering factors and use this approximation to predict total flexoelectric responses for a variety of insulators. Strategies to experimentally disentangle bulk and surface flexoelectric terms are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Frustration and Atomic Ordering in a Monolayer Semiconductor Alloy

Frustrated interactions can lead to short-range ordering arising from incompatible interactions of fundamental physical quantities with the underlying lattice. The simplest example is the triangular lattice of spins with antiferromagnetic interactions, where the nearest-neighbor spin-spin interactions cannot simultaneously be energy minimized. In this work, we show that engineering frustrated interactions is a possible route for controlling structural and electronic phenomena in semiconductor alloys. Using aberration-corrected scanning transmission electron microscopy in conjunction with density functional theory calculations, we demonstrate atomic ordering in a two-dimensional semiconductor alloy as a result of the competition between geometrical constraints and nearest-neighbor interactions. Statistical analyses uncover the presence of short-range ordering in the lattice. In addition, we show how the induced ordering can be used as another degree of freedom to considerably modify the band gap of monolayer semiconductor alloys.

2-dimensional systems↗
DOE OSTI2025

Accelerated Discovery of Cost-Effective Photoabsorber Materials for Near-Infrared (λ = 1600 nm) Photodetector Applications

Current infrared sensing devices are based on costly materials with relatively few viable alternatives known. To identify promising candidate materials for infrared photodetection, we have developed a high-throughput screening methodology based on high-accuracy r 2 SCAN and HSE calculations in density functional theory. Using this method, we identify ten already synthesized materials between the inverse perovskite family, the barium silver pnictide family, the alkaline pnictide family, and ZnSnAs 2 as top candidates. Among these, ZnSnAs 2 emerges as the most promising candidate due to its experimentally verified band gap of 0.74 eV at 0 K and its cost-effective synthesis through Bridgman growth. BaAgP also shows potential with an HSE-calculated band gap of 0.64 eV, although further experimental validation is required. Lastly, we discover an additional material, Ca 3 BiP, which has not been previously synthesized, but exhibits a promising optical spectra and a band gap of 0.56 eV. The method applied in this work is sufficiently general to screen wider bandgap materials in high-throughput and now extended to narrow-band gap materials.

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

The Bonding Nature and Adhesion of Polyacrylic Acid Coating on Li-Metal for Li Dendrite Prevention

The success of polyacrylic acid (PAA) to suppress Li dendrite growth suggests that the mechanical properties of polymer-based coatings, including the modulus, toughness, and interfacial adhesion are important design criteria. However, the measurement of the adhesion of thin PAA, as well as other polymer coatings to the reactive Li-metal anode surface is limited experimentally and challenging computationally. In this paper, a strategy was proposed to estimate the adhesion and delamination of the PAA(polymer)/Li interface, based on the bonding nature at the simpler PAA (oligomer)/Li interfaces using density functional theory calculations. It has been shown that the carboxylic acid groups in PAA reacted strongly with metallic Li, which significantly enhances the interfacial adhesion through the Li–O bonds formation, Li ionization and its incorporation into PAA, and –H or –OH termination of Li after decomposition of the COOH functional group. During delamination, it was found that the most likely PAA delamination route involved breaking partial Li–O bonds and lifting some ionized Li atoms from the Li-metal, especially for the Li atoms that showed a charge closer to +1 or are bonded with two O atoms from PAA. Based on the average bonding energies from PAA(oligomer)/Li interface delamination calculations, the work of separation, W sep , of the PAA(polymer)/Li interface was estimated to be ~1.0 (J/m 2 ). The high W sep of PAA (polymer)/Li was comparable with the Li 2 O/Li interface and higher than Li 2 CO 3 /Li and LiF/Li interfaces. This order correlated well with the areal density of Li–O bonds, which can serve as a descriptor for the interfacial adhesion. Furthermore, this computational approach can be applied to other interfaces with polymer-based coatings.

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

4 d element induced improvement of structural disorder and development of weakly reentrant spin-glass behavior in NiRuMnSn

The pursuit of efficient spin polarization in quaternary Heusler alloys with the general formula XX'YZ (where X, X', and Y are transition metals and Z is a p-block element) has been a subject of significant scientific interest. While previous studies showed that isoelectronic substitution of a 4d element in place of a 3d element in quaternary Heusler alloys improves the half-metallic ferromagnetic characteristics, our research here, on the quaternary Heusler alloy NiRuMnSn, reveals a strikingly different scenario. In this study, we present a detailed structural analysis of the material using x-ray absorption fine structure and neutron diffraction techniques which confirms the formation of a single-phase compound with 50:50 site disorder between Ni and Ru atoms at 4c and 4d sites. Contrary to expectations, our density functional theory calculations suggest a considerable decrease in spin polarization even in the ordered structure. Additionally, we report on the compound's exceptional behavior, displaying a rare reentrant spin-glass property below ~ 60 K, a unique and intriguing feature for quaternary Heusler-type compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic and optical properties of Zn-doped β-Ga 2 O 3 Czochralski single crystals

β-Ga 2 O 3 has several soluble deep acceptors that impart insulating behavior. In this study, we investigate Zn doping (0.25at.%) in bulk Czochralski and vertical gradient freeze β-Ga 2 O 3 . Representative crystals were assessed for orientation (electron backscatter diffraction and Raman spectroscopy), purity (glow discharge mass spectrometry and secondary ion mass spectrometry), optical properties (ultraviolet to near infrared absorption), and electrical properties (resistivity and current–voltage). Purity measurements indicate that Zn evaporation is insufficient to inhibit doping of Zn into β-Ga 2 O 3 . Hybrid functional calculations show Zn substitutes nearly equally on tetrahedral and octahedral sites, with less than ~0.1eV preference for the octahedral (Ga II ) site. Furthermore, calculations show that Zn Ga acts as a deep acceptor with trapping levels ~1.3 and ~0.9eV above the valence band for one and two holes, respectively. The solubility and electronic behavior of Zn dopants are consistent with measured concentrations >1×10 18 atoms/cm 3 and electrical measurements that show resistivity 10 11 –10 13 Ωcm, with no p-type conduction.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

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

Competing energetic states in γ – Fe 2 W O 6 with strong spin-charge-lattice coupling

We report magnetic and electronic properties of γ – Fe 2 W O 6 via neutron powder-diffraction measurements and first-principles density-function theory calculations. We reveal a magnetic ground state with two phases coexisting in which the minor phase is the same as an earlier report, despite the fact that both materials studied have the same space group but with slightly different atomic positions. Interestingly, both magnetic phases are well captured by first-principles calculations. Furthermore, we show that these two magnetic phases are correlated with electronic properties, with one being insulating and the other being metallic. These features suggest that γ – Fe 2 W O 6 exhibits competing energetic states in which spin, charge, and lattice degrees of freedom are strongly coupled to each other.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Understanding copper diffusion in CuInSe 2 with first-principles based atomistic and continuum models

We investigate the diffusion of copper in CuInSe 2 using thermodynamic and kinetic models based on density functional theory calculations, attempting to reconcile large differences in reported experimental diffusivities. We find that observations of rapid chemical diffusion can be explained by large thermodynamic factors, which we calculate using a compositionally constrained model of intrinsic point defect formation. We further characterize how copper diffusion coefficients depend on material synthesis conditions and exhibit their variation across the CuInSe 2 secondary phase diagram. In doing so, we identify stable off-stoichiometries that are dominated by either vacancy- or interstitial-mediated diffusion mechanisms. These results are employed in the development of a continuum reaction–diffusion model, which we use to simulate experimental depth profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Facile ab initio approach for self-localized polarons from canonical transformations

Electronic states in a crystal can localize due to strong electron-phonon (e-ph) interactions, forming so-called small polarons. Methods to predict the formation and energetics of small polarons are either computationally costly or not geared toward quantitative predictions. Here we show a formalism based on canonical transformations to compute the polaron formation energy and wave function using ab initio e-ph interactions. Comparison of the calculated polaron and band-edge energies allows us to determine whether charge carriers in a material favor a localized small polaron over a delocalized Bloch state. Due to its low computational cost, our approach enables efficient studies of the formation and energetics of small polarons, as we demonstrate by investigating electron and hole polaron formation in alkali halides and metal oxides and peroxides. We outline refinements of our scheme and extensions to compute transport in the polaron hopping regime.

36 MATERIALS SCIENCE↗