## Tunable Electronic Structure in Gallium Chalcogenide van der Waals Compounds

Brian Shevitski, 1, 2, ∗ Søren Ulstrup, 3, 4, † Roland J Koch, 3 Hui Cai, 5 Sefaattin Tongay, 5 Luca Moreschini, 3 Chris Jozwiak, 3 Aaron Bostwick, 3 Alex Zettl, 1 Eli Rotenberg, 3 and Shaul Aloni 2, ‡

1 Department of Physics, University of California at Berkeley, Berkeley, California 94720,USA 2 Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720,USA 3 Advanced Light Source, E. O. Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 4 Department of Physics and Astronomy, Interdisciplinary Nanoscience Center (iNANO), Aarhus University, 8000 Aarhus C, Denmark 5 School for Engineering of Matter, Transport and Energy, Arizona State University, Tempe, Arizona 85287, USA (Dated: September 23, 2019)

Transition metal monochalcogenides comprise a class of two-dimensional materials with electronic band gaps that are highly sensitive to material thickness and chemical composition. Here, we explore the tunability of the electronic excitation spectrum in GaSe using angle-resolved photoemission spectroscopy. The electronic structure of the material is modified by in-situ potassium deposition as well as by forming GaS x Se 1 -x alloy compounds. We find that potassium-dosed samples exhibit a substantial change of the dispersion around the valence band maximum (VBM). The observed band dispersion resembles that of a single tetralayer and is consistent with a transition from the direct gap character of the bulk to the indirect gap character expected for monolayer GaSe. Upon alloying with sulfur, we observe a phase transition from AB to AA ′ stacking. Alloying also results in a rigid energy shift of the VBM towards higher binding energies which correlates with a blue shift in the luminescence. The increase of the band gap upon sulfur alloying does not appear to change the dispersion or character of the VBM appreciably, implying that it is possible to engineer the gap of these materials while maintaining their salient electronic properties.

## I. INTRODUCTION

The ability to isolate monolayers and study their physical properties has led to a resurgence of interest in layered van der Waals (VDW) materials as two-dimensional systems. Transition-metal chalcogenides belong to this class of materials and demonstrate unique thermoelectric, photonic and electronic properties 1 . The composition, phase, and crystal structure of these materials can be tuned to display a wide range of electronic phases including metallic 2,3 , semiconducting 4-6 , superconducting 7-9 and charge-density wave 9-11 , inspiring utilization in a wide range of applications including electronic devices, sensors, and quantum devices 12 .

In their bulk form, monochalcogenides of Ga and In, such as GaSe, InSe, and GaS are van der Waals compounds. Gallium chalcogenides (GCs) are direct bandgap semiconductors that can be mechanically exfoliated from a bulk parent crystal. Their properties are highly sensitive to material thickness, chemical composition, and crystalline structure 13-15 . Bulk GaSe and GaS have direct band gaps of 2.10 eV and 3.05 eV, respectively, and transition to indirect gaps when the number of layers is reduced 16,17 . Monolayer field effect devices maintain mobilities on the order of of 0.1 cm V 2 -1 s -1 , on/off ratios of 10 4 -10 5 as well as quantum efficiencies that exceed those of graphene by several orders of magnitude when implemented in a device architecture 18-20 .

Alloys with mixed chalcogen content, GaS x Se 1 -x , allow for control of the bandgap, allowing engineering of the optical absorption and photoluminescence (PL) 21,22 . Interlayer coupling, defects, doping, and interaction with an underlying substrate may further modify the electronic structure of the GCs, as observed in exfoliated flakes of GaSe and GaTe on SiO2 and graphene 23-25 and epitaxially grown GaSe flakes on graphene 17,26,27 .

In this communication, we illustrate how the crystalline and electronic properties of GC materials are affected by changes in chemical composition, doping, and reduced dimensionality using high resolution scanning transmission electron microscopy (HRSTEM), angleresolved photoemission spectroscopy (ARPES), and PL.

HRSTEM reveals that GaSe tetralayers stack in the AB sequence ( ε phase) while sulfur alloys and GaS stack almost exclusively in the AA ′ sequence ( β phase). The complete E k ,k ,k ( x y z ) electronic structure of AB-stacked GaSe is determined using ARPES, allowing us to directly extract the band dispersion along high symmetry directions and determine valence band (VB) extrema and effective masses. We illustrate how the full energymomentum-resolved quasiparticle band structure of the GCs is modulated when the chemical composition and/or interaction between layers is changed. We explore the influence of potassium deposition on the band structure, finding it only causes a slight energy shift of the VBs and does not lead to degenerate doping which would enable ARPES visualization of the conduction band. However, potassium dosing does lead to a significant change of the dispersion around the top of the VB, implying a strong modification at the surface GaSe tetralayer compared to the bulk. Finally, we explore the effect of sulfur alloying on the band structure. We find that increasing sulfur content produces a rigid VB shift and results in a corresponding blue shift in the PL, suggesting that alloying is

