We report a lowering of work function for lanthanum hexaboride (LaB6) by monolayer hexagonal boron nitride (hBN) coating. Photoemission electron microcopy (PEEM) and thermionic emission electron microscopy (TEEM) both revealed that the hBN coated region of a LaB6 (100) single crystal has a lower work function compared to the bare (i.e., non-coated) and graphene coated regions. A broad and uniform brighter image of the hBN coated region in PEEM was quantitatively supported by a 0.4 eV decrease in the work function in photoelectron spectra compared to the bare region. TEEM results were consistent in that the hBN coated region exhibited thermionic emission at 905 °C, whereas the bare and graphene coated regions did not. A larger decrease in the work function for hBN coated LaB6 (100) compared to graphene coated LaB6 (100) was qualitatively supported by our density functional theory calculations. Adding an oxide layer in the calculations improved consistency between the calculation and experimental results. We followed up our calculations with synchrotron-radiation x-ray photoelectron spectroscopy and confirmed the presence of an oxide layer on our LaB6.

Free electrons are vital components for a wide range of systems from photodetectors and microscopes to accelerators. Multiple cathode types are capable of generating such free electrons (photo, thermionic, and field emission), and an intrinsic material property that determines their performance is the work function or the minimum energy required for electrons to escape from a material surface into a vacuum. For photocathodes, Spicer's three-step model1 identifies the work function as the main property that determines the escape rate of electrons into a vacuum when a material is excited with photons. Therefore, the work function is essential, because it determines the quantum efficiency (QE) of the photocathodes, where QE is a rate of free electrons generated per incident photons. Electrons overcome the work function via thermal energy to generate free electrons in thermionic cathodes. For field emission cathodes, free electrons are generated by electrons tunneling through the potential barrier of cathode surfaces via an electric field, where the barrier height is the work function. The challenge in improving the cathode performance lies in the preparation of robust low work function materials that allow cathodes to generate more free electrons at a given illumination wavelength, temperature, or electric field, at which the cathodes operate.

Here, we investigated the effects of a two-dimensional (2D) material coating on lanthanum hexaboride (LaB6) for photo- and thermionic-cathodes. 2D materials, such as hexagonal boron nitride (hBN) and graphene, are thermally stable up to at least 2200 °C, are chemically inert, and possess an excellent gas barrier property.2–7 In addition, there are experimental and theoretical reports on lowering the work function of metals and alkali antimonide semiconductors by using 2D material coating.8–12 Thus, if coating LaB6 with 2D materials lowers its work function, this could provide a breakthrough in achieving both robustness and a higher yield of free electron generation. In general, lowering the work function of a thermionic cathode by coating is challenging, because the coating materials need to withstand a high temperature of ∼1500 °C and should not increase the work function. This combination is extremely rare. For photocathodes, lowering of the work function can readily be achieved by well-established cesium (Cs) termination as demonstrated on gallium arsenide (GaAs) and other materials.13 The problem is that those surfaces are highly sensitive to environmental conditions; high chemical reactivity of Cs requires an ultrahigh vacuum (UHV) of 10−11 Torr for its operation, which severely limits operation environments. Even for the cathodes with Cs incorporation, in contrast to termination, such as the case of cesium antimonide (Cs3Sb) and cesium potassium antimonide (CsK2Sb), their operations still require 10−10 Torr. An exception is Cs2Te, which could operate in one or two orders of magnitude higher pressure; however, it requires 266 nm excitation due to its higher work function in contrast to 532 nm for other Cs terminated and incorporated photocathodes.

We chose the LaB6 cathode for our purpose, because it is a well-established thermionic cathode material that is commercially available for electron microscopes and spectrometers. Its clean surface has a low work function of ∼2.4 eV (Ref. 14) and a long lifetime of ∼5000 h due to a low evaporation rate at its operation temperature of ∼1500 °C. The low work function of LaB6 is also attractive as a photocathode—although the level of interest in this aspect is low compared to thermionic cathodes,15–20 interest in using LaB6 as a photocathode exists none-the-less due to its high temporal flexibility of generated free electrons compared to other types of cathodes.

