Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

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

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

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vacancy surface migration mechanisms in dilute nickel-chromium alloys

We investigate the unitary mechanisms related to the surface migration of vacancies in dilute Ni-Cr alloys via first-principle calculations. We survey a complete set of surface and sub-surface migration paths for vacancies near the (100) free surface and calculate the corresponding migration barriers. Our results show that a vacancy migrating towards the free surface will face lower energy barriers to migrate via atomic exchange with a neighboring Cr atom rather than with a Ni atom. Once a vacancy reaches the free surface, it will be trapped there. Our results also reveal that, when a Cr atom sits in the atomic plane just below the free surface, any in-plane vacancy hopping jump that would result in the vacancy sitting on top a subsurface Cr atom is energetically unfavorable. Taken together, these fundamental unitary surface migration mechanisms offer insights into the complex interactions between surface segregation and vacancy migration phenomena in Ni-Cr-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Topological antiferromagnetic semimetal for spintronics: A case study of a layered square-net system EuZnSb 2

Here, we use the first-principles and effective Hamiltonian methods to study the electronic structure and magnetic properties of a recently synthesized layered antiferromagnetic square-net topological semimetal EuZnSb2 [Phys. Rev. Res. 2, 033462 (2020)]. The main message of the paper is that effects of small changes in the band structure produced by the magnetic ordering and changes in the orientation of the Néel vector are amplified in such transport properties as the spin Hall conductivity. We predict that the effects of the broken symmetry introduced by the ordering of the Néel vector, being very weak in the bulk, are pronounced in the surface electronic dispersion, suggesting that surface probes may be more suited to measure them. The coexistence of the magnetism with many other competing phases make this material interesting and possibly useful for quantum spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Li interaction-induced phase transition from black to blue phosphorene

A comprehensive first-principle calculation has been carried out and revealed that at sufficiently high Li concentration and certain well-defined configurations a phase transition from black to blue phosphorene can take place. Blue phosphorene, a newly predicted allotrope of phosphorus, possesses unique crystalline and electronic structure and is a promising candidate, not only for fundamental research but also for electronic and optoelectronic applications. Methods to growth high quality blue phosphorene layers are highly desirable but challenging. Here, a novel kinetic pathway to grow blue phosphorene layers from black phosphorene layers via Li intercalation is proposed based on first principle study. This study pointed out that Li atoms intercalated in black phosphorene could act as ‘catalyst’ in the ‘reactive region’ of the lone pair of P atoms, leading to a P-P bond breaking and subsequently, a local structural transformation from orthorhombic lattice to an assembly of parallel narrow nanoribbons with rhombohedra-like symmetry. During Li deintercalation, these nanoribbons are self-mended and form blue phosphorene layers. The interlayer distance was found 4.60 Å for double layer with AA stacking and 4.13 Å for multilayer with ABC stacking, respectively, indicating a monolayer blue phosphorene can be mechanically exfoliated. Furthermore, this study also points out the possibility of new phases in other systems, where intercalation can lead to an unexpected structural phase transition and even a discovery of novel materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Deep-Level Defects and Impurities in InGaN Alloys

In this study, density functional theory calculations with a hybrid functional are used to examine the charge-state transition levels of native point defects and impurities in InGaN alloys, with the goal of identifying centers that play a role in defect-assisted recombination. Explicit alloy calculations are used to monitor the dependence of defect levels on indium content and distribution of In atoms. The relative shift (or lack thereof) of the charge-state transition levels of the different defects is explained by the atomic character of the defect state and whether it is derived from valence-band or conduction-band states of the host material or acts as an atomic-like impurity. The various possible atomic configurations of In and Ga cations for a given composition of InGaN lead to a distribution of charge-state transition levels. Defects on the nitrogen site lead to a larger spread in levels compared with defects on the cation site.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Machine learning guided discovery of superconducting calcium borocarbides

Pursuit of superconductivity in light-element systems at ambient pressure is of great experimental and theoretical interest. In this work, we combine a machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary Ca-B-C compounds. Three layered borocarbides (stable CaBC 5 and metastable Ca 2 BC 11 and CaB 3 C 3 ) are predicted to be phonon-mediated superconductors at ambient pressure. The stable CaBC 5 and the low-energy metastable Ca 2 BC 11 (with formation energy only 9.5 meV/atom above the convex hull) have a superconducting T c of 5.2 and 8.9 K, respectively. While the hexagonal CaB 3 C 3 possesses a T c of 26.1 K, it is metastable with formation energy of 153 meV/atom above the convex hull. The ML-guided approach opens up a way for greatly accelerating the discovery of new high-T c superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Split Ga vacancies and the unusually strong anisotropy of positron annihilation spectra in β–Ga 2 O 3

