Bidirectional heterostructures consisting of graphene and lateral MoS2/WS2 composites: a first-principles study.

First-principles calculations have been performed to explore the structural and electronic properties of bidirectional heterostructures composed of graphene and (MoS2) X /(WS2)4-X (X = 1, 2, 3) lateral composites and compare them with those of heterobilayers formed by graphene and pristine MS2 (M = Mo, W). The band gaps of the lateral heterostructures lie between those of pristine MoS2 and WS2. The weak coupling between the two layers can induce a tiny band-gap opening of graphene and formation of an n-type Schottky contact at the G-(MoS2) X /(WS2)4-X interface. Moreover, the combination ratio of MoS2/WS2 can control the electronic properties of G-(MoS2) X /(WS2)4-X . By applying external electric fields, the band gaps of (MoS2) X /(WS2)4-X (X = 0, 1, 2, 3, 4) monolayers undergo a direct-indirect transition, and semiconductor-metal transitions can be found in WS2. External electric fields can also be used effectively to tune the binding energies, charge transfers, and band structures (the types of Schottky and Ohmic contacts) of G-(MoS2) X /(WS2)4-X heterostructures. These findings suggest that G-(MoS2) X /(WS2)4-X heterostructures can serve as high-performance nano-electronic devices. This journal is © The Royal Society of Chemistry.

Introduction

Graphene has been used widely in various fields owing to its electronic, mechanical, and optical properties since it was discovered in 2004.[1,2] However, the zero gap of graphene limits its development and applications.[3] A vertical heterostructure formed by graphene and other two-dimensional matter, e.g., hexagonal boron nitride or transition-metal dichalcogenides (TMDs), via van der Waals forces is a promising method of solving this problem because it could lead to a band-gap opening while preserving the mechanical and electronic structures of graphene.[2-19] Most TMDs, of which the general chemical formula is MX2 (M = Mo, W, V, Co etc., X = S, Se, or Te, etc.), are promising candidates for electronic devices because they have a direct band gap, high carrier mobility, high on/off current ratios, and unique exaction properties.[2,6,12,16-24] Vertical heterostructures composed of a pristine TMDs monolayer and graphene have been the subject of extensive investigation.[2,10,12-16,19,25-28]

Lateral heterostructures have been fabricated within TMD monolayers (MoS2/WS2, WS2/WSe2, MoSe2/WSe2, etc.) and have opened up unprecedented opportunities to improve the performance of electronic devices with tunable electronic properties.[2,29-34] However, due to limitations of the experimental techniques, the electrochemical properties of such in-plane heterostructures have not been fully explored. Although theoretical studies have reported the electronic properties of some lateral heterostructures,[24,35-39] the effects of the component ratio and external electric fields on the properties of such hybrid systems remain unclear. We are left with an interesting question: would any peculiar properties appear in the vertical and lateral combined bidirectional heterostructures formed by graphene and in-plane TMDs composites in comparison with the heterobilayer from graphene and pristine TMDs?

Although extraordinary bidirectional heterostructures of γ-graphyne-MoSe2/WSe2 were investigated by Sun et al.,[35] the G-MoS2/WS2 system has not been studied. Inspired by the above factors, bidirectional heterostructures formed by graphene and lateral MoS2/WS2 composite are explored here on the basis of first-principles calculations. For comparison purposes, the heterobilayers from graphene and pristine MoS2 or WS2 are also considered. Because previous research has demonstrated that external electric fields can be used to tune the band gaps or other electronic structures of vertical heterostructures,[9,10,15,17-19,27,28,40-45] the structural and electronic properties of a (MoS2)/(WS2)4−, (X = 0, 1, 2, 3, 4) monolayer with different proportions and vertical heterobilayers with graphene under external electric fields are systematically investigated by density functional theory (DFT) calculations. The results indicate that external electric fields can induce consequential changes in the properties of a (MoS2)/(WS2)4− monolayer and G-(MoS2)/(WS2)4− heterostructures.

Computational methods

