The composition of III–V semiconductor alloys with multiple group V elements results from a complex interaction of each group V species with each other and with group IIIs. Molecular beam epitaxy growth conditions, such as the group V absolute fluxes and flux ratios, substrate temperature, group III growth rate, the presence of surfactants, and the alloy's strain state, all affect the composition of InAs1 – xSbx. These factors are codependent in a manner that is far from completely established. In this work, the authors examine how the sign and degree of strain affects the film’s composition. In this study, the authors show that InAs1 – xSbx alloys have some ability to resist the creation of strain by self-adjusting the incorporation of the group V elements that would otherwise be changed in an unfavorable way (by, for example, a substrate temperature change). This self-latching to the substrate lattice constant is beneficial, since it means that the fluxes may not have to be adjusted as precisely as otherwise needed.

InAs1 – xSbx has the smallest direct bandgap among the common III–V alloys. This is one of many desirable properties of InAs1 – xSbx, setting it apart as the most viable III–V challenger to HgCdTe for long-wavelength infrared (LWIR) detector applications.1 InAs1 – xSbx also has strong spin–orbit coupling and the largest g-factor among semiconductors that can interface well with superconductors. These properties drive recent interest in InAs1 – xSbx for studies of modern physics phenomena such as topologically protected states.2 

The growth of mixed group III alloys, such as InzGa1 – zAs, is generally very straightforward. When grown using stable fluxes and typical substrate temperatures (for which re-evaporation is negligible), the resulting films have mole fractions that are predictable and based directly on group III flux ratios.

In contrast, we cannot predict the constituent mole fractions in mixed group V alloys based only on the group V flux ratios, as they are also dependent on the interaction between the group V's and the group III's. For instance, the composition of InAs1 – xSbx results from a complex interaction of As and Sb with each other and with indium. Other factors affect the composition to varying degrees, such as the substrate temperature,3,4 the presence of surfactants,3,4 the growth rate,5 and, as we discuss in this article, the strain state. The composition of InAs1 – xSbx, whether grown as a bulk alloy or into a strained layer superlattice (SLS), is dependent on these factors, and the codependence of each of these parameters on each other is far from completely established. Despite these challenges, the quality of unstrained, unrelaxed InAs1 – xSbx bulk alloys, or InAs1 – xSbx /InAs(1 – y)Sby SLS structures grown on GaSb (001) substrates has greatly improved in recent years through highly tuned molecular beam epitaxy (MBE) growth techniques. We have repeatedly demonstrated thick, lattice-matched InAs1 – xSbx films of various compositions with dislocation densities below the transmission electron microscopy (TEM) detection limit and strong optical properties.6 

In this paper, we focus on strain, which influences compositions in both mixed group III and group V alloys, as reported in studies of In1 – xGaxN,7 and GaAs1 – xBix.8 We observe that the sign and degree of strain also influences incorporation in InAs1 – xSbx. Strained InAs1 – xSbx films tend to resist changes in the group V mole fractions despite growth conditions that should promote compositional change. For instance, a past study showed InAs1 – xSbx /InAs1 – ySby SLS structures have an overall maximum amount of Sb that will incorporate depending on its mismatch with the substrate. Changing the period or individual layer thickness, or the substrate temperature, did not substantially increase the incorporation.9 In this study, we find that InAs1 – xSbx grown close to lattice-matching conditions resists changing its group V mole fractions in a manner that would further increase its mismatch with the GaSb substrate. This effect of self-latching to the substrate lattice constant is beneficial since it means that the fluxes may not have to be adjusted as precisely as otherwise needed.

We grew the samples in a solid source Gen II MBE system equipped with valved, cracked Sb and As cells. The GaSb substrates were quartered 2 in. wafers, epi-ready, and exactly oriented along [001]. We examined all samples with high-resolution x-ray diffraction (XRD) and selected samples with TEM and energy dispersive spectroscopy (EDS). We collected four XRD scans, including two symmetric (004) scans measured before and after a 180° rotation in phi, and two asymmetric (115) scans. The first (115) scan had a high angle of instance, and the second had a low angle of incidence and a phi rotation of 180°. This standard method measures the in and out of plane lattice parameters, as described by Krost et al.,10 and allows us to extract the strain and composition.

