Koji
Okada
a,
Susumu
Fujii
bc,
Cédric
Tassel
*a,
Shenghan
Gao
a,
Hiroki
Ubukata
a,
Wenli
Pan
d,
Kentaro
Yamamoto
de,
Yoshiharu
Uchimoto
d,
Akihide
Kuwabara
*b and
Hiroshi
Kageyama
*a
aDepartment of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan
bNanostructures Research Laboratory, Japan Fine Ceramics Center, Nagoya 456-8587, Japan
cDivision of Materials and Manufacturing Science, Graduate School of Engineering, Osaka University, Osaka 565-0871, Japan
dGraduate School of Human and Environment Studies, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
eFaculty of Engineering, Nara Women's University, Kitauoyanishimachi, Nara 630-8506, Japan
First published on 26th May 2023
Unlike perovskite oxides, antiperovskites M3HCh and M3FCh (M = Li, Na; Ch = S, Se, Te) mostly retain their ideal cubic structure over a wide range of compositions owing to anionic size flexibility and low-energy phonon modes that promote their ionic conductivity. In this study, we show the synthesis of potassium-based antiperovskites K3HTe and K3FTe and discuss the structural features in comparison with lithium and sodium analogues. It is shown experimentally and theoretically that both compounds maintain a cubic symmetry and can be prepared at ambient pressure, in contrast to most of the reported M3HCh and M3FCh which require high pressure synthesis. A systematic comparison of a series of cubic M3HTe and M3FTe (M = Li, Na, K) revealed that telluride anions contract in the order of K, Na, Li, with a pronounced contraction in the Li system. This result can be understood in terms of the difference in charge density of alkali metal ions as well as the size flexibility of Ch anions, contributing to the stability of the cubic symmetry.
Subsequently, antiperovskites M3FCh (M = Li, Na; Ch = S, Se, Te) were synthesized mainly using the high pressure technique.11 Despite the use of the hard F anion instead of the soft H anion, M3FCh also adopts the cubic structure with the exception of Na3FS, and thus comparable ionic conduction is obtained due to the associated low energy FM6 octahedral rotation modes. These results, together with first-principles calculations, imply that highly polarizable Ch anions play a major role in stabilizing the cubic symmetry. Furthermore, we found a correlation between the soft phonon modes, the tolerance factor,12 and the activation energy for the mobile cations.
This work aims to expand the available compositional space of M3HCh and M3FCh by including potassium instead of Li and Na. We have successfully synthesized K3HTe and K3FTe under both high pressure and ambient pressure conditions. The preparation of a series of alkali metal-rich antiperovskites allows us to compare them systematically together with first-principles calculations. It is found that in M3BCh (M = Li, Na, K; B = H, F; Ch = S, Se, Te), not only the size-flexibility of the Ch-site anions but also the charge density of the M-site cations has an impact on their structures.
High-resolution synchrotron X-ray diffraction (SXRD) measurements were carried out at room temperature with a wavelength of λ = 0.42 Å at the BL02B2 beamline in SPring-8 (Japan), equipped with MYTHEN solid-state detectors. The K3HTe and K3FTe powder samples were sealed in borosilicate capillaries. Le Bail and Rietveld analyses of SXRD data were carried out using JANA2006 software,17 while the crystal structures were visualized using VESTA.18
Ionic conductivities of cold-pressed pellets of K3HTe and I-substituted K3HTe were measured by electrochemical impedance spectroscopy (EIS) with a constant voltage of 10 mV in the frequency range from 1 MHz to 0.1 Hz using an ECS ModuLab potentiostat/galvanostat. Around 0.2 g of the powder specimen was placed between two stainless-steel rods as an ion-blocking electrode in a custom-made Swagelok cell and pressed into a 10 mm diameter pellet using a hydraulic press at 360 MPa (thickness: 1–1.5 mm). The temperature was varied from 25 °C to 100 °C. Data at each temperature were collected after being held for 115 min to ensure temperature stabilization.
Ab initio calculations on K3HTe, K3FTe, KH, KF and K2Te were performed under 0 and 5 GPa using the plane-wave basis projector augmented wave (PAW) method19 implemented in the Vienna Ab initio Simulation Package (VASP).20,21 The generalized gradient approximation in the form of PBEsol22 was employed. The valence configurations in the PAW potentials used in this study are as follows: [3s2 3p6 4s1] for K, [1s1] for H, [2s2 2p5] for F, and [5s2 5p4] for Te. The remaining electrons were treated as frozen core electrons. The plane-wave energy cutoff was set to 650 eV. The first Brillouin zone was sampled using Γ-centered k-point grids with a spacing of 0.2 Å−1. The total energy convergence was set to 1.0 × 10−6 eV per atom, and atomic positions and lattice constants were relaxed until the residual forces of atoms were less than 1.0 × 10−3 eV Å−1.
Phonon band structures and vibrational free energies under harmonic approximation were calculated using the finite displacement method23 implemented in Phonopy.24 The displacement of atoms from their stationary positions to obtain second-order force constants was set to 0.01 Å. A 2 × 2 × 2 supercell with respect to the conventional unit cell was used to avoid the self-interaction between the displaced atoms. Non-analytical term correction25 was adopted using the Born effective charge tensors and dielectric constants calculated using density functional perturbation theory implemented in VASP.20,21
The nudged elastic band (NEB) method26 was used to calculate the energy barriers of K+ migration via vacancies and interstitial defects (dumbbell) in K3HTe and K3FTe. Charged supercells with the dimensions of 3 × 3 × 3 were used. The first Brillouin zone was sampled with one k-point (0.25, 0.25, 0.25) due to the large size of supercells. For NEB calculations, the total energy and force convergence were set to 1.0 × 10−6 eV per cell and 2.0 × 10−2 eV Å−1, respectively. Note that all the computational conditions employed here are consistent with our previous studies.10,11
Rietveld refinement of the SXRD data for K3HTe assuming the cubic antiperovskite structure (Pm−m) readily converged to GOF = 3.94, Rp = 5.44% and Rwp = 7.56% (Fig. 2a), with reasonable isotropic atomic displacement parameter, Uiso, values of 0.0385(2) Å2 and 0.0239(1) Å2 for K and Te, respectively. Refining occupancy factors, g, at K resulted in 0.988(1). Note that the Uiso of H and g of H and Te are not refined. Likewise, the Rietveld analysis of K3FTe converged successfully, yielding GOF = 5.52, Rp = 5.33% and Rwp = 7.31%, Uiso of 0.0235(2) Å2 for K, 0.0155(9) Å2 for F and 0.0165(1) Å2 for Te (Fig. 2b), and the occupancy of 0.990(1) for K and 0.960(5) for F. Thus, unless otherwise stated, we hereafter use the stoichiometric compositions of K3HTe and K3FTe.
Most of the alkali metal rich antiperovskites M3HCh and M3FCh (M = Li, Na) reported so far have been synthesized at high pressure,10,11 with the exception of Na3HTe and Na3FTe (unpublished data). First-principles calculations showed that the ideal cubic structure is dynamically stable at high pressure for most Na- and Li-based compounds, but at ambient pressure, many of them showed imaginary phonon modes appearing at the M and R points relevant to octahedral rotation. Here, we calculated the phonon band structures of K3HTe and K3FTe using the cubic structure at 0 and 5 GPa and found no imaginary phonon modes even under ambient pressure (Fig. 3a and b). The atom projected phonon density of states (DOS) for K3HTe indicates that the light H ion is responsible for high frequencies, whereas the heavy K and Te ions are for low frequencies. In K3FTe, this tendency is less pronounced. In addition, both compounds are lower in energy at ambient pressure than the starting materials K2Te and KH/KF, suggesting that they are thermodynamically stable phases (Fig. 3c and d). These results led us to synthesize K3HTe and K3FTe under ambient pressure, and we indeed obtained the XRD patterns with identical lattice constants for both cases (Fig. S2 and S3†).
![]() | ||
Fig. 3 Phonon band structures and partial phonon DOS of each ion of the cubic antiperovskites (a) K3HTe and (b) K3FTe at 0 GPa. The absence of imaginary phonon frequencies at the M and R points related to octahedral tilting, as found in, e.g., Na3HS and Na3FS,10,11 means that both compounds are dynamically stable. The Gibbs free energy ΔG (ΔG = ΔUinternal + PΔV − ΔFphonon, ΔFphonon = ΔEZPE − TΔSphonon) of (c) K3HTe relative to that of KH + K2Te and (d) K3FTe relative to that of KF + K2Te at 0 GPa, all calculated at ambient pressure. |
The discovery of K3HTe and K3FTe offers an opportunity to systematically compare the impact of the cations on stabilization in the two sets of cubic antiperovskites, M3HTe and M3FTe (M = Li, Na, K). Fig. 4a shows the lattice parameters of cubic M3HTe and M3FTe as a function of Shannon ionic radii of the alkali metal cations, rM.27 In both systems, the lattice expands linearly with increasing rM, but the slope is obviously different, with a larger value for M3HTe, which suggests that this phenomenon is related to the size flexibility of hydride anions (vs. F−).10,28 To obtain further insights into the structural features, the obtained lattice constants were used to estimate the ionic radius of H− anions, rH, where rH = a/2 − rM (rLi = 0.76 Å, rNa = 1.02 Å, rK = 1.38 Å), as employed previously (Fig. 4b).10 The expanded H− anion in Na3HTe compared to that in K3HTe looks reasonable because the smaller ionic radius of Na+ (vs. K+) is compensated by the larger H− anion. Following this trend, the H− radius of Li3HTe is expected to be larger than that of Na3HTe. However, the opposite result, in which H− contracts, was obtained, which implies that the hydride anion is not solely responsible for stabilizing and maintaining the cubic structure.
![]() | ||
Fig. 4 (a) The lattice parameters of (blue) M3HTe and (red) M3FTe (M = Li, Na, K).10,11 (b) The ionic radius of H− and F−, rH and rF, determined from the lattice constant, where rH = (a − 2rM)/2 and rF = (a − 2rM)/2. The fixed ionic radii of 0.76 Å, 1.02 Å and 1.38 Å were used for Li, Na and K, respectively.27 Error bars are within the circles in both figures. |
We hence took a theoretical approach to address this issue. We estimated the Bader volume29 for each ion using the relaxed structures at ambient pressure by first-principles calculations. In terms of the H− volume, its M dependence with a contraction at Li3HTe appears to be reproduced, as shown in Fig. 5. However, similar dependence is also found in the F− volume both qualitatively and quantitatively, with a smaller F− volume in Li3FTe than in Na3FTe. Instead, a more pronounced M-dependence is seen in the tellurium anion, where the Te anion decreases its size in the order of K, Na, and Li.
![]() | ||
Fig. 5 The Bader volume of the alkali metal cation (green), H−/F− anion (blue), and Te2− anion (yellow) in M3HTe (circles) and M3FTe (triangles). |
The large variation in the size of the Te anion as a function of the alkali metal cations or unexpected contractions of the B-site anions in Li systems is likely related to the surface charge density of the alkali metal which, when increasing from K to Na and Li, leads to stiffer anions. We have noticed similar behaviour in lattice constants or anionic sizes in binary compounds, such as rock-salt MH30–35 and antifluorite M2Te36–39 (Fig. S4†). Thus, the charge density of the M-site cation seems to predominantly modify the size of the Te anion to compensate for the size mismatch and preserve the cubic symmetry. Bearing in mind that Te has a high compressibility, the exceptionally small Te anion in Li3HTe is consistent with the fact that it can only be synthesized under high pressure.
Since M3HCh and M3FCh (M = Li, Na; Ch = S, Se, Te) exhibit excellent lithium or sodium ion conductivity,10,11 we performed nudged elastic band (NEB) calculations to estimate the migration energy barriers of K+ transport in K3HTe and K3FTe. As seen in Fig. 6, their energy barriers are comparable to those of Li-rich and Na-rich antiperovskites, as low as 265 meV and 352 meV, respectively, for the vacancy mechanism, and 76 meV and 102 meV for the interstitial dumbbell mechanism. The low migration energy barriers in K-rich antiperovskites indicate their potential as ionic conductors, once vacancies or interstitial defects are introduced by the aliovalent substitution approach, for example. K3HTe has a lower energy barrier than K3FTe, which is partly due to the softer phonon mode relevant to HK6 octahedral rotation (Fig. S5†). Attempts to measure the K+ ion conductivity of pristine and I-doped K3HTe resulted in high resistances (at least >106 Ω), which precluded proper estimation of bulk conductivities and activation energies (Fig. S6†). The large resistances may be attributed to the extremely high sensitivity of the sample surface to trace amounts of oxygen and water; thus, there is considerable room to adjust synthetic conditions in the future.
![]() | ||
Fig. 6 Migration energy barriers of (circles) M3HTe and (triangles) M3FTe (where M = Li, Na, K). The results of vacancy and interstitial dumbbell mechanisms are shown in red and blue, respectively. The schematics on the right side show the migration pathways of M+ ions around the HM6/FM6 octahedron. The migration energy barriers of Li- and Na-rich antiperovskites are taken from ref. 10 and 11. |
Footnote |
† Electronic supplementary information (ESI) available: Laboratory X-ray diffraction patterns of K3HTe and K3FTe, size-flexibility of soft anions in other compounds, migration barriers in K3HTe and K3FTe, and measured ionic conduction of pristine and iodine substituted K3HTe (PDF). See DOI: https://doi.org/10.1039/d3dt01039b |
This journal is © The Royal Society of Chemistry 2023 |