Long-range Li + dynamics in the lithium argyrodite Li 7 PSe 6 as probed by rotating-frame spin–lattice relaxation NMR †

Lithium-rich argyrodites belong to a relatively new group of fast ion conducting solids. They might serve as powerful electrolytes in all-solid-state lithium-ion batteries being, from a medium-term point of view, the key technology when safe energy storage systems have to be developed. Spin–lattice relaxation (SLR) nuclear magnetic resonance (NMR) measurements carried out in the rotating frame of reference turned out to be the method of choice to study Li dynamics in argyrodites. When plotted as a function of the inverse temperature, the SLR rates log 10 ( R 1 r ) reveal an asymmetric diﬀusion-induced rate peak. The rate peak contains information on the Li jump rate, the activation energy of the hopping process as well as correlation eﬀects. In particular, considering the high-temperature flank of the SLR NMR rate peak recorded in the rotating frame of reference, an activation energy of approximately 0.49 eV is found. This value represents long-range lithium jump diﬀusion in crystalline Li 7 PSe 6 . As an example, at 325 K the Li jump rate determined from SLR NMR is in the order of 1.4 (cid:2) 10 5 s (cid:3) 1 . The pronounced asymmetry of the rate peak R 1 r (1/ T ) points to correlated Li motion. It is comparable to that which is typically found for structurally disordered materials showing a broad range of correlation times.


I. Introduction
2][3][4][5][6][7] In particular, compared to liquid electrolytes, the use of solid ion conductors meets the challenges to design more powerful and highly safe lithium-ion batteries. 7,8These are expected to show a long lifetime which is, in particular, necessary for large-scale, decentralised stationary energy storage.Thus, research on materials for lithium-ion batteries is in line with the imperative goal to reduce our dependence on fossil fuels if we are to store electrical energy from renewable sources and, hence, to combat climate change.Finding suitable solid electrolytes will certainly provide a leap forward towards the development of safe all-solid-state lithium-ion batteries being useful for diverse applications of energy storage systems such as portable electronic devices or even (hybrid) electric vehicles.
During the last few years many Li containing solid electrolytes have been screened and also modified for such purposes.However, the number of suitable candidates turns out to be rather rare, see, e.g., ref. 8-14.Essentially, besides chemical inertness, which oxide materials usually show, they have to fulfill the necessary requirement to exhibit an ionic conductivity exceeding, in the ideal case, 10 À3 S cm À1 at room temperature.
Quite recently, a new class of Li-rich materials has been investigated, viz.crystalline lithium argyrodites, which exhibit such a sufficiently high Li + mobility. 9,156][17][18][19] Apart from first line shape analyses, 9,15 an in-depth Li nuclear magnetic resonance (NMR) spin-lattice relaxation (SLR) study is still missing.It turned out that the highest conductivity was found for the halides rather than for the pure sulfides or selenides. 151][22][23][24][25][26] Compared to conductivity measurements, diffusion-induced Li NMR SLR rates solely reflect Li + diffusivity rather than e À dynamics which might also contribute to the total conductivity probed.In the latter case, contributions arising from other charge carriers have to be carefully eliminated which in some cases might be fraught with great experimental difficulties.By using NMR spectroscopy instead, no special post-preparation of the material to be investigated is necessary and a powdered sample is sufficient to study Li dynamics in detail.Usually the powder is sealed in a small glass ampoule in an inert gas atmosphere to protect it from moisture.][29] In the present case, the low-temperature modification of Li 7 PSe 6 , whose structure has been well characterized by both X-ray powder diffraction and high-resolution NMR recently, 15 was chosen to study dynamic properties by temperaturevariable 7 Li NMR spectroscopy.The material served as a first model substance to enlighten Li ion transport properties of lithium-rich argyrodites by NMR SLR methods.Here, NMR SLR rate measurements performed in the so-called rotating frame of reference turned out to be the best choice to study Li selfdiffusion processes.The method takes advantage of an oscillating magnetic field B 1 (rather than the ordinary external magnetic field B 0 ), which is used to perturb the 7 Li spin system in the frame of an SLR NMR experiment. 30,31The associated radio (or locking) frequency o 1 is related to B 1 via gB 1 = o 1 where g is the magnetogyric ratio of the nucleus under investigation.Since o 1 is in the order of some tens of kHz and o 0 is in the upper MHz range, NMR SLR measurements performed in the rotating frame of reference are sensitive to Li jump rates with values in the order of 10 4 to 10 5 s À1 , thus, considerable longer correlation times can be probed than it is possible with SLR NMR measurements performed in the laboratory frame of reference being the more common technique to probe ion dynamics by NMR. 25,28n general, monitoring the complete diffusion-induced SLR NMR rate peak, which is obtained when the SLR rate R 1r is plotted vs. the inverse temperature 1/T, gives access to both the activation energy E a and the pre-exponential factor t 0 À1 determining the underlying Arrhenius relation (see below).Furthermore, it contains information about correlation effects such as structural disorder and Coulomb interactions. 32Therefore, it might be useful to roughly evaluate to what extent shortrange Li motions influence the R 1r rates measured and, thus, ion transport parameters.

