Kinetics of Silane Decomposition in High-Pressure ... - ACS Publications

Nov 20, 2017 - The resulting RMG-generated model reason- ably predicts silane concentration profiles at relatively low pressure for different temperat...
0 downloads 0 Views 1MB Size
Article Cite This: Ind. Eng. Chem. Res. 2017, 56, 14995−15000

pubs.acs.org/IECR

Kinetics of Silane Decomposition in High-Pressure Confined Chemical Vapor Deposition of Hydrogenated Amorphous Silicon Seyed Pouria Motevalian,*,† Stephen C. Aro,‡ Hiu Y. Cheng,‡ Todd D. Day,‡ Adri C. T. van Duin,§ John V. Badding,‡ and Ali Borhan† Department of Chemical Engineering, ‡Department of Chemistry, and §Department of Mechanical and Nuclear Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States

Ind. Eng. Chem. Res. 2017.56:14995-15000. Downloaded from pubs.acs.org by AUCKLAND UNIV OF TECHNOLOGY on 01/29/19. For personal use only.



S Supporting Information *

ABSTRACT: We study the kinetics of silane pyrolysis via confined high-pressure chemical vapor deposition (HPCVD) at pressures of 20−33 MPa in a microcapillary of 5.9 μm inner diameter. We find the growth rate to be first order with respect to silane concentration, with an activation energy of 53.7 ± 2.9 kcal/mol and a pre-exponential factor of 1.5 × 1010 m/s. The obtained activation energy is in the range of activation energies reported for hydrogen desorption from c-Si surfaces, suggesting that hydrogen desorption from the surface is the ratelimiting step in film growth. To further investigate this finding, reactive molecular dynamics simulations of thermal decomposition of silane on clean and hydrogen-passivated c-Si were performed. Homogeneous reactions were not observed in any of the simulations, in support of the hypothesis that heterogeneous silane decomposition on the silicon surface is the dominant mechanism for film deposition. In silane pyrolysis simulations on clean c-Si surfaces, almost all available silicon surface sites (i.e., dangling bonds) were occupied by silicon-hydrides (mostly tri- and dihydrides) upon exposure to gas-phase silane, whereas no reaction was observed during silane decomposition simulations on the hydrogen-passivated c-Si. Therefore, the results of the reactive molecular dynamics simulations indicate that the availability of dangling bonds resulting from hydrogen desorption from the surface is the rate-limiting step in film growth at high pressure.



INTRODUCTION Hydrogenated amorphous silicon (a-Si:H), one of the most technologically important semiconductors,1 is widely produced by the pyrolysis of silane (SiH4) at low pressure (mTorr) via chemical vapor deposition (CVD). The major challenge in synthesizing a-Si:H with CVD is to thermally decompose silane molecules at a low enough temperature to incorporate sufficient hydrogen to passivate the dangling bonds, while utilizing a high enough temperature to achieve rapid film growth. Because SiH4 precursor does not readily decompose until heated above 550 °C in a conventional CVD process, several methods have been developed to activate the precursor molecules to decompose at temperatures suitable for sufficient hydrogen incorporation.2 Examples include hot-wire CVD3−5 and plasma-enhanced CVD.6,7 Another method for accelerating reaction kinetics at low temperature is high-pressure CVD at 70−200 atm. This method allows for deposition rates of 0.3 Å/s at 450 °C, 30-fold higher than typical deposition rates of about 0.01 Å/s achieved in a conventional CVD reactor (with 13 Pa partial pressure of SiH4) at the same temperature.8,9 However, high pressure accelerates the homogeneous reactions that result in particle formation in the gas phase, leading to the deposition of lowquality a-Si:H films. Gas-phase reactions can be avoided by confining the reaction to microscale spaces, thereby resulting in the formation of smooth films.10 Recently, fabrication of singlecrystal silicon core fibers by laser crystallization of amorphous silicon deposited via high-pressure CVD has been reported.11 © 2017 American Chemical Society

The overall reaction for silane pyrolysis during confined HPCVD is

SiH4 → Si + 2H 2

(R1)

