Theoretical Investigation of Ta2O5, TaON, and ... - ACS Publications

Jan 17, 2017 - in the GW approximation can overcome the severe underestimation of the .... model combined with the GW corrections for absolute band...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIVERSITY OF SOUTH CAROLINA LIBRARIES

Article 2

5

3

5

Theoretical Investigation of TaO, TaON and TaN: Electronic Band Structures and Absolute Band Edges Zhi-Hao Cui, and Hong Jiang J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b12370 • Publication Date (Web): 17 Jan 2017 Downloaded from http://pubs.acs.org on January 30, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Theoretical Investigation of Ta2O5, TaON and Ta3N5: Electronic Band Structures and Absolute Band Edges Zhi-Hao Cui and Hong Jiang∗ Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Rare Earth Material Chemistry and Application, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

January 15, 2017 Abstract Early transition metal oxides, nitrides and oxynitrides have attracted a lot of interest due to their potential application in photovoltaics and photocatalysis. In this work, a systematic investigation is conducted to the electronic band structure of Ta2 O5 polymorphs, β-Ta3 N5 and β-TaON, which are crucial for the understanding of their photocatalytic properties, based on the state-of-the-art first-principles approaches. The calculated results imply that many-body perturbation theory in the GW approximation can well overcome the severe underestimation of the band gap caused by standard density functional theory (DFT) in the local and semilocal approximations, and lead to a quantitative agreement with experiment. The effects of the electron-phonon coupling on electronic band structure are considered by the Fr¨ olich model, and especially for ǫ-Ta2 O5 , a strong electron-phonon coupling is predicted as a result of small highfrequency dielectric constant and large effective masses. Based on the analysis in terms of the phenomenological ionic model, the band gap difference between three compounds can be physically attributed to not only the well-known energy difference between O-2p and N-2p orbitals, but also the influences of the Madelung potential on the conduction band energy. By comparing the calculated absolute band edge positions to the redox potentials for water reduction and oxidation, all of the three compounds are predicted to have the potential photocatalytic property for unassisted water splitting. In addition, we have also analyzed the stability and the band gaps of different Ta2 O5 polymorphs and found that β-Ta2 O5 , the phase commonly used in theoretical studies, is actually unstable and its unusually small band gap can be attributed to strong overlapping of neighboring atomic orbitals. On the other hand, ǫ-Ta2 O5 , which is much less well studied compared to β-Ta2 O5 , leads to calculated properties that are much more consistent with experimental findings of Ta2 O5 in general. Theoretical analysis and findings presented in this work can have general implications for the understanding of electronic band structure of other early transition metal compounds.

2

ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Since the discovery of photocatalytic water decomposition with the titanium dioxide electrode in 1972, 1 a lot of efforts have been made to find the potential semi-conductive materials for photocatalytic and photoelectrochemical applications. 2 Among all of the properties that are relevant to the photo-splitting of water, the band gap (Eg ) and the absolute band edge position (EVBM and ECBM ) with respect to the vacuum level are two crucial ones that determine whether a material is promising for unassisted water splitting. In particular, the band gap should be suitable for visible light absorption, and furthermore, ECBM should be more positive than the redox potential for water reduction (H2 /H2 O, -4.4 eV with respect to the vacuum level) while EVBM needs to be more negative than the redox potential for water oxidation (O2 /H2 O, -5.6 eV). 3–6 In the past decades, various transition metal oxide materials have been explored in order to satisfy the above two conditions. 7 However, only in the recent years, non-oxide materials, especially metal nitrides and oxynitrides, have attracted a lot of interest experimentally. 8–13 To fully realize the potential of these materials in practical applications, it is of great importance to establish a theoretical understanding of their electronic band structure, which is also necessary for rational design of new materials. First-principles electronic structure theory, at present mainly density functional theory (DFT) in the local density approximation (LDA) or various generalized gradient approximations (GGAs), has become an indispensable working tool in the interpretation and prediction of material properties. 14–16 In spite of its extensive applications, however, DFT in LDA/GGA actually faces severe problems for the theoretical description of the electronic band structure, including, in particular, the famous band gap problem, i.e. the systematically underestimation of band gaps of typical semiconductors by 30-100% . 17 The origin of this problem derives from the fact that even with the exact exchange-correlation functional (Exc [ρ]), the Kohn-Sham (KS) band gap εKS g , i.e. the energy difference between highest occupied and lowest unoccupied Kohn-Sham orbitals, often termed as the valence band maximum (VBM) and conduction band minimum (CBM) for solids, respectively, does 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

not correspond to the fundamental band gap Eg , due to the ignorance of the derivative discontinuity of Exc [ρ] (∆xc ). 18–20 Similar difficulties also exist for the description of the absolute band positions. 21 Currently many-body perturbation theory in the GW approximation (GWA) 17,22–24 is the state-of-the-art first-principles approach to describe the electronic band structures of semiconductors, 17,24–26 and therefore is promising for its application to solar energy conversion materials. 4,5,27–33 Moreover, the previous studies have shown the importance of GW correction to the accurate description of ionization potential with the uncertainty comparable with experiments. 21,28,34–37 In this work, we will mainly focus on three tantalum-based compounds, Ta2 O5 ,TaON and Ta3 N5 . The latter two compounds have been verified experimentally to have smaller band gaps with respect to common oxides, 3 and exhibit visible-light photocatalytic activity. 38–40 A series of studies on electronic structure of the three compounds based on DFT have been published. 41–49 Compared with pure metal oxides and nitrides, metal oxynitrides have the attractive feature that it can often combine the small band gap of nitrides with the good stability of oxides. 8–13,50,51 By considering systematically the oxide, nitride and oxynitride of the same transition metal element (i.e. Ta), our purpose in this work is to investigate the band gaps and the absolute band positions of these compounds using the state-of-the-art first-principles approaches, and further clarify the origin of the differences between them. The paper is organized as followings. The next section briefly introduces the methods to solve the band gap problem, including the hybrid functionals and the GW approach. The slab model combined with the GW corrections for absolute band position calculations is presented as well. In the third section, we present a systematic investigation of electronic band structure of Ta2 O5 , TaON and Ta3 N5 . The lattice properties after relaxation are presented and are compared with experimental studies. Then the results of band structures, absolute band positions of different methods (including electron-phonon coupling correction) are collected and the band gap difference is analyzed by the phenomenological ionic model. In addition, we also discuss the polymorphism of Ta2 O5 , and an analysis is made in order to

4

ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

interpret the gap difference of the multiple phases. The fourth section summarizes the main findings of this work.

Theory and Method The band gap problem and GW approximation As we have mentioned in the preceding section, conventional local or semi-local density functional approximations (LDA/GGA) are severely limited in accurate description of electronic band structure of materials, since the fundamental band gap does not correspond to the Kohn-Sham band gap, 18,52 Eg ≡ I − A = εKS g + ∆xc .

(1)

For the exchange-correlation energy in LDA/GGAs, which are explicit functionals of electron density, the derivative discontinuity ∆xc vanishes, which is the origin of band gap problem. In order to solve the problem, one approach is to mix a fraction of the Hartree-Fock exchange with LDA/GGAs exchange-correlation energy in the so-called hybrid functionals, which can recover the missing ∆xc to some extent. One typical hybrid functional is PBE0. 53 However, due to the long-range nature of the HF-type exchange, the computational cost of PBE0 is much larger than standard LDA/GGAs, especially for extended systems. 54 One way to improve the numerical efficiency is to use the screened hybrid functional, e.g. the HSE functional, 54,55 which has been proven to be able to describe the band gaps of many semiconducting materials very well. 56–58 Unlike density functional theory, many-body perturbation theory based on Green’s function provides a theoretically more rigorous framework for the treatment of quasi-particle excitations, as Green’s function itself describes the process of adding an electron to or removing an electron from the system intrinsically. 22,59 The electronic band structure can be

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 36

obtained in a theoretically rigorous way by solving the quasi-particle equation, 17,26 1 [− ∇2 + Vext (x) + VH (x)]Ψn (x) + 2

Z

dx′ Σxc (x, x′ ; En )Ψn (x′ ) = En Ψn (x).

(2)

