Abstract
Indium arsenide quantum dots are of great interest for next-generation telecom optoelectronics if their emission wavelength can be red shifted into the correct range. One method to achieve this is the deposition of a surfactant, such as bismuth, during quantum dot growth. Here, we present a series of indium arsenide quantum dot layers grown using several bismuth fluxes and two different growth temperatures. The effects of bismuth flux on quantum dot morphology and optical properties are studied by atomic force microscopy and photoluminescence measurements. Bimodal distributions of quantum dots are seen at low growth temperature, while at high temperature, a single dominant distribution is seen in most of the layers. A medium bismuth flux was seen to produce the highest integrated photoluminescence intensity at high growth temperature, whereas intensity saturates between medium and high fluxes at low growth temperatures. A significant increase in uncorrected aspect ratio seen for the layer grown with a low bismuth flux at high growth temperature presents a new opportunity for control of quantum dot morphology using bismuth.
1. Introduction
Since the 1980s there has been interest in using 3-dimensional quantum confinement to alter the optical and electrical properties of semiconductors [1, 2]. One method to achieve this is to embed nanoscale semiconductor blocks within a wider band gap semiconductor. These are referred to as quantum dots (QDs). Self-assembling indium arsenide (InAs) QDs grown on gallium arsenide (GaAs) have been studied for many years [3, 4] and are of key interest due to their potential in telecoms optoelectronics both as single photon emitters and laser diodes [5, 6]. There are various challenges to using InAs QDs in these applications, but most important is red shifting the emission wavelength into the telecoms O and C bands at 1300 and 1550 nm, respectively. This has been attempted by various means, including growing larger QDs, strain/band engineering in the capping layer, or using a surfactant during QD growth. Growing thicker InAs depositions to produce larger QDs can lead to the formation of dislocated InAs islands however and negatively affect the layer properties [7, 8]. Strain engineering using metamorphic buffers has yielded promising results [9–11] and growth on indium phosphide (InP) has produced devices operating at 1550 nm [12, 13].
High QD uniformity is also desirable in these structures, particularly for laser applications in order to reduce threshold currents [5]. Surfactant deposition during QD growth has been shown to have a significant impact on the size and uniformity of InAs QDs as well as impacting QD density [14–24].
Antimony (Sb) has been studied extensively as both a capping layer material, using GaAsSb [25–27], and as a surfactant during different stages of QD growth [23, 24, 28]. By capping with GaAsSb, dramatic red shifts to QD emission have been demonstrated [25, 26]. As a pregrowth layer Sb has been shown to increase QD density [29], when used as a surfactant during QD formation, however, Sb causes a blue shift in emission. This has been shown by Lu et al. on (001) GaAs substrates and was attributed to a reduction in QD size [23]. Matsuura et al. also observed this effect and ascribed it to an increase in the Stranski-Krastanov critical thickness [24].
Bismuth (Bi) can also be used as a surfactant during InAs QD growth. It is reported to lead to a reduction in QD density [16–18], and contrary to Sb, it has been shown to increase the average height of QDs [8, 16–18]. However, there are reports which have indicated that the effect of Bi is not straightforward. Dasika et al. [8] found that for InAs depositions of ≤2.3 monolayers (ML), Bi reduced QD density, but for depositions of ≥2.6 ML, a Bi flux increased QD density. It was also noted that Bi reduced the dislocated island density for all InAs depositions (2.3–3.3 ML), with the largest reduction for the thickest depositions. Chen et al. [19] found that a Bi flux reduced the density of QDs in the growth temperature range of 475–492°C but increased QD density for temperatures of 492–500°C. They also noted that the QDs became more uniform and the dislocated island density was reduced. Reyes et al. [17] found that low Bi ML/s slightly reduced photoluminescence (PL) intensity with small emission red shifts, but a flux of 0.06 ML/s caused a significant 165 meV red shift and a 68 times improvement of PL intensity.
As gallium arsenide bismide (GaAsBi) exhibits a large reduction in band gap per % Bi in GaAs [30], Bi is also a promising material for QD capping layers. Wang L et al. [20] studied the effect of Bi as a surfactant and constituent in the barrier and capping layers. They found that a 5% GaAsBi buffer reduced both the height and density of subsequent QDs when compared to GaAs. This was interpreted to indicate a delay in InAs QD nucleation caused by Bi segregation at the increased substrate temperature of 500°C for QD growth. When used in a capping layer grown at 280°C, 3% Bi induced a significant red shift of 163 meV in the emission wavelength. When the capping layer was grown at 500°C, thereby segregating all the supplied Bi; it was found to mitigate density loss during capping and improve QD uniformity.
In this work, we provide a fundamental study of the effects of Bi on InAs QDs using two layer series produced at growth temperatures of 380°C and 510°C each using several Bi fluxes. These temperatures were selected to document the effect of Bi in two regimes. The first, 380°C, is commonly used in growth of GaAsBi at which up to ~6% Bi can be incorporated [31]. The second, 510°C, is in the upper range of standard temperatures for InAs QD growth and will limit Bi to a surfactant role. Atomic force microscopy (AFM) and ensemble PL measurements are used to investigate the effects of different Bi fluxes on the morphology, density, and optical properties of the InAs QDs as well as their nucleation at low growth temperatures.
2. Methodology
The QD layers were grown in an Omicron MBE-STM system on semi-insulating (001) GaAs substrates measuring mm. The substrates were initially heated to ~620°C under an arsenic (As) overpressure to desorb the native oxide. Following this, a 300 nm GaAs buffer was deposited at 580°C using an As2:Ga atomic flux ratio of 1.7. Substrates were then cooled to the QD growth temperature of 380°C or 510°C and the As species adjusted to As4. The change in As species was done to widen the As flux window for efficient incorporation of Bi compared to As2, which has been discussed elsewhere [32]. Next, ML of InAs was deposited at 0.01 ML/s using an atomic flux ratio As4:In of 30 in the presence of an additional bismuth flux. After a 10 s pause to change the As flux via a needle valve, this QD layer was buried under a 60 nm GaAs spacer grown at the same temperature as the QDs, with an atomic flux ratio As4:Ga of 2.2. The spacer layer was added to facilitate optical measurements and prevent QD coupling [33]. Following this, an uncapped repeat QD layer was deposited to allow surface studies. The general structure of the layers is depicted in Figure 1, and the key growth conditions are included in Table 1. The GaAs growth rate was calibrated by observing RHEED oscillations and the InAs growth rate was calibrated by observing the RHEED transition at the onset of QD nucleation at 1.6 ML [34]. The As:III atomic flux ratios were calibrated using transitions between As-rich and Ga-rich reconstructions at the GaAs growth rate of 0.36 ML/s. The authors acknowledge that the conditions and structure displayed in Figure 1 do not represent the current state-of-the-art for InAs QD growth and therefore the InAs control layers grown without Bi will not display ideal characteristics. The reasoning for this basic design was to prevent ambiguity on the impact of the Bi flux on the QDs; as in the case of optimised conditions, further changes to the surface chemistry with the addition of a Bi flux would be expected to alter the previously tuned conditions and complicate any independent analysis. By instead starting with a traditional QD control layer, the full effect of the Bi flux can be observed unambiguously as has been performed on other GaAs orientations [18].

