Unlocking zeolite-like structures as a new family of interstitial oxide ion conductors: insights into carrier trapping, collective local distortion, and correlated disorder

Abstract

Zeolites have emerged as indispensable materials for applications ranging from catalysis and separation to adsorption and ion exchange, owing to their uniquely ordered porous architectures composed of well-defined channels and cavities. Inspired by the promising oxygen defect tolerance observed in various open structural frameworks, herein, we have developed a new zeolite-like feldspar structure, A₂B₂C₂O₈ with 4- and 8-membered rings, in the case of Sr_1−xLa_xGa₂Ge₂O_8+0.5x, as a new family of interstitial oxide ion conductors due to its open structural framework for accommodating and transporting interstitial oxide ions. Average structural analysis revealed that the interstitial oxygen occupies the centers of 4-membered rings, existing in a coordination equilibrium quasi-free state that confers high mobility; however, this contrasts with the experimentally observed low mobility. Local structural analysis further revealed that the local collective distortions in GaO₄ and GeO₄ tetrahedra, together with the correlated disorder of interstitial oxygen coordinated with Ge, hidden within the average structure, are critical to interstitial oxygen mobility from the 4-membered ring to the 8-membered ring. Our findings demonstrate zeolite-like structures as a new family of interstitial oxide ion conductors, offering new insights into the intricate interplay between oxide ion mobility, collective distortions, and correlated disorder at the local scale.

---

1. Introduction

Solid oxide fuel cells (SOFCs) are emerging as a transformative platform for electrochemical energy conversion, directly converting the chemical energy of fuels into electricity with high efficiency, fuel flexibility, and environmental friendliness.^1–3 Central to their functionality are oxide ion conductors, which act as electrolytes for transporting oxide ions in SOFCs and also have environmental applications in a wide range of devices such as oxygen separation membranes, pumps, and sensors.^4–7 In order to reduce the cost and enhance long-term stability for these practical applications, it is highly desirable to lower their operating temperatures in the future. The transport of oxide ions in oxide ion conductors, which is primarily governed by oxygen vacancies or interstitial defects that provide pathways for ionic migration, plays a critical role in determining the operating temperatures. However, developing high-performance oxide ion conductors operating at low temperature through tailoring the compositions, defects and structures remains a great challenge, particularly in balancing high oxide ion conductivity with chemical and mechanical stability under operating conditions.^8,9

To address these challenges, various structural prototypes have been developed successively as candidates for oxide ion conductors. These structures are either linked polyhedral networks or isolated polyhedral anion structures, which include the three-dimensional fluorite,¹⁰ pyrochlore,¹¹ perovskite,¹² langasite structures,¹³ as well as the low-dimensional structures containing tetrahedral units, such as LAMOX (La₂Mo₂O₉),¹⁴ scheelite,¹⁵ apatite,¹⁶ melilite-type structures,¹⁷ and borate-based oxide ion conductors, (Gd/Y)_1−xZn_xBO_3−0.5x.¹⁸ Summarizing the structural features and oxide ion conductivities of the above oxide ion conductors reveals that materials with open structural frameworks show remarkable potential for accommodating oxygen defects with high mobility. Among these open frameworks, the melilite and langasite families are notably two classical open structural types owing to their remarkable defect-tolerant properties.^13,17,19,20

The melilite family (A₂B(B′₂O₇)) features a two-dimensional (2D) layered structure composed of corner-sharing tetrahedral anionic layers (BO₄ with 4 bridging oxygens and B′O₄ with 3 bridging and 1 terminal oxygen) that form 5-membered rings accommodating the large A cations (Fig. 1a and b).²¹ By tuning the cationic composition and introducing excess positive charge, interstitial oxide ions can be stabilized within the pentagonal rings between A cations in melilite structures. Notably, the typical melilite La_1.54Sr_0.46Ga₃O_7.27 (ref. 17) demonstrates exceptionally high interstitial oxide ion conductivity (∼0.02 S cm⁻¹ at 600 °C), surpassing those of most commonly used electrolytes below 600 °C, with the interstitial oxygen defect migrating via a concerted “knock-on” mechanism involving interstitial and lattice oxygen.²² In contrast, the langasite family adopts a trigonal structure (A₃BC₃D₂O₁₄, Fig. 1c and d),^23–25 composed of layers of four-coordinated CO₄ tetrahedra (with 4 bridging oxygens) and three-coordinated DO₄ tetrahedra (with 3 bridging and 1 terminal oxygen), interconnected by BO₆ octahedra. This topology imparts exceptional compositional flexibility to the langasite structure. Rosseinsky et al.¹⁹ showed that the langasite structure can accommodate high concentrations of interstitial oxygen, although the interstitial oxygens are trapped within the 6-membered rings with limited mobility. These structural commonalities highlight the potential of open structural frameworks as a promising direction for developing advanced oxide ion conductors.¹³

Fig. 1 Typical open structural frameworks. (a and b) Melilite, (c and d) langasite and (e and f) zeolite-like feldspar structures viewed along the (a, b and e) [001] and (b, d and f) [100] directions.

Driven by the promise of open structural framework architectures, the zeolite-like family could emerge as a compelling candidate for advanced oxide ion conductors. Zeolites are inorganic crystalline materials constructed from corner-sharing TO₄ tetrahedra (T = Si, Al, P, etc.), forming precisely ordered microporous structures.²⁶ This unique architecture forms uniform channels and cavities that confer exceptional (hydro)thermal stability, a unique porous framework, and tunable acid sites. These features make zeolites indispensable in a wide range of applications, including catalysis, separation, adsorption, and ion exchange.²⁷ Beyond these traditional applications, zeolites have recently shown remarkable promise in energy storage systems. They were initially employed as functional layers on separators or electrodes and as additives in composite solid electrolytes for lithium-ion batteries. However, their potential in ion transport has not yet been adequately recognized.^28–30 More recently, Yu et al.³¹ reported a zeolite X (LiX) membrane serving as the solid electrolyte with superior environmental adaptability and electrochemical stability, addressing critical limitations of conventional solid electrolytes for Li–air batteries. The unique open structural features and outstanding properties in ion transport and chemical stability indicate that zeolite-based materials could be developed into oxide ion conducting electrolytes for advancing the SOFC technology. So far, the research into their role in oxide ion transport remains in its infancy.

Building on this concept, we explore a zeolite-like feldspar structure with the general formula AB₂C₂O₈, where A represents an alkaline earth metal (e.g., Ca²⁺ or Sr²⁺ or Ba²⁺), B is a trivalent metal (e.g., Al³⁺ or Ga³⁺), and C is a tetravalent metal (e.g., Si⁴⁺ or Ge⁴⁺) (Fig. 1e and f). This zeolite-like structure, composed of 4- and 8-membered rings,³² holds great promise for accommodating and transporting interstitial oxide ions. Herein, the zeolite-like case, strontium feldspar SrGa₂Ge₂O₈, was identified as a potential interstitial host due to its open structural framework and the variable coordination numbers of Ga and Ge that facilitate oxygen defect accommodation and migration. The interstitial oxygen defects were created by La-doping at the Sr site in SrGa₂Ge₂O₈ and the average structure of Sr_1−xLa_xGa₂Ge₂O_8+0.5x containing interstitial oxygen was determined by neutron powder diffraction (NPD). X-ray and neutron pair distribution function (PDF) analysis based on reverse Monte Carlo (RMC) modeling, combined with 1D and 2D ⁷¹Ga solid state NMR technique, extended X-ray absorption fine structure (EXAFS) analysis, and molecular dynamics (MD) simulations, further revealed that the local collective distortion and correlated disorder in the local structure, hidden within the average structure, are intricately associated with interstitial oxygen mobility. This study provides critical insights into the potential of the zeolite-like family as a new class of oxide ion conductors.

2. Experimental section

2.1 Synthesis method

