Abstract
A theoretical procedure for the creation of a stable drop-like static meniscus, appropriate for the growth of a single crystal tube, with a priori specified inner and outer radius, is presented. The method locates the controllable part of the pressure difference across the free surface. It consists in a set of calculus, which leads to the determination of the melt column height (between the horizontal crucible melt level and the shaper top level) in function of the pressure of the gas flow (introduced in the furnace for release the heat) in order to obtain the desired meniscus. The procedure is presented in general and is numerically illustrated for InSb tubes. The novelty is the algorithm for the exact determination of , which has to be used, the determination of the melt column height, and the evaluation of the effect of shaper radii. The setting of the thermal conditions, which assure that for the obtained static meniscus the solidification conditions are satisfied at the βrightβ places, is not considered here.
1. Introduction
The conventional melt growth techniques, as Bridgman growth [1β3] or Czochralski pulling [4β6] of single crystals, typically produce ingots of circular or square cross-sections which need to be cut in hundreds of slices to produce wafers. Using these processes, it is difficult to produce thin wafers from an ingot without wasting 40%β50% of material as kerfs during the cutting process. For this reason the E.F.G. technology can be more appropriate to produce single crystals with prescribed shapes and sizes which can be used without additional machining.
The growth of silicon tubes by E.F.G. process was first reported by Erris et al. [7]. In [7] a theory of tube growth by E.F.G. process is developed to show the dependence of tube wall thickness on the growth variables. The theory uses approximation reported in [8, 9], and it has been shown to be a useful tool understanding the feasible limits of the wall thickness control. A more accurate predictive model would require an increase of the acceptable tolerance range introduced by approximation.
Later, the heat flow in a tube growth system was analyzed in [10β19].
The state of the arts at the time 1993-1994, concerning the calculation of the meniscus shape in general in the case of the growth by E.F.G. method is summarized in [20]. According to [20], for the general differential equation describing the free surface of a liquid meniscus, possessing axial symmetry, there are no complete analysis and solution. For the general equation only numerical integrations were carried out for a number of process parameter values that were of practical interest at the moment. The authors of [21, 22] consider automated crystal growth processes based on weight sensors and computers. They give an expression for the weight of the meniscus, contacted with crystal and shaper of arbitrary shape, in which there are two terms related to the hydrodynamic factor.
In [23] it is shown that the hydrodynamic factor is too small to be considered in the automated crystal growth. In [24] a theoretical and numerical study of meniscus dynamics, under symmetric and asymmetric configurations, is presented. A meniscus dynamics model is developed to consider meniscus shape and its dynamics, heat and mass transfer around the die top and meniscus. Analysis reveals the correlations among tube thickness, effective melt height, pull rate, die top temperature, and crystal environmental temperature.
In [25] the effect of the controllable part of the pressure difference on the free surface shape of the static meniscus is analyzed for the tube growth by E.F.G. method for materials for which ; ; .
The present paper concerns also the shape and the stability of the free surface of a static meniscus (pulling rate equal to zero). More precisely, it is shown in which kind the explicit formulas reported in [25] can be combined in order to create a stable static meniscus having a free surface with prescribed size and shape, which is appropriate for the growth of a single crystal tube having a priori specified inner and outer radii. The free surface of a static meniscus is appropriate for the growth of a single crystal tube of constant inner radius and constant outer radius if the angle between the tangent lines to the free surface at the points (Figure 1 where the solidification conditions have to be assured) and the vertical is equal to the growth angle . Moreover, the function describing the free surface has to minimize the energy functional of the melt column (i.e., the meniscus has to be stable). In this paper we give a procedure for the choice of the melt column height, between the horizontal crucible melt level and shaper top level and of the pressure of the gas flow introduced in the furnace (for release the heat), in order to create a static meniscus of which free surface is appropriate for the growth of a single crystal tube of constant inner radius and outer radius .The thermal problem concerning the setting of the thermal conditions, which assure that for the obtained static meniscus at the level the solidification conditions are satisfied is not considered in this paper. The novelty consists in the fact that the free surface is not approximated by an arc with constant curvature, the computation takes into account the pressure of the gas flow, and the stability of the free surface is assured.

