Abstract
The diffusive behavior in a spinel-type Li+ ion battery material, Li[Ni1/2Mn3/2]O4, has been studied with positive and negative muon spin rotation and relaxation (μ ±SR) measurements in the temperature range between 200 and 400 K using a powder sample. The implanted μ + locates at an interstitial site near O2− ion so as to form a O–H like bond, while the implanted μ − is mainly captured by an oxygen nucleus, resulting in the formation of muonic oxygen. This means that local magnetic environments in Li[Ni1/2Mn3/2]O4 were investigated from the two different sites in the lattice, i.e., one is an interstitial site for μ +SR and the other is an oxygen site for μ −SR. Since both μ +SR and μ −SR detected an increase in the fluctuation rate of a nuclear magnetic field for temperatures above 200 K, the origin of this increase is clearly confirmed as Li diffusion. Assuming a random walk process with the hopping of thermally activated Li+ between a regular Li site and the nearest neighboring vacant octahedral sites, a self-diffusion coefficient of Li+ was found to range above 10−11 cm2/s at temperatures above 250 K with an activation energy of about 0.06 eV.
1 Introduction
In order to understand the operational processes of battery materials, it is crucial to understand the diffusive nature of ions in these materials, because the transport of ions is the fundamental mechanism behind their operation. Nuclear magnetic resonance (NMR) has been widely used for such a purpose [1], [2], [3], through the observation of nuclear magnetic fields in these materials. Recently, the positive muon spin rotation and relaxation (μ +SR) technique [4, 5] was found to be a good tool for detecting ion diffusion in solids [6]. Unlike in NMR, where the polarization is dependent on an externally applied high magnetic field, the initial muon spin polarization is aligned antiparallel to its momentum due to parity violation of a muon creation reaction [7]. As a result, μ +SR senses the internal magnetic field in materials without the need of any external magnetic field, i.e., even in zero field (ZF).
This leads to the following two important features: the first one is that μ +SR distinguishes the nuclear magnetic field from the electron magnetic field by simply applying a weak longitudinal field (LF) [8]. Here, “weak” means the field is comparable to the field distribution width formed by surrounding nuclear magnetic moments at the muon site, which typically ranges between 5 and 10 Oe, and “longitudinal field” means the magnetic field is parallel to the initial muon spin polarization. Such weak LFs decouple the contribution from a nuclear magnetic field on the μ + spin, while the effect of weak LF on the contribution from the electron magnetic field is negligibly small. Consequently, μ +SR provides information on ion diffusion even in materials that include paramagnetic ions [9, 10]. This is an attractive feature for studying Li+ [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], Na+ [10, 22, 23], and K+ [24] diffusion in cathode materials for an ion battery, because these materials always contain transition metal ions to keep charge neutrality during ion insertion and extraction reactions.
Another feature that makes μ +SR an attractive probe is its sensitivity, which in principle is temperature independent, since the muon spin polarization is an intrinsic feature of the technique, whereas the effect of an external field on aligning nuclear spins is suppressed by thermal excitation, resulting in the decrease of NMR sensitivity with temperature. Therefore, μ +SR is a very important tool to study ion dynamics at high temperatures.
Despite the above features, the information obtained with μ +SR is sometimes ambiguous. Whereas the implanted μ + is presumed to be static at low temperatures in its interstitial site (where its electrostatic potential is in a local minimum), in any solid there is a chance that the μ + will start to diffuse above certain temperatures. This means that the observed diffusive behavior could be caused by not only ion diffusion but also μ + diffusion. In order to confirm the diffusing species, μ + or ion, we have initiated series of μ −SR experiments on battery materials [25], [26], [27]. Here, μ −SR is the “anti-technique” of μ +SR by not using μ + but negative muons (μ −) [28], [29], [30], [31].
The implanted μ − is captured by a nucleus and forms a stable muonic atom. The muonic atoms, apart from the ones based on the diffusing species, is immobile until the decomposition temperature of the material. In other words, by means of μ ±SR, we are able to see the internal magnetic field from different viewpoints: at the atomic position in the lattice with μ −SR and at the interstitial sites with μ +SR. It should be noted that a μ −SR spectrum requires about 36 times higher counts than that of μ +SR to obtain the same statistical precision, since the spin polarization of μ − is down to about 1/6 of its initial polarization during the atomic capture process by a nucleus [32,33]. It is thus only reasonable to use the expensive and difficult μ −SR technique in order to determine the diffusing species, while the details of the diffusion processes are best determined by μ +SR.
A combined μ
±SR work has confirmed Li+ diffusion in an olivine-type LiMnPO4 above around 200 K [26] and in a spinel-type Li4Ti5O12
![Figure 1:
The crystal structure of cubic spinel Li[Ni1/2Mn3/2]O4 drawn with VESTA [41]. Li+ ions occupy the tetrahedral (8a) site, while Ni and Mn ions occupy the octahedral (16d) site. O2− ion is located in each corner of the tetrahedron and octahedron. Dark grey spheres represent a vacant 16c site.](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_001.jpg)
The crystal structure of cubic spinel Li[Ni1/2Mn3/2]O4 drawn with VESTA [41]. Li+ ions occupy the tetrahedral (8a) site, while Ni and Mn ions occupy the octahedral (16d) site. O2− ion is located in each corner of the tetrahedron and octahedron. Dark grey spheres represent a vacant 16c site.
2 Experimental
A powder sample of Li[Ni1/2Mn3/2]O4 was purchased from Toshima Manufacturing, Ltd. The μ −SR time spectra were measured on the decay muon beam line ARGUS [42, 43] at ISIS pulsed neutron and muon source of Rutherford Appleton Laboratory in UK. Approximately 22 g powder sample was placed in a copper container with a sample space of 30 mm diameter and 15 mm depth, made of 1 mm thick Cu plate [see Figure 2(a)]. A cap of the container was made of 3 mm thick Cu plate, which also works as a degrader to adjust the range of μ −s. The copper container was then set onto the end of the He-flow cryostat (Ministat, Oxford Instruments). The momentum of the μ − beam was adjusted to 70 MeV/c to maximize the number of μ − stopped in the sample, as shown in Figure 2(b). The μ −SR spectrum was measured at 200, 300, and 400 K with up to 50 Mevents for transverse field (TF) μ −SR and 200 Mevents for zero field (ZF) and 150 Mevents for longitudinal field (LF) μ −SR at a counting rate of 17 Mevents/h. Here TF [LF] means the applied magnetic field was perpendicular [parallel] to the initial μ − spin polarization.
![Figure 2:
(a) A copper container filled with a black powder sample of Li[Ni1/2Mn3/2]O4 for a μ
−SR experiment. The size of the sample space is 30 mm diameter and 15 mm depth. (b) The variation of the depth profile of the muons in the copper container filled with the sample by five different implantation momentum. The profile was simulated with SRIM [44].](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_002.jpg)
(a) A copper container filled with a black powder sample of Li[Ni1/2Mn3/2]O4 for a μ −SR experiment. The size of the sample space is 30 mm diameter and 15 mm depth. (b) The variation of the depth profile of the muons in the copper container filled with the sample by five different implantation momentum. The profile was simulated with SRIM [44].
The μ +SR time spectra were measured on the surface muon beam line EMU [43, 45] at ISIS. Approximately 2 g powder sample from the same batch as for the μ −SR measurements was packed into a standard ISIS titanium cell. The window of the cell was made of a titanium foil with 20 μm thickness. The cell was then mounted onto the Al plate of a closed cycle refrigerator-type cryo-furnace in the temperature range between 50 and 500 K. The counting rate was about 100 Mevents/h. The experimental techniques are described in more detail elsewhere [4, 5]. The obtained μ ±SR data was analyzed with musrfit [46].
3 Results
While the μ −SR study is the main focus of this paper, to put these results into context our discussion will start with the μ +SR measurements that were carried out on the same system.
3.1 μ +SR
Figure 3 shows the μ +SR time spectra recorded in wTF with H = 20 Oe at selected temperatures in order to understand the magnetic nature of the Li[Ni1/2Mn3/2]O4 sample. At 150 K, the spectrum exhibits a clear oscillation with a full asymmetry due to wTF, which means that the whole volume of the sample is in a paramagnetic state. As the temperature decreases from 150 K, the oscillating asymmetry decreases with decreasing temperature and then almost disappears below 100 K. In order to extract the μ +SR parameters, the wTF-μ +SR spectrum is fitted by a combination of an exponentially relaxing cosine oscillation due to wTF and an exponentially relaxing non-oscillatory signal, due to the 1/3 tail arising from the distribution of static electronic components in the powder sample when in an AF ordered state.
where A 0 is the initial asymmetry determined by the wTF measurements at high temperatures and P TF is the muon spin polarization function in wTF. A TF, f TF, ϕ TF and λ TF are the asymmetry, muon spin precession frequency, initial phase, and relaxation rate due to the applied TF. Here, f TF = γ μ /2π × 20 Oe = 13.554 kHz/Oe × 20 Oe ∼ 0.27 MHz. The A tail component corresponds to the signal from the magnetic ordered phase, in which the internal AF field is parallel to the initial muon spin polarization.
![Figure 3:
Weak transverse field (wTF) μ
+SR time spectra for Li[Ni1/2Mn3/2]O4 recorded at 50, 100, 110, 120, and 150 K. Solid lines represent the best fit using Eq. (1).](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_003.jpg)
Weak transverse field (wTF) μ +SR time spectra for Li[Ni1/2Mn3/2]O4 recorded at 50, 100, 110, 120, and 150 K. Solid lines represent the best fit using Eq. (1).
Figure 4 shows the temperature dependencies of the four wTF-μ
+SR parameters, i.e., A
TF, A
tail, λ
TF, and λ
tail. The A
TF(T) curve exhibits a step-like change at around 105 K, which indicates the appearance of a very large internal magnetic field compared with wTF at low temperatures. From the middle point of the step-like change in the A
TF(T) curve, T
N is determined as 104.2(5) K using a Sigmoid function. A relatively wide transition (ΔT
N = 41 K) indicates a random distribution of Ni and Mn at the octahedral site in the spinel lattice [40,47], where ΔT
N = T
90% − T
10% with
![Figure 4:
The temperature dependencies of (a) the weak transverse field asymmetry (A
TF) and tail asymmetry (A
tail) (b) their exponential relaxation rates (λ
TF and λ
tail) for Li[Ni1/2Mn3/2]O4. The data were obtained by fitting the wTF-μ
+SR spectrum with Eq. (1). A red solid line in (a) represents the best fit by a Sigmoid function and a blue dotted line indicates the predicted
A
tail
=
[
A
TF
(
200
K
)
−
A
TF
(
30
K
)
]
/
3
${A}_{\text{tail}}=[{A}_{\text{TF}}(200\text{\hspace{0.17em}K})-{A}_{\text{TF}}(30\text{\hspace{0.17em}K})]/3$
. The color of the symbol for the TF signal is slightly changed at temperatures below 90 K, since the wTF-μ
+SR signal comes from not the sample but from the Ti sample holder (see text).](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_004.jpg)
The temperature dependencies of (a) the weak transverse field asymmetry (A
TF) and tail asymmetry (A
tail) (b) their exponential relaxation rates (λ
TF and λ
tail) for Li[Ni1/2Mn3/2]O4. The data were obtained by fitting the wTF-μ
+SR spectrum with Eq. (1). A red solid line in (a) represents the best fit by a Sigmoid function and a blue dotted line indicates the predicted
As temperature decreases from 200 K, A
tail starts to increase below 130 K and approachs a predicted value for the 1/3 tail signal below 90 K, where A
tail well below T
N is predicted as
Figure 5 shows the ZF- and two LF-μ +SR spectra recorded at 150, 300, and 425 K, i.e., in a paramagnetic state. As temperature increases, the μ +SR spectrum indicates a clear change from a low-temperature static behavior to a high-temperature dynamic behavior. The μ +SR spectra were therefore fitted with a combination of an exponentially relaxing dynamic Kubo-Toyabe function (G DGKT) from the sample and a time independent background signal from the Ti sample holder:
where P LF is the muon spin polarization function in ZF and LF. A KT and A BG are asymmetries of the two signals, Δ is the field distribution width, ν is the field fluctuation rate, and λ KT is the exponential relaxation rate. For ν = 0 and H LF = 0 (i.e., ZF), G DGKT(t, Δ, ν, H LF) becomes a static Gaussian Kubo-Toyabe function (G zz KT) given by [48]:
which represents the time variation of the muon spin polarization due to the internal magnetic field formed by randomly oriented static nuclear dipoles with Gaussian distribution. Since the implanted muons feel both a nuclear magnetic field formed mainly by 6Li, 7Li, and 55Mn and a magnetic field caused by localized 3d electrons of Ni2+ (3d 8) and Mn4+ (3d 3) [36, 39, 49], an exponential is multiplied with G DGKT [50] in Eq. (2).

