Next Article in Journal
Thickness-Dependent Photoelectrochemical Water Splitting Properties of Self-Assembled Nanostructured LaFeO3 Perovskite Thin Films
Next Article in Special Issue
Electromechanical Characteristics by a Vertical Flip of C70 Fullerene Prolate Spheroid in a Single-Electron Transistor: Hybrid Density Functional Methods
Previous Article in Journal
Computational Study on Surface Bonding Based on Nanocone Arrays
Previous Article in Special Issue
Artificial Neural Network-Based Prediction of the Optical Properties of Spherical Core–Shell Plasmonic Metastructures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Study of Quenching Effects on Growth Processes and Size Distributions of Silicon Nanoparticles at a Thermal Plasma Tail

1
Graduate School of Engineering, Tohoku University, Sendai 980-8579, Japan
2
National Institute of Advanced Industrial Science and Technology, Nagoya 463-8560, Japan
3
Alma Mater Studiorum, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(6), 1370; https://doi.org/10.3390/nano11061370
Submission received: 28 April 2021 / Revised: 18 May 2021 / Accepted: 18 May 2021 / Published: 21 May 2021
(This article belongs to the Special Issue Computational Study of Nanomaterials)

Abstract

:
In this paper, quenching effects on silicon nanoparticle growth processes and size distributions at a typical range of cooling rates in a thermal plasma tail are investigated computationally. We used a nodal-type model that expresses a size distribution evolving temporally with simultaneous homogeneous nucleation, heterogeneous condensation, interparticle coagulation, and melting point depression. The numerically obtained size distributions exhibit similar size ranges and tendencies to those of experiment results obtained with and without quenching. In a highly supersaturated state, 40–50% of the vapor atoms are converted rapidly to nanoparticles. After most vapor atoms are consumed, the nanoparticles grow by coagulation, which occurs much more slowly than condensation. At higher cooling rates, one obtains greater total number density, smaller size, and smaller standard deviation. Quenching in thermal plasma fabrication is effectual, but it presents limitations for controlling nanoparticle characteristics.

1. Introduction

Nanoparticles offer unique electronic, optical, magnetic, and catalytic properties that differ from those of larger micrometer-size particles and those of bulk materials [1]. Nanoparticles are therefore highly attractive for use in applications expected to lead to breakthroughs in biomedical, environmental, and industrial fields. Even pure silicon nanoparticles are anticipated as promising materials for use in numerous applications such as lithography [2], electronic devices [3], electron transistors [4,5], floating gate memory devices [6,7,8,9], luminescent thermometers [10], and lithium ion battery electrodes [11,12,13]. Several Japanese groups have recently used thermal plasmas to achieve high-throughput production of silicon nanoparticles for lithium ion battery electrodes [14,15,16,17,18].
Thermal plasma, a kind of high-enthalpy plasma, is generated by high-power discharge at around atmospheric pressure. Thermal plasma is a unique fluid with high temperatures not only of electrons but also of heavy species, with high chemical reactivity, intensive light emission, and varying transport properties [19]. It is also an electrically conductive fluid. Therefore, it can be controlled by an electromagnetic field [20,21]. The high temperature of its heavy species, around 10,000 K, can vaporize even materials that have a high melting point or boiling point. In addition, the steep temperature-decrease gradient at its tail plays an important role in rapid conversion from material vapor to nanoparticles. By virtue of these two features, thermal plasma achieves one-step fabrication of nanoparticles with high yield [22]. Indeed, based on the use of inductively coupled thermal plasma (ICTP), Tanaka et al. [23] and Kodama et al. [24] synthesized titanium-dioxide nanoparticles at 6.7–12.3 g min−1 using a pulse modulation technique. Ohta et al. [15] produced silicon nanoparticles at 17 g min−1 using direct-current plasma jet assistance.
The number and size of nanoparticles are affected by temperature and flow fields in and around thermal plasma through growth processes [25]. Strong cross-correlations between the temperature at an upstream plasma fringe and nanoparticle concentration in a downstream region have been found for a transferred arc plasma system [26,27] and for a non-transferred arc plasma jet in two-dimensional [28] and three-dimensional fields [29] in recent numerical studies using a method suitable for simulating dynamic behaviors of thermal plasma–nonionized gas coexisting flows [30]. Actually, to control temperature and flow fields and to control nanoparticle growth processes in ICTP systems, several experimental studies adopted quenching strategies using a water-cooled coil [31], a water-cooled ball [15,32], direct plasma control by pulse modulation [33], counterflow injection [18,34], and radial gas injection [35,36]. Effects of those quenching methods were also investigated using two-dimensional simulations for practical conditions [37,38,39,40,41,42,43,44,45,46,47]. From a more fundamental perspective, specific effects of quenching on nanoparticle size distributions were investigated numerically using time-dependent models [48,49,50,51]. Although those works explicitly presented results obtained for spatial distributions or size distributions of nanoparticles, none clearly presented an implicit mechanism that dominates nanoparticle growth processes.
This study investigated quenching effects on the growth processes and size distributions of nanoparticles at a typical range of cooling rates at a thermal plasma tail. The investigative method is numerical, using a nodal-type model that expresses a size distribution evolving temporally with simultaneous homogeneous nucleation, heterogeneous condensation, interparticle coagulation, and melting point depression. To validate the present numerical calculations, silicon nanoparticle fabrication is demonstrated experimentally with and without quenching. Results of the study reveal the implicit mechanisms of processes related to conversion from vapor atoms to nanoparticles. Moreover, results clearly illustrate cooling rate effects on the total number, size, and dispersion of the produced nanoparticles.

2. Experiment Setup

Figure 1 portrays schematic illustrations of the nanoparticle fabrication systems (TP-40020NPS; JEOL Ltd., Tokyo, Japan) consisting of an ICTP torch and a chamber without or with a quenching portion. After argon gas (G1 grade, less than 0.1 ppm oxygen) was filled in the main chamber at 100 kPa, argon gas was injected continuously at 35 L min−1 from the torch top to sustain the plasma. Using a powder feeding system (TP-99010FDR; JEOL Ltd., Tokyo, Japan), coarse silicon powder (approx. 7 μm particle size, 99.99% purity; Kojundo Chemical Lab. Co., Ltd., Saitama, Japan) was then introduced through the feeder nozzle into the plasma torch at a 0.048 g min−1 feed rate, with carrier argon gas added at a flow rate of 3 L min−1. When the flow was quenched, additional argon gas was injected from the lower part of the torch toward the central axis at 80 L min−1. Details of the system geometries and discharge parameters were identical to those described for an earlier study [52,53].
Observation was achieved by using scanning electron microscopy (SEM: JSM-7800F; JEOL Co., Ltd., Tokyo, Japan). Size distributions for approximately 1200 nanoparticles were measured from SEM micrographs.

3. Numerical Model Description

3.1. Model Outline and Assumptions

