A Linear Microwave Plasma Source Using a Circular Waveguide Filled with a Relatively High-Permittivity Dielectric: Comparison with a Conventional Quasi-Coaxial Line Waveguide

: For a conventional linear microwave plasma source (LMPS) with a quasi-coaxial line transverse electromagnetic (TEM) waveguide, a linearly extended plasma is sustained by the surface wave outside the tube. Due to the characteristics of the quasi-coaxial line MPS, it is easy to generate a uniform plasma with radially omnidirectional surfaces, but it is difﬁcult to maximize the electron density in a curved selected region. For the purpose of concentrating the plasma density in the deposition area, a novel LMPS which is suitable for curved structure deposition has been developed and compared with the conventional LMPS. As the shape of a circular waveguide, it is ﬁlled with relatively high-permittivity dielectric instead of a quasi-coaxial line waveguide. Microwave power at 2.45 GHz is transferred to the plasma through the continuous cylindrical-slotted line antenna, and the radiated electric ﬁeld in the radial direction is made almost parallel to the tangential plane of the window surface. This research includes the advanced 3D numerical analysis and compares the results with the experiment. It shows that the electron density in the deposition area is higher than that of the conventional quasi-coaxial line plasma MPS.


Introduction
Microwave plasma sources (MPSs) are driven by applied frequencies of 1-100 GHz and lead to electron heating and excitation of a discharge from the electromagnetic waves that either penetrate the plasma or exist along the surface [1,2]. With weakly ionized MPS used for material processing, waves are important to transfer energy from the excited surface of the waveguiding structure to the bulk plasma, where the wave energy is absorbed [3]. The main element of the MPS, for material processing, is the microwave-to-plasma applicator because it determines the structure of the electromagnetic field in the plasma as well as the absorption energy efficiency of plasma [4]. The structure of the electromagnetic field in plasma is determined by the correlation of the internal plasma parameters (electron density, plasma frequency, etc.) and the external source parameters (waveguide structure, size of source, frequency, etc.). The energy absorption efficiency of MPS is determined by the wave penetration depth in the plasma, wave field strength at the plasma surface, and waveguiding structure.
Among the MPS for material processing, the sources generally have a form in which waves propagate along the transmission lines. The plasma is generated by the waves emitted from the dielectric surface which replaces a part of the waveguiding structure. Depending on the size and structure of the dielectric surface, plasma sources using the transmission line field can be classified in the following ways: (i) plasma is a part of the waveguiding structure (transmission line applicators) and (ii) plasma is sustained transmission line field can be classified in the following ways: (i) plasma is a part of the waveguiding structure (transmission line applicators) and (ii) plasma is sustained by the field leaking out of the structure (antenna applicators) [5]. For the transmission line applicator, the cross section of the plasma-dielectric contact is wide, whereas for the antenna applicator, the plasma-dielectric contact cross section is narrow. So, in the case of antenna applicator, the process stability is high since the wave propagation mode and the resonant cavity in the waveguiding structure are less affected by plasma generation. Hence, in material processing, which requires stable plasma operation, antenna applicator type MPSs are preferred.
Among the antenna apllicator type MPS, LMPSs are mainly used for large area processing [6]. Erstwhile LMPS models [7] were either in the form of discharge tubes placed inside the guiding structure or in the form of discharge tubes across the wide waveguide wall parallel to the electric field [8,9]. However, in terms of productivity improvement and large area material processing, the shape of the discharge tubes inside the guiding structure is not suited for the application, so the shape of the discharge outside the guiding structure is demanded. A SW sustained source called the duo-plasma line source, which uses the quasi-coaxial line TEM wave mode, was invented according to the shape of the discharge outside the tube with a waveguiding structure [10,11]. Due to the high scalability with a large processing area [12,13], duo-plasma line source has been applied to many plasma enhanced chemical vapor deposition (PECVD) processes such as hard coating, thin film, synthesis of graphene, and display panels, which require a large area, an operation pressure of 0.2~3Torr, and high speed processes [14][15][16][17].
In the quasi-coaxial line TEM wave mode, the isotropically emitted transverse wave field from the inner conductor is perpendicular to the tangential plane of the window surface in all radial directions, and the dielectric surface in all radial directions is exposed to interact with the plasma generation. Therefore, it is easy to generate a uniform plasma in the radially omnidirectional plasma-dielectric window surfaces. However, it is difficult to maximize the process productivity by concentrating the plasma density in a onedirectional selected deposition region. In this study, for the purpose of concentrating the plasma density in the selective area, a cylindrical slotted LMPS that had a high electron density in a specific deposition region was proposed and was compared with the quasicoaxial LMPS. The shape of the discharge outside the guiding structure [18][19][20], cylindrical-slotted LMPS was excited by feeding power in the circular TE11 mode.

Plasma Source Overview
Both quasi-coaxial line TEM and circular TE11 LMPSs were fabricated to compare the characteristics of plasma generation according to the waveguiding structure and the wave field direction. The circular TE11 wave mode was compared with quasi-coaxial line TEM mode, as shown in Figure 1. In the figure, the dashed and straight lines represent the magnetic and electric fields, respectively.  In the circular waveguide, with the diameter determined, the wave frequency must exceed the cut-off frequency of propagation. Among the lower order modes of circular waveguides, there are the TE 11 and TE 01 wave modes, in which the transverse electric field 3 of 20 direction is formed parallel to the quartz window surface. In the circular waveguide of the same specification, the cut-off frequency of TE 01 is twice more than that of the TE 11 wave mode. Therefore, it is easier to use the TE 11 mode with a lower cut-off frequency when the frequency is fixed at 2.45 GHz and the waveguide radius is low [1].

Experiment Experiment Setup
Based on the correlations with the plasma characteristics [21], the geometry parameters for the waveguiding structure were set to improve the power transfer efficiency from the wave generator to the plasma region. In the quasi-coaxial line TEM waveguiding structure, microwave power was delivered through a TE/TEM coupler which was bent at 90 • , as shown in Figure 2a.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 20 In the circular waveguide, with the diameter determined, the wave frequency must exceed the cut-off frequency of propagation. Among the lower order modes of circular waveguides, there are the TE11 and TE01 wave modes, in which the transverse electric field direction is formed parallel to the quartz window surface. In the circular waveguide of the same specification, the cut-off frequency of TE01 is twice more than that of the TE11 wave mode. Therefore, it is easier to use the TE11 mode with a lower cut-off frequency when the frequency is fixed at 2.45 GHz and the waveguide radius is low [1].