FIG. 1. (color online) Real space structure and HRSTEM images of the two most common GaS x Se 1 -x polytypes. Panels (a) and (b) show the crystal structure of the ε phase in the planes defined by the axes in the labels. The corresponding HRSTEM image in (c) shows a trigonal lattice, resulting from the projection of Bernal (AB) stacked honeycomb layers. The intensity cut in (d) was obtained from the red rectangle in (c) and shows bright atomic columns with approximately twice the intensity of dim atomic columns. Analogous images in (e)-(f) show the geometry of the β phase. The HRSTEM image in (g) shows a honeycomb mesh, with no intensity in the interstitial columns between hexagonal cells.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01003v2-figure-4fe47a07-6f51-461b-9d76-61029e8a96d4.jpg)

a viable method for tuning GC optical properties.

## II. EXPERIMENTAL

We study pure GaSe and GaS as well as three alloy compositions spanning the entire compositional range of GaS Se x 1 -x alloys. Crystals are grown using the modified Bridgman-Stockbarger method 28 . Briefly, stoichiometric amounts of Ga, S, and Se powder are weighed and mixed in sealed quartz ampoules to achieve the desired alloying ratio. The Ampoules are then heated in a zone furnace at 970 ◦ C for two weeks to grow single crystals 2-8 mm in diameter.

The HRSTEM samples are prepared by a two-step process. First, the crystals are cleaved using blue wafer dicing tape and transferred to Si substrates coated with 100 nm thick SiO . 2 Next, a Cu TEM grid with carbon film is adhered to flakes on the Si/SiO 2 wafer using a drop of isopropanol. Once the drop has dried, the SiO 2 is etched away using 1M NaOH. Lastly the grid is washed in deionized water and dried prior to imaging. HRSTEM is carried out using a double C s corrected FEI Titan 80300 operating at an accelerating voltage of 80 kV at the Molecular Foundry.

as defined by spatially scanning the sample with the synchrotron beam. The ARPES data are collected at the Microscopic and Electronic Structure Observatory (MAESTRO) at the Advanced Light Source (ALS) using the microARPES end-station equipped with a hemispherical Scienta R4000 analyzer. The beamline slit settings are adjusted so that the size of the beam is on the order of 20 µ m. VB and Ga 3d core level spectra are obtained primarily using a photon energy of 94 eV. S 2p and Se 3p core level data are collected using a photon energy of 300 eV. Photon energy scans of GaSe are acquired for photon energies between 21 eV and 140 eV. In order to relate the photoelectron kinetic energy, E kin , to the out-ofplane momentum, k z , we use the free-electron final state assumption where k 2 z = ( 2 m/h ¯ 2 ) (E kin +V ), where V 0 0 is the inner potential 29 . We find that V 0 = 10 2 eV . provides the best description of the data, given the outof-plane lattice constant c = 15.96 ˚ and the Brillouin A zone (BZ) periodicity of 2 π/ c for the measured polytypes of GaSe.

For ARPES measurements, large crystals of GCs are glued to a Cu sample holder using epoxy and cleaved in-situ in the ultra-high vacuum (UHV) chamber with a base pressure better than 5 × 10 -11 mbar. The samples are cooled to 85 K prior to cleaving and kept at this temperature during ARPES measurements. The cleaved single-crystal domain sizes are on the order of 100 µ m

The Fermi energy is determined on the clean Cu sample holder in contact with the crystals for all photon energies. The presence of surface photovoltage-induced energy shifts in the samples is checked by varying the photon flux. No such changes were observed, except on GaS, which excluded ARPES measurements on this material.

Potassium dosing experiments on GaSe are carried out in -situ using SAES getters mounted in the analysis chamber while the sample is kept at 85 K. Each dose takes 50 seconds, and between doses, an ( E,k ) spectrum of the VB and an energy spectrum over the Ga 3d/K 3p binding energy region are taken. The collection of these

FIG. 2. (color online) ARPES measurements at a photon energy of 94 eV of the electronic structure of pristine GaSe around the Γ point. (a) Sketch of the BZs and their orientation in the experiment with Γ, K and M symmetry points labeled.(b) Measured dispersion in the Γ -K -M direction obtained along the blue dashed line shown in (a). (c) Measured dispersion in the Γ -M -Γ direction obtained along the red dashed line shown in (a). The double headed arrows in (b)-(c) indicate the planes of constant energy from which the cuts in (d)-(h) are taken. The energy of each cut is also listed in bottom of each panel in units of electron volts (eV). High symmetry points Γ, K and M have been labeled in (e), (f) and (g), respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01003v2-figure-2b347d51-f724-4280-922c-e951bd4cb41d.jpg)

