Relationship between Carbon Nitride Structure and Exciton Binding

Oct 16, 2015 - Graphitic (g)-CxNyHz has become a popular element in photocatalytic water splitting cells. Recently, the crystal structures of complete...
6 downloads 9 Views 3MB Size
Article pubs.acs.org/JPCC

Relationship between Carbon Nitride Structure and Exciton Binding Energies: A DFT Perspective Sigismund Melissen,* Tangui Le Bahers, Stephan N. Steinmann, and Philippe Sautet Université de Lyon, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, Centre Nationale de Recherche Scientifique, 46 allée d’Italie, Lyon F-69007 Cedex, France S Supporting Information *

ABSTRACT: Graphitic (g)-CxNyHz has become a popular element in photocatalytic water splitting cells. Recently, the crystal structures of completely polymerized g-C3N4 and incompletely polymerized g-C6N9H3 crystals based on the triazine and heptazine monomers have been characterized. In this manuscript we evaluate the atomic and electronic nature of these structures using DFT. The study revealed strongly corrugated structures for the fully polymerized g-C3N4 and planar structures for the incompletely polymerized g-C6N9H3. The exciton binding energies of the heptazine-based structures are lower than that of their triazine-based analogues and lower for the completely polymerized structures than their incompletely polymerized analogues. The rather low dielectric constant and charge mobilities result in high exciton binding energies and hence low dissociation probabilities for these excitons. This confirms the necessity of a morphology inspired by bulk heterojunction architectures to ensure efficient charge carrier generation. The studied compounds can be considered intermediates between typical inorganic and organic semiconductors in terms of their photoabsorption properties.

1. INTRODUCTION The sun is a virtually inexhaustible natural source of energy, the total solar power over the earth’s surface amounting to 65 Petawatts at any instant, with roughly half of this energy falling within the visible light range.1 One major drawback of solar energy is its diurnal variation, making its direct storage into chemical bonds an attractive alternative to photovoltaics. In the search for such energy solutions, inspiration can be drawn from photosynthesis. In this biological process, the chloroplasts inside a plant’s leaves use photons to generate excited electrons and holes that are able, after several transfers, to initiate chemical reactions. Artificial photosynthesis reproduces all these steps in a manmade device, allowing the use of sunlight to induce chemical reactions. Among the possible chemical reactions, the splitting of water is the most interesting one since it leads to the production of hydrogen.2,3 The two redox half reactions are provided 2H+(aq) + 2e− → H 2(g)

(1a)

2H 2O(l) → O2 (g) + 4H+(aq) + 4e−

(1b)

synthesis history of these macromolecules contains elements of controversy, due to several insufficiently substantiated claims that a graphitic C 3N4 structure was found, and with confusionas terms such as polymer, graphitic, and graphene-like structures are used interchangedly. The main reasons for this confusion are that the obtained macromolecules are often both amorphous and insoluble in most organic solvents, making them difficult to characterize.8 Two forms of g-CxNyHz are described in the literature: triazine-based (leading to gt-CxNyHz) and heptazine-based (forming gh-CxNyHz). Their tectons are displayed in Figure 1. Although a gt-C3N4 structure was recently characterized,8,9 no definitive experimental evidence for crystalline gh-C3N4 has been put forward.10 Moreover, it is commonly believed that

Figure 1. g-CxNyHz tectons. (i) Melamine (R = NH2), based on triazine (R = H). (ii) Melem (R = NH2) based on heptazine (R = H).

Graphitic carbon nitride, with general formula g-CxNyHz, has been extensively studied as a light-absorbing element in photocatalytic water splitting cells.4 The existence of graphitic Si3N4-analogous C(IV) nitrides (g-C3N4) was postulated5,6 by Liu et al. 25 years ago; however, the notion of macromolecular C3N4-based structures dates back to the 19th century.7 The © 2015 American Chemical Society

Received: July 21, 2015 Revised: October 16, 2015 Published: October 16, 2015 25188

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196

Article

The Journal of Physical Chemistry C

Figure 2. 2D chemical structures of completely (left) and incompletely (right) polymerized gt- and gh-CxNy studied in this manuscript. Individual tectons are highlighted in blue.

influences the exciton binding energy Eb, and the charge carrier effective masses m* (invariably expressed in terms of the free electron mass me) involved in their mobilities μ. In Table 1, desired ranges of values for the electronic properties above are reported.24