2. The Free Surfaces Equations and the Pressure Difference Limits
For a single crystal tube growth by E.F.G. technique, in hydrostatic approximation, the outer free surface equation of the static meniscus is
and the inner free surface equation is Here, Ξ³ is the surface tension of the melt; Ο is the melt density; is the gravitational acceleration; are the coordinates with respect to the Oz axis, directed vertically upwards; is the radial coordinate with respect to the axis, oriented horizontal; are the outer and inner radius of the shaper, respectively; are the pressure difference across the outer and inner free surface, respectively:
In (2.3), denotes the hydrodynamic pressure in the meniscus melt due to the thermal and Marangoni convection; are the pressure of the gas flow introduced in the exterior and in the interior of the tube, respectively, for releasing the heat from the inner and outer side of the tube wall; denotes the melt column βheightβ between the horizontal crucible melt level, and the shaper top level (Figure 1). is positive when the crucible melt level is under the shaper top level and it is negative when the shaper top level is under the crucible melt level.
The solution of (2.1) has to satisfy the following conditions:
where is the tube outer radius; is the growth angle; is the contact angle between the outer free surface and the outer edge of the shaper top and (Figure 1).
Condition (2.4a) expresses that at the point (the left end of the outer free surface), where the solidification has to be realized, the angle between the tangent line to the free surface and the vertical is equal to the growth angle (i.e., the tangent to the tube outer wall is vertical).
Condition (2.4b) expresses that at the point (the right end of the outer free surface, where the free surface is attached to the shaper edge), the angle between the tangent line to the free surface and horizontal (i.e., the contact angle) is equal to .
Condition (2.4c) expresses that at the point the free surface is attached to the shaper edge.
Moreover, the solution has to minimize the energy functional of the melt column: The solution of (2.2) has to satisfy the following conditions: where: is the tube inner radius; is the growth angle; is the contact angle between the inner free surface and the inner edge of the shaper top and (Figure 1).
Condition (2.6a) expresses that at the point (the left end of the inner free surface, where the free surface is attached to the shaper edge), the angle between the tangent line to the free surface and horizontal (i.e., the contact angle) is equal to .
Condition (2.6b) expresses that at the point (the right end of the inner free surface), where the solidification has to be realized, the angle between the tangent line to the free surface and the vertical is equal to the growth angle (i.e., the tangent to the tube inner wall is vertical).
Condition (2.6c) expresses that at the point the free surface is attached to the shaper edge.
Moreover, the solution has to minimize the energy functional of the melt column:
Based on the mathematical theorems, rigorously proven in [25], the following statements, regarding the creation of an appropriate meniscus, can be formulated.
Statement 1
If the solution of the initial value problem (IVP)
is convex, then it does not represent the outer free surface of an appropriate drop-like static meniscus.
Comment 1. The above statement shows that the solutions of IVP (2.8), which represent the outer free surface of an appropriate drop-like static meniscus, can be obtained only for those values of for which the solution is not globally convex.
Statement 2
If , then the solution of the IVP (2.8) is globally convex, and it does not represent the outer free surface of an appropriate drop-like static meniscus.
Comment 2. The above statement locates a set of values which are not appropriate for the creation of a drop-like meniscus.
Statement 3
If the solution of the nonlinear boundary value problem (NLBVP) (2.1), (2.4a), (2.4b), and (2.4c) represents the outer free surface of an appropriate concave static meniscus on the closed interval , with , then the following inequalities hold:
Comment 3. The above statement shows that in order to obtain the outer free surface of a concave static meniscus, appropriate for the growth of a tube of outer radius , has to be searched in the range defined by the inequalities (2.9).
Consequence 1
If ( is the middle point of the interval ), then and inequalities (2.9) become
Comment 4. The above consequence shows that in order to obtain the outer free surface of a concave static meniscus, appropriate for the growth of a tube of outer radius , has to be searched in the range defined by the inequalities (2.10).
Consequence 2
If verifies the inequality
then there is no in the closed interval for which the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) possesses a concave solution.
Comment 5. The above consequence shows that if is in the range defined by the inequality (2.11), then it is impossible to obtain a static meniscus having concave outer free surface which is appropriate for the growth of a tube of outer radius situated in the range
Consequence 3
If , then and, according to (2.9), .
Comment 6. The above consequence shows that in order to obtain the outer free surface of a concave static meniscus, appropriate for the growth of a tube of outer radius , has to be very high.
Statement 4
If and verify the inequalities
then there exist and a concave solution of the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) on the interval .
Comment 7. The above consequence shows that if and verify (2.12), then a static meniscus having concave outer free surface, appropriate for the growth of a tube of outer radius situated in the range , is obtained.
Consequence 4
If for and the following inequalities hold:
then there exist in the interval and a concave solution of the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) on the interval .
Comment 8. The above consequence shows that if is in the range defined by the inequalities (2.13), then a static meniscus having concave outer free surface, appropriate for the growth of a tube of outer radius situated in the range , is obtained.
Statement 5
A concave solution of the NLBVP (2.1), (2.4a), (2.4b), and (2.4c) is a weak minimum of the energy functional of the melt column.
Comment 9. The above consequence shows that a static meniscus having a concave outer free surface, appropriate for the growth of a tube of outer radius situated in the range , is stable.
Statement 6
If and satisfy the inequalities
then the solution of the IVP (2.8) is concave on the interval where is the maximal interval of the existence of .
Comment 10. The above consequence shows that if is in the range defined by the inequalities (2.14), then eventually it can be used for the creation of a drop-like static meniscus, appropriate for the growth of a tube of outer radius in the range .
Statement 7
If and for a value , which satisfies a static meniscus appropriate for the growth of a tube of outer radius exists, then for the following inequalities hold:
Comment 11. The above statement locates those values of for which eventually nonglobally concave (drop-like) static meniscus, appropriate for the growth of a tube of outer radius , exists.
Statement 8
If the solution of the IVP
is convex, then it does not represent the inner free surface of an appropriate drop-like static meniscus.
Comment 12. The above statement shows that the inner free surface of an appropriate drop-like static meniscus can be obtained only for those values of for which the IVP (2.16) is not globally convex.
Statement 9
If satisfies the inequality
then the solution of the IVP (2.15) is convex on the interval , where is the maximal interval of the existence of .
Comment 13. The above statement shows that if is in the range defined by the inequality (2.17), then it is impossible to obtain a static meniscus, appropriate for the growth of a tube of inner radius situated in the range .
Statement 10
If the solution of the NLBVP (2.2) and (2.5) represents the inner free surface of an appropriate concave static meniscus on the closed interval with , then the following inequalities hold:
Comment 14. The above statement shows that in order to obtain the inner free surface of a concave static meniscus, appropriate for the growth of a tube of inner radius , has to be searched in the range defined by the inequalities (2.18).
Consequence 5
If , then (i.e., is the middle point of the interval ) and the inequalities (2.18) become:
Comment 15. The above consequence shows that in order to obtain the inner free surface of a concave static meniscus, appropriate for the growth of a tube of inner radius , has to be searched in the range defined by the inequalities (2.19).
Consequence 6
If verifies:
then there is no in the closed interval for which the NLBVP (2.2) and (2.5) possesses a concave solution.
Comment 16. The above consequence shows that if is in the range defined by the inequalities (2.19), then it is impossible to obtain a static meniscus, having concave inner free surface and appropriate for the growth of a tube of inner radius situated in the range .
Consequence 7
If , then and .
Comment 17. The above consequence shows that in order to obtain the inner free surface of a concave static meniscus appropriate for the growth of a tube of inner radius , has to be very high.
Statement 11
If and verify the inequalities
then there exist in the closed interval and a concave solution of the NLBVP (2.2) and (2.5) on the interval .
Comment 18. The above statement shows that if and verify (2.21), then a static meniscus having concave inner free surface and appropriate for the growth of a tube of inner radius in the range , is obtained.
Consequence 8
If for and the following inequalities hold:
then there exist in the interval and a concave solution of the NLBVP (2.2) and (2.5) on the interval .
Comment 19. The above consequence shows that if is in the range defined by the inequalities (2.22), then a static meniscus having a concave inner free surface and appropriate for the growth of a tube of inner radius , situated in the range , is obtained.
Statement 12
A concave solution of the NLBVP (2.2) and (2.5) is a weak minimum of the energy functional of the melt column.
Comment 20. The above statement shows that a static meniscus having a concave inner free surface appropriate for the growth of a tube of inner radius , situated in the range , is stable.
Statement 13
If satisfies and there exists such that on the interval an appropriate static meniscus exists, then verifies:
Comment 21. Inequality (2.23) locates the values of for which eventually nonglobally concave (convex-concave) static meniscus can be obtained.
Statement 14
If and there exists such that on an appropriate non-globally concave static meniscus exists, then satisfies the inequalities
Comment 22. Inequality (2.24) locates the values of for which eventually concave-convex static menisci can be obtained.
Statement 15
If and satisfy the inequalities
then the solution of the IVP (2.16) is concave on the interval where is the maximal interval of the existence of .
Comment 23. Inequality (2.25) locates the values of for which we can obtain eventually appropriate static menisci.
3. Creation of an Appropriate Drop-Like Static Meniscus for the Growth of a Tube of Inner Radius and Outer Radius
In this section it will be shown in which kind the explicit formulas presented in the above section can be used for the creation of an appropriate drop-like meniscus when are given a priori. In the same time the melt column βheightβ , which has to be used, is found in function of the pressure of the gas flow introduced in the furnace for release the heat.
3.1. Creation of the Outer Free Surface
For the creation of the appropriate outer free surface, the following limits were considered: for
These limits are represented for in Figure 2. The computations were performed in MathCAD V.13 using the following numerical data When , that is, , we have the following.

