Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Heterostructural interface engineering for ultrawide-gap nitrides from first principles: Ta C / Al N and Ta C / Ga N rocksalt-wurtzite interfaces

Epitaxial lattice matching is an important condition for the formation of coherent interfaces with low defect densities. However, lattice-matched substrates with the same crystal structure as the active layer are often not available, suggesting opportunities for utilizing heterostructural interfaces. For example, at high Al contents that are interesting for ultrawide-gap applications in power electronics, Al x ⁢Ga 1-x ⁢N semiconductor alloys in the (0001) orientation of the wurtzite (wz) structure become lattice-matched to (111)-oriented rocksalt (rs) Ta⁢C substrates. To predict the expected interface atomic structures under different synthesis conditions, we perform high-throughput density-functional-theory calculations, using an algorithm for systematic sampling of the possible stacking sequences of the atomic layers on the in-plane hexagonal lattice. The approach considers octahedral, tetrahedral, and prismatic coordination motifs, and is generally applicable for the modeling of commensurate rs/wz heterostructural interfaces. Our results provide guidance for synthesis control of substrate-film bonding and the polarity of ultrawide-gap Al x⁢ Ga 1-x⁢ N alloys on Ta⁢C substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio study on fcc Pr with correlation matrix renormalization theory

Here, we studied fcc praseodymium (Pr) with the ab initio correlation matrix renormalization theory (CMRT) explicitly calculating the Coulomb interactions among basis orbitals and without using any adjustable parameters to work for strongly correlated electron systems. We calculated its total energy in a paramagnetic ground state and studied the role of the correlated 4f electrons in the system. Good agreement was obtained between CMRT and experiments in the pressure volume dependence of the fcc phase. We also compared the CMRT results against other theoretical methods including local density approximation + dynamical mean-field theory and showed consistent results among them. Moreover, we found the normalized local charge fluctuation of the 4f electrons as the function of volume exhibits a clear slope change at the volume collapse region, indicating a switch in their correlation nature there.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Charge-transfer satellites and chemical bonding in photoemission and x-ray absorption of SrTi O 3 and rutile Ti O 2 : Experiment and first-principles theory with general application to spectroscopic analysis

First-principles, real-time-cumulant, and Bethe-Salpeter-equation calculations fully capture the detailed satellite structure that occurs in response to the sudden creation of the core hole in both photoemission and x-ray absorption spectra of the transition-metal compounds SrTi O 3 and rutile Ti O 2 . Analysis of the excited-state, real-space charge-density fluctuations betrays the physical nature of these many electron excitations that are shown to reflect the materials’ solid-state electronic structure and chemical bonding. This first-principles development of the cumulant-based core hole spectral function is generally applicable to other systems and should become a standard tool for all similar spectroscopic analysis going beyond the quasiparticle physics of the photoelectric effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Evidence of Ba-substitution induced spin-canting in the magnetic Weyl semimetal EuCd 2 As 2

Recently EuCd 2 As 2 was predicted to be a magnetic Weyl semimetal with a lone pair of Weyl nodes generated by A-type antiferromagnetism and protected by a rotational symmetry. However, it was soon discovered that the actual magnetic structure broke the rotational symmetry and internal pressure was later suggested as a route to stabilize the desired magnetic state. In this work we test this prediction by synthesizing a series of Eu 1-x BaxCd 2 As 2 single crystals and studying their structural, magnetic, and transport properties via both experimental techniques and first-principles calculations. We find that small concentrations of Ba (~3%–10%) lead to a small out-of-plane canting of the Eu moment. However, for higher concentrations this effect is suppressed and a nearly in-plane model is recovered. Studying the transport properties we find that all compositions show evidence of an anomalous Hall effect dominated by the intrinsic mechanism as well as large negative magnetoresistances in the longitudinal channel. A nonmonotonic evolution of the transport properties is seen across the series which correlates to the proposed canting suggesting canting may enhance the topological effects. Careful density functional theory calculations using an all-electron approach revise prior predictions finding a purely ferromagnetic ground state with in-plane moments for both the EuCd 2 As 2 and Eu 0.5 Ba 0.5 Cd 2 As 2 compounds, corroborating our experimental findings. This work suggests that Ba substitution can tune the magnetic properties in unexpected ways which correlate to changes in measures of topological properties, encouraging future work to locate the ideal Ba concentration for Eu moment canting.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Efficient simulations of charge density waves in the transition metal Dichalcogenide TiSe 2

