International Journal of Photoenergy
Volume 2015 (2015), Article ID 513639, 10 pages
http://dx.doi.org/10.1155/2015/513639

Research Article

Control of the Gas Flow in an Industrial Directional Solidification Furnace for Production of High Purity Multicrystalline Silicon Ingots

Lijun Liu,¹ Xiaofang Qi,¹ Wencheng Ma,¹ Zaoyang Li,¹ and Yunfeng Zhang²

¹Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China
²Yingli Green Energy Holding Co., Ltd., Baoding, Hebei 071051, China

Received 30 April 2015; Revised 2 July 2015; Accepted 6 July 2015

Academic Editor: Elias Stathatos

Copyright © 2015 Lijun Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A crucible cover was designed as gas guidance to control the gas flow in an industrial directional solidification furnace for producing high purity multicrystalline silicon. Three cover designs were compared to investigate their effect on impurity transport in the furnace and contamination of the silicon melt. Global simulations of coupled oxygen (O) and carbon (C) transport were carried out to predict the SiO and CO gases in the furnace as well as the O and C distributions in the silicon melt. Cases with and without chemical reaction on the cover surfaces were investigated. It was found that the cover design has little effect on the O concentration in the silicon melt; however, it significantly influences CO gas transport in the furnace chamber and C contamination in the melt. For covers made of metal or with a coating on their surfaces, an optimal cover design can produce a silicon melt free of C contamination. Even for a graphite cover without a coating, the carbon concentration in the silicon melt can be reduced by one order of magnitude. The simulation results demonstrate a method to control the contamination of C impurities in an industrial directional solidification furnace by crucible cover design.

1. Introduction

Multicrystalline silicon (mc-Si) produced by the directional solidification (DS) method has dominated the photovoltaic market in recent years. However, the concentrations of impurities in mc-Si ingots are higher than those in single crystals. This disadvantage leads to mc-Si solar cells having lower efficiency than single crystal solar cells. To improve the efficiency of solar cells, the impurities in mc-Si must be controlled. The main impurities in mc-Si ingots are oxygen (O) and carbon (C). Boron-oxygen complexes acting as recombination centers of photocarriers cause the serious problem of light-induced degradation of solar cells [1]. The C precipitates in a silicon ingot influence the wafer cutting process and reduce the solar cell efficiency [2]. Therefore, an effective way to control the O and C concentrations is essential to produce high purity mc-Si ingots for high efficiency solar cells.

There are various ways to control the concentrations of impurities to produce high purity mc-Si ingots in a DS furnace, such as process optimization and hot-zone design of the furnace. The former method mainly focuses on modifying the growth rate [3], silicon melt convection [4, 5], argon gas flow intensity [5–7], and furnace pressure [6–8] rather than influencing the pathway of impurity transport. The latter method redesigns the local configuration of the DS furnace, such as designing gas shields [9–11] and gas injector valve [12], so that the gas flow and the pathway of impurity transport can be controlled. It provides a new concept for impurity control in a DS furnace for mc-Si ingot growth. Recently, Liu et al. [13] pointed out that the silicon melt is severely contaminated by C impurities from the furnace chamber at the final stage of the melting process. To alleviate C contamination of the silicon melt before the solidification process, it is essential to control the transport pathway of the impurities in the furnace chamber and prevent them from being transported to the silicon melt free surface. For this purpose, in this study a crucible cover was designed as gas flow guidance in the GTS-DSS450 industrial furnace, as shown in Figure 1, which is different from previous furnaces such as GTS-DSS240 and GTS-HEM240 [10]. The cover can regulate the pathway of the argon gas flow in the furnace chamber. Gao et al. [9] numerically investigated the cover’s effect on O and C transport in a small laboratory-scale DS furnace. However, for an industrial DS furnace such as GTS-DSS450, which can produce a 450 kg ingot in one run, the cover’s effect on impurity transport is still unclear because of the large size and highly complex thermal interactions in the furnace. Very little research has been reported regarding the effect of crucible cover design on impurity transport in such an industrial furnace.

Figure 1: Configuration of an industrial directional solidification furnace.