spectra takes 70 seconds and 12 doses were applied for the results presented here, amounting to a total time of 24 minutes for collecting the dataset.

Judging from the change in photoemission intensity from the VB upon further dosing, we estimate a coverage rate of a complete monolayer after 12 doses. We observe that prolonged exposure to the synchrotron beam for the fully K-dosed sample leads to broadening and significant deterioration of the quality of the spectral features, preventing detailed photon energy or angle scans for the K-dosed samples.

The total energy and momentum resolution in the ARPES data are better than 20 meV and 0.01 ˚ A -1 , respectively.

## III. RESULTS AND DISCUSSION

## A. HRSTEM Measurements of Gallium Chalcogenides

The GaS Se x 1 -x system forms a continuous set of alloys for the entire compositional range (0 ≤ x ≤ 1) of the stoichiometry 30-32 . Like other layered VdW materials, the intralayer bonding is strong and mostly covalent while the interlayer bonding is weak and of VdW character. The unit cell consists of two tetrahedral pyramids stacked tip to tip with chalcogen atoms at the base and a Ga atom at the apex, forming a four-atom-thick tetralayer.

The layer stacking exhibits several different polytypes, depending on composition and growth conditions. In the ε phase (Fig. 1(a-b)), successive tetralayers are stacked in the sequence AB, resulting in Ga atoms from the A layer being vertically aligned with chalcogen atoms in the the B layer, and chalcogen atoms in the A layer being vertically aligned with the interstitial voids in the centers of the hexagons from the B layer, similar to Bernal stack-

ing in graphite. The inversion symmetric β phase (Fig. 1(e-f)) is obtained by stacking successive tetralayers in the sequence AA , resulting in Ga(chalcogen) atoms from ′ the A layer eclipsing chalcogen(Ga) atoms from the A ′ layer. The system can also form the rhombohedral γ phase (ABC stacking), but this stacking sequence is not observed experimentally in our samples.

Previous studies have explored the relationship between sulfur content and stacking phase 30-32 . GaS Se x 1 -x alloys with low sulfur content (x &lt; 0.15) stack in the ε phase while alloys with higher sulfur content (x &gt; 0.35) display the β stacking phase. Our experimental results (described below) also show this trend. We find that the ε phase occurs for x ≤ 0.3 and the β phase for x &gt; 0.3, The β and ε phases are thought to co-exist for intermediate sulfur content (0.15 &lt; x &lt; 0.35), phase segregated throughout the crystal, but our samples do not show experimental evidence for this theory.

Contrast in HRSTEM indicates differences in atomic number, density, and thickness throughout a specimen. Due to the specific stacking sequences exhibited by the polytypes of GaS x Se 1 -x alloys and the contrast mechanisms of STEM imaging, our imaging experiments (Fig. 1(c) and 1(g)) unambiguously reveal the stacking phases of these materials.

Figure 1(c) shows a typical HRSTEM image for an alloy with high sulfur content (x ∼ 0 75), . appearing as a trigonal lattice with alternating bright and dim atomic columns. Figure 1(d) shows a line plot of the intensity extracted from the region outlined by the red rectangle in Fig. 1(c). The heights of the peaks in the line trace indicate that there are approximately twice as many atoms in bright columns compared to the dim columns, consistent with the stacking geometry of the ε phase, shown in Fig. 1(b). Figure 1(g) shows an analogous image for pure GaSe which appears as a honeycomb mesh of atomic columns with no appreciable intensity in the interstitial voids between hexagons. The line profile (Fig. 1(h)) shows that the atomic columns have uniform intensity, consistent with AA ′ stacking exhibited by the β phase.

## B. ARPES Electronic Structure Measurements of GaSe

An overview of the electronic band structure of pristine GaSe is presented in Fig. 2. The data are acquired at a photon energy of 94 eV, which probes around the Γ point, as described by the BZ sketches in Fig. 2(a). Figures 2(b)-(c) show cuts along the Γ -K -M and Γ -M -Γ high symmetry directions as signified by the dashed lines in Fig. 2(a). The evolution of the dispersion with binding energy is presented via the constant binding energy cuts in Fig. 2(d)-(h), which reveal the trigonal symmetry of the dispersion as well as strong intensity variations between the bands in the BZs. The valence band maximum (VBM) is situated at Γ and is characterized by a parabolic-shaped lobe. The location of the VBM at Γ