We first coated the single crystal LaB6 (100) (AP-Tech) with hBN and graphene via an established wet-transfer method. Poly(methyl methacrylate) (PMMA) polymer support was coated onto hBN on copper (Cu) and graphene on Cu that were synthesized via chemical vapor deposition (CVD).21–26 Cu foils were then etched away using ammonium persulfate (APS) followed by transfers of PMMA supported hBN and graphene onto LaB6 (100) by exposure to water.27 As shown in Fig. 1(a), half of LaB6 (100) was first covered with graphene and dried in air. PMMA was then removed by immersing the entire assembly in an acetone bath for three cycles of 15 min each. Subsequently, a quarter of bare and graphene covered LaB6 (100) was covered by hBN, dried in air, and then PMMA was removed by the same protocol as that of graphene. The LaB6 (100) surface was flat with a roughness of below 1 nm by atomic force microscopy (AFM) [Fig. 1(b)]. We also performed AFM on transferred hBN and graphene and confirmed that they are continuous sub-nm films, thus, atomically thin [Figs. 1(c) and 1(d)]. Raman spectroscopy on hBN and graphene indicated that their crystal quality is high [Figs. 1(e) and 1(f)]. Specifically, the Raman spectrum of hBN had a signature E2g peak around 1370 cm−1,28 and that of graphene had signature G and 2D peaks around 1550 and 2700 cm−1 with no observable defect D peak around 1350 cm−1.29 

FIG. 1.

(a) Top view photograph of graphene and hBN coated LaB6 (100). LaB6 is fixed onto a nickel (Ni) substrate using silver paint. Regions of 2D material coatings are indicated by green (graphene) and red (hBN) dotted lines. (b) Atomic force microscopy (AFM) image of a bare LaB6 (100) surface. (c) AFM image of hBN wet-transferred onto a SiO2/silicon (Si) substrate. White arrows indicate a region of void to show a presence of the film. (d) AFM image of graphene wet-transferred onto a SiO2/Si substrate. The maximum on the height scale bar is 1 and 10 nm for (b) and (c)/(d), respectively. Lateral scale bars for (b)–(d) are 1 μm. (e) Raman spectrum of hBN wet-transferred onto the SiO2/Si substrate. The red arrow indicates the position of the E2g signature peak of hBN. (d) Raman spectrum of graphene wet-transferred onto the SiO2/Si substrate. The peak positions of G and 2D are labeled.

FIG. 1.

(a) Top view photograph of graphene and hBN coated LaB6 (100). LaB6 is fixed onto a nickel (Ni) substrate using silver paint. Regions of 2D material coatings are indicated by green (graphene) and red (hBN) dotted lines. (b) Atomic force microscopy (AFM) image of a bare LaB6 (100) surface. (c) AFM image of hBN wet-transferred onto a SiO2/silicon (Si) substrate. White arrows indicate a region of void to show a presence of the film. (d) AFM image of graphene wet-transferred onto a SiO2/Si substrate. The maximum on the height scale bar is 1 and 10 nm for (b) and (c)/(d), respectively. Lateral scale bars for (b)–(d) are 1 μm. (e) Raman spectrum of hBN wet-transferred onto the SiO2/Si substrate. The red arrow indicates the position of the E2g signature peak of hBN. (d) Raman spectrum of graphene wet-transferred onto the SiO2/Si substrate. The peak positions of G and 2D are labeled.

Close modal

The sample was then introduced into a system that can perform photoemission electron microcopy (PEEM), photoelectron spectroscopy (PES), and thermionic emission electron microscopy (TEEM) (Fig. 2). The sample temperature was gradually increased from room temperature up to 905 °C. TEEM was performed at elevated temperatures, and PEEM and PES were performed at room temperature after each annealing temperature to exclude thermal effects. In PEEM, ∼4.9 eV photons generated by a mercury (Hg) lamp were illuminated on the sample surface and photoelectrons with a spatial resolution of ∼20 nm to gain information on the spatial distribution of the work function. We followed our PEEM measurements with PES to quantify the difference in the work function for the regions of interest. For PES, the energy resolution was <0.2 eV, and the sample was biased to −3 V to efficiently collect photoelectrons. The Fermi level of the samples was determined using the Fermi edge of an evaporated gold (Au) thin film. We first used the helium (He) I line with an excitation energy of 21.22 eV to determine the Fermi edge of an evaporated Au film. We then used the Hg lamp on the same evaporated Au film and used the energy difference between He I and Hg lamp excitation to determine the work function of the samples. The purpose of performing TEEM was to correlate the work function observed in PEEM and PES with thermionic emission. We also performed synchrotron-radiation x-ray photoelectron spectroscopy (SR-XPS) separately to investigate the chemical state of our LaB6 (100) surfaces. Aside from experiments, we performed a density functional theory (DFT) calculation to theoretically investigate the effect of 2D material coating on the work function of LaB6 (100) surfaces.

