Magnetic nozzle radiofrequency plasma thruster approaching twenty percent thruster efficiency

Development of a magnetic nozzle radiofrequency (rf) plasma thruster has been one of challenging topics in space electric propulsion technologies. The thruster typically consists of an rf plasma source and a magnetic nozzle, where the plasma produced inside the source is transported along the magnetic field and expands in the magnetic nozzle. An imparted thrust is significantly affected by the rf power coupling for the plasma production, the plasma transport, the plasma loss to the wall, and the plasma acceleration process in the magnetic nozzle. The rf power transfer efficiency and the imparted thrust are assessed for two types of rf antennas exciting azimuthal mode number of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=+1$$\end{document}m=+1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=0$$\end{document}m=0, where propellant argon gas is introduced from the upstream of the thruster source tube. The rf power transfer efficiency and the density measured at the radial center for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=+1$$\end{document}m=+1 mode antenna are higher than those for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=0$$\end{document}m=0 mode antenna, while a larger thrust is obtained for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=0$$\end{document}m=0 mode antenna. Two-dimensional plume characterization suggests that the lowered performance for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=+1$$\end{document}m=+1 mode case is due to the plasma production at the radial center, where contribution on a thrust exerted to the magnetic nozzle is weak due to the absence of the radial magnetic field. Subsequently, the configuration is modified so as to introduce the propellant gas near the thruster exit for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=0$$\end{document}m=0 mode configuration and the thruster efficiency approaching twenty percent is successfully obtained, being highest to date in the kW-class magnetic nozzle rf plasma thrusters.