In Eq. 2, atomic unites are used, Vext (x) and VH (x) denote the external and Hartree potentials, and Σxc (x, x′ ; En ) denotes the exchange-correlation self-energy that accounts for all non-classical electron-electron interactions beyond the Hartree approximation. 17,23,26 In order to solve the quasi-particle energy En , Hedin derived a set of integrodifferential equations 22,23 and introduced the so-called GW approximation, in which Σxc in the space and time domain is approximated as the product of one-body Green’s function G and the screened Coulomb interaction W . In practice, the calculations can be further simplified by the G0 W0 approximation, 17 i Σxc (x, x ; ω) = 2π ′

Z



G0 (x, x′ ; ω ′ + ω)W0 (x′ , x; ω ′)eiηω dω ′

(3)

in which both G and W are calculated from Kohn-Sham orbitals in a given density functional approximation (often LDA or GGA), hence denoted as G0 and W0 , respectively. The quasiparticle energies are calculated by the first order perturbation theory based on Kohn-Sham orbital energies (ǫnk ) and wavefunctions(ψnk ),

Enk = ǫnk + Znk (ǫnk )R hψnk | Σxc (ǫnk ) − Vxc |ψnk i ,

(4)

where Znk is the quasi-particle renormalization factor, and Vxc is the Kohn-Sham exchangecorrelation potential.

Absolute band position with GW correction Absolute valence and conduction band edge positions with respect to the vacuum level, denoted as EVBM and ECBM , respectively, can be obtained from the slab model calculations. 21,60 6

ACS Paragon Plus Environment

Page 7 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In particular, EVBM can be determined by calculating the difference between the highest occupied KS orbital and the vacuum level in a slab model calculation, (s)

(s) EVBM,s = εVBM − Vvac .

(5)

Through the paper, we use “(s)” and “(b)” to denote the quantity calculated from the slab and the bulk system, respectively. However, in the previous work, 21 we have clarified the importance of considering the bulk correction to eliminate the finite-size effects and the possible influences of surface states in the direct slab model calculations (eq. 5), which leads to (b)

(b)

(s)

(s) EVBM,b = [εVBM − εref ] − [Vvac − εref ],

(6)

where εref is the macroscopically averaged electrostatic potential in the center of the slab (s)

(b)

(εref ) or in the bulk system (εref ). 21,61 As shown in our previous study, 21 when using LDA or GGAs, the ionization potential (I = −EVBM,b) calculated in terms of Eq. 6 is systematically underestimated with respect to experiment for typical semiconductors as a result of missing the contribution from the derivative discontinuity, and it is necessary to consider the GW quasi-particle correction, which, at the G0 W0 level, can be straightforwardly obtained by considering the GW correction to (GW )

the VBM of the corresponding bulk system ∆εVBM , 21 (GW )

(KS)

(GW )

EVBM = EVBM + ∆εVBM .

(7)

The energy of the conduction band edge with respect to the vacuum level (ECBM ) can then be easily calculated by adding the GW band gap of the bulk system to EVBM , (GW )

(GW )

ECBM = EVBM + Eg(GW ) .

7

ACS Paragon Plus Environment

(8)

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 36

Computational details All DFT and G0 W0 calculations are performed by using the Vienna ab-initio simulation package (VASP) 62,63 with the projector augmented wave (PAW) method. 64 For DFT calculations, both GGA and the hybrid functional are considered. The GGA functional in the PBEsol 65 parameterization is used for structural optimization since its good performance in the description of structural parameters and surface properties. 65 And the hybrid functional in the HSE06 form 54,55 is used for electronic band structure calculations to partly overcome the underestimation of band gaps by LDA/GGA functionals. 56,57 The crystal structures of all systems concerned are relaxed by using a Γ-centered k-mesh determined by requiring the number of k points × the length of lattice vectors ≈ 30 ˚ A and a large kinetic energy cutoff of 550 eV for plane waves. The latter is necessary to eliminate the possible influence of the Pulay stress on the accuracy of optimized lattice constants. 66 The energy convergence A. After the relaxation, criterion is 10−6 eV and the force convergence criterion is 10−4 eV/ ˚ the electronic structures are calculated with a kinetic energy cutoff of 450 eV. Ab initio molecular dynamics (AIMD) simulations within the canonical ensemble at 50 K, 200 K, 500 K and 1000 K are performed to check the stability of β-Ta2 O5 , in which the time step is chosen to be 0.25 fs, and the total simulation time is about 10 ps. For G0 W0 calculations in this work, we mainly use the KS orbital from PBEsol as the input, hence labeled as G0 W0 @PBEsol or simply G0 W0 for abbreviation. The band gap convergences with respect to the number of unoccupied states and k points have been carefully tested. The band structure diagrams of HSE06 and G0 W0 along high-symmetry k-path are obtained by using the Wannier interpolations. 67 The effects of the electron-phonon coupling (EPC) on electronic band structure are estimated based on the Fr¨olich model, 68,69 from which the renormalization energy of the band edges can be evaluated as, 70 ∆E ≈ −ωLO (α + 0.0159α2 + 0.000806α3),

8

ACS Paragon Plus Environment

(9)

Page 9 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

where α is the coupling constant,

α≡

e2 2m∗ ωLO 1 1 1 ( )2 ( − ). 2~ωLO ~ ε∞ ε0

(10)

We calculated the longitudinal optical phonon frequency ωLO at the Γ point (LO-TO splitting is considered, see Supporting Information for details), the high-frequency (i.e. electroniconly) dielectric constant ε∞ , the static dielectric constant ε0 by the density functional perturbation theory (DFPT) method, 71,72 and the effective mass m∗ at the VBM and CBM by the finite difference approach. For the slab model calculations of surface energies and absolute band positions, the thickness of the slab is as large as 13-25 ˚ A to ensure the center of the slab can be regarded as the bulk phase. We have also made a convergence test of surface energy with respect to slab thickness, see Table S1 in Supporting Information. The thickness of the vacuum is fixed as 15 ˚ A to eliminate the fictitious interaction between mirrors. The macroscopic electrostatic potential is calculated on a 4-times more dense Fourier grid, and then interpolated in spline form on a 10-times more dense grid to ensure the numerical accuracy.

Results and Discussions Crystal structures Ta2 O5 exists in various polymorphs, including the most commonly studied β-Ta2 O5 (space group Pmmm), 41,44,45,73–75 which is believed to be stable at low temperature and undergo a phase transition to α-Ta2 O5 at rather high temperature (∼1630 K). 73 Another typical phase is labeled by ǫ-Ta2 O5 (space group C2/c), 76 which can be prepared by both the hydrothermal approach 76 and high-pressure synthesis. 77 Besides, the X-ray powder diffraction measurements also indicate the existence of the polymorph with large number of oxygen vacancies, 78 hexagonal phase, 79,80 high-pressure prepared phase, 77 as well as the amorphous phase. 81 In

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1: Crystal structures of Ta2 O5 , TaON and Ta3 N5 . this work, we mainly focus on the first two representative phases, i.e. β-Ta2 O5 and ǫ-Ta2 O5 , since they are prepared in relatively mild conditions and are commonly considered to be stable enough for practical applications, especially for the β phase. We note in passing that our consideration of different polymorphic structures of Ta2 O5 is by no means exhaustive, and our purpose here is mainly to investigate the effects of changing crystal structures on electronic band structure. In contrast to the complicated polymorphism of Ta2 O5 , βTaON (space group P21 /c) and β-Ta3 N5 (space group Cmcm), in spite of the presence of other phases, 42,82 are considered to be the most important polymorphs in practice. All crystal structures considered in this work are illustrated in Figure 1. The lattice parameters optimized by PBEsol and those measured by experiment are collected in Table 1. For β-Ta2 O5 and ǫ-Ta2 O5 , the optimized lattice constants obtained from PBSsol calculations agree very well with experiment, which is consistent with the general good performance of the PBEsol GGA functional in predicting crystal structural

10

ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 1: Results of lattice parameters and relative energies of conventional unit cells of Ta2 O5 , TaON and Ta3 N5 from PBEsol. The quantities in the parentheses are obtained from the experimental literature cited in the Ref. column. ∆E denotes the relative energy with respect to β-Ta2 O5 . Species β-Ta2 O5 β ′ -Ta2 O5 ∗ ǫ-Ta2 O5 β-TaON β-Ta3 N5

