A Computational Study of Defects, Li-Ion Migration and Dopants in Li2ZnSiO4 Polymorphs

Lithium zinc silicate, Li2ZnSiO4, is a promising ceramic solid electrolyte material for Li-ion batteries. In this study, atomistic simulation techniques were employed to examine intrinsic defect processes; long range Li-ion migration paths, together with activation energies; and candidate substitutional dopants at the Zn and the Si sites in both monoclinic and orthorhombic Li2ZnSiO4 phases. The Li-Zn anti-site defect is the most energetically favourable defect in both phases, suggesting that a small amount of cation mixing would be observed. The Li Frenkel is the second lowest energy process. Long range Li-ion migration is observed in the ac plane in the monoclinic phase and the bc plane in the orthorhombic phase with activation energies of 0.88 eV and 0.90 eV, respectively, suggesting that Li-ion diffusivities in both phases are moderate. Furthermore, we show that Fe3+ is a promising dopant to increase Li vacancies required for vacancy-mediated Li-ion migration, and that Al3+ is the best dopant to introduce additional Li in the lattice required for increasing the capacity of this material. The favourable isovalent dopants are Fe2+ at the Zn site and Ge4+ at the Si site.


Introduction
The Li-Ion battery (LIB) is recognized as one of the most promising clean energy devices replacing the current, non-renewable fossil fuels. The development of new electrode and electrolyte materials with excellent electrochemical performance, low costs, environmentally benign natures and high abundance of constituent elements is a key feature in designing a promising Li-ion battery. Significant research activity has been devoted to preparing new classes of materials in the past decade to design high capacity batteries [1][2][3][4][5][6][7][8][9][10].
Silicon based materials are of great interest for designing electrode or electrolyte materials, as silicon is relatively safe, abundant and cheap. Furthermore, strong Si-O bonds in SiO 4 units provide great structural stability. A variety of silicate-based materials have been examined for use as cathode materials for LIBs either experimentally or theoretically [11][12][13][14].
Lithium metal silicates (Li 2 MSiO 4 (M = Fe, Mn, Co, Ni and so on)) are a class of materials in the orthosilicate family [4,[15][16][17][18][19][20][21]. These materials have gained the interest of the research community recently due to their good electrochemical performances and the existence of two Li per formula unit leading to the possibility of extracting more than one lithium. The theoretical capacity of silicate material is reported to be ca 166 mAh −1 with one Li removal and ca 300 mAh −1 for the two Li 2 of 13 removal [22]. Nyten et al. [4] first reported the orthorhombic crystal structure of Li 2 FeSiO 4 using X-ray diffraction. Thereafter, Nishimura et al. [16] synthesized a monoclinic structure of Li 2 FeSiO 4 with their high-resolution X-ray diffraction and transmission electron microscopy. Both crystal structures have been well studied theoretically and their Li-ion conductivities explained [23,24]. After the discovery of Li 2 FeSiO 4 , other lithium metal silicates, such as Li 2 MnSiO 4 , were discovered [18][19][20]. There was a significant debate on the crystal of Li 2 MnSiO 4 . In particular, Dominko et al. [18] reported an orthorhombic structure (space group Pmn2 1 ) which is similar to the orthorhombic Li 2 FeSiO 4 reported by Nyten et al. [4]. Politaev et al. [19] determined a monoclinic structure, which is isostructural with the structure reported by Nishimura et al. [16]. In previous studies [25,26], we examined defects and Li diffusion in polymorphs of Li 2 MnSiO 4 theoretically.
Lithium zinc silicate, Li 2 ZnSiO 4 , was studied experimentally as a candidate material for electrolyte and cathode material for LIBs and transparent glass ceramic material for broadband fibre amplifiers [27][28][29]. Adnan et al. [27] synthesized monoclinic Li 2 ZnSiO 4 ceramic powder using the sol gel method and determined the ionic conductivity at different temperatures. Recently, the orthorhombic structure of Li 2 ZnSiO 4 was synthesized by Sivagami et al. [28] and examined for cathode material for LIBs, though the crystal structure details are not reported. An electrochemical study shows 133 mAhg −1 discharge capacity at C/10 rate, suggesting that this material is comparable to other orthosilicate materials.
The theoretical modelling of materials is an active area of research as it can provide useful information about the crystal structures, defects and diffusion properties. In this study, we used atomistic simulation techniques to model both monoclinic and orthorhombic phases of Li 2 ZnSiO 4 , and examine defects, Li-ion diffusion and dopant properties in both phases. In previous studies, a variety of oxide materials, including Li, Na and Mg ion battery materials have been modelled using this technique [30][31][32][33][34][35][36][37][38].

Computational Methods
Classical, pair-wise potential-based calculations as implemented in the GULP code (version 3.4.1) provided by Curtin University of Technology, Perth, Australia [39] were employed. Here, ionic interactions were modelled using long-range and short-range forces. Long-range ionic interactions are coulombic attractions and short-range forces include electron-electron repulsion and van der Waals interactions (attractive). Short-range repulsive forces were modelled using the Buckingham potentials (refer to  Table S1). Geometry optimizations were carried out using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm [40]. In all optimized structures, the gradient norm was less than 0.001 eV/Å. Point defects were modelled using the Mott-Littleton method [41]. This methodology constructs two spherical regions around the defects. In the inner region, ions are relaxed more accurately. Though the enthalpies calculated in the simulation at hand were anticipated to be overestimates, as the ions are treated as spherical shapes with full charge models and dilute limits, the trend in relative energies is consistent.
Thermodynamically, the defect parameters (for example migration and formation energies) may be defined through the comparison of the real (i.e., defective) crystal to an isochoric or isobaric ideal (i.e., non-defective) crystal. These defect formation parameters can be interconnected via thermodynamic relations [42,43]. The present calculations correspond to the isobaric parameters for both the migration and formation processes [44][45][46][47].

Li 2 ZnSiO 4 Crystal Structures
Two different crystal structures of Li 2 ZnSiO 4 have been reported in the literature. The first structure reported by Yamaguchi et al. [48] crystallizes into a monoclinic phase (space group P2 1 /n), as shown in the Figure 1a. Both Zn 2+ and Si 4+ form corner-sharing tetrahedral units (ZnO 4 and SiO 4 ) Crystals 2019, 9, 563 3 of 13 pointing in opposite directions. Sivagami et al. [28] reported the synthesis of orthorhombic structure of Li 2 ZnSiO 4 (refer to Figure 1b) using a facile sol-gel technique. Although the X-ray diffraction study confirms the orthorhombic phase, crystal structure details, including lattice constants, are not available from that study. confirms the orthorhombic phase, crystal structure details, including lattice constants, are not available from that study. In this study, we examined both monoclinic and orthorhombic polymorphs of Li2ZnSiO4. The initial structure of orthorhombic phase was taken from the orthorhombic Li2FeSiO4 crystal structure reported by Nyten et al. [4]. Full geometry optimization (both ions positions and lattice constants) was performed on both structures. Calculations reproduced the experimental lattice parameters of monoclinic Li2ZnSiO4 well (refer to Table 1), validating the choice of potential parameters selected. Initial and final lattice parameters of orthorhombic phase are also reported, and the difference in the lattice parameters is less than 2.5%. This information would help in the future experimental characterization of orthorhombic Li2ZnSiO4 crystal structure. The present calculations find that the monoclinic crystal structure is 100 meV/formula unit lower in energy than the orthorhombic phase, showing the stability of monoclinic phase over the orthorhombic phase. Table 1. Comparison between the calculated structural parameters and corresponding experimental [48] or initial values [4] reported for monoclinic (P21/n) and orthorhombic (Pmn21) Li2ZnSiO4.  In this study, we examined both monoclinic and orthorhombic polymorphs of Li 2 ZnSiO 4 . The initial structure of orthorhombic phase was taken from the orthorhombic Li 2 FeSiO 4 crystal structure reported by Nyten et al. [4]. Full geometry optimization (both ions positions and lattice constants) was performed on both structures. Calculations reproduced the experimental lattice parameters of monoclinic Li 2 ZnSiO 4 well (refer to Table 1), validating the choice of potential parameters selected. Initial and final lattice parameters of orthorhombic phase are also reported, and the difference in the lattice parameters is less than 2.5%. This information would help in the future experimental characterization of orthorhombic Li 2 ZnSiO 4 crystal structure. The present calculations find that the monoclinic crystal structure is 100 meV/formula unit lower in energy than the orthorhombic phase, showing the stability of monoclinic phase over the orthorhombic phase. Table 1. Comparison between the calculated structural parameters and corresponding experimental [48] or initial values [4] reported for monoclinic (P2 1 /n) and orthorhombic (Pmn2 1 ) Li 2 ZnSiO 4 .

Intrinsic Defect Processes
Isolated point defects (vacancies and interstitials) are necessary to calculate the defect formation energies and Li-ion diffusion. The Frenkel and Schottky defect reaction energies were calculated by combining the point defect energies. The anti-site defect energy was calculated in the forms of isolated and cluster. In the case of the isolated form, substitutional energies of Li at the Zn site and Zn at the Li site were calculated separately, and then they were combined. In the cluster form, both individual impurity defects were present in the same supercell. Kröger-Vink notation [49] is used to represent these defect reactions as shown in Equations (1)- (8). Li vacancy is denoted by V Li and this defect has a negative charge, as the removal of a Li + ion introduces a negative charge in the lattice. The interstitial of Zn has a double positive charge and is denoted by Zn •• i . Li at the Zn site is denoted by Li Zn and this defect is negatively charged due to the sum of −2 charges introduced in the lattice by the Zn vacancy and +1 charge of Li. Superscript single dots and single apostrophes correspond to +1 and −1 charges, respectively.
Li Frenkel : Figure 2 summarises the reaction energies for defect processes. The Li-Zn anti-site defect energy (cluster form) is calculated to be the most dominant defect. This suggests that a small amount of cation inter mixing (Li at the Zn site (Li Zn ) and Zn at the Li site (Zn • Li )) would be present in these materials. The energy difference between the isolated and the cluster form in monoclinic structure is -0.46 eV (-0.35 eV for orthorhombic Li 2 ZnSiO 4 ). This energy is the binding energy between the isolated defects (Li Zn and Zn • Li ) to form the cluster (Li Zn : Zn • Li ) X . The exoergic binding energy suggests that isolated defects have the tendency to form clusters without energy cost. This defect has been observed in a variety of Li, Na and Mg ion battery materials either experimentally or theoretically [4,19,[36][37][38][50][51][52]. Experimental studies show that this defect is mainly due to synthesis conditions and cycling of the as-prepared material. In our previous modelling [25], the Li-Mn anti-site defect cluster was found to be the lowest energy defect process in Li 2 MnSiO 4 . However, the defect energy was higher than that found in Li 2 ZnSiO 4 due to the different charges of Zn and Mn. It is expected that cation exchange in Li 2 ZnSiO 4 (between Li + and Zn 2+ ) would be easier than in Li 2 MnSiO 4 (between Li + and Mn 2+ ) due to the smaller radius of Zn 2+ than that of Mn 2+ . The Li Frenkel is identified as the second most favourable defect energy process. Nevertheless, this process would require moderate temperature. High endothermic energies are observed for other Frenkel and Schottky processes, implying that they are highly unlikely to occur at normal temperatures. We examined the formation of Li 2 O in this material by considering Li 2 O Schottky-like reaction (Equation (6)). The formation energy is 2.35 eV per defect in monoclinic phase and 2.50 eV per defect in orthorhombic phase. The formation of Li 2 O would lead to the formation of V Li and V •• O in this material at high temperatures. Though there is a small energy difference in each reaction energy between both phases, the trend is observed to be the same.

Lithium Ion Diffusion
Long range diffusion of Li-ion with low activation energy is an essential condition for a promising Li-ion battery material. Experimental determination of diffusion pathways is generally challenging. Current methodology permitted us to calculate linear (direct) and non-linear (curved) pathways of various local Li hops together with activation energies. Fisher et al. [32] used the current methodology to examine the Li-ion diffusion pathway in LiFePO 4 . The calculated pathway was non-linear and in excellent agreement with the pathway determined in the neutron diffraction experiment [53].
Long range diffusion of Li-ion with low activation energy is an essential condition for a promising Li-ion battery material. Experimental determination of diffusion pathways is generally challenging. Current methodology permitted us to calculate linear (direct) and non-linear (curved) pathways of various local Li hops together with activation energies. Fisher et al. [32] used the current methodology to examine the Li-ion diffusion pathway in LiFePO4. The calculated pathway was nonlinear and in excellent agreement with the pathway determined in the neutron diffraction experiment [53]. In the monoclinic crystal structure, we identified four different local Li hops (refer to Figure 3) for the vacancy-assisted Li-ion migration. Table 2 lists the local Li hops and corresponding activation energies. Both direct and curved pathways were considered, as the later provided much more freedom to diffuse Li-ions in Li2ZnSiO4. In most of the cases, curved pathways yielded lower activation energies than those calculated in direct pathways (refer to Figure 4). Long range diffusion pathways constructed by connecting local Li hops show that Li-ions will more likely diffuse in the ac plane in a zig-zag pattern (A→D→A→D) as this route has the lowest overall activation energy of 0.88 eV. In Li2MnSiO4, the value was 0.94 eV due to the difference in ionic radii of Zn 2+ and Mn 2+ [25]. Table 3 reports two different long-range diffusion pathways considered, including a pathway (B→C→B→C) with high activation energy of 1.27 eV. In the monoclinic crystal structure, we identified four different local Li hops (refer to Figure 3) for the vacancy-assisted Li-ion migration. Table 2 lists the local Li hops and corresponding activation energies. Both direct and curved pathways were considered, as the later provided much more freedom to diffuse Li-ions in Li 2 ZnSiO 4 . In most of the cases, curved pathways yielded lower activation energies than those calculated in direct pathways (refer to Figure 4). Long range diffusion pathways constructed by connecting local Li hops show that Li-ions will more likely diffuse in the ac plane in a zig-zag pattern (A→D→A→D) as this route has the lowest overall activation energy of 0.88 eV. In Li 2 MnSiO 4 , the value was 0.94 eV due to the difference in ionic radii of Zn 2+ and Mn 2+ [25]. Table 3 reports two different long-range diffusion pathways considered, including a pathway (B→C→B→C) with high activation energy of 1.27 eV. Long range diffusion of Li-ion with low activation energy is an essential condition for a promising Li-ion battery material. Experimental determination of diffusion pathways is generally challenging. Current methodology permitted us to calculate linear (direct) and non-linear (curved) pathways of various local Li hops together with activation energies. Fisher et al. [32] used the current methodology to examine the Li-ion diffusion pathway in LiFePO4. The calculated pathway was nonlinear and in excellent agreement with the pathway determined in the neutron diffraction experiment [53]. In the monoclinic crystal structure, we identified four different local Li hops (refer to Figure 3) for the vacancy-assisted Li-ion migration. Table 2 lists the local Li hops and corresponding activation energies. Both direct and curved pathways were considered, as the later provided much more freedom to diffuse Li-ions in Li2ZnSiO4. In most of the cases, curved pathways yielded lower activation energies than those calculated in direct pathways (refer to Figure 4). Long range diffusion pathways constructed by connecting local Li hops show that Li-ions will more likely diffuse in the ac plane in a zig-zag pattern (A→D→A→D) as this route has the lowest overall activation energy of 0.88 eV. In Li2MnSiO4, the value was 0.94 eV due to the difference in ionic radii of Zn 2+ and Mn 2+ [25]. Table 3 reports two different long-range diffusion pathways considered, including a pathway (B→C→B→C) with high activation energy of 1.27 eV.   Figure 3).    Figure 3).   Three different local Li hops (P, Q and R) were identified in orthorhombic Li2ZnSiO4, as shown in Figure 5. Energy profile diagrams for local Li hops with activation energies are shown in Figure 6. In all cases, the activation energy calculated in a curved pathway is lower than that calculated in a direct pathway (refer to Table 4). A Li-ion diffuses in the bc plane in a zig-zag pattern (P→P→P→P) with the activation energy of 0.90 eV. The activation energy for the second pathway (Q→R→Q→R) is 1.07 eV, and in this pathway a Li-ion migrates in the ac plane (refer to Table 5). Similar pathways were calculated in Li2MnSiO4 and their corresponding activation energies were calculated to be 0.95 eV and 1.29 eV respectively [25], due to the larger radius of Mn 2+ than that of Zn 2+ .