There are several homogeneous and heterogeneous elementary reactions involved in this process. In the gas phase, monosilane (silane) is converted into other (poly) silanes (SinH2n+2), silylenes (SiH2), and disilenes (double-bonded isomers of silylene with the general chemical formula SinH2n).12 The pioneering work of Purnell and Walsh13 proposed the conversion of silane to monosilylene and hydrogen in the gas phase as the initial step of the decomposition, followed by conversion of monosilylene to higher silanes (disilane and trisilane) according to the following reactions: SiH4 → SiH 2 + H 2

(R2)

SiH4 + SiH 2 ↔ Si 2H6

(R3)

SiH 2 + Si 2H6 ↔ Si3H8

(R4)

They postulated that silicon film formation proceeds via decomposition of silylenes at the surface. Received: Revised: Accepted: Published: 14995

August 25, 2017 November 16, 2017 November 20, 2017 November 20, 2017 DOI: 10.1021/acs.iecr.7b03515 Ind. Eng. Chem. Res. 2017, 56, 14995−15000

Article

Industrial & Engineering Chemistry Research

silane adsorbing at the surface to form mono- and trihydrides that subsequently decamp to form silicon at the surface. Gates et al.21,22 studied the kinetics of surface reactions in low-pressure CVD using static secondary ion mass spectroscopy. They performed silane pyrolysis experiments at temperatures below 500 °C and at very low pressure to minimize the contribution of homogeneous reactions to film growth. They observed three different silicon-hydride species (monohydride, dihydride, and trihydride) at the surface, and proposed the reaction mechanism shown in Table 2. For experiments below

The simple reaction mechanism proposed by Purnell and Walsh was later modified by Coltrin et al.14 by adding several more reaction intermediates (mostly disilenes) to increase the number of reaction steps to 26 reversible reactions, with the rate constant for each reaction step calculated using Rice− Ramsperger−Kassel−Marcus (RRKM) theory. They assumed that the reaction of intermediates striking the surface results in the formation of the silicon film. The difference between the reactivity of different intermediate species at the surface was considered by assigning surface reaction probabilities to different gas-phase reaction byproducts. They concluded that in CVD with hydrogen and silane precursor, only monosilane and monosilylene from the gas phase contribute to film growth, with negligible contribution from the latter at low temperatures. Later, Ho et al.15 proposed a gas-phase reaction mechanism consisting of 10 elementary steps based on the study of Coltrin et al., and parametrized the temperature and pressure dependencies of the rate constant for each step in Troe form. Their predicted values of the density of silicon particles formed in the gas phase agreed well with laser-induced fluorescence measurements. The silane pyrolysis model of Ho et al. has been widely used in the literature,16−18 and details of their gas-phase mechanism are provided in Table 1.

Table 2. Surface Reaction Mechanism Proposed by Gates et al.21,22 k

1 SiH4 + 2* → SiH*3 + H*

k2

SiH*3 + * → H* + SiH*2 k3

2SiH*2 → H 2 + 2SiH* k4

2SiH* → H 2 + 2Si* k5

2H* → H 2 + 2* k6

Si* → film + *

Table 1. Gas-Phase Reaction Mechanism Proposed by Ho et al.15 SiH4 ↔ SiH2 + H2 Si2H6 ↔ SiH4 + SiH2 Si2H6 ↔ H3SiSiH + H2 Si3H8 ↔ SiH2 + Si2H6 Si3H8 ↔ SiH4 + H3SiSiH H3SiSiH ↔ H2SiSiH2 H3SiSiH + H2 ↔ SiH4 + SiH2 H3SiSiH + SiH4 ↔ Si2H6 + SiH2 H3SiSiH ↔ Si + SiH4 Si + Si2H6 ↔ SiH2 + H3SiSiH

(1) (2) (3) (4) (5) (6)

