Introduction

The improvement of Li-ion batteries (LIB) is focused on enhancing the specific power and energy density. To improve the battery energy density, it's necessary to utilize a high voltage positive electrode material. LiCoPO4 (lithium cobalt phosphate, LCP) is of particular interest due to a high theoretical capacity of around 167 mAh g−1, and high operating voltage (OV) ~ 4.8 V vs. Li+/Li1,2,3,4,5,6. LiCoPO4 belongs to the LiMPO4 olivine family (M = Co, Mn, Fe, and Ni) and demonstrates mediocre rate performance and cycle stability because of low diffusion of lithium ions in the lattice and poor electronic conductivity7,8,9. To increase the electrochemical performance of LiCoPO4, different strategies were used, such as coating of LCP with various carbon sources to increase electronic conductivity1,2,10,11,12,13,14 and nanosizing to improve ionic conductivity.

Moreover, in several works, cobalt ions were substituted with different transition metals such as Ni, Mn, and Fe, and this substitution caused significant improvement in the electronic and ionic conductivity of LiCoPO4, thus enhancing its electrochemical performance15,16,17,18,19. The improvement in the bulk conductivity is not due to antisite defects (Co or Fe on the Li site); but from the mobile polarons, and associated with lithium vacancies20,21. While the operating voltage was not noticeably lowered, the charge carrier transport in the crystal structure was improved, leading to a stabilization of the charging process and a decrease in the number of reactive species formed on the electrode surface.

During the electrochemical reactions, the electrode materials suffer from a significant structural change. In order to reliably track this structural evolution and relate it to the specific chemical transitions, it is necessary to use the operando approach. Operando XRD allows one to study full dynamics of the sample structure, track intermediate phases and avoid sample relaxation. This allows a deeper understanding of the internal mechanisms affecting the electrochemical characteristics of the material and its degradation, providing an understanding of the improvement of the properties of the material.

The changes in the crystallography that occur upon the LiCoPO4 to CoPO4 transition have been studied extensively6,14,22,23. Bramink and co-authors found that pure LiCoPO4 undergoes an electrochemical de-lithiation process by two subsequent two-phase steps, creating a Li2/3CoPO4 intermediate structure24. Strobridge and co-authors25 detected two intermediate phases in Fe-substituted LiCoPO4 upon de-lithiation. The 1st phase is Li2/3FexCo1−xPO4 for all concentrations, but the 2nd phase is Li1−x(Fe3+)x(Co2+)1−xPO4, for 0 < x < 1. Moreover, it was observed that for low lithium concentrations (0 < x < 0.5), de-lithiation between Li1−x(Fe3+)x(Co2+)1−xPO4 and Li2/3FexCo1−xPO4 occurs through a bulk single-phase mechanism25. Kang et al. have demonstrated by theoretical and experimental techniques that the doping by Fe2+ in the olivine structure reduces the Li-Co antisite mixing26. This effect refers to the preferred occupation of the cobalt site by iron, which leads to the expansion of an oxygen octahedron and improves the diffusion of Li+ in channels27. Furthermore, Allen and co-authors observed that Fe3+ has the ability to inhibit capacity fading, and the LiCoPO4 doped by Fe has better electronic conductivity than the pure LiCoPO428.

An in-depth understanding of the relationship among the structural characteristics and electrochemical routine is crucial for the rational and conscious development of doped LCPO materials; therefore, a combination of operando X-ray diffraction and absorption spectroscopy is useful to study the dynamics of the local atomic and electronic structure of electrode material during electrochemical transformations. In this paper, we employ such a combination at the synchrotron source to investigate the Fe doped LCP cathode material with a general formula LiFe0.5Co0.5PO4. We also use ab initio DFT calculations to address the peculiarities of redox processes in this dual-ion system.

Materials and methods

LCFP synthesis

The LiCo0.5Fe0.5PO4 was synthesized by mechanochemically assisted carbothermal reduction. The reagents of Li2CO3, Fe2O3 (a-modification), Co3O4, and (NH4)2HPO4 (qualification “pure for analysis”) were mixed with carbon “P 277” by a high-energy AGO-2 planetary mill under argon atmosphere for 5 min in a molar ratio (Co + Fe):P:Li = 1:1:3. The resulting product was annealed under argon flow for one hour at 750 °C. The final product was a single-phase LiCoFePO4 with an orthorhombic structure10.