II. Experiment
Synthesis and characterization of polycrystalline Li 7 PSe 6 with mm-sized crystallites (as well as the closely related Li 7 PS 6 ) are described in great detail by Kong et al. 15 Both Li 7 PS 6 and Li 7 PSe 6 exhibit two structural modifications with a fully reversible phase transition at 483 K and 437 K, respectively.While the high-temperature modification shows the typical argyrodite-type structure (space group F% 43m), the crystal structure of the low-temperature modification is closely related to that of the low-temperature form of a-Cu 7 PSe 6 crystallizing with orthorhombic symmetry. 15 7 i NMR measurements were conducted on a high performance digital Bruker Avance III spectrometer in connection with a shimmed cryomagnet of a nominal magnetic field of ca.7 T.This corresponds to a Larmor frequency of o 0 /2p = 116 MHz.A commercial high-temperature probe head (Bruker Biospin) was employed allowing p/2 pulse lengths of about 2.9 ms.The measurements were performed at temperatures ranging from 216 K to 416 K, which were monitored using a Eurotherm temperature controller in combination with a type T thermocouple.Measurements at temperatures down to 112 K were carried out using a cryo probe (Bruker) which was operated by evaporating liquid nitrogen.As compared to the standard probe, the p/2 pulse length of the cryo probe was ca.5.2 ms.Prior to the measurements the sample was vacuum dried and fire sealed in a glass NMR tube of 5 mm in diameter and 3 cm in length. 7Li NMR SLR rates 1/T 1 = R 1 were acquired by means of the saturation recovery pulse sequence 10 Â p/2 À t d À p/2 which is depicted in Fig. 1(A). 20,33Here, the initial ten p/2 pulses in rapid succession are used to destroy any longitudinal magnetization M z before its recovery is recorded as a function of the delay times t d . 7Li NMR SLR rates 1/T 1r = R 1r were recorded in the rotating frame of reference (SLRr).To this end the spin-locking technique [30][31][32][33][34][35][36] was employed (see Fig. 1(B)) where the relaxation of the equilibrium magnetization is probed at rather low magnetic field strengths B 1 .The corresponding locking frequency was chosen to be o 1 /2p E 45 kHz with the duration of the locking pulse t lock varying between 10 ms and 100 ms.The recycle delay for the SLRr experiments was set to at least 5 Â 1/R 1 .Both R 1 and R 1r were obtained by parameterizing the magnetic transients M z (t d ) and M r (t lock ), respectively, by stretched exponentials.In the case of R 1r the exponent g in M r (t lock ) p exp(À(t lock /T 1r ) g ) ranges from 0.4 to 0.9 (210 K r T r 430 K).For the sake of completeness, temperature-variable 7 Li NMR spin-spin relaxation (SSR) rates 1/T 2 = R 2 were monitored, which were recorded using a solid-echo pulse sequence 33 optimized for spin-3/2 nuclei (see Fig. 1(C)).t echo denotes the variable interpulse delay.The corresponding NMR spectra were either obtained after irradiation with a single p/2 pulse or using the solid-echo pulse sequence.In the latter case, echoes were recorded at fixed t echo of 20 ms; spectra were obtained by Fourier transformation beginning from the echo top.
][44][45][46][47][48] The technique gives access to single-spin two-time correlation functions 21,22,39 and takes advantage of the different electrical environments the Li spins sense during hopping from site to site.Here, so-called 7 Li sin-sin hopping correlation functions, 21,39 i.e., spin-alignment echo (SAE) decay curves, 24,25,43 were monitored at a fixed preparation time t p of 12 ms and a variable mixing time t m ranging from 10 ms to 100 s.Appropriate phase cycling is necessary to eliminate unwanted coherences. 39