tens of W to several kW; the electron thermal energy is often converted into the directed ion energy via electrostatic ion accelerations in a current-free double layer and an ambipolar electric field 17 , where the accelerated ions are neutralized by electrons overcoming the potential drop 18 . Fundamental studies on such an electronheated magnetic nozzle plasma thruster have shown that an internal azimuthal plasma current induces an axial Lorentz force in the magnetic nozzle, resulting in an increase in the thrust [19][20][21] . The azimuthal plasma current has been identified to be mainly driven by an electron diamagnetic drift, which is induced by a radial pressure gradient 22,23 . Therefore, the radial electron pressure is converted into the axial plasma momentum flux in the magnetic nozzle, where the electron temperature is decreased along the axis by losing their internal energy [24][25][26][27][28] . A two-dimensional magnetic nozzle thruster model in Ref. 29 can quantitatively explain the measured thrust and can be approximately rewritten by a one-dimensional model using a paraxial approximation, being similar to a physical nozzle model 30 . However, the thrust estimated by the paraxial approximation is a few tens of percent smaller than that by the two-dimensional model. It is expected that the two-dimensional structure of the plasma flow affects the thrust in the magnetic nozzle.
Fundamental laboratory experiments on the rf plasma sources have shown that a peripheral high density region is formed in an inductive mode, while a central density peak is formed simultaneously with an appearance of a helicon wave 2,31 . The excitation of the helicon wave often produces a high density plasma downstream of the source due to the wave propagating in the plasma column 32 . Several experiments in the expanding magnetic fields have shown that the high density conics are formed in the magnetic nozzles [33][34][35][36] . Two-dimensional measurement of an electron energy distribution has indicated that the high temperature electrons heated by the rf antenna field are transported along the expanding magnetic field 37 . Therefore, the structure of the electromagnetic fields would play an important role in the formation of the two-dimensional structures, which would affect the imparted thrust as discussed above. Numerical studies have implied that the rf power is absorbed in the magnetic nozzle region downstream of the source due to the wave propagation 38 as discussed in the experiment 39 . As the wave excitation and the rf power coupling are affected by the antenna structure, an azimuthal mode number m of the excitation mode, and the boundary conditions at the walls [40][41][42] , the optimization of the antenna structure for the thruster configuration is an important experimental issue to be further investigated. The above-mentioned experiments for the plasma source have often been performed with the m = 0 and +1 mode antennas. Investigations on the antenna effect on the thruster performance are required for the performance improvement since the thruster performance is not simply determined by the local plasma density, while different antennas have been tested for the different source geometries and magnetic fields, e.g., in Refs. [12][13][14][15] . Therefore, testing different antennas in the unchanged source geometry will provide an insight on the antenna optimization for the thruster performance.
In addition to the structure of the electromagnetic fields, analytical, numerical, and experimental studies have shown that a neutral density profile affects both the spatial and temporal behaviors of the plasma density, especially in highly ionized plasmas [43][44][45][46][47] . When the neutral gas is introduced from the upstream side of the source, the neutral density is depleted near the thruster exit due to the high ionization rate and the plasma density is lowered there as predicted in an analytical model 48 and observed in the experiment 49 . In such a density profile having the axial density peak in the upstream side, a non-negligible axial momentum flux is transferred to the radial wall, corresponding to the loss of the thrust to the wall 50,51 . One of the experiments has shown that the thrust loss can be inhibited by injecting the gas near the thruster exit and that the thrust can be increased by introducing the gas near the thruster exit 52 , in addition to the inhibition of the loss by the magnetic field 53 .
Several experiments have been performed to improve the performance of the magnetic nozzle rf plasma thruster, where the first direct thrust measurement showed the thruster efficiency less than a percent 54,55 . The subsequent experiment has shown that the thrust due to the magnetic nozzle can be increased by increasing the magnetic field, where the thrust approaches the theoretical limit assuming no plasma loss from the magnetic nozzle 22 . The upper limit of the thrust can be increased by enlarging the diameter of the source tube as predicted by combining the magnetic nozzle model with the global plasma production model 56 . The previously reported thruster efficiency calculated with the rf power is at most ∼ 10% for an rf power level of several kW 56 .
Here the investigation on the effect of the antenna structure is initiated from the comparison between the often-used two modes of m = 0 and m = +1 . In the present experiment, two different antennas exciting m = +1 and m = 0 modes shown in Fig. 1 are tested for the unchanged source tube in the magnetic nozzle rf plasma thruster which has the upstream gas injection port and is operated at the rf power of several kW, where the rf power transfer efficiency, the plasma density at the radial center, and the imparted thrust are compared between the two cases. A better power coupling and a higher plasma density at the radial center near the thruster exit are obtained for the m = +1 mode antenna, while an imparted thrust for the m = 0 mode antenna is larger than that for the m = +1 mode antenna. The configuration is subsequently modified so as to introduce the argon gas near the thruster exit for the m = 0 mode; the thruster efficiency estimated from the thrust, the mass flow rate, and the rf power, approaches twenty percent, being the highest efficiency to date in the magnetic nozzle rf plasma thruster operated with the rf power of the several kW level. The highest performance is resultantly originated from combining the insights of the effectiveness of the m = 0 mode antenna, the large diameter source tube, the sufficient magnetic field strength, and the gas injection near the thruster exit.

Experimental setup
Experiments are performed in a 1-m-diameter and 2-m-long cylindrical vacuum chamber, details of which are described in "Method" section. A pendulum thrust balance is installed inside the vacuum chamber and a thruster structure shown in Fig. 2a is attached to the balance. The thruster consists of a 9.5-cm-inner-diameter, 11-cm-outer-diameter, and 20-cm-long ceramic source tube and a solenoid located near the open source exit. The upstream side of the source tube is terminated by an insulator plate having a small center hole. Two  www.nature.com/scientificreports/ in Fig. 2, where the direction of the magnetic field at the radial center is indicated by a block arrow in Fig. 2a.
Either the m = +1 mode helical antenna (Fig. 1a) or the m = 0 mode loop antenna ( Fig. 1b) is wound around the source tube. The antenna is powered by a 13.56 MHz rf generator via an impedance matching box located outside of the chamber; the introduced argon gas is ionized, and the plasma is produced inside the source tube. Two capacitors in the matching box are tuned to minimize the reflected power and the power reflection is undetectable (less than 1 W) for all the conditions in the present experiment. The rf power transfer efficiency η p and the local plasma density n p are estimated by measuring the antenna current and by using a Langmuir probe. The equilibrium position of the pendulum balance is displaced when the plasma imparts a force to the thruster; the displacement induced only by the plasma production is measured by using a precise laser displacement sensor (with a resolution of 0.1 µm ). The absolute value of the force corresponding to the thrust F can be obtained by multiplying a calibration coefficient relating the displacement to the force. More detailed setup and experimental procedures can be found in "Method" section. Figure 3 shows the measured rf power transfer efficiency η p for various values of the solenoid current I B as a function of the rf power P rf , where the data in Fig. 3a,b are taken with the m = +1 and m = 0 mode antennas, respectively, and the upstream gas injector is used. As clearly seen in Fig. 3a, the power transfer efficiency η p discontinuously increases when increasing the power; the discontinuous change is considered to be the mode change from capacitive-to inductive-or wave-coupled modes. The rf power threshold of the transition from the low to high efficiencies is found to increase with an increase in the magnetic field strength. The similar feature has been observed in the helicon source experiments, e.g., in Ref. 57 . After the transition to the efficient coupling mode, which is probably the helicon wave discharge mode, the power transfer efficiency greater than 0.9 is obtained for all the magnetic field strength. On the other hand, the transition of the discharge mode is not observed, and the power transfer efficiency is nearly constant at about 0.85-0.9 for the m = 0 mode case as seen in Fig. 3b, which is lower than that for the m = +1 mode case.

