Next Article in Journal
Passive Daytime Radiative Cooling of Silica Aerogels
Next Article in Special Issue
MoS2 Nanoplatelets on Hybrid Core-Shell (HyCoS) AuPd NPs for Hybrid SERS Platform for Detection of R6G
Previous Article in Journal
Finite Element Analysis Model of Electronic Skin Based on Surface Acoustic Wave Sensor
Previous Article in Special Issue
Strain Relaxation of InAs Quantum Dots on Misoriented InAlAs(111) Metamorphic Substrates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Masked Droplet Deposition for Site-Controlled Ga Droplets

1
Center for Hybrid Nanostructures (CHyN), University of Hamburg, Luruper Chaussee 149, D-22761 Hamburg, Germany
2
Department of Physics, Paderborn University, Warburger Str. 100, D-33098 Paderborn, Germany
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(3), 466; https://doi.org/10.3390/nano13030466
Submission received: 22 December 2022 / Revised: 11 January 2023 / Accepted: 18 January 2023 / Published: 23 January 2023

Abstract

:
Site-controlled Ga droplets on AlGaAs substrates are fabricated using area-selective deposition of Ga through apertures in a mask during molecular beam epitaxy (MBE). The Ga droplets can be crystallized into GaAs quantum dots using a crystallization step under As flux. In order to model the complex process, including the masked deposition of the droplets and a reduction of their number during a thermal annealing step, a multiscale kinetic Monte Carlo (mkMC) simulation of self-assembled Ga droplet formation on AlGaAs is expanded for area-selective deposition. The simulation has only two free model parameters: the activation energy for surface diffusion and the activation energy for thermal escape of adatoms from a droplet. Simulated droplet numbers within the opening of the aperture agree quantitatively with the experimental results down to the perfect site-control, with one droplet per aperture. However, the model parameters are different compared to those of the self-assembled droplet growth. We attribute this to the presence of the mask in close proximity to the surface, which modifies the local process temperature and the As background. This approach also explains the dependence of the model parameters on the size of the aperture.

1. Introduction

Secure information and communication technologies can greatly benefit from quantum effects wherein semiconductor quantum dots (QDs) represent essential building blocks as deterministic sources for single photons and entangled photon pairs [1,2,3]. However, the often-used self-assembled epitaxial QDs [4] are usually located on a semiconductor surface without significant lateral ordering. On the other side, the site-controlled generation of QDs on a surface is essential for the realization of advanced optical devices like quantum photonic integrated circuits [5]. Several approaches have been demonstrated for site-control QDs, often based on molecular beam epitaxy (MBE) technology [6]. We note that a true site control of the QD position requires more technical effort compared to a regular lateral ordering, which can be achieved by optimized self organization strategies. Two general concepts are applicable for MBE-based site control on a surface. First, modifications of the substrate surface before MBE growth cause a rearrangement of the planarly deposited material for site control [7,8,9,10,11,12,13,14,15,16,17,18,19,20]. Second, the material beam fluxes to the substrate surface are modified for a area-selective deposition [21,22,23]. We consider the second approach to be more flexible, since a patterned substrate often only allows geometries where the distances are close to the natural diffusion length of the unpatterned surface. On the other hand, a modulation of the beam fluxes e.g., by shadow masks, requires additional equipment which must be compatible with the restrictive requirements of the MBE method regarding the ultra-high purity of the used environment.
The present manuscript discusses the generation of site-controlled Ga droplets by area-selective deposition through apertures in a mask. Here, a shadow mask restricts the arrival of new material to areas below the apertures. Figure 1 illustrates the difference between self-assembled formation of droplets without lateral order and area-selective deposition through apertures in a mask for site-controlled droplets. The design of the mask is described in ref. [24] and the process in ref. [23]. For QD generation, the Ga droplets can be crystallized into GaAs quantum dots using a crystallization step under As flux (droplet epitaxy [25]). The central topic of this manuscript is a multiscale kinetic Monte Carlo (mkMC) simulation of the area-selective formation of Ga droplets by deposition through an aperture. The simulated surface morphologies are compared with experimental results for a parameterization of the model and for a better understanding of the complex mechanism behind the area-selective deposition and the subsequent thermal annealing step.

2. Masked Droplet Deposition

