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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Thermodynamics of MgO Atomic Layer Deposition Surface Reactions

The realities of atomic layer deposition (ALD) surface reactions often deviate from the simple ligand exchange frequently used to illustrate the technique. A detailed understanding of these reactions is necessary to develop greater surface synthetic control including chemical selectivity for patterning or to target defects. Here, the thermodynamics of surface reactions relevant to MgO ALD were investigated using pyroelectric calorimetry to measure the time-resolved heat generation. These reactions show exothermic heat generation of 0.12 mJ/cm 2 and 0.15 mJ/cm 2 for alternating Mg(CpEt) 2 and H 2 O reactions, respectively. The total reaction heat closely matches the standard reaction enthalpy for bulk MgO, supplemented with first-principles molecular calculations. First-principles models further reveal that while simple ligand exchange is favorable during surface reactions, the required increase in Mg-coordination number from two in the precursor to six in bulk MgO requires additional coverage-dependent surface reactions, which depend on the availability and stability of proximal surface hydroxyls.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enhanced orbital anisotropy through the proximity to a SrTiO 3 layer in the perovskite iridate superlattices

We have used angle-dependent soft x-ray absorption spectroscopy (XAS) at the O K edge and first-principles calculations to investigate the electronic structures of iridate-based superlattices ( SrIr O 3 ) m / ( SrTi O 3 ) ( m = 1 , 2, 3, and ∞). We focus on the pre-edge Ir 5 d t 2 g – O 2 p orbital hybridization feature in the XAS spectra. By varying the measurement geometry relative to the incident photon polarization, we are able to extract the dichroic contrast and observe the systematic increase in the anisotropy of Ir 5 d orbitals as m decreases. First-principles calculations elucidate the orbital anisotropy coming mainly from the enhanced out-of-plane compression of Ir O 6 octahedra in the SrIr O 3 layers that are adjacent to the inserted SrTi O 3 layers. As m decreases, the increased volume fraction of these interfacial SrIr O 3 layers and their contact with the SrTi O 3 layers within the ( SrIr O 3 ) m / ( SrTi O 3 ) supercell lead to enhanced orbital anisotropy. Furthermore, the tilt and rotation of Ir O 6 octahedra are shown to be essential to understand the subtle orbital anisotropy in these superlattices, and constraining these degrees of freedom will give an incorrect trend. Overall, our results demonstrate that the structural constraint from the inserted SrTi O 3 layers, in addition to other electronic means such as polar interface and charge transfer, can serve as a knob to control the orbital degree of freedom in these iridate-based superlattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Machine Learning-Guided Discovery of Ternary Compounds Containing La, P, and Group 14 Elements

In this work, we integrate a deep machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary compounds. Using La–Si–P as a prototype system, we demonstrate that ML-guided first-principles calculations can efficiently explore crystal structures and their relative energetic stabilities, thus greatly accelerate the pace of material discovery. A number of new La–Si–P ternary compounds with formation energies less than 30 meV/atom above the known ternary convex hull are discovered. Among them, the formation energies of La 5 SiP 3 and La 2 SiP phases are only 2 and 10 meV/atom, respectively, above the convex hull. These two compounds are dynamically stable with no imaginary phonon modes. Moreover, by replacing Si with heavier-group 14 elements in the eight lowest-energy La–Si–P structures from our ML-guided predictions, a number of low-energy La–X–P phases (X = Ge, Sn, Pb) are predicted.

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

Biopolymer‐assisted Synthesis of P‐doped TiO 2 Nanoparticles for High‐performance Lithium‐ion Batteries: A Comprehensive Study

Abstract TiO 2 material has gained significant attention for large‐scale energy storage due to its abundant, low‐cost, and environmentally friendly properties, as well as the availability of various nanostructures. Phosphorus doping has been established as an effective technique for improving electronic conductivity and managing the slow ionic diffusion kinetics of TiO 2 . In this study, non‐doped and phosphorus doped TiO 2 materials were synthesized using sodium alginate biopolymer as chelating agent. The prepared materials were evaluated as anode materials for lithium‐ion batteries (LIBs). The electrodes exhibit remarkable electrochemical performance, including a high reversible capacity of 235 mAh g −1 at 0.1 C and excellent first coulombic efficiency of 99 %. An integrated approach, combining operando XRD and ex‐situ XAS, comprehensively investigates the relationship between phosphorus doping, material structure, and electrochemical performance, reinforced by analytical tools and first principles calculations. Furthermore, a full cell was designed using 2 %P‐doped TiO 2 anode and LiFePO 4 cathode. The output voltage was about 1.6 V with high initial specific capacity of 148 mAh g −1 , high rate‐capability of 120 mAh g −1 at 1 C, and high‐capacity retention of 96 % after 1000 cycles at 1 C.