500 °C, they reported an overall silane pyrolysis activation energy of 47 kcal/mol, a value that is close to the activation energy for hydrogen desorption from adsorbed monohydrides (reaction 4 in Table 2). Thus, they concluded that hydrogen desorption from the surface, resulting in the availability of new dangling bonds, is the rate-limiting step; this result is also supported by other studies.20,23−25 Whether the homogeneous or heterogeneous mechanism is predominantly responsible for silicon deposition is subject to debate. Some researchers believe homogeneous initiation of the growth process occurs wherein higher-order silanes (disilane and trisilane) and/or silylenes formed in the gas phase stick to the surface to deposit a semiconductor film.12−16,26 Others believe heterogeneous decomposition of silane at the surface occurs to deposit the film.20,23 The relative contributions of the homogeneous and heterogeneous reactions to film growth depend on deposition temperature, pressure, and the surfaceto-volume ratio of the reactor.27,28 In this study, we present the kinetics of film growth via confined HPCVD in a microcapillary, which results in the formation of optical fibers. We also show via reactive molecular dynamics simulations that in the confined HPCVD process, heterogeneous thermal decomposition of silane at the amorphous silicon surface is the dominant mechanism for deposition of the a-Si:H film due to the very high surface area to volume ratio present in the confined geometry. The kinetic data obtained in this study for a-Si:H film growth can be utilized to optimize the confined HPCVD process for fabrication of high-quality waveguide optical fibers with a plugged central channel over lengths of centimeters.9

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)

Odden et al.16 studied the kinetics of silane decomposition at relatively high temperature and pressure ranges of 690−830 K and 0.1−3.7 MPa, respectively. Using gas chromatography, they analyzed the residual gas after pyrolysis and found that it only contained monosilane and hydrogen, with traces of disilane in a few experiments. They adopted the reaction mechanism proposed by Ho et al., but added two more gas-phase reactions to regenerate monosilane from di- and trisilanes. They also modified the gas−surface interactions by considering only one dangling bond to be involved in the surface adsorption of silane in light of having a dense environment at high pressure. More recently, kinetic and thermodynamic data for silicon hydrides have been incorporated into Reaction Mechanism Generator (RMG), an open source software that generates kinetic models on the basis of a general understanding of how molecules react. The resulting RMG-generated model reasonably predicts silane concentration profiles at relatively low pressure for different temperatures.19 The heterogeneous decomposition of silane has been the subject of several studies. Robertson et al.20 performed several experiments at low pressure and concluded that silane pyrolysis on a silicon surface is completely surface-initiated. They also proposed that both hydrogen and silane adsorb at the surface by competing for free surface sites (i.e., dangling bonds), with



EXPERIMENTAL METHODS During each deposition experiment, a high-pressure (∼35 MPa) gas mixture of silane (less than 5%) and helium was introduced into bare-silica glass microcapillaries with 5.9 μm inner diameter. The microcapillaries were heated by a single furnace consisting of a ceramic casing wrapped in resistive heating wire and insulated with a thick (2−3 in.) layer of 14996

DOI: 10.1021/acs.iecr.7b03515 Ind. Eng. Chem. Res. 2017, 56, 14995−15000

Article

Industrial & Engineering Chemistry Research fiberglass. Wall temperatures throughout the furnace were measured via a k-type thermocouple housed in the ceramic casing using a hole of 1/8-in. outer diameter. This allowed for the thermocouple to have good thermal contact with the interior of the furnace at each measurement point. A second, identical thermocouple was used along with a PID feedback loop to control the temperature of the furnace via pulse-widthmodulation of the AC electric load (∼3A). A typical wall temperature profile is shown in Figure 1c. At the desired time

Cross sections of the sample at different axial positions along the microcapillary were mounted using standard electron microscopy techniques, and imaged using a FEI Nova NanoSEM 630 field-emission scanning electron microscope. Typical samples are shown in Figure 1a. This allowed measurement of the azimuthally averaged a-Si:H film thickness within each cross-section, which was then used to determine an axial profile for film thickness such as that shown in Figure 1b. Film thickness profiles were obtained by performing depositions at various precursor gas pressures and furnace temperatures.



RESULTS AND DISCUSSION The measured rate of a-Si:H film growth at different furnace temperatures and silane partial pressures is shown in Figure 2.

Figure 2. Film growth rate for different furnace temperatures and silane partial pressures.

The measured growth rates range from 0.008 to 0.272 Å/s, depending on the silane partial pressure and furnace temperature used in the experiment. Silicon film growth rate is an increasing function of both furnace temperature and silane partial pressure. For an axially uniform a-Si:H film, the time-evolution of the thickness of the deposited film (δ) can be related to the precursor pressure (P) and silane mole fraction (α) according to

