Abstract:Diffusion mechanisms of four intrinsic point defects in rutile TiO2, titanium interstitial (TiI), titanium vacancy (Vti), oxygen interstitial (OI) and oxygen vacancy (VO) are studied in the framework of density functional theory with quantum ESPRESSO suite. Diffusion processes are simulated by defect movement between two stable atomic configurations through using the climbing image nudged elastic band (CI-NEB) method.The initial and final atomic structure in the minimum energy path (MEP) are constructed with 3×3×4 perfect supercell matrix of 216 atoms. Considering that oxygen atoms build up TiO6 octahedron and half of the octahedral centers are occupied by Ti atoms in rutile, interstitial defect is constructed by adding one Ti or O atom to the empty oxygen octahedral center, and vacancy defect is constructed by removing one atom from crystal lattice grid. Structural relaxation is performed before performing the NEB calculation with gamma k point sampling in irreducible Brillouin zone with an energy cutoff of 650 eV. As rutile TiO2 has tetragonal symmetry (P42/mnm), the diffusion channel along the[100] direction is equivalent to the[010] direction. Then, the diffusion paths along the direction parallel to c axis ([001] direction) and perpendicular to the c axis ([100] or[110] direction) are chosen to find the minimum diffusion energy path of TiI and OI. As for VTi and VO, diffusion paths are established from the vacancy site to nearest lattice site of the same kind.Calculation results exhibit significant anisotropy of energy barrier and diffusion mechanism. Of all defect species, TiI diffusion along the[001] direction through interstitial mechanism has the lowest activation barrier of 0.5057 eV. In addition, diffusions along the[100] and[110] direction through kick-out mechanism show higher energy barriers of 1.0024 eV and 2.7758 eV, respectively. Compared with TiI, OI shows small barrier discrepancy between different diffusion directions, which is 0.859 eV along[001] and 0.902 eV along[100] direction. For vacancy defects, diffusion can occur only through the vacancy mechanism. The activation barrier energy of symmetrically inequivalent diffusion path of VO is 0.735 eV along the[110] direction, 1.747 eV along the[001] direction, and 1.119 eV from the TiO6 apex site to the equator site. On the other hand, VTi has two inequivalent paths with much larger diffusion energy barriers:2.375 eV along the[111] direction and 3.232 eV along the[001] direction. In summary, the TiI interstitial diffusion along the[001] direction (parallel to the c axis) has the lowest activation barrier in rutile TiO2, which is in excellent agreement with former experimental and theoretical data.

Threshold displacement energies and defect formation energies in Y2Ti2O7.

Threshold Displacement Energy of Lattice Atoms in Yttrium Ion Garnet, an Ab Initio Molecular Dynamics Study

Ab Initio Study of A Charged Vacancy in Yttrium Aluminum Garnet (Y3al5o12)

Ab initio Molecular Dynamics Simulation of Threshold Displacement Energies and Defect Formation Energies in Y4Zr3O12

An Ab Initio Molecular Dynamics Study on the Threshold Displacement Energies in Yttrium Aluminum Garnet

Impact of Point Defects on Electronic Structure in Y2Ti2O7

The formation and migration of point defects and the thermodynamic behavior of TiD 2 by first-principles calculations

Theoretical Study of Native Point Defects on Anatase TiO 2 (101) Surface

An Accurate Theoretical Study on Intrinsic Defect Energetics in Rutile Tio2

Ab Initio Molecular Dynamics Simulations of Ion–solid Interactions in Gd2Zr2O7 and Gd2Ti2O7

Ab Initio Molecular Dynamics Simulations of Threshold Displacement Energies in Srtio3

Theoretical investigation of structural, energetic and electronic properties of titanate pyrochlores

Ab Initio Study of Defect Properties in Ypo4

Intrinsic Defects in Yttrium Iron Garnet: A First-Principles Study

Chemical Descriptors of Yttria-Stabilized Zirconia at Low Defect Concentration: An ab Initio Study

Antisite Defect in Nonstoichiometric Yttrium Aluminum Garnet: Experimental and First-Principles Calculation

Point Defect Concentrations and Interactions in D019-Ti3Al from First-Principles Calculations

Threshold Displacement Energy In Gan: Ab Initio Molecular Dynamics Study

Ab Initiomolecular Dynamics Simulations of Low-Energy Recoil Events in ThO2, CeO2, and ZrO2

Density functional theory calculation of diffusion mechanism of intrinsic defects in rutile TiO2

Energetics of point defects in yttrium aluminum garnet doped with Mg and Si