Ultra-thin van der Waals crystals as semiconductor quantum wells.

Control over the quantization of electrons in quantum wells is at the heart of the functioning of modern advanced electronics; high electron mobility transistors, semiconductor and Capasso terahertz lasers, and many others. However, this avenue has not been explored in the case of 2D materials. Here we apply this concept to van der Waals heterostructures using the thickness of exfoliated crystals to control the quantum well dimensions in few-layer semiconductor InSe. This approach realizes precise control over the energy of the subbands and their uniformity guarantees extremely high quality electronic transport in these systems. Using tunnelling and light emitting devices, we reveal the full subband structure by studying resonance features in the tunnelling current, photoabsorption and light emission spectra. In the future, these systems could enable development of elementary blocks for atomically thin infrared and THz light sources based on intersubband optical transitions in few-layer van der Waals materials.

Introduction

Van der Waals crystals provide a platform for layer-by-layer material design which has expanded into a large multidisciplinary field during the last decade[1-8]. A broad variety of electronic and optoelectronic applications have been described and proof-of-concept devices have been reported[2,5]. Along with high quantum yield[3], compact design, and mechanical flexibility[5], one of the attractive features of two dimensional semiconductors (2DS) is their band-gap tunability. Due to the change in quantum confinement with thickness, the primary optical transition can be tuned over a broad energy range (>1 eV[4,6,8]), and further fine adjustment is available through changes in dielectric environment[7] and transverse electric fields[9].

While the majority of 2DS demonstrate primary optical transitions in visible and near-infrared (IR) ranges, their intersubband transitions hold a largely unexplored potential to expand their optical activity further into the IR and THz. Depending upon the choice of material, number of layers and doping (n- or p-type), these subband transitions densely populate a region from eV for bilayers to 0.05 eV for crystal thicknesess of 10 atomic layers and can be utilized for optoelectronic devices[10,11]. Nonetheless, experimental studies of higher energy subbands remained elusive until recently, when the first experimental observation of intersubband transitions in few-layer was demonstrated[12] using near-field THz absorption spectroscopy. While THz spectroscopy gives access to direct measurements of some intersubband transitions, it only works in a small energy window and requires specific doping of the materials studied. At the same time, angle-resolved photoemission spectroscopy (ARPES) only shows filled (valence) subbands and has diminished energy resolution when applied to small (micrometer-size) crystals[13].

Here, we present an experimental approach that enables a comprehensive experimental study of the full subband structure of atomically thin 2DS, namely resonant tunneling spectroscopy combined with photoluminescence excitation (PLE) measurements. When applied to InSe, this method enables us to map the subbands on both conduction and valence band sides of the spectrum in agreement with theoretical band structure modeling and available ARPES data, as well as to trace the analogy between exfoliated atomically thin crystals and quantum wells in conventional semiconductors[14].

Results and discussion

Experimental devices

Indium selenide is a new material in the 2DS family. It has a layer-dependent bandgap spanning from 1.2 eV for bulk to almost 3.0 eV for a monolayer[6,15,16] and features a crossover[13,17,18] from weakly indirect (in mono-, bi-, and tri-layer) to direct character (in crystals thicker than 4 layers). Its high crystal quality has recently been demonstrated, resulting in exceptionally high electron mobility[6] making it a promising avenue for developing atomically thin nanoelectronics[19]. For resonant tunneling measurements we build van der Waals heterostructures, where exfoliated few-layer -InSe is sandwiched between hexagonal boron nitride (hBN) crystals to provide tunneling barriers, and we employ graphene layers as source and drain contacts on the top and the bottom of the device as illustrated in Fig. 1a. The thickness of hBN was selected to be 4-3 layers to allow for large bias voltages without sample overheating. These structures have been produced using the dry-stacking method in an argon environment[20] to provide atomically clean and sharp interfaces, and to avoid InSe degradation[21,22]. The assembled stacks were deposited on an oxidized silicon wafer and contacts were defined using electron beam lithography (see Supplementary Figs. 1 and 2).

Fig. 1

InSe tunneling devices.

a Cross-sectional annular dark-field scanning transmission electron image of a typical LED device with matching layer schematic on the right. b Optical micrograph of a device used to study resonant tunneling through InSe subbands. Graphene is outlined in black, 2 layer thick (2L) InSe in blue, 4L in green and 5L in red. c Artificially colored electroluminescence (EL) map of the device shown in b: blue corresponding to EL around 1.98 eV, green 1.50 eV and red 1.35 eV with sampling bandwidth of 1 meV (detailed EL spectra can be found in SI). The scale bars are .