Polycrystalline samples of La_xSr_1−xGa₂Ge₂O_8+0.5x were synthesized using a conventional high-temperature solid-state reaction method. High-purity raw materials SrCO₃ (99.9%, Aladdin), La₂O₃ (99.99%, Aladdin), Ga₂O₃ (99.9%, Aladdin), and GeO₂ (99.99%, Aladdin) were weighed according to the desired stoichiometries and mixed thoroughly in ethanol using an agate mortar and pestle. The mixed powders were pressed uniaxially into pellets with a diameter of 10 mm and a thickness of approximately 1 mm. The pellets were sintered at 1100 °C for 20 h in air, with heating and cooling rates of 5 °C min⁻¹, obtaining a relative density of approximately 90%.

2.2 Characterization

Phase formation in the samples was examined by powder X-ray diffraction (XRD) using a Panalytical X'Pert PRO diffractometer with a PIXcel 1D detector using Cu Kα radiation. Variable-temperature XRD patterns were collected on this instrument equipped with an Anton Paar HTK 1200N high temperature attachment. Time-of-flight (TOF) neutron powder diffraction (NPD) data were collected from a general purpose powder diffractometer (GPPD) at the China Spallation Neutron Source (CSNS).³³ The total scattering neutron pair distribution function (nPDF) data were collected on the BL-1B NOMAD beamline at the Spallation Neutron Source (SNS), Oak Ridge National Laboratory. Average lattice fitting was conducted using PDFgui software.³⁴ Reverse Monte Carlo simulations for the nPDF data were carried out using RMCprofile software³⁵ combining the atomic movements and swapping and utilizing a 8 × 8 × 9 (77.19 Å × 73.63 Å × 77.21 Å) large supercell containing 30 528 atoms. The XRD and NPD data were analyzed by the Rietveld method, which was carried out using Topas-Academic software³⁶ and Jana 2006 software.³⁷

Impedance spectroscopy measurements were conducted using a Solartron 1260 frequency response analyzer over the 10⁻¹ to 10⁷ Hz frequency range within 300–1000 °C. Impedance data analysis was carried out using ZView software.³⁸ The oxygen transport number was determined by electromotive force (EMF) measurements on oxygen concentration cells of O₂//1%O₂ and O₂//5%H₂, within the temperature range of 500–1000 °C. The theoretical EMF values of the oxygen concentration cells were calculated using the Nernst equation.

Rotation electron diffraction (RED) data were collected on a JEOL-2100 electron microscope operating at 200 kV using a single-tilt tomography sample holder. The RED data were processed using REDp processing software.³⁹ X-ray photoelectron spectroscopy (XPS) measurements were performed on an AXIS Supra X-ray photoelectron spectrometer using 200 W monochromated Al Kα radiation. Morphological images and X-ray energy-dispersive spectroscopy (EDS) data were collected on a FEI Tecnai-G2 field emission transmission electron microscope (TEM) equipped with a Li/Si EDX detector. Ultrasoft X-ray absorption spectra (XAS) were collected on the 4B9B beamline of the photoelectron spectroscopy station at the Beijing Synchrotron Radiation Facility (BSRF).

Solid-state ⁷¹Ga nuclear magnetic resonance (NMR) spectra were acquired on a Bruker AVANCE NEO 850 (20 T) NMR spectrometer with a 1.3 mm magic angle sample spinning probe. All spectra were collected at a rotor frequency of 60 kHz. The 1D quantitative CT-selective ⁷¹Ga spectra were measured using a Hahn echo experiment with acquisition at the echo top using an interpulse delay of 20 μs. Triple-quantum ⁷¹Ga MQ-MAS 2D spectra were collected using a three-pulse sequence, a repetition delay of 0.5 s, and rotor-synchronized sampling of the indirect dimension ν_r = 240 kHz with 128 t₁ increments. The NMR parameters were fitted by means of the DMfit program.⁴⁰

2.2.1 Atomistic static lattice simulations. Atomistic static lattice simulations on defect formation in SrGa₂Ge₂O₈ were carried out using the General Utility Lattice Program (GULP) on the basis of the interatomic potential approach.⁴¹ The Buckingham potential function⁴² was used to model interaction between ions, with the shell model to describe the electronic polarizability for the structure modeling of SrGa₂Ge₂O₈. The potential parameters used for the atomistic simulation are listed in Table 1, among which the parameters for La–O, Sr–O, Ga–O, Ge–Ge, and O–O were taken from the previous studies.^22,43 The potential parameters were slightly modified through the relaxed fitting procedure implemented in the GULP package based on the initial parameters from the study of Tealdi et al.,²² to better reproduce the crystal structure of SrGa₂Ge₂O₈ via refining the A value while keeping the C value fixed. The optimized structural parameters showed good agreement with experimental data (Table S1†), validating the reliability of the interatomic potential model for subsequent static lattice and molecular dynamics simulations. The formation energies of interstitial oxide ions were calculated based on the appropriate combination of dopant, defect, and lattice energies of the binary oxides, as described by the following defect equation: .

Table 1 Buckingham interatomic potential and shell model parameters for the SrGa₂Ge₂O₈-based materials

Interaction	A (eV)	ρ (Å)	C (eV Å^6)	Y (e)	k (eV Å^−2)
La–O	4579.23	0.30437	0	3	99 999
Sr–O	1400.10150	0.357400	0	1.33	21.53
Ga–O	1241.55690	0.302428	3.323663	0	0
Ge–O	1482.37220	0.324420	10.64850	0	0
O–O	20 784.8710	0.122291	48.06810	−2.869	74.92

---

2.2.2 Molecular dynamics (MD) simulations. Interstitial oxygen migration in La-doped SrGa₂Ge₂O₈ was investigated using interatomic-potential-based MD simulations performed with the DL_POLY code in the NVT ensemble.⁴⁴ For the undoped composition, a 3 × 3 × 3 supercell model was constructed, comprising 1872 atoms (144 Sr, 288 Ga, 288 Ge, and 1152 O atoms). For the La-doped composition, a simulation box consisting of a 4 × 4 × 4 unit cell was generated, comprising a total of 3360 atoms: 64 La, 192 Sr, 512 Ga, 512 Ge, and 2080 O atoms, corresponding to the composition La_0.15Sr_0.85Ga₂Ge₂O_8.075. The systems were equilibrated at 1 atm and various temperatures ranging from 600 to 1200 °C for 135 ps in the NVT ensemble (constant number of particles N, volume V, and temperature T controlled using a thermostat). MD data analysis was carried out by using the Visual Molecular Dynamics (VMD) package,⁴⁵ and mean square displacement (MSD) data were exported through the nMoldyn code.⁴⁶ Oxygen diffusion coefficients were determined by linear fitting of the MSD versus time curves.

3. Results and discussion

3.1 Sr_1−xLa_xGa₂Ge₂O_8+x/2 solid solution

The XRD patterns of La-doped compositions Sr_1−xLa_xGa₂Ge₂O_8+x/2 (0 ≤ x ≤ 0.15, Fig. 2a) indicate the formation of a pure phase across the entire doping range, with no impurity phases detected. Rietveld refinement of the XRD patterns for the maximum La-doping concentration (x = 0.15, Fig. 2b) further confirms the pure phase formation. The refined cell parameters decrease gradually with increasing La content (x), attributed to the smaller ionic radius of La³⁺ (1.061 Å) compared to Sr²⁺ (1.13 Å), which causes contraction of the crystal lattice. Three-dimensional continuous rotation electron diffraction (3D cRED) patterns of Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 2d) yield a monoclinic unit cell with space group P2₁/c and cell parameters a = 9.02 Å, b = 9.31 Å, c = 8.63 (3) Å, β = 90.1°. The observed reflection conditions, including 0kl: k + l = 2n, h00: h = 2n; h0l: h + l = 2n; h00: h = 2n, indicate that the introduction of La³⁺ does not change the average crystal structure of Sr_0.85La_0.15Ga₂Ge₂O_8.075. Transmission electron microscopy (TEM) with elemental mapping (Fig. 2e) demonstrates the uniform distribution of all constituent elements in Sr_0.85La_0.15Ga₂Ge₂O_8.075, confirming the homogeneous incorporation of La³⁺ into the lattice.

