MATERIALS TRANSACTIONS
Online ISSN : 1347-5320
Print ISSN : 1345-9678
ISSN-L : 1345-9678
Materials Chemistry
Silicon Refining by Solidification from Liquid Si–Zn Alloy and Floating Zone Method
Yuanjia MaKouji YasudaAkifumi IdoTakeyuki ShimaoMing ZhongRika HagiwaraToshiyuki Nohira
Author information
JOURNAL FREE ACCESS FULL-TEXT HTML

2021 Volume 62 Issue 3 Pages 403-411

Details
Abstract

This study evaluated the refining ability of a proposed production process for solar-grade silicon utilizing the electrolytic reduction of SiO2 on a liquid zinc electrode in molten salt. The distribution behaviors of impurity elements during the precipitation of solid silicon from a liquid Si–Zn alloy were studied by thermodynamic calculations at 923 K. In the precipitation experiment, silicon granules were recovered from a liquid Si–Zn alloy, which was prepared from metallurgical-grade silicon. The impurity removal ratios exceeded 99% for C, Al, and Ca, and 90% for Fe. High removal ratios were attained for B and O as well. As the post-processing, a silicon ingot was produced from the precipitated silicon granules by the floating zone method. The Zn residue in the precipitated silicon was completely evaporated during the floating zone refining. The total content of metallic elements (Al, Ca, Fe, Ti, and Zn) was lower than 0.2 ppmw, even though metallurgical-grade silicon was used as the starting material.

1. Introduction

Photovoltaic (PV) power generation has attracted considerable attention as a clean and inexhaustible energy, which is expected to be a key technology to mitigate energy and environmental issues. Although PV power only accounted for 1.7% of the worldwide energy supply in 2017, according to World Energy Outlook 2018 by the International Energy Agency (IEA),1) it is expected to increase to 10% by 2040 if the cost of PV cell production is reduced.2) The global production of PV cells achieved 102.4 GW in 2018, approximately 350 times the value in 2000.3) Among the many types of solar cells, crystalline silicon solar cells accounted for 96.9% of the worldwide production in 2018.4) The global production of high-purity crystalline silicon also increased to 448,000 tons in 2018, by a factor of approximately 20 compared to that in 2000.5) Consequently, crystalline silicon solar cells are expected to remain mainstream in the PV industry in the long term.

High-purity silicon used in crystalline silicon solar cells is known as solar-grade silicon (SOG-Si), which requires a purity of 6–7 N. Currently, SOG-Si is produced by the following steps: (1) production of metallurgical-grade silicon (MG-Si) from SiO2 by carbothermal reduction, (2) production and purification of trichlorosilane, and (3) production of polycrystalline silicon from trichlorosilane gas by chemical vapor deposition (CVD).68) Steps (2) and (3) are referred to as the Siemens process. Although the purity of silicon produced by the Siemens process is very high, low productivity and high energy consumption hinder cost reduction. To develop a next-generation production process for SOG-Si, purification of MG-Si913) and metallothermic reduction of silicon halides by metal reductants1417) have been investigated.

Over the past two decades, we have been studying the electrolytic reduction of SiO2 to silicon in molten CaCl2 to develop a new production process for SOG-Si.1822) In molten CaCl2, electrochemical reduction of insulating SiO2 is realized by using a SiO2 contacting electrode, which provides the three-phase zone of conductor/SiO2/CaCl2. The produced silicon, which has high electrical conductivity at high temperature, becomes an additional electron pathway.   

\begin{equation} \text{SiO$_{2}$($s$)} + \text{4 e$^{-}$}\to \text{Si($s$)} + \text{2 O$^{2-}$} \end{equation} (1)
Several other research groups have also investigated the direct electrolytic reduction of SiO2.2333)

With the aim of efficient recovery of reduced silicon, we also proposed an electrolytic reduction process of SiO2 using a liquid zinc cathode in molten CaCl2.3436) Figure 1 schematically illustrates the proposed process.35,36) The overall process consists of three major steps: electrolysis, precipitation, and refining. In the electrolysis step, solid SiO2 is reduced to form a Si–Zn liquid alloy at the cathode. The separation of reduced silicon from the molten salt is simple because the reduced silicon forms a liquid Si–Zn alloy at the cathode.   

