Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A first-principles study on new high-pressure metastable polymorphs of MoO{sub 2}

The pressure-dependence of the stabilities of several MoO{sub 2} phases has been investigated by density-functional theory (GGA/PBE/PAW). Out of a set of 15 MX{sub 2} structures, the [SnO{sub 2}(II)], [α-PbO{sub 2}], and a modified rutile structure type were identified as possible metastable MoO{sub 2} polymorphs based on the analysis of thermodynamic properties and dynamic stability. High-pressure calculations suggest an orthorhombic TiO{sub 2} structure, dubbed [ortho-TiO{sub 2}], as a high-pressure polymorph at around 25 GPa. Furthermore, we find that the previously reported rutile-type MoO{sub 2} may be understood as a modified rutile type similar to the [VO{sub 2}] structure. - Graphical abstract: First-principles electronic structure and thermochemical calculations reveal four structure candidates as possible metastable polymorphs of MoO{sub 2}. Most promising is a distorted rutile-type similar to the known [VO{sub 2}] structure. An orthorhombic polymorph is proposed as a high-pressure polymorph. Display Omitted - Highlights: • Three possible metastable structure candidates for MoO{sub 2}. • Undistorted rutile type is improbable, a new distorted rutile-type MoO{sub 2} was suggested. • Orthorhombic phase of MoO{sub 2} (ortho-TiO{sub 2} type) should form at 25 GPa. • ab initio thermochemical data provided for MoO{sub 2}.

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

First-principles microkinetic study of methane and hydrogen sulfide catalytic conversion to methanethiol/dimethyl sulfide on Mo 6 S 8 clusters: activity/selectivity of different promoters

A large fraction of the global natural gas reserves is in the form of sour gas, i.e. contains hydrogen sulfide (H 2 S) and carbon dioxide (CO 2 ), and needs to be sweetened before utilization. The traditional amine-based separation process is energy-intensive, thereby lowering the value of the sour gas. Thus, there is a need to find alternative processes to remove, e.g. , hydrogen sulfide. Mo 6 S 8 clusters are promising candidates for transforming methane (CH 4 ) and hydrogen sulfide into methanethiol (CH 3 SH) and dimethyl sulfide (CH 3 SCH 3 ), which are high-value sulfur-containing products that can be further used in the chemical industry. Here first-principles microkinetics is used to investigate the activity and selectivity of bare and promoted (K, Ni, Cl) Mo 6 S 8 . The results show that methanethiol is produced via two different pathways (direct and stepwise), while dimethyl sulfide is formed via a competing pathway in the stepwise formation of methanethiol. Moreover, there is an increase in activity and a decrease in selectivity when adding an electropositive promoter (K), whereas the reverse behaviour is observed when adding an electronegative promoter (Cl). When adding Ni there is also a decrease in activity and an increase in selectivity; however, Ni is acting as an electron donor. The results provide insights and guidance as to what catalyst formulation is preferred for the removal of hydrogen sulfide in sour gas.

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

Degradation of Alkali-Based Photocathodes from Exposure to Residual Gases: A First-Principles Study

Photocathodes are a key component in the production of electron beams in systems such as X-ray free-electron lasers and X-ray energy-recovery linacs. Alkali-based materials display high quantum efficiency (QE), however, their QE undergoes degradation faster than metal photocathodes even in the high vacuum conditions where they operate. The high reactivity of alkali-based surfaces points to surface reactions with residual gases as one of the most important factors for the degradation of QE. In order to advance the understanding on the degradation of the QE, we investigated the surface reactivity of common residual gas molecules (e.g., O 2 , CO 2 , CO, H 2 O, N 2 , and H 2 ) on one of the best-known alkali-based photocathode materials, cesium antimonide (Cs 3 Sb), using first-principles calculations based on density functional theory. Furthermore, the reaction sites, adsorption energy, and effect in the local electronic structure upon reaction of these molecules on (001), (110), and (111) surfaces of Cs 3 Sb were computed and analyzed. The adsorption energy of these molecules on Cs3Sb follows the trend of O 2 (-4.5 eV) > CO 2 (-1.9 eV) > H 2 O (-1.0 eV) > CO (-0.8 eV) > N 2 (-0.3 eV) ≈ H 2 (-0.2 eV), which agrees with experimental data on the effect of these gases on the degradation of QE. The interaction strength is determined by the charge transfer from the surfaces to the molecules. The adsorption and dissociation of O containing molecules modify the surface chemistry such as the composition, structure, charge distribution, surface dipole, and work function of Cs 3 Sb, resulting in the degradation of QE with exposure to O 2 , CO 2 , H 2 O, and CO.

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