FIG. 2.

Schematic illustration of the PEEM/TEEM/PES measurement setup. An Hg lamp with ∼4.9 eV is used as an illumination source.

FIG. 2.

Schematic illustration of the PEEM/TEEM/PES measurement setup. An Hg lamp with ∼4.9 eV is used as an illumination source.

Close modal

Figure 3(a) is a schematic of 2D material coverage on the LaB6 (100) surface, which indicates graphene and hBN in red and green dashed lines, respectively. Figure 3(b) is a PEEM image of a sample at room temperature after being heated at 905 °C, which was the maximum temperature of our system. Interestingly, we observed a higher intensity of photoelectrons at the region corresponding to the hBN coating, which reveals that the work function is lower for the hBN coated LaB6 (100) surface compared to the non-coated and graphene-coated counterparts. Note that the hBN region is brighter than the graphene region, even where there is a graphene coating underneath [i.e., in between hBN and LaB6 (100)] confirming that hBN is playing a critical role in lowering the work function. To demonstrate that our concept is also applicable to thermionic cathodes, we performed TEEM. The regions of 2D material coverage on the LaB6 (100) surface are indicated in the same way as that of PEEM. 905 °C is the temperature that we started to observe thermionic emission as shown in Fig. 3(c). The results were consistent with PEEM and PES in that the brightest region aligned well with the hBN coated region, and the graphene coated region had brightness in between the hBN-coated region and the non-coated bare regions. The bright star shaped spot at the bottom left of the TEEM image is an artificial effect of the measurements. No observable signal was obtained for TEEM at 806 °C, indicating that thermionic emission at the temperature was below our detection limit [Fig. 3(d)].

FIG. 3.

(a) Top view schematic illustration of the LaB6 geometry. White, green, and red regions are bare LaB6, graphene covered LaB6, and hBN covered LaB6, respectively. (b) PEEM image taken at room temperature after heating up to 905 °C. (c) TEEM image taken at 905 °C. (d) TEEM image taken at 806 °C. Scale bars for (b)–(d) are 50 μm. Green and red dotted lines in (b)–(d) correspond to those in (a).

FIG. 3.

(a) Top view schematic illustration of the LaB6 geometry. White, green, and red regions are bare LaB6, graphene covered LaB6, and hBN covered LaB6, respectively. (b) PEEM image taken at room temperature after heating up to 905 °C. (c) TEEM image taken at 905 °C. (d) TEEM image taken at 806 °C. Scale bars for (b)–(d) are 50 μm. Green and red dotted lines in (b)–(d) correspond to those in (a).

Close modal

We followed up PEEM measurements with PES to quantify the brightness difference. As you can see from Fig. 4, the trend we observed in PEEM measurements was consistent with PES. The work function for non-coated and hBN coated regions was 4.2 and 3.8 eV, respectively. This demonstrates that as large as ∼0.4 eV of the work function can be lowered by simply coating LaB6 (100) with 2D materials such as hBN. hBN coated and hBN/graphene coated regions had similar work function; thus, it was difficult to distinguish them by the brightness in the PEEM image [Fig. 3(b)]. PES showed that the hBN/graphene coated region had a work function that is comparable to that of the hBN coated region (3.9 eV). The work function for the graphene coated region was 4.1 eV, slightly lower than the non-coated, bare surface. Our PEEM and PES results demonstrated that the work function lowering of LaB6 (100) that we achieved is directly applicable to photocathodes because both techniques generate free electrons via photoexcitation.

FIG. 4.

PES spectra of LaB6 coated with graphene and hBN taken at room temperature after heating up to 905 °C. The red arrow indicates a shift of low-energy cutoff of the spectra for hBN coated LaB6 compared to bare and graphene coated LaB6.

FIG. 4.

PES spectra of LaB6 coated with graphene and hBN taken at room temperature after heating up to 905 °C. The red arrow indicates a shift of low-energy cutoff of the spectra for hBN coated LaB6 compared to bare and graphene coated LaB6.

Close modal

