Search NASA⌕ Search

SEARCH · Search NASA

Results for “First Principles 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-Principles Calculations of Thermoelectric Transport Properties of Quaternary and Ternary Bulk Chalcogenide Crystals

Chalcogenide crystals have a wide range of applications, especially as thermoelectric materials for energy conversion. Thermoelectric materials can be used to generate an electric current from a temperature gradient based on the Seebeck effect and based on the Peltier effect, and they can be used in cooling applications. Using first-principles calculations and semiclassical Boltzmann theory, we have computed the Seebeck coefficient, electrical conductivity, electronic thermal conductivity, power factor, and figure of merit of 30 chalcogenide crystals. A Quantum Espresso package is used to calculate the electronic properties and locate the Fermi level. The transport properties are then calculated using the BoltzTraP code. The 30 crystals are divided into two groups. The first group has four crystals with quaternary composition (A 2 BCQ 4 ) (A = Tl; B = Cd, Hg; C = Si, Ge, Sn; Q = S, Se, Te). The second group contains 26 crystals with the ternary composition (A’B’Q 2 ) (A’ = Ag, Cu, Au, Na; B’ = B, Al, Ga, In; Q = S, Se, Te). Among these 30 chalcogenide crystals, the results for 11 crystals: Tl 2 CdGeSe 4 , Tl 2 CdSnSe 4 , Tl 2 HgSiSe 4 , Tl 2 HgSnS 4 , AuBSe 2 , AuBTe 2 , AuAlTe 2 , AuGaTe 2 , AuInTe 2 , AgAlSe 2 , and AgAlTe 2 are revealed for the first time. In addition, temperature-dependent transport properties of pure and doped AgSbSe 2 and AgSbTe 2 crystals with dopant compositions of AgSb 0.94 Cd 0.06 Te 2 and AgSbTe 1.85 Se 0.15 were explored. These results provide an excellent database for bulk chalcogenides crucial for a wide range of potential applications in renewable energy fields.

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

Tunable energy landscape of screw dislocation cores by compositional fluctuations in bcc high-entropy alloys from first-principles calculations

The energy landscape of screw dislocation cores plays a central role in dislocation-mediated deformation mechanisms in body-centered cubic (bcc) metals. In bcc high-entropy alloys (HEAs), this energy landscape is modulated by local compositional fluctuations, which has important implications for deformation processes in these materials. Through first-principles calculations, this study investigates high-symmetry screw dislocation core structures in NbTaMoW and NbTaTiHf bcc HEAs. The results show that alloying group IV transition metals lead to large local lattice distortions at dislocation cores, which is demonstrated to be an important factor governing fluctuations in core configurations along a dislocation line. Importantly, group IV elements near the core induce features in the energy landscape that are exclusive for HEAs, specifically lowering the energy of core configurations that are unstable in elemental bcc metals. A combined influence of these chemical effects with crystallographic details enables the activation of glide planes, a feature that has been linked to ductility improvements in bcc HEAs. These findings provide new insights into the atomic-scale mechanisms underlying dislocation mobility in bcc HEAs, offering a pathway for designing materials with tailored mechanical properties.

Borges, Pedro P P O↗
DOE OSTI2025

Chemical trends favoring interstitial cluster formation in bcc high-entropy alloys from first-principles calculations

Achieving high strength and ductility is a common goal in structural alloy design. Body-centered cubic high-entropy alloys (HEAs) commonly highlight the conflict between these properties, with stronger alloys being brittle and vice versa. Recent reports suggest interstitial solutes can be used to overcome this trade-off, in some cases providing both strength and ductility enhancements. This effect has been correlated with interstitial cluster formation, although the conditions favoring their formation remain incompletely understood. Using first-principles calculations of solution energies and diffusivities, we provide insights into thermodynamic and kinetic factors favoring interstitial solute clusters. Among C, N and O solutes, O interstitials display most desirable diffusion kinetics. Further, the results highlight the importance of local composition fluctuations in the HEAs to enable the formation of clusters of appreciable size. The results are explained in terms of bonding and distortion trends across solutes and HEA compositions to provide guidelines for alloy design.