Results
The plasma density n p as a function of the rf power P rf for the upstream gas injection is measured at the open source exit and the radial center, i.e., (z, r) = (0, 0) , for the m = +1 and m = 0 mode cases, as plotted in Fig. 4a,b, respectively. Simultaneously with the discontinuous change in the power transfer efficiency (see Fig. 3a), the density jump by an order of magnitude is detected as shown in Fig. 4a for the m = +1 mode case, where the maximum density is about 4 × 10 17 cm −3 . For the m = 0 mode case, no discontinuous change in the density is detected as well as the rf power transfer efficiency, where the maximum density at (z, r) = (0, 0) is about 3 × 10 17 cm −3 and slightly lower than that for the m = +1 mode.
For the same conditions as in Figs. 3 and 4, the thrust is measured by using the pendulum thrust balance for the m = +1 and m = 0 mode cases as shown in Fig. 5a,b, respectively, where the thrust for the low power transfer efficiency conditions in the m = +1 mode is very small due to the poor power coupling (Fig. 3) and the low plasma density (Fig. 4); hence only the data with the high power transfer efficiency mode are taken as in Fig. 5a. Both the data in Fig. 5a,b show that the thrust increases with an increase in the rf power P rf . As already described, the better power transfer efficiency is obtained for the m = +1 mode antenna; nevertheless the thrust for the m = 0 mode case (Fig. 5b) is significantly larger than that for the m = +1 mode case (Fig. 5a).
According to the simple two-dimensional model 22 , the thrust F imparted by the magnetic nozzle rf plasma thruster is given by  Typical error bar is about ± 5% . The discontinuous change in η p , which would be a discharge mode transition, is observed for the m = +1 case and the efficiency of η p > 0.9 is obtained for the efficient coupling conditions, while the efficiency for the m = 0 case is less than 0.9 with no mode transition. www.nature.com/scientificreports/ where p e0 , B r , B z , and p e are the maximum electron pressure in the source, the radial magnetic field, the axial magnetic field, and the electron pressure, respectively. In Eq. (1), the axial momentum flux lost to the radial source wall is neglected for simplicity, and the axisymmetric system is assumed. The first and second terms in the right-hand side (RHS) of Eq. (1) are the electron pressure force exerted to the back plate, and the Lorentz force arising from the azimuthal electron diamagnetic drift current and the radial magnetic field, respectively. Equation (1) indicates that the thrust could be affected by the two-dimensional profile of the electron pressure. Previous experiments have shown that the electron temperature is not uniform in the magnetically expanding rf plasmas; the precise identification of the electron pressure profile requires both the density and the electron temperature profiles as measured before 37 ; the full sweep of the I-V characteristics of the Langmuir probe have www.nature.com/scientificreports/ to be taken in the r-z plane. However it can be very qualitatively characterized by the ion saturation current I is of the Langmuir probe, since I is and p e are proportional to n p √ T e and n p T e , respectively. Figure 6 shows the two-dimensional profiles of the ion saturation current I is of the Langmuir probe taken at the rf power of P rf = 4 kW for the (a) m = +1 and (b) m = 0 mode cases. The contour lines in the magnetic nozzle region shows that the density or pressure decay along the axis for the m = +1 mode case (Fig. 6a) is conspicuous, compared with that for the m = 0 mode case, which qualitatively implies that the second term in Eq. (1) for the m = 0 mode case is larger than that for the m = +1 mode case.
To discuss in more detail, radial profiles of the ion saturation current I is at representative axial positions are shown in Fig. 7. For both the m = +1 and m = 0 mode cases, the peaks of I is are observed at the radially peripheral region of the plasma column, e.g., around r ∼ ±2.5 cm at z = 2 cm and r ∼ 6 cm at z = 10 cm as seen in Fig. 7a,b, respectively. Previous experiments have already shown that the peripheral high density and/ or electron temperature region is originated from the electrons heated by the rf antenna and transported along the expanding magnetic field lines 37 . Recent spatial measurement of the plasma momentum flux has also shown that the peripheral high density and temperature region significantly contributes to the thrust generation in the magnetic nozzle 58 . Furthermore, the plasma radius for the m = 0 mode case seems to be slightly larger than that for the m = +1 mode case, which would contribute to the increase in the thrust by the magnetic nozzle according to Eq. (1). For the m = +1 mode case, another peak of I is at the radial center can be observed as in Fig. 7b, which is probably due to the excitation of the helicon wave and the resultant ionization at the radial center as observed by Degeling et al. before 31 . However, the contribution of the high density plasma at the central region on the magnetic nozzle thrust term [the second term of the RHS in Eq. (1)] would be small due to the absence of B r at the radial center. Moreover, I is for the m = 0 mode at z = 20 cm is about 30% larger than that for the  Figure 6. Two-dimensional profiles of the ion saturation current I is of the Langmuir probe taken at the rf power of P rf = 4 kW for the (a) m = +1 and (b) m = 0 mode cases. For both the cases, the plasma expands in the magnetic nozzle and the density decay along the axis can be seen, where the density decay along the axis for the m = +1 mode case seems to be conspicuous, compared with the m = 0 mode case. www.nature.com/scientificreports/ m = +1 mode case as seen in Fig. 7c, despite the similar values of I is for the two antenna cases at z = 10 cm (see Fig. 7b). Although the mechanisms of the larger value of I is in Fig. 7c for the m = 0 mode case is still unclear, a possible reason would be a cross-field transport in the magnetic nozzle, which would be significantly affected by the pressure gradient and plasma instabilities. Summarizing the results in Figs. 3, 4, 5, 6 and 7, the large thrust can be obtained for the m = 0 mode antenna despite the poorer rf power transfer efficiency, compared with the m = +1 mode case, due to the larger contribution of the Lorentz force arising from the electron diamagnetic drift current and the radial magnetic field. The previous experiment has shown that the non-negligible thrust is lost to the radial wall when the propellant gas injected from the upstream side is highly ionized and the neutrals are depleted near the thruster exit 50 . Another experiment using a smaller diameter source tube than the present experiment has demonstrated that the thrust loss is inhibited by injecting the gas near the thruster exit 52 . For the m = 0 mode antenna configuration, which gives the better performance than the m = +1 mode antenna, the argon gas is introduced near the thruster exit by using the downstream gas injector as seen in Fig. 2a, where the gas flow rate is maintained at 70 sccm. Figure 8a shows the measured thrust as a function of the rf power P rf for the solenoid current of I B = 15 A and the large thrust is clearly obtained compared with the data in Fig. 5b, e.g., the thrust for the rf power of P rf ∼ 4 kW and the solenoid current of I B = 15 A is ∼ 40 mN and ∼ 55 mN in Figs. 5b and 8a, respectively.
The thruster efficiency η defined by  Fig. 6. The profile for the m = +1 mode antenna in (b) shows the center peak of I is in addition to the peripheral peak. The measured profiles imply that the plasma radius for the m = 0 mode case is slightly larger than that for the m = +1 mode case. Furthermore, the values of I is for the m = 0 mode case at z = 20 cm is about 30% larger than that for the m = +1 mode case as seen in (c). www.nature.com/scientificreports/ is used to assess the thruster performance, where F, ṁ , P rf are the thrust, the mass flow rate of the propellant, and the rf power. It should be noted that the rf power P rf is the output power from the rf generator; hence the dc-rf conversion efficiency in the generator and the electric power for the solenoid are not taken into account here. The thruster efficiency η calculated from the data in Fig. 8a is plotted in Fig. 8b; the thruster efficiency increases with an increase in the rf power and approaches about twenty percent, which is the highest values to data as far as the author knows. The presently described performance improvement also contributes to an active space debris removal technology by using the magnetic nozzle plasma thruster 59 . It should be noted that the power consumed for the solenoid is about 900 W, being 15% of the maximum rf power of 6 kW. Therefore, the power consumption for the solenoid would effectively lower the thruster efficiency by a factor of 1/1.15 ∼ 0.87 . The influence of the power loss at the solenoid would be lowered when operating the thruster at higher rf power level. Furthermore, replacement of the solenoid by permanent magnets suggested previously 60-62 will eliminate the power consumption by the solenoid, which is a further challenge for applying the sufficient magnetic field strength to the large diameter source tube.

