Revisiting the structural, electronic and photocatalytic properties of Ti and Zr based perovskites with meta-GGA functionals of DFT

The strongly constrained and appropriately normed (SCAN) functional of density functional theory (DFT) conforms to all possible exact constraints required of a meta-GGA functional and offers good approximations for structural and energetic properties of solids in comparison to experiments. However, SCAN is unable to fully overcome the underestimation of band gap for perovskite oxide materials suitable for photocatalysis. In the present work, we use a combination of meta-GGAs SCAN and modified Becke–Johnson local density approximation (mBJ-LDA) potential functional to accurately compute the structural, energetic, mechanical, vibrational and optoelectronic properties of Ti and Zr based ABO3 (A = Sr, Ba and B = Ti and Zr) perovskite oxides. In addition to evaluating their physical properties, the potential applications of these materials as photocatalyst operating in the UV region of the electromagnetic spectrum are also examined. We show that the structural, energetic, mechanical and vibrational properties calculated using SCAN are in better agreement with experimental data as compared to the commonly used semi-local functionals of DFT. However, the optoelectronic properties of the large band gap Ti and Zr based perovskite oxides are further improved if computed with the mBJ-LDA potential functional, whereby an even higher level of accuracy than with SCAN is achieved, with results that are comparable to the computationally expensive hybrid DFT functionals. On the whole, our DFT calculations indicate that a combination of SCAN and mBJ-LDA functionals for exploring the physical properties of large band gap perovskite oxides provide the means for identifying photocatalysts suitable for hydrogen production at low computational costs.


Introduction
In condensed matter physics and related sciences, Kohn-Sham (KS) density functional theory (DFT) has become a standard ab initio tool for exploring materials. For this reason DFT is being rigorously used for understanding the physical properties of materials and has proven to be particularly instrumental in predicting novel materials. [1][2][3][4][5] Since the precision of KS DFT calculations is strongly dependent on the selected approximation for computing exchange-correlation energy (i.e. E xc = E x + E c ), efforts for developing reliable functionals for determining E xc have seen a steady progress. 3,6,7 By examining self-consistent one-electron equations, one expects that a universal exchangecorrelation energy functional for KS DFT, wholly dependent on the ground-state electron density, r, can be obtained since E xc [r] is expressed by half of the Coulomb interaction in a double integral over space between electron and exchange-correlation hole. However, an approximate functional is always needed in real applications of KS DFT where semi-local approximations have traditionally allowed to formally compute E xc [r] from a single integral through a computationally efficient procedure.
Over the years, many DFT functionals have been developed which can be assigned to one of the four rungs of Jacob's ladder based on their level of sophistication and acquired accuracy. 8 Among these four rungs of Jacob's ladder of DFT functionals, the first three rungs are occupied by semi-local approximations. 8,9 A common feature of all the semi-local DFT functionals is that they can be expressed as single integral, a Computational Materials Modeling Laboratory, Department of Physics, E xc ¼ Ð e xc ðrÞd 3 r, where e xc is the exchange-correlation energy density per volume. For the first and second rungs, respectively, e xc is expressed in terms of electron density, r, alone (as in the case of local density approximation (LDA)) 10 and r and its first gradient, rr, (as in the case of generalized gradient approximation (GGA)). 11 On the other hand, combining electron density and its first derivative with kinetic-energy density and/or r 2 r gives us the so-called meta-GGA functionals, which occupy the third rung of Jacob's ladder. 7,9 Apart from the semi-local functionals, there also exist DFT functionals which allow computing E xc [r] through a double integral or Hartree-Fock (HF) exchange. Although accurate, these so-called hybrid functionals are computationally much more expensive than the semi-local functionals and thus occupy the fourth rung of the Jacob's ladder. 12 Recently, meta-GGA functionals have been increasingly employed to compute physical properties of complex oxide semiconductors. This can mainly be attributed to the fact that earlier theoretical studies clearly show that the structural and energetic properties of these materials in general, 9,13 and optoelectronic properties of large band gap semiconductors in particular, 14 are not accurately predicted by LDA and GGA functionals. The quest for improving accuracy in determining the physical properties of perovskite oxides can partly be ascribed to their potential applications in advance oxidation processes (AOPs) for which various research groups are making efforts in designing new and efficient photocatalysts. 15,16 In this context, perovskite oxides having the chemical composition ABO 3 constitute a potential class of ceramic materials which are candidates for numerous device applications. This can mainly be ascribed to the presence of BO 6 octahedra where the electronic states of the B and O atoms play dominant roles in the conduction and valence band edges, respectively. Consequently, this provides huge opportunities to considerably alter the physical properties of these materials through chemical as well as physical means. 17,18 In particular, strontium titanate (SrTiO 3 ), barium titanate (BaTiO 3 ), strontium zirconate (SrZrO 3 ) and barium zirconate (BaZrO 3 ) are important perovskite oxides in which the 3d and 4d electrons of Ti and Zr, respectively, are bonded with the 2p electrons of oxygen. The d electrons introduced by the presence of transition metals in the titanium and zirconium based perovskite oxides bring about interesting physical properties for these materials which substantiate the exceptional importance of these materials. These materials have therefore been thoroughly studied both experimentally and theoretically in the past for a wide range of applications. In their pristine forms, both the titanates and zirconates of Sr and Ba show an insulating nature under ambient conditions that persists even at very high temperatures despite structural phase transitions. Owing to their chemical and mechanical stabilities at high temperature, these compounds have found applications in hydrogen sensors, fuel cells and electronic devices. [19][20][21] Importantly, there are several reports available in the literature which place these materials among the most suitable candidates for photocatalysis applications. 16 Owing to these facts, it is inevitable that an accurate account of the structural, energetic and optoelectronic properties of ABO 3 (where A = Sr and Ba; B = Ti and Zr) perovskite oxides should be obtained using DFT calculations.
In 2015, Sun et al. 13 proposed the strongly constrained and appropriately normed (SCAN) meta-GGA functional, which conforms to all the known exact constraints for a meta-GGA. These include six constraints each for the exchange and correlation together with 5 constraints for the sum of the two. 13 In SCAN, the orbital kinetic energy densities of each spin are also added along with other ingredients (e.g. local spin density approximation (LSDA) and spin-density gradients) to obtain exchange and correlation energies, while the exchange and correlation potentials can still be computed using the Perdew-Burke-Ernzerhof (PBE) GGA. 11 As a result, structural and energetic properties computed using the SCAN functional show better accuracy for most types of bonds as compared to standard semi-local functionals of DFT. 22,23 In fact, the SCAN functional has proven to be more efficient than hybrid DFT for structural and energetic properties at only a fraction of the computational cost. 13,22,23 Despite these advantages, the underestimation of band gaps and incorrect optical properties for transition metal oxides computed using SCAN still requires one to resort to either semi-empirical methods 23 or hybrid DFT calculations. 24,25 In order to remain truly ab initio and within the regime of meta-GGA functionals of DFT with low computational cost, it is intuitive to compute the optoelectronic properties of large band gap Ti and Zr based perovskite oxides by employing modified Becke-Johnson local density approximation (mBJ-LDA) functional proposed by Tran and Blaha 14 using the ground state structural properties obtained from SCAN functional. Motivated by this, in the present study we use a combination of SCAN and mBJ-LDA meta-GGA functionals for computing the structural, energetic, mechanical, vibrational, electronic and optical properties of ABO 3 (where A = Sr and Ba; B = Ti and Zr) perovskite oxides using the all-electron full-potential linearized augmented plane wave method (FP-LAPW) of DFT. As the compounds under study are potential candidates for photocatalysis, 16 we examine their photocatalytic properties for producing hydrogen through water splitting using the electronic properties computed with SCAN and mBJ-LDA functional. We expect that the results provided in this study would prove helpful in designing and improving the photocatalytic activities of the ABO 3 perovskite oxides in future studies.