Migration Path Li-Li Separation (Å) Activation Energy (eV) Direct
The pathways examined in this study should help experimentalists in their future experimental work. Activation energies calculated in both structures show that Li-ion conductivity in Li2ZnSiO4 is moderate.  Three different local Li hops (P, Q and R) were identified in orthorhombic Li 2 ZnSiO 4 , as shown in Figure 5. Energy profile diagrams for local Li hops with activation energies are shown in Figure 6. In all cases, the activation energy calculated in a curved pathway is lower than that calculated in a direct pathway (refer to Table 4). A Li-ion diffuses in the bc plane in a zig-zag pattern (P→P→P→P) with the activation energy of 0.90 eV. The activation energy for the second pathway (Q→R→Q→R) is 1.07 eV, and in this pathway a Li-ion migrates in the ac plane (refer to Table 5). Similar pathways were calculated in Li 2 MnSiO 4 and their corresponding activation energies were calculated to be 0.95 eV and 1.29 eV respectively [25], due to the larger radius of Mn 2+ than that of Zn 2+ .
The pathways examined in this study should help experimentalists in their future experimental work. Activation energies calculated in both structures show that Li-ion conductivity in Li 2 ZnSiO 4 is moderate.       Figure 5).

Divalent Dopants
The Zn site was doped by some divalent dopants (M = Ni, Mg, Co, Fe, Ca, Sr and Ba) and the solution energy was calculated using the following equation.
The lowest solution energy is for Fe 2+ in both phases (refer to Figure 7). This suggests that the synthesis of Li 2 Zn 1-x Fe x SiO 4 (0.0 < x < 1.0) should be possible by experimentation. The second most favourable dopant is Mg 2+ in monoclinic phase and Ca 2+ in orthorhombic phase. The energy difference between the solution energies of Mg 2+ and Ca 2+ in orthorhombic phase is very small. This is reminiscent of the early experimental facts that in NaI, where the ionic radius r + of the host cation Na + differs only slightly from the ionic radius r 2+ of the divalent dopant Ca 2+ , the measurements show a small energy difference between the corresponding activation energies for the cation vacancy motion [54] in a similar fashion as observed in other alkali halides as well [55].  Table 4. Calculated Li-Li separations and corresponding activation energies in orthorhombic Li2ZnSiO4 (refer to Figure 5).