\begin{align} & \text{SiO$_{2}$($s$)} + \text{Si–Zn($l$, low Si conc.)} + \text{4 e$^{-}$}\\ &\quad \to \text{Si–Zn($l$, high Si conc.)} + \text{2 O$^{2-}$} \end{align} (2)
In the precipitation step, solid silicon is precipitated by decreasing the temperature of the liquid Si–Zn alloy.   
\begin{equation} \text{Si–Zn($l$, high Si conc.)}\to \text{Si($s$)} + \text{Si–Zn($l$, low Si conc.)} \end{equation} (3)
The liquid Si–Zn alloy with a low concentration of Si is reused as the cathode in the electrolysis step. The recovered silicon is then subjected to a refining step that comprises vacuum refining to remove residual Zn and directional solidification to manufacture SOG-Si ingots.

Fig. 1

A schematic drawing of the SOG-Si production process using electrochemical reduction of SiO2 powder on a liquid Si–Zn alloy cathode in molten CaCl2.35,36)

Silicon refining by solidification has been widely investigated not only for liquid silicon, but also for liquid Si–Al alloys,3740) Si–Cu alloys,41) Si–Sn alloys,42) and Si–Fe alloys.43,44) Although the refining of silicon using a liquid Si–Zn alloy is also considered to be possible, there are no reports so far, likely due to the experimental difficulty caused by the high vapor pressure of Zn. In our previous study,36) silicon granules were obtained by the electrolysis of SiO2 particles on a zinc cathode in molten CaCl2, and the concentrations of impurity elements (B, P, Al, Ca, Fe, Ti, and Zn) in the silicon granules were analyzed. However, the impurity behaviors during the precipitation and refining steps have not been quantitatively investigated. Through the previous experiment, we realized that the evaporation of zinc can be largely suppressed, even at 1123 K, the electrolysis temperature of SiO2, by submerging the Si–Zn alloy in molten CaCl2.35) The lower evaporation rate is explained by the low solubility of zinc in the molten salt. This method enables the distribution analysis of impurities during the refining of silicon by solidification from a liquid Si–Zn alloy.

In the present study, the distribution coefficients of impurity elements between solid Si and liquid Zn were calculated from thermodynamic data at the precipitation temperature of 923 K. Then, the behaviors of impurity elements in the precipitation step were experimentally evaluated using the liquid Si–Zn alloy, which was prepared from low-purity MG-Si. Finally, the precipitated silicon was melted and refined using the floating zone (FZ) method to analyze the behavior of impurities in the refining step.

2. Experimental

2.1 Silicon refining by solidification from liquid Si–Zn alloy

Figure 2 shows a schematic drawing of the experimental apparatus for the precipitation of solid silicon from a liquid Si–Zn alloy. 885 g of metallic zinc (FUJIFILM Wako Pure Chemical Corp., reagent grade, granule) and 197 g of CaCl2 (Kojundo Chemical Laboratory Co., Ltd., >99%) were charged in a graphite crucible (Toyo Tanso Co., Ltd., IG-110 grade, o.d. 80 mm × i.d. 70 mm × height 200 mm) and dried under vacuum at 773 K overnight. The evaporation of metallic zinc was suppressed by submerging the Si–Zn alloy in molten CaCl2. After the temperature was increased to 1123 K in a dry argon atmosphere, 20.0 g of MG-Si, which is equivalent to 5.0 at% of the prepared Si–Zn alloy, was charged through a carbon tube (Experiment 1) or a glass funnel (Experiment 2). The carbon tube had several holes cut at the face of the tube side to maintain both bath level and gas pressure equalization between the inside and outside of the tube. After the temperature was maintained at 1123 K for more than 24 h to sufficiently react with the silicon and zinc, the liquid Si–Zn alloy was cooled. Figure 3 shows the temperature changes during the cooling steps. In Experiment 1, the temperature was decreased from 1123 K to 723 K for 20 h, then maintained at 723 K for 5 h, and further reduced to 298 K for 20 h. In Experiment 2, the temperature was decreased from 1123 K to 773 K for 35 h, then maintained at 773 K for 10 h, and further decreased to 298 K for 5 h. A lump of zinc metal was recovered from the crucible after dissolution of CaCl2 with distilled water. The zinc metal lump was dissolved in a hydrochloric acid solution (20 mass%) to recover the precipitated silicon granules. The recovered silicon granules were then immersed in a hydrochloric acid solution (10 mass%) overnight.

