Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Efficient charge separation and visible-light response in bilayer HfS2-based van der Waals heterostructures

Biao Wanga, Xukai Luoa, Junli Changa, Xiaorui Chena, Hongkuan Yuana and Hong Chen*ab
aSchool of Physical Science and Technology, Southwest University, Chongqing 400715, People’s Republic of China. E-mail: chenh@swu.edu.cn
bKey Laboratory of Luminescent and Real-Time Analytical Chemistry, Ministry of Education, College of Chemistry and Chemical Engineering, Southwest University, Chongqing 400715, People’s Republic of China

Received 9th April 2018 , Accepted 27th April 2018

First published on 23rd May 2018


Abstract

Two-dimensional (2D) hafnium disulfide (HfS2) has been synthesized and is expected to be a promising candidate for photovoltaic applications, and at the same time the hexagonal BN sheet (h-BN) and graphene-like C3N4 sheet (g-C3N4) have also been fabricated and are expected to be applied in photocatalysis. In this work, we employ hybrid density functional theory to investigate HfS2-based van der Waals (vdW) heterojunctions for highly efficient photovoltaic and photocatalytic applications. HfS2/h-BN and HfS2/g-C3N4 heterostructures with direct bandgaps and efficient charge separation are both typical type-II semiconductors and have potential as photovoltaic structures for solar power. Moreover, compared with h-BN and g-C3N4 single-layers, HfS2/h-BN heterostructures with 6% tensile strain and HfS2/g-C3N4 heterostructures with 9% tensile strain have moderate bandgaps, whose optical absorption is obviously enhanced in the ultraviolet-visible (UV-VIS) light range and whose bandedges are suitable for photocatalytic water splitting. HfS2/h-BN heterostructures with 6% applied strain, being different from HfS2/g-C3N4 heterostructures with 9% strain, possess a direct bandgap and show complete separation of the photoinduced electron–hole pairs. Thus the HfS2/h-BN heterojunction with 6% strain has bright prospects for use in visible light photocatalytic water splitting to produce hydrogen.


1 Introduction

Since the discovery of graphene,1 ultrathin 2D nanomaterials2,3 have attracted widespread attention due to their peculiar properties and high utilization efficiency of specific surface areas, which are of significance for potential applications such as photocatalysis. Attributed to the steady improvements made in laboratories recently, many types of monolayer or few-layer transition metal dichalcogenides (TMDCs;4 e.g.WS2,5 MoS2,6 HfS2,7 ZrS2,8 etc.) have been successfully fabricated and have become the focus of researchers.

Because HfS2 is predominantly an ionic crystal with a moderate bandgap and a “two-dimensional” layered structure bonded by weak van der Waals forces,9 it has been inspiring researchers to study its properties from bulk to monolayer.10–12 Recently, few-layered hafnium disulfide nanosheets with a 1T-structure have been successfully made by chemical vapor deposition, which shows their bright prospects in the fields of photodetectors7 and field effect transistors13 due to their ultrafast photoresponse time, high photosensitivity and excellent field effect responses. Furthermore, theoretical calculations suggest that the monolayers of 1H- and 1T-HfS2 have potential applications as photocatalysts for water splitting.14 However, individual HfS2 nanosheets similar to graphitic carbon nitride (g-C3N4)15 are liable to rapidly recombine photogenerated electron–hole pairs and reduce solar conversion and photocatalytic efficiency.

For the sake of reducing bandgaps and enhancing the separation of the photoinduced electron–hole pairs,16 a number of vdW heterojunctions composed of different single-layer materials that are combined together by vdW forces have been manufactured and have been predicted to acquire more desirable properties.2,8,17–24 At present, 2D HfS2/phosphorene heterojunctions25 as field effect transistors and HfS2-based nanocomposites as Z-scheme type photocatalysts for hydrogen production26 have been studied, whose structures and photocatalytic mechanism are more complex than those of type-II heterojunctions. Because of the high cost and difficulty in manufacturing 2D heterojunctions in the laboratory, theoretical research is necessary to predict whether or not they can be used as catalysts for hydrogen evolution by water splitting.26 Moreover, the g-C3N4 and hexagonal boron nitride (h-BN) monolayers have become elementary slabs for vdW heterostructures.27 Here, we mainly investigate whether the bilayer vdW heterostructures with combined HfS2 and h-BN (g-C3N4) layers can form standard type-II heterojunctions, which are deemed to be able to enhance the separation of electron–hole pairs. In addition, their band gaps and band edge positions, the effect of the biaxial strain, the density of states and optical absorption spectra are calculated by density functional theory.

