Van der Waals metal-semiconductor junction: Weak Fermi level pinning enables effective tuning of Schottky barrier

See allHide authors and affiliations

Science Advances  22 Apr 2016:
Vol. 2, no. 4, e1600069
DOI: 10.1126/sciadv.1600069


Two-dimensional (2D) semiconductors have shown great potential for electronic and optoelectronic applications. However, their development is limited by a large Schottky barrier (SB) at the metal-semiconductor junction (MSJ), which is difficult to tune by using conventional metals because of the effect of strong Fermi level pinning (FLP). We show that this problem can be overcome by using 2D metals, which are bounded with 2D semiconductors through van der Waals (vdW) interactions. This success relies on a weak FLP at the vdW MSJ, which is attributed to the suppression of metal-induced gap states. Consequently, the SB becomes tunable and can vanish with proper 2D metals (for example, H-NbS2). This work not only offers new insights into the fundamental properties of heterojunctions but also uncovers the great potential of 2D metals for device applications.

  • Two-dimensional materials
  • heterojunctions


Metal-semiconductor junction (MSJ) is a critical component in electronic and optoelectronic devices. One of the most important parameters for the MSJ is the Schottky barrier (SB) height (Φ), an energy barrier height for charge carrier transport across the junction, which has a significant impact on device performance (1, 2). The importance of tuning Φ has recently been emphasized in devices based on two-dimensional (2D) semiconducting transition metal dichalcogenides (MX2, where M = Mo or W, and X = S, Se, or Te), which have attracted intense interest because of their unique properties and promising applications (35). Tuning Φ could enable a variety of improvements, such as reducing contact resistance, modulating carrier polarity in the channel for transistors, and enhancing selectivity of carrier extraction for photovoltaic cells. The first one is of particular importance in matching the performance with Si or III-V device analogs (6, 7).

However, it is rather difficult to tune Φ for 2D MX2 by using different common metals because of the effect of Fermi level pinning (FLP). The reason is that Φ is generally determined by the energy difference between the Fermi level (FL) and the semiconductor band edges in the junctionEmbedded Image(1)where Φe and Φh are the barrier heights for electrons and holes, respectively, EF is the Fermi energy, and EVBM and ECBM are the energy of the valence band maximum (VBM) and conduction band minimum (CBM), respectively, of the semiconductor in the junction. Neglecting the metal-semiconductor interaction, Φ should ideally follow the predictions of the Schottky-Mott modelEmbedded Image(2)where W is the work function of the metal, and Eea and Eip are the electron affinity and ionization potential of the semiconductor, respectively. EipEea = band gap. These quantities are the intrinsic properties of isolated materials before they form the junction. The Schottky-Mott model suggests that Φe (or Φh) is linearly dependent on the W of metals with a slope of ±1. However, in reality, Φ is usually insensitive to W, and the FL of the system is pinned to a fixed position in the semiconductor band gap, varying little with respect to the different metals used. The strength of FLP for a given semiconductor with a set of metals can be characterized by the linearly fitted slope of the Φ-versus-W plot (2)Embedded Image(3)If S = 1, the Schottky-Mott limit is recovered. Unfortunately, S is small for 3D metal–2D semiconductor junctions. Experiments measured S = 0.1 for 3D metal–MoS2 junctions (8), indicating a strong FLP.

In general, there exist many different models for explaining FLP [see the study by Tung (2) for a review], and no simple equations are applicable for all types of MSJs (2). Here, we focus on those that have been recognized to be relevant to 2D semiconductors. (i) Metal-induced gap states (MIGS) are formed in the semiconductor (912), as also observed in many 3D metal–3D semiconductor junctions (2, 1319). These states serve as a reservoir for electrons or holes and therefore pin the FL. (ii) An interface dipole formed by the charge redistribution at the interface can shift the electronic levels from their original positions, leading to a deviation from the Schottky-Mott limit (2, 9, 20, 21). (iii) Defects at the interface (created during material/device fabrication) could generate gap states that pin the FL (1, 2, 22), yet this can be neglected for a high-quality interface.

Here, we show that contrary to the conventional 3D metal–2D semiconductor junctions, FLP is weak for the MSJ formed by van der Waals (vdW) interactions, which is attributed to the suppression of MIGS in the semiconductor. This phenomenon allows for the tuning of Φ by using different 2D metals. On this basis, we identify the promising 2D metals that can form a low-Φ junction with 2D semiconductors.


We consider a wide range of 2D metals that have been experimentally realized, including the following: (i) triangular (T) and distorted triangular (T′) phases of MoX2 and WX2; (ii) group 5 MX2 (M = V, Nb, or Ta) and TiX2, which exhibit either hexagonal (H) or T phase at room temperature; and (iii) pristine and doped graphene. These materials have chemically saturated surfaces and bind with 2D semiconductors through vdW interactions.

