In Kim†
a,
Suji Choi†bc,
Ji-Hwan Kwonb,
Sang Jung Ahnb,
Min Sun Yeoma,
Ho Seong Leec,
Seong-Hoon Yic and
Young Heon Kim*d
aCenter for Supercomputing Applications, Korea Institute of Science and Technology Information, 245 Daehak-ro, Daejeon 34141, Republic of Korea
bKorea Research Institute of Standards and Science, 267 Gajeong-Ro, Yuseong-Gu, Daejeon 34113, Republic of Korea
cKyungpook National University, 80 Daehak-ro, Buk-gu, Daegu 41566, Republic of Korea
dGraduate School of Analytical Science and Technology (GRAST), Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, Republic of Korea. E-mail: y.h.kim@cnu.ac.kr
First published on 22nd December 2020
An in-depth understanding of thermal behavior and phase evolution is required to apply heterostructured nanowires (NWs) in real devices. The intermediate status during the vaporization process of InAs NWs in an Al2O3 shell was studied by conducting quenching during in situ heating experiments, using a transmission electron microscope. The formation of As clusters in the amorphous Al2O3 shell was confirmed by analyzing the high-angle annular dark field images and energy-dispersive X-ray spectra. The As clusters existed independently in the shell and were also observed at the end of the InAs pieces obtained after quenching. The formation process of the As clusters was demonstrated from a theoretical perspective. Moreover, an ab initio molecular dynamics simulation (AIMD) was conducted to study the atomic and molecular behaviors.
Thermal behavior and phase evolution in heterostructured nanomaterials have been studied recently using transmission electron microscopy (TEM).11–13 Recently well-established in situ heating TEM experiments can provide an alternative way to study the vaporization process for nanometer-sized materials by monitoring morphological and microstructural properties in real time with a high resolution of individual nanostructures. Moreover, the volumes of melting, vaporization, and condensation molecules can be quantified by considering the 3-dimensional morphology, and the exact sizes of these molecules, and the dynamic behaviors of the vaporization processes can be understood by analyzing the bright-field (BF) and high-resolution (HR) TEM images. Based on advanced TEM techniques, several studies have reported specific phase transition phenomena in nanomaterials by adopting various heterostructured nanomaterials.14–16
In this study, we report the observation of arsenic (As) clusters formed during the in situ heating TEM experiments of InAs NWs covered with an Al2O3 thin shell. Al2O3 has been widely used in shell structure engineering for producing nanotubes for drug delivery, and it has also been widely considered as a passivation layer in III-As compound semiconductors.15,16 In addition, an Al2O3 thin film imparts visibility to the phase evolution behavior of target materials during TEM imaging because of its amorphous structure. The atomistic evolution after the decomposition of indium and arsenic atoms is demonstrated by analyzing the results from molecular dynamics simulation.
We performed an ab initio molecular dynamics (AIMD) simulation to observe the dynamics of the surface atoms of InAs NWs. We adopted the projector-augmented wave method implemented in the Vienna ab initio simulation package.18 The flavor of Perdew–Berke–Ernzerhof exchange correlation function was adopted with a 400 eV plane wave basis energy cut-off.19 A slab model of In-terminated InAs(100) with a thickness of 20 Å (containing 64 atoms in total) was generated with a vacuum layer along the c-axis direction. The system was simulated based on the NVT ensemble for 40 ps at 1500 K using the Nosé-Hoover thermostat with 1.0 fs timestep. A set of 4 × 4 × 1 k-point grids was employed to sample the 1st Brillouin zone. To model the slab surface during the dynamics simulation, two In and As layers were fixed to mimic bulk states. Diffusion coefficients were calculated using TRAVIS, and the radial distribution functions were analyzed using VMD.20,21
Fig. 1 In situ TEM images during heat treatment indicating the morphological and microstructural evolutions of InAs/Al2O3 core/shell systems. (a) RT, (b) 600.0 °C, (c) 742.5 °C, and (d) 745.0 °C. |
Fig. 2 shows the HAADF-scanning TEM (STEM) images and the EDS results for the InAs/Al2O3 heterostructured NWs that remained localized after the completion of the in situ heating TEM experiments at 705 °C. In the HAADF-STEM image shown in Fig. 2(a), white bar-like structures appear to be present in the Al2O3 shell. However, the bar-like structure in region L1 in Fig. 2(a) is darker compared with other white bar-like structures. Moreover, it is slightly discontinuous and has non-uniform contrast. The EDS results obtained by the local elemental mapping method, which are depicted in Fig. 2(b)–(e), show the spatial distributions of the chemical elements In, As, O, and Al on the heterostructured NWs after the in situ heating TEM experiment. Region L1 was devoid of In atoms, although the other bar-like structures consisted of the chemical element In. On the other hand, the signal intensity for As atoms in region L1 was much stronger than the other areas. The EDS results indicated that region L1 was composed only of pure As atoms, whereas the other bar-like structures were composed of both In and As atoms. The signals for Al and O were intense in the shell area and almost same at the inside and outside of the shell.
Moreover, the shell area was completely devoid of In and As atoms. Therefore, we deduced that the chemical reaction between the InAs core and the Al2O3 shell did not occur during the in situ heating TEM experiment. Specifically, the strong signal intensity from As atoms was occasionally detected at the ends of the white bar-like structures, as indicated by the yellow arrows in Fig. 2(b), (c), and (e). The HAADF-STEM image in Fig. 2(f) and the EDS area mappings in Fig. 2(g) and (h) clearly show the localized As cluster at the end of the InAs core. The discontinuity of contrast was detected at the end of the InAs core, indicated by the yellow arrow in the HAADF-STEM image (Fig. 2(f)). The EDS mappings for In and As atoms clearly showed that As atoms were localized at the end of the InAs core, and the signal intensity was uniform in the area at the end of the InAs core. It is important to note that some of the vacant regions, e.g. the vacant area in region L3 in Fig. 2(a), are devoid of As clustering.
After the in situ heating process in the transmission electron microscope, HR TEM images were obtained for investigating the atomic structures. Fig. 3(a) shows the TEM image obtained from region L2 in Fig. 2(f), which was formed under the multi-beam condition. The interface between pre-existing InAs and the newly formed As cluster was abrupt, although a few thin steps with a thickness in the range of approximately two or three monolayers were observed at the interface. Even after the in situ heating TEM experiment, the thickness of the Al2O3 thin shell remained constant at 15 nm, and it was the same as that of the as-deposited film.
However, a slight change in contrast, as indicated by the blue arrow in the BFTEM image, occurred in the Al2O3 shell (Fig. 3(a)). Fig. 3(a) shows that the variation in the contrast in the Al2O3 shell is smooth in the left side, whereas it is slightly abrupt in the right side. The variation in the contrast was detected in the As cluster. Further, the variation in contrast in a TEM image can be caused by the differences in mass, thickness, crystallinity, and orientation of crystal structures in a specimen. In the case of the As cluster, the contrast change was predominantly caused by the crystallinity. The HRTEM image (Fig. 3(b)) obtained from region R1 in Fig. 3(a), in which bright-and-dark contrast is observed, shows the crystallized atomic structure in the As cluster.
On the other hand, the HRTEM image (Fig. 3(c)) obtained from region R2 in Fig. 3(a) indicates that the As cluster area in region R2 is devoid of crystallization. The inter-planar spacing (∼3.2 Å) observed in the As cluster in Fig. 3(b) is close to that of the {101} planes of the As crystal with a rhombohedral structure. Moreover, a Moiré pattern was observed near the middle of the cluster. From these observations, we deduced that As atoms might have been partially crystallized during the rapid cooling process after the in situ heating TEM experiment.
The existence of pores and/or pinholes in Al2O3 thin films deposited by the ALD method has been reported by a few research groups.22,23 The changes of atomic structures in Al2O3 deposited by ALD method might have been caused by the in situ heating experiments.24,25 Moreover, the accelerated crystallization of the amorphous Al2O3 shell by an electron beam at a low temperature (600 °C) has been reported by Huang et al.15 Therefore, we deduce that there are a few routes for the escape of atoms and/or molecules in the Al2O3 shells, probably formed during the deposition and the post in situ heating. In the actual experiment, the shape and/or morphology of the as-prepared Al2O3 thin shell gradually changed as the temperature increased; a shell shrink, a microcrack, and a small hole were finally generated at the local areas on the Al2O3 layer, indicated by the squares in Fig. 1(d). We therefore conclude that the decomposition proceeds with an increase in temperature, and the decomposed In and As atoms escape through the small paths in the thin shells.
In general, As vaporization occurs more rapidly than metal vaporization. Therefore, it causes the problem of leaving metal droplets on the surface while growing a thin film by using MBE at a high temperature.26–28
However, in our experiments, it was confirmed that a fertile region of As, not In, occurs in the Al2O3 thin shells. Our AIMD simulations show the distinct behavior of As atoms. At a highly elevated temperature of 1500 K (= 1226.85 °C), both In and As atoms rapidly fluctuate near the surface. Fig. 4(a) shows the radial distribution function (RDF) of In and As atoms near the surface plotted as a function of distance. For As atoms, the first peak followed by a local minimum is clearly shown, whereas such a peak is not observed for In atoms. The coordination number can be calculated from the integral of RDF up to the first minima, and the value of surface As atoms reads as ∼1.27.
Fig. 4 Radial distribution functions g(r) of As–As and In–In at the surface. The function of As atoms clearly shows a peak and the first local minimum indicating As–As pair, whereas In atoms do not exhibit such a behavior. (b) Evidence of existence of As2 (dimer), As3 (trimer) and As4 (tetramer) near the surface during the dynamics simulation. The red circle shows the evolution of such molecules. The atomic configurations are visualized using the VMD package.21 |
This result indicates that As atoms are more prone to form a molecule than In atoms. Indeed, we observed As dimers, trimers, and tetramers during the tens of ps molecular dynamics simulations, as shown in Fig. 4(b). According to our calculations, the size of an As trimer and tetramer is estimated as ∼6 Å, considering As–As bond lengths and the van der Waals radius of an As atom. We also calculated the diffusion coefficients from a linear regression of mean square displacements of surface atoms, and the value of In atoms is calculated as 4.38 × 10−9 m2 s−1, which is almost twice of that of As atoms, namely 2.52 × 10−9 m2 s−1. The relatively high diffusion coefficient value of In atoms strongly supports their relatively fast movements. In the InAs/Al2O3 core–shell structure, it is supposed that the diameter of the paths is generally smaller than the mean free path of gas molecules, resulting in the spreading of Knudsen diffusion. The Knudsen diffusion coefficient, DKA, in a cylindrical tube can be calculated by using the following equation:29,30
(1) |
At a fixed temperature (i.e., at 700 °C in our in situ heating experiment), the Knudsen diffusion coefficients of gas molecules were determined by the molecular weight because the radius, r, was fixed. The Knudsen diffusion coefficient of As atoms is larger than that of In atoms, but the Knudsen diffusion coefficients of As2 (dimer) and As4 (tetramer) molecules are smaller than the Knudsen diffusion coefficient of In atoms due to their heavy molecular weight. Therefore, we deduce that the In atoms escape more actively compared to As atoms through the InAs/Al2O3 core/shell system.
Footnote |
† These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2021 |