Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Many-body theory of phonon-induced spin relaxation and decoherence

First-principles calculations enable accurate predictions of electronic interactions and dynamics. However, computing the electron spin dynamics remains challenging. The spin-orbit interaction causes various dynamical phenomena that couple with phonons, such as spin precession and spin-flip e-ph scattering, which are difficult to describe with current first-principles calculations. Here, in this work, we show a rigorous framework to study phonon-induced spin relaxation and decoherence, by computing the spin-spin correlation function and its vertex corrections due to e-ph interactions. We apply this approach to a model system and develop corresponding first-principles calculations of spin relaxation in GaAs. Our vertex-correction formalism is shown to capture the Elliott-Yafet, Dyakonov-Perel, and strong-precession mechanisms—three independent spin decoherence regimes with distinct physical origins—thereby unifying their theoretical treatment and calculation. Our method is general and enables quantitative studies of spin relaxation, decoherence, and transport in a wide range of materials and devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Structural and vibrational properties of methane up to 71 GPa

Single-crystal synchrotron x-ray diffraction, Raman spectroscopy, and first principles calculations have been used to identify the structure of the high-pressure (HP) phase of molecular methane above 20 GPa up to 71 GPa at room temperature. The structure of the HP phase is trigonal R3, which can be represented as a distortion of the cubic phase B, previously documented at 7–15 GPa and confirmed here. The positions of hydrogen atoms in the HP phase have been obtained from first principles calculations, which also demonstrated the stability of this structure above 260 K at 25 GPa. The molecules occupy four different crystallographic sites in phase B and 11 sites in the HP phase, which result in splitting of molecular stretching modes detected in Raman spectroscopy and assigned here based on a good agreement with the Raman spectra calculated from the first principles. Furthermore, our study points out to a single HP phase unlike up to three HP phases proposed previously based on the Raman spectroscopy results only.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗
DOE OSTI2025

Mapping the three-dimensional fermiology of the triangular lattice magnet EuAg 4 Sb 2

In this paper, we report the temperature-field phase diagram as well as present a comprehensive study of the electronic structure and three-dimensional fermiology of the triangular-lattice magnet EuAg 4 ⁢Sb 2 , utilizing quantum oscillation measurements, angle-resolved photoemission spectroscopy, and first-principles calculations. The complex magnetic phase diagram of EuAg 4 ⁢Sb 2 highlights many transitions through nontrivial AFM states. Shubnikov-de Haas and de Haas-van Alphen oscillations were observed in the polarized ferromagnetic state of EuAg 4 ⁢Sb 2 , revealing three pairs of distinct spin-split frequency branches with small effective masses. A comparison of the angle-dependent oscillation data with first-principles calculations in the ferromagnetic state and angle-resolved photoemission spectra shows good agreement, identifying tubular hole pockets and hourglass-shaped hole pockets at the Brillouin zone center, as well as diamond-shaped electron pockets at the zone boundary. As the temperature increases, the frequency branches of the tiny hourglass pockets evolve into a more cylindrical shape, while the larger pockets remain unchanged. This highlights that variations in exchange splitting, driven by changes in the magnetic moment, primarily impact the small Fermi pockets without significantly altering the overall band structure. As a result, this is consistent with first-principles calculations, which show minimal changes near the Fermi level across ferromagnetic and simple antiferromagnetic states or under varying on-site Coulomb repulsion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of pressure on the electronic and magnetic properties of bulk NiI 2

Transition metal dihalides have recently garnered interest in the context of two-dimensional van der Waals magnets as their underlying geometrically frustrated triangular lattice leads to interesting competing exchange interactions. In particular, NiI 2 is a magnetic semiconductor that has been long known for its exotic helimagnetism in the bulk. Recent experiments have shown that the helimagnetic state survives down to the monolayer limit with a layer-dependent magnetic transition temperature that suggests a relevant role of the interlayer coupling. Here, we explore the effects of hydrostatic pressure as a means to enhance this interlayer exchange and ultimately tune the electronic and magnetic response of NiI 2 . We study first the evolution of the structural parameters as a function of external pressure using first-principles calculations combined with x-ray diffraction measurements. We then examine the evolution of the electronic structure and magnetic exchange interactions via first-principles calculations and Monte Carlo simulations. Further, we find that the leading interlayer coupling is an antiferromagnetic second-nearest-neighbor interaction that increases monotonically with pressure. The ratio between isotropic third- and first-nearest-neighbor intralayer exchanges, which controls the magnetic frustration and determines the magnetic propagation vector q of the helimagnetic ground state, is also enhanced by pressure. As a consequence, our Monte Carlo simulations show a monotonic increase in the magnetic transition temperature, indicating that pressure is an effective means to tune the magnetic response of NiI 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Evolution of lithium ordering with (de)-lithiation in β-LiVOPO 4 : insights through solid-state NMR and first principles DFT calculations

