Porous Graphitic Carbon Nitride: A Possible Metal-free Photocatalyst

Oct 22, 2014 - Hydrogen generation through photocatalytic water splitting with the aid of renewable solar energy is an important step toward the devel...
0 downloads 0 Views 860KB Size
Subscriber access provided by OHIOLINK ACCESS

Article

Porous Graphitic Carbon Nitride: A Possible Metal-free Photocatalyst for Water Splitting Kancharlapalli Srinivasu, Brindaban Modak, and Swapan K Ghosh J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/jp506538d • Publication Date (Web): 22 Oct 2014 Downloaded from http://pubs.acs.org on October 26, 2014

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 27

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

Porous Graphitic Carbon Nitride: A Possible Metalfree Photocatalyst for Water Splitting K. Srinivasu, Brindaban Modak and Swapan K. Ghosh* Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai – 400 085, India. Homi Bhabha National Institute and UM-DAE-Centre for Excellence in Basic Sciences, Mumbai, India

ACS Paragon Plus Environment

1

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 2 of 27

ABSTRACT. Hydrogen generation through photocatalytic water splitting with the aid of renewable solar energy is an important step towards the development of sustainable and alternative energy. In the present study, using the first principles calculations, we have explored the s-triazine based two-dimensional porous graphitic carbon nitride (g-CN) materials as potential photocatalyst for water splitting. For calculating the band structures more accurately, we have employed hybrid density functionals. The calculated band gap of the single layer g-CN is found to be 2.89 eV, which decreases to ~2.75 eV in multilayered structure. To improve the visible light activity, effect of doping with different non-metals on the electronic structure has been investigated. Among the different dopants studied, phosphorous is found to be more effective to reduce the band gap to 2.31 eV. The band edge potentials obtained from density functional calculations are corrected for vacuum potentials. The band alignments with respect to the water redox levels show that the thermodynamic criterion for the overall water splitting is satisfied. We have also carried out analogous studies on the heptazine based carbon nitride, g-C3N4 and the calculated band gaps as well as the position of the valance band maximum are consistent with the reported experimental results validating the computational method we have used.

Based on our

theoretical investigations, we can predict that the considered carbon nitride based materials should be a potential photocatalyst for water splitting under visible light.

KEYWORDS . Graphitic carbon nitride, Photocatalytic activity, Water splitting, Hybrid density functional, Substitutional doping.

ACS Paragon Plus Environment

2

Page 3 of 27

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. INTRODUCTION Search for sustainable, alternative energy sources has become an important area of research due to the world wide increasing energy demands in addition to the limited resources of the fossil fuels and their adverse environmental affects. Hydrogen is considered to be one of the best possible alternative energy carriers.1-2 However, the main challenges in developing hydrogen energy are its generation and storage. Hydrogen generation from water through photocatalytic water splitting using the solar energy is an ideal way which does not evolve any harmful byproducts.4 Following the development of photo electrochemical cell (PEC) by Fujishima and Honda5 using a semiconductor anode and a metal cathode, varieties of semiconductor based photocatalysts have been investigated for making an efficient catalyst to generate hydrogen from water. The mechanism of photocatalytic water splitting using semiconductor materials involves different steps, viz. generation of charge carriers (electron and holes) on light absorption, separation and diffusion of the photo generated charge carriers to the active sites on catalyst surface and further reaction (oxidation/reduction) with the adsorbed water to generate hydrogen and oxygen.3 Among them, TiO2 has been extensively investigated due to its potentiality to split water. However, its wide bandgap (~3.0 eV) limits its activity to the UV region of the solar spectrum.6 During past few years, extensive efforts have been dedicated to enhance the visible light activity of TiO2, mainly by doping with different metal of non-metal elements.6 Apart from TiO2, many other transition metal oxides, sulphides, selenides viz. ZnO, SrTiO3, NaTaO3 , ZnS, CdS, CdSe, etc. have been studied for photocatalytic water splitting.6-12

Other than these

conventional materials, polymeric semiconductor materials like graphitic carbon nitrides are also

ACS Paragon Plus Environment

3

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 4 of 27

found to have great potential as photocatalyst. Wang et al.13 have reported the melem-based graphitic carbon nitrides (g-C3N4) as a metal-free photocatalyst for visible-light driven water splitting. However, this material is reported to show very poor quantum yield which is attributed to the high recombination rate of the photo generated electron-hole pairs. To overcome this problem, many attempts have been made to tune the properties, viz. doping, metal decoration, introducing porosity, making graphene/g-C3N4 composites etc.14-17 Pan et al.18 have shown that the g-C3N4 nanotubes have better carrier mobilities as compared to the 2D sheets and further enhancement in photocatalytic activity upon functionalization of these tubes with metals like Pt and Pd is reported.