Conclusion
Performance of the magnetic nozzle rf plasma thruster is assessed for the two different rf antennas exciting the m = +1 and m = 0 mode electromagnetic fields, where the argon gas is introduced from the upstream side of the source tube. The larger thrust is obtained for the m = 0 mode antenna configuration due to the larger contribution of the diamagnetic thrust in the magnetic nozzle. For the m = 0 mode antenna configuration, the thruster efficiency approaching 20% is obtained by injecting the propellant gas near the thruster exit. The estimated thruster efficiency is highest to date in the magnetic nozzle rf plasma thrusters operated at the rf power of the several kW level.

Methods
Vacuum chamber and thruster settings. Figure 9 shows the schematic diagram of the experimental setup. The 1-m-diameter and 2-m-long vacuum chamber is evacuated by three 3000 Ls −1 turbomolecular pumps via gate valves to a base pressure less than 10 −4 Pa, where the effective total pumping speed for argon is about 4500 Ls −1 . The pendulum thrust balance is installed inside the vacuum chamber 63 , where a bottom aluminum plate is suspended from a top plate attached to the chamber by axially flexible metallic plates, allowing the bottom plate to be displaced by an axial force. The thruster structure including the source tube, the solenoid, and the back plate, is attached to the bottom plate. The rf antenna is wound around the source tube with no mechanical contact to ensure the pendulum motion of the balance. The antenna, which is made of a water-cooled copper tube, is covered by a thick insulator structure and is further shielded by a grounded metallic structure to suppress parasitic discharges outside the source tube 64 for both the m = +1 and m = 0 mode antennas. Argon gas is continuously introduced from either the upstream or downstream gas injector via a mass flow controller (with an accuracy of ±1 sccm), and the gas pressure is measured by using an ionization gauge connected to the chamber side port at z ∼ 60 cm, where the calibration coefficient for argon is used to obtain the absolute value of the pressure. The measured gas pressure is unchanged by the gas injection port and is about 25 mPa for 70 sccm argon gas flow rate. After introducing the gas, the DC solenoid current is turned on; further the rf power is fed to the antenna via a vacuum feedthrough and the impedance matching box; the plasma is produced inside the source tube and expands into the vacuum chamber. It is well known that the pressure inside the vacuum chamber and the resultant plasma behavior are significantly changed by the pumping speed. Previous experiment with the same vacuum chamber as the present www.nature.com/scientificreports/ experiment has tested the thrust measurement for different pumping speeds by closing two of three gate valves. As reported in 65 , no change in the thrust by the pumping speed has been detected for the gas flow rate similar to the present experiment.
Thrust measurement. The axial displacement induced by the plasma production, i.e., the thrust, is measured by a laser displacement sensor (with a resolution of 0.1 µm ). For the thrust assessment, the gas is continuously introduced to the thruster the chamber pressure reaches the steady state. The DC solenoid current is supplied; then it is observed that the equilibrium position of the pendulum is moved due to a magnetic force. At about 5-10 s after turning on the solenoid current, the rf power is turned on for about 10 s; then the equilibrium position is further moved due to the force exerted to the thruster by the plasma. From the difference in the equilibrium position with only the solenoid current and with both the solenoid current and the rf power, the displacement induced by the plasma can be obtained. It is noted that no displacement is detected for no gas flow rate (i.e., no plasma) even when the rf power is turned on. The absolute value of the force is obtained by multiplying the calibration coefficient relating the displacement to the force.