Following evaporation in the plasma, the material vapor is transported through the plasma to the chamber. Vapor reaching the plasma tail experiences a rapid temperature decrease. There, the vapor becomes supersaturated. Numerous nuclei are generated (homogeneous nucleation). Once nucleation occurs, the vapor condenses on the nuclei (heterogeneous condensation). Through these processes, nanoparticles grow. Simultaneously, the nanoparticles mutually collide and merge to form larger nanoparticles (interparticle coagulation) if they are still liquid and exist at temperatures higher than their melting points.
Those nanoparticle diameters vary widely from sub-nanometer scale to a few hundred nanometers, which correspond to monomer numbers in a single particle from 101 to 108. Furthermore, the spatial number density of the nanoparticles reaches 1018 m−3. Therefore, molecular dynamics calculations for every nanoparticle’s growth are impossible in a practical sense. Aerosol dynamics modelling makes it possible to calculate the collective and simultaneous growth processes of the numerous nanoparticles. Based on the aerosol representation, nanoparticle growth processes at thermal plasma tails have been predicted numerically using several computational models with different features summarized in Table 1. The models are classified in four types. Type A is a sophisticated model with a simple mathematical description treating nanoparticles’ growth by simultaneous homogeneous nucleation and heterogeneous condensation; however, it does not take interparticle coagulation into account [51,54,55,56]. It can express a size distribution with any shape so that it needs high costs of computational resources and time. Type B describes collective nanoparticles’ growth by simultaneous nucleation, condensation, and coagulation by a simple set of two aerosol equations and one vapor equation so that the computational costs are low [26,27,28,29,57,58]. However, it does not express a dispersed size distribution because it calculates physical quantities by averaged size approximation. Type C has been an often-used model which deals with collective nanoparticles’ growth by simultaneous nucleation, condensation and coagulation as well [39,40,42,43,44,45,46,47,59,60,61,62]. Although its computational costs are low, its mathematical formulation is more complex; nevertheless, it still requires an assumption of a lognormal size distribution. Type D can express a size distribution with any shape during collective nanoparticles’ growth through simultaneous processes of nucleation, condensation, and coagulation [37,48,49,53,63,64,65,66,67,68,69,70,71,72,73]. Its mathematical description is more complex; furthermore, its computational costs are higher because the equations as many as the nodes discretizing the size distribution.
This study adopts Type D with several assumptions: (i) spherical nanoparticles; (ii) negligible nanoparticle inertia; (iii) identical transport velocity and temperature of nanoparticles to those of the bulk gas flow; (iv) negligible heat generation caused by condensation; (v) negligible electric charge of nanoparticles; (vi) material vapor regarded as an ideal gas; and (vii) atmospheric pressure in the reaction chamber. Actually, Model type D discretizes the size distribution in the space of particle sizes expressed by particle volume. The discretized volumes are assigned as
v k + 1 = f q v k         ( k = 1 , 2 , . . . , k m a x )
where v represents the particle volume, k stands for the discrete point number for the particle size, and f q is a coefficient to determine the interval. f q = 1.16 and k m a x = 161 are set for this study. The minimum volume of a “particle” is chosen as v 2 = 10 v 1 , where v 1 denotes the volume of a monomer, which is equivalent to a vapor atom. In this size space, the net production rate of particles with volume v k is written as
Δ N 1 Δ t = N ˙ k ( n u c l ) + N ˙ k ( c o n d ) + N ˙ k ( c o a g )         ( k = 2 , . . . , k m a x )
where N k represents the number density of particles at node k, N ˙ signifies the increment per unit time, and superscripts (nucl), (cond), and (coag) stand for the contributions of homogeneous nucleation, heterogeneous condensation, and interparticle coagulation,   respectively .

3.2. Homogeneous Nucleation