High-accuracy first-principles-based rate coefficients for the reaction of OH and CH 3 OOH

Here, the ˙OH-initiated oxidation of methyl hydroperoxide, which plays an important role in the atmospheric chemistry of methane, was theoretically characterized using high-accuracy composite amHEAT-345(Q) coupled-cluster calculations followed by a two-dimensional E,J resolved master equation analysis. The reaction is found to proceed through two distinct hydrogen-bonded pre-reactive complexes leading to two product channels, in accord with the experimental observations: (i) ˙OH + CH 3 OOH → CH 3 OO˙ + H 2 O with a yield of 0.8 ± 0.1, and (ii) ˙OH + CH 3 OOH → HCHO + ˙OH + H 2 O with a yield of 0.2 ± 0.1. The calculated reaction enthalpies are within 0.2 kcal mol -1 of the benchmark ATcT values. Overall thermal rate coefficients obtained from first principles are found to be in the low-pressure limit at atmospheric pressure; the total rate coefficient can be expressed over the T = 200–450 K range as k(T) = 5.0 × 10 -12 × T -0.152 × exp(287/T) cm 3 s -1 , strongly supporting the experimental results of Vaghjiani and Ravishankara, with which this expression agrees within ca. 15%. The current results show that (i) is the principal reaction channel and support the view that, due to its inherently fast transformations, CH 3 OOH is an important redistribution species for HO x ˙ radicals in the Earth's atmosphere.

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

Low-index surface energies, cleavage energies, and surface relaxations for crystalline NiAl from first-principles calculations

NiAl surfaces frequently serve as a platform for studying a broad range of physical and chemical phenomena including chemisorption, catalysis, oxidation, alloy growth, and surface nanostructure formation. Knowledge of precise values for low-index surface energies of NiAl, the most fundamental quantities characterizing surface thermodynamics, is often invaluable for understanding of these phenomena. In all previous analyses for NiAl(100) and NiAl(111), Ni- and Al-terminations are not distinguished, and half of the cleavage energy (or equivalently, the average of the surface energies of two differently terminated surfaces) is always identified as the “surface energy”. No values are available for individual surface energies of Ni- or Al-terminated NiAl(100) or NiAl(111) surfaces, whereas knowledge of only cleavage energy is often insufficient for analyzing surface-associated behavior. As such, in this work we perform extensive first-principles density-functional-theory (DFT) calculations for surface energies and cleavage energies of NiAl(110), NiAl(100) and NiAl(111) by considering the chemical-potential-based formulations to clarify the ambiguity in their surface energies and cleavage energies. We obtain a surface-energy phase diagram for these three low-index surfaces versus the relevant chemical potential, as well as the chemical-potential-dependent Wulff plots for NiAl crystal equilibrium shapes. We also provide the surface-relaxation information from our DFT calculations for comparison with previous experimental data.

42 ENGINEERING↗
DOE OSTI2017

First-Principles Simulation of the (Li–Ni–Vacancy)O Phase Diagram and Its Relevance for the Surface Phases in Ni-Rich Li-Ion Cathode Materials

Despite several reports on the surface phase transformations from a layered to a disordered spinel and a rock-salt structure at the surface of the Ni-rich cathodes, the precise structures and compositions of these surface phases are unknown. The phenomenon, in itself, is complex and involves the participation of several contributing factors. Of these factors, transition metal (TM) ion migration toward the interior of the particle and hence formation of TM-densified surface layers, triggered by oxygen loss, is thermodynamically probable. Here in this study, we simulate the thermodynamic phase equilibria as a function of TM ion content in the cathode material in the context of lithium nickel oxides, using a combined approach of first-principles density functional calculations, the cluster expansion method, and grand canonical Monte Carlo simulations. We developed a unified lattice Hamiltonian that accommodates not only rock-salt like structures but also topologically different spinel-like structures. Also, our model provides a foundation to investigate metastable cation compositions and kinetics of the phase transformations. Our investigations predict the existence of several Ni-rich phases that were, to date, unknown in the scientific literature. Our simulated phase diagrams at finite temperature show a very low solubility range of the prototype spinel phase. Finally, we find a partially disordered spinel-like phase with far greater solubility that is expected to show very different Li diffusivity compared to that of the prototype spinel structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Combined First-Principles and Experimental Investigation into the Reactivity of Codeposited Chromium–Carbon under Pressure

