Local structures around the substituted elements in mixed layered oxides.

The chemical substitution of a transition metal (M) is an effective method to improve the functionality of a material, such as its electrochemical, magnetic, and dielectric properties. The substitution, however, causes local lattice distortion because the difference in the ionic radius (r) modifies the local interatomic distances. Here, we systematically investigated the local structures in the pure (x = 0.0) and mixed (x = 0.05 or 0.1) layered oxides, Na(M1-xM'x)O2 (M and M' are the majority and minority transition metals, respectively), by means of extended X-ray absorption fine structure (EXAFS) analysis. We found that the local interatomic distance (dM-O) around the minority element approaches that around the majority element to reduces the local lattice distortion. We further found that the valence of the minority Mn changes so that its ionic radius approaches that of the majority M.

The chemical substitution of a transition metal (M) is an effective method to improve the functionality of a material, such as its electrochemical123456 magnetic789 and dielectric10 properties. For example, NaFe0.5Co0.5O21 with an O3-type layered structure shows a high discharge capacity of 160 mAh/g and good cyclability, which is much higher than those of the parent O3-NaFeO2 and O3-NaCoO2. In addition, Na0.67Fe0.5Mn0.5O22 with a P2-type layered structure shows a high discharge capacity of 190 mAh/g, which is much higher than that of the parent P2-Na0.67MnO2. To thoroughly comprehend the effects of a partial substitution on the material properties, we first have to know the local structures around the respective Ms, e.g., the interatomic distance and its valence. The extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES) analyses around the K-edges are powerful techniques to determine the local structures around the respective Ms in the mixed compounds. Let us consider the situation in which the ionic radius (r) of the minority (substituted) M is larger than that of the majority (host) M. If the valence state of the minority M is extremely stable, the interatomic distance between the minority M and a neighboring atom is the same as in the parent compound that consists of the minority M. In this case, the resultant local lattice distortion around the minority M increases the Gibbs free energy. Hereafter, we will call this change in energy “distortion energy”. In contrast, if the distortion energy around the minority M is extremely high, M tends to take the same interatomic distance as in the majority M. In this case, the valence state of the minority M may change so that its ionic radius approaches that of the majority M. The actual systems are considered located between these two extreme cases. Thus, the local structure is governed by a subtle competition between the local distortion and the valences.

Among the transition metal compounds, the mixed layered oxides, Na(M1−M′)O2 (M and M′ are the majority and minority transition metals, respectively)1 are promising candidates for the cathode material for a sodium-ion battery11, (2) have large families with different combinations of transition metals, and (3) show simple crystal structure with alternating MO2 layers and Na sheets. The MO2 layer consists of edge-sharing MO6 octahedra formed by covalent bonding. In this sense, the layered oxides are suitable systems for the investigation of the local structures. The layered oxides have two crystal structures, i.e., the O3- and P2-types, depending on the stacking pattern of the MO2 layers and Na sheets. The oxygen sheets of the MO2 layers stack as AB|CA|BC (AB|BA) in the O3-type (P2-type) compounds, where the vertical line represents the Na sheet.

Here, we systematically investigated the interatomic distance around the respective Ms and their valence states in the mixed layered oxides. The systematic EXAFS analyses revealed that the local interatomic distance (d) around the minority M approaches that around the majority M. A similar effect is discernible in the interatomic distance (d) between the neighboring Ms. We further found that the valence of the minority Mn changes so that its ionic radius approaches that of the majority M.

EXAFS analysis