imperfections in the g-C xN yHz macrostructure do not necessarily hinder the photoactive process and possibly even constitute the catalytic sites and the ability of the fully polymerized photoactive material to perform the actual chemical conversions being generally contested.11 The overall photocatalytic performance of g-C3N4 is enhanced by adding cocatalysts through techniques such as nanoparticle decoration or by doping through elemental substitution.12 The general objective of nanoparticle decoration is the creation of catalytic sites for H2 and O2 evolution. For the first reaction, noble metals and their alloys are used,13,14 exploiting their possibilities to bond covalently to the NHx termini.15−17 For O2 evolution, transition metal oxides such as RuO2 are frequently used.18 Although the graphitic nature of the compound suggests a planar structure and some calculations are still performed on planar sheet geometries,10,19 several papers have stressed the importance of corrugation.20−23 In this study, four crystalline models of bulk g-CxNyHz will be studied using a theoretical approach. First, the fully polymerized and characterized g-C3N4 will be studied, starting from the results obtained earlier through molecular dynamics (MD) for the triazine-9 and heptazine-based20 variants. Then, the influence of the partial polymerization is investigated through g-C6N9H3 structures taken from the work by Schnick and co-workers. The proposed gh-C6N9H3 structure known as melon21 with P21/a symmetry and the fully characterized gtC6N9H3 structure22 from which LiCl was removed were considered as starting points. In Figure 2, the 2D chemical structures of the four systems are provided. Recently, the possibilities to calculate some of the fundamental properties governing a semiconductor’s ability to generate free charge carriers, marking the starting point for photocatalytic activity, using density functional theory (DFT) were pointed out24 and applied25,26 to several semiconductor families by some of us. These properties include the bandgap Eg (governing light absorption), the dielectric constant ϵr that

Table 1. Desired Ranges of Electronic Properties of Photoactive Semiconductors solar spectrum absorption

exciton dissociation

charge carrier diffusion

1.8 eV < Eg < 2.2 eV

Eb < 25 meV, ϵr > 10

m* < 0.5me

We have applied this approach9 to the g-CxNyHz structures and will detail our findings here as follows: the Computational Details section will outline the first-principles method employed here and briefly detail the semiconductor properties used to compare the studied structures. The comparison of experimentally obtained structures and our geometry optimizationsall of which are provided in the Supporting Informationis followed by a discussion of the influence of monomer nature and degree of polymerization on the material’s photoactivity. To the authors’ best knowledge, no systematic study into the ability of different g-CxNyHz structures to generate free charge carriers has been done, and the authors hope to draw some conclusions aiming to contribute to the discussion on structure−property relationships in g-CxNyHz-based photocatalytic cells.

2. COMPUTATIONAL DETAILS 2.1. DFT Calculations. All DFT calculations were performed employing the CRYSTAL14 suite, due to its natural treatment of systems with different types of periodicity and the efficiency of calculations performed with hybrid functionals using Gaussian type orbitals (GTOs).27,28 25189

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196

Article

The Journal of Physical Chemistry C

Figure 3. Optimized 3D-g-CxNyHz systems. (i) gt-C3N4, (ii) gh-C3N4*, (iii) gt-C6N9H3, and (iv) gh-C6N9H3. The unit cells are delineated in light gray. N atoms are in blue, C in gray, and H in white. dAr−Ar refers to the (average) interlayer distance. *In the Supporting Information, the geometry of gh-C3N4 is supplied in the axis system in which calculations were performed.

Gatti’s40 all-electron 6-31G(d,p) basis set41 was employed for H, C, and N. Convergence (10) but higher than a classical organic one (with values usually found between 2 and 4).71 Upon inspection of the band structures, the orbitals of the incompletely polymerized structures are more localized, as expected. Furthermore, strong band curvature is only obtained in the polymerization directions, with some interlayer mobility in the case of gt-C3N4. In the polymerization directions, the charge carrier effective masses are not strongly affected by the nature of the monomer, but a more complete polymerization allows for a better mobility of both the electrons and the holes. In the case of gt-C3N4, it must be remarked that due to orbital degeneracy a heavy and a light electron can be defined upon light absorption. Only the light electron is discussed here. For all systems, Eb is much greater than kBT, and Re is smaller than the tecton size, meaning that these excitons are localized and do not easily dissociate. The anisotropy of the different structures is expressed both in their dielectric constants, with ϵr,(100) > 3ϵr,(001), and effective masses, with m*(001) > 10m*(100), for both the hole and the electron. Thus, the exciton mobility is essentially limited to the two-dimensional molecular plane, accounting for the required factor 4 in eqs 4 and 5 (for the binding energies in 3D cf. Supporting Information).