a [˚ A] b [˚ A] c [˚ A] β[◦ ] ∆E [eV/f.u.] 6.43 (6.22) 3.67 (3.68) 7.72 (7.79) 0.00 6.36 3.75 7.60 -1.82 12.83 (12.79) 4.87 (4.85) 5.52 (5.53) 103.8 (104.3) -2.37 4.95 (4.96) 5.02 (5.03) 5.17 (5.18) 99.8 (99.6) 3.88 (3.89) 10.22 (10.21) 10.24 (10.26)

Ref. 74 77 83 84

The structure is optimized after a small displacement along the maximum imaginary phonon mode. ∗

properties, 65 and therefore the PBEsol optimized structure can be used as a good starting point for the following electronic structure calculations. Table 1 also collects the relative energy per formula unit (f.u.) of different polymorphic phases of Ta2 O5 , taking the energy of β-Ta2 O5 as zero. Surprisingly, β-Ta2 O5 is much less stable than ǫ-Ta2 O5 by about 2.37 eV per formula unit, which is a rather large energy difference, and the consideration of vibration contributions (zero-point energy and finite temperature effects) or high-order correlation effects beyond GGA (e.g. those accounted for by using adiabatic connection fluctuationdissipation theorem in the random phase approximation 85 ) is unlikely to change the relative stability of the two phases. The poor stability of β-Ta2 O5 compared to ǫ-Ta2 O5 is surprising since β-Ta2 O5 is often regarded as the stable phase at the ambient condition, 73,74 and is often used in theoretical study of Ta2 O5 . 44,45,75 It is therefore worth to have a closer look at the stability of βTa2 O5 . The phonon analysis indicates a large imaginary frequency (∼394 cm−1 ) of βTa2 O5 at the Γ point, consistent with the previous report. 41,75 In contrast, ǫ-Ta2 O5 shows no imaginary frequency from the DFPT phonon calculations at the Γ point. Furthermore, the AIMD simulations indicate that even at a low temperature of 50 K, the β phase would undergo a strong distortion after about 230 fs, which follows the vibration mode of the maximum imaginary frequency. By relaxing the symmetry of the β phase and slightly 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 36

distorting the structure in terms of the vibration mode of the maximum imaginary frequency, or taking a low-energy structure from the MD simulation trajectories, we can obtain a crystal structure that is similar to β phase, but without imaginary phonon frequencies, denoted as β ′Ta2 O5 henceforth (see Figure 1 (b)). Structurally β-Ta2 O5 and β ′ -Ta2 O5 differ mainly in two aspects: (i) the TaO6 octahedra in β ′ -Ta2 O5 are strongly distorted, and (ii) the connection between TaO6 octahedra also change significantly, which are clearly illustrated in Figure 2. We can therefore come to the conclusion that β-Ta2 O5 with the Pmmm symmetry is indeed unstable at low temperature and ambient pressure. In the next Section, we will further show that regarding electronic band structures, the results for β-Ta2 O5 are also inconsistent with experimental findings, which implies that β-Ta2 O5 may not be a realistic model to simulate Ta2 O5 commonly prepared experimentally.

Figure 2: (a) The top view from c direction of β-Ta2 O5 , and the maximum imaginary phonon mode is also labeled. (b) The top view from c direction of β ′ -Ta2 O5 . (c)The side view from b direction of β-Ta2 O5 . (d) The side view from b direction of β ′ -Ta2 O5 .

Electronic band structure In this section, we will present the results on the electronic band structure of the three compounds. The band structure diagrams from both PBEsol and G0 W0 @PBEsol are collected in Figure 3. It is obvious that G0 W0 predicts a larger band gap by pushing the VBM to lower energy and CBM to the opposite, with almost unchanged dispersion of band energies. For 12

ACS Paragon Plus Environment

Page 13 of 36

c)

6

6

b)

5

4

4

3

3

energy [eV]

5

2 1

2 1

0

0

-1

-1

-2 S 4

Γ

R

S

Γ

Y

X

Z

Γ

U

Z

d)

-2 Γ 4

3

3

2

2 energy [eV]

energy [eV]

a)

energy [eV]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1

0

-1

-1

Γ

Y

A

B

D

E

C

L

Z

Γ

X

Y

1

0

-2 Z

Y

-2 Γ

Y

T

Z

A

R

S

Γ

Figure 3: Band structures of different compounds: a) β-Ta2 O5 , b) ǫ-Ta2 O5 , c) β-TaON , d) β-Ta3 N5 . The black lines denote the results from PBEsol and the red dots denote the results from G0 W0 @PBEsol. In each case, the valence band maximum (VBM) from PBEsol is taken as the common energy zero (Fermi level) for both PBEsol and G0 W0 band structures.

Table 2: Results of the band gap (in eV) from different theoretical approaches and experiment. Note that the phase of Ta2 O5 used to measure the experimental band gap is undetermined. Species Ta2 O5 β-Ta2 O5 β ′ -Ta2 O5 ǫ-Ta2 O5 β-TaON β-Ta3 N5 a

Ref. 3;

b

PBEsol

HSE06 G0 W0

scGW

G0 W0 with EPC

Expt. 3.9 , 4.2b , 4.5c a

0.31 2.59 3.14 1.86 1.20

Ref. 86;

c

1.52 4.08 4.67 3.02 2.16

1.81 4.25 4.81 2.95 2.20

d

Ref. 39;

Ref. 87;

2.42 5.24 5.76 3.25 2.38

13

3.71 3.96 2.75 2.07

ACS Paragon Plus Environment

2.4 a , 2.8 2.1a

d

The Journal of Physical Chemistry

4.80

1.88

4.78

1.84 ε-Ta2O5

4.76 Eg [eV]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4.74 100 2.96

200

β-Ta2O5

1.80 1.76 100

300

Page 14 of 36

200

300

2.19

2.94

2.16 2.92

β-TaON

2.13

2.90 100 200 300 400 100 Emax [eV]

β-Ta3N5

200

300

Figure 4: The calculated G0 W0 band gaps fitted as a function of highest unoccupied state energy considered. clarity, we summarize the electronic band structure information from different approaches in Table 2. To obtain numerically converged G0 W0 band gaps, we fit the calculated band gap by (see Figure 4) 88 Eg (x) =

a + Eg∞ , x − ǫ0

(11)

where x denotes the energy of highest unoccupied state considered in GW calculations, a, ǫ0 and Eg∞ (i.e. the extrapolated band gap) are fitting parameters. It is interesting to note that for these materials with a 5d0 electronic configuration, the calculated G0 W0 band gap converges quite quickly with respect to the number of unoccupied states, in contrast to the situations observed in the systems with semicore d-states like ZnO. 88 In Table 2, we also present the band gaps of different compounds obtained from the quasi-particle self-consistent GW (scGW) scheme. 89,90 We can see that the band gaps from scGW are systematically larger than those from PBEsol-based G0 W0 , which is consistent in previous findings that the scGW approach without considering vertex corrections tend to significantly overestimate the band gaps of common semiconductors. 90,91 At the level of the GW approximation to the self-energy, the “best” results (i.e. in terms of the agreement with experiment) for normal systems with moderate electron correlation can often be obtained by G0 W0 or partially (energy-only) self-consistent GW0 with standard LDA or GGA Kohn-

14

ACS Paragon Plus Environment

Page 15 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Sham orbital energies and wave functions as the input. 33,88,90,92,93 Thus we mainly consider the G0 W0 @PBEsol results in this work. From Table 2, we can see that PBEsol-GGA severely underestimates the band gaps of all these compounds, and both the hybrid functional (HSE06) and G0 W0 @PBEsol increase the band gaps dramatically. It is noteworthy that the band gap of β-Ta2 O5 is unexpectedly small, much smaller than that of ǫ-Ta2 O5 and the experimental values for Ta2 O5 , and it is even smaller than that of TaON and Ta3 N5 calculated at the same level, which is contrary to the general trend that the band gap usually follows the order of oxide > oxynitride > nitride for compounds of the same metal. β−Ta2O5

Ta(6s) Ta(5d) O (2p)

β’-Ta2O5

ε−Ta2O5

0 -8

-6