The present experimental approach for the generation of site-controlled GaAs QDs by area-selective droplet deposition is described in earlier publications [23,24]. Therefore, we give here only a brief summary of the major points.
The samples were fabricated in a multi-chamber solid-source MBE system using semi-insulating (100) GaAs wafers as substrates [23]. In a first step, MBE growth was performed without a mask, and a 100 nm thick Al 0.3 Ga 0.7 As layer was deposited at usual MBE process parameters. Then, the As flux was minimized by closing the As shutter and valve, the substrate temperature was reduced to 350 °C, and a Ga pre-coverage of 1.2 monolayers (ML) was deposited at a rate of 0.7 ML/s to create a Ga-terminated surface. After that, the substrate was cooled down to 100 °C and transferred into the preparation chamber. There, a mask–substrate sandwich was created under ultrahigh vacuum (UHV) conditions employing a special manipulator in the preparation chamber [23]. The mask is a nano-patterned 100 nm thick Si 3 N 4 -membrane on a Si(100) support wafer [24]. Electron beam lithography and reactive ion etching produced circular holes in the membrane of the mask with varied diameters d from 140 nm to 5 µm. We note that the masks are fully compatible with the demanding ultra-high vacuum requirements of the MBE technology. The mask is installed in close contact with the substrate surface. However, we expect that there is still a small gap between mask and substrate (probably <1 µm). After transferring the mask–substrate sandwich back to the growth chamber, different amounts of Ga were deposited now through the mask for area-selective deposition at 100 °C and minimized As flux. Then, after a 60 s pause, the mask–substrate sandwich was transferred into the preparation chamber. After mask removal, the substrate was transferred back to the growth chamber and annealed there for 5 min at a temperature of 400 °C and at minimized As flux.

3. Multiscale Kinetic Monte Carlo Simulation

The simulation model for the masked droplet deposition expands a previous model for the self-assembled formation of Al and Ga droplets on a GaAs surface [26]. Reference [26] compares two approaches to model droplet formation; one uses mean-field rate equations and the other a multiscale kinetic Monte Carlo (mkMC) simulation. Obviously, a mean field model is not compatible with the present masked deposition, where the material flux to the surface is laterally modulated. Therefore, we expand the mkMC model to now include the effects of the mask.
The mkMC simulation reflects an atomistic picture and approximates the substrate surface by a square simulation field with m = m x × m y surface sites. On an (001) AlGaAs surface, the distance between two surface sites is d s = a / 2 0.5 0.40 nm, with the lattice constant a. The size of the simulation field m = 5 d / d s is adjusted according to the respective diameter d of the aperture (see Section 4). In contrast to the model for self-assembled droplet formation, where cyclic boundary conditions are assumed [26], we assume here that atoms diffusing out of the simulation field will not come back and, thus, are deleted. Considered objects on the surface are mobile atoms (monomers) and hemispherically-shaped droplets composed of s > 1 atoms. Possible processes are the arrival of atoms from the vapor beam on the surface, surface diffusion of mobile monomers by nearest-neighbor hopping, nucleation events by collisions between migrating monomers, attachment of mobile monomers to droplets, and escape of atoms from droplets. As an extension of the original mkMC model, in addition, the coalescence of droplets is considered due to higher coverage with droplet material.
In the following, the modeling of the respective processes is described in more detail. The flux of atoms from the vapor beam to the surface is site dependent, with F within the circular opening in the mask; otherwise, F = 0. An atom arriving on the surface can result either in the formation of a new monomer, a nucleation event by a direct hit to another monomer, or in droplet growth by a direct hit. Monomers hop to a nearest-neighbor surface site with rate D = ν exp E S / k B T , where T is the temperature, ν = 2 k B T / h is a vibrational frequency [27], k B Boltzmann’s constant, h Planck’s constant, and E S the activation energy for surface diffusion. This diffusion of monomers can cause a site-change, a nucleation event, or the attachment to a droplet. In the latter case, the increasing droplet volume can cause the coalescence of two droplets. Here, the touching droplets merge to form a new droplet. The new size s is the sum of the two merging droplets, and the new position is the center of mass between them. After each droplet coalescence, the appearance of further coalescence events due to increasing droplet volume is checked recursively. Finally, atoms escape from a droplet composed of s atoms with rate R E , s = 2 π r ζ ν exp E E / k B T , where E E is the activation energy for escape of monomers from droplets, ζ = exp ( r c / r ) describes the enhancement of the vapor pressure for small droplets due to the Gibbs-Thomson effect, r = 3 s / ( 2 π ) 3 is the droplet radius, r c = 2 γ V m o l / ( N A k B T ) , γ is the surface tension (0.67 N/m for Ga), V m o l the molar volume (11.8 × 10 6 m 3 /mol for Ga), and N A the Avogadro constant. An escape event from a droplet yields a new monomer at a random angle and distance r + 2 from the droplet center.
In ref. [26], it is described that an individual monomer performs several orders of magnitude more diffusion events compared to arrival plus escape. In order to speed-up the simulation and reduce the number of simulation steps, we replace monomer diffusion via a large number of nearest-neighbor hops by fewer jumps over longer distances [26]. This approximation is established by DeVita et al. as a multiscale kinetic Monte Carlo algorithm [28]. In detail, all diffusion events within the time interval τ I = 1 / ( m F + s n s R E , s ) , up to the occurrence of the next arrival or escape, are summarized. The monomer surface diffusion length is λ = τ D according to the Einstein relation, with the diffusion time τ . In the time interval τ I the diffusing monomer travels a distance d I = τ I D . Or, in other words, the rate R I = 1 / τ I = m F + s n s R E , s for traveling a distance d I by diffusion is given by the time interval up to the next arrival or escape. If there are other objects at a distance smaller than d I , the diffusion can result in either a displacement of the monomer, in nucleation by collision with another monomer, or in attachment to a droplet. The probability p for a collision with another object depends on the circular segment r / ( π d ) covered by the object, where d is the distance and r is the radius of the object. This gives the nucleation rate R N , i j = r D / ( π d i j 3 ) at which monomer i collides with a second monomer j, with the distance d i j between both. Accordingly, the rate of attachment to a droplet k is R A , i k = r D / ( π d i k 3 ) . In the mkMC model, the various rates sum up to a total activity rate
R t o t = m F + s n s R E , s + i R I + j i R N , i j + k R A , i k
The rate R t o t is used for the random selection of the next process in the simulation and for the calculation of the time interval d t = 1 / R t o t up to the next process.