determining the degree of polymerization. Overall, a good qualitative and quantitative agreement is obtained by comparing the geometries of the calculated structures to the experimentally characterized ones. Comment on the Phonon Band Spectrum of gt-C6N9H3. A measure for the interlayer interaction can be provided by phonon dispersion spectra. In particular, the splitting of the longitudinal and transverse acoustic (LA & TA) phonons starting at the Γ-point and determined at the A-point in the Brillouin Zone yields the required information. The validity of using an HSE06/D2 approach to compute such spectra can be seen by comparing experimentally and computationally obtained57,58 splitting values for hexagonal (h-) BN in the order of 100 cm−1. Graphite is known to exhibit a Γ-point splitting of 50 cm−1.59 Our calculation (cf. Table S1 and Figure S2 in the Supporting Information) of a splitting of ∼70 cm−1 confirms what can be hypothesized from both the geometry optimizations (with dAr−Ar equal to 3.14, 3.26, and 3.32 for hBN, gt-C6N9H3, and graphite, respectively) and the fact that in terms of electronegativity C lies between B and N. The relatively localized charges on h-BN allow for a stronger interlayer interaction than the delocalized electron cloud in the graphene layers in graphite, with gt-C6N9H3 lying between these two extremes. These spectra can be used further for a comprehension of properties such as the material’s thermal conductivity.60 3.2. Photoactivity Descriptors. The main numerical results of the photoactivity study are provided in Table 3, Table 3. Main Numerical Results of the Free Charge Carrier Generation Study gt-C3N4

gh-C3N4

gt-C6N9H3

gh-C6N9H3

gap nature

Γ→Γ

Γ→Γ

K→K

Γ→Γ

Eelg (eV) me*, mh* (me) ϵr {ϵ∞, ϵvib} Eb (meV) Eopt g (eV) Re (Å)

3.4 0.4, 1.4 6.8 {5.0, 1.8} 367 3.0 2.88

2.8 0.4, 1.4 7.2 {5.0, 2.3} 328 2.5 3.05

4.3 0.7, 1.9 4.6 {3.4, 1.2} 1350 3.0 1.16

3.4 0.9, 1.6 6.1 {4.1, 2.0} 840 2.6 1.38

while the band structures of the different studied systems are provided in Figure 5. To give the reader an idea of the influence of the choice of the PBE0 functional for the CPHF calculation,61 the band structures using the PBE0 functional were calculated and added to the Supporting Information, showing a similar band structure and a systematic increase (∼0.7 eV) in Eelg , as expected.24 A comparison between bandgap values obtained here and those found in the literature is provided in Figure 6. Note that in the reference study by Lotsch et al.21 for gh-C6N9H3 no photoabsorption study is carried out, and results by Niu et al.62 were used for comparison. All comparisons and particularly this last comparison demonstrate that the Wannier−Mott analysis presented in this paper acts as a first correction to determine Eopt g accurately. An advantage of this type of material is that for all structures the first electronic transition is direct. The first remarkable difference in their electronic properties is that the bandgap decreases markedly (∼0.6 eV) upon going from an incompletely to a completely polymerized system and upon going from a triazine- to a heptazine-based system. Upon

4. DISCUSSION The values for Eb obtained by Wei et al.19 using G0W0-BSE are between 1 and 2 eV. These values are slightly higher than the ones obtained by us using the Wannier−Mott model. The major reason for this discrepancy comes from the inclusion of 25192

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196

Article

The Journal of Physical Chemistry C

Figure 5. Band structures of different g-CxNyHz systems. The path corresponding to the [001] direction is highlighted in gray.

Figure 6. Comparison between the bandgap values obtained here and in the literature. Notes: §ref 9, ¶ref 10, ∥ref 62, †ref 63, ‡ref 64, *ref 65, and #ref 19.

Figure 7. g-CxNyHz valence and conduction band edge positions with respect to the water splitting half reactions for the optimized five-layer (5L) slabs. *NHE = Normal Hydrogen Electrode. Left-hand axis: electrochemical potential. Right-hand axis: absolute energy scale.

