Search NASA⌕ Search

SEARCH · Search NASA

Results for “HCP Mg”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A comprehensive diffusion mobility database comprising 23 elements for magnesium alloys

We report that reliable experimental diffusion coefficients of 10 key alloying elements in Mg obtained by the present authors together with experimental data in the literature enabled us to perform a systematic test of the reliability of diffusion coefficients obtained from DFT calculations. The computed activation energy values were found to be quite accurate (mostly within 0.2 eV) but the computed pre-factors were less reliable. Such insights allowed us to develop a practical and yet robust strategy to perform diffusion mobility assessments by adopting the computed activation energy while fitting only the pre-factor when available experimental data are limited to a narrow temperature range. The overall good agreement between the DFT data and experimental data also gave us the confidence to employ the computed data for those that were still missing or inaccessible from experimental measurements. A systematic assessment of both the measured and computed diffusion data in hcp Mg was performed using the above holistic approach to yield the most comprehensive open Mg mobility database to date, comprising 23 elements (Mg, Ag, Al, Be, Ca, Cd, Ce, Cu, Fe, Ga, Gd, In, La, Li, Mn, Nd, Ni, Pu, Sb, Sn, U, Y, Zn). This more reliable mobility database will contribute to future development of advanced Mg alloys. The holistic approach developed in this study will be very beneficial to the future establishment of reliable mobility databases for other alloy systems as well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Micromechanical Fields Associated with Irregular Deformation Twins in Magnesium

Understanding and controlling the development of deformation twins is paramount for engineering strong and stable hexagonal close-packed (HCP) Mg alloys. Actual twins are often irregular in boundary morphology and twin crystallography, deviating from the classical picture commonly used in theory and simulation. In this work, the elastic strains and stresses around irregular twins are examined both experimentally and computationally to gain insight into how twins develop and the microstructural features that influence their development. A nanoprecession electron diffraction (N-PED) technique is used to measure the elastic strains within and around a $\left\{ {10\overline{1}2} \right\}$ tensile twin in AZ31B Mg alloy with nm scale resolution. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) crystal plasticity model of the same sub-grain and irregular twin structure is employed to understand and interpret the measured elastic strain fields. The calculations predict spatially resolved elastic strain fields in good agreement with the measurement, as well as all the stress components and the dislocation density fields generated by the twin, which are not easily obtainable from the experiment. The model calculations find that neighboring twins, several twin thicknesses apart, have little influence on the twin-tip micromechanical fields. Furthermore, this work reveals that irregularity in the twin-tip shape has a negligible effect on the development of the elastic strains around and inside the twin. Importantly, the major contributor to these micromechanical fields is the alignment of the twinning shear direction with the twin boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shock melting and the hcp-bcc phase boundary of Mg under dynamic loading

The high-temperature and high-pressure response of magnesium has been investigated through shock-release experiments performed up to shock melting. The longitudinal and bulk sound speeds of Mg are reported along the Hugoniot from 25 to 56 GPa and used to determine the elastic properties and Grüneisen parameter. The hexagonal close-packed (hcp)–body-centered cubic (bcc) phase transition is marked by a reduction in the determined shear wave speed. Thermal softening is observed to begin between 40 and 44 GPa, with incipient melt at 55.5 GPa, in close agreement with previous diffraction measurements under laser shock loading. Examination of the release profiles showed that two different responses were observed, depending on the peak stress. When plotted in a stress-energy phase diagram, the two responses are observed to form separate lines that intersect the static hcp-bcc phase transition and incipient melt. The results indicate that the hcp-bcc transition occurs on the Hugoniot at 28.4 GPa and place the hcp-bcc-liquid triple point ~20 GPa. In conclusion, a multiphase equation of state is developed which places the melt boundary below the previously reported static measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Solutes that reduce yield strength anisotropies in magnesium from first principles