Charge density waves (CDWs) in transition metal dichalcogenides are the subject of growing scientific interest due to their rich interplay with exotic phases of matter and their potential technological applications. Here, using density functional theory with advanced meta-generalized gradient approximations (meta-GGAs) and linear response time-dependent density functional theory (TDDFT) with state-of-the-art exchange-correlation kernels, we investigate the electronic, vibrational, and optical properties in 1T-TiSe 2 with and without CDW. In both bulk and monolayer TiSe 2 , the electronic bands and phonon dispersions in either normal or CDW (semiconducting) phase are described well via meta-GGAs, which separate the valence and conduction bands just as HSE06 does but with significantly more computational feasibility. The experimentally observed humps of electron energy loss spectroscopy are successfully reproduced in TDDFT. Our work opens the door to simulating these complexities in CDW compounds from first principles by revealing meta-GGAs as an accurate low-cost alternative to HSE06.

Electronic properties and materials↗
DOE OSTI2026

Ferrimagnetic Order in Tetragonal Antiperovskite Mn3GeN

The crystal and magnetic structures of the nitride antiperovskite M⁢n3⁢GeN reveals ferrimagnetic order stemming from a distorted kagome-derived lattice of the Mn atoms. Polycrystalline M⁢n3⁢GeN was synthesized via a solid-state reaction and characterized using neutron powder diffraction, dc magnetometry, and first-principles calculations. Rietveld refinement reveals near-stoichiometric composition (M⁢n3⁢Ge⁢N0.992⁢(7)) adopting a tetragonal 𝐼⁢4/𝑚⁢𝑐⁢𝑚 structure at 𝑇=500K and below, featuring axially distorted and tilted [NM⁢n6] octahedra that result in a buckled Mn kagome lattice. On heating, the tetragonal distortion and octahedral tilt angle decrease continuously before transitioning to the cubic 𝑃⁢𝑚⁢3⁢𝑚 antiperovskite phase at 𝑇≈527K. Neutron diffraction and magnetometry together reveal noncollinear ferrimagnetic ordering. For 30K≤𝑇≤500K, the magnetic structure is described by magnetic space group 𝐼⁢𝑏⁢𝑎′⁢𝑚′ (72.544), with inequivalent Mn1 and Mn2 sublattices that couple antiferromagnetically to yield a net moment. Density-functional theory-based calculations show that the different local moments originate from the bandwidths associated with distinct Mn–N bond lengths. Temperature-dependent refinements reveal distinct differences in the thermal disordering profiles of the Mn1 and Mn2 sublattices. These findings reveal a subtlety in the magnetic and structural behavior of M⁢n3⁢GeN, highlighting the interplay between structural distortions, magnetic ordering, and electronic structure in kagome-derived antiperovskite materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Thermochemical Equilibria with Interacting Defects: Sr 1 − x Ce x Mn O 3 − δ Alloys for Water Splitting

Solar thermochemical hydrogen is one of the few potential routes towards direct fuel production from renewable energy sources, but the thermodynamic boundary conditions for efficient and economic energy conversion are challenging. Success or failure of a given oxide working material depends on the subtle balance between enthalpy and entropy contributions in the redox processes. Developing a mechanistic understanding of the behavior of materials on the basis of atomistic models and first-principles calculations is an important part of advancing the technology. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox-active defects start to interact with each other, thereby impeding the formation of additional defects. This problem is of more general importance to applications that rely on high levels of off-stoichiometry or doping, including, for example, batteries, thermoelectrics, and ceramic fuel cells. To account for such repulsive defect interactions, we introduce a statistical mechanics approach, defining an expression for the free energy of defect interaction based on limited sampling of defect configurations in density functional theory supercell calculations. The parameterization of this energy contribution as a function of defect concentration and temperature allows on-the-fly simulation of thermochemical equilibria. The approach consistently incorporates finite temperature effects by including the leading contributions to the temperature-dependent free energy for the case at hand, i.e., the ideal gas and configurational enthalpies and entropies. We demonstrate the capability and utility of the approach by simulating the water splitting redox processes for Sr 1 − x Ce x Mn O 3 − δ alloys. Published by the American Physical Society 2024

08 HYDROGEN↗
DOE OSTI2023

Study of the structure, structural transition, interface model, and magnetic moments of CrN grown on MgO(001) by molecular beam epitaxy