the vibrational term of the dielectric constant that reduces Eb by a factor ∼2. However, the qualitative conclusion remains similar: whatever model is used to calculate its energy, the excitons are not spontaneously dissociated at room temperature. This result confirms the similarities of these materials to classical semiconductor polymers that have exciton binding energies much higher than kBT. In the case of semiconductor polymers used in photovoltaics, the community developed the so-called “bulk heterojunction” architecture to allow the exciton

dissociation at the interface between the polymer and an electron-accepting material.72−74 For g-CxNyHz materials, the interface with the cocatalyst (Pt for instance) plays a similar role. The exciton is dissociated at the interface, and the electron is transferred to the co-catalyst in 25193

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196

Article

The Journal of Physical Chemistry C

Figure 8. Conceptualization of the free charge carrier generation in g-CxNyHz structures as proposed in this work.

the case of a hydrogen evolution reaction. As for the bulk heterojunction, the larger the interface, the better the exciton dissociation efficiency. This confirms why these materials have to be coated with highly dispersed nanoparticle cocatalysts,13,14,18 maximizing the uptake probability of the charge carriers before recombination. These ideas are visualized in Figure 8. The overall increases in charge carrier mobility and dielectric constant upon going from triazine- to heptazine-based systems and upon going from partially to completely polymerized systems are ultimately translated in the lowest Eb for gh-C3N4. Arguably, the use of incompletely polymerized g-C6N9H3 is still of intrinsic interest because of the ability of the NHx termini to bond to noble metal nanoparticles that perform the proton reduction reaction,15−17 but the experimental search for higher degrees of g-CxNyHz condensation in terms of tecton size and polymerization degree can be justified by the results above.

charge carrier uptake probability before recombination are vital for this type of material’s photocatalytic efficiency.



ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.5b07059. (1a) calibration of the optimization protocol; (1b) band structures calculated with PBE0; (2) crystal information: primary reference, number of atoms and Gibbs free energy per unit cell followed by structures in CIF-format; (3) exciton binding energies and radii within the framework of the 3D Wannier−Mott model; (4) further details on the phonon dispersion calculations (PDF)



AUTHOR INFORMATION

Corresponding Author

5. CONCLUSIONS A set of selected properties governing the capacity of a material to create free charge carriers upon illumination were computed for four graphitic (g)-CxNyHz structures by DFT using the range-separated hybrid functional HSE06. The effects of the nature of the monomer and degree of polymerization are studied. The incompletely polymerized g-C6N9H3 structures are planar and have relatively wide bandgaps. Excitons created in these structures have large binding energies; the advantage of this type of structure is the presence of primary and secondary amines that allow a cocatalyst such as Pd, Pt, or Au to bond to the material. The completely polymerized g-C3N4 structures show a marked corrugation of the molecular plane, and excitons generated in these materials have a lower binding energy, although still significantly higher than the thermal energy. The nonionicity and low polarizability of these materials is translated in a dielectric constant that is high compared to a typical polymeric semiconductor but low in comparison to the typical value for an inorganic semiconductor. This is subsequently translated in relatively localized excitons with a high binding energy, leading to a low dissociation probability at room temperature. It was found that heptazine-based materials and materials with higher degrees of polymerization should outperform triazine-based materials and those with a lower degree of polymerization. Bulk heterojunction-inspired interactions with dispersed nanoparticle cocatalysts maximizing

* E-ma il: [email protected]. Ph o n e : +33.4.7272.8640. Fax: +33.4.7272.8860. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Funding for this work was provided by the King Abdullah University of Science and Techonology (KAUST), within the framework of Special Academic Partnership Program “Water Splitting” (projects ENS 14.065 and KAUST 1974-02). The authors gratefully acknowledge Dr. A. Winfer and Dr. I. Shore (KAUST Supercomputing Laboratory), the computational resources provided by l’Institut du développement et des ressources en informatique scientifique (IDRIS, under project x2015080609) of the Centre Nationale de la Recherche Scientifique (CNRS), Prof. R. Orlando (University of Torino) for advice on the use of the CRYSTAL14 suite, and Prof. B. Lotsch (Max Planck Institute for Solid State Research), Dr. R. Kerber (Cambridge University), and Prof. K. Takanabe (KAUST) for valuable discussions on this topic.