Figure 1A shows a typical band structure of a 2D MSJ by using T-MoS2H-MoS2 as an example (T-MoS2 is a metal and H-MoS2 is a semiconductor; see computational details in the Supplementary Materials). The origin of each electronic state can be determined from its charge density distribution, based on which we identify the VBM and the CBM of the semiconductor in the junction (as marked by the blue and red dots, respectively, in Fig. 1A). We then plot the distribution of the states between the VBM and the CBM in the junction (Fig. 1A) and find that these states come from the metal rather than the semiconductor. Therefore, the MIGS in the semiconductor of the junction are negligible, which is further evidenced by projecting density of states (DOS) onto the semiconductor (Fig. 1A). This feature is in contrast to 3D metal–2D semiconductor junctions, where MIGS are significant (912). We attribute this phenomenon to the unique interaction between the metal and the semiconductor at the vdW MSJ, which is substantially weaker than the chemical bonding in other junctions (9, 12). Indeed, our calculations show that the distance between the 2D metal and the 2D semiconductor falls within 3 to 4 Å (depending on the computational methods and materials), a typical distance for vdW interactions. It is interesting to notice that the interaction can also be weakened by inserting a thin layer between the metal and the semiconductor to break their direct chemical bond (2325). However, in these cases, a proper separation material has to be chosen first. This need is avoided in our case, which also simplifies device fabrication.

Fig. 1 Weak FLP at the vdW MSJ.

(A) Electronic band structure and DOS of a typical vdW MSJ. T-MoS2H-MoS2 is used here as an example, where T-MoS2 is a metal and H-MoS2 is a semiconductor (SC). FL is set to zero. The VBM and CBM of the SC are marked by blue and red dots, respectively. The purple line shows the DOS projected on the semiconductor. Isosurfaces show the spatial distributions of the states. (B) SB heights (Φe for electrons and Φh for holes) between H-MoS2 and various 2D metals. The diagonal line shows the values predicted by the Schottky-Mott model. (C) Charge density change (averaged in the plane parallel to the interface) after the formation of the H-VS2H-MoS2 MSJ. Atom positions perpendicular to the basal plane are shown as distances relative to the center of the interface and marked by horizontal lines. a.u., arbitrary units.

The suppression of MIGS in the semiconductor can therefore lead to a weak FLP. To verify, Fig. 1B shows the Φ value (calculated by using Eq. 1) for H-MoS2 with various 2D metals. We choose to study T-MoS2, T-WS2, H-VS2, and T-VS2 because of their small lattice mismatch (<2%) with H-MoS2, which makes the computations efficient. Graphene (though strictly speaking a semimetal) is also selected because it has been used experimentally as an electrode for H-MoS2 (2630) (though a large supercell with commensurate structure is required for modeling; see fig. S1 for the structure model; commensurate structures are also used for other largely lattice-mismatched junctions, as shown later). Notably, we find that Φ largely follows the trend of the Schottky-Mott limit (the diagonal line), and its value can vary in a wide range, allowing the tuning of Φ by using different 2D metals. Although different computational methods (more specifically, density functionals) yield different S values, they all give S values higher than that of 3D metal–H-MoS2 junctions (see figs. S2 and S3 and related text in the Supplementary Materials), indicating a weaker FLP.

Although the vdW interaction is weak, it can redistribute the charge density at the interface and give rise to an interface dipole at the MSJ. This is shown in Fig. 1C, where the charge density (ρ) change after the formation of the junction is calculated for the example for H-VS2H-MoS2Embedded Image(4)We find an asymmetric charge accumulation/depletion across the interface, indicating the formation of an interface dipole. As discussed above, this interface dipole can shift the electronic levels from their original positions, leading to a deviation from the Schottky-Mott limit (2, 9, 20, 21).

Because FLP is weak for the vdW junction, Φh (or Φe) can be effectively reduced by using high-W (or low-W) metals as electrodes. To identify the promising candidates, we calculate the W of various 2D metals and compare them with the Eip and Eea values of 2D semiconductors, as shown in Fig. 2. These energies show a systematic variation as the cation or anion changes.

Fig. 2 Band alignment between 2D metals and semiconductors.

Left columns show the electron affinity and ionization potential of semiconductors. Right bars show the work function of metals. The phase is labeled in italics. C stands for pristine graphene, C20N is the N-doped graphene with the C/N ratio of 20:1, and similarly for C20B. For comparison, the work functions of some commonly used 3D metals are also shown.