Computational methods
In this work, the DFT calculations are performed using the WIEN2k package. 26 For the FP-LAPW method employed in the present work, the core and valence states are separated by muffin-tin spheres such that the electron wave function inside these spheres is augmented plane waves having spherical harmonics, while in the interstitial region the electron wave function is modeled using plane waves. The muffin-tin radii are chosen to be 2.1, 2.3, 1.8, 1.85. and 1.6 a.u. for Sr, Ba, Ti, Zr, and O atoms, respectively, while an energy separation of À6.0 Ry is used between the core and the valence states. With the present choice of muffin-tin radii and energy separation between the core and valence states, the Sr (4s, 4p, 5s), Ba (5s, 5p, 6s), Ti (3s, 3p, 3d, 4s), Zr (4s, 4p, 4d, 5s) and O (2s, 2p) states are treated as valence states. For the interstitial region, the plane wave cutoff R O K max (where R O is the muffin-tin radius of oxygen and K max is the largest K vector of the plane wave expansion) and the Fourier expansion of charge density vector G max are set to be 8.0 and 24 Bohr À1 , respectively. For the case of cubic perovskite unit cell, 84 k-points in the irreducible wedge of the first Brillouin zone (corresponding to a 14 Â 14 Â 14 k-mesh) are used and the energy and charge convergence criteria are set to 10 À5 Ry and 10 À3 e, respectively. During geometry optimization the atomic forces are always minimized to a value less than 1 mRy a.u. À1 using the SCAN functional for obtaining ground state atomic positions.
Since SCAN has previously been shown to perform better in calculating structural and energetic properties of band gap materials as compare to other functionals, 13 it has been used for computing all the ground state geometries and energetic properties in this work. On the other hand, the electronic and optical properties are computed using both SCAN and mBJ-LDA functionals using the SCAN optimized structures, where improved parameters of mBJ-LDA for semiconductors proposed by Koller et al. are utilized. 27 Thus, we combine the good structure from SCAN with the good electronic structure of mBJ-LDA (i.e. SCAN//mBJ-LDA simulations), since the latter is a potential functional and cannot be used for computing structural properties. 14 It is worth pointing out here that the correlation energy and correlation potential parts of mBJ-LDA can be computed using the PBE GGA. Since the SCAN functional uses PBE GGA for the correlation potential, 13,26 one might desire the same parameterization scheme and perform SCAN//mBJ-GGA simulation of the electronic properties. However, we note that the electronic band structure computed using mBJ-GGA yield band gaps that are farther from experiment, compared to mBJ-LDA. This can be ascribed to the fact that parametric optimization 27 for this potential functional has only been carried out for the default exchange-correlation energy obtained from LDA 26 and can lead to alightly larger deviation in band gap values when used for mBJ-GGA calculations.
For the present study, we consider the prototype cubic crystal (space group No. 221, Pm% 3m) structure of SrTiO 3 , BaTiO 3 , SrZrO 3 and BaZrO 3 perovskites and compute the structural properties using the SCAN functional. Since the energetic properties presented in this study also required computation of ground state properties of the oxygen molecule, elemental solids and binary oxides, all the structural parameters for these molecules and solids are also optimized using the SCAN functional. In the case of oxygen molecule, force minimization is performed for computing the optimized bond length. On the other hand, lattice parameters of elemental solids, binary oxides and ternary oxides are computed using self-consistent calculation for a series of fixed volumes of their respective unit cells. Where necessary, the c 0 /a 0 and b 0 /a 0 ratios have also been optimized. The SCAN optimized unit cells of all the ABO 3 compounds are subsequently used for computing the electronic and optical properties. For computing the mechanical properties, we have utilized the Elastic Constants for Cubic Cases (ELAS) program provided within the WIEN2k package. 26 For the computation of the vibrational properties of the ABO 3 compounds (using a 40-atom supercell), the frozen phonon approximation was employed using the phonopy code. 28 For all the calculations involving phonons, a force convergence criteria of 0.01 mRy a.u. À1 was used. In order to eliminate numerical noise in our phonon calculations, the charge density grid has been enlarged by a factor of 2, while symmetrization of the force constants was also taken into account during post processing of the forces with the phonopy code. 29

Structural and energetic properties
Since thermodynamic quantities, such as enthalpy of formation, require one to find the difference between the cumulative total energies of atomic species in their stable reference phases and the total energies of the compounds constituted by these atoms, it is important that the exchange-correlation functional is able to provide good approximation of the ground state energies. For this reason, we first examine the performance of the SCAN functional in determining the energetic properties of all the atoms constituting the ABO 3 (where A = Sr and Ba; B = Ti and Zr) compounds. Table 1 presents the structural properties and cohesive energies of O 2 dimer, fcc Sr, bcc Ba, hcp Ti and hcp Zr computed using the SCAN functional. In the present work, we have computed the cohesive energy, E c , as the difference between the total energy of the above-mentioned stable reference phases of an atom X and the energy of the same atom X placed inside a large supercell, i.e. Table 1 Structural properties (Å) and cohesive energies (eV per atom) of oxygen molecule, face centered cubic unit cell of strontium, body centered cubic unit cell of barium, hexagonal closed packed unit cell of titanium and hexagonal closed packed unit cell of zirconium computed using the SCAN functional. For the sake of comparison, results of the GGA functional reported in previous DFT studies and experimental data are also presented In the case of oxygen, E stable-phase X in eqn (1) is the total energy of a SCAN optimized O 2 dimer, which is simulated in a large asymmetric orthorhombic cell using only the G-point. On the other hand, E stable-phase X values in the case of Sr, Ba, Ti and Zr are computed from their respective SCAN optimized crystal unit cells. E atom X in eqn (1) is the energy of an isolated single atom which has been computed by placing a single atom in an asymmetric orthorhombic unit cell using the G-point. For a correct determination of the E stable-phase X and E atom X energies, spin-polarized calculations are adopted which ensure proper occupation.
It is evident from Tables 1 and 2 that the lattice parameters for elemental as well as binary solids obtained using the SCAN functional are in better agreement with experimental data compared to the GGA functional. [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48] Importantly, the SCAN calculated values of E c for Sr, Ba, Zr and Ti also a show significantly smaller deviation from experimental data which is usually much larger for semi-local DFT functionals. 13 Importantly, it is evident from our calculated values of bond length and cohesive energy for O 2 molecules that the SCAN functional does not suffer from the over-binding issue faced in other semi-local exchange-correlation energy functionals, thus providing an accuracy comparable to hybrid DFT calculations. 13,40,41 On the whole, better agreement of the structural properties and cohesive energies of O 2 and elemental metals computed using SCAN with the experimental data ( Table 1) clearly indicates that this meta-GGA exchange-correlation functional of DFT facilitates accurate determination of the thermodynamic properties of ABO 3 perovskite oxides.
As the primary result of a DFT calculation for solids is the determination of total energy, an important thermodynamic quantity that can be derived from the calculated total energy is the reaction energy, E React , which plays an important role for identifying the application of the subject material in various technologies. 42 For the present study, E React is the energy of perovskite oxides relative to cumulative energies of the binary oxides corresponding to the chemical reaction AO + BO 2 -ABO 3 . For this reason, we can compute E React by subtracting the sum of total energies of its binary oxides AO (E AO t ) and BO 2 E BO 2 t from the total energy of ABO 3 E ABO 3 t according to the relation 42 However, since the total energies of ABO 3 , AO and BO 2 compounds are also related to their respective enthalpies of formation by the formulas and respectively, we can see that E React can also be expressed as In the above equations, E t is the total energy of the optimized unit cells of the concerned oxides. It is evident from eqn (2) to (6) that only a good determination of the ground state total energies of AO, BO 2 and ABO 3 would ensure accurate computation of E React . For the present study, we have computed the ground state lattice parameters and enthalpies of formation of the relevant AO and BO 2 binary oxides using the SCAN meta-GGA functional, which are presented in Table 2. A comparison of the lattice parameters presented in Table 2 clearly shows that the SCAN functional gives a much better account of the structural properties of SrO, BaO, TiO 2 and ZrO 2 compared to the common semi-local functionals belonging to the first two rungs of the Jacob's ladder. 8 At the same time, the calculated enthalpies of the formation for the four AO and BO 2 oxides clearly show that large overestimation commonly seen in the GGA calculations is eliminated by the employment of the SCAN functional. This can mainly be ascribed to the fact that the SCAN functional does not suffer from the overbinding of oxygen bonds, which allows it to perform better for both molecules and solids. 13,23 The lattice parameters of the ABO 3 (where A = Sr and Ba; B = Ti and Zr) compounds under consideration computed using the SCAN meta-GGA functional are presented in Table 3, while the total energy vs volume and pressure vs volume data are shown in Fig. 1. When compared to experimental data, the SCAN calculated lattice parameters of SrTiO 3 and BaZrO 3 are only slightly larger (o0.5%), while the lattice parameters for BaTiO 3 and SrZrO 3 are slightly smaller (o0.12%). It is also  Table 3 that the enthalpies of formation computed using the SCAN functional for ABO 3 compounds are in better agreement with experimental values compared to GGA. These results further confirm the superior performance of this meta-GGA functional in determining the thermodynamics of complex oxides. From Fig. 1, it can also be seen that the pressure vs. volume data computed using the Birch-Murnaghan's equations of state give a linear variation of pressure with volume for the SCAN calculations. Since the mechanical properties of perovskite materials (discussed later) are strongly dependent on smooth variation of pressure vs volume of the unit cell, Fig. 1 already indicates good performance of the SCAN meta-GGA functional. 13 It is also worth pointing out here that the experimental enthalpies of the formation for the ABO 3 compounds are less than the sum of enthalpies of the formation of the concerned binary oxides (i.e. DH AO f þ DH BO 2 f ), which indicates that the structure of the cubic perovskite oxides is energetically stable. 49 Since experimental studies have shown the existence of BaTiO 3 in the cubic phase, 51,52 our calculated values of the enthalpies of the formation of BaTiO 3 , BaO and , showing that the SCAN meta-GGA functional has good agreement with experimental observations which have earlier been satisfied only through the use of the semi-empirically corrected GGA+U method. 42 As mentioned earlier, determination of E React for the chemical reaction AO + BO 2 -ABO 3 is important for determining the suitability of a perovskite oxide for various applications. 42 In addition, E React computed using DFT have also been used to examine the suitability of one-step and two-step techniques for synthesizing novel ternary compounds. 5 In the case of compounds crystallizing in the perovskite structure, one can relate E React with the dimensionless Goldschmit's tolerance factor, t. 56 As the value of Goldschmit's tolerance factor indicates the stability of the perovskite structure, one expects an increase in t with increase in the values of E React . 57 In Fig. 2 we present our SCAN calculated E React values as a function of increasing t of the four perovskite oxides. The comparison of our calculated results with the previous GGA+U calculations 42 and experimental data 58 as shown in Fig. 2 makes it evident that the SCAN meta-GGA functional is efficient in predicting reaction energies. Importantly, the E React values obtained for SrTiO 3 , BaTiO 3 and BaZrO 3 with SCAN are found to be in excellent agreement with experimental ones.

Mechanical properties
Since the mechanical response of a material to applied forces can be derived from its elastic constants, DFT calculations are commonly employed to supplement experimental results 59 or to predict mechanical properties 60 of ABO 3 (where A = Sr and Ba; B = Ti and Zr) perovskite oxides. In the case of cubic crystals the mechanical properties can be completely defined in terms of three independent elastic constants C 11 , C 12 and C 44 , which correspond to longitudinal compression, transverse expansion and shear elastic modulus, respectively. 61 Since the three elastic constants for cubic structure can be computed from straight forward total energy calculations, in the present study, we have used the SCAN meta-GGA functional to compute the elastic constants of the cubic perovskite structure of alkaline earth metal titanates and zirconates, which are presented in Table 4. The comparison of our SCAN results with previous GGA 59 calculations for SrTiO 3 reveals that both these functionals show small deviations from experimental data. 62 However, larger deviations of both SCAN and GGA results from the experimental data are evident in the case of BaTiO 3 . 51 Since the large deviation of our calculated values of elastic constants for BaTiO 3 can be due to the fact that this compound generally adopts a distorted perovskite structure, it is difficult to judge the performance of SCAN functional based solely on the comparison between the data presented in Table 4.
Although various important mechanical parameters (e.g. shear modulus, Young's modules, Zener anisotropy factor, and Kleinman parameter) can be derived from C 11 , C 12 and C 44 using simple mathematical formulas, 61 we note that only experimental data of elastic constants for SrTiO 3 and BaTiO 3 are available in the literature. 51,62 However, Yamanak et al. 63,64 have carried out studies for measuring the thermophysical properties of SrZrO 3 and BaZrO 3 in which they report the shear modulus, G (GPa), and Young's modulus, Y (GPa), of these perovskite oxides. For this reason, we first compute the bulk modulus, B (GPa), Voigt's shear modulus, S V (GPa), and Reuss's shear modulus, S R (GPa) for the four ABO 3 compounds using the formulas and respectively. The values computed using eqn (7) through (9) are subsequently used for computing G and Y from and respectively. On comparing our calculated values of G and Y with GGA and experimental data presented in Fig. 3, it is clear that the SCAN meta-GGA functional performs better than the semi-local GGA functionals for predicting the mechanical properties of perovskites considered in this work.

Vibrational properties
Vibrational properties of solids can be evaluated in terms of phonons that represent the collective dynamic excitations in a periodic lattice. Apart from being significantly important in indicating phase stability and structural transitions, phonons are crucial in modulating the transport properties of charge carriers in bulk [65][66][67] and nano-materials. 68,69 Since lattice vibrations can provide additional momentum to charge carriers of indirect band gap materials, 70 phonons also play a key role in interband transitions. For this reason, a reliable knowledge of the phonon modes in a perovskite photocatalyst can prove useful in designing improved absorption of electromagnetic radiations 71,72 as well as controlling the electron-hole recombination rates. 73,74 For the five atom cubic unit cell of the perovskites considered in this work, there are five triply degenerate phonons at the Brillouin zone center (G-point).
Among these phonons, one of the triply degenerate phonons is acoustic while the remaining four are optical. The coherent movement of charges in the acoustic phonons gives them zero frequency at the G-point, while the four triply degenerate optical phonons in a cubic perovskite structure have non-zero frequencies at the G-point. 75 Three of these four triply degenerate optical phonons are infrared active and are designated as  75 On the other hand, the fourth triply degenerate (designated as G 25 ) phonon mode is optically inactive and corresponds to a B-O 3 torsion mode. 75 It is well know that the dielectric and ferroelectric properties of the cubic ABO 3 perovskite oxides are controlled by phonons, 51,53,55,76 which can be modified by varying the stress or temperature. It is therefore intuitive to assess the performance of a DFT functional in predicting experimentally observed phonon modes and instabilities (e.g. ferroelectric at low temperatures). Unfortunately, a correct description of the vibrational properties from the frozen phonon approximation using semi-local as well as advanced functionals of DFT at 0 K has remained a challenge. 50,[77][78][79][80][81] For instance, plane wave (PW) based DFT calculations by Wahl et al. 50 utilizing both semi-local and Heyd-Scuseria-Ernzerhof (HSE06) schemes only confirm the experimentally observed Sr-TiO 3 ferroelectric soft mode of SrTiO 3 to be real within the LDA approximation (see Table 5). On the contrary, all of LDA, GGA and HSE06 confirm the imaginary Ba-TiO 3 ferroelectric soft mode of BaTiO 3 . 50 However, recent employment of the SCAN functional for computing phonon modes of SrTiO 3 and BaTiO 3 at the G-point using the projector augmented wave (PAW) method has shown better performance of this meta-GGA exchange-correlation parameterization scheme. 79 Since the PW/PAW based DFT calculations of phonons for the considered class of compounds are found to be sensitive to the choice of pseudopotentials as well as the treatment of semi-core states as valence or core states, 79,80,82 an accurate all-electron account of vibrational frequencies at the G-point through FP-LAPW method of DFT using the SCAN meta-GGA functional is therefore mandated.
Since performance evaluation of various DFT functionals in predicting the vibrational properties of ABO 3 perovskites can be carried out on the basis of optical phonon frequencies at the G symmetry point, 50,[75][76][77][78] in Table 5 we present our SCAN results together with the available earlier theoretical and experimental values. It is worth pointing out here that unlike the usual space group dependent unit cell calculations performed within the WIEN2k code for structural and electronic properties, the phonopy code requires forces caused by finite atomic displacements to be computed within supercells which are deprived of space group. 28 For this reason, in this work, the atomic forces for the 40-atom supercells have been computed using G-point sampling. From Table 5 it is evident that the G-point phonon frequencies computed using SCAN in the present FP-LAPW calculations are in excellent agreement with the available experimental data. [83][84][85] One can clearly see that both SrTiO 3 and BaZrO 3 show no tendency for ferroelectric distortion at low temperatures. On the other hand, both BaTiO 3 and SaZrO 3 show a tendency to ferroelectric distortions at low temperatures, which has also been observed in earlier studies. 50,82,84 It is particularly important to notice that our calculated values for SrTiO 3 and BaTiO 3 are in better agreement with the experimental data compared to HSE06. 50 Moreover, our FP-LAPW results of the A-BO 3 ferroelectric soft mode for SrTiO 3 and BaZrO 3 computed using SCAN are the closest to the experimental data when compared with earlier theoretical results. 79,82,83,85 The better performance of SCAN can mainly be attributed to the fact that an accurate determination of the vibrational properties using the frozen phonon approximation requires very precise computation of the curvature of the total energy surface with respect to the atomic displacements introduced in the solid. As SCAN already shows good performance in determining the total energy dependent physical properties of atoms, molecules and solids, 13 from Table 5 it is evident that it also performs better than the computationally expensive HSE06 functional in determining the vibrational properties of ABO 3 perovskites. 50,79

Electronic properties
Using the lattice parameters optimized with the SCAN functional, we have computed the electronic density of states (DOS) and band structure diagrams of the selected alkaline earth metal titanates and zirconates, which are presented in Fig. 4 and 5, respectively. For these properties both SCAN and mBJ-LDA functionals are employed since the former DFT functional is known to underestimate the band gaps of solids despite its good performance for structural and energetic properties. 13,23 By comparing Fig. 4 and 5, one can see that the overall features of the electronic properties for SrTiO 3 , BaTiO 3 , SrZrO 3 and BaZrO 3 computed using the two meta-GGA functionals are similar. However, both the DOS and band structure diagrams show that the positioning of the electronic states belonging to conduction and valence bands computed using SCAN//mBJ-LDA simulations are modified, which results in increase of the band gaps. From Fig. 4 we can see that the upper valence bands for SrTiO 3 and BaTiO 3 are mainly composed of the O-2p states along with a minor contribution coming from the Ti-3d states. In a similar fashion, Fig. 5 shows that the upper valence bands for SrZrO 3 and BaZrO 3 are mainly composed of the O-2p states along with a minor contribution coming from the Zr-4d states. Irrespective of the choice of the functional, it is also evident that for all the ABO 3 perovskites considered in this study the conduction band edge is mainly composed of the B atom's d states along with a minor contribution coming from the O-2p states. Moreover, the contribution of Sr-5s and Ba-6s states in the conduction band minimum is negligible; however, these states dominate the conduction band at higher energies. This is understandable since the ABO 3 perovskite oxides are composed of alternating layers of AO and BO 2 , such that stronger covalent bonding between B and O atoms dominates the conduction and valence band edges, respectively. 32,35 From the electronic properties of ABO 3 (where A = Sr and Ba; B = Ti and Zr) shown in Fig. 4 and 5, it is evident that the fundamental energy band gap between the O-2p states at the valence band edge and the B atom's d states at the conduction band edge is along the R -G direction. On the other hand, the optical band gap for all these compounds is slightly larger than the fundamental band gap, which is along the G -G direction. The SCAN and mBJ-LDA calculated band gaps for ABO 3 (where A = Sr and Ba; B = Ti and Zr) along the R -G and G -G directions are presented in Table 6 and compared with earlier theoretical results computed using GGA and HSE06 functionals as well as with experimental data. [86][87][88][89][90][91] It is clear that the SCAN meta-GGA functional employed in this study underestimates band gaps for the perovskite oxides considered here, which are only slightly larger than the band gaps predicted by the GGA functional. 23,86,92 The comparison of the fundamental band gaps computed using SCAN//mBJ-LDA simulations and the ones computed using hybrid DFT in ref. 86 shows that a relatively smaller shift of the B atom's d states making up the conduction band edge at the G symmetry point is achieved with the meta-GGA (i.e. mBJ-LDA) functional employed here. Since the band structures shown in Fig. 4 and 5 reveal that the application of mBJ-LDA does not shift the O-2p states located at the R symmetry point, our calculated mBJ-LDA R -G band Table 5 The G-point optical phonon frequencies (THz) for the cubic phase of ABO 3 (where A = Sr and Ba; B = Ti and Zr) compounds computed using the SCAN meta-GGA functional in the present study compared with the available theoretical results computed using the SCAN, LDA, GGA and HSE06 functionals and experimental data. Imaginary frequencies are indicated by imaginary unit, i, after the value gaps are smaller than the ones obtained using the HSE06 functional. 86 On the other hand, Fig. 4 and 5 clearly show that the O-2p states located at the G symmetry point below the Fermi level are pushed towards more negative energies with the application of the mBJ-LDA functional. As a consequence, the optical band gaps of SrTiO 3 and SrZrO 3 computed using the SCAN//mBJ-LDA simulations (Table 6) are in better agreement with the experimental data. [87][88][89][90][91] On the other hand, optical band gaps of cubic BaTiO 3 and BaZrO 3 are slightly smaller when compared with the hybrid DFT results. 86 Despite the fact that different experimental values for the optical band gaps of BaZrO 3 are reported in the literature, 90,91 comparison of our calculated band gap with earlier works suggests that the performance of the mBJ-LDA functional for barium compounds is also satisfactory. In particular, the small difference in the fundamental and optical band gap of cubic BaTiO 3 (0.11 eV) is in good agreement with the previous experimental reports which reported difficulty in determining the direct or indirect nature of the band gap for this material. 88 Thus, our results allow us to claim that the SCAN//mBJ-LDA simulations are a better choice than HSE06 for reliably evaluating the optical properties and photocatalytic performance of ABO 3 (where A = Sr and Ba; B = Ti and Zr) perovskite oxides.

Optical properties
For practical purposes, the application of a material in optical devices can be directly related to its absorption properties. Since the optical band gaps of ABO 3 (where A = Sr and Ba; B = Ti and Zr) ceramics make them suitable for applications in devices operating in the 220 nm-360 nm wavelength range, it is intuitive to examine their optical properties. As a long range of absorption appears to be plausible for transition between states having varying degrees of spatial overlaps, the band gap alone does not give complete information of absorption properties.
To address this issue, optical absorption in a wide range of wavelength has been obtained for SrTiO 3 , BaTiO 3 , SrZrO 3 and BaZrO 3 perovskite oxides by computing the real and imaginary parts of the complex dielectric function. 93,94 The optical absorption plots calculated using the SCAN and mBJ-LDA meta-GGA functionals are presented in Fig. 6. Irrespective of the choice of the functional used, it is evident that the optical absorption is maximum in the far-UV region (o200 nm) of the electromagnetic spectrum. Since the SCAN functional is unable to accurately reproduce the experimental band gaps of the compounds under study, one can see that the threshold of the optical absorption for this functional is in the visible region for the titanates and UV region for the zirconates. However, the results computed using SCAN//mBJ-LDA simulations show that the threshold of optical absorption for all the compounds under consideration is in the UV region, in agreement with the experimental observations. 87-91

Photocatalytic performance
Despite the fact that the optical absorption plots shown in Fig. 6 indicate that SrTiO 3 , BaTiO 3 , SrZrO 3 and BaZrO 3 are suitable for applications in optical devices operating in the UV region of the electromagnetic spectrum, the application of these perovskite oxides for hydrogen production through photocatalysis requires that the band edge potentials are appropriate with reference to water oxidation and reduction potentials. 16 It is well known that for successfully producing hydrogen by water splitting, a material's conduction band minima (CBM) should be more negative than the reduction potential of H + /H 2 , while the valence band maxima (VBM) should be more positive than the oxidation potential of H 2 O/O 2 . 16 In order to estimate positions of the conduction and valence band edges of the compounds under study, we have computed band edge potential using the procedure detailed in ref. [95][96][97][98] for the SCAN, mBJ-LDA and experimental optical band gaps presented in Table 6. The calculated results are presented in Fig. 7, where one can see that the edge potential of CBM computed using   Fig. 7 again makes it evident that the mBJ-LDA functional is able to provide better agreement with experiment as compared to the SCAN functional. An important feature evident from our calculated mBJ-LDA results and experimental data presented in Fig. 7 is that the Ti based compounds have both their CBM and VBM edge potentials located close to the H + /H 2 and H 2 O/O 2 potentials, respectively, when compared to the Zr based compounds. Since reports from 1980s already show insignificant photocatalytic activity of pristine SrTiO 3 for splitting water, 99,100 one might suspect that pristine perovskite oxides with band gap 4 4 eV (e.g. LiTaO 3 101 and BaZrO 3 102 ) are more suited for hydrogen production. However, this is not the case, because SrTiO 3 loaded with a co-catalyst show an improved photocatalytic activity. This is due to the lower rate of recombination of the photogenerated electron and holes achieved in Rh loaded SrTiO 3 , which results in better hydrogen evolution (628 mmol h À1 ) 103 when compared to pristine BaZrO 3 (104.5 mmol h À1 ). 102 Therefore, recombination is an important factor influencing the photocatalytic performance of a material 104 that can be evaluated in terms of the transfer rate of photogenerated electrons and holes in the reduction and oxidation processes. 105 Unfortunately, the transfer rate of photogenerated electrons and  holes cannot be directly measured from DFT calculations, however, it can indirectly be evaluated from the absolute ratio of the effective masses of the photogenerated holes and electrons in valence and conduction bands, respectively, according to the formula 105 A large value of D would indicate a better transfer rate of photogenerated electrons and holes, hence lower rate of recombination. On the other hand, smaller values of D would indicate higher rate of recombination. In eqn (12), E k is the energy corresponding to wavevector k. Since the fundamental band gaps of Ti and Zr based perovskite oxides considered in this work are along the R -G symmetry points, we have considered the curvatures of the bands at G point for computing the effective mass of electron and the curvatures of the bands at R point for computing the effective mass of holes. 105 A comparison of the results computed using SCAN and mBJ-LDA (Table 7) clearly indicates that effective masses obtained from the former meta-GGA functional are slightly smaller; indicating that SCAN gives a larger curvature for both the conduction and valence band edges. This inability of the SCAN functional in predicting the accurate band structures of materials having band gaps 4 1 eV is analogous to the results obtained with a GGA functional 23 which also results in an underestimation of the effective masses. 106 On the contrary, application of mBJ-LDA causes narrowing of the band dispersion that results in large effective masses, which are overestimated as compared to experimental data. 106,107 Nevertheless, the values of D computed using both the SCAN and mBJ-LDA functionals clearly show that an improved separation of charges should be expected from Zr based perovskite oxides. These results are in accordance with the photocatalytic hydrogen evolution under UV irradiation of pristine SrTiO 3 99,100 and BaZrO 3 . 102 We therefore conclude that the employment of SCAN//mBJ-LDA simulations for investigating oxide semiconductors using DFT allows a reliable account of the structural, electronic and photocatalytic properties which are usually possible through computationally expensive hybrid functionals or GW methods. 13,14 Moreover, our results clearly suggest that the employment of a combination of SCAN and mBJ-LDA meta-GGA functionals for examining the modulation of band edge potentials of large band gap ceramics by means of doping can be important in predicting new photocatalysts.

Conclusions
In conclusions, we have performed a systematic analysis of combined application of SCAN and mBJ-LDA meta-GGA functionals for computing the physical properties of alkaline earth metal titanate and zirconate in cubic perovskite structure using the FP-LAPW method of DFT. The potential of these perovskite oxides for photocatalytic hydrogen production by splitting water is also examined. Our results clearly show that the experimentally observed structural, energetic, mechanical and vibrational properties of large band gap perovskite oxides are accurately predicted by the SCAN functional. In particular, the accurate prediction of the structural and energetic properties of atomic species and the binary oxides involved in synthesis of the concerned perovskite oxides shows that SCAN meta-GGA functional is able to provide an accuracy comparable to that of hybrid DFT calculations. The comparison of SCAN calculated reaction energies for SrTiO 3 , BaTiO 3 , SrZrO 3 and BaZrO 3 shows good agreement with experimental observations, which have previously been achieved through GGA+U calculations. By computing the elastic constants and frequencies of the optical phonons at the Brillouin zone center of the concerned perovskite oxides, we also show that SCAN functional performs better than the semi-local GGA as well as HSE06 functional for accurately predicting the mechanical and vibrational properties of these materials. Despite its good performance for the above-mentioned physical properties of perovskite oxides, SCAN functional is found to fall short of reproducing the experimentally observed band gaps of alkaline earth metal titanates and zirconates. To address this issue, we show that the application of SCAN//mBJ-LDA simulations allows one to achieve better accuracy in determining the optoelectronic properties of SrTiO 3 , BaTiO 3 , SrZrO 3 and BaZrO 3 . In agreement with experiments, all four compounds are found to be indirect band gap materials showing absorption thresholds in the UV region of the electromagnetic spectrum. Using the electronic band gaps predicted by mBJ-LDA meta-GGA functional we also show that our ab initio predicted band edge potentials and effective masses of photogenerated electrons and holes are consistent with photocatalytic hydrogen evolution observed in SrTiO 3 and BaZrO 3 . Since ab initio calculations provide the means to simulate changes in material's properties, our results clearly demonstrate that combining SCAN and mBJ-LDA functionals of DFT would prove useful as a computationally economical and robust technique for designing new photocatalysts based on Ti and Zr based perovskite oxides. This can particularly prove useful in future DFT studies for exploring the

Conflicts of interest
There are no conflicts to declare.