In Series A and B, we alter the composition and strain by changing the substrate temperature. In series C, we alter the strain and composition by changing the Sb/As flux ratio, which causes differing degrees of compressive and tensile lattice mismatch in stacked triple-layers. Section III details the growth parameters with the data in chart form for each series.

Table I summarizes all data and growth conditions for series A and B. The InAs1 – xSbx samples have two different Sb/As flux ratios for a large range of substrate temperatures as measured by pyrometer. Series A (samples A1–A11) was grown with an Sb/As flux ratio of 1.05 and temperatures ranging from 320 to 490 °C. For our typical growth temperature of 415 °C, we would expect that an InAs1 – xSbx film grown with Sb/As of 1.05 would have an Sb mole fraction of ∼0.48. Series B (samples B1–B11) had a flux ratio of 0.07 and a 340–490 °C temperatures range. Under this flux condition, we would expect InAs1 – xSbx to have an Sb mole fraction near lattice match (0.089) at our typical growth temperature of 415 °C.

TABLE I.

Growth conditions and data for samples in series A (flux ratio = 1.05) and B (flux ratio = 0.07). Sb/As refers to the flux ratio, T refers to the substrate temperature during the InAs1 – xSbx deposition, t is the film thickness, hc is the Matthews–Blakeslee critical thickness, f refers to the lattice mismatch, and Rel refers to the film relaxation.

SampleSb/AsT
(°C)
Sb mole fractiont
(Å)
hc
(Å)
fRel %
A1 1.05 320 0.650 7500 11.0 −0.037 100 
A2 1.05 335 0.580 7500 15.1 −0.033 100 
A3 1.05 350 0.559 7500 16.5 −0.031 100 
A4 1.05 370 0.528 7500 18.8 −0.029 100 
A5 1.05 385 0.498 7500 20.6 −0.028 100 
A6 1.05 400 0.507 7500 21.4 −0.028 100 
A7 1.05 415 0.484 7500 22.8 −0.027 99 
A8 1.05 425 0.445 7500 27.2 −0.024 97 
A9 1.05 450 0.436 7500 28.3 −0.023 97 
A10 1.05 475 0.386 7500 35.9 −0.020 92 
A11 1.05 490 0.388 6000 36.3 −0.020 96 
B1 0.07 340 0.102 2500 2011 0.001 100 
B2 0.07 355 0.101 2500 2209 0.001 100 
B3 0.07 370 0.092 2500 13732 0.000 77 
B4 0.07 385 0.089 2500 ∞ 0.000 
B5 0.07 400 0.089 2500 ∞ 0.000 
B6 0.07 415 0.089 2500 ∞ 0.000 
B7 0.07 430 0.082 2500 4048 0.001 48 
B8 0.07 445 0.068 2500 1137 0.002 85 
B9 0.07 460 0.061 6000 813 0.002 98 
B10 0.07 475 0.056 6000 656 0.002 100 
B11 0.07 490 0.053 6000 593 0.003 100 
SampleSb/AsT
(°C)
Sb mole fractiont
(Å)
hc
(Å)
fRel %
A1 1.05 320 0.650 7500 11.0 −0.037 100 
A2 1.05 335 0.580 7500 15.1 −0.033 100 
A3 1.05 350 0.559 7500 16.5 −0.031 100 
A4 1.05 370 0.528 7500 18.8 −0.029 100 
A5 1.05 385 0.498 7500 20.6 −0.028 100 
A6 1.05 400 0.507 7500 21.4 −0.028 100 
A7 1.05 415 0.484 7500 22.8 −0.027 99 
A8 1.05 425 0.445 7500 27.2 −0.024 97 
A9 1.05 450 0.436 7500 28.3 −0.023 97 
A10 1.05 475 0.386 7500 35.9 −0.020 92 
A11 1.05 490 0.388 6000 36.3 −0.020 96 
B1 0.07 340 0.102 2500 2011 0.001 100 
B2 0.07 355 0.101 2500 2209 0.001 100 
B3 0.07 370 0.092 2500 13732 0.000 77 
B4 0.07 385 0.089 2500 ∞ 0.000 
B5 0.07 400 0.089 2500 ∞ 0.000 
B6 0.07 415 0.089 2500 ∞ 0.000 
B7 0.07 430 0.082 2500 4048 0.001 48 
B8 0.07 445 0.068 2500 1137 0.002 85 
B9 0.07 460 0.061 6000 813 0.002 98 
B10 0.07 475 0.056 6000 656 0.002 100 
B11 0.07 490 0.053 6000 593 0.003 100 

