Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Controlling topology through targeted composite symmetry manipulation in magnetic systems

The possibility of selecting magnetic space groups by orienting the magnetization direction or tuning magnetic orders offers a vast playground for engineering symmetry-protected topological phases in magnetic materials. In this work, we study how selective tuning of symmetry and magnetism can influence and control the resulting topology in a two-dimensional magnetic system, and we illustrate such a procedure in the ferromagnetic monolayer MnPSe 3 . Density functional theory calculations reveal a symmetry-protected accidental semimetallic (SM) phase for out-of-plane magnetization, which becomes an insulator when the magnetization is tilted in-plane, reaching band-gap values close to 100 meV. We identify an order-2 composite antiunitary symmetry and threefold rotational symmetry that induce the band crossing, and we classify the possible topological phases using symmetry analysis, which we support with tight-binding and k · p models. Breaking of inversion symmetry opens a gap in the SM phase, giving rise to a Chern insulator. We demonstrate this explicitly in the isostructural Janus compound Mn 2 ⁢P 2 ⁢S 3 ⁢Se 3 , which naturally exhibits Rashba spin-orbit coupling that breaks inversion symmetry. Our results map out the phase space of topological properties of ferromagnetic transition-metal phosphorus trichalcogenides, and they demonstrate the potential of the magnetization-dependent metal-to-insulator transition as a spin switch in integrated two-dimensional electronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
DOE OSTI2026

Experimental confirmation of the magnetic ordering transition induced by an electronic structure change in the metallic triangular antiferromagnet Co 1/3⁢ TaS 2

We report angle-resolved photoemission spectroscopy (ARPES) studies combined with DFT+DMFT calculations to confirm that the magnetic ordering vector transition from Q=(1/2,0,0) to Q=(1/3,0,0) in the metallic triangular antiferromagnets Co 1/3±𝜀⁢ TaS 2 (𝜀 ≈ 0.007) is induced by the electronic structure change in the system. The ARPES-measured Fermi surface (FS) maps of Co 0.325 ⁢TaS 2 show two hexagonal and one circular holelike FSs around Γ, which matches well with the triple- Q state by taking into account the contribution of nesting vectors occurring between Co 3⁢𝑑 and Ta 5⁢𝑑 orbitals. In the case of Co 0.340 ⁢TaS 2 , a new electron pocket around K appears, and the FS geometry changes as a result of the correlation effect of Co 4 ⁢S 18 tripods forming in the system. The magnetic susceptibility calculations based on the charge-self-consistent DFT+DMFT band structures and the random phase approximation indicate that the most stable magnetic ordering vector (1/2,0,0) split into (1/6,0,0) and (1/2,0,0), which is consistent with the magnetic phase transition around 𝑥=1/3 in Co 𝑥 ⁢TaS 2 .

Fermi surface↗
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 OSTI2020

Photoluminescence, photophysics, and photochemistry of the V B - defect in hexagonal boron nitride

Extensive photochemical and spectroscopic properties of the V B - defect in hexagonal boron nitride are calculated, concluding that the observed photoemission associated with recently observed optically detected magnetic resonance is most likely of (1) 3 E'' → (1) 3 A' 2 origin. Rapid intersystem crossing from the defect's triplet to singlet manifolds explains the observed short excited-state lifetime and very low quantum yield. New experimental results reveal smaller intrinsic spectral bandwidths than previously recognized, interpreted in terms of spectral narrowing and zero-phonon-line shifting induced by the Jahn-Teller effect. Different types of computational methods are applied to map out the complex triplet and singlet defect manifolds, including the doubly ionized formulation of the equation-of-motion coupled-cluster theory that is designed to deal with the open-shell nature of defect states, and mixed quantum-mechanics/molecular-mechanics schemes enabling 5763-atom simulations. Finally, two other energetically feasible spectral assignments from among the singlet and triplet manifolds are considered, but ruled out based on inappropriate photochemical properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. 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. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
DOE OSTI2024

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermally induced structural competitiveness and metastability of body-centered-cubic iron under nonequilibrium conditions

The structure and stability of iron near melting at multi-megabar pressures are of significant interest in high pressure physics and earth and planetary sciences. While the body-centered cubic (BCC) phase is generally recognized as unstable at lower temperatures, its stability relative to the hexagonal close-packed (HCP) phase at high temperatures (approximately 0.5 eV) in the Earth’s inner core (IC) remains a topic of ongoing theoretical and experimental debate. Motivated by a series of studies attempting to understand the state of the IC and resulting in contradictory findings, we aim to explore the physical factors that may contribute to this confusion in a broader context. Our ab initio calculations show a significant drop in energy, the emergence of a plateau and a local minimum in the potential energy surface, and stabilization of all phonon modes of BCC iron at elevated electron temperatures (>1–1.5 eV). These effects increase the competition among the BCC, HCP, and the face-centered cubic (FCC) phases and lead to the metastability of the BCC structure. Furthermore, at pressure-temperature conditions relevant to the IC, we can expect substantial coexistence of FCC and HCP iron due to their similar free energies and low transition barrier (∼0.2 eV/atom), and the thermodynamic stability of BCC iron is enhanced by its substantial lattice vibration entropy. In conclusion, these findings provide a clear theoretical framework for understanding iron phase relations and solidification processes, both experimentally and in the IC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Theoretical study of the crystal and electronic properties of α - RuI 3