Thrust balance calibration.
After installing all the thruster components (including the source tube, the solenoid, the power cable for the solenoid, the gas injector, the antenna) and before pumping the vacuum chamber, a calibration procedure is performed by applying known axial forces to the balance and measuring the displacements. To apply the known axial force, known mass pieces are put on a calibration basket attached to a horizontal thread connected to the thruster, where a second thread is further attached to the support mounted on the chamber. When the angle θ , the horizontal distance L, and the height H, are given as indicated in Fig. 9, the axial force F cal exerted to the thruster structure can be calculated from the force balance in a simple mechanics as The distance L and the height H are measured by a straightedge with the error of 0.5% , and the mass of the piece measured by an electronic mass balance (with a resolution of 0.001 g) is 0.3234 ± 0.005 g. The error σ F in the force can be estimated as 1.7% from the propagation law of the error given by where σ m , σ L , and σ H are the errors of m, L, and H, respectively. Figure 10 shows the typical relation between the measured displacement and the applied force (open squares). The characteristic can be fitted by a linear line as drawn by a solid line in Fig. 10; hence the linearly is well maintained over the force range including the thrust detected in the present experiment. The fitted line gives the calibration coefficient as ∼ 0.475 mN/µm , which is confirmed to be unchanged after venting the chamber.
rf power transfer efficiency. The rf power transfer efficiency η p is defined by a ratio of the rf power absorbed by the plasma to the net power from the rf generator. Assuming that the rf power is consumed by the plasma and the rf circuit including the antenna, the power transfer efficiency can be written as www.nature.com/scientificreports/ where R p , R vac , R total are the resistances of the plasma, the rf circuit, and the total load (the plasma and the circuit) during the plasma production. R total and R vac are estimated from the rf antenna current measured by a current probe (with an accuracy of ± 1% ), and the net rf powers with the plasma and with no plasma (i.e., no gas), respectively, where the impedance matching circuit is tuned so as to minimize the reflected rf power. The method to obtain the rf power transfer efficiency has been widely used 66,67 and has also shown good agreement with a theoretical model 68 . Since the rf antenna is water cooled to being maintained at a constant temperature, the resistance of the rf antenna is considered to be unchanged during the experiment.
Langmuir probe. The Langmuir probe used in the present experiment has a radially facing 3-mm-diameter planar tip, where the one side of the tip surface is covered by a ceramic paste. A bipolar voltage power supply is connected to the probe via a resistor. By sweeping the bias voltage, the current-voltage characteristic can be obtained, where the current signal is obtained by measuring the voltage across the resistor via an precise isolation amplifier. The tangential line of the electron current in a semi-logarithm plot of the current-voltage characteristic can give an electron temperature. The measured electron temperature T e at (r, z) = (0, 0) is 5.3 ± 0.5 eV for the m = +1 mode antenna and 5.4 ± 0.6 eV for the m = 0 mode antenna. The plasma density n p can be estimated from the ion saturation current I is of the Langmuir probe, which is measured by biasing the probe to −70 V and given by where e, u B , and S are the elementary charge, the Bohm velocity given by u B = √ k B T e /m i with the Boltzmann constant k B and the ion mass m i , and the collecting surface area of the Langmuir probe, respectively. By mounting the Langmuir probe on the movable stage installed inside the vacuum chamber, the two-dimensional profile of the ion saturation current I is can be measured.

Data availability
The data presented in this study are available from the author (K.T.) upon reasonable request.