2 Computational method

Our first-principles calculations are carried out using the Vienna ab initio simulation package (VASP)28,29 with the projector augmented wave (PAW) method.30 Electronic exchange and correlation effects are described by the generalized gradient approximation (GGA)31 in the parametrization of Perdew–Burke–Ernzerhof (PBE).32 Since it is well known that PBE generally underestimates the bandgaps of semiconductors, the Heyd–Scuseria–Ernzerhof hybrid functional (HSE06),33 which is considered a more precise and widely applied method, has been chosen to calculate their electronic and optical properties. The vdW interactions are taken into account in these heterojunctions by employing Grimme’s DFT-D3 method.34 A default mixing parameter (α = 0.25) is adopted to compute the exchange correlation energy in our work. In the following computations, a Monkhorst–Pack mesh35 for k-point sampling is set as 7 × 7 × 1. A cutoff energy of 550 eV is chosen. The convergence criteria of the geometry optimizations are set to less than 10−5 eV for the total energy and 0.01 eV Å−1 for the force on each ion. A vacuum space of 20 Å is inserted to avoid the interactions between adjacent heterojunctions. The vacuum space is used to evaluate the absolute position of the band compared with that of the standard hydrogen electrode potential. Therefore, the absolute bandedge positions are computed via subtracting the vacuum level by means of averaging the LOCPOT file.

3 Results and discussion

3.1 Structural stability

HfS2 monolayers in a 1T-phase have been theoretically predicted to be globally stable.14 Moreover, 1T-phase few-layered HfS2 nanosheets have been successfully fabricated in experiments.7 Therefore, we adopt the 1T-phase of 2D HfS2 in the following calculations. Our calculated optimized lattice constant of the HfS2 monolayer is 3.61 Å, which matches well with the previous experimental value of 3.62 Å.10 As shown in Fig. 1, the geometric configurations of these heterojunctions are composed of a 2 × 2 HfS2 (a = 3.61 Å) single-layer supercell, and a unit cell of g-C3N4 (a = 7.13 Å)36 nanosheet or a 3 × 3 h-BN (a = 2.50 Å)37 single-layer supercell (Fig. 1). There exists lattice mismatch, which is illustrated by this formula: ε = (lotherslHfS2)/lHfS2, where lothers represents the unit cell of the g-C3N4 single-layer (a = 7.13 Å) and the 3 × 3 h-BN supercell (a = 7.50 Å), and lHfS2 corresponds to the 2 × 2 HfS2 supercell (a = 7.22 Å). As listed in Table 1, the lattice mismatches are −1.16% and 3.89% for g-C3N4 and h-BN, which indicates that the HfS2 nanosheet meets the geometric lattice-matching demand to form the vdW heterojunctions.
image file: c8ra03047b-f1.tif
Fig. 1 The optimized geometry of the vdW heterojunctions. The pink, blue, gray, yellow and sky blue balls indicate B, N, C, S and Hf atoms, respectively. Top: top view of the heterojunctions, bottom: side view of the heterojunctions. Left: HfS2/h-BN heterojunction, right: HfS2/g-C3N4 heterojunction.
Table 1 The geometric parameters, binding energies and work functions of the heterojunctions
Structure ε (%) R (Å) Eb (eV) Φ (eV)
HfS2/h-BN 3.89 3.40 −0.11 6.09
HfS2/g-C3N4 −1.16 3.35 −0.58 6.74


