Rapid Low-Dimensional Li+ Ion Hopping Processes in Synthetic Hectorite-Type Li0.5[Mg2.5Li0.5]Si4O10F2.

Understanding the origins of fast ion transport in solids is important to develop new ionic conductors for batteries and sensors. Nature offers a rich assortment of rather inspiring structures to elucidate these origins. In particular, layer-structured materials are prone to show facile Li+ transport along their inner surfaces. Here, synthetic hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2, being a phyllosilicate, served as a model substance to investigate Li+ translational ion dynamics by both broadband conductivity spectroscopy and diffusion-induced 7Li nuclear magnetic resonance (NMR) spin-lattice relaxation experiments. It turned out that conductivity spectroscopy, electric modulus data, and NMR are indeed able to detect a rapid 2D Li+ exchange process governed by an activation energy as low as 0.35 eV. At room temperature, the bulk conductivity turned out to be in the order of 0.1 mS cm-1. Thus, the silicate represents a promising starting point for further improvements by crystal chemical engineering. To the best of our knowledge, such a high Li+ ionic conductivity has not been observed for any silicate yet.

Introduction

The diffusion of small cations and anions plays an important role in many devices such as sensors[1−3] and batteries.[4−11] Although for some applications, e.g., in the semiconductor industry[12] or in the development of breeding materials[13−15] for fusion reactor blankets, diffusion is unwanted, in other branches, materials with extremely high diffusion coefficients are desired.[16−21] Finding the right material and tailoring its dynamic properties further requires an in-depth understanding of the origins that determine fast ion dynamics.[22−26]

In many cases, layer-structured materials[27−30] are known to offer fast diffusion pathways along the buried interfaces or between larger gaps inside the crystal structure.[31] Commonly, powder samples are synthesized and investigated; however, polycrystalline samples do not allow orientation-dependent measurements of ionic conductivity. Fortunately, anisotropic properties of ionic transport in polycrystalline samples can be probed by nuclear magnetic resonance (NMR) spin–lattice relaxation experiments.[5,30,32−35] The most prominent examples, whose Li+ diffusion properties were studied in this way,[36] include graphite[37,38] or transition metal chalcogenides such as TiS2,[39−43] NbS2,[44,45] and SnS2.[46] Recently, two-dimensional (2D) Li+ diffusion has also been determined in hexagonal LiBH4.[29,47] Fast fluoride, F–, diffusion is observed in MeSnF4 (Me = Pb, Ba) and, as has been shown quite recently, also in layer-structured RbSn2F5.[48] In these materials, spatial constraints guide the ions over long distances. The principle of guided ions has also been found in Li12Si7.[49,50] In binary silicide, the Li+ ions are subjected to a fast one-dimensional diffusion process along the surface of a virtual pipe formed by the stacked Si5 rings.[49] In channel-structured materials with the ions diffusing inside the channels, their diffusion pathways may, however, easily be blocked by foreign, immobile ions.[51] The same effect might, of course, also influence 2D translational ion dynamics but to a lesser degree.

To understand Li+ diffusion in structures offering 2D diffusion pathways, we chose hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2 as a model system, see Figure , and studied the Li+ self-diffusion properties and electrical ionic transport. While the latter is investigated by broadband conductivity spectroscopy,[52] we took advantage of 7Li NMR spin–lattice relaxation measurements[32,53] to shed light on Li+ translational dynamics.

Figure 1

(a) Idealized crystal structure of hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2 (drawn according to the space group C2/m, the real structure is clearly turbostratic) with two distinct crystallographic sites occupied by the Li ions. Those between the Mg-rich layers (see (b)), which are here shown following a homogenous charge distribution, are expected to have access to fast 2D diffusion pathways.

Hectorite is, as montmorillonite, a clay mineral of the smectite group (gr. σμηκτις, soil with the ability to clean, to soak something up, e.g., swelling through water uptake). Such silicates belong to the family of phyllosilicates (gr. φυλλο, leaf). Fluorohectorite is a 2D host material with a rigid 2:1 sandwich-like structure;[54] each lamella is composed of two silicate and one Mg2+ rich layer where the earth alkaline cations are octahedrally coordinated by oxygen and fluorine anions. An idealized structure of synthetic (turbostratic) Li-bearing fluorohectorite [Li0.5]inter[Mg2.5Li0.5]octa[Si4]tetraO10F2 is shown in Figure . Natural hectorite of the composition Na0.3(Mg, Li)3Si4O10(OH)2, including also minor amounts of Ca2+, K+, and Al3+, was named for its occurrence 5 km south of Hector (San Bernardino County, California, USA).[55] It was first described in 1941 and found in a bentonite deposit, altered from clinoptilolite, a natural zeolite, derived from volcanic tuff and ash with a high glass content.[55] Cation exchange strategies allow one to effectively replace the Na ions in hectorite by Li ions; also larger inorganic or organic molecules (pillars) can be introduced, resulting in porous so-called pillared clays.[54] Exemplarily, the Khan group successfully used hectorite-based materials as passive and active filler materials to prepare nanocomposite polymer-based (gel) electrolytes as well as to develop composite LiCoO2-based cathode materials.[56−59] Here, our hypothesis is that the interlayer gap in such host structures, particularly that of hectorite, offers indeed fast diffusion pathways for small charge carriers such as Li+ and Na+ ions.

Verifying the hypothesis of 2D transport is, however, challenging if only powder samples are available at hand. Fortunately, 7Li NMR spin–lattice relaxation measurements,[49,53] also successfully used to characterize electrolytes for batteries,[60,61] represent a unique tool to study such anisotropic properties even for powdered samples, since the spectral density functions J governing the NMR spin–lattice relaxation rates 1/T1 (∝ J) possess specific features for 1, 2, and 3D ion transport.[30,62] As intimated above, low-dimensional ionic transport has, so far, been probed only for a limited number of materials by NMR relaxation techniques.[36] The present study contributes to this research field and is aimed at answering the question whether 2D silicate structures offer an assortment of materials with enhanced ion diffusion properties.

