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 145 records · Page 8

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↗

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↗

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↗

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↗

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↗

Samarium: from a distorted-fcc phase to melting under dynamic compression using in-situ x-ray diffraction

Abstract Lattice and electronic structure interactions for f -electrons are fundamental challenges for lanthanide equation of state development. Difficulties in first-principles calculations, such as density functional theory (DFT), emphasize the need for well-characterized experimental data. Here, we measure in-situ x-ray diffraction of shocked samarium (Sm) and temperature along the Hugoniot for the first time, providing direct evidence for phase transitions. We report direct evidence of a distorted fcc (dfcc) phase at 23 GPa. Shocked samarium melts from the dfcc phase starting at 33 GPa (1333 K), with complete melt at 40 GPa (1468 K). Previous work indicated shock melt at 27 GPa (1200 K), underscoring the significance of x-ray measurements for detecting phase transitions. Interestingly, our observed melting is in sharp contrast with the melting reported by a diamond anvil cell study. These experimental data can tightly constrain first principles calculations and serve as key touchstones for equation of state modeling.

36 MATERIALS SCIENCE↗

HfTi_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-titanium (Hf-Ti). 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 Ti. 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↗

NbTi_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 niobium-titanium (Nb-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP), https://github.com/ORNL/ScalableWorkflow_VASP_Calculations. 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 Nb and Ti. 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↗

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗

Search for active and inactive ion insertion sites in organic crystalline materials

The position of mobile active and inactive ions, specifically ion insertion sites, within organic crystals, significantly affects the properties of organic materials used for energy storage and ionic transport. Identifying the positions of these atomic (and ionic) sites in an organic crystal is challenging, especially when the element has low X-ray scattering power, such as lithium (Li) and hydrogen, which are difficult to detect by powder X-ray diffraction. First-principles calculations, exemplified by density functional theory (DFT), are very practical for confirming the relative stability of ion positions in materials. However, the lack of effective strategies to identify ion sites in these organic crystalline frameworks renders this task extremely challenging. This work presents two algorithms: the (i) efficient location of ion insertion sites from extrema in electrostatic local potential and charge density (ELIISE), and the (ii) ElectRostatic InsertioN (ERIN), which leverage charge density and electrostatic potential fields accessed from first-principles calculations, combined with the Simultaneous Ion Insertion and Evaluation (SIIE) workflow. SIIE inserts all ions simultaneously—to determine ion positions in organic crystals. We demonstrate that these methods accurately reproduce known ion positions in 16 organic materials and identify previously overlooked low-energy sites in tetralithium 2,6-naphthalenedicarboxylate (Li4NDC), an organic electrode material, highlighting the importance of inserting all ions simultaneously, as in the SIIE workflow. These algorithms are also integrated with off-the-shelf machine learning potentials, yielding promising results comparable to first-principles findings.