As shown in Fig. 1, the equilibrium distances of the fully relaxed nanocomposites, marked as R in Table 1, are 3.40 Å and 3.35 Å for HfS2/h-BN and HfS2/g-C3N4 heterojunctions, respectively, and represent the space between the S atom in the HfS2 layer and the other layer. The equilibrium distances of the studied heterostructures are close to those of other vdW heterostructures,8,24 which are typical distances of vdW heterostructures. Therefore, the vdW correction of interactions has been applied when optimizing the structure and calculations. The binding energy (Eb) is defined by the equation: Eb = Ex/HfS2 – (Ex + EHfS2), where Ex/HfS2, Ex and EHfS2 represent the total energy of the nanocomposites, the independent single-layers (x = h-BN or g-C3N4) and the isolated HfS2 monolayer, respectively. The computed binding energies of the HfS2/h-BN and HfS2/g-C3N4 heterojunctions (Table 1) are −0.11 eV and −0.58 eV respectively, which means that these nanocomposites are stable.

3.2 Band structure and density of states

Because Hf is a heavy element, the spin-orbital coupling (SOC) correction has been investigated in previous research and it was discovered that the SOC effect in the HfS2 single-layer is small.26 The band structure of a HfS2 monolayer with SOC correction is calculated and presented in Fig. 2. Furthermore, we examine the effect of SOC on the HfS2 monolayer and related heterojunctions, and find that its influences on their bandgaps are all less than 0.1 eV. Therefore, in consideration of the heavy price of calculations, the SOC correction isn’t employed in the following calculations. Using the HSE06 functional, the calculated bandgap of the HfS2 monolayer is 1.99 eV, which approximates to the previous computed value of 2.06 eV.12,26 The bandedge of the HfS2 monolayer can’t straddle the water reduction potential shown in Fig. 2, which means that the HfS2 single-layer isn’t able to be directly used for photocatalytic hydrogen production from water splitting. The calculated bandgaps of the g-C3N4 and h-BN monolayers are 2.76 eV and 5.69 eV, respectively. Moreover, the catalytic activity depends on the related work function (Φ), which represents the minimum quantity of work required to remove an electron from the interior of materials to the vacuum level. Work function is defined by this equation: Φ = EvEF, where Ev and EF refer to the vacuum energy and the Fermi energy. The calculated values of the work function for HfS2/h-BN and HfS2/g-C3N4 heterojunctions are 6.09 and 6.74 eV, respectively. The values of the work function and the photocatalytic activity have an inverse relationship,14 which means that the catalytic activity of the HfS2/h-BN heterostructure is better than that of the HfS2/g-C3N4 heterojunction.
image file: c8ra03047b-f2.tif
Fig. 2 The band structure of a HfS2 monolayer with SOC correction.

HfS2/h-BN and HfS2/g-C3N4 heterostructures are both direct bandgap semiconductors, whose bandgaps have been observed to reduce to 1.78 eV and 1.29 eV, as presented in Table 2. Because of these suitable bandgaps which are less than 1.8 eV, these heterostructures can effectively take advantage of solar power. It’s well known that the redox potentials of water are connected to the pH value.38 The reduction potential for H+/H2 is EH+/H2 = −4.44 eV + pH × 0.059 eV, and the oxidation potential for O2/H2O is EO2/H2O = −5.67 eV + pH × 0.059 eV. As shown in Table 2, Eg, EVBM and ECBM denote the bandgap and the energy levels of the valence band maximum (VBM) and the conduction band minimum (CBM) compared to those of the vacuum level, respectively. Unfortunately, the CBMs of these heterojunctions shown in Fig. 3 aren’t higher than the reduction potential of water at pH = 0, which indicates that these heterostructures may not be used for water splitting to produce hydrogen. However, the light-generated electrons and holes are completely separated in HfS2/h-BN and HfS2/g-C3N4 heterostructures. Fig. 4 shows the total density of states (TDOS) and partial density of states (PDOS). The CB in the HfS2/h-BN heterojunction is composed of Hf 5d and S 3p states that originate from the HfS2 slab, while the VB of the heterostructure largely consists of the N 2p state which stems from the h-BN layer. Thus, the CB and VB in the heterojunction are seated in two different slabs. Similarly, the CB in HfS2/g-C3N4 is composed of Hf 5d and S 3p states located in the HfS2 layer, while the VB is composed of C 2p and N 2p states located in the other layer. Therefore, both of the heterostructures form typical type-II semiconductors, which can improve the efficiency of utilizing sunlight. With moderate bandgaps and full separation of photogenerated electron–hole pairs, HfS2/h-BN (g-C3N4) heterostructures have profound development potential as photovoltaic materials for solar power, except for in water splitting.