-4

-2

0 2 4 energy [eV]

6

8

10

Figure 5: DOS of ǫ-Ta2 O5 , β ′ -Ta2 O5 and β-Ta2 O5 by PBEsol. All three DOS are shifted so that the Ta’s inner core 5s states are aligned. The dashed line indicates the VBM of ǫ-Ta2 O5 . The dramatic difference between the band gaps of β-Ta2 O5 and ǫ-Ta2 O5 is remarkable, and calls for an in-depth understanding of its physical origin. To this end, we use the phenomenological ionic model, 31,94,95 in which the band gap of a typical ionic insulating system can be decomposed as 1 a c Eg = (Aa − I c ) + e(|VMa | + |VMc |) − (∆Epol + ∆Epol ) − (WVB + WCB ), 2

15

ACS Paragon Plus Environment

(12)

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

where “a” denotes the anionic species whose states dominate the VBM, and “c” indicates the cationic species whose states dominate the CBM. In the first term (free ions contribution), A and I are the atomic electron affinity and ionization energy, respectively, which are certainly the same for different polymorphic structures of Ta2 O5 . The difference of the third term (the polarization contribution) between different compounds is expected to be small because of their qualitatively similar chemical bonding. We have also confirmed that the contribution of the second term (the Madelung potential VM ), calculated in terms of Bader effective charges, is nearly identical in different phases of Ta2 O5 . Thus the only remaining factor is the band broadening, which can be clearly illustrated by the density of states (DOS) shown in Figure 5. It is clear that both the valence band and conduction band are much wider in β-Ta2 O5 than in β ′ -Ta2 O5 and ǫ-Ta2 O5 , which is therefore the main origin of the smaller band gap in the β phase. The strong band broadening can be further attributed to the crystal structure features of the β phase. From Figure 1, we can see that the TaO6 octahedra in the β phase are very regular with only little distortion, compared with those in the β ′ or ǫ phase. This feature will significantly enhance the overlap between neighboring atomic orbitals and therefore lead to a large band broadening. In addition, the Ta-O bond length in the β phase is slightly shorter (about 1.9 ˚ A), which can also enhance the effective overlap between atomic orbitals on the ligand O atoms and the central Ta atom. From that point of view, one can also understand naturally that the structurally more strongly distorted β ′ -Ta2 O5 has a significantly larger band gap due to relatively weaker inter-atomic orbital overlap (see also the band structure of β ′-Ta2 O5 in Figure S1 in Supporting Information). Another evidence of strong band broadening is the crystal-field splitting of Ta-5d states. In all of the three phases, Ta and O form the octahedral coordination (although the β ′ phase is distorted to a large extent). But from Figure 5, only ǫ-Ta2 O5 and β ′ -Ta2 O5 show obvious splitting between t2g and eg band, and β-Ta2 O5 , with the almost regular octahedral coordination, shows only a flat, continuous d band. The phenomenon derives from the strong band broadening in the β phase, which makes the gap between t2g and eg “squeezed out”. To summarize, the

16

ACS Paragon Plus Environment

Page 16 of 36

Page 17 of 36

reason why the β phase has such a small band gap is the strong overlap between orbitals derived from the nearly regular octahedral coordination. Together with the instability of the β phase discussed in the previous section, we can plausibly argue that β-Ta2 O5 with its currently accepted structure can not be regarded as a realistic representation of Ta2 O5 that is used in practical applications. Considering its remarkable stability compared to β-Ta2 O5 and β ′ -Ta2 O5 , we think ǫ-Ta2 O5 is more suitable to serve as the representative phase of Ta2 O5 in the ambient condition, or at least as the main component of amorphous Ta2 O5 samples that are usually investigated in experiment. In the following, we will mainly consider theoretical results from ǫ-Ta2 O5 unless stated otherwise. Ta(6s) Ta(5d) O(2p) N(2p)

Ta3N5

PDOS [arb. unit]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TaON

Ta2O5

-8

-6

-4

-2

0

2

energy [eV]

4

6

8

10

Figure 6: DOS of ǫ-Ta2 O5 , β-TaON and β-Ta3 N5 by PBEsol. All three DOS are shifted so that the Ta’s inner core 5s states are aligned. The dashed line indicates the VBM of ǫ-Ta2 O5 . Another obvious feature from Table 2 is the trend that the band gap decreases in the order of oxide > oxynitride > nitride, which is consistent with experimental findings. 3 To obtain an in-depth understanding of underlying physics, we plot the DOS of three compounds, shown in Figure 6. When aligned in term of the semicore Ta-5s states that are present in all three compounds, the VBM of Ta3 N5 and TaON, which are both dominated by N-2p characters, is about 0.5 eV higher than that of ǫ-Ta2 O5 , which is consistent with the energetic order of O- and N-2p atomic orbital. Another obvious feature is that conduction bands of TaON and 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 36

Ta3 N5 move towards lower energy. This can be attributed to the difference of the electrostatic potential at the Ta site, since the CBM is mainly of Ta-5d characters. In particular, the calculated Madelung potential at the position of Ta, decreases from oxide to nitride (29.77 V in Ta2 O5 , 28.80 V in TaON , and 27.11 V in Ta3 N5 ), which indicates that the CBM state is more stable in Ta3 N5 than in Ta2 O5 from a purely electrostatic point of view. Figure 6 also shows that the band widths of the three compounds are similar, which is consistent with nearly same Ta-O(N) bond lengths (all being about 2.0 ˚ A) in their structures. We can therefore conclude that both the orbital characters and the electrostatic factor are essential to understand the band gap difference between the three compounds: higher-energy of N2p orbital pushes up the VBM while the less repulsive electrostatic potential experienced by electrons occupying CBM pushes downs the conduction band energy in oxynitride and nitride. Table 3: Results of high-frequency and static dielectric constant, longitudinal optical phonon frequency at Γ point (in 2π·THz), effective mass at VBM and CBM (in m0 ), the corresponding coupling constant and the renormalization energy (in eV) by PBEsol. The dielectric constant and the effective mass are arithmetic trace-averaged and geometric trace-averaged respectively. Species β ′ -Ta2 O5 ǫ-Ta2 O5 β-TaON β-Ta3 N5

ε∞ 5.49 5.68 8.89 10.97

ε0 32.67 32.45 26.09 42.85

ωLO m∗VBM 29.48 3.26 40.87 5.64 32.13 1.22 30.24 0.52

m∗CBM 0.71 1.71 0.86 0.64

αVBM 2.89 3.09 0.83 0.51

αCBM 1.35 1.70 0.70 0.73

∆EVBM +0.37 +0.55 +0.11 +0.06

∆ECBM -0.17 -0.30 -0.09 -0.07

We further compare our theoretically calculated band gaps to the experiment values. Since the results from HSE06 and G0 W0 are close to each other, and the latter is in general more accurate than the former, we will mainly consider the G0 W0 results in the following discussion. It should be noted that experimental band gaps often suffer from several uncertainty factors including the sample quality, the particular experimental techniques, and the way to extract the band gap from experiment, to name a few only. 88 As a result, the experimental values for the band gap of a particular compound usually exhibit significant scattering. It is therefore necessary to exercise some caution when comparing theory and experiment as far 18

ACS Paragon Plus Environment

Page 19 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