Structural phase transition is studied in high quality CrN thin films grown by molecular beam epitaxy on MgO(001) substrates. Cross-sectional transmission electron microscopy and x-ray diffraction reveal that the epitaxial relationship between CrN film and MgO substrate is [100] CrN /[100] MgO , [110] CrN /[110] MgO , and [001] CrN /[001] MgO . The films show tensile strain/compression at the CrN/MgO(001) interface, which relaxes gradually with the film growth. Temperature dependent x-ray diffraction measurements show a first-order structural phase transition. In addition to the experimental measurements, first-principles theoretical calculations have been carried out for finding a stable model for the CrN/MgO interface. Furthermore, these calculations determine two possible models for the interface, where a monolayer of chromium oxide is formed between the CrN and MgO layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electronic structure of van der Waals ferromagnet CrI 3 from self-consistent vertex corrected GW approaches

Electronic structure of the layered van der Waals ferromagnet CrI 3 is studied with self consistent diagrammatic approaches beyond the GW approximation. Considerable improvement in the calculated band gap as compared to the non-self-consistent G0W0 results has been found. Certain spectral features in the valence bands discovered recently by the angle resolved photoemission spectroscopy, are reproduced better when we use full frequency dependent self energy. Density functional theory and the quasiparticle self-consistent GW method which are based on the frequency-independent self energy are unable to resolve these features. The non-locality effects in the diagrams beyond the GW approximation are large for both polarizability and self energy. This finding can potentially have an impact on the development of methods like the GW+DMFT (dynamical mean field theory).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Influence of crystalline order and defects on the absolute work functions and electron affinities of TiO2- and SrO-terminated n-SrTiO3(001)

We establish how the detailed structure and composition of the terminal layer of n-SrTiO3 (001) determine key surface electronic properties. We have measured and calculated from first principles the absolute work functions and electron affinities of bulk SrNb0.01Ti0.99O3(001) with both TiO2 and nominally SrO terminal planes. The match between theory and experiment is excellent at 4.6(1) eV for the TiO2 termination if an ideal, bulk-truncated surface structure is assumed. In contrast, the ideal SrO-termination leads to a calculated work function that is ~1.3 eV lower than the experimental value of 3.6(1) eV. We show that this discrepancy is associated with defects on the SrO surface that serve as electron scavengers. These defects deplete the concentration of itinerant electrons in the subsurface region and increase the negative charge density on the surface, thus increasing the work function. Several different surface defect configurations were modeled; the ones that yield the best agreement with experiment involve Sr vacancies in the terminal layer along with OH–/H– pairs or O22– occupying anion sites adjacent to the Sr vacancies.

Chambers, Scott A.↗
DOE OSTI2024

Extreme compression of planetary gases: High-accuracy pressure-density measurements of hydrogen-helium mixtures above fourfold compression

Hydrogen (H 2 ) and helium (He), the most abundant elements in the universe, pose a unique challenge in measuring the equation of state of the mixture, owing to their differing physical properties. There remains a need for data with high enough precision to discriminate between existing equation of state (EOS) mix models in order to understand the internal structure of gas-giant planets. Here, we have measured the EOS of precompressed H 2 - He mixtures at conditions directly relevant to the planetary interiors using hypervelocity gas guns and Sandia’s Z machine with less than 10% uncertainty in density, enabling validation of mixture models. We precompressed 50:50 molar mixtures of H 2 -He to 0.1–0.2 GPa and directly measured particle velocity (in gas-gun experiments) and shock velocities (in Z-machine experiments). To complement the experimental efforts, we also computed the Hugoniots of precompressed H 2 -He mixtures using density-functional-theory-based molecular dynamics. Furthermore, we observe approximately 3- to 4.3-fold density compression at pressures up to 44 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Assessment of empirical interatomic potential to predict thermal conductivity in ThO 2 and UO 2

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. In this study we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the U x Th 1- x O 2 system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO 2 involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO 2 , to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO 2 ) or experimental measurements (UO 2 ). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Surface intermediates identified using IR spectroscopy during the catalytic oxidation of C 2 H 6 on IrO 2 (1 1 0) at near-ambient pressures