Experiment Setup
Based on the correlations with the plasma characteristics [21], the geometry parameters for the waveguiding structure were set to improve the power transfer efficiency from the wave generator to the plasma region. In the quasi-coaxial line TEM waveguiding structure, microwave power was delivered through a TE/TEM coupler which was bent at 90°, as shown in Figure 2a. To maximize the power transfer ratio when converting from a rectangular TE10 to a quasi-coaxial line TEM waveguiding structure, a rod-shaped antenna was placed in the field concentration region. The rod-shaped antenna was made of copper as an inner conductor. The diameter of the inner conductor was 9 mm, and the axial length was 600 mm. The microwave propagated mainly in the space between the inner rod and the quartz tube. The space between the inner conductor and the quartz tube was filled with atmospheric pressure for the device cooling by the flow of compressed air [10]. Inside the plasma chamber, quartz window was placed to separate the field generator and plasma excitation region which is filled with argon gas at low pressure. The quartz window = 4 was designed with an inner diameter of 32 mm, a thickness of 2 mm, and a total length of 550 mm in the axial direction. The compressive strength of the quartz window was 117,000 kg/cm 2 , the Young modulus was 73,900 kg/mm 2 , and the tensile strength was 510 kg/cm 2 , which could sustain the mechanical load without deformation in the corresponding axial length. The plasma was ignited by the excitation of a discharge from the electromagnetic waves that penetrate or exist along the quartz window surface.
The propagation wave mode was converted from the rectangular TE10 wave mode to the circular TE11 wave mode as shown in Figure 2b. In the circular TE11 waveguiding structure, microwave power was transferred to the plasma through the continuous cylindrical slot antenna along the wave propagation direction. The circular waveguide was made of copper material with an inner diameter of 30 mm and a thickness of 1 mm. In the cylindrical slotted waveguiding structure, the cut surface of the circular TE11 waveguide was cut to a radial thickness of 10 mm along with an axial length of 500 mm from the entrance of the plasma chamber to the end of the continuous axis. The quartz To maximize the power transfer ratio when converting from a rectangular TE 10 to a quasi-coaxial line TEM waveguiding structure, a rod-shaped antenna was placed in the field concentration region. The rod-shaped antenna was made of copper as an inner conductor. The diameter of the inner conductor was 9 mm, and the axial length was 600 mm. The microwave propagated mainly in the space between the inner rod and the quartz tube. The space between the inner conductor and the quartz tube was filled with atmospheric pressure for the device cooling by the flow of compressed air [10]. Inside the plasma chamber, quartz window was placed to separate the field generator and plasma excitation region which is filled with argon gas at low pressure. The quartz window r = 4 was designed with an inner diameter of 32 mm, a thickness of 2 mm, and a total length of 550 mm in the axial direction. The compressive strength of the quartz window was 117,000 kg/cm 2 , the Young modulus was 73,900 kg/mm 2 , and the tensile strength was 510 kg/cm 2 , which could sustain the mechanical load without deformation in the corresponding axial length. The plasma was ignited by the excitation of a discharge from the electromagnetic waves that penetrate or exist along the quartz window surface.
The propagation wave mode was converted from the rectangular TE 10 wave mode to the circular TE 11 wave mode as shown in Figure 2b. In the circular TE 11 waveguiding structure, microwave power was transferred to the plasma through the continuous cylindrical slot antenna along the wave propagation direction. The circular waveguide was made of copper material with an inner diameter of 30 mm and a thickness of 1 mm. In the cylindrical slotted waveguiding structure, the cut surface of the circular TE 11 waveguide was cut to a radial thickness of 10 mm along with an axial length of 500 mm from the entrance of the plasma chamber to the end of the continuous axis. The quartz window and plasma generation process were set to be the same as in Figure 2a. When the inner diameter of the circular waveguide was 30 mm and was filled with free space material characteristics, the cut-off frequency in the circular waveguide was approximately 5.9 GHz. For the circular waveguide of 30 mm diameter, it was necessary to lower the cut-off fre-Appl. Sci. 2021, 11, 5358 4 of 20 quency below 2.45 GHz for the wave propagation in the applied frequency. Accordingly, a rod type (diameter = 30 mm, r = 9) alumina was inserted inside the circular waveguiding structure, and the cut-off frequency was reduced to a third of 5.9 GHz.
To compare the source characteristics, experiments were conducted under 500 mTorr Ar plasma condition where 0.6 kW microwave power was applied. The experimental system used for this study is shown in Figure 3. Continuous wave (CW) power at 2.45 GHz was transmitted from a magnetron (ASTEX MKS Fl20608-3kW Mag Head, ASTEX Smart Match Fl20606-3kW) to the waveguide mode converter by a rectangular TE 10 wave mode through the WR340 waveguide. The power supply (MKS Smart Power generator Fl20160-1-3kW) was connected to a three-phase line and outputted 300-3000 W. As shown in Figure 3, the microwave transmission line from the generator to the plasma includes a WR340 3-Stub tuner, dummy load (Fl20164), and circulator (AX3161) as the product of National Electronics.
window and plasma generation process were set to be the same as in Figure 2a. When the inner diameter of the circular waveguide was 30 mm and was filled with free space material characteristics, the cut-off frequency in the circular waveguide was approximately 5.9 GHz. For the circular waveguide of 30 mm diameter, it was necessary to lower the cut-off frequency below 2.45 GHz for the wave propagation in the applied frequency. Accordingly, a rod type (diameter = 30 mm, = 9) alumina was inserted inside the circular waveguiding structure, and the cut-off frequency was reduced to a third of 5.9 GHz.
To compare the source characteristics, experiments were conducted under 500mTorr Ar plasma condition where 0.6 kW microwave power was applied. The experimental system used for this study is shown in Figure 3. Continuous wave (CW) power at 2.45 GHz was transmitted from a magnetron (ASTEX MKS Fl20608-3kW Mag Head, ASTEX Smart Match Fl20606-3kW) to the waveguide mode converter by a rectangular TE10 wave mode through the WR340 waveguide. The power supply (MKS Smart Power generator Fl20160-1-3kW) was connected to a three-phase line and outputted 300-3000 W. As shown in Figure 3, the microwave transmission line from the generator to the plasma includes a WR340 3-Stub tuner, dummy load (Fl20164), and circulator (AX3161) as the product of National Electronics. In the circular TE11 LMPS, the microwave transmission line from the power generator to the waveguide mode converter was same as that of the quasi-coaxial waveguiding structure, as shown in Figure 3. In the circular waveguiding structure, WR340 twisted waveguide was used to set the electric field in the required direction because the plasma interacting surface typically faced the deposition region, and the waveguiding structure was changed from a rectangular WR340 waveguide to a circular waveguiding structure through a homemade rectangular to a circular waveguide converter.
For the experiment with the linear MWP source, a cylindrical chamber with a diameter of 500 mm and a height of 150 mm was fabricated, as shown in Figure 4a,b. The microwave power was designed to propagate along the straight line in the quartz tube. To increase the spatial uniformity of the gas flow, a showerhead type gas inlet structure was used, which was located opposite to the deposition region, 40 mm away from the quartz window surface. The images of the plasma generated near the window surface were taken, as shown in Figure 4b. It represents the images of the MWP plasma optical emission patterns along the linear direction from the front side and right side views. The In the circular TE 11 LMPS, the microwave transmission line from the power generator to the waveguide mode converter was same as that of the quasi-coaxial waveguiding structure, as shown in Figure 3. In the circular waveguiding structure, WR340 twisted waveguide was used to set the electric field in the required direction because the plasma interacting surface typically faced the deposition region, and the waveguiding structure was changed from a rectangular WR340 waveguide to a circular waveguiding structure through a homemade rectangular to a circular waveguide converter.
For the experiment with the linear MWP source, a cylindrical chamber with a diameter of 500 mm and a height of 150 mm was fabricated, as shown in Figure 4a,b. The microwave power was designed to propagate along the straight line in the quartz tube. To increase the spatial uniformity of the gas flow, a showerhead type gas inlet structure was used, which was located opposite to the deposition region, 40 mm away from the quartz window surface. The images of the plasma generated near the window surface were taken, as shown in Figure 4b. It represents the images of the MWP plasma optical emission patterns along the linear direction from the front side and right side views. The mounting position of the camera in the plasma chamber is shown in Figure 4b as the front and side views. The equipment used to record the shape of plasma was a cylindrical camera (model: Microsoft Life Cam Studio) suitable for the rounded type chamber flange. The brightness of the emitted light was slightly distorted because the plasma shape was captured by a camera equipped with a neutral density(ND) filter (product of HORUSBENNU, 49 mm), which was used to reduce the brightness to 1/400. mounting position of the camera in the plasma chamber is shown in Figure 4b as the front and side views. The equipment used to record the shape of plasma was a cylindrical camera (model: Microsoft Life Cam Studio) suitable for the rounded type chamber flange. The brightness of the emitted light was slightly distorted because the plasma shape was captured by a camera equipped with a neutral density(ND) filter (product of HORUSBENNU, 49 mm), which was used to reduce the brightness to 1/400.  Figure 5 shows the plasma diagnostics using a homemade single Langmuir probe(LP) in the deposition region. The deposition process is usually carried out 20-25 mm away from the window surface. Hence, the diagnosis was performed in the deposition area at five points (20-25 mm away from the window surface) at 20 mm intervals from the center, as shown in Figure 5. The tip of the LP was made of a 100µm diameter tungsten wire (10 mm in length), which was designed to satisfy the collisionless thin sheath condition [22]. The electron energy distribution function (EEDF) was obtained by calculating the second derivative of the measured I-V curve [23,24]. The second derivative of the measured current with respect to the probe bias voltage, " , was proportional to the electron energy probability function (EEPF) ( ), which was related to the EEDF as ( ) = / ( ). The electron density is calculated by integrating the EEDF, and the electron temperature was measured as an effective electron temperature , which corresponds to a mean electron energy from the integrals of the EEDF [23,24].  Figure 5 shows the plasma diagnostics using a homemade single Langmuir probe(LP) in the deposition region. The deposition process is usually carried out 20-25 mm away from the window surface. Hence, the diagnosis was performed in the deposition area at five points (20-25 mm away from the window surface) at 20 mm intervals from the center, as shown in Figure 5. The tip of the LP was made of a 100 µm diameter tungsten wire (10 mm in length), which was designed to satisfy the collisionless thin sheath condition [22]. mounting position of the camera in the plasma chamber is shown in Figure 4b as the front and side views. The equipment used to record the shape of plasma was a cylindrical camera (model: Microsoft Life Cam Studio) suitable for the rounded type chamber flange. The brightness of the emitted light was slightly distorted because the plasma shape was captured by a camera equipped with a neutral density(ND) filter (product of HORUSBENNU, 49 mm), which was used to reduce the brightness to 1/400.  Figure 5 shows the plasma diagnostics using a homemade single Langmuir probe(LP) in the deposition region. The deposition process is usually carried out 20-25 mm away from the window surface. Hence, the diagnosis was performed in the deposition area at five points (20-25 mm away from the window surface) at 20 mm intervals from the center, as shown in Figure 5. The tip of the LP was made of a 100µm diameter tungsten wire (10 mm in length), which was designed to satisfy the collisionless thin sheath condition [22]. The electron energy distribution function (EEDF) was obtained by calculating the second derivative of the measured I-V curve [23,24]. The second derivative of the measured current with respect to the probe bias voltage, " , was proportional to the electron energy probability function (EEPF) ( ), which was related to the EEDF as ( ) = / ( ). The electron density is calculated by integrating the EEDF, and the electron temperature was measured as an effective electron temperature , which corresponds to a mean electron energy from the integrals of the EEDF [23,24]. The electron energy distribution function (EEDF) was obtained by calculating the second derivative of the measured I-V curve [23,24]. The second derivative of the measured current with respect to the probe bias voltage, I e , was proportional to the electron energy probability function (EEPF) f e (ε), which was related to the EEDF as g e (ε) = ε 1/2 f e (ε). The electron density is calculated by integrating the EEDF, and the electron temperature was measured as an effective electron temperature T e f f , which corresponds to a mean electron energy from the integrals of the EEDF [23,24].