Wang et al.19 have observed enhanced photocatalytic activity for the

nonocomposites made of g-C3N4 and MoS2 in comparison to pure g-C3N4. Very recently, Shi et al.20 have synthesised the composites of mesoporous carbon and g-C3N4 and observed significant improvement in photocatalytic activity.

Apart from the g-C3N4, another polymeric carbon

nitride, the poly (triazine imide), is also found to be a potential photocatalyst for water splitting.21-23 In the present study, through first-principles based calculations, we have explored another form of graphitic carbon nitride with CN stoichiometry (g-CN) as a photocatalyst for water splitting under visible light. Though the g-C3N4 form of carbon nitride has been studied extensively, the number of studies on g-CN is very less. Cao et al.24-25 synthesized the g-CN nanotubes through the reaction of cyanuric chloride (C3N3Cl3) with sodium using NiCl2 as a catalyst. Different forms of g-CN, viz. nanotubes, nanoribbons and microspheres have been synthesized and their photoluminescence spectra have shown peaks centered around 430-450 nm, due to the band gap emission. Because of the porous nature of carbon nitride and their bio-compatibility, different carbon nitride systems are also studied for their application in biotechnology.26-27 Here, we have

ACS Paragon Plus Environment

4

Page 5 of 27

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

considered the free standing single layer as well as multi layers of g-CN. To further enhance the visible activity of g-CN, we have also studied the effect of substitutional doping with different non-metal elements like boron, oxygen, phosphorous etc. Analogous studies on g-C3N4 are also being carried out to compare the computational results with the existing experimental results. To the best of our knowledge, this is the first study to explore g-CN as a possible active photocatalyst for hydrogen generation through water splitting.

2. COMPUTATIONAL DETAILS: All our spin-polarized periodic density functional theory (DFT) based electronic structure calculations have been carried out using the projector augmented wave (PAW)28-29 potentials as implemented in the first principles based Vienna ab initio Simulation Package (VASP).30-31 Plane-wave basis sets with a kinetic energy cutoff of 550 eV have been used. The exchangecorrelation energy density functional, Exc [ρ] has been treated through the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE)32. An energy cutoff of 1 x 10-6 eV is set as the convergence criteria for the electronic self-consistent field iterations. The atomic positions are optimized until the maximum Hellmann−Feynman force on each atom is less than 0.01 eV Å-1. Automatically generated Gamma centered 9 × 9 × 1 Monkhorst−Pack set of kpoints were used to sample the Brillouin zone.33 As it is known that the band gaps calculated from PBE functional are underestimated, we have also used the more accurate hybrid functional, developed by Heyd, Scuseria, and Ernzerhof (HSE03 and HSE06)34 using the PBE optimized geometries to improve the accuracy of band structure calculations. All the initial geometries and the reported figures are generated using the graphical software XCrySDen and VESTA.35-36

ACS Paragon Plus Environment

5

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 27

3. RESULTS AND DISCUSSION: At first, we will look into the electronic structure of the single layer two-dimensional g-CN studied through the density functional theory. The optimized unit cell structures of the twodimensional g-CN considered is shown in Figure 1a. The optimized cell parameter from PBE method is found to be 7.127 Å with the C-C and C-N bond distances of 1.510 Å and 1.341 Å respectively. From the PBE calculations, g-CN is found to have a direct band gap of 1.54 eV. However, the band gap of single layer g-CN calculated through the HSE03 functional is found to be 2.89 eV which is more close to the experimental value25 as compared to the corresponding gap from HSE06 functional (3.18 eV). To verify the accuracy of the method, we have also calculated the band gap of g-C3N4 using both HSE03 and HSE06 methods and the band gap is found to be 2.82 eV and 3.0 eV respectively indicating that the HSE03 result is closer to the experimental value of 2.73 eV.37 Hence, here onwards we will discuss only the HSE03 results and the corresponding PBE and HSE06 results are given in Table 1. The band dispersion plot of g-CN single layer calculated through the HSE03 method is shown in Figures 1b. From the band structure, it is clear that both the valance and conduction bands are well dispersed and no localized states are present to act as recombination centers for the photo generated electron hole pairs. This is in contrast to the case of g-C3N4 for which localised band edges are reported in the band structure and claimed to be the reason for its poor efficiency observed.18 Hence the life time of the photo-generated charge carriers in g-CN is expected to be longer which can lead to better photocatalytic efficiency. We have also calculated the band decomposed charge density of the valance band maximum (VBM) and conduction band minimum (CBM) and the corresponding iso-surface plots are shown in Figures 1c and 1d respectively. From the isosurface plots, it is clear that the VBM mostly originates from the non-bonding nitrogen 2p states