Figure 1 plots the Sb mole fraction as determined by x-ray as a function of growth temperature. The horizontal dashed line denotes the lattice-matched Sb mole fraction of 0.089.

FIG. 1.

Antimony mole fraction as a function of substrate temperature for samples grown with Sb/As flux ratios of 1.05 (A) and 0.07 (B).

FIG. 1.

Antimony mole fraction as a function of substrate temperature for samples grown with Sb/As flux ratios of 1.05 (A) and 0.07 (B).

Close modal

Overall, the Sb mole fraction decreases with increasing growth temperature for the substrate temperatures examined in this study. We discussed growth temperature effects on Sb incorporation in our previous study.3,4 Haugan et al. more recently observed a similar decrease in Sb in InAs1 – xSbx for increasing growth temperatures over the range of 400–440 °C.11 In Fig. 1, we notice a striking plateau in this trend over a ∼50 °C temperature range for compositions near lattice match in series B.

Figure 2 illustrates this trend clearly by plotting the percentage change in Sb mole fraction versus the percentage change in substrate temperature. The baseline is the Sb mole fraction for the two samples grown at 415 °C. Negative changes in mole fraction percentages refer to a decrease in the Sb mole fraction relative to the sample grown at 415 °C, and negative changes in temperature percentages refer to temperatures lower than 415 °C. For visual clarity, we only labeled the endpoints of each series. Samples A1–A11 are completely relaxed, and the compositions are not affected by strain; therefore, we can easily fit a linear relationship for this temperature range. The behavior of the B series as shown in Fig. 2 indicates a resistance to growth temperature-induced composition changes for films grown near lattice-matched conditions.

FIG. 2.

Percent change of mole fraction vs substrate temperature as compared with samples grown at 415 °C.

FIG. 2.

Percent change of mole fraction vs substrate temperature as compared with samples grown at 415 °C.

Close modal

All samples except B1–B8 have thicknesses that are orders of magnitude above the Matthews–Blakeslee12 critical thickness for InAsSb. Although some films retain strain well above critical thickness,13 it is reasonable to expect samples that are magnitudes of order thicker to be completely relaxed. This assumption proves reasonable in the XRD measurements. Figure 3 plots the XRD measured relaxations as a function of (a) mole fraction and (b) growth temperature. The A series films are highly mismatched to the substrate, and are almost completely relaxed, as shown in Figs. 1 and 3. Therefore, we do not expect strain influences in the group V mole fractions for A1–A11. The linear composition/substrate temperature relationship shown in Figs. 1 and 2 for series A affirms this expectation.

FIG. 3.

InAs1 – xSbx film relaxation as a function of (a) Sb mole fraction and (b) substrate temperature. For visual clarity, the data points were not all labeled in (a), but the sample name increments numerically from left to right for both series A and B.

FIG. 3.

InAs1 – xSbx film relaxation as a function of (a) Sb mole fraction and (b) substrate temperature. For visual clarity, the data points were not all labeled in (a), but the sample name increments numerically from left to right for both series A and B.

Close modal

For the B series, the Sb mole fraction plateaus for a temperature range of around 50°. These mole fractions are at or extremely close (within 0.0025) of the lattice match mole fraction. It implies a “latching effect” where the film has a limited but noticeable ability to incorporate selectively the group V species needed to minimize strain. Again, we chose this Sb/As ratio because it reliably produces lattice-matched films for grown at 415 °C and measured by XRD at room temperature. We notice in Fig. 3 that the “latching effect” seems more pronounced on the “colder” side of this temperature.

Figure 4 plots relaxation as a function of lattice mismatch for series B. Relaxation for compressive and tensile films is not symmetric about the lattice-matched position. This is not surprising since there are numerous reports showing a larger critical thickness for tensile strained versus compressively strained films.14,15

FIG. 4.

Film relaxation as a function of lattice mismatch for series B.

FIG. 4.

Film relaxation as a function of lattice mismatch for series B.

Close modal

