Abstract

Single-crystalline InGaAs/GaAs semiconductor micro/nanotubes have been obtained by the strain-driven self-rolling mechanism. This approach combines the advantages of bottom-up (epitaxial growth) and top-down (postgrowth processing) techniques, offering an exceptional opportunity to realize complex three-dimensional nanoarchitectures by using conventional photolithography and wet-etching processes. The method employed to obtain micro/nanotubes with selected orientation and length is described in detail. By means of high-resolution scanning electron microscopy characterization, we show a clear shape difference between single-wall and multiwalls tubes and we discuss it on the basis of strain release, taking into account also possible shape deformations induced during micro/nanotubes drying. We analyse the In-segregation profile in the nominal In0.20Ga0.80As/GaAs bilayer and we show its effect on the actual diameter of the tubes, concluding that a more accurate description of the structure should consider an In0.20Ga0.80As/In0.10Ga0.90As/GaAs trilayer. This work will be useful to set up reliable methodologies for the realization of strain-driven micro/nanotubes with controlled properties, necessary for their implementation in a large number of application fields.

1. Introduction

The miniaturization of components and devices is nowadays an essential requirement for almost all application fields. Scaling down of current technologies to obtain a real three-dimensional (3D) definition at the nanoscale is complicated and expensive, conditions that strongly limit the achievement of mass-production. Almost two decades ago, Prinz and coworkers [1] proposed a method to define micro/nanoarchitectures through the controlled release of strained membranes [2]. Since then, several authors studied and implemented this technique, highlighting the possibility of obtaining miniaturized devices without the need of complex procedures [35]. This technique combines bottom-up and top-down methodologies. Bottom-up methods are used to create a two-dimensional (2D) system containing a certain amount of strain energy, as in the epitaxial growth of mismatched heterostructures. The top-down method consists of a photolithographic technique used to define a pattern on the structure. The strain of the system can be released by removing a sacrificial layer using wet-etching procedures. This results in a shape transition of the freestanding 2D system into a 3D structure.

Several 3D shapes have been reported [4] pointing out the unique feature of this approach: the combination of nanoscale dimensions given by the epitaxial growth and microscale size defined by lithography. One of the most studied 3D shape is based on tube formation as a result of strain-driven rolling-up of 2D membranes [3]. These micro/nanotubes (μtubes) have been proposed for several applications: tubular-based photonic integrated components (filters, modulators, directional couplers, and lasers) obtained by embedding active layers such as quantum dots and quantum wells in the heterostructure [6], novel sensors directly integrated into microfluidic channels (lab-in-a-tube) [7], and also energy-storage microdevices [8].

A deeper understanding of the relationship between structural parameters of the starting 2D system and the properties of strain-driven 3D structures is clearly needed to implement such applications. To this aim, we focus on the study of μtubes formed starting from InGaAs/GaAs bilayers grown by Molecular Beam Epitaxy (MBE). MBE is a technique that allows a high degree of control on the deposition of heterostructures and, then, is particularly suited for this work [9]. The semiconductor structure is patterned using conventional optical lithography and the bilayer is released by a wet-etching. This last step allows the study of the very first stages of the rolling-up process which, as we observe, shows deviations from the ideal rolled-up μtube. The actual structure of the bilayer is investigated through the comparison of photoluminescence (PL) results and energy levels calculation of the as-grown heterostructure, based on structure parameters derived by X-ray diffraction (XRD). The results of this study indicate the presence of an In-segregation profile, typical of In-containing MBE-grown structures [10]. Due to this effect, we conclude that the actual structure should be better described by an In0.20Ga0.80As/In0.10Ga0.90As/GaAs trilayer, rather than the nominal In0.20Ga0.80As/GaAs bilayer. Using the trilayer parameters, we calculate the μtube diameter and we find a substantial agreement with the experimental diameter measured by scanning electron microscopy (SEM).

Our study highlights the relevance of 2D structure deposition procedures and parameters to control the features of the final 3D architectures. This is of interest to correctly design rolled-up micro/nanostructures with desired properties for direct implementation into practical devices.

2. Materials and Methods