In this paper, the industrial DS furnace GTS-DSS450 was used to numerically investigate the effect of the cover design on the argon gas flow and coupled oxygen/carbon transport in the furnace chamber as well as the contamination of impurities in the silicon melt. All thermal interactions, such as melt convection, argon gas flow, solid conduction, thermal radiation, and phase change, were included in the global simulations. The main chemical reactions of O and C in the DS furnace were considered. This study provides a possible way to produce high purity mc-Si ingots through the optimal design of the crucible cover in an industrial DS furnace.

2. Model Description

2.1. Models of Global Heat Transfer and Coupled Oxygen/Carbon Transport

Global modeling, including melt convection, argon flow, solid thermal conduction, and thermal radiation, was first carried out for the DS furnace. The configuration of the DS furnace to produce mc-Si ingots is shown in Figure 1. The furnace is subdivided into a number of subblocks. The structured/unstructured combined mesh scheme is applied. The major assumptions in the model are as follows: (1) the quasi-steady-state assumption is applied, (2) the geometry of the furnace configuration is axisymmetric, (3) radiative surfaces are diffuse-gray, (4) the melt flow is incompressible and the Boussinesq assumption is applicable, and (5) the argon gas in the furnace chamber is ideal and completely transparent. The algorithms for the global modeling of heat transfer in a DS furnace have been described in our previous paper [14]. The low Mach approximation and the ideal gas law are applied to the argon gas. The governing equations for argon gas flow are as follows:

where is the argon gas density, is the gas velocity, is the gas pressure, is the dynamic viscosity, is the strain rate tensor, is the gravity acceleration vector, is the heat capacity, is the temperature, is the heat conductivity, and are the reference density and reference pressure, respectively, and is the specific gas constant, which is 208.1 J/(kgK) for argon gas. In the current simulation, the flow rate and inlet temperature of the argon gas are 20 L/min and 300 K, respectively. The furnace pressure is 0.6 bar. The temperature of the outer wall of the chamber is taken to be 300 K.

Based on the global thermal and flow fields, the coupled oxygen/carbon transport was simulated. A generally accepted pathway for O and C transport in a silicon solidification furnace is as follows [6, 7]. (1) O dissolves from the quartz crucible wall into the silicon melt through a Si₃N₄ coating. (2) The dissolved O atoms in the silicon melt combine with a silicon atom and evaporate into the argon gas flow in the form of SiO gas at the gas/melt interface. (3) SiO is then carried away by the argon flow to all of the graphite fixture’ surfaces and reacts with them to form CO gas. (4) The CO gas is transported back with the argon gas flow to the gas/melt interface and then dissolves into the silicon melt in the form of C and O atoms. (5) Finally, the C and O atoms in the silicon melt are segregated into the growing crystal. The chemical reaction equations corresponding to the above five reactions have been described in [7].

The effects of O and C species on the density of the silicon melt were neglected because of their low concentrations. Thus, the governing equations of O and C transport in the silicon melt are expressed as

where and are the mass fractions of O and C atoms in the silicon melt, respectively, is the density of melt, is the velocity of melt, and and are the diffusivities of O and C atoms in the melt, respectively. Their values are taken to be 5.0 × 10⁻⁸ m²/s [15].

The governing equations for the transport of SiO and CO species in the argon gas are as follows:

where and are the mass fractions of SiO and CO in argon gas, respectively, is the density of argon gas, and is the argon gas flow velocity. and are the diffusivities of SiO and CO in the argon gas, respectively, which are expressed as [2]

The boundary conditions for the impurity species transport in the silicon and argon gas regions are as follows:

(a) On the inner wall of the quartz crucible, the equilibrium concentration of O atoms can be expressed as [16]

For the C atoms at the inner wall of the crucible, a zero flux boundary condition is applied.

(b) At the melt-gas interface, O, C, SiO, and CO coexist. Equilibrium relationships for the same element between the silicon melt and argon gas are given as follows [2]:

where and are the molar concentrations of SiO and CO in the argon gas, respectively, and are the molar concentrations of O and C in the silicon melt, respectively, and is the universal gas constant, which is equal to 8.314 J/(molK).

Two more boundary conditions are needed to solve for the four unknown variables at the melt-gas interface according to the law of mass conservation:

(c) SiO reacts with carbon at the surfaces of the hot graphite fixtures to generate CO. The reaction is assumed to be reversible and the Gibbs free energy of the reaction is expressed as follows [2]:

The equilibrium constant of the reaction is . The coupled boundary conditions at the graphite fixture’ surfaces can be expressed as

(d) At the argon gas inlet, the concentrations of SiO and CO are set to zero. At other solid component surfaces, zero fluxes of SiO and CO are applied. At the argon gas outlet, zero gradients of SiO and CO are used. Zero radial gradients for all the species are applied along the centerline of the fluid regions in the furnace.

The relationship between the mass fraction and the molar concentration in the aforementioned equations and boundary conditions is , where is the atomic weight of the impurity species. For the current simulations, the fully molten state (i.e., only the silicon melt exists in the crucible and crystallization has not started) was chosen. The first reason is that it is convenient to analyze the effect of the crucible cover design on the impurity contamination of the silicon melt before the solidification process. Another reason is that there is always a stabilization period in a DS process after the silicon feedstock is melted and before solidification starts in which the thermal system can be treated as quasi-steady state.

2.2. Crucible Cover Designs

Three crucible cover designs were selected and their effects on the transport of oxygen and carbon impurities in an industrial DS furnace for growing mc-Si were compared. The first design is called the no-cover design. That is, there is no cover on the crucible in the GTS-DSS450 furnace. The second design is the GTS-DSS450 furnace with a flat cover, as shown in Figure 1. The third design is shown in Figure 2, in which the flat cover is replaced by a U-shaped cover (U cover). There are two parameters in the third design: the gap between the melt surface and the bottom surface of the cover and the gap between the side surface of the cover and the inside surface of the crucible . For the current designs under consideration, is 40 mm and is 20 mm.

Figure 2: Design of the crucible U cover.

To investigate the effect of the crucible cover design on the contamination of impurities in the silicon melt by the environment in the furnace chamber, the content of impurities in the silicon feedstock was not considered in the simulations. For covers made of molybdenum or graphite with a perfect coating, no chemical reaction occurs on the cover’s surface. The hot heaters are the unique C source in the simulation. If the cover is made of graphite without a coating, SiO will be consumed on the cover’s surface to produce CO and the cover will be another C source. Both of the above cases will be discussed in the following sections. For all of the simulations, the quartz crucible is the only O source.

3. Results and Discussion

3.1. Comparison of Impurity Transport Pathways

The argon flow pattern and impurity transport pathway in the furnace chamber for the different cover designs were first investigated. Figure 3 shows the stream lines of argon flow above the silicon melt free surface for the different designs. In Figure 3(a), a strong gas recirculation developed in the space above the silicon melt. The impurities evaporated from the melt free surface can be easily taken away, while CO produced at the hot heaters’ surfaces can be transported back at the same time. In Figure 3(b), a flat cover made of molybdenum or tungsten was used and gas recirculation was limited to a smaller area above the melt free surface than with the no-cover design. It is difficult for CO to be transported back to the gas/melt interface by convection and diffusion. However, gas recirculation under the flat cover may also prohibit impurities that come out of the melt surface from being taken away. In Figure 3(c), gas recirculation in the space between the cover and the melt free surface disappears, leaving only a small backflow eddy. The argon flow can easily take away impurities that evaporate from the melt free surface and prevent CO from being transported back.

(a)

(b)

(c)

Figure 3: Stream lines above the melt free surface for different cover designs: (a) no cover, (b) flat cover, and (c) U cover.

In the above discussion, molybdenum or other materials with high melting points and low chemical activity were used to make the cover an ideal insulator. That is, SiO does not react with the cover to produce CO. However, a graphite board with or without a coating is usually used in the actual production process considering the production cost. Even for such a graphite cover with a coating, the coating can be damaged through repeated use. Therefore, cases with chemical reaction occurring on the cover’s surface should also be investigated. For the flat-cover design, CO produced at the lower surface of the cover can be transported to the melt free surface by argon gas recirculation. For the U-cover design, there is no gas recirculation above the melt free surface to take back the CO, but the distance between the melt free surface and the cover is small, which makes it easy for CO to diffuse from the cover to the melt surface. Therefore, the effects of the different designs on the production of high purity mc-Si are not clear.

3.2. Oxygen and Carbon Distributions