as the band gap is concerned. For example, the band gap of Ta2 O5 , measured by several groups, 3,86,87 often without characterizing the crystal structure simultaneously, scatter in the range of 3.9-4.5 eV. With these uncertainties taken into account, we can see that for Ta3 N5 the G0 W0 band gap agrees very well with experiment, but for TaON and Ta2 O5 (using the theoretical result for ǫ-Ta2 O5 ), the G0 W0 band gaps show a certain overestimation, especially for Ta2 O5 . Such over-estimation of the band gap by LDA or GGA-based G0 W0 has also been observed in other early transition metal oxides with apparent d0 electronic configuration such as TiO2 , 28,29,96 SrTiO3 , 96,97 ATaO3 (A=Li, Na, K), 31,97 LiNbO3 , 98 and WO3 , 69 to name a few typical cases only. For such compounds, there is increasing evidence indicating that the electron-phonon coupling can play a significant role in reducing the band gap. 28,69,98 Following Refs. 69, we estimate the effects of EPC on electronic band structure by calculating the valence and conduction band renormalization energy ∆E based on the simplified Fr¨olich model (see eq. 9), and the related quantities for β ′ -Ta2 O5 , ǫ-Ta2 O5 , βTa3 N5 , and β-TaON are collected in Table 3. We can see that compared with Ta-oxynitride and nitride, ǫ-Ta2 O5 shows a significantly smaller electronic dielectric constant (i.e. ε∞ ) and larger effective mass (m∗ ), which can be attributed to its larger band gap and flatter valence and conduction band edges, as shown in Figure 3. As a result, the consideration of EPC reduces the band gap of ǫ-Ta2 O5 by as much as 0.85 eV. A slightly weaker decreasing is also observed for β ′ -Ta2 O5 . For Ta3 N5 and TaON , the effects of EPC are much weaker, and the band gap decreases only by 0.1-0.2 eV. We note in passing that our result for Ta3 N5 is in agreement with the previous calculations. 99 Having considered the EPC effects, the G0 W0 band gaps of all three compounds show a rather good agreement with experiment. Since experimental band gaps of the compounds considered in this work are all obtained from optical absorption, excitonic effects that describe electron-hole interaction may also need to be considered. However, in 3D bulk materials, the excitonic effects on the band gap are expected to be weak. 100 In particular, Morbec et al 99 reported, based on the GW plus Bethe-Salpeter equation (BSE) approach, that the exciton binding energy in Ta3 N5 is only

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 36

about 0.05 eV. We can therefore expect that considering excitonic effects plays a marginal role for the band gaps of Ta-oxides and oxynitrides as well.

Absolute band edge positions We further consider the energy of VBM and CBM with respect to the vacuum level, the socalled absolute band edge positions, using the approach described in the preceding section. Since the absolute band energies is highly surface dependent, it is important to search for the most stable exposed surface of concerned systems. The relative stability of different exposed surfaces of β-TaON 101 and β-Ta3 N5 102 has been theoretically investigated previously. In contrast, the surface stability of Ta2 O5 has been rarely studied previously to our best knowledge, probably due to the complex polymorphism of Ta2 O5 . It was found that in both β-TaONand β-Ta3 N5 , the (100) surface is the most stable, while in TaON, the (111) surface is only slightly less stable than the (100) one. 101 Based on the previous study, for TaON we consider only its (100) and (111) surfaces with either O- or N-termination, for βTa3 N5 we consider its (100) surface, and for the surfaces of Ta2 O5 , we have considered several low-index surfaces of ǫ-Ta2 O5 to the find the most stable exposed surface. The calculated surface energies are collected in Table 4. The detailed surface structures and corresponding macroscopic electrostatic potentials can be found in Figure S2 in Supporting Information. Among the low-index surfaces of ǫ-Ta2 O5 we have considered, the (100) surface has the lowest surface energy, and therefore will be used for the calculation of absolute band positions. For TaON, we find that on both (100) and (111) surfaces, the O-termination is more stable than the N-termination. The calculated surface energies are quite close to those reported in Ref. 101, except that the difference between (100) and (111) surfaces are slightly larger. Both (100) and (111) surfaces can be important in reality as a result of the closeness of their surface energies. For β-Ta3 N5 , our calculated surface energy, 1.41 J/m2 , for the (100) surface is larger than that in Ref. 102 by about 0.3 J/m2 , possibly due to the use of different computational approaches. 20

ACS Paragon Plus Environment

Page 21 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 4: Results of surface stability. The value in parentheses are from the previous work indicated in Ref. column. System ǫ-Ta2 O5

β-TaON

β-Ta3 N5 a

from Ref. 101;

b

surface slab thickness [˚ A] (100) 23.8 (010)-O 14.1 (010)-Ta 14.5 (001) 13.8 (100)-O 18.9 (100)-N 18.9 (111)-O 16.0 (111)-N 16.2 (100) 13.6

Esurf [J/m2 ] 0.92 1.13 2.28 1.71 1.44 (1.65a ) 3.47 1.53 (1.66a ) 2.32 1.41 (1.1b)

from Ref. 102.

Table 5: Results of absolute band positions from PBEsol and G0 W0 for the surfaces of Ta2 O5 , TaON , and Ta3 N5 compared to available experimental values. Surface PBEsol ǫ-Ta2 O5 (100) (010)-O β-TaON (100)-O (111)-O β-Ta3 N5 (100) a b

-7.22 -7.66 -7.04 -6.52 -5.42

EVBM [eV] G0 W0 Expt. -7.83a, -7.93b -7.94 -8.37 -6.6b -7.67 -7.15 -6.02a,b -6.04

PBEsol -4.90 -5.34 -5.38 -4.86 -4.38

ECBM [eV] G0 W0 Expt. -3.93a , -4.03b -3.98 -4.41 -4.1b -4.92 -4.40 -4.14a , -3.92b -3.99

by electrochemical measurement from Ref. 3; by ultra-violet photoemission spectroscopy measurement from Ref. 3

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 36

We calculated absolute band positions with or without the GW correction and the results are collected in Table 5.

The related experimental results obtained from photo-

electrochemical (PEC) measurement or ultra-violet photoemission spectroscopy (UPS) are also collected for comparison. We can see the GW correction is pretty significant to push EVBM to lower energy and ECBM to higher energy. The G0 W0 results for ǫ-Ta2 O5 obtained from the (100) surface are nearly identical to experimentally measured values, and those from the (010)-O terminated surface show a systematic shift towards lower energy by about 0.4 eV. For β-TaON, the calculated G0 W0 VBM and CBM energies from the (100)-O surfaces are systematically lower by nearly 1 eV than experimental values, but the discrepancy is significantly reduced when considering the results obtained from the (111)-O surface, which could imply that the (111)-O surface of β-TaON may play a more important role than (100) surface in determining the absolute band positions of this compound. For β-Ta3 N5 , the G0 W0 results again agree very well with experimental values. Considering the significant error bar in experimental values, the overall agreement between the G0 W0 prediction and experiment is remarkable. Our investigations based on the G0 W0 approach have therefore confirmed, from a theoretical perspective, that in all of the three compounds, EVBM is more negative than -5.6 eV, and ECBM is more positive than -4.4 eV, which indicates their potential application for unassisted water splitting, and visible light photocatalytic activity for TaON and Ta3 N5 , in particular. It is important to note that our discussion on the application of these materials in photo-splitting of water is based on the results obtained from vacuum calculations, similar to a lot of previous studies. 28,32,33,37,69 However, an explicit consideration of the interface with water can be important for a more realistic and quantitative treatment of the issue, and has been addressed by several recent studies. 103–105

22

ACS Paragon Plus Environment

Page 23 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Conclusions To summarize, we have investigated electronic band structures and absolute band positions of Ta2 O5 , TaON and Ta3 N5 by different first-principles approaches, including GGA (PBEsol), the hybrid functional (HSE06) and the GW approach (G0 W0 @PBEsol). Both the hybrid functional and the GW approach can well reproduce experimental band gaps, especially after the effects of the electron-phonon coupling on electronic band structures have been taken into account based on the simplified Fr¨olich model. The latter are particularly significant for ǫ-Ta2 O5 as a result of its small high-frequency dielectric constant and large effective mass. We also analyzed the origin of the band gap difference among the oxide, oxynitride and nitride, and we found that both the energy difference between oxygen and nitrogen 2p-orbitals, and the difference in the electrostatic (Madelung) potential in these compounds play important roles in determining their electronic band structure. We have investigated the relative stability and electronic band structures of several polymorphic phases of Ta2 O5 , and we found that the commonly studied β-Ta2 O5 is structurally unstable, and leads to a significantly distorted structure once the crystal symmetry is relaxed; ǫ-Ta2 O5 is energetically much more stable, and its band gap obtained from the G0 W0 approach with the EPC correction agrees very well with commonly measured experimental values, implying that ǫ-Ta2 O5 is more suitable than β-Ta2 O5 as a realistic representation of Ta2 O5 that is measured experimentally. We further show that electronic band structure of Ta2 O5 in different polymorphic structures is strongly affected by different inter-atomic orbital overlapping, which can be further related to inter- and intra-TaO6 octahedra structural distortion. We have also evaluated the absolute valence and conduction band positions with respect to the vacuum level in Ta2 O5 , TaON and Ta3 N5 by combining the PBEsol slab model calculations with the G0 W0 correction to the bulk valence band maximum, which can well reproduce experimental findings, and has confirmed that TaON and Ta3 N5 are visible light active and have appropriate band edge positions for photo-catalytic splitting of water to hydrogen and oxygen. 23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Acknowledgments This work is supported by National Natural Science Foundation of China (Project Numbers 21373017, 21673005 and 21321001). ZHC acknowledges the support of the Hui-Chun Chin and Tsung-Dao Lee Chinese Undergraduate Research Endowment (CURE).