FIG. 3. (color online) ARPES measurements of the k z -dispersion of GaSe. (a) Sketch of bulk BZ highlighting the Γ -A -H -K plane. (b) Dispersion along the Γ -A direction, as marked by the blue dashed line in (a). Fits of the position of the VBM intensity are shown via open circles. The solid curve following the fitted dispersion is a guide for the eye. The BZ size of 2 π/ c in the z -direction is marked by a double-headed arrow and two vertical dashed lines. (c)-(d) Dispersion extracted at the given k z values, corresponding to the Γ -K direction (red dashed line). (e) Sketch of 1st and 2nd BZs with an outline of the Γ -A -L -M plane. (f) Constant energy contour of the plane sketched in (e) extracted at a binding energy of 0.2 eV with high symmetry points annotated by tick marks. (g) Dispersion along the Γ -M -Γ direction obtained at a k z value of 4.75 ˚ A -1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01003v2-figure-5cd01ef8-cb83-48bc-86ab-40842bca7522.jpg)

is consistent with theoretical predictions that GaSe is a direct-gap semiconductor with the conduction band minimum (CBM) also located at Γ 33 . Additional sub bands dispersing towards M and K appear above binding energies of 1.5 eV. The material remains rather featureless away from Γ before the onset of this higher binding energy range, implying most of the optical properties of GaSe will derive from the central parabolic lobe. The ARPES data discussed here are fully consistent with the electronic structure data presented in Ref. 34.

We investigate the k z -dispersion in the three dimensional BZ of GaSe shown in Fig. 3(a). By measuring the photon energy dependence of the dispersion at normal emission we are able to trace the band structure along the Γ -A direction, as seen in Fig. 3(b). A strong k z -dispersion is visible in the center of the parabolic lobe, which is further highlighted via the fits (open circles) and guide line in Fig. 3(b), whereas the bands at binding energies higher than 1.5 eV have a much more flat k z -

dispersion. The fitted dispersion reveals an apparent BZ doubling effect; however, this occurs due to a suppression of the photoemission intensity in adjacent BZs 34 , similar to observations in photoemission experiments on graphite 35,36 . The band structure along the Γ -K line is shown for two neighboring BZs in the k z direction in Fig. 3(c-d). The main intensity moves from the bottom corners of the parabolic lobe in Fig. 3(c) to the VBM in Fig. 3(d). Similar intensity variations are observed in neighboring in-plane BZs, as sketched in Fig. 3(e). The ( k || , k z )-contour at a binding energy of 0.2 eV shown in Fig. 3(f) shows a cut around the VBM for the Γ -A -L -M plane (see Fig. 3(e)). The intensity maxima are seen to shift between the 1st and 2nd BZs and the absolute intensity levels also vary strongly within the same BZ. A cut along the Γ -M line at k z = 4.75 ˚ A -1 in Fig. 3(g) shows the upper VB parabolic lobe in two neighboring zones. In the first BZ, the intensity is concentrated in the bottom corner while in the 2nd BZ the intensity is concentrated around the VBM.

This behavior can be reconciled with the orbital character of the bands, previously studied with detailed calculations 13 . Since the states near the VBM are mainly of Se p z character, the interlayer interactions will cause intensity modulations reflecting the coupling of the GaSe tetralayers, leading to a dispersion resembling a single tetralayer for the k z value used for the slice shown in Fig. 3(c), while the k z value for Fig. 3(d) shows a dispersion with a more bulk-like character 13 . The subbands just below the VBM exhibit less dramatic behavior due to the p x y ( ) character of these states which are less sensitive to the adjacent GaSe tetralayers 37,38 .

## C. Potassium Deposition of GaSe

The effect of K dosing of bulk GaSe is explored in Fig. 4. Alkali metal deposition offers a possibility to chemically donate electrons near the surface of GaSe and thereby dope the material. The ARPES measurements along the Γ -K -M direction in Fig. 4(a-b) before and after complete K dosing show that a significant band structure change occurs, rather than strong electron doping. Amaximum energy shift of 0.1 eV of the VB towards higher binding energies is observed, which may originate from surface band bending due to the added charge or a slight doping of the material. After further dosing, we observe that the spectral weight shifts from the top (see Fig. 4(a)) to the bottom (see Fig. 4(b)) of the VB lobe. In Fig. 4(c), the VBM dispersion (open circles), extracted via fits to energy distribution curves (EDCs), has been fitted with 2nd and 4th order polynomial functions in the clean and K doped cases, respectively. In the pristine GaSe, a global VBM is observed at Γ and at a binding energy of (0.09 ± 0.03) eV. The effective mass around the top of the parabola is m ∗ = (1.1 ± 0.2) m 0 , where m 0 is the free electron mass. In K dosed GaSe, the VBM splits off into two maxima located at ± 0.3 ˚ A -1