ACS Paragon Plus Environment

6

Page 7 of 27

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

whereas the CBM is originating from the C=C anti-bonding π states (carbon pz states) and nitrogen pz states. As the measured band gap is 2.89 eV, it cannot utilize the major part of the solar spectrum and hence, it is required to engineer the band gap to enhance the visible light absorption efficiency of these materials. In one of the recent studies on g-C3N4, through the first principle calculations, Wu et al.38 have shown that the bilayer has much better visible light adsorption efficiency as compared to the single layer which is attributed to the interlayer coupling. To verify this effect on the g-CN, we have also studied the electronic structure of the bilayer as well as the trilayer of g-CN. For the bilayer, we have considered two different stacking modes, AA and AB. From the optimized structures, the AB stacking mode as shown in Figure 2a is found to be energetically more favourable as compared to the AA mode. The optimized inter layer separation from the PBE method is found to be 3.578 Å. However, as it is known that the van der Waals interactions are important in this kind of systems, we have reoptimized the structure using the dispersion corrected PBE (PBE-D2)39 method and the interlayer separation is observed to be 3.088 Å indicating the strong coupling between the two g-CN layers. The interlayer interaction energy per unit cell is found to be -0.57 eV. The calculated band structure plot of the bilayer is given in Figure 2b and the corresponding band gap is found to be 2.81 eV which is 0.08 eV less as compared to that of the single layer. This indicates that the band gap is reduced due to the interlayer coupling in g-CN. To see any further effect of inter-layer coupling on the electronic structure, we have also studied the electronic structure of trilayer g-CN. For trilayer, we have added one more g-CN sheet above the bilayer in two different ways, ABA and ABC. It has been found that the ABC mode of stacking as shown in Figure 2c is energetically more favourable with interlayer

ACS Paragon Plus Environment

7

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 27

separation of 3.029 Å and the corresponding inter-layer interaction energy per unit cell is calculated to be -1.19 eV. The calculated band dispersion plot is reported in Figure 2d and the corresponding band gap is found to be 2.75 eV which is 0.14 eV less as compared to the single layer band gap. To mimic the multi layer structure more accurately, we have considered another system with six layers of g-CN stacked in ABCABC manner. The calculated band gap of this system is found to be 2.72 eV which is very close to that of tri-layer g-CN (2.75 eV) indicating that it is nearly approaching the multilayer nature. These results indicate that though the band gap is getting lowered on forming multiple layers, the reduction is not to the extent to achieve appreciable visible light activity. Doping with foreign elements has been shown to be an effective method for tuning the electronic band structure of semiconductor materials due to the difference in the energy levels of the dopant element as compared to that of the host element.40-41 Here, to further reduce the band gap of g-CN, so that the material can absorb the major part of the solar spectrum, we have studied the effect of doping with different non-metal elements like boron, oxygen, sulphur, and phosphorous. In the case of boron doping, we have considered two possible structures such that in one of the structures boron is substituted in place of nitrogen while in the other, boron replaces the carbon atom. From the optimized energies, it is found that the boron substitution in place of carbon is energetically more favourable as compared to the case of replacement of nitrogen by boron. The band gap of the boron doped system is found to be 2.84 eV which is comparable to that of the undoped system and hence boron doping is not helpful in reducing the band gap. In the cases of oxygen, sulphur and phosphorous doping, substitution of nitrogen as shown in Figure 3a is found to be energetically more preferred as compared to replacement of the carbon. The calculated band gaps in oxygen and sulphur doped g-CN are found to be 2.81 eV and

ACS Paragon Plus Environment

8

Page 9 of 27

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