Cleaved pieces of each layer were analysed by AFM using a nominally 7 nm wide tip in 3-4 different locations with each scan covering 1 μm2. This was to ensure a suitably large QD population was identified for each layer to perform statistical analysis on. This AFM data is presented as recorded without corrections. It is important to note that AFM is not a direct representation of the buried QDs due to compositional changes which occur during overgrowth [35, 36]. Room temperature PL measurements were taken using a liquid nitrogen cooled Ge photodetector, Horiba monochromator, and a chopped 532 nm laser operated at excitation powers of 90-900 mW. Analysis of the AFM data was performed using custom, homemade MATLAB software designed to identify individual QDs and extract their morphology characteristics. Large QDs which had minor overlap, <2 nm, due to close proximity were manually separated using this software and characterised independently. QDs with overlap larger than this which were beginning to coalesce were treated a single QD. QD populations identified from each scan were then combined into a single population for each layer. QD densities were calculated from these combined populations based on the total area covered, and height distributions were calculated as histograms indicating the probability of a QD being at a given height.
3. Results and Discussion
AFM images of the layers are shown in Figure 4. It can be seen for layer L0 that at low growth temperature and without Bi, there is negligible QD formation. There appear to be some QDs starting to form, but these are at a low density and embedded in the rough wetting layer. With the addition of a Bi flux, there is a dramatic change in surface morphology where a high number of QDs form for L1. This effect is reminiscent of recent results by Lewis et al. who observed abrupt QD formation on GaAs (110) surfaces with the addition of a Bi flux [18]. With further increasing Bi flux, taller QDs are formed in layers L2 and L3. At high growth temperature, QDs are formed in all layers. The QDs grown at high temperature are also visibly more uniform compared to their low temperature counterparts.
Figure 2 displays the AFM QD height distributions for all of the samples. As seen in the raw AFM images, as the Bi flux is increased, the peak QD height increases at both growth temperatures. The layers grown at 380°C all display bimodal distributions, with a large number of short QDs between 1 and 4 nm in height, distinct from the population of taller QDs which are generally 4–13 nm tall. The exception to this is L0 which did not produce any tall QDs. The distribution of short QDs is only seen in the layers grown at 380°C and is believed to originate from 3D growth transitions below 1.6 ML which have been seen to occur at low growth temperatures [34–39]. At high growth temperature, most of the layers have single distributions and random spreads of QDs at shorter heights. H1, however, does appear display bimodal characteristics which has been observed in QD morphology comparisons and will be discussed in Figure 3. Common to all the high temperature layers are peak QD heights which are consistently higher than the low growth temperature equivalents.