The growth of the InGaAs/GaAs bilayer structure was carried out in a solid-source Varian Gen II MBE system, on (100)-oriented GaAs semi-insulating substrate (see Figure 1(a)). After the desorption of the native oxide from the GaAs substrate surface, a 5 nm thick AlAs buffer layer was grown at 600°C; the thickness of AlAs was chosen to optimize the postgrowth processing of the structure [11]. The substrate temperature was then lowered to 500°C during a growth interruption of 60 s under Arsenic flux to grow the InGaAs/GaAs bilayer. This lower temperature is instrumental to minimize Indium-segregation and desorption during InGaAs deposition and alloy intermixing at InGaAs-GaAs interface. The InGaAs thickness is well below the critical value predicted by the Matthews–Blakeslee model for the formation of misfit dislocations [12]. The evolution of the reflection high-energy electron diffraction (RHEED) pattern during MBE deposition indicated that the InGaAs was pseudomorphic to GaAs, thus assuring a high crystal quality of the structure. The strain state of InGaAs/GaAs bilayer was confirmed by double crystal XRD measurements in conjunction with detailed rocking curve analysis. By XRD analysis we were able to obtain a reliable estimation of the thicknesses and compositions of layers in the structure [13], as discussed in Results and Discussion.

Continuous-wave photoluminescence (PL) was performed at 10 K on the as-grown structure using the 532 nm line of a DPSS laser as excitation, with a maximum power density of about 5 W/cm2. The PL spectra were measured by a Fast-Fourier Transform spectrometer with a 1 meV resolution and a cooled Ge detector.

XRD values of thicknesses and compositions were used to calculate the quantum confined energy levels of the planar structure (Figure 1(a)), in order to compare the expected emission energies with PL characterization. Calculation of the quantum system composed by the AlAs, InGaAs, and GaAs layers was performed by means of the software Tibercad, following the procedure described in depth in [14]. We considered a conduction-band offset of 0.52 (inclusive of strain effects) and a value of the exciton binding energies of 10 meV [15].

The design parameters of the epitaxial bilayer determine the μtube diameters and the μtube nanowall thicknesses, since these features are controlled by the value of the lattice mismatch and by the relative thickness of InGaAs and GaAs layers. On the contrary, the final shape and the positioning of 3D microstructures on the substrate are controlled by the photolithographic patterning. To define the pattern, we used conventional photolithography and two consecutive wet-chemical etching processes with different etching solutions.

First, on the surface of epitaxial structure we prepared a positive photoresist micropattern consisting of an array of rectangular shapes with the same width (40 μm) and different lengths, ranging from 50 to 3200 μm (50, 100, 200, 400, 800, 1600, and 3200 μm). The edges of rectangles were precisely aligned along the directions that are known to be the preferred rolling directions on (100) GaAs, as a consequence of Young’s modulus anisotropy in cubic semiconductor crystals [16, 17].

Then, we performed the first wet-chemical etching (H2SO4 : H2O2 : H2O = 1 : 8 : 80 for 12 s) to transfer the geometric pattern down to the AlAs sacrificial buffer layer. The etched profiles obtained were investigated by Atomic Force Microscopy (AFM). We measured a uniform etching depth of  nm, deep enough to uncover the AlAs surrounding the rectangular mesa structures.

Finally as shown in Figures 1(b) and 1(c), to release the bilayer from the substrate, AlAs sacrificial layer was underetched in a high selective diluted HF solution (HF : H2O = 1 : 2) [18]. The addition of a surfactant (benzalkonium chloride) and the gentle stirring of HF solution during the etching process were required to optimize the etching efficiency and to improve the spatial homogeneity of the 3D InGaAs/GaAs μtubes [19]. Both the wet-chemical etchings were performed at room temperature. The structures were rinsed in two steps: first in deionized (DI) water for 5 minutes and then in methanol for the following 5 minutes. The structures were dried on a hot-plate for 10 s at 70°C. The use of a low surface tension solvent (0.0227 N/m for methanol against 0.0728 N/m for water, both at 25°C) allows preventing the tubes from collapsing during the evaporation of the rinsing solution [20, 21]. In addition, we followed the work by Seleznev et al. [22] to calculate the collapsing pressure and the critical pressure (see and of [22], resp.) relative to our structures. When exceeds , the tubes are expected to collapse. Considering methanol as rinsing liquid, we obtain bar and = 1.39 bar, thus showing that our structures are not expected to collapse.

An optical microscope image of the μtube array obtained by longer etching time (45 s) is shown in Figure 2.