An optical micrograph of a typical device consisting of three InSe regions with number of layers 2, 4, and 5 shown in Fig. 1b. Each of these terraces has separate graphene contacts allowing them to have bias applied individually and control the position of the chemical potentials and in the source and drain graphene layers, respectively. To understand vertical charge transport in this system, we have reconstructed its band diagram using a tight-binding-calculated InSe band-structure[6,18] and the relative band alignment measured in our recent ARPES studies[13] and presented in Supplementary Fig. 4 using InSe exfoliated from the same bulk ingot. Here, the graphene Dirac point was near the NL-InSe (N layer-InSe) conduction band edge. The resulting band diagram is shown in Fig. 2 for a 4L-LED device in an unbiased (a), weakly biased (b), strongly biased in forward (c), and strongly biased in reverse configuration (d) 4L-InSe LED device (for the full band structure of InSe, see Supplementary Fig. 8).

Fig. 2

Resonant tunneling spectroscopy of few-layer InSe crystals at T = 4.2K.

Band alignment in 4L InSe device, for unbiased (a), weakly biased (b) and strongly biased in forward (c) and reversed (d) configurations. (blue) and density of states (DoS) of the conduction bands from tight binding model (purple, arb. u.) for 4L (e) and 5L (f) InSe film, with the number of layers established by AFM topography. Energy scale along the top axis was found using the EL onset values. Calculated evolution of the conduction subbands with bias for 4L (g) and 5L (h) InSe (higher energy bands are also shown). Symbols indicate the bias at which the steps were measured in (e) and (f) for forward (brown) and reverse (orange) directions. Evolution of PL and EL with bias for 4L (i) and 5L (j) InSe. The orange dotted line indicates electrostatic calculations of the band gap reduction with electric field, non-dispersive line near 1.5 eV is due to hydrocarbon contamination located outside the sample.

Resonant tunneling spectroscopy

All these anticipated vertical transport regimes can be recognized in the experimentally measured tunneling current derivative () curves shown in Fig. 2e, f (for curves see Supplementary Fig. 3). As soon as a small bias is applied, shows a sharp increase (for 4L-InSe (e) and 5L-InSe (f)), indicating that crosses the edge of the first InSe conduction subband, in agreement with the band diagram in Fig. 2b. Upon further increase in , displays a series of steps. We attribute this behavior to the step-like increase in the density of states of InSe when crosses the conduction subbands . Eventually, of the drain graphene electrode reaches the top subband () of the InSe valence band Fig. 2c, so that electrons and holes injected into the InSe from opposite sides can recombine, producing a bright electroluminescence (EL) signal, shown in the intensity map Fig. 1c[23]. At this point the tunneling current density reaches , and becomes unstable due to the contribution from the recombination process.

Using bias values at the onset of EL, , (indicated as gray shaded areas in Fig. 2e, f) we can determine the rate of change of the InSe potential with bias, as , where is the offset between the Dirac point in graphene and the InSe valence band edge at zero bias[13]. The EL signal is observed to have a slight asymmetry for forward and reverse bias voltages due to the different hBN barrier thicknesses (4L and 3L) on either side: e.g. the thresholds of EL for 4L InSe LED are 2.83 V and . Therefore, the relation between and energy in the band structure of the crystal is slightly different for forward and reverse bias (top axis in Fig. 2e, f). Using the experimentally determined values of the parameter for both devices, we are able to compare the experimentally measured with the theoretically calculated density of electron states (dashed purple line). We note that the linear density of states (DoS) in graphene makes the assignment of observed steps unambiguous, as its only feature, the Dirac point, is located near zero bias where the tunneling current is not measurable. This enables us to determine the energies of subband edges experimentally.

For completeness, in the theoretical analysis we took into account the shifts of both conduction and valence subbands caused by the out-of-plane electric field due to the applied bias as,where and are the respective total thicknesses of the InSe and hBN in the stack, and is the ratio between the out-of-plane dielectric constants of InSe and hBN. This leads to a clearly pronounced quadratic red shift of the lowest subbands energy seen in Fig. 2g, h, similar to the previously observed quantum confined Stark effect in monolayers[24]. At the same time, subbands in the middle of the spectrum ( for 4L and for 5L) remain almost constant, whereas the highest subband is blue shifted. The DoS calculations shown in Fig. 2e, f take into account the bias dependence of the band structure and the resulting subband onsets match closely (within 40 meV) with the step-like features observed in our experiment (shown as gray dashed lines at the highest slope).

The calculated subband shifts for both conduction and valence bands agree well with the shift in the photoluminescence (PL) energy shown in Fig. 2i, j. The observed evolution of EL with bias also accurately follows the trend predicted for PL and at 3 V is shifted by almost 200 meV since the and bands move in opposite directions. It is also clearly non-symmetric, with an almost double EL linewidth for the negative biases, which can be explained with partial filling of the subband when the hBN barrier of the drain is thicker than that of the source. This also agrees with a much faster increase in EL intensity for negative , indicating that the recombination rate is primarily determined by a steeper rise in the DoS of the flatter valence band.