Simulation Methods
Conventionally, analyses of long linear microwave sources reveal limitations to solving the transmission model self-consistently [25][26][27]. Self-consistency requires that the wave field and continuity equations be solved simultaneously to account for the effects of the plasma radial and axial inhomogeneities. The plasma spatial inhomogeneity can be calculated by using numerical methods [28][29][30][31][32][33][34].
Based on numerical methods, in this study, a fluid plasma simulation model(using COMSOL Multiphysics) in the 3-D domain was presented to solve the plasma and wave equations self-consistently [35]. The spatially self-consistent description of the transmis-sion line could be achieved because the plasma fluid simulation could solve the wave and plasma equations simultaneously in 3-D space. The correlation between the plasma generation, longitudinal power absorption, and sustaining the electric field generated from electromagnetic waves could be calculated.
In Figure 6, the simulation geometry was designed in the same way as in the Figure 2 (experiment). The height and the diameter of the chamber were designed the same as the experiment. The parts of the power applicator and structures in the simulation were designed, as shown in the upper parts or panels of the figure. In the case of the quasi-coaxial TEM LMPS, shown in Figure 6a, wave excitation was applied at one side (power feeding port) of the WR340 waveguide structure as the rectangular TE 10 mode. The rectangular TE 10 wave was converted to a quasi-coaxial TEM mode through a rod-shaped antenna at the location 3/4 × λ wr340 (λ wr340 : WR340 guide wavelength, 14.36 cm in 2.45 GHz) away from the power input port where the electric field was maximized. On the other hand, in the circular TE 11 mode, shown in Figure 6b, the WR340 waveguide and the converter of the rectangular to circular waveguiding structure was omitted to reduce the computational loads. As the wave mode was conserved from the rectangular TE 10 to the circular TE 11 [1], the power applicator was designed to be simplified by applying a cylindrical TE 11 wave power directly to one side of the cylindrical alumina surface (power feeding port), as shown in Figure 6b. In the electromagnetic wave model in COMSOL Multiphysics [35], power application method could be selected so that the input power of the active port is fully deposited on the conductive material [36]. In both the cases of this study, this power application model was selected to reduce the calculation time. 500 mTorr Ar plasma with 2.45 GHz and 0.6 kW microwave power was applied. The background gas temperature was set to 300 K, and the neutral gas density was assumed to be uniform and was calculated under the assumption of an ideal gas. The velocity of the electron was assumed to have a Maxwellian distribution. solving the transmission model self-consistently [25][26][27]. Self-consistency requires that the wave field and continuity equations be solved simultaneously to account for the effects of the plasma radial and axial inhomogeneities. The plasma spatial inhomogeneity can be calculated by using numerical methods [28][29][30][31][32][33][34].
Based on numerical methods, in this study, a fluid plasma simulation model(using COMSOL Multiphysics) in the 3-D domain was presented to solve the plasma and wave equations self-consistently [35]. The spatially self-consistent description of the transmission line could be achieved because the plasma fluid simulation could solve the wave and plasma equations simultaneously in 3-D space. The correlation between the plasma generation, longitudinal power absorption, and sustaining the electric field generated from electromagnetic waves could be calculated.
In Figure 6, the simulation geometry was designed in the same way as in the Figure  2 (experiment). The height and the diameter of the chamber were designed the same as the experiment. The parts of the power applicator and structures in the simulation were designed, as shown in the upper parts or panels of the figure. In the case of the quasicoaxial TEM LMPS, shown in Figure 6a, wave excitation was applied at one side (power feeding port) of the WR340 waveguide structure as the rectangular TE10 mode. The rectangular TE10 wave was converted to a quasi-coaxial TEM mode through a rod-shaped antenna at the location 3 4 ⁄ × λ (λ : WR340 guide wavelength, 14.36 cm in 2.45 GHz) away from the power input port where the electric field was maximized. On the other hand, in the circular TE11 mode, shown in Figure 6b, the WR340 waveguide and the converter of the rectangular to circular waveguiding structure was omitted to reduce the computational loads. As the wave mode was conserved from the rectangular TE10 to the circular TE11 [1], the power applicator was designed to be simplified by applying a cylindrical TE11 wave power directly to one side of the cylindrical alumina surface (power feeding port), as shown in Figure 6b. In the electromagnetic wave model in COMSOL Multiphysics [35], power application method could be selected so that the input power of the active port is fully deposited on the conductive material [36]. In both the cases of this study, this power application model was selected to reduce the calculation time. 500 mTorr Ar plasma with 2.45 GHz and 0.6 kW microwave power was applied. The background gas temperature was set to 300 K, and the neutral gas density was assumed to be uniform and was calculated under the assumption of an ideal gas. The velocity of the electron was assumed to have a Maxwellian distribution.  Table 1 lists the Ar chemical reaction set used in the simulation. In this simulation, electrons, Ar + ions, and Ar * (first Ar metastable species) were considered [37][38][39][40]. In  Table 1 lists the Ar chemical reaction set used in the simulation. In this simulation, electrons, Ar + ions, and Ar* (first Ar metastable species) were considered [37][38][39][40]. In particular, the reaction set containing the Ar excited states and step-ionization was selected to increase the accuracy of the calculations in the pressure regime [36][37][38][39][40][41][42]. The cross sections are given for the charged particle reaction sets, and the rate coefficients are given for the reactions of quenching and ionization by atomic collisions. Table 1. Selected reaction for Argon discharges (∆ε: threshold energy, σ c : cross-section, k f : rate coefficient (m 3 /s)) [38][39][40][41].