Supporting Information Information on the details for phonon calculations, the band structure of β ′ -Ta2 O5 , surface structures and the corresponding macroscopic electrostatic potentials for different surfaces considered in this work, and the convergence of the calculated surface energies with respect to the slab thickness. This material is available free of charge via the Internet at http://pubs.acs.org.

References (1) Fujishima, A.; Honda, K. Electrochemical Photolysis of Water at a Semiconductor Electrode. Nature 1972, 238, 37. (2) Hu, S.; Shaner, M. R.; Beardslee, J. A.; Lichterman, M.; Brunschwig, B. S.; Lewis, N. S. Amorphous TiO2 Coatings Stabilize Si, GaAs, and GaP Photoanodes for Efficient Water Oxidation. Science 2014, 344, 1005–1009. (3) Chun, W.-J.; Ishikawa, A.; Fujisawa, H.; Takata, T.; Kondo, J. N.; Hara, M.; Kawai, M.; Matsumoto, Y.; Domen, K. Conduction and Valence Band Positions of Ta2O5, TaON, and Ta3N5 by UPS and Electrochemical Methods. J. Phys. Chem. B 2003, 107, 1798–1803. (4) Liao, P.; Carter, E. A. New Concepts and Modeling Strategies to Design and Evaluate

24

ACS Paragon Plus Environment

Page 24 of 36

Page 25 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Photo-electro-catalysts Based on Transition Metal Oxides. Chem. Soc. Rev. 2013, 42, 2401 –2422. (5) Ping, Y.; Rocca, D.; Galli, G. Electronic Excitations in Light Absorbers for Photoelectrochemical Energy Conversion: First Principles Calculations Based on Many Body Perturbation Theory. Chem. Soc. Rev. 2013, 42, 2437 –2469. (6) Chen, X.; Shen, S.; Guo, L.; Mao, S. S. Semiconductor-based Photocatalytic Hydrogen Generation. Chem. Rev. 2010, 110, 6503. (7) Kudo, A.; Miseki, Y. Heterogeneous Photocatalyst Materials for Water Splitting. Chem. Soc. Rev. 2009, 38, 253–278. (8) Ahmed, M.; Guo, X. A Review of Metal Oxynitrides for Photocatalysis. Inorg. Chem. Front. 2016, 3, 578 – 590. (9) Ebbinghaus, S. G.; Abicht, H.-P.; Dronskowski, R.; Mller, T.; Reller, A.; Weidenkaff, A. Perovskite-related Oxynitrides - Recent Developments in Synthesis, Characterisation and Investigations of Physical Properties. Prog. Solid State Chem. 2009, 37, 173–205. (10) Fuertes, A. Metal Oxynitrides as Emerging Materials with Photocatalytic and Electronic Properties. Mater. Horiz. 2015, 2, 453 – 461. (11) Fuertes, A. Chemistry and Applications of Oxynitride Perovskites. J. Mater. Chem. A 2012, 22, 3293. (12) Moriya, Y.; Takata, T.; Domen, K. Recent Progress in the Development of (Oxy)nitride Photocatalysts for Water Splitting Under Visible-light Irradiation. Coord. Chem. Rev. 2013, 257, 1957 – 1969.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 36

(13) Maeda, K. (Oxy)nitrides with d0-electronic Configuration as Photocatalysts and Photoanodes that Operate Under a Wide Range of Visible Light for Overall Water Splitting. Phys. Chem. Chem. Phys. 2013, 15, 10537 – 10548. (14) Martin, R. M. Electronic Structure: Basic Theory and Practical Methods; Cambridge University Press: Cambridge, UK, 2004. (15) Carter, E. A. Challenges in Modeling Materials Properties Without Experimental Input. Science 2008, 321, 800 – 803. (16) Bhatt, M. D.; Lee, J. S. Recent Theoretical Progress in the Development of Photoanode Materials for Solar Water Splitting Photoelectrochemical Cells. J. Mater. Chem. A 2015, 3, 10632–10659. (17) Aryasetiawan, F.; Gunnarsson, O. The GW Method. Rep. Prog. Phys. 1998, 61, 237. (18) Perdew, J. P.; Parr, R. G.; Levy, M.; Balduz, J. L. Density-Functional Theory for Fractional Particle Number: Derivative Discontinuities of the Energy. Phys. Rev. Lett. 1982, 49, 1691. (19) Sham, L. J.; Schl¨ uter, M. Density-Functional Theory of the Energy Gap. Phys. Rev. Lett. 1983, 51, 1888. (20) Perdew, J. P.; Levy, M. Physical Content of the Exact Kohn-Sham Orbital Energies: Band Gaps and Derivative Discontinuities. Phys. Rev. Lett. 1983, 51, 1884. (21) Jiang, H.; Shen, Y.-C. Ionization Potentials of Semiconductors from First-principles. J. Chem. Phys. 2013, 139, 164114. (22) Hedin, L. New Method for Calculating the One-Particle Green’s Function with Application to the Electron-Gas Problem. Phys. Rev. 1965, 139, A796. (23) Hedin, L.; Lundqvist, B. I. Effects of Electron-electron and Electron-phonon Interactions on the One-electron States of Solids. Solid State Phys. 1970, 23, 1. 26

ACS Paragon Plus Environment

Page 27 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(24) W. G. Aulbur, L. J.; Wilkins, J. W. Quasiparticle Calculations in Solids. Solid State Phys. 1999, 54, 1. (25) Rinke, P.; Qteish, A.; Neugebauer, J.; Scheffler, M. Exciting Prospects for Solids: Exact-exchange Based Functionals Meet Quasiparticle Energy Calculations. Phys. Stat. Sol. (B) 2008, 245, 929. (26) Jiang, H. The GW method: Basic Principles, Latest Developments and Its Applications for d- and f-electron Systems. Acta. Phys. -Chim. Sin. 2010, 26, 1017. (27) Vidal, J.; Botti, S.; Olsson, P.; Guillemoles, J.-F.; Reining, L. Strong Interplay between Structure and Electronic Properties in CuIn(S,Se)2: A First-Principles Study. Phys. Rev. Lett. 2010, 104, 056401. (28) Kang, W.; Hybertsen, M. S. Quasiparticle and Optical Properties of Rutile and Anatase TiO2. Phys. Rev. B 2010, 82, 085203. (29) Chiodo, L.; Garcia-Lastra, J. M.; Iacomino, A.; Ossicini, S.; Zhao, J.; Petek, H.; Rubio, A. Self-energy and Excitonic Effects in the Electronic and Optical Properties of TiO2 Crystalline Phases. Phys. Rev B 2010, 82, 045207. (30) Christensen, N. E.; Svane, A.; Laskowski, R.; Palanivel, B.; Modak, P.; Chantis, A. N.; van Schilfgaarde, M.; Kotani, T. Electronic Properties of 3R-CuAlO2 Under Pressure: Three Theoretical Approaches. Phys. Rev. B 2010, 81, 045203. (31) Wang, H.; Wu, F.; Jiang, H. Electronic Band Structures of ATaO3(A = Li, Na, and K) from First-Principles Many-Body Perturbation Theory. J. Phys. Chem. C 2011, 115, 16180. (32) Jiang, H. Structural and Electronic Properties of ZrX2 and HfX2 (X = S and Se) from First Principles Calculations. J. Chem.Phys. 2011, 134, 204705.

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(33) Jiang, H. Electronic Band Structures of Molybdenum and Tungsten Dichalcogenides by the GW Approach. J. Phys. Chem. C 2012, 116, 7664 – 7671. (34) Toroker, M. C.; Kanan, D. K.; Alidoust, N.; Isseroff, L. Y.; Liao, P.; Carter, E. A. First Principles Scheme to Evaluate Band Edge Positions in Potential Transition Metal Oxide Photocatalysts and Photoelectrodes. Phys. Chem. Chem. Phys. 2011, 13, 16644. (35) Hinuma, Y.; Gr¨ uneis, A.; Kresse, G.; Oba, F. Band Alignment of Semiconductors from Density-functional Theory and Many-body Perturbation Theory. Phys. Rev. B 2014, 90, 155405. (36) Gr¨ uneis, A.; Kresse, G.; Hinuma, Y.; Oba, F. Ionization Potentials of Solids: The Importance of Vertex Corrections. Phys. Rev. Lett. 2014, 112, 096401. (37) Stevanovic, V.; Lany, S.; Ginley, D. S.; Tumas, W.; Zunger, A. Assessing Capability of Semiconductors to Split Water Using Ionization Potentials and Electron Affinities Only. Phys. Chem. Chem. Phys. 2014, 16, 3706 – 3714. (38) Abe, R.; Takata, T.; Sugihara, H.; Domen, K. Photocatalytic Overall Water Splitting Under Visible Light by TaON and WO3 with an IO3-/I- Shuttle Redox Mediator. Chem. Comm. 2005, 3829 – 3831. (39) Abe, R.; Higashi, M.; Domen, K. Facile Fabrication of an Efficient Oxynitride TaON Photoanode for Overall Water Splitting into H2 and O2 under Visible Light Irradiation. J. Am. Chem. Soc. 2010, 132, 11828–11829. (40) Ho, C.-T.; Low, K.-B.; Jash, P.; Shen, H.; Snee, P. T.; Meyer, R. J. Formation of Sol Gel-Derived TaOxNy Photocatalysts. Chem. Mater. 2011, 23, 4721 – 4725. (41) Andreoni, W.; Pignedoli, C. A. Ta2O5 Polymorphs: Structural Motifs and Dielectric Constant from First Principles. Appl. Phys. Lett. 2010, 96, 062901.