The work function of 4.2 eV we observed for non-coated LaB6 (100) was much higher than ∼2.4 eV reported for its clean surfaces. This is probably due to the presence of an oxide layer.14,30 It is known that LaB6 surfaces, including the (100) form native oxide layer in air, and annealing at ∼1500 °C under 1 × 10−7 Torr is necessary for its complete removal. Our annealing temperature of 905 °C was well below ∼1500 °C. Moreover, our observed work function of 4.2 eV is close to the reported value of ∼4.0 eV for the LaB6 (100) surface with a saturated oxide layer.14 We performed SR-XPS at the surface chemistry end-station (BL23SU) of SPring-8 in Japan using soft x-rays (711 eV) to gain information on the oxidation states of our LaB6 (100) surfaces. As expected, we observed lanthanum oxide contribution in the O 1s peak even after annealing at 905 °C (Fig. 5). Entire peak positions and shapes as well as peak positions of deconvoluted peaks were consistent with those reported for La2O3 prior to and after annealing at 800 °C.31,32 Specifically, while there was a dominant broad peak around 531.47 eV prior to annealing, the peak shape changed drastically after the annealing and peaks formed at 531.03 and 531.85 eV. The peak at lower binding energy has been assigned to lattice oxygen. There is ongoing discussion on the exact assignment of the peak at the higher binding energy due to the complex structure, which contains a mixture of both carbonates and hydroxyls. The results suggest that the annealing of our air exposed LaB6 (100) surfaces at 905 °C forms La2O3. We observed the B 1s peak around 187.5 eV for wide scans prior to and after annealing, indicating that the thickness of our oxide layer is below a few nm as expected for a native oxide layer.

FIG. 5.

XPS O 1s spectra of non-coated LaB6 prior to [room temperature (RT)] and after the annealing at 905 °C. The intensity is normalized to peak intensities, and the peaks are deconvoluted.

FIG. 5.

XPS O 1s spectra of non-coated LaB6 prior to [room temperature (RT)] and after the annealing at 905 °C. The intensity is normalized to peak intensities, and the peaks are deconvoluted.

Close modal

We performed density functional theory (DFT) calculations on the work function of LaB6 (100) with 2D material coatings to gain insight into the mechanism of observed work function lowering. Specifically, we calculated the work function of oxidized LaB6 (100) surfaces and these surfaces after coating with hBN and graphene monolayers. We also calculated the case for pristine LaB6 (100) surfaces as a reference. Our calculations were performed with the use of DFT and the projector augmented-wave (PAW) method33 as implemented in the Vienna ab initio simulation package (VASP).34 The generalized gradient approximation (GGA) of the Perdew–Burke–Ernzerhof (PBE)35 functional was used to represent the exchange–correlation interaction. The substrates were represented with slab models with a thickness of ∼12 Å and a vacuum gap in the direction normal to the surface to eliminate the interactions between the replicas due to the periodic boundary conditions. The surfaces of the slabs are terminated with La atoms based on previous experiments. The size of the primitive cell of the LaB6 slab was 4.155 × 4.155 Å2, and the size of a rectangular cell of graphene/hBN was 2.467 × 4.274 Å2. Therefore, a 3 × 1 supercell of the LaB6 slab (12.465 × 4.155 Å2) and a 5 × 1 supercell of graphene/hBN (12.335 × 4.274 A2) had the lattice mismatch of ∼3%.

We started our calculation on a pristine LaB6 (100) surface as a reference to ensure that our calculated value was consistent with experiment values. The calculated work function was 2.20 eV, which was indeed in good agreement with an experimental value of 2.3 eV (Ref. 36) [Fig. 6(a), Table I]. After coating with the hBN monolayer, its work function decreased to 1.87 eV. However, the calculated work function increased to 3.44 eV with graphene coating [Fig. 6(b), Table I]. This modification originates from the interfacial dipoles induced by the coating as discussed later. These calculations qualitatively support the work function modifications of LaB6 that we observed experimentally by the 2D material coatings.

FIG. 6.

Schematics of the atomic arrangement used for DFT calculations. Large green, small green, brown, and red spheres represent lanthanum (La), boron (B), carbon (C), and oxygen (O) atoms, respectively. (a) LaB6 (100), (b) top (left) and side (right) views of graphene coated LaB6 (100), and (c) top (left) and side (right) views of oxidized LaB6 (100).

FIG. 6.

Schematics of the atomic arrangement used for DFT calculations. Large green, small green, brown, and red spheres represent lanthanum (La), boron (B), carbon (C), and oxygen (O) atoms, respectively. (a) LaB6 (100), (b) top (left) and side (right) views of graphene coated LaB6 (100), and (c) top (left) and side (right) views of oxidized LaB6 (100).