Fig. 2

A schematic drawing of the experimental apparatus for the precipitation of solid silicon from a liquid Si–Zn alloy. (a) Thermocouple, (b) graphite crucible, (c) carbon tube (Exp. (1)), (d) funnel (Exp. (2)), (e) molten CaCl2, (f) liquid zinc, and (g) MG-Si.

Fig. 3

Temperature changes during the cooling step.

2.2 Silicon refinement and zinc evaporation by the FZ method

A portion of the silicon granules recovered in the precipitation experiment (Experiment 2) was crushed into powder in an agate mortar. Silicon powder with sizes below 0.1 mm was selected by a sieve and molded into a cylinder with a diameter of approximately 7 mm by isostatic pressing under a pressure of 60 MPa for 1 min. The molded silicon cylinder was hung with a platinum wire and inserted into a quartz tube of a floating zone furnace (Canon Machinery Inc., Desktop type A) as a feed silicon, as shown in Fig. 4. A 10N-purity poly-silicon rod (Furuuchi Chemical Corp., diameter 5 mm) was used as the seed silicon. The zone melting was conducted twice to evaporate the residual Zn to obtain a silicon ingot in an inverse argon stream with a uniform speed of 20–30 mm h−1 and a lamp voltage of 55–65 V. Silicon wafers were cut out from the prepared silicon ingot using a diamond wheel saw (South Bay Technology Inc., SBT650).

Fig. 4

A schematic drawing of the floating zone (FZ) method.

2.3 Impurity analysis of silicon samples

The analyses of impurity elements in the silicon samples were conducted using the high-frequency combustion–infrared absorption method (LECO, CS-LS600) for carbon analysis, inert gas fusion–infrared absorption method (LECO, TC600) for oxygen analysis, and glow discharge mass spectrometry (GD-MS; VG Elemental, VG 9000) for the other elements. Silicon samples were immersed in a hydrofluoric acid solution (20–25 mass%) for 20 min and then washed with distilled water and acetone before GD-MS analysis.

3. Thermodynamic Calculations

In order to predict the purification ability of the precipitation step, the distribution behaviors of the impurity elements were evaluated. When the temperature was decreased to 923 K from the electrolysis temperature of 1123 K, solid silicon corresponding to 5 at% was recovered. Although the liquid Zn equilibrated with solid Si at 923 K contains 1 at% Si, the distribution behavior of impurities was evaluated assuming that the liquid Zn contained no Si. That is, the distribution coefficients of impurities between solid Si and liquid Zn were calculated at infinite dilution.

When the standard state of impurity element A is defined as its pure stable state at 923 K, the activities of A in an equilibrium state are identical in both phases.   

\begin{equation} a_{\text{A(Zn(${l}$))}} = a_{\text{A(Si(${s}$))}} \end{equation} (4)
The activity is the product of activity coefficient and mole fraction as follows:   
\begin{equation} \gamma_{\text{A(Zn(${l}$))}}x_{\text{A(Zn(${l}$))}} = \gamma_{\text{A(Si(${s}$))}}x_{\text{A(Si(${s}$))}} \end{equation} (5)
Here, γA(Zn(l)) and γA(Si(s)) are the activity coefficients, and xA(Zn(l)) and xA(Si(s)) are the mole fractions of A in liquid Zn and solid Si, respectively. The activity coefficient, γA, is identical to the activity coefficient at infinite dilution, γ°A, when the concentration is low enough and the Henry’s law can be applied. Thus, the distribution coefficient of A at infinite dilution, k°A, is defined as the ratio of the mole fraction of solid Si phase to the liquid Zn phase as follows:   
\begin{equation} k^{\circ}{}_{\text{A}} = \frac{x_{\text{A(Si(${s}$))}}}{x_{\text{A(Zn(${l}$))}}} \end{equation} (6)
  