More precisely, we have integrated the following system: for , and different . The obtained outer radii versus are represented in Figure 3, which shows that the desired outer radius is obtained for .

Actually, as it can be seen on the same figure, for we can obtain also a second outer radius , which is not anymore in the desired range . Moreover, the outer free surface of this meniscus is not globally concave; it is a convex-concave meniscus (Figure 4).

(a)

(b)
Taking into account [7, 23, 24], the melt column height in this case is , where is the pressure of the gas flow (introduced in the furnace for release the heat from the outer side of the tube wall). When [7, 24], then is negative, ; that is, the crucible melt level has to be with above the shaper top level. When , then ; that is, the crucible melt level has to be with above the shaper top level.
Additional remarks are as follows.
(i)For in the range , according to Statement 2, the outer free surface of the meniscus is globally convex, and, according to Statement 1, such a meniscus is not appropriate for the growth of a tube of outer radius . (ii)For in the range , according to Consequence 2, it is impossible to obtain concave outer free surface, which is appropriate for the growth of a tube of outer radius .(iii)For in the range , according to the Statement 3, it is impossible to obtain concave outer free surface, which is appropriate for the growth of a tube of outer radius .(iv)For in the range , according to the Statement 4, the outer free surface is concave and is appropriate for the growth of a tube which outer radius is higher than .3.2. Creation of an Appropriate Inner Free Surface
For the creation of the appropriate inner free surface, the following limits were considered for
These limits are represented for in Figure 5. The computations were performed in MathCAD V.13 using the same numerical data as for the creation of the appropriate outer free surface.