The Computational Mesh & Geometry
To generate the 3-D meshes based on the advancing-front method (AFM) meshing algorithm, tetrahedral meshes were used for the discretization in 3-D geometry models. The mesh in the quasi-coaxial TEM model consists of approximately 8.7 × 10 5 elements and the mesh in the circular TE 11 model consists of approximately 7.5 × 10 5 domain elements. The meshes used in these models were set finer (0.1~1.2 mm) near the quartz window surface than those in the rest of the discharge chamber (10 mm in maximum).

Governing Equations
The 3-D simulation model in this work coupled two physical models, namely the electromagnetic field and plasma models. The continuous 2.45 GHz electromagnetic wave field (sinusoidal) was considered time-harmonic with amplitude E [34]. The frequency analysis calculated the sinusoidal response for one period and deduced average rms loads. It applied those to the transients in the plasma model with boundary conditions depending on the time t.

The Electromagnetic Field Model
The electric field in the microwave model satisfied the following equation deduced from the Maxwell's equations. The wave equation for this case assumed that the time variation of the electric and magnetic field components was harmonic and the medium was non-magnetic [43].
This Helmholtz equation for the electric field was solved simultaneously in the waveguiding structure and in the plasma chamber. Here, E denotes the electric field(HF field), k 0 is the wave number in free space (k 0 = ω √ µ 0 ε 0 ), ω is the angular frequency (2π× 2.45 GHz), ε 0 and µ 0 are the free space permittivity 8.854 × 10 −12 [F/m] and free space permeability (4π × 10 −7 N/A 2 ), j is imaginary unit, σ is conductivity, respectively. ε r is the relative permittivity depending on the property of the materials. As the material property, the relative permittivity of the quartz window is ε q r = 4. The alumina used in TE 11 mode had a relative permittivity of ε al r = 9, the relative permittivity of copper is ε co r = 1, the relative permittivity of air is ε air r = 1, and the relative permittivity of plasma is expressed as ε p r in the following equations.
where q is the elementary charge and m e is the mass of an electron. ω pe denotes the plasma frequency that depends on the electron density n e , and the v en is the electron-neutral collision frequency for momentum transfer [44].
v en = k tot Ar (ε)·n Ar (4) v en is expressed as a function of mean electron energy ε and the number of neutral particles densities [34,45]. The expression k tot i is the total collision rate of neutral species i. By calculating the mean electron energy (ε) with a Maxwell distribution function and substituting it into the first moment Boltzmann equation, σ p (plasma conductivity) can be obtained as a function of the plasma parameter.
The unmagnetized plasma conductivity σ p is expressed as a function of n e , v en , and ω. The conductivity of copper is 6 × 10 7 S/m. Both quartz and air are set to non-conductive materials (σ q = σ air = 0) and non-magnetic materials (µ q r = µ air r = 1). For non-magnetized high frequency plasmas, the relative permeability µ p r = 1. The density of the abosrbed power from the electromagnetic is calculated as in the following equation [1,[46][47][48][49].
Q h [44] is the density of the absorbed microwave power, Re denotes the real part of the corresponding expression, and E * is the complex conjugate of E.

Plasma Model
Normally, it is known that the gas flow and the temperature highly affect the spatial spreading of neutral gases [43]. However, in this model, the neutral fluid velocity was set to 0 and the neutral gas temperature was set to 300 K to simplify the model in the 3-D domain. Hence, the equations used in this simulation and the terms related to the neutral fluid velocity were omitted. ∂ ∂t (n e ) + ∇·Γ e = R e (7) µ e = q m e ν en , D e = k B T e m e ν en The continuity equation for the electron density is given as above.
Here, E denotes the electric field(DC field), n e is the electron density and R e is electron production rate. In this simulation, the mobility of electrons was calculated as a function of the electron energy by using BOLSIG+ [50], and the diffusivity was determined from Einstein's relation. In this model, the electron mobility used to calculate the electron transport properties was taken from the electron impact reactions set to v en . Here, µ e is the electron mobility, and D e is the electron diffusivity [49][50][51][52][53].

∂ ∂t
(n ε ) + ∇·Γ ε + E·Γ e = S en +Q h /q (10) Appl. Sci. 2021, 11, 5358 9 of 20 The conservation equation for the electron energy density is given by the aforementioned equations, where, n ε is the electron energy density, S en is the energy loss/gain from the collisions, and Q h is the heat source for plasma model absorbed from the electromagnetic wave power. In addition, T e is the electron temperature, µ ε is the electron energy mobility, and D ε is the electron energy diffusivity [51][52][53][54]. Both µ ε and D ε are calculated from the electron mobility µ e , using the Einstein relation.
As aforementioned, for non-electron species (neutral, excited species, and ions), the equation is solved for the mass fraction of each species. This is the simplified form of the Maxwell-Stefan equation, which makes it computationally less expensive in solving the transport equations. Here, w k is the mass fraction of the kth species, j k is the diffusive flux vector, and R k is the rate coefficient for species k. The rate coefficients of electron-heavy particle reactions are calculated from the cross-section of the corresponding reaction σ as a mean value [34]. In addition, ρ is the density of the mixture, v k is the diffusion velocity for species k, D k is the diffusion coefficient (diffusion coefficient of Ar: D Ar is 0.01 m 2 /s, which is same for Ar * and Ar + ), z k is the charge number of species k, and µ k is the mobility for species k [55][56][57][58][59].
As aforementioned, the electrostatic field can be calculated using Poisson's equation. The space charge density, ρ, is computed based on the plasma chemistry (electron density and ion density) [35]. V is the electric potential and n k is the number density of species k.
The boundary conditions of the waveguiding structure power input ports were used to supply electromagnetic waves (input power) in the active ports. At metallic boundaries, the perfect electric conductor (PEC) boundary condition was used, except for the copper material region where the antenna rod was in the quasi-coaxial MWP and the cylindrical guiding structure was in the circular guiding structure MWP. The PEC boundary condition sets the tangential component of the electric field to 0 [35]. The wall boundary condition used to calculate the transport equation was the inhomogeneous Neumann boundary condition. The value of the normal flux to the boundary is described by following equation [34]: where n denotes the unit vector normal to the boundary and v k is the velocity magnitude (thermal velocity for neutrals and Bohm velocity for ions), n k is the density of the k species, and γ k is the sticking coefficient. In this study, the sticking coefficient of each species at the wall boundary was set to 1.
The wall boundary condition [60] in the electron energy equation, which is the homogeneous Neumann condition, were used to set the normal flux of the energy density to 0. The electron energy settings in the metallic material wall boundary conditions were set to 0, the electron reflection coefficient at the wall was set to 0, and the walls in contact with the plasma were set to have a ground potential (V = 0).
In most cases, when calculating the microwave plasma by using the numerical methods, numerical problems arose with regions where n e (electron density) values were close to n cr (critical density). To avoid this, COMSOL Multiphysics provides a method using a "Doppler broadening parameter" which allows for convergence [49]. However, in this model, the use of this method was limited to prevent the deterioration of the computational accuracy.

Results and Discussion
Normally, the generation of optical emission is associated with the Ar* excitation state species, which means electron heating occurs. The electron heating mechanism is based on the electric field profiles, and in this study, the wave field was fed by a wave launcher from the end of the tube. In the linear MWP, which applied power to one side port, electron density attenuation appeared along the linear axis due to the difference in the power absorption along the wave propagation direction, which could result in a specific electron density distribution. Moreover, if the reflection of the waves was made at the opposite boundary, the formation of a standing wave occurred. A standing wave formed in the antenna structure and made inhomogeneous line strengths of an electric field, which caused a spatial difference in the plasma power absorption along the axial direction. Figure 7 shows that the microwave power is applied from the left side in the front view and the right side in the side view. For TE 11 wave mode, as shown in Figure 7a, the brightness of the emitted light is higher at the beginning of wave propagation, and the faint pattern becomes clearer towards the end. From Figure 7b, we see that from the right-side view the optical emission pattern can be observed more clearly. In the quasi-coaxial line TEM wave mode, the same optical emission pattern and the decrease in brightness along the axial directions are observed as in the TE 11 mode. On the other hand, emitted light pattern in the TEM structure was different in the intensity of the emitted light brightness along the axial line and the wavelength of the emitted light wave pattern. In the TEM structure, emitted light brightness concentrated at the beginning of the axial line was relatively lesser and was evenly distributed, as shown in Figure 7c. The wavelength of the emitted light pattern was longer than that of the TE 11 mode. The half wavelength observed on the optical emission pattern was measured to be 23 mm in TE 11 and 47 mm in TEM. set to 0, the electron reflection coefficient at the wall was set to 0, and the walls in contact with the plasma were set to have a ground potential (V = 0).
In most cases, when calculating the microwave plasma by using the numerical methods, numerical problems arose with regions where (electron density) values were close to (critical density). To avoid this, COMSOL Multiphysics provides a method using a "Doppler broadening parameter" which allows for convergence [49]. However, in this model, the use of this method was limited to prevent the deterioration of the computational accuracy.

Results and Discussion
Normally, the generation of optical emission is associated with the Ar* excitation state species, which means electron heating occurs. The electron heating mechanism is based on the electric field profiles, and in this study, the wave field was fed by a wave launcher from the end of the tube. In the linear MWP, which applied power to one side port, electron density attenuation appeared along the linear axis due to the difference in the power absorption along the wave propagation direction, which could result in a specific electron density distribution. Moreover, if the reflection of the waves was made at the opposite boundary, the formation of a standing wave occurred. A standing wave formed in the antenna structure and made inhomogeneous line strengths of an electric field, which caused a spatial difference in the plasma power absorption along the axial direction. Figure 7 shows that the microwave power is applied from the left side in the front view and the right side in the side view. For TE11 wave mode, as shown in Figure 7a, the brightness of the emitted light is higher at the beginning of wave propagation, and the faint pattern becomes clearer towards the end. From Figure 7b, we see that from the rightside view the optical emission pattern can be observed more clearly. In the quasi-coaxial line TEM wave mode, the same optical emission pattern and the decrease in brightness along the axial directions are observed as in the TE11 mode. On the other hand, emitted light pattern in the TEM structure was different in the intensity of the emitted light brightness along the axial line and the wavelength of the emitted light wave pattern. In the TEM structure, emitted light brightness concentrated at the beginning of the axial line was relatively lesser and was evenly distributed, as shown in Figure 7c. The wavelength of the emitted light pattern was longer than that of the TE11 mode. The half wavelength observed on the optical emission pattern was measured to be 23 mm in TE11 and 47 mm in TEM.     -10 show the simulation results of the spatial profiles of the electric field, electron density, and electron temperature, respectively, for each source along the axial and radial directions. Figure 8b shows that the coordinate direction is based on the center of the waveguide, and the X-axis is directed perpendicular (radial direction, θ = 0 • ) from the waveguide to the plasma column; the Y-axis (radial direction, θ = 90 • ) is the direction of the height, and the Z-axis is the wave propagation direction (axial direction) in the waveguiding structure. waveguiding structure.
In Figure 8, the E-field profile (norm | | (E-field amplitude, + + )) is calculated in the 3-D domain. In Figure 8, at TEM wave mode, the direction of the radially emitted wave electric field is perpendicular to the tangential plane of the quartz window surface in the XY plane. In the radial direction, the electric field strength was maximum at the copper antenna surface (inner conductor) and decreased as it was emitted in all directions of plasma-dielectric contact window surfaces [10].   waveguiding structure.
In Figure 8, the E-field profile (norm | | (E-field amplitude, + + )) is calculated in the 3-D domain. In Figure 8, at TEM wave mode, the direction of the radially emitted wave electric field is perpendicular to the tangential plane of the quartz window surface in the XY plane. In the radial direction, the electric field strength was maximum at the copper antenna surface (inner conductor) and decreased as it was emitted in all directions of plasma-dielectric contact window surfaces [10].   On the other hand, in the TE11 mode, the electric field radiated locally in the radial direction from the continuous line slot antenna, and the radially radiated electric field was parallel to the tangential plane of the quartz window surface. In the radial direction, the electric field strength was maximal at the dielectric surface (alumina), which contacted a slot surface of the cylindrical waveguiding structure and it decreased with radial emission.
In the standing waveform, the half-wavelength of the propagating quasi-coaxial line TEM mode through the coaxial line was approximately 48 mm. In the TE11 mode, the half In Figure 8, the E-field profile (norm |E| (E-field amplitude, E 2 x + E 2 y + E 2 z )) is calculated in the 3-D domain. In Figure 8, at TEM wave mode, the direction of the radially emitted wave electric field is perpendicular to the tangential plane of the quartz window surface in the XY plane. In the radial direction, the electric field strength was maximum at the copper antenna surface (inner conductor) and decreased as it was emitted in all directions of plasma-dielectric contact window surfaces [10].
On the other hand, in the TE 11 mode, the electric field radiated locally in the radial direction from the continuous line slot antenna, and the radially radiated electric field was parallel to the tangential plane of the quartz window surface. In the radial direction, the electric field strength was maximal at the dielectric surface (alumina), which contacted a slot surface of the cylindrical waveguiding structure and it decreased with radial emission.
In the standing waveform, the half-wavelength of the propagating quasi-coaxial line TEM mode through the coaxial line was approximately 48 mm. In the TE 11 mode, the half wavelength of the wave propagating through the dielectric was 22 mm, which was approximately 2.1 times smaller. The calculated wavelength by the simulation corresponded to the experimental results, which affected the electron density distribution along the Z-direction.
The electric field attenuation in the axial direction had a considerable influence on the axial distribution of the inhomogeneous plasma. Figure 11 shows the simulation of the attenuation of the normalized E-field amplitude along the wave propagation direction. In general, the rate of attenuation of the wave propagating to the lossy media (X-direction) is expressed as a value proportional to e −αx (α: attenuation coefficient) [1]. On the other hand, in the TE11 mode, the electric field radiated locally in the radial direction from the continuous line slot antenna, and the radially radiated electric field was parallel to the tangential plane of the quartz window surface. In the radial direction, the electric field strength was maximal at the dielectric surface (alumina), which contacted a slot surface of the cylindrical waveguiding structure and it decreased with radial emission.
In the standing waveform, the half-wavelength of the propagating quasi-coaxial line TEM mode through the coaxial line was approximately 48 mm. In the TE11 mode, the half wavelength of the wave propagating through the dielectric was 22 mm, which was approximately 2.1 times smaller. The calculated wavelength by the simulation corresponded to the experimental results, which affected the electron density distribution along the Z-direction.
The electric field attenuation in the axial direction had a considerable influence on the axial distribution of the inhomogeneous plasma. Figure 11 shows the simulation of the attenuation of the normalized E-field amplitude along the wave propagation direction. In general, the rate of attenuation of the wave propagating to the lossy media (X-direction) is expressed as a value proportional to e (α : attenuation coefficient) [1]. The attenuation rate of the electric field intensity (standing wave amplitude) along the wave propagation (Z-direction) could also be represented by the attenuation coefficient α. In the presence of the lossy media, the power absorption was concentrated at the beginning of the applied power [25]. The attenuation coefficient α was affected by the wavelength of the standing wave. In this simulation, the value of the coefficient was Figure 11. Attenuation of electric field norm |E| along Z-direction.
The attenuation rate of the electric field intensity (standing wave amplitude) along the wave propagation (Z-direction) could also be represented by the attenuation coefficient α. In the presence of the lossy media, the power absorption was concentrated at the beginning of the applied power [25]. The attenuation coefficient α was affected by the wavelength of the standing wave. In this simulation, the value of the coefficient was represented by 5 in the quasi-coaxial line TEM and 10 in the circular TE 11 . The value was inversely proportional to the wavelength. In the circular TE 11 wave mode with a short wavelength (which propagated in dielectric space), the power absorption was higher than the quasi-coaxial line TEM at the beginning of wave propagation.
In the linear microwave source, from a macroscopic point of view, two factors affected the spatial distribution of the electron density significantly, namely, the attenuation of the field strength along the wave propagation direction (axial direction) and the coupled field profile at the window surface (radial direction).
First, the electron density distribution along the wave propagation direction was proportional to the attenuation of the electric field strength. As mentioned in the previous section, power absorption mechanism in the plasma was related to the heat source from the electromagnetic wave power. The absorbed power was represented as a function of the variable E. Thus, the Re(J·E) of the power dissipation in the plasma chamber followed the E field profile. In both the cases of LMPS, power dissipation in the axial direction was concentrated at the surface of dielectric near field region.
However, the power dissipation in the XY plane appeared differently depending on the shape of the emitted field direction. In Figures 12 and 13, at a distance of 5 mm, 15 mm, and 25 mm from the quartz window surface (XZ plane at Y = 0), the electron density along the Z-axis is attenuated in proportion to the decrease in the electric field strength, as shown in Figure 11. In the circular TE 11 mode, the attenuation coefficient of the normalized electron density was five times higher than that of the quasi-coaxial line TEM mode, and the electron generation was concentrated at the initial 100 mm of the power application. Near the window surface, the same fluctuations in the electron density appeared as a change in the intensity of the amplitude with a standing wave effect. At a distance of 5 mm, 15 mm, and 25 mm from the quartz window surface along the X-axis, in the case of quasi-coaxial line TEM mode, whose wavelength was longer than that of the circular TE 11 , the tendency of oscillation was higher than that of the circular TE 11 .The short wavelength of circular TE 11 mode was less affected and showed a relatively uniform electron density spatial distribution at 25 mm away from the window surface (with smoothing effect) [46]. mm, and 25 mm from the quartz window surface (XZ plane at Y = 0), the electron density along the Z-axis is attenuated in proportion to the decrease in the electric field strength, as shown in Figure 11. In the circular TE11 mode, the attenuation coefficient of the normalized electron density was five times higher than that of the quasi-coaxial line TEM mode, and the electron generation was concentrated at the initial 100 mm of the power application. Near the window surface, the same fluctuations in the electron density appeared as a change in the intensity of the amplitude with a standing wave effect. At a distance of 5 mm, 15 mm, and 25 mm from the quartz window surface along the X-axis, in the case of quasi-coaxial line TEM mode, whose wavelength was longer than that of the circular TE11, the tendency of oscillation was higher than that of the circular TE11.The short wavelength of circular TE11 mode was less affected and showed a relatively uniform electron density spatial distribution at 25 mm away from the window surface (with smoothing effect) [46].
Although the generation of electrons in the circular TE11 mode was concentrated at the initial 100 mm of wave propagation, the electron density (X = 5, 15, and 25 mm) at Z > 100 mm was 1.5 to 2 times higher than the quasi-coaxial TEM. In the TE11 mode, the electron density spatial distribution was formed close to the planar type. The waveguiding structure and the emitted field direction affected the attenuation of the electron density along the radial distance, the ratio of the electron density of 5-25 mm in the X-axis (at Z = 0) was 20% in the quasi-coaxial TEM mode and 33% in the circular TE11 mode.  When the same power per unit volume was applied to the plasma generation in the simulation, the number of electrons generated in the total volume was 2.8 × 10 15 m −3 in the quasi-coaxial line TEM mode and 2.3 × 10 15 m −3 in the circular TE11 mode. Although the number of electrons produced in the entire volume was 20% higher in the quasi-coaxial Figure 13. Attenuation of electron density along the axial direction (Z-axis) in the XZ plane (circular TE 11 waveguiding structure), distance (5 mm, 15 mm, 25 mm for black, red, and blue lines, respectively) from the window surface along the X-axis, normalized norm |E| on the cut-surface of the cylindrical waveguiding structure (surface of the inner dielectric rod shown by the green line), and the window surface (purple line) are to compare the electron density distribution along the E-fied.
Although the generation of electrons in the circular TE 11 mode was concentrated at the initial 100 mm of wave propagation, the electron density (X = 5, 15, and 25 mm) at Z > 100 mm was 1.5 to 2 times higher than the quasi-coaxial TEM. In the TE 11 mode, the electron density spatial distribution was formed close to the planar type. The waveguiding structure and the emitted field direction affected the attenuation of the electron density along the radial distance, the ratio of the electron density of 5-25 mm in the X-axis (at Z = 0) was 20% in the quasi-coaxial TEM mode and 33% in the circular TE 11 mode.
When the same power per unit volume was applied to the plasma generation in the simulation, the number of electrons generated in the total volume was 2.8 × 10 15 m −3 in the quasi-coaxial line TEM mode and 2.3 × 10 15 m −3 in the circular TE 11 mode. Although the number of electrons produced in the entire volume was 20% higher in the quasi-coaxial line TEM mode than that in the circular TE 11 , the number of electrons generated in the actual deposition region was 50% higher in the circular TE 11 mode than that in the quasi-coaxial TEM mode.
In Figure 10, the spatial profiles of the electron temperature for each source are represented in the axial and radial directions. For the axial direction (XZ plane, Y = 0), the electron temperature was nearly constant along the axial distance. For radial distances in the XZ plane (Y = 0), near the window surface, the electron temperature was relatively high with a value of 3.5 eV while in the case of radial distances > 10 mm from the quartz window surface, the electron temperature with the value of 0.8~1.3 eV weakly depended on the radial distance. In Figure 10, the electron temperature distribution of the XY plane in each source varied depending on the shape of the emitted field profile such as the electron density distribution. The electron density and the electron temperature near the substrate surface (z = 20 and 25 mm) measured by a single LP are compared with the simulation results in Figures 14 and 15. The results of the experimental diagnosis show that the tendency is consistent with the simulation results. The experimental results in Figure 14a shows that the tendency of the electron density distribution along the Z-axis is consistent with the simulation result in Figure 14b. The electron density decreased along the wave propagation direction, and the attenuation ratio was higher in the circular TE 11 mode. The electron density fluctuations due to the standing wave effect were prominent in the quasi-coaxial TEM mode, and the effect decreased along the X-axis. When comparing the electron density according to the wave mode, in both cases, the electron density near the substrate was measured to be approximately 1.5~2 times higher in the circular TE 11 mode than that in the quasi-coaxial line TEM mode. In the experiment, the rate of increase in density was 2, which was more than the simulation value of 1.5. This is because the electron loss from the quartz window surface along the direction of the electric field was not considered in detail in this simulation model. because the electron loss from the quartz window surface along the direction of the electric field was not considered in detail in this simulation model. In the experiment, the energy distribution function of electrons of radial distance ≥ 20 mm followed the Maxwell distribution function. The electron temperature profile shown in Figure 15a is consistent with the simulation results in Figure 15b, which shows that when the radial distances > 20 mm from the quartz window surface, the electron temperature is weakly dependent on the radial distance.
To consider each source and the plasma characteristics more deeply, it was necessary to consider the effects by resolving the electric field into the radial and axial direction electric field components. The material properties of plasma act like a lossy medium, which cause not only the transverse field components but also the axial field components to exist [16,17]. The presence of an axial component of the wave electric field in the quasicoaxial line source was detailed by H. Nowakowska et al. [61]. In the presence of the axial component of the wave field, in the high electron density region (ne > 10 18 ), internal wave was represented by a wave that propagated near the window surfaces along the axial direction. Power dissipation in the plasma chamber is concentrated near the dielectric window surface region, so the presence of an internal wave is important as a component that affects the electron heating in the axial direction.   In Figure 12, a normalized norm | | of the inner conductor (green line) and the window surface (purple line) in the quasi-coaxial TEM source are set to compare with the electron density distribution along the axial direction (Z-axis) in the XZ plane. Normalized norm | | was expressed as the RSS (Root Sum Square: ∑ ) of the axial and radial In the experiment, the energy distribution function of electrons of radial distance ≥20 mm followed the Maxwell distribution function. The electron temperature profile shown in Figure 15a is consistent with the simulation results in Figure 15b, which shows that when the radial distances >20 mm from the quartz window surface, the electron temperature is weakly dependent on the radial distance.
To consider each source and the plasma characteristics more deeply, it was necessary to consider the effects by resolving the electric field into the radial and axial direction electric field components. The material properties of plasma act like a lossy medium, which cause not only the transverse field components but also the axial field components to exist [16,17]. The presence of an axial component of the wave electric field in the quasi-coaxial line source was detailed by H. Nowakowska et al. [61]. In the presence of the axial component of the wave field, in the high electron density region (ne > 10 18 ), internal wave was represented by a wave that propagated near the window surfaces along the axial direction. Power dissipation in the plasma chamber is concentrated near the dielectric window surface region, so the presence of an internal wave is important as a component that affects the electron heating in the axial direction.
In Figure 12, a normalized norm |E| of the inner conductor (green line) and the window surface (purple line) in the quasi-coaxial TEM source are set to compare with the electron density distribution along the axial direction (Z-axis) in the XZ plane. Normalized norm |E| was expressed as the RSS (Root Sum Square: i ) of the axial and radial components of the electric field. Figure 12 shows that the electron density distribution along the axial direction in the XZ plane follows the norm |E| profile in the window surface rather than the profile of the norm |E| in the inner conductor surface. This result is due to the axial field formed near the window surface region, and it can be observed in more detail in Figure 16. Figure 16 shows the distribution of the power dissipation, and the electron density distributions according to the electric field components in the axial and radial directions.
conductive material, including the inner copper antenna, was less than 0.1 W. Therefore, when the total power of 600 W was applied to generate the plasma, power dissipation by the axial field was doubled compared to the radial field. Thus, in the quasi-coaxial TEM MWP sources, it was observed that axial field components in a horizontal form on the window surface compared to radial field components in a perpendicular form on the window surface (electron losses occur at the wall boundary) dominantly affected the electron heating in the plasma discharge. Figure 16. Electron density and power dissipation Q h,z (=Re ( 1 2 σ·E z ·E * z )) distributions according to the |E z | profile, electron density and power dissipation Q h,x (=Re ( 1 2 σ·E x ·E * x )) distributions according to the |E x | profile in the XZ plane of quasi co-axial line TEM waveguiding structure. (a) Electron density distribution according to the profile of |E z |, (b) the distribution of power dissipation value Q h,z according to the profile of |E z |, (c) electron density distribution according to the profile of |E x |, and (d) the distribution of power dissipation value Q h,x according to the |E x | profile.
In the TEM wave mode, because the radially radiated electric field was directed in all directions with the same value, E x (= E y ) represented the component of the radial field, and E z represented the component of the axial field. In the case of the axial electric field E z , the field strength is the highest at 4 × 10 3 V/m for the window surface, and for the radial electric field E x , the field strength is the highest at 3 × 10 4 V/m for the inner conductor surface. In Figure 16a,c, it is possible to compare the spatial distribution of the electron density according to each electric field component. Similar to Figure 12, it can be seen that the electron density distribution on the XZ plane follows the axial field component and not the radial field component. In Figure 16b,d, the spatial distribution characteristics of the power dissipation deposited in the plasma can be seen along each component of the electric field. In the figures, the maximum strength of the radial field is 3 × 10 4 V/m and has a value greater than 4 × 10 3 V/m of the axial field, while the maximum value of the power dissipation along the radial field near the surface region is 1 × 10 5 W/m 3 , which is lower than 4 × 10 5 W/m 3 of the axial field. When the power dissipation value deposited in the plasma discharge in the plasma chamber was volume-integrated according to the radial and axial fields, in the quasi-coaxial TEM MWP, the power dissipation value by the axial field ≈ 400 W, and the power consumption by the radial field ≈ 200 W. In this simulation model, the power application method that inputted 600 W was fully deposited to the conductive material, and the plasma was applied. The power consumed by the conductive material, including the inner copper antenna, was less than 0.1 W. Therefore, when the total power of 600 W was applied to generate the plasma, power dissipation by the axial field was doubled compared to the radial field. Thus, in the quasi-coaxial TEM MWP sources, it was observed that axial field components in a horizontal form on the window surface compared to radial field components in a perpendicular form on the window surface (electron losses occur at the wall boundary) dominantly affected the electron heating in the plasma discharge.
In Figure 13, normalized norm |E| on the cut surface of the cylindrical waveguiding structure (green line) and window surface (purple line) in the circular TE 11 wave MWP source are set to compare with the electron density distribution along the axial direction (Z-axis) in the XZ plane. Unlike the result of Figure 12, Figure 13 shows that electron density distribution along the axial direction in the XZ plane simultaneously follows the norm |E| profile in the cut surface of the cylindrical waveguiding structure and also in the window surface. This was because the axial field had a small effect on the electron generation in the spatial distribution. It represented that electron generation was mainly affected by the distribution of the radial direction the electric field component. When the volume-integrated power dissipation value deposited on the plasma discharge in the plasma chamber was expressed as components of the radial and axial fields, in the circular TE 11 wave MWP, the power dissipation value by the axial field ≈ 40 W, and the power dissipation by the radial field ≈ 560 W. In this numerical model, the power application method in the circular TE 11 source was applied in the same way as in the quasi-coaxial TEM source, and the power consumed by the conductive material was also less than 0.1 W. Therefore, it can be seen that when a total of 600 W was applied to generate the plasma, the dissipation of power by the radial field was more dominant than that of the axial field component.
In Figure 17, the spatial distribution of the power dissipation and electron density distributions are represented according to the electric field components, which are indicated by resolving the radial fields into the E x and E y components. In the case of E x , the field strength was the highest for 6 × 10 3 V/m near the window surface, and the field direction was perpendicular to the window surface. Similarly, for the case of E y , the field strength was the highest for 1.5 × 10 4 V/m, and the field direction was parallel to the window surface. In both these cases, as shown in Figure 17, the spatial distributions of the electron generation and power dissipation coincided with the electric field concentrated region. When the volume-integrated power dissipation value deposited for plasma generation was calculated by resolving into the E x and E y field components, the power dissipation value E x ≈ 300 W and E y ≈ 260 W. Among the radial field components involved in the plasma discharge, a dominant component was E x , which was due to the field emitted geometry characteristics of the waveguiding structure.

Conclusions
A conventional quasi co-axial TEM waveguiding structure and a newly developed circular TE11 waveguiding structure sources are presented to compare the plasma characteristics along the waveguiding structure. To compare the source characteristics, experiments were conducted under 500mTorr Ar plasma condition, where 0.6 kW microwave power was applied. The characteristics of the plasma was diagnosed by a Figure 17. Electron density and power dissipation Q h,x (=Re ( 1 2 σ·E x ·E * x )) distributions according to the |E x | profile, electron density and power dissipation Q h,y (=Re ( 1 2 σ·E y ·E * y )) distributions according to the E y profile in the XZ plane of circular TE 11 waveguiding structure. (a) Electron density distribution according to the profile of |E x |, (b) the distribution of power dissipation value Q h,x according to the profile of |E x |, (c) electron density distribution according to the profile of E y , and (d) the distribution of power dissipation value Q h,y according to the E y profile.

Conclusions
A conventional quasi co-axial TEM waveguiding structure and a newly developed circular TE 11 waveguiding structure sources are presented to compare the plasma characteristics along the waveguiding structure. To compare the source characteristics, experiments were conducted under 500mTorr Ar plasma condition, where 0.6 kW microwave power was applied. The characteristics of the plasma was diagnosed by a single LP and analyzed by the numerical 3-D model using plasma fluid simulation. Attenuation of the electric field strength along the wave propagation direction and the coupled field profile at the window surface affected the spatial distribution of the electron density significantly. In the TE 11 mode, the attenuation coefficient of the normalized electron density was 4.8 times higher than that of the TEM mode. The field direction affected the attenuation of the electron density along the radial distance, and the electron-density reduction rate along the radial direction was higher in the TEM mode than that in the TE 11 mode. The plasma characteristics of each plasma source were considered by resolving the wave electric field into the radial and axial field components. The influence of the axial electric field component was more dominant in the quasi-coaxial TEM MWP source. The number of electrons produced in the entire volume was 20% higher in the TEM mode than that in the TE 11 mode, but the number of electrons generated in the actual deposition region was 50% higher in the TE 11 mode than that in the TEM mode.