Long-Range Path
Overall Activation Energy (eV) P→P→P→P 0.90 Q→R→Q→ R 1.07

Divalent Dopants
The Zn site was doped by some divalent dopants (M = Ni, Mg , Co, Fe, Ca, Sr and Ba) and the solution energy was calculated using the following equation.
The lowest solution energy is for Fe 2+ in both phases (refer to Figure 7). This suggests that the synthesis of Li2Zn1-xFexSiO4 (0.0 < x < 1.0) should be possible by experimentation. The second most favourable dopant is Mg 2+ in monoclinic phase and Ca 2+ in orthorhombic phase. The energy difference between the solution energies of Mg 2+ and Ca 2+ in orthorhombic phase is very small. This is reminiscent of the early experimental facts that in NaI, where the ionic radius r + of the host cation Na + differs only slightly from the ionic radius r 2+ of the divalent dopant Ca 2+ , the measurements show a small energy difference between the corresponding activation energies for the cation vacancy motion [54] in a similar fashion as observed in other alkali halides as well [55].

Trivalent Doping
Here, we consider some trivalent dopants (M = Al, Ga, Fe, In, Sc, Y, Gd and La) on Zn and Si sites. Trivalent dopants at the Zn site would introduce Li vacancies that facilitate Li self-diffusion in Li2ZnSiO4. The solution of X2O3 was considered via the following reaction equation, as written using the Kröger-Vink notation.