FIG. 4. (color online) Effect of K dosing on the GaSe electronic structure. (a)-(b) ARPES Measurements of the dispersion in the Γ -K -M direction (a) before and (b) after K dosing. The dashed curves present the fitted location of the VBM peak intensity. (c) Location of the VBM peak intensity (open circles) fitted with a 2nd order and a 4th order polynomial for the clean and K dosed samples, respectively. (d) EDCs (open circles) and results of Voigt function fits (curves) in the clean and K dosed cases. The EDCs were extracted at Γ along the lines shown in (a)-(b). The position of the top-most peak is plotted via dashed lines in (a)-(b) and open circles in (c). (e) Evolution of the photoemission intensity at Γ as a function of K dosing steps. The open circles and guidelines mark the fitted peak positions obtained by performing the EDC analysis shown in (d) in each K dosing step. (f)-(g) Ga 3d core level measurements (f) before and (g) after K dosing. The data (thick red curves) have been fitted (dashed curves) with Doniach-Sunjic line shapes (filled peaks) for each of the observed components.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01003v2-figure-4cd32891-4af3-4516-94fc-df658bf2e475.jpg)

with a binding energy of (1.30 ± 0.04) eV and an effective mass of m ∗ = (1.7 ± 0.2) m 0 .

The EDCs extracted at normal emission in Fig. 4(d) trace the peak positions at each K dosing step. Voigt line shapes have been fitted to each of the peaks, which have been plotted together with the photoemission intensity at normal emission at each K dose in Fig. 4(e), revealing a gradual shift of the ARPES intensity from the pristine bulk dispersion seen in Fig. 4(a) to the modified dispersion shown in Fig. 4(b). The modified dispersion strongly resembles the 'bow-shape' or'inverted sombrero' dispersion of single-layer GaSe, where the inversion of the VB is characterized by the energy difference between the VBM and the local energy minimum at Γ 13,22,39-41 . Here, we find that the band inversion is (48 ± 12) meV, which compares well with a value of 80 meV obtained from recent DFT calculations 40 . The value is substantially smaller than 120 meV found for single-layer GaSe grown on graphene on silicon carbide 41

and than the value of (150 +/- 10) meV that we determine for the similar dispersion at k z = 3 15 ˚ . A -1 in Fig. 3(c). We expect that interlayer interactions between the Se p z orbitals play a significant role for these dispersion changes around the VBM 13 and that the gradual change observed in Fig. 4(e) signals a modification of the topmost GaSe layer, leading to a weaker coupling to the underlying bulk. This may be facilitated by K intercalation in the van der Waals gap between the layers, as observed in metallic 42 , and semiconducting 43 transition metal dichalcogenides, or by a chemical interaction on the surface of the crystal, which significantly changes the coupling between individual layers as seen in graphene multilayers 44,45 . In these situations one may also expect a change of the interlayer separation of tetralayers as well as a shift of inner potential, which also contribute to the spectral changes we observe. The possibility of chemical changes near the surface is further supported by the observation of an additional core level component in the binding energy range between the Ga 3d and K 3p core levels for the K-dosed sample, which is demonstrated by the comparison between the clean GaSe core level spectrum in Fig. 4(f) and the spectrum for the fully K-dosed sample in Fig. 4(g).

FIG. 5. (color online) Composition of GaS x Se 1 -x alloys from core level data. (a)-(c) Core level measurements (red thick curves) in the Se 3p and S 2p binding energy region for (a) GaSe and (b)-(c) GaS x Se 1 -x alloys. Fit results (black curves) to Doniach-Sunjic line shapes (filled peaks) are provided in each panel. The composition x of the alloys determined via the fitted peak components is stated in (b)-(c). Note the intensity scale varies between the panels.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01003v2-figure-9c48b4cf-ec2d-42d4-b2a8-df00ee12f0b8.jpg)

## D. ARPES and PL Measurements of GaS Se x 1 -x

GaS Se x 1 -x alloys offer a promising route to tune the optical gap and luminescent properties of III-VI semiconductors. However, it is unknown how robust the electronic structure remains between different alloys. Here we compare pristine GaSe with two different GaS x Se 1 -x alloys (alloys I and II). The stoichiometry x is checked via core level measurements of the Se 3p and S 2p binding energy regions, shown in Fig. 5. By comparing the spectral weight of the Se 3p core levels, we are able to estimate that the composition of alloy I is x = 0.27 ± 0.05 and the composition of alloy II is x = 0.61 ± 0.05. These values are consistent when we perform the same analysis on the Se 3d core levels (not shown). The core level data is in agreement with the stoichiometry of precursor added to the growth ampoule.

Figure 6(a)-(c) present ARPES measurements of the VBs for the three systems. The overall band dispersion is found to be very similar with a constant effective mass of m ∗ = (1.1 ± 0.2) m 0 . The binding energy position of the VBM changes significantly between the three systems, indicating a change of the electronic band gap. It is also possible that Fermi level pinning plays a strong