To illustrate the behavior of the subbands in the context of few-layer InSe as a quantum well, we show in Fig. 3d the dispersion of the bulk crystal together with the subband energies for 5-layer InSe at the in-plane point. The 5-layer subband energies can be qualitatively matched with the bulk bands, quantized at 5 discrete momenta as is familiar in the analysis of subband spectra in semiconductor quantum wells[25]. We further illustrate the quantum well behavior of the subbands in Fig. 3c by plotting conduction band wavefunction coefficients for orbitals on the (antisymmetric) indium pair in each layer, showing how the subband with the lowest effective has no nodes, with the number of nodes increasing with .

Fig. 3

PLE of 5L and 7L InSe at T = 4.2K.

a, b Left - 2D band structure of 5L and 7L InSe around computed using DFT-parametrized tight-binding dispersion with spin-orbit interaction taken into account. Right - PL (red), PLE measured at the black arrow (blue) and projected optical DoS (purple) of lamellae for 5 and 7L InSe. Relevant optical transitions from to are labeled as i, j. Orange dotted line shows ARPES intensity around -point[13], plotted to match - optical transition in PLE with the first valence subband in ARPES. Inset shows scanning electron micrograph of the lamella overlaid with its PL map sampled at 1.30(2) eV, scale bar is . c Amplitude , of In -orbital contribution to the subbands across the layers for . Band edges in bulk InSe are at Z points of the 3D Brillouin zone, hence the number of nodes in the layer-dependent microscopic wavefunction counts down from 4 for the lowest energy subband to 0 for the highest. d Solid lines - out-of-plane dispersion of bulk InSe at in-plane point. Orange dots - -point subbands for 5L InSe at given by quantization of momentum in quantum well picture (dashed lines). Insets: center - 3D Brillouin zone of InSe, bottom left - amplitude , of Se -orbital contribution to top valence subband. e Exciton binding energy () extracted from PLE (see inset) as a function of number of layers in the film. The error bars correspond to the width of continuum absorption step. Colors indicate different samples. The dashed curve is displayed as a guide to the eye.

Probing the valence subbands by PLE

Finally, we probed the subband states of InSe films using PLE, where the energy of the excitation light was continuously varied up to 2.7 eV while the photoluminescence intensity of the lowest energy exciton was recorded. The resulting spectra shown in blue Fig. 3a, b for and 7 layers with the corresponding detection energy marked by the black arrow on the PL curve (red). Here, InSe films were encapsulated in thick () hBN crystals without graphene contacts. The optical excitation from the lowest conduction band to the highest valence subband is mostly coupled to light polarized perpendicular to the plane of 2D crystal due to the intralayer electric dipole moment[18]. To provide both in- and out-of-plane polarization of photons for the excitation, we have extracted a cross-sectional block (lamella), with an hBN/InSe/hBN area of  ×  and a depth of using a focused ion beam instrument (see Supplementary Figs. 6 and 7) and rotated it such that the viewing direction is along the basal plane of the InSe crystal (for PLE spectra with in-plane polarized detection of N = 2 to N = 7 see Supplementary Fig. 5).