Trivalent Doping
Here, we consider some trivalent dopants (M = Al, Ga, Fe, In, Sc, Y, Gd and La) on Zn and Si sites. Trivalent dopants at the Zn site would introduce Li vacancies that facilitate Li self-diffusion in Li 2 ZnSiO 4 . The solution of X 2 O 3 was considered via the following reaction equation, as written using the Kröger-Vink notation.
Solution enthalpies are reported in Figure 8a. Calculations reveal that Fe 3+ is the most favourable dopant in both the monoclinic and orthorhombic structures. Endoergic solution enthalpies suggest that this process can take place at high temperatures. Additionally, it should be stressed that calculations in this study were carried out at 0 K. Therefore, calculations with high temperature would make that process more favourable. Solution enthalpies are reported in Figure 8a. Calculations reveal that Fe 3+ is the most favourable dopant in both the monoclinic and orthorhombic structures. Endoergic solution enthalpies suggest that this process can take place at high temperatures. Additionally, it should be stressed that calculations in this study were carried out at 0 K. Therefore, calculations with high temperature would make that process more favourable. Next, the same dopants were considered at the Si site. Negative charges introduced by substitutional dopants in the lattice were compensated by positively charges Li interstitials. This defect-engineering strategy would increase the capacity of the as-prepared Li2ZnSiO4 material. The following equation was used to calculate the solution energy for this process. Figure 8b reports the solution enthalpies. The most favourable dopant was found to be Al 3+ . The possible composition would be Li2+xZnSi1-xAlxO4 (0.0 < x< 1.0). The exact amount of x should be determined experimentally. A similar strategy was applied to Li2MnSiO4 and Al 3+ was found to be the promising dopant at the Mn site [25]. In all cases, the solution enthalpies calculated for the Next, the same dopants were considered at the Si site. Negative charges introduced by substitutional dopants in the lattice were compensated by positively charges Li interstitials. This defect-engineering strategy would increase the capacity of the as-prepared Li 2 ZnSiO 4 material. The following equation was used to calculate the solution energy for this process. Figure 8b reports the solution enthalpies. The most favourable dopant was found to be Al 3+ . The possible composition would be Li 2+x ZnSi 1-x Al x O 4 (0.0 < x < 1.0). The exact amount of x should be determined experimentally. A similar strategy was applied to Li 2 MnSiO 4 and Al 3+ was found to be the promising dopant at the Mn site [25]. In all cases, the solution enthalpies calculated for the monoclinic phase are slightly lower in energy than those calculated for the orthorhombic phase. The second most favourable dopant is Fe 3+ . Other dopants exhibit high solution enthalpies, meaning they are highly unlikely to take place at normal temperatures.