2.49eV respectively. However, both oxygen and sulphur doping is found to introduce isolated donor states below the conduction band which can act as recombination centre and hence can decrease the efficiency of the catalyst. However, in the case of phosphorous doped g-CN, the band gap is found to be reduced considerably to 2.31 eV without introducing any undesirable midgap states in the forbidden region and the corresponding band structure plot is shown in Figure 3b. These results on non-metal doping reveal that the phosphorous doping will be more suitable to bring down the band gap of g-CN without creating any undesirable midgap states. To verify the analogous doping effects on g-C3N4, we have studied the electronic structure of similar non-metal doped g-C3N4. We have verified the different possible sites for substitutional doping viz. three different nitrogen sites, one two coordinated nitrogen and two different three coordinated nitrogen sites as well as one carbon site which is three coordinated. In the cases of boron and phosphorous, substitution in place of carbon is found to be more preferred as compared to all three different nitrogen sites, whereas in the cases of sulphur and oxygen, nitrogen which is coordinated to two carbons (pyridine like nitrogen) is found to be energetically preferred site for substitution and these results are consistent with the earlier results on sulphur and phosphorous doping.16 The optimized geometry of g-C3N4 and calculated density of states for g-C3N4 along with all the doped systems are reported in Figures S1 and S2 (supporting information). The calculated band gap in the cases of boron doped system is found to be 3.24 eV, and it is found to create midgap states as shown in Figure S2 which can trap the photo-generated carriers thereby decreasing the efficiency. However, it can be observed that the conduction band becomes highly dispersed as compared to that of the undoped system. The calculated band gap in oxygen and sulphur doped g-C3N4 is found to be 2.98 eV and 3.15 eV respectively and the increase in bandgap upon sulphur doping is consistent with the experimental

ACS Paragon Plus Environment

9

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 10 of 27

results.37 As both O and S are n-type dopants, occupied donor states are found to be created close to the bottom of conduction band as it is shown in Figure S2. In the case of phosphorous doped g-C3N4, the calculated band gap is found to be 2.70 eV which is 0.12 eV less as compared to the undoped g-C3N4 and this result is consistent with the recent experimental observations of Hu et al.42 However, in this case also dopant donor states are fond below the conduction band. These results indicate that, except phosphorous all other doped systems are found to have wider band gap as compared to the pristine g-C3N4. To have a look at the optical absorption characteristics of these g-CN systems considered, we have calculated the complex frequency dependent dielectric ε (ω ) = ε1 (ω ) + iε 2 (ω ), function,

where ε1 and ε2 represent the real and imaginary parts of the dielectric function,

respectively. The absorption coefficient α(ω) has been obtained from the relation43

α (ω ) = 2ω

(

2

2

ε1 (ω ) + ε 2 (ω ) − ε1 (ω )

)

1

2

.

The calculated absorption spectra for the pure single layer g-CN system along with its boron and phosphorous substituted systems are shown in Figure 4. The undoped g-CN is found to have a strong absorption peak around 260-300 nm which is the characteristic of π to π* transition observed in the s-triazine compounds and the result is consistent with the experimental observations (260-280 nm).25 This indicates that the absorption of pure g-CN is limited to the UV region of the spectrum. In the case of boron doped g-CN, though the peak is broadened slightly as compared to pure g-CN, still it is in the UV region only. However, in the case of phosphorous doped g-CN, the absorption curve is considerably extended towards the visible region indicating its improved visible light activity as compared to the undoped and boron doped systems.

ACS Paragon Plus Environment

10

Page 11 of 27

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

Though improving the visible light activity is a necessary condition for photocatalytic water splitting using solar light, it is not sufficient. To generate hydrogen and oxygen from the water splitting, the band edges should be positioned appropriately with respect to the redox levels of water. As the position of the band edges obtained directly from DFT are known to be shifted when a pseudopotental based method is used along with the periodic boundary condition, we have corrected the band edges using standard method.44-45 We have calculated the averaged local potential along the z-direction (perpendicular to the g-CN layer) and the plot of the electrostatic potential energy as a function of vacuum thickness (distance) is shown in Figure 5. The plot shows that the energy increases and gets saturated which can be taken as the electrostatic potential energy in vacuumed and all the Eigen values are corrected with respect to this potential. To assure the accuracy of the method, we have first calculated the band edge potentials of the gC3N4 for which the experimental results are available. The calculated valence band maximum (VBM) for g-C3N4 is found to be 2.06 eV below the water reduction level which is very close to the experimental reported value (2.12 eV below the water reduction level) indicating the accuracy of the method. The conduction band minimum (CBM) of g-C3N4 is found to be located 0.76 eV above the water reduction potential as it has been shown in the band alignment plot in Figure 6. The corresponding band alignments for boron and phosphorous doped g-C3N4 are also shown in Figure 6 and it shows that the boron doped system is found to have more elevation of conduction band which is 1.19 eV above the water reduction. From these results on g-C3N4 it is clear that though the band gap is increased on B-doping, its conduction band becomes more dispersed which can result in better carrier mobilities and also its reduction capacity increases due to the considerable elevation of the conduction band. As shown in Figure 6, for the single layer g-CN, the CBM and VBM are found to be located at -4.24 eV and -7.13 eV respectively.