Figure 1. (a) FESEM images of the a-Si:H film cross-section in a 5.9 μm-diameter microcapillary. (b) The axial profile of film thickness determined from FESEM images of the a-Si:H film cross-section after 24 h for a 2% silane mixture at 30 MPa. (c) The axial profile of furnace temperature.

⎛M P⎞ dδ = kα⎜ w ⎟ dt ⎝ ρRT ⎠

(1)

where t, k, R, T, Mw, and ρ denote the time, reaction rate constant, universal gas constant, temperature, molecular weight, and density of the a-Si:H film, respectively. Here, we have assumed the precursor to be an ideal gas mixture, and the overall silane decomposition reaction to be first-order with respect to silane concentration. Equation 1 was used to calculate the reaction rate constant (k) at different temperatures from the corresponding experimental measurements of film growth rate shown in Figure 2. Assuming an Arrhenius temperature-dependence for the reaction rate constant, the activation energy (E0) and pre-exponential factor (A0) were then obtained from a semilog plot of the reaction rate constant versus 1/T. The Arrhenius plot of the surface reaction rate constant versus 1/T is shown in Figure 3. The vertical error bars represent the combination of errors involved in the calculation

point, the deposition was quenched by removing the microcapillary from the furnace, venting it, and flushing it with high-pressure helium to purge the silane. The high-pressure helium and silane mixtures were prepared by first evacuating a stainless steel high-pressure vessel under high vacuum (9−10 Torr) via a turbomolecular pump (Pfeiffer HiPace 300) backed by a scroll pump (Edwards XDS series). Once evacuated, the desired amount of silane gas was cryogenically condensed into the stainless steel vessel from a pressure-monitored antechamber of known volume. Upon warming to room temperature, the silane pressure within the vessel was measured. The vessel was then pressurized to 35 MPa with helium carrier gas using a pneumatically actuated diaphragm gas compressor. Film deposition was subsequently carried out for sufficiently short times to avoid mass transfer limitations. 14997

DOI: 10.1021/acs.iecr.7b03515 Ind. Eng. Chem. Res. 2017, 56, 14995−15000

Article

Industrial & Engineering Chemistry Research

lations of thermal decomposition of silane at a silicon surface using ReaxFF.31−35 In particular, ReaxFF simulations were performed to provide additional support for the notion that hydrogen desorption from the surface is the rate-limiting step for film growth at high pressures in confined microcapillaries. ReaxFF is a bond-order-based potential that includes a geometry-dependent, polarizable charge as well as van der Waals interactions, which together enable the method to simulate both the covalent and the nonbonded interactions involved in a system comprised of a silicon surface and silane molecules. Additional details about ReaxFF and force-field testing are provided in the Supporting Information. Dynamic simulations of silane decomposition at a silicon surface were performed for a system comprised of 15 silane molecules (at 7 MPa pressure and 400 °C temperature in a 22 × 22 × 50 Å box) and a c-Si slab with 256 silicon atoms (as shown in the upper left-hand side panel embedded in Figure 4). Details of silane energy minimization and c-Si slab preparation via NVT simulation are provided in the Supporting Information. Starting with the initial configuration shown in the upper left-hand corner of Figure 4, canonical ensemble simulations of silane decomposition on clean and hydrogenpassivated c-Si surfaces were performed, and the number of unreacted silane molecules was monitored throughout the simulations. Homogeneous reactions were not observed during ReaxFF simulations of high-pressure silane pyrolysis on either clean or hydrogen-passivated surfaces, supporting the notion that, unlike silicon film growth in reactor geometries with very small surface-to-volume ratios, high-pressure CVD of silicon in confined domains is predominantly due to a heterogeneous reaction mechanism. Figure 4 shows the comparison between the simulation results for clean and hydrogen-passivated c-Si surfaces in terms of the percentage of unreacted silane molecules in the gas phase as a function of time. Starting from the initial configuration shown in the upper left-hand side panel, almost all silane molecules decompose at the clean c-Si surface to form siliconhydrides (mostly di- and trihydrides) within 10 ps of simulation time, as shown by the final configuration in the lower righthand side image in Figure 4. However, no reaction was observed on the hydrogen-passivated c-Si surface during NVT