We report a systematic first-principles study on positron annihilation parameters in the β–Ga 2 O 3 lattice and Ga monovacancy defects complemented with orientation-dependent experiments of the Doppler broadening of the positron-electron annihilation. We find that both the β–Ga 2 O 3 lattice and the considered defects exhibit unusually strong anisotropy in their Doppler broadening signals. This anisotropy is associated with low symmetry of the β–Ga 2 O 3 crystal structure that leads to unusual kind of one-dimensional confinement of positrons even in the delocalized state in the lattice. In particular, the split Ga vacancies recently observed by scanning transmission electron microscopy produce unusually anisotropic positron annihilation signals. Furthermore, we show that in experiments, the positron annihilation signals in β–Ga 2 O 3 samples seem to be often dominated by split Ga vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Equation of state, phase transitions, and band-gap closure in Pb Cl 2 and Sn Cl 2

The equations of state and band-gap closures for Pb Cl 2 and Sn Cl 2 were studied using both experimental and theoretical methods. Here we measured the volume of both materials to a maximum pressure of 70 GPa using synchrotron-based angle-dispersive powder x-ray diffraction. The lattice parameters for both compounds showed anomalous changes between 16–32 GPa, providing evidence of a phase transition from the cotunnite structure to the related Co 2 Si structure, in contrast to the postcotunnite structure as previously suggested. First-principles calculations confirm this finding and predict a second phase transition to a Co 2 Si-like structure between 75– 110 GPa in Pb Cl 2 and 60–75 GPa in Sn Cl 2 . Band gaps were measured under compression to ~70 GPa for Pb Cl 2 and ~66 GPa for Sn Cl 2 and calculated up to 200 GPa for Pb CI 2 and 120 GPa for Sn Cl 2 . We find an excellent agreement between our experimental and theoretical results when using the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional, which suggests that this functional could reliably be used to calculate the band gap of similar AX 2 compounds. Experimental and calculated band-gap results show discontinuous decreases in the band gap corresponding to phase changes to higher-coordinated crystal structures, giving insight into the relationship between interatomic geometry and metallicity.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Magnetoelectric control of topological phases in graphene

Topological antiferromagnetic (AFM) spintronics is an emerging field of research, which involves the topological electronic states coupled to the AFM order parameter known as the Néel vector. The control of these states is envisioned through manipulation of the Néel vector by spin-orbit torques driven by electric currents. In this work, we present a different approach favorable for low-power AFM spintronics, where the control of the topological states in a two-dimensional material, such as graphene, is performed via the proximity effect by the voltage induced switching of the Néel vector in an adjacent magnetoelectric AFM insulator, such as chromia. Mediated by the symmetry protected boundary magnetization and the induced Rashba-type spin-orbit coupling at the interface between graphene and chromia, the emergent topological phases in graphene can be controlled by the Néel vector. Using density functional theory and tight-binding Hamiltonian approaches, we model a graphene / Cr 2 O 3 (0001) interface and demonstrate nontrivial band gap openings in the graphene Dirac bands asymmetric between the K and K ' valleys. This gives rise to an unconventional quantum anomalous Hall effect (QAHE) with a quantized value of 2 e 2 / h and an additional steplike feature at a value close to e 2 / 2 h , and the emergence of the spin-polarized valley Hall effect (VHE). Additionally, depending on the Néel vector orientation, we predict the appearance and transformation of different topological phases in graphene across the 180° AFM domain wall, involving the QAHE, the valley-polarized QAHE, and the quantum VHE, and the emergence of the chiral edge states along the domain wall. These topological properties are controlled by voltage through magnetoelectric switching of the AFM insulator with no need for spin-orbit torques.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrogen Dissociation Controls 1-Hexyne Selective Hydrogenation on Dilute Pd-in-Au Catalysts