Table 2 Bandgaps and bandedges of the heterojunctions at pH = 0
Structure Eg (eV) EVBM (eV) ECBM (eV) Bandgap Straddling water redox potentials
HfS2/h-BN 1.78 −6.51 −4.74 Direct No
HfS2/h-BN with −6% strain 0.62 −5.77 −5.14 Direct No
HfS2/h-BN with −3% strain 1.27 −6.21 −4.94 Direct No
HfS2/h-BN with 3% strain 2.03 −6.63 −4.60 Direct No
HfS2/h-BN with 6% strain 2.28 −6.71 −4.43 Direct Yes
HfS2/h-BN with 9% strain 2.52 −6.80 −4.28 Direct Yes
HfS2/g-C3N4 1.29 −7.19 −5.90 Direct No
HfS2/g-C3N4 with −6% strain 0.01 −5.34 −5.33 Direct No
HfS2/g-C3N4 with −3% strain 0.38 −5.60 −5.22 Direct No
HfS2/g-C3N4 with 3% strain 1.92 −6.65 −4.73 Indirect No
HfS2/g-C3N4 with 6% strain 2.29 −6.88 −4.59 Direct No
HfS2/g-C3N4 with 9% strain 2.72 −7.15 −4.43 Indirect Yes



image file: c8ra03047b-f3.tif
Fig. 3 Band structures of HfS2/h-BN (a) and HfS2/g-C3N4 (b) heterojunctions.

image file: c8ra03047b-f4.tif
Fig. 4 TDOS and PDOS for HfS2/h-BN (a) and HfS2/g-C3N4 (b) heterojunctions.

In previous studies it was confirmed that the in-plane biaxial strains can effectively adjust the optical absorption and electronic properties of vdW heterostructures.8,24 In order to explore their potential application in the field of photocatalytic water splitting, HfS2/h-BN and HfS2/g-C3N4 heterostructures with biaxial strains have been comprehensively investigated. During relaxation of the geometric structure, the lattice constant lx is fixed to imitate the strain on the heterostructures. Strains are defined as follows: η = (lxl0)/l0, where l0 and lx represent the lattice of the original heterostructures and strained heterojunctions. Employing the in-plane biaxial compressive and tensile strain on these nanocomposites, the value of η ranges from −6% to 9%, as shown in Fig. 5, where the bandedge positions are presented as a function of the applied strains. With increasing tensile strains, the lattices of the heterostructures become larger, which indicates that the distance of the nearest atoms is enlarged. On the basis of tight-binding theory,39 their energy-bands will narrow down and the bandgaps will get bigger. Consequently, the enlargement of the bandgaps is accompanied by an increase in the biaxial strain as shown in Fig. 5. Due to their porosity in the g-C3N4 layer,40 the bandedge variation with strain in the HfS2/g-C3N4 heterojunction is different from that in the HfS2/h-BN heterojunction. Applying 6% or 9% tensile strain, the bandedge of the HfS2/h-BN heterostructure can straddle water redox potentials. Similarly, that of the HfS2/g-C3N4 heterostructure with 9% tensile strain can straddle water redox potentials as well. However, HfS2/g-C3N4 with increasing pressure, whose bandgap approaches to zero, can transform from semiconductor to metal. With the aim of meeting the requirements for photocatalytic water splitting, a HfS2/h-BN heterostructure with 6% strain and HfS2/g-C3N4 heterostructure with 9% strain were investigated using the following calculations.


image file: c8ra03047b-f5.tif
Fig. 5 The bandedge positions of the HfS2/h-BN heterojunction (a) and HfS2/g-C3N4 heterojunction (b) as a function of biaxial strain.

