1. Introduction
Two-dimensional (2D) nanomaterials, layered transition metal chalcogenide family, recently have attracted much attention in nanoelectronics and optoelectronic devices, owing to the distinctive electronic and optical properties [1-5]. Bulk ZrS2, is held together with strong covalent in-plane bonds and weak van der Waals forces between the layers, and each layer consists of three atomic layers, forming the hexagonal structure. It is similar to graphene, easily using mechanical or chemical exfoliation to form monolayer and few layer sheets [6, 7]. Recently Zeng groups founded an effective approach to fabricate monolayer ZrS2 nanosheet [8]. Thus monolayer ZrS2 based field effect transistor (FET) devices are fabricated, which have higher electric sensitivity, and provide excellent semiconducting electronic properties [9]. Particularly, focused on the optical properties, two-dimensional (2D) ZrS2 are also suitable for solar cell applications, and present outstanding properties [10]. For the practical applications in optoelectronics, tunable electronic properties of 2D transition metal disulfide (TMD) materials are especially crucial. Recently Yan Li etc. reported many properties of monolayer ZrS2, such as strain effects on the monolayer band structure [11]. External electric field provides, besides strain and doping, another general route to tune the bandgap and crystal structure of the material, thereby inducing a novel physical phenomena. Current studies have confirmed that the electronic structure of graphene and MoS2 bilayer can be tuned in a large range by the applied external electric field [12-14]. Qihang Liu et al. systematically demonstrated the effect of vertical electric field on the electronic structure of MoS2 bilayer [14]. Recently further theoretical studies have been reported, such as applying an external electric field to a rippled MoS2 monolayer [15]. Current studies have confirmed that MoS2 bilayer undergoes semiconductor-metal transition by the applied external electric field [14]. However, about the bilayer ZrS2, the properties on applying external electric field, are rarely reported up to now. Another way of tuning the electronic structure could be by application of feasible normal compressive strain in experiment. Due to relatively weaker Van der Waals interlayer interactions normal compressive strain has been achieved as shown for other bilayer materials [16, 17]. In addition, ZrS2 thin sheets are easily grown with vertically aligned layers, so the vertical pressure is easy to realize in experiments. However, regarding the bilayer ZrS2 the result on applying normal strain has not been reported. The responses of electronic properties under electric field and external strain are interesting issues for discussion in band-gap engineering. In the following, the electronic structure of bilayer ZrS2 with applying the electric field and the compressive strain normal to the bilayer are respectively discussed using first-principles calculations. The physical mechanisms of the band gap decrease are discussed, increasing the applied electric field and strain. The application of electric field and strain on few-layered ZrS2 are expected to tune its electronic transport properties in nanoelectronics and nanophotonics.
2. Computational methods
The first-principles calculations were performed within the density-functional theory (DFT) using the generalized gradient approximation (GGA) proposed by Perdew, Burke, and Ernzerhof for the exchange and correlation energy [18]. Numerical calculation was implemented by the Vienna ab initio simulation package (VASP). Van der Waals interactions in bilayer sheets, which are important in determining the electronic properties with the different interlayer distance, were included in the Grimme's DFT-D2 method [19]. On applying certain compressive strain, the distances between the layers were fixed using selective dynamics in Cartesian coordinates, and the atomic positions and lattice constants under each strain were optimized by using the conjugate gradient method. A plane-wave basis set with kinetic energy cutoff of 600 eV is used. Brillouin zone (BZ) sampling with Monkhorst-Pack(MP) method of 24×24×1 k-points were chosen. The convergence for energy was chosen as 10-6 eV/atom between two steps, and the maximum Hellmann-Feynman atomic force was less than 0.01 eV/Å upon ionic relaxation. The direction of the external electric field is normal to the plane of bilayers, and in VASP, the external uniform electric field is treated with adding an artificial dipole sheet (i.e., dipole correction) in the unit cell [20]. To avoid the interaction between periodic images of slabs in the z-direction, a vacuum of 20 Å is sufficient.
3. Results and discussion
3.1 Electronic properties of bilayer ZrS2
ZrS2 compounds crystallize in the simplest 1T-CdI2 polytype, as shown in Fig. 1(a), and each layer consists of a chalcogen-metal-chalcogen sandwich structure. For ZrS2 bilayer, which can be constructed by pairing two ZrS2 monolayers together, as shown in Fig. 1(b), here we only investigate the AA stacking, which is the most stable configuration and it can be easily obtained by the mechanical exfoliation of bulk ZrS2. For bulk ZrS2, our results about the optimized lattice parameter
3.2 Electronic properties of bilayer ZrS2 with external electric field
The vertical electric field is a powerful tool to achieve tunable electric properties in band gap engineering [21, 22]. To exploit the tunable electric properties of monolayer and bilayer ZrS2 with the external electric field, we have applied electric field perpendicular to the ZrS2 layers along the
In contrast, ZrS2 bilayer is sensitive to the external electric field. With the electric field increasing, the band gap of ZrS2 bilayer can be flexibly decreased, as shown in Fig. 4. Up to 1.4 V/Å, ZrS2 bilayer undergoes a semiconducting-metallic transition. Similar to MoS2, carbon nano-tubes [23] and boron nitride (BN) sheets [24], Stark effect, result in a shift of the bands and change the band structure when applying an external electric field, especially in the conduction region. As shown in Fig. 4, the energy bands such as band A and band B located at point M, are separated from each other entirely. The stronger the electric field is, the larger the band splitting is, as shown in Fig. 4(b), thus the band gap is smaller. With the electric field increasing, at a critical electric field, semiconducting metallic transition occurs.
3.3 Electronic properties of bilayer ZrS2 with vertical compressive strain
Strain provides another general route to tune the bandgap and crystal structure of the material. For ZrS2 monolayer, tuning electronic properties with applying tensile strain has been discussed, with the small tensile strain, its bandgap increases [11, 25]. On the other hand, by application of compressive strain, it has been experimentally and theoretically demonstrated that the band structure of bulk TMDs can be modified. Thus, we also studied the effect of vertical strain on the electronic properties of ZrS2 bilayer. Here we applied compressive strain on the ZrS2 bilayer along the
4. Conclusion
By means of density functional theory with van der Waals correction, the influence of vertical electric field and compressive strain on the electronic properties of bilayer ZrS2 has been systematically investigated. The band gaps of the ZrS2 bilayer decrease monotonically with increasing electric field. At critical electric fields, semiconducting metallic transition occurs. We have also observed that the vertical compressive strain has a significant impact on the electronic properties of ZrS2 bilayer. Overall, the tuning of electronic properties due to external electric field and compressive pressure would promote their potential applications in electronics and optoelectronics.