ACS Paragon Plus Environment

11

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 27

The CBM is located just (0.26 eV) above the water reduction level whereas the VBM is found to be 1.41 eV below the water oxidation potentials. In the case of multilayer systems, the VBM is located at around -7.2 eV and the CBM is around -4.3 to -4.4 eV. In the case of boron doped gCN, though the band gap does not alter much, both the conduction and valence band positions are found to be lifted as compared to the pure g-CN indicating improved reduction capacity. In the case of phosphorous doping, the change in CBM is found to be negligible but the VBM shift upward is significant as compared to the undoped system. Thus, in all the systems the CBM is located above the water reduction potential and the VBM is well below the water reduction potential. In both the cases of g-CN and g-C3N4, boron doping is found to lift the CBM thereby improving the reduction potential of the catalyst. These band alignments show that both the oxidation and reduction reactions of water are thermodynamically feasible in both g-C3N4 and gCN. Due to the better dispersion of band edges, the g-CN, especially its phosphorous doped counterparts with lower band gap of 2.31 eV can be expected to be a good photocatalyst material for hydrogen generation under visible light. As the doping with non-metal elements has been well explored in graphene and other forms of carbon nitride (g-C3N4), it is experimentally quite possible to synthesize this photocatalyst systems.37, 42, 46

4. CONCLUSIONS In summary, we have explored the possible photocatalytic activity of s-triazine based twodimensional graphitic carbon nitrides toward the water splitting under visible light. Single layer 2D g-CN system is found to have a direct band gap of 2.89 eV, while the multilayer g-CN materials are found to have slightly improved visible light absorption due to the interlayer coupling. Absence of localized states in the band dispersion plot indicates better charge carrier mobilities. The calculated band gap and band edge potentials of g-C3N4 are very close to the

ACS Paragon Plus Environment

12

Page 13 of 27

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

experimentally reported results indicating the accuracy of the method used. Effect of doping with non-metal elements like boron, oxygen, phosphorous, etc. on the electronic band structure of g-CN as well as g-C3N4 has been investigated. Substitutional doping with phosphorous in gCN is found to improve the visible light absorption significantly by bringing down the band gap to 2.31 eV without introducing any midgap states. Our results on sulphur and phosphorous doped g-C3N4 are matching well with the reported experimental results. The calculated optical absorption spectra show significant enhancement in visible light absorption efficiency on phosphorous doping. The positions of band edges in pure g-CN as well as phosphorous doped counterpart with respect to water redox levels are found to satisfy the thermodynamic criteria for overall water splitting. ASSOCIATED CONTENT

Supporting Information. Optimized geometry if g-C3N4, DOS plots for g-C3N4 and its nonmetal doped counter parts are reported. This material is available free of charge via the Internet at http://pubs.acs.org.

Corresponding Author *Email: [email protected] Phone: 91-22-25595092

ACKNOWLEDGMENT We thank the BARC computer center for providing the high performance parallel computing facility. We also thank Dr B.N. Jagatap for his encouragement and support. The work of S.K.G. is supported through Sir J.C. Bose Fellowship from the Department of Science and Technology, India.

ACS Paragon Plus Environment

13

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 14 of 27

REFERENCES (1) Lewis, N. S.; Nocera, D. G. Powering the Planet: Chemical Challenges in Solar Energy Utilization. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 15729–15735. (2) Züttel, A.; Remhof, A.; Borgschulte, A.; Friedrichs, O. Hydrogen: The Future Energy Carrier. Philos. Trans. R. Soc., A 2010, 368, 3329−3342. (3) Yerga, R. M. N.; Galván, M. C. A.; del Valle, F.; de la Mano, J. A. V.; Fierro, J. L. G.; Water

Splitting

on

Semiconductor

Catalysts

under

Visible-Light

Irradiation.