ABBREVIATIONS CB, Conduction Band; CPHF, Coupled-Perturbed Hartree− Fock; DFT, Density Functional Theory; GTOs, Gaussian Type Orbitals; SCF, Self-Consistent Field; VB, Valence Band 25194

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196

Article

The Journal of Physical Chemistry C



(22) Wirnhier, E.; Döblinger, M.; Gunzelmann, D.; Senker, J.; Lotsch, B. V.; Schnick, W. Poly(Triazine Imide) with Intercalation of Lithium and Chloride Ions [(C3N3)2(NH(x)Li(1‑x))LiCl]: a Crystalline 2D Carbon Nitride Network. Chem. - Eur. J. 2011, 17, 3213−21. (23) Sehnert, J.; Baerwinkel, K.; Senker, J. Ab Initio Calculation of Solid-State NMR Spectra for Different Triazine and Heptazine Based Structure Proposals of g-C3N4. J. Phys. Chem. B 2007, 111, 10671− 10680. (24) Le Bahers, T.; Rérat, M.; Sautet, P. Semiconductors Used in Photovoltaic and Photocatalytic Devices: Assessing Fundamental Properties from DFT. J. Phys. Chem. C 2014, 118, 5997−6008. (25) Melissen, S.; Labat, F.; Sautet, P.; Le Bahers, T. Electronic Properties of PbX3CH3NH3 (X = Cl, Br, I) Compounds for Photovoltaic and Photocatalytic Applications. Phys. Chem. Chem. Phys. 2015, 17, 2199−2209. (26) Kandiel, T. A.; Anjum, D. H.; Sautet, P.; Le Bahers, T.; Takanabe, K. Electronic Structure and Photocatalytic Activity of Wurtzite CuGaAs Nanocrystals and Their Zn Substitution. J. Mater. Chem. A 2015, 3, 8896−8904. (27) Dovesi, R.; Orlando, R.; Erba, A.; Zicovich-Wilson, C. M.; Civalleri, B.; Casassa, S.; Maschio, L.; Ferrabone, M.; De La Pierre, M.; D’Arco, P.; et al. CRYSTAL14: A Program for the Ab Initio Investigation of Crystalline Solids. Int. J. Quantum Chem. 2014, 114, 1287−1317. (28) Dovesi, R.; Saunders, V. R.; Roetti, C.; Orlando, R.; ZicovichWilson, C. M.; Pascale, B.; Civalleri, B.; Doll, K.; Harrison, N. M.; Bush, I. J. et al. CRYSTAL14 User’s Manual; University of Torino: Italy, 2014. (29) Ferrero, M.; Rérat, M.; Orlando, R.; Dovesi, R. Coupled Perturbed Hartree-Fock for Periodic Systems: The Role of Symmetry and Related Computational Aspects. J. Chem. Phys. 2008, 128, 014110. (30) Ferrero, M.; Rérat, M.; Orlando, R.; Dovesi, R. The Calculation of Static Polarizabilities of 1−3D Periodic Compounds. The Implementation in the CRYSTAL Code. J. Comput. Chem. 2007, 29, 1450−1459. (31) Ferrero, M.; Rérat, M.; Kirtman, B.; Dovesi, R. The Calculation of Static Polarizabilities of 1−3D Periodic Compounds. The Implementation in the CRYSTAL Code. J. Chem. Phys. 2008, 129, 244110. (32) Noel, Y.; Zicovich-Wilson, C. M.; Civalleri, B.; D’Arco, P.; Dovesi, R. Polarization Properties of ZnO and BeO: An Ab Initio Study Through the Berry Phase and Wannier Functions Approaches. Phys. Rev. B: Condens. Matter Mater. Phys. 2001, 65, 14111. (33) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207−8215. (34) 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. (35) Krukau, A.; Vydrov, O.; Izmaylov, A.; Scuseria, G. Influence of the Exchange Screening Parameter on the Performance of Screened Hybrid Functionals. J. Chem. Phys. 2006, 125, 224106. (36) Barone, V.; Hod, O.; Peralta, J. E.; Scuseria, G. E. Accurate Prediction of the Electronic Properties of Low-Dimensional Graphene Derivatives Using a Screened Hybrid Density Functional. Acc. Chem. Res. 2011, 44, 269−279. (37) Henderson, T. M.; Paier, J.; Scuseria, G. E. Accurate Treatment of Solids with the HSE Screened Hybrid. Phys. Status Solidi B 2011, 248, 767−774. (38) Grimme, S. Semiempirical GGA-type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787−1799. (39) Wang, V.; Liu, R. J.; He, H. P.; Yang, C. M.; Ma, L. Hybrid Functional with Semi-Empirical Van der Waals Study of Native Defects in Hexagonal BN. Solid State Commun. 2014, 177, 74−79. (40) Gatti, C.; Saunders, V. R.; Roetti, C. Crystal Field Effects on the Topological Properties of the Electron Density in Molecular Crystals: The Case of Urea. J. Chem. Phys. 1994, 101, 10686−10696.

