First-principles investigation of thermodynamics and electronic transitions in vacancy-ordered rare-earth perovskite nickelates

Abstract

Controlled introduction of oxygen vacancies offers an effective route to induce metal-to-insulator transition in strongly correlated rare-earth nickelates ($R$NiO$_3$) at room temperature. However, the role played by the rare-earth cations on the structure, thermodynamic stability, and electronic properties of oxygen-deficient nickelates remains unclear. Here, we employ density functional theory calculations with Hubbard corrections (DFT + $U$) to investigate the whole family of $R$NiO$_{2.5}$ ($R$ = Pr-Er) compounds in two commonly observed oxygen-vacancy ordered configurations, namely brownmillerite, and square planar. We find that square planar polymorph is always more stable ($\sim$0.4 eV/u.f) than the brownmillerite for all rare-earth cations, owing to the exceedingly low volumetric strains (< 1\%). Formation energy of $R$NiO$_{2.5}$ gradually increases with decreasing size of $R$ owing to stronger Ni-O covalent interactions in pristine $R$NiO$_3$ with small $R^{3+}$ cations. This necessitates more oxygen-lean environments for synthesis of $R$NiO$_{2.5}$ with smaller $R^{3+}$ cations. Analysis of the density of states and band structures reveals that electronic structure of $R$NiO$_{2.5}$ is governed by two factors: (a) localization of electron on NiO$_6$ octahedra yielding a Mott insulating state with strong correlations as Ni $e_g$ is half filled, and (b) crystal field splitting in the NiO$_4$ tetrahedra/square planar polyhedra. Brownmillerite $R$NiO$_{2.5}$ is metallic, while square planar $R$NiO$_{2.5}$ is an insulator with a predicted gap of $\sim$ 0.2-0.3 eV, depending on the $R^{3+}$ cation. Crystal orbital Hamilton population (COHP) analysis indicates that the Ni-O bond belonging to square-planar NiO$_4$ polyhedra exhibit much greater covalent character than those in NiO$_6$ octahedra in square planar $R$NiO$_{2.5}$.