Characterization using electrochemistry

Electrodes for the electrochemical characterization were made through mixing around 80 wt% of LCFP active component with 15 wt% of carbon (Timcal Super P Conductive) and 5 wt% of PVDF/NMP binder. The resulting slurry was cast on the aluminum foil and dried at 150 °C under a vacuum overnight. The sample mass loading was around 2–3 mg per cm2 after drying, and an electrode with a diameter of 10 mm was employed. The cells were assembled in an argon glovebox with a lithium foil as an anode, a glass fiber filter (Whatman, Grade GF/C) as a separator, and 10 μl of 1 M LiPF6 solution in 1:1 (by weight) EC:DMC (Sigma Aldrich) as an electrolyte. The assembled cells were cycled in a galvanostatic mode from 3.0 to 5.0 V vs. Li/Li+ at RT with a charge/discharge rate of C/10.

Operando synchrotron-based XAS and XRD characterization

X-ray absorption spectra and simultaneous diffraction patterns were obtained in the operando mode at the European Synchrotron Radiation Facility ESRF's Swiss-Norwegian Beamline BM31 (Grenoble, France). Home made electrochemical cells with X-ray transparent windows made from glassy carbon were used for cycling29. The resulting cathode powder was placed in the cell with a mass load of 13 mg/cm2. The cells were cycled galvanostatically using a BAT-SMALL battery analyzer (Astrol Electronic AG, Switzerland) in the region from 3 to 5 V vs Li/Li+ with a cycling rate of C/10. Diffraction patterns were recorded each 5.5 min, while the Fe and Co K-edge absorption spectra—for 9.5 min, resulting in a total measurement time of 15 min per point.

2D diffraction images were performed with 5 s acquisition time by the Dexela 2D CMOS detector (PerkinElmer, USA). The wavelength λ = 0.49796 Å and the sample-to-detector distance of 289.55 mm were calibrated using silicon powder and LaB6 standards.

Fe and Co K-edge X-ray absorption spectra were measured in a single continuous scan in transmission mode employing a Si(111) monochromator. Ionization chambers were used to monitor intensity before and after assembling, and a Gd2O3 reference was measured simultaneously with the sample for energy alignment.The observed spectra were aligned, flattened, and normalized using the Demeter package's Athena function30,31,32. Then all spectra were interpolated on the same set of energy points and submitted as the input for the principle component analysis (PCA) via open-source PyFitIt software as described precisely in the Mathematical details of the PCA procedure are presented below, and were not included in the main text, but were appropriately referenced:

Any XANES spectrum can be projected onto an orthogonal basis set of functions. The set of projections over each basis function then describes the entire spectrum. We construct basis set functions in the procedure of principal component analysis applied to the theoretical training set. A theoretical XANES dataset \(\mathbf{X}\) is decomposed in the following way:

$$\mathbf{X}={\mathbf{U}}_{(\mathrm{m}\times \mathrm{m})}{{\varvec{\Sigma}}}_{(\mathrm{m}\times \mathrm{n})}{{\mathbf{V}}^{\mathrm{T}}}_{(\mathrm{n}\times \mathrm{n})}$$
(1)

where \(\mathbf{X}\) can be considered as a matrix of dimensions \((m\times n)\), where m is the number of XANES energy points while n is the number of theoretical spectra constituting it, \(\mathbf{U}\) and \(\mathbf{V}\) are two square unitary matrices, \({\varvec{\Sigma}}\) is a diagonal rectangular matrix while \(\mathrm{T}\) denotes the transpose operator. The diagonal elements of \({\varvec{\Sigma}}\) are referred to as the singular dataset values, whose magnitudes are proportional to the amount of variance of the related component. Columns in matrix \(\mathbf{U}\) have the dimensionality of a XANES spectrum and are referred to as abstract mathematical components (i.e. they do not look like line spectra, but their proper linear combination does). Matrix \({\varvec{\Sigma}}\mathbf{V}\) provides the weights which need to be employed to reconstruct each spectrum of \(\mathbf{X}\) from \(\mathbf{U}\). Considering Eq. (1), the ith XANES spectrum \({{\varvec{\upmu}}}_{\mathrm{i}}\) of \(\mathbf{X}\), can be rewritten as31:

$${{\varvec{\upmu}}}_{\mathrm{i}}\left(\mathbf{E},\mathbf{p}\right)=\sum_{\mathrm{j}=1}^{\mathrm{m}}{\mathrm{h}}_{\mathrm{ij}}(\mathbf{p}){\mathbf{u}}_{\mathrm{j}}(\mathbf{E})$$
(2)

where \(\mathbf{E}=({\mathrm{E}}_{1},\dots {\mathrm{E}}_{\mathrm{m}})\) is the set of XANES energy points, \({\mathbf{u}}_{\mathrm{j}}\) represents the jth column vector of \(\mathbf{U}\) while \({h}_{\mathrm{ij}}\) is the fraction of the jth component in the ith spectrum provided by matrix \({\varvec{\Sigma}}\mathbf{V}\)32.

Ab initio DFT calculations

The VASP software was used to perform ab initio DFT calculations33,34 using GGA-PBE projector-augmented wave pseudopotentials with spin polarization. GGA + U correction was 5.05 and 3.84 for Co and Fe, respectively6. A two-step geometry optimization incorporating cell volume and then atomic locations was done for each alkali metal concentration. For ionic relaxation and electronic minimization, the conjugate gradient method and the specific Davidson block iteration scheme techniques were utilized. Following that, optimized structures were employed for precise single-point computations of total energy, charge densities, and atomic charges using the Bader technique35,36. To determine Bader charge values for reference compounds including Fe and Co with known oxidation states the same computational techniques were used.

Result and discussion

According to synchrotron powder XRD (Fig. 1), the as-synthesized LiFe0.5Co0.5PO4 is characterized as a single orthorhombic phase without any carbon peaks according to the open crystallography database reference #9004884. Rietveld refinement of the XRD profile suggests that it has a Pbnm space group with the lattice parameters of a = 4.69 Å, b = 10.25 Å, and c = 5.96 Å.

Figure 1
figure 1

Rietveld refinement of the synchrotron XRD pattern of as-synthesized LiCo0.5Fe0.5PO4.

Figure 2 displays the first two charges/discharge profiles of LiCo0.5Fe0.5PO4 cathode. Two voltage plateaus were observed, the first one at ~ 3.5 V corresponds to the Fe2+/Fe3+ redox couple, and the second one at ~ 4.9 V—to the Co2+/Co3+ redox couple. This material demonstrates the capacity of around 89 mAh/g and 93 mAh/g for the first and the second discharges, respectively. The 1st specific charge capacity is slightly higher than the theoretical capacity of 167 mAh/g, which indicates that some side reactions occur, which might include the CEI (cathode EI) formation and electrolyte decomposition at high voltage. However, in the second cycle, the charge-specific capacity decreases.

Figure 2
figure 2

First two charge/discharge profiles of LiCo0.5Fe0.5PO4.

To reveal the structural transformations upon lithiation/delithiation, we carried out operando synchrotron diffraction experiment. As reported in the previous work, ex situ powder XRD data for several cobalt phosphate compositions (LiCoPO4, Li0.7CoPO4, and CoPO4) exhibit reflections from two olivine phases with various unit cell volume and a substantial contribution from the fraction of an amorphous phase37. We assume that the cathode undergoes partial amorphization during lithium extraction.

Operando X-ray diffraction