The cover is first assumed to be made of molybdenum or other materials with high melting points and low chemical activity. Therefore, no SiO is consumed and no CO is produced on the cover’s surface. The effects of the different cover designs on the transport of SiO and CO in the furnace chamber are compared. Figure 4 shows the distributions of SiO and CO above the melt free surface for the different cover designs. Comparing the SiO distributions in Figures 4(a) and 4(b), the concentrations of SiO for the no-cover and flat-cover designs are almost the same above the melt free surface. The value is about 1.0 × 10⁻⁹ mol/cm³. That is, because the flat-cover design only slightly changes the argon flow pattern compared with the no-cover design, as shown in Figures 3(a) and 3(b), it has little effect on the evaporation and transport of SiO. However, for the U-cover design, the concentration of SiO ranges from 1.0 × 10⁻¹¹ to 1.0 × 10⁻⁹ mol/cm³ in the space between the cover and the melt free surface, as shown in Figure 4(c). The reason for this is that there is no gas recirculation above the melt free surface for such a hot-zone design and no SiO is carried back. The maximum concentration of SiO is located at the melt/gas/crucible junction point for all three designs. It is the configuration of the crucible that prevents the argon gas from reaching the junction point, and SiO evaporated here cannot be easily taken away.

(a)

(b)

(c)

Figure 4: Distributions of SiO (left) and CO (right) concentration above the melt free surface for different cover designs without surface reaction: (a) no cover, (b) flat cover, and (c) U cover.

Comparing the CO concentrations in Figure 4(a) with those in Figure 4(b), the cover design has a significant influence on the CO distributions above the melt free surface. For the no-cover and flat-cover designs, the CO concentrations are in the order of 10⁻¹⁰ and 10⁻¹¹ mol/cm³, respectively. The flat-cover design can reduce the CO concentration by one order of magnitude compared with the no-cover design by preventing CO produced at the heaters’ surfaces from being transported back to the melt free surface. In Figure 4(b), the largest gradient of CO is located in the space between the crucible top surface and the cover bottom surface, where the distance between the crucible and the cover may play an important role in controlling CO transport. In Figure 4(c), the concentration of CO at the melt free surface can be reduced to 10⁻²⁹ mol/cm³ for the U-cover design, which means that it is almost free of carbon. The reasons for this are that the silicon feedstock is assumed to be pure and the hot heaters’ surfaces are the only carbon source in the current simulation. In addition, the distance between the inner surface of the crucible and the outside surface of the cover is very small and no gas recirculation occurs in the area between the cover and the melt free surface, both of which make it difficult for CO produced at the heaters’ surfaces to be transported back by gas convection or diffusion.

After investigating the effect of the cover design on the distributions of SiO and CO in the DS furnace chamber, the concentrations of O and C atoms in the silicon melt were compared for the different designs. Figure 5 shows the distributions of O and C in the silicon melt for the different cover designs under the same conditions as in Figure 4. The O concentration in the silicon melt is almost the same for the three designs, and the values range from 10¹⁶ to 10¹⁷ atoms/cm³. The almost identical values of the O concentration for the no-cover and flat-cover designs correspond to the similar argon flow patterns and SiO distributions above the melt free surface, as shown in Figures 4(a) and 4(b). However, for the U-cover design, the O concentration in Figure 5(c) is close to those in Figures 5(a) and 5(b), while the patterns of argon flow and the SiO distribution above the melt free surface are quite different to those of the other two designs, as shown in Figures 3 and 4. This is because of the fact that it is the melt flow rather than the argon gas flow that dominates the transport of O atoms in the melt. As a result, the cover design has little effect on the O concentration in the silicon melt. Because the quartz crucible is the only O source, the maximum concentrations of O in the silicon melt for the three designs are almost the same. Compared with the effect of the cover design on the O concentration, the effect on the C concentration in the silicon melt is significant. The C concentrations shown in Figures 5(a) and 5(b) are 10¹⁹ and 10¹⁷ atoms/cm³, respectively. The flat-cover design can reduce the C concentration in the silicon melt by two orders of magnitude compared with the no-cover design. For the U-cover design, as shown in Figure 5(c), the silicon melt is almost free of C contamination. Therefore, we can conclude that the U-cover design is the best choice to reduce C contamination in the silicon melt.

(a)

(b)

(c)

Figure 5: Distributions of O (left) and C (right) concentration in the silicon melt for different cover designs without surface reaction: (a) no cover, (b) flat cover, and (c) U cover.