To interpret the PLE data we note that the material studied was lightly p-doped, so that emission, used for detection involved equilibrium (rather than photoexcited) holes in the subband (dominated by the P orbitals of Se[18]). Due to the p-doping, recombination is promoted, making PLE particularity sensitive to electron excitations () to the subband. As to transitions, the photo-excited electrons have to undergo intersubband relaxation () before recombining radiatively with an equilibrium hole, suppressing their visibility in PLE due to the competition with non-radiative processes. The theoretically calculated projected optical DoS is plotted in Fig. 3a, b (purple dashed lines). It accounts for absorption of both out- and in-plane polarized photons: including contributions from spin-flip interband transitions generated by spin-orbit mixing of P orbitals of Se with the P orbital (marked in green in Fig. 3a, inset). While such mixing is weak for the states in the highest subband , it is enhanced for the deeper subbands, which are closer to the P band on the energy scale. The photon energy thresholds for each transition , determined from the spectra of standing waves subbands in the film, Fig. 3c, are indicated (as “n,0") in Fig. 3a, b, and coincides with the features in the observed PLE spectra. For comparison, on the same plot we show micro-ARPES spectra measured on 5L and 7L films of InSe, which reflect the subband structure of the valence band in these films with noticeably broader features.

Also, on all samples studied, detailed inspection of the low energy onset of the PLE spectra have shown a noticeable excitonic peak followed by a plateau due to the free particle absorption. Measuring the difference between the peak energy (red line in the inset of Fig. 3e) and the onset of continuum absorption (black line)[26], we estimate the layer-dependent exciton binding energies as shown in Fig. 3e. The sharp increase of the binding energy for thin samples agrees with those typically seen in other 2DS[27], while for thick flakes () it closely matches (within 1 meV) that of bulk InSe[28] indicated by dashed horizontal line in Fig. 3e.

Discussion

In summary, by the combined resonant tunneling and photoluminescence excitation spectroscopy approach we have achieved a comprehensive understanding of the entire subband spectrum of few-layer atomically thin films of InSe, including its evolution due to the quantum confined Stark effect. Having been exemplified here for InSe, this approach can be extended to a variety of 2D systems, including hetero- and homo-bilayers of transition metal dichalcogenides with a view to detect features of moiré superlattice minibands[29]. For optoelectronic applications, electrical injection of carriers in vertical devices is a possible route towards creating infrared and terahertz light sources that employ intersubband transitions. The above results demonstrate that, given appropriate doping to provide electrons or holes for recombination, one can control carrier injection into any designated subband of an atomically thin crystal to enable such emitters.

Methods

Device fabrication

For all devices Bridgman-grown bulk rhombohedral -InSe[13] crystals were mechanically exfoliated onto PPC-coated (polypropylene carbonate) silicon wafers and crystals displaying large uniform terraces were identified using optical microscopy. These flakes were picked up with graphene/hBN stack carried on a PMMA membrane[20] followed by a second thin hBN crystal and a graphene monolayer. Finally, the assembled stack was released onto a thick (50 nm) hBN exfoliated on an oxidized silicon wafer. The exfoliation and transfer of thin InSe took place in an argon filled glovebox to protect it from degradation[19]. Both graphene crystals were then contacted using electron beam lithography followed by CHF etching and Cr/Au (1.5 nm/50 nm) deposition. Using an additional e-beam lithography step, the top graphene layer was divided into separate electrodes following outlines of the InSe terraces using Ar/O plasma etching.

For the measurements of PL and PLE, devices without graphene have been made to prevent the emission quenching[30,31]. As before, InSe was encapsulated between hBN crystals and then covered with an additional thick hBN layer (>200 nm) followed by Cr/Au (1 nm/80 nm) to prevent Ga ion damage of the InSe during the milling process. Using FEI Helios dual-beam focused ion beam (FIB)/scanning electron microscope (SEM), an additional layer of platinum () was deposited over the selected region to provide further protection. Ga+ FIB milling was used to define the lamella at 30 kV, with current 7 nA to remove bulk of the material and current 1 nA to trim the lamella to its correct size. An OmniProbe micromanipulator was used to extract the lamella, rotate it by , and position it on an OMICRON transmission electron microscopy (TEM) grid. After milling, the damaged edges of the 2D stack were removed by FIB polishing using decreasing acceleration voltages (5 kV, 47 pA and 2 kV, 24 pA). The final thickness of the specimen was .

Optical measurements

All optical measurements have been performed at 4 K using an AttoDry100 cryostat and a AttoCFM I inset from Attocube. The illumination power on the sample was kept ~0.4 mW and the spot size was . For the excitation, a supercontinuum white light laser (Fianium WhiteLase WL-SC-400-15-PP) combined with a Contrast filter (LLTF SR-VIS-HP8) have been used resulting in accurate wavelength control between 400 nm and 1000 nm with a linewidth of  nm. For PLE measurements, the laser power was stabilized before the sample using an RD40-UV laser power controller from Brockton Electro-Optics Corp. The PL/PLE spectra were recorder using spectrometer diffraction grating (Princeton Instruments Acton Spectrapro SP-2500i with 300 g/mm) and nitrogen-cooled CCD camera (Princeton Instruments Pylon PIX-100BRX).

Tunneling spectroscopy

The tunneling measurements were carried out at 4.2 K by applying a DC bias and measuring the current between the source and drain graphene electrodes using a Keithley 2614B source-meter.

Tight-binding model

Details concerning the spin-orbit coupling, the calculation of the density of states and the absorption spectra can be found in the Supplementary Fig. 9.

Cross-sectional STEM

For high resolution scanning transmission electron microscope (STEM) imaging the FIB cross section received further thinning with a Helios Dual beam FIB-SEM, using sequentially lower milling currents of 30 kV, 16 kV, 5 kV, and 2 kV for ion beam milling and polishing. A probe side aberration corrected FEI Titan G2 80-200 kV was used with a probe convergence angle of 21 mrad, a HAADF inner angle of 48 mrad and a probe current of  pA. To ensure the electron probe was parallel to the basal planes, the cross-sectional FIB sample was aligned to the relevant Kikuchi bands of the Si substrate and the 2D crystals.