Tetravalent Doping
Finally, tetravalent dopants (M = Ge, Ti, Sn, Zr and Ce) were considered at the Si site. The following reaction was constructed to calculate the solution enthalpy. As the charges of the dopants and the Si were the same, no charge compensation was needed. Figure 9 shows the solution enthalpies calculated. In both structures, the trend in solution enthalpies is the same and there is only a very small difference in solution enthalpies. The promising dopant for this process is found to be the Ge 4+ . The preference of Ge 4+ is due to the small difference between the ionic radii of Si 4+ (0.26 Å) and Ge 4+ (0.39 Å). The positive solution enthalpy of GeO 2 indicates that energy should be provided in the form of heat for this process. This is due to the stronger Si-O bonds present in the SiO 4 unit than the Ge-O bonds present the in GeO 4 unit. High solution enthalpies are observed for other dopants. Solution enthalpies for TiO 2 and CeO 2 are approximately 5 eV and 6 eV respectively, suggesting that they are highly unlikely to occur. second most favourable dopant is Fe 3+ . Other dopants exhibit high solution enthalpies, meaning they are highly unlikely to take place at normal temperatures.

Tetravalent Doping
Finally, tetravalent dopants (M = Ge, Ti, Sn, Zr and Ce) were considered at the Si site. The following reaction was constructed to calculate the solution enthalpy. As the charges of the dopants and the Si were the same, no charge compensation was needed.
(12) Figure 9 shows the solution enthalpies calculated. In both structures, the trend in solution enthalpies is the same and there is only a very small difference in solution enthalpies. The promising dopant for this process is found to be the Ge 4+ . The preference of Ge 4+ is due to the small difference between the ionic radii of Si 4+ (0.26 Å) and Ge 4+ (0.39 Å). The positive solution enthalpy of GeO2 indicates that energy should be provided in the form of heat for this process. This is due to the stronger Si-O bonds present in the SiO4 unit than the Ge-O bonds present the in GeO4 unit. High solution enthalpies are observed for other dopants. Solution enthalpies for TiO2 and CeO2 are approximately 5 eV and 6 eV respectively, suggesting that they are highly unlikely to occur.

