Screening of Potential TM-Olivine Candidates As Cathodes for Mg-Ion Batteries
Atomic scale calculations of the olivine structure, MgTMSiO4, with varying transition metals (TM), were performed with the Vienna ab initio simulation package (VASP). This is a density functional theory (DFT) code and the PBE generalized gradient approximation (GGA) was chosen. The unit cell was allowed to relax during ionic optimisations, and single point energy calculations on the final structure were performed for the final electronic relaxation. Hybrid calculations with 25% Hartree-Fock exchange were performed as benchmark single point calculations for some structures, as well as Hubbard GGA+U calculations to investigate the effect of localized strongly correlated d electrons.
We will present results on a number of TMs in order to screen promising candidates from a thermodynamic and kinetic point of view. Some preliminary results are shown in Figure 1. Figure 1a shows that for MgTMSiO4 with TM = V, Mn and Fe, the formation of an intermediary MgxTMSiO4, x=0.25, 0.5 or 0.75, is stable with respect to the extremes with full or no magnesiation. Unit cell expansion during the uptake and release of Mg may be of crucial importance for long term stability of the final cathode. Figure 1b) shows that the maximum unit cell expansion from TMSiO4 to MgTMSiO4 is much more prominent for TM= Mn than TM = V or Fe.
Nudged elastic band (NEB) calculations for the diffusion of Mg between different sites in MgMnSiO4 are shown in Figure 1c and d. Figure 1c shows that moving Mg from to a neighbouring (empty) Mg site has a 1 eV barrier. It should be noted that although the unit cell contains 4 formula units, removing either of the four Mg atoms gives 4 equivalent Mg0.75MnSiO4. This is reflected in Figure 1c where the start and end points of the NEB path are at the same energy level. In the perfect MgMnSiO4 crystal Mg occupies sites denoted 4a, and Mn the sites denoted 4c. However, scrambling is possible, and Figure 1d shows the NEB path for Mg diffusing to an empty Mn (4c) site in the case of MgMn0.75SiO4. The 0.8 eV energy difference between end and start point indicates that empty 4c sites are thermodynamically unfavourable for MgMnSiO4 and that defects arising from Mn vacancies are easily reoccupied by Mg.
Figure 1 Preliminary results of calculations on MgTMSiO4 olivine materials as potential candidates for cathodes in Mg-ion batteries. (a) Reaction energies for x MgTMSiO4 + (1-x) TMSiO4 → MgxTMSiO4. (b) Maximum unit cell expansion during magnesiation. NEB-path with diffusion energies when Mg diffuses from (c) one 4a to a neighbouring 4a position in Mg0.75MnSiO4 (d) from a 4a to a 4c site in MgMn0.75SiO4.