Experimental Section

Preparation and Characterization of Li0.5[Mg2.5Li0.5]Si4O10F2

Lithium fluorohectorite with the nominal composition of Li0.5[Mg2.5Li0.5]<Si4 > O10F2 was prepared by a two-step ion exchange procedure starting from melt synthesized Na-fluorohectorite (Na0.5[Mg2.5Li0.5]<Si4 > O10F2). The Na-fluorohectorite was synthesized according to a preparation route published earlier.[54,63] To ensure the complete exchange of Na+ by Li+, Na-fluorohectorite was first treated with n-butylammonium chloride (C4H12ClN (C4), 2 M) using a 100-fold excess of its cation exchange capacity (CEC) to obtain the so-called C4-fluorohectorite. Typically, 500 mg of Na-fluorohectorite was exchanged overnight five times with 40 mL of C4 solution at 80 °C. The resulting product was washed five times with an ethanol–water mixture (1:1) and once more with pure ethanol. Atomic absorption spectroscopy (AAS), see below, was used to check the completeness of the ion exchange reaction. After this, the dried C4-fluorohectorite was treated with a 100-fold excess of the CEC of LiOH (2 M, 4 times 40 mL, 6 h, room temperature). While removing excessive LiOH by washing 6 times with 40 mL of deionized water, the ionic strength decreased to a value lower than 0.02 M; furthermore, complete delamination by repulsive osmotic swelling was observed.[63−66] Finally, the gel of delaminated Li-fluorohectorite was freeze-dried followed by vacuum drying at 160 °C.

AAS

To verify complete Na+-exchange, AAS was used. For this purpose, approximately 20 mg of the samples were weighed into clean Teflon-flasks of 15 mL volume. After the addition of 1.5 mL of 30 wt % HCl (Merck), 0.5 mL of 85 wt % H3PO4 (Merck), 0.5 mL of 65% HNO3 (Merck), and 1 mL of 48 %HBF4 (Merck), the samples were digested in a 1200 Mega microwave digestion apparatus (MLS, Mikrowellen-Labor-Systeme) for 6.5 min and heated at 600 W. Afterward, the closed sample container was cooled to room temperature; the clear solution was diluted to 100 mL in a volumetric flask and analyzed with a Varian AA100-spectrometer. To determine the Na+-content by AAS, a five-point calibration curve (0.0 to 0.2 mg/100 mL) was measured by diluting a Na+ standard AAS-solution (1000 mg/mL) in nitric acid.

CHN-Analysis

A Elementar Unicode equipped with a combustion tube filled with tungsten(VI)-oxide-granules was used to analyze the elements C, H, and N at a combustion temperature of 1050 °C. Samples were dried at 120 °C prior to the measurement.

Physisorption Measurements

Argon adsorption measurements were performed using a Quantachrome Autosorb at 87.35 K using samples that were dried at 120 °C for 24 h in a high vacuum. The pore sizes and volumes were calculated using a nonlocal DFT model (software version 2.11, Ar on zeolite/silica, cylindrical pores, equilibrium model).

Conductivity Spectroscopy

To perform broadband conductivity measurements, the hectorite sample was first ground in a mortar before the sample was uniaxially pressed (0.5 tons) into a cylindrical pellet with a diameter of 5 mm and a thickness of approximately 1 mm; the exact thickness was determined with a vernier calliper. The pellet was equipped on both sides with Au electrodes that blocked Li+ ion transport. For this purpose, we used an EM SCD 050 sputter coater from LEICA. The thickness (100 nm) was controlled with a quartz thickness monitoring system. The pellet was then placed inside an active ZGS cell (Novocontrol), which was connected to a Novocontrol Concept 80 broadband dielectric spectrometer. We measured the complex impedance, conductivity, and permittivity as a function of temperature and in a frequency range of 0.01 Hz to 10 MHz. The ZGS sample holder was placed in a cryostat, allowing us to measure conductivity isotherms from 133 to 433 K in steps of 20 K. Temperature regulation and monitoring was carried out with a QUATRO controller (Novocontrol). In order to eliminate residual moisture, the pellet was pre-dried in the impedance cell at 433 K (160 °C) for 4 h. During the entire measurement, the impedance cell was permanently flushed with dry, freshly evaporated nitrogen gas.

Magic Angle Spinning (MAS) NMR

6Li (73.6 MHz, 120 W power amplifier, 1024 scans, recycle delay 5T1 (see below), pulse length 3 μs) and 19F (470.6 MHz, 50 W, 32 scans, recycle delay 1 s, pulse length 2.1 μs) MAS one-pulse NMR spectra were recorded with a 500-MHz Avance spectrometer (Bruker). The measurements were carried out at a rotation speed of 25 kHz using 2.5 mm rotors. The spectra were recorded using one-pulse sequences with ambient bearing gas. Crystalline LiCH3COO and liquid CFCl3 served as reference materials to determine the isotropic chemical shifts δiso.

Time-Domain 7Li NMR Measurements