ChemSusChem 2009, 2, 471 – 485. (4) Bard, A. J.; Fox, M. A. Artificial photosynthesis: Solar Splitting of Water to Hydrogen and Oxygen. Acc. Chem. Res. 1995, 28, 141-145. (5) Fujishima, A.; Honda, K. Electrochemical Photolysis of Water at a Semiconductor Electrode. Nature 1972, 238, 37-38. (6) K. Hashimoto, H. Irie, A. Fujishima, TiO2 Photocatalysis: A Historical Overview and Future Prospects. Jpn. J. Appl. Phys. Part 2005, 44, 8269-8285. (7) Chen, X.; Shen, S.; Guo, L.; Mao, S. S. Semiconductor-based Photocatalytic Hydrogen Generation. Chem. Rev. 2010, 110, 6503–6570. (8) Kudo, A.; Miseki, Y. Heterogeneous Photocatalyst Materials for Water Splitting. Chem. Soc. Rev., 2009, 38, 253–278. (9) Maeda, K.; Domen, K. Photocatalytic Water Splitting: Recent Progress and Future Challenges. J. Phys. Chem. Lett. 2010, 1, 2655–2661. (10)

Maeda, K.; Domen, K. New Non-Oxide Photocatalysts Designed for Overall

Water Splitting under Visible Light. J. Phys. Chem. C 2007, 111, 7851-7861.

ACS Paragon Plus Environment

14

Page 15 of 27

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

(11) Modak, B.; Srinivasu, K.; Ghosh, S. K. Photocatalytic Activity of NaTaO3 Doped with N, Mo, and (N,Mo): A Hybrid Density Functional Study. J. Phys. Chem. C, 2014, 118, 10711–10719. (12) Liu, P.; Nisar, J.; Ahuja, R.; Pathak, B. Layered Perovskite Sr2Ta2O7 for Visible Light Photocatalysis: A First Principles Study. J. Phys. Chem. C, 2013, 117, 5043–5050. (13) Wang, X.; Maeda, K.; Thomas, A.; Takanabe, K.; Xin, G.; Carlsson, J. M.; Domen, K.; Antonietti, M. A Metal-free polymeric Photocatalyst for Hydrogen Production from Water under Visible Light. Nat. Mater. 2008, 8, 76-80. (14) Zhang, J.; Chen, X.; Takanabe, K.; Maeda, K.; Domen, K.; Epping, J. D.; Fu, X.; Antonietti, M.; Wang, X. Synthesis of a Carbon Nitride Structure for Visible-Light Catalysis by Copolymerization. Angew. Chem. Int. Ed. 2010, 49, 441 –444. (15) Zhang, Y.; Mori, T.; Niu, L.; Ye, J.; Non-covalent Doping of Graphitic Carbon Nitride Polymer with Graphene: Controlled Electronic Structure and Enhanced Optoelectronic Conversion. Energy Environ. Sci., 2011, 4, 4517-4521. (16) Ma, X.; Lv, Y.; Xu, J.; Liu, Y.; Zhang, R.; Zhu, Y. A Strategy of Enhancing the Photoactivity of g‑C3N4 via Doping of Nonmetal Elements: A First-Principles Study. J. Phys. Chem. C, 2012, 116, 23485–23493. (17) Zheng, Y.; Liu, J.; Liang, J.; Jaroniec, M.; Qiao, S. Z. Graphitic Carbon Nitride Materials: Controllable Synthesis and Applications in Fuel Cells and Photocatalysis. Energy Environ. Sci., 2012, 5, 6717-6731. (18) Pan, H.; Zhang, Y. W.; Shenoy, V. B.; Gao, H. Ab Initio Study on a Novel Photocatalyst: Functionalized Graphitic Carbon Nitride Nanotube. ACS Catal. 2011, 1, 99–104.

ACS Paragon Plus Environment

15

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 16 of 27