The lithium-ion battery cathode material β-VOPO 4 is capable of intercalating more than one Li ion per transition metal ion due to the accessibility of both the V 5+ /V 4+ and V 4+ /V 3+ redox couples at ~4.5 V and ~2.3 V vs. Li, respectively, giving a theoretical capacity greater than ~300 mA h g -1 . The ability to perform full and reversible two Li-ion intercalation in this material, however, has been a matter of debate and the poor crystallinity of the fully lithiated phase has thus far precluded its complete structural characterisation by conventional diffraction-based methods. In this work, 7 Li and 31 P NMR spectroscopy, in combination with first principles DFT calculations, indicate that chemical lithiation results in a single phase β-Li 2 VOPO 4 exhibiting a complex Li ordering scheme with lithium ions occupying multiple disordered environments. 2D NMR 7 Li correlation experiments were used to deduce the most likely Li ordering for the β-Li 2 VOPO 4 phase from amongst several DFT optimised structures. In contrast, electrochemically lithiated β-Li 2- x VOPO 4 discharged to 1.6 V exhibits, in addition to β-Li 2 VOPO 4 , a β-Li 1.5 VOPO 4 phase. Additionally, the existence of β-Li 1.5 VOPO 4 is not reflected in the flat galvanostatic charge and discharge curves nor is evident from diffraction-based methods due to the very close structural similarity between the β-Li 1.5 VOPO 4 phase and β-Li 2 VOPO 4 phases. We demonstrate that solid state NMR spectroscopy, in combination with DFT results, provides a powerful tool for identifying intermediate states formed during charge/discharge of these complex phosphates as these phases can be distinguished from the end member phases primarily by the nature of the lithium ordering.

25 ENERGY STORAGE↗
DOE OSTI2021

Disorder and Itinerant Magnetism in Full Heusler Pd 2 TiIn

We report electronic and magnetic properties of full Heusler Pd 2 TiIn based on first principles calculations. This compound has been variously characterized as magnetic or non-magnetic. We use first principles calculations with accurate methods to reexamine this issue. We find that ideal ordered Heusler Pd 2 TiIn remains non-magnetic, in accord with prior work. Furthermore, we do find that it is possible to explain the magnetism seen in experiments through disorder and in particular we find that site disorder can lead to moment formation in this compound. In addition, we find an alternative low energy cubic crystal structure, which will be of interest to explore experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers

Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Electronic and magnetic properties of the topological semimetal SmMg 2 Bi 2

Dirac semimetals show nontrivial physical properties and can host exotic quantum states like Weyl semimetals and topological insulators under suitable external conditions. Here, by combining angle-resolved photoemission spectroscopy measurements (ARPES) and first-principle calculations, we demonstrate that the Zintl-phase compound SmMg 2 Bi 2 is in close proximity to a topological Dirac semimetallic state. ARPES results show a Dirac-like band crossing at the zone center near the Fermi level (E F ), which is further confirmed by first-principle calculations. Theoretical studies also reveal that SmMg 2 Bi 2 belongs to a Z 2 topological class and hosts spin-polarized states around the E F . Zintl's theory predicts that the valence state of Sm in this material should be Sm 2+ , however, we detect many Sm-4f multiplet states (flat-bands) whose energy positions and relative intensities suggest the presence of both dominant Sm 2+ and minor Sm 3+ . The small concentration (2.5%) of Sm 3+ in the bulk of a crystal is inferred to arise from Sm vacancies in the crystal. It is also evident that these flat bands and other dispersive states are strongly hybridized when they cross each other. Due to the presence of Sm 3+ ions, the temperature dependence of the magnetic susceptibility χ (T) shows a Curie-Weiss-like contribution in the low-temperature region, in addition to the Van Vleck-like behavior expected for the Sm 2+ ions. The present study will help to better understand the electronic structure, magnetism, and transport properties of related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Fermiology and transport properties of the candidate topological crystalline insulator SrAg 4 ⁢Sb 2

Compared to time-reversal symmetry-protected ℤ 2 topological insulators and Dirac/Weyl semimetals, there are significantly fewer candidates for topological crystalline insulators. SrAg 4 ⁢Sb 2 is predicted to exhibit topological crystalline insulator behavior when considering spin-orbit coupling. In this study, we systematically investigate single crystals of SrAg 4 ⁢Sb 2 using electrical transport and magnetic torque measurements, along with first-principles calculations. Our transport data reveals its compensated semimetal nature with a magnetoresistance up to around 700% at 2 K and 9 T. Analysis of de Haas–van Alphen oscillations uncovers a Fermi surface consisting of three distinct Fermi pockets with light effective masses. Comparison between the three-dimensional fermiology obtained from our oscillation data and the first-principles calculations demonstrates excellent agreement. This confirms the accuracy of the calculations, which indicate a band inversion centered at the 𝑇 point and identify the existence of nontrivial tube and needle hole Fermi pockets at Γ, alongside one trivial diamond electron pocket at the 𝐹 point in the Brillouin zone. Furthermore, symmetry and topology analysis results in two potential sets of topological invariants, suggesting the emergence of two-dimensional gapless Dirac surface states either on the 𝑎⁢𝑏 planes or on both the 𝑎⁢𝑏 planes and mirror planes, protected by crystal symmetries. Therefore, SrAg 4 ⁢Sb 2 emerges as a promising candidate topological crystalline insulator.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Is Ba 3 In 2 O 6 a high- T c superconductor?