Increasing the selectivity of the catalytic hydrogenation of alkynes to alkenes is of major importance for the processing of petrochemicals and the production of fine chemicals. Achieving high selectivity for alkene formation at high conversions, however, remains a long-standing challenge in heterogeneous catalysis. Here, the mechanism and origin of the high selectivity of dilute Pd-in-Au catalysts has been studied by a combination of first-principles calculations, microkinetic simulations, and isotopic exchange hydrogenation experiments using Pd 0.04 Au 0.96 nanoparticles embedded in raspberry colloid-templated silica. The Pd is predominantly in the form of isolated atoms, only surrounded by Au atoms, based on prior studies. The simulations indicate that the rate-limiting process for 1-hexyne hydrogenation on Pd monomers in Au(111) is H 2 dissociation, which has a large free energy barrier of 0.86 eV at 363 K and 0.2 bar of H 2 . The C–H bond formation steps, on the other hand, proceed with lower barriers, which contrasts with previous studies of extended Pd catalysts. The microkinetic simulations identify the sizable H 2 dissociation barrier and the small barrier for the hydrogenation of 1-hexyne as key factors that lead to a high selectivity for the production of 1-hexene from 1-hexyne, even at high conversion. The unconventional H 2 dissociation limiting process in combination with the low coverage of weakly bound hydrocarbon intermediates explains the near-zero order of 1-hexyne found experimentally. Furthermore, the partial hydrogenation of 1-hexyne to form 1-hexene is shown to be an irreversible process from our isotopic exchange hydrogenation experiments and is explained by the strongly exothermic nature of the reaction. Diluting active species, like Pd, in a less active host metal, like Au, hence appears promising as a means of tuning the binding energy of reactants and altering reaction profiles, leading to distinct kinetic behavior for an optimal catalytic activity and selectivity. Furthermore, the combination of microkinetic modeling, density functional theory calculations, and isotopic exchange experiments is thus demonstrated to be an effective approach to modeling important catalytic phenomena.

1-hexyne↗
DOE OSTI2025

HfNb_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-niobium (Hf-Nb). 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 Nb. 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,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetoelastic coupling and spin contributions to entropy and thermal transport in biferroic yttrium orthochromite

Direct engineering of material properties through exploitation of spin, phonon, and charge-coupled degrees of freedom is an active area of development in materials science. However, the relative contribution of the competing orders to controlling the desired behavior is challenging to decipher. In particular, the independent role of phonons, magnons, and electrons, quasiparticle coupling, and relative contributions to the phase transition free energy largely remain unexplored, especially for magnetic phase transitions. Here, we study the lattice and magnetic dynamics of biferroic yttrium orthochromite using Raman, infrared, and inelastic neutron spectroscopy techniques, supporting our experimental results with first-principles lattice dynamics and spin-wave simulations across the antiferromagnetic transition at T N ~ 138 K. Spectroscopy data and simulations together with the heat capacity (C p ) measurements, allow us to quantify individual entropic contributions from phonons (0.01 ± 0.01k B atom -1 ), dilational (0.03 ± 0.01k B atom -1 ), and magnons (0.11 ± 0.01k B atom -1 ) across TN. High-resolution phonon measurements conducted in a magnetic field show that anomalous T-dependence of phonon energies across TN originates from magnetoelastic coupling. Additionally, phonon scattering is primarily governed by the phonon–phonon coupling, with little contribution from magnon–phonon coupling, short-range spin correlations, or magnetostriction effects; a conclusion further supported by our thermal conductivity measurements conducted up to 14 T, and phenomenological modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local positional and spin symmetry breaking as a source of magnetism and insulation in paramagnetic EuTi O 3

