Lithium aluminate, LiAlO 2 , is a material that is presently being considered as a tritium breeder material in fusion reactors and coating material in Li-conducting electrodes. Here, we employ atomistic simulation techniques to show that the lowest energy intrinsic defect process is the cation anti-site defect (1.10 eV per defect). This was followed closely by the lithium Frenkel defect (1.44 eV per defect), which ensures a high lithium content in the material and inclination for lithium diffusion from formation of vacancies. Li self-diffusion is three dimensional and exhibits a curved pathway with a migration barrier of 0.53 eV. We considered a variety of dopants with charges +1 (Na, K and Rb), +2 (Mg, Ca, Sr and Ba), +3 (Ga, Fe, Co, Ni, Mn, Sc, Y and La) and +4 (Si, Ge, Ti, Zr and Ce) on the Al site. Dopants Mg 2+ and Ge 4+ can facilitate the formation of Li interstitials and Li vacancies, respectively. Trivalent dopants Fe 3+ , Ni 3+ and Mn 3+ prefer to occupy the Al site with exoergic solution energies meaning that they are candidate dopants for the synthesis of Li (Al, M) O 2 (M = Fe, Ni and Mn) compounds.Energies 2019, 12, 2895 2 of 10
Computational MethodsThe present computational study was based on classical pair-wise potential calculations to describe LiAlO 2 via the General Utility Lattice Program (GULP) code [23,24]. In the present approach, the total energy (lattice energy) is determined by long range (i.e., Coulombic) and short range [i.e., electron-electron repulsive and attractive intermolecular forces (van der Waals forces)]. The latter were modelled using Buckingham potentials (refer to the supplementary information). The van der Waals forces arising from the spontaneous formation of instantaneous dipoles are very important as the formation energy results are sensitive to those forces. The present modelling approach takes into the account of Van der Waals forces as a function of the interatomic distance (r) [23,24]. Two-body Buckingham potential mentioned in the supplementary information consists of two parts. The first part of the equation represents the Pauli repulsion (electron-electron) and the second part denotes the van der Waals interaction. Ionic positions and lattice parameters were relaxed using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm [24]. Convergence criteria dictated that in relaxed configurations, forces on each atom was <0.001 eV/Å. To introduce point defects in the lattice we used the Mott-Littleton methodology [25] similarly to recent work [26][27][28]. It is established that although the present methodology may overestimate the defect formation enthalpies at the dilute limit, the trends will be unchanged [29]. In a thermodynamic perspective defect parameters can be described by comparing the real (i.e., defective) crystal to an isobaric or isochoric ideal (i.e., non-defective) crystal. Defect formation parameters can be interconnected through thermodynamic relations [30,31], with the present atomistic simulations corresponding to the isobaric parameters for the mig...