REFERENCES

(1) Lewis, N.; Crabtree, G. Basic Research Needs for Solar Energy Utilization; Office of Science US Dept of Energy, 2005. (2) McEvoy, J.; Brudvig, G. Water-Splitting Chemistry of Photosystem II. Chem. Rev. 2006, 106, 4455−4483. (3) Pinaud, B. A.; Benck, J. D.; Seitz, L. C.; Forman, A. J.; Chen, Z.; Deutsch, T. G.; James, B. D.; Baum, K. N.; Baum, G. N.; Ardo, S.; et al. Technical and Economic Feasibility of Centralized Facilities for Solar Hydrogen Production via Photocatalysis and Photoelectrochemistry. Energy Environ. Sci. 2013, 6, 1983−2002. (4) 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. 2009, 8, 76−80. (5) Liu, A. Y.; Cohen, M. L. Prediction of New Low Compressibility Solids. Science 1989, 245, 841−842. (6) Liu, A. Y.; Cohen, M. L. Structural Properties and Electronic Structure of Low-Compressibility Materials: β-Si3N4 and hypothetical β-C3N4. Phys. Rev. B: Condens. Matter Mater. Phys. 1990, 41, 10727− 10734. (7) Liebig, J. About Some Nitrogen Compounds. Ann. Pharm. 1834, 10, 10. (8) Kroke, E. gt-C3N4 - The First Stable Binary Carbon(IV) Nitride. Angew. Chem., Int. Ed. 2014, 53, 11134−11136. (9) Algara-Siller, G.; Severin, N.; Chong, S. Y.; Björkman, T.; Palgrave, R. G.; Laybourn, A.; Antonietti, M.; Khimyak, Y. Z.; Krasheninnikov, A. V.; Rabe, J. P.; et al. Triazine-Based Graphitic Carbon Nitride: a Two-Dimensional Semiconductor. Angew. Chem. 2014, 126, 7580−7585. (10) Zuluaga, S.; Liu, L.-H.; Shafiq, N.; Rupich, S. M.; Veyan, J.-F.; Chabal, Y. J.; Thonhauser, T. Structural Band-Gap Tuning in g-C3N4. Phys. Chem. Chem. Phys. 2015, 17, 957−962. (11) Gillan, E. G. Synthesis of Nitrogen-Rich Carbon Nitride Networks from an Energetic Molecular Azide Precursor. Chem. Mater. 2000, 12, 3906−3912. (12) Cao, S.; Yu, J. g-C3N4-Based Photocatalysts for Hydrogen Generation. J. Phys. Chem. Lett. 2014, 5, 2101−2107. (13) Wang, Y.; Wang, X.; Antonietti, M. Polymeric Graphitic Carbon Nitride as a Heterogeneous Organocatalyst: from Photochemistry to Multipurpose Catalysis to Sustainable Chemistry. Angew. Chem., Int. Ed. 2012, 51, 68−89. (14) Di, Y.; Wang, X.; Thomas, A.; Antonietti, M. Making MetalCarbon Nitride Heterojunctions for Improved Photocatalytic Hydrogen Evolution with Visible Light. ChemCatChem 2010, 2, 834−838. (15) He, D.; Zeng, C.; Xu, C.; Cheng, N.; Li, H.; Mu, S.; Pan, M. Polyaniline-Functionalized Carbon Nanotube Supported Platinum Catalysts. Langmuir 2011, 27, 5582−5588. (16) Cheng, H.; Chen, L.; Cooper, A. C.; Sha, X.; Pez, G. P. Hydrogen Spillover in the Context of Hydrogen Storage Using SolidState Materials. Energy Environ. Sci. 2008, 1, 338−345. (17) Tao, S.; Yu, L.-J.; Pang, R.; Huang, Y.-F.; Wu, D.-Y.; Tian, Z.-Q. Binding Interaction and Raman Spectra of π-Conjugated Molecules Containing CH2/NH2 Groups Adsorbed on Silver Surfaces: A DFT Study of Wagging Modes. J. Phys. Chem. C 2013, 117, 18891−18903. (18) Maeda, K.; Wang, X.; Nishihara, Y.; Lu, D.; Antonietti, M.; Domen, K. Photocatalytic Activities of Graphitic Carbon Nitride Powder for Water Reduction and Oxidation under Visible Light. J. Phys. Chem. C 2009, 113, 4940−4947. (19) Wei, W.; Jacob, T. Strong Excitonic Effects in the Optical Properties of Graphitic Carbon Nitride g-C3N4. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 085202. (20) Gracia, J.; Kroll, P. Corrugated Layered Heptazine-Based Carbon Nitride: the Lowest Energy Modifications of C3N4 Ground State. J. Mater. Chem. 2009, 19, 3013−3019. (21) Lotsch, B. V.; Döblinger, M.; Sehnert, J.; Seyfarth, L.; Senker, J.; Oeckler, O.; Schnick, W. Unmasking Melon by a Complementary Approach Employing Electron Diffraction, Solid-State NMR Spectroscopy, and Theoretical Calculations-Structural Characterization of a Carbon Nitride Polymer. Chem. - Eur. J. 2007, 13, 4969−80. 25195

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196