Close modal
TABLE I.

Summary of the DFT calculated and experimentally obtained work function of bare and 2D material coated LaB6 (100). The numbers in ()s for oxidized LaB6 (100) indicate the modification induced by 2D material coatings.

Work function (eV) LaB6 (100) hBN\LaB6 (100) Gr\LaB6 (100) Oxidized LaB6 (100) hBN\oxidized LaB6 (100) Gr\oxidized LaB6 (100)
DFT calculations  2.20  1.87  3.44  3.12 (0)  2.99 (−0.13)  3.88 (+0.76) 
Experiments in this study  2.336   NA  NA  4.2 (0)  3.8 (−0.4)  4.1 (−0.1) 
Work function (eV) LaB6 (100) hBN\LaB6 (100) Gr\LaB6 (100) Oxidized LaB6 (100) hBN\oxidized LaB6 (100) Gr\oxidized LaB6 (100)
DFT calculations  2.20  1.87  3.44  3.12 (0)  2.99 (−0.13)  3.88 (+0.76) 
Experiments in this study  2.336   NA  NA  4.2 (0)  3.8 (−0.4)  4.1 (−0.1) 

After confirming that our calculation values for the clean surface were consistent with the experimental values, we moved on to a more representative case for our experiments with an oxide layer. The calculated work function of the oxidized LaB6 (100) surface was dependent on the coverage of oxygen (O) as reported experimentally. For the surface with O atoms at the equivalent metal bridge sites (with O coverage of 1/2 on these available sites30), the calculated work function was 3.60 eV. This is close to ∼4.0 eV reported experimentally for the oxygen saturated surface.14 However, a large supercell was required to match this surface structure with graphene and hBN coatings, which means that the calculations involving 1/2 O coverage are computationally expensive. Therefore, we used an oxidized LaB6 (100) surface with 1/3 O coverage instead in our study to represent an oxidized LaB6 surface [Fig. 6(c), Table I]. Its calculated work function was 3.12 eV, slightly lower than the case of 1/2 O coverage as expected. When the surface was covered with an hBN monolayer, the work function decreased to 2.99 eV. This is qualitatively consistent with the case of clean LaB6 (100), indicating that the calculations predict lowering of the work function by hBN coating on both the clean and oxide surfaces. The work function increased to 3.88 eV after graphene coating, and this was also qualitatively consistent with the case of clean surfaces. Band diagrams illustrating the work function change of oxidized LaB6 (100) by graphene and hBN coating are shown in Fig. 7. For the semi-metallic graphene, the charge transfer from LaB6 to the conduction bands of the coating materials induces inward pointing dipoles, which increase the work function [Figs. 7(a) and 7(b)]. The wave-like curves in Fig. 7(c) are the planar-averaged electron transfer along the direction perpendicular to the surface. The inset of Fig. 7(c) shows the corresponding electron transfers at the interface of oxidized LaB6 (100) surface and graphene. Blue and yellow regions indicate electron depletion and accumulation, respectively. On the other hand, for monolayer hBN, because of its large bandgap, an outward dipole is formed at the interface originating from exchange repulsion, which decreases the work function [Fig. 7(d) and 7(e)]. Figure 7(f) in comparison to Fig. 7(c) clearly shows that the surface dipole has different directions for graphene and BN coated surfaces. For both the hBN and graphene coating, the degree of work function modification on oxidized surfaces was smaller compared to the clean surfaces, indicating that the oxide layer reduces the change transfers between LaB6 and 2D materials.

FIG. 7.

Schematic band alignment of graphene with the oxidized LaB6 (100) surface: (a) before contact and (b) after contact. (c) Charge redistribution and interfacial dipole moment. Schematic band alignment of hBN with the oxidized LaB6 (100) surface: (d) before contact and (e) after contact. (f) Charge redistribution and interfacial dipole moment. Large green, small green, brown, red, and white spheres in insets of (c) and (f) represent lanthanum (La), boron (B), carbon (C), oxygen (O), and hydrogen (H) atoms, respectively. Blue and yellow regions show the electron transfer at the interface of oxidized LaB6 (100) surface and 2D materials.

FIG. 7.

