Improving Photocatalytic Properties of SrTiO3 Through SB, NB Doping DFT
Improving Photocatalytic Properties of SrTiO3 Through SB, NB Doping DFT
A systematic study using hybrid density functional theory has been carried out to investigate the synergistic
effect of Sb and N doping on the photocatalytic properties of SrTiO3 under visible light. The calculated band
gap (3.19 eV) for SrTiO3 with the Heyd, Scuseria, and Ernzerhof hybrid functional is found to be very close to
the experimentally observed value of 3.2 eV. Although doping with N is able to enhance the visible light
activity by reducing the effective band gap to 2.31 eV, the localized occupied and unoccupied states in
the forbidden region may affect the photocatalytic activity. However, the presence of Sb not only
passivates those unoccupied states completely, but also shifts the localized occupied states near the
valence band to form a continuum band structure. The introduction of N into the SrTiO3 crystal
structure is favored by the presence of Sb. In the codoped system charge compensation is established,
thereby unwanted vacancy formation will be minimized. The absorption curve for the (Sb, N)-codoped
SrTiO3 is found to shift towards the visible region due to reduction in band gap to 2.66 eV. Moreover, the
Received 18th July 2014
Accepted 10th September 2014
band alignment shows that the (Sb, N)-codoping makes SrTiO3 thermodynamically more suitable for
hydrogen production as compared to the undoped system. Based on the present study, we can propose
DOI: 10.1039/c4ra07289h
(Sb, N)-codoping as one of the effective approaches to improve the photocatalytic activity of SrTiO3 for
www.rsc.org/advances water splitting under visible light irradiation.
This journal is © The Royal Society of Chemistry 2014 RSC Adv., 2014, 4, 45703–45709 | 45703
View Article Online
synthesized (La, N)-codoped SrTiO3, which shows improved pho- 12.5% and N concentration of 4.17%. The exchange and corre-
tocatalytic activity under visible light. However, this photoactivity lation energy density functionals were dened by generalized
has been found to be limited to the decomposition of organic gradient approximation (GGA) with Perdew–Burke–Ernzerhof
compounds. Although, it has been shown in the theoretical study (PBE) scheme41 during geometry optimization. K-point mesh of 8
of Wei et al. that codoping of either nonmetal (H, F, Cl, Br, I) or 8 8 was found to be sufficient to attain convergence. The
metal (V, Nb, Ta, Sc, Y, La) is able to passivates the N-induced total energy convergence tolerance was set to 106 eV per atom
Published on 10 September 2014. Downloaded by Christian Albrechts Universitat zu Kiel on 23/10/2014 22:40:24.
discrete states, their applicability for the photo-splitting of water for the self consistent iteration. Electronic structure calculations
is still to be explored.26 Recently, Yu et al. synthesized (Cr, N)- have been carried out by employing the HSE hybrid functional,
codoped SrTiO3, which has been shown to generate only H2 where the short-range interactions are calculated with exact
during water splitting under visible light.27 In the present study, we exchange mixing in both HF and DFT. The exchange-correlation
employ Sb as codopant to improve the photoactivity of N- energy is expressed as
doped SrTiO3 under visible light. The advantages for choosing Sb
as codopant are: (a) being very similarity in ionic radius (Sb5+: Ref ¼ EHSE
XC ¼ aEX (m) + (1 a)EX
SR PBE,SR
(m) + EPBE,LR
X (m) + EPBE
C (1)
0.60 Å; Ti4+: Ref ¼ 0.605 Å),28 Sb can be easily tted at the Ti lattice
site without causing major lattice distortion in the host crystal where, ‘a’ stands for the mixing coefficient. The screening
structure; (b) introduction of Sb into the N-doped system is parameter m denes the short ranged (SR) and long ranged (LR)
expected to maintain charge balance in the codoped SrTiO3 (Ti4+ + part of the interaction. In the present study, we use the mixing
O2 ¼ Sb5+ + N3); (c) use of a non transition metal (Sb) may be exchange parameter of 28% and standard screening parameter
more preferable over transition metal to avoid localized ‘d’- of 0.2 Å1 to reproduce the experimental band gap of SrTiO3.11
states in the forbidden region. It has been successfully shown in K-point mesh of 3 3 3 was set for the hybrid functional
the earlier reports that introduction of Sb as a codopant calculations. The valence states considered during the calcula-
signicantly increases the photocatalytic property of the mon- tions are: Sr (4s24p65s2), Ti (4s23d2), Sb (5s25p3), O (2s22p4), and
odoped entity.16,29–36 As for example, codoping of Sb into the Rh- N (2s22p3). The energy cutoff of 600 eV has been chosen for the
doped SrTiO3 effectively reduces the electron–hole recombi- plane wave basis sets. To investigate the absorption behavior
nation rate by xing Rh oxidation state to +3 forming a charge frequency-dependent dielectric function calculation has been
compensated system.16 Similarly, codoping of Sb into the Cr- carried out for the codoped and undoped SrTiO3.
doped SrTiO3 has been found to enhance the photoactivity for
H2 evolution by suppressing the formation of hexavalent 3. Results and discussion
chromium (Cr6+) as well as vacancy, which are known to
diminish the photoactivity of SrTiO3.29 However, detailed Salient features of the electronic structure of undoped SrTiO3,
theoretical studies exploring the effect of Sb on the electronic monodoped SrTiO3 using N and Sb as dopants, and (Sb, N)-
structure of SrTiO3 are rare. Here, we present a systematic study codoped SrTiO3 are discussed below.
to investigate the synergistic effect of both N (as a substituent of
O) and Sb (as a substituent of Ti) as codopants, and compare the 3.1. SrTiO3
results with that of the N-doped, Sb-doped and undoped
Cubic SrTiO3 is described by a framework of TiO6 perfect
SrTiO3, using density functional theory as a tool. Since the use
octahedron, where the center position is occupied by Ti.
of Heyd, Scuseria, and Ernzerhof (HSE)37 functional is known to
Fig. 1a shows the band structure plot for the undoped SrTiO3
successfully reproduce the experimental band gap for undoped
along high symmetric k-path of the Brillouin zone (the hori-
SrTiO3, we employ the same functional in our calculations.
zontal dashed line represents the Fermi level). The calculated
Conclusions have been drawn by analyzing the band structure,
band gap (3.19 eV) shows excellent agreement with the
partial density of states (PDOS), optical spectrum, and band
experimental value (3.2 eV).5 Analysis of the projected density
alignment with respect to the water redox levels.
of states (Fig. 2) indicates that the valence band maximum
(VBM) is dominated by O 2p states while conduction band
2. Computational methods minimum (CBM) is composed of Ti 3d states. We will now
discuss the effect of doping and codoping on the electronic
Spin-polarized DFT calculations have been carried out using
structure of SrTiO3.
Vienna ab initio simulation package (VASP)38 electronic structure
code. The electron-ion interaction has been treated by projector
augmented wave (PAW) potential,39 which is widely used for the 3.2. N-doped SrTiO3
studies of electronic structure of semiconductor materials. For Since, doping of nitrogen in the oxide based material mainly
the Brillouin zone integration we employ Monkhorst and Pack involves replacement of oxygen by nitrogen, we consider only
scheme to generate gamma-centered k-point sets.40 The mono- substitutional doping. From the optimized cell structure,
doped system is modeled by replacing either one of the Ti or O obtained in our calculation, it is found that N doping does not
atoms by Sb or N, respectively from a 2 2 2 supercell (40 make any signicant change in the SrTiO3 crystal structure. The
atoms) of the cubic (space group: Pm3m) SrTiO3 crystal structure. Ti–N bond length (1.976 Å) in N-doped SrTiO3 is slightly higher
In the case of codoping, we simultaneously introduce both Sb than the Ti–O bond length (1.974 Å) in undoped crystal. This
and N into SrTiO3. This corresponds to Sb concentration of leads to increase in lattice parameter by only a small extent.
45704 | RSC Adv., 2014, 4, 45703–45709 This journal is © The Royal Society of Chemistry 2014
View Article Online
Fig. 2 DOS and PDOS plots for undoped SrTiO3. The vertical dashed
lines indicate the Fermi level.
This journal is © The Royal Society of Chemistry 2014 RSC Adv., 2014, 4, 45703–45709 | 45705
View Article Online
Bader charges on Ti and O centers are 2.04|e| and 1.23|e|, formation energy in each case.45,46 The defect formation energy
respectively, while, for Sb-doped SrTiO3, the values on the Sb and relates chemical potential of the elements as
O centers are 2.80|e| and 1.29|e|, respectively.
DHf ¼ Edoped + nOmO + nTimTi ESrTiO3 nNmN nSbmSb (3)
3.4. (Sb, N)-codoped SrTiO3 where, Edoped represents total energy for the doped/codoped
Published on 10 September 2014. Downloaded by Christian Albrechts Universitat zu Kiel on 23/10/2014 22:40:24.
where, EN-SrTiO3, ESb-SrTiO3, E(Sb,N)-SrTiO3 and ESrTiO3 represent the mSr + mTi + 3mO ¼ mSrTiO3(bulk) (4)
energy of the N-doped, Sb-doped, (Sb, N)-codoped, and undo-
The chemical potential of the element cannot exceed that of
ped SrTiO3 supercell, respectively. The positive value of the
the respect bulk (Sr/Ti) or gaseous (O) state, i.e. mSr # mSr(bulk),
defect pair binding energy (1.58 eV per supercell) indicates that
mTi # mTi(bulk) and mO # mO(gas).
the codoped system is sufficiently stable.
The expression for the heat of formation SrTiO3 is given by
Let us discuss the synergistic effect of both N and Sb on the
electronic structure of SrTiO3. The band structure plot (Fig. 1d)
D ¼ mSrTiO3(bulk) mSr(bulk) mTi(bulk) 3mO(gas) (5)
for the codoped SrTiO3 shows that the Fermi level resides above
the VBM, similar to that of an intrinsic semiconductor. It is Since the heat of formation for SrTiO3 is a negative quantity
interesting to observe that the band gap reduces signicantly to we can write
2.66 eV, which ensures the absorption of visible light. This is the
consequence of elevation of VB edge associated with the mSr(bulk) + D # mSr # mSr(bulk) (6a)
formation of (N 2p, O 2p) hybridized state (Fig. 4). The discrete
acceptor states as appearing in the case of N-doped SrTiO3 are mTi(bulk) + D # mTi # mTi(bulk) (6b)
completely passivated, resulting into a continuum band struc-
ture. This is due to the presence of Sb which compensates for 3mO(gas) + D # 3mO # 3mO(gas) (6c)
the electron deciency introduced through N-doping. It has
The chemical potential, mSr(bulk), mTi(bulk), and mSb have been
been quantied by Bader charge analysis, which indicates a
calculated from the energy of an atom in the respective bulk
higher electronic charge on the N centre in the (Sb, N)-codoped
structure. On the other hand, mO(gas) and mN have been calcu-
SrTiO3 (1.36|e|) as compared to that in N-doped SrTiO3
lated from the energy of an oxygen atom (1/2EO2) and nitrogen
(1.15|e|). On the other hand, the Ti 3d defect states arising in
atom (1/2EN2) in the corresponding gaseous molecule conned
the case of Sb-doped SrTiO3 is no longer present in the (Sb, N)-
at the centre of a 20 20 20 Å3 cubic box.
codoped SrTiO3. Presence of both Sb and N leads to formation
The calculation of the formation energy (Fig. 5a and b) indi-
of charge compensated system, which will minimize formation
cates that in both the cases substitutional doping is energetically
of undesired defects and thus reduce charge carrier loss.
more favored in comparison to the interstitial doping. Fig. 5a
shows the variation of formation energy for N-doped SrTiO3 as a
3.5. Defect formation energy function of oxygen chemical potential (mO mO(gas)). Under O-
rich condition the formation energy has a large positive value,
In order to nd favorable growth condition for doping of the
which shows a decreasing trend as we proceed towards O-poor
individual element as well as both, we calculate defect
condition. As shown in Fig. 5b, that the formation energy for
Sb-doping is highly positive under Ti-rich condition. This indi-
cates that the doping with individual element is highly unfa-
vorable under host rich condition, may be due to requirement of
large energy for vacancy formation in the substitution process.
However, this is found to be feasible in presence of both the
dopant elements as shown in Fig. 5c. The calculation indicates
that the O-rich condition is more suitable than the Ti-rich
condition for the growth of the (Sb, N)-codoped SrTiO3.
45706 | RSC Adv., 2014, 4, 45703–45709 This journal is © The Royal Society of Chemistry 2014
View Article Online
Fig. 5 The variation of defect formation energy with the chemical potential of O (mO mO-gas) and Ti (mTi mTi-bulk) for (a) N-doped SrTiO3 (b) Sb-
doped SrTiO3 (c) (Sb, N)-codoped SrTiO3. The color lines (5c) corresponds to different formation energies for the (Sb, N)-codoped SrTiO3. The
codoped system cannot be formed in the lower half of the diagonal axis (large white region) in the case of 5c because the chemical potential of Sr
(mSr) is exceeded (eqn (6a)).
This journal is © The Royal Society of Chemistry 2014 RSC Adv., 2014, 4, 45703–45709 | 45707
View Article Online
seen from Fig. 7, that the VBM of (Sb, N)-codoped SrTiO3 is 2 2 2 (Fig. 5c) and 2 2 3 supercell (Fig. 8) has been
located 0.20 eV below the water oxidation level and the CBM found to be almost similar. This justies the choice of supercell
position is 1.23 eV above the water reduction level. Hence, size used in the calculation.
overall water splitting is thermodynamically feasible in the
(Sb, N)-codoped SrTiO3. Interestingly, the CBM of the (Sb, N)-
SrTiO3 is located even 0.43 eV above the CBM of the undoped 4. Conclusion
Published on 10 September 2014. Downloaded by Christian Albrechts Universitat zu Kiel on 23/10/2014 22:40:24.
SrTiO3. This implies that the modied system should be more The hybrid DFT calculations reported here reveal the signi-
feasible for H2 evolution due to enhancement of reducing cance of codoping of Sb and N into SrTiO3 for improving the
power with respect to the undoped SrTiO3. photocatalytic activity under visible light. The positive defect
pair binding energy indicates that the (Sb, N)-codoped SrTiO3 is
3.8. Effect of supercell size sufficiently stable. The band gap is reduced signicantly
without introducing any isolated midgap states. As the codoped
To investigate the effect of supercell size, we have also consid- system is charge compensated, it is expected to reduce unde-
ered larger supercell (2 2 3). The formation energy has been sirable vacancy formation, and consequently defect assisted
calculated following the same procedure as discussed in the carrier loss can be avoided. Presence of Sb also facilitates the
case of 2 2 2 supercell. Fig. 8 shows the formation energy introduction of N into the SrTiO3 crystal structure by reducing
variation as a function of chemical potential of Ti and O. The formation energy. The band alignment for the codoped SrTiO3
calculated formation energy for (Sb, N)-codoped SrTiO3 using is such that both the photo-oxidation and photo-reduction
processes associated with water splitting are thermodynami-
cally feasible. Moreover, the reducing power of the modied
SrTiO3 is higher than that of the undoped material, making it
more suitable for H2 generation. Hence we can propose that the
codoping of Sb and N into the SrTiO3 crystal is one of the
promising approaches to enhance the efficiency for photo-
splitting of water using visible light.
Acknowledgements
We thank the BARC computer centre for providing the high
performance parallel computing facility. The authors like to
acknowledge Dr A. K. Samanta, Dr C. Majumder and Dr D.
Karmakar for valuable discussions. We also thank Dr B. N.
Jagatap for his encouragement and support. The work of Prof.
Swapan K. Ghosh is supported through Sir J. C. Bose Fellowship
Fig. 7 The band alignment of the undoped, and codoped SrTiO3 with from the Department of Science and Technology, India.
respect to the water redox levels.
References
1 A. Kudo and Y. Miseki, Chem. Soc. Rev., 2009, 38, 253–278.
2 X. Chen, S. Shen, L. Guo and S. S. Mao, Chem. Rev., 2010, 110,
6503–6570.
3 K. Maeda and K. Domen, J. Phys. Chem. Lett., 2010, 1, 2655–
2661.
4 T. Hisatomi, J. Kubota and K. Domen, Chem. Soc. Rev., 2014,
DOI: 10.1039/c3cs60378d.
5 M. Cardona, Phys. Rev., 1965, 140, A651.
6 H. Irie, Y. Maruyama and K. Hashimoto, J. Phys. Chem. C,
2007, 111, 1847–1852.
7 J. W. Liu, G. Chen, Z. H. Li and Z. G. Zhang, J. Solid State
Chem., 2006, 179, 3704–3708.
Fig. 8 The variation of defect formation energy with the chemical 8 D. Wang, J. Ye, T. Kako and T. Kimura, J. Phys. Chem. B, 2006,
potential of O (mO mO-gas) and Ti (mTi mTi-bulk) for (Sb, N)-codoped 110, 15824–15830.
SrTiO3 using 2 2 3 supercell. The color lines corresponds to
9 W. Wei, Y. Dai, H. Jin and B. Huang, J. Phys. D: Appl. Phys.,
different formation energies for the (Sb, N)-codoped SrTiO3. The
codoped system cannot be formed in the lower half of the diagonal 2009, 42, 055401.
axis (large white region) because the chemical potential of Sr (mSr) is 10 H. Yu, S. Ouyang, S. Yan, Z. Li, T. Yu and Z. Zou, J. Mater.
exceeded (eqn (6a)). Chem., 2011, 21, 11347–11351.
45708 | RSC Adv., 2014, 4, 45703–45709 This journal is © The Royal Society of Chemistry 2014
View Article Online
11 P. Reunchan, N. Umezawa, S. Ouyang and J. Ye, Phys. Chem. 30 R. Niishiro, R. Konta, H. Kato, W.-J. Chun, K. Asakura and
Chem. Phys., 2012, 14, 1876–1880. A. Kudo, J. Phys. Chem. C, 2007, 111, 17420–17426.
12 R. Konta, T. Ishii, H. Kato and A. Kudo, J. Phys. Chem. B, 31 T. Ikeda, T. Nomoto, K. Eda, Y. Mizutani, H. Kato, A. Kudo
2004, 108, 8992–8995. and H. Onishi, J. Phys. Chem. C, 2008, 112, 1167–1173.
13 S. W. Bae, P. H. Borse and J. S. Lee, Appl. Phys. Lett., 2008, 92, 32 C. Di Valentin, G. Pacchioni, H. Onishi and A. Kudo, Chem.
104107. Phys. Lett., 2009, 469, 166–171.
Published on 10 September 2014. Downloaded by Christian Albrechts Universitat zu Kiel on 23/10/2014 22:40:24.
14 H.-C. Chen, C.-W. Huang, J. C. S. Wu and S.-T. Lin, J. Phys. 33 M. Niu, W. Xu, X. Shao and D. Cheng, Appl. Phys. Lett., 2011,
Chem. C, 2012, 116, 7897–7903. 99, 203111–203113.
15 S. Kawasaki, K. Akagi, K. Nakatsuji, S. Yamamoto, 34 P. Reunchan, S. Ouyang, N. Umezawa, H. Xu, Y. Zhang and
I. Matsuda, Y. Harada, J. Yoshinobu, F. Komori, J. Ye, J. Mater. Chem. A, 2013, 1, 4221–4227.
R. Takahashi, M. Lippmaa, C. Sakai, H. Niwa, M. Oshima, 35 R. Asai, H. Nemoto, Q. Jia, K. Saito, A. Iwase and A. Kudo,
K. Iwashina and A. Kudo, J. Phys. Chem. C, 2012, 116, Chem. Commun., 2014, 50, 2543–2546.
24445–24448. 36 R. Niishiro, S. Tanaka and A. Kudo, Appl. Catal., B, 2014,
16 K. Furuhashi, Q. Jia, A. Kudo and H. Onishi, J. Phys. Chem. C, 150–151, 187–196.
2013, 117, 19101–19106. 37 J. Paier, M. Marsman, K. Hummer, G. Kresse, I. C. Gerber
17 J. Wang, S. Yin, M. Komatsu, Q. Zhang, F. Saito and T. Sato, J. and J. G. Ángyán, J. Chem. Phys., 2006, 124, 154709.
Mater. Chem., 2003, 13, 2348–2352. 38 G. Kresse and D. Joubert, Phys. Rev. B: Condens. Matter Mater.
18 H. W. Kang and S. B. Park, Chem. Eng. Sci., 2013, 100, 384– Phys., 1999, 59, 1758.
391. 39 P. E. Blöchl, Phys. Rev. B: Condens. Matter Mater. Phys., 1994,
19 N. Li and K. L. Yao, AIP Adv., 2012, 2, 032135. 50, 17953.
20 C. Zhang, Y. Jia, Y. Jing, Y. Yao, J. Ma and J. Sun, Comput. 40 H. J. Monkhorst and J. D. Pack, Phys. Rev. B: Solid State, 1976,
Mater. Sci., 2013, 79, 69–74. 13, 5188.
21 J. Wang, S. Yin, M. Komatsu, Q. Zhang, F. Saito and T. Sato, J. 41 J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson,
Photochem. Photobiol., A, 2004, 165, 149–156. M. R. Pederson, D. J. Singh and C. Fiolhais, Phys. Rev. B:
22 Y. Y. Mi, S. J. Wang, J. W. Chai, J. S. Pan, C. H. A. Huan, Condens. Matter Mater. Phys., 1992, 46, 6671.
Y. P. Feng and C. K. Ong, Appl. Phys. Lett., 2006, 89, 231922. 42 H. H. Wang, F. Chen, S. Y. Dai, T. Zhao, H. B. Lu, D. F. Cui,
23 P. Liu, J. Nisar, B. Pathak and R. Ahuja, Int. J. Hydrogen Y. L. Zhou, Z. H. Chen and G. Z. Yang, Appl. Phys. Lett., 2001,
Energy, 2012, 37, 11611–11617. 78, 1676–1678.
24 M. Miyauchi, M. Takashio and H. Tobimatsu, Langmuir, 43 G. Henkelman, A. Arnaldsson and H. Jónsson, Comput.
2004, 20, 232–236. Mater. Sci., 2006, 36, 254–360.
25 J. Wang, S. Yin, M. Komatsu and T. Sato, J. Eur. Ceram. Soc., 44 W.-J. Yin, S.-H. Wei, M. M. Al-Jassim and Y. Yan, Phys. Rev.
2005, 25, 3207–3212. Lett., 2011, 106, 066801.
26 W. Wei, Y. Dai, M. Guo, L. Yu, H. Jin, S. Han and B. Huang, 45 R. Long and N. J. English, Chem. Mater., 2010, 22, 1616–1623.
Phys. Chem. Chem. Phys., 2010, 12, 7612–7619. 46 W.-J. Shi and S.-J. Xiong, Phys. Rev. B: Condens. Matter Mater.
27 H. Yu, S. Yan, Z. Li, T. Yu and Z. Zou, Int. J. Hydrogen Energy, Phys., 2011, 84, 205210.
2012, 37, 12120–12127. 47 M. Fox, Optical Properties of Solids, Oxford University Press,
28 R. D. Shannon, Acta Crystallogr., Sect. A: Cryst. Phys., Diffr., USA, 2002.
Theor. Gen. Crystallogr., 1976, 32, 751–767. 48 Y. Xu and M. A. A. Schoonen, Am. Mineral., 2000, 85, 543.
29 H. Kato and A. Kudo, J. Phys. Chem. B, 2002, 106, 5029–5034.
This journal is © The Royal Society of Chemistry 2014 RSC Adv., 2014, 4, 45703–45709 | 45709