For 2D semiconductors, we find that (i) as the atomic number of X increases (from S to Se to Te), the CBM and VBM energies increase (except that the CBM of WTe2 is lower than that of WSe2) and band gap decreases; and (ii) moreover, for the common X system, the CBM and the VBM energies of WX2 are higher than those of MoX2. These observations agree with literature results calculated using different methods (31, 32). Therefore, MoS2 has the lowest CBM and VBM, WTe2 has the highest VBM, and WSe2 has the highest CBM. These results explain why MoS2 is usually n-type in experiments but WSe2 is p-type: compared with MoS2, WSe2 has higher CBM and VBM and therefore has higher Φe but lower Φh for given metals, resulting in easier hole injection but more difficult electron injection (33).

For 2D metals, we find that (i) for the same M, the W of MX2 increases as the X atomic number decreases, and (ii) for the same X, the W of MX2 increases as the M atomic number (in the same group) decreases (except for VX2, which is magnetic, in which case the spin polarization changes the order between VX2 and NbX2). These trends originate from the coupling between the M d orbitals and the X p orbitals. As shown by the projected DOS (Fig. 3A), the states near the FL are contributed by the M d states and X p states, and their spatial distribution suggests that they are antibonding-coupled (Fig. 3B). For common M systems, because M d states have higher energy than X p states (31) (Fig. 3C), the p-d coupling weakens and the repulsion level reduces as the p energy decreases from Te to Se to S. This decreases the energy of the antibonding states, and thus the FL moves down and W increases. For common X systems, the M d energy decreases, and the Coulomb repulsion between the d electrons increases when M moves up in the same group (34); thus, the FL decreases and W increases. We also find that for different phases of MX2, the FL decreases as the stability increases. For example, T′-MoX2 and WX2 are more stable than the corresponding T phases (35) and have lower FLs; the situation is similar for group 5 H-MX2 compared with their T phases. Notably, we find that some of the metallic MX2 (group 5 H-MX2 and TiS2) have higher W than Pt, which has the highest W among elemental metals. Of these, H-NbS2 has the highest W, suggesting that it could be a promising electrode for achieving low-Φh contact with the semiconductor, as shown later.

Fig. 3 Origin of work function variation in 2D metals.

(A) DOS of metallic MX2 (T-MoS2 is used here as an example). FL is set to zero. The black line shows the total DOS, and the other lines are the projected DOS on the orbitals of M and X. (B) Charge density distribution of the states in the electron volt range (−0.025, 0.025) (top and side views). (C) Schematic of coupling between M d states and X p states.

We also explore the potential of using doped graphene as electrodes. B has fewer electrons than C, thus B doping increases the W of graphene; inversely, N has more electrons than C, and its doping decreases the W of graphene. In particular, at a high N-doping concentration (C/N, 20:1), W is as low as that of Sc and below all 2D MX2 metals, implying that it could be a promising electrode for electron injection.

After identifying the promising metal electrodes based on their W, we calculate their Φ with the 2D semiconductors, as shown in Fig. 4 (see fig. S4 for representative band structures and table S1 for a comparison of the results calculated with different methods). Indeed, H-NbS2 can form a low-Φh junction. In particular, Φh is negative when H-NbS2 contacts H-WTe2, MoTe2, WSe2, or MoSe2, indicating a spontaneous electron transfer to H-NbS2 and a hole injection into the semiconductor upon contact. Note that this strong charge transfer shifts the FL and the VBM from their original positions, leading to a deviation from the ideal Schottky-Mott model. Similarly, C20N can form a low-Φe contact. However, as opposed to H-NbS2, which has been experimentally synthesized, an efficient N-doping approach up to high concentrations has yet to be developed.

Fig. 4 SB heights for H-NbS2 and C20N metals, with various 2D semiconductors.

C20N with tellurides are not calculated because of the extremely large supercell required for modeling to reduce the lattice mismatch. The red bars indicate the CBM, the blue bars represent the VBM that are set to zero, and the black dots denote the position of FL in the corresponding MSJ.


Using a 2D metal as an electrode has other benefits. First, its transparency and flexibility are compatible with similar features of the 2D semiconductor channel, and it can be integrated for transparent and flexible electronics (28, 30). Second, a 2D metal has limited electronic DOS, and therefore low quantum capacitance. Thus, when charge is accumulated by applying a dielectric-mediated voltage, its W changes markedly compared with conventional metals, which have high DOS. This unique feature of a 2D metal leads to a gate-tunable W and therefore Φ, as has been observed in experiments (26, 3639). Third, it has been shown that a flat interface between the metal and the semiconductor could facilitate carrier transport (40). The atomically flat interface is difficult to achieve by using conventional metals but can be easily realized at the vdW contact. In addition, using a 2D metal as an electrode allows full encapsulation of the semiconductor, which prevents contamination from the environment (29, 39) or stabilizes the reactive semiconductor such as phospherene (41). Moreover, suppression of MIGS reduces the electron-hole recombination at the interface, leading to a higher energy conversion efficiency for optoelectronic devices.