Using Labusch-type solid solution strengthening models parameterized with DFT-computed solute-dislocation interaction energies, we perform a computational search for 63 solutes across the periodic table to find those that lower anisotropy ratios (non-basal to basal CRSS) of magnesium potentially increasing its ductility per the von Mises criterion. For this purpose, we compute changes in strength for solutes as a function of composition and temperature, and compute anisotropy ratios for solutes that include both rare earth and non-rare earth elements. Here we specifically focus on solute-dislocation interaction energies in the following DFT-optimized dislocations as representative of three non-basal plastic deformation modes: $\langle c + a \rangle$ edge, (10$\bar1$2) tension twinning edge, and the (10$\bar{1}$1) compression twinning edge. We find that solute-induced changes in non-basal deformation modes can be approximated using a second-order polynomial in the size misfit of the solutes, which permits rapid screening of solutes. Our approach to identify solutes known to improve strengthening incorporates solute solubility, and suggests other solutes that not have been previously explored for strengthening. The 8 rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios at their optimal concentrations, are: Gd, Tb, Dy, Nd, Ho, Er, Tm, and Yb. The 12 non-rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios, are: Y, Mn, Sc, Pb, Ca, Ag, Bi, Tl, Zn, Li, Ga, and Al. Of these, Gd, Nd, Er, Yb, Y, Mn, Ca, Zn, Li, and Al are used in commercial Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Geometric compatibility measure m' for twin transmission: A predictor or descriptor?

In this work, the geometric compatibility factor m' is critically analyzed to assess whether it can be used to interpret/predict twin transmission (TT) across grain boundaries (GBs). This geometric measure is widely used to relate the likelihood of TT to the misalignment of both the shear and plane-normal directions within a twin set (i.e., incoming and outgoing twin). Here, using a large set of electron back scattering diffraction (EBSD) data, a detailed statistical analysis of twin-GB interactions is performed for {${1\bar{01}}2$} tensile twins in hexagonal close-packed (HCP) metals Mg, Zr, and Ti at different strain levels. In addition, a full-field crystal plasticity model is employed to quantify the role of local stresses and the applicability of m' as a criterion for the TT process. This combined study addresses the following three main questions: (i) What is the fidelity of m' in describing experimentally observed TTs? (ii) Can m' be used as a metric to predict/anticipate TT? (iii) Does m' naturally capture local stress effects? As a descriptor, m' cannot rationalize ~25% of TT events observed in Mg or more than 50% of TT events in Zr and Ti. As a predictor, the m'-measure does not predict TT events in over ~50% of twin-GB interactions analyzed. Further, the applicability of m' to describe and predict TT events decreases with an increase in elastic anisotropy, plastic anisotropy, and macroscopic strain levels. Finally, the twinning simulations reveal that m' does not capture the key effects of local stresses on variant selection upon twin transmission. The local stress induced by the twinning shear transformation plays a dominant role in driving the TT process compared to the geometric alignment of the constituting twins, i.e., m'.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Grain boundary structure search by using an evolutionary algorithm with effective mutation methods

Grain boundaries (GBs) accommodate the misorientation between adjacent grains in a polycrystalline material. GBs are geometrically described by the macroscopic and microscopic degrees of freedom. Besides, at the atomistic level, GBs exhibit complicated behaviors under varying thermodynamic conditions. The complexity of atomistic GB structures demands stochastic searching for possible states. The effectiveness of stochastic search methods relies on techniques to recreate and select atomistic structures. In this work, we developed a new mutation operator that can induce direct and collective atomistic structure changes to boost the search efficiency of exploring GB structures with evolutionary algorithms (EA). We implemented the mutation methods along with innovative selection, crossover, boundary condition preprocessing methods to form an EA-based package to explore GB structures in grand canonical ensembles with atomistic simulations. We used this package to study the [001] symmetric tilt grain boundaries (STGBs) in FCC copper (Cu), the [110] STGBs in BCC tungsten (W), and the $[1\bar{2}10]$ STGBs in HCP magnesium (Mg). The results show that our design and implementation based on new mutation procedures, selection, and boundary conditions provide a high-quality search of atomistic GB structures in the grand canonical ensemble for different crystal lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High‐Throughput Evaluation of Hardening Coefficients of Eight Alloying Elements in Magnesium