It has been suggested that Ba 3 In 2 O 6 might be a high-Tc superconductor. Experimental investigation of the properties of Ba 3 In 2 O 6 was long inhibited by its instability in air. Recently epitaxial Ba 3 In 2 O 6 with a protective capping layer was demonstrated, which finally allows its electronic characterization. The optical bandgap of Ba3In2O6 is determined to be 2.99 eV in-the (001) plane and 2.83 eV along the c-axis direction by spectroscopic ellipsometry. First-principles calculations were carried out, yielding a result in good agreement with the experimental value. Various dopants were explored to induce (super-)conductivity in this otherwise insulating material. Neither A- nor B-site doping proved successful. The underlying reason is predominately the formation of oxygen interstitials as revealed by scanning transmission electron microscopy and first-principles calculations. Additional efforts to induce superconductivity were investigated, including surface alkali doping, optical pumping, and hydrogen reduction. To probe liquid-ion gating, Ba 3 In 2 O 6 was successfully grown epitaxially on an epitaxial SrRuO 3 bottom electrode. So far none of these efforts induced superconductivity in Ba 3 In 2 O 6 , leaving the answer to the initial question of whether Ba 3 In 2 O 6 is a high-T c superconductor to be 'no' thus far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Topological Circular Dichroism in Chiral Multifold Semimetals

Uncovering the physical contents of the nontrivial topology of quantum states is a critical problem in condensed matter physics. Here, we study the topological circular dichroism in chiral semimetals using linear response theory and first-principles calculations. We show that, when the low-energy spectrum respects emergent SO(3) rotational symmetry, topological circular dichroism is forbidden for Weyl fermions, and thus is unique to chiral multifold fermions. This is a result of the selection rule that is imposed by the emergent symmetry under the combination of particle-hole conjugation and spatial inversion. Using first-principles calculations, we predict that topological circular dichroism occurs in CoSi for photon energy below about 0.2 eV. In conclusion, our Letter demonstrates the existence of a response property of unconventional fermions that is fundamentally different from the response of Dirac and Weyl fermions, motivating further study to uncover other unique responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Distinguishing Thermal Fluctuations from Polaron Formation in Halide Perovskites

Recent angle-resolved photoelectron spectroscopy (ARPES) measurements of the hole effective mass in CsPbBr$_3$ revealed an enhancement of $\sim$50 % compared to the bare mass computed from first principles for CsPbBr$_3$ at $T = 0 K$. This large enhancement was interpreted as evidence of polaron formation. Employing accurate finite-temperature first-principles calculations, we show that the calculated hole effective mass of CsPbBr$_3$ at $T = 300 K$ can explain experimental results without invoking polarons. Thermal fluctuations are particularly strong in halide perovskites compared to conventional semiconductors such as Si and GaAs, and cannot be ignored when comparing with experiment. We not only resolve the debate on polaron formation in halide perovskites, but also demonstrate the general importance of including thermal fluctuations in first-principles calculations for strongly anharmonic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Population Balance Models for Catalytic Depolymerization: From Elementary Steps to Multiphase Reactors

Here, the ongoing accumulation of plastic waste in landfills and in the environment is driving research on chemical processes and catalysts to recycle polymers. Traditional modeling strategies are not applicable to these processes because they involve too many reactants and intermediates, one for each molecular weight and each functionalization. To model the kinetics, we have developed population balance models (PBMs) that account for macromolecular reactants in the bulk and macromolecular catalytic intermediates. These PBMs couple to each other through polymer adsorption and desorption models and to traditional rate equations for small molecule products and co-reactants (like hydrogen or ethylene). The models, in combination with experimental data, are being used in many ways: (i) to test mechanistic hypotheses, (ii) to extract rate parameters, (iii) to quantitatively compare catalyst activities, (iv) to account for mass transfer and vapor–liquid partitioning in two-phase reactors, and (v) to design novel support architectures and catalysts that mimic the processive action of natural depolymerization enzymes. Some key theoretical advances allow PBMs to be constructed from elementary rates and mechanisms, as opposed to traditional formulations with pseudoelementary rate parameters invoked as fitting parameters. We discuss ways to build these models “bottom-up” from first-principles calculations and ways to extract model parameters from “top down” analyses of rate data. The combination provides a quantitative bridge between first-principles calculations and the kinetics of complex macromolecular transformations for polymer upcycling and beyond.

Manis, Lela K. [University of Illinois at Urbana-C↗