FIG. 6. (color online) VBM and PL peaks positions for GaS Se x 1 -x alloys. (a)-(c) ARPES measurements of the upper VB region for (a) GaSe, (b) GaS 0 27 . Se 0 73 . and (c) GaS 0 61 . Se 0 39 . . The dotted line presents the fitted position of the VBM and the dashed curves exhibit the fitted dispersion around the VBM in each case. The corresponding VBM binding energy positions are given in units of eV with an error bar of ± 0 03 eV. (d) PL intensity for the same samples as in . (a)-(c) and Fig. 5. (e)-(f) Positions of the (e) PL peaks and (f) VBM as a function of composition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01003v2-figure-eebbb039-0226-464c-b567-91a66b561f42.jpg)

role for the band positions, as the defect concentration is likely different in the alloys. PL measurements, shown in Fig. 6(d), reveal a change of the optical gap. The PL peak positions plotted in Fig. 6(e) show that the optical gap increases as more S is added in the alloy. The same trend is seen for the VBM binding energy position in Fig. 6(f), which fits with a simultaneous increase of the quasiparticle band gap. The trend determined here is fully consistent with optical absorption measurements of GaS Se x 1 -x 22 . Note that ARPES measurements of pure GaS were attempted but the sample was found to charge too severely, which is likely due to the larger band gap of the material compared to the other GaS x Se 1 -x compounds investigated here.

The ε and β stacking phases observed in HRSTEM are not expected to give rise to striking differences in the electronic structure 37-39 . However, detailed inspection of the measurements in Fig.. 6(a)-(c) does show some variation in the intensity and structure of the dispersive bands in the higher binding energy region above the central VB lobe. While different stacking phases may affect the photoemission intensity, the change in chemical environment caused by the mixed chalcogen content directly influences the in-plane p x y ( ) orbitals that the higher binding energy bands are composed of, thereby leading to the possibility of a modified dispersion.

## IV. CONCLUSION

We have measured the bulk electronic structure of GaSe using ARPES and thereby identified the bulk band dispersion, the structure of the VBM, and the complex photoemission intensity variations associated with the tetralayer unit cell. Potassium deposition of GaSe leads to a dramatic change of the VBM dispersion from a parabolic shape characteristic for the bulk to an inverted bow-like shape that is usually associated with single-layer samples. Our data are consistent with an increased band

- ∗ Authors contributed equally.
- † ulstrup@phys.au.dk; Authors contributed equally.
- ‡ saloni@lbl.gov
- 1 M. Wasala, H. I. Sirikumara, Y. Raj Sapkota, S. Hofer, D. Mazumdar, T. Jayasekera, and S. Talapatra, J. Mater. Chem. C 5 , 11214 (2017).
- 2 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306 , 666 (2004).
- 3 J. Feng, X. Sun, C. Wu, L. Peng, C. Lin, S. Hu, J. Yang, and Y. Xie, Journal of the American Chemical Society 133 , 17832 (2011).
- 4 K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, Phys. Rev. Lett. 105 , 136805 (2010).
- 5 A. Splendiani, L. Sun, Y. Zhang, T. Li, J. Kim, C.-Y. Chim, G. Galli, and F. Wang, Nano Lett. 10 , 1271 (2010).

gap and a transition from a direct to an indirect band gap semiconductor, which is caused by a strong modification of the top-most GaSe tetralayer due to the interaction with the deposited K atoms. We have investigated the effect of sulfur alloying on the crystalline and electronic structure and used HRSTEM to observe a transition from the ε phase to the β phase for sulfur content above ∼ 30 percent. We have also shown that alloying causes a rigid shift of the VBM binding energy position without a significant change in the actual band dispersion around the VBM. This rigid shift is consistent with an increase in the optical band gap measured by PL.

## ACKNOWLEDGMENTS

This work was supported primarily by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, under Contract No. DE-AC02-05CH11231, within the van der Waals Heterostructures Program (KCWF16), which provided for HRSTEM and ARPES support. This work was also supported in part by the National Science Foundation under Grant No. DMR-1807233, which provided for PL measurements and preliminary exploration of sample growth parameters. S.U. acknowledges financial support from VILLUM FONDEN under the Young Investigator Program (Grant No. 15375). R.J.K. was supported by a fellowship within the Postdoc-Program of the German Academic Exchange Service (DAAD). The Advanced Light Source is supported by the Director, Office of Science, Office of Basic Energy Sciences, of the U.S. Department of Energy under Contract No. DEAC02-05CH11231. Work at the Molecular Foundry was supported by the Office of Science, Office of Basic Energy Sciences, of the US Department of Energy under Contract No. DE-AC02-05CH11231. B.S. acknowledges support from the NSF LSAMP BD graduate fellowship.