\begin{equation} = \frac{\gamma^{\circ}{}_{\text{A(Zn(${l}$))}}}{\gamma^{\circ}{}_{\text{A(Si(${s}$))}}} \end{equation} (7)
where, γ°A(Zn(l)) and γ°A(Si(s)) are the activity coefficients of A at infinite dilution in liquid Zn and solid Si, respectively. Although the equilibrium expressed in eq. (4) is attained for the A–Zn–Si ternary system, the value of γ°A(Zn(l)) and γ°A(Si(s)) can be calculated from A–Zn and A–Si binary systems with the assumption that interaction between impurity and Zn in solid Si and that between impurity and Si in liquid Zn can be ignored.

Typical impurity elements of B, C, O, P, Al, Ca, Fe, and Ti were selected for the present evaluation. Since the temperature dependences of γ°A have been already reported for B–Si,45) Al–Si,46) and Al–Zn47) binary systems, the values at 923 K were calculated based on these reports. For other systems, γ°A was evaluated from the phase diagrams and thermodynamic data. Table 1 summarizes the equilibrated phases with saturated solid Si or liquid Zn at 923 K, which are given by the binary phase diagrams for the A–Si and A–Zn systems.48,49) Furthermore, their standard Gibbs energy of formation, ΔG°f, at 923 K is also given in Table 1.5065)

Table 1 Equilibrated phases with saturated solid Si or liquid Zn at 923 K and their standard Gibbs energy of formation at 923 K.

For some cases, the equilibrium between solid Si and compound AySiz is considered, where y and z are the stoichiometry of the compound. In an equilibrium state, the activities of A in the Si phase, aA(Si(s, sat)), and in the AySiz phase, $a_{\text{A}(\text{A}_{y}\text{Si}_{z})}$, are identical. The same is applicable to the activities of silicon aSi(Si(s, sat)) and $a_{\text{Si}(\text{A}_{y}\text{Si}_{z})}$.   

\begin{equation} a_{\text{A(Si(${s}$, sat))}} = a_{\text{A(A${_{y}}$Si${_{z}}$)}} \end{equation} (8)
  
\begin{equation} a_{\text{Si(Si(${s}$, sat))}} = a_{\text{Si(A${_{y}}$Si${_{z}}$)}} \end{equation} (9)
The Gibbs energy of formation of AySiz, ΔG°f(AySiz), is given by the following equation:   
\begin{equation} \Delta G^{\circ}{}_{\text{f}}(\text{A$_{y}$Si$_{z}$}) = RT\ln\left(\frac{a_{\text{A(A${_{y}}$Si${_{z}}$)}}^{y}a_{\text{Si(A${_{y}}$Si${_{z}}$)}}^{z}}{a_{\text{A${_{y}}$Si${_{z}}$(A${_{y}}$Si${_{z}}$)}}}\right) \end{equation} (10)
Here, R is the gas constant (8.3145 J mol−1 K−1), T is the absolute temperature, and $a_{\text{A}_{y}\text{Si}_{z}(\text{A}_{y}\text{Si}_{\text{z}})}$ is the activity of AySiz in the AySiz compound. Here, the standard state of AySiz is its pure stable state at 923 K. If there is no solubility of AySiz in other phases, $a_{\text{A}_{y}\text{Si}_{z}(\text{A}_{y}\text{Si}_{\text{z}})}$ is regarded as unity.   
\begin{equation} a_{\text{A${_{y}}$Si${_{z}}$(A${_{y}}$Si${_{z}}$)}} = 1 \end{equation} (11)
The activity of Si in the solid Si phase is also regarded as unity because the amount of impurities in the Si phase is very small.   
\begin{equation} a_{\text{Si(Si(${s}$, sat))}} = 1 \end{equation} (12)
From these equations, the activity of A in the solid silicon phase is obtained from the following equation:   
\begin{equation} a_{\text{A(Si(${s}$, sat))}} = \exp\left(\frac{\Delta G^{\circ}{}_{\text{f}}(\text{A$_{y}$Si$_{z}$})}{yRT}\right) \end{equation} (13)
Then, γ°A(Si(s)) is calculated from the solubility of A in the solid Si phase, xA(Si(s, sat)), as follows:   
\begin{equation} \gamma^{\circ}{}_{\text{A(Si(${s}$))}} = \frac{a_{\text{A(Si(${s}$, sat))}}}{x_{\text{A(Si(${s}$, sat))}}} \end{equation} (14)
  