(19) Wang, J.; Guan, Z.; Huang, J.; Li, Q.; Yang, J. Enhanced Photocatalytic Mechanism for the Hybrid g-C3N4/MoS2 Nanocomposite. J. Mater. Chem. A, 2014, 2, 7960–7966. (20) Shi, L.; Liang, L.; Ma, J.; Wang, F.; Sun, J. Remarkably Enhanced Photocatalytic Activity of Ordered Mesoporous Carbon/g-C3N4 Composite Photocatalysts Under Visible Light. Dalton Trans., 2014, 43, 7236–724. (21) Schwinghammer, K.; Tuffy, B.; Mesch, M. B.; Wirnhier, E.; Martineau, C.; Taulelle, F.; Schnick, W.; Senker, J.; Lotsch, B. V. Triazine-based Carbon Nitrides for Visible-LightDriven Hydrogen Evolution. Angew.Chem.Int.Ed. 2013, 52, 2435 –2439. (22) McDermott, E. J.; Wirnhier, E.; Schnick, W.; Virdi, K. S.; Scheu, C.; Kauffmann, Y.; Kaplan, W. D.; Kurmaev, E. Z.; Moewes, A. Band Gap Tuning in Poly(triazine imide), a Nonmetallic Photocatalyst. J. Phys. Chem. C, 2013, 117, 8806–8812. (23) Schwinghammer, K.; Mesch, M. B.; Duppel, V.; Ziegler, C.; Senker, J.; Lotsch, B. V. Crystalline Carbon Nitride Nanosheets for Improved Visible-Light Hydrogen Evolution. J. Am. Chem. Soc., 2014, 136, 1730–1733. (24) Cao, C.; Huang, F.; Cao, C.; Li, J.; Zhu, H. Synthesis of Carbon Nitride Nanotubes via a Catalytic-Assembly Solvothermal Route. Chem. Mater. 2004, 16, 5213−5215. (25) Li, J.; Cao, C.; Hao, J.; Qiu, H.; Xu, Y.; Zhu, H. Self-assembled One-dimensional Carbon Nitride Architectures. Diamond Relat. Mater. 2006, 15, 1593−1600. (26) Li, X.; Zhou, J.; Wang, Q.; Kawazoe, Y.; Jena, P. Patterning Graphitic C−N Sheets into a Kagome Lattice for Magnetic Materials. J. Phys. Chem. Lett. 2013, 4, 259−263. (27) Li, X.; Zhang, S.; Wang, Q. Stability and Physical Properties of a Tri-ring Based Porous g-C4N3 sheet. Phys. Chem. Chem. Phys., 2013, 15, 7142--7146.

ACS Paragon Plus Environment

16

Page 17 of 27

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

(28) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B. 1994, 50, 17953– 17979. (29) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector AugmentedWave Method. Phys. Rev. B. 1999, 59, 1758–1775. (30) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B. 1996, 54, 11169-11186. (31) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mat. Sci. 1996, 6, 15-50. (32) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (33) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-zone Integrations. Phys. Rev. B 1976, 13, 5188-5192. (34) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207. (35) Kokalj, A. XCrySDen—a New Program for Displaying Crystalline Structures and Electron Densities. J. Mol. Graph. Model. 1999, 17, 176-179 (36) Momma, K.; Izumi, F. VESTA: A Three-dimensional Visualization System for Electronic and Structural Analysis. J. Appl. Cryst. 2008, 41, 653-658 (37) Liu, G.; Niu, P.; Sun, C.; Smith, S. C.; Chen, Z.; Lu, G. Q.; Cheng, H. Unique Electronic Structure Induced High Photoreactivity of Sulfur-Doped Graphitic C3N4. J. Am. Chem. Soc. 2010, 132, 11642–11648.

ACS Paragon Plus Environment

17

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 27

(38) Wu, F.; Liu, Y.; Yu, G.; Shen, D.; Wang, Y.; Kan, E. Visible-Light-Absorption in Graphitic C3N4 Bilayer: Enhanced by Interlayer Coupling. J. Phys. Chem. Lett. 2012, 3, 3330−3334. (39) Grimme, S. Semiempirical GGA-type Density Functional Constructed with a LongRange Dispersion Correction. J. Comput. Chem., 2006, 27, 1787-1799. (40) Modak, B.; Srinivasu, K.; Ghosh, S. K. Band Gap Engineering of NaTaO3 Using Density Functional Theory: A Charge Compensated Codoping Strategy. Phys. Chem. Chem. Phys., 2014, 16, 17116-17124. (41) Modak, B.; Srinivasu, K.; Ghosh, S. K. Improving Photocatalytic Properties of SrTiO3 Through (Sb, N) Codoping: A Hybrid Density Functional Study. RSC Adv., 2014, 4, 45703–45709. (42) Hu, S.; Ma, L.; You, J.; Li, F.; Fan, Z.; Wang, F.; Liu, D.; Gui, J. A Simple and Efficient Method to Prepare Phosphorus Modified g-C3N4 Visible Light Photocatalyst. RSC Adv., 2014, 4, 21657–21663. (43) Fox, M. Optical Properties of Solids, Oxford University Press, USA, 2002. (44) 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–16654. (45) Kang, J.; Tongay, S.; Zhou, J.; Li, J.; Wu, J. Band Offsets and Heterostructures of Two-dimensional Semiconductors. Appl. Phys. Lett. 2013, 102, 012111. (46) Choi, C. H.; Park, S. H.; Woo, S. I. Binary and Ternary Doping of Nitrogen, Boron, and Phosphorus into Carbon for Enhancing Electrochemical Oxygen Reduction Activity. ACS Nano 2012, 6, 7084−7091.