28

ACS Paragon Plus Environment

Page 28 of 36

Page 29 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(42) Bredow, T.; Lumey, M.-W.; Dronskowski, R.; Schilling, H.; Pickardt, J.; Lerch, M. Structure and Stability of TaON Polymorphs. Z. Anorg. Allg. Chem. 2006, 632, 1157– 1162. (43) Chen, S.; Wang, L.-W. Intrinsic Defects and Electronic Conductivity of TaON: Firstprinciples Insights. Appl. Phys. Lett. 2011, 99, 222103. (44) Sahu, B. R.; Kleinman, L. Theoretical Study of Structural and Electronic Properties of β-Ta2O5 and δ-Ta2O5. Phys. Rev. B 2004, 69, 165202. (45) Nashed, R.; I.Hassan, W. M.; Ismail, Y.; K.Allam, N. Unravelling the Interplay of Crystal Structure and Electronic Band Structure of Tantalum Oxide (Ta2O5). Phys. Chem. Chem. Phys. 2013, 15, 1352–1357. (46) Yashima, M.; Lee, Y.; Domen, K. Crystal Structure and Electron Density of Tantalum Oxynitride, a Visible Light Responsive Photocatalyst. Chem. Mater. 2007, 19, 588– 593. (47) M. Fang, C.; Orhan, E.; Wijs, G. A. d.; T.Hintzen, H.; Groot, R. A. d.; Marchand, R.; Saillard, J.-Y.; With, G. d. The Electronic Structure of Tantalum (Oxy)nitrides TaON and Ta3N5. J. Mater. Chem. 2001, 11, 1248–1252. (48) Li, P.; Fan, W.; Li, Y.; Sun, H.; Cheng, X.; Zhao, X.; Jiang, M. First-Principles Study of the Electronic, Optical Properties and Lattice Dynamics of Tantalum Oxynitride. Inorg. Chem. 2010, 49, 6917–6924. (49) Harb, M.; Sautet, P.; Nurlaela, E.; Raybaud, P.; Cavallo, L.; Domen, K.; Basset, J.-M.; Takanabe, K. Tuning the Properties of Visible-light-responsive Tantalum (Oxy)nitride Photocatalysts by Non-stoichiometric Compositions: A First-principles Viewpoint. Phys. Chem. Chem. Phys. 2014, 16, 20548–20560.

29

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(50) Maeda, K.; Domen, K. New Non-Oxide Photocatalysts Designed for Overall Water Splitting under Visible Light. J. Phys. Chem. C 2007, 111, 7851–7861. (51) Wang, S.; Wang, L.-W. Atomic and Electronic Structures of GaN / ZnO Alloys. Phys. Rev. Lett. 2010, 104, 065501. (52) Godby, R. W.; Schl¨ uter, M.; Sham, L. J. Self-energy Operators and Exchangecorrelation Potentials in Semiconductors. Phys. Rev. B 1988, 37, 10159. (53) Burke, K.; Ernzerhof, M.; Perdew, J. P. The Adiabatic Connection Method: A Nonempirical Hybrid. Chem. Phys. Lett. 1997, 265, 115–120. (54) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207. (55) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Erratum: “Hybrid Functionals Based on a Screened Coulomb Potential” [J. Chem. Phys.118, 8207 (2003)]. J. Chem. Phys. 2006, 124, 219906. (56) Heyd, J.; Peralta, J. E.; Scuseria, G. E.; Martin, R. L. Energy Band Gaps and Lattice Parameters Evaluated with the Heyd-Scuseria-Ernzerhof Screened Hybrid Functional. J. Chem. Phys. 2005, 123, 174101. (57) Marsman, M.; Paier, J.; Stroppa, A.; Kresse, G. Hybrid Functionals Applied to Extended Systems. J. Phys. : Condens. Matter 2008, 20, 064201. (58) Garza, A. J.; Scuseria, G. E. Predicting Band Gaps with Hybrid Density Functionals. J. Phys. Chem. Lett. 2016, 7, 4165 – 4170. (59) Fetter, A. L.; Walecka, J. D. Quantum Theory of Many-particle Systems; McGrawHill, New York, 1971. (60) Bechstedt, F. Prinples of Surface Physics; Springer: Berlin, 2003. 30

ACS Paragon Plus Environment

Page 30 of 36

Page 31 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(61) Fall, C. J.; Binggeli, N.; Baldereschi, A. Deriving Accurate Work Functions from Thin-slab Calculations. J. Phys.: Condens. Matter 1999, 11, 2689. (62) Kresse, G.; Furthm¨ uller, J. Efficient Iterative Schemes for Ab Initio Total-energy Calculations Using a Plane-wave Basis Set. Phys. Rev. B 1996, 54, 11169. (63) Shishkin, M.; Kresse, G. Implementation and Performance of the Frequency-dependent GW Method with the PAW Framework. Phys. Rev. B 2006, 74, 035101. (64) Bl¨ochl, P. E. Projector Augmented-wave Method. Phys. Rev. B 1994, 50, 17953. (65) Perdew, J. P.; Ruzsinszky, A.; Csonka, G. I.; Vydrov, O. A.; Scuseria, G. E.; Constantin, L. A.; Zhou, X.; Burke, K. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Phys. Rev. Lett. 2008, 100, 136406. (66) Kresse, G.; Marsman, M.; Furthm¨ uller, J. VASP the Guide. 2016; Sec. 7.6. (67) Mostofi, A. A.; Yates, J. R.; Lee, Y.-S.; Souza, I.; Vanderbilt, D.; Marzari, N. Wannier90: A Tool for Obtaining Maximally-localised Wannier Functions. Comput. Phys. Commun. 2008, 178, 685–699. (68) Mitra, T. K.; Chatterjee, A.; Mukhopadhyay, S. Polarons. Phys. Rep. 1987, 153, 91–207. (69) Ping, Y.; Rocca, D.; Galli, G. Optical Properties of Tungsten Trioxide from Firstprinciples Calculations. Phys. Rev. B 2013, 87, 165203. (70) Smondyrev, M. A. Diagrams in the Polaron Model. Theor. Math. Phys. 1986, 68, 653–664. (71) Gajdoˇs, M.; Hummer, K.; Kresse, G.; Furthm¨ uller, J.; Bechstedt, F. Linear Optical Properties in the Projector-augmented Wave Methodology. Phys. Rev. B 2006, 73, 045112. 31

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 36