The HfS2/h-BN heterojunction with 6% applied tensile strain is a direct bandgap semiconductor, whose bandgap is calculated to be 2.28 eV. The VBM and CBM of the heterojunction are both at the G point, as plotted in Fig. 6. Meanwhile the HfS2/g-C3N4 heterostructure with 9% applied tensile strain, whose VBM is at the G point and CBM is between the G and K points, is an indirect bandgap semiconductor with a bandgap of about 2.72 eV. Considering the moderate and direct bandgap, the photocatalytic performance of the HfS2/h-BN heterojunction with 6% strain is superior to that of the other heterostructures.


image file: c8ra03047b-f6.tif
Fig. 6 The band structures of the HfS2/h-BN heterostructure with 6% strain (a) and HfS2/g-C3N4 heterostructure with 9% strain (b) compared to the vacuum level.

The TDOS and PDOS of the HfS2/h-BN heterostructure with 6% strain and HfS2/g-C3N4 heterostructure with 9% strain are presented in Fig. 7. The CB of the HfS2/h-BN heterostructure with 6% strain is mainly composed of Hf 5d, Hf 6s and S 3p orbitals derived from the HfS2 layer, while the VB of the heterojunction is primarily composed of the N 2p state stemming from the h-BN layer. Obviously, this heterojunction is a classic type-II semiconductor, where the CB and VB are located on different layers in favor of complete separation of the electron–hole pairs and separate production of hydrogen and oxygen. For the HfS2/g-C3N4 heterostructure with 9% strain, the VB of the heterojunction is composed of N 2p and S 3p orbitals that originate from both layers, and the CB is composed of the Hf 5d, Hf 6s and S 3p orbitals. As a result, the HfS2/g-C3N4 heterostructure with 9% strain is incapable of completely separating the photoinduced carriers. Therefore, considering the efficient division of the electron–hole pairs, the HfS2/h-BN heterostructure with 6% strain is more suitable for photocatalysis than the HfS2/g-C3N4 heterostructure with 9% strain.


image file: c8ra03047b-f7.tif
Fig. 7 The TDOS and PDOS for the HfS2/h-BN heterostructure with 6% strain (a) and HfS2/g-C3N4 heterostructure with 9% strain (b).

3.3 Charge transfer analysis

Since strengthened photocatalytic performance arises from the effective separation and transfer of the photogenerated carriers,41 the three dimensional (3D) charge density differences at the interfaces of the heterojunctions are depicted in Fig. 8. This can be calculated using the following equation: ρdif = ρHfS2/YρHfS2ρY, where Y and HfS2/Y denote h-BN or g-C3N4 and HfS2-based heterojunctions, respectively. As shown in the figure, a great deal of electrons are redistributed in the interfaces of the heterojunctions, which leads to built-in electric fields and facilitates separation of the photogenerated electron–hole pairs. For the sake of quantitative analysis of the charge variation at these interfaces, Bader charge analysis42 has been applied to these heterojunctions. Based on the calculations, the Hf, B and C atoms are inclined to deplete electrons, while the S and N atoms are inclined to aggregate electrons. Moreover, the calculated quantity of the charge transfer from the h-BN and g-C3N4 layer (the top layer in Fig. 8) to the HfS2/layer (the bottom layer in Fig. 8) is 0.0056 e and 0.0165 e, respectively. The calculated quantity of the charge transfer is in the same order of magnitude as that for typical vdW heterostructures.43 Obviously, a great deal of charge accumulation is observed in the interface of the lower layer shown in Fig. 8, while the charge depletion appears on the upper layer. Therefore, our work reveals that these HfS2-based heterojunctions can effectively enhance the separation and transfer of electron–hole pairs and then improve the photocatalytic performance.
image file: c8ra03047b-f8.tif
Fig. 8 3D charge density differences in the HfS2/h-BN heterojunction with 6% applied strain (a) and HfS2/g-C3N4 heterojunction with 9% applied strain (b) with an isovalue of 0.00012 e Å−3. The yellow and sea-green isosurfaces represent charge aggregation and depletion, respectively.

3.4 Optical properties