El Halya, Nabil↗
DOE OSTI2023

Molecular beam epitaxy and crystal structure of majority a -plane-oriented and substrate-strained Mn 3 Sn thin films grown directly on sapphire (0001)

The Kagome antiferromagnet Mn 3 Sn has garnered much attention due to the presence of exciting properties such as anomalous Hall and Nernst effects. This paper discusses the synthesis of crystalline Mn 3 Sn thin films, prepared on Al 2 O 3 (0001) substrates at 453 ± 5°C using molecular beam epitaxy. The growth is monitored in situ using reflection high energy electron diffraction and measured ex situ using x-ray diffraction, Rutherford back-scattering, and cross-sectional scanning transmission electron microscopy. Our analysis shows the in-plane lattice constants of a 1,M = 4.117 ± 0.027 Å and a 2,M = 4.943 ± 0.033 Å, which is a very unexpected result when compared to the bulk a-plane Mn 3 Sn. This indicates a strain in the film and makes it challenging to provide a straightforward explanation. In an effort to explain our results, we discuss two possible orientation relationships between the Mn 3 Sn films and the sapphire substrates. Samples prepared under these conditions appear to have smooth surfaces locally, but overall the film has a 3D island morphology. First-principles calculations provide atomic models of the Mn 3 Sn (112¯0) lattice on Al 2 O 3 (0001) high symmetry sites, indicating that the L3-R90° is the most stable configuration. A detailed discussion of the experimental data and theoretical results, as well as strain effects, is provided.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Cation and anion topotactic transformations in cobaltite thin films leading to Ruddlesden-Popper phases.

Topotactic transformations involve structural changes between related crystal structures due to a loss or gain of material while retaining a crystallographic relationship. The perovskite oxide La0.7Sr0.3CoO3 (LSCO) is an ideal system for investigating phase transformations due to its high oxygen vacancy conductivity, relatively low oxygen vacancy formation energy, and strong coupling of the magnetic and electronic properties to the oxygen stoichiometry. While the transition between cobaltite perovskite and brownmillerite (BM) phases has been widely reported, further reduction beyond the BM phase lacks systematic studies. In this paper, we study the evolution of the physical properties of LSCO thin films upon exposure to highly reducing environments. We observe the rarely reported crystalline Ruddlesden-Popper phase, which involves the loss of both oxygen anions and cobalt cations upon annealing where the cobalt is found as isolated Co ions or Co nanoparticles. First-principles calculations confirm that the concurrent loss of oxygen and cobalt ions is thermodynamically possible through an intermediary BM phase. The strong correlation of the magnetic and electronic properties to the crystal structure highlights the potential of utilizing ion migration as a basis for emerging applications such as neuromorphic computing.

Chiu, I-Ting↗
NASA NTRSDate not supplied

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗
DOE OSTI2022

CASM — A software package for first-principles based study of multicomponent crystalline solids

CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fast and accurate machine learning prediction of phonon scattering rates and lattice thermal conductivity

Abstract Lattice thermal conductivity is important for many applications, but experimental measurements or first principles calculations including three-phonon and four-phonon scattering are expensive or even unaffordable. Machine learning approaches that can achieve similar accuracy have been a long-standing open question. Despite recent progress, machine learning models using structural information as descriptors fall short of experimental or first principles accuracy. This study presents a machine learning approach that predicts phonon scattering rates and thermal conductivity with experimental and first principles accuracy. The success of our approach is enabled by mitigating computational challenges associated with the high skewness of phonon scattering rates and their complex contributions to the total thermal resistance. Transfer learning between different orders of phonon scattering can further improve the model performance. Our surrogates offer up to two orders of magnitude acceleration compared to first principles calculations and would enable large-scale thermal transport informatics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