The temperature variation of the ZF- and two LF-μ +SR spectra recorded at (a) 150 K, (b) 300 K, and (c) 425 K. The applied LF’s were 5 and 10 Oe. Green solid lines represent the best fit using Eq. (2).
In the paramagnetic state, the TF-μ
+SR spectrum is also fitted with an Abragam equation [51, 52] using Δ+ and ν
+ and exponential relaxation function
Combining with Eq. (2), the TF-, ZF-, and LF-μ
+SR spectra were fitted using common A
KT, A
BG, and A
TF at temperatures between 150 and 425 K, while Δ+, ν
+, and
Figure 6 shows the temperature dependencies of the μ +SR parameters obtained by fitting the TF-, ZF-, and two LF-μ +SR spectra with Eqs. (4) and (8). Here, a superscript of each symbol, i.e., “+” means that the data obtained from μ +SR measurements. The field distribution width (Δ+), which roughly corresponds to the spin-spin relaxation rate (1/T 2), is found to be almost temperature independent.
![Figure 6:
The temperature dependencies of the μ
+SR parameters for Li[Ni1/2Mn3/2]O4: (a) the field distribution width (Δ+), (b) the field fluctuation rate (ν
+) and the exponential relaxation rate (
λ
KT
+
${\lambda }_{\text{KT}}^{+}$
), and (c) the three asymmetries (
A
TF
+
${A}_{\text{TF}}^{+}$
,
A
KT
+
${A}_{\text{KT}}^{+}$
, and
A
BG
+
${A}_{\text{BG}}^{+}$
). “+” of each symbol means that the data was obtained from μ
+SR measurements. Δ+, ν
+,
A
KT
+
${A}_{\text{KT}}^{+}$
, and
A
BG
+
${A}_{\text{BG}}^{+}$
were obtained by fitting the TF-, ZF- and LF-μ
+SR spectra with Eq. (2). Note that Δ+ corresponds roughly to a spin-spin relaxation rate (1/T
2), while ν
+ corresponds roughly to a spin-lattice relaxation rate (1/T
1). More correctly, 1/T
1 would be the relaxation rate of the nuclear moments that contribute to Δ+, not the relaxation rate of the μ
+ spin.](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_006.jpg)
The temperature dependencies of the μ
+SR parameters for Li[Ni1/2Mn3/2]O4: (a) the field distribution width (Δ+), (b) the field fluctuation rate (ν
+) and the exponential relaxation rate (
The field fluctuation rate, ν
+, which corresponds roughly to the spin-lattice relaxation rate, increases almost monotonically with increasing temperature, being consistent with Li diffusion. Further confirmation for Li diffusion is given in Section 3.2. On the contrary,
In the above fit, each asymmetry is assumed to be temperature independent, because the fractions of muons stopped in the sample and the Ti cell are, in principle, independent of temperature.
3.2 μ −SR
When considering the positive muon results, we cannot rule out the possibility that the muon may start to diffuse in the sample at some temperature. To be certain we are studying the change in the nuclear magnetic field from a “fixed view point” we have also measured the μ −SR spectrum for Li[Ni1/2Mn3/2]O4 using the same sample to that for μ +SR. Figure 7 shows the time histogram of the TF-μ −SR spectrum with TF = 50 Oe recorded at 300 K. Here, forward [backward] means upstream [downstream] detector positions with respect to the sample position and the incoming muon beam direction. Since the decay asymmetry is very small and the lifetime of the μ − depends upon the nucleus on which it is captured by, the histogram of μ −SR was fitted by a combination of four different lifetimes (τ i=1−4). Note that, for μ +SR, the time histogram includes only one component, i.e., μ + with τ = 2.19703 µs:
where N
i
is a normalization constant at t = 0 for the i-th muon lifetime (τ
i
), i.e., muonic O (μ
−O) with τ
1 = τ
O = 1.7954 μs [53], muonic Li (μ
−Li) with τ
2 = τ
Li = 2.188 μs [53], muonic Mn (μ
−Mn) with τ
3 = τ
Mn = 0.2325 μs [53], and an unknown component with τ
4 = 26.8(6) μs, A
i
is the muon decay asymmetry for each process, ω
i
= 2πf
i
is the angular frequency of the μ
− spin precession caused by the applied TF, and ϕ
i
is the initial phase. Since
![Figure 7:
The time histogram of the TF-μ
−SR spectrum in Li[Ni1/2Mn3/2]O4 for the (a) forward counter and (b) backward counter. Red open circles represent the experimental data, green solid lines represent the fit result using Eq. (5), and blue solid lines represent the histograms of the four muonic-atom lifetimes.](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_007.jpg)
The time histogram of the TF-μ −SR spectrum in Li[Ni1/2Mn3/2]O4 for the (a) forward counter and (b) backward counter. Red open circles represent the experimental data, green solid lines represent the fit result using Eq. (5), and blue solid lines represent the histograms of the four muonic-atom lifetimes.
As a first approximation, A
i
is non-zero only for the main component, i.e., μ
−O and negligible for the other elements. This is because both Li and Mn have nuclear spins and will therefore form hyperfine-coupled states with the μ
−, leading to a precession at a dramatically different frequencies compared with the applied value. Furthermore, A
i
is less than 1/6 of A
0 for μ
+SR (typically 0.24) due to the reduction of spin polarization,
The relative asymmetric counting efficiency can be estimated for each element and are listed in Table 1. It is noted that both
Parameters of the TF-μ −SR histograms for the Li[Ni1/2Mn3/2]O4 sample obtained by fitting with Eq. (5), where τ O, τ Li, and τ Mn were fixed to the values reported in Ref. [53]. Since it is extremely difficult to determine short τ components using data at early times because of the difficulty for determining a correct time zero (t 0) of the pulsed μ ± beam, the histograms were fitted only in the time domain between 0.4 and 16 μs.
n | Component | τ (μs) |
|
|
|
α i/α O |
---|---|---|---|---|---|---|
1 | O | τ O = 1.7954 [53] | 1 | 1 | α O = 1.062(5) | 1 |
2 | Li | τ Li = 2.188 [53] | 0.014(2) | 0.010(2) | α Li = 0.8(3) | 0.7(2) |
3 | Mn | τ Mn = 0.1569 [53] | 1.05(2) | 1.33(4) | α Mn = 1.34(6) | 1.26(6) |
4 | ? | τ 4 = 26.8(6) | 0.00110(2) | 0.00073(2) | α 4 = 0.71(3) | 0.67(3) |
In fact, the TF-, ZF-, and two LF-μ
−SR spectra for Li[Ni1/2Mn3/2]O4 exhibit a sudden decrease with time particularly below around 8 µs, regardless of external magnetic field (see Figure 8). This is well explained by the effect of the τ
4 component on the μ
−SR spectra, because the τ
4 component with
![Figure 8:
The TF-, ZF-, and two LF-μ
−SR spectra for Li[Ni1/2Mn3/2]O4 recorded at 300 K. The applied TF was 50 Oe, and the applied LF were 20 and 50 Oe.](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_008.jpg)
The TF-, ZF-, and two LF-μ −SR spectra for Li[Ni1/2Mn3/2]O4 recorded at 300 K. The applied TF was 50 Oe, and the applied LF were 20 and 50 Oe.
Therefore, as seen in Figure 9, the TF-, ZF-, and two LF-μ −SR spectra were fitted in the time domain between 0.8 and 6 μs with a combination of an exponentially relaxing cosine oscillation and an exponentially relaxing non-oscillatory signal for the TF-spectrum, and with a combination of an exponentially relaxing dynamic Kubo–Toyabe signal and an exponentially relaxing non-oscillatory signal for the ZF- and LF-spectra. In addition, the following three parameters, i.e., A TF, A TFe, and λ TFe, are common for the TF-, ZF-, and LF-μ −SR spectra at each temperature, and Δ and ν are common for the ZF- and LF-μ −SR spectra at each temperature:
where
![Figure 9:
The TF-, ZF-, and two LF-μ
−SR spectra for Li[Ni1/2Mn3/2]O4 recorded at (a) 200 K, (b) 300 K, and (c) 400 K in the time domain between 0 and 7 μs. The applied TF was 50 Oe, and the applied LF were 20 and 50 Oe. (b) is the same to Figure 8.](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_009.jpg)
The TF-, ZF-, and two LF-μ −SR spectra for Li[Ni1/2Mn3/2]O4 recorded at (a) 200 K, (b) 300 K, and (c) 400 K in the time domain between 0 and 7 μs. The applied TF was 50 Oe, and the applied LF were 20 and 50 Oe. (b) is the same to Figure 8.
The exponentially relaxing non-oscillatory signal (A TFe and A LFe), which is common for the TF-, ZF-, and LF-μ −SR spectra, is most likely originating from different lifetimes of μ − between μ −Li and μ −O [26], but both μ − spins are also depolarized by the localized 3d moments. The other possibility for the origin of the A TFe signal is the formation of μ −Al in the surroundings of the sample. The lifetime difference between μ −Al and μ −O could provide a comparable relaxation to the experiment.
Figure 10 shows the temperature dependencies of Δ−, ν
−, and λ
KT
− for Li[Ni1/2Mn3/2]O4. A superscript of each symbol, i.e., “-” means that the data obtained from μ
−SR measurements. Although the observed Δ− slightly increases with temperature, Δ− is approximately temperature independent and very consistent with the predicted value at the oxygen site in the Li[Ni1/2Mn3/2]O4 lattice. Here, the predicted value (
![Figure 10:
The temperature dependencies of the μ
−SR parameters for Li[Ni1/2Mn3/2]O4, (a) the field distribution width (Δ−), (b) the field fluctuation rate (ν
−), and (c) the exponential relaxation rate (
λ
KT
−
${\lambda }_{\text{KT}}^{-}$
). “-” of each symbol means that the data was obtained from μ
−SR measurements. Δ−c, ν
−c, and
λ
KT
−
c
${\lambda }_{\text{KT}}^{-\text{c}}$
are the fitting results with a common Δ− (
Δ
-
,
c
=
0
.
30
(
3
)
μ
s
-
1
${{\Delta}}^{-,\text{c}}=0.30(3)\hspace{0.17em}\mu {\text{s}}^{-1}$
) in the whole temperature range measured. In (a), a dotted line represents Δ−,calc at the oxygen site predicted with dipole field calculations. The μ
+SR results, i.e., Δ+, ν
+, and
λ
KT
+
${\lambda }_{\text{KT}}^{+}$
are also plotted for comparison.](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_010.jpg)
The temperature dependencies of the μ
−SR parameters for Li[Ni1/2Mn3/2]O4, (a) the field distribution width (Δ−), (b) the field fluctuation rate (ν
−), and (c) the exponential relaxation rate (
The value of ν
− determined from fitting with Δ− as a free parameter increases with temperature with a slope (dν
−/dT) that is steeper than that of ν
+. This is likely to be due to a correlation of the fitted parameters. However, when the spectra are fitted using Δ−c instead, the increase in ν
−c is comparable to that of ν
+. Finally,
4 Discussion
Here, we attempt to estimate the self-diffusion coefficient of Li+ (
![Figure 11:
The relationship between ν and inverse temperature. Solid lines represent the best fit using a thermal activation process:
ν
=
ν
0
exp
[
−
E
a
/
(
k
B
T
)
]
$\nu ={\nu }_{0}\mathrm{exp}[-{E}_{a}/({k}_{\text{B}}T)]$
, in the temperature range between 200 and 400 K, where E
a
is an activation energy and k
B is the Boltzmann constant. The right vertical axis is converted into
D
Li
J
${D}_{\text{Li}}^{\text{J}}$
using Eq. (8).](/document/doi/10.1515/zpch-2021-3102/asset/graphic/j_zpch-2021-3102_fig_011.jpg)
The relationship between ν and inverse temperature. Solid lines represent the best fit using a thermal activation process:
In order to estimate
If we assume that the field fluctuation rate ν
± corresponds to a hopping rate, i.e., the jump frequency of Li+ ions between the 8a and 16c sites,
where N
i
is the number of Li sites in the i-th jump path, Z
v,i
is the vacancy fraction, and s
i
is the jump distance. Based on the above consideration, n = 1, N
1 = 4, Z
1 = 1, and s
1 = 1.77 Å. Since
5 Summary and future
We have measured not only μ
+SR but also μ
−SR spectra for a spinel-type battery material, Li[Ni1/2Mn3/2]O4, in order to investigate the diffusing species above around 250 K. The observed dynamic behavior both with μ
+SR and μ
−SR evidenced Li diffusion. There exists a wide range of cathode materials for a Li-ion battery, which are mainly classified within the following three groups; layered rocksalt-, spinel-, and olivine-type compounds. We have already confirmed Li diffusion with μ
±SR in the spinel Li4Ti5O12 [34] and Li[Ni1/2Mn3/2]O4, and the olivine LiMnPO4 [26]. Moreover, for a typical layered rocksalt-type compound, Li
x
CoO2, the estimated
Funding source: Swedish Research Council http://dx.doi.org/10.13039/501100004359
Award Identifier / Grant number: Dnr. 2016-06955
Funding source: Japan Society for the Promotion Science http://dx.doi.org/10.13039/501100001691
Award Identifier / Grant number: JP18H01863
Award Identifier / Grant number: JP20K21149
Acknowledgments
We thank the staff of ISIS, in particular David Hunt for making the sample cell, for help with the μ ±SR experiments (RB1720324 [58] and RB2010138 [59]).
-
Author contributions: J.S. conceived the experiments. K.O., O.K.F., M.M., and J.S. performed the μ ±SR experiments in ISIS together with K.I. and the local ISIS support from S.P.C. and A.D.H. J.S. analyzed the data and created the first draft. All authors review the manuscript in several steps and have accepted responsibility for the entire content of this submitted manuscript and approved submission.
-
Research funding: M.M. and O.K.F. were partly supported by the Swedish Research Council (VR) through a neutron project grant (BIFROST, Dnr. 2016-06955). This work was supported by Japan Society for the Promotion Science (JSPS) KAKENHI Grant No. JP18H01863 and JP20K21149.
-
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
References
1. Heitjans, P., Indris, S. J. Phys. Condens. Matter 2003, 15, R1257; https://doi.org/10.1088/0953-8984/15/30/202.Search in Google Scholar
2. Grey, C. P., Dupr, N. Chem. Rev. 2004, 104, 4493; https://doi.org/10.1021/cr020734p.Search in Google Scholar
3. Chandran, C. V., Heitjans, P. Chapter One - Solid-State NMR Studies of Lithium Ion Dynamics Across Materials Classes, Vol. 89 of Annual Reports on NMR Spectroscopy; Academic Press, 2016; pp. 1–102.10.1016/bs.arnmr.2016.03.001Search in Google Scholar
4. Kalvius, G. M., Noakes, D. R., Hartmann, O. Handbook on the Physics and Chemistry of Rare Earths, Vol. 32, chap. 206; North-Holland: Amsterdam, 2001; pp. 55–451; https://doi.org/10.1016/s0168-1273(01)32005-6.Search in Google Scholar
5. Yaouanc, A., de Réotier, P. D. Muon Spin Rotation, Relaxation, and Resonance, Application to Condensed Matter, Oxford, New York, 2011.Search in Google Scholar
6. Sugiyama, J., Mukai, K., Ikedo, Y., Nozaki, H., Månsson, M., Watanabe, I. Phys. Rev. Lett. 2009, 103, 147601; https://doi.org/10.1103/physrevlett.103.147601.Search in Google Scholar PubMed
7. Wu, C. S., Ambler, E., Hayward, R. W., Hoppes, D. D., Hudson, R. P. Phys. Rev. 1957, 105, 1413; https://doi.org/10.1103/physrev.105.1413.Search in Google Scholar
8. Hayano, R. S., Uemura, Y. J., Imazato, J., Nishida, N., Yamazaki, T., Kubo, R. Phys. Rev. B 1979, 20, 850; https://doi.org/10.1103/physrevb.20.850.Search in Google Scholar
9. Sugiyama, J. J. Phys. Soc. Jpn. 2013, 82, SA023; https://doi.org/10.7566/jpsjs.82sa.sa023.Search in Google Scholar
10. Månsson, M., Sugiyama, J. Phys. Scripta 2013, 88, 068509; https://doi.org/10.1088/0031-8949/88/06/068509.Search in Google Scholar
11. Sugiyama, J., Ikedo, Y., Mukai, K., Nozaki, H., Månsson, M., Ofer, O., Harada, M., Kamazawa, K., Miyake, Y., Brewer, J. H., Ansaldo, E. J., Chow, K.H., Watanabe, I., Ohzuku, T. Phys. Rev. B 2010, 82, 224412; https://doi.org/10.1103/physrevb.82.224412.Search in Google Scholar
12. Sugiyama, J., Nozaki, H., Harada, M., Kamazawa, K., Ofer, O., Månsson, M., Brewer, J. H., Ansaldo, E. J., Chow, K. H., Ikedo, Y., Miyake, Y., Ohishi, K., Watanabe, I., Kobayashi, G., Kanno, R. Phys. Rev. B 2011, 84, 054430; https://doi.org/10.1103/physrevb.84.054430.Search in Google Scholar
13. Sugiyama, J., Nozaki, H., Harada, M., Kamazawa, K., Ikedo, Y., Miyake, Y., Ofer, O., Månsson, M., Ansaldo, E. J., Chow, K. H., Kobayashi, G., Kanno, R. Phys. Rev. B 2012, 85, 054111; https://doi.org/10.1103/physrevb.85.054111.Search in Google Scholar
14. Sugiyama, J., Mukai, K., Nozaki, H., Harada, M., Månsson, M., Kamazawa, K., Andreica, D., Amato, A., Hillier, A. D. Phys. Rev. B 2013, 87, 024409; https://doi.org/10.1103/physrevb.87.024409.Search in Google Scholar
15. Sugiyama, J., Mukai, K., Harada, M., Nozaki, H., Miwa, K., Shiotsuki, T., Shindo, Y., Giblin, S. R., Lord, J. Phys. Chem. Chem. Phys. 2013, 15, 10402; https://doi.org/10.1039/c3cp51662h.Search in Google Scholar PubMed
16. Sugiyama, J., Nozaki, H., Umegaki, I., Mukai, K., Miwa, K., Shiraki, S., Hitosugi, T., Suter, A., Prokscha, T., Salman, Z., Lord, J.S., Månsson, M. Phys. Rev. B 2015, 92, 014417; https://doi.org/10.1103/physrevb.92.014417.Search in Google Scholar
17. Umegaki, I., Kawauchi, S., Sawada, H., Nozaki, H., Higuchi, Y., Miwa, K., Kondo, Y., Månsson, M., Telling, M., Coomer, F. C., Cottrell, S. P., Sasaki, T., Kobayashi, T., Sugiyama, J. Phys. Chem. Chem. Phys. 2017, 19, 19058; https://doi.org/10.1039/c7cp02047c.Search in Google Scholar PubMed
18. Powell, A. S., Lord, J. S., Gregory, D. H., Titman, J. J. J. Phys. Chem. C 2009, 113, 20758; https://doi.org/10.1021/jp9091249.Search in Google Scholar
19. Baker, P. J., Franke, I., Pratt, F. L., Lancaster, T., Prabhakaran, D., Hayes, W., Blundell, S. J. Phys. Rev. B 2011, 84, 174403; https://doi.org/10.1103/physrevb.84.174403.Search in Google Scholar
20. Powell, A. S., Stoeva, Z., Lord, J. S., Smith, R. I., Gregory, D. H., Titman, J. J. Phys. Chem. Chem. Phys. 2013, 15, 816; https://doi.org/10.1039/c2cp43318d.Search in Google Scholar PubMed
21. Johnson, I. D., Ashton, T. E., Blagovidova, E., Smales, G. J., Lübke, M., Baker, P. J., Corr, S. A., Darr, J. A. Sci. Rep. 2018, 8, 4114; https://doi.org/10.1038/s41598-018-22435-1.Search in Google Scholar PubMed PubMed Central
22. Sugiyama, J., Nozaki, H., Umegaki, I., Harada, M., Higuchi, Y., Ansaldo, E. J., Brewer, J. H., Miyake, Y., Kobayashi, G., Kanno, R. J. Phys. Conf. 2014, 551, 012012; https://doi.org/10.1088/1742-6596/551/1/012012.Search in Google Scholar
23. Jensen, A. C. S., Olsson, E., Au, H., Alptekin, H., Yang, Z., Cottrell, S., Yokoyama, K., Cai, Q., Titirici, M.-M., Drew, A. J. J. Mater. Chem. A 2020, 8, 743; https://doi.org/10.1039/c9ta10113f.Search in Google Scholar
24. Matsubara, N., Nocerino, E., Forslund, O. K., Zubayer, A., Papadopoulos, K., Andreica, D., Sugiyama, J., Palm, R., Guguchia, Z., Cottrell, S. P., Kamiyama, T., Saito, T., Kalaboukhov, A., Sassa, Y., Masese, T., Månsson, M., Månsson. Sci. Rep. 2020, 10, 18305; https://doi.org/10.1038/s41598-020-75251-x.Search in Google Scholar PubMed PubMed Central
25. Sugiyama, J., Umegaki, I., Nozaki, H., Higemoto, W., Hamada, K., Takeshita, S., Koda, A., Shimomura, K., Ninomiya, K., Kubo, M. K. Phys. Rev. Lett. 2018, 121, 087202; https://doi.org/10.1103/physrevlett.121.087202.Search in Google Scholar PubMed
26. Sugiyama, J., Forslund, O. K., Nocerino, E., Matsubara, N., Papadopoulos, K., Sassa, Y., Cottrell, S. P., Hillier, A. D., Ishida, K., Månsson, M., Brewer, J. H. Phys. Rev. Res. 2020, 2, 033161; https://doi.org/10.1103/physrevresearch.2.033161.Search in Google Scholar
27. Sugiyama, J., Umegaki, I., Takeshita, S., Sakurai, H., Nishimura, S., Forslund, O. K., Nocerino, E., Matsubara, N., Månsson, M., Nakano, T., Yamauchi, I., Ninomiya, K., Kubo, M. K., Shimomura, K. Phys. Rev. B 2020, 102, 144431; https://doi.org/10.1103/physrevb.102.144431.Search in Google Scholar
28. Nagamiya, S., Nagamine, K., Hashimoto, O., Yamazaki, T. Phys. Rev. Lett. 1975, 35, 308; https://doi.org/10.1103/physrevlett.35.308.Search in Google Scholar
29. Nagamine, K., Nagamiya, S., Hashimoto, O., Kobayashi, S., Yamazaki, T. Hyperfine Interact. 1976, 2, 407; https://doi.org/10.1007/bf01021193.Search in Google Scholar
30. Brewer, J. H. Hyperfine Interact. 1984, 19, 873; https://doi.org/10.1007/bf02066132.Search in Google Scholar
31. Nishida, N. Hyperfine Interact. 1993, 79, 823; https://doi.org/10.1007/bf00567614.Search in Google Scholar
32. Winston, R. Phys. Rev. 1963, 129, 2766; https://doi.org/10.1103/physrev.129.2766.Search in Google Scholar
33. Buckle, D. C., Kane, J. R., Siegel, R. T., Wetmore, R. J. Phys. Rev. Lett. 1968, 20, 705; https://doi.org/10.1103/physrevlett.20.705.Search in Google Scholar
34. Sugiyama, J. Unpublished work in J-PARC; 2021.Search in Google Scholar
35. Ariyoshi, K., Iwakoshi, Y., Nakayama, N., Ohzuku, T. J. Electrochem. Soc. 2004, 151, A296; https://doi.org/10.1149/1.1639162.Search in Google Scholar
36. Amdouni, N., Zaghib, K., Gendron, F., Mauger, A., Julien, C. J. Magn. Magn Mater. 2007, 309, 100; https://doi.org/10.1016/j.jmmm.2006.06.018.Search in Google Scholar
37. Goodenough, J. B., Kim, Y. Chem. Mater. 2010, 22, 587; https://doi.org/10.1021/cm901452z.Search in Google Scholar
38. Shimoda, K., Murakami, M., Komatsu, H., Arai, H., Uchimoto, Y., Ogumi, Z. J. Phys. Chem. C 2015, 119, 13472; https://doi.org/10.1021/acs.jpcc.5b03273.Search in Google Scholar
39. Mukai, K., Sugiyama, J. J. Electrochem. Soc. 2010, 157, A672; https://doi.org/10.1149/1.3374338.Search in Google Scholar
40. Sugiyama, J., Nozaki, H., Umegaki, I., Mukai, K., Cottrell, S. P., Shiraki, S., Hitosugi, T., Sassa, Y., Suter, A., Salman, Z., Prokscha, T., Månsson, M. JPS Conf. Proc. 2018, 21, 011015.Search in Google Scholar
41. Momma, K., Izumi, F. J. Appl. Crystallogr. 2008, 41, 653; https://doi.org/10.1107/s0021889808012016.Search in Google Scholar
42. Matsuzaki, T., Ishida, K., Nagamine, K., Watanabe, I., Eaton, G., Williams, W. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 2001, 465, 365; https://doi.org/10.1016/s0168-9002(01)00694-5.Search in Google Scholar
43. Hillier, A. D., Lord, J. S., Ishida, K., Rogers, C. Phil. Trans. Math. Phys. Eng. Sci. 2019, 377, 20180064; https://doi.org/10.1098/rsta.2018.0064.Search in Google Scholar PubMed PubMed Central
44. Ziegler, J. F., Ziegler, M., Biersack, J. Nuclear instruments and methods in physics research section B: beam interactions with materials and atoms 268. In 19th International Conference on Ion Beam Analysis, 2010, p. 1818; https://doi.org/10.1016/j.nimb.2010.02.091.Search in Google Scholar
45. Giblin, S. R., Cottrell, S. P., King, P. J. C., Tomlinson, S., Jago, S. J. S., Randall, L. J., Roberts, M. J., Norris, J., Howarth, S., Mutamba, Q. B., Rhodes, J., Akeroyd, F. A. Nucl. Instrum. Methods Phys. Res. 2014, 751, 70; https://doi.org/10.1016/j.nima.2014.03.010.Search in Google Scholar
46. Suter, A., Wojek, B. Phys. Procedia 2012, 30, 69; https://doi.org/10.1016/j.phpro.2012.04.042.Search in Google Scholar
47. Mukai, K., Ikedo, Y., Kamazawa, K., Brewer, J. H., Ansaldo, E. J., Chow, K. H., Månsson, M., Sugiyama, J. RSC Adv. 2013, 3, 11634; https://doi.org/10.1039/c3ra40878g.Search in Google Scholar
48. Kubo, R., Toyabe, T. Magnetic Resonance and Relaxation; North-Holland: Amsterdam, 1996.Search in Google Scholar
49. Mukai, K., Sugiyama, J. Solid State Commun. 2010b, 150, 906; https://doi.org/10.1016/j.ssc.2010.03.003.Search in Google Scholar
50. Matsuzaki, T., Nishiyama, K., Nagamine, K., Yamazaki, T., Senba, M., Bailey, J. M., Brewer, J. H. Phys. Lett. A 1989, 123, 91.10.1016/0375-9601(87)90664-5Search in Google Scholar
51. Keren, A. Phys. Rev. B 1994, 50, 10039; https://doi.org/10.1103/physrevb.50.10039.Search in Google Scholar PubMed
52. Abragam, A. Principles of Nuclear Magnetism, International Series of Monographs on Physics; Oxford University Press: Oxford, 1961.10.1063/1.3057238Search in Google Scholar
53. Suzuki, T., Measday, D. F., Roalsvig, J. P. Phys. Rev. C 1987, 35, 2212; https://doi.org/10.1103/physrevc.35.2212.Search in Google Scholar PubMed
54. Kojima, K. M., Yamanobe, J., Eisaki, H., Uchida, S., Fudamoto, Y., Gat, I. M., Larkin, M. I., Savici, A., Uemura, Y. J., Kyriakou, P. P., Rovers, M. T., Luke, G. M. Phys. Rev. B 2004, 70, 094402; https://doi.org/10.1103/physrevb.70.094402.Search in Google Scholar
55. Borg, R. J., Dienes, G. J. An Introduction to Solid State Diffusion; Academic Press: San Diego, 1988; pp 53–77, chap. 3; https://doi.org/10.1016/b978-0-12-118425-4.50008-2.Search in Google Scholar
56. der Ven, A. V., Ceder, G. Electrochem. Solid State Lett. 2000, 3, 301.10.1149/1.1391130Search in Google Scholar
57. Nakamura, K., Ohno, H., Okamura, K., Michihiro, Y., Moriga, T., Nakabayashi, I., Kanashiro, T. Solid State Ionics 2006, 177, 821; https://doi.org/10.1016/j.ssi.2006.02.021.Search in Google Scholar
58. Sugiyama, J. How Muons See the Motion of Negative Hydrogen Ions in Solids; STFC ISIS facility, 2020.Search in Google Scholar
59. Sugiyama, J. Negative muon spin rotation and relaxaion on spinel Li[Ni1/2Mn3/2]O4, STFC ISIS facility, 2021. https://doi.org/10.5286/ISIS.E.RB2010138.Search in Google Scholar
© 2021 Jun Sugiyama et al., published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.
Articles in the same Issue
- Frontmatter
- Preface
- Special issue on the occasion of the 75th birthday of Paul Heitjans
- Contribution to Special Issue dedicated to Paul Heitjans
- Unusual cation coordination in nanostructured mullites
- A novel high entropy spinel-type aluminate MAl2O4 (M = Zn, Mg, Cu, Co) and its lithiated oxyfluoride and oxychloride derivatives prepared by one-step mechanosynthesis
- Two new quaternary copper bismuth sulfide halides: CuBi2S3Cl and CuBi2S3Br as candidates for copper ion conductivity
- Sintering behavior and ionic conductivity of Li1.5Al0.5Ti1.5(PO4)3 synthesized with different precursors
- Status and progress of ion-implanted βNMR at TRIUMF
- How Li diffusion in spinel Li[Ni1/2Mn3/2]O4 is seen with μ ±SR
- Nuclear magnetic resonance (NMR) studies of sintering effects on the lithium ion dynamics in Li1.5Al0.5Ti1.5(PO4)3
- Anion reorientations and cation diffusion in a carbon-substituted sodium nido-borate Na-7,9-C2B9H12: 1H and 23Na NMR studies
- Site preferences and ion dynamics in lithium chalcohalide solid solutions with argyrodite structure: I. A multinuclear solid state NMR study of the system Li6PS5-xSexI and of Li6AsS5I
- Site preferences and ion dynamics in lithium chalcohalide solid solutions with argyrodite structure: II. Multinuclear solid state NMR of the systems Li6PS5−x Se x Cl and Li6PS5−x Se x Br
- Independent component analysis combined with Laplace inversion of spectrally resolved spin-alignment echo/T 1 3D 7Li NMR of superionic Li10GeP2S12
- How the cation size impacts on the relaxational and diffusional dynamics of supercooled butylammonium-based ionic liquids: DPEBA–TFSI versus BTMA–TFSI
- Solid-state NMR studies of non-ionic surfactants confined in mesoporous silica
- Inorganic-organic hybrid materials based on the intercalation of radical cations: 2-(4-N-methylpyridinium)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazol-1-oxyl-3-N-oxide in fluoromica clay
- Lithium tracer diffusion in near stoichiometric LiNi0.5Mn1.5O4 cathode material for lithium-ion batteries
- On the CaF2-BaF2 interface
- The ionic conductivity of alkali aluminum germanium phosphate glasses – comparison of Plasma CAIT with two electrode DC measurements
- Thin-film chemical expansion of ceria based solid solutions: laser vibrometry study
- Predicting conductivities of alkali borophosphate glasses based on site energy distributions derived from network former unit concentrations
- Ionic transport in K2Ti6O13
- F anion transport in nanocrystalline SmF3 and in mechanosynthesized, vacancy-rich Sm1—x BaxF3—x
- An overview of thermotransport in fluorite-related ionic oxides
Articles in the same Issue
- Frontmatter
- Preface
- Special issue on the occasion of the 75th birthday of Paul Heitjans
- Contribution to Special Issue dedicated to Paul Heitjans
- Unusual cation coordination in nanostructured mullites
- A novel high entropy spinel-type aluminate MAl2O4 (M = Zn, Mg, Cu, Co) and its lithiated oxyfluoride and oxychloride derivatives prepared by one-step mechanosynthesis
- Two new quaternary copper bismuth sulfide halides: CuBi2S3Cl and CuBi2S3Br as candidates for copper ion conductivity
- Sintering behavior and ionic conductivity of Li1.5Al0.5Ti1.5(PO4)3 synthesized with different precursors
- Status and progress of ion-implanted βNMR at TRIUMF
- How Li diffusion in spinel Li[Ni1/2Mn3/2]O4 is seen with μ ±SR
- Nuclear magnetic resonance (NMR) studies of sintering effects on the lithium ion dynamics in Li1.5Al0.5Ti1.5(PO4)3
- Anion reorientations and cation diffusion in a carbon-substituted sodium nido-borate Na-7,9-C2B9H12: 1H and 23Na NMR studies
- Site preferences and ion dynamics in lithium chalcohalide solid solutions with argyrodite structure: I. A multinuclear solid state NMR study of the system Li6PS5-xSexI and of Li6AsS5I
- Site preferences and ion dynamics in lithium chalcohalide solid solutions with argyrodite structure: II. Multinuclear solid state NMR of the systems Li6PS5−x Se x Cl and Li6PS5−x Se x Br
- Independent component analysis combined with Laplace inversion of spectrally resolved spin-alignment echo/T 1 3D 7Li NMR of superionic Li10GeP2S12
- How the cation size impacts on the relaxational and diffusional dynamics of supercooled butylammonium-based ionic liquids: DPEBA–TFSI versus BTMA–TFSI
- Solid-state NMR studies of non-ionic surfactants confined in mesoporous silica
- Inorganic-organic hybrid materials based on the intercalation of radical cations: 2-(4-N-methylpyridinium)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazol-1-oxyl-3-N-oxide in fluoromica clay
- Lithium tracer diffusion in near stoichiometric LiNi0.5Mn1.5O4 cathode material for lithium-ion batteries
- On the CaF2-BaF2 interface
- The ionic conductivity of alkali aluminum germanium phosphate glasses – comparison of Plasma CAIT with two electrode DC measurements
- Thin-film chemical expansion of ceria based solid solutions: laser vibrometry study
- Predicting conductivities of alkali borophosphate glasses based on site energy distributions derived from network former unit concentrations
- Ionic transport in K2Ti6O13
- F anion transport in nanocrystalline SmF3 and in mechanosynthesized, vacancy-rich Sm1—x BaxF3—x
- An overview of thermotransport in fluorite-related ionic oxides