The material α-RuCl 3 , with a two-dimensional Ru honeycomb sublattice, has attracted considerable attention because it may be a realization of the Kitaev quantum spin liquid. Recently, a new honeycomb material, α-RuCl 3 , was prepared under moderately high pressure, and it is stable under ambient conditions. However, different from α-RuCl 3 , α-RuI 3 was reported to be a paramagnetic metal without long-range magnetic order down to 0.35 K. Here, the structural and electronic properties of the quasi-two-dimensional α-RuI 3 are theoretically studied. First, based on first-principles density functional theory calculations, the ABC stacking honeycomb-layer R$\bar{3}$ (No. 148) structure is found to be the most likely stacking order for α-RuI 3 along the c axis. Furthermore, both R$\bar{3}$ and P$\bar{3}$1c are dynamically stable because no imaginary frequency modes were obtained in the phononic dispersion spectrum without Hubbard U. Moreover, the different physical behavior of α-RuI 3 compared to α-RuCl 3 can be understood naturally. The strong hybridization between Ru 4d and I 5p orbitals decreases the “effective” atomic Hubbard repulsion, leading the electrons of RuI3 to be less localized than in RuCl 3 . As a consequence, the effective electronic correlation is reduced from Cl to I, leading to the metallic nature of α-RuI 3 . Based on the DFT+U (Ueff=2 eV) plus spin-orbital coupling, we obtained a spin-orbit Mott insulating behavior for α-RuCl 3 and, with the same procedure, a metallic behavior for α-RuI 3 , in good agreement with experimental results. Furthermore, when introducing large (unrealistic) U eff =6 eV, the spin-orbit Mott gap opens in α-RuI 3 as well, supporting the physical picture we are proposing. Our results provide guidance to experimentalists and theorists working on two-dimensional transition metal tri-iodide layered materials.

2-dimensional systems↗
DOE OSTI2025

Weyl semimetal phases and intrinsic spin-Hall conductivity in SbAs ordered alloys

Here, using density functional theory calculations, we investigated possible Weyl semimetal (WSM) phases in antimony arsenide ordered alloys Sb 1-x ⁢As x (x=0, 1/6, 1/3, 1/2, 2/3, 5/6, 1). We find WSM phases for all As compositions of Sb 1-x ⁢As x with broken inversion symmetry, in contrast to Bi 1-x⁢ Sb x where only compositions x=1/2 and 5/6 were predicted to exhibit WSM phases. The WSM phases in Sb 1-x ⁢As x are characterized by the presence of 12 Weyl points, located within 55 meV from the Fermi level in the case of x = 1/2. The robust spin-orbit coupling strength and Berry curvature in these alloys produce large spin-Hall conductivity in the range of 176–602 ($\hslash$/e)(S/cm), comparable to that in the BiSb alloys. Finally, Sb 0.5 ⁢As 0.5 is predicted to be almost lattice matched to GaAs(111), with the Fermi level within the gap of the semiconductor, facilitating growth and characterization, and thus, offering promising integration with conventional semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Superconductivity in CH 4 and BH – 4 containing compounds derived from the high-pressure superhydrides

Inspired by the synthesis of the high-pressure Fm3m LaH 10 superconducting superhydride, systematic density functional theory (DFT) calculations are performed to study ternaries that could be derived from it by replacing two of the hydrogen atoms with boron or carbon and varying the identity of the electropositive element. Though many of the resulting alkali-metal and alkaline-earth MC 2 H 8 phases are predicted to be dynamically stable at mild pressures, their superconducting critical temperatures (T c s) are low because their metallicity results from the filling of an electride-like band. Substitution with a trivalent element leads to phases with substantial metal d- character at the Fermi level whose T c s are typically above 40 K. Here, among the MB 2 H 8 phases examined, KB 2 H 8 , RbB 2 H 8 and CsB 2 H 8 are predicted to be dynamically stable at very mild pressures, and their stability is rationalized by a DFT-Chemical Pressure analysis that elucidates the role of the M atom size. Quantum anharmonic effects strongly affect the properties of KB 2 H 8 , the highest predicted T c compound, near 10 GPa, but molecular dynamics simulations reveal it would decompose below its T c at this pressure. Nonetheless, at ca. 50 GPa KB 2 H 8 is predicted to be thermally stable with a superconducting figure of merit surpassing that of the recently synthesized LaBeH 8 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Anisotropy, frustration, and saddle point in the twisted kagome compound ErPdPb