Borges, Pedro P P O↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2024

First-principles calculation of lattice distortion, electronic structure, and bonding properties of GeTe-based and PbSe-based high-entropy chalcogenides

The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.

Physics↗
DOE OSTI2022

A thermochemical database from high-throughput first-principles calculations and its application to analyzing phase evolution in AM-fabricated IN718

A comprehensive thermochemical database is constructed based on high–throughput first-principles phonon calculations of over 3000 atomic structures in limited concentrations in Ni, Fe, and Co alloys involving a total of 26 elements including Al, B, C, Cr, Cu, Hf, La, Mn, Mo, N, Nb, O, P, Re, Ru, S, Si, Ta, Ti, V, W, Y, and Zr, providing thermochemical data largely unavailable from existing experiments. Here, the database can be employed to predict the equilibrium phase compositions and fractions directly from first-principles by minimizing the chemical potential of a multicomponent system with a fixed overall chemical composition and a fixed temperature. It is applied to the additively manufactured nickel-based IN718 superalloy to analyze the phase evolution with temperature. IN718 is known for its great performance in tensile, fatigue, creep, and rupture strength, combined with easy fabrication and corrosion resistance. In particular, we successfully predicted the formation of L1 0 -FeNi, γ’-Ni 3 (Fe,Al), α-Cr, δ-Ni 3 (Nb,Mo), γ”-Ni 3 Nb, and η-Ni 3 Ti at low temperatures (below 680 K), γ’-Ni 3 Al, δ-Ni 3 Nb, γ”-Ni 3 Nb, α-Cr, and γ-Ni(Fe,Cr,Mo) at intermediate temperatures (between 680 and 1140 K), and δ-Ni 3 Nb and γ-Ni(Fe,Cr,Mo) at high temperatures (above 1140 K) in IN718. These predictions are validated by EDS mapping of compositional distributions and corresponding identifications of phase distributions. The database is expected to be a valuable source for future thermodynamic analysis and microstructure prediction of alloys involving the 26 elements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗
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 OSTI2021

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗
DOE OSTI2023

Computational Materials Design for Ceramic Nuclear Waste Forms Using Machine Learning, First-Principles Calculations, and Kinetics Rate Theory

Ceramic waste forms are designed to immobilize radionuclides for permanent disposal in geological repositories. One of the principal criteria for the effective incorporation of waste elements is their compatibility with the host material. In terms of performance under environmental conditions, the resistance of the waste forms to degradation over long periods of time is a critical concern when they are exposed to natural environments. Due to their unique crystallographic features and behavior in nature environment as exemplified by their natural analogues, ceramic waste forms are capable of incorporating problematic nuclear waste elements while showing promising chemical durability in aqueous environments. Recent studies of apatite- and hollandite-structured waste forms demonstrated an approach that can predict the compositions of ceramic waste forms and their long-term dissolution rate by a combination of computational techniques including machine learning, first-principles thermodynamics calculations, and modeling using kinetic rate equations based on critical laboratory experiments. By integrating the predictions of elemental incorporation and degradation kinetics in a holistic framework, the approach could be promising for the design of advanced ceramic waste forms with optimized incorporation capacity and environmental degradation performance. Such an approach could provide a path for accelerated ceramic waste form development and performance prediction for problematic nuclear waste elements.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

First-principles calculations of metal surfaces. I. Slab-consistent bulk reference for convergent surface properties

The first-principles computation of the surfaces of metals is typically accomplished through slab calculations of finite thickness. The extraction of a convergent surface formation energy from slab calculations is dependent upon defining an appropriate bulk reference energy. I describe a method for an independently computed, slab-consistent bulk reference that leads to convergent surface formation energies from slab calculations that also provides realistic uncertainties for the magnitude of unavoidable nonlinear divergence in the surface formation energy with slab thickness. In conclusion, the accuracy is demonstrated on relaxed, unreconstructed low-index aluminum surfaces with slabs with up to 35 layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First-Principles Calculation of Electroweak Box Diagrams from Lattice QCD