Since the “latching effect” occurs for substrate temperatures ranging from 385 to 415 °C, it might be worth changing our standard growth temperature for InAs1 – xSbx from 415 °C to the middle of this range (∼400 °C). We have not yet optimized the growth temperature for InAs1 – xSbx. In a future work, we will determine if the optical and crystalline properties of InAs1 – xSbx grown at temperatures closer to the middle of the “latching window” improve relative to films grown at 415 °C. The findings in this study indicate that InAs1 – xSbx may be tolerant to small deviations in flux or substrate temperature. It may be beneficial to grow materials in the middle of the temperature range that enhances this tolerance.

Series C was initially grown to expand our dataset for establishing the nonlinear relationship of the Sb/As flux ratio to the InAsSb composition. Multilayer stacks allow multiple data points from a single growth, saving substrates and growth time. We also wanted to examine the compressive versus tensile defect morphology in InAsSb.

We grew two samples with three InAs1 – xSbx layers at a substrate temperature of 415 °C. The targeted Sb/As flux ratios were 0.09, 0.18, and 0.27. Figure 5 shows a schematic of the triple layer structures C1 and C2. For Sample C1, we began by depositing the layer that we expected to be closest to (although slightly in excess of) the lattice match mole fraction and then increased the Sb/As ratio in two steps, resulting in three 0.75 μm thick layers. We reversed the layer order for Sample C2, etc., i.e., we deposited the most highly mismatched, higher Sb content first and then reduced the Sb concentration for the next two layers.

FIG. 5.

Schematics of the triple layer structures of set B, with the Sb/As flux ratio listed for each layer. L1 refers to layer 1, etc.

FIG. 5.

Schematics of the triple layer structures of set B, with the Sb/As flux ratio listed for each layer. L1 refers to layer 1, etc.

Close modal

The group III cell temperatures and the group V valve position for C1 and C2, along with the shutter sequence (which is the same for both samples), are shown in Fig. 6.

FIG. 6.

Cell temperatures for group III's and valve positions for the group V elements for (a) sample C1 and (b) sample C2. A valve position of 300 is “fully opened.” The shutter sequence (c) was identical for both samples.

FIG. 6.

Cell temperatures for group III's and valve positions for the group V elements for (a) sample C1 and (b) sample C2. A valve position of 300 is “fully opened.” The shutter sequence (c) was identical for both samples.

Close modal

Figure 7 shows the recorded chamber and beam flux monitor pressure throughout the growth. The pressure gauges logged the same overall system pressures versus time (which are mostly driven by As) corresponding to the same Sb/As ratios for samples C1 and C2. From this and the plot of the growth logs in Fig. 6, we can see that the source cell temperatures and the flux conditions were the same in layer 1, sample C1 and in layer 3, sample C2, just as layer 3 in C1 had the same growth conditions as layer 1 in C2. The second layer in both samples had the same growth conditions. Therefore, the only difference between the layers in C1 and C2 that had the same growth conditions was its lattice mismatch with the underlying layer.

FIG. 7.

Chamber pressure throughout the growth for samples C1 and C2.

FIG. 7.

Chamber pressure throughout the growth for samples C1 and C2.

Close modal

Figure 8 shows the (004) x-ray rocking curves for samples C1 and C2. At first glance, it seems that although exposed to the same flux and cell temperature conditions, layer 1 in C1 and layer 3 in C2 have different mole fractions, as listed in the schematic of Fig. 9. Layer 3 in sample C2 is on the right hand side of the substrate peak, whereas layer 1 in sample C1 is just to the left of the substrate peak and near lattice match. The panalytical epitaxy software finds that layer 3, sample C2 has a mole fraction of 0.02 Sb vs 0.09 for layer 1, sample C1. All of the other layers in the two samples give nominally the same mole fraction for layers grown with the same Sb/As ratio. The relaxations listed are relative to the GaSb substrate (Fig. 10).

FIG. 8.

(004) rocking curves for samples C1 and C2.

FIG. 8.

(004) rocking curves for samples C1 and C2.

Close modal
FIG. 9.

Schematic with relaxation and mole fraction of InAs1 – xSbx layers according to XRD.

FIG. 9.

Schematic with relaxation and mole fraction of InAs1 – xSbx layers according to XRD.

Close modal
FIG. 10.

EDS profile of InAs1 – xSbx layers for sample C2. The left side of the profile with 100% Sb is from the GaSb substrate, followed by the three InAsSb layers of decreasing Sb composition, denoted by L1, L2, and L3.

