School of Applied Science, Taiyuan University of Science and Technology, Taiyuan 030024, China
2.
Department of Chemistry and Chemical Engineering, Jinzhong University, Jinzhong, Shanxi 030619, China
3.
Shanxi Key Laboratory of Metal Forming Theory and Technology, School of Material Science and Engineering, Taiyuan University of Science and Technology, Taiyuan 030024, China
Received Date:
07 January 2021 Revised Date:
17 June 2021 Available Online:
10 September 2021
Abstract:
The first-principles method based on density functional theory was performed to investigate the electronic structure and optical properties of the armchair nanoribbons (AC-g-C3N4NRs) and zigzag g-C3N4 nanoribbons (ZZ-g-C3N4Rs). The results show that the edge H atoms of AC-g-C3N4NRs and ZZ-g-C3N4NRs can exist stably. The valence band maximums (VBMs) of AC-g-C3N4NRs are mainly contributed by most of N atoms, whereas the VBMs of ZZ-g-C3N4NRs are contributed by the N atoms near the CH edge. The conduction band minimums (CBMs) of AC-g-C3N4Rs mainly belong to C and N atoms near the one edge or two edges of AC-g-C3N4NRs, while the CBMs of ZZ-g-C3N4Rs mainly belong to C and N atoms near the NH edge of ZZ-g-C3N4NRs. The absorption coefficient and the reflectivity of AC-g-C3N4NR or ZZ-g-C3N4NR increased with the increasing width of the corresponding nanoribbon. An obvious blueshift phenomenon of the absorption coefficient could be generated in the low-energy range as the width of AC-g-C3N4NR increased.
As a low-cost, high-stable, nontoxic and visiblelight-responsive photocatalyst, graphitic carbon nitride (g-C3N4) has a multifunctional application in photocatalytic hydrogen evolution[1-3], CO2 reduction[4-5] and photo-catalytic degradation of pollutants[6-7]. In order to further enhance the photocatalytic performance of g-C3N4, bulk g-C3N4 is usually cut into nano-sized g-C3N4, such as g-C3N4 nanosheets[3, 8-11] and g-C3N4 nanoribbons (g-C3N4 NRs)[12-13], which possesses favorable photocatalytic activity because of the larger specific surface area with abundant active sites and short diffusion distance of photogenerated charge carriers.
Due to the quantum size effect, the properties of g-C3N4 NRs are considerably different from those of g-C3N4 nanosheets. In experiments, Wang et al. have successfully synthesized the Mn-doped g-C3N4 NR catalyst by a two-step calcination method[13]. Zhao et al. have fabricated g-C3N4 NR on graphene sheets by using a simple one-step hydrothermal method[12]. However, litttle attention is paid to theoretical research on the armchair nanoribbons (AC-g-C3N4NRs) or zigzag g-C3N4 nanoribbons (ZZ-g-C3N4NRs) which can be obtained by cutting nanosheets along specific directions, and thus the electronic and optical properties of these above two g-C3N4 NRs as function of widths are not clear. Exposure of these properties will help the design and fabrication of g-C3N4 NRs-based electronic and optical devices in experiment.
In this work, the rest of the paper is organized as follows: the computational method and models of AC-g-C3N4NRs and ZZ-g-C3N4NRs with the width of 7 (i.e. the number of the atom chains) are given in Section 2. The binding energies, band gaps, band structure, partial charge density, partial density of states and optical properties of g-C3N4 NRs are analyzed and the conclusions are drawn in the last section.
1.
Computational method and models
The computation of electronic and optical properties was performed within the framework of the density functional theory (DFT) implemented in the Vienna ab-initio Simulation Package (VASP)[14-18]. The electronionic core interactions were treated by the projected augmented wave (PAW) potentials[19]. The Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional[20] within the generalized-gradient approximation (GGA) was used in order to yield the correct ground-state structure of the systems. An efficient Broyden/Pulay mixing scheme[21-22] was used for the mixing of the charge density. The cut-off energy of the plane waves was 500 eV. The energies and the forces on each ion were converged to less than 10-5 eV·atom-1 and 0.1 eV·nm-1, respectively. The Gaussian smearing broadening was chosen as 0.05 eV. The Brillouin zones were sampled by 1×1×11 and 1×1×25 K-point meshes according to Monkhorst-Pack scheme[23] for calculations of electronic structures and optical properties, respectively. The number of energy bands was increased to 300 when the optical properties were computed. The absorption coefficient and reflectivity can be derived from the computational dielectric function. Compared with the experiment value 2.7 eV, the computed band gap of bulk g-C3N4 was 1.28 eV and the difference in band gap was 1.42 eV. The same conclusion is also drawn by Wu et al[24]. However, the PBE computations can still reveal the variation tendency of the electronic structures and optical properties.
The width of AC-g-C3N4NR (or ZZ-g-C3N4NRs) was classified by the number of the atom chains NA (or NZ) across the ribbon width and denoted as NA-AC-g-C3N4NR (or NZ-ZZ-g-C3N4NR). In the present work, we focus our attention on the structures of AC-g-C3N4NR and ZZ-g-C3N4NRs with the width NA or NZ=4~10. As examples, geometry structures of 7-AC-g-C3N4NR and 7-ZZ-g-C3N4NR are shown in Fig. 1a and 1b, respectively. g-C3N4NRs periodically extend along the z direction. To avoid the interactions between the neighboring ribbons, g-C3N4NRs were separated from each other by the vacuum region with 1.6 nm in both edge-to-edge (i.e. x direction) and layer-to-layer (i.e. y direction).
Figure 1
Figure 1.
Geometry structures of monolayer g-C3N4NRs with H termination: (a) 7-AC-g-C3N4NR; (b) 7-ZZ-g-C3N4NR
A periodic unit cell is presented between the two red solid lines; Green straight line or zigzag curve shows atom chain; Total number of the atom chains across AC-g-C3N4NRs and ZZ-g-C3N4NRs width is also described by NA and NZ, respectively; Gray, blue and white balls represent the C, N and H atoms, respectively
The binding energies Eb are usually adopted to investigate the binding strength between H atoms and g-C3N4NRs. Eb was computed according to Eb=(Epassivated-Ebare
-nEH)/n, where n is the number of H atoms, Epassivated and Ebare refer to the total energy of g-C3N4NR with H termination and g-C3N4NR without H termination, respectively. EH is the half of the total energy of a H2 molecule at 0 K and 0 GPa and was computed by equation EH
=1/2EH2. The evolution of the binding energies of AC -and ZZ-g-C3N4NRs as a function of the width NA or NZ is shown in Fig. 2a. It can be seen that the binding energies of these fourteen g-C3N4NRs are all negative values, indicating that the edge H atoms can exist stably. The black curve and red curve represent that the binding energy of the AC-g-C3N4NR has a small fluctuation, whereas that of ZZ-g-C3N4NR has a large fluctuation with the increasing number of atom chains. In addition, the smaller binding energies of ZZ-g-C3N4NRs with the width NZ
=4, 6, 8, 10 indicate that H atoms are preferred to adsorb on the edges of these systems. Fig. 2b displays the evolution of the band gaps of AC-and ZZ-g -C3N4NRs as a function of the width NA or NZ. It can be seen from Fig. 2b that the band gaps of the remaining eleven semiconductor g-C3N4NRs present a small fluctuation except for ZZ-g-C3N4NR with NZ=5, 7, 9 which are metallic systems with no band gaps. Table 1 shows the molar ratio (ε) of C and N atoms of the AC-and ZZ-g-C3N4NRs. Table 1 and Fig. 2 indicate that the binding energies and band gaps of AC-g-C3N4NRs have no obvious change rule, while ZZ-g-C3N4NRs have the smaller binding energies and the larger band gaps when the ε of C and N atoms is 3/4.
Figure 2
Figure 2.
Evolution of the binding energies (a) and band gaps (b) of AC-g-C3N4NRs and ZZ-g-C3N4NRs as a function of the width NA or NZ
In the following sections, the aforementioned stable semiconductor AC-g-C3N4NRs with the width NA
=4, 7, 10 and semiconductor ZZ-g-C3N4NRs with the width NZ=4, 10 will be taken as examples for further investigating the electronic and optical properties of g-C3N4NRs.
2.2
Band structure, partial charge density and partial density of states
Fig. 3 shows band structures of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10. Fig. 3a and 3d reflect that AC-g-C3N4NR and ZZ-g-C3N4NR with NA or NZ=4 are both indirect band gap semiconductors because of the valence band maximum (VBM) and the conduction band minimum (CBM) locating at different points of K-space, whereas Fig. 3b, 3c and 3e are all the direct band gap systems due to the VBM and CBM located at Γ point of K-space.
Figure 3
Figure 3.
Band structures of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10
Red and green bands denote the top valence band and the bottom conduction band, respectively; Blue circles shown on the red and green bands indicate the VBM and CBM, respectively; Fermi level (EF) was set to 0 eV, and it is the energy level at which the Fermi-Dirac distribution function of an assembly of fermions is equal to one-half
In order to intuitively recognize compositions of the VBMs and CBMs, Fig. 4 shows partial charge densities for VBMs and CBMs of AC-g-C3N4NRs with (a) NA= 4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ= 4, (e) NZ=10. It can be seen that the VBMs of AC-g-C3N4NRs (the red isosurface in Fig. 4a~4c are mainly distributed on most of N atoms, whereas the VBMs of ZZ-g-C3N4NRs (the red isosurface in Fig. 4d and 4e are contributed by the N atoms near the CH edge. Moreover, these charge distributions contributing to VBMs are all within the g-C3N4 plane. However, the CBMs of AC-g-C3N4NRs (the green isosurface in Fig. 4a~4c) mainly belong to C and N atoms near the one edge or two edges of AC-g-C3N4NRs, while CBMs of ZZ-g-C3N4NRs (the green isosurface in Fig. 4d and 4e) mainly belong to C and N atoms near the NH edge of ZZ-g-C3N4NRs. Interestingly, these charge distributions are perpendicular to the g-C3N4 plane.
Figure 4
Figure 4.
Partial charge densities for VBMs (red isosurface) and CBMs (green isosurface) of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10
Furthermore, we computed the partial density of states of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA
=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10. As shown in Fig. 5, the VBMs of AC-and ZZ-g-C3N4NRs are composed of the 2p states of N atoms, while the CBMs are mainly contributed by the 2p states of C and N atoms of the systems. In terms of element compositions, this is in good agreement with the above analysis of the partial charge densities.
Figure 5
Figure 5.
Partial density of states of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10
In the present work, for comparison, the optical properties of g-C3N4NRs, bulk g-C3N4 and the armchair graphene nanoribbons (AC-GNRs) are characterized by absorption coefficient (α), reflectivity (R) and energy loss (L), which are defined as follows[25-26]:
Where the real part ε1 and imaginary part ε2 of the complex dielectric function ε = ε1+jε2 can be obtained based on the computational electronic states, and ω is the circular frequency.
Fig. 6 shows the computed optical absorption coefficients of (a) AC-g-C3N4NRs with NA=4, 7, 10, (b) ZZ-g -C3N4NRs with NZ=4, 10, (c) AC-GNRs with NA=4, 7, 10 and (d) bulk g-C3N4 as a function of energy. In the work, ZZ-GNRs are not investigated because the computed band gaps are 0. Seen from the Fig. 6a and 6b, the absorption coefficient of AC-g-C3N4NR or ZZ-g-C3N4NR increased with the increasing width of the corresponding nanoribbon. In Fig. 6a, these absorption peaks of 4-AC-g-C3N4NR, 7-AC-g-C3N4NR and 10-AC-g-C3N4NR (i. e. the black, red and green absorption peaks) in the low-energy range of 0~5 eV are specifically located at 3.58, 3.84, and 4.31 eV, respectively. Therefore, an obvious blueshift phenomenon can be generated in the low-energy range as the width increase of AC-g-C3N4NR. In Fig. 6b, the black and green absorption peaks are respectively located at 4.28 and 4.32 eV in the low energy region, indicating that there is only a very weak blueshift. By contrast, Fig. 6c displays an obvious redshift phenomenon in the low-energy range as the width increase of AC-GNR. Compared with bulk g-C3N4 (shown in Fig. 6d), the AC-and ZZ-g-C3N4NRs have the smaller absorption coefficient.
Figure 6
Figure 6.
Computed absorption coefficients of (a) AC-g-C3N4NRs with NA=4, 7, 10, (b) ZZ-g-C3N4NRs with NZ=4, 10,
(c) AC-GNRs with NA=4, 7, 10 and (d) bulk g-C3N4
Fig. 7 shows the computational reflectivity of (a) AC-g-C3N4NRs with NA=4, 7, 10 and (b) ZZ-g-C3N4NRs with NZ=4, 10, (c) AC-GNRs with NA=4, 7, 10 and (d) bulk g-C3N4 as a function of energy. It also can be seen from Fig. 7a and 7b that the reflectivity is increasing with the increased width increase of AC-or ZZ-g-C3N4NR, and the strongest reflectivities of these NRs are all in the low energy range of 2~5 eV. By contrast, AC-GNRs can draw the same conclusion (Fig. 7c). In addition, bulk g-C3N4 has a larger reflectivity compared with g-C3N4NRs.
Figure 7
Figure 7.
Computed reflectivity of (a) AC-g-C3N4NRs with NA=4, 7, 10, (b) ZZ-g-C3N4NRs with NZ=4, 10,
(c) AC-GNRs with NA=4, 7, 10 and (d) bulk g-C3N4
The electronic and optical properties of g-C3N4 nanoribbons have been investigated by using the first-principles calculations. The results are as follows:
(1) The edge H atoms of AC-and ZZ-g-C3N4NRs studied in the present work can exist stably. AC-g-C3N4NRs with NA=4~10 and ZZ-g-C3N4NRs with NZ=4, 6, 8, 10 are semiconductors, whereas ZZ-g-C3N4NR with NZ=5, 7, 9 are metallic systems with no band gaps.
(2) The VBMs of AC-g-C3N4NRs are mainly distributed on most of N atoms, whereas the VBMs of ZZ-g -C3N4NRs are contributed by the N atoms near the CH edge. The CBMs of AC-g-C3N4NRs mainly belong to C and N atoms near the one edge or two edges of AC-g-C3N4NRs, while the CBMs of ZZ-g-C3N4NRs mainly belong to C and N atoms near the NH edge of ZZ-g-C3N4NRs. The VBMs of AC-and ZZ-g-C3N4NRs are composed of the 2p states of N atoms, while the CBMs are mainly contributed by the 2p states of C and N atoms of the systems.
(3) The absorption coefficient and the reflectivity of AC-g-C3N4NR or ZZ-g-C3N4NR are increased with the width increase of the corresponding nanoribbon. For the absorption coefficient, an obvious blueshift phenomenon can be generated in the low-energy range as the width increase of AC-g-C3N4NR.
Acknowledgements:
The authors would like to acknowledge the support from the National Natural Science Foundation of China (Grants No.11705124, 51871158, 11704274).
[1]
Li Y R, Wang Z W, Xia T, Ju H X, Zhang K, Long R, Xu Q, Wang C M, Song L, Zhu J F, Jiang J, Xiong Y J. Adv. Mater., 2016, 28(32): 6959-6965 doi: 10.1002/adma.201601960
[2]
Cao S W, Jiang J, Zhu B C, Yu J G. Phys. Chem. Chem. Phys., 2016, 18(28): 19457-19463 doi: 10.1039/C6CP02832B
[3]
Sun S D, Li J, Cui J, Gou X F, Yang Q, Liang S H, Yang Z M, Zhang J M. Inorg. Chem. Front., 2018, 5(7): 1721-1727 doi: 10.1039/C8QI00242H
[4]
Xia P F, Zhu B C, Yu J G, Cao S W, Jaroniec M. J. Mater. Chem. A, 2017, 5(7): 3230-3238 doi: 10.1039/C6TA08310B
[5]
Gao G P, Jiao Y, Waclawik E R, Du A J. J. Am. Chem. Soc., 2016, 138(19): 6292-6297 doi: 10.1021/jacs.6b02692
Ong W J, Tan L L, Ng Y H, Yong S T, Chai S P. Chem. Rev., 2016, 116: 7159-7329 doi: 10.1021/acs.chemrev.6b00075
[8]
Teng Z Y, Lv H Y, Wang C Y, Xue H G, Peng H, Wang G X. Carbon, 2017, 113: 63-75 doi: 10.1016/j.carbon.2016.11.030
[9]
Wang Y L, Tian Y, Yan L K, Su Z M. J. Phys. Chem. C, 2018, 122(14): 7712-7719 doi: 10.1021/acs.jpcc.8b00098
[10]
Wu H P, Liu Y Z, Kan E J, Ma Y M, Xu W J, Li J, Yan M C, Lu R F, Wei J F, Qian Y. J. Phys. D: Appl. Phys., 2016, 49(29): 295301 doi: 10.1088/0022-3727/49/29/295301
[11]
Li H S, Hu H Q, Bao C J, Hua J, Zhou H C, Liu X B, Liu X D, Zhao M W. Phys. Chem. Chem. Phys., 2015, 17(8): 6028-6035 doi: 10.1039/C4CP05560H
[12]
Zhao Y, Zhao F, Wang X P, Xu C Y, Zhang Z P, Shi G Q, Qu L T. Angew. Chem. Int. Ed., 2014, 53(50): 13934-13939 doi: 10.1002/anie.201409080
[13]
Wang J C, Cui C X, Kong Q Q, Ren C Y, Li Z, Qu L, Zhang Y, Jiang K. ACS Sustainable Chem. Eng., 2018, 6(7): 8754-8761 doi: 10.1021/acssuschemeng.8b01093
Wu H Z, Zhong Q H, Bandaru S, Liu J, Lau W M, Li L L, Wang Z L. J. Phys. : Condens. Matter, 2018, 30: 155303
[25]
Zhu B C, Zhang J F, Jiang C J, Cheng B, Yu J G. Appl. Catal. B, 2017, 207: 27-34
[26]
Yuan Y, Huang Y H, Ma F, Zhang Z Q, Wei X M. J. Mater. Sci., 2017, 52(14): 8546-8555
Figure 1
Geometry structures of monolayer g-C3N4NRs with H termination: (a) 7-AC-g-C3N4NR; (b) 7-ZZ-g-C3N4NR
A periodic unit cell is presented between the two red solid lines; Green straight line or zigzag curve shows atom chain; Total number of the atom chains across AC-g-C3N4NRs and ZZ-g-C3N4NRs width is also described by NA and NZ, respectively; Gray, blue and white balls represent the C, N and H atoms, respectively
Figure 3
Band structures of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10
Red and green bands denote the top valence band and the bottom conduction band, respectively; Blue circles shown on the red and green bands indicate the VBM and CBM, respectively; Fermi level (EF) was set to 0 eV, and it is the energy level at which the Fermi-Dirac distribution function of an assembly of fermions is equal to one-half
Figure 4
Partial charge densities for VBMs (red isosurface) and CBMs (green isosurface) of AC-g-C3N4NRs with (a) NA=4, (b) NA=7, (c) NA=10 and ZZ-g-C3N4NRs with (d) NZ=4, (e) NZ=10