To prepare the sample for NMR measurements, the hectorite powder sample was ground using a mortar and pestle under an Ar atmosphere in a glove box. After this, it was dried in a vacuum at 160 °C before being fire-sealed in Duran ampoules under vacuum. We recorded both 7Li NMR spectra and 7Li NMR spin–lattice relaxation rates as a function of temperature (173 to 433 K). For this purpose, we used a Bruker 300-MHz spectrometer that was connected to a Bruker cryomagnet with a nominal magnetic field of 7 T. This field corresponds to a nominal 7Li Larmor frequency of ω0/2π = 116 MHz. While laboratory frame spin–lattice relaxation rates (1/T1) were recorded with the saturation recovery pulse sequence,[67] the spin-lock technique[67,68] was applied to measure transversal magnetization transients leading to the rate 1/T1ρ, which characterizes spin–lattice relaxation in the so-called rotating frame of reference. The 90° pulse length (tp) was 2 μs (200 W broadband amplifier) and showed only a slight dependence on temperature. To measure 1/T1, we used a comb of 10 closely spaced (80 μs) radio frequency pulses to destroy any longitudinal magnetization M. The subsequent recovery of M as a function of the delay time td (4 scans) was followed with a 90° detection pulse; the delay between each scan, although not needed, was set to 1 s. In general, M(td) follows a stretched exponential behavior according to M(td) = M0(1 – exp[−(td/T1)γ′]), 0<γ′≤1. M0 denotes the equilibrium magnetization that is reached at td → ∞. Here, M(td) follows a bi-exponential behavior if the full area under the free induction decays (FIDs) is used to construct M(td), see the Supporting Information, Figure S2. In this case, we used a sum of stretched exponentials to parameterize the transient.[31,69,70] Alternatively, the two components of the FIDs were separately analyzed, as explained below. For the corresponding spin-lock Mρ(td) transients (see also Figure S2), we used a single, stretched exponential function to analyze the decay functions: Mρ(tlock) ∝ exp[−(tlock/T1ρ)γ]), tlock denotes the duration of the locking pulse. We varied tlock from 30 μs to 30 ms (8 scans per locking time) and used a power level that results in a locking frequency ω1/2π of ca. 20 kHz. The recycle delay between each scan was at least 5T1. This value is orders of magnitude lower than ω0/2π used for the laboratory-frame measurements; hence, spin-lock NMR is sensitive to much slower Li ion dynamics compared to probing via T1 measurements.

Results and Discussion

Synthesis of the Li-Fluorohectorites

Direct melt synthesis of Li-fluorohectorites is troublesome. In general, melt-synthesized Li-fluorohectorites may contain a range of impurities and could suffer from isomorphous substitution effects. Thus, the charge density of the silicate layers turned out to be heterogeneous as described earlier.[71] In contrast to the Li-bearing compound, melt synthesis of Na-fluorohectorite followed by long-term annealing yields phase pure materials with a homogenous charge density and, thus, a uniform intracrystalline reactivity.[54,63] This uniform intracrystalline reactivity is especially important for a uniform cation distribution in the interlayer space, which in turn is expected to have a considerably high impact on cation transport properties. Therefore, the Li-fluorohectorite used for the present study was prepared by taking advantage of a complete ion exchange of Na-fluorohectorite. Unfortunately, the affinity of Li+ for the interlayer space of hectorites is lower than that of Na+.[72] Consequently, complete ion exchange could not be accomplished directly, even at very high incoming Li+-concentrations. Therefore, Na+ was fully exchanged with C4 taking advantage of the much higher selectivity of organic cations for the interlayer space.[73] The completeness of ion exchange at this stage was verified by AAS.

In the next step, Li exchange to Li-fluorohectorite could be driven by deprotonation of the intercalated C4. Elemental analysis proved complete replacement of C4 within experimental errors. At low ionic strength, i.e., at values being smaller than 0.02 M as achieved during the washing procedure,[63] Li-hectorite delaminates spontaneously in aqueous suspension by repulsive osmotic swelling whereupon a liquid crystalline gel is obtained. Freeze drying of this gel leads to a voluminous, spongelike, and hierarchically porous material (Figure b).

Figure 2

(a) Powder X-ray diffraction pattern of dry Li-fluorohectorite Li0.5[Mg2.5Li0.5]Si4O10F2 with marked reflections of the 00 L series and the hk bands; the inset shows the representative bandlike structure of the hectorite. (b) Scanning electron micrograph of freeze-dried and partially restacked Li-fluorohectorite.

The powder X-ray diffraction pattern of dry Li-fluorohectorite is shown in Figure a and shows the characteristic reflections of the 00 L-series with a basal spacing of 10.2 Å. Missing hkl reflections and occurrence of the λ-shape of 02/11-, 13/20-, 04/22-, 15/24-, and 06-bands indicated a turbostratic type of stacking of adjacent silicate layers. Upon drying, partial restacking is forced with large overlapping areas between adjacent individual layers (20 μm diameter) and leads to bandlike aggregate structures (Figure b). Consequently, Li+ cations are residing at external basal planes and in the interlayer space. The ratio of the two Li sites can only be estimated. The thickness of the restacked tactoids cannot be determined by analyzing the full width at half maximum (fwhm) of basal reflections via the Scherrer equation because the broadening is largely determined by random interstratification of the slightly varying stacking vector. Because of turbostratic disorder, the PXRD pattern contains little structural information beyond basal spacing, which is obvious when comparing a simulated XRD pattern (see Figure S1) for a hypothetical, fully 3D ordered Li-hectorite structure (C2/m (no. 12), a = 5.25(2) Å, b = 9.08(3) Å, c = 10.37(4) Å, β = 96.51(8)°). Since the hk-bands are located approximately at the position of hk0 reflections, this unit cell could be refined by applying the experimental turbostratic pattern, although with some uncertainty in the ab-dimensions. As ion exchange is a topotactic reaction, the structure of the silicate layers is fully preserved as indicated by ab-dimensions being similar to dimensions found for a 3D ordered member of the hectorite family (a = 5.2434(10) Å; b = 9.0891(18) Å).[74] The latter, refined structure of the silicate layer could, therefore, be safely applied to obtain Figure and to draw the scheme shown in the inset of Figure a.

