Ultra-high-purity (UHP) gases are widely used in semiconductor fabrication plants (fabs), both as reactants in some processes like deposition and etching   , and as inert gases for purging various parts of the system. For example, UHP nitrogen is used in very large quantities in semiconductor fabs as diluent and purge medium. The UHP gases are typically delivered through pipes from gas storage facilities outside the fab to the point of use (POU). UHP gases, such as nitrogen, are usually stored in the cryogenic tanks where gas is produced by controlled vaporization of the stored liquid and transported through delivery pipes to the fab. Most semiconductor processes are very sensitive to impurities and require a well-controlled quality   . Failure to maintain a stable purity level at the (POU) can lower the performance of process tools and even slow or shut down the processing line  .
Moisture, selected as the model impurity species in this study, is one of the most common impurities. While affecting many fabrication processes even at the parts per billion (ppb) levels, it is also one of the most challenging impurities to remove and control due to its strong adsorption on the surfaces of transport line, components, and reactors    . Additionally, when gases are supplied from cryogenic sources, there is usually some drift and variation in the moisture level at the POU due to two major factors  :
1) Source drift: This is caused by the build-up of moisture in the cryogenic tanks due to the selective retention of moisture in the liquid phase during vaporization. As impurity level goes up with time and gas usage, and the gas purity can no longer satisfy the process requirements (usually after 70% - 80% of a tank is used), the remaining content is rejected and wasted, causing financial and environmental losses.
2) Temperature-induced variations: This is caused by adsorption (capture) and desorption (release) of moisture on the surfaces of delivery pipes and flow-control devices. The delivery pipes in most fabs are long and therefore a major sink for moisture adsorption. Moreover, these pipes, running from the storage area to the POU are usually in the open environment. Therefore, ambient temperature variations would have a large effect on the moisture adsorption and desorption, leading to drifts and spikes in impurity at the POU.
Single tank systems are widely used as cryogenic sources for many gases in industry today. However, this configuration results in impurity drift due to selective retention and accumulation of impurities in one phase during vaporization. For example, moisture is selectively retained in the liquid phase of most inert and process gas tanks, resulting in impurity drift upward, as liquid is vaporized and gas is withdrawn from the cryogenic tank. The drift necessitates shutting down the supply when the impurity level reaches the highest acceptable level. At this time, the partially used tank needs to be returned to the supplier, often with large amount of gas remaining unused in the tank. This increases the cost and the negative environmental impact of the operation due to the partial usage of tank content.
In recent years, some modifications have been introduced in cryogenic tank design to reduce the effect of source drift  . These changes increase the cost and the complications of supply tanks significantly. However, there has been no systematic study of disturbances caused by pipelines or developments of techniques to mitigate them. Purifiers can be used to remove moisture at the POU  ; however, they are expensive and often unreliable under transient flow conditions. A previous study by our group focused on some preliminary study of the effect of pipeline geometry and temperature variations on the moisture level at the POU  . The objective of the current work has been to develop a method for mitigating the effects of the variations in impurity levels, caused by the above two factors and assuring a stable and well-controlled purity level at the POU.
2. Method of Approach
The proposed concept and the associated system in this study are schematically shown in Figure 1. The system, chosen to illustrate the method, consists of two cryogenic tanks, two mass flow controllers (MFC), a gas delivery system, a real-time sensor for on-line measurement of impurities, and the auxiliary electronics for running the process simulator, data acquisition, and MFCs programing and control. The flow streams from the two tanks are controlled by two MFCs, then mixed and transported to the POU through a transport pipeline. The impurity level at the POU is monitored by an online sensor. The process simulator provides input to the MFC controller unit that adjusts the flow out of each MFC. Overall, the set-up provides a dynamic flow-mixing scheme run by the output from the process simulator. The system properties, operating conditions and purity requirements are all inputs to the process simulator. The details of the process simulator and its functions are described in the following two sections: first one for cancelling the drift in the cryogenic tanks and the second one for mitigating the effects of the pipeline temperature variations.
2.1. Multiple Tanks to Resolve the Source Drift Problem
As pointed out in the Introduction section, the single-tank configuration results in impurity drift due to the accumulation of impurities like moisture in the liquid phase during vaporization. This leads to partial usage of tank content and wasting of gas. To resolve this drift issue, a dual tank system is proposed, where flow rates from two tanks are varied and mixed so that a stable and well-controlled purity level is obtained.
Assuming that gas and liquid in the tank are in equilibrium, the following equations determine the operation of tank 1:
Figure 1. Schematic of proposed system.
where , , and are the total tank volume, the gas phase volume, and the liquid phase volume in the tank, respectively. and are the impurity concentrations in the vapor and liquid phases in the tank, is the total amount of impurity in the tank, Q is the withdrawn flow rate, and are the mole fractions of impurity in the vapor and the liquid phases, H is Henry’s law constant, and α is the ratio of liquid to vapor densities. Equations (2) and (3) together are the mass balance for the impurity in the tank; Equation (4) shows the liquid-vapor equilibrium for the impurity, and Equation (5) is the overall balance for gaseous content. Solving these five equations simultaneously gives:
The total flow rate, Q, which is the sum of the flow rates from the tanks is given by the process demand at the POU. The flow rate from the first tank, , as a fraction of Q, is given by:
Equations (6) and (7) can be modified as
Similarly, the equations in tank 2 are
The impurity concentration in the mixed flow from two tanks is
Given any desired function, the mixing function, , is given by:
2.2. Programmed Mixing to Mitigate the Temperature-Induced Variations
The analysis of the processes taking place in the pipeline is through solving the conservation equations of transport and reactions in the pipe. However, application of the conventional numerical approach in this particular case is complicated since the inlet concentration is not known; the known boundary condition is the outlet and not the inlet concentration. The proposed solution to this numerical issue is the use of an equivalent system model consisting of a series of well-mixed cell, shown in Figure 2, where the number of cells is selected to represent the extent of dispersion in the pipe. This modeling approach is widely used in the analysis of non-ideal tubular reactors  .
The mass balance for impurity in the gas phase for the first cell is given by
The mass balance for any intermediate cell i in the chain is
and the corresponding equation for the last cell is
The mass balance equations for the impurity adsorbed on the pipe surface are
for any cell i in the chain (19)
where and are the inlet and the outlet impurity concentrations of the first cell; and are the inlet and the outlet impurity concentrations of cell i; similarly, and are the inlet and outlet impurity concentrations of the last cell. , , and are the adsorbed concentrations in these cells; is the adsorption rate constant, is the desorption rate constant, and is the surface site density in each cell; S is the surface area of each cell, V is the cell volume, and Q is the total volumetric flow rate.
The initial conditions of Equations (15) to (20) are
Figure 2. Schematic of mixing cells model.
at , (22)
where is the initial adsorbed concentration, is the initial gas phase concentration, and is the adsorbed concentration in equilibrium with .
The adsorption rate constant, , and the desorption rate constant, , are functions of temperature, given by
where and are the activation energies for adsorption and desorption, R is gas constant, and is ambient temperature in the pipe.
When the mixed flow enters the pipe, its temperature rises from tank temperature, , to the pipe temperature which is assumed to be the ambient temperature, . Therefore, the mixed concentration, , is related to the inlet concentration of cell 1, , by
and the change in the outlet concentration of the last cell, , with temperature is given by
where is the initial temperature in the pipe.
For a given desired , the function , which is the ratio of flow rate from tank 1 to the total flow rate, is determined by the simultaneous solution to the above equations.
3. Results and Discussion
3.1. Two-Tank Configuration for Canceling the Impurity Drift
Nitrogen and moisture are selected as the model main gas and impurity, respectively in this example. However, the methodology and the process simulator tool are not limited to these two compounds and can be applied to other gases as well. In this example, two tanks are used simultaneously as the cryogenic source for nitrogen. Initially, tank 1 is a partially used and tank 2 is a new full tank. Since the typical gas usage in semiconductor processing tools is transient and periodic, a cyclic flow rate is used to represent the total gas supplied by the system.
In general, the present analysis can be used to determine the mixing function needed for obtaining any desired impurity concentration in the mixed flow from the cryogenic tanks, . In this example, the desired concentration, , is assumed to be constant. The parameters used in this example are given in Table 1. Henry’s law constant of 0.39 was calculated
Table 1. Parameters in two-tank model study.
from experimental data presented in previous study  . It is assumed that the flow rate needed at the POU alternates between a low value and a high value at regular intervals. This kind of variation is typical for most semiconductor fabrication applications and tools. Using the analysis method of Section 2.1, the mixing function, , that gives a constant mixed concentration , is determined; the results are shown in Figure 3.
The mixing function profile shows that initially , indicating that the mixed flow is fully supplied from tank 1. As time goes on, decreases because the moisture level in tank 1 increases and more flow from tank 2 is required to compensate for this increase. After about 44 hours, is nearly zero, indicating that the mixed flow is almost provided from tank 2. At this time, moisture level in tank 1 is too high and the desired mixed concentration cannot be obtained by mixing the flow from these two tanks. Therefore, a tank change is needed at this time to replace the exhausted tank 1 with a new tank.
Figure 4 shows liquid volume in two tanks during this 44-hour cycle before tank replacement. At the end of this cycle, the liquid volume in tank 2 decreases from 30 m3 to 12.1 m3, which is almost the same volume as the initial liquid volume (12 m3) in tank 1. The liquid volume in tank 1 goes from 12 m3 to 2.4 m3 after 44 hours. The 44-hour cycle represents the lifetime of tank 1 before replacing it with a new tank. At this time, tank 2 is in conditions similar to that of tank 1 at the beginning of the cycle. Therefore, a new cycle starts and repeats the same changes. After each cycle, the tank which runs out of liquid nitrogen will be taken out for recharge or replacement.
To illustrate the advantages of the proposed system configuration in comparison with a typical conventional supply system, the time profiles of impurity in both cases are shown in Figure 5. In the conventional single-tank system, a stabilized moisture level is not achievable, and the impurity concentration goes up with
Figure 3. Temporal profile of the mixing function, .
Figure 4. Temporal profiles of the liquid content of the two tanks.
Figure 5. Comparison of the impurity profiles in two-tank and single-tank configurations.
time; after around 28 hours, the moisture level exceeds 35 ppb. However, using the proposed controlled mixing with the appropriate mixing function, , a stream with constant concentration of 35 ppb is obtained and maintained for the entire cycle time. Moreover, the results show that the two-tank configuration with the proposed mixing technique results in 92% usage of the liquid nitrogen in each tank compared with less than 80% usage of liquid nitrogen in a conventional single-tank system.
3.2. Application of the Proposed Scheme to Cancel the Combined Effects of Tank Drift and Pipe Temperature Variations
In this example, a more complete system is analyzed, where the delivery tanks are followed by a pipeline that takes the gas from the tanks to the POU. Therefore, the transport processes, the adsorption and desorption effects, and the temperature changes in the ambient around the pipe all contribute to the impurity variations at the POU. For illustration, a desired constant moisture level of 40 ppb is assumed at the POU, which is the pipe outlet in this case. A cyclic flow rate between a high value and a low value is used to represent the POU requirement. In one cycle, lasts 20 minutes, and lasts 60 minutes. The cyclic ambient temperature (diurnal effect) is assumed to follow:
The initial moisture concentration C0 is 7.79E−6 mol/m3; therefore, the desired concentration at POU is:
The mixing function is calculated in two steps. In the first step, the inlet concentration, , is determined using the mixing cells model described in Section 2.2. The parameters of the system used in this example are listed in Table 2.
In the first step, the pipe inlet concentration, , is calculated from the mixing cells model of Section 2.2. starts as and then changes due to temperature variations, going down as temperature goes up. The comparison of results for 3 and 4 tanks in Figure 6 shows that there is not a significant difference in the final results between these two cases.
In the second step, the procedure described in Section 2.1 is used to find from , calculated in step 1. Due to the presence of the pipeline, the system configuration in this case is different from that of example in Section 3.1; therefore, the resulting mixing function, , will also be different.
In this example, the cycle starts with one tank initially full and one tank with 1/3 liquid nitrogen remaining. The used tank is called tank 1, and the full tank is labeled tank 2. Tank 1 starts with the moisture level of 40 ppb, which is the same as the initial moisture level in the delivery pipe. The flow from tank 2 has the initial moisture level of 20 ppb. Using , obtained from the analysis of the pipe
Table 2. Parameters in combination of models.
Figure 6. Comparison of Cin profiles obtained by assuming 3 and 4 cells in modeling the pipeline.
model, is calculated using Equation (25). The calculated mixing function, , and the liquid volume in tanks are shown in Figure 7 and Figure 8. The results show that at the end of 16 days, b is nearly zero. The liquid volume in tank 2 drops from 30 m3 to 9.77 m3, which is close to the initial liquid volume (9.65 m3) in tank 1. During this time, the liquid volume in tank 1 goes from 9.65 m3 to 0.80 m3 after 16 days. Therefore, 16 days represents one cycle and the time for replacement of the exhaust tank. Around 97% of liquid nitrogen will be used after one cycle, which represents a much better usage efficiency than that from
Figure 7. Temporal profile of the mixing function, .
Figure 8. Temporal profiles of the liquid content of the two tanks.
the conventional single tank configuration.
3.3. Parametric Study of Temperature
An important application of the methodology developed in this study is in parametric study for optimizing the operation of an existing system or designing a new system. Using the proposed approach, the analysis of a system over a wide range of configurations and operating conditions is relatively fast and inexpensive compared to running experiments on a physical set up. For example, to see the effect of temperature variations, three types of temperature amplitudes (10, 20, and 30 degrees centigrade) were compared, using the same temperature functionality given by Equation (27). The time profiles of inlet concentration, , mixing function, , and liquid volumes, and are shown in Figures 9-11, respectively. Results in Figure 9 show that a larger temperature variation results in a larger variation in the impurity concentration, primarily due to a larger effect of adsorption and desorption; it will also increase the extent of
Figure 9. Temporal profiles of Cin under various temperature amplitudes.
Figure 10. Temporal profiles of mixing function under various temperature amplitudes.
Figure 11. Temporal profiles of liquid content under various temperature amplitudes.
change in the mixing function, .
Most semiconductor fabrication processes are very sensitive to impurities in purge and process gases. Therefore, maintaining a stable purity level in gas supply systems is of critical importance for this industry. Moisture is one of the most common and also most problematic impurities in the ultra-high-purity gases. Most UHP gases are supplied from cryogenic storage tanks and transported to the point of use (POU) by pipelines that are exposed to ambient environment. This arrangement results in variations of moisture level delivered to the POU due to two primary factors: The first is the drift due to the phase segregation of impurities during liquid vaporization in the cryogenic tanks. Moisture in most cryogenic systems segregates favorably in the liquid phase. This results in the retention and accumulation of moisture in the liquid phase and a gradual increase in the impurity in the gas delivered from these sources. The second factor is the variation of the impurity due to the adsorption and desorption on the surfaces of transport pipelines exposed to the varying outside temperature.
In this study, a new approach is proposed and analyzed that will cancel the variations caused by both of these two factors. In this approach, the stability and high purity at the POU are achieved through controlled mixing of gases from two cryogenic tanks. The mixing ratio is a dynamic function of time that is found through the application of a process simulator developed in this work. The simulator includes both the system geometry and configuration as well as the operating conditions and variations of external factors such as the ambient temperature. The proposed approach is expected to result in significant savings in the usage of gases. For example, compared to the conventional single-tank method, which can only use less than 80% of liquid nitrogen from a liquid nitrogen tank, the proposed method increases the liquid nitrogen utilization to over 90%.
The formulation of the simulator is generally valid for most cryogenic gas supply systems. However, the continuation of this work should include application to the control of impurities in reactive gases such as moisture in ammonia, which is another key issue in the semiconductor fabs.
Financial and technical supports from Intel Corporation and the SRC Engineering Research Center for Environmentally Benign Semiconductor Manufacturing are acknowledged.
: ratio of liquid to vapor densities;
: mixing function;
: initial gas phase concentration in the pipeline, (mol/m3);
: outlet impurity concentration of the first cell, (mol/m3);
: inlet impurity concentration of the last cell, (mol/m3);
: outlet impurity concentration of the last cell, (mol/m3);
: inlet impurity concentration of cell i, (mol/m3);
: outlet impurity concentration of cell i, (mol/m3);
: inlet impurity concentration of the first cell, (mol/m3);
: impurity concentration of the liquid phase in tank 1, (mol/m3);
: adsorbed concentration in the first cell, (mol/m2);
: adsorbed concentration in the last cell, (mol/m2);
: adsorbed concentration in cell i, (mol/m2);
: initial adsorbed concentration, (mol/m2);
: mixed concentration, (mol/m3);
: impurity concentration of the vapor phase in tank 1, (mol/m3);
: initial impurity concentration of the vapor phase in tank 1, (mol/m3);
: impurity concentration of the vapor phase in tank 2, (mol/m3);
: initial impurity concentration of the vapor phase in tank 2, (mol/m3);
d: pipe diameter, (m);
: activation energy for the adsorption, (kcal/mol);
: activation energy for the desorption, (kcal/mol);
H: Henry’s law constant;
: adsorption rate constant, (m3/mol・s);
: adsorption rate constant at 298 K, (m3/mol・s);
: desorption rate constant, (1/s);
: desorption rate constant at 298 K, (1/s);
L: pipe length, (m);
: total amount of impurity in tank 1, (mol);
P: pressure, (Pa);
Q: volumetric flow rate, (m3/s);
: flow rate from tank 1, (m3/s);
: high flow rate, (m3/s);
: low flow rate, (m3/s);
R: gas constant, (J/mol・K);
S: surface area of each cell, (m2);
: surface site density in each cell, (mol/m2);
: running time for high flow rate, (min);
: running time for low flow rate, (min);
: initial temperature in the pipe, (K);
: ambient temperature in the pipe, (K);
V: cell volume, (m3);
: gas phase volume in tank 1, (m3);
: initial gas volume in tank 1, (m3);
: gas phase volume in tank 2, (m3);
: initial gas volume in tank 2, (m3);
: liquid phase volume in tank 1, (m3);
: liquid phase volume in tank 2, (m3);
: total tank volume, (m3);
: total volume in tank 1, (m3);
: mole fraction of impurity in the vapor phase in tank 1;
: mole fraction of impurity in the liquid phases in tank 1.
 Watanabe, T., Funke, H.H., Torres, R., Raynor, M.W., Vininski, J. and Houlding, V.H. (2003) Contamination Control in Gas Delivery Systems for MOCVD. Journal of Crystal Growth, 248, 67-71. https://doi.org/10.1016/S0022-0248(02)01889-4
 Fine, S.M., McGuire, J.T., Choi, B.S., Bzik, T., Crofton, K.T., Melnyk, A., Perez, M. and Sheriff, D.P. (1996) Optimizing the UHP Gas Distribution System for a Plasma Etch Tool. Solid State Technology, 39, 71-77.
 Haider, A.M. and Shadman, F. (1991) Desorption of Moisture from Stainless Steel Tubes and Alumina Filters in High Purity Gas Distribution Systems. IEEE Transactions on Components, Hybrids, and Manufacturing Technology, 14, 507-511. https://doi.org/10.1109/33.83935
 Krishnan, S. and Laparra, O. (1997) Contamination Issues in Gas Delivery for Semiconductor Processing. IEEE Transactions on Semiconductor Manufacturing, 10, 273-278. https://doi.org/10.1109/66.572082
 Yao, J., Wang, H., Dittler, R., Geisert, C. and Shadman, F. (2010) Application of Pressure-Cycle Purge (PCP) in Dry-Down of Ultra-High-Purity Gas Distribution Systems. Chemical Engineering Science, 65, 5041-5050. https://doi.org/10.1016/j.ces.2010.06.002
 Dittler, R., Shadman, F., Guzman, R. and Muscat, A. (2014) Reducing Ultra-High-Purity (UHP) Gas Consumption by Characterization of Trace Contaminant Kinetic and Transport Behavior in UHP Fabrication Environments, ProQuest Dissertations and Theses.
 Vininski, J., Yilcelen, B., Torres, R. and Houlding, V. Bulk Specialty Gas Systems for Ammonia Gas Phase Delivery vs. Liquid Phase Delivery. http://www.mathesongas.com/electronics/technical-library
 Wu, J. and Shadman, F. (2018) Impurity Drift and Variations in High-Purity Gas Delivery Systems. International Journal of Engineering and Mathematical Modelling, 2018, 1-9. http://www.orb-academic.org/index.php/journal-of-engineering-modelling/issue/view/44
 Waldram, S.P. (1985) Non-Ideal Flow in Chemical Reactors. In: Waldram, S.P., Ed., Comprehensive Chemical Kinetics, Elsevier, Amsterdam, Vol. 23, 223-281. https://doi.org/10.1016/S0069-8040(08)70148-1