We report upon a comprehensive investigation of the subthreshold characteristics of the ballistic electron emission microscopy (BEEM) current in ballistic electron emission spectroscopy. Starting from the Bell-Kaiser model, we derive an analytical equation to describe the subthreshold behavior of the BEEM current. It is found that the BEEM current in this region should exhibit a subthreshold swing of ∼60 mV/decade at room temperature, which we experimentally verified. This finding provides a rule of thumb for the detectability of the subthreshold behavior in a spectrum. For spectra where the subthreshold behavior is discernible above the signal noise, it is demonstrated that significant deviations in the near-threshold region can occur when fitting with a simple quadratic model that ignores the subthreshold behavior. To take the subthreshold behavior into account, a simple analytical model is proposed. This model not only fits significantly better in the near threshold region than the square model, but also gives a barrier height closer to the one extracted from the Bell-Kaiser model. More significantly, this model provides a quick method to estimate the subthreshold BEEM current amplitude based on the BEEM current above the barrier height.

1.
W. J.
Kaiser
and
L. D.
Bell
,
Phys. Rev. Lett.
60
,
1406
(
1988
).
2.
L. D.
Bell
and
W. J.
Kaiser
,
Phys. Rev. Lett.
61
,
2368
(
1988
).
3.
W.
Yi
,
A. J.
Stollenwerk
, and
V.
Narayanamurti
,
Surf. Sci. Rep.
64
,
169
(
2009
).
4.
Here, threshold is defined as the voltage corresponding to the extracted barrier height value.
5.
C.
Tivarus
, Ph. D. thesis (
The Ohio State University
, Columbus,
2005
).
6.
M.
Razavy
,
Quantum Theory of Tunneling
(
World Scientific
,
Singapore
,
2003
).
7.
C. J.
Chen
,
Introduction to Scanning Tunneling Microscopy
, 2nd ed. (
Oxford University Press
,
New York
,
2007
).
8.
The extracted barrier heights from the Bell-Kaiser fitting are rather insensitive to the choice of the average workfunction in the range of 3–4.7 eV and the tunneling gap distance of 7–15 Å; this has also been pointed out in Ref. 9.
9.
G. N.
Henderson
,
P. N.
First
,
T. K.
Gaylord
, and
E. N.
Glytsis
,
Phys. Rev. Lett.
71
,
2999
(
1993
).
10.
C.
Tivarus
,
J. P.
Pelz
,
M. K.
Hudait
, and
S. A.
Ringel
,
Phys. Rev. Lett.
94
,
206803
1
(
2005
).
11.
M.
Prietsch
and
R.
Ludeke
,
Phys. Rev. Lett.
66
,
2511
(
1991
).
12.
A.
Bannani
,
C.
Bobisch
, and
R.
Möller
,
Science
315
,
1824
(
2007
).
13.
In order to improve the fitting in the near threshold region, one can introduce a factor, f1.2, to compensate for the approximations made in the model by changing (eV-ΦBH)/(kT) to (eV-ΦBH)/(fkT). We find that introducing this factor can improve the fitting in the near threshold region, and the extracted barrier height value does not significantly change (∼0.01 eV or less).
You do not currently have access to this content.