- 6 K. Watanabe, T. Taniguchi, and H. Kanda, Nat. Mater. 3 , 404 (2004).
- 7 I. Guillam´ on, H. Suderow, S. Vieira, L. Cario, P. Diener, and P. Rodi` ere, Phys. Rev. Lett. 101 , 166407 (2008).
- 8 E. Navarro-Moratalla, J. O. Island, S. Ma˜ nas-Valero, E. Pinilla-Cienfuegos, A. Castellanos-Gomez, J. Quereda, G. Rubio-Bollinger, L. Chirolli, J. A. Silva-Guill´n, e N. Agra¨ ıt, G. A. Steele, F. Guinea, H. S. J. van der Zant, and E. Coronado, Nat. Commun. 7 , 11043 (2016).
- 9 M. M. Ugeda, A. J. Bradley, Y. Zhang, S. Onishi, Y. Chen, W. Ruan, C. Ojeda-Aristizabal, H. Ryu, M. T. Edmonds, H.-Z. Tsai, A. Riss, S.-K. Mo, D. Lee, A. Zettl, Z. Hussain, Z.-X. Shen, and M. F. Crommie, Nat. Phys. 12 , 92 (2015).
- 10 G. Liu, S. Rumyantsev, M. A. Bloodgood, T. T. Salguero, and A. A. Balandin, Nano Lett. (2018).
- 11 X. Xi, L. Zhao, Z. Wang, H. Berger, L. Forr´, J. Shan, o and K. F. Mak, Nat. Nanotechnol. 10 , 765 (2015).

- 12 W. Choi, N. Choudhary, G. H. Han, J. Park, D. Akinwande, and Y. H. Lee, Materials Today 20 , 116 (2017).
- 13 D. V. Rybkovskiy, A. V. Osadchy,
- and E. D. Obraztsova, Phys. Rev. B Condens. Matter 90 , 235302 (2014).
- 14 H. Cai, E. Soignard, C. Ataca, B. Chen, C. Ko, T. Aoki, A. Pant, X. Meng, S. Yang, J. Grossman, F. D. Ogletree, and S. Tongay, Advanced Materials 28 (2016).
- 15 J. Robertson, J. Phys. C: Solid State Phys. 12 (2001).
- 16 S. Lei, L. Ge, Z. Liu, S. Najmaei, G. Shi, G. You, J. Lou, R. Vajtai, and P. M. Ajayan, Nano Letters 13 , 2777 (2013).
- 17 Z. Ben Aziza, D. Pierucci, H. Henck, M. G. Silly, C. David, M. Yoon, F. Sirotti, K. Xiao, M. Eddrief, J.-C. Girard, and A. Ouerghi, Phys. Rev. B 96 , 035407 (2017).
- 18 P. Hu, Z. Wen, L. Wang, P. Tan, and K. Xiao, ACS Nano 6 , 5988 (2012).
- 19 Y. Zhou, Y. Nie, Y. Liu, K. Yan, J. Hong, C. Jin, Y. Zhou, J. Yin, Z. Liu, and H. Peng, ACS Nano 8 , 1485 (2014).
- 20 D. J. Late, B. Liu, J. Luo, A. Yan, H. S. S. R. Matte, M. Grayson, C. N. R. Rao, and V. P. Dravid, Advanced Materials 24 , 3549 (2012).
- 21 H. Serizawa, Y. Sasaki, and Y. Nishina, Journal of the Physical Society of Japan 48 , 490 (1980).
- 22 C. S. Jung, F. Shojaei, K. Park, J. Y. Oh, H. S. Im, D. M. Jang, J. Park, and H. S. Kang, ACS Nano (2015), 10.1021/acsnano.5b04876.
- 23 H. Chen, Y. Li, L. Huang, and J. Li, Rsc Adv 5 , 50883 (2015).
- 24 O. Del Pozo-Zamudio, S. Schwarz, M. Sich, I. A. Akimov, M. Bayer, R. C. Schofield, E. A. Chekhovich, B. J. Robinson, N. D. Kay, O. V. Kolosov, A. I. Dmitriev, G. V. Lashkarev, D. N. Borisenko, N. N. Kolesnikov, and A. I. Tartakovskii, 2D Materials 2 (2015).
- 25 W. Kim, C. Li, F. A. Chaves, D. Jimnez, R. D. Rodriguez, J. Susoma, M. A. Fenner, H. Lipsanen, and J. Riikonen, Advanced Materials 28 , 1845 (2016).
- 26 Z. Ben Aziza, H. Henck, D. Pierucci, M. G. Silly, E. Lhuillier, G. Patriarche, F. Sirotti, M. Eddrief, and A. Ouerghi, ACS Nano (2016).
- 27 X. Li, L. Basile, B. Huang, C. Ma, J. Lee, I. Vlassiouk, A. Puretzky, M. Lin, M. Yoon, M. Chi, J. Idrobo, C. Rouleau, B. Sumpter, D. Geohegan, and K. Xiao, Acs Nano 9 , 8078 (2015).
- 28 Y. Ni, H. Wu, C. Huang, M. Mao, Z. Wang, and X. Cheng, Journal of Crystal Growth 381 , 10 (2013).
- 29 S. H¨ ufner, Photoelectron Spectroscopy , 3rd ed. (Springer, Berlin, 2003).
- 30 C. Whitehouse and A. Balchin, Journal of Materials Science 13 , 2394 (1978).
- 31 J. Terhell, V. Brabers, and G. van Egmond, Journal of Solid State Chemistry 41 , 97 (1982).
- 32 G. Arancia, M. Grandolfo, C. Manfredotti, and A. Rizzo, physica status solidi (a) 33 , 563 (1976).
- 33 D. Olgu´ ın, A. Rubio-Ponce, and A. Cantarero, Eur. Phys. J. B 86 , 350 (2013).
- 34 L. Plucinski, R. Johnson, B. Kowalski, K. Kopalko, B. Orlowski, Z. Kovalyuk, and G. Lashkarev, Physical Review B 68 , 125304 (2003).
- 35 E. L. Shirley, L. J. Terminello, A. Santoni, and F. J. Himpsel, Physical Review B 51 , 13614 (1995).
- 36 M. Mucha-Kruczy´ski, O. Tsyplyatyev, A. Grishin, E. Mcn Cann, V. I. Fal'ko, A. Bostwick, and E. Rotenberg, Physical Review B 77 , 195403 (2008).
- 37 M. Camara, A. Mauger, and I. Devos, Physical Review B 65 , 125206 (2002).
- 38 S. Nagel, A. Baldereschi, and K. Maschke, Journal of Physics C: Solid State Physics 12 , 1625 (1979).
- 39 D. Rybkovskiy, N. Arutyunyan, A. Orekhov, I. Gromchenko, I. Vorobiev, A. Osadchy, E. Y. Salaev, T. Baykara, K. Allakhverdiev, and E. Obraztsova, Physical Review B 84 , 085314 (2011).
- 40 Z. Ben Aziza, V. Zlyomi, H. Henck, D. Pierucci, M. G. Silly, J. Avila, S. J. Magorrian, J. Chaste, C. Chen, M. Yoon, K. Xiao, F. Sirotti, M. C. Asensio, E. Lhuillier, M. Eddrief, V. I. Fal'ko, and A. Ouerghi, Physical Review B 98 , 115405 (2018).
- 41 M.-W. Chen, H. Kim, D. Ovchinnikov, A. Kuc, T. Heine, O. Renault, and A. Kis, npj 2D Materials and Applications , 1 (2018).
- 42 K. Rossnagel, E. Rotenberg, H. Koh, N. V. Smith, and L. Kipp, Physical Review Letters 95 , 126403 (2005).
- 43 T. Eknapakul, P. D. C. King, M. Asakawa, P. Buaphet, R.H. He, S.-K. Mo, H. Takagi, K. M. Shen, F. Baumberger, T. Sasagawa, S. Jungthawan, and W. Meevasana, Nano Letters 14 , 1312 (2014).
- 44 T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Rotenberg, Science 313 , 951 (2006).
- 45 T. Ohta, A. Bostwick, J. L. McChesney, T. Seyller, K. Horn, and E. Rotenberg, Physical Review Letters 98 , 206802 (2007).
Brian Shevitski Søren Ulstrup Roland J. Koch Hui Cai Sefaattin Tongay Luca Moreschini Chris Jozwiak Aaron Bostwick Alex Zettl Eli Rotenberg Shaul Aloni
Tunable Electronic Structure in Gallium Chalcogenide van der Waals Compounds
Transition metal monochalcogenides comprise a class of two-dimensional
materials with electronic band gaps that are highly sensitive to material
thickness and chemical composition. Here, we explore the tunability of the
electronic excitation spectrum in GaSe using angle-resolved photoemission
spectroscopy. The electronic structure of the material is modified by
$\textit{in-situ}$ potassium deposition as well as by forming
GaS$_{x}$Se$_{1-x}$ alloy compounds. We find that potassium decouples the
top-most tetra-layer of the GaSe unit cell, leading to a substantial change of
the dispersion around the valence band maximum (VBM). The observed band
dispersion of a single tetralayer is consistent with a transition from the
direct gap character of the bulk to the indirect gap character expected for
monolayer GaSe. Upon alloying with sulfur, we observe a phase transition from
AB to $\text{AA}^{\prime}$ stacking. Alloying also results in a rigid energy
shift of the VBM towards higher binding energies which correlates with a blue
shift in the luminescence. The increase of the band gap upon sulfur alloying
does not appear to change the dispersion or character of the VBM appreciably,
implying that it is possible to engineer the gap of these materials while
maintaining their salient electronic properties.