Liquid–solid diffusion couples (LSDCs) are employed to generate a composition gradient in the single‐phase hexagonal closed‐packed (hcp) solid solution with compositions up to the solubility limit of various solutes in Mg. Nanoindentation scanning across the composition gradient in LSDCs allows effective evaluation of composition‐dependent hardness of eight alloying elements (Al, Ca, Ce, Gd, Li, Sn, Y, and Zn) in the hcp Mg phase. The hardening coefficients, an indicator of the potency of solid‐solution hardening, are evaluated from the measured composition‐hardness data and correlated with various materials properties such as atomic radius, shear modulus, and elastic modulus of the solutes. The rank of hardening potency of Al, Gd, Sn, Y, and Zn measured by nanoindentation is in good agreement with that measured by microindentation reported in the literature. The hardening coefficient (potency) from the strongest to the weakest is Ce > Ca > Y ≈ Gd > Zn > Al ≈ Sn > Li in Mg‐based hcp binary solid solutions. The hardening coefficient is found to be closely correlated with the strengthening potency.

36 MATERIALS SCIENCE↗
DOE OSTI2023

White Laue and powder diffraction studies to reveal mechanisms of HCP-to-BCC phase transformation in single crystals of Mg under high pressure

Mechanisms of hexagonal close-packed (HCP) to body-centered cubic (BCC) phase transformation in Mg single crystals are observed using a combination of polychromatic beam Laue diffraction and monochromatic beam powder diffraction techniques under quasi-hydrostatic pressures of up to 58 ± 2 GPa at ambient temperature. Although experiments were performed with both He and Ne pressure media, crystals inevitably undergo plastic deformation upon loading to 40–44 GPa. The plasticity is accommodated by dislocation glide causing local misorientations of up to 1°–2°. The selected crystals are tracked by mapping Laue diffraction spots up to the onset of the HCP to BCC transformation, which is determined to be at a pressure of 56.6 ± 2 GPa. Intensity of the Laue reflections from HCP crystals rapidly decrease but no reflections from crystalline BCC phase are observed with a further increase of pressure. Nevertheless, the powder diffraction shows the formation of 110 BCC peak at 56.6 GPa. The peak intensity increases at 59.7 GPa. Upon the full transformation, a powder-like BCC aggregate is formed revealing the destructive nature of the HCP to BCC transformation in single crystals of Mg.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The {10$\overline{1}$2} non-cozone twin-twin interactions in Mg: A stability and mobility study using 3-D atomistic simulations

Due to easy activation of tensile twinning in Mg, multiple {10$\overline{1}$2} twin variants can interact with each other and form twin-twin boundaries over the course of plastic deformation. Previous studies using 2-D settings provide only a partial understanding of these interactions, especially the non-cozone ones. Here, atomistic simulations are used to study the 3-D structural characteristic and evolution of the non-cozone {10$\overline{1}$2} twin-twin junctions. The study reveals the existence of new twin-twin boundaries (TTBs) such as TTB BP and TTB K2 , formed after the interaction between the basal prismatic and conjugate twin interfaces with the coherent twin boundary. Here, for both non-cozone twin-twin interactions, the {10$\overline{1}$2} TTB and its associated twin-twin junctions are found to play a major role in the {10$\overline{1}$2} twin's stability and mobility. Specifically, they promote the growth of the 3-D twin in both the normal and forward directions during the interaction and hinder the detwinning process upon unloading.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Critical Shape for the Growth of Grain Boundary Twin Embryos in Mg and Mg Alloys: Crystal Plasticity Modeling

Application of polycrystalline hexagonal close packed (HCP) metals in engineering designs has been constrained by their anisotropic responses due to twinning and limited plasticity. In deformation, twins most often initiate at grain boundaries (GBs), and thicken and propagate across the grain. In this work, the GB twin embryos in Mg and Mg alloys, and the conditions that influence their propagation are investigated. Using a micromechanical crystal plasticity model, the role of embryo shape on the driving forces prevailing at the embryo boundaries that could support its expansion is studied. The modeled embryos are either planar, extending more in the shear direction than normal to the twin plane, or equiaxed. Results show that the thinner the embryo, the greater the driving forces for both thickening and forward propagation. Alloys with low prismatic-to-basal critical resolved shear stress (CRSS) ratios promote embryo thickening and large CRSS values for the slip mode that primarily accommodates the twin shear encourage propagation. The neighboring grains with orientations that enable local accommodation of the embryo twin shear by pyramidal slip promote forward propagation but have little effect on thickening. When two like embryos lie along the same GB, their paired interaction promotes forward propagation but hinders thickening.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dataset of simulated vibrational density of states and X-ray diffraction profiles of mechanically deformed and disordered atomic structures in Gold, Iron, Magnesium, and Silicon

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2021