Schematic band alignment of graphene with the oxidized LaB6 (100) surface: (a) before contact and (b) after contact. (c) Charge redistribution and interfacial dipole moment. Schematic band alignment of hBN with the oxidized LaB6 (100) surface: (d) before contact and (e) after contact. (f) Charge redistribution and interfacial dipole moment. Large green, small green, brown, red, and white spheres in insets of (c) and (f) represent lanthanum (La), boron (B), carbon (C), oxygen (O), and hydrogen (H) atoms, respectively. Blue and yellow regions show the electron transfer at the interface of oxidized LaB6 (100) surface and 2D materials.

Close modal

In summary, lowering of the work function for LaB6 by 2D material coating is demonstrated. Specifically, PEEM and TEEM both revealed that an hBN coated region of the LaB6 (100) single crystal has a lower work function compared to bare (i.e., non-coated) and graphene coated regions. A broad and uniform brighter image of an hBN coated region in PEEM was quantitatively supported by a 0.4 eV decrease in the work function in photoelectron spectra compared to the bare region. TEEM results were consistent in that the hBN coated region exhibited thermionic emission at 905 °C, whereas the bare and graphene coated regions did not. A larger decrease in the work function for hBN coated LaB6 (100) compared to graphene coated LaB6 (100) was qualitatively supported by our density functional theory (DFT) calculations. Adding a presence of an oxide layer on LaB6 in the calculations improved consistency between the calculations and experimental results. We followed up our calculations with SR-XPS and confirmed the presence of an oxide layer on our LaB6 (100) single crystal.

This work was supported by the U.S. Department of Energy (DOE) Office of Science U.S.–Japan Science and Technology Cooperation Program in High Energy Physics (HEP) and the JSPS KAKENHI (Grant No. JP17KK0125). The XPS measurements were performed using synchrotron radiation at Beamline BL23SU of SPring-8 (Project Nos. 2021A-E14, 2021B-E22, 2022A-E20, and 2022B-E13) with the approval of the Japan Synchrotron Radiation Research Institute (JASRI) (Proposal Nos. 2021A3831, 2021A3801, 2021B3831, 2021B3801, 2022A3831, 2022A3801, 2022B3831, and 2022B3801). A Part of this work was supported by the JAEA Advanced Characterization Nanotechnology Platform under the remit of “Nanotechnology Platform” of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan (Grant Nos. JPMXP09A21AE0013 and JPMXP09A21AE0042). A part of this work was supported by “Advanced Research Infrastructure for Materials and Nanotechnology in Japan (ARIM)” of the MEXT, Japan (Proposal Nos. JPMXP1222AE0019 and JPMXP1222AE0029). The work was also supported, in part, by the Los Alamos National Laboratory (LANL) Laboratory Directed Research and Development (LDRD) Program through Directed Research (DR) “Cathodes And Rf Interactions In Extremes (CARIE)” (Project No. 20230011DR). Studies were performed, in part, at the Center for Integrated Nanotechnologies, an Office of Science User Facility, operated for the U.S. Department of Energy (DOE), Office of Science. Los Alamos National Laboratory, an affirmative action equal opportunity employer, is managed by Triad National Security, LLC for the U.S. Department of Energy’s NNSA, under contract 89233218CNA000001. This research used resources provided by the Los Alamos National Laboratory Institutional Computing Program, which is supported by the U.S. Department of Energy National Nuclear Security Administration under Contract No. 89233218CNA000001.

N.A.M. owns the patent US 8,823,259 on a concept of graphene protection of photocathodes and the patent US 10,679,816 on a concept of graphene protection of thermionic cathodes. The other authors declare no competing financial interest.

