Kit
McColl
and
Furio
Corà
*
Department of Chemistry, University College London, London, UK. E-mail: f.cora@ucl.ac.uk
First published on 30th January 2019
Anatase TiO2 is a candidate high-power electrode material in Li-ion and Na-ion batteries and has been explored as a Mg battery cathode material, although Mg capacity in undoped anatase is limited. Here we use hybrid-exchange density functional theory calculations to investigate the underlying factors affecting Mg intercalation and mobility in anatase. We find that at the dilute limit, Mg ions have 5-fold coordinated insertion sites, and activation barriers for migration are a surprisingly low 537 meV. As the concentration of Mg inserted into the structure is increased, a cooperative distortion of the lattice occurs, contracting the c lattice parameter. The distortion results in stable orderings of Mg ions in sites which are 6-fold coordinated, but also results in migration barriers that exceed 1500 meV in Mg0.5TiO2 due to a collective relaxation of the host lattice. The total increase in barrier is predominantly a result of the stabilisation of the insertion sites, as opposed to a destabilisation of the activated sites along the migration pathway. The insertion sites in the dilute limit can be described as frustrated, and it is this unfavourable insertion geometry under dilute conditions that allows the Mg ions to migrate with low activation barriers. The limited performance for Mg2+ storage can therefore be attributed to the loss of frustrated coordination at high Mg concentration, which restricts Mg mobility and therefore capacity. Strategies to enhance the capacity of Mg in anatase should therefore aim to inhibit the c lattice parameter contraction or otherwise destabilise stable orderings of Mg in Mg0.5TiO2 to retain the frustrated coordination of Mg ions at high Mg concentrations.
Mg is heavier and less electropositive than Li, but these deficits are compensated by its divalency and the reduced weight of a battery design that avoids the anode intercalation host, meaning specific capacities of Mg battery systems can approach or exceed those of Li-ion batteries using a graphite anode.2 Magnesium is also inherently safer than lithium in a battery system, and is far more abundant and cheaper than lithium, making it a more scalable technology for future energy storage solutions. Issues finding Mg battery electrolytes compatible with both electrodes remain,3 however a key bottleneck in the development of Mg batteries is the search for a suitable cathode material.2 Intercalation of divalent ions presents a number of challenges compared to monovalent ions. Multivalent cathode materials are more prone to conversion reactions, and this is a particular issue for Mg due to the high thermodynamic stability of MgO.2 Divalent Mg ions also experience a stronger electrostatic interaction with the anionic framework of the cathode host than Li ions, so despite their similar size, Mg ions typically experience significantly higher activation barriers for movement, leading to poor kinetics.6 Suitable cathode materials are therefore limited.
The pioneering work of Aurbach first demonstrated that Mg battery chemistry was possible, yet the Chevrel phase Mo6S8 cathode used offered a limited voltage.7 Since this work, extensive experimental and computational work has investigated materials that may be suitable for Mg intercalation cathodes, often motivated by their known performance as Li-ion intercalation hosts.2 One of the materials that has been the subject of numerous studies as Li intercalation host is TiO2.8–10 Despite not offering high operating voltages, the anatase11 and TiO2 (B)12 polymorphs of titania show good rate performance for Li intercalation, making them attractive for high power applications. Anatase TiO2 has also been explored as a possible cathode material for Na-ion batteries.13 As a Mg-battery cathode, intercalation to Mg0.5TiO2 provides a high theoretical capacity of 294.5 mA h g−1. Furthermore, anatase TiO2 displays excellent stability, and the low toxicity and high abundance of titanium make TiO2 an attractive sustainable battery material.
Comparatively, there has been little research on TiO2 polymorphs as a Mg battery material. Early work by Sheha achieved only limited capacities using an anatase cathode and a Mg anode.14 Su et al. used a LiBH4/Mg (BH4)2/tetraglyme electrolyte and a Mg metal anode to achieve a reversible low rate (0.2C) capacity of 145 mA h g−1, and 85 mA h g−1 at rates of 2C, with an anatase cathode demonstrating co-intercalation of Mg and Li.15 A similar hybrid Li+/Mg2+ electrochemical setup by the same group achieved pseudocapacitive storage of 115 mA h g−1 at 2C using a TiO2 (B) cathode.16 Meng et al. demonstrated a reversible capacity of 35 mA h g−1 in TiO2 (B) nanoflakes via a double-layer capacitive process in a Mg only arrangement. In contrast, a Li+/Mg2+ hybrid arrangement provided a capacity of 114 mA h g−1 at 1 A g−1, but predominantly via Li intercalation and surface pseudocapacitive charge storage, rather than Mg-intercalation.17 Zhang et al. found that in a Mg only battery design, electrochemical insertion of Mg was only possible up to a stoichiometry of Mg0.1TiO2 in commercially available anatase particles (<20 nm), giving a reversible capacity of 40 mA h g−1 after 50 cycles.18 Koketsu et al. reported that F− and OH− incorporation into anatase nanoparticles (∼5 nm) could enhance voltage, rate performance and capacity for Mg2+ insertion (∼165 mA h g−1 at a rate of 20 mA g−1), and also allow the reversible insertion of Al3+ ions.19 The anatase was determined to have a high concentration (∼22%) of Ti4+ vacancies to charge balance substitution of F− into O2− sites, and these sites were identified as the origin of the improved capacity, whilst also offering a vacancy-mediated diffusion mechanism, and increased insertion energies.20 Li-insertion rates into this Ti-deficient system were also shown to be enhanced in a separate work by the same group.21 Recently, Sheng et al. demonstrated that V-doped anatase nanoparticles (∼10 nm) could sustain a capacity of 121.9 mA h g−1 at 20 mA g−1. The introduction of V was predicted by calculations to provide enhanced voltages compared to undoped TiO2. However the charge storage was predominantly of a faradaic pseudocapacitive type, rather than a diffusion-controlled intercalation mechanism.22
The enhanced capacity reported by Koketsu et al. suggests that there may be doping strategies that can make anatase a viable Mg battery cathode material. Conveniently, doped anatase has been studied as a Li-ion battery electrode material in previous experimental literature, and with improvements and changes in performance noted for a range of dopants including Ni,23 Nb,24,25 Sn,26 Fe,27 and Mo,28 providing an indication of routes to explore. However for doping strategies to be rationally designed, the origin of the limited performance for undoped anatase should be fully understood.
Here we employ hybrid density functional theory calculations to address this problem. We investigate the insertion geometry of Mg, voltage and ion mobility at dilute concentrations. We then assess how these properties change as further Mg is added to the structure. Finally, we discuss doping strategies that may be suitable to enhance the limited experimentally observed capacity.
The local geometry of the Mg ion is indicated in Fig. 2. Mg–O bond distances for Mg located in the geometric centre of the octahedral hole are 1.90 Å in the xy plane and 2.75 Å in the z direction. The displacement from the centre of the octahedral hole site results in a 5-fold distorted square pyramidal coordination, with Mg–O distances of 2.055 Å to the peak of the pyramid, and four further bonds of 1.959 Å and 1.953 Å to ions we label O(a), and 1.939 Å and 1.935 Å to O(b). The Mg–O distance across the octahedral centre is 3.18 Å.
The calculated intercalation energy for a single Mg ion is −1.13 eV. This is in reasonable agreement with Koketsu et al. who find a value of −1.02 eV using GGA + U.19 Legrain et al. find the insertion energy to be −1.74 eV using GGA, which however misrepresents electronic localisation on reduced Ti3+.37
The insertion of an Mg2+ ion introduces two electrons for charge balance, which localise on two separate Ti ions (Fig. S2†), reducing them formally to Ti3+, according to the following equation (in Kröger–Vink notation):
(1) |
The Ti–O bond lengths of these reduced Ti ions are modified (Table S2†) by a combination of electrostatic distortion due to the insertion of the positively charged Mg ion, and the electrons localised on each Ti3+. The electronic density of states of TiO2 and Mg0.03TiO2 are indicated in Fig. S1 and S3.† The two Ti ions reduced upon insertion of Mg result in a split double defect state ∼0.8 eV below the conduction band edge.
Upon insertion of Mg to form Mg0.03TiO2, the lattice parameters of anatase change with an expansion of the a and b parameters by +1.4% and +0.4% respectively (where the a direction is along the Mg–O(a) bonds). The c parameter contracts by −0.8%.
Tunnels extending in the x and y direction intersect at the octahedral insertion sites. However since Mg ions in these sites are displaced from the centre, this means that under dilute conditions they only occupy a single tunnel extending in one direction (x or y) at any moment. To move between tunnels, Mg ions must jump across the centre of the insertion site (Fig. 3B) which we denote as pathway 2. The hop distance for pathway 2 is 1.01 Å, which corresponds to twice the distance that the Mg ions relax away from the centre of the site (Fig. 2). Moving between tunnels allow ions to migrate on a 2D [001] plane. Using an alternating combination of pathways 1 & 2, Mg ions can achieve mobility in the z direction, and thus 3D mobility is possible. As we shall discuss later, structural changes reduce the length of pathway 2 to zero at high Mg2+ concentrations, and 3D mobility can be achieved directly by hops between octahedral insertion sites, along what is effectively only pathway 1. The activation barrier for pathway 2 (221 meV) is lower than for pathway 1 (537 meV) under dilute Mg conditions, it is therefore not rate limiting and for these reasons we will not discuss pathway 2 further.
Fig. 4 Local geometry and energetics of the Mg ion along the migration profile under a dilute concentration of Mg. (a) Geometry of the Mg ion in the stable insertion site. (b) Energetics of the Mg ion along the migration pathway, indicating the geometry in the stable and metastable sites. (c) Geometry of the Mg ion in the metastable site at 50% path length. A more detailed description of the local geometry of the Mg ion in the activated and metastable sites is provided in the ESI.† |
In the metastable intermediate site the Mg ion is 6-fold coordinated, in a distorted octahedral environment (Fig. S5†). In the activated site the Mg ion has a similar coordination (Fig. S6†), but with a non-linear O–Mg–O axial bond angle of 170°, and a closer proximity to the two O(a) ions it was initially coordinated to. The difference between the geometry of the two sites is relatively small however, resulting in the small energy difference (31 meV). If viewed normal to the direction of migration, the pathway the Mg ion takes in the structure is sinusoidal, rather than linear, and thus does not strictly follow the line of the shared octahedral edge. Initially the Mg ion moves towards the O-ion at the end of the Mg–O pyramid, in the c direction, compressing the Mg–O bond to 1.914 Å (from 2.054 Å), but retaining a five-fold coordination with a bond of 2.06 Å to the O(b) ion along the b axis opposite to the direction of motion. Crossing the activation barrier, the Mg ion moves in the b direction, breaking the bond to the more distant O(b) ion, but forming new bonds with two O(a) ions to achieve a 6-fold coordination. In the plateau region, the Mg ion has 6-fold coordination, and the second half of the pathway is a reflection of the first. The change in coordination number along the pathways at the dilute limit is therefore 5 → 6 → 5.
The barrier of 537 meV is relatively low for Mg-ion mobility in a cathode material, and is similar to that of Li mobility in anatase (500 meV) under dilute conditions.44 Generally activation barriers predicted for the migration of multivalent ions are significantly higher than for Li migration, due to the stronger electrostatic interactions between the multivalent ion and the anionic framework of the cathode.2 The barrier of 537 meV corresponds to a diffusivity (D) of ∼4.5 × 10−12 cm2 s−1, assuming an Arrhenius-like activated process, D = va2exp(−Em/kT), where v is the jump frequency (≈1012 s−1), a is the jump distance (2.36 Å for anatase) and Em is the activation barrier.6 We note that such an assumption neglects any contribution from Coulombic interactions between Mg ions and electrons localised upon Ti3+ ions (small polarons), which are known to contribute to limited diffusion for Li-ion cathode materials45 and may affect the mobility of Mg ions.46,47 For comparison with anatase, the Mg diffusion barriers in the α and δ polymorphs of V2O5, which is considered one of the most promising oxides for Mg ion battery cathodes2 are ∼1100 and ∼600 meV respectively.48,49 Similarly low barriers to that of anatase, of ∼600–800 meV have been identified in spinel-structured oxides.50 Interestingly, the migration barrier we find in anatase TiO2 is similar in magnitude to the lowest barrier found for Mg mobility in layered α-MoO3 (580 meV). This migration pathway in α-MoO3 is not between the layers, but rather an intra-layer pathway that takes a similar route along the line of shared octahedral edges, between 5-coordinated stable sites.51 It has been estimated by Rong et al. that for adequate battery operation, migration barriers for ion mobility must be below ∼525 meV for micron-sized particles, and ∼650 meV for nanoparticles.6 The barrier for Mg diffusion in anatase under dilute conditions falls below this limit for nano-size particles.
Since the thermodynamics allow additional Mg insertion, the limited capacity observed experimentally may therefore have a kinetic origin. In the anatase LixTiO2 system, insertion of Li leads to a stable ordering at Li0.5TiO2 and a tetragonal to orthorhombic distortion with increasing Li content.54 The orthorhombic distortion yields significantly reduced Li mobility,44 and insertion of Li into bulk anatase particles is often limited to Li0.5TiO2–Li0.6TiO2, with phase segregation into Li-rich and Li-poor regions.38,39,44 We therefore continue our investigation, considering whether similar structural changes occur in the MgxTiO2 system.
Fig. 6 Calculated formation energies of MgxTiO2 structures. The dotted line links the ground state configurations for TiO2 and Mg0.5TiO2. |
Fig. 7 Structural changes between TiO2 and Mg0.5TiO2, illustrated with a (2 × 2 × 1) unit cell. The Mg0.5TiO2 ordering shown is the most stable calculated (Fig. S13b†). |
Fig. 8 Coordination of Mg ions in the lowest energy configuration of Mg0.5TiO2 illustrated in Fig. S13b.† Mg–O bond lengths are indicated by arrows pointing to the bonds. |
As noted earlier, not only does the structure undergo an orthorhombic distortion, but there is also a significant contraction in the c lattice parameter by −10.7% (Table S3†) as the Ti3+ ions adopt a more regular octahedral coordination. It is this c axis contraction that allows the Mg ions to coordinate to both O-ions along the c axis of the octahedral hole in Mg0.5TiO2. This coordination is not possible in Mg0.03TiO2, where the distances from the centre of the octahedral interstitial site to the O-ions along the c axis are too long for efficient contact for ions the size of Mg2+ and Li+. Under those conditions, the Mg2+ ion relaxes towards one O-ion and is thus 5-fold coordinated (Fig. 2). A cooperative effect at higher concentrations of Mg in the structure is required to contract the lattice along the c axis enough for each inserted Mg ion to be able to achieve 6-fold coordination.
As it has been noted earlier, and indicated in Fig. 5, the energy gained upon insertion of each Mg ion increases as the concentration of Mg is increased. We can therefore conclude that the 6-fold coordination of the Mg ions is more favourable than the 5-fold coordination at the dilute limit, and it is the structural distortion upon increasing Mg concentration that stabilises further the inserted Mg ions. Specifically it is the contraction of the c parameter that allows the more stable 6-fold coordination to be achieved for each Mg ion.
Fig. 9 indicates the activation barriers for each of these pathways, compared with the barrier for mobility at the dilute limit. The barriers in Mg0.5TiO2 are 2120, 1710 and 1639 meV, for pathways (i), (ii) and (iii) respectively, the lowest of which corresponds to ∼17 orders of magnitude reduction in Mg diffusivity at room temperature, compared to Mg0.03TiO2 (537 meV), if Arrhenius-like behaviour is assumed. These significantly increased barriers are likely to be the origin of the low Mg capacities observed experimentally. We note that ions moving along pathway (i) hop to sites that are similar in energy to the initial Mg position (+51 meV), whereas Mg ions moving along pathways (ii) and (iii) reach sites that are 660 and 880 meV higher in energy respectively than their starting positions. This mean that Mg ions moving from less favourable sites to the stable sites will have lower barriers to overcome (∼980 meV and ∼650 meV for the reverse of pathways (ii) and (iii) respectively). This behaviour will tend to favour the formation of stable Mg-rich ordered regions, whereby Mg ions will be pinned to stable lattice configurations and it will be difficult for Mg ions to be extracted from these regions. These observations complement the results of the formation energy calculations (Fig. 6a), which indicates a preference for phase segregation into TiO2 and Mg0.5TiO2 regions.
Fig. 9 Activation barriers for Mg moving along three possible symmetry inequivalent pathways in Mg0.5TiO2. The y axes are plotted indicating the activation barrier with respect to the most stable Mg site in Mg0.5TiO2 (left), and the insertion energy per Mg ion (right). Pathways (i)–(iii) are defined in Fig. S14 of the ESI.† Also plotted for reference (against insertion energy) is the barrier for Mg moving along pathway 1 under dilute conditions (Mg0.03TiO2), smoothed for clarity. |
As indicated in Fig. 5, Mg ions become more stabilised (i.e. their insertion energies become more negative) as the concentration of Mg in the structure is increased, and this can be attributed to a more favourable 6 fold coordination for the Mg ions in Mg0.5TiO2. As the Mg ions become more stabilised, the energy difference between the insertion sites and the activated sites increases, and so does the activation barrier.
It is instructive to consider the energetics of the Mg ions relative to metallic Mg to provide a consistent scale of their stability (Fig. 9, right hand y axis). In this energy frame, the transition state energy in the Mg0.03TiO2 and Mg0.5TiO2 phases is similar. The change in activation barrier (∼1100 meV) between pathway 1 in Mg0.03TiO2 and pathways (ii) and (iii) in Mg0.5TiO2 is due mainly to the Mg stabilisation in Mg0.5TiO2 as measured by the increased insertion energy per ion (∼900 meV) between these two concentrations. The energy of the activated site varies by only ∼250 meV. Thus, for pathways (ii) and (iii) it is the stabilisation of the insertion sites, rather than the destabilisation of the activated sites at high Mg concentrations that make the most substantial contribution to the increase in activation barriers.
We also note that pathways (ii) and (iii) connect intercalation sites of non-equal energy in Mg0.5TiO2, i.e. stable 6-coordinated sites and less stable 5-fold coordinated sites.
The existence of stabilised and destabilised sites adjacent to each other suggests that the distortions that result in stabilisation of certain sites are fairly localised. This result indicates that if cooperative distortions of the anatase lattice that yield 6-fold coordinated sites can be prevented, migration barriers for Mg ions in the Mg0.5TiO2 lattice could be reduced significantly. This would result in higher Mg mobility at high Mg concentration in anatase and allow for greater capacity.
Whilst the stabilisation of Li insertion sites has not been explicitly addressed as the origin for the substantial increases in activation barriers for Li moving in LixTiO2 at high Li concentrations,44 the displacement of Li ions from the centre of the octahedral insertion site, structural distortions and increase in insertion energy per Li between 0 < x < 0.5 follow the same trends as we describe here for Mg,41,54 and we therefore suggest the behaviour is equivalent. It follows from this that strategies to improve Mg capacity in anatase may also be applicable for improving Li storage capacity.
The use of structures in which Mg ions have ‘frustrated’ or ‘non-ideal’ coordination has been proposed as a design principle for achieving high Mg ion mobility.6 Frustrated coordination has to date predominantly been identified as a property of metastable polymorphs,47 such as Mn2O4 spinel.50 A further example is the metastable tunnel bronze ζ-V2O5 phase, which has been demonstrated as a viable Mg ion cathode material,58 and is predicted by DFT calculations to have low activation barriers (620–860 meV) for Mg mobility.46 In fact it has been suggested that insertion of multivalent ions into frustrated coordination sites ‘almost certainly’ requires a structure that is thermodynamically metastable.6 Whilst metastable phases like the V2O5 tunnel bronzes may be favourable for high ion mobility, if they are significantly higher in energy than the ground state polymorphs (α-V2O5), they may be prone to phase transitions under repeated electrochemical cycling, or conversion reactions.2
The identification of frustrated coordination for Mg ions in anatase, and the indication of similar effect for Li-ions, thus indicates that frustration is not a property unique to materials that would generally be classed as metastable phases. Rutile is the ground state polymorph of TiO2 under ambient conditions, but the energy difference between anatase and rutile is small.59,60 Anatase TiO2 has high thermal and chemical stability, does not readily undergo transformation to rutile, and is not generally considered ‘metastable’ in the same way as, for example, chemically de-intercalated tunnel bronze structures.47
Anatase TiO2 therefore provides a stable structure in which frustrated coordination of Mg (and Li-ions) occurs under dilute conditions and this leads to low migration barriers for the ions. The limitation of the anatase structure is therefore not its stability under cycling, but the loss of frustrated coordination at high Mg (or Li) concentrations.
Strategies to increase the electrochemical performance of anatase as a Mg (or Li) ion cathode material should therefore aim to retain the frustrated coordination of the Mg (Li) ions, by inhibiting the c parameter contraction that results in stable Mg (Li) orderings.
One route to achieve this effect could be via isovalent substitution of larger ions than Ti4+ (i.e. Zr4+, r = 0.72 Å or Hf4+, r = 0.71 Å), which would expand the anatase lattice parameters.61,62 Incorporation of Nb5+ (r = 0.64 Å) is also known to cause expansion of the anatase lattice,63 and improvements have been demonstrated for Li ion storage in Nb-doped anatase.24,64
The isovalent substitution of V4+ (r = 0.58 Å) for Ti4+ by Sheng et al. was not found to improve intercalation capacity for Mg in anatase, but instead offered improved pseudocapacitive storage.22 In light of the current discussion, we speculate that the smaller size of V4+ compared to Ti4+ does not cause lattice expansion, and is unlikely to inhibit the distortions discussed above.
An alternative strategy is the creation of Ti4+ vacancies that can provide insertion sites for Mg.19,20 By providing vacancies into which Mg can intercalate, the cation deficient structure may be resistant against the structural distortions (i.e. c parameter contraction) that we have identified as stabilising Mg and thus limiting mobility and capacity. The distortion arises from Mg2+ filling the interstitial sites, whereas insertion into Ti vacancy sites is not expected to create the same distortion. The authors report the creation of Ti vacancies, achieved by an anionic doping strategy, substituting F− and OH− into oxide ion sites. An alternative method proposed to create Ti vacancies in anatase is doping with Mo6+.28,65 Higher-valent cation doping is a general strategy to induce structural cation vacancies,66 and may be effective in anatase.
It is worth nothing however that the chemistry taking place during heterovalent substitutions (Nb5+, Mo6+ into Ti4+ sites) is not always well understood, with a complex interplay of charge compensating defects occurring in the lattice, with possible clustering of vacancies in various combinations which are often challenging to probe experimentally or represent in models.
A rule of thumb to judge the effectiveness of doped anatase phases towards Mg intercalation is to quantify the c/a ratio, with high values expected to yield better resistance to loss of structure.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8ta09939a |
This journal is © The Royal Society of Chemistry 2019 |