4. Results and Discussion

For the area-selective deposition, an amount of θ = 40 ML of Ga is deposited through apertures in a mask at a flux of 0.35 ML/s and a low temperature, nominally of T = 100 °C (see Section 2). We note that measurements of such low temperatures are not very precise. Apertures of varied diameter d = 140 nm, 230 nm, 340 nm, and 600 nm are used. The following thermal annealing step of 5 min is performed at T = 400 °C. The average number of Ga droplets in the area of a mask aperture is taken from scanning electron microscopy (SEM) images. Examples are shown in Figure 2 where typically, data from 20 apertures are used to determine the average droplet number. The data in Table 1 clearly demonstrate that small apertures and a thermal annealing step are essential to achieve perfect site control with one droplet per aperture.
The average droplet number is also obtained from mkMC simulated surface morphologies and compared with the experimental values. Here, the size of the simulation field m x = m y = 5 d / d s is adjusted according to the respective diameter d of the aperture, with the distance d s between two surface lattice sites. For the two activation energies in the mkMC model, we start with values obtained in previous simulations [26] for self-assembled Ga droplet nucleation without a mask. For T 300 °C the values are E S = 0.115 eV for monomer diffusion and E E = 1.24 eV for the escape of atoms from droplets. For T > 300 °C, the value of E E = 1.24 eV + 0.06 (T[°C] − 300)/100 eV is temperature dependent, which indicates that here, additional processes that modify the binding energy of adatoms to the droplets become relevant.

4.1. Simulation of the Area-Selective Deposition

For the area-selective deposition step at θ = 40 ML and T = 100 °C, simulations using the parameters for self-assembled Ga droplet formation ( E S = 0.115 eV, E E = 1.24 eV) yield only a poor reproduction of the experimental droplet numbers in Table 1. To explain this discrepancy, we assume that the presence of a mask at a very small distance to the surface and with only small apertures modifies the local process parameters in the experiments. This can be the substrate temperature, the As background, or a combination of both. Furthermore, the usage of a different MBE chamber in the present experiments can also modify the As background. A modification of the As background is expected to influence the activation energies E S and E E in the simulation. To evaluate the respective influence, we have performed numerous simulation runs.
Figure 3 shows pairs of E E , E S where the simulated droplet numbers agree with the experimental values after the deposition step. Two temperatures 50 °C and 100 °C as well as three aperture diameters d = 240 nm, 340 nm, and d = 600 nm are considered. As a general trend, a reduced E S is compensated by an increase in E E . However, the data in Figure 3a indicate a saturation for E E above about 1.45 eV, where the escape of atoms from droplets becomes negligible small. Due to the saturation, there is no matching value of E E for the value of E S = 0.115 eV taken from the self-assembled droplet formation. For T = 100 °C, there are fewer simulated values, since the computation time increases exponentially with the temperature. Therefore, a linear regression is used to extrapolate for a smaller E E .
To reduce the number of free model parameters, we assume in the following that only E S is sensitive to the mask, whereas E E = 1.24 eV agrees with the value of the self-assembled process. This approach reflects the known influence of an As background on the average activation energy for surface diffusion during GaAs epitaxy [29]. Here, a higher E S corresponds to a higher As background. In the present experiments, the mask can induce a higher As background and, thus, an increasing E S at a smaller d. The second column in Table 2 gives values of E S which provide agreement with the experiments as a function of d for constant E E = 1.24 eV and T = 100 °C.
In a second approach to explaining the influence of the aperture size, the local temperature T can depend on d. Figure 4 shows pairs of E E and E S with agreement between simulation and experiment for adjusted values of T. Interestingly, again due to a saturation, the value of E S is always above that for the self-assembled process. This indicates an elevated As background even in the case of a d-dependent temperature. Assuming a constant E E = 1.24 eV and E S = 0.46 eV, the adjusted temperatures are given in the third column of Table 2. A comparison between experimental droplet morphologies and the simulation results obtained using this approach are shown in Figure 2.
To summarize this part, the area-selective deposition related results indicate that the activation energy E S for surface diffusion is higher compared to the self-assembled droplet formation. This is probably related to a higher As background which can be caused by the usage of a different MBE chamber or by the presence of the mask. The reason for the additional dependence of the droplet number on the aperture size is not unequivocal, and the effect can be explained by an aperture-size dependence either of the temperature or of the As background and, thus, of E S and possibly also of E E .