Fig. 2 Solid solution of Sr_1-xLa_xGa₂Ge₂O_8+x/2. (a) XRD patterns for compositions with 0 ≤ x ≤ 0.15. (b) Rietveld refinement plot of Sr_0.85La_0.15Ga₂Ge₂O_8.075. (c) Refined cell parameters as a function of x. (d) 3D ED reconstructed reciprocal lattice of Sr_0.85La_0.15Ga₂Ge₂O_8.075 and the (0kl), (h0l) and (hk0) diffraction planes. (e) TEM elemental mapping of Sr_0.85La_0.15Ga₂Ge₂O_8.075.

3.2 Electrical and stability performances

La-doping at the Sr site in SrGa₂Ge₂O₈ can be described using the following defect equation: . Here, the substitution of La³⁺ for Sr²⁺ introduces interstitial oxygen . Alternating current (AC) impedance spectroscopy was used to study the mobility of the introduced interstitial oxygen. The typical complex impedance plots of Sr_0.85La_0.15Ga₂Ge₂O_8.075 at different temperatures (Fig. 3a and S1†) show a remarkable Warburg electrode response, indicative of ionic migration. This ionic migration behavior is also observed for the composition with x = 0.05 (Fig. S2 and S3†). This electrode response is highly sensitive to oxygen partial pressure, as evidenced in Fig. 3b, further confirming oxide ion migration within the structural framework of Sr_0.85La_0.15Ga₂Ge₂O_8.075. The oxygen transport numbers of Sr_0.85La_0.15Ga₂Ge₂O_8.075 were determined using an oxygen concentration cell method. Electromotive force measurements were conducted in oxygen concentration cells (O₂//1%O₂ and O₂//5%H₂, Fig. 3c) over a temperature range of 500–1000 °C, yielding an oxygen transport number of ∼80% with minimal electronic conduction. The conductivities of Sr_0.85La_0.15Ga₂Ge₂O_8.075 under different oxygen partial pressures (Fig. 3d) show minor p-type electronic conduction in the high oxygen partial pressure range of 1–10⁻⁴ atm. The bulk conductivities of Sr_1−xLa_xGa₂Ge₂O_8+x/2 (0 ≤ x ≤ 0.15, Fig. 3e) show that increasing the concentration of interstitial oxygen enhances oxide ion conductivity by approximately two orders of magnitude compared with the parent composition, while the activation energy decreases from 2.42(1) eV to 1.35(4) eV (Fig. 3f). Among the compositions, the x = 0.15 composition, Sr_0.85La_0.15Ga₂Ge₂O_8.075, exhibits the highest oxide ion conductivity (1.72 × 10⁻⁴ S cm⁻¹) with an activation energy of 1.35 eV.

Fig. 3 Electrical and stability performance of Sr_1–xLa_xGa₂Ge₂O_8+x/2. (a and b) Typical complex impedance plots of x = 0.15 at (a) 750 °C in air and (b) 900 °C under different oxygen partial pressures. (c) Oxygen transport numbers of x = 0.15 determined using O₂//1%O₂ and O₂//5%H₂ cells. (d) Conductivities for x = 0.15 in an oxygen partial pressure range of 1–10⁻³⁰ atm. (e) Bulk conductivities in air and (f) derived activation energies. (g) VT-XRD patterns of x = 0.15 and (h) the refined cell parameters as a function of temperature. (i) TG-DSC curves of x = 0.15.

Variable temperature XRD (VT-XRD) patterns of Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 3g, and S4†) demonstrate excellent phase stability across the temperature range of 25–1000 °C, further supported by the TG-DSC data (Fig. 3i). Notably, the phase stability is also maintained under reducing conditions (Fig. S5†). The volume thermal coefficient (TEC) value derived from VT-XRD (Fig. 3h) is α_V = 2.154 × 10⁻⁵ K⁻¹, which closely matches the TEC values of widely used cobalt-containing cathodes (2.0–2.5 × 10⁻⁵ K⁻¹).^47–49 These thermal stability and compatibility features highlight the potential of the zeolite-like feldspar structure as a promising candidate for SOFC applications.

3.3 Average structure

Combining Rietveld analysis with the charge-flipping method applied to the NPD data of Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 4a), using the SrGa₂Ge₂O₈ structural model in monoclinic symmetry (P2₁/c) with its Ga and Ge ordering arrangement,⁵⁰ revealed the presence of interstitial oxygen O_i at the (0 0 0) site (Table 2). In contrast to the parent structure (Fig. 4b), the interstitial oxygen was located at the center of 4-membered rings rather than the larger 8-membered rings, as marked by black circles in Fig. 4c. The interstitial oxygen, coordinated with two Ga and two Ge atoms, exhibits Ga–O_i and Ge-O_i bond lengths of 2.292 Å and 2.060 Å, respectively, indicative of underbonding (BVS = 1.30, Table 2) and approximate ionic bonding character (Fig. 4d). This configuration allows the oxygen to remain in a quasi-free state of coordination equilibrium, enabling high mobility. Notably, Ga and Ge adopt (4 + 1) GaO₄₊₁ and GeO₄₊₁ coordination environments. The introduced interstitial oxygen induces shifts in the binding energies of Ga-3d and Ge-3d orbitals (Fig. 4e) towards higher energy regions, reflecting the influence of the additional oxygen coordination.

Fig. 4 Average structure of Sr_0.85La_0.15Ga₂Ge₂O_8.075. (a) Rietveld refinement plot of the NPD data. (b) Parent structure in comparison with (c) crystal structure from NPD data refinement. (d) Enlarged structure around interstitial oxygen in the square ring. (e) Binding energies of Ga-3d and Ge-3d orbitals. (f) Interpretation of the oxygen K-edge XAS spectrum of the 3d transition metal oxide, and the ratio of t_2g and e_g states is 6 : 4. (g) Ultrasoft O K-edge XAS spectra in comparison with the parent composition. (h) Ga K-edge EXAFS (points) and curve fit (line), shown in k³ weighted k-space. (i) Ga K-edge EXAFS (points) and curve fit (line), shown in R-space (FT magnitude and imaginary component). (j) Ga K-edge wavelet transform EXAFS contour plot.

Table 2 The final refined structural parameters of Sr_0.85La_0.15Ga₂Ge₂O_8.075 from NPD data^a

Atom	Wyckoff site	x	y	z	Occ.	Uiso (Å^2)	BVS
a Space group: P21/c (14), a = 9.2044(3) Å, b = 9.6497(2) Å, c = 8.5795 (3) Å, β = 90.396(3)°, V = 762.01 (4) Å^3, and Z = 4.
Sr1	4e	0.8847(5)	0.4135(6)	0.246(1)	0.87(5)	0.040(1)	1.87
La1	4e	0.8847(5)	0.4135(6)	0.246(1)	0.13(5)	0.040(1)	2.16
Ga1	4e	0.0510(7)	0.1926(8)	0.5691(7)	1	0.033(1)	3.30
Ga2	4e	0.2334(7)	0.9240(9)	0.9330(7)	1	0.019(1)	3.65
Ge1	4e	0.0555(6)	0.2001(8)	0.9395(6)	1	0.041(1)	3.93
Ge2	4e	0.2463(6)	0.9197(9)	0.5567(6)	1	0.025(1)	3.75
O1	4e	0.1967(9)	0.091(1)	0.015(1)	1	0.041(1)	2.10
O2	4e	0.1908(9)	0.086(1)	0.481(1)	1	0.034(1)	2.19
O3	4e	0.120(1)	0.374(1)	0.9337(9)	1	0.031(1)	1.85
O4	4e	0.114(1)	0.371(1)	0.5652(9)	1	0.023(4)	1.88
O5	4e	0.9110(9)	0.204(1)	1.0760(9)	1	0.040(1)	2.25
O6	4e	0.9042(9)	0.193(1)	0.4310(9)	1	0.023(3)	2.21
O7	4e	−0.0014(8)	0.1359(6)	0.760(1)	1	0.035(6)	1.75
O8	4e	0.332(1)	0.9245(7)	0.743(1)	1	0.031(3)	1.75
Oi	2a	0	0	0	0.10(3)	0.038 (1)	1.30