High-pressure synthesis in the diamond anvil cell suffers from the lack of a general approach for the control of precursor stoichiometry and homogeneity. Here, we present results from a new method we have developed that uses magnetron cosputtering to prepare stoichiometrically precise and atomically mixed amorphous films of Cr:C. Laser-heated diamond anvil cell experiments carried out on a flake of this sample at pressures between 13.5 and 24.3 GPa lead to the observation of Cr 3 C (Pnma) over the entire pressure range–in good agreement with our in-house theoretical predictions–but also reveal two other metastable phases that were not expected: a novel monoclinic chromium carbide phase and the NaCl-type CrC (Fm3̅m) phase. The unexpected stability of CrC is investigated by using first-principles methods, revealing a large stabilizing effect tied to substoichiometry at the carbon site. These results offer an important case study into the current limitations of crystal structure prediction methods with regard to phase complexity and bolster the growing need for advanced theoretical approaches that can more completely survey experimentally unexplored phase space.

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

A method of orbital analysis for large-scale first-principles simulations

An efficient method of calculating the natural bond orbitals (NBOs) based on a truncation of the entire density matrix of a whole system is presented for large-scale density functional theory calculations. The method recovers an orbital picture for O(N) electronic structure methods which directly evaluate the density matrix without using Kohn-Sham orbitals, thus enabling quantitative analysis of chemical reactions in large-scale systems in the language of localized Lewis-type chemical bonds. With the density matrix calculated by either an exact diagonalization or O(N) method, the computational cost is O(1) for the calculation of NBOs associated with a local region where a chemical reaction takes place. As an illustration of the method, we demonstrate how an electronic structure in a local region of interest can be analyzed by NBOs in a large-scale first-principles molecular dynamics simulation for a liquid electrolyte bulk model (propylene carbonate + LiBF{sub 4})

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

Structure, Magnetism, and First-Principles Modeling of the Na 0.5 La 0.5 RuO 3 Perovskite

High-purity polycrystalline Na 0.5 La 0.5 RuO 3 was synthesized by a solid-state method, and its properties were studied by magnetic susceptibility, heat capacity, and resistivity measurements. In this study, we find it to be an orthorhombic perovskite, in contrast to an earlier report, with random La/Na mixing. With a Curie–Weiss temperature of –231 K and an effective moment of 2.74 μB/mol Ru, there is no magnetic ordering down to 1.8 K. A broad hump at 1.4 K in the heat capacity, however, indicates the presence of a glassy magnetic transition, which we attribute to the influence of the random distribution of Na and La on the perovskite A-sites. Comparison to CaRuO 3 , a structurally ordered ruthenate perovskite with a similar A-site ionic radius and magnetic properties, is presented. First-principles calculations indicate that the Na–La distribution determines the local magnetic exchange interactions between Ru ions, favoring either antiferromagnetic or ferromagnetic coupling when the local environment is Na- or La-rich. Thus, our data and analysis suggest that mixing cations with different charges and sizes on the A-site in this perovskite results in magnetic frustration through a balance of local magnetic exchange interactions.

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

Metastable Li 1+δ Mn 2 O 4 (0 ≤ δ ≤ 1) Spinel Phases Revealed by in Operando Neutron Diffraction and First-Principles Calculations

Here, we examine the reaction mechanisms driven by the lithiation of Li 1+δ Mn 2 O 4 (0 ≤ δ ≤ 1) spinels via in operando neutron powder diffraction (NPD). New reaction mechanisms are proposed involving solid-solution regions within both cubic and tetragonal spinel phases in addition to a continuous phase transition between them. In operando NPD is an ideal tool to follow the light elements such as lithium and oxygen in cathode materials which are often the key to fully understand their structural evolutions. Here, we report a novel methodology to prepare an extremely thick electrode with ~378 mg·cm –2 loading density suitable for the in operando NPD studies. Enabled by such thick electrode, we find that the metal oxygen M–O (M = Li and Mn) bond lengths in both end members LiMn 2 O 4 and Li 2 Mn 2 O 4 experience pronounced changes larger than dictated by the change in lattice parameters because of the locally formed Jahn–Teller distorted Mn 3+ . First-principles density functional theory calculations confirm these metastable intermediates and further propose atomistic reaction pathways for the phase transition by coupling a global structure search algorithm. In conclusion, these findings redefine the conventional understandings on two-phase reactions of this spinel material.

36 MATERIALS SCIENCE↗
DOE OSTI2016

First-principles study of inert gas incorporation and migration in zirconium nitride