Figure 5 shows the areal densities of the QDs. As the layers grown at 380°C displayed a bimodal QD height distribution, the densities of the two QD distributions have been separately recorded and were calculated using the Gaussian fits shown in Figure 2. Due to both distributions for H1 being much taller than the small dots observed in the low growth temperature layers, these have been combined in Figure 5. For layer L1, the density of tall QDs is almost identical to that of the high temperature equivalent, H1. As the Bi flux is increased at low temperature, the density of tall QDs decreases and the number of shorter QDs increases, resulting in a higher density of short QDs than tall QDs for layer L3. The impact of Bi on the QD density of the high growth temperature layers is simpler; at a low Bi flux, there is no significant change in QD density. However, as the Bi flux is increased for layers H2 and H3, there is a large decrease in QD density to 60–70% of the H1 density. This is similar to the results seen by Fan et al. [16] who noted a 50% decrease in QD density when using a Bi surfactant.
Using the area covered by each QD as identified in the MATLAB software and assuming the QDs were circular in base, an indicative aspect ratio (AR) for each QD was calculated from the uncorrected AFM data by dividing the QD height by estimated diameter. These are plotted for each layer in Figure 3 as a scatter of AR against the QD height. L0 has been omitted due to large uncertainty and the small sample size of its QD population. A dashed grey line has been added to each layer using a fit to layer H0 which had the tightest data spread and most closely represents standard InAs QD growth in this series. Looking at the low growth temperature layers, which displayed bimodal distributions in QD height, it can be seen that both height distributions follow a similar evolution in AR with height as shown by the H0 fit line. The separation of the shorter and taller QD distributions driven by increasing Bi flux can be also seen, as in Figure 2. The higher growth temperature layers generally show a tighter spread of AR compared to the low growth temperature layers, except for H1 which has two distinct distributions of QDs, which was not distinctly clear from Figure 2. Whilst these two distributions overlap in height, the taller distribution has a significantly larger aspect ratio, which is not seen in any of the other layers. The average aspect ratio for these taller QDs is approximately 0.31. A similar shift in AR to a higher value has been seen in this temperature range by Saito et al. [40]. They concluded that the driving force behind this shape transition was a combination of QD height and base strain leading to a higher aspect ratio becoming more energetically favourable. As this is only observed for the lowest Bi flux layer at this temperature and not at higher fluxes or without it, it appears that it may be possible to use Bi to control this transition.
Figure 6 shows the room temperature PL spectra taken for the layers. All of the layers display a low energy emission between 0.7 and 0.9 eV, referred to as a low energy feature (LEF); the source of this is considered to be from QD clusters and chains which have been observed by TEM elsewhere [41]. These have been seen on the surface by AFM, and due to QD dissociation during GaAs overgrowth [35, 36, 42], their densities are expected to be larger in the buried QD layers. It is noteworthy that in [41], this emission was only seen in layers grown at 300 and 450°C and not in a layer grown at 520°C. It is present here in all of the 510°C layers shown in Figure 6, however, and the intensity changes little across the range of Bi fluxes, suggesting that a Bi surfactant has little effect on the formation of these clusters. As expected, due to the associated reduction of antisites and other defects, the higher growth temperature layers produced stronger QD emission in all cases. The lower growth temperature layers displayed longer wavelength emission however, with L2 and L3 having a ground state peaks at 0.968 eV. It is possible that this comparative red shift is due to a change in the QD In content between the two growth temperatures. However, as the growth rate is low, 0.01 ML/s, this is not believed to be the case as other reports have found that low growth rates reduced In-Ga intermixing [43]. Another potential cause is the incorporation of a small percentage of Bi into the QDs grown at 380°C although it is still possible that the difference in emission wavelength is produced by morphology changes during the capping process which have been demonstrated to vary dramatically over a similar temperature range [36]. The longest wavelength achieved from the layers grown at 510°C came from H1 with a ground state peak at 0.987 eV. This layer did not have the tallest QDs from the 510°C series nor is it expected that the capping process would be any different for this layer compared to the others grown at this temperature. Also the solubility of Bi at this temperature prohibits incorporation. It is therefore expected that this red shift is due to emission originating from the QDs of large aspect ratio, seen in Figure 3. Whilst the statistics in Figure 3 were derived from the surface QD population, it is expected that the optically active, buried QDs will have retained a bimodal distribution despite compositional and morphology changes which occur during the capping process. The large emission red shift caused by a change in aspect ratio has been modelled by Usman et al. [44] and may be expected to be larger than the blue shift from reduced QD height. Figure 7 shows how the integrated PL (IPL) intensity of the layers varies with bismuth flux at each growth temperature. It is clear that the use of any bismuth flux produces an increase in IPL intensity at both growth temperatures. Comparing the IPL from the LEF in the layers grown at 380°C shows no clear trend; however, in the layers grown at 510°C, the IPL of the LEF increases with a decrease in dot density. The authors note that these trends in IPL intensity for both the QDs and LEF closely mimic the peak PL intensity trends for all the layers.
4. Conclusions
In summary, a systematic series of quantum dot layers grown under various bismuth fluxes and growth temperatures has been presented. Quantum dot nucleation at low growth temperatures was induced through the use of a bismuth flux as demonstrated by layers L0 and L1. In all layers which formed quantum dots, there was a consistent increase in quantum dot height with increasing bismuth flux and a decrease in the density of large quantum dots at medium and high bismuth fluxes for both growth temperatures. Integrated photoluminescence intensity was improved compared to growth without Bi but appears to saturate between the medium and high Bi flux. The change in surface quantum dot aspect ratio observed at higher growth temperatures in layer H1 suggests that bismuth could be used to control quantum dot shape. If replicable at lower growth temperatures, where bismuth may incorporate into the dot structure, this could provide an extra pathway for MBE growth to produce quantum dots emitting at even longer wavelengths.
Data Availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declare no conflicts of interest in this work.
Authors’ Contributions
NJB wrote the manuscript, recorded PL data, and analysed the PL and AFM data. NJB and MRC synthesised the layers. RDR and JPRD oversaw the project.
Acknowledgments
The authors would like to gratefully acknowledge Dr. Deborah Hammond and Prof. Graham Leggett from the Sheffield Surface Analysis Centre (SSAC) for the AFM measurements and Dr. Edmund Clarke for useful discussions on quantum dot growth. This work has been supported by an EPSRC grant (EP/S036792/1). The work of RDR was supported by the Royal Academy of Engineering under the Research Fellowships scheme.