Conclusions
In this study, atomistic-scale simulations based on classical potentials were employed to study the defect chemistry, Li ion diffusion and doping behaviour in monoclinic and orthorhombic crystal structures of Li2ZnSiO4. The Li-Zn anti-site defect is the dominant defect that would be observed in this material. The Li Frenkel is the second lowest defect energy process, ensuring that Li vacancies and Li interstitials would be present at equilibrium at slightly high temperatures. In both phases, long-range Li-ion diffusion is two dimensional with a moderate diffusion rate and the activation energy of ~ 0.90 eV. Here, we show that Li vacancies required for the vacancy-assisted Li-ion selfdiffusion can be introduced in the lattice by doping Fe 3+ at the Zn site. Doping of Al 3+ at the Si site would introduce extra Li that could enhance the capacity of this material. Promising isovalent

Conclusions
In this study, atomistic-scale simulations based on classical potentials were employed to study the defect chemistry, Li ion diffusion and doping behaviour in monoclinic and orthorhombic crystal structures of Li 2 ZnSiO 4 . The Li-Zn anti-site defect is the dominant defect that would be observed in this material. The Li Frenkel is the second lowest defect energy process, ensuring that Li vacancies and Li interstitials would be present at equilibrium at slightly high temperatures. In both phases, long-range Li-ion diffusion is two dimensional with a moderate diffusion rate and the activation energy of~0.90 eV. Here, we show that Li vacancies required for the vacancy-assisted Li-ion self-diffusion can be introduced in the lattice by doping Fe 3+ at the Zn site. Doping of Al 3+ at the Si site would introduce extra Li that could enhance the capacity of this material. Promising isovalent dopants at the Zn and the Si sites are Fe 2+ and Ge 4+ respectively. Though there is a small difference in calculated defect energies between both phases, the trend was observed to be the same.
Funding: This research was financially supported by the European Union's H2020 Programme under Grant Agreement number 824072-HARVESTORE.