Characterizing surface intermediates during catalysis is essential for developing mechanistic models of catalytic reactions and for identifying rate-controlling steps. Here, in this work, we used polarization-dependent reflection absorption infrared spectroscopy (PD-RAIRS) to investigate surface intermediates that form on IrO 2 (1 1 0) during the catalytic oxidation of C 2 H 6 at pressures near 0.5 Torr. Our results show that adsorbed CO and HCO 2 are the only intermediates detectable with PD-RAIRS across a wide range of temperatures and reactant compositions during complete C 2 H 6 oxidation, giving rise to dominant peaks at about 2093 and 1310 cm −1 , respectively. These assignments were corroborated by RAIRS measurements following CO and HCOOH adsorption in ultrahigh vacuum, RAIRS with C 2 H 6 – 18 O 2 mixtures during catalysis and vibrational frequencies and IR intensities computed using density functional theory. We find that increasing the C 2 H 6 mole fraction in the reaction mixture raises the CO surface coverage relative to HCO 2 , correlating with increasing catalytic oxidation rates. The formation of CO and HCO 2 surface species at temperatures well below catalytic ignition indicates that the overall rate of C 2 H 6 oxidation on IrO 2 (1 1 0) is governed by steps late in the reaction sequence, such as H 2 O formation, CO oxidation and HCO 2 dehydrogenation. These findings provide key insight into the mechanism of alkane oxidation on IrO 2 (1 1 0) and valuable input for first-principles microkinetic modeling.

Pope, Connor [Univ. of Florida, Gainesville, FL (U↗
DOE OSTI2025

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

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

Mitigation of the internal p-n junction in CoS 2 -contacted FeS 2 single crystals: Accessing bulk semiconducting transport

Pyrite FeS 2 is an outstanding candidate for a low-cost, nontoxic, sustainable photovoltaic material, but efficient pyrite-based solar cells are yet to materialize. Recent studies of single crystals have shed much light on this by uncovering a p-type surface inversion layer on n-type (S-vacancy doped) crystals, and the resulting internal p-n junction. This leaky internal junction likely plays a key role in limiting efficiency in pyrite-based photovoltaic devices, also obscuring the true bulk semiconducting transport properties of pyrite crystals. Here, we demonstrate complete mitigation of the internal p-n junction in FeS 2 crystals by fabricating metallic CoS 2 contacts via a process that simultaneously diffuses Co (a shallow donor) into the crystal, the resulting heavy n doping yielding direct Ohmic contact to the interior. Low-temperature bulk transport studies of controllably Co- and S-vacancy doped semiconducting crystals then enable a host of previously inaccessible observations and measurements, including determination of donor activation energies (which are as low as 5 meV for Co), observation of an unexpected second activated transport regime, realization of electron mobility up to 2100 cm 2 V –1 s –1 , elucidation of very different mobilities in Co- and S-vacancy-doped cases, and observation of an abrupt temperaturedependent crossover to bulk Efros-Shklovskii variable-range hopping, accompanied by an unusual form of nonlinear Hall effect. Aspects of the results are interpreted with the aid of first-principles electronic structure calculations on both Co- and S-vacancy-doped FeS 2 . Furthermore, this work thus demonstrates unequivocal mitigation of the internal p-n junction in pyrite single crystals, with important implications for both future fundamental studies and photovoltaic devices

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

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗
DOE OSTI2021

Large-Scale, GPU-Enhanced DFTB Approaches for Probing Multi-Component Alloys

The major goals of the project are to develop and apply GPU-enhanced density functional theory tight binding (DFTB) calculations for efficient simulations of structural materials of alloy systems. The figure below depicts the GPU algorithmic developments on the left that will be carried out, and the solid-state structures on the right are representative of the complex alloys that can be computed with the GPU-enhanced DFTB approach. The methodology and tools developed in this project encompass accurate intermolecular potentials and GPU enhancements to the density functional tight binding (DFTB) approach for high-throughput ab initio molecular dynamics calculations of multi-component alloys at elevated temperatures. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of large, multi-component alloys at the predictive quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum nature of chemical systems; however, these methods cannot tackle the large sizes relevant to these multi-component systems. The DFTB-based ab initio molecular dynamics approach (coupled with our in-house GPU capabilities for enhanced speed) used in this project provides a viable approach for probing these large systems at a quantum level of detail that is significantly faster than current first-principle methods Collectively, the capabilities developed in this project directly respond to DOE HBCU-OMI, AOI 2-1 initiatives by (1) enabling accurate and efficient predictions and (2) bringing a fundamental understanding of structural interactions in these complex systems at elevated temperatures.

20 FOSSIL-FUELED POWER PLANTS↗