Figure 3. Arrhenius plot for the surface reaction rate constant versus furnace temperature.

of the reaction rate constant (i.e., the errors associated with pressure and temperature variations, as well as azimuthal variations of the a-Si:H film thickness measured by fieldemission scanning electron microscopy), while the horizontal error bar represents temperature controller error of 0.5 °C. The slope of the linear regression fit in Figure 3 is proportional to the activation energy (E0) and yields a value of E0 = 53.7 ± 2.9 kcal/mol, while the intercept corresponds to a pre-exponential factor of 1.5 × 1010 m/s. The error bar for the activation energy is calculated on the basis of a 99% confidence interval. The R2 value of 0.99 for the linear regression shows that the experimentally measured reaction rate is well described by first-order kinetics. The predicted activation energy falls in the range of activation energies reported for hydrogen desorption reactions in heterogeneous silane pyrolysis (e.g., 45, 47, and 55 kcal/mol for reactions 3, 4, and 5 in Table 2, respectively).29,30 This suggests that hydrogen desorption from the surface to produce new dangling bonds is the rate-limiting step in HPCVD of a-Si:H within micrometer-sized confinements. This is consistent with the finding of previous kinetic studies at pressures below 1300 Pa, where hydrogen desorption was reported as the rate-limiting step for silicon film growth.20 Additional insight into the silane pyrolysis reaction was obtained by performing reactive molecular dynamics simu-

Figure 4. Percent of unreacted silane molecules in the gas phase as a function of time for ReaxFF simulations on clean and hydrogen-passivated c-Si surfaces. The initial configuration of the ReaxFF system is shown in the upper left-hand side panel. The final configurations for simulation with hydrogen-passivated and clean c-Si slabs are shown in the upper and lower right-hand side panels, respectively. 14998

DOI: 10.1021/acs.iecr.7b03515 Ind. Eng. Chem. Res. 2017, 56, 14995−15000

Article

Industrial & Engineering Chemistry Research simulations, as shown by the upper right-hand side panel in Figure 4. Therefore, the results of these ReaxFF simulations suggest that silane decomposition is mainly controlled by the availability of surface reaction sites (i.e., dangling bonds), which further supports the notion that hydrogen desorption from the surface to produce new dangling bonds is the rate-limiting step in confined HPCVD of a-Si:H. The experimental study of Odden et al.16 on silane pyrolysis at pressures of 0.1−3.7 MPa showed that concentrations of the reactive gas-phase intermediates, that is, Si2H6, Si3H8, and SiH2, were 2, 4, and 12 orders of magnitude smaller than silane concentration, respectively. The reaction rate constants for adsorption of Si2H6 and Si3H8 at the surface (reactions 2 and 3 in Table 2, respectively) are only 1 order of magnitude larger than that for adsorption of silane.27 Hence, assuming first-order kinetics, the rate of silane adsorption at the surface was 10 and 1000 times the adsorption rates of di- and trisilane, respectively, in the experiments of Odden et al.16 Therefore, heterogeneous decomposition of silane was the dominant source of the deposited a-Si:H film in their experiments, consistent with the results of our ReaxFF simulations. Geometric confinement can play an important role in determining the dominant film growth mechanism. Previous studies have demonstrated that decomposition of SiH4 carried out in reactors with centimeter dimensions leads to the formation of silicon particles in the gas phase that subsequently deposit at the substrate surface.16,36,37 It is possible, however, to avoid particle formation in the gas phase and promote heterogeneous film growth by confining the precursor in extreme aspect-ratio reactors, such as silica fiber capillaries, as reported previously.9,38,39 Therefore, confining HPCVD reactors to capillaries with diameters on the order of 100 μm or smaller can lead to the predominance of heterogeneous pathway in a-Si:H film growth.