\begin{equation} = \cfrac{\exp\biggl(\cfrac{\Delta G^{\circ}{}_{\text{f}}(\text{A$_{y}$Si$_{z}$})}{yRT}\biggr)}{x_{\text{A(Si(${s}$, sat))}}} \end{equation} (15)

In a similar manner, γ°A(Zn(l)) is obtained by replacing Si(s) with Zn(l) in the above discussion. Here, the activity of Zn in the Zn(l) phase, aZn(Zn(l, sat)), needs further consideration. As shown in the phase diagrams, the solubility of impurities in the liquid Zn phase can reach several atomic percent at 923 K. Therefore, the value of aZn(Zn(l, sat)) cannot be regarded as unity where the standard state is pure liquid Zn, but rather regarded as the mole fraction of Zn, xZn(Zn(l, sat)), assuming Raoult’s law.   

\begin{equation} a_{\text{Zn(Zn(${l}$, sat))}} = x_{\text{Zn(Zn(${l}$, sat))}} \end{equation} (16)
  
\begin{equation} a_{\text{Zn(A${_{y}}$Zn${_{z}}$)}} = a_{\text{Zn(Zn(${l}$, sat))}} = 1 - x_{\text{A(Zn(${l}$, sat))}} \end{equation} (17)
Consequently, the value of γ°A(Zn(l)) is obtained as follows:   
\begin{equation} a_{\text{A(Zn(${l}$, sat))}} = (1 - x_{\text{A(Zn(${l}$, sat))}})^{-\frac{z}{y}}\exp\left(\frac{\Delta G^{\circ}{}_{\text{f}}(\text{A$_{y}$Zn$_{z}$})}{yRT}\right) \end{equation} (18)
  
\begin{equation} \gamma^{\circ}{}_{\text{A(Zn(${l}$))}} = \cfrac{(1 - x_{\text{A(Zn(${l}$, sat))}})^{-\frac{z}{y}}\exp\biggl(\cfrac{\Delta G^{\circ}{}_{\text{f}}(\text{A$_{y}$Zn$_{z}$})}{yRT}\biggr)}{x_{\text{A(Zn(${l}$, sat))}}} \end{equation} (19)

Table 2 summarizes the solubilities of the impurity elements in solid Si and liquid Zn phases at 923 K.46,48,51,60,6676) The solubilities of P(Zn(l)), Al(Zn(l)), Ca(Zn(l)), Fe(Zn(l)), and Ti(Zn(l)) were read from the phase diagrams.48,76) Most of the other solubility values were calculated from the reported functions of solubility with regards to temperature.46,48,51,60,6676) As for B(Zn(l)), the solubility was calculated from the interaction parameter of liquid B–Zn under the assumption of regular solution model, ΩB,Zn, 2.98 × 104 J mol−1 at 923 K, reported by Chen et al.68) In this calculation, the standard state of B is not a pure solid B, which is the stable phase at 923 K, but pure liquid B. When solid B and B-saturated liquid B–Zn are in an equilibrium state, the activity of B in liquid B–Zn with liquid B as the standard state, aB(l, Zn(l, sat)), is expressed by the Gibbs energy of fusion of B, ΔGfus(B).   