Figure 3 displays the evolution of the selected reflections in the operando synchrotron XRD patterns during the first two lithiation/de-lithiation cycles of a cell at the C/10 rate. No amorphization of the active material was observed. At several planes (“0 2 0”, “0 2 1”, and “1 3 0”), this observation is in agreement with Strobridge et al.25. During the first charge, around 0.5 Li per formula unit leaves the structure of LiCo0.5Fe0.5PO4, accompanied by an increase in the lattice parameter “a” to 4.72 Å and a decrease of b and c to 9.92 Å and 5.82 Å, respectively, as shown in Fig. 4a. Since the lithium loss was approximately 0.3 per formula unit in the first discharge, this phase partially recovered to the original LiCo0.5Fe0.5PO4 state, with the following lattice parameters a = 4.715 Å, b = 10.005 Å, and c = 5.855 Å These results are similar to those previously reported for Pnma-Li0.7CoPO4 structure with a = 10.070 Å, b = 5.851 Å and c = 4.717 Å38. It is clear that upon charging, the lattice parameter a is expanding, while b and c are shrinking. At the end of charge, the b and c cell parameters were higher than expected for Fe0.5Co0.5PO4, but c was smaller, revealing that the full conversion was not achieved during the operando experiment. Moreover, upon lithiation and delithation processes, the profiles of the potential remain stable, which proposed small or insignificant changes in the mechanism of the reaction, as demonstrated in Fig. 4b.

Figure 3
figure 3

Operando synchrotron XRD patterns of LiCo0.5Fe0.5PO4 upon the first two charge–discharge cycles in the 2θ range: (a) from 5.3 to 8°, indicating the shift to the lower 2θ values at this region, (b) from 9.5 to 12° indicating the shift to lower 2θ values at this region. The red line represents the charging process, and the blue line is the discharge process.

Figure 4
figure 4

Evolution of (a) lattice parameters and (b) unit cell volume obtained from the operando XRD compared to the cell potential and Li loss upon the cycling of LiCo0.5Fe0.5PO4.

In situ synchrotron X-ray absorption measurements

Figure 5 demonstrates the evolution of operando Fe and Co K-edge XANES spectra of LiCo0.5Fe0.5PO4 cathode upon cycling. Though for Fe, there is a very notable chemical shift of the absorption edge, indicating the Fe2+/Fe3+ redox reaction (Fig. 5a), for Co K-edge, this shift is hardly noticeable (Fig. 5b).

Figure 5
figure 5

A series of operando XANES spectra for the Fe (a) and Co (b) K-edge in the LiCo0.5Fe0.5PO4 cathode material obtained during the first cycle (the red line is a charging process, and the blue line represents the discharge process).

This set of spectra was further processed using PCA. Because the XAS spectrum of a multiphase bulk mixture is a superposition of XAS spectra for each phase, PCA can be used to mathematically decompose the sequence of spectra for the electrochemical phase transition process to get the pristine spectra of the relevant phases and their concentrations. Initially, all spectra were interpolated to a single energy range. Singular value decomposition was then used to obtain a certain number of components from the sequence of spectra. In our situation, two abstract spectra were more transformed into chemically meaningful information via target transformation. The following physical constraints were expected: spectra should be normalized, the values of concentrations should be positive, and at least one of the components under consideration should represent either a fully reduced or fully oxidized phase. In this way, two components were successfully identified for both Co and Fe, as shown in Fig. 6.

Figure 6
figure 6

PCA components extracted from the series of operando XANES spectra for the Fe compared to the experimental spectra for LiFePO4 and FePO4 (ref.40) (a) and Co (b) K-edge in the LiCo0.5Fe0.5PO4 cathode material, compared to the experimental spectra of the reference compounds (LiCoPO4@UiO-66)7. (c) Cell voltage of Fe during the first cycle dots on a black voltage profile mark the start time for measuring each successive XAS spectrum. (d) PCA phase concentration of the Fe2+ and Fe3+ components compared to the cell potential. (e) Cell Voltage of Co during the first cycle dots on a black voltage profile marks the start time for measuring each successive XAS spectrum. (f) PCA phase concentration of the Co2+ and Co3+ components compared to the cell potential.

There is a significant change between the first and the second components because the first component corresponds to Fe2+ and the second—to Fe3+ in CoFePO4, as demonstrated by the comparison with the experimental XAS spectra for LiFePO4 and FePO4 in Fig. 6a. For the Co K-edge, in a similar way, the first component can be attributed to the Co2+ in LiCo0.5Fe0.5PO4 and the second one—to the Co3+ in CoFePO4 due to a rare study of in situ XAS of LiCoPO4, the current result compared to our previously reported work LiCoPO4 coating by UiO-66 meatlorganic framework7, as shown in Fig. 6b.