III. Results and discussion
Temperature-variable 7 Li NMR spectra of crystalline Li 7 PSe 6 yield first insights into the Li dynamics of the argyrodite.In Fig. 2 various NMR spectra are shown which were recorded at temperatures ranging from 216 to 417 K.They were acquired using a single pulse sequence, thus, quadrupole intensities (see below) are largely suppressed.Typically, at low temperatures, i.e. in the so-called rigid-lattice regime, a dipolar broadened NMR (central) line is observed which can be approximated with a Gaussian.Due to the increasing motion of the Li spins at elevated temperatures dipole-dipole interactions are continuously averaged.Hence, the central NMR line width obeys a motional narrowing (MN) from which the Li jump rate can be roughly estimated.In the present case, even at temperatures as low as 200 K MN has already set in.This can be clearly seen when Fig. 3 is considered where the NMR line width (full width at half maximum), directly read out from the NMR spectra recorded, is plotted vs. T.The dashed line is simply drawn to guide the eye; the rigid-lattice regime is reached at temperatures lower than 150 K.
Whereas at the highest temperatures (T > 350 K, that is, the regime of extreme narrowing) the NMR central line, corresponding to the spin-quantum transition I = AE1/2 $ I = 81/2, can be well approximated with a Lorentzian, the NMR line in the intermediate temperature range can be represented by the sum of a Lorentzian line and a Gaussian line.However, a strongly heterogeneous narrowing of the NMR central line, which in many cases is reflected by a pronounced difference in line width of the two components, 50,51 seems to be absent in the present case.Thus, the lack of clearly recognizable twophase spectra hampers a quantitative analysis according to that, e.g., presented in ref. 52.Nevertheless, the shape of the central transitions shown in Fig. 2(a) points to a distribution of Li jump rates and, therefore, also to a distribution of Li migration pathways in the Li-argyrodite studied.
Since the MN curve shown in Fig. 3 is obtained by plotting the overall line widths of the spectra recorded, it is only of limited use for the extraction of dynamic jump parameters.Significant MN is expected when the jump rate t À1 reaches the order of the rigid lattice line width.For example, one might take use of the inflexion point of the curve (T infl E 235 K) to estimate t À1 .Using n rl E 6.8 kHz for the NMR line width measured in the rigid lattice regime, the corresponding Li jump rate t À1 is approximately given by t Moreover, from the so-called onset temperature (T onset E 180 K) an activation energy E a,MN can be roughly estimated according to the expression of Waugh and Fedin, 53 E a,MN /meV = 1.62 Â T onset /K.In the present case this yields E a,MN E 0.3 eV.However, considering the analysis of Faske et al., 51 in the case of a distribution of lithium jump rates, this approach is at risk to underestimate the activation energy obtained and to overestimate the jump rate deduced from MN, respectively.In particular, this is the case when motion-induced line narrowing spans a broad temperature range, here this is more than 100 K, and T onset is mainly governed by the fraction of fast ions which are presumably subject to, e.g., short-range (or spatially localized) jump processes. 52n the present case, a more precise access to Li jump rates is given by SLR NMR.Before discussing the corresponding results, it is useful to have a closer look on the features of the NMR spectra recorded in the rigid-lattice regime and in the limit of extreme narrowing, i.e., at the highest temperatures where the line width is mainly governed by inhomogeneities of the external magnetic field.Interestingly, in the latter regime a well-defined quadrupole powder pattern shows up.Such electric quadrupole intensities are due to the interaction of the quadrupole moment of the 3/2-spin nucleus, which is a result of the non-spherical charge distribution of the nucleus, being immersed in an electronic environment characterized by a non-vanishing electric field gradient (EFG).The latter is produced by the electric charge distribution in the direct neighbourhood of the nucleus under observation.This additional interaction alters the Zeeman levels and thus the associated (angular) Zeeman frequency o 0 towards o 0 AE o Q . 54The corresponding quadrupole frequency is given by with the quadrupole coupling constant d = e 2 qQ/h.Here, e and eq are the proton charge and the principle component of the electric field gradient tensor, respectively.Q is the electric quadrupole moment of the 7 Li nucleus and h denotes Planck's constant.The angles Y and F specify the orientation of the external field B 0 in the principle axis system of the EFG tensor at the site of the 7 Li spin.The parameter Z Q reflects any asymmetries of the EFG tensor.Here, the 7 Li NMR line recorded at 417 K can be well represented by a quadrupole powder pattern with d = 14.6 kHz and Z Q = 0.6.Certainly, the values should be regarded as mean values characterizing the EFG tensor the Li spins are exposed to due to rapid hopping at elevated temperatures.In the low-temperature modification of Li 7 PSe 6 the Li ions occupy electrically (and magnetically) inequivalent crystallographic sites.In general, the absence of a complete averaging of electric quadrupole interactions at high T has been found for many Li ion conductors and might be regarded as a universal feature. 47,55,56It is referred to as anisotropic (or inhomogeneous) diffusion which is certainly traced back to the complex crystal structure of the lithiumcontaining argyrodite under investigation, cf. also Fig. 7 in ref. 16.
In the absence of lithium motions on the kHz scale, quadrupole intensities in the rigid-lattice regime can be made visible by recording NMR solid echoes.Owing to receiver dead-time effects, fast decaying components of an FID are difficult to detect.In Fig. 4 an NMR spectrum recorded in the rigid-lattice regime is shown which has been acquired at 190 K with the help of a two-pulse solid-echo sequence.Besides the Gaussian shaped central line, whose area fraction is about 36% of the total area, a broad quadrupole foot (64%) becomes apparent.The width (fwhm) of the latter component is about 54 kHz and points to a mean rigid-lattice quadrupole coupling constant in the order of 100 kHz.This value is comparable with those of other Li ion conductors studied.With increasing temperature, the rigid-lattice quadrupole interactions are partly averaged due to Li motions with rates as high as 10 5 s À1 .Finally, the well-defined powder pattern, as characterized above, shows up at higher T.
][27][28]37,57,58 In Fig. 5 the R 1r rates measured in both the laboratory and rotating frames of reference are shown in an Arrhenius plot.The R 1 rates reveal a  diffusion-induced low-temperature flank pointing to a relatively small activation energy of only 0.28 (1) eV.The so-called low-T regime of NMR SLR is sensitive to short-range rather than to long-range lithium diffusion. 20Interestingly, the SLR rates probed in the laboratory frame of reference do not exhibit a R 1 (1/T) rate peak; above 300 K the rates are independent of temperature and amount to be about 10 s À1 .Note that at T o 235 K the rates are increasingly governed by non-diffusive and weakly temperature-dependent background effects which, most likely, can be ascribed to coupling of the Li spins with paramagnetic impurities.However, in the case of NMR SLR measurements in the rotating frame of reference performed at a locking frequency in the order of 45 kHz the R 1r rates pass through a well-defined maximum showing up at T max = 325 K.This allows the direct determination of the Li jump rate t À1 which can be estimated using the maximum condition o 1 t E 0.5.Note that this estimation is based on an exponential correlation function G(t) leading to symmetric NMR rate peaks (see below).Here, at T max the rate t À1 is in the order of 1.4 Â 10 5 s À1 .Compared to the result from NMR motional narrowing it becomes evident that the rate t MN À1 seems to overestimate the jump process because of the reasons discussed above.
A more detailed analysis is achieved by using an appropriate relaxation model 20,59 to evaluate the SLR NMR rates R 1r (1/T) measured.The dashed line shown in Fig. 5 represents a complete relaxation fit according to ( Note that the correlation rate t c À1 is in the order of the jump rate t À1 .In general, J 3D denotes the spectral density function characterizing electric and magnetic field fluctuations due to three-dimensional ion hopping.Here, the expression of the spectral density function, to which R 1r (o 1 , T) is proportional, is based on the relaxation model of Bloembergen, Purcell and Pound (BPP) introduced for 3D, isotropic diffusion. 60riginally, the exponent b is given by b = 2, which originates from the assumption of an underlying single exponential motional correlation function G(t).J 3D is obtained from G(t) by Fourier transformation.In the case of b = 2 this results in a symmetric (Lorentzian shaped) spectral density function J 3D . 20alues smaller than two (for b the range 1 o b r 2 is valid) as introduced by Ku ¨chler et al. 61 lead to asymmetric NMR SLR rate peaks whereby the slope of the low-temperature flank, that is the regime where o 1 t c 1 holds, is smaller than that of the corresponding high-temperature side characterized by the limit o 1 t { 1.In particular, asymmetric NMR SLR rate peaks have been observed for structurally disordered ion conductors.3][64] A pronounced heterogeneous dynamics, being closely related to structural disorder which can be represented by an irregular formed potential landscape, is expected to include both long-range dynamics as well as localized hopping processes.
In the present case, the R 1r (o 1 ,T) data can be best represented by a fit with b = 1.41 (2).This result certainly points to a non-exponential motional correlation function.It might directly  The upper figure shows the temperature dependence of g used to parameterize the corresponding magnetization transients.For comparison, the temperature behaviour of the R 2 SSR NMR rates and the decay rates deduced from stimulated echo NMR (see the inset of Fig. 7) are also shown.See the text for further explanation.
This journal is c the Owner Societies 2013 reflect a heterogeneous potential landscape the Li ions are exposed to.In Li 7 PSe 6 the Li ions occupy multiple, distinctly connected sites as revealed by high-resolution 6 Li NMR spectroscopy. 15A relatively low b value, as found here, is comparable to that recently reported for Li ion hopping in Li 7 La 3 Zr 2 O 12 (LLZ); LLZ also exhibits a highly complex crystallographic structure where the Li ions are distributed among various inequivalent positions.Let us mention that b o 2 minimally influences the maximum condition used above to estimate the jump rate at the temperature where the rate peak shows up.Taking into account b = 1.41(2), the exact maximum condition is given by it turned out that D sd is 1.3(8) Â 10 À15 m 2 s À1 at T = 300 K. Here, the jump distance a was estimated to be a E 2.9 Å which is based on the structural data given in ref. 15.With D sd the corresponding ion conductivity can also be determined via the Nernst-Einstein equation.According to this relation the so-called solid-state diffusion coefficient D NE is obtained which is given by N denotes the charge carrier density and q represents the charge of the Li ions.s dc is the ionic dc conductivity.Note that D NE is directly linked with D sd according to Here, H R is the Haven ratio and f denotes the so-called correlation factor.In the present case we assumed f/H R E 1 to estimate the ionic conductivity s dc .Once again, N has been estimated based on the crystallographic data recently obtained for Li 7 PSe 6 . 15ith D sd given above, the ionic conductivity turns out to be in the order of 2 Â 10 À6 S cm À1 at 300 K.This is in reasonably good agreement with preliminary conductivity measurements on similar samples. 16It should be noted that even higher ion conductivities have been probed for some halide-containing Li-argyrodites. 9,16,17dentifying E a = 0.49(2) eV (see above) with the activation energy reflecting long-range Li diffusion in Li 7 PSe 6 , the slope of the low-T flank of the rate peak of Fig. Indeed, the crystallographic structure of Li-containing argyrodites contains structural motifs, viz.Li ions residing in close-by tetrahedral voids, which might help explain the differences between localized, i.e., short-range Li dynamics on the one hand and long-range Li diffusion on the other hand. 15This is also supported by recent theoretical investigations on halidecontaining Li-argyrodites. 16,17or comparison with the SLR NMR rates, the g exponents describing the deviation of the magnetization transients M r (t lock ), see Fig. 6, from single-exponential time behaviour are also included in Fig. 5. Whereas a straightforward interpretation of g is generally not available, one might read out two features showing up in the present case.At very low temperatures the exponent is approximately 0.5 which might point to (diffusion-induced) spin-lattice relaxation due to coupling of the 7 Li spins primarily to paramagnetic impurities.At even lower temperatures homonuclear spin-diffusion might play an additional role.With increasing diffusivity, g initially increases before a minimum value is attained at T max where the corresponding maximum peak R 1r (1/T) shows up.Interestingly, when T is further increased and Li diffusion becomes rapid the exponent g steadily increases towards g = 1, i.e., singleexponential time behaviour.Note that at very high temperatures the various electrical interactions the Li spins are exposed to are greatly averaged resulting in a quadrupole powder pattern which can be simulated with a single set of (mean) parameters.Thus, when considering the underlying motional correlation function, being at least partly reflected by the exponent g, it might also change as a function of temperature and seems to tend to become a single exponential in the limit o 1 t { 1.Such an observation is in agreement with the predictions of Ngai and co-workers, 66 who have analyzed the shape of motional correlation functions from both NMR and impedance spectroscopy measurements.Dashed and solid lines represent fits according to M r (t lock ) p exp(À(t lock /T 1 ) g ).At high temperatures, i.e. in the limit o 1 t c 1, the stretching exponent g increases towards single exponential time behaviour.The decay curves recorded above For the sake of completeness, R 2 SSR NMR rates were recorded at temperatures down to 190 K.In the rigid-lattice regime, T 2 is approximately 70 ms.This is in good agreement with the second moment, i.e., the line width of the central transition in this temperature range (6.8 kHz), see above.It is worth mentioning that, as expected, the 7 Li NMR solid echoes are composed of a fast and slowly decaying contribution which mirror quadrupolar and dipolar contributions.Irrespective of what is analyzed in detail, i.e., (i) the complete area under the echo, (ii) the amplitude of the quadrupole part or (iii) the area of the slowly decaying dipolar part, at very low T the corresponding transients follow single-exponential time behaviour (see Fig. S1 and S2, ESI †).However, the situation changes upon increasing T.Then, a fast decaying and temperature independent contribution with R 2 0 = 30 ms shows up (see Fig. S2, ESI †).These rates might reflect a very slow dynamic process.Interestingly, the stimulated echoes, which are presented below, also seem to point to a subset of lithium ions with a very low diffusivity.
For comparison, in Fig. 5 the temperature-dependent R 2 rates of the slowly decaying part of the transients are included.The rates start to decrease at T = 260 K.At this temperature the associated jump rate is in the order of 10 4 s À1 , which is in good agreement with the Arrhenius parameters deduced from the R 1r (1/T) peak.However, the corresponding activation energy, which can be estimated from the linearly decaying part of the R 2 curve, turns out to be only 0.4 eV (see also Fig. S2, ESI †).Moreover, the maximum R 1r rate exceeds R 2 .Similarly, this has been observed by, e.g., Bertermann et al. 56 who recorded R 2 rates using a Hahn echo pulse sequence.In doing so, Hahn echo decay seems to be determined by dipolar rather than by quadrupolar interactions as expected.
The present difference between R 1r and R 2 also indicates that the two rates probed are not necessarily influenced by exactly the same relaxation mechanism and/or the same subset of (differently coupled) Li spins.Considering the pre-factor obtained from the R 1r (1/T) fit (R 1r = C 0 Â J 3D (2o 1 )), which amounts to C 0 = 2.95(9) Â 10 9 s À2 , thus, being in the order of the square of the rigid-lattice quadrupole coupling constant (see above), we suppose that, besides dipolar ones, fluctuating quadrupolar interactions also play a large role in influencing R 1r .However, a conclusive answer to what extent dipolar interactions affect the SLR NMR rates probed might only be possible by future 6 Li NMR experiments carried out at resonance frequencies comparable to those used for the 7 Li NMR measurements.
In addition to NMR SLR measurements in the rotating frame of reference, we have recorded 7 Li spin-alignment echo NMR decay curves which are sensitive to slow Li exchange processes between electrically inequivalent Li sites in the argyrodite.Thus, the method takes advantage of the distinct quadrupolar interactions of the jumping Li ions.The decay curves shown in Fig. 7 were obtained by plotting the intensity of the spinalignment echo (read out at the echo maximum, i.e., at t = t p ) as a function of mixing time.The S 2 curves reflect the probability to find an Li nuclear spin, which initially (t m = 0) resides on a specific site, on an electrically equivalent site at a later time (t = t m ).
Interestingly, echo damping proceeds in two steps. 22,26,67ost likely, the first decay step directly reflects Li hopping; the second step, however, might be induced by (quadrupolar) spinlattice relaxation and/or spin-diffusion effects.The latter are expected to be characterized by a decay rate t 00 being in the order of seconds.The full decay curves were parameterized with a combination of two stretched exponentials according to Below 270 K, the stretching parameter g SAE is approximately 0.5.It decreases with increasing T reaching g SAE E 0.25 at ambient temperature.The fits yield the decay rates t SAE À1 shown in the inset of Fig. 7. Obviously, the first decay step points to Li jumps which are thermally activated similar to those being responsible for motion-induced SLR NMR in the limit o 0 t c 1 (see above).9][70] In the examples studied so far, the corresponding activation energy E a,SAE turned out to be very similar to the value E a,dc deduced from dc-conductivity measurements.However, in the present case, low-temperature echo damping seems to be mainly influenced by a diffusion process with a much lower thermal activation.This behaviour might change when 6 Li (spin-1 nucleus) stimulated echoes are recorded instead of 7 Li (spin-3/2 nucleus) echoes.This is because the latter are exposed to rather large homonuclear dipole-dipole interactions which might strongly influence the generation of a pure spin-alignment state.
Irrespective of this possible influence on the 7 Li NMR S 2 curves, even at temperatures as low as 220 K the two-step decay behaviour remains.For example, the (normalized) correlation Fig. 7 7 Li NMR (sin-sin) two-time correlation functions recorded using the function recorded at 241 K is characterized by a so-called plateau value S N = B/(A + B) (see eqn (6)) which is approximately 0.12.In general, 1/S N is related to the number of electrically inequivalent Li sites the ion visits during self diffusion. 39The inverse, 1/S N , equals the number of available Li sites provided these are equally populated and visited by the spins with the same probability.S N = 0.12 leads, at first glance, to approximately eight different Li sites.In fact, there are at least seven electrically inequivalent Li sites occupied in Li 7 PSe 6 .It is worth mentioning that, once again, 7 Li- 7 Li dipole-dipole interactions might influence 1/S N which usually point to a larger number of sites than actually present in the crystal structure.Such a situation has also been found for SAE NMR measurements on b-eucryptite, LiAlSiO 4 , single crystals. 71With increasing temperature T, the value S N decreases which might be attributed to the increasing averaging of electric quadrupolar as well as magnetic dipolar interactions due to sufficiently fast Li hopping.
Alternatively, the emergence of a distinct plateau might also be interpreted in terms of a small and decoupled sub-ensemble of very slow Li ions that show up at large mixing times and low temperature.However, since the corresponding decay rates associated with the second step (1/t 00 ) do not reveal a significant temperature dependence, the above-given explanation, pointing to the averaging of quadrupole interactions, seems to be more suitable to understand the NMR response.Moreover, at higher temperatures the jump process occurring during the evolution period determined by t p might also affect the correlation functions probed.This phenomenon is called motional phase averaging, see ref. 21 for further details.
In contrast to the diffusion process being responsible for the maximum of R 1r (1/T), the first decay step of the SAE NMR decay curves reflects Li spins diffusing on a rather long time scale.In this context it is worth mentioning once again that, in general, SAE NMR is insensitive to Li + jumps between electrically inequivalent sites.This also holds for sufficiently fast exchange processes for which echo formation is not possible or for hopping processes leading to elimination of fluctuating EFGs.The latter means fast jump processes involving a small number of preferred but electrically inequivalent sites.Such a process would lead to an averaged but mixing-time independent EFG and is expected to contribute to SLR NMR but it would be invisible for SAE NMR.This view might help understand the difference found between t À1 from R 1r (1/T) and t SAE

À1
. Such a scenario of fast and slowly diffusing spins, 16 which allow exchange between each other, is also in agreement with a relatively large final state amplitude for the slow process most likely involving the majority of the Li sites available.As mentioned above, exchange between closely adjacent voids in Li 7 PSe 6 , see ref. 15 and 16, might represent such a preferred (spatially localized and fast) hopping process to which SAE NMR is not sensitive, cf. also the illustration of Fig. 7 in ref. 16.The relatively low activation energy deduced from SAE NMR might point to the conclusion that at higher T further hopping processes come into play which are finally detected by R 1r (1/T) in the high-temperature limit.
For comparison with a recent study on Li 4 SnS 4 , let us note that Dehnen and co-workers 72 reported on SAE data which reveal low-temperature echo decay damping pointing to an extremely slow Li motional process.This process is reported to be thermally activated by 0.05 eV only.Very recently, by our own group a similar motion process with an activation energy of 0.09 eV was observed for a structurally disordered ion conductor.In contrast to the SAE NMR response probed here, but in analogy to the well-known nearly constant loss (NCL) phenomenon in conductivity spectroscopy, in ref. 73 we interpreted this feature as a result from strictly localized (or caged) motions, see also ref. 74.
Finally, instead of a single (Arrhenius-type) diffusion process we found indications for various motional mechanisms present in Li 7 PSe 6 .This also includes the possibility that diffusion mechanisms change as a function of temperature and dominate others in certain ranges.It also helps understand why it is hardly possible to perform a so-called joint fit, see, e.g., ref. 28, of the temperature-variable and frequency-dependent SLR NMR rates R 1r and R 1 .Here, the rates R 1 might be interpreted as the result of a superposition of different processes leading to 0.28 eV (Fig. 5) in the end.In contrast, long-range Li + ion dynamics at higher T seem to be determined by an activation energy of 0.49 eV which can be deduced from the R 1r (1/T) peak in the relaxation limit o 1 t c { 1.

IV. Conclusion and outlook
In summary, local electric structures and Li self-diffusion in a polycrystalline sample of Li 7 PSe 6 have been examined by temperature-variable 7 Li NMR spectroscopy.Line shape as well as SLR measurements have been used to get insights into translational Li dynamics present.While the first seem to be mainly governed by local exchange processes, with the help of SLR measurements short-as well as long-range diffusion parameters are accessible.In particular, it turned out that SLR measurements performed in the rotating frame of reference are highly suited to characterize Li hopping motions taking place on different length scales in the Li-containing argyrodite.From the diffusion-induced R 1r SLR rate peak probed, both the activation energy (0.49(2) eV) as well as the pre-exponential factor (1.6(9) Â 10 13 s À1 ) of the underlying Arrhenius relation characterizing long-range dynamics (at higher T) could be probed.At ambient temperature, the Li self-diffusion coefficient is approximately 1.3(8) Â 10 À15 m 2 s À1 .This points to an Li ion conductivity in the order of 2 Â 10 À6 S cm À1 at 300 K.In contrast, at lower T Li dynamics seem to be characterized by activation energies ranging from 0.21 eV (R 1r ) to 0.28 eV (R 1 ).Most likely, the relevant diffusion processes being valid in this T range also entail short-range dynamics.
9]75 Currently, chlorine-containing and bromine-containing argyrodites are investigated in our laboratory.In addition to the first 7 Li SAE NMR measurements presented here, the analysis of This journal is c the Owner Societies 2013 Phys.Chem.Chem.Phys., 2013, 15, 7123--7132 7131 6 Li NMR hopping correlation functions might be very useful to shed more light on the dynamic processes probed via the generation of stimulated echoes.

Fig. 1 (
Fig.1 (A)The saturation recovery pulse sequence used to probe the SLR NMR rates R 1 : a comb of closely spaced rf-pulses (p/2) x 0 (n E 10) destroys the (initial) longitudinal magnetization (M z = M eq ) so that M z = 0 at t = 0.The subsequent recovery of M z (t) is probed with a single p/2-pulse which is sent after a variable relaxation delay t = t d .(B) The pulse sequence of the spin-lock technique to record the rotating-frame NMR rates R 1r : immediately after the p/2-pulse the magnetization M pointing along the (Ày 0 ) is locked by the field B 1 .The decay of M r(Ày 0 ) is then probed by plotting the height (or area) of the free induction decay (FID) as a function of the locking pulse length t lock .See also ref.37.In a simplified manner (C) and (D) depict the pulse sequences to record SSR NMR rates and stimulated echo decay curves.See the text for further explanation.

Fig. 2
Fig.2(a) and (b)7 Li NMR spectra of ternary Li 7 PSe 6 recorded at the temperatures indicated.The line width steadily decreases reflecting a pronounced motional narrowing which indicates rapid Li exchange on the NMR time scale defined by the rigid lattice line width of Li 7 PSe 6 .Pronounced and well-defined quadrupole powder patterns emerge at temperatures chosen to be higher than 300 K.It is fully developed at 417 K. (c) Magnification of the quadrupole powder pattern of the7 Li NMR recorded at 417 K.The spectrum shown at the bottom displays a simulation of the measured one with the help of WinSolids software 49 using a single set of fitting parameters describing the NMR line shape.See the text for further details.

Fig. 3
Fig.3Motional narrowing of the7 Li NMR linewidth of polycrystalline Li 7 PSe 6 .While filled symbols (K) represent data points recorded using a standard probe, boxes display the line widths (fwhm) of NMR spectra recorded using a cryo probe.The rigid lattice value, n rl E 6.8 kHz, is reached at temperatures below 150 K.The onset of MN is expected when the Li jump rate reaches values in the order of 10 3 s À1 .See the text for further explanation.

Fig. 4 7
Fig. 4 7 Li solid-echo (first-order quadrupolar) NMR spectrum recorded using a two-pulse sequence at 116 MHz and 190 K.The area fractions of the dipolar broadened central and quadrupole contribution are in agreement with the theoretical ratio of 4 : 6. Upon increasing temperature, the central line and the broad quadrupole signal narrow due to rapid Li exchange.

Fig. 5
Fig. 5 Arrhenius plot of the 7 Li NMR SLR rates measured at a Larmor frequency o 0 /2p of 116 MHz (R 1 ) and a locking frequency o 1 /2p of 45(2) kHz (R 1r ).The dashed line represents a modified BPP fit taking into account a sub-quadratic frequency dependence of the R 1r rates (R 1r p o 1 Àb , b o 2) in the limit o 1 t c 1.
5 might be related to Li jump processes taking place on a much shorter length scale.This is expressed by o 1 t c 1 being valid in this temperature range.In the present case, the low-T flank is characterized by E a,low = 0.21(1) eV, thus, fulfilling the relation E a,low = (b À 1)E a .It directly links the frequency dependence of R 1r p o 1 Àb to the asymmetry of the diffusion-induced relaxation rate peak R 1r (1/T).

Fig. 6 7
Fig. 6 7 Li NMR magnetization transients M r (t lock ) of Li 7 PSe 6 recorded at the temperatures indicated and at o 0 /2p = 116 MHz and o 1 /2p = n 1 = 45(2) kHz.Dashed and solid lines represent fits according to M r (t lock ) p exp(À(t lock /T 1 ) g ).At high temperatures, i.e. in the limit o 1 t c 1, the stretching exponent g increases towards single exponential time behaviour.The decay curves recorded above 360 K represent magnetization transients of the high-temperature flank of the R 1r (1/T) peak shown in Fig.5.
B denotes Boltzmann's constant.At ambient temperature, this yields t À1 E 9.5 Â 10 4 s À1 .Using the Einstein-Smoluchowski equation, which relates t À1 with the self-diffusion coefficient D sd via