4.2. Simulation of the Thermal Annealing Step

As a final process step, the samples are annealed without a mask for 5 min at T = 400 °C. This step reduces the droplet number substantially, as is visible in Table 1. Figure 5 shows examples of area-selective Ga droplets after annealing. Clearly visible is the strong influence of the aperture size on the droplet number.
The simulations of the annealing step start from the aperture-size dependent droplet morphologies which are described in Figure 2. Only the activation energy for the escape of atoms from droplets is modified at the high temperature, according to E E = 1.24 eV + 0.06 (T[°C] − 300)/100 eV (see Section 3). A comparison between experimental and simulated droplet morphologies in Figure 5 indicates that the model nicely reproduces the experimental trend. However, the simulation produces slightly higher droplet numbers. This deviation may be caused by the idealized model assumptions. As described in Section 3, the model considers an abrupt border between the area below the aperture with a constant impinging flux F and the area shadowed by the masked with F = 0. In the experiments, the geometry of the evaporation cell in combination with the distance between mask and substrate will lead to an additional transition region with a nonuniform flux and an enlarged diameter of the deposition area. Furthermore, the experiments indicate that a significant fraction of the deposited Ga is not found in the droplets (50% and more) [23]. This loss of the droplet material is attributed to an outdiffusion under the area shadowed by the mask. The simulations also indicate a material loss, but only up to 10%, which agrees with the higher droplet numbers in comparison to the experiments.
For a demonstration of perfect site control with one droplet per aperture, a further sample was fabricated using a smaller aperture diameter of d = 140 nm. The other parameters are as in Figure 2 and Figure 5. Figure 6 demonstrates the creation of a single site-controlled Ga droplet using this small aperture. In the simulations, a temperature of T = 50 °C is assumed during deposition, which reflects the value for d = 230 nm (see Section 4.1). The diameter of the droplet taken from SEM is 100 nm, and the simulated droplet (Figure 6) has a very similar diameter of 98 nm.
In order to study the annealing process in more detail, we have performed simulation runs with varied annealing temperatures T and annealing times t a . Results are plotted in Figure 7. Clearly visible is the decreasing droplet number at a higher T and a longer annealing time. This reduction is attributed to two effects. First, the droplet coarsening by Qstwald ripening is well known to reduce the droplet number [30]. Second, the diffusion of material to areas below the mask far away from the aperture will also reduce the droplet number. We note that these results confirm the choice of the experimental annealing parameters T = 400 °C, t a = 300 s, and d 230 nm as useful for the realization of perfect site control with one Ga droplet per aperture. Using such parameters, the fabrication of a perfectly ordered array of Ga droplets by deposition through a mask was demonstrated [23].

5. Conclusions

A multiscale kinetice Monte Carlo simulation [26] of self-assembled Ga or Al droplet formation on AlGaAs surfaces is expanded for area-selected deposition. Model results are compared with experimental data on the area-selective deposition [23] of Ga droplets through apertures in a mask [24]. In the experiments, a thermal annealing step after the area-selective deposition is very important for the reduction of droplet density down to one site-controlled droplet per aperture. The model is able to quantitatively reproduce the experimental behavior, including the area-selective deposition and the annealing step. However, the model parameters need to be changed in comparison to the self-assembled droplet growth. We assume here that the mask in close proximity to the surface modifies the local process temperature and the As background. A higher As background is expected to yield a higher activation energy for surface diffusion [29]. This approach also explains the dependence of the model parameters on the size of the aperture.
For a functionalization, the studied site-controlled Ga droplets can be crystallized into GaAs quantum dots under an arsenic flux. This method is established as droplet epitaxy (DE) [25,31]. We assume that the position and number of the Ga droplets will not change during crystallization, so that the site-control of the initial droplets is also maintained for the GaAs quantum dots. Droplet epitaxial quantum dots are promising building blocks for quantum photonic devices and quantum networks [25]. In an alternative approach, we expect that the masked deposition technique can be used also for the fabrication of site-controlled Al droplets. Al droplets allow the generation of various types of quantum structures through the local droplet etching (LDE) technique [32,33]. LDE takes place at much higher temperatures compared to DE, where the gradient of the As concentration causes a substantial diffusion of As from the AlGaAs substrate into the Al droplets. This yields the formation of nanoholes in the AlGaAs surface, which can be filled, e.g., with GaAs for the generation of quantum structures. Both concepts, DE and LDE, would be greatly improved by a site control of the resulting quantum structures. We expect that here, the present simulation model will be helpful for optimizing the process parameters.