Here we consider theoretically the paramagnetic phases of EuTiO 3 that represent configurations created by two sets of microscopic degrees of freedom (m-DOFs): positional symmetry breaking due to octahedral rotations and magnetic symmetry breaking due to spin disorder. The effect of these sets of m-DOFs on the electronic structure and properties of the para phases is assessed by considering sufficiently large (super) cells with the required nominal global average symmetry, allowing, however, the $\mathcal{local}$ positional and magnetic symmetries to be lowered. We find that tendencies for local symmetry breaking can be monitored by following total energy lowering in mean-fieldlike density-functional theory, without recourse for strong correlation effects. While most nominally cubic ABO 3 perovskites are known for their symmetry breaking due to the B-atom sublattice, the case of f-electron magnetism in EuTiO 3 is associated with A-sublattice symmetry breaking and its coupling to structural distortions. We find that (i) paramagnetic $\mathcal{cubic}$ EuTiO 3 has an intrinsic tendency for both magnetic and positional symmetry breaking, while paramagnetic $\mathcal{tetragonal}$ EuTiO 3 has only magnetic symmetry lowering and no noticeable positional symmetry lowering with respect to low-temperature antiferromagnetic tetragonal phase. (ii) Properly modeled paramagnetic tetragonal and cubic EuTiO 3 have a nonzero local magnetic moment on each Eu ion, consistent with the experimental observations of local magnetism in the para phases of EuTiO 3 significantly above the Néel temperature. Interestingly, (iii) the local positional distortion modes in the short-range ordered para phases are inherited from the long-range ordered low-temperature antiferromagnetic ground-state phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting the Activity and Selectivity of Bimetallic Metal Catalysts for Ethanol Reforming using Machine Learning

Machine learning is ideally suited for the pattern detection in large uniform datasets, but consistent experimental datasets on catalyst studies are often small. Here we demonstrate how a combination of machine learning and first-principles calculations can be used to extract knowledge from a relatively small set of experimental data. The approach is based on combining a complex machine-learning model trained on an extensive computational library of transition-state energies with simple linear regression models of experimental catalytic activities and selectivities from the literature. Using the combined model, we identify the key C–C bond scission reactions involved in ethanol reforming and perform a computational screening for ethanol reforming on monolayer bimetallic catalysts with architectures TM-Pt-Pt(111) and Pt-TM-Pt(111) (TM = 3d transition metals). The model also predicts four promising catalyst compositions for future experimental studies. In conclusion, the approach is not limited to ethanol reforming but is of general use for the interpretation of experimental observations as well as for the computational discovery of novel catalytic materials.

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

First-principles thermodynamic assessment of Sr-containing secondary phase formation in strontium-substituted lanthanum manganites for solid oxide cell applications

Sr-secondary phase formation is a potentially significant degradation mode with direct impact upon solid-oxide cell (SOC) commercial viability. A first-principles based thermodynamic study was performed for La 1−x Sr x MnO 3±δ (LSM) perovskites to assess their stability against formation of different Sr-secondary phases, including SrO, SrCrO 4 , SrSO 4 , SrCO 3 , and Sr(OH) 2 , for SOC applications. The Sr-secondary phase formation reaction free energies were determined via a thermodynamic model by combining ab initio lattice dynamics calculations for the solid phases and ab initio thermodynamic data for the gas phases. The current approach expands the previously reported thermodynamic modeling studies by integrating first-principles based point defect equilibria into the thermodynamic analysis. The modeling results obtained using this new approach indicate an increased tendency to form the SrO oxide upon decreasing the oxygen partial pressure. Additionally, the enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with the contaminant gas species as a function of temperature and pressure.

defect thermodynamics modeling↗
DOE OSTI2023

Nanoscale imaging of He-ion irradiation effects on amorphous TaO x toward electroforming-free neuromorphic functions

Resistive switching in thin films has been widely studied in a broad range of materials. Yet, the mechanisms behind electroresistive switching have been persistently difficult to decipher and control, in part due to their non-equilibrium nature. Here, we demonstrate new experimental approaches that can probe resistive switching phenomena, utilizing amorphous TaO x as a model material system. Specifically, we applied scanning microwave impedance microscopy and cathodoluminescence (CL) microscopy as direct probes of conductance and electronic structure, respectively. These methods provide direct evidence of the electronic state of TaO x despite its amorphous nature. For example, CL identifies characteristic impurity levels in TaO x , in agreement with first principles calculations. We applied these methods to investigate He-ion-beam irradiation as a path to activate conductivity of materials and enable electroforming-free control over resistive switching. However, we find that even though He-ions begin to modify the nature of bonds even at the lowest doses, the films' conductive properties exhibit remarkable stability with large displacement damage and they are driven to metallic states only at the limit of structural decomposition. Finally, we show that electroforming in a nanoscale junction can be carried out with a dissipated power of <20 nW, a much smaller value compared to earlier studies and one that minimizes irreversible structural modifications of the films. Importantly, the multimodal approach described here provides a new framework toward the theory/experiment guided design and optimization of electroresistive materials.

36 MATERIALS SCIENCE↗