When , that is, , we obtain the following.
(i)If there exists a concave inner free surface for the growth of a tube of inner radius , then, according to the Statement 10, this can be obtained for a value of which is in the range .Taking into account the above fact, in order to create a concave inner free surface, appropriate for the growth of a tube which inner radius is equal to , we have solved the IVP (2.16) for different values of in the range . More precisely, we have integrated the following system: for , and different . The obtained inner radii versus are represented in Figure 6, which shows that the desired inner radius is obtained for .

Taking [7, 23, 24], the melt column height in this case is , where is the pressure of the gas flow (introduced in the furnace for release the heat from the inner side of the tube wall). When [7, 24], then is positive, ; that is, the crucible melt level has to be with under the shaper top level. When , then ; that is, the crucible melt level has to be with above the shaper top level.
Additional remarks are as follows.
(i)For in the range , according to Statement 9, the meniscus inner free surface is convex on the interval , and according to Statement 8, such a meniscus is not appropriate for the growth of a tube of inner radius .(ii)For in the range if an appropriate static meniscus exists, then the inner free surface of the meniscus is not globally concave.(iii)If is in the range and an appropriate nonglobally concave inner free surface exists, then is less than (iv)For in the range the obtained inner free surface is appropriate for the growth of a tube of inner radius which is less than the desired radius .3.3. Creation of a Concave Meniscus, Appropriate for the Growth of a Tube of Inner Radius and Outer Radius
For creating a concave meniscus, appropriate for the growth of a tube with outer radius β[m] and inner radius β[m] the melt column heights (with respect to the crucible melt level) have to be
In the following we will present three different procedures to create such a meniscus.
Procedure A [; ]
When the shaper outer top is at the same level as the shaper inner top, with respect to the crucible melt level, then the relation holds and hence: . It follows that the pressure of the gas flow, introduced in the furnace for releasing the heat from the inner wall of the tube, , has to be higher than the pressure of the gas flow introduced in the furnace for release the heat from the outer wall of the tube, ; .The results of the integration of the system (3.3) and (3.5) for and when the shaper outer top is at the same level as the shaper inner top, with respect to the crucible melt level, are represented in Figure 7.Figure 7 shows also that the inner free surface is higher than the outer free surface:

Procedure B [; ]
When the pressure of the gas flow introduced in the furnace for releasing the heat from the outer wall of the tube, , is equal to the pressure of the gas flow introduced in the furnace for release the heat from the inner wall of the tube , then the equalities and imply the equality
In other words, when , then the level difference between the shaper inner top and outer top has to be . In this condition the inner surface top level with respect to the shaper outer top level is (see Figure 8).

Hence the level difference between the inner free surface top and outer free surface top in this case is .
Due to the difficulties which can appear in this case in the creation of an appropriate thermal field, which assures the solidification in the right places, it is more appropriate to use and gas flows having the property (Procedure A) . For this choice the difference between the inner free surface height and outer free surface height with respect to the shaper top level is only , and the creation of the thermal field, which assures the solidification in the right places, can be easier.
Procedure C [; ]
Finally, to create an appropriate concave meniscus for which the inner meniscus is equal to the outer meniscus height, we will proceed by changing both the pressures , of the gas flows and the inner and outer top levels of the shaper.
For example, if we take for creating an appropriate outer free surface, then the crucible melt level has to be with above the shaper outer top level.
Now we take . After that, in order to satisfy , we take the pressure equal to . For this choice, we will have an appropriate drop-like meniscus for which , as it is shown in Figure 9, obtained by integration of systems (3.3) and (3.5) for the given numerical data.

4. The Effect of the Choice of the Shaper Radii
In this sequence we intend to evaluate numerically the effect of the choice of the shaper radii for the growth of a tube which outer radius is equal to , and inner radius is equal to .
For this purpose, the procedure described in sequence 3 will be applied for the growth of a tube having the same outer radius β[m] and inner radius β[m], but the shaper radii changed: the inner radius of the shaper now is (instead of ), and the outer radius of the shaper now is (instead of ).
Applying the procedure in this case we find that to create a static meniscus which outer free surface is appropriate for the growth, has to be , and to create a static meniscus which inner free surface is appropriate for the growth, has to be . The so-obtained outer and inner free surfaces are stable. The inner free surface height is , and the outer free surface height is .
When we take the shaper inner top at the same level as the shaper outer top, then the pressure difference between the pressure of the inner gas flow and the outer gas flow has to be and the level difference between the inner and outer free surface top is equal to .
When we take the pressure of the inner gas flow equal to the pressure of the outer gas flow , then the level difference between the shaper inner top and the shaper outer top has to be and the level difference between the inner and outer free surface top is equal to the difference
5. Conclusions
(i)Stable and drop-like static meniscus, appropriate for the growth of a single-crystal cylindrical tube with a priori specified inner and outer radii, can be created by the choice of the difference between the inner top level and outer top level of the shaper, or by the choice of the pressure of the gas flows, introduced inside and outside of the tube for release the heat from the tube walls.(ii)The values of the differences between the inner and outer top level of the shaper, or the pressure in the inner and outer gas flows, depend on the shaper radii. This dependence can be relevant.(iii)The realization of the thermal conditions, which assure that for the obtained static meniscus the solidification conditions are satisfied at the right places, is easier when we take the shaper inner top at the same level as the shaper outer top, and we use inner and outer gas flows having different pressures (). The setting of the thermal conditions is not discussed in this paper.Acknowledgments
The authors thank Professor Koichi Kakimoto for the useful discussions and suggestions concerning the tube growth process. They are grateful to the Romanian National Authority for Research for supporting the research under the Grant ID 354 no. 7/2007.