Author Contributions

Conceptualization, C.H. and D.R.; methodology, S.F., V.Z. and A.A.; software, S.F. and C.H.; validation, C.H. and D.R.; investigation, S.F., V.Z. and A.A.; resources, C.H. and D.R.; data curation, S.F., V.Z. and A.A.; writing—original draft preparation, C.H.; writing—review and editing, C.H., D.R. and V.Z.; supervision, C.H. and D.R.; project administration, C.H. adn D.R.; funding acquisition, C.H. and D.R. All authors have read and agreed to the published version of the manuscript.

Funding

S.F., A.A. and C.H. received funding from the DFG via HE 2466/2-1 and from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 721394. V.Z. and D.R. received funding from the DFG via the collaborative research center SFB-TRR 142 and from the BMBF via the Q.com project No. 16KISO114.

Data Availability Statement

Not applicable.

Acknowledgments

S.F., A.A. and C.H. would like to thank the DFG for financial support via via HE 2466/2-1 and the European Union’s Horizon 2020 research and innovation programme for support under the Marie Skłodowska-Curie grant agreement No 721394. V.Z. and D.R. would like to thank the the DFG for financial support via the collaborative research center SFB-TRR 142 and the BMBF for financial support via the Q.com project No. 16KISO114.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nowak, A.K.; Portalupi, S.L.; Giesz, V.; Gazzano, O.; Savio, C.D.; Braun, P.F.; Karrai, K.; Arnold, C.; Lanco, L.; Sagnes, I.; et al. Deterministic and electrically tunable bright single-photon source. Nat. Commun. 2014, 5, 3240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Müller, M.; Bounouar, S.; Jöns, K.D.; Glässl, M.; Michler, P. On-demand generation of indistinguishable polarization-entangled photon pairs. Nat. Photonics 2014, 8, 224–228. [Google Scholar] [CrossRef] [Green Version]
  3. Huber, D.; Reindl, M.; Huo, Y.; Huang, H.; Wildmann, J.S.; Schmidt, O.G.; Rastelli, A.; Trotta, R. Highly indistinguishable and strongly entangled photons from symmetric GaAs quantum dots. Nat. Commun. 2017, 8, 15506. [Google Scholar] [CrossRef] [Green Version]
  4. Henini, M. (Ed.) Handbook of Self Assembled Semiconductor Nanostructures for Novel Devices in Photonics and Electronics, 1st ed.; Elsevier Science: Amsterdam, The Netherlands, 2008. [Google Scholar]
  5. Khasminskaya, S.; Pyatkov, F.; Słowik, K.; Ferrari, S.; Kahl, O.; Kovalyuk, V.; Rath, P.; Vetter, A.; Hennrich, F.; Kappes, M.M.; et al. Fully integrated quantum photonic circuit with an electrically driven light source. Nat. Photonics 2016, 10, 727–732. [Google Scholar] [CrossRef]
  6. Schmidt, O. Lateral Alignment of Epitaxial Quantum Dots; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  7. Kohmoto, S.; Nakamura, H.; Ishikawa, T.; Nishikawa, S.; Nishimura, T.; Asakawa, K. Site-controlled self-organization of InAs quantum dots. Mater. Sci. Eng. B 2002, 88, 292–297. [Google Scholar] [CrossRef]
  8. Mehta, M. Focused ion beam implantation induced site-selective growth of InAs quantum dots. Appl. Phys. Lett. 2007, 91, 123108. [Google Scholar] [CrossRef]
  9. Atkinson, P.; Schmidt, O.G.; Bremner, S.P.; Ritchie, D.A. Formation and ordering of epitaxial quantum dots. Comptes Rendus Phys. 2008, 9, 788–803. [Google Scholar] [CrossRef]
  10. Sünner, T.; Schneider, C.; Strauß, M.; Huggenberger, A.; Wiener, D.; Höfling, S.; Kamp, M.; Forchel, A. Scalable fabrication of optical resonators with embedded site-controlled quantum dots. Opt. Lett. 2008, 33, 1759–1761. [Google Scholar] [CrossRef]
  11. Martín-Sánchez, J.; Alonso-González, P.; Herranz, J.; González, Y.; González, L. Site-controlled lateral arrangements of InAs quantum dots grown on GaAs(001) patterned substrates by atomic force microscopy local oxidation nanolithography. Nanotechnology 2009, 20, 125302. [Google Scholar] [CrossRef]
  12. Schneider, C.; Huggenberger, A.; Sünner, T.; Heindel, T.; Strauß, M.; Göpfert, S.; Weinmann, P.; Reitzenstein, S.; Worschech, L.; Kamp, M.; et al. Single site-controlled In(Ga)As/GaAs quantum dots: Growth, properties and device integration. Nanotechnology 2009, 20, 434012. [Google Scholar] [CrossRef]
  13. Hakkarainen, T.V.; Tommila, J.; Schramm, A.; Tukiainen, A.; Ahorinta, R.; Dumitrescu, M.; Guina, M. Structural characterization of InAs quantum dot chains grown by molecular beam epitaxy on nanoimprint lithography patterned GaAs(100). Nanotechnology 2011, 22, 295604. [Google Scholar] [CrossRef]
  14. Skiba-Szymanska, J.; Jamil, A.; Farrer, I.; Ward, M.B.; Nicoll, C.A.; Ellis, D.J.P.; Griffiths, J.P.; Anderson, D.; Jones, G.A.C.; Ritchie, D.A.; et al. Narrow emission linewidths of positioned InAs quantum dots grown on pre-patterned GaAs(100) substrates. Nanotechnology 2011, 22, 065302. [Google Scholar] [CrossRef] [PubMed]
  15. Unrau, W.; Quandt, D.; Schulze, J.H.; Heindel, T.; Germann, T.D.; Hitzemann, O.; Bimberg, D. Electrically driven single photon source based on a site-controlled quantum dot with self-aligned current injection. Appl. Phys. Lett. 2012, 101, 211119. [Google Scholar] [CrossRef]
  16. Jöns, K.D.; Atkinson, P.; Müller, M.; Heldmaier, M.; Ulrich, S.M.; Schmidt, O.G.; Michler, P. Triggered Indistinguishable Single Photons with Narrow Line Widths from Site-Controlled Quantum Dots. Nano Lett. 2013, 13, 126–130. [Google Scholar] [CrossRef]
  17. Jamil, A. On-chip generation and guiding of quantum light from a site-controlled quantum dot. Appl. Phys. Lett. 2014, 104, 101108. [Google Scholar] [CrossRef] [Green Version]
  18. Sala, E.M.; Bollani, M.; Bietti, S.; Fedorov, A.; Esposito, L.; Sanguinetti, S. Ordered array of Ga droplets on GaAs(001) by local anodic oxidation. J. Vac. Sci. Technol. Nanotechnol. Microelectron. Mater. Process. Meas. Phenom. 2014, 32, 061206. [Google Scholar] [CrossRef]
  19. Herranz, J.; Wewior, L.; Alén, B.; Fuster, D.; González, L.; González, Y. Role of re-growth interface preparation process for spectral line-width reduction of single InAs site-controlled quantum dots. Nanotechnology 2015, 26, 195301. [Google Scholar] [CrossRef] [Green Version]
  20. Yang, L.; Carter, S.G.; Bracker, A.S.; Yakes, M.K.; Kim, M.; Kim, C.S.; Vora, P.M.; Gammon, D. Optical spectroscopy of site-controlled quantum dots in a Schottky diode. Appl. Phys. Lett. 2016, 108, 233102. [Google Scholar] [CrossRef]
  21. Ohkouchi, S.; Nakamura, Y.; Nakamura, H.; Asakawa, K. Indium nano-dot arrays formed by field-induced deposition with a Nano-Jet Probe for site-controlled InAs/GaAs quantum dots. Thin Solid Film. 2004, 464, 233–236. [Google Scholar] [CrossRef]
  22. Ohkouchi, S.; Nakamura, Y.; Ikeda, N.; Sugimoto, Y.; Asakawa, K. In situ mask designed for selective growth of InAs quantum dots in narrow regions developed for molecular beam epitaxy system. Rev. Sci. Instrum. 2007, 78, 073908. [Google Scholar] [CrossRef] [Green Version]
  23. Zolatanosha, V.; Reuter, D. Site-controlled droplet epitaxy of GaAs quantum dots by deposition through shadow masks. J. Vac. Sci. Technol. Nanotechnol. Microelectron. Mater. Process. Meas. Phenom. 2018, 36, 02D105. [Google Scholar] [CrossRef]
  24. Zolatanosha, V.; Reuter, D. Robust Si3N4 masks for 100nm selective area epitaxy of GaAs-based nanostructures. Microelectron. Eng. 2017, 180, 35–39. [Google Scholar] [CrossRef]
  25. Gurioli, M.; Wang, Z.; Rastelli, A.; Kuroda, T.; Sanguinetti, S. Droplet epitaxy of semiconductor nanostructures for quantum photonic devices. Nat. Mater. 2019, 18, 799–810. [Google Scholar] [CrossRef] [PubMed]
  26. Heyn, C.; Feddersen, S. Modeling of Al and Ga Droplet Nucleation during Droplet Epitaxy or Droplet Etching. Nanomaterials 2021, 11, 468. [Google Scholar] [CrossRef]
  27. Clarke, S.; Vvedensky, D.D. Origin of Reflection High-Energy Electron-Diffraction Intensity Oscillations during Molecular-Beam Epitaxy: A Computational Modeling Approach. Phys. Rev. Lett. 1987, 58, 2235–2238. [Google Scholar] [CrossRef] [PubMed]
  28. DeVita, J.P.; Sander, L.M.; Smereka, P. Multiscale kinetic Monte Carlo algorithm for simulating epitaxial growth. Phys. Rev. B 2005, 72, 205421. [Google Scholar] [CrossRef] [Green Version]
  29. Heyn, C.; Harsdorff, M. Simulation of GaAs growth and surface recovery with respect to gallium and arsenic surface kinetics. Phys. Rev. 1997, 55, 7034–7038. [Google Scholar] [CrossRef]
  30. Heyn, C.; Stemmann, A.; Hansen, W. Dynamics of self-assembled droplet etching. Appl. Phys. Lett. 2009, 95, 173110–173112. [Google Scholar] [CrossRef]
  31. Abbarchi, M.; Mastrandrea, C.A.; Kuroda, T.; Mano, T.; Sakoda, K.; Koguchi, N.; Sanguinetti, S.; Vinattieri, A.; Gurioli, M. Exciton fine structure in strain-free GaAs/Al0.3Ga0.7As quantum dots: Extrinsic effects. Phys. Rev. B 2008, 78, 125321–125324. [Google Scholar] [CrossRef]
  32. Heyn, C.; Stemmann, A.; Köppen, T.; Strelow, C.; Kipp, T.; Grave, M.; Mendach, S.; Hansen, W. Highly uniform and strain-free GaAs quantum dots fabricated by filling of self-assembled nanoholes. Appl. Phys. Lett. 2009, 94, 183113–183115. [Google Scholar] [CrossRef]
  33. Heyn, C.; Küster, A.; Ungeheuer, A.; Gräfenstein, A.; Hansen, W. Excited-state indirect excitons in GaAs quantum dot molecules. Phys. Rev. B 2017, 96, 085408. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of self-assembled droplet formation with planar deposition of the droplet material on a substrate, random nucleation, and formation of metal droplets without lateral order. (b) Schematic of an area-selective deposition process where Ga droplets are site controlled on an Al 0.3 Ga 0.7 As surface by apertures in a mask which consists of a 100 nm thick Si 3 N 4 -membrane on a Si(100) support wafer.
