Ahmed H. Bibya,
Basant A. Alib and
Charles B. Musgrave*bcd
aMaterials Science and Engineering Program, University of Colorado Boulder, Boulder, CO 80309, USA
bDepartment of Chemical and Biological Engineering, University of Colorado, Boulder, Colorado 80303, USA. E-mail: charles.musgrave@utah.edu
cDepartment of Chemical Engineering, University of Utah, Salt Lake City, Utah 84112, USA
dDepartment of Materials Science and Engineering, University of Utah, Salt Lake City, Utah 84112, USA
First published on 22nd May 2025
Perovskite-type electrolytes (PTEs) exhibit high room-temperature ionic conductivity (∼1 mS cm−1) but suffer from poor stability against metal anodes, limiting their application in metal all-solid-state batteries. We investigated alternatives to the prototypical PTE, lithium lanthanum titanate (LLTO), to identify PTE compositions with improved stability against metal-anodes and electrolyte performance. We evaluated composition–stability relationships of PTEs using DFT calculations of thermodynamic and electrochemical stabilities, diffusion barriers, and mechanical properties. Substituting Ti with Zr (LLZrO) and Hf (LLHfO) resulted in promising metastability against metallic Li anodes, although their diffusion barriers (0.188–0.576 eV) generally exceed those of LLTO. Unlike LLZrO, which lacks mechanical stability, LLHfO is predicted to be mechanically stable with balanced elasticities, high hardness, good ductility, low anisotropy, and the potential to suppress dendrite formation. These characteristics make LLHfO a promising Li-metal-compatible alternative to LLTO with strong overall performance metrics. Analogous Na-PTEs are stable against Na-metal; however, their high Na+ diffusion barriers (>0.85 eV) must be lowered to achieve practical ionic conductivity.
Broader contextThe research situates itself within the effort to enhance the performance of metal all-solid-state batteries by addressing the stability limitations of perovskite-type electrolytes (PTEs). While conventional electrolytes like lithium lanthanum titanate (LLTO) offer high ionic conductivity at room temperature, their instability against metal anodes restricts their wider application. By leveraging density functional theory (DFT), the study explores alternative PTE compositions that could offer improved compatibility with lithium and sodium metal anodes. In particular, substituting titanium with zirconium (LLZrO) and hafnium (LLHfO) reveals that although these alternatives exhibit higher diffusion barriers compared to LLTO, LLHfO shows promising mechanical stability—combining high hardness, good ductility, and low anisotropy—to potentially suppress dendrite formation. For sodium-based batteries, however, the high diffusion barriers remain a critical challenge, indicating the need for further compositional or defect engineering. Overall, this systematic investigation not only identifies LLHfO as a promising candidate for lithium-metal batteries but also lays the groundwork for future modifications aimed at optimizing PTEs for high-performance energy storage applications. |
Since its first synthesis by Inaguma et al. in 1993, lithium lanthanum titanate (LLTO) has become the archetypal perovskite-type electrolyte (PTE).13,16 Despite its high room-temperature bulk ionic conductivity (∼1 mS cm−1),17 LLTO suffers from poor stability against Li metal anodes, limiting its use as an SSE in MASSBs.18 The general formula of LLTO is Li3xLa2/3−x□1/3−2xTiO3, where □ represents A-site vacancies. PTEs are typically aliovalent substitutional solid solutions with A-site deficiencies, formed by substituting a lower-valent A-site cation (e.g., Sr2+) with a higher-valent one (e.g., La3+). The resulting A-site vacancies preserve charge neutrality via ionic compensation.19 To enable vacancy-mediated ion diffusion, a fraction of the A-site vacancies accommodates Li+, while the remainder serve as hopping sites for Li+.7 The transition state (TS) for Li+ diffusion is the 3c site, which is located between neighboring A-sites and consists of four oxygen ions in the perovskite octahedral channel.8 The composition of PTE takes the general formula A′mA′′n□oBX3, where A′, A′′, □, B, and X are the working-ion, large-ionic-radius cation, A-site vacancy, small-ionic-radius cation, and anion, respectively.21 m, n, and o represent the concentrations of the working-ion, the large-ionic-radius cation, and the A-site vacancy, respectively. For a monovalent working-ion, compositions can be categorized as working-ion-poor or working-ion-rich. Working-ion-poor compositions typically exhibit lower and more anisotropic ionic conductivity, while working-ion-rich compositions provide higher, isotropic ionic conductivity.20
Each crystallographic site in the SSE, along with the overall composition, plays a distinct role in performance by influencing geometric constraints, ion redox activity, the diffusion energy landscape of the working-ion, and the local bonding environment.4,22 Achieving electrochemical stability in PTEs requires avoiding multivalent elements that are susceptible to reduction or oxidation.17 Enhancing ionic conductivity in PTEs is guided by two key design strategies: first, the working-ion octahedral diffusion channels can be expanded by incorporating large-radius cations, ensuring that the A′′-site cation is larger than the B-site cation, which is critical for maintaining perovskite structural stability.4,20 Second, interactions at the hopping TS can be reduced by strengthening the B–X bond, which weakens the A′–X interaction,20,24 and selecting X-site elements with high polarizability.4 However, because highly polarizable X-site elements are often multivalent, their use must be carefully considered. To improve PTE performance, previous work leveraged several strategies including tuning Li+ concentrations in LLTO,17 incorporating polymers,17,24 and doping or alloying elements at the A, B, and X sites in both LLTO and other PTEs.17,20 However, these efforts improved ionic conductivity, but to date, no PTE has demonstrated stability against metal anodes. Notably, the search for pristine, metal-stable alternatives to LLTO has largely been neglected. Instead, most studies have focused on directly doping or alloying non-LLTO PTEs, often without a comprehensive evaluation of SSE performance of their pristine compositions.
In this work, we aim to identify pristine PTE compositions that are stable against metal anodes and to elucidate key factors governing SSE performance. To this end, we studied how PTE composition influences stability against metal anodes. To benchmark performance, we selected PTE compositions composed of elements with relatively stable oxidation states, similar to those found in LLTO. For the A′-site, we primarily focused on Li+ and extended our investigation to Na+ in cases where the composition showed stability with Li+. For the A′′-site, we selected La3+, the largest trivalent ion with a relatively stable oxidation state. For the B-site, we examined two groups of elements in addition to the benchmark Ti case. First, we investigated Zr and Hf, which are in the same group as Ti in the periodic table and share similar characteristics. Notably, Hf4+ and Zr4+ exhibit more stable oxidation states, larger ionic radii, and lower electronegativities than Ti, all of which are anticipated to strengthen the B–X bonds.25–27 Second, we examined Sn and Pb, both p-block (post-transition) metals with comparatively large ionic radii. Sn has a relatively stable 4+ oxidation state,25–27 while Pb is a well-known B-site element for a variety of applications.28 We considered the chalcogens O, S, Se, and Te as potential X-site anions, assuming each to adopt the 2− oxidation states. Although heavier chalcogens (S, Se, and Te) can exhibit multiple oxidation states, they may enhance ionic conductivity if their 2− oxidation state is maintained.4 Sulfur, in particular, enables the high-energy-density anion-redox chemistry that underpins low cost metal–sulfur batteries.29
For the studied PTE compositions, we first evaluated a broad set of compositions (Fig. 1) to determine their structural, thermodynamic, and electrochemical stabilities. We also calculated electronic bandgaps to assess electronic conductivity. Finally, we conducted an in-depth analysis (Fig. 1) by computing the diffusion barriers, mechanical properties, and Projected Crystal Orbital Hamilton Population (pCOHP) of promising candidates. Detailed structural models and computational details are provided in the ESI.† We analyzed two distinct regimes: a computationally efficient low working-ion concentration (lc) model (A′0.125A′′0.625□0.25B1X3) for compositional screening, and a high working-ion concentration (hc) model (A′0.3125A′′0.5625□0.125B1X3) for in-depth analysis.
We evaluated the perovskite structural stability using the Goldschmidt tolerance factor (t, eqn (S1)†), the octahedral factor (μ, eqn (S2)†), and the machine-learned tolerance factor τ (eqn (S3)†), which all utilize Shannon ionic radii. Additionally, we used the weighted average effective ionic A-site radius (, eqn (S4)†) to assess A-site size compatibility. Stability was defined within the ranges 0.75 < t < 1, μ > 0.414, and τ < 4.18.23,25,30–32 Table 1 shows that all studied compositions were structurally stable, except those containing heavier chalcogens (S, Se, and Te) at the X-site. This instability is attributed to their large ionic radii, which lower μ below 0.414 and raise τ above 4.18, thereby destabilizing the octahedral framework.
Structure | t | μ | τ | Ef (eV per atom) | ED (eV per atom) | Ew (eV) | Eg (eV) | Comment |
---|---|---|---|---|---|---|---|---|
a Denotes metastable electrochemical stability, indicating that the reduction potential onset is metastable. | ||||||||
Step 1: Evaluating B-site elements | ||||||||
lc-LLTO | 0.98 | 0.43 | 3.86 | −3.50 | −5.50 | [1.53, 3.45] | 3.69 | Benchmark |
lc-LLZrO | 0.92 | 0.51 | 3.87 | −3.71 | −5.69 | [0, 3.19]a | 5.27 | Passed |
lc-LLHfO | 0.93 | 0.51 | 3.85 | −3.80 | −5.92 | [0, 3.28]a | 5.68 | Passed |
lc-LLSnO | 0.94 | 0.49 | 3.83 | −2.51 | −4.58 | Unstable | 3.06 | Excluded |
lc-LLPbO | 0.90 | 0.55 | 4.03 | −2.00 | −4.61 | [2.95, 3.20] | 0.73 | Excluded |
Step 2: Evaluating X-site elements | ||||||||
lc-LLTS | 0.93 | 0.33 | 4.59 | −1.78 | −4.84 | Unstable | 1.33 | Excluded |
lc-LLZrS | 0.89 | 0.39 | 4.48 | −1.99 | −5.03 | Unstable | 2.23 | Excluded |
lc-LLHfS | 0.89 | 0.39 | 4.47 | −2.00 | −5.31 | Unstable | 2.45 | Excluded |
lc-LLTSe | 0.92 | 0.31 | 4.82 | −1.59 | −4.45 | Unstable | 0.93 | Excluded |
lc-LLTTe | 0.90 | 0.27 | 5.20 | −1.12 | −4.19 | Unstable | 0.01 | Excluded |
Step 3: Evaluating Na+ as a working-ion (A′) | ||||||||
lc-NLTO | 0.98 | 0.43 | 3.84 | −3.50 | −5.49 | [0, 3.03]a | 3.55 | Passed |
lc-NLZrO | 0.93 | 0.51 | 3.83 | −3.72 | −5.69 | [0, 2.68] | 5.23 | Passed |
lc-NLHfO | 0.93 | 0.51 | 3.81 | −3.80 | −5.94 | [0, 2.26] | 5.63 | Passed |
Due to the disordered solid-solution nature of PTEs, we employed Pymatgen's OrderDisorderedStructureTransformation module to generate symmetry distinct atomic configurations within the lc and hc models.33 We applied the Ewald summation (Coulomb) method to calculate electrostatic energies to rank configurations.34 DFT calculations were performed on the ten lowest-energy configurations (LECs) to identify the most probable configuration for the lc and hc compositions. Formation and decomposition energies were then calculated using the DFT energies of the LEC's, Pymatgen's PhaseDiagram module,33 and data from the materials project (MP) repository.35 Results are reported in Table 1. All studied compositions exhibited negative decomposition and formation energies, indicating thermodynamic stability against decomposition into competing compounds and elemental phases.
We defined the electrochemical stability window (Ew) as the potential range over which a given phase has decomposition energies ≥−50 meV per atom, as described in eqn (S8)† and illustrated in Fig. 2(a). This threshold accounts for both the DFT margin of error affecting calculated formation energies and the metastability arising from kinetic stabilization effects.36 Ideally, the electrochemical stability window of a given composition should span a wide potential range relative to the metal anode reference. It should start with a reduction potential onset at 0 V to ensure stability against metal anodes and extend to a high oxidation potential onset (i.e., >3 V) to support high-voltage cathode chemistries. Under grand canonical ensemble conditions relevant to battery environments, the computed Ew values in Table 1 indicate that six compositions are largely unstable, exhibiting low thermodynamic favorability relative to decomposition into competing phases. Notably, LLZrO, LLHfO, and NLTO (with Na+ as the working-ion) exhibit reasonable metastability against metal anodes along with high oxidation potential onset, as shown in Fig. 2(a). NLZrO and NLHfO also showed stability against metal anodes, as shown in Fig. S1(b).† To assess the electronic insulation, a key requirement of SSEs, we used the HSE06 hybrid DFT functional to compute electronic bandgaps.37 Table 1 shows that only four compositions had bandgaps <2 eV, with LLTTe exhibiting a near-zero bandgap, indicative of poor electronic insulation. During compositional screening, B-site substitution with Zr and Hf combined with O on the X-site yielded promising results for both Li+ and Na+ compositions, as summarized in Table 1. The hc composition results closely mirror those of the lc variants during compositional screening (see Tables 1 and S1†).
We calculated the working-ion diffusion barriers for the stable lc and hc compositions using the climbing-image solid-state nudged elastic band (CI SS-NEB) method (see ESI†) and report the results in Fig. S2(a–d).† As a solid solution, PTE exhibits configurational disorder. Consequently, working-ions in PTEs diffuse along site-to-site trajectories with asymmetric activation barriers that strongly depend on the local atomic arrangement surrounding the diffusing ion. This leads to direction-dependent hopping barriers. To account for this variability, we employed the kinetically-resolved activation (KRA) barrier, which is defined as the average of all the diffusion barriers referenced to the end points of a given ion hop.38,39 The results are summarized in Fig. 2(b) and Table S2.† Generally, for disordered materials, diffusion barriers from a single configuration do not yield quantitative predictions of ionic conductivity, which requires extensive sampling and statistical mechanics analysis.38 Furthermore, the uncertainty of the DFT-calculated diffusion barriers is ∼0.1 eV.40 Nonetheless, diffusion barrier calculations provide valuable information about the relative ion transport kinetics in PTEs. In addition to their promising metastability against metallic Li anodes, LLZrO and LLHfO exhibited Li+ diffusion barriers between 0.188–0.576 eV, only moderately higher than those of LLTO (0.040–0.311 eV). This implies a lower room-temperature ionic conductivity, potentially by an order of magnitude, assuming Arrhenius behavior. In contrast, the analogous diffusion barriers for Na+ exceed 0.85 eV. Therefore, despite their favorable metastability against Na anodes, these PTEs are unlikely to be viable for Na all-solid-state batteries without further modifications to lower their diffusion barriers. The high barrier is attributed to the large ionic diameter of Na+ (2 Å), which is nearly twice the size of the steric opening at the TS (bottleneck) in both lc and hc structures (see Table S2†).25
Mechanical properties are crucial in SSEs for safe, durable and reliable battery handling and operation. The stress–strain method, based on the generalized Hooke's law, was employed to calculate the elastic constants using VASPKIT.41 The LECs for LLTO and LLHfO were mechanically stable, satisfying all the Born stability criteria, as described in the ESI.† However, for LLZrO, the LEC showed a negative C55 and one negative elastic tensor eigenvalue, suggesting either a structural distortion or configurational shift. We instead evaluated the second most stable configuration that has an energy 3 meV per atom above the LEC. This was more tractable than attempting to predict a mechanically stable distorted phase and found that it satisfied the mechanical stability criteria. The mechanical properties of the three studied compositions, reported in Table 2, were used to assess the elastic stability of PTEs. High bulk modulus (B), Young's modulus (E), and Debye temperature increase SSE stiffness and reduce compressibility, though excessively high values may deteriorate ionic conductivity.42–44 Accordingly, LLHfO and LLZrO exhibited more favorable elastic properties than LLTO due to their moderate and balanced stiffness.45,46 According to the Pugh criterion (B/G > 1.75, where G is the shear modulus), LLTO, LLZrO, and LLHfO are all ductile,47 making them robust under battery operating conditions and manufacturing-friendly. According to Monroe and Newman's linear elasticity model, G of an SSE must be at least twice that of lithium metal to suppress dendrite propagation.48,49 All three compositions have higher computed values of G than that required to suppress dendrite formation, as discussed in the ESI.†50 The computed Poisson ratios (v) predict that all three compositions are highly rigid. Vickers hardness reflects the strength of the bonds of a material and serves as an indicator of its structural integrity under mechanical stress. We found that LLHfO and LLZrO have higher Vickers hardness than LLTO and are thus likely to be more mechanically stable against stress during battery operation. In SSEs, anisotropy affects ionic conductivity and mechanical stability. The universal anisotropy number (UAN) quantifies the degree of anisotropy, where a UAN greater than 1 indicates significant directional dependence.51 LLTO and LLHfO's computed UANs of 0.14 and 0.11 suggest more uniform mechanical behavior and ionic conductivity than LLZrO, which has a UAN of 0.99. LLHfO shows the best combination of isotropy and stiffness, whereas LLZrO, while softer (lower E, B) than LLTO, has a significantly higher anisotropy, which could influence crack propagation or ionic percolation in a real polycrystal.
Mechanical property, unit | LLTO | LLZrO | LLHfO |
---|---|---|---|
Bulk modulus (B), GPa | 277.54 | 115.60 | 135.92 |
Young's modulus (E), GPa | 251.46 | 150.21 | 186.44 |
Shear modulus (G), GPa | 93.20 | 58.52 | 73.32 |
Poisson's ratio (v) | 0.35 | 0.28 | 0.27 |
Pugh's ratio (B/G) | 2.98 | 1.98 | 1.85 |
Vickers hardness, GPa | 4.88 | 6.70 | 8.91 |
Universal anisotropy number (UAN) | 0.14 | 0.99 | 0.11 |
Debye temperature, K | 609.40 | 438.80 | 427.40 |
We used Lobster52 to perform pCOHP analysis on the Li–O and transition metal (TM)–O bonds in the studied compositions to examine the interaction between the working-ion (A′) and the X-site anion. The analysis was performed on both the initial (Li@A-site) and TS (Li@TS) configurations, extracted from the NEB images shown in Fig. S3.† The Li–O peaks are considerably weaker than the TM–O peaks, confirming their greater ionicity and weaker bonding. Notably, the Li–O anti-bonding peak for LLTO lies higher than those of LLZrO and LLHfO, indicating stronger Li–O bonding in LLTO than in LLZrO and LLHfO. For Zr–O bonds, the presence of anti-bonding interactions below the Fermi level highlights their tendency to dissociate under stress or elevated temperatures. This observation aligns with the results from mechanical stability analyses. Conversely, the small bonding peaks in the anti-bonding contribution of Li–O bonds in LLTO are negligible and positioned far above the Fermi level, indicating minimal influence on its bonding nature. Additionally, the Integrated Crystal Orbital Hamilton Population (ICOHP) was used to quantify the Li–O and TM–O bond strengths, with higher (−ICOHP) values indicating stronger bonds. As shown in Fig. 3(a), the results indicate that Zr and Hf form stronger bonds with oxygen due to their lower electronegativity relative to Ti, which slightly weakens the Li–O bonds in LLZrO and LLHfO relative to LLTO. To examine the Li–O bonding changes upon diffusion, the average −ICOHP values for Li–O bonds of the diffusing Li+ were calculated at the Li@TS configurations. As shown in Fig. 3(b), the Li–O bonds of LLZrO and LLHfO are slightly weaker than those in LLTO, reducing the interaction of oxygen with the diffusing Li+, as intended. Nonetheless, the persistence of high diffusion barriers in LLZrO and LLHfO despite weaker Li–O bonding (A′–X interaction) indicates that factors such as the lattice deformation energy may dominate diffusion kinetics, suggesting that this design principle may need to be revisited in future studies.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5eb00089k |
This journal is © The Royal Society of Chemistry 2025 |