The calculated optical absorption spectra of the g-C3N4, h-BN and HfS2 single-layers and related heterojunctions as a function of wavelength are illustrated in Fig. 9. Compared with the optical absorption of the h-BN and g-C3N4 monolayer, the optical absorption intensity of the corresponding HfS2-based heterojunctions is dramatically strengthened in the ultraviolet-visible (UV-VIS) light region, which is in accordance with reduction of the band gap. Whether the HfS2/h-BN heterostructures have applied strain or not, the optical absorption edges all shift to the visible light range (more than 500 nm), as plotted in Fig. 9.
image file: c8ra03047b-f9.tif
Fig. 9 Optical properties of HfS2/h-BN heterojunctions (a) and HfS2/g-C3N4 heterojunctions (b).

Moreover, the optical absorption spectra of the HfS2/h-BN heterojunction (with 6% strain) depicted in Fig. 9(a) are obviously enhanced in the ultraviolet range in contrast to those of the HfS2 single-layer. Compared with the optical absorption spectrum of the h-BN monolayer, it can only be speculated that the significant red shift of the HfS2/h-BN heterojunction’s absorption spectrum (with 6% strain) is mainly induced by a transition from the N 2p to Hf 5d orbital, as shown in Fig. 5(a) (Fig. 7(a)), which is important for strengthening the heterojunction’s photocatalytic activity. Consequently, with an ideal band gap, complete charge separation and visible-light response, the HfS2/h-BN heterojunction can be used in the field of photovoltaics. On account of its moderate bandgap, well positioned band edge, effective charge separation and reinforced UV-VIS light absorption, the HfS2/h-BN heterojunction with 6% strain can be used for photocatalytic water splitting.

As shown in Fig. 9(b), the optical absorption edge of the HfS2/g-C3N4 heterojunction (with 9% strain) shifts to a longer wavelength range. In particular, the HfS2/g-C3N4 heterojunction takes full advantage of visible light, which is similarly derived from a transition from the N 2p to Hf 5d orbital, as shown in Fig. 5(b). Because of its relatively narrow bandgap, efficient charge separation and visible-light response, the HfS2/g-C3N4 heterojunction is a promising material for photovoltaic applications. With an indirect and moderate bandgap, the HfS2/g-C3N4 heterojunction with 9% applied strain has a well positioned band edge but isn’t able to completely separate charge carriers. Therefore, the photocatalytic performance of the HfS2/h-BN heterojunction with 6% strain is superior to that of the HfS2/g-C3N4 heterojunction with 9% strain.

4 Conclusions

In conclusion, we have systematically studied the HfS2/h-BN and HfS2/g-C3N4 heterojunctions by using the HSE06 hybrid functional. We find that these heterostructures are both direct bandgap semiconductors, whose bandgaps have been observed to reduce to 1.78 eV and 1.29 eV, respectively. The CB and VB of the heterostructures are both derived from different layers, effectively avoiding recombination of the photoinduced electron–hole pairs. Owing to an ideal band gap, complete separation of the photoinduced charge carriers and a visible-light response, these heterojunctions can be widely used in the photovoltaic field. However, the bandedges of these heterojunctions are incapable of straddling water redox potentials.

In addition, applying biaxial strain can effectively adjust the electronic properties of HfS2/h-BN and HfS2/g-C3N4 heterojunctions. With pressure of a certain extent, the bandgap of the HfS2/g-C3N4 heterojunction can reduce to 0 eV, which may indicate that the heterojunction transforms from semiconductor to metal. The band edges of HfS2/h-BN composites with 6% tensile strain and HfS2/g-C3N4 heterojunctions with 9% tensile strain are both able to straddle water redox potentials. The photocatalytic performance of a HfS2/h-BN heterojunction with 6% strain, which is a typical type-II semiconductor, has advantages over that of a HfS2/g-C3N4 heterojunction with 9% strain. Because of the moderate band gap, effective charge separation, well positioned band edge and reinforced UV-VIS light absorption, HfS2/h-BN heterojunctions with applied strains are a promising photocatalyst for water splitting. These calculations pave the way for exploiting high-performance HfS2-based photocatalysts.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