---

To further study the oxygen coordination environments and the associated variations in electronic structure, ultrasoft O K-edge X-ray absorption spectroscopy (XAS) was employed, with Au as the reference (Fig. S6†). The schematic diagram in Fig. 4f illustrates the oxygen K-edge XAS spectrum of the 3d transition metal oxide, with the oxygen 1s core state (shown in green) at ∼534 eV binding energy. Typically, O 1s electrons are excited to occupy the e_g and t_2g orbitals. However, the presence of excess interstitial oxygen contributes to increased t_2g orbital occupancy, as indicated by a reduction in the relative peak intensities from ΔL to approximately 0.9ΔL (Fig. 4g). As a result, the composition Sr_0.85La_0.15Ga₂Ge₂O_8.075 shows a relatively higher t_2g orbital occupancy compared to the parent structure, further confirming the impact of the interstitial oxygen in modifying the electrical properties at the atomic scale.

To further investigate the structural variations associated with coordination numbers of B-site Ga, Ga K-edge extended X-ray absorption fine structure (EXAFS) measurements were performed. For the parent SrGa₂Ge₂O₈, the fitting of the Fourier-transformed (FT) k³-weighted EXAFS spectra in both k-space (Fig. S7a†) and R-space (Fig. S7b†) was carried out without including the scattering path of the interstitial oxygen, leading to an R factor of 0.9% (Table S2†). In contrast, the FT k³-weighted EXAFS spectra of Sr_0.85La_0.15Ga₂Ge₂O_8.075, shown in Fig. 3h and 4i for k-space and R-space, respectively, were fitted with the inclusion of the scattering path for the interstitial oxygen, giving an R factor of 1.7% (Table S3†). Furthermore, wavelet transform analysis of the EXAFS spectra for both the parent SrGa₂Ge₂O₈ (Fig. S7c†) and Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 4j) reveals a weak but discernible signal corresponding to the scattering path of interstitial oxygen. This observation confirms the formation of (4 + 1)-coordinated GaO₄₊₁ environments around Ga atoms, induced by the interstitial oxygen, and highlights the role of interstitial oxygen in modulating the local atomic structure.

3.4 Local structure

⁷¹Ga solid-state nuclear magnetic resonance (NMR) measurements were performed to study the local structure of Sr_1−xLa_xGa₂Ge₂O_8+x/2. The structure of SrGa₂Ge₂O₈ contains two Ga crystallographic sites with a multiplicity ratio of 1 : 1. ⁷¹Ga MAS NMR spectrum at 20 T and a MAS rate of 60 kHz of pure SrGa₂Ge₂O₈ (Fig. 5a) shows one signal with a complex line shape. In the ⁷¹Ga MQMAS experiment (Fig. 5b, c and Table 3), two signals were revealed at 197.8 and 204.0 ppm. The simulation of the 1D ⁷¹Ga spectrum yielded an integral intensity ratio of 1 : 1, in good agreement with the crystal structure. Compared with the parent composition, the La-doped compositions displayed essentially the same ⁷¹Ga resonance signals except for one additional quadrupolar signal at 40.7 ppm with a C_Q value of 8.4 MHz and an η_Q value of 0 (Fig. 5a and c). Based on literature data we assigned this resonance to six-fold coordinated gallium in gallium oxide.^51,52 Ga₂O₃ contains two Ga crystallographic sites with a multiplicity ratio of 1 : 1. The second signal from the four-coordinated gallium of Ga₂O₃ oxide was not observed due to its low-intensity, second-order quadrupolar enlargement, and the overlap of the intense SrGa₂Ge₂O₈ signals. It should be noted that the spectrum of parent SrGa₂Ge₂O₈ contains the trace of this signal from the 6-coordinate Ga. The existence of Ga₂O₃ in the samples was not detected by XRD, indicating the higher sensitivity of NMR to traces of secondary phases.

Fig. 5 Local structures of Sr_1-xLa_xGa₂Ge₂O_8+x/2. (a) ⁷¹Ga solid-state NMR spectra. (b and c) 2D MQMAS ⁷¹Ga solid-state NMR spectra of (b) parent SrGa₂Ge₂O₈ and (c) Sr_0.85La_0.15Ga₂Ge₂O_8.075. (d) Geometry optimized structure containing interstitial oxygen at (e) Ga1Ge2 and (f) Ga2Ge1 layers. (g) Small box fitting of neutron PDF data of Sr_0.85La_0.15Ga₂Ge₂O_8.075 with (h) the short-range fitting in the range of 1–10 Å, where the red arrows indicate mismatches.

Table 3 ⁷¹Ga isotropic chemical shifts (δ_iso), quadrupolar constants (C_Q), asymmetry parameters (η_Q), and integral intensities (I) obtained from the simulation of the ⁷¹Ga MAS and MQMAS spectra at 20 T of Sr_1−xLa_xGa₂Ge₂O_8+x/2

Atom	δiso, ppm (±0.2 ppm)	CQ, MHz (±0.1 MHz)	ηQ (±0.05)
Ga1	197.8	6.7	0.39
Ga2	204.0	9.8	0.27
Ga2O3	40.7	8.4	0
	200.0	11.0	0.85

---

The nearly unchanged ⁷¹Ga resonance signals in the NMR spectra of the La-doped compositions compared with the parent one suggest that locally the interstitial oxide ions do not form strong bonding with Ga. This is consistent with the NPD structural analysis (Fig. 4d), which revealed that Ga forms 4 covalent bonds (1.828 Å) with lattice oxygen and one ionic bond (2.292 Å) with interstitial oxygen. The coordination of gallium close to oxygen interstitials in La-doped compositions can thus be described as 4 + 1. Unfortunately, for this reason, probing structural variations arising from interstitial oxygen in the ⁷¹Ga NMR spectra could be hindered given the low concentration of interstitial oxide ions and the much weaker ionic bonding nature of Ga–O_i. Our conclusion is confirmed by observations in the case of La₂Ga₃O_7.5 in which the Ga1 site has 4 + 1 coordination with a chemical shift (δ_iso = 187.0 ppm) in the range corresponding to four-coordinated gallium.⁵³

Geometry optimization of the Sr_0.85La_0.15Ga₂Ge₂O_8.075 structure revealed structural features consistent with the average structure (Fig. 5d–f) and yielded a defect formation energy of ∼1.80 eV, indicating a thermodynamically favorable environment for accommodating interstitial defects within the zeolite-like feldspar structure. Although the defect formation energy is significantly lower than that (2.67 eV) observed for the typical melilite La_1.50Sr_0.50Ga₃O_7.25, this feldspar phase has a much lower capability for accommodating interstitial oxide ion defects than the melilite. Therefore, the defect formation energies cannot be directly compared to quantitatively assess the defect formation capability particularly when derived from calculations based on interatomic potential methods using empirical parameters. The interstitial defect content in the feldspar phase could be experimentally modulated further using different synthetic routes, as indicated by the melilite gallate and aluminates.^53,54