Figure 6d shows that at the first charge, the Fe2+ ions are almost completely converted to Fe3+, thus, about 50% of Li ions should be extracted from the cathode material. Figure 6e,f demonstrate that the electrode surface remains predominately Co2+ in LiCoPO4, irrespective of the charge process. However, upon discharge process has noted a shift in the center of peaks to the higher energy suggesting the growth of some Co3+, with additional loss of various characteristics in the bulk of the material. The determination of several features missing in Co3+ compounds suggests that (~ 60)Co2+ was not fully oxidized39. There is also a peculiar disturbance in Co phase concentration in the region corresponding to the Fe2+ redox activity. Such a reaction of Co on the change in the local atomic and electronic coordination of Fe might be indicative of the homogeneous distribution of d-metals in the structure of cathode material and the homogeneous Li insertion/extraction mechanism. X-ray diffraction is sensitive to the long-range order, and the conclusion about the solid solution mechanism is made based on this experimental technique. On the contrary, X-ray absorption fine structure is sensitive to the local atomic structure within 5–7 Å around the absorber. Therefore the “components” retrieved from PCA or linear combination fits are related to the local order around 3d metals but not to crystalline phases.

To understand both cobalt and iron performance during the charge/discharge process of LiCoFePO4, ab initio DFT calculations were carried out to estimate the local atomic charges according to the Bader method for different Li stoichiometry in LixCoFePO4 (x = 0, 0.25, 0.5, 0.75 and 1). Figure 7 demonstrates the unit cell averaged atomic charges for cobalt, iron, and oxygen compared to the reference values. The calculation suggests that upon the lithiation, Fe2+ should fully convert to Fe3+, which is in agreement with the PCA results shown above. However, the suggested Co3+/Co2+ conversion rate during the lithiation is only slightly above 50%, similar to what was observed in the operando XAS measurements. However, for the oxygen, there is also a substantial change in the Bader charge upon the delitiation, which suggest that there is a reversible anion redox reaction on oxygen, which takes part in the total charge compensation mechanism to make up for the deficiency in Co oxidation.

Figure 7
figure 7

Unit cell averaged local atomic (Bader) charges for Fe (a), Co (b) and O (c) as a function of x in LixCo0.5Fe0.5PO4 (x = 0, 0.25, 0.5, 0.75, and 1). The values for each species were averaged over non-equivalent atomic positions. Reference charge values (dashed lines) were obtained from calculations for reference compounds with different Co and Fe oxidation states.

Conclusion

The current work reports operando characterization of lithium cobalt iron phosphate charge/discharge process as a positive electrode in Li-ion half-cell by synchrotron diffraction and absorption in terms of Co and Fe k-edge. 115 patterns and 200 spectra of synchrotron diffraction and absorption obtained over 36 h of electrochemical cycling exhibited the reaction of reversible changes in redox, structural, and electronic properties for Co and Fe sites during the first cycle. The electrochemical performance was tested. Cycling curves exhibited two voltage plateaus, one corresponding to the Fe2+/Fe3+ redox at 3.5 V, and the second one—to the Co2+/Co3+ redox at 4.9 V. The structural evolution of LiCo0.5Fe0.5PO4 during cycling was tracked by operando synchrotron XRD. Upon the first charge, the lattice parameter “a” increases; however, the “b, c” decrease. On the other hand, a decreases, and b, c increase during the discharge process. Simultaneously, the electronic and local atomic structure of the material was examined using operando synchrotron XAS. The results proved that the Fe2+ ions are fully converted to Fe3+, however, Co2 is only partially oxidized to Co3+. The percentage of 3d metals that undergo redox reactions is determined by quantitative analysis of spectroscopic data using PCA. This value is then compared to the number of electrons transported through the cell to confirm the electrochemically active phase. DFT calculations are performed using experimental data as a benchmark. When the theoretical model reproduces the observed values, it may be expanded to include parameters that the experimental setup does not allow for. In this way, we examined the LixCoFePO4 performance in the whole range of lithium stoichiometry x = 0 . . . 1. When x is smaller than 0.75 in LixCoFePO4, the anionic redox process associated with oxygen atoms in the lattice, which accompanies iron charge shifts and dominates at high voltages, is predicted using DFT and Bader charge analysis.