Surface reaction of the crucible cover was not taken into account in the above discussion. However, it is not always consistent with the actual production process, in which the crucible cover is usually made of graphite and SiO is consumed on the cover’s surface to produce CO. In the following discussion, chemical reaction on the cover’s surface is considered. Figure 6 shows the distributions of SiO and CO above the melt free surface for the different cover designs that take part in the chemical reaction. The SiO concentrations above the melt free surface in Figures 6(a) and 6(b) are quite different. In Figure 6(a), the SiO concentration is about 6.0 × 10⁻¹⁰ mol/cm³ for the flat cover. However, for the U-cover design, the SiO concentration ranges from 1.0 × 10⁻¹² to 1.0 × 10⁻⁹ mol/cm³ in the space between the cover and the melt free surface, as shown in Figure 6(b). This indicates that the U-cover design can still reduce the SiO concentration near the melt free surface compared with the flat-cover design, even though chemical reaction on the crucible cover’s surface exists. Comparing the SiO concentration in Figure 6(a) with that in Figure 4(b), the SiO concentration (about 6.0 × 10⁻¹⁰ mol/cm³) above the melt free surface in Figure 6(a) is slightly lower than that in Figure 4(b) (about 1.0 × 10⁻⁹ mol/cm³) and the distribution patterns are similar. The reason for this is that the crucible cover chemically reacts and consumes the SiO gas that evaporates from the melt free surface when the cover is made of carbon material. Therefore, the SiO concentration shown in Figure 6 is reduced.

(a)

(b)

Figure 6: Distributions of SiO (left) and CO (right) concentration above the melt free surface for different cover designs with surface reaction: (a) flat cover and (b) U cover.

For the CO distribution, the U-cover design can still reduce the CO concentration in the space between the melt free surface and the cover compared with the flat-cover design even with chemical reaction on the cover’s surface, as shown in Figure 6. However, the CO concentrations in both Figures 6(a) and 6(b) are higher than those in Figures 4(b) and 4(c), respectively, because the crucible cover is another carbon source in the furnace chamber.

Finally, the O and C atom distributions in the silicon melt are compared for the different cover designs with chemical reaction on the cover’s surface. The concentration of O is about 10¹⁶ to 10¹⁷ atoms/cm³ in Figures 7(a) and 7(b). This indicates that the cover design has little effect on the oxygen concentration in the silicon melt. However, for the C distribution, the U-cover design can reduce the concentration by one order of magnitude compared with the flat cover, as shown in Figure 7. Comparing the C concentration in Figure 7(a) with that in Figure 5(a), the graphite flat cover without a coating cannot reduce the C concentration in the silicon melt. However, for the graphite U cover without a coating, it can still reduce the C concentration by one order of magnitude compared with the no-cover design by comparison of Figure 7(b) with Figure 5(a).

(a)

(b)

Figure 7: Distributions of O (left) and C (right) concentration in the silicon melt for different cover designs with surface reaction: (a) flat cover and (b) U cover.

4. Conclusions

Global simulations were carried out including all types of heat transfer and coupled oxygen/carbon transport in an industrial DS furnace. Three types of cover design, namely, no-cover, flat-cover, and U-cover designs, were compared to investigate their effect on oxygen and carbon transport during the DS process. The U-cover design can significantly change the argon flow pattern and SiO distribution above the melt free surface compared with the other two designs, while it had little effect on the O concentration in the silicon melt. For the carbon impurity, the cover design can significantly influence the CO distribution above the melt free surface as well as the carbon concentration in the silicon melt. The flat-cover design without chemical reaction on its surface can reduce the C concentration in the silicon melt by two orders of magnitude compared with the no-cover design. For the U-cover design without chemical reaction on its surface, the silicon melt is almost free of C contamination. For the graphite covers without a coating, the flat-cover design has little influence on the C concentration and the U-cover design can reduce the C concentration in the silicon melt by one order of magnitude compared with the no-cover design. This study may provide a possible way to produce high purity mc-Si for solar cells.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 51176148 and 51406156) and Fundamental Research Funds for the Central Universities of China.

References

