{'fraction_non_alphanumeric': 0.02347013366750209, 'fraction_numerical': 0.019227756892230576, 'mean_word_length': 4.815607682380627, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The role of metal substitutions in the development of Li batteries, part II: solid electrolytes

Of the three main components in state-of-the-art Li-ion batteries (cathode, anode and electrolyte), metal substitutions into known structures have primarily served a pivotal role in developing the cathode and solid electrolytes only. Previously, we reviewed the role of substitutions in optimizing cathodes, and these included multiple classes of materials that can only be commercialized with improved electrolytes able to operate at higher potentials than liquid carbonate-based electrolytes. Solid electrolytes are considered by many to be the path forward in this respect. In this review, we will discuss the design principles established for solid electrolytes, as well as both the successes and limitations of using substitutions to improve the performance metrics of highest import for critical applications such as electric vehicles and grid storage. A massive research push is underway in order to bring all-solid-state Li batteries to widespread use, particularly for electric vehicles, and metal substitutions are playing a critical role in developing the solid electrolytes to the point where they will, hopefully, out-perform the established Li-ion batteries with liquid electrolytes. This review further emphasizes metrics that are badly needed in further studies in order to systematically optimize solid electrolytes in a way that correlates well to better battery performance. Although the focus of this review article is on metal substitutions, we also discuss a few cases where non-metallic substitutions have been needed to bring a class of materials into viability where further metal substitutions become significant.


Introduction
There is a global need for a safer, higher energy and power density Li rechargeable battery due to the projected massive boost in electric vehicle demand. 1 If one also considers the need for grid storage, the requirement for better batteries becomes extremely critical. All-solid-state Li batteries (ASSLB) have the potential to fulfill these demands by overcoming the limitations of the commercial Li-ion battery that relies on liquid electrolytes. In an ASSLB, the conventional flammable organic liquid electrolyte is replaced by a nonflammable solid electrolyte. Inorganic solid electrolyte can potentially overperform liquid electrolyte by facilitating the use of high potential cathode, which were discussed in part I of this review, such as LiNi 0.5 Mn 1.5 O 4 (5 V) and LiCoPO 4 (4.8 V), and they also may permit the use of Li metal anodes that represent important safety hazards in traditional Li-ion batteries. Due to the narrow electrochemical stability window of current liquid electrolytes, high potential cathodes can oxidize and consume it increasing the internal resistance; 2 moreover, some cathodes tend to degrade by dissolution because of the liquid nature of the electrolyte. 3 In contrast, inorganic solid electrolytes with wider electrochemical stability windows would allow the integration of high potential cathodes and inhibit the degradation due to dissolution. 4 Furthermore, it aids in developing new chemistry cathodes such as Li-S batteries, 5 where solid electrolytes prevent polysulfide dissolution into the electrolyte. Solid electrolytes are also being explored for use in Li-O 2 batteries. 6 Although replacing the conventional carbon anode with lithium metal increases the energy density, but dendrites can form. The dendrite growth can be hindered by the high bulk modulus of the inorganic solid electrolyte (ISE) at high current densities, thus achieving high energy and power density. 7 Besides these advantages, bipolar multicell stack configuration becomes possible with ISE, as well as operation at lower temperature where liquid electrolyte freezes. 8 From the above brief discussion of the required role of ISEs, it is clear that a number of metrics are needed to determine the viability of a particular solid electrolyte, as illustrated in the radar plots in Fig. 1. The requirements for a good ISE are to have high ionic conductivity at a wide range of temperatures, including room temperature, low electronic conductivity to prevent short circuit and limit dendrite growth within the electrolyte, wide electrochemical stability window and high electrode wettability that leads to a low interfacial area specific resistance (ASR). Other properties such as mechanical properties (some flexibility is required to accommodate volume changes in the electrodes), and thermal stability are also important but are not often considered in the literature and so will not be discussed extensively in this review. The methods used to determine the key metrics will also be reviewed later in this article. It should be noted that the radar plot for a particular family (say oxides) in Fig. 1 is a generalization, and metal substitutions can have significant consequences on multiple metrics in each radar plot. As such, the landscape for optimizing solid electrolytes proves to be very complex, and the need for systematic studies that consider numerous metrics is extremely high.
One way to achieve better performance ISE is by tuning the structures of already known ISEs with cationic/anionic substitutions. Some of these structures are oxides such as perovskite, argyrodite, Li superionic conductor LISICON-type, Na superionic conductor NASICON-type, anti-perovskite or sulfides such as Thio-LISICON, and Li 10 GeP 2 S 12 (LGPS) like structures. Such substitutions can enhance the bulk ionic conductivity by increasing the number of mobile Li ions 9 and tuning the potential energy landscape of Li ions. There is a long list of mechanisms involved in conductivity enhancement by substitution including inducing Li sublattice disorder, 10 lowering bottleneck energy, 11 or lowering the phase transition temperature to a superionic conductor phase. 9 Substitution can additionally influence the grain boundary conductivity and sinterability (densification) of the ISE. 12 The effect of substitution is not limited to these properties and may often have unexpected detrimental impacts on other properties such as the electronic conductivity and electrochemical stability window of the ISE. It is therefore essential to do systematic and proper screening of multiple key ISE metrics in order to fully determine whether a substitution is beneficial or not. Understanding the exact mechanism of enhancement can give us hints for engineering better ISE such as multi-substitution.
In this review, we start by describing the different synthesis methods of ISE. We then consider in detail the experimental methods to be used to determine the key metrics needed to screen solid electrolytes. Finally, we survey the role that metal substitutions have played in the development of the various families of known ISE structures that are still being studied extensively for all-solid-state Li batteries. For each family, the influence of the substitutions on the different properties of the ISE and the underlying mechanisms are discussed. Additionally, the different substitution strategies used to optimize specific electrochemical properties are summarized.

Synthesis of solid electrolytes
ISEs are synthesized by various methods and can be made either into bulk powders or thin films. The methods used to make bulk powders include solid-state synthesis and wet-chemical methods such as co-precipitation, sol-gel, solvothermal, and hydrothermal. In the solid-state synthesis method, the reagents are ground together by ball milling or mortar and pestle and then pressed into green pellets for heat treatment. The drawback of the solid-state method is that high temperatures and long sintering times are required because the synthesis reactions happen at solid-solid interfaces, where energy is needed to overcome lattice energy and diffuse anions/cations. 14 Usually intermediate grinding/repelleting is mandatory to form fresh interphases to yield pure phase materials. For instance, one calcination step is not enough to completely react all the raw materials in Al 3+ substituted Li 7 La 3 Zr 2 O 12 (LLZO) prepared by the solid-state method. 15 Furthermore, it is not possible to produce a dense solid electrolyte using the solid-state method in a single heat step specially if non-oxide reactants are used such as carbonate based reagents, because cracks and voids form inside the bulk of the pellet during calcination; therefore, the calcined pellets must be re-ground and repressed into pellets for the sintering step. Thus, the calcination temperature and the property of the calcined powder formed after grinding influence the sintering temperature and the densification of the final ISE. For example, varying the calcination temperature of Al 3+ substituted LiTi 2 (PO 4 ) 3 LTP electrolyte from 750 to 900 1C with the same sintering temperature increases the relative density of the final pellet from 85% to 98%. 16 Sometimes, these solid mixtures are heated under a specific atmosphere; for example, heating the unsubstituted LLZO reagents mixture under CO 2 can mitigate the lithium loss, thus stabilizing the cubic phase. 17 On the other hand, if volatile or airsensitive reagents such as sulfur are used, the mixing/pelletizing must be done in an argon-filled glovebox. The mixture can then be sealed in a tube usually made of quartz for the heat treatment. A carbon-coated tube is used to avoid contamination from the tube, or the samples are placed in a carbon crucible.
The alternative methods to solid-state synthesis are all solution based. In the co-precipitation method, reagents are dissolved into aqueous solutions mixed and then a precipitating agent is added such as ammonium hydroxide. After filtering the precipitate, it is dried, preheated, and then pressed into pellets. The co-precipitation is less time consuming than the other methods. However, the drawback is that all the cations should be soluble in the aqueous mixture and do not precipitate before adding the precipitator. The sol-gel method by contrast proves to be far more versatile. The main types of sol-gel synthesis are the citrate (Pechini method) and polymerization. The citrate method involves adding citrate as a chelating agent to an aqueous solution containing the reactants. Once the solvent volume starts to reduce, a gel forms and total removal produces solid powder which then goes through the same high temperature heat treatment as the methods mentioned above. The formation of citrate complexes facilitates atomic scale mixing of the cations thus achieving a high homogeneity compared to solid-state synthesis. While the polymerization involves the hydrolysis of alkoxide reagents followed by the condensation forming a gel that is then dried, calcined, and pressed into pellets. The sol-gel method yields small size powders lowering the sintering temperature compared to the solid-state method.
In the synthesis approaches mentioned above, the last step is the sintering and many techniques have been successfully used to perform the sintering, such as furnace heating, floating zone (single crystal), 18 spark plasma, 19 and microwave assisted. 20 The advantage of spark plasma and microwave assisted sintering over furnace heating is that less high temperature dwell time is needed to yield a dense ISE but with less grain size growth. 21 Some synthesis methods used for a specific family of materials such as the melt-quench method for the glass-ceramic type solid electrolytes do not need a sintering step where the sample is heated until melt and then quenched to form glass after it is annealed to crystalize it into glassceramic material. 22 The sintering process can be done under a specific atmosphere for instant O 2 -assisted sintering of a particular garnet solid electrolyte increased its relative density from 95% to 97% compared to the sintering conducted under the same conditions but in air. 23 Usually, the sintering step requires high temperatures (as high as 1400 1C is not uncommon for ceramic oxides for example) where Li volatilization from the sample and contaminants diffusing from the substrate become critical issues to be mindful of. For example, during the sintering of Li 3x La (2/3)Àx TiO 3 (LLTO) at 1250 1C for 5 h the Li content drops from 0.5 to 0.3 per formula unit. 24 Another concern is side reactions with the substrate. For example, sintering LLZO pellet at 1220 1C for 12 h in an alumina crucible causes the side in contact with the alumina to have a 16 at% unintentional Al content. 25 At high temperatures, total prevention of lithium loss is not possible; thus, an excess amount of lithium (10%) is often added to the starting stoichiometry. The contamination from the substrate occurs typically with alumina substrates and can be avoided by using platinum or zirconia substrates or by using a sacrificial powder. 25,26 The sacrificial powder approach involves covering the sample pellet with powder made of the same composition as the sample (it also sits on a bed of this powder). This sacrificial powder will then interact with the substrate (preventing sample/substrate contact) and lose more Li than the sample, thereby yielding a pellet with the desired composition and no contaminant. The composition of the prepared ISE must then be confirmed by elemental analysis, such as inductively coupled plasma (ICP).
The synthesis methods of ISE mentioned above are used to prepare thick electrolytes (pellets) typically used for screening of electrolytes for bulk-type all-solid-state batteries. However, when making full solid batteries it is advantageous to make this layer as thin as possible to maximize energy density, and a few such cases will be discussed herein. Numerous thin-film ISE deposition techniques have been utilized to make full solid batteries. These are divided into two categories: chemical vapor deposition (CVD) (e.g. atomic layer deposition ALD 27,28 ) and physical vapor deposition (e.g. pulse laser deposition PLD). 29 CVD is quite common due to the versatility of this method. In CVD, the film is formed through the reaction of the gaseous reactants. This reaction can be driven by (1) increasing the diffusion of the gaseous reactants via reducing the pressure (low-pressure CVD), 30 (2) plasma (plasma enhanced CVD), 31 (3) the heat induced by laser (laser CVD). 32,33 The other more frequently used CVD method is the metalorganic chemical vapor deposition (MOCVD) in which organometallic (metal alkyl) reagents are used. 34 The choice of specific CVD technique depends on the nature of the materials and the required thickness and adhesion. An advantage of these deposition techniques is that it allows to the construction of complicated 3D hierarchical designs. The other common method used to form a thin layer ISE is by spin coating the precursors (sol-gel) and then heat treating. 35 Other than thin film and bulk forms of ISE, nanostructured ISE such as nanowires and nanoparticles are being explored. These nanowires can be fabricated by electrospinning. 36 The nanostructured ISE have been used in a polymer-ISE composite electrolyte where the ISE can act as an ionic channel and a filler lowering polymer glass transition temperature, thus enhancing the overall room temperature conductivity. 37 However, typically, most of the ISE screening is done on bulk pellets and will therefore be the primary focus of this review.

Screening methods for solid electrolytes
The continued improvement of the performance of ISE is essential for developing an ASSLB functioning at room temperature and out-classing the established Li-ion batteries. The inorganic solid electrolyte in ASSLB contributes to the total internal ionic resistance by the conductivity of its bulk, grain boundaries, and the electrolyte/electrode interfaces. There are several key properties that ISE must fulfill, the most fundamental of which are: high ionic conductivity, low electronic conductivity, and high stability with the electrodes under operating conditions. This section discusses the various methods used to determine these properties and presents some methods used to further understand the diffusion mechanisms at play.

Ionic conductivity
Firstly, an ISE should be a good ionic conductor (410 À4 S cm À1 ) at room temperature and an electronic insulator (o10 À10 S cm À1 ). It should be noted that liquid electrolytes have ionic conductivities on the order of 10 À2 S cm À1 and many researchers consider 10 À3 S cm À1 to be the minimum ionic conductivity required of solid electrolytes to be truly competitive. Only the ion of interest (lithium here) should be mobile with a transport number close to one because any small electronic conductivity can cause dendrite nucleation from the bulk of the electrolyte, especially at high current densities (s electronic 4 10 À10 S cm À1 at 1 mA cm À2 ). 38 The total electrical (ionic + electronic) conductivity of ISE is measured using electrochemical impedance spectroscopy by scanning MHz to mHz frequency range and applying 10-100 mV potential after depositing electrical contacts on the two faces of a polished ISE pellet. This contact can be deposited by thermal evaporation, magnetron sputtering, e-beam evaporation, or applying a colloidal gold paste and then annealing it. 39 Other contacts such as silver and platinum can also be used. The acquired impedance spectrum is usually represented by a Nyquist plot. A typical Nyquist plot of an ISE consists of highfrequency region semicircles and a low-frequency spike as shown in Fig. 2(a). The spectrum is fitted to an equivalent electronic circuit. The semicircles are modeled by a resistor and capacitor in parallel and the spike as a capacitor as shown in Fig. 2(a). The capacitor is usually replaced by a constant phase element to fit the non-ideal capacitance of the electrochemical system components. The low-frequency spike is assigned to the capacitance of the blocking electrode, while the semicircles are assigned to the bulk and grain boundary. The capacitance of the semicircles helps to differentiate between the bulk and grain boundary. A typical bulk and grain boundary capacitances are 10 À12 F and 10 À11 -10 À8 F respectively which are calculated from a standard permittivity of B10 using parallel plate capacitor equation. 40 For instance, the grain boundary conductivity in LLTO is much lower than the bulk due to the orientation of the domains blocking lithium transport from one domain to the other. 41 In contrast, in LLZO some researchers found grain boundary conductivity higher than the bulk. 42 The bulk and grain boundary are not always resolved as two semicircles. However, sometimes only one semicircle is observed, and that is due to the property of the material or the frequency range used for the measurement. If only one semicircle is observed (total conductivity) with the other one out of the measurement frequency region, then measuring the impedance at low temperatures (À100 to À130 1C) can separate the bulk and grain boundary conductivities. 43,44 Prior to studying the impact of substitutions, it is critical to correctly assign the contribution of the bulk and grain boundary to the total conductivity to understand the effect of specific substitutions.
It is also important to recognize the usefulness of performing the conductivity measurements as a function of temperature. The key transitions required to enable conduction (e.g. lithium sublattice order-disorder, defects formation) can occur at various temperature regimes. 45 For example exploring these transitions in intrinsic ion conductors which conductivity depends on the intrinsic defects that increase with temperature. Likewise studying the temperature dependent conduction mechanisms of extrinsic ion conductors which conductivity is due to substitution like doping in extrinsic semiconductor. Typically, such data is plotted as log(sT) vs. 1/T to investigate whether a specific ion conductor exhibits Arrhenius or non-Arrhenius type ion conduction. The activation energy and preexponential factor extracted from the Arrhenius plot are used for comparing performances and developing solid electrolyte by helping to understand the underlying mechanisms causing the conduction differences. This is especially useful when conductivities are so small that they cannot be measured at room temperature. In such cases, the impact of substitutions represents important information to develop the electrolyte but cannot be observed unless impedance measurements are performed at higher temperatures, and conductivities can then be extrapolated to room temperature. Next, the total ionic conductivity is calculated by measuring the electronic conductivity and subtracting it from the total electrical conductivity.

Electronic conductivity
The Hebb-Wagner method is used to measure the electronic conductivity by applying DC polarization on an asymmetric (À)Li/SSE/Cu(+) cell with a reversible and an irreversible (ion blocking) electrode such as lithium and copper respectively. 46,47 The electronic conductivity is calculated from the voltage vs. steady-state current plot of the chronoamperometry measurements at different voltages as shown in Fig. 2(b). The lithium transference number can then be calculated from the ionic contribution to the DC total conductivity obtained at t = 0 s prior to steady-state or total conductivity from the electrochemical impedance spectroscopy (EIS) measurement. In case the investigated ISE is not stable against lithium electrode the upper limit of electronic conductivity can be measured using two blocking electrodes instead. 38 Thus, it is important for comparing the electronic conductivity to have similar cell configuration moreover the same applied potential because for some ISEs the electronic conductivity is not ohmic at all applied potentials. 48 Electronic conductivity below 10 À10 S cm À1 is needed for dendrite free ASSLB operating at 1 mA cm À2 . 38

Electrochemical stability at the anode
Besides appropriate conductivities, solid electrolytes require a wide electrochemical stability window. Alternatively, they may form a stable ionically conductive solid electrolyte interphase (SEI) with the electrode. The stability against the anode is often done vs. lithium as this is considered to be key in order for solid batteries to out-perform Li-ion batteries in terms of energy density. The instability against lithium metal is due to the lithium insertion reaction forming a reduced decomposition phase as shown in Fig. 3(b). Stability against lithium electrode is best probed by measuring the impedance (using EIS) of a Li/ISE/Li symmetric cell at different time intervals during electrochemical cycling. Fig. 4 shows the interfacial resistance between Nb 5+ substituted LLZO and metal electrode increases by time, indicating the formation of an electronic conductive phase propagating through the bulk of the electrolyte. 49 The same symmetric cell is also used to study the interfacial resistance between the lithium electrode and the ISE conditioned (sanding dry/wet, 50 acid etching, 51 applying a coating 52 ) surface signifying lithium wettability. Furthermore, the galvanostatic cycling (plating/stripping) of this symmetric cell at different current densities can demonstrate in progressive cycles if any dendrites will form causing shortage or if the internal resistance will continuously increase (measured voltage at same current density increase by each cycle), which could be due to bad electrode/electrolyte contact or electrolyte decomposition into ion resistive phase.

Electrochemical stability with the cathode
The upper limit of the electrochemical stability window is the tendency for lithium extraction forming a lithium deficient oxidized interphase, as shown in Fig. 3(a). This limit has been explored in the literature using a variety of methods. The most common is to use a Li/ISE/Au cell using cyclic voltammetry, as illustrated in Fig. 5(a). However, this method has been shown to systematically overestimate the electrochemical stability window due to the limited electrical contact between the current collector and the ISE (i.e. the poor electronic conductivity of the ISE prevents the reactions that would otherwise take place at the high potentials). 53 The alternative method that should be used is to make a Li/TE/ISE+CA/Au (where TE is a functional electrolyte, ISE is the electrolyte under investigation and CA is a conductive additive). The CA is used to increase the amount of electronic paths such that poor electronic conductivity does not prevent the parasitic reactions from taking place. The test electrolyte (TE) can either be a liquid organic electrolyte 54 or the ISE under investigation as long as its ionic conductivity is adequate. This technique systematically shows a more representative narrower electrochemical stability window which also fits better the computational predictions. For example, the cyclic voltammetry of Li 10 GeP 2 S 12 (LGPS) with added Pt black shows a lower oxidation limit of 2.5 V compared to stability up to 4 V as determined by the Li/ISE/Au CV method as shown in Fig. 5. 53 No stabilities determined by the Li/ISE/Au CV method will be reported in this review article as they are deemed inaccurate. It is worth to mention that electrochemical stability window is also investigated computationally using different approaches such as calculating the electronic bandgap and band position which only assumes electron transfer. The other two methods, the stoichiometry stability method and phase stability method, present the process of stability against the electrodes (lithium extraction/insertion) as shown in Fig. 3 more accurately by including the chemical potential of the lithium which contains the potentials of the electrodes. 55 Finally, the solid electrolyte should have optimal mechanical properties; for instance, it should be soft and deformable to wet the electrode better and be able to accommodate the volume expansion of the electrode; otherwise, if brittle, it cracks, but at the same time, it should have enough hardness to prevent the dendrite growth. Although the mechanical properties and thermal stability of great importance for the ISE performance, will not be discussed in this review.

Characterization of structures and mechanism determination
The structure is most commonly identified by powder X-ray diffraction. In a new crystal structure, synchrotron X-ray diffraction or neutron diffraction is required to resolve the new structure. The mechanism underlying lithium transport is predicted by conducting neutron diffraction and relaxation nuclear magnetic resonance at different temperatures to locate the exact positions of the lithium ion and reveal their local dynamics. Moreover, from the neutron diffraction data, the Maximum Entropy Method (MEM) gives information about the Li + diffusion pathways as shown in Fig. 6(b). 56 MEM is complemented with the Bond Valence (BV) method that can provide the path that minimizes the empirically determined valence  mismatch for lithium ion diffusion. BV calculations are extremely easy/cheap to perform and yield iso-valence mismatch surface maps as shown in Fig. 6(a). From BV calculations, the bond valence landscape (BVL) can be extracted which shows the relative energy of the barriers that can then be compared to the energy barriers calculated from computational simulation and from the lithium nuclear densities. Moreover, molecular dynamic calculations are utilized to predict the possible mechanisms and pathways of lithium diffusion. 57 All these techniques are complemented with each other to produce a more accurate description of lithium migration.

Impact of substitutions on solid electrolyte metrics
A great number of different substitutions into ISEs have been performed over the years. These have been done for a variety of reasons, we discuss the main driving forces for substitutions into ISEs in this section. First and foremost, substitutions are performed to improve the bulk ionic conductivity of the ceramic structure. The bulk ionic conductivity is governed by the composition and crystal structure of the host framework, which determines the potential surface landscape of the lithium diffusion pathway. The pathways typically include at least one bottleneck through which the lithium ion diffusion is limited. The energy barrier of the bottleneck can be lowered by introducing larger ionic radius cations on neighboring sites in order to increase the lattice size and thereby enlarge the bottleneck making more room for Li to pass through. 58 On the other hand, inducing lithium sublattice disorder (partially occupied sites) makes the lithium sites energetically equivalent reducing the defect formation energy contribution to the activation energy. Therefore, different lithium and vacancy distribution control the number of mobile charge carriers leading to higher or lower conductivity. Substitutions can therefore impact the bulk conductivity via a wide variety of mechanisms. We discuss these in detail in later sections focusing on each class of electrolyte, but it is worth illustrating a few such mechanisms here as they do occur regardless of the electrolyte class. For example, the creation or destruction of vacancies on the lithium site is a common means to tune conductivity. This is well illustrated in stoichiometric tetragonal LLZO garnet where lithium ions fully occupy the lowest energy tetragonal site making the diffusion of lithium ions very difficult but introducing vacancies, for instance, by super-valent substitution such as the partial substitution of Li + by Ga 3+ perturbs the lithium distribution and makes diffusion more facile. 59 The result is two partially occupied tetrahedral and octahedral sites, which increase the conductivity by three orders of magnitude. 60 Another common approach to tune conductivity is to make sub-valent substitution in order to add interstitials to the structure to create local distortion and flatten the energy landscape of the lithium pathway. 61 Substitutions are also often made in an attempt to improve the grain boundary ionic conductivity. In cases where the grain boundary conductivity is low, tuning the microstructure becomes essential. It is also important to recognize that the presence of a secondary phase at the grain boundaries also contributes to the total ionic conductivity (either positively or negatively depending on this phases impact on sintering properties). 62 The presence of voids between the grains drastically lowers the ionic conductivity; hence measuring the relative density of the materials is crucial for comparing the performance of ionic conductors. 63 The relative density of ISE also in great measure dictates its resistance against dendrite penetration at high current densities. 64 The relative density/porosity changes by the synthesis method, sintering method/temperatures/time, and introduction of substitutions.
Substitutions into the solid electrolyte is also often used as a means to tune the interfacial area-specific resistance (ASR) with the electrodes that include wettability and interphase stability. For example, fluorination of LLZO facilitates homogenous lithium plating and stripping cycles at different current densities as shown in Fig. 7. 65 It should be mentioned that substitutions into the bulk of the material in order to impact the surface resistance may have detrimental impacts on the bulk properties. Often times, this approach is in competition with the use of a buffer layer (i.e. a thin layer of a different solid electrolyte at the surface of the cathode/anode). For example, coating lithium cathode with lithium niobate 66 and coating the anode with alumina layer 67 can mitigate the electrolyte/electrode instability hence keeping the interfacial impedance low.  is the most promising in the lithium-stuffed garnet-type solid electrolytes. LLZO has the basis chemical formula Li 7 La 3 Zr 2 O 12 but further substitutions must be made to make this a viable solid electrolyte. LLZO has two polymorphs: tetragonal (t-LLZO, I41/acd space group) and cubic (c-LLZO, Ia % 3d space group). The host framework of t-LLZO as shown in Fig. 8(a) consists of La 3+ at two dodecahedral 8b and 16e sites and Zr 4+ at octahedral 16c site. The LaO 8 and ZrO 6 polyhedra are connected by edge-sharing while the lithium ions are ordered by fully occupying tetrahedral 8a site and two distorted octahedral 16f and 32g sites. The host framework of c-LLZO shown in Fig. 8(c) consists of La 3+ at the dodecahedral 24c sites and Zr 4+ at octahedral 16a sites. The LaO 8 and ZrO 6 polyhedra are connected by edge-sharing while the framework interstitial space is accommodated by the lithium ions and vacancies. In the c-LLZO the lithium ions only partially occupy tetragonal 24d sites and distorted octahedral 96h site as shown in Fig. 8(c). The LiO 4 tetrahedra of c-LLZO shares face with the LiO 6 distorted octahedra through which lithium ions hope (this is the bottleneck). c-LLZO has higher conductivity than t-LLZO thus, stabilizing the cubic phase is a strategy used to enhance conductivity. LLZO has a good electrochemical stability window. Its stability against the lithium anode is due to the formation of a solid electrolyte interphase (SEI), which increases the interfacial resistance slightly but also passivates the electrolyte decomposition. 49 A disadvantage of the LLZO is its air sensitivity which forms Li 2 CO 3 increasing the interfacial resistance against the electrodes. 70 4.1.2 Impact of substitutions. Substitutions on both the TM and Li sites have been widely used to enhance the ionic conductivity by stabilizing the c-LLZO and optimizing lithium concentration and distribution. The maximum amount of lithium that cubic garnet structure can have is 7.5 per unit formula with long-range vacancies ordering and short-range ordering at 6.5 lithium per unite formula thus highest conductivity can be achieves at 6.4 AE 0.1. 71 Supervalent substitution introduces vacancies, which stabilizes the c-LLZO by lowering the order-disorder transition temperature. Fig. 9(b) shows the different substitutions done on the LLZO system, while panel (a) shows all the substitutions that are expected to be stable based on computation. The colored border indicates  the different substitution sites, and the fill corresponds to the total conductivity contour. This figure clearly illustrates numerous other substitutions can still be considered for the LLZO garnet. Tables 1 and 2 summarize the results for key properties for metal substitutions into the garnet structure. These tables provide summaries of the impact of substitutions, and also dramatically illustrate properties that have to date been underexplored (e.g. the electronic conductivity in multiple substitutions has not been widely studied despite the necessity to ensure high electronic resistance to limit dendrite formation within the electrolyte). In the rest of this section, we highlight particular results from these tables and discuss their broader significance for the design of advanced solid electrolytes.

Impact of substitutions into specific classes of ISE
Partial substitution of the Li in LLZO can stabilize the c-LLZO, such substitutions include Al 3+ , Ga 3+ , Zn 2+ , Fe 3+ and these result in enhanced ionic conductivities as shown in Fig. 9(b). At high sintering temperatures, unintentional Al 3+ substitution from alumina substrate can occur. Al-LLZO has been reported with ionic conductivities as high as 3.1 Â 10 À4 S cm À1 at its solubility limit of x = 0.25 in Li 7À3x Al x La 3 Zr 2 O 12 whereas the ionic conductivity of unsubstituted LLZO is 2.0 Â 10 À7 S cm À1 . 72 Gallium substitution also stabilizes the c-LLZO by substituting at the 24d tetrahedral lithium site and producing vacancies to preserve charge neutrality. As the Ga 3+ content increases in Li 7Àx Ga x La 3 Zr 2 O 12 the lattice size and bottleneck size do not change; however, due to the Li-Ga coulombic repulsion, the vacancies around Ga 3+ occupied sites are inaccessible (trapped), causing the lithium conductivity to decrease in the range 0 r x r 0.10. On the other hand, for 0.10 o x r 0.30, the conductivity decrease trend becomes relatively flat because the extra vacancies confine percolation pathways for lithium ions motion as shown in Fig. 10. 73 In Li 7À3x Ga x La 3 Zr 2 O 12 , the highest conductivity of 1.46 Â 10 À3 S cm À1 is achieved at x = 0.25. 59 Some reports show that Ga occupies the 96h site and stabilizes the c-LLZO with I% 4% 3d space group. Accordingly, the reason the higher ionic with their corresponding total ionic conductivity (color map), and relative density (RD%). 59,63,81,99,[105][106][107][108][109][110][111][112][113][114][115][116][117][118]    conductivity of Ga-LLZO compared to Al-LLZO is due to the different site occupation where 96h site occupation does not block lithium long range path as it does the 24d site. 74 Compared to Ga 3+ substitution Zn 2+ likewise stabilizes the cubic phase. In both systems, part of the compensating vacancies are added to the tetragonal site while the rest added to the octahedral site. The amount of vacancies created by Zn 2+ is not enough to introduce vacancies to a non-vicinal octahedral site as in Ga-LLZO to facilitate lithium mobility leading to more than an order of magnitude lower conductivity. The negative aspect of Li site substitution is that it hinders the lithium path while deactivating the nearest neighboring octahedral site vacancies forming defect clusters; thus, high content of substituent at this site can block the lithium path. 75 Some substituents have a dual role of stabilizing c-LLZO and enhancing sinterability. For instance, the lower sintering temperature of 1000 1C for Ga-LLZO vs. 1100 1C for Al-LZZO and the associate higher densification of the Ga-LLZO are both attributed to the low eutectic point of Li-Ga-O compounds which is present as a liquid phase at the grain-boundaries during synthesis. 76 Dense ISEs will have enhanced ionic conductivity and they also help impede dendrite formation. Fe 3+ substituted LLZO has high ionic conductivity, but it is unstable against the lithium anode because it gets reduced to Fe 2+ yielding an electronically conductive phase that propagates through the bulk and creates a short circuit. 77 This again illustrates the critical importance of screening all key metrics in substitution studies of solid electrolytes.
Partial substitution of Zr 4+ site similarly stabilizes c-LLZO, and its advantage over Li site substitution is that it does not hinder the lithium paths. Such substitutions include Ge 4+ , Nb 5+ , Ta 5+ , Sc 3+ . A small amount of alkaline earth metal substitution around 0.05 likewise substitutes the Zr 4+ site. It expands the lattice size because alkaline earth metals have larger ionic radii than Zr, which causes the bottleneck size to increase, thus lowering the motional energy. Thus as shown in Fig. 11 larger the alkaline earth metal is higher ionic conductivity is attained. 78 Alkaline earth metal substitution also improves the sinterability by its lower melting point compared to ZrO 2 . 12 On the other hand large amount of alkaline earth metal substituents around 0.5 substitutes the La site. 79 Ta 5+ substitution of LLZO stabilizes the c-LLZO and at optimal lithium concentration of 0.6 per unit formula with the presence of a small amount of Al from the substrate distributed at the grain boundaries shows highest conductivity of 10 À3 S cm À1 . 80 When Ta-LLZO is co-substituted with Ga or Al the conductivity decreases because they occupy Li sublattice hindering their diffusion. 81 Alkaline-earth elements substitution of LLZO with contents above 0.05 substitute the La 3+ enhancing the sinterability. The density of alkaline earth element substituted LLZO increases with the ionic radius of the alkaline earth elements.
The interphase resistance of substituted LLZO with lithium anode depends on the nature of the substituents; for instance, in Al, Ta, Nb substituted LLZO, when in contact with lithium metal, Zr is reduced forming an oxygen deficient layer increasing the interfacial resistance. The amount of oxygen deficient layer forms in the substituted LLZO is as follows Ta 4+ o Nb 5+ o Al 3+ . This interphase was explored by measuring the change in interfacial impedance of Li/substituted-LLZO/Li at different time intervals and complemented by XPS measurement to explore the reduced species. Although in Nb-LLZO similar amount of oxygen deficient interphase forms but its reduction leads to the formation of a conductive phase that propagates through the bulk of ISE. While in Al substituted LLZO Zr 4+ reduces into protecting oxygen deficient layer. 49 The lithium stability of the high ionic conductive Nb-LLZO is improved by co-substituting with Ga. 82 Multiple substitutions further improves the ionic conductivity as shown in Table 1, such as substituting Li + by Ga 3+ to stabilize the cubic phase and Zr 4+ by Sc 3+ to increase the number of lithium ions and improve the density. Additionally, Sc 3+ causes local disorder of the Li network, enhancing the local mobility. 83 Multi-substitution schemes are also employed, such as Li 6.65 Ga 0.05 La 2.95 Ba 0.05 Zr 1.75 Ta 0.25 O 12 that includes supervalent substitution of Li + and Zr 4+ by Ga 4+ and Ta 5+ respectively and this stabilizes the c-LLZO by introducing vacancies besides Ga 4+ has a role of enhancing the sinterability. 84 Simultaneously, subvalent substitution La 3+ by Ba 2+ increases the lithium concentration and

Materials Advances
Review the bottleneck size due to its larger ionic radius. Additionally, multi-substituion decreases oxygen deficiency. These all enhance the lithium ion conductivity. 85 The electronic conductivities of substituted LLZO are surveyed in Table 2. The electronic conductivity of single crystal Ga-LLZO is around 10 À10 S cm À1 while for the polycrystalline material it is 10 À8 S cm À1 . The higher electronic conductivity in the polycrystalline sample is due the grain boundary and not of the intrinsic bulk, which implies that researchers must be mindful of the electronic conductivity in the final device and not solely in a perfect crystal. 86 Further studies demonstrate that coating Ca-Nb-LLZO with LiAlO 2 before sintering lowers the electronic conductivity from 3.59 to 1.01 Â 10 À8 S cm À1 . 87 This illustrates (as in part I of this review) that coatings can often be used to enhance surface properties without excess substitution into the bulk of the materials, such that coating methods are often in competition with bulk substitutions in real applications. An important current drawback of LLZO is that Li 2 CO 3 forms when LLZO is exposed to air and this reduces the ionic conductivity. 88 Though the focus of this review is the impact of substitutions, it is worth briefly summarizing here the state-of-the-art for each class of ISE in order to consider what future role substitutions may play. The interfacial resistance between the garnet electrolyte and the electrodes currently limits the performance of garnet based all-solid-state batteries. Different strategies have been used to reduce garnet/Li interface impedance. Some examples include using molten Li-C, 89 coating garnet with C, surface acid treatment, 51 wet sanding, 50 thin layer of Ag, Sn or Ge, 90-92 solid polymer electrolyte interface, 93 and Al 2 O 3 layer deposition. 94 Similarly, the interfacial resistance between garnet and cathode is high because of the rigid nature of the garnet materials and the high cathode/garnet interfacial instability. One way to reduce the interfacial resistance is by co-sintering the composite cathode/electrolyte, but the electrolyte and electrode should not react at that temperature; otherwise, a more resistive interface can form. For instance, when Al-LLZO is co-sintered with LiCoO 2 (LCO), a resistive interface forms due to Al diffusion at high temperature forming t-LLZO at the interface. 95 This can be avoided by using Ta-LLZO instead. A discharge capacity of 110 mA h g À1 with 50 mA cm À2 current density at 50 1C was achieved in a LiCoO 2 / Ta-LLZO/Li cell configuration with a composite cathode composed of LiCoO 2 and Ta-LLZO co-sintered at 1050 1C. This battery showed good performance for 100 cycles as shown in the Fig. 12(a). 96 Although the areal discharge capacity dropped significantly at higher discharge current densities as shown in Fig. 12(b) but it performed better than other ASSB cell configurations at these current densities. Other reagents can be added to composite cathode to form an interfacial layer between the cathode active materials and the solid electrolyte during co-sintering. For example, formation of either a spinel layer on a layered oxide cathode 97 or Li 2.3Àx C 0.7+x B 0.3Àx O 3 at the LLZO/LCO interface 98 reduces the interfacial resistance. Low melting point sintering aids, such as Li 3 BO 3 , can also be added to lower the required sintering temperature and thereby avoid elemental diffusion and reaction between electrolyte/electrode that can occur at high temperatures. Another method such as depositing a thin film of the cathode active material on top of the electrolyte directly 99,100 or with interface layer within between, such as a thin layer of Nb is also explored but this method has limited capacity. 101 Despite all these developments garnet based ASSLBs continue to struggle to achieve high rate capacity and long cycle life. Continued research and development is required to further improve these ASSLB.

Perovskite
Perovskite materials have the generic formula ABX 3 and are made of corner shared BX 6 octahedra and A located at dodecahedral sites as shown in Fig. 13. Lithium lanthanum titanates (LLTO) with composition Li 3x La (2/3)Àx TiO 3 (1/25 o x o 1/6) are the most promising electrolytes in the perovskite family because of their high bulk conductivities reaching B10 À3 S cm À1 . 140 Lithium and lanthanum share occupation of the A sites while titanium fully occupies the B sites. The ionic conductivity of perovskite ISE is strongly dominated by the grain boundary conductivity that is very poor B10 À4 -10 À5 S cm À1 ; however, substitutions can reduce the substantial grain boundary resistance. Table 3 summarizes the impact of substitutions made into LLTO and the resulting effect on total conductivity. The lattice size changes by introducing different size substituents to the A-site such as subvalent partial substitution of La 3+ by larger ionic radius Sr 2+ expands the lattice size leading to lower barrier energy and increases the lithium concentration by charge compensation. These two factors enhance the bulk conductivity. 141 Likewise, partial substitution on the b-site affects the ionic conductivity. Similar to the other ISE families, substituting Ti 4+ with Ge 4+ increases the density of the solid electrolyte, which in turn improves the total conductivity by improving the transport at the grain boundaries. 142 On the other hand, although partial substitution of Ti 4+ by Nb 5+ lowers the bulk conductivity, it also yields a higher grain boundary conductivity, which is the limiting factor of LLTO total conductivity such that Nb-substitution does improve the overall conductivity. 143 A disadvantage of LLTO electrolyte is that it gets lithiated by the metallic lithium anode converting into conductive phase thereby preventing its integration into ASSLB with Li metal anode.

NASICON-type
NASICON-type ISEs are lithium ion conductors with the structure of a sodium super ionic conductor with chemical formula LiM 2 (XO 4 ) 3 . The structure, shown in Fig. 14, is rhombohedral and can be described in the R% 3c space group. It is made of XO 4 tetrahedra at 18e site and MO 6 octahedra at 12c sites joined by corner sharing while the lithium ions fully occupy antiprismatic 6b sites (Li1). LiM 2 (XO 4 ) 3 has many polymorphs, but the rhombohedral structure yields the highest ionic conductivity. Fig. 15 shows the various substitutions attempted to improve this class of electrolytes. These substitutions are on two sites: the X sites (P, Si) and the M sites (Ge, Ti, Sc, Zr, Hf). For the M sites, the activation energy decreases as the ionic radius increases due to the larger lattice with expanded bottleneck size as shown in Fig. 16. However, when we resort to larger radii cations such as Zr, Hf, and Sc, the expansion of the lattice causes a structural distortion wherein Li ions move to a more stable site (M2) and transforms the structure into monoclinic/ triniclic phases that are poorly conductive. The optimum M substitution in terms of ionic conductivity is therefore considered to be Ge as illustrated in Fig. 15. Three classes of NASICON-type materials have been studied with parent structures with M = Ti (LTP), Ge (LGP), and Zr (LZP). The rest of this section discusses progress in these families. Although LiGe 2 (PO 4 ) 3 LGP can be synthesized in the rhombohedral structure without any substituents, the conductivity is quite small (B10 À6 S cm À1 ). All the Li ions occupy the 6b site with the 18e site vacant. Fig. 15 shows that subvalent partial substitutions of Ge 4+ with Cr 3+ and Al 3+ dramatically improves the conductivity. This is attributed to lithium ions now occupying both the 6b and 18e sites. Subvalent substitution of LGP by Al 3+ introduces additional lithium ions to the Li3 site and due  to the repulsion with the Li1 site, the lithium ions redistribute between Li1 and Li3. This repulsion causes the volume of the LiO 6 octahedra to increase, expanding c of the lattice, but the size of the oxygen triangle windows T1 and T2 shown in Fig. 14 does not expand. The coulombic interaction between Li1 and Li3 pushes Li1 into Li3*, which acts as a spring flattening Li ion energy landscape, decreasing the activation energy. Therefore the mechanism behind the conductivity enhancement of Al- LGP is that by increasing the energy of lithium, the local minima become less profound, thereby decreasing the migration barriers. 146 Optimal Al 3+ substitution content of 0.4 mol per unit formula of LGP achieves room temperature total conductivity of 1.22 Â 10 À3 S cm À1 which is three order of magnitude higher than of unsubstituted LGP and is certainly of interest in real battery applications. 147 Higher substituent content increases activation energy due to the smaller bottleneck size. 148 Further co-substituting with Y tweaks the microstructure and grain morphology producing a denser product with lower intergrain resistance. 149 Other co-substitution explored on Al-LGP is Sr where due to its large ionic radius has limited solubility of 0.17 per unit formula above which secondary phase (ScPO 4 ) forms at the grain boundary lowering intergrain lithium diffusion. Higher Sr contents also distorts the cubic phase as mentioned above. The optimal composition of Li 1.5 Al 0.33 Sc 0.17 Ge 1.5 (PO 4 ) 3 shows higher conductivity than LGP due to the larger Li diffusion channel size and facilitating new diffusion channels direct 6b-36f path through which the Li ions move in concerted fashion instead moving through 18e by isolated single-ion diffusion. 150 One drawback of LGP is that it gets reduced when in contact with metallic lithium anodes yielding an electronically conductive material. 151 LiZr 2 (PO 4 ) 3 (LZP) has four polymorphs which are presented in Fig. 17. At synthesis temperature above 1100 1C, a 0 triclinic  phase (P% 1) forms which transforms into a rhombohedral phase (R% 3c) at around 40 1C. Whereas the b 0 monoclinic phase (P21/n) forms at synthesis temperatures below 1100 1C and transforms into b orthorhombic phase (Pbna) at around 300 1C. Both a and b has ionic conductivity but a phase has higher. Thus, lowering the transition temperature to stabilize room temperature a phase is a promising strategy to achieve high conductivity. Such stabilization of the R% 3c phase is accomplished by subvalent substitution of LZP with La 3+ such substitutions are shown in Fig. 15. Moreover, the lithium concentration is increased by neutrality preservation and due to its larger ionic radius, the lattice size and triangle bottleneck window size expands leading to lower migration energy. The sum effect is enhanced bulk conductivity. It worth mentioning that different synthesis methods can also lower the phase transition temperature. For example, the sol-gel method can stabilize the rhombohedral without any substitution by reducing the a 0 to a LZP phase transition temperature from 40 to 15 1C, yielding total and bulk conductivities of 1.5 Â 10 À6 and 1.5 Â 10 À4 S cm À1 , respectively. 152 LTP conductivity is enhanced by subvalent substitution as shown in Fig. 15 such as with Al 3+ , Cr 3+ . 153 The mechanism of the ionic conductivity enhancement is the same as in the subvalent substitution of LGP. The lithium in Al-LTP moves through the Li1-Li3-Li3-Li1 pathway which is explored by maximum-entropy method (MEM) 154 Introducing Al to LTP in the range of 0.1-0.5 mol per unit formula increases the bulk conductivity three orders of magnitude to reach B10 À3 S cm À1 which was measured on a single crystal. 155 Low temperature EIS was used to separate the bulk and grain conductivities in polycrystalline samples of LATP as shown in the Arrhenius plots in Fig. 18. This reveals poor grain boundary conductivity with a high activation energy. Furthermore, Fig. 18 also shows that LATP has a decreasing bulk conductivity when exposed to air and the process can be reversed with heating. 44 Ti 4+ -Based NASICON-type LTP gets reduced by Li metal into an electronically conductive phase. Thus, a buffer layer is currently required (e.g. LIPON layer 156 or polymer electrolyte 157 ) in order to protect the LATP.
Still of interest for further study, this class of material may yet be optimized with further substitutions. It should however be noted that these materials are extremely rigid such that assembly in a ASSLB is difficult and maintaining the battery structure during operation is very difficult unless low volume expansion electrodes are used and this greatly limits the energy density. 158

LISICON
Lithium super ionic conductor (LISICON) is a family of materials with the structure of Li 4 SiO 4 or g-Li 3 PO 4 with monoclinic and orthorhombic structures, respectively, as shown in Fig. 19. Both are made of XO 4 tetrahedra (X = P, Si), but with different orientations. The three lithium ions per unit formula of g-Li 3 PO 4 are located at the tetrahedral 8d and 4c sites. The XO 4 tetrahedral site can be partially substituted by different valence cations changing the charge carrier concentration and lattice size, hence tuning the ionic conductivity. Table 4 shows the resulting conductivities for the various substitutions attempted into these structures. Here, we focus on the one material that has a conductivity about 5 Â 10 À5 S cm À1 (W substituted Li-Ge-O) and refer a reader interested in this class of materials to a review article explores this class more. 175 Li 3.70 Ge 0.85 W 0.15 O 4 has a dramatically improved conductivity compared to the parent Li 4 GeO 4 due to the supervalent W 6+ substitution creating Li vacancies, and enlarging the oxygen windows such that the Li-diffusion energy barrier is lower. 176 The drawback of LISICON is that they have high conductivity only at high temperatures; for instance Li 3.5 Si 0.5 P 0.5 O 4 has a room temperature conductivity of B10 À7 S cm À1 and B10 À3 S cm À1  at 300 1C. At low temperature ordering of the interstitial lithium ions decreases the conductivity. However, at high temperature the two octahedral interstitial sites become partially occupied with weaker interaction with the oxygen at high temperature that lowers the activation energy. The ionic conductivity of LICISON is modest and still needs further enhancement.

Thio-LISICON
The anionic substitution of O 2À by more polarizable S 2À in the LISICON materials discussed in the previous section leads to orders of magnitude higher ionic conductivity ISE called thio-LISICON. The structure of thio-LISICON is the same structure as the LISICON only O 2À is replaced by a more polarizable S 2À anion. Li 4 GeS 4 has b-Li 3 PS 4 thio-LISICON structure and can be represented in the Pnma space group. When Ge is substituted by Sn forms a solid solution Li 4 Ge 1Àx Sn x S 4 . Although the introduction of larger ionic radius Sn increases the triangular bottleneck window but does not affect the ionic conductivity because it enhances only short-range diffusion. However, the Li 4 SnS 4 end of the solid solution shows an order of magnitude higher ionic conductivity 1.4 Â 10 À6 S cm À1 due to the enhanced short range diffusion combined with the reposition/redistribution of lithium ions which decreases the Li1-Li2 distance shown in Fig. 20(c) converting their configuration from edge shared to face shared facilitating long range diffusion chain along the b-axis. Moreover, the introduction of new partially occupied Li4 sites enables 3D conduction by connecting the linear diffusion paths as shown in Fig. 20(d). 185 The low ionic conductivity of Li 4 SnS 4 can be enhanced by substituting it with As forming Li 4Àx Sn 1Àx As x S 4 (x = 0 to 0.250) solid solution with the highest ionic conductivity being 1.39 Â 10 À3 S cm À1 for Li 3.833 Sn 0.833 As 0.166 S 4 . 186 The advantage of the thio-LISICON mentioned above compared to thiophosphate based ISEs that are discussed below is their air stability because thiophosphates react with moister to release H 2 S gas. The ionic conductivity of parent b-Li 3 PS 4 increases 3 order of magnitude by substituting P with Si. It forms a very limited solid solution Li 3+x [Si x P 1Àx ]S 4 (0.15 o x o 0.33) where the b-Li 3 PS 4 structure is entropy stabilized by splitting lithium site causing disorder which flattens the energy landscape thus reaching ionic conductivity of 1.22 Â 10 À3 S cm À1 for Li 3.25 [Si 0.25 P 0.75 ]S 4 . 187 Furthermore, the Li 4 GeS 4 -Li 3 PS 4 solid solution is divided into three regions (I, II, and III), as shown in Fig. 21, according to their superstructures caused by the different cation ordering. In this solid solution, Li 3.25 Ge 0.25 P 0.75 S 4 has the highest conductivity. The Li 2 GeS 4 -Li 2 PS 4 phase diagram is shown in Fig. 22, which reveals the regions where each phase the thio-LISICON and LGPS (discussed separately later in this review) exists. The downside of Li 3.25 Ge 0.25 P 0.75 S 4 is its instability against lithium. Substitution of Ge by Al into Li (4À1/3) Al 1/3 P 2/3 S 4 enhances stability against lithium but lowers the conductivity to 8.02 Â 10 À4 S cm À1 at 25 1C. 188 This impressive class of solid electrolytes is receiving a great deal of interest in all solid batteries, so we now briefly summarize the state-of-the-art in thio-LISICON ASSLB. Alloy anodes are widely used to integrate sulfide electrolytes in all-solid-state   batteries since they are unstable against lithium anode.
Although the potential of alloy anodes still are outside the electrochemical stability window of sulfide electrolytes, its reduction thermodynamic driving force is lower compared to lithium anode. Such alloying metals are Al, Ga, In, Sn, or Sb; from this list In is the most promising because of its ductility. 189 Moreover, sulfide-based electrolytes have a low upper electrochemical stability window such that a buffer layer between the cathode and electrolyte is needed. On the contrary, it is compatible with the sulfur cathode, making it a good candidate for Li-S batteries. Recently, an impressive areal capacity of 7.8 mA h cm À2 was obtained for a Li-S cell utilizing Li 3 PS 4 as the electrolyte (this is well above the roughly 3 mA h cm À2 achieved in cathodes for state-of-the-art Li-ion batteries discussed in part I of this review). It should be noted that this performance was obtained at a slow cycling of about C/65 (the current density was 0.12 mA cm À2 ) and the capacity faded rapidly within a few cycles, indicating that much work is needed to further improve the performance of all-solid batteries. 190 In Fig. 23(a) the electrochemical performance of 7.7 mg cm À2 loading at different cycles is show. At higher loadings as shown in Fig. 23(c) the gravimetric capacity dropped due the higher polarization of thick electrode. Although as shown in Fig. 23(b) the coulombic efficiency of high loading cell is almost 100% for the first 10 cycles but the charge discharge capacity dropped. Furthermore, the cell configuration had an average voltage in the range 1.2-1.5 V, thereby making their energy about 80% that of current Li-ion batteries operating at average voltage of about 3.8-3.9 V (as discussed in part I of this review). Nonetheless, these full solidbattery results are extremely promising and will undoubtedly encourage a great deal of further interest in ASSLB designs with thio-LISICON electrolytes.

Argyrodite
Argyrodite has Li 7 PS 6 chemical formula with two polymorphs. Fig. 24 shows both structures: the low ionic conductivity room temperature tetragonal phase (Pna21) and the high conductivity high temperature cubic phase (F% 43m). The framework of the cubic phase is made of PS 4 tetrahedra located at the 4b site and S 2À atoms located at the 4a, 4d, and 16e site while the lithium ions are found at the interstitial sites. The cubic structure is stabilized at room temperature by anionic or cationic substitutions. Anionic partial substitution such as substituting S 2À by halide leads to Li 6 PS 5 X (X = Cl, Br, I). 193 The introduction of more polarizable halides softens the lattice decreasing the activation energy. The softer lattice of Li 6 PS 5 Br compared to Li 6 PS 5 Cl causes a decrease of the pre-exponential factor. In Br À and Cl À substitution, the structure is disordered where the 4d S 2À site is mixed occupied S 2À /X À leading to enhanced conductivity. On the other hand, Li 6 PS 5 I has very low ionic conductivity because of the ordered structure due to the difference in the ionic size of S 2À and I À . There are three possible lithium migration paths divided into intercage and intracage as shown in Fig. 24. S 2À /I À disorder facilitates the intracage transport, thus enhancing the ionic conductivity. 194 In contrast, partial subvalent cation substitution of P 5+ at the 4b site by for example, Si 4+ or Ge 4+ produces solid solution Li 7+x M x P 1Àx S 6 (M = Si, Ge) where the conductivity is improved by stabilizing the cubic phase. The ionic conductivity in these solid solutions increases to a specific substituent content, after which decreases. 195 Both Ge 4+ and Si 4+ substitution regardless of their radius cause S 2À /I À disorder and increase the lithium occupancy at both the 24g and 48h sites which lead to stronger Li-Li coulombic interactions causing an increase in the double jump path and wider jump through a trigonal planer window and intercage jump through a polyhedral window. In other words, it flattens the potential surface; hence the activation energy decreases. 196 Moreover, Si 4+ substitution expands the  electrochemical stability window to reach 6 V. 197 Although Li 6 PS 5 I has low conductivity as mentioned above, iodine based interface helps to form uniform plating and stripping against the lithium anode threby reducing dendrite formation.
Cationic substitution of P 5+ by Sn 4+ in Li 7 PS 6 is not possible but becomes possible in Li 6 PS 5 I since it has a larger lattice. The larger lattice size caused by larger ionic radius of Sn 4+ compared to P 5+ and higher lithium content in Li 6+x P 1Àx Sn x S 5 I enhances conductivity. Furthermore its air stability improved because Sn is soft acid and prefers to bond to soft base S 2À in contrast P 5+ is prone to oxidation where PS 4 reacts with water forming H 2 S. 198 Analogous to Li 6+x M x P 1Àx S 5 I (M = Si, Ge, Sn), thioantimonate argyrodite Li 6+x M x Sb 1Àx S 5 I (M = Si, Ge, Sn) yields an impressive ionic conductivity of 2.4 Â 10 À2 S cm À1 for Li 6.6 Si 0.6 Sb 0.4 S 5 I which is among the highest reported solid electrolyte ionic conductivities. This high ionic conductivity is due to minor S 2À /I À disorder and significant lithium cation disorder. 61

LGPS-like
Li 10 GeP 2 S 12 (LGPS) structure as shown in Fig. 25(b) consists of two tetrahedral sites where the 4d site is mixed occupied by Ge/ P and the 2b site is fully occupied by P. There are four lithium  sites two of each is fully occupied and the other two 16h and 8f sites are partially occupied making a chain along c-axis. The fully occupied lithium octahedral 4d site connects the XS 4 tetrahedra along the a-axis through edge sharing with the tetrahedra at 4d site and corner sharing with PS 4 . While the other 4c fully occupied lithium site connects the PS 4 tetrahedra along the a-axis through edge sharing with tetrahedra at 4d site and corner sharing with XS 4 . The lithium transport pathway is through the chain of interstitials at the 16h and 8f sites along the a-axis. 69 These 1D lithium diffusion paths are interconnected by 8f-4d-8f and 16h-4c-16h lithium paths enabling 3D lithium diffusion as shown in Fig. 25(c). 56 Although substituting of Ge 4+ by Sn 4+ expands the unit cell, the bottleneck size S3-S2 shown Fig. 26 becomes smaller, leading to higher migration energy. 202 Moreover, this substitution causes the lattice to soften and the coulombic interaction between S 2À and Li + to increase because of the longer Sn-Li bond and lower electronegativity making the S more electron dense which all likewise increases the activation energy. 11 Table 5 summarizes the substitutions in LGPS like structure. The drawback of LGPS family ISE their very narrow electrochemical stability window because sulfur and phosphorus are prone to oxidation and reduction, respectively. When LGPS is used with metal anode, it degrades into Li 2 S Li 3 P and Ge-Li alloy leading to interfacial resistance growth. The ionic conductivity of Li 2 S dictates this resistance growth. The decomposition rate of LGPS against lithium anode extrapolates to a 4.6 kO cm 2 interfacial resistance within a year increasing the cell overpotential (internal resistance). 203 Its narrow electrochemical stability window can be resolved by adding buffer layer such as coating the LCO cathode with LiNbO 3 and in situ formation of LiH 2 PO 3 on the lithium anode. 204 Substitution in LGPS has been investigated computationally for Li 10AE1 MP 2 X 12 (M: Ge, Si, Sn, Al or P, and X: O, S or Se). 205 The anionic substitution of S 2À by O 2À decreases the conductivity because of its smaller ionic radius and polarizability compared to sulfur. This substitution narrows the lithium diffusion path, however, moving to Se doesn't enhance the conductivity further which means an optimal channel size is obtained with S 2À . By contrast, the effect of cationic substitution on the ionic conductivity is small, as shown in the    Table 6. 205 As an effort to enhance the narrow electrochemical stability window of LGPS, partial substitutions such as Li + /Ba 2+ and S 2À /O 2À have been performed, where the ionic conductivity remained high and a higher electrochemical stability was reported. 206,207 These electrochemical stability windows are overestimated because they were measured on a Li/SE/Au cell that is flawed as discussed earlier, but the fact that the substitutions yield improvements is sound.
The state-of-the-art for LGPS ASSLB is impressive. Fig. 25(a) shows the Arrhenius plot of LGPS-like structure ISE with different substitutions/compositions accomplished in the literature. The optimal chemical formula Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3 with LGPS structure holds the record with an ionic conductivity of 2.5 Â 10 À2 S cm À1 . 208 In an all solid battery with the LiCoO 2 + LGPS/LGPS/Li configuration, the instability of LGPS against the cathode and anode was resolved by buffer layers on each electrode, thereby achieving a capacity of 113.7 mA h g À1 at the 500th cycle at a cycling rate of 0.1C. 204 Despite the challenges with air stability, this class of ISEs will continue to receive a great deal of attention due to the unmatched ionic conductivities.

Anti-perovskite
While all the previously discussed classes of structures have been studied for a significant period of time that has allowed time for significant exploration of the impact of substitution, anti-perovskites represent a newer class of ISEs where the impact of substitutions has not been explored so thoroughly. Anti-perovskite has the same structure as perovskite with inverted charges where oxygen occupies the octahedral B site OLi 6 with the lithium ions at the vertices, and the dodecahedral A site is occupied by a halide or hydroxide with conductivity reaching around B10 À3 S cm À1 . The performance of lithium rich anti-perovskites Li 3 OA is enhanced by tuning the content of A site. The conductivity enhancement moving from I À to Cl À results from the optimal radius size ratio of Cl À to O and Li expanding the Li ion transport channel reaching room temperature conductivity of 8.5 Â 10 À4 S cm À1 . 210 Partial substitution of the Li can tune the charge carriers concentration. Perfect lithium anti-perovskite crystal is not conducive due to the absence of lithium vacancies. The presence of point defects Frenkel or Schottky (Li A deficiency) improves its conductivity. The conductivity of anti-perovskite can be enhanced by anionic mixing and cationic partial substitution. 141 Partial substitution of Li by divalent cation such as Ba 2+ incorporates vacancies into the structure, improving the performance of the glassy ISE. 211 Hydrogen substitution to obtain Li 3Àx OH x Cl lowers energy barriers by facilitating Schottky defects formation. 212 A site mixing can generate lattice mismatch which improves lithium diffusion for instance, Li 3 OCl 1Àx Br x has conductivity higher than the end members. 213 One the other B site partial substitution of O 2À by S 2À produces double perovskite with lower defect formation energy, hence improving room temperature conductivity. 214 It worth mentioning that in all these antiperovskite structures, slight off stoichiometry can produce vacancy facilitating the lithium transport through vacancy hopping route. The ordered vacancy in Li 2 OHCl orthorhombic structure at high temperature converts into cubic phase shown in Fig. 27(a) with disordered vacancy enhancing the conductivity. The partial substitution of OH À by F À stabilizes the cubic phase at room temperature and reduces the number of À OH groups which hydrogen hinters lithium motion; thus, the conductivity increases. 215 On the other hand the substitutions that introduce octahedra tilting or and distortion as shown in Fig. 27(b) decreases the ionic conductivity. 216 The disadvantage Table 6 The ionic activation energy and conductivity calculated computationally for different possible cationic and anionic substitutions in LGPS structure. 205 Reproduced with permission from ref. 205 of anti-perovskite ISE is their low oxidation potential which was predicted by DFT calculation to be between 2.5-3.0 V, 217,218 which currently limits their application to low voltage battery. We consider this class of ISE to be significantly improved by further substitution studies, with a particular focus on the stability against the cathode.

Design principles for optimizing solid electrolytes
The impact of metal substitutions into inorganic solid electrolytes on electrochemical properties and the underlying mechanism were reviewed. Substitutions in ISEs can enhance bulk ionic conductivity by increasing the Li carrier concentration and decreasing the bottleneck energy. There are a number of important design principles that were discussed throughout this review. These can be summarized as follows.
In order to use metal substitutions to enhance the total ionic conductivity, this can be done either by tuning: (1) The bulk conductivity (structure tuning): Tuning lattice volume by isovalent substitution to lower the bottleneck (barrier) energy.
Induce lattice distortion. Increase Li sublattice disorder through aliovalent substitution by creating Li vacancies.
Stabilize the high-temperature ion conductive phase by substitution.
Create defects by aliovalent substitution through charge compensation via introducing or increasing the concentration of the vacancies or interstitial Li ions.
Tweak the potential energy surface of the Li ion by inserting Li ions into high-energy sites to activate concerted ion migration through cation substitution.
The inductive property of the cation substituent can change the bond energy of Li-anion, thus lithium energy potential landscape.
(2) The grain boundary contribution to the total conductivity (microstructure tuning): Sintering additive to increase the density of the solid electrolyte, thus enhancing the grain boundary ionic conductivity.
Engineer the orientation of the domains and control the size of the grains.
It is now essential to look at the impact of metal substitutions on properties other than ionic conductivity. In particular, it is vital to obtain accurate values for electronic conductivity, and stability at both high and low electrochemical potential. We note in particular a lack of electronic conductivity measurements for the substituted ISE in many cases and encourage the ISE research community to measure this important metric systematically. Similarly, the use of a Li/electrolyte/inert metal cell measured by the cyclic voltammetry method is deemed to be an unreliable method of measuring the electrochemical stability window as it has been demonstrated to give far better stability results than in real battery operation (and computational studies). The systematic reporting of all important metrics for solid electrolytes is therefore deemed critical so that metal substitution studies continue to play an important role in developing full-solid batteries.

Current challenges and future perspectives
The ionic conductivities of a number of ISE electrolytes have reached values close to and even higher than those of liquid electrolytes. However, their performance in full batteries remains far behind the conventional lithium-ion battery due to both the intrinsic property limitations and its interface nature with the electrodes. For example, thiophosphate based electrolytes (argyrodite, LGPS-like, thio-LISICON) have good mechanical properties because of their formability, but their stability against air (O 2 and H 2 O) must be improved. Moreover, their electrochemical stability window is very narrow and they are chemically unstable against lithium. Further studies on the impact of metal substitutions must focus more strongly on these properties. On the other hand, oxide base electrolytes are more stable in air. However, their brittle nature makes them mechanically incompatible for direct contact with the electrode since they cannot accommodate the expansion/contraction of the electrode. Additionally, their synthesis requires high temperatures. It is unlikely that metal substitutions can play a role in improving the mechanical properties, instead composite electrolytes are expected to be required.
Importantly, despite the high ionic conductivity of the present solid-state electrolytes, the electrolyte/electrode interface behaves as a bottleneck, increasing the internal resistance and lowering the critical current density. Most of the all-solid batteries currently operate at low current densities and have low areal capacity compared to liquid-based lithium-ion ones. The challenge lies in understanding these interfaces better to improve them. Results to date suggest that innovative architectures may be required to overcome current limitations. Engineering better interfaces with both improved physical contact (wettability and loss of contact by progressive cycles) and chemical/electrochemical stability of the electrode/electrolyte interface will be required to achieve higher areal capacities at higher C-rates. Metal substitutions should continue to play an important role here, especially in improving the stability of the electrolyte, as long as substitution studies start to place a higher importance on reliably screening these properties as highlighted many times in this review.

Conflicts of interest
There are no conflicts to declare.