\begin{equation} \ln\left(\frac{a_{\text{B(${l}$, Zn(${l}$, sat))}}}{a_{\text{B(${l}$)}}}\right) = -\frac{\Delta G_{\text{fus}}(\text{B})}{RT} \end{equation} (20)
Since its standard state is pure liquid (aB(l) = 1),   
\begin{equation} \ln a_{\text{B(${l}$, Zn(${l}$, sat))}} = -\frac{\Delta G_{\text{fus}}(\text{B})}{RT} \end{equation} (21)
is given. Then, aB(l, Zn(l, sat)) is calculated to be 2.29 × 10−2 at 923 K from the reported ΔGfus(B) value of 2.90 × 104 J mol−1.51,60,61) In the equilibrium state, the following equation holds.   
\begin{align} RT\ln a_{\text{B(${l}$, Zn(${l}$, sat))}} &= \varOmega_{\text{B,Zn}}(1 - x_{\text{B(Zn(${l}$, sat))}})^{2} \\ &\quad + RT\ln x_{\text{B(Zn(${l}$, sat))}} \end{align} (22)
By substituting aB(l, Zn(l, sat)) = 2.29 × 10−2 and ΩB,Zn = 2.98 × 104, the value of xB(Zn(l, sat)) was calculated to be 4.7 × 10−4.

Table 2 Solubility of each impurity element in solid Si and liquid Zn at 923 K.

The calculated activity coefficients of impurity elements at infinite dilution in solid Si and liquid Zn are summarized in Table 3. Then, the distribution coefficients can be calculated using eq. (4), which are also listed in Table 3. The distribution coefficient of O (1.2 × 104) is substantially large, indicating that almost all the O impurity in the liquid Si–Zn alloy transfers to solid Si. On the other hand, the distribution coefficients of C, Al, Ca, Fe, and Ti are 5.1 × 10−3, 7.9 × 10−5, 2.3 × 10−3, 3.6 × 10−8, and 2.4 × 10−9, respectively. Thus, these impurities are considered to remain in the liquid Si–Zn alloy during the precipitation step. The removal of B and P is known to be difficult in simple solidification refining from liquid Si because of the large distribution coefficients, that is, 0.8 for B and 0.35 for P.77) In the present calculation, the distribution coefficients are 0.15 for B and 0.26 for P. The obtained smaller values compared to those for the Si(s)/Si(l) system demonstrate that the effective removal of B and P is possible by utilizing the liquid Si–Zn alloy system. All of these calculations indicate that the electrolytic process using a liquid Si–Zn alloy as a cathode has the advantages of high refining ability and high resistance against impurity inclusion. Even when carbon is deposited onto the cathode by the reduction of CO32− ions produced at a carbon anode, it is not contained in the solid Si phase in the precipitation step. The applicability of graphite anodes in the electrolysis step is a significant advantage for practical applications.

Table 3 Activity coefficient of each impurity element at infinite dilution in solid Si and liquid Zn with the standard states of their stable phase and distribution coefficient between solid Si and liquid Zn at 923 K.

For reference, the distribution coefficients of B for the liquid Si–Al alloy were reported to be 0.49 at 1473 K and 0.22 at 1273 K,45) and those of P were 0.12 at 1373 K and 0.061 at 1173 K.78) The distribution coefficients of both elements become smaller at lower temperatures. Although 923 K is even lower than the cases of the liquid Si–Al alloy, the distribution coefficients of the liquid Si–Zn alloy, 0.15 for B, and 0.26 for P, are higher than those for the Si–Al alloy. This is explained by the fact that the affinity of Zn to these impurities is lower than that of Al.

4. Results and Discussion

4.1 Silicon refining by solidification from liquid Si–Zn alloy and the FZ method

Figure 5(a) shows an optical image of the top view of the zinc ingot obtained from the precipitation experiment (1). The top portion consists of a metallic colored main body and gray–blue parts. Although the bottom part cannot be seen from the figure, the entire bottom was a metallic color. Figure 5(b) shows a cross-sectional SEM image of the zinc ingot. The top portion contains many 100–1000-µm particles. From EDX analysis, only Si was detected for the deposited particles, while only Zn was detected in the main body. Most deposited particles were observed at the top of the ingot, and a small amount of silicon particles were also found at the edge of the middle part of the lump. Considering the densities of liquid Zn, ρZn(l) = 6.56 g cm−3,79) and solid Si, ρSi(s) = 2.32 g cm−3,80) at 723 K, silicon particles floated at the top of the liquid zinc during the precipitation step. Figure 5(c) shows an optical image of the silicon granules recovered after acid leaching of the zinc lump. The granules have a size of approximately 1 mm and exhibit a metallic luster.