Uranium-Molybdenum alloy based fuels are the most promising fuel for the future high performance research reactors using low enriched uranium (LEU). However, when used in dispersion with aluminum, significant amount of interaction products form at the interface between the UMo fuel particle and the surrounding Al during in-pile tests which leads to non-acceptable swelling. To reduce the swelling in U-Mo/Al dispersion fuel at high burnups, different diffusion barrier coatings, e.g., ZrN, were deposited on the U-Mo particles to avoid the interaction between fuel particles and the Al matrix. These diffusion barrier coatings are very effective in reducing the formation of interaction layer and therefore the fuel swelling. So far, ZrN is selected as the candidate coating material for the European Dispersion Fuel program due to the superior thermophysical properties of ZrN. The thermodynamic and kinetic stability of ZrN coating with respect to Al and U-Mo has been studied by the current authors using density functional theory calculations. So far there are very few studies of the diffusion behaviors of inert gases in ZrN. A recent publication on irradiated (Pu,Zr)N reveals that the majority of fission-induced Xe is retained in the material, while He is release to the fuel pin to a large extent. This raises the questions on the underlying differences of inert gas behavior in nitride fuel. Therefore, the purpose of this study is to understand the incorporation and migration behavior of inert gases He, Kr and Xe in ZrN using density functional theory (DFT) based first-principles calculations. Using density-functional theory calculations we provided a systematic study of the incorporation and migration of fission gases in ZrN. We investigated the formation energies of several defects (vacancies, interstitials, divacancies, Frenkel pairs, and Schottky defects) and the incorporation energies of inert gases in these defects. Results show that the formation energies of point defects are strongly affected by the chemical potentials of N and Zr species. The energetically most favorable incorporation site for fission gases is found to be the first nearest neighbor Schottky defect, in which the fission gas atom sits at somewhere between the Zr and N vacancies. The predicted solution energies of rare gases, i.e, He, Kr and Xe, in ZrN, is almost linearly dependent on the atom size of inert gases. We compared the solution energy Xe in ZrN with several other nuclear materials with the same B1 structure (MX), e.g., TiN, UN and UC. The solution energies of Xe in MX are closely related to the size of vacancy sites. The potential migration mechanisms of Xe in ZrN include interstitial and vacancy-assisted diffusion. The calculations show that vacancies essentially act as Xe trapping sites. The migration behavior of Xe in N and Zr divacancy and Schottky defects were compared with that through tetrahedral interstitials. It is found that the migration barrier is much lower with the interstitial mechanism than those with the vacancy mechanism. Therefore, at low temperatures Xe diffusion in ZrN is mostly through the interstitial mechanism or existing grain boundaries while at high temperature volume diffusion could be preferred, i.e., diffusion with Zr and N vacancies. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2021

A first principles framework to predict the transient performance of latent heat thermal energy storage

Thermal energy storage (TES) is increasingly recognized as an essential component of efficient Combined Heat and Power (CHP), Concentrated Solar Power (CSP), Heating Ventilation and Air Conditioning (HVAC), and refrigeration as it reduces peak demand while helping to manage intermittent availability of energy (e.g., from solar or wind). Latent Heat Thermal Energy Storage (LHTES) is a viable option because of its high energy storage density. Parametric analysis of LHTES in terms of dimensionless numbers is highly desired as a tool to model LHTES systems. One approach is to develop a model equation so as to minimize the error between the model and data obtained from experiments or simulations. While this approach can produce an accurate correlation applicable within the range of data used for its creation, it does not provide physical understanding of the rate-limiting process controlling the transient behavior of the device. In this paper we present an alternative approach whereby the potential rate-limiting processes are identified from first principles and then the key process is determined as a function of time as a LHTES device is charged. For example, in a simple geometry, the melt-fraction can be expected to vary linearly in time if the heat transfer rate is limited by natural convection of the phase changing material and we show it scales with the PCM Grashof number as $Gr^1_p$ and PCM Prandtl number as $Pr_p^{(1/3)}$. On the other hand, if surface area of solid PCM limits the heat transfer rate, the melt fraction increases asymptotically to reach full melting. The existence of these linear and asymptotic regions and the $Gr^1_pP r^{1/3}_p$ shape of the melt fraction curve is verified using our database of 64 simulations. Of practical importance in designing LHTES devices is the melt fraction at which the heat transfer rate ceases to be limited by convection, after which the heat storage rate deteriorates. For our geometry, this is found to be about 90%. This test case of our methodology shows the value of our approach, that predicting heat storage rate based on the rate-limiting physical phenomenon as a function of time is an effective approach to modeling LHTES devices.

25 ENERGY STORAGE↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

First-Principles Study of Surface Properties of U{sub 3}Si{sub 2} for Accident Tolerant Fuel

