When capacitively coupled plasma (CCP) is used to generate large-area plasma, the standing wave effect becomes significant, which results in the hindering of the uniform plasma process such as in a plasma etcher or plasma chemical vapor deposition. In this study, the transmission line modeling method is applied to calculate the voltage distribution over atmospheric-pressure CCP electrodes with the size of 1 m × 0.2 m. The measured plasma impedance in our previous study was used in the present calculation. The results of the calculations clearly showed the effects of excitation frequency and the impedance of the plasma on the form of the voltage distribution caused by the standing wave effect. In the case of 150 MHz frequency, the standing wave effect causes a drastic change in the voltage distribution via plasma ignition; however, the change is small for 13.56 MHz. It was also clarified that the power application position is important for obtaining a uniform voltage distribution.
Keywords:CCP; Atmospheric-pressure plasma CVD; TLM; Voltage distribution
Plasma-enhanced chemical vapor deposition (PECVD) is an important and widely used process for forming various kinds of thin films in the electronics industry to fabricate, for example, very-large-scale integration and solar cells. For PECVD, capacitively coupled plasma (CCP) has the advantage of generating the large-area plasma necessary to process large substrates. However, when the electrodes become large relative to the wavelength of the electromagnetic wave used to generate the plasma, the standing wave effect will become significant, deteriorating the uniformity of the film thickness obtained [1-5].
It is considered that the voltage distribution over the CCP electrode greatly affects not only the distribution of plasma characteristics, such as plasma density and electron temperature, but also the deposited film thickness uniformity, especially in the case of PECVD. In this study, the voltage distribution over the CCP electrode was calculated by the transmission line modeling (TLM) method  for 13.56 and 150 MHz and with different power application positions. A two-dimensional (1 m × 0.2 m) plane electrode was modeled, and the impedances of the atmospheric-pressure plasma obtained from IV (current and voltage) measurements and analysis  were used for the calculation.
A one-dimensional model of electrodes and plasma (including sheath) is shown in Figure 1. Radio-frequency voltage is applied to the upper electrode, and the lower electrode is grounded. We assume that only the upper electrode has resistance and inductance for simplicity. This simplified model is useful enough to calculate a relative voltage between two electrodes, because only relative voltage is important for plasma.
Figure 1. One-dimensional model of electrodes and plasma (including sheath).
Plasma will be generated in the space between the upper and lower electrodes. In this model, electrodes (upper and lower) and plasma are divided into small elements of length ΔX. The voltage U is assumed to be constant within the elements. Symbol δ is the thickness of current flow (skin depth). The currents flowing into and out from the element are shown by the arrows in Figure 1. The plasma is assumed to be able to be represented by the parallel connection of the capacitance Cp and the conductance Gp.
One can derive a one-dimensional wave equation from the above mentioned one-dimensional model and extend it to the following two-dimensional wave equation. In the case of a two-dimensional model, the electrode will be divided in two directions, and the widths of the element will be ΔX and ΔY. For simplicity, the element widths ΔX and ΔY were set to be equal.
Here, L and R are inductance and resistance per unit length (in current flow direction) of the electrodes of element width (ΔX or ΔY), and Cp and Gp are the capacitance and conductance of plasma per unit length of element width, respectively. F(x,y,t) is the external force (causes voltage to change) applied to the upper electrode at position (x,y).
Electrode resistance R and inductance L
When radio-frequency power is applied to the electrodes, the current will flow only on the surface of the electrodes owing to the skin effect. The effective electrode resistance per unit length R (of width w) is determined by the following equation :
where σ is the conductivity of the electrode material, δ is the skin depth, and w is the width of the current flow. The skin depth δ is determined by
where ω is the angular frequency, and μ is the magnetic permeability of the electrode material.
The inductance of a pair of two parallel plates (electrodes) per unit length (of width w) can be calculated using 
where d is the distance between the upper and lower electrodes, and w is the width of the current flow. When aluminum is used as the electrode material, the conductivity σ is 0.33 × 108 ohm−1 m−1.
Table 1 shows a summary of these values calculated for the cases of 150 MHz and 13.56 MHz.
Plasma conductance Gp and capacitance Cp
In the case of atmospheric-pressure plasma, since the gap between the electrodes is usually too narrow (≤1 mm) to perform Langmuir probe analysis, we performed plasma impedance analysis in our previous study . A combination of the measurement of the current and voltage waveforms outside of the apparatus and calculation using the electrically equivalent circuit model enabled us to derive the impedance Zp of the plasma-filled capacitor. Figure 2 shows the measured impedance of atmospheric-pressure helium plasma (real (Figure 2a) and imaginary (Figure 2b) parts of Zp) as a function of applied power density, for 150 MHz and 13.56 MHz excitations using a metal electrode with a diameter of 10 mm and a gap of 1 mm. As shown in Figure 2, the plasma impedance Zp changes depending on the applied power; this is known as a nonlinear characteristic of the plasma. However, it is also shown that the impedance becomes constant (the system is linear) in a considerably wide power range when sufficiently high power is applied to the plasma. Although taking the nonlinear characteristic of plasma into account will give more exact results, we consider that it is still meaningful to calculate the voltage distribution on the assumption that the plasma impedance is constant, since plasma equipment is often used in such a saturated area.
Figure 2. Real (a) and imaginary (b) parts of plasma impedance vs. applied power density. Electrode diameter, 1 cm; electrode gap, 1 mm.
The plasma conductance Gp and the susceptance Bp per unit length of element width are calculated from a given plasma impedance Zp (Zp = Rp' − Xpj) using
Then the plasma (parallel) capacitance Cp per unit length of element width at a particular frequency ω (shown in Figure 3) can be calculated from plasma susceptance Bp, as
Figure 3. Conversion of plasma impedance (left) to admittance (right).
Wavelength and phase velocity in the electrodes
The propagation constant γ ≡ α + βj of the solution of Equation 1 is
Its real part α (attenuation coefficient) and imaginary part β (phase propagation constant) are described as
The phase velocity v of the electromagnetic wave propagating in the system described by Equation 1 is
The wavelength λ is calculated using
From these equations, it is clear that the wavelength on the electrode is governed not only by the electrode configuration but also the impedance of plasma. Both the attenuation coefficient α and the wavelength λ greatly affect how a standing wave is formed on the electrode.
Results and discussion
Equation 1 can be numerically solved by a finite differential method. Calculation was performed for an electrode size of 0.2 × 1 m2, with the edges of the electrodes assumed to be open. Usually, plasma equipment is designed so that the edge of the electrode is not exposed to the plasma. Sometimes, the edges of the electrode will be supported by dielectric materials such as quartz and ceramics, in which case the edges are terminated by the capacitance formed by the dielectrics. In such a case, in order to minimize the power loss, the electrode supporting system will be designed so that the capacitance becomes as small as possible, in which case the impedance is close to that of the open case. The electrode was divided into small elements of which the size is 0.01 × 0.01 m (ΔX = ΔY = 0.01 m). Both Cp and Gp are assumed to stay constant with relatively small variation in the electrode voltage. Cp and Gp values were calculated from the measured impedance of atmospheric-pressure helium plasma (Zp) shown in Figure 2.
Table 2 shows the plasma impedance Zp, admittance Yp, and (parallel) capacitance Cp used for the calculations. The propagation constant γ and the wavelength λ are also shown. It is seen that the wavelength λ on the electrode is considerably shorter than that in free space.
Figure 4 shows the calculated two-dimensional distribution of the voltage amplitude at each point on the electrode during plasma generation. The power was applied at the center of the electrode.
Figure 4. Two-dimensional distribution of voltage amplitude on the electrode during plasma generation. Power was applied at the center of the electrode. (a) 150 MHz and (b) 13.56 MHz.
The central cross-sectional distributions of the plots in Figure 4 are shown in Figure 5, where voltage distribution is along the central cross-sectional line in the direction of electrode length. Voltages oscillate between their maximum and minimum with the driving frequency. Dotted lines in Figure 5 show instantaneous voltage profiles at elapsed times of 9.35 and 181.77 ns for 150 and 13.56 MHz, respectively. They always remain between the maximum voltage (upper solid line) and the minimum voltage (lower solid line). It is clearly seen that voltage variation is considerably larger for 150 MHz than for 13.56 MHz. The voltage variation over the electrode is approximately 58% and 12% for 150 and 13.56 MHz, respectively.
Figure 5. Voltage distributions along the central cross-sectional line on the electrode. Power was applied at the center of the electrode. (a) 150 MHz and (b) 13.56 MHz.
The spatial differentiation of the instantaneous voltage gives the electric field in the horizontal direction at the differentiation point. As shown in Figure 5, the gradient of the instantaneous voltage is largest at the driving point. According to the calculation, the largest gradient of the instantaneous voltage in 150 MHz case was approximately 0.45 V/m, while the average electric field across the electrodes was 5,000 V/m. This means that the current flowing in the horizontal direction is small enough compared with that flowing in the vertical direction. Since the difference was even larger in the 13.56 MHz case, the current flowing in the horizontal direction can be neglected.
Very different voltage distribution profiles are obtained when radio-frequency power is applied on both ends of the electrode, as shown in Figure 6. The phase of radio frequency was set to be the same. The voltage variations over the electrode are approximately 39% and 11% for 150 and 13.56 MHz, respectively. Therefore, this type of power application would be more advantageous for obtaining more uniform plasma over the electrode.
Figure 6. Voltage distributions along the central cross-sectional line on the electrode during plasma generation. Power was applied on both ends of the electrode with the same phase. (a) 150 MHz and (b) 13.56 MHz.
Figure 7 shows the results of the calculations of voltage distribution before plasma ignition. When there is no plasma between the electrodes, the conductance G is zero and the capacitance C is determined by
where ϵ0 is the permittivity of vacuum. S and d are the electrode area and the distance between the upper and lower electrodes, respectively.
Figure 7. Voltage distribution on the electrode before plasma ignition. Power was applied at the center of the electrode. (a) 150 MHz and (b) 13.56 MHz.
Comparing Figure 7 with Figure 5, a slight difference is seen in the case of 13.56 MHz. When 150 MHz is applied, however, the voltage distribution before plasma ignition is considerably different from that after plasma ignition.
From the attenuation coefficient α shown in Table 2, the resistive loss in the 150 MHz case is larger than that in the 13.56 MHz case. However, the resistive loss only causes a monotonic decay in voltage amplitude from the driving point along the wave-propagation direction. Since Figure 5 does not show a monotonic decay in voltage from the driving point, the drastic change in the voltage pattern in the 150 MHz case is considered to be caused mainly by the standing wave effect.
The interference pattern may change sensitively with the changes in various parameters (e.g. electrode shape, setup, and plasma parameters) in the case of 150 MHz. It can be said that in the case of 13.56 MHz, the expected or measured voltage distribution before plasma ignition is useful for designing the electrode setup. However, in the case of 150 MHz, careful design of the electrode setup should be required to obtain stable and uniform plasma generation.
A mathematical model for calculating the voltage distribution over a parallel-plate plasma electrode was constructed by the TLM method. In this study, driving frequencies of 150 MHz and 13.56 MHz were compared. Actually measured atmospheric-pressure helium plasma impedance was used for these calculations. In the case of 150 MHz frequency, the standing wave effect caused a drastic change in the voltage distribution on the electrode by plasma ignition; however, the change was small for 13.56 MHz. Thus, in the case of 13.56 MHz, the expected or measured voltage distribution before plasma ignition is useful for designing the electrode setup. However, in the case of 150 MHz, careful design of the electrode setup should be required to obtain stable and uniform plasma generation. It was also shown that the power application position is important for obtaining uniform voltage distribution. It is considered that the voltage distribution will greatly affect the plasma density distribution and therefore film thickness uniformity in the case of plasma CVD. The TLM method is applicable to circular electrodes as well, and not only to atmospheric-pressure plasma but also to low-pressure plasma. The simulation by the TLM method will be useful in optimizing the configurations of parallel-plate plasma systems.
The authors declare that they have no competing interests.
MS developed the finite differential element code, ran the simulations, and wrote the manuscript. FT, HO, HK, and KY assisted in designing the work, discussed the results, and proofread the manuscript. All the authors read and approved of the final manuscript.
This work was supported in part by Grants-in-Aid for Scientific Research [nos. 20676003, 21656039, 22246017, and Global COE Program (H08)] from the Ministry of Education, Culture, Sports, Science and Technology, Japan.
Plasma Sources Sci Technol 1997, 6:170-178. Publisher Full Text
Plasma Sources Sci Technol 2004, 13:436-445. Publisher Full Text
Thin Solid Films 2012, 520:4219-4236. Publisher Full Text
Hiroaki K, Hiromasa O, Kiyoshi Y: High-rate and low-temperature film growth technology using stable glow plasma at atmospheric pressure. In Materials Science Research Trends. Edited by Olivante LV. New York: Nova; 2008::197.