Patent Publication Number: US-2022223906-A1

Title: Lithium and sodium superionic conductors

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
     This application is a divisional of U.S. application Ser. No. 16/346,144, filed Apr. 30, 2019 entitled “LITHIUM AND SODIUM SUPERIONIC CONDUCTORS” which is a National Stage of PCT/US17/59340, filed Oct. 31, 2017 which claims the benefit of priority of U.S. Provisional Patent Application Ser. No. 62/415,167, filed Oct. 31, 2016, entitled “NOVEL LITHIUM SUPERIONIC CONDUCTOR SOLID ELECTROLYTES”, and U.S. Provisional Patent Application Ser. No. 62/547,478, filed Aug. 18, 2017, entitled “NOVEL SODIUM SUPERIONIC CONDUCTOR SOLID ELECTROLYTE” both owned by the assignee of the present application and herein incorporated by reference in their entirety. 
    
    
     GOVERNMENT FUNDING 
     This invention was made with government support under DE-SC0012118 awarded by the Department of Energy and under CMMI1436976 awarded by the National Science Foundation The government has certain rights in the invention 
    
    
     BACKGROUND 
     All-solid-state rechargeable lithium-ion batteries (ss-LIBs) are a revolutionary architecture for energy storage that holds the promise to be both more energy dense as well as safer than traditional organic-liquid-electrolyte-based lithium-ion batteries. The crucial enabling component in ss-LIBs is the non-flammable lithium superionic conductor solid electrolyte, which must possess a challenging suite of properties, including high Li +  conductivity, excellent phase stability, electrochemical and mechanical compatibility with the electrodes, and ideally, stability under ambient air. 
     Thiophosphate lithium superionic conductors have emerged as one of the leading contenders for ss-LIB solid electrolyte applications. For instance, the Li 10 GeP 2 S 12  (LGPS) family of materials and the Li 7 P 3 S 11  glass-ceramic have ionic conductivities exceeding 10 mS/cm, on par or even exceeding that of organic solvent electrolytes. More recently, Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3 , which also has the LGPS structure, was reported to have an ionic conductivity as high as 25 mS/cm at room temperature. Besides typically higher ionic conductivities, sulfides have the advantage of being softer than oxides, allowing low porosity and intimate electrode-electrolyte contact to be achieved via cold-pressing techniques. 
     Despite these advances, the number of known lithium superionic conductors remains relatively few, and existing materials all suffer from various limitations. For instance, Li 10 GeP 2 S 12  is expensive due to the use of Ge, while the Sn and Si substituted analogues tend to form conducting, non-passivating phases upon contact with the Li anode. Other candidates, such as Li 7 P 3 S 11 , are metastable and tends to form non-passivating layers with the typical LiCoO 2  cathode used in today&#39;s LIBs. There is thus a need for the development of new lithium superionic conductors that can potentially mitigate some of these limitations. 
     SUMMARY 
     Systems and methods according to present principles meet the needs of the above in several ways. 
     Presented are new, earth-abundant lithium superionic conductors, Li 3 Y(PS 4 ) 2  and Li 5 PS 4 Cl 2 , that emerged from a comprehensive screening of the Li—P—S and Li-M-P—S chemical spaces. Both candidates are derived from the relatively unexplored quaternary silver thiophosphates. One key enabler of this discovery is the development of a first-of-its-kind high-throughput first principles screening approach that can exclude candidates unlikely to satisfy the stringent Li +  conductivity requirements using a minimum of computational resources. Both candidates are predicted to be synthesizable, and are electronically insulating. Systems and methods according to present principles enable new, all-solid-state rechargeable lithium-ion batteries. 
     Li 3 Y(PS 4 ) 2  may in some cases have superior properties. Compared to the existing lithium superionic conductors, Li 3 Y(PS 4 ) 2 , the more promising among the two, is expected to present an overall better balance of properties as a lithium superionic conductor electrolyte for all-solid-state battery applications. It has better predicted phase and electrochemical stability. Its Li +  conductivity is sufficiently high (&gt;2 mS/cm) such that it is not likely a limiting factor, and can be increased multi-fold to ˜7 mS/cm via aliovalent doping strategies. 
     Regarding the new screening approach for superionic conductor electrolytes, its validity has been confirmed and the same has been shown to efficiently extract promising candidates with a minimal cost of computational resources. 
     In one aspect, the invention is directed towards a high-throughput screening method for superionic conductor, comprising: determining an initial pool of candidate structures from existing Li—P—S, Ag—P—S ternary and Li-M-P—S, Ag-M-P—S quaternary structures; filtering out unstable materials; and performing diffusivity screening on remaining materials. 
     Implementations of the invention may include one or more of the following. The filtering may be performed by phase stability analysis. The diffusivity screening may be performed by a three step approach. The three steps may include topological analysis, quick diffusivity estimation, and long AEID simulations. The long AIMD simulations may be performed at multiple temperatures for a converged diffusivity of the most promising candidates. The method may further include performing dopant and composition optimization. 
     In another aspect, the invention is directed to a superionic conductor, found by the method described above. 
     In yet another aspect, the invention is directed to a superionic conductor, having the structures of Li 3 Y(PS 4 ) 2  or Li 5 PS 4 C 12 . 
     Advantages of certain implementations of the invention may include one or more of the following. Li superionic conductors as solid electrolytes have several advantages, including overcoming the potential leakage and flammability problems arising from use of traditional organic liquid electrolytes. The lithium superionic conductors according to present principles further exhibit high conductivity (&gt;2 mS/cm), good phase and electrochemical stability and low electronic conductivity, which are highly promising to boost the electrochemical performance of all-solid-state battery applications. 
     In another aspect, the new screening approach for superionic conductor electrolytes has been employed to a predict a new Na superionic conductor (Na 3 Y(PS 4 ) 2 ), which has excellent phase stability at 0 K and high conductivity (&gt;10 mS/cm). It is the first sodium conductor from the computational or experimental side that exhibits both good phase stability and very high conductivity. 
     Other advantages will be understood from the description that follows, including the figures and claims. 
     This Summary is provided to introduce a selection of concepts in a simplified form. The concepts are further described in the Detailed Description section. Elements or steps other than those described in this Summary are possible, and no element or step is necessarily required. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended for use as an aid in determining the scope of the claimed subject matter. The claimed subject matter is not limited to implementations that solve any or all disadvantages noted in any part of this disclosure. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIG. 1  is a flowchart of the screening procedure for new lithium superionic conductors, where the relative sizes of the arrows indicate the number of candidates remaining after the current screening step. 
         FIG. 2  is a plot of MSD 1200K  versus MSD 800K  for a wide range of known alkali conductors and the new candidates described herein. 
         FIG. 3 a    shows the crystal structure (conventional cell) of Li 3 Y(PS 4 ) 2  and  FIG. 3 b    shows the crystal structure (conventional cell) of Li 5 PS 4 C 12 . 
         FIG. 4  shows Arrhenius plots of promising lithium superionic conductor candidates from AEID simulations, where solid lines indicate the initial Li 3 Y(PS 4 ) 2  and Li 5 PS 4 Cl 2  candidates and dashed lines indicate the further optimized candidates obtained via aliovalent doping or substitution of Li 3 Y(PS 4 ) 2 , (Li 3 La(PS 4 ) 2 , Li 3.125 Y 0.875 Ca 0.125 (PS 4 ) 2  and Li 20.875 Y 0.875 Zr 0.125 (PS 4 ) 2 ). 
         FIG. 5A  and  FIG. 5B  are isosurfaces of Li +  probability density distribution P for Li 3 Y(PS 4 ) 2  and Li 5 PS 4 Cl 2  from AEID simulations at 800K with P=0.0001 a 0   −3  (a 0  is the Bohr radius). 
         FIG. 6( a )  shows the investigated Li vacancy diffusion paths in Li 3 Y(PS 4 ) 2  viewed along the b direction and  FIGS. 6( b )-6( d )  show the calculated CI-NEB migration barriers for selected percolating paths. 
         FIG. 7 a    shows the investigated Li vacancy diffusion paths in Li 5 PS 4 Cl 2  viewed along the a direction and  FIG. 7 b    shows the investigated Li vacancy diffusion paths along the b direction and  FIGS. 7( c )-7( e )  show calculated CI-NEB migration barriers for selected paths. 
         FIG. 8 a    shows the calculated element-projected density of states for Li 3 Y(PS 4 ) 2  (band gap Eg=3.41 eV) and  FIG. 8( b )  shows the calculated element-projected density of states for Li 5 PS 4 Cl 2  (Eg=3.57 eV) using the HSE screened hybrid functional. 
         FIG. 9  is a schematic diagram of an all-solid-state lithium battery 
         FIG. 10  is a flowchart of the screening procedure for new sodium superionic conductors, 
         FIG. 11  is schematic diagram of an all-solid-state Na-ion battery. 
         FIG. 12  shows the computed Na 2 S—P 2 S 5 —Y 2 S 3  phase diagram. 
         FIG. 13  shows XRD patterns of Ag 3 Y(PS 4 ) 2 . 
         FIG. 14  shows preliminary XRD result of Na 3 Y(PS 4 ) 2 . 
     
    
    
     DETAILED DESCRIPTION 
     Introduction 
     When surveying the space of known lithium thiophosphates, an interesting observation is that many of them have analogues in Ag thiophosphates. For example, Li 7 P 3 S 11  and Li 3 PS 4  bear remarkable structural similarity to Ag 7 P 3 S 11  and Ag 3 PS 4  respectively. The highly interesting Li argyrodite superionic conductors with formula Li 6 PS 5 X (X═Cl, Br, I) derive their name from the mineral argyrodite (Ag 8 GeS 6 ), and show promising Li conductivities of &gt;1 mS/cm for X═Cl and Br. Many Ag compounds are also known to exhibit extraordinarily high ionic conductivities; for example, α-AgI is perhaps the best known, and one of the first superionic conductors ever discovered. 
     Inspired by this observation, we have performed a comprehensive screening of the ternary Li—P—S and quaternary Li-M-P—S(where M is a non-redox active element) chemical spaces for new lithium superionic conductors using an efficient screening approach based on high-throughput density functional theory (DFT) calculations. The scope of this work extends beyond the known Li thiophosphates and includes new candidates obtained from Ag for Li substitution of Ag thiophosphates. The screening yielded two highly promising candidates, Li 3 Y(PS 4 ) 2  and Li 5 PS 4 Cl 2 , which are predicted to satisfy the necessary combination of excellent phase and electrochemical stability, high Li +  conductivity, and low electronic conductivity. We also show that the conductivity of the more promising Li 3 Y(PS 4 ) 2  material can be further enhanced multi-fold via aliovalent doping. Finally, we will discuss the relative merits of this new superionic conductor compared to current state-of-the-art superionic conductors. 
     Initial Candidate Selection 
     The initial pool of candidate lithium superionic conductors was constructed from the following: 
     1. All known ordered Li—P—S and Li-M-P—S structures from the 2015 version of the In-organic Crystal Structure Database (ICSD). Only non-redox-active elements were allowed for M. 
     2. Substitution of Ag with Li on all known ordered Ag—P—S and Ag-M-P—S structures from the ICSD. 
     Unique structures were identified from the pooled candidates using an in-house structure matching algorithm implemented in the Python Materials Genomics (pymatgen) materials analysis library. 
     DFT Calculations 
     All DFT calculations were performed using the Vienna Ab initio Simulation Package (VASP) within the projector augmented-wave approach. The exchange-correlation functional and calculation parameters were carefully selected to achieve a balance between computational accuracy and cost for the different types of calculations. 
     Spin-polarized calculations using the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA) functional was used for all structural relaxations. The convergence parameters, e.g. k-point density of at least 1000/(number of atoms in the unit cell) and energy cutoff of 520 eV, were similar to those used in the Materials Project (MP), which have been tested extensively over a broad range of chemistries. 
     All structures were fully relaxed using parameters similar to those used in the Materials Project (MP), which has been extensively tested over a broad range of chemistries and mate-rials. 
     All calculations were spin-polarized and performed using the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA) functional. A k-point density of at least 1000/(number of atoms in the unit cell) and an energy cutoff of 520 eV was used. Where available, pre-relaxed structures were first obtained from the MP using the Materials Application Programming Interface (API) to reduce computational cost. 
     The phase stability of a compound was estimated by determining its energy above the convex hull E hull  in the relevant Li—P—S and Li-M-P—S phase diagrams. Stable compounds have an E hull  of 0, and the higher the value, the more unstable the compound is at 0 K. Apart from the compounds of primary interest in this work, the energies of existing compounds were extracted from the MP database using the Materials Application Programming Interface (API). To account for overbinding of sulfur in DFT calculations, an energy correction of −0.66 eV per S atom for sulfides was applied. 
     The phase stability of all compounds of interest were estimated by constructing the relevant Li—P—S and Li-M-P—S phase diagrams using the convex hull construction. The energy above hull E hull  is then used as an estimate of thermodynamic stability. Stable compounds have an E hull  of 0, and the higher the value, the more unstable the compound is at 0 K. To account for overbinding of sulfur in PBE, an energy correction for sulfides was applied. 
     The electrochemical stability was assessed using the lithium grand potential phase diagram approach. In this approximation, Li is treated as the main mobile species and the solid electrolyte/electrode interface can be modeled as an open system with respect to Li. The relevant thermodynamic potential is therefore the grand potential, which can be approximated as φ≈E−μ Li N Li  in which E, N Li  and μ Li  are DFT total energy, number of lithium atoms in the open system, and lithium chemical potential, respectively. The phase equilibria at the anode and charged cathode can be approximated as the lithium superionic conductor composition at high μ Li =μ ∘   Li  and low μ Li =μ ∘   Li −5 eV (μ ∘   Li  is the chemical potential of metallic Li), respectively. 
     Automated non-spin-polarized Ab initio molecular dynamics (AIMD) simulations were performed in an NV T ensemble at elevated temperatures with a Nose-Hoover thermostat. A smaller plane-wave energy cutoff of 280 eV, a minimal Γ-centered 1×1×1 k-point mesh, and a time step of 2 fs were adopted. The simulation supercell sizes were at least 9 Å along each lattice direction. In line with previous studies, the simulation cell parameters were fixed at the fully relaxed cell parameters at 0 K. The Li +  diffusivity was obtained via a linear fit of the mean square displacement (MSD) with time, and Arrhenius plots were constructed from simulations at multiple temperatures to obtain the activation energy E a  and extrapolated room-temperature self-diffusivity D 300K  and conductivity σ 300K . 
     Climbing image nudged elastic band (CI-NEB) calculations were performed to determine the vacancy migration barriers for the most promising candidates. Overall charge neutrality was achieved via adding a positive background charge. The forces were converged to within 0.05 eV/Å. 
     Regarding the electronic structure band gap calculations were performed using the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional, due to the well-known underestimation of band gaps by semi-local functionals. 
     Parameterization of Screening Criteria 
       FIG. 1  is a flowchart showing a schematic of one example of the high-throughput (HT) screening framework described herein, which is tiered based on considerations of the relative importance of each property for lithium solid electrolyte applications and the computational cost required to calculate it. 
     First and foremost, all technologically relevant materials must be synthesizable, i.e., exhibit good phase stability. In this work, we have adopted a cutoff of E hull &lt;30 meV/atom, which is based on similar cutoffs adopted in previous HT computational materials screening efforts as well as the fact that Li 7 P 3 S 11  and Li 10 GeP 2 S 12 , both well-known superionic conductors, have been predicted to have an E hull  of 21-25 meV/atom. Of course, depending on circumstances, different cutoffs may be employed as well. 
     Second, a lithium superionic conductor must have a high Li +  conductivity at room temperature (σ 300K ). Due to the near unity transference number of lithium superionic conductors, σ 300K  exceeding 0.1 mS/cm should suffice for comparable performance with organic liquid electrolytes, though σ 300K &gt;1 mS/cm is preferred. However, obtaining converged diffusivity and conductivity numbers from AEID simulations is a highly computationally demanding process, usually requiring at least hundreds of picoseconds of simulation time (˜50,000-100,000 time steps) at multiple temperatures. Because we are interested only in superionic conductors with extremely high diffusivity, we have adopted a three-step diffusivity screening that includes a topological screening step, a quick estimation step, and a converged screening step. 
     The topological screening, which is the first step in the screening, is based purely on topological considerations. Only materials exhibiting &gt;1D diffusion networks with a minimum bottleneck size r C  of 1.75 Å, are considered as suitable candidates for lithium superionic conductors. This cutoff is slightly smaller than the channel size for the Li 10 GeP 2 S 12  superionic conductor (1.84 Å). A looser cutoff is used to avoid screening out too many candidates in the first screening step. The topological evaluation was carried out using the open source software Zeo ++ . 
     Regarding the quick estimation step, quick estimates of the diffusivity and activation energy were obtained using the mean square displacements obtained from short AIMD simulations of 50 ps at 800 K (MSD 800K ) and 1200 K (MSD 1200K ).  FIG. 2  shows a log-log plot of the MSD 1200K  versus MSD 800K  for a wide range of known compounds as well as the materials screened in this work. Square markers indicate known materials, which include well-established superionic conductors as well as relatively poorer conductors. Circle markers indicate new candidates screened. 
     The known superionic conductors evaluated include an approximate ordered model (Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5 , see Supplementary Information for details) for the recently reported Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3  superionic conductor, which has the LGPS structure and an extraordinarily high conductivity of 25 mS/cm. 
     The Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3  superionic conductor reported recently has an extraordinarily high ionic conductivity of 25 mS/cm, and has the same framework as the Li 10 GeP 2 S 12  (LGPS) that reported earlier. To estimate its diffusion characteristics for comparison with our proposed candidates, we first constructed a model based on an approximate composition of Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5 . Starting from the conventional cell of LGPS with formula Li 20 Ge 2 P 4 S 24 , all Ge were replaced with Si, one P atom was replaced with Si, and one S atom was replaced with Cl, yielding a cell formula of Li 20 Si 3 P 3 S 23 Cl, which reduces to Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5 . An enumeration was performed using the algorithm of Hart et al. was performed to yield all symmetrically distinct orderings of Si/P and S/Cl, and all structures were fully relaxed using DFT calculations employing the same parameters as outlined above. The ordering with the lowest energy structure was then used for subsequent investigations, e.g., AIMD, stability analyses, etc. It should be noted that the experimental structure is a disordered one, but we do not expect the diffusion characteristics to be significantly affected by the choice of the starting structure. 
     We may observe that all known superionic conductors fall within the white region bounded by MSD 800K &gt;5 Å 2  and MSD 1200K /MSD 800K &lt;7, and have therefore used these criteria in our screening process. The former criterion ensures a minimum baseline diffusivity, while the second criterion ensures that the activation energy is below ˜400 meV. 
     Non-spin-polarized ab initio molecular dynamics (AIMD) simulations were performed in an NV T ensemble at elevated temperatures with a Nose-Hoover thermostat. A smaller plane-wave energy cutoff of 280 eV, a minimal Γ-centered 1×1×1 k-point mesh, and a time step of 2 fs were adopted. The volume (V) was fixed at the relaxed 0 K volume for AEID simulations at elevated temperatures, in line with the usual approximations used in previous works. The simulation supercell sizes were at least 9 Å along each lattice direction. All calculations were automated by an in-house automated AEID workflow. 
     From the AEID simulations, the Li +  self-diffusivity can be obtained via the following expression: 
     
       
         
           
             
               
                 
                   
                     D 
                     = 
                     
                       
                         
                           
                             〈 
                             
                               [ 
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 
                                   r 
                                   ⁡ 
                                   
                                     ( 
                                     t 
                                     ) 
                                   
                                 
                               
                               ] 
                             
                             〉 
                           
                           2 
                         
                         
                           2 
                           ⁢ 
                           d 
                           ⁢ 
                           t 
                         
                       
                       = 
                       
                         
                           M 
                           ⁢ 
                           S 
                           ⁢ 
                           D 
                         
                         
                           2 
                           ⁢ 
                           d 
                           ⁢ 
                           t 
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
     where d is the dimensionality factor that equals 3 for 3D crystal structure, and {[Δr(t)]} is the average Li +  mean square displacement (MSD) over a time duration t. The self-diffusivity was obtained via a linear fit of the MSD vs 2 dt. The Arrhenius plot was constructed from diffusivities at multiple temperatures to obtain the activation energy (E a ) and the extrapolated room-temperature self-diffusivity (D 300K ). 
     The room-temperature Li +  conductivity was then estimated via Nemst-Einstein relation: 
       σ 300K =(ρ z   2   F   2   /RT )× D   300K ,   (2)
 
     where ρ, R and F are the molar density of Li +  in the unit cell, gas constant and Faraday&#39;s constant, respectively, and T=300 K and z=+1 were used in the expression. 
     Short AEID simulations of 60 ps were performed for the quick screening step in order to derive the mean square displacement cutoffs. The first 10 ps (˜5,000 time steps) were used for heating up as well as for equilibration, and the trajectories from 10 ps to 60 ps were used to estimate the MSD. Based on our previous AEID calculations, the diffusivities for most superionic conductors are found on par or beyond the magnitude of 10 −6  cm 2 /s at 800 K. By combining Eqn. (1) and benchmarking results shown in  FIG. 2 , we set MSD 800K &gt;5 Å 2  as the baseline diffusivity criterion. 
     Assuming that the diffusivity follows an Arrhenius relationship, we can also write the diffusivity as: 
     
       
         
           
             
               
                 
                   
                     D 
                     = 
                     
                       
                         D 
                         0 
                       
                       ⁢ 
                       
                         e 
                         
                           
                             - 
                             
                               E 
                               a 
                             
                           
                           
                             k 
                             ⁢ 
                             T 
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
     where E a  is the activation barrier and k is Boltzmann&#39;s constant. 
     Combining Eqn. (1) and (3), we can write: 
     
       
         
           
             
               
                 
                   
                     M 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     S 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     D 
                   
                   = 
                   
                     2 
                     ⁢ 
                     
                       D 
                       0 
                     
                     ⁢ 
                     d 
                     ⁢ 
                     t 
                     ⁢ 
                     
                       D 
                       0 
                     
                     ⁢ 
                     
                       e 
                       
                         
                           - 
                           
                             E 
                             a 
                           
                         
                         
                           R 
                           ⁢ 
                           T 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
           
         
       
     
     Let us consider the ratio of MSD at 1200 K and 800 K for the same simulation time period t. 
     
       
         
           
             
               
                 
                   
                     
                       M 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       S 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         D 
                         
                           1200 
                           ⁢ 
                           K 
                         
                       
                     
                     
                       M 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       S 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         D 
                         
                           800 
                           ⁢ 
                           K 
                         
                       
                     
                   
                   = 
                   
                     e 
                     
                       
                         
                           - 
                           
                             
                               E 
                               a 
                             
                             
                               1200 
                               ⁢ 
                               k 
                             
                           
                         
                         + 
                         
                           
                             E 
                             a 
                           
                           
                             800 
                             ⁢ 
                             k 
                           
                         
                       
                       , 
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
             
               
                 
                   
                     E 
                     a 
                   
                   = 
                   
                     
                       2400 
                       · 
                       k 
                     
                     ⁢ 
                     
                       
                           
                       
                       ⁢ 
                       
                           
                       
                     
                     ⁢ 
                     
                       
                         ln 
                         ⁡ 
                         
                           ( 
                           
                             
                               M 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               S 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 D 
                                 
                                   1200 
                                   ⁢ 
                                   K 
                                 
                               
                             
                             
                               M 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               S 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               
                                 D 
                                 
                                   800 
                                   ⁢ 
                                   K 
                                 
                               
                             
                           
                           ) 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       For 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         
                           M 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           S 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           
                             D 
                             
                               1200 
                               ⁢ 
                               K 
                             
                           
                         
                         
                           M 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           S 
                           ⁢ 
                           
                               
                           
                           ⁢ 
                           
                             D 
                             
                               800 
                               ⁢ 
                               K 
                             
                           
                         
                       
                     
                     &lt; 
                     7 
                   
                   , 
                   
                     
                       E 
                       a 
                     
                     &lt; 
                     
                       402 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         meV 
                         . 
                       
                     
                   
                 
               
               
                 
                     
                 
               
             
           
         
       
     
     The relative trends observed in  FIG. 2  are consistent with the known properties of the conductors, with faster conductors such as Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5  exhibiting larger MSD 800K  (higher diffusivity at 800 K) and lower MSD 1200K /MSD 800K  (lower activation barriers) than slightly poorer conductors such as Li 10 GeP 2 S 12 . It should be noted that though Li 7 P 3 S 11  was experimentally reported to have an ionic conductivity of 17 mS/cm, recent computational work by the authors of this work suggests that the intrinsic conductivity of this material may be as high as 60 mS/cm, which accounts for its high MSD 800K . Also, though Li 3 OCl x Br 1-x  and Na 3 PS 4  were initially reported to be superionic conductors, the pristine crystalline phases have been shown to have low room-temperature conductivities, and higher ionic conductivities are obtained only with the introduction of a large number of defects/dopants (e.g., doped-Na 3 PS 4 ) and/or via amorphization. 
     Finally, in the converged screening step, longer AEID simulations at six temperatures were performed on the materials that pass the first two screening steps to obtain converged diffusivities (and conductivities) and activation barriers. 
     The above three-step screening process allows us to rapidly eliminate poor candidates with a minimum amount of computational resources, and devote expensive AIMD simulations to obtain converged diffusivity statistics on the most promising materials. 
     Besides excellent Li +  conductivity, a solid electrolyte for all-solid-state rechargeable 
     lithium-ion batteries must also be electronically insulating and exhibit good electrochemical stability against the electrodes. An assessment of these properties were carried out for the most promising candidates. 
     Identification of Potential Candidates 
     Table 1 summarizes the phase stability, topological parameters and rapid AIMD screening results of all new Li—P—S and Li-M-P—S candidates. The rapid AIMD screening results are also presented in  FIG. 2  for comparison with known superionic conductors. We note that though we did perform Ag for Li substitution of Ag—P—S compounds, the derived compounds are all well-known compounds such as Li 7 P 3 S 11  and Li 3 PS 4  that have already been extensively explored as lithium superionic conductors, and thus will not be considered further here. 
     
       
         
           
               
               
               
               
               
               
             
               
                 TABLE 1 
               
               
                   
               
               
                   
                 Source  
                 E hull   
                 r c   
                 MSD 800K   
                 MSD 800K   
               
               
                 Compound 
                 (ICSD number) 
                 (meV/atom) 
                 (Å) 
                 (Å 2 ) 
                 MSD 800K   
               
               
                   
               
             
            
               
                   
               
            
           
           
               
            
               
                 Promising candidates 
               
            
           
           
               
               
               
               
               
               
            
               
                 Li 3 Y(PS 4 ) 2  (C2/c) 
                 Ag 3 Y(PS 4 ) 2  (417658) 
                 2 
                 1.88 
                 65.1 
                 4.5 
               
               
                 Li 3 PS 4 Cl 2  (C2mm) 
                 Ag 3 PS 4 Cl 2   (416587) 
                 17 
                 1.76 
                 77.9 
                 3.1 
               
            
           
           
               
            
               
                 Candidates failing at least one screening criteria 
               
            
           
           
               
               
               
               
               
               
            
               
                 Li 15 P 4 S 16 Cl 3  (I 4 3d) 
                 Ag 18 P 4 S 16 Cl 3  (416586) 
                 8 
                 1.76 
                 2.2 
                 145.1 
               
               
                 LiZnPS 4  (I 4 ) 
                 ICSD (95785) 
                 0 
                 1.83 
                 24.7 
                 13.9 
               
               
                 LiAl(PS 3 ) 2  (C2/c) 
                 ICSD (425979) 
                 0 
                 1.82 
                 27.8 
                 16.5 
               
               
                 Li 2 Zn(PS 3 ) 2  (C2/c) 
                 Ag 2 Zn(PS 3 ) 2  (72719) 
                 33 
                 1.84 
                 — 
                 — 
               
               
                 LiIn(PS 3 ) 2  (P 3 1c) 
                 AgIn(PS 3 ) 2  (202185) 
                 0 
                 0 
                 — 
                 — 
               
               
                 LiZnPS 4  (Pna2 1 ) 
                 AgZnPS 4  (48197) 
                 16 
                 0 
                 — 
                 — 
               
               
                   
               
            
           
         
       
     
     Among the new quaternary compounds, only Li 3 Y(PS 4 ) 2  (LYPS) and Li 5 PS 4 Cl 2  (LP-SCl) satisfy all the initial screening criteria: low E hull , r C &gt;1.75 Å, MSD 800K &gt;5 Å and MSD 1200K /MSD 800K &lt;7. Their MSD 800K  are on par with that of the leading LGPS-based candidate, Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5 , but their MSD 1200K /MSD 800K  ratios are slightly higher. Though Li 15 P 4 S 16 Cl 3 , LiZnPS 4  and LiAl(PS 3 ) 2  are also predicted to have fairly low E hull  and reasonably high MSD 800K , their MSD 1200K /MSD 800K  are far too high, indicating high activation barriers. The remaining candidates do not pass either the phase stability criterion or the topological screening. During the preparation of this manuscript, it has come to our attention that the LiZnPS 4  candidate in Table 1 has been investigated as a superionic conductor. Our screening calculations show that the stoichiometric LiZnPS 4  compound fails the MSD ratio cutoff by a factor of 2, which is consistent with the high activation barriers reported for the stoichiometric compound reported in Richards et al.&#39;s work. A more in-depth comparison of our proposed candidates with known superionic conductors is provided below. 
     The crystal structures of LYPS and LPSCl are shown in  FIG. 3 , and the relaxed lattice parameters are given in Table 2. LYPS and LPSCl belong to the monoclinic C2/c and orthorhombic C2mm space groups, respectively. Both LYPS and LPSCl structures have two symmetrically distinct Li sites, which are labeled as Li1 and Li2 in  FIGS. 6 and 7 . 
     
       
         
           
               
               
               
               
               
               
               
               
             
               
                 TABLE 2 
               
               
                   
               
               
                   
                   
                 a 
                 b 
                 c 
                 α 
                 β 
                 γ 
               
               
                 Compound 
                 atoms/cell 
                 (Å) 
                 (Å) 
                 (Å) 
                 (°) 
                 (°) 
                 (°) 
               
               
                   
               
             
            
               
                   
               
            
           
           
               
               
               
               
               
               
               
               
            
               
                 Li 3 Y(PS 4 ) 2   
                 56 
                 17.122 
                 9.290 
                 9.137 
                 90.0 
                 122.3 
                 90.0 
               
               
                 Li 5 PS 4 Cl 2   
                 24 
                 7.212 
                 10.494 
                 6.024 
                 90.0 
                 90.0 
                 90.0 
               
               
                   
               
            
           
         
       
     
     Li +  Conductivities and Mechanisms 
     Long AEID simulations of at least 200 ps at multiple temperatures were performed on the promising LYPS and LPSCl candidates.  FIG. 4  shows the Arrhenius plot of the diffusivity versus 1000/temperature for the two candidates, and Table 3 summarizes the key conductivity properties extracted. In  FIG. 4 , dashed lines indicate the further optimized candidates obtained via aliovalent doping or substitution of Li 3 Y(PS 4 ) 2  (Li 3 La(PS 4 ) 2 , Li 3.125 Y 0.875 Ca 0.125 (PS 4 ) 2  and Li 2.875 Y 0.875 Zr 0.125 (PS 4 ) 2 ). 
     
       
         
           
               
               
               
               
               
             
               
                 TABLE 3 
               
               
                   
               
               
                   
                 σ 300K   
                 Error range of σ 300K   
                 E α   
                 D 300K   
               
               
                 Formula 
                 (mS/cm) 
                 (mS/cm) 
                 (meV) 
                 (cm 2 /s) 
               
               
                   
               
             
            
               
                 Li 3 Y(PS 4 ) 2   
                 2.16 
                 [1.46, 3.19] 
                 278 
                 3.56 × 10 −8   
               
               
                 Li 5 PS 4 Cl 2   
                 1.85 
                 [1.38, 2.47] 
                 304 
                 1.36 × 10 −8   
               
               
                   
               
            
           
         
       
     
     The estimated activation energies E a  for LYPS and LPSCl are 278 meV and 304 meV, respectively. The extrapolated room-temperature conductivities are 2.16 mS/cm for LYPS and 1.85 mS/cm for LPSCl, i.e., both candidates are indeed verified to be lithium superionic conductors. 
     To further understand the atomistic diffusion mechanisms and pathways, the Li +  probability density function (PDF) was calculated from the AIMD simulations of the two candidates at 800 K, and CI-NEB calculations were performed to calculate the vacancy migration barriers in the identified pathways.  FIGS. 5A and 5B  shows isosurfaces of Li +  probability density distribution P (light blue) for Li 3 Y(PS 4 ) 2  and Li 5 PS 4 Cl 2  from AIMD simulations at 800 K with P=0.0001 a 0   −3  (a 0  is the Bohr radius). From  FIGS. 5A and 5B , we observe that both candidates have 3D diffusion networks. Schematics of the identified diffusion pathways and the corresponding calculated CI-NEB migration barriers are given in  FIGS. 6 and 7  for LYPS and LPSCl, respectively.  FIG. 6( a )  shows the investigated Li vacancy diffusion paths in Li 3 Y(PS 4 ) 2  viewed along b direction.  FIGS. 6( b )-6( d )  show the calculated CI-NEB migration barriers for selected percolating paths.  FIG. 7  shows the investigated Li vacancy diffusion paths in Li 5 PS 4 Cl 2  viewed along the a direction ( FIG. 7 a   ), the b direction ( FIG. 7 b   ).  FIGS. 7( c )-7( e )  show calculated CI-NEB migration barriers for selected paths. 
     For LYPS, there are five symmetrically distinct hops between neighboring Li sites, namely, A→B, B→F, B→C, C→E and C→H (see  FIG. 6 ). We may observe that the lowest barrier percolating pathway in LYPS is A→B→C→D (all Li1) with a barrier of 210 meV. To enable &gt;1D diffusion, the next lowest barrier pathway is G→C→H (all Li2) along the c direction with an overall barrier of 268 meV. The E→B→F pathway has a much higher overall barrier of 406 meV, which suggests it is more likely for a vacancy at E to diffuse to B via the lower barrier E→G→C→B path. The estimated effective 3D vacancy migration barrier (268 meV) is therefore in good agreement with the activation energy from AIMD simulations (278 meV). 
     For LPSCl, we may observe that the crystal structure comprises layers of Li1 and Li2 (see  FIG. 7 ) in the b-c plane stacked along the a direction. Within each Li1 plane, the 
     vacancy migration barriers are relatively low (A→B→C→D with overall barrier of 166 meV, 
     and B→E with barrier of 217 meV). For 3D diffusion, the H→G→F Li1-only path connecting 
     different Li1 layers has the lowest overall barrier of 321 meV. This is again in reasonably good agreement with the AIMD activation energy of 304 meV. All other paths involving vacancy hops between Li1 and Li2 sites have significantly higher barriers (&gt;380 meV). 
     Electronic Band Gap 
       FIG. 8  shows the calculated HSE densities of states for LYPS and LPSCl.  FIG. 8( a )  shows the density of states for Li 3 Y(PS 4 ) 2  (band gap E g =3.41 eV) and  FIG. 8( b )  shows the density of states for Li 5 PS 4 Cl 2  (E g =3.57 eV) using the HSE screened hybrid functional. Both candidates are large band gap (&gt;3 eV) insulators. In both candidates, the valence band maximum is dominated by anion S-p character, whereas the conduction band minimum has contributions from the framework cations (Y and P) as well as anion S. 
     The band gap is also an upper limit for the intrinsic stability of the material against 
     reduction (acceptance of an electron) and oxidation (loss of an electron). Similar to other sulfide-based solid-electrolytes, the intrinsic electrochemical stability of the two candidates are limited to ˜3.5 eV. 
     Electrochemical Stability 
     Better estimates of the electrochemical stabilities of LYPS and LPSCl were obtained using the lithium grand potential approach. Table 4 summarizes the predicted phase equilibria at the solid electrolyte/anode (metallic Li) interface and solid electrolyte/charged 5V cathode interface. The dominant product at the anode is Li 2 S in all cases, which is a good electronic insulator and reasonable Li conductor, especially as an amorphous interphase. The other products at the anode are YP, a semiconductor with band gap of ˜1 eV, 55  and Li 3 P. On the solid electrolyte/charged cathode, P 2 S 5  is always predicted to be one of the products. However, the presence of S 2 Cl 2  and PCl 5  at the LPSCl/cathode interface may prove problematic in real-world applications as they undergo hydrolysis readily to form HCl. 
     
       
         
           
               
               
               
             
               
                 TABLE 4 
               
               
                   
               
               
                   
                 Phase equilibria at 5 V cathode 
                 Phase equilibria at anode 
               
               
                 Electrolyte 
                 μLi = μºLi eV 
                 μLi = (μºLi − 5) eV 
               
               
                   
               
             
            
               
                   
               
            
           
           
               
            
               
                 Promising candidates 
               
            
           
           
               
               
               
            
               
                 Li 3 Y(PS 4 ) 2   
                 YPS 4  + 0.5 P 2 S 5  + 1.5S 
                 YP + Li 3 P + 8 Li 2 S 
               
               
                 Li 5 PS 4 Cl 2    
                 0.067 PCl 5  + 0.833S 2 Cl 2  + 0.467 P 2 S 5   
                 Li 3 P + 4 Li 2 S + 2 LiCl 
               
            
           
           
               
            
               
                 Known superionic conductors 
               
            
           
           
               
               
               
            
               
                 Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5   
                 0.75 P 2 S 7  + 0.25 S 2 Cl 2  + 
                 1.5 Li 3 P + 0.3 Li 21 Si 5  + 0.5 
               
               
                 Li 10 GeP 2 S 12   
                 1.5 SiS 2  + 2.75 S 
                 LiCl + 11.5 Li 2 S 2 Li 3 P + 0.25 
               
               
                 Li 7 P 3 S 11   
                 P 2 S 5  + GeS 2  + 5 S  
                 Li 15 Ge 4  + 12 Li 2 S 3  Li 3 P + 11 
               
               
                   
                 1.5 P 2 S 5  + 3.5 S 
                 Li 2 S 
               
               
                   
               
            
           
         
       
     
     For comparison, Table 4 also presents the predicted phase equilibria for the Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5  model of the recently reported Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3  superionic conductor as well as Li 10 GeP 2 S 12  and Li 7 P 3 S 11 . Similar to the candidates identified in this work, Li 2 S is predicted to be the dominant product at the anode/electrolyte interface in all instances, with the small band gap Li 3 P comprising a relatively small fraction. For Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5  and Li 10 GeP 2 S 12 , there is an additional Li 21 Si 5  or Li 15 Ge 4  phase, which also have a small band gap, consistent with previous experimental studies. At the cathode/Li 10 S 1.5 P 1.5 S 11.5 Cl 0.5  interface, S 2 Cl 2  is predicted to be one of the products, though the proportion is much less in comparison to LPSCl due to the much lower content of Cl. 
     Analysis 
     From the results in the preceding sections, Li 3 Y(PS 4 ) 2  (LYPS) and Li 5 PS 4 Cl 2  (LPSCl) have emerged as promising new lithium superionic conductors based on a comprehensive screening of the Li—P—S and Li-M-P—S chemical spaces. Both candidates exhibit good phase stability (low E hull ) and excellent topological characteristics (&gt;1 D large conduction channels), and are predicted to be electronic insulators with high Li +  conductivities (exceeding 1 mS/cm). The Li +  conduction mechanisms and migration barriers were elucidated using CI-NEB calculations, and the results further confirm the predictions from the AIMD simulations. In addition, both candidates comprise entirely of earth-abundant elements, making them practical from a cost perspective. 
     Both candidates are derived from the replacement of Ag with Li in known quaternary Ag thiophosphates in the ICSD. We see this as a further positive attribute of the two candidates as ion exchange from the known Ag-based compounds is therefore a potential initial synthesis route that can be explored. For example, ion exchange has similarly been used to synthesize the well-known Li 7 P 3 S 11  superionic conductor from Ag 7 P 3 S 11 . We speculate that due to the significantly larger ionic radii of Ag (129 pm) compared to Li (90 pm), Li-substituted Ag compounds may present large percolating voids conducive to fast 3D Li mobility. Such a strategy can certainly be expanded beyond just the thiophosphate chemistries that are the focus of this work. However, we would point out that the large ionic radii difference between Ag +  and Li +  can potentially lead to incompatibility of Li with the Ag-based host framework, which is why a computational assessment of phase stability is a critical first step to determine the likelihood of synthesis. Also, not all Ag compounds have percolating 3D diffusion networks of sufficient channel size. Here again, the efficient tiered screening approach outlined in this work based on inexpensive topological analysis followed by more computationally intensive first principles calculations can provide useful guidelines. 
     Between the two candidates, LYPS is believed to be the more promising one. Not only is LYPS predicted to have a marginally higher Li +  conductivity than LPSCl in AIMD 
     simulations, it is also predicted to be significantly more stable (E hull =2 meV/atom) and its lack of Cl means that there is likely to be fewer issues with reaction products at higher voltages. We will note that like all sulfide-based materials, air and moisture stability may be a potential area of concern, though this limitation has not prevented the development of prototype all-solid-state rechargeable lithium-ion batteries utilizing other sulfide solid electrolytes. Like other sulfides, both materials are predicted to be relatively soft, which should make it easier to achieve low porosity using cold-pressing methods. 
     Further Optimization of LYPS 
     To explore if further enhancement of the conductivity of LYPS is possible, we performed isolvalent substitutions and aliovalent doping of LYPS. La 3+  was examined as a potential substitute for Y 3+  due to its slightly larger ionic radii (117 pm compared to 104 pm for Y 3+ ). The computed E hull  for Li 3 La(PS 4 ) 2  is 20 meV, significantly higher than LYPS, and its ionic conductivity is only slightly higher at 3.27 mS/cm with a slightly lower activation energy of 263 meV (see  FIG. 4 ), which can be attributed to the ˜5% volume expansion and increase in channel size (1.92 Å). 
     Unlike isovalent substitutions, aliovalent doping can have the additional effect of introducing Li +  vacancies or interstitials. Both Ca 2+  or Zr 4+  dopants that have comparable ionic radii to Y 3+  were explored using a 1×1×2 supercell of LYPS, with the introduction of Li +  interstitials and vacancies, respectively. Table 5 summarizes the dopant formation energies and room-temperature Li +  conductivities for the doped structures. Both Ca 2+  and Zr 4+  were found to have reasonably low dopant formation energies of 0.63 eV and 0.26 eV, respectively. From AIMD simulations, we find that aliovalent doping of LYPS with the introduction of either vacancies or interstitials can lead to multi-fold increases in its ionic conductivity. Substitution of 12.5% of Y 3+  with Ca 2+  and Zr 4+  leads to extrapolated room temperature conductivities of 7.14 mS/cm and 5.25 mS/cm, respectively, with corresponding decreases in activation energies to 231 meV and 241 meV, respectively (see  FIG. 4 ). 
     
       
         
           
               
               
               
               
               
               
               
             
               
                 TABLE 5 
               
               
                   
               
               
                   
                   
                   
                   
                   
                 error range of 
                   
               
               
                   
                   
                 E f   
                 E hull   
                 σ 300K   
                 
                   
                 
                 E a   
               
               
                 Dopant 
                 Formula 
                 (eV) 
                 (meV/atom) 
                 (mS/cm) 
                 (mS/cm) 
                 (meV) 
               
               
                   
               
             
            
               
                 Ca 
                 Li 3.125 Y 0.875 Ca 0.125 (PS 4 ) 2   
                 0.63 
                 6 
                 7.14 
                 [4.67, 10.92] 
                 231 
               
               
                 Zr 
                 Li 2.875 Y 0.875 Zr 0.125 (PS 4 ) 2   
                 0.26 
                 4 
                 5.25 
                 [3.77, 7.31]  
                 241 
               
               
                   
               
               
                     indicates data missing or illegible when filed 
               
            
           
         
       
     
     Due to computational cost considerations, our explorations of dopant optimization is limited by the size of the supercell accessible within AIMD simulations. Nevertheless, the doping results are a proof of concept that there is significant scope for further fine-tuning of dopant and Li concentration in LYPS to achieve even higher conductivities, a claim that we hope will be verified by experimental efforts at synthesizing undoped and doped LYPS. 
     Comparison with Other State-of-the-Art Superionic Conductors 
     In comparison with state-of-the-art sulfide superionic conductors such as Li 7 P 3 S 11 , the LGPS family (Li 10 GeP 2 S 12  and Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3 , LYPS (undoped or doped) has slightly lower Li +  conductivity. However, bulk ionic conductivity is no longer the critical factor in all-solid-state battery performance beyond 1 mS/cm. Indeed, other properties such as interfacial stability play a far more critical role. For instance, though the Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3  superionic conductor recently reported has an extraordinarily high room temperature ionic conductivity of 25 mS/cm, its interfacial stability is much poorer than the Li 9.6 P 3 S 12  composition in the same structure, which has a lower conductivity of ˜1 mS/cm. The result is that Li 4 Ti 5 O 12 , which has a voltage of 1.5 V against Li/Li + , had to be used as the anode with Li 9.54 Si 1.74 P 1.44 S 11.7 Cl 0.3 , lowering achievable energy densities due to the low overall operating voltage of ˜2.5 V. In contrast, full cell performance at a relatively high operating voltage of up to 4.2 V was demonstrated for Li 9.6 P 3 S 12  with standard graphitic anodes. 
     LYPS compares favorably to these known superionic conductors in terms of both phase and electrochemical stability. The calculated E hull  of LYPS is only 2 meV/atom, sub-stantially lower than that of Li 7 P 3 S 11  (21 meV/atom), Li 10 GeP 2 S 12  (25 meV/atom) and Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5  (30 meV/atom). Recently, computational evidence has been reported of extraordinarily high Li conductivities exceeding 50 mS/cm in the Li 1+2x Zn 1−x PS 4  solid solution, a compound that was also considered in our screening. However, these high conductivities were obtained only with the introduction of a large number of defects, requiring high predicted synthesis temperatures exceeding 950 K. In comparison, doped LYPS with conductivities of up to 7 mS/cm still maintains a relatively low E hull  and small dopant formation energies. 
     In terms of interfacial stability, there are no reaction products of major concern at the cathode/LYPS interface, unlike Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5  where the presence of Cl is predicted to result in the formation of S 2 Cl 2 . On the anode/LYPS interface, the presence of the semiconducting YP phase may be of potential concern, though its band gap (˜1 eV) is still higher than the Li—Si alloys (e.g., band gaps of 0.6 eV for Li 12 Si 7  and 0.08 eV for Li 7 Si 3 ) predicted to form at the anode/Li 10 Si 1.5 P 1.5 S 11.5 Cl 0.5  interface. An electrically insulating interface is desired for passivation to avoid further propagation of the reaction front. Furthermore, Li—Si alloys are also well known to undergo significant volume expansion (in excess of 300%) at high lithiation, which may be detrimental to maintaining intimate electrode/electrolyte contact. 
     In summary, LYPS may present an overall better balance of properties as a lithium superionic conductor solid electrolyte for all-solid-state battery applications. It has clearly better predicted phase stability, and likely better interfacial stability based on the predicted phase equilibria at the electrode/electrolyte interface. Its conductivity, though somewhat lower than some of the state-of-the-art candidates, is sufficiently high that it is not likely to be a limiting factor, and can potentially be further improved with the demonstrated doping strategies. 
     Applications 
     Lithium superionic conductor electrolytes identified in accordance with the techniques described herein may be used in a variety of applications. For instance, they may be employed in a solid state battery.  FIG. 9  is a schematic diagram of an all-solid-state lithium battery during a charge process. The battery contains three components: a metallic lithium anode, a charged cathode, and a lithium superionic conductor electrolyte. During the charge process, Li +  ions are shuttled from cathode to anode through the solid electrolyte. The reverse procedure takes place during a discharge process. 
     Other Superionic Conductor Materials 
     Inspired by the similarity of the chemical space of Li—P—S and the recently identified superionic conductor Na—P—S as well as the prediction of Li 3 Y(PS 4 ) 2  as described above as a new lithium superionic conductor electrolyte (SICE), we have conducted a high-throughput (HT) screening in Na-M-P—S chemical spaces for potential sodium SICEs. To increase the coverage of these chemical spaces, additional compounds substituted from existing A-M-P—X (A=Li, Na, Ag, K; X═S, O; M is non-redox-active element) chemical spaces are also included. From our calculations, we have identified one highly promising new sodium superionic conductor Na 3 Y(PS 4 ) 2  for Na-ion battery applications. 
     Our system and method according to present principles can potentially address the safety issues caused by the use of organic liquid electrolytes in sodium ion batteries. The material is completely new in that no such structure is known in the literature. We predicted the new sodium superionic conductor Na 3 Y(PS 4 ) 2  from density functional theory (DFT) calculations, and it has been subsequently synthesized and confirmed by X-ray diffraction (XRD) measurements. Based on our prediction, it should exhibit good phase stability, high Na +  conductivity, low electronic conductivity and good electrochemical stability. 
     The flow chart describing one example of a procedure of our first-principles high-throughput (HT) screening for superionic conductors is given in  FIG. 10 , which is similar to the flowchart shown in Figure for 1 for the procedure for identifying lithium superionic conductors. 
     Our initial candidate pool comprised two groups: (i) existing Na-M-P—S(M is non-redox-active elements) compounds from Inorganic Crystal Structure Databse (ICSD) 2016 version and (ii) A/Na, X/S substitution from A-M-P—X (A=Li, Na, Ag, K; X═S, O) chemical space. We first filtered out highly unstable materials (measured by E hull ; the higher the value, the more unstable a compound is) with E hull &gt;30 meV/atom. This is followed by a three-step diffusivity screening process. This screening procedure involves a topological analysis excluding those have too small or only 1D Li diffusion channels, a new short AIMD estimation process from two aspects of consideration (baseline diffusivity and activation energy) and a converged long-term AIMID simulation at multiple temperatures. The detailed cutoffs of the criteria are modified from our previous published work for Li superionic conductors ( Chem. Mater.  29, 2474-2484 (2017)) due to the Li and Na ionic radius difference. 
       FIG. 11  is the working schematics of a Na-ion battery. Anode, cathode and solid electrolyte are three components of the all-solid-state battery. During charging process, Na +  ions are transported from cathode to anode through electrolytes. The reverse procedure happens at discharge process. 
     Because of the rareness of sodium superionic conductors and more requirements not only limited to conductivity, it is important to search for new sodium solid electrolytes with a combination of good phase stability, electronic conductivity and electrochemical stability. This motivates our first-principles guided new sodium superionic conductors investigation. There are three main stages during the development of our system and method according to present principles. 
     (i) Concept stage: We used our modified HT screening strategy to search for good sodium superionic conductors and predicted new Na 3 Y(PS 4 ) 2  materials with good phase stability and excellent Na +  conductivity. This computational work was followed by experimental synthesis attempts, which can verify the theoretical prediction and form a complete computational guided new material discovery chain. 
     (ii) Computer modeling simulation stage: The first-principles high throughput screening approach shown in  FIG. 10  is developed to search for the promising sodium conductors. Taking the great advantages of the robust computational method, we have investigated and developed several chemical spaces which are less studied before and  FIG. 12  gives the phase diagram that we newly contributed to estimate the phase stability of Na 3 Y(PS 4 ) 2  (the circular dots are all stable phases in the Na 2 S—P 2 S 5 —Y 2 S 3  phase diagram). 
     (iii) Experimental data stage: We have successfully reproduced the pure Ag 3 Y(PS 4 ) 2  phase (XRD patterns shown in  FIG. 13 ). New Na 3 Y(PS 4 ) 2  material has already been synthesized and XRD results show the signal peaks of this new compound ( FIG. 14 ). 
     Rechargeable Na-ion batteries have enjoyed significant attention and interest in recent years due to the abundant sodium resources and more possible new chemistries. As the key component of all-solid-state batteries, sodium solid electrolytes overcome the leakage and flammable problem caused by organic liquid electrolytes. For commercial application, a combination of multiple properties should be taken into consideration. The promising material we predicted (Na 3 Y(PS 4 ) 2 ) has sufficiently high conductivity (&gt;10 mS/cm) and good phase stability as well as electrochemical stability. Our computational guided new material discovery strategy can greatly speed up the procedure of searching for promising candidates and reduce the efforts in experiments. Experimental achievement also verifies the validation of our HT screening method.