Article

The Journal of Physical Chemistry C

Graphene Oxide Composite: The Role of Oxygen. Chem. Mater. 2015, 27, 1612−1621. (64) Meek, G. A.; Baczewski, A. D.; Little, D. J.; Levine, B. G. Polaronic Relaxation by Three-Electron Bond Formation in Graphitic Carbon Nitrides. J. Phys. Chem. C 2014, 118, 4023−4032. (65) Xu, Y.; Gao, S.-P. Band Gap of C3N4 in the GW Approximation. Int. J. Hydrogen Energy 2012, 37, 11072−11080. (66) Komsa, H. P.; Krasheninnikov, A. V. Effects of Confinement and Environment on the Electronic Structure and Exciton Binding Energy of MoS2 from First Principles. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 1−6. (67) 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. (68) Butchosa, C.; Guiglion, P.; Zwijnenburg, M. A. Carbon Nitride Photocatalysts for Water Splitting: A Computational Perspective. J. Phys. Chem. C 2014, 118, 24833−24842. (69) Guiglion, P.; Butchosa, C.; Zwijnenburg, M. A. Polymeric Watersplitting Photocatalysts; a Computational Perspective on the Water Oxidation Conundrum. J. Mater. Chem. A 2014, 2, 11996. (70) Trasatti, S. The Absolute Electrode Potential: an Explanatory Note. Pure Appl. Chem. 1986, 58, 955−966. (71) Blom, P. W. M.; Mihailetchi, V. D.; Koster, L. J. A.; Markov, D. E. Device Physics of Polymer: Fullerene Bulk Heterojunction Solar Cells. Adv. Mater. 2007, 19, 1551−1566. (72) Yu, G.; Gao, J.; Hummelen, J. C.; Wudl, F.; Heeger, A. Polymer Photovoltaic Cells - Enhanced Efficiencies via a Network of Internal Donor-Acceptor Heterojunctions. Science 1995, 270, 1789−1791. (73) Gunes, S.; Neugebauer, H.; Sariciftci, N. Conjugated PolymerBased Organic Solar Cells. Chem. Rev. 2007, 107, 1324−1338. (74) Coakley, K. M.; McGehee, M. D. Conjugated Polymer Photovoltaic Cells. Chem. Mater. 2004, 16, 4533−4542.