Figure 2a portrays a schematic image showing calculations for homogeneous nucleation. When material vapor is transported to the plasma tail, its temperature decreases along with decreased saturation pressure of the material. Consequently, the material vapor becomes supersaturated, which engenders homogeneous nucleation. Therein, the production rate of particles having volume v k by homogeneous nucleation is written as
N ˙ k ( n u c l ) = I ξ k ( n u c l )         ( k = 2 , . . . , k m a x )
here, I is the homogeneous nucleation rate [74] obtained as
I = N s 2 S 2 σ π m 1 exp [ Θ 4 Θ 3 27 ( ln S ) 2 ]
where N s signifies the number density of the vapor atoms in the equilibrium state, σ stands for the surface tension, and m 1 denotes the mass of a monomer. Additionally, S expresses the supersaturation ratio defined as
S = N 1 N S = p 1 p S
where N 1 and p 1 respectively represent the number density and pressure of vapor monomers, and p S denotes the saturation pressure. Furthermore, Θ stands for the dimensionless surface tension, which is given as
Θ = σ s 1 k B T
where s 1 represents the monomer surface area, k B is Boltzmann’s constant, and T stands for the temperature. In Equation (3), ξ k ( n u c l ) represents the partition coefficient for discretized size nodes written as
ξ k ( n u c l ) = { v * v k ( v k v * < v k + 1 ) v * v 1 ( v * < v 1 ) 0 ( o t h e r )
here, v * represents the volume of a nucleus that is regarded as a particle having the critical diameter given as described in an earlier report [74],
v * = π 6 ( 4 σ v 1 k B T ln S ) 3
Therefore, the critical diameter is equivalent to the diameter of a nucleus generated by homogeneous nucleation in this model.

3.3. Heterogeneous Condensation

Figure 2b represents a schematic image of calculation for heterogeneous condensation. When nuclei are generated by homogeneous nucleation, the material vapor is still supersaturated, which causes heterogeneous condensation. Therein, the net production rate of particles having volume v k by heterogeneous condensation can be written as
N ˙ k ( c o n d ) = i = 2 k m a x ( ξ i , k ( c o n d ) δ i , k ) N i Δ t         ( k = 2 , . . . , k m a x )
here, ξ i , k ( c o n d ) is the partition coefficient for heterogeneous condensation, which is given as
ξ i , k ( c o n d ) = { v k + 1 ( v i + Δ v i ) v k + 1 v k ( v k v i + Δ v i < v k + 1 ) ( v i + Δ v i ) v k 1 v k v k 1 ( v k 1 v i + Δ v i < v k ) 0 ( o t h e r )
In Equation (9), δ i , k is Kronecker’s delta, and
δ i , k = { 1 ( i = k ) 0 ( i k )
The growth rate of particles having diameter d i is given as the following equation [64] modified from a general description as [75]
Δ v i Δ t = 2 π d i D 1 v 1 ( N 1 N S i ) [ 0.75 α ( 1 + Kn i ) ( 0.75 α ) + 0.2832 α Kn i + Kn i + Kn i 2 ]
where D 1 is the diffusion coefficient of monomers [76], Kn is the Knudsen number defined as the ratio of the mean free path of gas to the particle radius, and α is the accommodation coefficient, which is set as 1.0 for this study. Additionally, N S i denotes the number density of equilibrium vapor atoms corrected considering the surface curvature for node i as [77]
N S i = N S exp ( 4 σ v 1 d i k B T )
Particles with volume v i gain volume Δ v i . Consequently, they form new particles having volume v i + Δ v i . The new particles are split into adjacent nodes under the mass conserving condition by Equation (10).

3.4. Interparticle Coagulation

Figure 2c portrays a schematic image of calculation for interparticle coagulation. The growing nanoparticles mutually collide and coagulate. When their temperature becomes higher than their melting point, the particles merge and form larger particles, thereby decreasing their own numbers. The net production rate of the particles at node k by interparticle coagulation can therefore be expressed as
N ˙ k ( c o a g ) = 1 2   i = 2 k m a x j = 2 k m a x ξ i , j , k ( c o a g ) β i , j N i N j N k i = 2 k m a x β i , k N i         ( k = 2 , . . . , k m a x )
where β i , j represents the collision frequency function between two particles having volumes v i and v j covering a wide size range from the free molecular regime to the continuum regime [64]. The first term and the second term of the right-hand side in Equation (14) express the gain of node k by the collision between the particles at the other nodes i and j, and the loss of node k by the collision between the particles at node k and the particles at the other nodes, respectively. Particles with volume v i and v j collide and coagulate. They consequently form new particles with volume v i + v j . Figure 2c shows that the new particles are split into the adjacent nodes under the mass conserving condition presented in Figure 2. The partition coefficient for coagulation of size ξ ( c o a g ) is given as
ξ i , j , k ( c o a g ) = { v k + 1 ( v i + v j ) v k + 1 v k ( v k v i + v j < v k + 1 ) ( v i + v j ) v k 1 v k v k 1 ( v k 1 v i + v j < v k ) 0 ( o t h e r )

3.5. Vapor Consumption

The total number of material atoms is conserved through conversion from vapor to particles. The population balance equation of the material vapor is calculated simultaneously considering vapor atom consumption by nucleation and condensation:
Δ N 1 Δ t = k = 2 k m a x I ξ k ( n u c l ) n * i = 2 k m a x N i Δ v i v 1 Δ t
In Equation (16), n * represents the number of monomers composing a thermodynamically neutrally stable nucleus. The number density of vapor atoms directly affects the rates of nucleation and condensation.

3.6. Melting Point Depression

Melting points of nanoscale particles are lower than those of bulk materials [78]. This melting point depression is considered in the present model. The melting point T m e l t ,   k of particles with diameter d k can be estimated as
T m e l t ,   k = T m e l t ,   b u l k × ( 1 ε d k )         ( k = 2 , . . . , k m a x )
Therein, T m e l t ,   b u l k denotes the melting point of a bulk material. Parameter ε is a characteristic property determined by the solid and liquid surface energies and the bulk melting enthalpy. This study adopts the value of ε = 1.16 nm for silicon spherical particles [78]. Figure 3 depicts the melting points of bulk silicon and nanoscale spherical silicon estimated by Equation (17). The present model treats a particle with temperature lower than this depressed melting point as a solid body and does not calculate the interparticle coagulation if both particles are solid in computation.

3.7. Computational Conditions

Computations were performed under typical conditions at a small-type ICTP’s tail where the temperature decreases monotonically at 18 cooling rates ranging from 2.0 × 103 K s−1 to 6.4 × 104 K s−1. In response to these cooling rates, the time increment Δ t for computations was set as 5.0 μs to 0.15625 μs, which gave sufficient resolution for the processes. Corresponding to the silicon feed rate of the experiment, it is assumed that, after complete evaporation by the plasma, the silicon vapor atoms are transported with the argon carrier gas. Therein, the molar fraction of silicon to argon is set to 1.257% at 2700 K as the starting condition of the computations. The material properties of silicon were referred from a textbook used for this field [79].

4. Results and Discussion

4.1. Experiment Results and Model Validation

Figure 4 depicts SEM images of silicon particles obtained with and without quenching. The primary particles are spherical particles having nanometer-range diameters. Apparently, the particles obtained with quenching are much smaller than those without quenching. Figure 5 shows the primary particle size distributions, which quantitatively indicate the same tendency in the experiments and the computations. With quenching, the size distributions exhibit smaller and narrower ranges. The experiment without quenching obtained mean diameter of 150 nm and standard deviation of 82 nm. The experiment with quenching obtained the mean diameter of 27 nm and standard deviation of 9 nm. On the other hand, computation with the cooling rate of 3.1 × 103 K s−1 obtained mean diameter of 149.5 nm and standard deviation of 39.7 nm. Computation with the cooling rate of 6.2 × 104 K s−1 obtained mean diameter of 26.8 nm and standard deviation of 9.98 nm. Those computation results agree to an acceptable degree with the experimentally obtained results. Therefore, in this report, the processes with the cooling rates of 6.2 × 104 K s−1 and 3.1 × 103 K s−1 are discussed as conditions with and without quenching, respectively.
In spite of that agreement, the results include some differences among them. The causes of the differences are discussed hereinafter. In general, experiments of practical material processing by thermal plasmas inherently include errors because of temporal fluctuations and spatial non-uniformity of temperature and flow fields, as reported from experimental studies [80,81] and earlier computational studies [26,30,82,83,84,85]. Such turbulent-like fluid-dynamic behaviors affect not only the evaporation of precursory particles [85,86] but also nanoparticle formation and transport [21,28,29]. Indeed, in the experiment reported herein, some large unvaporized precursor particles were found. On the other hand, the present numerical demonstrations were conducted under ideal conditions assuming complete evaporation of the precursor particles and nanoparticle growth at linear cooling rates. The electric charge might also affect the nanoparticle growth, even in the region downstream from the plasma tail [73]. However, that effect was neglected in the present computation. Recent studies conducted by computation [87,88] and experimentation [89,90] have demonstrated that ion transport by cataphoresis affected the material vapor distribution in and around plasma.
In addition, the surface tension of bulk silicon was used to calculate homogeneous nucleation by Equation (4) and the critical size by Equation (8). The diameters of major nuclei generated in this study were estimated as below a few nanometers, for which the surface tension of the bulk material might have not been valid. Validated data for such small sizes were unavailable. Therefore, bulk silicon data were used. It is noteworthy that numerical calculations of homogeneous nucleation rates and critical sizes inherently lead to errors because of this limitation. Nevertheless, Equations (4) and (8) are still accepted as the most reliable equations available.
It is noteworthy that nodal discretization of the size space introduces errors. To reduce error accumulation, the large number of nodes with small geometric spacing factor fq in Equation (1) should be used, as discussed in an earlier report of the relevant literature [66]. In the present study, several combinations of the node number and geometric spacing factor were tested. For instance, double the number of nodes led to changes of less than 1.5% in the particle size distribution. The set of an adequately large number of nodes and small geometric spacing factor was chosen to restrain error accumulation.
Heat generation by condensation was neglected. If all the vapor would be condensed, then the heat generation would be estimated as 1.1 W, which would heat the surrounding argon gas by approximately 24 K. Condensation occurred at 2500–2600 K. Therefore, the temperature increase effect was negligible. Even if no heat were transferred from silicon nanoparticles to argon gas, the temperature increase could be considerable and the temperature difference between them could be large. Actually, the condensation heat would be transferred to both and both temperatures would increase. If their temperature difference was large, then it would affect the growth rate. An additional model should be implemented to elucidate that point.
It might be significant to discuss about agglomerates as secondary larger particles as shown in SEM images of Figure 4. The rietveld analysis was performed from the X-ray diffraction (XRD) results. Calculated by the Halder-Wagner plot, the crystallite sizes were 73 nm with quenching and 74 nm without quenching. None of them were significantly smaller than the particle size obtained from the images. With quenching, the crystallite size is considered to be larger than the particle size due to the presence of unvaporized precursor particles. Meawhile, without quenching, the crystallite size is about half that of the particle size, but the effect of twinning during the solidification process from droplets cannot be ruled out. Therefore, in order to confirm these things experimentally, it is necessary to analyze the microstructure by Electron Back Scattered Diffraction (EBSD) measurements. A computational work treating aggregation was also reported using a different-type model [46], whereas the main scope of this study is more fundamental clarification of the vapor-to-particle conversion stage in quenching processes. The implementation of an aggregation model would be addressed in a further investigation.
The present model still has room for improvement. Nevertheless, the discussion presented above reinforces the assertion that the present computation validly predicts silicon nanoparticle growth processes in an ICTP system.

4.2. Implicit Mechanism of Collective Nanoparticle Growth

Figure 6 presents evolutions of molar densities of silicon vapor atoms and of the saturated state with and without quenching. The horizontal axes showing temperature are reversed corresponding to time progress with the temperature decrease. Both with and without quenching, the molar densities of the vapor atoms increase with the decrease in temperature. Then they exceed that of the saturated state. When the supersaturation ratios defined as Equation (5) reach the maximum values of 1.63 at 2586 K without quenching and 1.87 at 2567 K with quenching, the vapor atoms begin to decrease drastically. These decreases are caused by rapid conversion of the vapor atoms to nanoparticles.
Figure 7 portrays evolutions of the conversion ratios, which represent how many vapor atoms have changed into nanoparticles. Approximately 40% and 50% of the vapor atoms are converted to nanoparticles at those rapid conversions. When the molar densities of the vapor atoms become slightly higher than that of the saturated state, as depicted in Figure 6, the conversions proceed moderately. For both conditions, 99% of the vapor atoms complete the conversions at 2100 K.
Figure 8 portrays the evolutions of vapor consumption rates by homogeneous nucleation and heterogeneous condensation. Without and with quenching, the rates of nucleation and condensation increase remarkably below 2600 K. When the number of nuclei increases and consequently the total surface area increases, vapor consumption by condensation becomes overwhelming. The consumption rate of condensation is approximately 4.0 × 10 5 times as large as that of nucleation without quenching. On the other hand, that of condensation is approximately 340 times as large as that of nucleation with quenching. After the number densities of the vapor atoms decrease to that of the saturated state, the gradients of the consumption rates of condensation become more moderate, thereby leading to the moderately increasing curves in Figure 7.
Figure 9 depicts the evolutions of the homogeneous nucleation rates. As the vapor atoms become highly supersaturated, the nucleation rates increase and reach the highest values of 1.4 × 10 17 m−3s−1 at 2586 K without quenching and 4.8 × 10 21 m−3s−1 at 2567 K with quenching. After these moments, the nucleation rates decrease drastically.
Figure 10 displays the nucleus diameter evolution. The nucleus diameters take minimum values of 1.15 nm without quenching and 0.916 nm with quenching. The temperatures of those minimum values coincide with temperatures at which the homogeneous nucleation rates reach the highest values. Those results of Figure 9 and Figure 10 indicate that a larger number of smaller nuclei are generated under the quenching condition where the cooling rate is simply higher. The total numbers of the vapor atoms are the same with and without quenching. Then they are consumed mainly by heterogeneous condensation on the nuclei as presented in Figure 8. Therefore, because a larger number of the nuclei share the vapor atoms with quenching, fewer vapor atoms condense on a nucleus. As a result, the quenching condition will yield a larger number of smaller nanoparticles.
Figure 11 portrays the instantaneous growth rates by condensation for 1 K decrease for particle diameters around the nucleus diameters when the homogeneous nucleation rates reach the highest values. The growth rate without quenching is 0.16–4.15 nm K−1, which is much higher than that with quenching of 0.03–0.33 nm K−1. Moreover, the condition without quenching takes 323 μs, which is much longer than that with quenching of 16.1 μs for a 1 K decrease. Therefore, under the condition without quenching, the particles grow larger, whereas smaller nuclei are generated continuously over a longer period. Therefore, the condition without quenching produces a more dispersed size distribution with a larger standard deviation.

4.3. Evolution of Particle Size Distribution and Cooling Rate Dependency

The evolutions of particle size distributions with and without quenching are presented from Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 for representative temperatures. The corresponding movies are also presented in Video S1 as supplementary data. Therein, τ denotes the time elapsed after the homogeneous nucleation rate exceeded 1.0 × 10 6 m−3s−1. It is noteworthy that one text [77] states that “particle formation can be conveniently observed experimentally” for this value. The conditions with and without quenching have different τ at the same temperature because the rates of temperature decrease differ.
Figure 12a shows that, at 2586 K, when the condition without quenching exhibits the highest nucleation rate, particles smaller than 5 nm present the largest number, whereas particles larger than 100 nm have already grown. On the other hand, the condition with quenching has only particles smaller than 5 nm as presented in Figure 12b.
Figure 13a shows that, at 2583 K, when the homogeneous nucleation rate is decreasing drastically, the particles have a bimodal distribution with two peaks under the condition without quenching because the particle growth changes from the nucleation dominant mode to the condensation dominant mode as indicated in Figure 8 and Figure 9. For Figure 13b, the particle growth is still in the nucleation dominant mode under the condition with quenching.
Figure 14b shows that, at 2567 K, when the condition with quenching offers the highest nucleation rate, the particles are still smaller than 30 nm. Particles smaller than 5 nm have the largest number. On the other hand, in Figure 14a, most particles have grown to larger particles under the condition without quenching. Particles smaller than 5 nm remain because the growth rate by heterogeneous condensation is lower for smaller particles because of their larger Knudsen numbers as indicated by Equation (12) and depicted in Figure 11. Those smaller nanoparticles will collide and merge with larger nanoparticles more frequently, thereby decreasing their numbers slowly.
Figure 15a,b show unimodal size distributions without and with quenching at τ = 163.75 ms and τ = 8.1873 ms, respectively. At these moments, the temperatures are 2100 K identically and 99% of the vapor atoms have already been converted into nanoparticles as shown in Figure 7. Therefore, after these moments, the nanoparticles grow by interparticle coagulation.
Figure 16a,b show size distributions obtained without and with quenching at 1700 K which is the melting point of bulk silicon. As depicted in Figure 3, the melting point of silicon nanoparticles are lower than bulk silicon. Therefore, the nanoparticles are still liquid and keep growing by interparticle coagulation. Without quenching, the product consists of particles large than 40 nm as shown in Figure 16a. Those nanoparticles will solidify and finish coagulation growth soon because the melting points for large nanoparticles are slightly lower than bulk silicon. On the other hand, with quenching, the product still includes many particles smaller than 40 nm as shown in Figure 16b. Due to the considerable depression of the melting points, their coagulation growth will progress longer. In addition, Figure 3 tells that the growing particles will start to solidify from large ones to small ones and eventually reach their final states. It is noteworthy that these size distributions in Figure 16 at 1700 K are not much different from those in Figure 15 at 2100 K, which means that growth by interparticle coagulation is much slower than that by heterogeneous condensation.
Figure 17, Figure 18 and Figure 19 depict the cooling rate dependence on the total number density, the number mean diameter, and the standard deviation of the finally obtained silicon nanoparticles, respectively. A higher cooling rate can lead to higher total number density, smaller size, and smaller standard deviation. In an actual fabrication process, attention must be devoted to the fact that cooling rate control changes those obtained quantities simultaneously. It is therefore interesting that those quantities have no linear dependence on the cooling rate. These results suggest that quenching in thermal plasma fabrication is effectual, but it has limitations for controlling nanoparticle characteristics.

5. Conclusions

Quenching effects on the growth process and size distribution of silicon nanoparticles in the typical range of cooling rates at a thermal plasma tail were investigated computationally using a nodal-type model. Results suggest that the size distribution evolves temporally with simultaneous homogeneous nucleation, heterogeneous condensation, interparticle coagulation, and melting point depression. Furthermore, silicon nanoparticle fabrication was demonstrated experimentally with and without quenching. Those experiment results support the validity of the numerically obtained size distributions.
The implicit mechanism of the collective growth processes of silicon nanoparticles was revealed. When the vapor atoms became highly supersaturated, they decreased drastically by rapid conversion of the vapor atoms to nanoparticles. At these rapid conversions, approximately 40% and 50% of the vapor atoms were converted to nanoparticles. When the number of nuclei increased by nucleation and the total surface area consequently increased, vapor consumption by condensation became overwhelming. A larger number of smaller nuclei were generated under the quenching condition where the cooling rate was higher. Therein, fewer vapor atoms condensed on a nucleus because more nuclei shared the vapor atoms. As a result, the quenching condition yielded a larger number of smaller nanoparticles. After most of the vapor atoms were converted into particles, the particles grew by interparticle coagulation much more slowly than heterogeneous condensation.
As the cooling rate was increased, the total number density increased, the size decreased, and the standard deviation increased. Control of the cooling rate changed those quantities simultaneously. Furthermore, those quantities were found to have no linear dependence with the cooling rate. These results suggest that quenching in thermal plasma fabrication was effectual but that it had limitations for controlling nanoparticle characteristics. A sensitivity analysis for the individual effects on those quantities in respect to the cooling rate would be further exploration of the observed “no linear dependence”.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nano11061370/s1, Video S1: Evolutions of size distributions of silicon nanoparticles (a) without and (b) with quenching. t0 denotes the time at which the homogeneous nucleation rate exceeded 1.0 × 10 6 m−3s−1: t0 = 29.8032 ms without quenching and t0 = 1.4902 ms with quenching. τ is defined as the time elapsed after t0.

Author Contributions

M.S. developed the numerical model and algorithm and designed the computational investigation. Y.H. conducted the experiment. E.G. analyzed the data and discussed the results. M.S. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Scientific Research (B) (KAKENHI: Grant No. 18H01371).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.W. Synthesis and properties of nanophase materials. Mater. Sci. Eng. A 1993, 168, 189–197. [Google Scholar] [CrossRef] [Green Version]
  2. Bapat, A.; Perrey, C.R.; Campbell, S.A.; Carter, C.B.; Kortshagen, U. Synthesis of highly oriented, single-crystal silicon nanoparticles in a low-pressure, inductively coupled plasma. J. Appl. Phys. 2003, 94, 1969–1974. [Google Scholar] [CrossRef]
  3. Bapat, A.; Anderson, C.; Perrey, C.R.; Carter, C.B.; Campbell, S.A.; Kortshagen, U. Plasma synthesis of single-crystal silicon nanoparticles for novel electronic device applications. Plasma Phys. Control. Fusion 2004, 46, B97–B109. [Google Scholar] [CrossRef] [Green Version]
  4. Fu, Y.; Dutta, A.; Willander, M.; Oda, S. Carrier conduction in a Si-nanocrystal-based single-electron transistor-I. Effect of gate bias. Superlattices Microstruct. 2000, 28, 177–187. [Google Scholar] [CrossRef]
  5. Nishiguchi, K.; Oda, S. Electron transport in a single silicon quantum structure using a vertical silicon probe. J. Appl. Phys. 2000, 88, 4186–4190. [Google Scholar] [CrossRef] [Green Version]
  6. Banerjee, S.; Huang, S.; Yamanaka, T.; Oda, S. Evidence of storing and erasing of electrons in a nanocrystalline-Si based memory device at 77 K. J. Vac. Sci. Technol. B Microelectron. Nanometer Struct. Process. Meas. Phenom. 2002, 20, 1135–1138. [Google Scholar] [CrossRef]
  7. Ostraat, M.L.; De Blauwe, J.W.; Green, M.L.; Bell, L.D.; Brongersma, M.L.; Casperson, J.; Flagan, R.C.; Atwater, H.A. Synthesis and characterization of aerosol silicon nanocrystal nonvolatile floating-gate memory devices. Appl. Phys. Lett. 2001, 79, 433–435. [Google Scholar] [CrossRef]
  8. Tiwari, S.; Rana, F.; Hanafi, H.; Hartstein, A.; Crabb’e, E.F.; Chan, K. A silicon nanocrystals based memory. Appl. Phys. Lett. 1996, 68, 1377–1379. [Google Scholar] [CrossRef] [Green Version]
  9. Tiwari, S.; Rana, F.; Chan, K.; Shi, L.; Hanafi, H. Single charge and confinement effects in nano-crystal memories. Appl. Phys. Lett. 1996, 69, 1232–1234. [Google Scholar] [CrossRef]
  10. Botas, A.M.P.; Brites, C.D.S.; Wu, J.; Kortshagen, U.; Pereira, R.N.; Carlos, L.D.; Ferreira, R.A.S. A new generation of primary luminescent thermometers based on silicon nanoparticles and operating in different media. Part. Part. Syst. Charact. 2016, 33, 740–748. [Google Scholar] [CrossRef]
  11. Larcher, D.; Beattie, S.; Morcrette, M.; Edström, K.; Jumasc, J.-C.; Tarascon, J.-M. Recent findings and prospects in the field of pure metals as negative electrodes for Li-ion batteries. J. Mater. Chem. 2007, 17, 3759–3772. [Google Scholar] [CrossRef]
  12. Liu, X.H.; Zhong, L.; Huang, S.; Mao, S.X.; Zhu, T.; Huang, J.Y. Size-dependent fracture of silicon nanoparticles during lithiation. ACS Nano 2012, 6, 1522–1531. [Google Scholar] [CrossRef]
  13. Zuo, X.; Zhu, J.; Müller-Buschbaum, P.; Cheng, Y.-J. Silicon based lithium-ion battery anodes: A chronicle perspective review. Nano Energy 2017, 31, 113–143. [Google Scholar] [CrossRef]
  14. Kambara, M.; Kitayama, A.; Homma, K.; Hideshima, T.; Kaga, M.; Sheem, K.-Y.; Ishida, S.; Yoshida, T. Nano-composite Si particle formation by plasma spraying for negative electrode of Li ion batteries. J. Appl. Phys. 2014, 115, 143302. [Google Scholar] [CrossRef]
  15. Ohta, R.; Fukada, K.; Tashiro, T.; Dougakiuchi, M.; Kambara, M. Effect of PS-PVD production throughput on Si nanoparticles for negative electrode of lithium ion batteries. J. Phys. D Appl. Phys. 2018, 51, 105501. [Google Scholar] [CrossRef]
  16. Tanaka, Y.; Shimizu, K.; Akashi, K.; Onda, K.; Uesugi, Y.; Ishijima, T.; Watanabe, S.; Sueyasu, S.; Nakamura, K. High rate synthesis of graphene-encapsulated silicon nanoparticles using pulse-modulated induction thermal plasmas with intermittent feedstock feeding. Jpn. J. Appl. Phys. 2020, 59, SHHE07. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, X.; Hayashida, R.; Tanaka, M.; Watanabe, T. Synthesis of carbon-coated silicon nanoparticles by induction thermal plasma for lithium ion battery. Powder Technol. 2020, 371, 26–36. [Google Scholar] [CrossRef]
  18. Zhang, X.; Liu, Z.; Tanaka, M.; Watanabe, T. Formation mechanism of amorphous silicon nanoparticles with additional counter-flow quenching gas by induction thermal plasma. Chem. Eng. Sci. 2020, 230, 116217. [Google Scholar] [CrossRef]
  19. Boulos, M.I.; Fauchais, P.; Pfender, E. Thermal Plasmas Fundamentals and Applications; Plenum Press: New York, NY, USA, 1994; Volume 1, Available online: https://www.springer.com/gp/book/9780306446078 (accessed on 20 May 2021).
  20. Sato, T.; Shigeta, M.; Kato, D.; Nishiyama, H. Mixing and magnetic effects on a nonequilibrium argon plasma jet. Int. J. Therm. Sci. 2001, 40, 273–278. [Google Scholar] [CrossRef]
  21. Shigeta, M. Numerical study of axial magnetic effects on a turbulent thermal plasma jet for nanopowder production using 3D time-dependent simulation. J. Flow Control Meas. Vis. 2018, 6, 107–123. [Google Scholar] [CrossRef] [Green Version]
  22. Shigeta, M.; Murphy, A.B. Thermal plasmas for nanofabrication. J. Phys. D Appl. Phys. 2011, 44, 174025. [Google Scholar] [CrossRef]
  23. Tanaka, Y.; Tsuke, T.; Guo, W.; Uesugi, Y.; Ishijima, T.; Watanabe, S.; Nakamura, K. A large amount synthesis of nanopowder using modulated induction thermal plasmas synchronized with intermittent feeding of raw materials. J. Phys. Conf. Ser. 2012, 406, 012001. [Google Scholar] [CrossRef] [Green Version]
  24. Kodama, N.; Tanaka, Y.; Kita, K.; Uesugi, Y.; Ishijima, T.; Watanabe, S.; Nakamura, K. A method for large-scale synthesis of Al-doped TiO2 nanopowder using pulse-modulated induction thermal plasmas with time-controlled feedstock feeding. J. Phys. D Appl. Phys. 2014, 47, 195304. [Google Scholar] [CrossRef] [Green Version]
  25. Shigeta, M.; Nishiyama, H. Numerical analysis of metallic nanoparticle synthesis using RF inductively coupled plasma flows. Trans. ASME J. Heat Transf. 2005, 127, 1222–1230. [Google Scholar] [CrossRef]
  26. Shigeta, M.; Tanaka, M.; Ghedini, E. Numerical analysis of correlation between arc plasma fluctuation and nanoparticle growth-transport under atmospheric pressure. Nanomaterials 2019, 9, 1736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Asai, S.; Tanaka, M.; Yamane, S.; Miyasaka, F.; Shigeta, M.; Nomura, K.; Ogino, Y. Recent progresses of welding and joining engineering. J. Jpn. Weld. Soc. 2020, 89, 322–335. (In Japanese) [Google Scholar] [CrossRef]
  28. Shigeta, M. Modeling and simulation of a turbulent-like thermal plasma jet for nanopowder production. IEEJ Trans. Electr. Electron. Eng. 2019, 14, 16–28. [Google Scholar] [CrossRef] [Green Version]
  29. Shigeta, M. Simulating turbulent thermal plasma flows for nanopowder fabrication. Plasma Chem. Plasma Process. 2020, 40, 775–794. [Google Scholar] [CrossRef]
  30. Shigeta, M. Turbulence modelling of thermal plasma flows. J. Phys. D Appl. Phys. 2016, 49, 493001. [Google Scholar] [CrossRef]
  31. Watanabe, T.; Okumiya, H. Formation mechanism of silicide nanoparticles by induction thermal plasmas. Sci. Technol. Adv. Mater. 2004, 5, 639–646. [Google Scholar] [CrossRef]
  32. Kambara, M.; Hideshima, T.; Kaga, M.; Yoshida, T. Powders: Plasma Spray PVD for High-Throughput Production, Encyclopedia of Plasma Technology, 1st ed.; Shohet, J.L., Ed.; Taylor & Francis Group: Boca Raton, FL, USA, 2017; pp. 1176–1190. Available online: https://www.taylorfrancis.com/books/mono/10.1081/E-EPLT/encyclopedia-plasma-technology-leon-shohet (accessed on 20 May 2021).
  33. Ishisaka, Y.; Kodama, N.; Kita, K.; Tanaka, T.; Uesugi, Y.; Ishijima, T.; Sueyasu, S.; Watanabe, S.; Nakamura, K. High-rate synthesis of Si nanowires using modulated induction thermal plasmas. Appl. Phys. Express 2017, 10, 096201. [Google Scholar] [CrossRef]
  34. Kong, P.C.; Lau, Y.C. Plasma synthesis of ceramic powders. Pure Appl. Chem. 1990, 62, 1809–1816. [Google Scholar] [CrossRef]
  35. Leparoux, M.; Schreuders, C.; Fauchais, P. Improved plasma synthesis of Si-nanopowders by quenching. Adv. Eng. Mater. 2008, 10, 1147–1150. [Google Scholar] [CrossRef]
  36. Li, J.G.; Wang, X.; Tang, C.; Ishigaki, T.; Tanaka, S. Energy transfer enables 1.53 μm photoluminescence from erbium-doped TiO2 semiconductor nanocrystals synthesized by Ar/O2 radio-frequency thermal plasma. J. Am. Ceram. Soc. 2008, 91, 2032–2035. [Google Scholar] [CrossRef]
  37. Proulx, P.; Bilodeau, J.-F. A model for ultrafine powders production in a thermal plasma reactor. Plasma Chem. Plasma Process. 1991, 11, 371–386. [Google Scholar] [CrossRef]
  38. Girshick, S.L.; Chiu, C.-P.; Muno, R.; Wu, C.Y.; Yang, L.; Singh, S.K.; McMurry, P.H. Thermal plasma synthesis of ultrafine iron particles. J. Aerosol Sci. 1993, 24, 367–382. [Google Scholar] [CrossRef]
  39. Bilodeau, J.F.; Proulx, P. A mathematical model for ultrafine iron powder growth in thermal plasma. Aerosol Sci. Technol. 1996, 24, 175–189. [Google Scholar] [CrossRef]
  40. Désilets, M.; Bilodeau, J.F.; Proulx, P. Modelling of the reactive synthesis of ultra-fine powders in a thermal plasma reactor. J. Phys. D Appl. Phys. 1997, 30, 1951–1960. [Google Scholar] [CrossRef]
  41. Cruz, A.C.D.; Munz, R.J. Vapor phase synthesis of fine particles. IEEE Trans. Plasma Sci. 1997, 25, 1008–1016. [Google Scholar] [CrossRef]
  42. Aristizabal, F.; Munz, R.J.; Berk, D. Modeling of the production of ultra fine aluminium particles in rapid quenching turbulent flow. J. Aerosol Sci. 2006, 37, 162–186. [Google Scholar] [CrossRef]
  43. Mendoza-Gonzalez, N.Y.; Goortani, B.M.; Proulx, P. Numerical simulation of silica nanoparticles production in an RF plasma reactor: Effect of quench. Mater. Sci. Eng. C 2007, 27, 1265–1269. [Google Scholar] [CrossRef]
  44. Shigeta, M.; Watanabe, T. Numerical investigation of cooling effect on platinum nanoparticle formation in inductively coupled thermal plasmas. J. Appl. Phys. 2008, 103, 074903. [Google Scholar] [CrossRef]
  45. Shigeta, M.; Watanabe, T. Two-dimensional analysis of nanoparticle formation in induction thermal plasmas with counterflow cooling. Thin Solid Films 2008, 516, 4415–4422. [Google Scholar] [CrossRef]
  46. Mendoza-Gonzalez, N.Y.; Morsli, M.E.; Proulx, P. Production of nanoparticles in thermal plasmas: A model including evaporation, nucleation, condensation, and fractal aggregation. J. Therm. Spray Technol. 2008, 17, 533–550. [Google Scholar] [CrossRef]
  47. Shigeta, M.; Watanabe, T. Model integration for metal nanoparticle synthesis by an rf thermal plasma flow with counterflow cooling. Trans. Jpn. Soc. Mech. Eng. Part B 2009, 75, 2019–2028. (In Japanese) [Google Scholar] [CrossRef]
  48. Girshick, S.L.; Chiu, C.-P.; McMurry, P.H. Modelling particle formation and growth in a plasma synthesis reactor. Plasma Chem. Plasma Process. 1988, 8, 145–157. [Google Scholar] [CrossRef]
  49. Girshick, S.L.; Chiu, C.-P. Numerical study of MgO powder synthesis by thermal plasma. J. Aerosol Sci. 1990, 21, 641–650. [Google Scholar] [CrossRef]
  50. Joshi, S.V.; Liang, Q.; Park, J.Y.; Batdorf, J.A. Effect of quenching conditions on particle formation and growth in thermal plasma synthesis of fine powders. Plasma Chem. Plasma Process. 1990, 10, 339–358. [Google Scholar] [CrossRef]
  51. Shigeta, M.; Watanabe, T.; Nishiyama, H. Numerical investigation for nano-particle synthesis in an RF inductively coupled plasma. Thin Solid Films 2004, 457, 192–200. [Google Scholar] [CrossRef]
  52. Hirayama, Y.; Suzuki, K.; Yamaguchi, W.; Takagi, K. Cold welding behavior of fine bare aluminum powders prepared by new low oxygen induction thermal plasma system. J. Alloys Compd. 2018, 768, 608–612. [Google Scholar] [CrossRef]
  53. Hirayama, Y.; Shigeta, M.; Liu, Z.; Yodoshi, N.; Hosokawa, A.; Takagi, K. Anisotropic Nd-Fe ultrafine particles with stable and metastable phases prepared by induction thermal plasma. J. Alloys Compd. 2021, 873, 159724. [Google Scholar] [CrossRef]
  54. Shigeta, M.; Watanabe, T. Numerical analysis for preparation of silicon-based intermetallic nano-particles in induction thermal plasma flow systems. JSME Int. J. Ser. B Fluids Therm. Eng. 2005, 48, 425–431. [Google Scholar] [CrossRef] [Green Version]
  55. Shigeta, M.; Watanabe, T. Numerical analysis for co-condensation processes in silicide nanoparticle synthesis using induction thermal plasma at atmospheric pressure conditions. J. Mater. Res. 2005, 20, 2801–2811. [Google Scholar] [CrossRef]
  56. Shigeta, M.; Watanabe, T. Multi-component co-condensation model of Ti-based boride/silicide nanoparticle growth in induction thermal plasmas. Thin Solid Films 2007, 515, 4217–4227. [Google Scholar] [CrossRef]
  57. Nemchinsky, V.A.; Shigeta, M. Simple equations to describe aerosol growth. Model. Simul. Mater. Sci. Eng. 2012, 20, 045017. [Google Scholar] [CrossRef]
  58. Shigeta, M. Simple nonequilibrium model of collective growth and transport of metal nanomist in a thermal plasma process. Theor. Appl. Mech. Jpn. 2015, 63, 147–154. [Google Scholar] [CrossRef]
  59. Murphy, A.B. Formation of titanium nanoparticles from a titanium tetrachloride plasma. J. Phys. D Appl. Phys. 2004, 37, 2841–2847. [Google Scholar] [CrossRef]
  60. Watanabe, T.; Shigeta, M. Nanoparticle synthesis by thermal plasmas. In Nanoparticles: Properties, Classification, Characterization, and Fabrication; Kestell, A.E., DeLorey, G.T., Eds.; Nova Science Publishers: New York, NY, USA, 2010; pp. 149–212. Available online: http://www.novapublishers.org/catalog/product_info.php?products_id=18755 (accessed on 20 May 2021).
  61. Watanabe, T.; Shigeta, M. Nanoparticle synthesis by thermal plasmas. In Nanomaterials: Properties, Preparation and Processes; Cabral, V., Silva, R., Eds.; Nova Science Publishers: New York, NY, USA, 2011; pp. 115–176. Available online: http://www.novapublishers.org/catalog/product_info.php?products_id=11125 (accessed on 20 May 2021).
  62. Tanaka, Y.; Onda, K.; Akashi, K.; Furukawa, R.; Nakano, Y.; Ishijima, T.; Uesugi, Y.; Sueyasu, S.; Watanabe, S.; Nakamura, K. A numerical study on nanoparticle synthesis in pulse-modulated induction thermal plasmas with intermittent feedstock powder feeding by method of moment. In Proceedings of the Gaseous Electronics Conference, Online, 6 October 2020. GT3.00006. [Google Scholar]
  63. Rao, N.; Girshick, S.L.; Heberlein, J.; McMurry, P.; Jones, S.; Hansen, D.; Micheel, B. Nanoparticle formation using a plasma expansion process. Plasma Chem. Plasma Process. 1995, 15, 581–606. [Google Scholar] [CrossRef]
  64. Shigeta, M.; Watanabe, T. Growth mechanism of silicon-based functional nanoparticles fabricated by inductively coupled thermal plasmas. J. Phys. D Appl. Phys. 2007, 40, 2407–2419. [Google Scholar] [CrossRef]
  65. Shigeta, M.; Watanabe, T. Two-directional nodal model for co-condensation growth of multi-component nanoparticles in thermal plasma processing. J. Therm. Spray Technol. 2009, 18, 1022–1037. [Google Scholar] [CrossRef]
  66. Shigeta, M.; Watanabe, T. Growth model of binary alloy nanopowders for thermal plasma synthesis. J. Appl. Phys. 2010, 108, 043306. [Google Scholar] [CrossRef]
  67. Cheng, Y.; Shigeta, M.; Choi, S.; Watanabe, T. Formation mechanism of titanium boride nanoparticles by RF induction thermal plasma. Chem. Eng. J. 2012, 183, 483–491. [Google Scholar] [CrossRef]
  68. Shigeta, M.; Miyake, M.; Tanaka, M. Numerical analysis of collective growth of primary nanoparticles in arc welding. Q. J. Jpn. Weld. Soc. 2015, 33, 365–375. (In Japanese) [Google Scholar] [CrossRef] [Green Version]
  69. Shigeta, M.; Watanabe, T. Effect of precursor fraction on silicide nanopowder growth under thermal plasma conditions: A computational study. Powder Technol. 2016, 288, 191–201. [Google Scholar] [CrossRef]
  70. Shigeta, M.; Watanabe, T. Effect of saturation pressure difference on metal-silicide nanopowder formation in thermal plasma fabrication. Nanomaterials 2016, 6, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Colombo, V.; Ghedini, E.; Gherardi, M.; Sanibondi, P.; Shigeta, M. A two-dimensional nodal model with turbulent effects for the synthesis of Si nano-particles by inductively coupled thermal plasmas. Plasma Sources Sci. Technol. 2012, 21, 025001. [Google Scholar] [CrossRef]
  72. Bora, B.; Saikia, B.J.; Borgohain, C.; Kakati, M.; Das, A.K. Numerical investigation of nanoparticle synthesis in supersonic thermal plasma expansion. Vacuum 2010, 85, 283–289. [Google Scholar] [CrossRef]
  73. Shigeta, M.; Minami, S.; Tanaka, M. Modeling for collective growth of fume primary particles with charge effect in arc welding. Weld. World 2018, 62, 203–213. [Google Scholar] [CrossRef]
  74. Girshick, S.L.; Chiu, P.C.; McMurry, P.H. Time-dependent aerosol models and homogeneous nucleation rates. Aerosol Sci. Technol. 1990, 13, 465–477. [Google Scholar] [CrossRef]
  75. Seinfeld, J.H.; Pandis, S.N. Atmospheric Chemistry and Physics, From Air Pollution to Climate Change; Wiley: New York, NY, USA, 1998. [Google Scholar]
  76. Hirschfelder, J.O.; Curtiss, C.F.; Bird, R.B. Molecular Theory of Gases and Liquids; Wiley: New York, NY, USA, 1954. [Google Scholar]
  77. Friedlander, S.K. Smoke, Dust and Haze: Fundamentals of Aerosol Dynamics, 2nd ed.; Oxford University Press: New York, NY, USA, 2000. [Google Scholar]
  78. Guisbiers, G.; Kazan, M.; Overschelde, O.V.; Wautelet, M.; Pereira, S. Mechanical and thermal properties of metallic and semiconductive nanostructures. J. Phys. Chem. 2008, 112, 4097–4103. [Google Scholar] [CrossRef]
  79. The Japan Institute of Metals and Materials. Metal Data Book; Maruzen: Tokyo, Japan, 1993. (In Japanese) [Google Scholar]
  80. Kawajiri, K.; Sato, T.; Nishiyama, H. Experimental analysis of a DC–RF hybrid plasma flow. Surf. Coat. Technol. 2003, 171, 134–139. [Google Scholar] [CrossRef]
  81. Takana, H.; Jang, J.; Igawa, J.; Nakajima, T.; Solonenko, O.P.; Nishiyama, H. Improvement of in-flight alumina spheroidization process using a small power argon DC-RF hybrid plasma flow system by helium mixture. J. Therm. Spray Technol. 2011, 20, 432–439. [Google Scholar] [CrossRef]
  82. Shigeta, M. Time-dependent 3-D simulation of an argon RF inductively coupled thermal plasma. Plasma Sources Sci. Technol. 2021, 21, 055029. [Google Scholar] [CrossRef]
  83. Shigeta, M. Numerical investigation of nano-material processing by thermal plasma flows. In Sustained Simulation Performance 2012; Resch, M.M., Wang, X., Bez, W., Focht, E., Kobayashi, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 169–182. [Google Scholar] [CrossRef]
  84. Shigeta, M. Three-dimensional flow dynamics of an argon RF plasma with dc jet assistance: A numerical study. J. Phys. D Appl. Phys. 2013, 46, 015401. [Google Scholar] [CrossRef]
  85. Shigeta, M.; Sato, T.; Nishiyama, H. Numerical simulation of a potassium-seeded turbulent RF inductively coupled plasma with particles. Thin Solid Films 2003, 435, 5–12. [Google Scholar] [CrossRef]
  86. Shigeta, M.; Sato, T.; Nishiyama, H. Computational simulation of a particle-laden RF inductively coupled plasma with seeded potassium. Int. J. Heat Mass Transf. 2004, 47, 707–716. [Google Scholar] [CrossRef]
  87. Xiang, J.; Park, H.; Tanaka, K.; Shigeta, M.; Tanaka, M.; Murphy, A.B. Numerical study of the effects and transport mechanisms of iron vapour in tungsten inert-gas welding in argon. J. Phys. D Appl. Phys. 2020, 53, 044004. [Google Scholar] [CrossRef]
  88. Xiang, J.; Chen, F.F.; Park, H.; Tanaka, K.; Shigeta, M.; Tanaka, M.; Murphy, A.B. Numerical study of the metal vapour transport in tungsten inert-gas welding in argon for stainless steel. Appl. Math. Model. 2020, 79, 713–728. [Google Scholar] [CrossRef]
  89. Tanaka, K.; Shigeta, M.; Tanaka, M.; Murphy, A.B. Imaging spectroscopy for transient transport of chromium vapor during helium TIG welding. Q. J. Jpn. Weld. Soc. 2020, 38, 21s–24s. [Google Scholar] [CrossRef]
  90. Tanaka, K.; Shigeta, M.; Tanaka, M.; Murphy, A.B. Investigation of transient metal vapour transport processes in helium arc welding by imaging spectroscopy. J. Phys. D Appl. Phys. 2020, 53, 425202. [Google Scholar] [CrossRef]
Figure 1. Schematic illustrations of nanoparticle fabrication in ICTP systems (a) without and (b) with quenching.
Figure 1. Schematic illustrations of nanoparticle fabrication in ICTP systems (a) without and (b) with quenching.
Nanomaterials 11 01370 g001
Figure 2. Illustrations depicting computations of (a) homogeneous nucleation, (b) heterogeneous condensation, and (c) interparticle coagulation.
Figure 2. Illustrations depicting computations of (a) homogeneous nucleation, (b) heterogeneous condensation, and (c) interparticle coagulation.
Nanomaterials 11 01370 g002
Figure 3. Melting points of bulk silicon and nano-size spherical silicon estimated by Equation (17).
Figure 3. Melting points of bulk silicon and nano-size spherical silicon estimated by Equation (17).
Nanomaterials 11 01370 g003
Figure 4. SEM images of silicon nanoparticles obtained through experiments (a) without and (b) with quenching.
Figure 4. SEM images of silicon nanoparticles obtained through experiments (a) without and (b) with quenching.
Nanomaterials 11 01370 g004
Figure 5. Size distributions of silicon nanoparticles obtained from (a) experiments and (b) computations.
Figure 5. Size distributions of silicon nanoparticles obtained from (a) experiments and (b) computations.
Nanomaterials 11 01370 g005
Figure 6. Evolutions of molar densities of vapor atoms and saturated states.
Figure 6. Evolutions of molar densities of vapor atoms and saturated states.
Nanomaterials 11 01370 g006
Figure 7. Evolutions of conversion ratios.
Figure 7. Evolutions of conversion ratios.
Nanomaterials 11 01370 g007
Figure 8. Evolutions of vapor consumption rates by homogeneous nucleation and heterogeneous condensation.
Figure 8. Evolutions of vapor consumption rates by homogeneous nucleation and heterogeneous condensation.
Nanomaterials 11 01370 g008
Figure 9. Evolutions of homogeneous nucleation rates.
Figure 9. Evolutions of homogeneous nucleation rates.
Nanomaterials 11 01370 g009
Figure 10. Evolutions of nucleus diameters.
Figure 10. Evolutions of nucleus diameters.
Nanomaterials 11 01370 g010
Figure 11. Instantaneous growth rates by heterogeneous condensation at temperatures of the highest homogeneous nucleation rate (a) without and (b) with quenching.
Figure 11. Instantaneous growth rates by heterogeneous condensation at temperatures of the highest homogeneous nucleation rate (a) without and (b) with quenching.
Nanomaterials 11 01370 g011
Figure 12. Size distributions of silicon nanoparticles at 2586 K (a) without and (b) with quenching.
Figure 12. Size distributions of silicon nanoparticles at 2586 K (a) without and (b) with quenching.
Nanomaterials 11 01370 g012
Figure 13. Size distributions of silicon nanoparticles at 2583 K (a) without and (b) with quenching.
Figure 13. Size distributions of silicon nanoparticles at 2583 K (a) without and (b) with quenching.
Nanomaterials 11 01370 g013
Figure 14. Size distributions of silicon nanoparticles at 2567 K (a) without and (b) with quenching.
Figure 14. Size distributions of silicon nanoparticles at 2567 K (a) without and (b) with quenching.
Nanomaterials 11 01370 g014
Figure 15. Size distributions of silicon nanoparticles at 2100 K (a) without and (b) with quenching.
Figure 15. Size distributions of silicon nanoparticles at 2100 K (a) without and (b) with quenching.
Nanomaterials 11 01370 g015
Figure 16. Size distributions of silicon nanoparticles at 1700 K (a) without and (b) with quenching.
Figure 16. Size distributions of silicon nanoparticles at 1700 K (a) without and (b) with quenching.
Nanomaterials 11 01370 g016
Figure 17. Cooling rate dependence on total number density of silicon nanoparticles.
Figure 17. Cooling rate dependence on total number density of silicon nanoparticles.
Nanomaterials 11 01370 g017
Figure 18. Cooling rate dependence on number mean diameter of silicon nanoparticles.
Figure 18. Cooling rate dependence on number mean diameter of silicon nanoparticles.
Nanomaterials 11 01370 g018
Figure 19. Cooling rate dependence on standard deviation of silicon nanoparticles.
Figure 19. Cooling rate dependence on standard deviation of silicon nanoparticles.
Nanomaterials 11 01370 g019
Table 1. Computational growth models based on aerosol representation.
Table 1. Computational growth models based on aerosol representation.
Model Type A
[51,54,55,56]
Model Type B
[26,27,28,29,57,58]
Model Type C
[39,40,42,43,44,45,46,47,59,60,61,62]
Model Type D
[37,48,49,53,63,64,65,66,67,68,69,70,71,72,73]
Nucleationconsideredconsideredconsideredconsidered
Condensationconsideredconsideredconsideredconsidered
Coagulationnot consideredconsideredconsideredconsidered
Size distributionanymono-disperselognormalany
Mathematical descriptionsimplersimplermore complexmore complex
Computational costshigherlowerlowerhigher
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shigeta, M.; Hirayama, Y.; Ghedini, E. Computational Study of Quenching Effects on Growth Processes and Size Distributions of Silicon Nanoparticles at a Thermal Plasma Tail. Nanomaterials 2021, 11, 1370. https://doi.org/10.3390/nano11061370

AMA Style

Shigeta M, Hirayama Y, Ghedini E. Computational Study of Quenching Effects on Growth Processes and Size Distributions of Silicon Nanoparticles at a Thermal Plasma Tail. Nanomaterials. 2021; 11(6):1370. https://doi.org/10.3390/nano11061370

Chicago/Turabian Style

Shigeta, Masaya, Yusuke Hirayama, and Emanuele Ghedini. 2021. "Computational Study of Quenching Effects on Growth Processes and Size Distributions of Silicon Nanoparticles at a Thermal Plasma Tail" Nanomaterials 11, no. 6: 1370. https://doi.org/10.3390/nano11061370

APA Style

Shigeta, M., Hirayama, Y., & Ghedini, E. (2021). Computational Study of Quenching Effects on Growth Processes and Size Distributions of Silicon Nanoparticles at a Thermal Plasma Tail. Nanomaterials, 11(6), 1370. https://doi.org/10.3390/nano11061370

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