Tackling structural complexity in Li 2 S-P 2 S 5 Solid-State Electrolytes using Machine Learning Potentials
The lithium thiophosphate (LPS) material class provides promising candidates for solid-state electrolytes (SSEs) in lithium ion batteries due to high lithium ion conductivities, non-critical elements, and low material cost. LPS materials are characterized by complex thiophosphate microchemistry and structural disorder influencing the material performance. To overcome the length and time scale restrictions of ab initio calculations to industrially applicable LPS materials, we develop a near-universal machine-learning interatomic potential for the LPS material class.The trained Gaussian Approximation Potential (GAP) can likewise describe crystal and glassy materials and different P-S connectivities P m S n . We apply the GAP surrogate model to probe lithium ion conductivity and the influence of thiophosphate subunits on the latter. The materials studied are crystals (modifications of Li 3 PS 4 and Li 7 P 3 S 11 ), and glasses of the x Li 2 S-(100 -x )P 2 S 5 type (x = 67, 70 and 75).The obtained material properties are well aligned with experimental findings and we underscore the role of anion dynamics on lithium ion conductivity in glassy LPS.The GAP surrogate approach allows for a variety of extensions and transferability to other SSEs.
I. INTRODUCTION<!>A. Lithium Ion Mobility<!>III. CONCLUSION
<p>While lithium-ion batteries with liquid electrolytes entered the market in 1991, all-solidstate lithium-ion batteries (ASS-LIBs), although investigated for decades, are still not widely in use. [1][2][3][4][5][6] ASS-LIBs promise several advantages in comparison to liquid electrolyte batteries: higher power density, minimized safety and fire hazards, longer cycle lifetimes, more comprehensive temperature ranges, and enhanced energy density by potential usage of Li metal anodes. 1,7,8 Solid electrolytes of the Li 2 S-P 2 S 5 material class have gained substantial attention due to their favorable properties. 6 First, they possess high conductivities of up to 10 −2 Scm −1 , which ranks them among the most conductive solid electrolytes such as Li 10 GeP 2 S 12 or Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 . 9,10 Secondly, they are composed of the earth-abundant elements sulfur and phosphorous enabling sustainable applications at large scales.</p><p>However, this material class's potential is hampered by the poor understanding of the relevant structure-property relations. This manifests itself in huge deviations in Li-ion conductivity between theory and experiment. As such, β-Li 3 PS 4 serves as an illustrating example.</p><p>Experimental studies report a lithium ion conductivity of approximately 10 −7 Scm −1 , making the material unsuitable for industrial battery applications. 11 In contrast, an ab initio study predicts a conductivity of 10 −1 Scm −1 ; a six orders of magnitude deviation from experiment that would make the material the new record holder in solid-state lithium ion conduction. 12 Such huge discrepancies often arise from computational limitations that constrain tractable system sizes and sampling times. In the LPS case, high-resolution TEM images for instance revealed the presence of crystalline nanoparticles in otherwise amorphous regions, highlighting that conductivity calculations of ideal crystals are too short-sighted for this materials class. 13 The problem is further accentuated by the complex chemical structure of LPS. 14,15 A large structural variability at the molecular level, more precisely different thiophosphate poly-anions, are found in all crystalline and amorphous materials. 6,16 For a detailed description of the lithium ion conductivity in LPS we thus need to tackle these challenges: structural complexity of LPS glass compounds, chemical reactivity of thiophosphates, and the influence of anion composition on the lithium ion conductivity.</p><p>Here we tackle these challenges by replacing the computationally demanding direct firstprinciples calculations with a surrogate machine-learning (ML) model. Once trained, this Gaussian Approximation Potential (GAP) model allows for an upscaling of both time and length scale: molecular dynamics (MD) simulations covering up to several nanoseconds and system sizes of several thousand atoms become feasible. Furthermore, the flexibility offered by the ML approach allows one to implement a GAP model that is more versatile and can better represent the crucial complex chemistry than a classical force field. 17 We present a data-efficient iterative training approach to extend an earlier ML force field to yield a near-universal description of the LPS material class. 18,19 In the first part of this work we present our data-efficient training protocol and evaluate the GAP model on a.) its predictive accuracy for lithium ion conductivity and b.) its ability to reproduce two known phase transitions in crystalline Li 3 PS 4 . The second part focuses on the influence of anion composition on the lithium ion conductivity of different LPS glass compounds.</p><!><p>We obtain the reactive GAP model used to describe the LPS class by fitting to densityfunctional theory (DFT) training data computed with the FHI-aims full-potential package (see SI for full details of the computational settings). 20 The underlying approach is based on three consecutive steps: defining the anion lattice, sampling of Li-sites, and fine-tuning the materials density. In the first place, only the dominant anion species (e.g. PS 3− 4 and P 2 S 4− 7 ) are taken into consideration and utilized in a ratio that represents the desired stoichiometry correctly. For a data-efficient sampling of lithium sites, we sample Li-ion distributions on stable and meta-stable Li sites in a quasi-Monte Carlo like fashion. The materials density is obtained by an iterative compression scheme. Convergence, a detailed step-by-step description of the underlying algorithm, and numerical error assessments of the training procedure are given in the SI. The benefit of this approach is that it allows the free tuning of stoichiometries and polyanion ratios. Note that we use a charge free approach in this work. We previously showed that neglecting long-range electrostatics is unproblematic for describing lithium ion transport in isotropic bulk-like systems. 18 As a first validation of our GAP model, we turn to the Li-ion conductivity of crystalline LPS materials (α, β, γ Li 3 PS 4 and Li 7 P 3 S 11 ) at finite temperature, predicted from MD sim-ulations via the Nernst -Einstein equation (see the SI for details). Using the GAP model we evaluate the ionic conductivity from the mean-square-displacement (MSD) sampled during 2 ns MD simulations for every crystalline compound at various temperatures between 400 and 800 K. Room temperature (RT) conductivities are extrapolated from a linear fit. Note that for crystalline Li 7 P 3 S 11 we required longer simulation times of up to 13 ns to reach converged conductivities, i.e. time scales that would essentially be prohibitive for direct ab initio MD. While Li-ion conductivity in LPS is usually dominated by diffusion of defects (Li + vacancies), Li 7 P 3 S 11 exhibits a more collective Li + motion yielding the observed high conductivity. 18,21,22 As seen in Figure 1, a broad range of Li-ion conductivities are exhibited in LPS. While high RT conductivities of up to 3.6 10 −3 Scm −1 are found for α-Li 3 PS 4 and Li 7 P 3 S 11 , β and γ-Li 3 PS 4 exhibit poor RT conductivities of 10 −5 to 10 −7 Scm −1 . These crystalline RT conductivities are in good agreement with experimental literature, although the extrapolated RT conductivity of β-Li 3 PS 4 is somewhat overestimated. 22 For each specific temperature and stoichiometry 20 glass structures are generated to enable an ensemble averaging. We consider three different stoichiometries in the analysis that cover the range from fully tetrahedral (Li 3 PS 4 ) via mixed (Li 7 P 3 S 11 ) to fully bridged tetrahedral (Li 4 P 2 S 7 ) thiophosphate moieties. These three dominant anion subunits are depicted in Figure 3    We can probe the phase transition quantitatively by studying the radial distribution functions (RDFs) of the sulfur sublattice as a function of simulation temperature (Figure 2). The β-and γ-phase share a HCP (hexagonal close-packed) sulfur sublattice, which is transformed to a BCC (body centered cubic) lattice in the α-phase. 17 For both sublattices, the S-S RDF displays a distinct peak at 3.4 Å, attributed to the intramolecular S-S distance.</p><p>In the HCP sublattice, a second distinct peak at 4.3 Å is observed. The latter is missing in the BCC structure. Both β-and γ-phase show the characteristic double-peak in the lowtemperature RDF, while the second peak vanishes for temperatures above 650 K. This same phase transition as also been observed in experimental studies and ab initio simulations. 23,24 The here obtained temperature between 600 and 700K for the phase transition again matches fairly well with the experimentally reported 746 K. 24 C. The Role of Anion Composition in Li 2 S-P 2 S 5 Glasses</p><p>Concluding that we can correctly describe the lithium ion dynamics and structural changes in crystalline Li 3 PS 4 we now turn to the influence of anion composition on the lithium ion conductivity in LPS glasses. As shown above, the RT conductivity generally increases with Li 2 S content of the glass material. The increasing Li-ion conductivity is partly attributable to the lithium mass percentage increase at equal densities. This larger concentration of charge carriers yields higher conductivities for similar diffusion coefficients, accounting for an increase in conductivity of ∼ 30%. As this is much less than the above described rough tripling of the conductivity, we suspect the different anion compositions in the sulfur sub-lattice to be another, dominant factor.</p><p>Existing data on the origin of ion conductivity suppression by the anion lattice is quite ambiguous. For example, experimental studies report a strong conductivity suppression by P 2 S 4− 6 , attributed to meso-scale precipitation of the non-conducting Li 2 P 2 S 6 phase. 25,26 On the contrary, density of state calculations report that P 2 S 4− 7 should suppress ion conduction at the atomic scale. 27 The charge transfer along the covalent bond between the P and the bridging S lowers the positive partial charge of the P centers, which supposedly attracts Li + ions to the P 2 S 4− 7 anions more strongly than the other thiophosphate anions. These are just two illustrative examples discussed as possible origins of ion conductivity suppression by the anion lattice.</p><p>First, we analyse the anion composition at different temperatures for all three stoichiometries. Violin plots depicting the building block distributions at different temperatures within the structure ensembles are displayed in Figure 3. For the Li 3 PS 4 glass, the simple PS 3− perature range. Hypo-thiodiphosphate P 2 S 4− 6 occurs only in small concentrations ≤ 10 at.% and shows no strong temperature dependence. Up to 25 at.% of pyro-thiodiphosphate P 2 S 4− 7 occur at the lower temperature but gradully disappear between 600 and 700 K. These anion ratios are in agreement with experimental ratios found for Li 3 PS 4 . 28 In both, the Li 7 P 3 S 11 and Li 4 P 2 S 7 glasses, the P 2 S 4− 6 content instead increases between 400 and 700 K, even though the increase is not too pronounced in comparison to the width of the distribution in the ensemble. The found P 2 S 4− 6 contents in Li 7 P 3 S 11 and Li 4 P 2 S 7 are slightly higher compared to experimental data.   In order to explore the effect of the anion lattice motion on the Li-ion conductivity we either constrain the sulfur positions, or the phosphorous positions, and compare the Li ion conductivity obtained within MD simulations with these two frozen lattices to the unconstrained Li-ion conductivity. As seen in the right panel of Figure 4, the Li-ion conductivity decreases for all three glass stoichioemetries for both frozen lattices. However, while in the case of frozen phosphorus we observe only a slight decrease in conductivity, freezing the sulfur degrees of freedom reduces the conductivity by approximately two orders of magnitude. We observe the largest decrease of the conductivity for the stoichiometry consisting of the highest PS 3− 4 (Li 3 PS 4 ) content, and the smallest change for the lowest PS 3− 4 content (Li 4 P 2 S 7 ). This suggests that the motion of sulfur throughout the Li-ion conduction plays a significant role.</p><p>Smith and Siegel showed that in glassy Li 3 PS 4 , lithium migration occurs via a mechanism that combines a concerted motion of lithium ions with re-orientations of PS 3− 4 anions. 29 This effect, known as the 'paddlewheel' mechanism, can directly attribute the increasing Li-ion conductivity with increasing PS 3− 4 content. So far, the paddlewheel effect has only been shown in Li 3 PS 4 , but our results confirm this effect occurs as long as PS 3− 4 is present.</p><p>Hence, the conductivity of Li 4 P 2 S 6 and Li 4 P 2 S 7 decrease as well when the sulfur lattice is frozen, but the effect is not as pronounced as in Li 3 PS 4 . As the Li-ion conductivity of all three stoichiometries is almost identical when the sulfur lattice is frozen, this actually suggests that the higher accessible volume of Li in Li 3 PS 4 arises from re-orientations of PS 3− 4 anions. Hence, both effects can not be decoupled, but rather the re-orientation of PS 3− 4 anions generates new Li sites. Together with the increased overall Li content, this thus fully rationalizes why the Li-ion conductivity increases with higher Li 2 S content.</p><!><p>All of the herein described effects, collective Li-ion motion of crystalline Li 7 P 3 S 11 , phase transitions of crystalline Li 3 PS 4 , and the conductivity/anion-composition relation in glassy LPS, could not be studied before by a single interatomic potential, preventing the relative identification of trends and common origins. While not only this can now be achieved by our machine learning surrogate model, the general structure of the training protocol furthermore allows for a variety of extensions, including additional selection criteria, 19,30 using an electrostatic baseline in the model, 31 doping with transition metals, and modeling of solid/solid interfaces. 32,33 We correspondingly see much prospects in the use of ML potentials to further elucidate atomic scale processes in complex battery materials.</p>