We present the first realistic lattice QCD calculation of the $γ W$-box diagrams relevant for beta decays. The nonperturbative low-momentum integral of the $γ W$ loop is calculated using a lattice QCD simulation, complemented by the perturbative QCD result at high momenta. Using the pion semileptonic decay as an example, we demonstrate the feasibility of the method. By using domain wall fermions at the physical pion mass with multiple lattice spacings and volumes, we obtain the axial γ W -box correction to the semileptonic pion decay, $\square$$\frac{VA}{γW}|_π$ = 2.830 (11) stat (26 )syst × 10 –3 , with the total uncertainty controlled at the level of ~1% . This study sheds light on the first-principles computation of the $γ W$-box correction to the neutron decay, which plays a decisive role in the determination of $|V_ud|$ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Database and deep-learning scalability of anharmonic phonon properties by automated brute-force first-principles calculations

Understanding the anharmonic phonon properties of crystal compounds—such as phonon lifetimes and thermal conductivities—is essential for investigating and optimizing their thermal transport behaviors. These properties also impact optical, electronic, and magnetic characteristics through interactions between phonons and other quasiparticles and fields. In this study, we develop an automated first-principles workflow to calculate anharmonic phonon properties and build a comprehensive database encompassing more than 6500 inorganic compounds. Utilizing this dataset, we train a graph neural network model to predict thermal conductivity values and spectra from structural parameters, demonstrating a scaling law in which prediction accuracy improves with increasing training data size. High-throughput screening with the model enables the identification of materials exhibiting extreme thermal conductivities—both high and low. The resulting database offers valuable insights into the anharmonic behavior of phonons, thereby accelerating the design and development of advanced functional materials.

Ohnishi, Masato [University of Tokyo (Japan); Inst↗
DOE OSTI2021

Dy adsorption on and intercalation under graphene on 6 H -SiC(0001) surface from first-principles calculations

Previous experimental observations motivate clarification of configuration stabilities and kinetic processes for intercalation of guest atoms into a layered van der Waals material such as a graphene-SiC system. From our first-principles density functional theory (DFT) calculations, we analyze Dy 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 intercalation of a single Dy atom into the gallery between TLG and BLG is more favorable than adsorption on TLG but that intercalation into the gallery underneath BLG is highly unfavorable. We obtain diffusion barriers of ~0.45 and 0.54 eV for a Dy atom diffusing on and under TLG, respectively. We find that the direct penetration of a Dy atom from the graphene top into the gallery under TLG is almost inhibited below a temperature of ~1400 K due to a large global barrier of at least ~3.5 eV. Instead, we find that a single Dy atom on TLG can easily intercalate by crossing a TLG step (e.g., a zigzag step presaturated by a Dy chain or a reconstructed zigzag step zz57). Additionally, we also perform DFT calculations for different Dy coverages to demonstrate how the favorability of Dy intercalation, as well as the corresponding interlayer spacings, depend on the coverage. Consequently, we can provide general insight and guidance for extensively studied systems involving intercalation of foreign atoms into graphene on a SiC substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First-principles calculation of spin and orbital contributions to magnetically ordered moments in Sr 2 IrO 4

Here, we show how an accurate first-principles treatment of the canted-anti-ferromagnetic ground state of Sr 2 IrO 4 , a prototypical 5d correlated spin-orbit coupled material, can be obtained without invoking any free parameters, such as the Hubbard $\textit{U}$ or tuning the spin-orbit coupling strength. Our theoretically predicted iridium magnetic moment of 0.250$μ_B$, canted by 12.6° off the $\textit{a}$ axis, is in accord with experimental results. By resolving the magnetic moments into their spin and orbital components, we show that our theoretically obtained variation of the magnetic scattering amplitude $\langle M_m \rangle$ as a function of the polarization angle is consistent with recent nonresonant magnetic x-ray scattering measurements. The computed value of the band gap (55 meV) is also in line with the corresponding experimental values. A comparison of the band structure to that of the cuprates suggests the presence of incommensurate charge-density wave phases in Sr 2 IrO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