contingent on the availability of surface dangling bonds resulting from hydrogen desorption from the surface.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.iecr.7b03515. Reactive molecular dynamic simulations using ReaxFF (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: seyedpouria.motevalian@pfizer.com. ORCID

Seyed Pouria Motevalian: 0000-0003-3768-789X John V. Badding: 0000-0002-4517-830X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the Penn State Materials Research Science and Engineering Center (NSF DMR-1420620). S.P.M. wishes to thank Mohammad Haghighi Parchini for insightful comments and valuable suggestions.



REFERENCES

(1) Street, R. A. Hydrogenated Amorphous Silicon; Cambridge University Press: New York, 1991. (2) Jasinski, J. M.; Gates, S. M. Silicon chemical vapor deposition one step at a time: fundamental studies of silicon hydride chemistry. Acc. Chem. Res. 1991, 24 (1), 9−15. (3) Wiesmann, H.; Ghosh, A.; McMahon, T.; Strongin, M. a-Si: H produced by high-temperature thermal decomposition of silane. J. Appl. Phys. 1979, 50 (5), 3752−3754. (4) Jadkar, S.; Sali, J.; Funde, A.; Bakr, N. A.; Vidyasagar, P.; Hawaldar, R.; Amalnerkar, D. Deposition of hydrogenated amorphous silicon (a-Si: H) films by hot-wire chemical vapor deposition (HWCVD) method: Role of substrate temperature. Sol. Energy Mater. Sol. Cells 2007, 91 (8), 714−720. (5) Soni, S.; Phatak, A.; Dusane, R. High deposition rate device quality a-Si: H films at low substrate temperature by HWCVD technique. Sol. Energy Mater. Sol. Cells 2010, 94 (9), 1512−1515. (6) Cocorullo, G.; Minarini, C.; Rubino, A.; Della Corte, F.; Rendina, I.; Terzini, E. Amorphous silicon waveguides and light modulators for integrated photonics realized by low-temperature plasma-enhanced chemical-vapor deposition. Opt. Lett. 1996, 21 (24), 2002−2004. (7) Scott, B. A.; Reimer, J. A.; Longeway, P. A. Growth and defect chemistry of amorphous hydrogenated silicon. J. Appl. Phys. 1983, 54 (12), 6853−6863. (8) Sazio, P. J.; Amezcua-Correa, A.; Finlayson, C. E.; Hayes, J. R.; Scheidemantel, T. J.; Baril, N. F.; Jackson, B. R.; Won, D.-J.; Zhang, F.; Margine, E. R. Microstructured optical fibers as high-pressure microfluidic reactors. Science 2006, 311 (5767), 1583−1586. (9) Baril, N. F.; Keshavarzi, B.; Sparks, J. R.; Krishnamurthi, M.; Temnykh, I.; Sazio, P. J.; Peacock, A. C.; Borhan, A.; Gopalan, V.; Badding, J. V. High-Pressure Chemical Deposition for Void-Free Filling of Extreme Aspect Ratio Templates. Adv. Mater. 2010, 22 (41), 4605−4611. (10) He, R.; Day, T. D.; Sparks, J. R.; Sullivan, N. F.; Badding, J. V. High Pressure Chemical Vapor Deposition of Hydrogenated Amorphous Silicon Films and Solar Cells. Adv. Mater. 2016, 28 (28), 5939−5942. (11) Ji, X.; Lei, S.; Yu, S.-Y.; Cheng, H. Y.; Liu, W.; Poilvert, N.; Xiong, Y.; Dabo, I.; Mohney, S. E.; Badding, J. V.; Gopalan, V. Single-



CONCLUSION We studied the kinetics of silane pyrolysis at pressures of 20− 33 MPa within the extreme aspect-ratio confinement of a 5.9 μm-diameter microcapillary reactor with elevated wall temperature. Silane flowing in the reactor decomposed at the microcapillary wall to deposit a thin film of hydrogenated amorphous silicon. Film growth rates at different temperatures were determined by measuring the thickness of the deposited films using scanning electron microscopy. The overall decomposition reaction was found to be first-order with respect to silane concentration, with activation energy and pre-exponential factor of 53.7 ± 2.9 kcal/mol and 1.5 × 1010 m/s, respectively. The resulting activation energy falls in the range of activation energies reported for hydrogen desorption reactions at c-Si surface, suggesting that hydrogen desorption from the surface is the rate-limiting step in film growth. To further investigate this observation, reactive molecular dynamics simulations of the thermal decomposition of silane on silicon surfaces were performed at 7 MPa pressure and 400 °C temperature using ReaxFF. No gas-phase reactions were observed during the ReaxFF simulations of high-pressure silane pyrolysis, in agreement with previously reported experimental results highlighting the effects of confinement. Furthermore, upon silane exposure of the simulated silicon surface, almost all of the available silicon surface sites (i.e., dangling bonds) were rapidly occupied with silicon-hydrides (mostly tri- and dihydrides), with further progress in film growth being 14999

DOI: 10.1021/acs.iecr.7b03515 Ind. Eng. Chem. Res. 2017, 56, 14995−15000

Article

Industrial & Engineering Chemistry Research

reactive molecular dynamics study, Part I. J. Phys. Chem. C 2012, 116 (30), 16111−16121. (33) Senftle, T. P.; Hong, S.; Islam, M. M.; Kylasa, S. B.; Zheng, Y.; Shin, Y. K.; Junkermeier, C.; Engel-Herbert, R.; Janik, M. J.; Aktulga, H. M. The ReaxFF reactive force-field: development, applications and future directions. npj Comput. Mater. 2016, 2, 15011. (34) Van Duin, A. C.; Dasgupta, S.; Lorant, F.; Goddard, W. A. ReaxFF: a reactive force field for hydrocarbons. J. Phys. Chem. A 2001, 105 (41), 9396−9409. (35) Van Duin, A. C.; Strachan, A.; Stewman, S.; Zhang, Q.; Xu, X.; Goddard, W. A. ReaxFFSiO reactive force field for silicon and silicon oxide systems. J. Phys. Chem. A 2003, 107 (19), 3803−3811. (36) Odden, J.; Egeberg, P.; Kjekshus, A. From monosilane to crystalline silicon, Part I: Decomposition of monosilane at 690−830K and initial pressures 0.1−6.6 MPa in a free-space reactor. Sol. Energy Mater. Sol. Cells 2005, 86 (2), 165−176. (37) Odden, J. O.; Egeberg, P. K.; Kjekshus, A. From monosilane to crystalline silicon. Part III. Characterization of amorphous, hydrogencontaining silicon products. J. Non-Cryst. Solids 2005, 351 (14), 1317− 1327. (38) Sazio, P. J. A.; Amezcua-Correa, A.; Finlayson, C. E.; Hayes, J. R.; Scheidemantel, T. J.; Baril, N. F.; Jackson, B. R.; Won, D.-J.; Zhang, F.; Margine, E. R.; Gopalan, V.; Crespi, V. H.; Badding, J. V. Microstructured Optical Fibers as High-Pressure Microfluidic Reactors. Science 2006, 311 (5767), 1583. (39) Baril, N. F.; He, R.; Day, T. D.; Sparks, J. R.; Keshavarzi, B.; Krishnamurthi, M.; Borhan, A.; Gopalan, V.; Peacock, A. C.; Healy, N.; Sazio, P. J. A.; Badding, J. V. Confined High-Pressure Chemical Deposition of Hydrogenated Amorphous Silicon. J. Am. Chem. Soc. 2012, 134 (1), 19−22.

Crystal Silicon Optical Fiber by Direct Laser Crystallization. ACS Photonics 2017, 4 (1), 85−92. (12) Giunta, C. J.; McCurdy, R. J.; Chapple-Sokol, J. D.; Gordon, R. G. Gas-phase kinetics in the atmospheric pressure chemical vapor deposition of silicon from silane and disilane. J. Appl. Phys. 1990, 67 (2), 1062−1075. (13) Purnell, J.-H.; Walsh, R. In The pyrolysis of monosilane. Proc. R. Soc. London, Ser. A 1966, 293, 543−561. (14) Coltrin, M. E.; Kee, R. J.; Evans, G. H. A Mathematical Model of the Fluid Mechanics and Gas-Phase Chemistry in a Rotating Disk Chemical Vapor Deposition Reactor. J. Electrochem. Soc. 1989, 136 (3), 819−829. (15) Ho, P.; Coltrin, M. E.; Breiland, W. G. Laser-induced fluorescence measurements and kinetic analysis of Si atom formation in a rotating disk chemical vapor deposition reactor. J. Phys. Chem. 1994, 98 (40), 10138−10147. (16) Odden, J.; Egeberg, P.; Kjekshus, A. From monosilane to crystalline silicon, Part II: Kinetic considerations on thermal decomposition of pressurized monosilane. Int. J. Chem. Kinet. 2006, 38 (5), 309−321. (17) Girshick, S.; Swihart, M.; Suh, S. M.; Mahajan, M.; Nijhawan, S. Numerical Modeling of Gas-Phase Nucleation and Particle Growth during Chemical Vapor Deposition of Silicon. J. Electrochem. Soc. 2000, 147 (6), 2303−2311. (18) Liu, S.-S.; Xiao, W.-D. Numerical simulations of particle growth in a silicon-CVD fluidized bed reactor via a CFD−PBM coupled model. Chem. Eng. Sci. 2014, 111, 112−125. (19) Slakman, B. L.; Simka, H.; Reddy, H.; West, R. H. Extending Reaction Mechanism Generator to Silicon Hydride Chemistry. Ind. Eng. Chem. Res. 2016, 55 (49), 12507−12515. (20) Robertson, R.; Gallagher, A. Reaction mechanism and kinetics of silane pyrolysis on a hydrogenated amorphous silicon surface. J. Chem. Phys. 1986, 85 (6), 3623−3630. (21) Gates, S.; Greenlief, C.; Kulkarni, S.; Sawin, H. Surface reactions in Si chemical vapor deposition from silane. J. Vac. Sci. Technol., A 1990, 8 (3), 2965−2969. (22) Gates, S.; Kulkarni, S. Kinetics of surface reactions in very lowpressure chemical vapor deposition of Si from SiH4. Appl. Phys. Lett. 1991, 58 (25), 2963−2965. (23) Claassen, W.; Bloem, J.; Valkenburg, W.; Van Den Brekel, C. The deposition of silicon from silane in a low-pressure hot-wall system. J. Cryst. Growth 1982, 57 (2), 259−266. (24) Oura, K.; Lifshits, V. G.; Saranin, A. A.; Zotov, A. V.; Katayama, M. Hydrogen Interaction with Clear and Modified Silicon Surfaces. Surf. Sci. Rep. 1999, 35 (1), 1−69. (25) Jasinski, J. M.; Gates, S. M. Silicon Chemical Vapor Deposition One Step at a Time. Acc. Chem. Res. 1991, 24 (1), 9−15. (26) Hsin, W. C.; Tsai, D. S.; Shimogaki, Y. Surface reaction probabilities of silicon hydride radicals in SiH4/H2 thermal chemical vapor deposition. Ind. Eng. Chem. Res. 2002, 41 (9), 2129−2135. (27) Buss, R. J.; Ho, P.; Breiland, W. G.; Coltrin, M. E. Reactive sticking coefficients for silane and disilane on polycrystalline silicon. J. Appl. Phys. 1988, 63 (8), 2808−2819. (28) Weerts, W.; de Croon, M.; Marin, G. The adsorption of silane, disilane and trisilane on polycrystalline silicon: a transient kinetic study. Surf. Sci. 1996, 367 (3), 321−339. (29) Greenlief, C. M.; Gates, S. M.; Holbert, P. A. Reaction kinetics of surface silicon hydrides. J. Vac. Sci. Technol., A 1989, 7 (3), 1845− 1849. (30) Sinniah, K.; Sherman, M. G.; Lewis, L. B.; Weinberg, W. H.; Yates, J. T., Jr; Janda, K. C. New mechanism for hydrogen desorption from covalent surfaces: The monohydride phase on Si (100). Phys. Rev. Lett. 1989, 62 (5), 567. (31) Buehler, M. J.; van Duin, A. C.; Goddard, W. A., III Multiparadigm modeling of dynamical crack propagation in silicon using a reactive force field. Phys. Rev. Lett. 2006, 96 (9), 095505. (32) Newsome, D. A.; Sengupta, D.; Foroutan, H.; Russo, M. F.; van Duin, A. C. Oxidation of silicon carbide by O2 and H2O: a ReaxFF 15000

DOI: 10.1021/acs.iecr.7b03515 Ind. Eng. Chem. Res. 2017, 56, 14995−15000