Field-Emission SEM was used to investigate the morphology of the μtubes in 54°-tilted configurations. The 54°-tilted images are automatically corrected for vertical parallax to allow length measurements along the growth direction. SEM images were obtained with 5 keV acceleration voltage of the electron beam and 5 mm working distance.

3. Results and Discussion

As shown by Chun et al. [23, 24], μtubes grown on GaAs (001) substrates preferentially roll-up along the directions. Additionally, the actual rolling direction of μtubes is determined by geometrical conditions. Given a rectangular mesa shape oriented along the directions, experimental results and theoretical calculations demonstrate that long-side rolling is energetically favoured. This is shown in the optical microscope image of Figure 2, where μtubes are oriented along the direction corresponding to the longer side of the rectangular mesas.

In order to study the features of the rolling process, we analysed the shape of μtubes as a function of the etching time () of the underlying sacrificial layer. By increasing , a larger fraction of the InGaAs/GaAs bilayer detaches from the underlying substrate and is free to bend. In Figure 3 we plot the rolling distances () experimentally measured from several SEM images of the μtubes as functions of the etching times. Assuming that = 0 when = 0, from the plot we extrapolate a very low value for the initial etching velocity (). This may stem from the existence of an oxidized AlAs layer that constitutes a barrier for the etching. As soon as the oxide is etched away, the bilayer starts detaching with a = 0.68 μm/s constant velocity proving that the procedure we chose grants a uniform etching of the sacrificial layer.

By reducing to a few tens of seconds, we can analyse the very first stages of rolling-up. With reference to Figure 4, we distinguish three regimes:(i) = 20 s, = μm: this is the initial stage of rolling-up where the InGaAs/GaAs freestanding layer starts bending (Figure 4(a));(ii) = 30 s, = μm: the bilayer rolls up into the tube. SEM magnification shows an irregular shape of the tube section (Figure 4(b));(iii) = 45 s, = μm: this stage corresponds to the formation of μtubes with several foldings. SEM image shows a perfectly circular shape of the tube section (Figure 4(c)).SEM images in Figure 4 show that only the μtubes obtained with a 45 s etching possess the expected circular sections while the structures resulting from shorter etching times show an irregular shape. Usually, shape deformations in single-wall μtubes are attributed to partial or complete collapse during the drying process [22]. In the previous section we showed that our specific conditions for μtube drying prevent the collapse of μtubes. Anyhow, the model proposed by Seleznev et al. [22] does not take into account capillary forces between the substrate and the walls that could induce an additional shape deformation. Our experimental study allows studying and comparing the very first stages of the rolling-up process and gives the possibility of accurately analysing different shape effects. While nonhomogeneous curvature radius of the bilayer in the proximity of the substrate can be influenced by capillary forces, Figures 4(a) and 4(b) show a change of curvature radius also in the proximity of the μtube free edge. A reason for this behaviour could be that an additional lateral strain is released in the proximity of the free edge of the bilayer, while this is not possible in the inner region of the freestanding bilayer. As a result, the curvature radius becomes larger as the distance from the edge decreases, as indeed observed in Figure 4(a). Figure 4(b), corresponding to a complete μtube, still maintains an irregular shape. As clearly visible from the SEM image, a strong curvature change of the μtube section takes place when the outer folding overlaps with the underlying free edge of the bilayer. By further increasing the number of foldings (Figure 4(c)), the shape becomes regular with a homogeneous curvature radius. We explain this behaviour due to the balance between different curvature radii of the overlapping bilayer: when the number of foldings increases, the free edge is flattened and is forced to assume the curvature of the inner regions of the bilayer. The dependence of the lattice parameter of the tubes on the proximity to the edge bound to the substrate has been shown by Krause et al. [25] through X-ray microdiffraction analysis of several types of tubes. This experimental evidence supports the idea of the existence of regions of the freestanding bilayer subjected to different strain.

On the = 45 s μtube (Figure 4(c)) we measured an internal diameter of 1.582 μm. Following the work by Deneke et al. [19] the diameter of rolled-up μtubes can be related to the thickness and composition of the materials that constitute the bilayer. As is known, the epitaxial growth of In-containing compounds on GaAs is affected by segregation, a phenomenon that may strongly change the effective composition profile of the structure [10]. Through X-ray measurements and calculations of the energy-band structure, we already showed that PL emission of buried In(Ga)As/GaAs quantum wells is strongly dependent on their effective In-content profile [15]. To take into account such effects, we analysed the PL emission of the as-grown structure and modeled the band structure in order to compare the experimental and calculated emission energies; these study is performed starting from X-ray measured thicknesses and compositions (Figure 5). As a result of the segregation, the nominal In0.2Ga0.8As (5 nm)/GaAs (11.5 nm) bilayer turns out to be better described by an In0.2Ga0.8As (3.3 nm)/In0.1Ga0.9As (4 nm)/GaAs (11.5 nm) trilayer. As a further confirmation, the thickness of the trilayer structure (18.8 nm) well compares to the thickness of the μtube walls measured from high-magnification SEM images ( nm), as shown in Figure 6.

In order to calculate the expected value of μtube diameter by using the actual parameters of the structure, we followed the approach proposed by Nishidate and Nikishkov in [26] for the calculation of generalized plane strain deformation of multilayer structures. On the basis of plane strain elasticity theory, authors derive curvature estimation for strained multilayer structure with an initial strain due to lattice mismatch. In our study the strained structure consists of three layers of thicknesses = 3.3 nm, = 4 nm and = 11.5 nm. The curvature of a trilayer can be calculated by setting in of [26], where values are defined as , , , and and , , and , with , , and given by the lattice constants of In0.2Ga0.8As, In0.1Ga0.9As, and GaAs, respectively. Following this notation, , , and are the initial strains in the lower-In0.2Ga0.8As, middle-In0.1Ga0.9As, and upper-GaAs layer. Exact values of Young modulus and Poisson ratio for the different materials were considered to solve the generalized multilayer plane strain equation using Matlab. A value of the μtube diameter of 1.764 μm was obtained, against the experimental value of μm obtained from the statistical analysis of μtubes as the one shown in Figure 4(c). The discrepancy between theoretical and experimental values could be due to additional strains on the structure surface and/or to the consumption of the top GaAs oxide during wet-chemical process in HF solution, as suggested by Deneke et al. [19]. The typical thickness of GaAs oxide in these structures is in the order of a few monolayers (ML) [27]. By taking into account a GaAs thickness reduction of 4 ML, we calculated a diameter of 1.608 μm in excellent agreement with the experimental value.

4. Conclusions

Here we have presented the fabrication of InGaAs/GaAs μtubes with nanometer-sized walls, thanks to the strain-induced self-rolling of MBE-grown bilayers. Lengths and location of μtubes on the substrate are controlled by mesa structures defined by conventional photolithographic processing, while micro/nanotube diameters depend on the lattice mismatch and on the thickness of the InGaAs/GaAs structure.

We have discussed in detail all the technological steps needed to obtain such 3D microstructures: (i) MBE growth of starting structure, consisting of layers of AlAs, GaAs, and pseudomorphic InGaAs, (ii) XRD and PL characterization of the epitaxial structure, (iii) photolithography, and (iv) wet-chemical etching.

SEM characterization of rolled-up structures has allowed us to understand the mechanisms of the self-rolling process. In particular, we observe an evident difference between single-wall and multiwalls tubes that we ascribe in part to capillary forces between the substrate and the μtube walls but also to a different strain status of the rolled-up bilayer. When the outer edge of the bilayer is free to elastically relax an increase in the curvature radius is observed in comparison with the situation when the edge is constrained to bend under other foldings of the bilayer.

Structural parameters of the epitaxial layer were obtained by XRD characterization and validated by comparing PL emission energies with model calculations of quantum levels related to the AlAs/InGaAs/GaAs quantum well. We showed that the structure is actually an In0.2Ga0.8As/In0.1Ga0.9As/GaAs trilayer, rather than an In0.2Ga0.8As/GaAs bilayer. This picture is confirmed by SEM high-magnification of walls and by theoretical calculation of curvature radius.

A trilayer is expected to have rather different properties than a bilayer such as band profiles, strain status, and refraction index mismatch with air. All of these are crucial parameters for applications of nanotubes as nanoelectronic elements, as systems for strain engineering of embedded nanostructures (such as quantum dots), as optical and chemical sensors.

Hence, we believe that the results we have shown here will contribute to a deeper understanding of the properties of these 3D nanostructures.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work has been partially supported by Fondazione Monte Parma. The authors acknowledge E. Buffagni, F. Rossi, and C. Ferrari for XRD characterization.