EXAFS analysis is one of the most powerful methods to determine the local interatomic distance around the respective Ms in these mixed compounds. Figure 1 shows a prototypical Fourier transformation of the χ(k)k3 − k plots around the majority Mn and minority Fe in P2-NaMn0.9Fe0.1O2. χ and k are the oscillatory components of the normalized absorption and angular wavenumber, respectively. Both the plots show a two-peak feature, which is attributable to the paths to the first- (O) and second- (TM) nearest neighbor elements, respectively. The peak positions (1.50 and 2.47 Å) in the Mn plot are slightly smaller than those (1.55 and 2.52 Å) in the Fe plot. This indicates that the interatomic distances around the majority Mn are shorter than those around the minority Fe. To evaluate precisely the interatomic distances, we performed least-squares fitting of the FT[χ(k)k3] − R plots with the EXAFS equation (see Method section). In the O3 compounds, we included the contributions from the first- (O), second- (M), and third- (Na) nearest neighbor elements. In the P2 compounds, we included the contributions from the first- (O), third- (M), and fifth- (O) nearest neighbor elements. We omitted the contribution from the second- (Na2) and fourth- (Na1) nearest neighbor elements because their occupancies (g = 0.21 – 0.41) are small. The interatomic distances thus obtained, i.e., d and d, are listed in Table 1.

Figure 1

FT[χ(k)k3]−R plots of the majority Mn and minority Fe in P2-NaMn0.9Fe0.1O2.

The curves were obtained by Fourier transformation of the χ(k)k3−k plots in the k-range from 2 Å−1 to 11.2 Å−1, where χ and k are the oscillatory components of the normalized absorption and angular wavenumber, respectively.

Table 1

Interatomic distances in layered oxides.

Compound	Majority dM-O (Å)	Minority dM-O (Å)	Majority dM-M (Å)	Minority dM-M (Å)
O3-NaFeO2	2.041 (11)	—	3.039 (8)	—
O3-NaFe0.9Mn0.1O2	2.033 (8)	1.930 (9)	3.022 (6)	2.997 (8)
O3-NaFe0.9Co0.1O2	2.037 (8)	1.886 (9)	3.031 (6)	2.858 (12)
O3-NaCoO2	1.905 (6)	—	2.851 (6)	—
O3-NaCo0.95Fe0.05O2	1.910 (6)	2.009 (10)	2.858 (12)	2.915 (9)
P2-NaMnO2	1.907 (8)	—	2.898 (9)	—
P2-NaMn0.9Fe0.1O2	1.908 (6)	1.998 (12)	2.901 (6)	2.905 (11)
P2-NaMn0.9C0.1O2	1.909 (7)	1.914 (9)	2.880 (7)	2.848 (8)

The values were obtained using EXAFS analyses.

Interatomic distance between M and O

In Fig. 2, we plotted the interatomic distances (d) between M and O in the mixed layered oxides against those (d0) in the corresponding pure layered oxides, i.e., O3-NaFeO2, O3-NaCoO2, and P2-NaMnO2. In the majority M [Fig. 2(a)], the d values coincide with the d0 values within experimental error, indicating that the partial (<10%) substitution does not influence the local structure around the majority M.

Figure 2

Interatomic distances (d) between M and O in the mixed layered oxides around (a) majority and (b) minority Ms. Horizontal axes are the M-O distances (d0) in the corresponding pure layered oxides, i.e., O3-NaFeO2, O3-NaCoO2, and P2-NaMnO2.

In the minority M [Fig. 2(b)], however, considerable differences are observed between d and d0. In the Fe-substituted compounds, i.e., O3-NaCo0.95Fe0.05O2 and P2-NaMn0.9Fe0.1O2, dFe-O (=2.009 and 1.998 Å) is shorter than d0Fe-O (=2.041 Å) by 1.6–2.1%. We note that the d of the minority Fe is much longer than d0Co-O (=1.905 Å) and d0Mn-O (=1.907 Å). In short, dFe-O of the minority Fe approaches d0Co-O and d0Mn-O. A similar trend is observed in the Mn-substituted compound. In O3-NaFe0.9Mn0.1O2, dMn-O (=1.930 Å) is longer than d0Mn-O (=1.907 Å) by 1.2% and approaches d0Fe-O (=2.041 Å). We note that the same trend is observed even if we only include the contributions from the nearest neighboring O and M in the EXAFS analyses (Fig. S12).

Interatomic distance between M and M

In Fig. 3, we plotted the interatomic distances (d) between the neighboring Ms in the mixed layered oxides against the mean-field values (dMF). The mean-field values are expressed, for example, in NaFe1-yCoyO2, as dMFCo- = (1 − y)d0Fe-Fe + yd0Fe-Co around Fe and dMFM-M = (1 − y)d0Co-Fe + yd0Co-Co around Co, where d0Fe-Fe and d0Co-Co are the values in O3-NaFeO2 and O3-NaCoO2, respectively. The value d0Fe-Co (=d0Co-Fe) is defined by (d0Fe-Fe + d0Co-Co)/2. In the majority M [Fig. 3(a)], d coincides with the mean-field values within experimental error.

Figure 3

Interatomic distances (d) between the neighboring Ms in the mixed layered oxides around (a) majority and (b) minority Ms. Horizontal axes are the mean-field values (dMF) evaluated from the values of the pure layered oxides, i.e., O3-NaFeO2, O3-NaCoO2, and P2-NaMnO2.

In the minority M [Fig. 3(b)], however, considerable differences are observed between d and dMF. In the Fe-substituted compounds, i.e., O3-NaCo0.95Fe0.05O2 and P2-NaMn0.9Fe0.1O2, dFe- is shorter than the mean-field value, reflecting the compressed dFe-O (<d0Fe-O). A similar trend is observed in the Mn-substituted compound. In O3-NaFe0.9Mn0.1O2, dMn- is longer than the mean-field value, reflecting the elongated dMn-O (>d0Mn-O). We note that the same trend is observed even if we only include the contributions from the nearest neighboring O and M in the EXAFS analyses (Fig. S13).

XANES spectra

Now, let us discuss the valence states of the minority Ms and their relation to the interatomic distance, i.e., d and d. Figure 4 shows the XANES spectra of the pure and mixed layered oxides around the (a) Mn K-edge, (b) Fe K-edge, and (c) Co K-edge. The black curve represents the spectra of the pure layered oxides. The blue and red curves correspond to the spectra of the majority and minority Ms in the mixed layered oxides, respectively. The main peak is attributable to the 1s–4p transition. Their peak shifts are crude measures of the valence change: the valence increases or decreases if the peak shows a blue or red shift, respectively.

Figure 4

XANES spectra of pure and mixed layered oxides around the (a) Mn K-edge, (b) Fe K-edge, and (c) Co K-edge. The black curve represents the spectra of pure layered oxides. The blue and red curves correspond to the spectra of the majority and minority Ms in the mixed layered oxides, respectively. The insets show magnified spectra in the pre-edge region.

In the Mn K-edge spectra [Fig. 4(a)], the main peak of P2-NaMnO2 (black curve) is located at 6556.0 eV. We observed no detectable shift in the majority Mn spectrum (blue curves) in P2-NaMn0.9Fe0.1O2, indicating that the valence state of Mn is unchanged. In the minority Mn spectrum (red curve) in O3-NaFe0.9Mn0.1O2, the peak (=6554.0 eV) shows a significant red shift, indicating a decrease in the Mn valence from Mn3+ (t2g3eg1) with high-spin (HS) configuration. A trace of the valence change is also discernible in the pre-edge spectra [inset of Fig. 4(a)], which is attributable to the 1s–3d transition. The P2-NaMnO2 (black curve) spectrum shows a two-peak feature, whose higher- and lower-lying components are attributable to the 1s–3dt2g and 1s–3 deg transitions, respectively. The relative intensity of the higher peak decreases in O3-NaFe0.9Mn0.1O2 (red curve), suggesting the partial formation of the eg-electron. Thus, the XANES spectrum indicates that the valence of the minority Mn approaches the divalent HS configuration (t2g3eg2). Such a valence change reduces the local lattice distortion around Mn, because the r value increases from r = 0.550 Å for Mn3+ to 0.780 Å for Mn2+. Here, we consider the Na deficiency effect on the Mn valance in O3-NaFe0.9Mn0.1O2, whose actual chemical formula is Na0.86Fe0.9Mn0.1O2 (Table 2). The Na deficiency tends to increases the Mn valence to compensate the decrease in Na+. This trend is opposite to reduction of the Mn valence as suggested by the red shift of the main peak.

Table 2

Lattice constants (a and c), z coordinate (z) of O, and Na percentage (x) of layered oxides.

Compound	a (Å)	c (Å)	z	x
O3-NaFeO2	3.02157 (5)	16.07402 (31)	0.23298 (12)	0.9831 (45)
O3-NaFe0.9Mn0.1O2	3.01862 (12)	16.15361 (65)	0.23249 (21)	0.8632 (77)
O3-NaFe0.9Co0.1O2	3.01954 (7)	16.08690 (40)	0.23275 (14)	0.9771 (53)
O3-NaCoO2	2.89022 (4)	15.60949 (24)	0.22948 (11)	0.9624 (43)
O3-NaCo0.95Fe0.05O2	2.89368 (9)	15.63520 (66)	0.22982 (21)	0.9572 (71)
P2-NaMnO2	2.87337 (8)	11.16667 (57)	0.08771 (37)	0.6734 (169)
P2-NaMn0.9Fe0.1O2	2.87730 (7)	11.19772 (44)	0.08559 (27)	0.7218 (141)
P2-NaMn0.9C0.1O2	2.86177 (13)	11.19690 (84)	0.07971 (34)	0.4443 (163)

In the O3-type compounds, the atomic coordinates were Na (0,0,0), M (0,0,1/2), and O (0,0, z).

The x value is the occupancy at the Na sites. In the P2-type compounds, the atomic coordinates were Na1 (1/3,2/3,3/4), Na2 (0,0,1/4), M (0,0,0), and O (1/3,2/3, z). The x value is the sum of the occupancies at the Na1 and Na2 sites.

In the Fe K-edge spectra [Fig. 4(b)], the main peak of O3-NaFeO2 (black curve) is located at 7126.5 eV. We observed no detectable shift of the majority Fe spectrum (blue curves) in O3-NaFe0.9Co0.1O2. The peak (=7127.5 eV) in the minority Fe spectrum (red curve) in O3-NaCo0.95Fe0.05O2 shows a significant blue shift, suggesting that the Fe valence increases from Fe3+ (t2g3eg2). The pre-edge spectrum [black curve in the inset of Fig. 4(b)] of O3-NaFeO2 is rather broad, suggesting a two-peak feature consistent with HS Fe3+ state (t2g3eg2)12. The spectral weight at the high-energy side increases slightly in O3-NaCo0.95Fe0.05O2 (red curve), implying the partial formation of the eg-hole. Thus, the XANES spectrum implies that the valence of the minority Fe approaches the tetravalent HS configuration. Such a valence change reduces the local lattice distortion around Fe, because the r value decreases from r = 0.645 Å for Fe3+ to 0.585 Å for Fe4+. On the other hand, Ménétrier et al.13 performed 57Fe Mössbauer spectroscopy in Li-overstoichiometric Li1.1Co1FeO2 (y < 0.08), in which Fe occupies the square-pyramidal site. This spectroscopy revealed that Fe at the square-pyramidal site takes the trivalent HS configuration (t2g3eg2). Another possible explanation for the unexpected XANES spectrum is that the spectral change is originated in the impurity phase. A close inspection of the XRD pattern [Fig. S2(b)] reveals weak additional reflections that are unaccounted by the trigonal model. We performed two-phase Rietveld analysis with adding spinel phase (Fd-3m; Z = 8). The additional reflections are well reproduced by 5 wt% spinel phase with a = 8.0805(2) Å (Fig. S14). Judging from the cell parameter, the spinel phase is considered to be Fe-substituted Co3O4. Then, the tetrahedrally-coordinated Fe in the impurity phase may be responsible for the XANES spectrum. We further observed a significant blue shift in the P2-NaMn0.9Fe0.1O2 spectrum (broken red curve), suggesting increases in the Fe valence. We note that P2-NaMnO2 does not show any cooperative Jahn-Teller distortion of HS Mn3+. The valence change, however, may be due to the Na deficiency effect in P2-NaMn0.9Fe0.1O2 (g = 0.31).

Thus, the minority elements adjust their valences so that their ionic radii approach that of the majority M. This empirical relation is effective in the design of functional materials with such a partial substitution because we can predict the valence change of the minority Ms. Specifically, the valence state of the substituted M has a determining effect on the physical properties of the oxide, such as the redox voltage in the battery material, exchange interaction in magnetic material, and electron-lattice interaction in the dielectric material.

In the Co K-edge spectra [Fig. 4(c)], the main peak of O3-NaCoO2 (black curve) is located at 7725.5 eV. We note that the peak position of P2-NaCoO2 is almost the same as that of O3-NaCoO2. We observed no detectable peak shift of the majority Co spectrum (blue curve) in O3-NaCo0.95Fe0.05O2. The peak (=7724.0 eV) in the minority Co (red curve) in O3-NaFe0.9Co0.1O2 shows a significant red shift, indicating that the Co valence decreases from Co3+ (t2g6). The pre-edge spectra of O3-NaFe0.9Co0.1O2 [red curve in the inset of Fig. 4(c)], however, remain with a single peak with no lower-lying absorption corresponding to the 1s–3dt2g transition. In addition, the dCo- value of O3-NaFe0.9Co0.1O2 is much smaller than the mean-field value [Fig. 3(b)], which is curious because Co3+ (t2g6) has the smallest ionic radius (=0.550 Å). These observations imply that Co3+ takes the ligand-hole state (t2g6eg1, where means the ligand hole) rather than the trivalent LS state (t2g6eg0). Due to the complicated ligand-hole state of the Co ion, it seems to be difficult to put forth a simplified argument based only on the ionic radius.

Summary

The systematic EXAFS analyses revealed that the local interatomic distance (d) around the minority M approaches that around the majority M to reduce the local lattice distortion. The valence of the minority Mn changes so that its ionic radius approaches that of the majority M. This empirical relation is effective in the design of functional materials with such a partial substitution.

Method

Sample preparation

Three pure and five mixed layered oxides were prepared from Na2O2, Na2CO3, Co3O4, Fe2O3, and Mn3O4 by solid-state reaction. In Fe-based O3-type compounds, Na2O2, Fe2O3, and the minority M sources were mixed in a 1.2: 1 − x: x atomic ratio and calcined at 943 K for 24 h in air. In the Co-based O3-type compounds, Na2O2, Co3O4, and the minority M sources were mixed in a 1.2: 1 − x: x atomic ratio and calcined at 823 K for 16 h in O2. In the Mn-based P2-type compounds, Na2CO3, Mn3O4, and the minority M sources were mixed in a 0.7: 1 − x: x atomic ratio and calcined at 1273 K for 24 h in air. The laboratory X-ray diffraction patterns show O3- or P2-type structure without detectable impurities.

Crystal structure

To determine the crystal structure and Na percentage (x), we performed synchrotron-radiation X-ray powder diffraction measurements at the BL8A beamline of the Photon Factory, KEK. The samples were finely ground and placed in 0.3 mmϕ glass capillaries. The capillaries were sealed and mounted on the Debye-Scherrer camera. The powder diffraction patterns were detected with an imaging plate. The exposure time was 5–10 minutes. The wavelength (=0.689033 Å) of the X-ray was calibrated by the lattice constant of standard CeO2 powders. The diffraction patterns of the O3-type compounds were analyzed by the Rietveld method (Rietan-FP14: Figs S1 and S2) with a trigonal model (R-m; Z = 3, hexagonal setting). The diffraction patterns of the P2-type compounds were analyzed by the Rietveld method (Fig. S3) with a hexagonal model (P63/mmc; Z = 2). We observed neither traces of impurities nor secondary phases. The obtained parameters, i.e., lattice constant (a and c), z coordinate (z) of O, and Na percentage (x), are listed in Table 2.