The interstitial oxygen in the average structure is in a quasi-free state of coordination equilibrium, which presumably could promote its high mobility. However, this contrasts with the experimentally observed limited interstitial oxygen mobility in Sr_0.85La_0.15Ga₂Ge₂O_8.075, which is closely related to local structural variations that are not captured in the average structure. Although spectroscopic techniques such as O-XAS, ⁷¹Ga-EXAFS, and solid-state NMR revealed some information on the local structure, to gain deeper insights into the local structures, the neutron partial distribution function (nPDF) technique was employed on Sr_0.85La_0.15Ga₂Ge₂O_8.075. Indeed, small-box fitting of the nPDF data for Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 5g) and short-range fitting within the range of 0–10 Å (Fig. 5h) revealed mismatches, indicated by red arrows in Fig. 5h. These mismatches indicate the presence of local structural disorder/distortions in Sr_0.85La_0.15Ga₂Ge₂O_8.075, which could have correlation with the interstitial oxygen migration behavior. These local structural disorder/distortions were further investigated through reverse Monte Carlo (RMC) simulations using large supercell models fitted to the PDF data.

3.5 Collective local distortion and correlated disorder in Sr_1−xLa_xGa₂Ge₂O_8+x/2

To further elucidate the local structure and its correlation with interstitial oxygen migration behavior in Sr_0.85La_0.15Ga₂Ge₂O_8.075, a large-scale big-box model was constructed using an 8 × 8 × 9 supercell (77.19 Å × 73.63 Å × 77.21 Å, Fig. 6a). The reverse Monte Carlo (RMC) method was used to analyse the spatial variation of local atomic distributions in Sr_0.85La_0.15Ga₂Ge₂O_8.075. The big box refinement of neutron PDF data of Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 6b) reveals various local structures (Fig. 6c–g) responsible for stabilizing interstitial oxygen. By tracking interstitial oxygens, we employed statistical analysis of the RMC-derived big-box model to elucidate the stabilization mechanism of interstitial oxygens. The statistical analysis reveals that all interstitial oxygens are stabilized through collective tetrahedral distortions (Fig. 6d₂–g₂) that facilitate coordination with Ge atoms (Fig. 6c and d₁–g₁), revealing the correlated disorder of interstitial oxygen coordinated with Ge. Specifically, these interstitial oxygens predominantly coordinate with one Ge atom (31.86%, Fig. 6d₁) or two Ge atoms (57.48%, Fig. 6e₁), with additional configurations involving one Ga and one Ge atom (3.15%, Fig. 6f₁) or two Ge atoms and one Ga atom (5.51%, Fig. 6g₁). In other words, the local collective tetrahedral distortions arise from the preferential coordination of interstitial oxygen with Ge atoms, accompanied by cooperative local relaxations of the surrounding atoms to accommodate the interstitial oxygen and maintain its coordination with Ge.

Fig. 6 Local structure of Sr_0.85La_0.15Ga₂Ge₂O_8.075 revealed by RMC simulations on neutron PDF data. (a) 3D atomic configuration with an approximately 77.19 Å × 73.63 Å × 77.21 Å cell for RMC modeling. (b) Big box refinement of neutron PDF data. (c) The proportion of interstitial oxygen in different coordination environments from the statistics of the RMC-fitted model. (d–f) Coordination environments of interstitial oxygen: coordinated with (d) one Ge, (e) two Ge, (f) one Ga and one Ge, and (g) one Ga and two Ge.

This behavior arises from the differences in electronegativity (χ) and ionic radius (r) between Ge (χ = 2.01, r = 0.39 Å) and Ga (χ = 1.81, r = 0.47 Å). The higher electronegativity and smaller ionic radius of Ge favor stronger covalent bonding with interstitial oxygen, inducing a confinement effect that traps interstitial oxygens within the 4-membered rings of the zeolite-like feldspar framework. This similar preferential coordination of interstitial oxygen with Ge is also observed in the B-site mixed melilite structure.⁵⁵ In other words, while local collective distortions facilitate the accommodation of interstitial oxygens, they simultaneously impose geometric constraints that limit long-range ion transport. These findings provide direct evidence of the structural variations induced by the incorporation of interstitial oxygen and its role in modulating local ionic mobility: the mobility of the interstitial oxygen ions is stabilized and facilitated by local collective distortions within the GaO₄ and GeO₄ tetrahedra, but is trapped by coordination with Ge atoms within 4-membered rings. This disrupts the coordination equilibrium quasi-free state of interstitial oxygen determined in the average structure, highlighting the crucial role of local structural variations in correlation with the experimentally observed low oxide ion mobility in Sr_0.85La_0.15Ga₂Ge₂O_8.075.

3.6 Interstitial oxygen migration mechanism

To explore the underlying mechanism of interstitial oxygen transport in Sr_0.85La_0.15Ga₂Ge₂O_8.075, molecular dynamics (MD) simulations based on the interatomic potential method were employed. The mean square displacement (MSD) values for cations and oxide ions in the parent SrGa₂Ge₂O₈ (Fig. S8†) show only lattice vibrations within the temperature range of 873–1473 K, with no evidence of long-range ionic migration. In contrast, MSD analysis for Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 7a and b) reveals remarkable long-range migration of interstitial oxygens, accompanied by localized but non-negligible motions of Ge cations, while all other cations exhibit only lattice vibrations. These small but non-negligible motions of Ge atoms (Fig. 7a) indicate that effective interstitial ion transport likely requires local motions of Ge atoms, coupled to the dynamics of the surrounding Ge/GaO₄ frameworks. Notably, the interstitial oxygens exhibit long-range migration with an activation energy of 0.82 eV and oxide ion diffusion coefficients ranging from 5.8 × 10⁻⁹ to 3.9 × 10⁻⁸ cm² s⁻¹ over the temperature range of 873–1473 K (Fig. 7c), highlighting the significant enhancement in oxide ion mobility introduced by the La doping.

Fig. 7 Interstitial oxygen migration in Sr_0.85La_0.15Ga₂Ge₂O₈. (a) Mean square displacement (MSD) values of Sr, La, Ga, Ge, and O atoms. (b) MSD values for oxygen atoms at different lattice sites. (c) Arrhenius plot of oxide ion diffusion coefficients. (d) Scattering plots of La and Sr cations. (e) Scattering plots viewed along the [010] direction. (f and h) Scattering plots of the Ga2Ge1 layer and (g) the Ga1Ge2 layer, extracted from (e). (i) Oxide ion migration pathways in the Ga1Ge2 layer and (j) Ga2Ge1 layer. Deep blue and light blue tetrahedra represent Ga1O₄ and Ga2O₄ tetrahedra, respectively; deep green and light green tetrahedra correspond to Ge2O₄ and Ge1O₄ tetrahedra. Oxygen atoms are color-coded as follows: red/yellow (O3/O4), purple (O1/O2), brown (O5/O6), pink (O7/O8), and orange (O9).

The scatter plots of cationic and oxide ion positions for Sr_0.85La_0.15Ga₂Ge₂O_8.075 (Fig. 7d and e) reveal that all oxygen atoms are involved in oxygen exchange without any cationic exchange, consistent with the above MSD results. To analyse the oxide ion migration pathways in more detail, the scatter plots for the Ga1Ge2 layer and Ga2Ge1 layer were extracted and are shown in Fig. 7g and h, respectively. In the Ga1Ge2 layer (Fig. 7i), it is observed that the interstitial oxygen is initially stabilized within the 4-membered ring formed by two GaO₄ and two GeO₄ tetrahedra (Fig. 7i₁). During the simulation, the interstitial oxygen transitions between coordination states associated with its migration. The interstitial oxygen then coordinates with two GaO₄ and one GeO₄ tetrahedra (Fig. 7i₂), followed by coordination with one GaO₄ and one GeO₄ tetrahedra (Fig. 7i_3–4). Eventually, the interstitial oxygen coordinates with a single GaO₄ tetrahedron within a larger 8-membered ring (Fig. 7i₅), realizing the long-range migration of interstitial oxide ion from the 4-membered ring to the 8-membered ring in the zeolite-like structure. This migration is driven by local structural distortions in the GaO₄ and GeO₄ tetrahedra, which create dynamic pathways for oxide ion migration. The local structural distortions in the GaO₄ and GeO₄ tetrahedra disrupt the local lattice symmetry and facilitate oxide ion mobility. The same migration mechanism is also observed in the Ga2Ge1 layer (Fig. 7j and j_1–6) in Sr_0.85La_0.15Ga₂Ge₂O_8.075. These MD simulation results are consistent with the local structural variations derived from RMC modeling.