(41) Ditchfield, R.; Hehre, W. J.; Pople, J. A. Self-Consistent Molecular-Orbital Methods. IX. An Extended Gaussian-Type Basis for Molecular-Orbital Studies of Organic Molecules. J. Chem. Phys. 1971, 54, 724−728. (42) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188−5192. (43) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods Without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158−6170. (44) Gilat, G. Analysis of Methods for Calculating Spectral Properties in Solids. J. Comput. Phys. 1972, 10, 432−465. (45) Gilat, G.; Raubenheimer, L. J. Accurate Numerical Method for Calculating Frequency-Distribution Functions in Solids. Phys. Rev. 1966, 144, 390−395. (46) Gilat, G.; Raubenheimer, L. J. Erratum: Accurate Numerical Method for Calculating Frequency-Distribution Functions in Solids. Phys. Rev. 1966, 147, 670. (47) Pascale, F.; Zicovich-Wilson, C. M.; López Gejo, F.; Civalleri, B.; Orlando, R.; Dovesi, R. The Calculation of the Vibrational Frequencies of Crystalline Compounds and its Implementation in the CRYSTAL Code. J. Comput. Chem. 2004, 25, 888−897. (48) Zicovich-Wilson, C. M.; Pascale, F.; Roetti, C.; Saunders, V. R.; Orlando, R.; Dovesi, R. Calculation of the Vibration Frequencies of αQuartz: The Effect of Hamiltonian and Basis Set. J. Comput. Chem. 2004, 25, 1873−1881. (49) Panda, A.; Renshaw, C. K.; Oskooi, A.; Lee, K.; Forrest, S. R. Excited State and Charge Dynamics of Hybrid Organic/Inorganic Heterojunctions. II. Experiment. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 045303. (50) He, X.-F. Excitons in Anisotropic Solids: The Model of Fractional-Dimensional Space. Phys. Rev. B: Condens. Matter Mater. Phys. 1991, 43, 2063−2069. (51) Zhu, X.-Y.; Yang, Q.; Muntwiler, M. Charge-Transfer Excitons at Organic Semiconductor Surfaces and Interfaces. Acc. Chem. Res. 2009, 42, 1779−1787. (52) Gregg, B. A. The Photoconversion Mechanism of Excitonic Solar Cells. MRS Bull. 2005, 30, 20−22. (53) Björkman, T. Van der Waals Density Functional for Solids. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 165109. (54) Mayo, S. L.; Olafson, B. D.; Goddard, W. A. DREIDING: a Generic Force Field for Molecular Simulations. J. Phys. Chem. 1990, 94, 8897−8909. (55) Stokes, H. T.; Hatch, D. M. FINDSYM: Program for Identifying the Space-Group Symmetry of a Crystal. J. Appl. Crystallogr. 2005, 38, 237−238. (56) Burke, K.; Perdew, J. P.; Wang, Y. In Electronic Density Functional Theory; Dobson, J. F., Vignale, G., Das, M. P., Eds.; Springer US: New York, 1998; DOI: 10.1007/978-1-4899-0316-7_7, pp 81− 111. (57) Dai, J.; Wu, X.; Yang, J.; Zeng, X. C. Unusual Metallic Microporous Boron Nitride Networks. J. Phys. Chem. Lett. 2013, 4, 3484−3488. (58) Geick, R.; Perry, C. H.; Rupprecht, G. Normal Modes in Hexagonal Boron Nitride. Phys. Rev. 1966, 146, 543−547. (59) Nicklow, R.; Wakabayashi, N.; Smith, H. G. Lattice Dynamics of Pyrolytic Graphite. Phys. Rev. B 1972, 5, 4951−4962. (60) Wu, Y. H.; Yu, T.; Shen, Z. X. Two-Dimensional Carbon Nanostructures: Fundamental Properties, Synthesis, Characterization, and Potential Applications. J. Appl. Phys. 2010, 108, 071301. (61) Czaja, W.; Gränacher, I. A Relation Between Gap Energy and Dielectric Constant in Semiconductors of Diamond Structure. Helv. Phys. Acta 1963, 36, 1073−1080. (62) Niu, P.; Yin, L.-C.; Yang, Y.-Q.; Liu, G.; Cheng, H.-M. Increasing the Visible Light Absorption of Graphitic Carbon Nitride (Melon) Photocatalysts by Homogeneous Self-Modification with Nitrogen Vacancies. Adv. Mater. 2014, 26, 8046−8052. (63) Xu, L.; Huang, W.-Q.; Wang, L.-L.; Tian, Z.-A.; Hu, W.; Ma, Y.; Wang, X.; Pan, A.; Huang, G.-F. Insights into Enhanced Visible-Light Photocatalytic Hydrogen Evolution of g-C3N4 and Highly Reduced 25196

DOI: 10.1021/acs.jpcc.5b07059 J. Phys. Chem. C 2015, 119, 25188−25196