Appl. Phys. Lett., Vol. 75, No. 1, 5 July 1999
Miller et al.
47
FIG. 2. In the abrupt displacement approximation, the continuous bending
of lattice planes is approximated as a step function.
To do this, image simulations were performed by inte-
grating the two-beam Howie–Whelan equations7,9 using the
MATHEMATICA software package. These are a system of
coupled first-order ordinary differential equations that can be
used to simulate contrast due to strained crystalline speci-
mens if a strain model is available.
FIG. 3. Simulated image contrast can be used to improve the accuracy of
strain measurements. As the slope of ␣(r)ϭ␣maxr/r0 increases, the extrema
move closer to the center of the island. Here, ␣max/(/2) varies from 5
͑lightest͒ to 6 ͑darkest͒. By analyzing the positions of the extrema in experi-
mental images, accurate strain measurements are possible.
Unfortunately, no simple analytical model exists for the
strain field due to isolated Ge islands on Si͑001͒. Other re-
searchers have shown that image simulations using finite el-
ement ͑FE͒ strain fields yield contrast that is qualitatively
similar to experimental images,11,12 however such simula-
tions have not been used for quantitative measurements. It is
difficult to measure sample properties such as strain using FE
simulations since it is difficult to adjust simulation param-
eters: one FE simulation is needed to generate each image
simulation.
Surprisingly, results can be well explained by using a
greatly simplified strain model. We find that when the dis-
placement field decays rapidly compared with the length
scale relevant to imaging, i.e., the electron extinction dis-
tance ͑typically ϳ100 nm͒, qualitative features can be pre-
dicted using a very simple approximation which we call the
abrupt displacement approximation ͑ADA͒. In this approxi-
mation, the continuous displacement field is replaced by a
step function a short distance into the substrate ͑Fig. 2͒. The
size of the step increases linearly from the center to the edge
of the island.
In other words, provided that the displacement field de-
cays rapidly enough—within 10–20 nm of the strain
source—the ADA allows us to replace any continuous dis-
placement field with a stacking fault-like9 abrupt displace-
ment which is at a fixed depth and has a variable magnitude.
This approximation is useful since TEM image contrast due
to stacking faults is well understood. Comparison of this
model with image simulations using full strain fields shows
that the ADA works well.
Properties of image contrast depend on the scalar prod-
uct of the diffraction vector used for imaging, g, and the
displacement field, u, via the quantity ␣ϭ2 g.u. If u in-
creases linearly parallel to g, then ␣ϭ2g(⑀r), for rϽr0
where r is parallel to g, r0 is the radius of the island, and ⑀
ϭ⑀rr , the radial strain in the substrate. Using these assump-
tions with the ADA, we find that image contrast is periodic
with ␣, nearly following a sinusoid. By counting extrema
that lie inside the radius of the island the maximum value of
the phase ␣max can be bounded. If n is the number of extrema
between the center of the island and the island’s edge, then
More precise measurements are made by analyzing the posi-
tions of intensity extrema inside the island’s radius.
Line profiles across islands imaged with a gϭ(400) dif-
fraction vector show three extrema, allowing us to place
bounds on ␣max. Analyzing the positions of extrema closest
to the particle edges improves the measurement. The average
position of the outermost extrema for three islands was found
to be 0.887 r0 , which, when compared with a contrast table
͑Fig. 3͒, leads to a value ␣m(4a0x0)ϭ8.98 (ϭ5.72/2). By
simple algebra, the maximum lattice plane displacement is
⑀r0ϭ0.19 nm and using a weak-beam TEM measurement of
the average island radius, we find that ⑀ϭ0.19 nm/22.5 nm
ϭ0.84Ϯ0.17%. We expect that by improving experimental
procedures, measurement error less than 10% will be readily
achievable.
Our measurement agrees well with FE simulations and
with other measurements of the strain. FE results show a
maximum displacement ⑀r0ϭ0.216 nm. Converting our
strain measurement to a stress ͑using 170 GPa for the Si
Young’s modulus and 0.262 for the Poisson ratio͒, we find
rrϭ1.6 GPa, a value which lies at the upper end of a re-
cently measured range of stress values for pure Ge on
Si͑001͒ during the ‘‘nucleation and growth’’ phase of evolu-
tion, prior to dislocation introduction.13 A combined x-ray
scanning tunneling microscopy ͑STM͒ measurement of
strain in an earlier growth phase, when the presence of
‘‘huts’’ or ‘‘pyramids’’ dominates, led to a strain measure-
ment at the edge of the island equal to 0.5%.14 Strain has also
been measured in islands composed of dilute SiGe alloys on
Si͑001͒, and values similar to those measured here have been
predicted and measured.6,15
We offer the following observations about the measured
strain value: first, since the measured maximum displace-
ment matches the results of FE simulations for pure Ge, this
suggests that Si has not diffused into the Ge to relax strain.
This stands in contrast to electron diffraction strain measure-
ments of pure Ge islands deposited and capped with Si at a
much higher temperature ͑750 °C͒. In this case the islands
were observed to incorporate nearly 40% Si into the Ge.16
Second, we note that the maximum displacement of lat-
tice planes, 0.19 nm, is equal to the Burgers vector for an
2nϪ1͒ Ͻ␣maxϽ 2nϩ1͒
.
2
2
edge dislocation in Si. Displacement of lattice planes greater
136.165.238.131 On: Fri, 19 Dec 2014 03:14:40