FIG. 10.

EDS profile of InAs1 – xSbx layers for sample C2. The left side of the profile with 100% Sb is from the GaSb substrate, followed by the three InAsSb layers of decreasing Sb composition, denoted by L1, L2, and L3.

Close modal

The 0.02 Sb mole fraction in layer 3, sample C2 does not seem plausible. We examined sample C2 with EDS. Figure 10 shows the mole fraction profile, and the XRD and EDS data are plotted together in Fig. 11.

FIG. 11.

Mole fraction of InAs1 – xSbx layers (L1 = layer 1, etc.) for samples B1 and B2 for according to XRD and EDS.

FIG. 11.

Mole fraction of InAs1 – xSbx layers (L1 = layer 1, etc.) for samples B1 and B2 for according to XRD and EDS.

Close modal
FIG. 12.

(a) (002) DF and (b) (220) DF images of sample C1. (c) HRTEM of layer 3.

FIG. 12.

(a) (002) DF and (b) (220) DF images of sample C1. (c) HRTEM of layer 3.

Close modal

We notice in Fig. 11 that nearly fully relaxed layers exhibit the same mole fraction in C1 and C2 for the same Sb/As ratio. The EDS and XRD data for sample C2 also agree for layers 2 and 3.

There is a large disparity between the EDS data and the XRD data for sample C2, layer 3. TEM results, which we will discuss later, show cracks, and non-normal relaxation in layer 3, which is likely affecting the XRD results. Since the XRD and EDS are in reasonable agreement for the two other layers in the sample, and the EDS value of 0.13 for layer 3 is more in line with expectations and would be unaffected by unusual diffraction due to tensile strain relaxation, we believe that the EDS value is much more reliable.

Layer 1 in sample C1 is the only layer that we expected to show strain effects. For an Sb/As ratio of 0.09, we expect a mole fraction of ∼0.125. This is close enough to lattice match that it is unlikely to relax fully. The XRD indicates that it is 13% relaxed, with a mole fraction of 0.09, as shown in Fig. 11. The film has incorporated Sb nearer to lattice match conditions than the Sb/As flux ratio would suggest. Compare this to the mostly relaxed layer 3 in sample C2, which has a 0.13 Sb mole fraction (which is closer to that expected for Sb/As = 0.09) despite being grown with the same Sb/As ratio as sample C1, layer 1.

As with the samples in the first two series (A and B), we see that strain affects incorporation near lattice-matched conditions (layer 1, sample C1), whereas the other layers (relaxed) in samples C1 and C2 are unaffected by strain and have nearly the same mole fraction for the same Sb/As ratio. For nearly fully relaxed films, the sign of the mismatch (compressive versus tensile) seems to have no effect on mole fraction.

Figure 12 shows (a) (002) and (b) (220) dark field (DF) TEM images of sample C1. The defect morphology is consistent with that expected for compressive strain relaxation.16 The curved lines move with the beam and are thickness fringes. Figure 12(c) is an HRTEM image of layer 3, showing no unusual strain contrast, with flat, undistorted lattice fringes measuring as expected for a relaxed layer InAs0.74Sb0.26. Since tetragonal distortion and dislocation relaxation accommodate the strain relaxation, simple XRD analysis of the four collected scans gives accurate values for the strain and composition.

Figure 13 shows (002) DF images of sample C2. Layers 2 and 3 are tensile, and the sample relieves some of the strain by cracking, particularly in layer 3, and with surface undulations (as seen, for example, at the top of the crack on the left hand side). Simple XRD analysis measuring tetragonal distortions is, therefore, not appropriate for calculating the strain and the relaxation in these films.

FIG. 13.

(002) Bright Field (BF) image for sample C2.

FIG. 13.

(002) Bright Field (BF) image for sample C2.

Close modal

InAs1 – xSbx is of interest because of its potential as an LWIR detector and for studies of modern physics phenomena such as topologically protected states. Through well controlled MBE growth and structural/chemical analysis, we have improved the quality of unstrained, unrelaxed InAs1 – xSbx across a wide mole fraction spectrum. This required a detailed study of mixed group V alloying, which is not straightforward due to the complex interaction of As and Sb with each other and with In.