T. Saitoh, X. Wang, H. Hashigami et al., “Suppression of light degradation of carrier lifetimes in low-resistivity CZ-Si solar cells,” Solar Energy Materials and Solar Cells, vol. 65, no. 1, pp. 277–285, 2001.

D. E. Bornside, R. A. Brown, T. Fujiwara, H. Fujiwara, and T. Kubo, “Effects of gas-phase convection on carbon contamination of Czochralski-grown silicon,” Journal of the Electrochemical Society, vol. 142, no. 8, pp. 2790–2804, 1995.

M. Trempa, C. Reimann, J. Friedrich, and G. Müller, “The influence of growth rate on the formation and avoidance of C and N related precipitates during directional solidification of multi crystalline silicon,” Journal of Crystal Growth, vol. 312, no. 9, pp. 1517–1524, 2010.

C. Reimann, M. Trempa, T. Jung, J. Friedrich, and G. Müller, “Modeling of incorporation of O, N, C and formation of related precipitates during directional solidification of silicon under consideration of variable processing parameters,” Journal of Crystal Growth, vol. 312, no. 7, pp. 878–885, 2010.

B. Gao, S. Nakano, and K. Kakimoto, “Reduction of oxygen impurity in multicrystalline silicon production,” International Journal of Photoenergy, vol. 2013, Article ID 908786, 6 pages, 2013.

B. Gao, S. Nakano, and K. Kakimoto, “Global simulation of coupled carbon and oxygen transport in a unidirectional solidification furnace for solar cells,” Journal of the Electrochemical Society, vol. 157, no. 2, pp. H153–H159, 2010.

Z. Y. Li, L. J. Liu, W. C. Ma, and K. Kakimoto, “Effects of argon flow on impurities transport in a directional solidification furnace for silicon solar cells,” Journal of Crystal Growth, vol. 318, no. 1, pp. 304–312, 2011.

Y.-Y. Teng, J.-C. Chen, C.-W. Lu, H.-I. Chen, C. Hsu, and C.-Y. Chen, “Effects of the furnace pressure on oxygen and silicon oxide distributions during the growth of multicrystalline silicon ingots by the directional solidification process,” Journal of Crystal Growth, vol. 318, no. 1, pp. 224–229, 2011.

B. Gao, X. J. Chen, S. Nakano, and K. Kakimoto, “Crystal growth of high-purity multicrystalline silicon using a unidirectional solidification furnace for solar cells,” Journal of Crystal Growth, vol. 312, no. 9, pp. 1572–1576, 2010.

X. Ma, L. Zheng, H. Zhang, B. Zhao, C. Wang, and F. Xu, “Thermal system design and optimization of an industrial silicon directional solidification system,” Journal of Crystal Growth, vol. 318, no. 1, pp. 288–292, 2011.

Y.-Y. Teng, J.-C. Chen, C.-W. Lu, and C.-Y. Chen, “Numerical investigation of oxygen impurity distribution during multicrystalline silicon crystal growth using a gas flow guidance device,” Journal of Crystal Growth, vol. 360, no. 1, pp. 12–17, 2012.

M. P. Bellmann, D. Lindholm, and M. M'Hamdi, “A novel method for gas flow and impurity control in directional solidification of multi-crystalline silicon,” Journal of Crystal Growth, vol. 399, pp. 33–38, 2014.

X. Liu, B. Gao, and K. Kakimoto, “Numerical investigation of carbon contamination during the melting process of Czochralski silicon crystal growth,” Journal of Crystal Growth, vol. 417, pp. 58–64, 2014.

L. J. Liu and K. Kakimoto, “Partly three-dimensional global modeling of a silicon Czochralski furnace. I. Principles, formulation and implementation of the model,” International Journal of Heat and Mass Transfer, vol. 48, no. 21-22, pp. 4481–4491, 2005.

L. J. Liu, S. Nakano, and K. Kakimoto, “Investigation of oxygen distribution in electromagnetic CZ–Si melts with a transverse magnetic field using 3D global modeling,” Journal of Crystal Growth, vol. 299, no. 1, pp. 48–58, 2007.

H. Matsuo, R. B. Ganesh, S. Nakano et al., “Thermodynamical analysis of oxygen incorporation from a quartz crucible during solidification of multicrystalline silicon for solar cell,” Journal of Crystal Growth, vol. 310, no. 22, pp. 4666–4671, 2008.