X-ray absorption spectroscopy

The X-ray absorption spectroscopy (XAS) measurements were conducted at BL-9C of the Photon Factory, KEK. The powder was finely ground, and mixed with BN, and pressed into pellets with 5 mm in diameter. The weight of the powder in the pellets was optimized by making the edge jump Δμt ≈ 1 for majority elements and 0.32–0.77 for minority elements. The storage ring was operated at 2.5 GeV and 450 mA, and the maximum photon flux was 1 × 1011 photons s−1. The XAS were recorded in a transmission mode with a Si(111) double-crystal monochromator at 300 K. N2 (100%) gas-filled ion camber was used to monitor the incident X-ray intensity and Ar (15%) and N2 (85%) gas-filled ion camber was used to monitor the transmitted intensity. The higher harmonics were eliminated by a detune to 60% of incident X-ray. The ratio of higher harmonics contamination estimated by the detection efficiency of the ion chambers was less than 2 × 10−4. The beam size at the sample position was 0.8 × 0.6  mm (h × w). The wavelengths of the monochromator were calibrated with the absorption edge of Fe, Mn, and Co foil. In X-ray near edge structure (XANES) analyses, the background subtraction and normalization were performed using the ATHENA program15.

EXAFS analysis

The extended X-ray fine structure (EXAFS) analyses were performed using the ATHENA and ARTEMIS programs15. First, the oscillatory components were extracted using the ATHENA program after background subtraction and normalization of the absorption spectra. Thus, we obtained χ(k)k3−k plots, where χ and k are the oscillatory components of the normalized absorption and angular wavenumber, respectively. Then, we performed Fourier transformation of χ(k)k3 in the k-range from 2 Å−1 to 11.2 Å−1.

The plots thus obtained FT[χ(k)k3]−R were analyzed using the ATHENA program (Figs S4–S11). In the plane wave and single-scattering approximation, χ(k) around the K-edge is expressed by the EXAFS equation as:

where S0, Nj, Rj, Fj, σj2, and φj are the passive electron reduction factor, degeneracy of path, path length, effective scattering amplitude, mean square displacement, and effective scattering phase shift of the jth atom, respectively. The angular wavelength k is defined by where m, E, and E0 are the electron mass, energy of the incident X-ray, and energy shift, respectively. In the usual EXAFS analysis, S0, Nj, σj2, E0, and Rj are adjusted to minimize the difference between the experimental data and EXAFS equation. In the O3 compound, we included the contribution from the first- (O), second- (M), and third- (Na) nearest neighbor elements. In the P2 compounds, we included the contributions from the first- (O), third- (M), and fifth- (O) nearest neighbor elements. We omitted the contribution from the second- (Na2) and fourth- (Na1) nearest neighbor elements because the occupancy (g = 0.21–0.41) is small. In both cases, we fixed Nj for the three elements at the crystallographic value (=6). We used the same S0 and E0 values for the three elements. The least-squares fittings were performed in the R range from 1 Å to 3 Å.

The details of the analytical results are shown in Figs S4–S11. We found that the R values of the third paths, i.e. the M–Na distances in the O3 compounds and the longer M - O distances of the P2 compounds, are seriously deviated from the crystallographic values. In addition, the σj2 values for these paths are significantly large (~0.01–0.03 Å2) as compared with those of the first and second paths. In this sense, the EXAFS analysis model including the third paths remains ambiguous. We performed another EXAFS analysis with omitting the third path to validate the discussion above. We observed the same trend of the M – O and M – M distances between the compounds even if we omitted the third path. (Figs S12 and S13) Thus, we confirmed that the R value of the M – O and M – M distances are reliable even though those of the third path remains ambiguous.

Additional Information

How to cite this article: Akama, S. et al. Local structures around the substituted elements in mixed layered oxides. Sci. Rep. 7, 43791; doi: 10.1038/srep43791 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.