Although this work focuses on the vdW junction, the band alignment shown in Fig. 2 could also provide guidelines for designing a low-resistance “edge” contact junction, where 2D metals and semiconductors are chemically bonded in the same basal plane (42). In addition, Fig. 2 also offers guidelines for designing photovoltaic cells based on 2D materials; in this case, it is critical to select electrodes with an FL close to the VBM or CBM of the semiconductor to obtain a high open-circuit voltage.

We discuss several practical issues for NbS2, which is predicted to be the most promising electrode. In reality, it is easier to obtain a few-layer material than a monolayer; hence, it is interesting to investigate how the electronic structure changes as a function of thickness [we use 3R stacking for multilayer calculations because it is the most common one for NbS2 (43); note that “3R” only denotes the interlayer stacking pattern, and each layer still has the same structure as the monolayer considered above, that is, H-NbS2]. We find that NbS2 remains metallic independent of the thickness, in agreement with calculations in the literature (44). A recent experimental study shows that, in contrast to bulk NbS2, which is a metal, thin NbS2 nanoflakes demonstrate semiconducting properties; this is attributed to the presence of surface defects and/or the nonstoichiometry of the samples that can trap the electrons and quench the carriers (43). The activation energy is found to be very small (1.4 to 1.5 meV) (43). This suggests that it is possible to recover the metallic nature of thin NbS2 by improving the sample quality. We also find that the difference between the W of a monolayer and that of a few-layer (or bulk) material is small (tens of millielectron volts). Therefore, the low Φh should be expected for multilayer NbS2 as well. Note that Φ would be further lowered by the image charge effect (1) when the carriers approach the MSJ during the device operation. The 2D metals may be subject to oxidation when directly exposed to ambient air. This can be overcome by using encapsulation techniques (for example, seal the device with BN) (29, 39, 41), which have been used to stabilize reactive 2D materials such as black phosphorus (which otherwise would quickly degrade in air) (41).

Notice that a recent experimental study shows the presence of interface states at the bilayer-monolayer lateral heterojunction (45), which would trap/scatter the charge carriers transporting along the basal plane. However, these would not affect the carrier injection, which occurs at the bilayer region away from the bilayer-monolayer interface. Finally, we notice a very recent experimental study demonstrating that a low contact resistance for hole transport is achieved by using Nb-doped 2D WSe2 or MoS2 as contact material (46). This further supports our prediction that NbS2 would be a promising candidate.

In summary, a large and untunable SB is overcome by replacing conventional metals with 2D metals, which form a vdW MSJ with 2D semiconductors. This type of MSJ exhibits a weak FLP, which is attributed to the suppression of MIGS. This work not only advances the understanding of fundamental properties of heterojunctions but also shows novel functions of 2D metals in device applications.


Figures in the main text were obtained by using density functional theory with projector-augmented wave pseudopotentials (47, 48) and the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional (49), as implemented in the Vienna ab initio simulation package (50, 51). We also used a local-density approximation functional (52), as well as PBE with vdW corrections (D3) (53), in certain cases to check the robustness of the conclusions, as shown in the Supplementary Materials. Monkhorst-Pack sampling (54) of sufficient k-points was used. The vacuum distance between periodic images was kept to >15 Å. All of the structures were relaxed until the final force on each atom was <0.01 eV/Å. Eea, Eip, and W were calculated by subtracting the corresponding electronic levels with respect to the vacuum level in the supercell.


Supplementary material for this article is available at

fig. S1. Supercell used for modeling the graphene–H-MoS2 junction.

fig. S2. Comparison of pinning factor between vdW and chemically bonded junctions.

fig. S3. Φ between 2D H-MoS2 and various 2D metals calculated by using the PBE-D3 method.

fig. S4. Band structure of H-NbS2H-WSe2 and C20N–H-MoS2 junctions.

table S1. Comparison of Φ between the PBE and the PBE-D3 methods.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Y. L. acknowledges discussions with B. Yakobson and B. Wood on NbS2. Funding: This research was funded by the U.S. Department of Energy (DOE) under contract no. DE-AC36-08GO28308 with the National Renewable Energy Laboratory (NREL) and supported by the DOE Energy Efficiency and Renewable Energy (EERE) DE-261 EE0006336 FPACE-II. This work used computational resources at the NREL (sponsored by the DOE EERE) and the National Energy Research Scientific Computing Center (supported by DOE Office of Science DE-AC02-05CH11231). Author contributions: Y. L. conceived the idea and performed the calculations. Y. L., P. S., and S.-H.W. analyzed the data. Y. L. and S.-H.W. wrote the paper. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Navigate This Article