The DMol3 code was used for geometric optimization and electronic properties calculations.[46,47] The generalized gradient approximation (GGA) with a Perdew–Burke–Ernzerhof (PBE) functional and the DFT semi-core pseudo-potentials (DSPP) double numerical atomic basis set plus polarization (DNP) were employed.[48] To show the expected differences with the different function, the on-site Coulomb interaction (GGA + U approach) was adopted to study the electric properties using Vienna ab initio simulation package program (VASP) package.[49,50] Due to the poor description of the weak van der Waals interactions of the popular PBE functional,[56] an empirical dispersion-corrected density functional theory approach proposed by Grimme was used.[51] The optimal geometric convergence criteria of energy, force, and displacement were 1.0 × 10−5 Ha, 0.002 Ha Å−1, and 0.005 Å, respectively. The K-point was set to 9 × 9 × 1 for structural optimization and to 15 × 15 × 1 to calculate electronic properties. The vacuum thickness was set to more than 20 Å to avoid interactions between the neighboring layers. A vertical external electric field in the range of 0 to 1.0 V Å−1 was applied to (MoS2)/(WS2)4− lateral heterostructures and −1.0 to 1.0 V Å−1 to the G-(MoS2)/(WS2)4− heterobilayer, with a positive electric field pointing from (MoS2)/(WS2)4− to graphene and a negative field pointing in the reverse d direction.[39,41] Electron density differences were computed by CASTEP code.[52]