The kagome lattice, with its inherent geometric frustration, provides a rich platform for exploring intriguing magnetic phenomena and topological electronic structures. In reduced-symmetry structures, such as twisted kagome systems involving rare earth elements, additional anisotropy can arise, enabling intriguing properties including spin-ice states, magnetocaloric effects, noncollinear magnetic ordering, and the anomalous Hall effect. Here, we report the synthesis of single crystals of ErPdPb, which features a twisted kagome lattice net of Er atoms within the hexagonal ZrNiAl-type structure, and we investigate its magnetic, electronic, and thermal properties. The material exhibits a highly anisotropic, correlated magnetic state below 2.7 K, as evidenced by magnetic, transport, and heat capacity measurements. Density functional theory (DFT) calculations indicate strong easy-axis anisotropy, consistent with experiment and crystal-field expectations, as well as quasi-one-dimensional bands and a spin-split saddle point at the zone center. Here, the coexistence of competing magnetic interactions along different crystallographic directions suggests an inherent degree of frustration. ErPdPb thus provides a promising platform for exploring the interplay of frustration, anisotropy, and electronic structure in a twisted kagome lattice.

Antiferromagnets↗
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

A Panoramic View of MXenes via an Atomic Coordination‐Based Design Strategy

Two‐dimensional (2D) transition metal carbides and nitrides, known as MXenes, possess unique physical and chemical properties, enabling diverse applications in fields ranging from energy storage to communication, catalysis, sensing, healthcare, and beyond. Despite extensive research and notable advancements, a fundamental understanding of MXenes’ phase diversity and its connection to their hierarchical precursors, including the intermediate MAX phases and the ancestral bulk phases, remains limited. Here, in this study, it is hypothesized that the atomic coordination environments adopted by transition metal and nonmetallic atoms in their three‐dimensional (3D) bulk precursors may persist in 2D MXenes to govern their phase diversity. Using high‐throughput modeling based on first‐principles density functional theory, a wide range of MXene phases is unveiled and comprehensively evaluate their relative stabilities across a large chemical space. The key to the approach lies in considering various atomic coordination environments drawn from four types of ancestral bulk phases. Through this comprehensive structural library of MXenes, general guiding principles are uncovered, such as a close alignment between the phase stability of MXenes and that of their 3D precursors. These findings introduce a new design strategy in which the atomic coordination environments in bulk phases can serve as reliable predictors for accessing the diverse structural landscape of MXenes.

MXenes↗
DOE OSTI2022

Electrical Manipulation of Topological Phases in a Quantum Anomalous Hall Insulator

Quantum anomalous Hall phases arising from the inverted band topology in magnetically doped topological insulators have emerged as an important subject of research for quantization at zero magnetic fields. Though necessary for practical implementation, sophisticated electrical control of molecular beam epitaxy (MBE)-grown quantum anomalous Hall matter have been stymied by growth and fabrication challenges. Here, in this study, a novel procedure is demonstrated, employing a combination of thin-film deposition and 2D material stacking techniques, to create dual-gated devices of the MBE-grown quantum anomalous Hall insulator, Cr-doped (Bi,Sb) 2 Te 3 . In these devices, orthogonal control over the field-induced charge density and the electric displacement field is demonstrated. A thorough examination of material responses to tuning along each control axis is presented, realizing magnetic property control along the former and a novel capability to manipulate the surface exchange gap along the latter. Through electrically addressing the exchange gap, the capabilities to either strengthen the quantum anomalous Hall state or suppress it entirely and drive a topological phase transition to a trivial state are demonstrated. The experimental result is explained using first principle theoretical calculations, and establishes a practical route for in situ control of quantum anomalous Hall states and topology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Design of a Stable Heusler Alloy with Switchable Metal-to-Half-Metal Transition at Finite Temperature

We report the electronic structure and metallicity of Si-doped Mn2VGe Heusler alloys within the first-principles density-functional theory framework are discussed. Mn2VGe is found to assume two stable structures at (cubic) lattice constants of 5.7 and 6.05 Å, which are well-separated in energy and correspond to half-metallic (low-spin) and metallic (high-spin) phases, respectively. Substitution of Ge by Si reduces the energy difference between these two phases, which become nearly degenerate at an Si concentration of ≈3.125 at% at zero temperature. The analysis shows that the switching between the high- and low-spin phases for this composition can be triggered via a pressure of 2.8 GPa at 300 K. Si-doped Mn 2 VGe can thus provide an ultrafast, low power, cost-effective materials platform for spintronics applications.

36 MATERIALS SCIENCE↗