In this study, we demonstrate that InAs1 – xSbx grown close to lattice matching resists changing its mole fraction in a manner that would further increase its mismatch with the GaSb substrate. We also determine that relaxation is not symmetric for tensile and compressive strains in InAs1 – xSbx. Layers that are nearly fully relaxed are largely unaffected by strain in terms of incorporation for a given Sb/As ratio, whether the film is compressive or tensile relative to the underlying layer.

It may be possible to optimize future growths by better defining the center of the growth condition window where the “latching effect” occurs, which would make the material more tolerant of normal deviations in flux and substrate temperature that occur over the time period required for MBE growth.

Stony Brook University acknowledges support of the Army Research Office (Award No. W911NF-16-2-0053).

1.
A. J
Ciani
,
C. H.
Grein
,
W. L.
Sarney
,
S. P.
Svensson
,
D. V.
Donetsky
, and
G. L.
Belenky
, “
A comparison of indium arsenide antimonide and mercury cadmium telluride as long wavelength infrared detector materials
,”
J. Electron. Matter.
(submitted).
2.
W.
Mayer
 et al., arXiv:1909.12571 (
2019
).
3.
W. L.
Sarney
,
S. P.
Svensson
,
E. M.
Anderson
,
A. M.
Lundquist
,
C.
Pearson
, and
J. M.
Millunchick
,
J Cryst. Growth
406
,
8
(
2014
).
4.
S. P.
Svensson
,
W. L.
Sarney
,
B. C.
Connelly
,
E. M.
Anderson
, and
J. M.
Millunchick
,
J. Cryst. Growth
425
,
234
(
2015
)
5.
W. L.
Sarney
and
S. P.
Svensson
,
J. Vac. Sci. Technol. B
33
,
060604
(
2015
).
6.
G.
Belenky
,
D.
Donetsky
,
G.
Kipshidze
,
D.
Wang
,
L.
Shterengas
,
W. L.
Sarney
, and
S. P.
Svensson
,
Appl. Phys. Lett.
99
,
141116
(
2011
).
7.
M. V.
Durnev
,
A. V.
Omelchenko
,
E. V.
Yakovlev
,
I. Yu.
Evstratov
, and
S. Yu.
Karpov
,
Phys. Status Solidi A
208
,
2671
(
2011
).
8.
L.
Nattermann
,
P.
Ludewig
,
E.
Sterzer
, and
K.
Volz
,
J. Cryst. Growth
470
,
15
(
2017
).
9.
W. L.
Sarney
,
S. P.
Svensson
,
M. K.
Yakes
,
Y.
Xu
,
D.
Donetsky
, and
G.
Belenky
,
J. Appl. Phys.
124
,
35304
(
2018
).
10.
A.
Krost
,
G.
Bauer
, and
Joachim
Woitok
, “High resolution x-ray diffraction,” in Optical Characterization of Epitaxial Semiconductor Layers, edited by G. Bauer and W. Richter (Springer, Berlin, 1996), Chap. 6, pp. 286–391.
11.
H. J.
Haugan
,
G. J.
Brown
, and
J. A.
Peoples
,
J. Vac. Sci. Technol. B
35
,
02B107
(
2017
).
12.
J. W.
Matthews
and
A. E.
Blakeslee
,
J. Cryst. Growth
27
,
118
(
1974
).
13.
W. L.
Sarney
,
L.
Salamanca-Riba
,
J.
Bruno
, and
R.
Tober
,
Solid State Electron.
46
,
1643
(
2002
).
14.
G.
Salviati
,
C.
Ferrari
,
L.
Lazzarini
,
L.
Nasi
,
A. V.
Drigo
,
M.
Berti
,
D.
De Salvador
,
M.
Natali
, and
M.
Mazzer
,
Appl. Surf. Sci.
188
,
36
(
2002
).
15.
P.
Maigne
,
M.
Gendry
,
T.
Venet
,
Y.
Tahri
, and
G.
Hollinger
,
Appl. Phys. Lett.
69
,
682
(
1996
).
16.
P. M. J.
Marée
,
J. C.
Barbour
,
J. F.
van der Veen
,
K. L.
Kavanagh
,
C. W. T.
Bulle-Lieuwma
, and
M. P. A.
Viegers
,
J. Appl. Phys.
62
,
4413
(
1987
).