Hisato Yamaguchi: Conceptualization (lead); Funding acquisition (lead); Investigation (lead); Writing – original draft (lead). Shuichi Ogawa: Data curation (lead); Supervision (lead). Ryunosuke Yusa: Data curation (supporting). Gaoxue Wang: Data curation (lead). Michael Thompson Pettes: Data curation (supporting). Fangze Liu: Data curation (supporting). Yasutaka Tsuda: Data curation (supporting). Akitaka Yoshigoe: Resources (supporting); Supervision (supporting). Tadashi Abukawa: Resources (supporting); Supervision (supporting). Nathan A. Moody: Resources (supporting).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
W. E.
Spicer
,
Phys. Rev.
112
(
1
),
114
122
(
1958
).
2.
A. K.
Geim
and
I. V.
Grigorieva
,
Nature
499
(
7459
),
419
425
(
2013
).
3.
K. S.
Novoselov
,
A.
Mishchenko
,
A.
Carvalho
, and
A. H.
Castro Neto
,
Science
353
(
6298
),
aac9439
(
2016
).
4.
H.
Yamaguchi
,
J.
Granstrom
,
W.
Nie
,
H.
Sojoudi
,
T.
Fujita
,
D.
Voiry
,
M.
Chen
,
G.
Gupta
,
A. D.
Mohite
,
S.
Graham
, and
M.
Chhowalla
,
Adv. Energy Mater.
4
(
4
),
1300986
(
2014
).
5.
F.
Liu
,
N. A.
Moody
,
K. L.
Jensen
,
V.
Pavlenko
,
C. W. N.
Villarrubia
,
A. D.
Mohite
, and
G.
Gupta
,
Appl. Phys. Lett.
110
(
4
),
041607
(
2017
).
6.
C. S.
Rout
,
P. D.
Joshi
,
R. V.
Kashid
,
D. S.
Joag
,
M. A.
More
,
A. J.
Simbeck
,
M.
Washington
,
S. K.
Nayak
, and
D. J.
Late
,
Sci. Rep.
3
(
1
),
3282
(
2013
).
7.
U. V.
Patil
,
A. S.
Pawbake
,
L. G. B.
Machuno
,
R. V.
Gelamo
,
S. R.
Jadkar
,
C. S.
Rout
, and
D. J.
Late
,
RSC Adv.
6
(
54
),
48843
48850
(
2016
).
8.
E.
Starodub
,
N. C.
Bartelt
, and
K. F.
McCarty
,
Appl. Phys. Lett.
100
(
18
),
181604
(
2012
).
9.
P. A.
Khomyakov
,
G.
Giovannetti
,
P. C.
Rusu
,
G.
Brocks
,
J.
van den Brink
, and
P. J.
Kelly
,
Phys. Rev. B
79
(
19
),
195425
(
2009
).
10.
G.
Wang
,
P.
Yang
,
N. A.
Moody
, and
E. R.
Batista
,
NPJ 2D Mater. Applicat.
2
(
1
),
17
(
2018
).
11.
G.
Wang
,
P.
Yang
, and
E. R.
Batista
,
Phys. Rev. Mater.
4
(
2
),
024001
(
2020
).
12.
G.
Giovannetti
,
P. A.
Khomyakov
,
G.
Brocks
,
V. M.
Karpan
,
J.
van den Brink
, and
P. J.
Kelly
,
Phys. Rev. Lett.
101
(
2
),
026803
(
2008
).
13.
A. H.
Sommer
,
Photoemissive Materials
(
Robert E. Krieger Publishing Company
,
1980
).
14.
R.
Nishitani
,
C.
Oshima
,
M.
Aono
,
T.
Tanaka
,
S.
Kawai
,
H.
Iwasaki
, and
S.
Nakamura
,
Surf. Sci.
115
(
1
),
48
60
(
1982
).
15.
P. E.
Oettinger
,
Appl. Phys. Lett.
56
(
4
),
333
334
(
1990
).
16.
D. J.
Bamford
,
M. H.
Bakshi
, and
D. A. G.
Deacon
,
Nucl. Instrum. Methods Phys. Res. Sect. A: Accel., Spectrom., Detectors Assoc. Equip.
318
(
1
),
377
380
(
1992
).
17.
S. H.
Kong
,
J.
Kinross-Wright
,
D. C.
Nguyen
, and
R. L.
Sheffield
,
Nucl. Instrum. Methods Phys. Res. Sect. A: Accel., Spectrom., Detectors Assoc. Equip.
358
(
1
),
272
275
(
1995
).
18.
M. E.
Conde
,
S. I.
Kwon
,
A. T.
Young
,
K. N.
Leung
, and
K. J.
Kim
,
Appl. Phys. Lett.
65
(
21
),
2660
2662
(
1994
).
19.
P. G.
O'Shea
,
J. A.
Lancaster
,
R.
Sachtschale
, and
C. R.
Jones
,
Appl. Phys. Lett.
73
(
3
),
411
413
(
1998
).
20.
K.
Torgasin
,
K.
Morita
,
H.
Zen
,
K.
Masuda
,
T.
Katsurayama
,
T.
Murata
,
S.
Suphakul
,
H.
Yamashita
,
T.
Nogi
,
T.
Kii
,
K.
Nagasaki
, and
H.
Ohgaki
,
Phys. Rev. Accel. Beams
20
(
7
),
073401
(
2017
).
21.
X.
Li
,
W.
Cai
,
J.
An
,
S.
Kim
,
J.
Nah
,
D.
Yang
,
R.
Piner
,
A.
Velamakanni
,
I.
Jung
,
E.
Tutuc
,
S. K.
Banerjee
,
L.
Colombo
, and
R. S.
Ruoff
,
Science
324
(
5932
),
1312
1314
(
2009
).
22.
Y.
Uchida
,
K.
Kawahara
,
S.
Fukamachi
, and
H.
Ago
,
ACS Appl. Electron. Mater.
2
(
10
),
3270
3278
(
2020
).
23.
Y.
Uchida
,
T.
Iwaizako
,
S.
Mizuno
,
M.
Tsuji
, and
H.
Ago
,
Phys. Chem. Chem. Phys.
19
(
12
),
8230
8235
(
2017
).
24.
Y.
Ogawa
,
B.
Hu
,
C. M.
Orofeo
,
M.
Tsuji
,
K.-I.
Ikeda
,
S.
Mizuno
,
H.
Hibino
, and
H.
Ago
,
J. Phys. Chem. Lett.
3
(
2
),
219
226
(
2012
).
25.
H.
Ago
,
Y.
Ito
,
N.
Mizuta
,
K.
Yoshida
,
B.
Hu
,
C. M.
Orofeo
,
M.
Tsuji
,
K.-I.
Ikeda
, and
S.
Mizuno
,
ACS Nano
4
(
12
),
7407
7414
(
2010
).
26.
H.
Ago
,
Y.
Ogawa
,
M.
Tsuji
,
S.
Mizuno
, and
H.
Hibino
,
J. Phys. Chem. Lett.
3
(
16
),
2228
2236
(
2012
).
27.
A.
Reina
,
H.
Son
,
L.
Jiao
,
B.
Fan
,
M. S.
Dresselhaus
,
Z.
Liu
, and
J.
Kong
,
J. Phys. Chem. C
112
(
46
),
17741
17744
(
2008
).
28.
R. V.
Gorbachev
,
I.
Riaz
,
R. R.
Nair
,
R.
Jalil
,
L.
Britnell
,
B. D.
Belle
,
E. W.
Hill
,
K. S.
Novoselov
,
K.
Watanabe
,
T.
Taniguchi
,
A. K.
Geim
, and
P.
Blake
,
Small
7
(
4
),
465
468
(
2011
).
29.
A. C.
Ferrari
,
J. C.
Meyer
,
V.
Scardaci
,
C.
Casiraghi
,
M.
Lazzeri
,
F.
Mauri
,
S.
Piscanec
,
D.
Jiang
,
K. S.
Novoselov
,
S.
Roth
, and
A. K.
Geim
,
Phys. Rev. Lett.
97
(
18
),
187401
(
2006
).
30.
N.
Yamamoto
,
E.
Rokuta
,
Y.
Hasegawa
,
T.
Nagao
,
M.
Trenary
,
C.
Oshima
, and
S.
Otani
,
Surf. Sci.
348
(
1
),
133
142
(
1996
).
31.
J. P. H.
Li
,
X.
Zhou
,
Y.
Pang
,
L.
Zhu
,
E. I.
Vovk
,
L.
Cong
,
A. P.
van Bavel
,
S.
Li
, and
Y.
Yang
,
Phys. Chem. Chem. Phys.
21
(
40
),
22351
22358
(
2019
).
32.
I. Z.
Mitrovic
,
M.
Althobaiti
,
A. D.
Weerakkody
,
V. R.
Dhanak
,
W. M.
Linhart
,
T. D.
Veal
,
N.
Sedghi
,
S.
Hall
,
P. R.
Chalker
,
D.
Tsoutsou
, and
A.
Dimoulas
,
J. Appl. Phys.
115
(
11
),
114102
(
2014
).
33.
G.
Kresse
and
D.
Joubert
,
Phys. Rev. B
59
(
3
),
1758
1775
(
1999
).
34.
G.
Kresse
and
J.
Furthmüller
,
Comput. Mater. Sci.
6
(
1
),
15
50
(
1996
).
35.
J. P.
Perdew
,
K.
Burke
, and
M.
Ernzerhof
,
Phys. Rev. Lett.
77
(
18
),
3865
3868
(
1996
).
36.
H.
Yamauchi
,
K.
Takagi
,
I.
Yuito
, and
U.
Kawabe
,
Appl. Phys. Lett.
29
(
10
),
638
640
(
1976
).
Published open access through an agreement with Alamos National Laboratory