(72) Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and Related Crystal Properties from Density-functional Perturbation Theory. Rev. Mod. Phys. 2001, 73, 515. (73) Lagergren, S.; Magn´eli, A. On the Tantalum - Oxygen System. Acta. Chem. Scand. 1952, 6, 444–446. (74) Aleshina, L. A.; Loginova, S. V. Rietveld Analysis of X-ray Diffraction Pattern from β-Ta2O5 Oxide. Crystallogr. Rep. 2002, 47, 415–419. (75) Wu, Y.-N.; Li, L.; Cheng, H.-P. First-principles Studies of Ta2O5 Polymorphs. Phys. Rev. B 2011, 83, 144105. (76) Izumi, F.; Kodama, H. A New Modification of Tantalum(V) Oxide. Journal of the Less Common Metals 1979, 63, 305–307. (77) Zibrov, I. P.; Filonenko, V. P.; Sundberg, M.; Werner, P.-E. Structures and Phase Transitions of B -Ta2O5 and Z -Ta2O5: Two High-pressure Forms of Ta2O5. Acta Crystallogr., Sect. B: Struct. Sci 2000, 56, 659–665. (78) Stephenson, N. C.; Roth, R. S. Structural Systematics in the Binary System Ta2O5WO3. V. The Structure of the Low-temperature Form of Tantalum Oxide L-Ta2O5. Acta Crystallogr. Sect. B: Struct. Crystallogr. Cryst. Chem. 1971, 27(5), 1037–1044. (79) Terao, N. Structure des Oxides de Tantale. Jpn. J. Appl. Phys. 1967, 6, 21. (80) Fukumoto, A.; Miwa, K. Prediction of Hexagonal Ta2O5 Structure by First-principles Calculations. Phys. Rev. B 1997, 55, 11155. (81) Wu, X. M.; Soss, S. R.; Rymaszewski, E. J.; Lu, T. M. Dielectric Constant Dependence of Poole-Frenkel Potential in Tantalum Oxide Thin Films. Mater. Chem. Phys. 1994, 38, 297–300.

32

ACS Paragon Plus Environment

Page 33 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(82) Schilling, H.; Stork, A.; Irran, E.; Wolff, H.; Bredow, T.; Dronskowski, R.; Lerch, M. γ-TaON: A Metastable Polymorph of Tantalum Oxynitride. Angew. Chem., Int. Ed. 2007, 46, 2931–2934. (83) Armytage, D.; Fender, B. E. F. Anion Ordering in TaON: A Powder Neutrondiffraction Investigation. Acta Crystallogr., Sect. B: Struct. Crystallogr. and Cryst. Chem. 1974, 30, 809–812. (84) Brese, N. E.; O’Keeffe, M.; Rauch, P.; DiSalvo, F. J. Structure of Ta3N5 at 16 K by Time-of-flight Neutron Diffraction. Acta Crystallogr. Sect. C Cryst. Struct. Commun. 1991, 47, 2291–2294. (85) Cui, Z.-H.; Wu, F.; Jiang, H. First-principles Study of Relative Stability of Rutile and Anatase TiO2 Using the Random Phase Approximation. Phys. Chem. Chem. Phys. 2016, 18, 29914 – 29922. (86) Knausenberger, W.; Tauber, R. N. Selected Properties of Pyrolytic Ta2O5 Films. J. Electrochem. Soc. 1973, 120, 927. (87) Chaneliere, C.; Autran, J. L.; Devine, R. A. B.; Balland, B. Tantalum Pentoxide (Ta2O5) Thin Films for Advanced Dielectric Applications. Mater. Sci. Eng., R 1998, 22, 269–322. (88) Jiang, H.; Blaha, P. GW with Linearized Augmented Planewaves Extended by Highenergy Local Orbitals. Phys. Rev. B 2016, 93, 115203. (89) Faleev, S. V.; van Schilfgaarde, M.; Kotani, T. All Electron Self-consistent GW Approximation: Application to Si, MnO, and NiO. Phys. Rev. Lett. 2004, 93, 126406. (90) Shishkin, M.; Marsman, M.; Kresse, G. Accurate Quasiparticle Spectra from SelfConsistent GW Calculations with Vertex Corrections. Phys. Rev. Lett. 2007, 99, 246403. 33

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 36

(91) An, W.; Wu, F.; Jiang, H.; Tian, G.-S.; Li, X.-Z. Systematic Investigation on Topological Properties of Layered GaS and GaSe under Strain. J. Chem. Phys. 2014, 141, 084701. (92) Shishkin, M.; Kresse, G. Self-consistent GW Calculations for Semiconductors and Insulators. Phys. Rev. B 2007, 75, 235102. (93) Jiang, H.; Gomez-Abal, R.; Rinke, P.; Scheffler, M. Electronic Band Structure of Zirconia and Hafnia Polymorphs from the GW Perspective. Phys. Rev. B 2010, 81, 085119. (94) Cox, P. A. Transition Metal Oxides: An Introduction to their Electronic Structure and Properties; Clarendon Press, 1992. (95) Woodward, P. M.; Mizoguchi, H.; Kim, Y. I.; Stoltzfus, M. W. In Metal Oxides: Chemistry and Applications; Fierro, J. L. G., Ed.; CRC Press, 2006; Chapter 6, p 133. (96) Kang, G.; Kang, Y.; Han, S. Influence of Wave-function Updates in GW Calculations on Titanates. Phys. Rev. B 2015, 91, 155141. (97) Benrekia, A. R.; Benkhettou, N.; Nassour, A.; Driz, M.; Sahnoun, M.; Lebegue, S. Structural, Electronic and Optical Properties of Cubic SrTiO3 and KTaO3: Ab Initio and GW Calculations. Physica B 2012, 407, 2632 – 2636. (98) Schmidt, W. G.; Albrecht, M.; Wippermann, S.; Blankenburg, S.; Rauls, E.; Fuchs, F.; Rdl, C.; Furthmller, J.; Hermann, A. LiNbO3 Ground- and Excited-state Properties from First-principles Calculations. Phys. Rev. B 2008, 77, 035106. (99) Morbec, J. M.; Narkeviciute, I.; Jaramillo, T. F.; Galli, G. Optoelectronic Properties of Ta3N5: A Joint Theoretical and Experimental Study. Phys. Rev. B 2014, 90, 155204.

34

ACS Paragon Plus Environment

Page 35 of 36

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(100) Onida, G.; Reining, L.; Rubio, A. Electronic Excitations: Density Functional Versus Many-body Green’s Function Approaches. Rev. Mod. Phys. 2002, 74, 601. (101) Esch, T. R.; Vidic, J.; Bredow, T. TaON Surface Energetics and Optical Spectra. Phys. Status Solidi RRL 2014, 8, 587–591. (102) Watanabe, E.; Ushiyama, H.; Yamashita, K. Theoretical Studies on the Stabilities and Reactivities of Ta3N5 (100) Surfaces. Chem. Phys. Lett. 2013, 561562, 57–62. (103) Wu, Y.; Chan, M. K. Y.; Ceder, G. Prediction of Semiconductor Band Edge Positions in Aqueous Environments from First Principles. Phys. Rev. B 2011, 83, 235301. (104) Ping, Y.; Sundararaman, R.; Goddard III, W. A. Solvation Effects on the Band Edge Positions of Photocatalysts from First Principles. Phys. Chem. Chem. Phys. 2015, 17, 30499 – 30509. (105) Pham, T. A.; Ping, Y.; Galli, G. Modelling Heterogeneous Interfaces for Solar Water Splitting. Nat. Mater. 2017, Advance Online Publication,DOI: 10.1038/NMAT4803.

35

ACS Paragon Plus Environment

The Journal of Physical Chemistry

TOC Graphic 0 -2

Eabs [eV]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 36

-4 -6

H / H2

VBM

H2O/O2

GW Expt.

-8 -10 -12

Ta2O5

TaON

36

+

CBM

Ta3N5

ACS Paragon Plus Environment