Gopidi, Harshan Reddy [Argonne National Laboratory↗

A Polar Magnetic and Insulating Double Corundum Oxide: Mn2MnSbO6 with Ordered Mn(II) and Mn(III) Ions

A new magnetic insulator Mn2MnSbO6 with a polar crystal structure and an ordered Mn2+ and Mn3+ arrangement was synthesized under a high pressure of 7.5 GPa and 1300 degrees C. The crystal structure of Mn2MnSbO6, investigated by synchrotron powder X-ray diffraction, was found to be isomorphous with that of Ni3TeO6-type, space group R3. The non-centrosymmetric structure was confirmed by the second-harmonic generation measurements. The X-ray absorption near-edge spectroscopy measurement confirmed the nominal oxidation states of Mn22+Mn3+SbO6. Magnetic measurements indicate that Mn2MnSbO6 orders antiferromagnetically below 44 K and undergoes a field-induced spin-flop transition at 5 K. First-principles calculations indicate an antiferromagnetic ground state with up/up/up/down/down/down (uuuddd) spin configuration of the six crystallographically unique Mn ions in the c-axis doubled magnetic structure. The density functional theory calculations also substantiate the experimentally observed charge ordering of the Mn2+/Mn3+ ions and the insulating behavior due to a bandgap of 0.52 eV. To the best of our knowledge, this is the first double corundum oxide containing Jahn-Teller active Mn3+ ions.

Feng, Hai L.↗

Origins of heat transport anisotropy in MoTe2 and other bulk van der Waals materials

Knowledge of how heat flows anisotropically in van der Waals (vdW) materials is crucial for thermal management of emerging 2D materials devices and design of novel anisotropic thermoelectric materials. Despite the importance, anisotropic heat transport in vdW materials is yet to be systematically studied and is often presumably attributed to anisotropic speeds of sound in vdW materials due to soft interlayer bonding relative to covalent in-plane networks of atoms. In this work, we investigate the origins of the anisotropic heat transport in vdW materials, through time-domain thermoreflectance (TDTR) measurements and first-principles calculations of anisotropic thermal conductivity of three different phases of MoTe 2 . MoTe 2 is ideal for the study due to its weak anisotropy in the speeds of sound. We find that even when the speeds of sound are roughly isotropic, the measured thermal conductivity of MoTe 2 along the c-axis is 5–8 times lower than that along the in-plane axes. We derive meaningful characteristic heat capacity, phonon group velocity, and relaxation times from our first principles calculations for selected vdW materials (MoTe 2 , BP, h-BN, and MoS 2 ), to assess the contributions of these factors to the anisotropic heat transport. Interestingly, we find that the main contributor to the heat transport anisotropy in vdW materials is anisotropy in heat capacity of the dominant heat-carrying phonon modes in different directions, which originates from anisotropic optical phonon dispersion and disparity in the frequency of heat-carrying phonons in different directions. Further, the discrepancy in frequency of the heat-carrying phonons also leads to ~2 times larger average relaxation times in the cross-plane direction, and partially explains the apparent dependence of the anisotropic heat transport on the anisotropic speeds of sound. This work provides insight into understanding of the anisotropic heat transport in vdW materials.

36 MATERIALS SCIENCE↗

MnSi 2 Te 4 : A van der Waals Antiferromagnetic Semiconductor with Large Negative Magnetoresistance

Magnetism in van der Waals semiconductors offers significant potential for fundamental research on low-dimensional magnetism and the development of high-performance two-dimensional spintronic devices. Here, we report the growth, physical properties, and first-principles calculations of a new dual-octahedral transition metal chalcogenide (DTMC) MnSi 2 Te 4 . MnSi 2 Te 4 features a layered structure with an intralayer heterostructure, where the metal octahedra and nonmetal dimeric octahedra form zigzag chains alternately. Property characterization reveals that MnSi 2 Te 4 is a collinear G-type antiferromagnetic semiconductor, with a Néel temperature T N of 18.6 K and a significant unsaturated negative magnetoresistance (NMR) reaching −42.5% at 9 T and 100 K. First-principles calculations on the electronic band structure demonstrate that the large NMR primarily originates from the spin splitting due to parity-time symmetry breaking. This study not only discovers a new member of DTMCs with a unique crystal structure and large NMR, but also establishes a promising platform for investigating next-generation spintronic devices.

Liao, Ke [Chinese Academy of Sciences (CAS), Beiji↗

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗

Selective Antisite Defect Formation in WS 2 Monolayers via Reactive Growth on Dilute W-Au Alloy Substrates

Defects are ubiquitous in 2D materials and can affect the structure and properties of the materials and also introduce new functionalities. Methods to adjust the structure and density of defects during bottom-up synthesis are required to control the growth of 2D materials with tailored optical and electronic properties. In this study, the authors present an Au-assisted chemical vapor deposition approach to selectively form S W and S2 W antisite defects, whereby one or two sulfur atoms substitute for a tungsten atom in WS 2 monolayers. Guided by first-principles calculations, they describe a new method that can maintain tungsten-poor growth conditions relative to sulfur via the low solubility of W atoms in a gold/W alloy, thereby significantly reducing the formation energy of the antisite defects during the growth of WS 2 . The atomic structure and composition of the antisite defects are unambiguously identified by Z-contrast scanning transmission electron microscopy and electron energy-loss spectroscopy, and their total concentration is statistically determined, with levels up to ≈5.0%. Furthermore, scanning tunneling microscopy/spectroscopy measurements and first-principles calculations also verified these antisite defects and revealed the localized defect states in the bandgap of WS 2 monolayers. This bottom-up synthesis method to form antisite defects should apply in the synthesis of other 2D materials.

2D materials↗

Atomic structures of interfacial solute gateways to θ' precipitates in Al-Cu alloys

Many materials employed in critical structural applications depend upon metastable strengthening precipitates that transform or dissolve at elevated temperatures. Herein, aberration-corrected scanning transmission electron microscopy and first-principles calculations are used to accurately determine the atomic structure of the highly mobile, semi-coherent precipitate interfaces that control this process in the classic θ' (Al 2 Cu) precipitate in the Al-Cu system. Semi-coherent {110} interfaces are found to be composed of an array of unexpected misfit dislocations that are arranged in two different structural units. Dislocations accommodate nearly all of the misfit between the Al matrix and strengthening phase. Cu is observed to segregate to the compressed edge of the dislocation cores at specific sites in this interface. First-principles calculations revealed the energetic landscape that facilitates these sites to become entry and exit gateways of Cu atoms in this semi-coherent interface. In conclusion, this investigation reveals critical features within semi-coherent interfaces that determine the thermal stability of precipitation-hardened alloys.

36 MATERIALS SCIENCE↗

High-Pressure Synthesis of Bulk Cobalt Cementite, Co 3 C

Transition metal carbides find use in a wide range of advanced high-resilience applications including high-strength steels, heat shields, and deep-earth drills. However, carbides of the mid-to-late transition metals remain difficult to isolate and characterize on account of their metastability, which precludes the preparation of high-quality bulk single crystal samples using traditional solid-state methods. Herein, we report a combined computational and experimental survey of the cobalt-carbon binary system under high pressures and demonstrate that pressure offers a route toward the bulk synthesis of the metastable cementite-type cobalt carbide, Co 3 C, which under ambient conditions can only be prepared in low-dimensional thin film or nanoparticle forms. First-principles calculations reveal two competitive low-energy stoichiometric phases under ambient pressures - Pnnm-Co 2 C (Fe 2 C-type) and Pnma-Co 3 C (Fe 3 Ctype) - consistent with the known low-dimensional phases that have been studied for their promising magnetic properties. However, the calculated formation enthalpy of Pnma-Co 3 C decreases steadily with the applied pressure, while that of Pnnm-Co 2 C increases. Here, we pursue these results using high-pressure laser-heated synthesis methods coupled with in situ X-ray diffraction and observe the formation of Pnma-Co 3 C above 4.8 GPa. We determine the experimental bulk modulus of Co 3 C to be K o = 237 GPa (K p = 4.0). First-principles calculations of the phonon modes in Co 3 C reveal dynamical instabilities at ambient pressure that are absent under compression. These results offer a promising new route for the synthesis of rare-earth-free magnets.

36 MATERIALS SCIENCE↗

Machine learning of 27Al NMR electric field gradient tensors for crystalline structures from DFT

NMR crystallography has emerged as a promising technique for the determination and refinement of atomic coordinates in crystal structures. The crystal structure of compounds containing quadrupolar nuclei, such as 27Al, can be improved by directly comparing solid-state NMR measurements to DFT computations of the electric field gradient (EFG) tensor. The non-negligible computational cost of these first-principles calculations limits the applicability of this method to all but the most well-defined structures. We developed a fast, low-cost machine learning model to predict EFG parameters based on local structural motifs and elemental parameters. We computed 8081 EFG tensors from 1681 27Al crystalline solids using DFT and benchmarked them against 105 experimentally measured 27Al sites. Surprisingly, simple local geometric features dominate the predictive performance of the resulting random-forest model, yielding an R2 value of 0.98 and an RMSE of 0.61 MHz for CQ, the quadrupolar coupling constant. This model accuracy should enable pre-refining future structural assignments before finally validating with first-principles calculations. Such a catalogue of 27Al NMR tensors can serve as a tool for researchers assigning complex NMR spectra influenced by the nuclear electric quadrupole interaction.

Sun, He↗