The effects of free surfaces on deformation twinning in HCP metals

Deformation twinning is a predominant mode of plastic deformation for hexagonal close packed metals, like Mg and Ti. The heterogenous microstructure and the local stresses associated with twinning play a key role in their mechanical response and fracture. Surface analyses, like electron microscopy, are frequently employed to spatially map microstructure and micromechanical fields in order to study twinning behavior. However, these measurements are inherently influenced by the vicinity of the free surface. In this work, an elasto-visco-plastic fast-Fourier-transform (EVP-FFT) polycrystal modeling approach is employed to investigate the effects of free surfaces on twin development before and after loading. We compare calculated micromechanical fields on free surfaces with those calculated inside the bulk and, in some cases, experimental surface measurements. The results indicate that the creation of free surfaces can promote twin propagation and growth and can influence twin morphology by causing a twin lamella to become larger, more blunted and irregular. The structure along the twin boundaries are also affected, due to the higher driving stresses that extend prismatic-basal and basal-prismatic facets. Furthermore, free surfaces invoke different slip activities in the twin and the surrounding parent crystal by enhancing basal, prismatic and pyramidal slip in some localized regions, while reducing slip in others. We demonstrate that the simulated free-surface effects lead to better qualitative and quantitative agreement with experimental measurements from scanning electron microscopy and digital image correlation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase-field model of deformation twin-grain boundary interactions in hexagonal systems

Here, a phase-field model for describing dynamic interactions between deformation twins and grain boundaries in a hexagonal close packed (HCP) metal is established. It is applied to simulating the coupled evolution mechanisms of grains and twins in magnesium (Mg) by probing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries. We analyze the effect of the strain relaxation near a grain boundary, which is related to the geometric compatibility arising from the misorientation angle between adjoining grains. Phase-field simulations demonstrate that twin transmission across a grain boundary into a neighboring grain leads to grain boundary migration towards the neighboring grain with a reduced grain boundary width. The preferred nucleation site of new twin variant in the neighboring grain is related to the elastic interaction energy distribution. These predicted transmission behaviors agree well with existing experimental observations and molecular dynamics simulations. By analyzing set of systematic phase-field simulation results, we establish a rotation-angle-related twin variant selection rule for analyzing the transmission of $\{10\bar{1}2\}$ deformation twins across grain boundaries in Mg and its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fundamental Thermodynamic, Kinetic, and Mechanical Properties of Lithium and Its Alloys

Lithium alloying reactions are beneficial in promoting uniform plating and stripping of lithium metal in all-solid-state batteries. First-principles calculations are performed to predict thermodynamic, kinetic, and mechanical properties of lithium and several important Li–M alloys (M = Mg, Ag, Zn, Al, Ga, In, Sn, Sb, and Bi). While the Li–Mg binary system forms a solid solution, most other lithium–metal alloys prefer stoichiometric intermetallic compounds with common local motifs that enable fast Li diffusion. Lithium and Li-rich alloys exhibit an unusually flat energy landscape along paths that connect BCC to close-packed structures like FCC and HCP, with important implications for mechanical properties. Very low migration barriers for Li diffusion that rival those of superion conductors are predicted, both in pure Li and in Li–M intermetallics. However, vacancy concentration, which is crucial for substitutional diffusion, is predicted to be low in metallic Li and most Li–M intermetallics. Compounds such as B32 LiAl and LiGa as well as D0 3 Li 3 Sb and Li 3 Bi exhibit structural vacancies at higher ends of their voltage windows, which together with low migration barriers leads to exceptionally high Li mobilities. In the Li–Mg solid solution, the addition of Mg is found to decrease the vacancy tracer diffusion coefficient by an order of magnitude.

36 MATERIALS SCIENCE↗