Due to the high thermal conductivity and uranium density compared to UO{sub 2} fuel used in current Light Water Reactors (LWRs), U-Si fuels, such as the U{sub 3}Si{sub 2} and U{sub 3}Si{sub 5} compounds, are currently under investigation for the development of accident tolerant fuels (ATF) in response to the Fukushima Daiichi accident in Japan. In order to perform meaningful mesoscale and macroscale fuel performance simulations, fundamental material properties of nuclear fuels should be known. One example is the simulations of fission gas bubble evolution and gas bubble-induced swelling in nuclear fuels using mesoscale methods, such as phase field and rate theory. To properly describe the gas bubble behavior in the fuel, it is crucial to know the surface energies of the fuel. Unfortunately, there is no study of the surface properties of U{sub 3}Si{sub 2} from both experiment and theory. To this end, we systematically investigate the surface properties of U{sub 3}Si{sub 2} using density functional theory (DFT) based first-principles calculations. Of several low index facets, we studied the (001) and (100) surfaces of tetragonal U{sub 3}Si{sub 2} with all the possible terminations. Since all the studied surfaces are not stoichiometric, we also investigated the effect of elemental chemical potentials on the formation of surface. We systematically investigate the surface properties of U{sub 3}Si{sub 2} using density functional theory calculations. Of several low-index facets, the (001) and (100) surfaces show similar averaged surface energies calculated using even number of layers. At the high Si chemical potential, the Si-terminated U{sub 3}Si{sub 2}(100) and U{sub 3}Si{sub 2}(001) surfaces show lower energies than those of the U-terminated surfaces, with U,Si-terminated U{sub 3}Si{sub 2}(001) as the most stable surface. The surface energies of all non-stoichiometric U{sub 3}Si{sub 2} are linearly dependent on the Si chemical potential. The surface energies of U-terminated surfaces decrease as Si potential decreases. The U{sub 2}-terminated U{sub 3}Si{sub 2}(100) surface becomes that the most stable surface at U-rich side. Based on the current results, an averaged surface energy of U{sub 3}Si{sub 2}, i.e., 1.45 J/m{sup 2}, is recommended for the description of the rounded facets of gas bubbles in U{sub 3}Si{sub 2} by mesoscale simulations. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

First-principles calculation of lattice distortion, electronic structure, and bonding properties of GeTe-based and PbSe-based high-entropy chalcogenides

The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.

Physics↗
DOE OSTI2018

First-principles study of electronic structure and Fermi surface in semimetallic YAs

In the course of searching for new systems, which exhibit nonsaturating and extremely large positive magnetoresistance, electronic structure, Fermi surface, and de Haas-van Alphen characteristics of the semimetallic YAs compound were studied using the all-electron full-potential linearized augmented-plane wave (FP–LAPW) approach in the framework of the generalized gradient approximation (GGA). In the scalar-relativistic calculation, the cubic symmetry splits fivefold degenerate Y- d orbital into low-energy threefold-degenerate and twofold degenerate doublet states at point around the Fermi energy. Furthermore one of them, together with the threefold degenerate character of As-p orbital, render the YAs semimetal with a topologically trivial band order and fairly low density of states at the Fermi level. Including spin–orbit (SO) coupling into the calculation leads to pronounced splitting of the state and shifting the bands in the energy scale. Consequently, the determined four different 3-dimensional Fermi surface sheets of YAs consists of three concentric hole-like bands at and one ellipsoidal electron-like sheet centred at the X points. In full accordance with the previous first-principles calculations for isostructural YSb and YBi, the calculated Fermi surface of YAs originates from fairly compensated multi-band electronic structures.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Mixed-Domain Charge Transport in the S-Se System from First Principles

Lithium-sulfur (Li-S) batteries are emerging systems heralded for the inexpensive cathode component (in S)and higher energy density than that of typical Li-ion batteries [1–5]. For the latter reason in particular, they are drawing high levels of interest for long-range electric air-craft applications where high gravimetric energy density is of utmost importance. Unfortunately, S by itself has a much too low electrical conductivity to be useful asa stand-alone cathode material. One effort to alleviate these shortcomings is to alloy sulfur with selenium [6–13],which is expected to improve the paltry conductivity of S as Se, depending on the phase, is 10∼30 times more conductive. This, however, comes at the cost of the heavier mass of Se that works against the gravimetric energy density. The charge transport behavior in the S1−xSex alloy system, not yet fully known or understood, is necessary for one to navigate the trade-offs between conductivity and other variables such as energy density en route to ultimately determining the optimum Se content, on demand. Charting carrier mobilities and conductivities in the S-Se system throughout the compositional spectrum from first principles at a predictive level is therefore critical.

Junsoo Park↗
DOE OSTI2022

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