Fig. 5

(a) An optical image of a top view and (b) a SEM image of a cross section of the zinc ingot obtained after the precipitation experiment (1) in molten CaCl2, and (c) an optical image of the silicon granules recovered after acid leaching of the zinc ingot.

Figure 6 shows a cross-sectional SEM image of the zinc ingot obtained from the precipitation experiment (2). Similar to the zinc ingot obtained from experiment (1), many Si particles were deposited in the upper part. On the other hand, the image shows much larger size particles of 300–2000 µm and more inclusion of zinc into the gaps of silicon granules, compared with experiment (1). This result is attributed to the progressive crystal growth of silicon particles during precipitation at a lower cooling rate (experiment (1): 20 K h−1, (2): 10 K h−1).

Fig. 6

A SEM image of a cross section of the zinc ingot obtained after the precipitation experiment (2) in molten CaCl2.

Figure 7(a) shows an optical image of the silicon powder obtained by crushing the silicon granules prepared in the precipitation experiment (2). From this silicon powder, a silicon ingot was obtained by FZ melting, as shown in Fig. 7(b). The ingot has a gray–blue metallic luster and a diameter of 3–5 mm and length of 25 mm. The non-uniform patterns observed on the surface suggest that the silicon ingot is polycrystalline. In general, the crystallinity can be improved by slowing the moving speed during FZ melting. The ingot was cut with a diamond wheel saw at the broken lines indicated in Fig. 7(b). The sliced wafers were used to analyze the impurity contents. Specifically, parts A, B, and C were used for GD-MS analysis, carbon analysis, and oxygen analysis, respectively.

Fig. 7

Optical images of (a) silicon powder for silicon ingot and (b) silicon ingot obtained by FZ melting.

4.2 Impurity analysis of silicon samples

Table 4 lists the impurity contents in the silicon granules obtained in the precipitation experiments (1) and (2) and the silicon wafers obtained by FZ melting. The values for the original MG-Si are also listed for reference. The removal fractions in the precipitation experiments, rA, which are calculated from the contents in the original MG-Si, xA(MG-Si), and the precipitated silicon, xA(prec Si), are shown in parentheses.   

\begin{equation} r_{\text{A}} = \left(1-\frac{x_{\text{A(prec Si)}}}{x_{\text{A(MG-Si)}}}\right)\times 100 \end{equation} (23)
For the precipitated silicon granules, the removal fractions of Al, Ca, and Ti exceeded 99%, and the fractions of B, C, and Fe exceeded 90%. The high removal fractions agree with the low distribution coefficients calculated in Section 3. It is noteworthy that the removal fraction of B was above 95%. The removal of B is known to be difficult in the solidification refining from liquid silicon because the distribution coefficient between Si(s)/Si(l) is as high as 0.8.77) The removal of P has also been achieved at fractions of 54.2% and 66.7%. However, the removal fractions are smaller than expected, which might be due to the formation of a ZnSiP2 compound with a melting point of 1643 K.81,82) Concerning the removal fraction of O, high values of 80.8% and 90.0% were obtained, which is inconsistent with the calculated distribution coefficient of 1.2 × 104. As the concentration of O in MG-Si is obviously higher than the solubility limit in the solid Si phase,83) the presence of some oxides, except for SiO2, such as CaO, in MG-Si is indicated. Thus, the dissolution of oxides into molten CaCl2 contributed to the removal of O.

Table 4 Impurity contents in silicon samples with removal fractions and acceptable impurity levels for SOG-Si.