This work was financially supported by the National Natural Science Foundation of China (Grant No. 11645002 and 10904125), the Natural Science Foundation of Chongqing (Grant No. CSTC-2011BA6004 and CSTC-2008BB4253) and the Fundamental Research Funds for the Central Universities (Grant No. XDJK2017C062).

References

  1. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva and A. A. Firsov, Science, 2004, 306, 666–669 CrossRef PubMed.
  2. F. Bonaccorso, L. Colombo, G. H. Yu, M. Stoller, V. Tozzini, A. C. Ferrari, R. S. Ruoff and V. Pellegrini, Science, 2015, 347, 41–51 CrossRef PubMed.
  3. H. Zhang, ACS Nano, 2015, 9, 9451–9469 CrossRef PubMed.
  4. Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman and M. S. Strano, Nat. Nanotechnol., 2012, 7, 699–712 CrossRef PubMed.
  5. J. N. Coleman, M. Lotya and V. Nicolosi, et al., Science, 2011, 331, 568–571 CrossRef PubMed.
  6. X. L. Fan, Y. Yang, P. Xiao and W. M. Lau, J. Mater. Chem. A, 2014, 2, 20545–20551 Search PubMed.
  7. B. J. Zheng, Y. F. Chen, Z. G. Wang, F. Qi, Z. S. Huang, X. Hao and P. J. Li, 2D Mater., 2016, 3, 035024 CrossRef.
  8. X. R. Zhang, Z. S. Meng, D. W. Rao, Y. H. Wang, Q. Shi, Y. Z. Liu, H. P. Wu, K. M. Deng, H. Y. Liu and R. F. Lu, Energy Environ. Sci., 2016, 9, 841–849 Search PubMed.
  9. F. Mehran, T. D. Schultz, M. W. Shafer and W. J. Fitzpatrick, Phys. Rev. Lett., 1974, 33, 1085–1088 CrossRef.
  10. D. L. Greenaway and R. Nitsche, J. Phys. Chem. Solids, 1965, 26, 1445–1458 CrossRef.
  11. C. Kreis, S. Werth, R. Adelung, L. Kipp, M. Skibowski, E. E. Krasovskii and W. Schattke, Phys. Rev. B: Condens. Matter Mater. Phys., 2003, 68, 235331 CrossRef.
  12. J. Kang, H. Sahin and F. M. Peeters, Phys. Chem. Chem. Phys., 2015, 17, 27742 RSC.
  13. S. H. Chae, Y. Jin, T. S. Kim, D. S. Chung, H. Na, H. Nam, H. Kim, D. J. Perello, H. Y. Jeong, T. H. Ly and Y. H. Lee, ACS Nano, 2016, 10, 1309–1316 CrossRef PubMed.
  14. D. Singh, S. K. Gupta, Y. Sonvane, A. Kumarc and R. Ahujad, Catal. Sci. Technol., 2016, 6, 6605–6614 Search PubMed.
  15. Z. W. Zhao, Y. J. Sun and F. Dong, Nanoscale, 2015, 7, 15–37 RSC.
  16. H. L. Wang, L. S. Zhang, Z. G. Chen, J. Q. Hu, S. J. Li, Z. H. Wang, J. S. Liu and X. C. Wang, Chem. Soc. Rev., 2014, 43, 5234–5244 RSC.
  17. J. E. Padilha, A. Fazzio and A. J. R. da Silva, Phys. Rev. Lett., 2015, 114, 066803 CrossRef PubMed.
  18. K. S. Novoselov, A. Mishchenko, A. Carvalho and A. H. Neto, Science, 2016, 353, 461–473 CrossRef PubMed.
  19. D. H. Deng, K. S. Novoselov, Q. Fu, N. F. Zheng, Z. Q. Tian and X. H. Bao, Nat. Nanotechnol., 2016, 11, 218–230 CrossRef PubMed.
  20. J. Y. Liu, Q. Q. Huang, K. Zhang, Y. Y. Xu, M. Z. Guo, Y. Q. Qian, Z. G. Huang, F. C. Lai and L. M. Lin, Nanoscale Res. Lett., 2017, 12, 259 CrossRef PubMed.
  21. F. P. Li, W. Wei, P. Zhao, B. B. Huang and Y. Dai, J. Phys. Chem. Lett., 2017, 8, 5959–5965 CrossRef PubMed.
  22. N. Wang, D. Cao, J. Wang, P. Liang, X. S. Chen and H. B. Shu, J. Mater. Chem. C, 2017, 5, 9687–9693 RSC.
  23. X. K. Luo, G. Z. Wang, Y. H. Huang, B. Wang, H. K. Yuan and H. Chen, Phys. Chem. Chem. Phys., 2017, 19, 28216–28224 RSC.
  24. B. Wang, A. L. Kuang, X. K. Luo, G. Z. Wang, H. K. Yuan and H. Chen, Appl. Surf. Sci., 2018, 439, 374–379 CrossRef.
  25. J. Chang, J. Appl. Phys., 2015, 117, 214502 CrossRef.
  26. C. F. Fu, Q. Q. Luo, X. X. Liao and J. L. Yang, J. Mater. Chem. A, 2016, 4, 18892–18898 Search PubMed.
  27. A. K. Geim and I. V. Grigorieva, Nature, 2013, 499, 419–425 CrossRef PubMed.
  28. G. Kresse and J. Furthmüller, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 11169–11186 CrossRef.
  29. G. Kresse and D. Joubert, Phys. Rev. B: Condens. Matter Mater. Phys., 1999, 59, 1758–1775 CrossRef.
  30. P. E. Blöchl, Phys. Rev. B: Condens. Matter Mater. Phys., 1994, 540, 17953–17979 CrossRef.
  31. J. A. White and D. M. Bird, Phys. Rev. B: Condens. Matter Mater. Phys., 1994, 50, 4954–4957 CrossRef.
  32. M. Ernzerhof and G. E. Scuseria, J. Chem. Phys., 1999, 110, 5029–5036 CrossRef.
  33. J. Heyd and G. E. Scuseria, J. Chem. Phys., 2003, 118, 8207–8215 CrossRef.
  34. S. Grimme, J. Antony, S. Ehrlich and S. Krieg, J. Chem. Phys., 2010, 132, 154104 CrossRef PubMed.
  35. H. J. Monkhorst and J. D. Pack, Phys. Rev. B: Solid State, 1976, 13, 5188–5192 CrossRef.
  36. X. C. Wang, K. Maeda, A. Thomas, K. Takanabe, G. Xin, J. M. Carlsson, K. Domen and M. Antonietti, Nat. Mater., 2009, 8, 76C–80C CrossRef PubMed.
  37. O. O. Kurakevych and V. L. Solozhenko, Acta Crystallogr., Sect. C: Cryst. Struct. Commun., 2007, 63, 80–82 Search PubMed.
  38. V. Chakrapani, J. C. Angus, A. B. Anderson, S. D. Wolter, B. R. Stoner and G. U. Sumanasekera, Science, 2007, 318, 1424–1430 CrossRef PubMed.
  39. W. A. Harrison, Solid State Commun., 2002, 124, 443–447 CrossRef.
  40. D. M. Liang, T. Jing, Y. C. Ma, J. X. Hao, G. Y. Sun and M. S. Deng, J. Phys. Chem. C, 2016, 120, 24023–24029 Search PubMed.
  41. S. S. Ding, W. Q. Huang, Y. C. Yang, B. X. Zhou, W. Y. Hu, M. Q. Long, P. Peng and G. F. Huang, J. Appl. Phys., 2016, 119, 205704 CrossRef.
  42. W. Tang, E. Sanville and G. Henkelman, J. Phys.: Condens. Matter, 2009, 21, 084204 CrossRef PubMed.
  43. L. Xu, W. Q. Huang, L. L. Wang, Z. A. Tian, W. Y. Hu, Y. M. Ma, X. Wang, A. L. Pan and G. F. Huang, Chem. Mater., 2015, 27, 1612–1621 CrossRef.

This journal is © The Royal Society of Chemistry 2018
Click here to see how this site uses Cookies. View our privacy policy here.