The radial distribution functions (RDFs) of Ga/Ge–O interactions (Fig. S9–S16†), derived from MD simulations, provide valuable insight into the local structural disorder in Sr_0.85La_0.15Ga₂Ge₂O_8.075. In the RDFs of Ga2–O interactions over the simulation time range of 0–120 ps, it is observed that the originally coordinated oxygen atoms (O1, O3, O5, and O8, Fig. S9†) surrounding Ga2, are gradually joined by additional oxygen atoms (O2, O6, and O9, Fig. S9†). This behavior is also evident in other Ga–O and Ge–O (Fig. S10–S16†) interactions, demonstrating that all oxygen sites participate in the long-range oxide ion migration.

Zeolite and zeolite-like feldspar-type frameworks can be transformed into oxide-ion conductors through targeted local structural modifications. The targeted aliovalent La³⁺ → Sr²⁺ substitution employed here is one example of such a strategy, readily extendable to other open-framework systems, including aluminogermanates, aluminophosphates, and gallium phosphates. In addition, the selection of B-site framework cations with larger ionic radii and variable coordination numbers (e.g., Ga or Ge instead of Al) can expand tunnel dimensions and promote coordination flexibility, both of which are critical for enhancing ionic mobility. In contrast, Si-based frameworks may be less suitable due to the rigid tetrahedral coordination preference of Si, which limits both structural flexibility and ionic mobility.^56–58 We further propose that zeolite-like frameworks incorporating compositionally pure B-site cations capable of variable coordination represent a promising and underexplored class of oxide-ion conductors. Taken together, these insights establish that oxide-ion transport in zeolite-like structures is governed by the interplay between framework openness, B-site coordination flexibility, and aliovalent doping, providing a promising space for discovering related materials across the broader zeolite family.

The findings highlight the exceptional potential of the zeolite-like feldspar structure for developing interstitial oxide ion conductors for SOFCs, extending its utility beyond traditional applications in building materials and dielectric ceramics.⁵⁹ Although the interstitial oxygen in zeolite-like case Sr_0.85La_0.15Ga₂Ge₂O_8.075 exhibits limited mobility due to higher migration energies induced by the correlated disorder of interstitial oxygen coordinationed with Ge, its favorable ability to accommodate of interstitial oxygen defects suggests promising applications in catalytic processes where ion mobility is less critical. Moreover, zeolite-like structures featuring interstitial oxygen defects offer additional advantages, including increased active sites, tunable electronic structures to enhance catalytic activity, optimized redox cycles, and improved reaction selectivity. These characteristics render zeolite-like structures attractive for applications in electrocatalysis, water splitting, CO₂ reduction, methane oxidation, and exhaust treatment. Furthermore, in many cases, the information derived from the averaged structure often masks critical local structural details, making it challenging to correlate structural features with the ionic transport behaviors in materials. Therefore, employing research methods that focus on local structural investigations has become increasingly essential. In this study, NMR, EXAFS, and neutron PDF based on RMC modeling offer the advantages of uncovering structural features hidden within the average structure. Notably, we identify the correlated disorder of interstitial oxygen species coordinated with Ge atoms within 4-membered rings, accompanied by local collective distortions of neighboring GaO₄ and GeO₄ tetrahedra, which act as carrier traps during oxide ion migration. These features are distinct from both uncorrelated simple local distortions, referring to isolated displacements of individual atoms or small clusters from their ideal lattice positions, and conventional periodic supercell modulations.^60,61 By bridging the gap between local structures and ionic migration behavior, this study uncovers fundamental mechanisms hidden within the average structure and offers a pathway for designing and optimizing ion-conducting materials through tailored modifications of local structures in zeolite-like frameworks.

4. Conclusions

In summary, the zeolite-like structural family, exemplified by strontium feldspar Sr_1−xLa_xGa₂Ge₂O_8+x/2, has been demonstrated as a new structure type for accommodating and transporting interstitial oxide ions. NPD analysis revealed that interstitial oxygen is located at the center of 4-membered rings, coordinated to two Ga and two Ge atoms in a coordination equilibrium quasi-free state, presumably with high mobility within the average structural framework, which however contrasts with the experimentally observed low mobility. Multiple complementary local structural analysis techniques, including 1D and 2D ⁷¹Ga NMR, EXAFS, and neutron PDF based on RMC modeling, demonstrated that the stabilization of interstitial oxygen is facilitated by local collective distortions within the GaO₄ and GeO₄ tetrahedra, but its migration is hindered by the correlated disorder of interstitial oxygen coordinated with Ge atoms within 4-membered rings. The resulting material, Sr_0.85La_0.15Ga₂Ge₂O_8.075, demonstrates excellent thermo-mechanical stability and a suitable TEC value well-matched to that of widely used cobalt-containing cathodes. Although its moderate oxide-ion conductivity and non-negligible electronic conductivity constrain its applicability as a SOFC electrolyte, the combination of a zeolite-like framework capable of accommodating oxygen defects and high thermal stability and compatibility indicates that the zeolite-like family holds strong potential not only for the development of advanced oxide-ion conductors, if the oxide ion mobility can be enhanced, but also for the catalytic applications associated with breaking and reforming of metal–oxygen bonds assisted by and correlated with oxide ion mobility. This finding not only demonstrates the zeolite-like feldspar structure as a new family of oxide ion conductors but also bridges the gap between local structural variations and ionic migration behavior, providing a solid foundation for designing and optimizing ion-conducting materials through local structural modifications in zeolite-like frameworks.

Data availability

The authors confirm that the data supporting the findings of this study are available within the article and its ESI.†

Author contributions

X. K. and X. L. supervised the experimental design and led the overall study. X. W. performed the majority of the experiments as well as data analysis. A. R. assisted with solid-state NMR measurements and contributed to the discussion and interpretation of NMR data. C. L. helped with the neutron PDF data collection. S. D. and L. H. help with the NPD data collection. C. L., H. X., K. L., Q. L., J. S., X. X., participated in the careful discussion on the local structure of zeolite-like structures. The corresponding authors of this manuscript are X. K. and X. L. All authors discussed the results and commented on the manuscript.

Conflicts of interest

The authors declare no competing interests.

Acknowledgements

This research was supported by the Guangxi Natural Science foundation (No. 2025GXNSFBA069588), the Guangxi Key Development Program (No. GuikeAB25069467), the National Natural Science Foundation of China (No. 22205017, 22090043, and 22090042), the National Key R&D Program of China (2020YFA0406202), and the Guilin University of Technology Research Startup Project (No. RD2400002912). The authors are grateful for the Guangxi BaGui Scholars Special Funding. Financial support from the IR-RMN-THC Fr3050 CNRS for conducting the research is gratefully acknowledged.

References