TiV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys titanium-vanadium (Ti-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ti and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of the third dimension in searching for Majorana fermions in α-RuCl 3 via phonons

Understanding phonons in α-RuCl 3 is critical to analyze the controversy around the observation of the half-integer thermal quantum Hall effect. While many studies have focused on the magnetic excitations in α-RuCl 3 , its vibrational excitation spectrum has remained relatively unexplored. We investigate the phonon structure of α-RuCl 3 via inelastic neutron-scattering experiments and density-functional-theory calculations. Our results show excellent agreement between experiment and first-principles calculations. After validating our theoretical model, we extrapolate the low-energy phonon properties. We find that the phonons in α-RuCl 3 that either propagate or vibrate in the out-of-plane direction have significantly reduced velocities and therefore have the potential to dominate the observability of the elusive half-integer plateaus in the thermal Hall conductance. In addition, we use low-energy interlayer phonons to resolve the low-temperature stacking structure of our large crystal of α-RuCl 3 , which we find to be consistent with that of the R$\bar{3}$ space group, in agreement with neutron diffraction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Machine Learning-Guided Exploration of Ternary Metal Borohydrides

We employ deep machine learning (ML) combined with first-principles calculations to explore energetically favorable ternary metal borohydrides. Using La–B–H as a prototype system, we demonstrate that iteratively trained ML models can efficiently screen hundreds of thousands of hypothetical structures and accurately select a small fraction of promising structures and compositions for further studies by first-principles calculations. Such an ML-guided approach dramatically accelerates the pace of materials discovery. A number of new La–B–H ternary compounds with formation energies within 100 meV/atom above the known ternary convex hull are discovered, including a known stable La(BH 4 ) 3 phase. Moreover, by replacing La with Group 1, 2, 3, 13, and 14 elements in the four lowest-energy La–B–H structures from our ML-guided predictions, several low-energy X–B-H (X = Mg, Ca, Sr, Ba, Sc, Y, Ac, Al, Ga, In, Si, Ge, Sn, Pb) compounds are predicted.

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

Influence of Annealing and Composition on the Crystal Structure of Mixed-Halide, Ruddlesden–Popper Perovskites

Mixed-halide two-dimensional (2D) hybrid organic–inorganic perovskites offer an important opportunity to control the band gap for applications in optoelectronic devices. In this study, we focus on phenethylammonium lead halide [(PEA) 2 Pb(I 1– x Br x ) 4 ] films, the pure iodide form of which is one of the most widely studied optically active 2D perovskite systems. Resonant infrared, matrix-assisted pulsed laser evaporation is used to grow films to explore the effects of post-growth annealing and mixed-halide composition. The composition, crystal structure, and optical properties are studied for as-grown and annealed films of (PEA) 2 Pb(I 1– x Br x ) 4 for x = 0, 0.25, 0.5, 0.75, and 1. First-principles calculations are used in conjunction with the experimental data to explain the mixed-halide behavior that does not trend monotonically with the bromide content. Important results of this work are as follows: (i) X-ray diffraction reveals evidence for halide phase separation around x = 0.25, consistent with first-principles calculations, whereas no phase separation is observed for x = 0.5 and above and (ii) a unique photoluminescence (PL) peak splitting is observed for x = 0.75, a composition for which no compositional phase separation is observed. The PL splitting is tentatively explained by the coexistence of two distinct types of halide short-range ordering at x = 0.75. Overall, this study demonstrates that mixed-halide n = 1 Ruddlesden–Popper perovskites are not simple random alloys but that instead, they display distinct sites and ordering preferences of the different halide anions. These preferences are critical to understand and rationally tune the properties of the materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The structure of sc16 GaP obtained at 17.5 GPa and 1400 K

We report using a laser-heated diamond anvil cell, a powdered sample of gallium phosphide was compressed to 17.5 GPa and heated up to 1400 K. The material obtained was characterized at room temperature using synchrotron x-ray diffraction. Experimental results were compared with first-principles calculations. The polymorph observed assumes the simple cubic structure (sc16) initially reported for GaAs at high pressures. Microdiffraction mapping showed variable grain sizes of the synthesized phase, with the largest grains located in the middle of the heated spot. Structural refinements were performed on selected grains. The structure predicted based on first-principles calculations is in close agreement with the experiments. The two Ga–P bonds show similar lengths in sc16-GaP; however, the bond angles differ, resulting in a distorted tetrahedral coordination geometry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Engineering Anomalously Large Electron Transport in Topological Semimetals

Abstract Anomalous transport of topological semimetals has generated significant interest for applications in optoelectronics, nanoscale devices, and interconnects. Understanding the origin of novel transport is crucial to engineering the desired material properties, yet their orders of magnitude higher transport than single‐particle mobilities remain unexplained. This work demonstrates the dramatic mobility enhancements result from phonons primarily returning momentum to electrons due to phonon‐electron dominating over phonon–phonon scattering. Proving this idea, proposed by Peierls in 1932, requires tuning electron and phonon dispersions without changing symmetry, topology, or disorder. This is achieved by combining de Haas ‐ van Alphen (dHvA), electron transport, Raman scattering, and first‐principles calculations in the topological semimetals MX 2 (M = Nb, Ta and X = Ge, Si). Replacing Ge with Si brings the transport mobilities from an order magnitude larger than single particle ones to nearly balanced. This occurs without changing the crystal structure or topology and with small differences in disorder or Fermi surface. Simultaneously, Raman scattering and first‐principles calculations establish phonon–electron dominated scattering only in the MGe 2 compounds. Thus, this study proves that phonon‐drag is crucial to the transport properties of topological semimetals and provides insight to engineer these materials further.

36 MATERIALS SCIENCE↗