ACS Paragon Plus Environment

18

Page 19 of 27

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

Figure Captions Figure 1. (a) Optimized geometries of the g-CN along with its (b) band structure plot using HSE03 methods and the band decomposed charge density plot of (c) highest valance and (d) lowest conduction states

Figure 2. Optimized supercell structures of (a) double layer g-CN with AB stacking and (c) triple layer with ABC stacking along with the band structure of (b) double layer and (d) triple layer calculated through HSE03 functionals.

Figure 3. (a) Optimized unit cell structure and (b) band structure plot (from HSE03) of phosphorous doped g-CN

Figure 4. Calculated optical absorption spectra of different g-CN systems considered

Figure 5: Plot of averaged electrostatic potential as a function of vacuum thickness for the g-CN (a) single layer, (b) bilayer and (c) trilayer.

Figure 6. Calculated VBM and CBM positions of g-CN and its boron and phosphorous doped counter parts.

ACS Paragon Plus Environment

19

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 27

Table 1. The calculated band gap results of different g-CN systems using PBE, HSE03 and HSE06 methods.

System

Band gap (eV) PBE

HSE03

HSE06

Single layer

1.54

2.89

3.18

Bilayer

1.49

2.81

3.10

Trilayer

1.45

2.75

3.03

six-layer

1.41

2.72

2.99

g-CN-B

1.71

2.84

3.11

g-CN-P

1.15

2.31

2.52

ACS Paragon Plus Environment

20

Page 21 of 27

(a)

(b)

8 6

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

4 2 0 -2 -4

Κ

Γ

(c)

Μ

Γ

(d)

Figure 1. (a) Optimized geometries of the g-CN along with its (b) band structure plot using HSE03 methods and the band decomposed charge density plot of (c) highest valance and (d) lowest conduction states.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

(c)

(a)

8

(b)

Energy (eV)

6

4 2 0 -2 -4Γ Γ

(d)

8

6

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

Page 22 of 27

4 2 0

-2

Κ

Μ

Γ

-4Γ Γ

Κ

Μ

Γ

Figure 2. Optimized supercell structures of (a) double layer g-CN with AB stacking and (c) triple layer with ABC stacking along with the band structure of (b) double layer and (d) triple layer calculated through HSE03 functionals.

ACS Paragon Plus Environment

22

Page 23 of 27

(a)

(b)

8 6

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

4 2 0 -2 -4

Γ

Κ

Μ

Γ

Figure 3. (a) Optimized unit cell structure and (b) band structure plot (from HSE03) of phosphorous doped g-CN

ACS Paragon Plus Environment

23

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

Absorp. coeff (arb.unit)

The Journal of Physical Chemistry

200

Page 24 of 27

g-CN B@g-CN P@g-CN

300

400

500

600

700

Wavelength (nm) Figure 4. Calculated optical absorption spectra of different g-CN systems considered

ACS Paragon Plus Environment

24

Page 25 of 27

5

0

Energy (eV)

0

-5

(a) -10

-15

-5

(b)

-10 -15

-20

-20

0

5

10

15

20

25

0

5

0

10

15

20

25

30

35

0

Distance (A )

Distance (A )

5 0

Energy (eV)

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

-5

(c)

-10 -15 -20 0

5

10

15

20

25

30

35

0

Distance (A )

Figure 5: Plot of averaged electrostatic potential as a function of vacuum thickness for the g-CN (a) single layer, (b) bilayer and (c) trilayer.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

-3.0

B @ g -C 3 N 4 P @ g-C 3 N 4

g -C 3 N 4

g -C N

-4.5

B @ g -C N P @ g -C N

2.31 eV

2.85 eV

2.89 eV

3.24 eV

2.70 eV

+

2.82 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 26 of 27

H /H 2 O 2 /H 2 O

-6.0

-7.5

Figure 6. Calculated VBM and CBM positions of g-C3N4 and g-CN along with their boron and phosphorous doped counter parts.

ACS Paragon Plus Environment

26

Page 27 of 27

TOC Absorp. coeff (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

200

g-CN B@g-CN P@g-CN

300

400

500

600

700

Wavelength (nm)

ACS Paragon Plus Environment

27