Figure 1. (a) Schematic of self-assembled droplet formation with planar deposition of the droplet material on a substrate, random nucleation, and formation of metal droplets without lateral order. (b) Schematic of an area-selective deposition process where Ga droplets are site controlled on an Al 0.3 Ga 0.7 As surface by apertures in a mask which consists of a 100 nm thick Si 3 N 4 -membrane on a Si(100) support wafer.
Nanomaterials 13 00466 g001
Figure 2. SEM images (left) of Ga droplets after area-selective deposition without annealing for varied diameters d of the mask aperture together with corresponding simulated surface morphologies (right). The red circles in the simulated surfaces indicate the size of the aperture. The experimental parameters are F = 0.35 ML/s, θ = 40 ML, and T = 100 °C. The simulations are performed using E S = 0.46 eV, E E = 1.24 eV, as well as an aperture size-dependent temperature of 50 °C at d = 230 nm, 65 °C at d = 340 nm, and 100 °C at d = 600 nm (see text).
Figure 2. SEM images (left) of Ga droplets after area-selective deposition without annealing for varied diameters d of the mask aperture together with corresponding simulated surface morphologies (right). The red circles in the simulated surfaces indicate the size of the aperture. The experimental parameters are F = 0.35 ML/s, θ = 40 ML, and T = 100 °C. The simulations are performed using E S = 0.46 eV, E E = 1.24 eV, as well as an aperture size-dependent temperature of 50 °C at d = 230 nm, 65 °C at d = 340 nm, and 100 °C at d = 600 nm (see text).
Nanomaterials 13 00466 g002
Figure 3. To parameterize the area-selective deposition step, pairs of E E , E S are determined which yield simulated droplet numbers that are in agreement with the experimental values for (a) T = 50 °C and (b) T = 100 °C. Three aperture diameters d are considered as indicated. The continuous lines in (b) show results of linear fits. The dashed blue lines indicate the activation energies obtained for self-assembled droplet formation.
Figure 3. To parameterize the area-selective deposition step, pairs of E E , E S are determined which yield simulated droplet numbers that are in agreement with the experimental values for (a) T = 50 °C and (b) T = 100 °C. Three aperture diameters d are considered as indicated. The continuous lines in (b) show results of linear fits. The dashed blue lines indicate the activation energies obtained for self-assembled droplet formation.
Nanomaterials 13 00466 g003
Figure 4. Pairs of E E and E S which yield simulated droplet numbers in agreement with the experimental values after deposition. The temperature is assumed to depend on the aperture diameter d as indicated. The dashed blue lines indicate the activation energies obtained for self-assembled droplet formation.
Figure 4. Pairs of E E and E S which yield simulated droplet numbers in agreement with the experimental values after deposition. The temperature is assumed to depend on the aperture diameter d as indicated. The dashed blue lines indicate the activation energies obtained for self-assembled droplet formation.
Nanomaterials 13 00466 g004
Figure 5. SEM images (left) of Ga droplets after area-selective deposition and thermal annealing for varied diameters d of the mask aperture, together with corresponding simulated surface morphologies (right). The red circles indicate the size of the aperture. The experimental and simulation parameters for the area-selective deposition step are as in Figure 2. Thermal annealing is performed for 5 min at T = 400 °C. The simulations of the annealing step use a T-dependent E E , as is described in the text.
Figure 5. SEM images (left) of Ga droplets after area-selective deposition and thermal annealing for varied diameters d of the mask aperture, together with corresponding simulated surface morphologies (right). The red circles indicate the size of the aperture. The experimental and simulation parameters for the area-selective deposition step are as in Figure 2. Thermal annealing is performed for 5 min at T = 400 °C. The simulations of the annealing step use a T-dependent E E , as is described in the text.
Nanomaterials 13 00466 g005
Figure 6. SEM image (left) of a single site-controlled Ga droplet for d = 140 nm together with a corresponding simulated surface morphology (right). The red circles indicate the size of the aperture. The experimental and simulation parameters are as in Figure 5. In the simulations, during deposition T = 50 °C is assumed (see Section 4.1).
Figure 6. SEM image (left) of a single site-controlled Ga droplet for d = 140 nm together with a corresponding simulated surface morphology (right). The red circles indicate the size of the aperture. The experimental and simulation parameters are as in Figure 5. In the simulations, during deposition T = 50 °C is assumed (see Section 4.1).
Nanomaterials 13 00466 g006
Figure 7. Simulated droplet number after annealing (a) as a function of annealing temperature T and (b) as a function of the annealing time t a for different aperture diameters d.
Figure 7. Simulated droplet number after annealing (a) as a function of annealing temperature T and (b) as a function of the annealing time t a for different aperture diameters d.
Nanomaterials 13 00466 g007
Table 1. Experimental average number of Ga droplets in the area of mask apertures of varied diameters d, determined using SEM.
Table 1. Experimental average number of Ga droplets in the area of mask apertures of varied diameters d, determined using SEM.
d [nm]Droplets after DepositionDroplets after Annealing
23020.8 ± 1.221.1 ± 0.07
34034.3 ± 1.362.4 ± 0.21
60052.0 ± 1.9010.8 ± 0.41
Table 2. Two approaches for simulation parameters which reproduce the experimental droplet numbers after area-selective deposition for different aperture diameters d. E E = 1.24 eV is assumed to be constant.
Table 2. Two approaches for simulation parameters which reproduce the experimental droplet numbers after area-selective deposition for different aperture diameters d. E E = 1.24 eV is assumed to be constant.
d [nm] E S [eV] at T = 100 °CT [°C] at E S = 0.46 eV
2300.6650
3400.5965
6000.46100
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Feddersen, S.; Zolatanosha, V.; Alshaikh, A.; Reuter, D.; Heyn, C. Modeling of Masked Droplet Deposition for Site-Controlled Ga Droplets. Nanomaterials 2023, 13, 466. https://doi.org/10.3390/nano13030466

AMA Style

Feddersen S, Zolatanosha V, Alshaikh A, Reuter D, Heyn C. Modeling of Masked Droplet Deposition for Site-Controlled Ga Droplets. Nanomaterials. 2023; 13(3):466. https://doi.org/10.3390/nano13030466

Chicago/Turabian Style

Feddersen, Stefan, Viktoryia Zolatanosha, Ahmed Alshaikh, Dirk Reuter, and Christian Heyn. 2023. "Modeling of Masked Droplet Deposition for Site-Controlled Ga Droplets" Nanomaterials 13, no. 3: 466. https://doi.org/10.3390/nano13030466

APA Style

Feddersen, S., Zolatanosha, V., Alshaikh, A., Reuter, D., & Heyn, C. (2023). Modeling of Masked Droplet Deposition for Site-Controlled Ga Droplets. Nanomaterials, 13(3), 466. https://doi.org/10.3390/nano13030466

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop