Richard S.
Schäufele
ab,
Miguel
Vazquez-Pufleau
b,
Afshin
Pendashteh
b and
Juan J.
Vilatela
*b
aDepartment of Applied Physics, Universidad Autónoma de Madrid, Cantoblanco 28049, Madrid, Spain
bIMDEA Materials, Madrid 28049, Spain. E-mail: juanjose.vilatela@imdea.org; Tel: +34 915493422
First published on 10th December 2021
Synthesis of inorganic nanowires/nanotubes suspended in the gas through floating catalyst chemical vapour deposition (FCCVD) produces exceptional growth rates of 5–1000 micron per second, several orders of magnitude faster than conventional substrate processes. It leads to nanowire lengths >100 microns and thus to the possibility of direct assembly into freestanding macroscopic networks as a continuous process. This work studies the different reaction paths controlling conversion and selectivity in FCCVD applied to the synthesis of silicon nanowires (SiNWs) from silane, grown through an aerosol of gold catalyst nanoparticles. There are two main competing reactions: catalysed growth of SiNWs and non-catalysed formation of amorphous Si nanoparticles. The mass fraction of the two populations can be precisely determined by XRD and Raman spectroscopy, enabling high-throughput screening of reaction parameter space. The experimental data and accompanying analytical model show that selectivity is kinetically controlled by the ratio of precursor/hydrogen carrier gas, through its inhibition of the pyrolisis of silane into silylene. In contrast, the rate of SiNW growth is largely unaffected by hydrogen and not limited by precursor availability. These results provide a framework to describe the kinetics of nanomaterials growth by FCCVD.
The most common methods for nanowire (NW) synthesis include wet-chemical methods, such as supercritical fluid–liquid–solid (SFLS) growth,6 top-down methods, such as electrochemical etching,7,8 and gas-phase processess based on the vapor–liquid–solid (VLS) growth mode.9 These methods have different merits according to their target application, and different challenges for scale-up. Top-down methods, for example, are constrained by the size and shape of the crystalline monolith from which the excess material is removed. Wet-chemical methods are batch processes, with throughput of NWs per batch in the range of around 100 mg, and require numerous additional processing steps to harvest the 1D nanostructures.10 The resulting NWs produced by wet-chemical methods often have moderate purity,11 with a relatively high surface contamination reminiscent of the carbonaceous precursor species and solvents used in the reaction.12
Gas-phase processes for synthesis of nanomaterials are highly scalable,13 enable large control over their chemistry and morphology, and are applicable to an enormous library of different nanomaterials.9 For nanowires and nanotubes (NT), their synthesis is commonly conducted through thermo-catalytic decomposition of precursors and incorporation into a catalyst resting on a substrate, generally following the well-known VLS mechanism.14
In a particularly interesting synthesis mode, the catalyst for NW growth is not supported on a substrate but floating in the gas stream, as demonstrated for carbon nanotubes (CNTs),15 boron nitride nanotubes (BNNT),16 III–V semiconductor nanowires,17 and more recently silicon nanowires (SiNWs).18 Growth through floating catalyst chemical vapour deposition (FCCVD) eliminates the need for substrates, enabling a continuous process where nanowires are directly assembled into macroscopic solids from the gas phase.15 The solids produced by this method have shown bulk mechanical, electrical and thermal properties in the high-performance range and above reference engineering materials.19
Very importantly, growth of nanowires through FCCVD occurs under conditions that are significantly different from substrate-based CVD. For instance, the FCCVD process is in a flow-through reactor (Fig. 1a), instead of a batch reactor, thus, mixing of the aerosol of catalyst nanoparticles and the precursor gas is essentially instantaneous and homogeneous, avoiding precursor depletion zones. As a consequence, at similar process conditions, floating catalyst particles are exposed to higher effective precursor concentrations resulting in increased catalytic activity and lower diffusion bottlenecks. These are further enhanced by unlocking higher growth temperatures, which are not inhibited anymore by substantial precursor consumption from undesired thin film growth on the substrate.17 Overall, atmospheric FCCVD leads to 20–1000 times faster growth rates than substrate processess,18,20 and produces high aspect ratios for nanowires (>200, e.g.Fig. 1b and c) and for nanotubes. This enables short residence times only in the order of <15 seconds, in contrast with common reaction times of >20 minutes for substrate-based processes.
Of particular interest is to determine the factors that control selectivity and conversion in the growth of NW by FCCVD. For CNTs, BNNTs, and III–V semiconductor NWs, high-purity materials are typically produced under very dilute synthesis conditions and low conversion, whereas reaction conditions for high conversion and/or high throughput often produce formation of a large fraction of soot-like nanoparticles and other impurities. Under conditions of high selectivity (>95% nanowire volume fraction) in GaAsNW synthesis, for example, conversion is estimated at around 0.25%.21–23 Conversion in FCCVD synthesis of carbon nanotubes ranges between 0.07–10%,24 but at high conversions typically with more than 10% volume fraction of carbonaceous impurities.
This paper sets out to study growth of 1D nanostructures by FCCVD using SiNWs as a model system. Through direct measurements of reaction throughput and sample purity-crystallinity, the results lead to maps of conversion and selectivity. These clarify the role of hydrogen, which besides acting as a “carrier” gas suppresses non-catalysed growth of solid silicon from late derivates in silane pryolisis. The experimental results are rationalised with a kinetic model describing the two dominant reactions.
This work focuses on the role of the carrier gas in the conversion and selectivity in the decomposition of SiH4 in the presence of Au. We experimented with reaction conditions that featured SiH4 molar fractions of 0–3%, H2 molar fractions of 28–72% and the remaining part being N2. Total gas flow rate ranged from 4slm to 8slm and were controlled with OMEGA FMA5500A-type mass flow controllers and Bronkhorst Coriolis flow controllers. Other synthesis parameters were fixed: reactor pressure at 900 mbar, reaction zone temperature of 650 °C, and Au catalyst nanoparticles average diameter of 22.3 nm ± 8.5 nm. Residence times are calculated from the total gas flow traversing the hot zone of the reactor tube, taken as 30 cm. They were in the range of 5–13 seconds for the experiments in this study.
SiNW samples were studied with a scanning electron microscope (SEM; FEI Helios NanoLab 600i) in order to determine average diameter of SiNWs and AuNPs, average length and average aspect ratio. These were obtained from image analysis of SEM micrographs of over 168 different SiNWs (see ESI†), utilizing the software ImageJ. High resolution transmission electron microscopy (HRTEM) images, EDX profiles and SAED patterns were taken in a FEI Talos F200X operating at 80 kV. Raman spectroscopy was carried out in a Renishaw inVia micro-Raman spectrometer with a laser wavelength of 532 nm (2.33 eV) using 50× objectives and a low power configuration in order to avoid heating effects and sample damage. Samples were weighed using a Precisa ES 125SM analytical balance. These samples were freestanding sheets directly removed from the paper filter. X-ray diffraction data were collected at the NCD-SWEET, the wide-angle/small-angle X-ray scattering (WAXS/SAXS) beamline of the Spanish synchrotron ALBA. Scattering of the samples was collected at a radiation wavelength of λ = 1.0 Å. Prior to scattering collection, the sample holder position was calibrated using silver behenate (AgBh) for SAXS and chromium oxide (Cr2O3) for WAXS. The patterns were first corrected for the background scattering and then analyzed using DAWN software (v. 2.20), obtaining radial profiles after azimuthal integration.
Mathematical fitting of experimental data using the proposed kinetic models was done by iteration, using a supervised minimum square methodology. Details are included in ESI.†
The SiNW morphology and overall purity of materials synthesised in this FCCVD mode depends on various synthesis parameters, but in our exploration of parameter space so far, we have observed a strong dependence on the molar fraction of H2, N2 and SiH4. As an example, in Fig. 2e–h we show compositional information of a sample produced with a higher SiH4 molar ratio and a lower H2 molar ratio (SiH4:
H2
:
N2 molar ratio of approximately 2.7
:
54
:
43.3), but otherwise under nominally identical synthesis conditions in terms of catalyst concentration and size, reaction temperature and pressure. Residence time is also similar, at 6.9 s and 8.4 s for the samples in Fig. 2a and e, respectively. The second sample has a dark brown colour and contains a high degree of impurities in addition to SiNWs. These impurities are irregular-shaped particulates, either aggregated as clusters trapped between the SiNWs network or coating the SiNW surface (Fig. 2e). Indeed, we observe a substantial increase in apparent nanowire diameter (44.0 nm ± 15.1 nm) in this sample, while average size of Au seed nanoparticles remains at approximately 24.4 nm ± 6.8 nm. Via EDX analysis and selected area electron diffraction (SAED) we find the impurities to be amorphous silicon (a-Si) (Fig. 2f–g).
![]() | ||
Fig. 3 Crystallinity is plotted against H2 molar fraction and SiH4 molar fraction for various FCCVD-grown samples. Colour coding of data points is maintained throughout Fig. 3, 4 and Fig. S4† in order to simplify tracking of individual experiments on the various presented plot formats. Insets highlight respective sample colours. |
The compositional difference between SiNW samples are clearly visible in their Raman spectra. The spectrum for the high purity sample (Fig. 2d) has a sharp peak corresponding to the 1st order TO phonon of crystalline Si (c-Si), a small peak at 483 cm−1 found in SiNWs27 (Fig. S1†), and a very weak peak at 496 cm−1 corresponding to the 1st order TO phonon of a-Si (Fig. S1†). The same three peaks are also present in the low purity sample (Fig. 2h), however, the a-Si peak is significantly more pronounced. The higher abundance of a-Si is in line with SEM and TEM analysis.
We found a direct correlation between the relative intensity of the Raman peaks for a-Si and c-Si, and the overall sample purity determined by electron microscopy (Fig. S2†). Furthermore, because the impurities produced in this temperature range are amorphous, a Raman-derived nanowire purity can be obtained which is directly correlated with crystalline mass fraction, (Fig. S6†).28 Hence, sample crystallinity (purity) can be quantitatively determined from the Raman spectra, calculated for example as:
![]() | (1) |
Very importantly, this crystallinity parameter is indicative of the selectivity of the FCCVD reaction. Since the weighted Raman intensities are proportional to the mass fraction of the different phases (Fig. S6†), C is equivalent to the fraction of material synthesised through growth catalysed by Au.‡
The approach is demonstrated in Fig. 3. It shows the crystallinity and colour of SiNW samples produced at different SiH4:
H2
:
N2 molar fractions plotted against C. Crystallinity increases with a higher content of H2 and a lower content of SiH4, reaching a plateau at approximately 98% C. This is consistent with SEM analysis of the samples, presented in Fig. S2,† although above
= 93% a-Si cannot be detected by SEM.
A large space of synthesis conditions was found to lead to successful growth of samples with predominantly SiNWs. The most relevant to understand the relation between selectivity and conversion are summarised in Fig. 4. The colour codes represent different experimental sets, for example, at constant total gas flow rate (yellow spheres), constant flow rate of hydrogen and nitrogen gas (green spheres), or constant molar ratio of SiH4/H2 (red spheres).
Fig. 4a presents a plot of crystallinity for different molar fractions of SiH4 and H2 (the remaining fraction is N2). It shows that conditions leading to high crystallinity correspond to a low SiH4 concentration and/or high H2 content. For instance, high crystallinity of above 95% is only achieved for SiH4 molar fractions below 2%, and requiring between 34% and 64% H2. In other words, selectivity is favoured under dilute conditions in hydrogen.
Fig. 4b presents the conversion data for the same reaction conditions. The data follow the opposite behaviour. Conversion increases with increasing SiH4 molar fraction, which is expected, but H2 has also a strong effect. The yellow datapoints, for example, show that replacement of H2 with N2 at constant SiH4 molar fraction increases conversion from 0.13% to 0.56%.
The balance between crystallinity and conversion is more clearly observed in Fig. 4c. It shows a clear drop in crystallinity with increasing conversion, with the data falling on an exponential decay curve. The curve represents the limits of conversion and selectivity for different combinations of SiH4:
H2
:
N2 fractions. Its absolute values depend on specific details of the synthesis reactor, but the trend is inherent to the FCCVD reaction.
Given the observed dependence of the reaction on the fractions of both SiH4 and H2, it is insightful to analyse a plot of conversion against the molar ratio of SiH4/H2, as shown in Fig. 4d. In this representation the values of crystallinity define areas in the two-dimensional graph. For example, the region of crystallinity above 88% stretches up to a conversion of approximately 0.25% and corresponds to a low SiH4/H2 ratio. Conversely, a high SiH4/H2 ratio increases conversion almost linearly to 0.6%, but reduces crystallinity below 70%.
Solid silicon is a stable product of the FCCVD reaction; it has a lower chemical potential μsolid than silicon as a gas and than the eutectic with Au.30 This thermodynamic product can be formed via two different mechanisms: VLS or vapor–solid (VS) growth. The VLS route is catalysed by gold and leads to growth of crystalline NWs, with the lower activation energy (reported to be between 17.5–22 kJ mol−1).31,32
The main competing mechanism is non-catalyzed formation of solid Si through direct deposition as particulates or a coating on nascent nanowires. The corresponding activation energy is reported to be 29–38 kJ mol−1 and thus higher than what is required in the first step of VLS growth. The selectivity of the process could potentially be affected by hydrogenation/passivation of active surfaces at SiNWs and a-Si particles in the presence of H2 gas, which could change the activation energy for the catalysed and non-catalysed reactions33,34 and thus favour one of the products. This is unlikely though. The small SiH4/H2 molar ratios used in our study imply that both the catalyst and nascent NWs are essentially in a H2-rich atmosphere for all reaction conditions.
Instead, the selectivity of the process is controlled by the effect of hydrogen on the kinetics of the two reactions. And its effect is particularly prominent because of the high collision rate between SiHx molecules compared to the rate of arrival on Au catalyst nanoparticles.
For both the VLS and VS routes, the conversion of SiH4 into solid Si occurs through silane derivates. As a product of the reaction, hydrogen is known to affect the reversible decomposition reactions of SiH4 into silane derivates:35
![]() | (2) |
We expect that early derivates can take part in the VLS route, for example via further decomposition at the Au catalyst particle. In contrast, non-catalysed formation of solid Si from silane pyrolisis is likely to require late derivates. They form via silane pyrolysis by undergoing a series of polymerization reactions29,36,37 which release hydrogen. Therefore, based on Le Chatelier's principle, higher hydrogen availabilty is expected to reduce reaction rate. Indeed, the kinetic constants reported for non-catalised formation of nanoparticles are a factor of 2 lower for early reaction times (conversions between 3–35%) than for later reaction conditions (conversion between 39–75%).35 Thus, the role of hydrogen in the selectivity of the process is in avoiding the VS route through suppression of extended SiH4 decomposition. This is schematically shown in reaction (3), where multiple lines indicate enhanced reaction pathways:
![]() | (3) |
If selectivity is governed by the suppression of a-Si, this implies that the reaction conditions leading to an increase in conversion simply produce more a-Si in addition to SiNWs. This is indeed the case. From electron microscopy analysis of samples covering a wide range of crystallinity and reaction conversion values we find the number density of SiNWs to remain fairly constant (Fig. S4d†). This is a consequence of the catalyst activity being close to 100% and largely independent of SiH4/H2 ratio (see ESI†).
In order to simplify the problem and constrain its mathematical complexity the following assumptions are made:
(a) The experimental reactor is well described by the plug flow reactor model, i.e. there is perfect mixing in the radial direction but no mixing along the axial direction.
(b) The reactor is considered isothermal, equal to the setpoint of the furnace.
(c) The reaction goes to completion until exhaustion of silane. This is supported by kinetic estimation showing that even in the case of inhibited reaction, more than 99% percent of the reaction is over after the first 0.5 s. In contrast, the lower bound of residence time in the reactor is about 5 s.
(d) The losses to the reactor walls or formation of side-products, such as sub-2 nm clusters,37 do not affect the ratio of the two primary measurable products i.e. nanowires and amorphous nanoparticles which are described by eqn (1) and expressed as a function of mass in eqn (5)
Under these assumptions, the rate of disappearance of silane is equal to the sum of the rate of formation of nanoparticles and nanowires:
![]() | (4) |
Then crystallinity C can be obtained by solving the following equation:
![]() | (5) |
The upper bound t of the integral is until completion based on the assumption in (c).
A very approximate expression to describe the growth rate of silicon nanowires is reported in the literature for synthesis catalysed by Au nanoparticles supported on a substrate.38 It relates longitudinal growth velocity (vnw, in nm) to silane molar fraction (xSiH4,percentage), total pressure (Ptot, in mbar) and temperature (T, in K) by:
![]() | (6) |
It follows that the mass production rate of nanowires (ṁnw) is given by multiplying its linear growth velocity vnw times its cross sectional area, times the total number of nanowires Nonw and density of silicon ρ:
![]() | (7) |
However, eqn (6) and (7) predict indefinite growth with increasing precursor availability, which is adequate to describe growth for substrate CVD under low pressure (<1.6 mbar, equivalent to around 0.18 silane molar ratio). In the synthesis of 1D nanomaterials by FCCVD growth rate is exceptionally fast, between approx. 1 μs−1 for nanowires and 100000 μs−1 for nanotubes, and may thus be limited by other factors. Indeed, we find that SiNW aspect ratio is very high for all synthesis conditions and has a very weak dependence on hydrogen concentration (Fig. S4†). Therefore, eqn (7) is valid up to a silane saturation concentration [SiH4]sat.
The nanoparticle formation rate is taken as the kinetics of homogeneous nucleation from Hogness et al.,39 who derived this kinetic expression and constant based on pressure changes in a reactor fed with high purity silane.
![]() | (8) |
![]() | (9) |
Eqn (4) and (9) have been succesfully used in studies of competing reactions during nanoparticle growth,34 but not under variable hydrogen concentration. This can be included in the model by considering a simple silicon nanoparticles formation mechanism following the reaction:
![]() | (10) |
And thus
![]() | (11) |
However, this expression is a function of the concentration of the silylene radical , which is not known and very difficult to measure due to its high reactivity and transient nature. Therefore, we reformulate the rate expression of silicon nanoparticle synthesis as a function of more easily measurable variables. The first step is to formulate the rate of formation of the silylene intermediate
, which is given by:
![]() | (12) |
Using the pseudo steady state hypothesis, which assumes that the intermediate radical is consumed as quickly as it is generated and therefore its concentration remains constant during the reaction, we obtain:
![]() | (13) |
![]() | (14) |
![]() | (15) |
When no hydrogen is present this expression reduces the case where the rate limiting step is the formation of silylene.29,40
Through an iterative process, we use the growth rate eqn (15), (7) and (4) to fit the experimental data for selectivity as a function of precursor and hydrogen concentration, with three unknown parameters: the factor related to supression of nanoparticle formation (k3/k2), the growth enhancement parameter for FCCVD growth (ηFCCVD), and the silane saturation concentration above which it does not limit growth rate([SiH4sat]). As shown in Fig. 5, there is good agreement between the experimental and calculated data (see supplementary video of 3D plot). The mathematical description of inhibition of nanoparticles formation at higher hydrogen concentration is well described by the homogeneous nucleation inhibition factor. Similarly, the decline of nanowire presence at higher silane concentration is in agreement with the proposed saturation effect.
From the fitting, we obtain that the a-Si nanoparticle suppression term is k3/k2 = 1.4 × 107 cm3 g−1. The only report found on the kinetics of silane decomposition for nanoparticle growth by White et al.35 has a large scatter in kinetic constants, preventing a comparison with the value obtained in this work but in qualitative agreement with our results. The nanowire growth enhancement factor ηFCCVD = 25, which is a factor of 25 times higher than proposed for substrate-grown SiNWs.38 The saturation precursor concentration [SiH4sat] gives 0.05% molar ratio, which is equivalent to a collision rate with catalyst particles of 109 Si-atoms per s. This is 1/100th relative to the diffusion of Si in a Si–Au eutectic liquid particle, thus, it is unlikely that Si–Au diffusion is the rate limiting step in NW growth, which is in agreement with observations on substrate-grown SiNWs31 Similarly, we note that under conditions of nearly 100% growth of SiNWs, the ratio SiH4/hydrogen fraction does not produce significant changes in aspect ratio. This implies that the collision rate of silane, or its derivates, on the catalyst, is not the main factor limiting SiNW length. In view of these results, we expect that the limiting processes are either precursor decomposition (including potential adsorption and reaction of precursor at the AuNP surface, and desorption of byproducts) or the precipitation of the SiNW.
Lastly, this work helps rationalise growth of highly pure SiNWs by FCCVD and the competing processes to increase throughput and sample purity. This aspect is important for envisaged applications of SiNW networks, particularly as high-performance lithium-ion batteries anode materials41–43 and sensors.44
Footnotes |
† Electronic supplementary information (ESI) available: Propagation of relevant Raman peaks at varying gas mixture ratios, correlation of SEM analysis and crystallinity, c-Si contamination on SiNWs at high temperatures, SiNW info on diameter, length, aspect ratio and number concentration, morphology of sample grown without hydrogen, methodology of determining crystallinity, explanation of approach of modelization of experimental data. See DOI: 10.1039/d1nr07261g |
‡ This assignation is valid for all synthesis experiments shown in the manuscript, performed at 650 °C and 900 mbar. Impurities in these samples – if present – were always a-Si of similar morphology to that shown in Fig. 2e–h. At around 720 °C and 900 mbar, the reaction also produced non-catalysed crystalline particulates (see ESI Fig. S3†), in agreement with the threshold found in particle formation studies on silane pyrolisis.29 |
This journal is © The Royal Society of Chemistry 2022 |