As the interlayer surface area in the collapsed dry state is not accessible to probe gas molecules, the degree of restacking can, however, be deduced from the specific surface area determined by physisorption. Given the large diameter, the contribution of edges to the total surface area can be safely neglected. The specific surface area per gram for an individual delaminated silicate layer was calculated from the known ab dimension and the formula weight. The ratio of this calculated surface area (759.368 m2/g) and the surface area measured by Brunauer–Emmett–Teller analysis (5.174 m2/g) yields the average number nl of layers in the stack, nl = 147. This estimation indicates that most of the Li+ ions are indeed located in the interlayer space.

Ionic Transport and Li+ Diffusivity

In Figure a, the broadband conductivity isotherms of layer-structured Li0.5[Mg2.5Li0.5]Si4O10F2 are shown. The isotherms, which show the real part of the complex conductivity σ’ as a function of frequency ν, were recorded over a dynamic range of 9 decades.

Figure 3

(a) Conductivity isotherms of layer-structured, hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2 recorded from T = 133 K to T = 433 K; the measurements span a frequency range of 9 decades. In regime I, the isotherms are governed by the grain boundary (g.b.) response, while in regime II, bulk properties dominate the shape of the curves. See text for further details. (b) Arrhenius plot showing both log10(σ’DCT/(S cm–1 K)) and log10(1/τM / s–1) vs. the inverse temperature 1000/T. The left axis refers to σ’DC showing either the temperature behavior of the bulk ion conductivity or that governed by the highly resistive g.b. regions. The latter denotes the total conductivity seen in regime I of the isotherms shown in (a). 1/τM refers to the two maxima (i = 1, 2) seen in the electric modulus curves M”(ν), see Figure b, whereas 1/τM1 refers to the g. b. response and 1/τM1 denotes the characteristic electrical relaxation rates of the M”(ν)-peak representing the bulk response. Lines, either dashed or dotted, show Arrhenius fits with the activation energies indicated. At room temperature, σ’DC, bulk ≈ 0.14 mS cm–1 is expected through extrapolation of the Arrhenius line shown. For comparison, the rates 1/τρ from electrical resistivity measurements carried out at a fixed frequency (1.2 MHz, 10 MHz, see Figure b) are included as well.

Figure 4

(a) Electrical modulus curves M”(ν) of hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2, which were recorded at the temperatures indicated. (b) Same curves as in (a) but using a double logarithmic plot to visualize also those peaks that are governed by a higher permittivity (and capacitance), see arrow pointing at the peak that appears at 0.03 Hz (233 K). The latter reflects g.b. responses, while the main modulus peaks are produced by bulk electrical relaxation processes characterized by a much lower permittivity, see the arrow pointing at the peak at ν = 8 × 105 Hz. The dashed lines show the response of a measurement on a second pellet prepared, indicating the good reproducibility of the analysis.

Figure 5

(a) Change of the electrical permittivity of layer-structured Li0.5[Mg2.5Li0.5]Si4O10F2 as a function of frequency and temperature. While the bulk region (ε’b(0) = 10) is hardly seen in this representation, the g.b. response (ε’gb(0) = 360) and electrode polarization (ε’p(0)>106) dominate the isotherms, which were recorded from 133 to 433 K in steps of 20 K. (b) Temperature behavior of the electric resistivity is expressed as M”/ω and measured at two different frequencies ω/2π = 1.2 and 10.0 MHz. The fact that the flanks at high temperatures do not coincide is an indication of low-dimensional ion transport; 3D ion dynamics would yield the behavior indicated by the dashed-dotted lines. From the peak maxima, two electrical relaxation rates can be deduced that are also included in Figure b. See text for further explanations.

The conductivity curves can be subdivided into several parts, labelled A, B, and C (Figure a). To explain these features, we first consider the isotherm recorded at the highest temperature, that is, at T = 433 K (regime I). At low frequencies ν, electrode polarization appears, which manifests itself as a drop in conductivity (A). With increasing frequency, σ’ enters the so-called direct current (DC) plateau that corresponds to long-range ion transport in Li0.5[Mg2.5Li0.5]Si4O10F2 (B). The associated capacitance C of this distinct plateau is in the order of 300 pF (= Cgb). Without any doubt, such a high value corresponds to an electrical process that is dominated by g.b. regions.[75] Increasing the frequency further leads to the associated dispersive part (C) of the DC plateau, which follows a σ’ ∝ νp = 1 behavior as indicated by the dashed line in Figure a.

With decreasing temperature, σ’DC, gb decreases and the curves are shifted toward lower frequencies. At sufficiently low T, a second DC plateau could be revealed (see regime II). This plateau does only appear as a curvature of the isotherms shown in regime II of Figure a. Its dispersive regime is only hardly visible; at very low temperatures and high frequency, the exponent p in the Jonscher-type relation σ’ ∝ νp approaches 1, indicating the beginning of a nearly constant loss behavior.[76−81] Most importantly, while total conductivities, dominated by the g.b. response, can easily be read off from the distinct plateaus in regime I (see Figure a), σ’DC, b values characterizing bulk ion transport are only accessible by careful evaluation of the curvature seen in regime II. Filled symbols in Figure a mark estimations of these values.

The temperature dependence of both values, σ’DC, b and σ’DC, gb, is analyzed using the Arrhenius plot shown in Figure b. Conductivities dominated by ion-blocking grain boundaries obey the Arrhenius law with an activation energy Ea gb of 0.58 eV. Lines in Figure b refer to σ’DCT = σ0 exp. (−Ea/(kBT)) with σ0 being the pre-exponential factor and kB denoting Boltzmann’s constant. Compared to σ’DC, gb the bulk conductivity values σ’DC, b turned out to be higher by 6 orders of magnitude, and bulk ion transport has to be characterized by an activation energy of Ea, b = 0.36 eV. Extrapolating the estimated values σ’DC, b toward room temperature yields 0.14 mS cm–1 at 293 K. This value is comparable to those reported for garnet-type oxides such as cubic-Li7La3Zr2O12 being considered as powerful electrolytes for ceramic electrochemical energy storage systems.[10]

To support our claim that the curvatures in regime II of Figure a are suitable to estimate bulk ion conductivities of hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2, we analyzed electric modulus curves, M’’(ν), which we recorded at a fixed temperature and variable frequencies.[48] The corresponding curves are shown in Figure a by using a half-logarithmic plot and in Figure b by taking advantage of the double-logarithmic scaling to illustrate the change of M’’ with frequency.

In general, the complex modulus is given by the inverse complex permittivity, M* = 1/ε*. The real part of ε*, denoted as ε’, is shown in Figure a. As seen for σ’, also ε’(ν) passes through the same electrical regimes viz. electrode polarization, g.b., and bulk response. The latter is only slightly seen for the isotherms recorded at very low temperatures. The permittivity value ε’b(ν → 0) = 10 (see Figure a) corresponds to Cb< 10 pF if we use the relation for a plate capacitor, C = ε0ε’bA/d, to estimate the associated bulk capacitance. Here, ε0 is the vacuum permittivity, A denotes the area of the pellet and d its thickness. Importantly, the imaginary part of M*, M’’, is dominated by bulk effects as the amplitude of M’’ is proportional to the inverse of Cb, (gb), M’’ ∝ 1/C. Hence, we expect a clearly visible bulk M’’(ν) peak accompanied by a second one, the g.b. response, being drastically reduced in magnitude. As compared to the main peak, the g.b. peak (denoted as peak 1 in the following, see Figure b) is expected to be shifted toward lower frequencies by at least 6 orders of magnitude, as σ’DC, gb/σ’DC, b ≈ 10.[6] While in the semilog plot, see Figure a, only the main peak (peak 2) is seen, in the double-logarithmic representation of the data, two modulus peaks are indeed recognizable, see the arrows in Figure b that exemplarily point to these two peaks that belong to the curve measured at 233 K. The ratio of the peak amplitudes is given by M’’max(peak 2)/M’’max(peak 1), which, in the present case, amounts to 30–40. Of course, this value reflects the ratio ε’gb/ε’b (see Figure a) and the ratio Cgb/Cb as well. Note that these values are rough estimations, as they assume Debye-like responses. In general, deviations from Debye-like behavior are seen if correlated motion or, as argued in the following, low-dimensional transport governs the electrical relaxation processes.[82,83]

Here, the most important parameter, which we deduce from the two peaks of the modulus curves shown in Figure b, is the characteristic relaxation rate τM–1. It corresponds to the frequency at which the peak appears, i.e., at M’’ = M’’max. In Figure b, log10(τM1–1) and log10(τM2–1) of the two peaks 1 and 2 are plotted vs 1000/T. τM2–1(1/T) reveals almost the same activation energy (0.35 eV) as σ’DC, b, thus supporting our interpretation of σ’DC, b to represent a parameter characterizing a bulk response. The difference in activation energies is less than 0.015 eV, which is within the error range of the values. The same holds for the activation energies obtained for τM1–1 and σ’DC, gb (0.60 eV vs 0.58 eV). Note that the rates τM–1 shown in Figure b refer to technical frequencies νmax; using angular frequencies instead (ωmax = 2π νmax) would simply shift the Arrhenius line by an additive constant (log10(2π) ≈ 0.8) upward on the log10 scale.

In addition to the possibility to extract characteristic relaxation rates from M’’(ν) peaks, we measured the quantity M’’/ω as a function of the inverse temperature at a fixed frequency.[52,84,85]M’’/ω corresponds to the resistivity ρ’ and is expected to pass through an electrical relaxation peak that should be parameterizable with a Lorentzian shaped function: ρ’ = M’’/ω = τρ/(1 + (τρω)β) with 1<β≤2. Here, τρ denotes the electrical relaxation time, and β is a parameter that expresses the deviation of the peak, log10(ρ’) vs 1/T, from symmetric behavior. Symmetric rate peaks are obtained for β = 2; β<2, however, indicates correlated motion affecting the electrical process probed. In the case of 3D motions, the latter case would result in asymmetric peaks with the so-called low-temperature flank, showing a lower slope than the high-temperature flank. If measured at different frequencies, the rates in the latter regime are expected to coincide. This feature is indicated in Figure b by the dashed-dotted lines. The activation energy of the high-T flanks should correspond to Ea sensed by σDC b and τM2–1, which are characterized by 0.36 eV. Astonishingly, in the present case, Ea in the high-T limit ranges from only 0.06 to 0.08 eV. Moreover, the high-temperature flanks of the two peaks clearly do not coincide. As an example, at 330 K, the corresponding values M’’/ω differ by more than one order of magnitude revealing a dispersive behavior of the underlying electrical relaxation function. We attribute this deviation seen by the resistivity measurements to the low-dimension, i.e., 2D ionic transport in the Li-bearing phyllosilicate Li0.5[Mg2.5Li0.5]Si4O10F2. This assignment is considered to be in analogy to the interpretation of diffusion-induced NMR spin–lattice relaxation rate (1/T1) peaks as a function of the inverse temperature at various frequencies. In the case of, e.g., the 2D Li-ion conductor LiTiS2, the dispersive behavior of the NMR correlation function in the high-T, that is, the low-frequency, limit was in quantitative agreement with the NMR theory.[40] Although, in principle, NMR and conductivity relaxation are governed by different correlation functions,[83,86] the conclusion that the high-T flank is sensitive to the dimensionality of diffusion holds qualitatively also for the resistivity.

At the peak maximum, i.e., at T = Tmax, the condition τρω ≈ 1 is valid. With ν = ω/2π for each frequency (ν = 1.2 MHz and ν = 10.0 MHz), a rate τρ–1(1/Tmax) can be estimated, which we included in Figure b as well. τρ–1 agrees with τM2–1deduced from the M’’(ν) peaks if we extrapolate the corresponding Arrhenius line toward higher temperatures. Hence, long-range ion transport in the bulk of Li0.5[Mg2.5Li0.5]Si4O10F2 is given by an activation energy of 0.35 eV. This value is comparable to those energies ordinarily seen for Li+ ionic transport in, e.g., garnet-type oxides.[10,87,88] Identifying Ea = 0.35 eV with the activation energy, which should actually be seen in the high-T regime of the ρ’(1/T) peak, the parameter β would take a value of approximately 1.29, as β is given by 0.10 eV/0.35 eV = β – 1, if we adopt the analogy between conductivity and NMR relaxation more quantitatively.[5,30] In this relation, the value of 0.10 eV represents the activation energy in the low-T regime of the ρ’(1/T) peak, see Figure b. β = 1.29 indicates rather strong ion–ion correlation effects. While such effects reduce the slope of the low-T flank, dimensionality effects will only influence the rates in the high-T regime. Semi-empirical models for 2D diffusion predict that the slope of the high-T flank should take a value that is approximately 75% of the slope in the low-T limit.[62] Here, the combination of ion–ion correlation and 2D diffusion results in apparently symmetric peak shapes even for low-dimensional dynamics. A very similar behavior has been observed earlier by Li+ nuclear spin-relaxation in LiTiS2.[40]

In conclusion, we found evidence that a 2D ionic transport is present in the phyllosilicate under study, which is governed by a rather low activation energy of only 0.35 eV. From Figure b, we see that the ratios τM2–1/τM1–1 (= τM1/τM2) and σ’DC, b/σ’DC, gb differ by at least one order of magnitude. This observation is independent of T, as the differences in Ea for the two processes show only little change with temperature. Hence, we might have underestimated rather than overestimated the bulk ion conductivity, which we read off from the σ’(ν) curvatures seen in regime I (Figure a) by one order of magnitude. Such a rapid ionic transport process is expected to be measurable also in 7Li nuclear spin relaxation. Indeed, both NMR line-shape measurements and spin-lock NMR spin–lattice relaxation rates point to temporary spin fluctuations being extremely rapid on the NMR time scale.

In Figure , 7Li NMR line shapes are shown, which we recorded at temperatures ranging from 173 to 433 K. As can be seen in Figure b, the line consists of two components. A broader line with a width of ca. 12 kHz is superimposed by a narrower one whose width narrows with temperature because of increasing dipole–dipole averaging. Such averaging originates from motional processes rendering the homonuclear Li–Li dipolar interactions time dependent. The broad line can either be interpreted as a result of electric quadrupole interactions of the Li spins with nonvanishing electric field gradients or as a signal representing a group of Li ions with much lower diffusivity. As we do not recognize the emergence of a sharp quadrupole powder pattern at high temperatures, which would be in line with the universal temperature behavior of 7Li NMR line shapes affected by averaging processes,[89,90] we assume that the broad line originates from a slow Li+ subensemble. Most likely, this ensemble represents those ions sharing the same crystallographic position as Mg2+. We assume that the Li ions located at 2a, i.e., between the layers, give rise to the narrow line observed.

Figure 6

(a) 7Li NMR line shapes recorded at the temperatures indicated. At 173 K, the fwhm of the superimposed complete line is given by 2.08 kHz. It reduces to 516 Hz at 413 K. This overall line is composed of two contributions; a deconvolution is seen in (b). In (b), the change of fwhm, either determined by using the full line shape or by analyzing the narrow component only, is shown as a function of temperature. Inset: deconvoluted 7Li NMR spectra with the area of the Lorentzian (narrowed line) taking 28%, i.e., ca. 30%, of the full area under the line. See text for further explanation.

To analyze the change in line width, we read off the fwhm of the total line; the line width as a function of temperature is shown in Figure . Starting from ca. 2 kHz, it reduces to values of ca. 500 Hz at 413 K, which is the so-called regime of extreme narrowing. We observe that the line width sharply increases at low temperature reaching 0.75 kHz at 250 K, at even higher temperatures, a shallower decrease is seen, which is still not fully completed at T = 450 K. To find out whether this is an exclusive feature of the overall line, we deconvoluted the lines with the help of a sum of a Lorentzian and a Gaussian function. Here, the narrow, Lorentzian-shaped line reveals the same temperature dependence as the overall line. The mainly homonuclear interactions of this line are averaged when T = 250 K is reached. Further averaging of residual couplings, likely between Li(2a) and Li(4 h), is seen at higher temperatures. Importantly, even at T = 173 K, the line has not reached its temperature-independent rigid-lattice value, i.e., rather rapid Li+ exchange processes are present in Li0.5[Mg2.5Li0.5]Si4O10F2. We estimate that at 225 K, the associated average jump rate 1/τ should be in the order of 10 kHz or higher, as full averaging will be achieved if 1/τ greatly exceeds the spectral width of the NMR signal in the rigid lattice regime.

The deconvolution of the 7Li NMR spectra tells us that approximately 30% (see the deconvoluted spectra shown in the inset of Figure b) of the total number of Li ions in Li0.5[Mg2.5Li0.5]Si4O10F2 have access to fast diffusion pathways, see the inset of Figure b. This value might be even higher, as we cannot exclude that the broad component includes some quadrupole satellite intensities associated with the narrow line. In such a case, the value might reach 50%, which would be in line with what the chemical formula Li0.5[Mg2.5Li0.5]Si4O10F2 suggests. To shed light on this assumption, we recorded high-resolution 6Li MAS NMR spectra, see Figure , for which first-order quadrupolar interactions and dipole–dipole broadening, both interactions being smaller for 6Li anyway, are artificially averaged out through fast sample rotation. In agreement with the static line shapes, the 6Li MAS line can be best approximated with a superposition of two distinct lines whereby the area under the narrow one again amounts to approximately 30%. A certain extent of Li–Mg exchange could serve as an explanation to understand the reduced spin density between the layers in Li0.5[Mg2.5Li0.5]Si4O10F2. For the sake of completeness, in Figure , the corresponding 19F MAS NMR spectrum is shown, which reveals the two magnetically inequivalent F sites in Li0.5[Mg2.5Li0.5]Si4O10F2.

Figure 7

(a) 6Li and 19F MAS NMR spectra of Li0.5[Mg2.5Li0.5]Si4O10F2 recorded under ambient bearing gas conditions (300 K). The 6Li MAS NMR spectrum can be best approximated with a combination of two spectral components as indicated. A deconvolution was possible with a sum of a Lorentzian line (narrow line) and a Gaussian one (broad line). Values denote line widths. Isotropic shifts were allowed to float. Approximately 32% of the Li ions, see also Figure b, are responsible for the narrow line. (b) In 19F MAS NMR, two lines with equal intensities are seen reflecting the two magnetically inequivalent F sites in Li0.5[Mg2.5Li0.5]Si4O10F2. Asterisks mark spinning sidebands appearing at an interval of 25 kHz for each F signal.

The two-component features seen in both 7Li and 6Li NMR also affect the evolution of the magnetization transients, recorded in the former case, to extract diffusion-induced 7Li NMR spin–lattice relaxation rates. As mentioned above, M(td) can only be parameterized satisfactorily with a single, stretched exponential at sufficiently low temperatures (T<225 K). Above this threshold, which coincides with the beginning of the shallowly decaying part of the NMR line width (see Figure b), the transients start to follow a two-exponential time behavior, see also Figure S2. Therefore, we used a sum of two stretched exponentials, see the Supporting Information, to parameterize the overall transients; the same procedure, using the same terminology, is described in detail elsewhere.[31,70] The two rates 1/T1,slow and 1/T1,fast are shown in the Arrhenius plot of Figure a. The corresponding stretching exponents are displayed in the upper part of Figure . While 1/T1,slow does not depend on temperature within the error limit of the analysis, 1/T1,fast reveals a shallow diffusion-induced peak at Tmax = 385 K. The low-T flank points to a very low activation energy of only 0.06 eV, which is likely to be dominated by strictly localized, i.e., short-range, Li+ hopping processes and/or correlation effects.[30,32] Almost the same value was probed by the resistivity measurements, see Figure b.

Figure 8

(a) 7Li NMR spin–lattice relaxation rates of Li0.5[Mg2.5Li0.5]Si4O10F2 recorded in the laboratory frame of reference at a Larmor frequency of 116 MHz. The lower part shows an Arrhenius plot. The two-component 7Li NMR line shape produces bi-exponential magnetization transients M(td) leading to two rates 1/T1fast and 1/T1slow that reflect the fast and slow ions in Li-bearing fluorohectorite, also seen in 7Li NMR lines. Alternatively, the stepwise analysis of the underlying FIDs also yields two relaxation rates. Stretching exponents γ are shown in the upper part of the figure. (b) Comparison of the spin-lock 7Li NMR spin–lattice relaxation rates (1/T1ρ, 20 kHz) with those measured in the laboratory frame of reference. The rate 1/T1ρ passes through a diffusion-induced maximum at approximately 205 K. At this temperature, the corresponding γ passes through a minimum. Lines (dashed and solid) are to guide the eyes. Linear parts were analyzed with the Arrhenius law to extract the activation energies (in eV) as indicated.

Alternatively, we separated the two spectral contributions directly in the time domain of the underlying FIDs by individual analysis of the corresponding components that reveal quite different effective spin–spin relaxation rates. This procedure is described in detail elsewhere.[31,69,70,91] The FIDs are composed of a sharply (FIDfast) and a slowly decaying part (FIDslow), leading to two distinct rates, see Figure a. While the first decay reflects the broad spectral component, the latter decay, FIDslow, arises from the motionally narrowed component that represents fast Li+ ions. At least for high temperatures, the rates associated with FIDslow agree with those obtained when analyzing the full M(td) curves. Again, a relatively low activation energy of 0.12 eV can be estimated from the corresponding low-T flank. This value resembles that also seen in the resistivity peaks ρ’(1/T) shown in Figure b, cf. the data recorded at 1.2 MHz. For the separated rates of the slow subensemble, a temperature independent behavior is seen for T>330 K.

To further characterize Li+ ion hopping with methods that are also able to probe long-range ion dynamics, we performed variable-temperature spin-lock NMR measurements.[92] The corresponding 1/T1ρ rates are shown in Figure b. The magnetization transients (Figure S2) do not allow for a separate investigation as it has been carried out for the 1/T1 measurements. Hence, they are to be regarded as averaged values. Likely, the fast Li+ ions will dominate the overall nuclear spin response. Coming from low temperatures, the rates experience a sharp increase with temperature and pass through a maximum at T = 210 K; the maximum in1/T1ρ corresponds to a minimum in the stretching factor, see arrow in the upper graph of Figure b. The slope of the low-T flank of this prominent peak yields an activation energy of 0.39 eV, which is very similar to that seen by both conductivity spectroscopy and the modulus analysis, cf. the temperature dependence of σ’DC, bT, and τM2–1 (Figure b). Thus, we found a strong indication that the peak probed by 1/T1ρ, which appeared at 205 K, reflects the same temporary fluctuations as detected by electrical relaxation. Interestingly, the activation energy of its high-T flank is comparable with those also characterizing the modulus peaks in Figure b. A value of 0.08 eV is also very similar to the estimated ones from 1/T1 NMR measurements (0.06 and 0.12 eV, see Figure a). Assuming a symmetric 1/T1(1/T) peak, the activation energies in the high-T limit would be characterized by the same activation. In general, this similarity would be expected if the two methods probe the same (low-dimensional) diffusion process.[53] Importantly, dimensionality effects affect the slope of the NMR rate peaks in the high-T regime and yield apparent, reduced values.[32,43]

In summary, electrical modulus measurements and spin-lock 7Li NMR support the finding that long-range ion transport in hectorite-type Li0.5[Mg2.5Li0.5]Si4O10F2 is characterized by an activation energy of approximately 0.35 eV (see Figure b). In contrast to long-range transport, barriers with values of approximately 0.1 eV govern the local hopping processes. However, the values from high-T flanks, either seen in 1/T1ρ or in M’’(ν), are affected by dimensionality effects. As also recognized for the M’’(ν) peaks, we observe that the rates 1/T1ρ of the high-T flank do not coincide (see arrow in Figure b), with those expected from 1/T1 as is indicated by the dashed line in Figure b. This line extrapolates the high-T flank to even higher temperatures. The same feature is seen in Figure b. This dispersive behavior strongly supports our idea about 2D diffusion in the synthetic hectorite studied. Finally, 1/T1ρ(1/T) revealed a second peak appearing at Tmax = 350 K (marked by a star in Figure b), which, most likely, is to be assigned to the low translational ion mobility of the Li+ ions sharing sites with Mg2+ in the layers of the phyllosilicate.

To directly compare the different responses seen by NMR and electrical relaxation, we finally estimated Li+ jump rates from the maxima seen in 1/T1 and 1/T1ρ at Tmax = 385 K and Tmax = 205 K, respectively, by using the maximum condition τNMRω0(1) ≈ 1. The rates obtained are included in Figure . In the present case, they agree very well with those from electrical measurements. At 293 K, the Arrhenius line yields a mean jump rate 1/τ of approximately 3 × 107 s–1, which translates into a 2D diffusion coefficient D2D = a2/(4τ) of 4.7 × 10–13 m2 s–1 if we use a ≈ 2.5 Å. The crystallographic Li–Li distance is ca. 5 Å, which we think is too large for a single Li+ jump. Assuming temporarily occupied interstitial positions connecting two or more regular sites in a distance of 2.5 Å seems to be more reasonable. D = 4.7 × 10–13 m2 s–1 would indeed correspond to an estimated ionic conductivity σest., almost reaching the order of 0.1 mS cm–1 if we use the (idealized) cell volume of 490 Å3 to estimate a charge carrier concentration of ca. N1 = 2 × 1027 m–3. σest. is related to DNE via the Nernst-Einstein relation, which reads in general DNE = σkBT/(N1q), with q being the charge of the Li+ ions. We assumed that D2D and DNE are connected by D2D = (HR/f)DNE with HR being the Haven ratio and f denoting the correlation factor. Our estimation anticipated HR/f being of the order of unity. Taking into account that only 30% of the ions participate in Li+ long-range ion transport, D2D = (HR/f)DNE is fulfilled if we chose HR/f ≈ 0.3. The relatively high value of D2D obtained clearly renders hectorite as a fast ion conductor.

Figure 9

Arrhenius plot of the Li+ hopping rates 1/τ deduced from both nuclear spin relaxation and electrical relaxation taking advantage of the modulus representation and resistivity measurements, see Figure . Data from NMR, if restricted to the maxima of the spin–lattice relaxation rate peaks, agree well with the results from broadband electrical characterization. The dashed line represents a linear fit yielding 0.35(1) eV. For comparison, jump rates measured for LiTiS2 (x = 0.7, 0.41 eV)[43] and LiC6 (x = 1, 0.55 eV)[38] are also shown, see dotted lines.

Rapid ion exchange in Li0.5[Mg2.5Li0.5]Si4O10F2 benefits from a wide interlayer gap giving the ions the necessary space to freely diffuse in the phyllosilicate. To the best of our knowledge, this is the first time that such high cation diffusivities were probed in a synthetic clay mineral. For comparison, in Figure , we included jump rates of two other 2D materials, LiTiS2 (0.41 eV)[43] and LiC6 (0.55 eV),[38] studied by NMR (and other methods) earlier. Inspired by the present results, one might think about crystal-chemical modifications that could lead to even higher diffusivities. Nature may offer a range of further silicate structures, including also layered ones, from which the necessary inspiration can be drawn to develop both powerful electrolytes and active materials for, e.g., lithium-ion batteries. Exemplarily, the Khan group studied hectorite-containing composite materials for this purpose.[56−59]

Conclusions

Li+ ion diffusion and electrical transport in the hectorite-type phyllosilicate Li0.5[Mg2.5Li0.5]Si4O10F2 was studied by broadband conductivity spectroscopy, modulus analysis, and 7Li NMR spin–lattice relaxation measurements. Electrical and nuclear spin relaxation confirmed our hypothesis of rapid interlayer (2D) Li+ exchange processes in the silicate. This process is characterized by an overall activation energy of approximately 0.35 eV. Under ambient conditions, conductivity spectroscopy points to an ionic conductivity as high as 0.14 mS cm–1, representing a favorable starting point for further improvements by crystal-chemical engineering. In addition, we derived a consistent picture of Li+ ion dynamics and showed that the electric and magnetic fluctuations probed originate from the same translational process. Our study supports the general idea that spatial confinement, able to guide the charge carriers over long distances, is helpful in enabling fast ion transport.