1. B. C. H. Steele and A. Heinzel, Materials for fuel-cell technologies, Nature, 2001, 414, 345–352.
2. T. Famprikis, P. Canepa, J. A. Dawson, M. S. Islam and C. Masquelier, Fundamentals of inorganic solid-state electrolytes for batteries, Nat. Mater., 2019, 18, 1278–1291.
3. X. Yang, A. J. Fernandez-Carrion and X. Kuang, Oxide Ion-Conducting Materials Containing Tetrahedral Moieties: Structures and Conduction Mechanisms, Chem. Rev., 2023, 123, 9356–9396.
4. D. Schweke, Y. Mordehovitz, M. Halabi, L. Shelly and S. Hayun, Defect Chemistry of Oxides for Energy Applications, Adv. Mater., 2018, 30, 1706300.
5. Z. Zou, Y. Li, Z. Lu, D. Wang, Y. Cui, B. Guo, Y. Li, X. Liang, J. Feng, H. Li, C. W. Nan, M. Armand, L. Chen, K. Xu and S. Shi, Mobile Ions in Composite Solids, Chem. Rev., 2020, 120, 4169–4221.
6. C. Zhao, Y. Li, W. Zhang, Y. Zheng, X. Lou, B. Yu, J. Chen, Y. Chen, M. Liu and J. Wang, Heterointerface engineering for enhancing the electrochemical performance of solid oxide cells, Energy Environ. Sci., 2020, 13, 53–85.
7. J. H. Myung, D. Neagu, D. N. Miller and J. T. Irvine, Switching on electrocatalytic activity in solid oxide cells, Nature, 2016, 537, 528–531.
8. E. D. Wachsman and K. T. Lee, Lowering the Temperature of Solid Oxide Fuel Cells, Science, 2011, 334, 935–939.
9. Y. Zhang, B. Chen, D. Guan, M. Xu, R. Ran, M. Ni, W. Zhou, R. O'Hayre and Z. Shao, Thermal-expansion offset for high-performance fuel cell cathodes, Nature, 2021, 591, 246–251.
10. T. Kudo and H. Obayashi, Oxygen Ion Conduction of the Fluorite-Type Ce_1-xLn_xO_2-x/2 (Ln = Lanthanoid Element), J. Electrochem. Soc., 1975, 122, 142–147.
11. J. Shamblin, M. Feygenson, J. Neuefeind, C. L. Tracy, F. Zhang, S. Finkeldei, D. Bosbach, H. Zhou, R. C. Ewing and M. Lang, Probing disorder in isometric pyrochlore and related complex oxides, Nat. Mater., 2016, 15, 507–511.
12. M. Yashima, T. Tsujiguchi, Y. Sakuda, Y. Yasui, Y. Zhou, K. Fujii, S. Torii, T. Kamiyama and S. J. Skinner, High oxide-ion conductivity through the interstitial oxygen site in Ba₇Nb₄MoO₂₀-based hexagonal perovskite related oxides, Nat. Commun., 2021, 12, 556.
13. Z. J. Luo, X. H. Li, X. G. Wang, S. H. Deng, L. H. He, K. Lin, Q. Li, X. R. Xing and X. J. Kuang, The Langasite Family for the Development of Oxygen-Vacancy-Mediated Oxide Ion Conductors, Chem. Mater., 2024, 36, 2835–2845.
14. P. Lacorre, F. Goutenoire, O. Bohnke, R. Retoux and Y. Laligant, Designing fast oxide-ion conductors based on La₂Mo₂O₉, Nature, 2000, 404, 856–858.
15. X. Y. Yang, A. J. Fernandez-Carrion, J. H. Wang, F. Porcher, F. Fayon, M. Allix and X. J. Kuang, Cooperative mechanisms of oxygen vacancy stabilization and migration in the isolated tetrahedral anion Scheelite structure, Nat. Commun., 2018, 9, 4484.
16. M. L. Tate, D. A. Blom, M. Avdeev, H. E. A. Brand, G. J. McIntyre, T. Vogt and I. R. Evans, New Apatite-Type Oxide Ion Conductor, Bi₂La₈[(GeO₄)₆]O₃: Structure, Properties, and Direct Imaging of Low-Level Interstitial Oxygen Atoms Using Aberration-Corrected Scanning Transmission Electron Microscopy, Adv. Funct. Mater., 2017, 27, 1605625.
17. X. J. Kuang, M. A. Green, H. Niu, P. Zajdel, C. Dickinson, J. B. Claridge, L. Jantsky and M. J. Rosseinsky, Interstitial oxide ion conductivity in the layered tetrahedral network melilite structure, Nat. Mater., 2008, 7, 498–504.
18. X. H. Li, L. Q. Yang, Z. Y. Zhu, X. G. Wang, P. H. Chen, S. C. Huang, X. Y. Wei, G. H. Cai, P. Manuel, S. H. Yang, J. H. Lin, X. J. Kuang and J. L. Sun, Borates as a new direction in the design of oxide ion conductors, Sci. China Mater., 2022, 62, 2737–2745.
19. M. Diaz-Lopez, J. F. Shin, M. Li, M. S. Dyer, M. J. Pitcher, J. B. Claridge, F. Blanc and M. J. Rosseinsky, Interstitial Oxide Ion Conductivity in the Langasite Structure: Carrier Trapping by Formation of (Ga,Ge)₂O₈ Units in La₃Ga_5-xGe_1+xO_14+x/2 (0 < x ≤ 1.5), Chem. Mater., 2019, 31, 5742–5758.
20. L. Corti, D. Iuga, J. B. Claridge, M. J. Rosseinsky and F. Blanc, Disorder and Oxide Ion Diffusion Mechanism in La_1.54Sr_0.46Ga₃O_7.27 Melilite from Nuclear Magnetic Resonance, J. Am. Chem. Soc., 2023, 145, 21817–21831.
21. J. M. S. Skakle and R. Herd, Crystal chemistry of (RE,A)₂M₃O₇ compounds (RE=Y, lanthanide; A=Ba, Sr, Ca; M=Al, Ga), Powder Diffr., 1999, 14, 195–202.
22. C. Tealdi, P. Mustarelli and M. S. Islam, Layered LaSrGa₃O₇-Based Oxide-Ion Conductors: Cooperative Transport Mechanisms and Flexible Structures, Adv. Funct. Mater., 2010, 20, 3874–3880.
23. D. V. Roshchupkin, D. V. Irzhak, O. A. Plotitsyna, R. R. Fakhrtdinov, O. A. Buzanov and A. P. Sergeev, Piezoelectric crystal Ca₃TaGa₃Si₂O₁₄: synthesis, structural perfection, piezoelectric and acoustic properties, J. Surf. Invest., 2012, 6, 490–493.
24. N. Lee, Y. J. Choi and S.-W. Cheong, Magnetic control of ferroelectric polarization in a self-formed single magnetoelectric domain of multiferroic Ba₃NbFe₃Si₂O₁₄, Appl. Phys. Lett., 2014, 104, 072904.
25. S. Zhang, Y. Zheng, H. Kong, J. Xin, E. Frantz and T. R. Shrout, Characterization of high temperature piezoelectric crystals with an ordered langasite structure, J. Appl. Phys., 2009, 105, 114107.
26. Y. Li and J. H. Yu, Emerging applications of zeolites in catalysis, separation and host–guest assembly, Nat. Rev. Mater., 2021, 6, 1156–1174.
27. Y. Li and J. H. Yu, New Stories of Zeolite Structures: Their Descriptions, Determinations, Predictions, and Evaluations, Chem. Rev., 2014, 114, 7268–7316.
28. J. C. Barbosa, R. Gonçalves, C. M. Costa, V. de Zea Bermudez, A. Fidalgo-Marijuan, Q. Zhang and S. Lanceros-Méndez, Metal–organic frameworks and zeolite materials as active fillers for lithium-ion battery solid polymer electrolytes, Mater. Adv., 2021, 2, 3790–3805.
29. Y. Li, X. Wang, J. Liang, K. Wu, L. Xu and J. Wang, Design of A High Performance Zeolite/Polyimide Composite Separator for Lithium-Ion Batteries, Polymers, 2020, 12, 764.
30. D. Chen, X. Wang, J. Liang, Z. Zhang and W. Chen, A Novel Electrospinning Polyacrylonitrile Separator with Dip-Coating of Zeolite and Phenoxy Resin for Li-ion Batteries, Membranes, 2021, 11, 267.
31. X. Chi, M. Li, J. Di, P. Bai, L. Song, X. Wang, F. Li, S. Liang, J. Xu and J. Yu, A highly stable and flexible zeolite electrolyte solid-state Li-air battery, Nature, 2021, 592, 551–557.
32. M. Calleri and G. Gazzoni, The Structures of (Sr,Ba) [(AI,Ga)₂(Si,Ge)₂O₈]. III. The Crystal Structures of the ParacelsianLike Modifications of Synthetic SrGa₂Ge2O₈ and BaGa₂Ge₂O₈, Acta Cryst. v, 1976, 32, 1196–1205.
33. J. P. Xu, Y. G. Xia, Z. D. d. Li, H. C. Chen, X. L. Wang, Z. Z. Sun and W. Yin, Multi-physics instrument: total scattering neutron time-of-flight diffractometer at China Spallation Neutron Source, Nucl. Instrum. Methods Phys. Res., Sect. A, 2021, 1013, 165642.
34. C. L. Farrow, P. Juhas, J. W. Liu, D. Bryndin, E. S. Božin, J. Bloch, T. Proffen and S. J. L. Billinge, PDFfit2 and PDFgui: computer programs for studying nanostructure in crystals, J. Phys.: Condens. Matter, 2007, 19, 335219.
35. M. G. Tucker, D. A. Keen, M. T. Dove, A. L. Goodwin and Q. Hui, RMCProfile: reverse Monte Carlo for polycrystalline materials, J. Phys.: Condens. Matter, 2007, 19, 335218.
36. A. Bruker, TOPAS V3: General profile and structure analysis software for powder diffraction data, Bruker AXS, Karlsruhe, Germany, 2005.
37. M. Dušek, V. Petříček, L. Palatinus, J. Rohlíček, K. Fejfarová, V. Eigner, M. Kučeráková, M. Poupon, M. Henriques, J. Plášil and E. Skořepová, Jana 2006 Cookbook, Czech Republic, Institute of Physics ASCR, Na Slovance 2, Prague, 2019.
38. D. Johnson, Zview for Windows, Impedance/Gain Phase Graphing and Analysis Software, Scribner Associates, North Carolina, 2001.
39. W. Wan, J. Sun, J. Su, S. Hovmoller and X. Zou, Three-dimensional rotation electron diffraction: software RED for automated data collection and data processing, J. Appl. Crystallogr., 2013, 46, 1863–1873.
40. D. Massiot, F. Fayon, M. Capron, I. King, S. Le Calvé, B. Alonso, J. O. Durand, B. Bujoli, Z. Gan and G. Hoatson, Modelling one- and two-dimensional solid-state NMR spectra, Magn. Reson. Chem., 2001, 40, 70–76.
41. J. D. Gale and A. L. Rohl, The General Utility Lattice Program (GULP), Mol. Simul., 2003, 29, 291–341.
42. C. Catlow, Mott-Littleton calculations in solid-state chemistry and physics, J. Chem. Soc., Faraday Trans. 2, 1989,(85), 335–340.
43. G. Sastre and J. D. Gale, Derivation of an Interatomic Potential for Germaniumand Silicon-Containing Zeolites and Its Application to the Study of the Structures of Octadecasil, ASU-7, and ASU-9 Materials, Chem. Mater., 2003, 15, 1788–1796.
44. S. Nosé, A unified formulation of the constant temperature molecular dynamics methods, J. Chem. Phys., 1984, 81, 511–519.
45. W. Humphrey, A. Dalke and K. Schulten, VMD: Visual Molecular Dynamics, J. Mol. Graphics, 1996, 14, 33–38.
46. T. Rog, K. Murzyn, K. K. Hinsen and G. R. Kneller, Software News and Updates nMoldyn: A Program Package for a Neutron Scattering Oriented Analysis of Molecular Dynamics Simulations, J. Comput. Chem., 2002, 24, 658–666.
47. Y. Song, Y. Chen, M. Xu, W. Wang, Y. Zhang, G. Yang, R. Ran, W. Zhou and Z. P. Shao, A Cobalt Free Multi Phase Nanocomposite as Near Ideal Cathode of Intermediate Temperature Solid Oxide Fuel Cells Developed by Smart Self Assembly, Adv. Mater., 2020, 32, 1906979.
48. X. Ding, C. Cui and L. Guo, Thermal expansion and electrochemical performance of La_0.7Sr_0.3CuO_3-δ-Sm_0.2Ce_0.8O_2-δ composite cathode for IT-SOFCs, J. Alloys Compd., 2009, 481, 845–850.
49. W. Zhou, Z. Shao, R. Ran and R. Cai, Novel SrSc_0.2Co_0.8O_3-δ as a cathode material for low temperature solid-oxide fuel cell, Electrochem. Commun., 2008, 10, 1647–1651.
50. M. W. Phillips, H. Kroll and H. Pentinghaus, The Structures of Synthetic Paracelsian Analogs, SrGa₂Si₂O₈ and SrGa₂Ge₂O₈, Am. Mineral., 1975, 60, 659–666.
51. H. K. C. Timken and E. Oldfield, Solid-State Gallium-69 and Gallium-71 Nuclear Magnetic Resonance Spectroscopic Studies of Gallium Analogue Zeolites and Related Systems, J. Am. Chem. Soc., 1987, 109, 7669–7673.
52. D. Massiot, I. Ian Farnan, N. Gautier, D. Trumeau, A. Trokiner and J. P. Coutures, ⁷¹Ga and ⁶⁹Ga nuclear magnetic resonance study of β-Ga₂O₃: resolution of four- and six-fold coordinated Ga sites in static conditions, Solid State Nucl. Magn. Reson., 1995, 4, 241–248.
53. J. Fan, V. Sarou-Kanian, X. Yang, M. Diaz-Lopez, F. Fayon, X. Kuang, M. J. Pitcher and M. Allix, La₂Ga₃O_7.5: A Metastable Ternary Melilite with a Super-Excess of Interstitial Oxide Ions Synthesized by Direct Crystallization of the Melt, Chem. Mater., 2020, 32, 9016–9025.
54. J. G. Xu, J. H. Wang, A. Rakhmatullin, S. Ory, A. J. Fernández-Carrión, H. B. Yi, X. J. Kuang and M. Allix, Interstitial Oxide Ion Migration Mechanism in Aluminate Melilite La_1+xCa_1–xAl₃O_7+0.5x Ceramics Synthesized by Glass Crystallization, ACS Appl. Energy Mater., 2019, 2, 2878–2888.
55. X. Li, X. Wei, X. Wang, C. Lou, W. Zhang, J. Xu, Z. Luo, M. Tang, S. Deng, L. He, X. Xing, J. Sun and X. Kuang, B-site mixed cationic tetrahedral layer confined the concentration and mobility of interstitial oxygen in mellite family, J. Mater. Chem. A, 2023, 11, 5615–5626.
56. I. R. Evans, J. S. O. Evans, H. G. Davies, A. R. Haworth and M. L. Tate, On Sr_1–xNa_xSiO_3–0.5x New Superior Fast Ion Conductors, Chem. Mater., 2014, 26, 5187–5189.
57. P. H. Chien, Y. Jee, C. Huang, R. Dervisoglu, I. Hung, Z. Gan, K. Huang and Y. Y. Hu, On the origin of high ionic conductivity in Na-doped SrSiO(3), Chem. Sci., 2016, 7, 3667–3675.
58. R. D. Bayliss, S. N. Cook, S. Fearn, J. A. Kilner, C. Greaves and S. J. Skinner, On the oxide ion conductivity of potassium doped strontium silicates, Energy Environ. Sci., 2014, 7, 2999–3005.
59. V. Fuertes, J. J. Reinosa, J. F. Fernández and E. Enríquez, Engineered feldspar-based ceramics: a review of their potential in ceramic industry, J. Eur. Ceram. Soc., 2022, 42, 307–326.
60. A. J. Garcia-Adeva, D. R. Conradson, P. Villella and S. D. Conradson, Signatures of Collective Local and Nanoscale Distortions in Diffraction Experiments, J. Phys. Chem. B, 2003, 107, 6704–6716.
61. T. R. Welberry and B. D. Butler, Diffuse X-ray Scattering from Disordered Crystals, Chem. Rev., 1995, 95, 2369–2403.

---

Footnote

† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5sc02898a

---