Comparing the impurity contents in the silicon granules obtained from the precipitation experiments (1) and (2), larger amount of Zn was detected in experiment (2). This is explained by the larger amount of Zn left in the gaps of the granules, as observed in the cross-sectional SEM image. For other impurity elements, no significant difference was found. Summarizing the above, high processing rate and small inclusion of zinc can be achieved at a high cooling rate. On the other hand, large silicon granules are obtained at a low cooling rate, which will be advantageous for silicon recovery from liquid Si–Zn alloy by filtering.

For the silicon wafer obtained by FZ melting, metal elements (Al, Ca, Fe, Ti, and Zn) were further removed from the precipitated silicon granules; the total concentration was less than 0.2 ppmw. It should be mentioned that Zn was effectively removed to less than 0.01 ppmw owing to its high vapor pressure,84) even though a relatively large amount of 5900 ppmw remained in the silicon granules. In this respect, this contrasts with the previous solidification refining from Si–M (M = Al, Cu, Fe, Sn, etc.) alloys. The residue of the alloying element in silicon is known as one of the major challenges of solidification refining. For the solidification refining from the Si–Al alloy, approximately 400–600 ppmw of Al remained in the precipitated silicon.85) As the distribution coefficient of Al in Si is 2 × 10−3,86) at least twice the directional solidification refining is required for Al removal. Thus, the extremely easy removal of Zn is considered to be a great advantage as an alloying element. In the silicon wafer, P and O were also further removed to 13 ppmw and 10 ppmw, respectively. In contrast, B and C increased compared with the precipitated silicon granules, perhaps because carbon is considered to be contaminated by a diamond wheel saw. The increase in B needs further study.

Compared with the acceptable level for SOG-Si shown in Table 4, the contents of Al and Ca in the silicon wafer meet the level in spite of their high contents in MG-Si. It is noteworthy that the concentration of O is lower than that in actual silicon wafers used for solar cells, that is, approximately 20 ppmw,83) which is produced by the CZ method. Although the contents of Fe, Ti, B, and P are higher than the acceptable level, this is mainly due to the high impurity contents in the starting material of MG-Si. When electrochemically reduced silicon, which has lower impurity contents than MG-Si, is used as the starting material, lower concentrations of Fe, Ti, B, and P are expected to be achieved.

5. Conclusion

The distribution behaviors of impurity elements between solid Si and liquid Zn were studied by thermodynamic calculations and experiments. The distribution coefficients of impurity elements were calculated from the reported data at 923 K. The considerably low calculated coefficients suggested that almost all of the C, Al, Ca, and Fe would remain in the liquid Si–Zn alloy during the silicon precipitation. As the distribution coefficients of B and P were calculated to be 0.15 and 0.26, respectively, their efficient removals were expected. Precipitation experiments were successfully performed by submerging the liquid Si–Zn alloy in molten CaCl2. From the impurity analysis of the obtained silicon granules, high removal fractions were achieved for C, Al, Ca, and Fe, which were consistent with the calculated distribution coefficients. The high removal fraction of B was also confirmed when compared with the solidification from liquid silicon. Contrary to the prediction from the calculation, O was removed from the starting material of MG-Si, which was explained by the dissolution of residual oxides such as CaO into molten CaCl2. The obtained silicon granules were transformed into wafers using the FZ method. For the silicon wafers, the total content of metal impurities (Al, Ca, Fe, Ti, and Zn) was less than 0.2 ppmw, wherein Al and Ca met the acceptable levels for SOG-Si. In the silicon wafers, the concentrations of P and O also decreased. In particular, the concentration of Zn decreased from 5900 ppmw in the granules to less than 0.01 ppmw in the wafers, which demonstrates the great advantage of Zn as an alloying element.

Acknowledgement

This study was partially supported by Core Research for Evolutionary Science and Technology (CREST) from the Japan Science and Technology Agency (JST); Grant-in-Aid for Scientific Research A, Grant Number 16H02410, from the Japan Society for the Promotion of Science (JSPS), the Japan Prize Foundation, the Kato Foundation for Promotion of Science, and the Joint Usage/Research Program on Zero-Emission Energy Research, Institute of Advanced Energy, Kyoto University (ZE28A12, ZE29A01).

REFERENCES
 
© 2021 The Mining and Materials Processing Institute of Japan
feedback
Top