The optimized lattice constants of graphene, MoS2, and WS2 were 2.472 Å, 3.214 Å, and 3.206 Å, respectively. In line with the modeling methods of graphene-MoS2,[16,45] silicene-MoS2,[44] and germanene-MoS2 heterostructures,[43,53] a 4 × 4 substrate layer of (MoS2)/(WS2)4− was chosen to match the 5 × 5 graphene monolayer, producing a tolerable lattice mismatch for this system (∼4.0%). To quantitatively describe the interaction strength between the graphene sheet and the (MoS2)/(WS2)4− monolayer, the binding energies were defined aswhere Etotal, E(MoS, and Egraphene are the energies of the G-(MoS2)/(WS2)4− heterostructures, isolated (MoS2)/(WS2)4−, and graphene monolayers, respectively, and N is the number of carbon atoms in the supercell.

Results and discussion

(MoS2)/(WS2)4− monolayer

The electronic properties of pristine MoS2 or WS2 and the (MoS2)/(WS2)4− lateral heterostructures systems illustrated in Fig. 1, i.e., (MoS2)1/(WS2)3, (MoS2)2/(WS2)2, and (MoS2)3/(WS2)1, were examined first. Fig. 2 presents the calculated electronic band structures of G-(MoS2)/(WS2)4− heterostructures. As shown in Fig. 2, pristine MoS2 and WS2 (x = 0 and 4) have respective direct band gaps at the K point of 1.59 eV and 1.81 eV, which are almost identical to previously reported values.[22,36,37,45] The band gaps of the hybrid lateral heterostructures with different combination MoS2/WS2 ratios lie between those of pristine MoS2 and WS2, and they all maintain the direct band gaps of a pristine TMDs monolayer.

Fig. 1

Top and side views of stable configurations of WS2 (a), (MoS2)1/(WS2)3 (b), (MoS2)2/(WS2)2 (c), (MoS2)3/(WS2)1 (d), MoS2 (e), respectively. The yellow, green, and red balls denote S, Mo and W atoms, respectively.

Fig. 2

Band structures of WS2 (a), (MoS2)1/(WS2)3 (b), (MoS2)2/(WS2)2 (c), (MoS2)3/(WS2)1 (d), and MoS2 (e), respectively. The horizontal red lines represent the Fermi level. (f) Calculated band gap of (MoS2)/(WS2)4− monolayer as a function of electric field intensity.

Although some theoretical studies of MoS2/WS2 lateral heterostructures have been reported in recent years,[24,35-38] the effects of external electric fields on the electronic properties of this lateral hybrid monolayer have been ignored. As the external electric field varies from 0 to 0.6 V Å−1, the bond distance of WS2 monolayer increases with the electric field strength. And the bond distance decreases when the electric field strength more than 0.7 V Å−1. The bond distance of MoS2 was also increases with the electric field strength, but the change is weak (<0.1 Å). The change of bond distance of the lateral hybrid (MoS2)/(WS2)4− (X = 1, 2, 3) heterostructures are not obvious. As shown in Fig. 2, the applied electric fields are associated with a decreasing trend in band gaps for the above systems because field-induced repulsion among the electronic levels results in a shift in conduction band minimum (CBM) toward the Fermi level (EF), consequently, a reduction in band gaps (Fig. S1† and 2).

Furthermore, a direct–indirect band gap transition of (MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) under vertical electric fields has also been observed. The pristine MoS2 and (MoS2)/(WS2)4− lateral heterostructures maintain semiconductor characteristics under an external electric field (F ≤ 1.0 V Å−1). While for WS2, the conduction band overlaps the valence band across the Fermi level at 0.6 V Å−1 to 0.7 V Å−1 (Fig. S1†), and the band gap of WS2 shows a semiconductor–metal transition, which is consistent with a report by Affandi.[54] As shown in Fig. S1,† the band structures of MoS2 are more stable than that of WS2 under the external electronic fields with different intensity. We suspect that the 5d electrons of W may cross the Fermi level under the strong electric field, which would cause the materials to be metallic. This is in good agreement with previous studies.[20,25,54] As a result, the electric properties of (MoS2)/(WS2)4− monolayer with high W concentration become not stable under the electric field and present different trend.

G-(MoS2)/(WS2)4− heterostructures

Vertical hybrid structures of graphene stacked on a (MoS2)/(WS2)4− monolayer were constructed and calculated. A supercell is formed by alternate stacking of graphene sheets on a (MoS2)/(WS2)4− monolayer, with 16 metal (Mo or W) atoms, 32 S atoms, and 50 C atoms. To analyze the optimal stacking patterns, two representative approaching modes were considered. The C atom was placed on top of the metal atoms (W or Mo) or was located above the top of the S atoms (Fig. S2†). The binding energies for different stacking patterns differ by less than 0.1 meV per C atom. This indicates that interlayer bonding takes place irrespective of stacking pattern, which is in agreement with G-MoS2 and Si-MoS2.[16,44] Only the most stable configuration shown in Fig. 3 is considered in the following discussions. The interlayer distances between a (MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) monolayer and graphene are ∼3.2 Å, which suggests that the coupling between them is dominated by van der Waals forces. Graphene exhibits a small wrinkle (∼0.1 Å) in its planar atomic framework for all G-(MoS2)/(WS2)4− heterostructures, similar to the previous calculation.[16,26,55]

Fig. 3

Top and side views of stable configurations of G-WS2 (a), G-(MoS2)1/(WS2)3 (b), G-(MoS2)2/(WS2)2 (c), G-(MoS2)3/(WS2)1 (d), G-MoS2 (e), respectively. The gray balls denote C atom.

Fig. 4 shows the calculated Eb values of G-(MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) heterobilayers. The calculated Eb of G-WS2 is −26.9 meV per C atom, which is nearly identical to previously reported values based on the GGA-PBE function (−26.8 meV and −29 meV).[15,55] The Eb of G-MoS2 is −10.2 meV per C atom, which is smaller than the value (−23 meV) reported by Ma et al. based on local density approximation (LDA) calculations,[16] which is not surprising given that the LDA method overestimates the calculated results. The binding energies of graphene with (MoS2)1/(WS2)3, (MoS2)2/(WS2)2, and (MoS2)3/(WS2)1 lateral heterostructures were −24.1, −18.9, and −14.4 meV, respectively. This clearly shows that the Eb between graphene and the lateral hybrid MoS2/WS2 with different mixing proportions lies between those of graphene with pristine MoS2 or WS2. As expected, binding strength increases with proportion of WS2 units in (MoS2)/(WS2)4−, in the order of G-WS2 > G-(MoS2)1/(WS2)3 > G-(MoS2)2/(WS2)2 > G-(MoS2)3/(WS2)1 > G-MoS2. These Eb values confirm weak van der Waals interactions between the graphene sheet and the (MoS2)/(WS2)4− monolayer, including pristine MX2 and hybrid lateral heterostructures. As shown in the Fig. S3,† the d band of MoS2 has much sharper peak than WS2 around the Fermi level, indicating that the 4d electrons of Mo atoms are more localized than the 5d electrons of W atoms. The more delocalized 5d orbitals leads to the enhancement of interlayer coupling with the concentration of W increases, and consequently results in the strengthen of binding energy.

Fig. 4

Binding energies (Eb) as a function of the electric field intensity for graphene-(MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) heterobilayers.

To probe the effect of external electric fields on the interaction strength of G-(MoS2)/(WS2)4− heterostructures, vertical electric fields in the range of −1.0 V Å−1 to 1.0 V Å−1 were applied to the heterobilayer. The computed Eb as a function of electric fields for G-(MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) is plotted in Fig. 4. The stability order of the heterostructures was maintained regardless of electric field intensity. An applied positive or negative electric field can result in a more negative Eb for the above systems, which suggests that the external electric field can intensify the interaction between the layers and enhance heterobilayer stability. Moreover, the direction of the electric field has a weak influence on Eb, with the effect of the positive field being slightly more pronounced. For G-(MoS2)/(WS2)4− heterobilayers, the change of bond length and the interlayer distance are both less than 0.1 Å under the external electric field. So, there is no obvious structural changes in the heterobilayers under the external electric field. Enhancement of interlayer binding in the heterostructures can be attributed to built-in electric fields that are asymmetrical in the sense that graphene and the (MoS2)/(WS2)4− monolayer have different electron affinities. Field-driven charge transfer among the component layers then enhances the electrostatic interaction between heterobilayers.

To further explore the bonding nature of graphene and the (MoS2)/(WS2)4− monolayer, the charge density difference of the G-(MoS2)/(WS2)4− hybrid systems, indicated in Fig. 5(a–e), was used to visualize electron transfer upon formation of the interlayer. This can be defined asΔwhere ρtotal, ρ(MoS, and ρgraphene are the charge densities for G-(MoS2)/(WS2)4− heterostructures, the isolated (MoS2)/(WS2)4− monolayer, and graphene, respectively. The pink and blue areas denote electron accumulation and depletion, respectively, which clearly demonstrate the charge transfer direction between graphene and (MoS2)/(WS2)4− in the heterobilayers. In these heterostructures, electrons are lost on the graphene side, while they accumulate mainly on the (MoS2)/(WS2)4− side, and a small amount of electrons accumulate on the graphene sheet. The electrons spontaneously transfer from (MoS2)/(WS2)4− monolayer to the graphene sheet. An interface dipole is defined by the potential step, forming in ther interface, is due to the difference in the work function of the MoS2/WS2 and graphene. Although such a value of charge transfer is quite small, it leads to an interface dipoles that cannot be ignored. This is consistent with experimental phenomena.[27,28] As supported by Mulliken population analysis (Fig. 5f), the charge accumulated on the graphene sheet increases with proportion of WS2 units in (MoS2)/(WS2)4−, although the change is small. This small charge transfer suggests a weak van der Waals force between the interlayers, which is consistent with the Eb results.

Fig. 5

Differential electron density of G-WS2 (a), G-(MoS2)1/(WS2)3 (b), G-(MoS2)2/(WS2)2 (c), G-(MoS2)3/(WS2)1 (d), and G-MoS2 (e), respectively. Pink and blue areas correspond to the accumulation and depletion of electronic densities. (f) Mulliken populations of per C atom in graphene for (MoS2)/(WS2)4− as a function of the electric field strength.

The effect of the external electric field on charge redistribution among the G-(MoS2)/(WS2)4− interlayers was further examined. Fig. 5f displays the average Mulliken population per C atom in graphene as a function of electric field strength. The extent of charge transfer is in the order of G-WS2 < G-(MoS2)1/(WS2)3 < G-(MoS2)2/(WS2)2 < G-(MoS2)3/(WS2)1 < G-MoS2 regardless of field strength and direction. The negative electric fields cause more electrons to flow from (MoS2)/(WS2)4− to graphene; thus, the graphene monolayer is more negatively charged with increasing field intensity. When the electric field reaches −1.0 V Å−1, the graphene in G-WS2, G-(MoS2)1/(WS2)3, G-(MoS2)2/(WS2)2, G-(MoS2)3/(WS2)1, and G-MoS2 gains 0.0128 e, 0.0127 e, 0.0124 e, 0.0123 e, and 0.0121 e per C atom, respectively. When the external electric field is applied in the positive direction, the heterostructures exhibit the opposite response: the transferred electrons of G-WS2 gradually decrease to zero at 0.2–0.3 V Å−1, while other G-(MoS2)/(WS2)4− hybrid systems reach zero at 0.1–0.2 V Å−1. A further increase in the strength of the electric field alters the transferred direction of electrons from graphene to the (MoS2)/(WS2)4− monolayer, with graphene being more positively charged. In brief, graphene sheets lose electrons as the electric field changes from negative to positive, reaching 0.0113 e, 0.0120 e, 0.0124 e, 0.0129 e, and 0.0134 e per C atom at F = 1.0 V Å−1 for G-WS2, G-(MoS2)1/(WS2)3, G-(MoS2)2/(WS2)2, G-(MoS2)3/(WS2)1, and G-MoS2, respectively.

The degree of charge transfer between heterobilayers can reflect variation in the interface dipole and the magnitude of the band gap opening. To analyze the impact of the weak van der Waals forces between the interlayer on the band gaps of the hybrid G-(MoS2)/(WS2)4− heterostructures, the calculated gaps are presented in Fig. 6(a–e). Small band gap openings of 4.6, 4.3, 4.3, 4.4, and 2.7 meV for respective stable heterostructures of G-WS2, G-(MoS2)1/(WS2)3, G-(MoS2)2/(WS2)2, G-(MoS2)3/(WS2)1, and G-MoS2 were observed. The derived gaps of G-WS2 and G-MoS2 are consistent with previously reported results (3.9 meV for WS2, 2 meV for MoS2),[16,26] indicating the plausibility of our calculated results. As shown in Fig. S4,† the band structure which used the on-site Coulomb interaction (GGA + U approach) is same to our result in current research. And the difference of band gap is small (2 meV). So, our results are credible. Fig. 6(a–e) shows that the Dirac point of graphene around the Fermi level is well preserved in the band structure of all G-(MoS2)/(WS2)4− heterostructures, which suggests a high carrier mobility of the hybrid systems. As shown in Fig. 6f, when an external electric field is applied, the band gaps at the Dirac point of graphene change only slightly (<0.5 meV) in those heterobilayers, regardless of the field direction.

Fig. 6

Band structures of G-WS2 (a), G-(MoS2)1/(WS2)3 (b), G-(MoS2)2/(WS2)2 (c), G-(MoS2)3/(WS2)1 (d), and G-MoS2 (e), respectively. The insets in figure are a magnification of the area near the linear dispersion band at the K point. The horizontal red lines represent the Fermi level. (f) Band gaps of G-(MoS2)/(WS2)4− heterostructures as a function of the electric field intensity.

External electric fields have a large effect on the Schottky barrier of the G-(MoS2)/(WS2)4− heterostructures in our calculations. When the Fermi level is close to the CBM, known as the n-type Schottky barrier (ΦBn), it is defined as ΦBn = EC − EF. Similarly, the p-type Schottky contact corresponds to the Fermi level close to the valence band maximum (VBM), and the Schottky barrier (ΦBp) is defined as ΦBp = EF − EV, where EV and EF are the CBM and VBM, respectively. The Fermi level was referenced at 0 eV. To quantitatively characterize regulation of the external electric field on the Schottky barrier of the above heterobilayers, the evolutions of Schottky contact type and Schottky barrier height of G-(MoS2)/(WS2)4− heterostructures as a function of electric field intensity are supplied in Fig. 7(f. The corresponding band structures of G-(MoS2)/(WS2)4− heterostructures under external electric fields at 1.0, −0.5, 0.5 and 1.0 V Å−1 are also shown.

Fig. 7

Band structures of G-WS2 (a1−4), G-(MoS2)1/(WS2)3 (b1−4), G-(MoS2)2/(WS2)2 (c1−4), G-(MoS2)3/(WS2)1 (d1−4) and G-MoS2 (e1−4) under electric fields at F = −1.0, −0.5, 0.5 and 1.0 V Å−1 (1–4). The horizontal red lines represent the Fermi level. (fa–fe) Schottky barrier in the G-(MoS2)/(WS2)4− heterostructures as a function of the electric field intensity.

Combining Fig. 6(a–e) with Fig. 7(f, it becomes clear that, without the external electric field, coupling of graphene with the (MoS2)/(WS2)4− monolayer causes the heterostructures to form an n-type Schottky contact. In addition, as the X value increases in G-(MoS2)/(WS2)4− heterostructures, the position of the Fermi level shifts toward the CBM of the (MoS2)/(WS2)4− monolayer and, consequently, the n-type Schottky barrier height decreases. Thus, different ratios of Mo and W in lateral composites can effectively control the Schottky barrier of G-(MoS2)/(WS2)4− heterostructures.

When positive electric fields are applied, the CBM of the (MoS2)/(WS2)4− monolayer shifts down to the Fermi level to reduce the n-type Schottky barrier height gradually as the intensity increases. The ΦBn of the G-(MoS2)/(WS2)4− heterostructures disappears, along with a transition to Ohmic contact to present metallic characteristics, when the electric field intensity increases to 0.7, 0.6, 0.6, 0.5, and 0.5 V Å−1 for G-(MoS2)/(WS2)4− heterostructures (X = 0, 1, 2, 3, 4), respectively. Moreover, the Fermi level falls to below the Dirac point of graphene, leading to p-type (hole) doping in graphene. The trend in variation of the Schottky barrier height of G-MoS2 under positive electric fields is consistent with results reported by Hu et al.[41] As the electric field varies from 0 to −0.1 V Å−1, the n-type Schottky contact transforms to p-type Schottky contact for the G-WS2, G-(MoS2)1/(WS2)3, and G-(MoS2)2/(WS2)2 interfaces. While the applied negative electric field strength is more than −0.1 V Å−1, the VBM of WS2, (MoS2)1/(WS2)3, and (MoS2)2/(WS2)2 above the Fermi level declines, and the ΦBp decreases gradually. Finally, the corresponding heterostructures change into the Ohmic contact when the respective field strength approaches −0.6, −0.7, and −0.7 V Å−1. The Dirac point of graphene decreases to below the Fermi level, resulting in n-type (hole) doping in graphene. For the G-(MoS2)3/(WS2)1 and G-MoS2 interfaces, at a strength less than −0.1 V Å−1, they retain an n-type Schottky contact. Variation of the electric field from −0.1 to −0.2 V Å−1 can induce a transition from n-type to p-type Schottky contacts. Under applied negative electric fields strength greater than −0.2 V Å−1, the VBM of (MoS2)3/(WS2)1 and MoS2 below the Fermi level decreases gradually to reduce ΦBp. When the strength increases to −0.7 and −0.8 V Å−1, conversion from the p-type Schottky contact to the Ohmic contact occurs at the G-(MoS2)3/(WS2)1 and G-MoS2 interface. The graphene's Dirac point crosses the Fermi level and moves into the valence band of (MoS2)3/(WS2)1 and MoS2, which is below the Fermi level, resulting in n-type (hole) doping in graphene. In summary, electric fields can induce an electrostatic potential difference at the G-(MoS2)/(WS2)4− interface and affect the charge transfer process between graphene and the (MoS2)/(WS2)4− monolayer. Applied external electric fields can therefore control contact formation (Schottky and Ohmic contacts) in G-(MoS2)/(WS2)4− heterostructures and tune the type of doping in graphene.

Summary and conclusions

We performed density functional theory calculations to systematically study the structural and electronic properties of (MoS2)/(WS2)4− monolayers and G-(MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) heterobilayers. The impact of external electric fields was also explored. We found that the band gaps of lateral hybrid (MoS2)/(WS2)4− (X = 1, 2, 3) heterostructures lie between those of the pristine MoS2 and WS2. Applied external vertical electric fields can gradually decrease the band gaps and produce a direct–indirect transition. External electric fields can also induce a semiconductor–metal transition in the WS2 monolayer when the field strength is more than 0.6 V Å−1. Interlayer coupling in graphene and the (MoS2)/(WS2)4− (X = 0, 1, 2, 3, 4) monolayer is dominated by weak van der Waals interactions. A tiny band gap opening (<5 meV) at the Dirac point of graphene suggests such hybrid materials have a high carrier mobility. The Eb, charge transfers, and ΦBn of G-(MoS2)/(WS2)4− bidirectional heterostructures increase with proportion of WS2 units in (MoS2)/(WS2)4−, while they all lie between those of G-WS2 and G-MoS2. Furthermore, external electric fields can greatly enhance coupling strength, efficiently tuning the amount and direction of charge transfer between the interlayer. External electric fields can therefore not only control the Schottky barrier height, but also the type of Schottky barrier (n-type and p-type) and contact formation type (Schottky and Ohmic contacts) at the G-(MoS2)/(WS2)4− interface. External electric fields can tune the position of the Dirac point relative to the Fermi level to achieve effective doping (p-type and n-type) in graphene. Due to their electronic properties, G-(MoS2)/(WS2)4− heterostructures provide a significant opportunity to design high-performance nano-electronic devices, such as tunable G-(MoS2)/(WS2)4−-based Schottky barrier diodes and field-effect transistors.

Conflicts of interest

There are no conflicts to declare.