<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Phys.</journal-id>
<journal-title>Frontiers in Physics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Phys.</abbrev-journal-title>
<issn pub-type="epub">2296-424X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fphy.2019.00227</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Fast and Automatic Control of a Frequency-Tuned Radiofrequency Plasma Source</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Takahashi</surname> <given-names>Kazunori</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/152749/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Hanaoka</surname> <given-names>Kengo</given-names></name>
</contrib>
<contrib contrib-type="author">
<name><surname>Ando</surname> <given-names>Akira</given-names></name>
</contrib>
</contrib-group>
<aff><institution>Department of Electrical Engineering, Tohoku University</institution>, <addr-line>Sendai</addr-line>, <country>Japan</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Agnes Granier, Centre National de la Recherche Scientifique (CNRS), France</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Alan A. Howling, &#x000C9;cole Polytechnique F&#x000E9;d&#x000E9;rale de Lausanne, Switzerland; Saikat Chakraborty Thakur, University of California, San Diego, United States</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Kazunori Takahashi <email>kazunori&#x00040;ecei.tohoku.ac.jp</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Plasma Physics, a section of the journal Frontiers in Physics</p></fn></author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>01</month>
<year>2020</year>
</pub-date>
<pub-date pub-type="collection">
<year>2019</year>
</pub-date>
<volume>7</volume>
<elocation-id>227</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>10</month>
<year>2019</year>
</date>
<date date-type="accepted">
<day>10</day>
<month>12</month>
<year>2019</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2020 Takahashi, Hanaoka and Ando.</copyright-statement>
<copyright-year>2020</copyright-year>
<copyright-holder>Takahashi, Hanaoka and Ando</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>A frequency-tuned radiofrequency (rf) plasma source is automatically controlled to minimize the rf power reflection and to maintain the constant net rf power corresponding to the forward power minus the reflected power. The experiment is performed with a power amplifier operational for the frequency of 37&#x02013;43 MHz, a compact helicon source consisting of a loop antenna, a solenoid, and an insulator tube. The rf power is supplied to the rf antenna via an impedance matching circuit consisting of only fixed small ceramic capacitors. It is demonstrated that the reflection coefficient of the rf power is minimized within &#x0007E;10 ms, and the net power is successfully maintained at a constant level.</p></abstract>
<kwd-group>
<kwd>rf plasma source</kwd>
<kwd>impedance matching</kwd>
<kwd>rf generator</kwd>
<kwd>low-pressure plasma</kwd>
<kwd>compact plasma source</kwd>
</kwd-group>
<counts>
<fig-count count="9"/>
<table-count count="0"/>
<equation-count count="10"/>
<ref-count count="25"/>
<page-count count="8"/>
<word-count count="4946"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Radiofrequency (rf) plasma source operated in the frequency of 1&#x02013;100 MHz has been widely used for various applications such as industrial plasmas [<xref ref-type="bibr" rid="B1">1</xref>, <xref ref-type="bibr" rid="B2">2</xref>], space propulsion devices [<xref ref-type="bibr" rid="B3">3</xref>, <xref ref-type="bibr" rid="B4">4</xref>], and fusion plasmas [<xref ref-type="bibr" rid="B5">5</xref>]. To transfer the rf power to a load, the load impedance has to be matched to the output impedance of the generator and the characteristic impedance of a power transmission line, being typically 50 &#x003A9;. The impedance matching is generally performed by inserting two variable capacitors and an inductor (if necessary), where the capacitances are adjusted by mechanical motion.</p>
<p>Fundamental studies on and development of space propulsion devices utilizing the rf plasma sources have vigorously progressed in recent years [<xref ref-type="bibr" rid="B6">6</xref>&#x02013;<xref ref-type="bibr" rid="B11">11</xref>]. Since the weight and size of the rf system mountable on the spacecraft is very limited, development of the compact rf system is required for the rf plasma thruster. Further addressing space propulsion, the time lag of the communication between the spacecraft and the ground control station would be a serious problem for the impedance matching. Therefore, automatic and fast control of the impedance tuning is useful for further development of the rf plasma thruster.</p>
<p>The compact, automatically controlled, and fast-controlled rf plasma system is also useful for terrestrial industrial plasma devices such as plasma sputtering and etching reactors in a recently developed new concept of the semiconductor manufacturing system, being called &#x0201C;Minimal Fab&#x0201D; [<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B13">13</xref>]. In the traditional semiconductor manufacturing system, the diameter of the silicon wafer has been increased to yield sufficient economic benefit by low-mix high-volume production in the large-scale manufacturing system (called &#x0201C;Mega Fab&#x0201D;). On the other hand, Minimal Fab consists of various manufacturing tools in processing a half-inch wafer. Since the wafer is transferred between the processing tools in a localized clean wafer transfer system, no clean room is required for the device production; the facility and production costs can be reduced for the high-mix and low-volume products. The frame and wafer transfer structures are strictly standardized to connect the processing tools. Therefore, all of the components required for each tool, e.g., a plasma reactor, a pumping system, an rf generator, the wafer transfer mechanics, and the programmable logic controller, and so on, have to be contained within the compact frame. Furthermore, the precise control of the pulsed plasma would also be required for some plasma processing methods, e.g., plasma etching by Bosch process [<xref ref-type="bibr" rid="B14">14</xref>], where gas species are alternatively switched for etching and passivation.</p>
<p>As described above, the available space for the rf system is strictly limited in the applications; hence, the development of the compact rf system for the plasma source is a significantly important technology. To overcome this issue, a frequency-variable tuning has been tested in several laboratory experiments. Charles et al. has used the frequency-variable power generator operated at around 13.56 &#x000B1; 0.678 MHz with the fixed small ceramic capacitors for impedance tuning [<xref ref-type="bibr" rid="B15">15</xref>]. One of the authors has used an amplifier operational for the broader range of the frequency to match the impedance for the industrial plasma source [<xref ref-type="bibr" rid="B16">16</xref>], where the frequency was manually adjusted for the initial laboratory test. These experiments have been performed with the linear rf amplifier classified as class AB, while class D or class E switching amplifiers have also been used for the frequency-tunable matching [<xref ref-type="bibr" rid="B17">17</xref>&#x02013;<xref ref-type="bibr" rid="B20">20</xref>]. The traditional impedance matching technique can perform nearly perfect impedance matching by adjusting two capacitors, while it would be difficult to obtain the zero power reflection for the frequency tuning since only one control parameter (frequency) is adjustable. Since the primary parameter deciding the plasma characteristics in terms of the rf system would be the net power given by the forward power minus the reflected power, it would be favorable to maintain the constant net power during plasma production rather than the forward power when the perfect impedance matching is difficult.</p>
<p>The aim of the present study is to develop a compact, automatic, and fast-controlled rf system for future applications of the rf plasma source to the terrestrial plasma processing tools in &#x0201C;Minimal Fab&#x0201D; and the space propulsion device called a helicon thruster. By mounting an autonomously operational controller in the rf system, frequency tuning for the impedance matching and power control for maintaining the constant net rf power, which corresponds to the forward power minus the reflected power, are simultaneously performed. Details of the rf system providing the impedance tuning within about 10 ms is described here.</p></sec>
<sec id="s2">
<title>Experimental Setup</title>
<p><xref ref-type="fig" rid="F1">Figure 1A</xref> shows the schematic diagram of the plasma reactor used for the present experiment, which is very similar to that of the previous experiment [<xref ref-type="bibr" rid="B16">16</xref>]. The experiment is performed with a weakly magnetized rf plasma source attached to a 114 mm<sup>3</sup> vacuum chamber, which is evacuated by a 70-l/s turbomolecular pump to a base pressure of about 10<sup>&#x02212;4</sup> Pa. The source consists of a 50-mm-outer-diameter, 42-mm-inner-diameter, and 75-mm-long ceramic source tube wound by a double-turn rf loop antenna and a solenoid providing a static magnetic field. The solenoid is situated at <italic>z</italic> = &#x02212;15 mm, and dc current of <italic>I</italic><sub>B</sub> = 4.5 A is supplied to the solenoid from a dc power supply, providing a magnetic field strength of about 45 G at the solenoid center, where the open source exit is defined as <italic>z</italic> = 0. The rf loop antenna consisting of a 4-mm-diameter copper tube is wound at <italic>z</italic> = &#x02212;50 &#x000B1; 10 mm, where the antenna radius is about 27 mm. The upstream side of the ceramic tube is terminated at <italic>z</italic> = &#x02212;75 mm by a grounded vacuum flange having a gas injection port and a small viewport. Argon gas is introduced from the gas injection port via a mass flow controller, and the gas flow rate is maintained at 10 sccm, giving the chamber pressure of about 0.65 Pa, which is measured by an ionization gauge attached to the side flange of the chamber. The rf antenna is powered by a broadband class AB rf amplifier via an impedance matching circuit consisting of only fixed ceramic capacitors (<italic>C</italic><sub>1</sub> and <italic>C</italic><sub>2</sub>) as drawn in <xref ref-type="fig" rid="F1">Figure 1B</xref>; the argon plasma is produced inside the source tube. <italic>R</italic><sub>p</sub>, <italic>L</italic><sub>p</sub>, <italic>R</italic><sub>a</sub>, and <italic>L</italic><sub>a</sub> in <xref ref-type="fig" rid="F1">Figure 1B</xref> are the plasma resistance, plasma inductance, antenna resistance, and antenna inductance, respectively, which will be discussed for the matching network design in &#x0201C;Design of Matching Network.&#x0201D; <xref ref-type="fig" rid="F2">Figure 2</xref> shows the photograph of the source, including the impedance matching circuit; the matching circuit can be made very compact as shown in the photograph of <xref ref-type="fig" rid="F2">Figure 2</xref> as it does not contain the variable capacitors. To reduce the power loss at the capacitors, high-Q and ultra-low equivalent series resistance capacitors (HQCC series of AVX Corporation) are mounted on the circuit. The plasma source structure and the matching circuit are air-cooled by using a 40-mm-diameter small fan; it is confirmed that the plasma source can be operated in steady state (at least for 20&#x02013;30 min in the laboratory experiment) for 200 W rf power. It should be noted that a similar previous experiment has shown plasma density of about 10<sup>18</sup> m<sup>&#x02212;3</sup> for the rf power of about 100 W [<xref ref-type="bibr" rid="B16">16</xref>]. An optical emission intensity <italic>I</italic><sub>opt</sub> is measured by a photodetector via an optical fiber connected to the viewport of the upstream flange to very roughly know the relative plasma density.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>(A)</bold> Schematic diagram of a weakly magnetized rf plasma source attached to the 114 mm<sup>3</sup> vacuum chamber. <bold>(B)</bold> Schematic of an impedance matching circuit containing a plasma resistance <italic>R</italic><sub>p</sub> and a plasma inductance <italic>L</italic><sub>a</sub>, where <italic>C</italic><sub>1</sub> and <italic>C</italic><sub>2</sub> are chosen as 390 and 73 pF, respectively. <italic>L</italic><sub>a</sub> and <italic>R</italic><sub>a</sub> are the antenna inductance and resistance, respectively.</p></caption>
<graphic xlink:href="fphy-07-00227-g0001.tif"/>
</fig>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Photograph of the plasma source and the impedance matching circuit consisting of the fixed ceramic capacitors.</p></caption>
<graphic xlink:href="fphy-07-00227-g0002.tif"/>
</fig>
<p><xref ref-type="fig" rid="F3">Figure 3</xref> shows the schematic diagram of the presently developed rf system. An rf signal from a voltage-controlled oscillator (VCO) is input into a voltage variable attenuator (VVA) and further input into a 200-W class AB rf power amplifier (AMP) operational between 37 and 43 MHz. Both the output impedance of the amplifier and the coaxial cables are 50 &#x003A9;. It should be noted that the dc&#x02013;rf conversion efficiency of the power amplifier is about 75%; the bandwidth of the amplifier is limited to 37&#x02013;43 MHz to ensure the conversion efficiency. The amplified rf power is transferred to a load via a bi-directional coupler (50 &#x000B1; 0.5 dB coupling and 20 dB directivity). The rf signals coupled with the forward and reflected powers are converted into dc voltages by using zero-bias schottky diode rf detectors and low-pass filters (LPFs); then, the absolute value of the reflection coefficient |&#x00393;<sub>r</sub>| given by</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mo>&#x00393;</mml:mo></mml:mrow><mml:mrow><mml:mtext>r</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mtext>r</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mtext>f</mml:mtext></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>can be calculated. The driving frequency and the output power can be controlled by the voltage signals to the VCO and the VVA. To autonomously control both the frequency and the rf power, a CompactRIO Single-Board Controller (National Instruments, sbRIO-9637) is mounted on the system. The controller includes the Linux Real-Time operating system (OS), field-programmable gate array (FPGA), analog&#x02013;digital (AD) converters, and digital&#x02013;analog (DA) converters; it has ports of the digital input/output (DIO), analog input (AI), and analog output (AO). The voltage signal to the VCO and the VVA can be controlled by the AO channels, and the dc voltage corresponding to the forward and reflected powers are monitored by the AI channels, where the DIO channels are used to control the RF on/off and to monitor some signals from the thermal interlock. Since the CompactRIO Single-Board Controller has the OS, it can autonomously operate even if a personal computer (PC) is disconnected. The detailed procedure of the autonomous control is described in the &#x0201C;Results.&#x0201D;</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Schematic diagram of the rf system containing a board controller including the FPGA and the AD/DA converters.</p></caption>
<graphic xlink:href="fphy-07-00227-g0003.tif"/>
</fig></sec>
<sec id="s3">
<title>Design of Matching Network</title>
<p>Since it is difficult to perfectly match the load impedance to the output impedance of the rf amplifier for the frequency tuning system as described in &#x0201C;Introduction,&#x0201D; it is significantly important to select proper capacitors in the matching network. Many modeling studies have been performed to characterize the electrical property of the rf plasma source over the wide range of parameters [<xref ref-type="bibr" rid="B21">21</xref>]. However, the very simple equivalent circuit shown in <xref ref-type="fig" rid="F1">Figure 1B</xref> is used here to explore the suitable capacitances for the matching network since the detailed model strongly depends on the specification of the plasma source. The total impedance <italic>Z</italic><sub>L</sub> of the matching circuit including the plasma in <xref ref-type="fig" rid="F1">Figure 1B</xref> is given by</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mo>&#x003C9;</mml:mo><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:msub><mml:mrow><mml:mi>Z</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>R</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mi>j</mml:mi><mml:mo>&#x003C9;</mml:mo><mml:mi>L</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mo>&#x003C9;</mml:mo><mml:msub><mml:mrow><mml:mi>C</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E4"><label>(4)</label><mml:math id="M4"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mi>L</mml:mi><mml:mo>=</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msub><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>The reflection coefficient given by Equation (1) can be rewritten as</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M5"><mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:mrow><mml:msub><mml:mo>&#x00393;</mml:mo><mml:mtext>r</mml:mtext></mml:msub></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>Z</mml:mi><mml:mtext>L</mml:mtext></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mi>Z</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>Z</mml:mi><mml:mtext>L</mml:mtext></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>Z</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mtext>&#x000A0;</mml:mtext><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
<p>where <italic>Z</italic><sub>0</sub> is the output impedance (50 &#x003A9;) of the amplifier being equal to the characteristic impedance of the coaxial cable connecting the amplifier and the matching network. One of the experiments measured the resistance and the reactance of the load including the rf antenna and the plasma for the inductively coupled plasma source [<xref ref-type="bibr" rid="B22">22</xref>], where the equivalent circuit is very similar to <xref ref-type="fig" rid="F1">Figure 1B</xref>. The results showed the plasma resistance of about a few ohms, which is higher than the typical circuit resistance of an order of 10<sup>&#x02212;1</sup> &#x003A9;. Furthermore, the reactance is mostly unchanged or decreases a fraction with an increase in the discharge power; it can be considered that the plasma reactance is much smaller than the antenna reactance (<italic>L</italic> &#x0007E; <italic>L</italic><sub>a</sub>). The inductance of the <italic>N</italic>-turn antenna having the radius <italic>a</italic> and the axial length <italic>l</italic> is generally given by</p>
<disp-formula id="E6"><label>(6)</label><mml:math id="M6"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mtext>a</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mi>&#x003C0;</mml:mi><mml:msup><mml:mrow><mml:mi>a</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mi>l</mml:mi></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003BC;<sub>0</sub> is the vacuum permeability. <italic>K</italic><sub>n</sub> is the Nagaoka coefficient given by [<xref ref-type="bibr" rid="B23">23</xref>]</p>
<disp-formula id="E7"><label>(7)</label><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mtext>n</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>4</mml:mn></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:msqrt><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:msqrt></mml:mrow></mml:mfrac><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mi>K</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mfrac><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mn>2</mml:mn><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mi>E</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E8"><label>(8)</label><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:msup><mml:mrow><mml:mi>l</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>/</mml:mo><mml:mn>4</mml:mn><mml:msup><mml:mrow><mml:mi>a</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>K</italic>(<italic>k</italic>) and <italic>E</italic>(<italic>k</italic>) are the complete elliptic integrals of the first and second kind, respectively. The numerically calculated inductance for the given antenna geometry of <italic>l</italic> = 20 mm, <italic>a</italic> = 27 mm, and <italic>N</italic> = 2 is <italic>L</italic><sub>a</sub> &#x0007E; 0.235 &#x003BC;H. The two capacitors in <xref ref-type="fig" rid="F1">Figure 1B</xref> are chosen so as to have a minimum value of |&#x00393;<sub>r</sub>| within the frequency range of 37&#x02013;43 MHz and nearly zero value of |&#x00393;<sub>r</sub>| for <italic>R</italic> of a few ohms.</p>
<p><xref ref-type="fig" rid="F4">Figure 4A</xref> shows the typical calculation of |&#x00393;<sub>r</sub>| as a function of the frequency for various values of <italic>R</italic>, where <italic>C</italic><sub>1</sub> = 390 pF, <italic>C</italic><sub>2</sub> = 73 pF, and <italic>L</italic> = 0.235 &#x003BC;H are used. It is found that the nearly zero value of |&#x00393;<sub>r</sub>| can be obtained for <italic>R</italic> = 2 &#x003A9; at the frequency of &#x0007E;40 MHz. It should be noted that a slight decrease in |&#x00393;<sub>r</sub>| at &#x0007E; 40 MHz is still obtained for <italic>R</italic> = 0.25 &#x003A9; case, which would be enough to ignite the discharge.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p><bold>(A)</bold> Reflection coefficient |<italic>&#x00393;</italic><sub>r</sub>| calculated by Equation (5) as a function of the frequency for various values of the resistance <italic>R</italic>, where <italic>C</italic><sub>1</sub> = 390 pF, <italic>C</italic><sub>2</sub> = 73 pF, and the antenna inductance <italic>L</italic><sub>a</sub> &#x0007E; 0.235 &#x003BC;H calculated by Equation (6) are used for the calculation. <bold>(B)</bold> Reflection coefficient for <italic>R</italic> = 2 &#x003A9; when changing the antenna inductance by &#x000B1; 10 % of the <italic>L</italic><sub>a</sub> used in <bold>(A)</bold>.</p></caption>
<graphic xlink:href="fphy-07-00227-g0004.tif"/>
</fig>
<p>For the high frequency above 10 MHz, the reactance is significantly changed by the geometrical error of the antenna, the parasitic capacitance and inductance in the circuit. <xref ref-type="fig" rid="F4">Figure 4B</xref> shows the calculated |&#x00393;<sub>r</sub>| for the <italic>R</italic> = 2 &#x003A9; cases with the inductances of <italic>L</italic><sub>a</sub> and <italic>L</italic><sub>a</sub> &#x000B1; 10 %. The frequency giving the minimum |&#x00393;<sub>r</sub>| is found to be shifted by the change of the inductance. Therefore, the measurement of the reflection coefficient by a vector network analyzer is required to be performed, whereby the values of the capacitance should be adjusted. However, the above-mentioned design procedure is useful to roughly estimate the capacitors in the matching network.</p></sec>
<sec sec-type="results" id="s4">
<title>Results</title>
<p>Before showing the automatic and fast control of the frequency tuning rf system, the basic characteristics of the loads without and with the plasma are shown. <xref ref-type="fig" rid="F5">Figure 5A</xref> shows the reflection coefficient |&#x00393;<sub>r</sub>| as a function of the rf frequency. The data plotted by crosses and filled circles are taken with no gas (no plasma) and with argon flow rate of <italic>C</italic><sub>Ar</sub> = 10 sccm and solenoid current of <italic>I</italic><sub>B</sub> = 4.5 A, respectively. It is noted that the voltage signals to the VVA are maintained at a constant level, and the data are taken by changing the voltage signal to the VCO, i.e., the rf frequency. The optical emission intensity <italic>I</italic><sub>opt</sub>, forward rf power <italic>P</italic><sub>f</sub>, and reflected rf power <italic>P</italic><sub>r</sub> measured for <italic>C</italic><sub>Ar</sub> = 10 sccm and <italic>I</italic><sub>B</sub> = 4.5 A are plotted by filled circles in <xref ref-type="fig" rid="F5">Figure 5B</xref> and by open circles and filled squares in <xref ref-type="fig" rid="F5">Figure 5C</xref>, respectively. The net power <italic>P</italic><sub>net</sub> is simply calculated by <italic>P</italic><sub>net</sub> = <italic>P</italic><sub>f</sub> &#x02013; <italic>P</italic><sub>r</sub> and plotted by crosses in <xref ref-type="fig" rid="F5">Figure 5C</xref>. The optical emission intensity shows that plasma can be produced for the frequency range of 39.5&#x02013;42 MHz.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p><bold>(A)</bold> Measured reflection coefficient |<italic>&#x00393;</italic><sub>r</sub>| as a function of the rf frequency for <italic>C</italic><sub>Ar</sub> = 0 sccm and <italic>I</italic><sub>B</sub> = 0 A case (<italic>crosses</italic>) and for <italic>C</italic><sub>Ar</sub> = 10 sccm and <italic>I</italic><sub>B</sub> = 4.5 A case (<italic>filled circles</italic>). <bold>(B)</bold> The optical emission intensity <italic>I</italic><sub>opt</sub> measured from the viewport on the upstream flange. <bold>(C)</bold> The forward (<italic>P</italic><sub>f</sub>, <italic>open circles</italic>), reflected (<italic>P</italic><sub>r</sub>, <italic>filled squares</italic>), and net (<italic>P</italic><sub>net</sub>, <italic>crosses</italic>) rf powers.</p></caption>
<graphic xlink:href="fphy-07-00227-g0005.tif"/>
</fig>
<p>For the no plasma case (<italic>C</italic><sub>Ar</sub> = 0 sccm and <italic>I</italic><sub>B</sub> = 0 A), the reflection coefficient |&#x00393;<sub>r</sub>| is found to be minimized at the frequency of &#x0007E;40.8 MHz as seen by the crosses in <xref ref-type="fig" rid="F5">Figure 5A</xref>, while the minimum |&#x00393;<sub>r</sub>| is obtained at &#x0007E;41.2 MHz for the case of <italic>C</italic><sub>Ar</sub> = 10 sccm and <italic>I</italic><sub>B</sub> = 4.5 A. Since the frequency giving the minimum |&#x00393;<sub>r</sub>| is changed by the presence of the plasma and the variation of the impedance is expected to be non-linear, frequency tuning is required to minimize the power reflection when using only the fixed capacitors in the impedance matching circuit. It is clearly observed that the optical emission intensity <italic>I</italic><sub>opt</sub> increases when the power reflection is reduced (as shown in <xref ref-type="fig" rid="F5">Figures 5A,B</xref>). Although the gain (attenuation) of the VVA is maintained at a constant level in <xref ref-type="fig" rid="F5">Figure 5</xref>, the output power is fairly changed by the frequency due to the not constant gain of the amplifier and the impedance mismatch. Therefore, the voltage signal to the VVA has to be adjusted to control the output power from the amplifier. Furthermore, the optical emission intensity <italic>I</italic><sub>opt</sub> seems to fit the net power <italic>P</italic><sub>net</sub>, implying the importance of maintaining not the forward power <italic>P</italic><sub>f</sub> but the net power <italic>P</italic><sub>net</sub>.</p>
<p>In order to automatically control the frequency and the gain of the VVA to minimize the reflection coefficient |&#x00393;<sub>r</sub>| and to adjust the net rf power <italic>P</italic><sub>net</sub> close to the target value, the CompactRIO Single-Board Controller is controlled by a Labview program. <xref ref-type="fig" rid="F6">Figure 6</xref> shows the brief flow chart of the controller. Before turning on the rf power, the variation step of the voltage signal to the VCO (i.e., frequency step, &#x00394;<italic>f</italic>), the variation step of the voltage signal to the VVA (i.e., the gain step, &#x00394;<italic>G</italic>), the target value of the net power <italic>P</italic><sub>set</sub>, the initial values of the frequency <italic>f</italic><sub>0</sub>, and the initial gain of the VVA <italic>G</italic><sub>0</sub> are given as the external parameters. During the period of &#x0201C;RF ON,&#x0201D; the forward and reflected powers (<italic>P</italic><sub>f</sub>, <italic>P</italic><sub>r</sub>) are measured, and the net power <italic>P</italic><sub>net</sub> and the reflection coefficient |&#x00393;<sub>r</sub>| are calculated on the CompactRIO Single-Board Controller. Based on the measurement of <italic>P</italic><sub>net</sub> and |&#x00393;<sub>r</sub>|, the frequency of the VCO and the gain of the VVA are changed according to the equations of</p>
<disp-formula id="E9"><label>(9)</label><mml:math id="M9"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mo>=</mml:mo><mml:mi>f</mml:mi><mml:mo>-</mml:mo><mml:mtext>sign&#x000A0;</mml:mtext><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mo>&#x00393;</mml:mo></mml:mrow><mml:mrow><mml:mtext>r</mml:mtext></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>f</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mo>&#x00394;</mml:mo><mml:mi>f</mml:mi><mml:mtext>&#x000A0;</mml:mtext><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mrow><mml:mo>&#x00393;</mml:mo></mml:mrow><mml:mrow><mml:mtext>r</mml:mtext></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E10"><label>(10)</label><mml:math id="M10"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:mi>G</mml:mi><mml:mo>=</mml:mo><mml:mi>G</mml:mi><mml:mo>-</mml:mo><mml:mtext>sign&#x000A0;</mml:mtext><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mtext>net</mml:mtext></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mtext>set</mml:mtext></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x00394;</mml:mo><mml:mi>G</mml:mi><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>By multiplying the signs of d|&#x00393;<sub>r</sub>|/d<italic>f</italic> and (<italic>P</italic><sub>net</sub> &#x02013; <italic>P</italic><sub>set</sub>) as in Equations (9) and (10), the voltage signals to the VCO and VVA can be varied toward the minimum reflection coefficient and the target value of the net power, respectively. Around the frequency giving the reduced reflection coefficient, the variation step of the frequency is reduced by multiplying |&#x00393;<sub>r</sub>| in Equation (9), providing the stable operation around the frequency giving a small value of |&#x00393;<sub>r</sub>|.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Flow chart of the automatic, fast, and autonomous control of the rf system. &#x00394;<italic>f</italic>, &#x00394;<italic>G, P</italic><sub>set</sub>, <italic>f</italic><sub>0</sub>, and <italic>G</italic><sub>0</sub> are the variation step of the frequency, the variation step of the gain of the VVA, the target value of the net rf power, the initial frequency, and the initial gain of the VVA, respectively.</p></caption>
<graphic xlink:href="fphy-07-00227-g0006.tif"/>
</fig>
<p>Typical temporal evolutions of the frequency <italic>f</italic>, the reflection coefficient |&#x00393;<sub>r</sub>|, the optical emission intensity <italic>I</italic><sub>opt</sub>, the rf powers of <italic>P</italic><sub>f</sub>, <italic>P</italic><sub>r</sub>, and <italic>P</italic><sub>net</sub> are shown in <xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F8">8</xref>, where the initial frequency is set at <italic>f</italic><sub>0</sub> = 37.6 MHz in <xref ref-type="fig" rid="F7">Figure 7</xref> and 40.6 MHz in <xref ref-type="fig" rid="F8">Figure 8</xref>, respectively. The target value of the net power is chosen as <italic>P</italic><sub>set</sub> = 40 W here, and the rf power is turned on for <italic>t</italic> = 0&#x02013;100 ms for both cases. As can be seen in both figures, the reflection coefficient is rapidly decreased by automatic frequency tuning; then, plasma is produced and optical emission is detected. For the case of <xref ref-type="fig" rid="F7">Figure 7</xref>, longer time is required to reduce the power reflection since the initial frequency for this case is more distant from the matching point than that in the case of <xref ref-type="fig" rid="F8">Figure 8</xref>. However, stable plasma production can be obtained at about <italic>t</italic> &#x0007E; 10 ms. When choosing the initial frequency properly, stable plasma production can be obtained within a few of milliseconds as seen in <xref ref-type="fig" rid="F8">Figure 8</xref>. Furthermore, it is demonstrated that the net power <italic>P</italic><sub>net</sub> equal to the target value <italic>P</italic><sub>set</sub> is successfully sustained during the discharge, as seen in both <xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F8">8</xref>, providing the fairly constant optical emission intensity <italic>I</italic><sub>opt</sub>. These data demonstrate that the automatic and fast control of the impedance and power tuning in the presently developed system can provide stable plasma production during pulse, which can provide high-speed minimization of power reflection.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Temporal evolutions of <bold>(A)</bold> the rf frequency <italic>f</italic>, the optical emission intensity <italic>I</italic><sub>opt</sub>, and the reflection coefficient |<italic>&#x00393;</italic><sub>r</sub>|, and <bold>(B)</bold> the rf powers <italic>P</italic><sub>f</sub>, <italic>P</italic><sub>r</sub>, and <italic>P</italic><sub>net</sub>. The initial frequency and the target value of the net power are chosen as <italic>f</italic><sub>0</sub> = 37.6 MHz and <italic>P</italic><sub>set</sub> = 40 W, respectively.</p></caption>
<graphic xlink:href="fphy-07-00227-g0007.tif"/>
</fig>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p><bold>(A,B)</bold> The same data set as in <xref ref-type="fig" rid="F7">Figures 7A,B</xref> for <italic>f</italic><sub>0</sub> = 40.6 MHz and <italic>P</italic><sub>set</sub> = 40 W.</p></caption>
<graphic xlink:href="fphy-07-00227-g0008.tif"/>
</fig>
<p>One of the applications of the presently developed rf system is industrial plasma processing tool such as a plasma etcher operated in Bosch process, where etching and passivation gases are alternatively switched, e.g., as in Chang et al. [<xref ref-type="bibr" rid="B24">24</xref>]. Therefore, impedance tuning has to be finished for a much shorter time than the discharge pulse. Thus, high-speed impedance tuning would be an important technique to stably sustain the discharge when shortening the periods of each process. The plasma etcher utilizing the inductively coupled or helicon plasma source generally uses two rf amplifier for plasma production and substrate bias voltage [<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B25">25</xref>]. Moreover, the available space for the rf system in the recently proposed &#x0201C;Minimal Fab&#x0201D; (see &#x0201C;Introduction&#x0201D;) is significantly limited [<xref ref-type="bibr" rid="B12">12</xref>]. Therefore, the system having two rf generators has to be made compact to be fitted to the frame designed for Minimal Fab. <xref ref-type="fig" rid="F9">Figure 9</xref> shows the photograph of the presently developed rf system including the two rf amplifiers and the CompactRIO Single-Board Controller. It can be successfully made within the size of 220 mm in width, 135 mm in height, and 260 mm in depth and would be contained in the frame designed for Minimal Fab. The target value of the net power and the initial frequency can be set by putting the analog voltage signal, which can facilitate the initial test of the plasma source.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p>Photograph of the rf system containing the two rf generators, the CompactRIO Single-Board Controller, and the dc power supply. The total size is within 220 mm in width, 135 mm in height, and 260 mm in depth.</p></caption>
<graphic xlink:href="fphy-07-00227-g0009.tif"/>
</fig></sec>
<sec sec-type="conclusions" id="s5">
<title>Conclusion</title>
<p>An automatic, fast, and autonomously controlled rf plasma source operational at variable frequency is developed. The rf system can be operated for the frequency range of 37&#x02013;43 MHz, and the maximum rf power is about 200 W, where the impedance matching circuit consists of only the fixed ceramic capacitors, leading to a compact plasma source. Frequency tuning is accomplished within about 10 ms, and the net power corresponding to the forward power minus the reflected power is successfully maintained at a constant level during the pulse, which stabilizes the plasma source operation. For future application to the plasma etcher in Minimal Fab, the rf system containing the two rf generators is made compact, being contained within the size of 220 mm in width, 135 mm in height, and 260 mm in depth. The system will also be useful for space propulsion device utilizing the rf plasma source. The actual process of plasma etching and application to propulsion device remain further as challenging issues.</p></sec>
<sec sec-type="data-availability-statement" id="s6">
<title>Data Availability Statement</title>
<p>All datasets generated for this study are included in the article/supplementary material.</p></sec>
<sec id="s7">
<title>Author Contributions</title>
<p>The rf generator was designed and assembled by KT. The controller program was written by KT. KH installed the generator to the plasma source designed and assembled by KT. The data were taken and analyzed by KT and KH. The results were discussed by KT, KH, and AA. KT wrote the manuscript. KH and AA reviewed the manuscript.</p>
<sec>
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec>
</sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hopwood</surname> <given-names>J</given-names></name></person-group>. <article-title>Review of inductively coupled plasmas for plasma processing</article-title>. <source>Plasma Sources Sci Technol.</source> (<year>1992</year>) <volume>1</volume>:<fpage>109</fpage>. <pub-id pub-id-type="doi">10.1088/0963-0252/1/2/006</pub-id></citation></ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Donnelly</surname> <given-names>VM</given-names></name> <name><surname>Kornblit</surname> <given-names>A</given-names></name></person-group>. <article-title>Plasma etching: yesterday, today, and tomorrow</article-title>. <source>J Vac Sci Technol A.</source> (<year>2013</year>) <volume>31</volume>:<fpage>050825</fpage>. <pub-id pub-id-type="doi">10.1116/1.4819316</pub-id></citation></ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Charles</surname> <given-names>C</given-names></name></person-group>. <article-title>Plasmas for spacecraft propulsion</article-title>. <source>J Phys D Appl Phys.</source> (<year>2009</year>) <volume>42</volume>:<fpage>163001</fpage>. <pub-id pub-id-type="doi">10.1088/0022-3727/42/16/163001</pub-id></citation></ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takahashi</surname> <given-names>K</given-names></name></person-group>. <article-title>Helicon-type radiofrequency plasma thrusters and magnetic plasma nozzles</article-title>. <source>Rev Mod Plasma Phys.</source> (<year>2019</year>) <volume>3</volume>:<fpage>3</fpage>. <pub-id pub-id-type="doi">10.1007/s41614-019-0024-2</pub-id></citation></ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kraus</surname> <given-names>W</given-names></name> <name><surname>Schiesko</surname> <given-names>L</given-names></name> <name><surname>Bonomo</surname> <given-names>F</given-names></name> <name><surname>Fantz</surname> <given-names>U</given-names></name> <name><surname>Heinemann</surname> <given-names>B</given-names></name> <name><surname>Hurlbatt</surname> <given-names>A</given-names></name> <etal/></person-group>. <article-title>First beam extraction experiments at BATMAN Upgrade</article-title>. <source>AIP Conf Proc.</source> (<year>2018</year>) <volume>2052</volume>:<fpage>040004</fpage>. <pub-id pub-id-type="doi">10.1063/1.5083738</pub-id></citation></ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Olsen</surname> <given-names>CS</given-names></name> <name><surname>Maxwell</surname> <given-names>GB</given-names></name> <name><surname>Carter</surname> <given-names>MD</given-names></name> <name><surname>Chang-D&#x000ED;az</surname> <given-names>FR</given-names></name> <name><surname>Giambusso</surname> <given-names>M</given-names></name> <name><surname>Glover</surname> <given-names>TW</given-names></name> <etal/></person-group>. <article-title>Investigation of plasma detachment from a magnetic nozzle in the plume of the VX-200 magnetoplasma thruster</article-title>. <source>IEEE Trans Plasma Sci.</source> (<year>2015</year>) <volume>43</volume>:<fpage>252</fpage>&#x02013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1109/TPS.2014.2321257</pub-id></citation></ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takahashi</surname> <given-names>K</given-names></name> <name><surname>Charles</surname> <given-names>C</given-names></name> <name><surname>Boswell</surname> <given-names>RW</given-names></name></person-group>. <article-title>Approaching the theoretical limit of diamagnetic-induced momentum in a rapidly diverging magnetic nozzle</article-title>. <source>Phys Rev Lett.</source> (<year>2013</year>) <volume>110</volume>:<fpage>195003</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.110.195003</pub-id><pub-id pub-id-type="pmid">23705713</pub-id></citation></ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takahashi</surname> <given-names>K</given-names></name> <name><surname>Ando</surname> <given-names>A</given-names></name></person-group>. <article-title>Laboratory observation of a plasma-flow-state transition from diverging to stretching a magnetic nozzle</article-title>. <source>Phys Rev Lett.</source> (<year>2017</year>) <volume>118</volume>:<fpage>225002</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.118.225002</pub-id><pub-id pub-id-type="pmid">28621999</pub-id></citation></ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Charles</surname> <given-names>C</given-names></name> <name><surname>Boswell</surname> <given-names>RW</given-names></name> <name><surname>Bish</surname> <given-names>A</given-names></name></person-group>. <article-title>Low-weight fixed ceramic capacitor impedance matching system for an electrothermal plasma microthruster</article-title>. <source>J Propul Power.</source> <volume>30</volume>:<fpage>1117</fpage>&#x02013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.2514/1.B35119</pub-id></citation></ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aanesland</surname> <given-names>A</given-names></name> <name><surname>Rafalskyi</surname> <given-names>D</given-names></name> <name><surname>Bredin</surname> <given-names>J</given-names></name> <name><surname>Grondein</surname> <given-names>P</given-names></name> <name><surname>Oudini</surname> <given-names>N</given-names></name> <name><surname>Chabert</surname> <given-names>P</given-names></name> <etal/></person-group>. <article-title>The PEGASES gridded ion&#x02013;ion thruster performance and predictions</article-title>. <source>IEEE Trans Plasma Sci.</source> (<year>2015</year>) <volume>43</volume>:<fpage>321</fpage>&#x02013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/TPS.2014.2369534</pub-id></citation></ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Goebel</surname> <given-names>DM</given-names></name></person-group>. <article-title>Analytical discharge model for RF ion thrusters</article-title>. <source>IEEE Trans Plasma Sci.</source> (<year>2008</year>) <volume>36</volume>:<fpage>2111</fpage>&#x02013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.1109/TPS.2008.2004232</pub-id></citation></ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Khumpuang</surname> <given-names>S</given-names></name> <name><surname>Hara</surname> <given-names>S</given-names></name></person-group>. <article-title>A MOSFET fabrication using a maskless lithography system in clean-localized environment of Minimal Fab</article-title>. <source>IEEE Trans Semicond Manuf.</source> (<year>2015</year>) <volume>28</volume>:<fpage>393</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1109/TSM.2015.2429572</pub-id></citation></ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Khumpuang</surname> <given-names>S</given-names></name> <name><surname>Imura</surname> <given-names>F</given-names></name> <name><surname>Hara</surname> <given-names>S</given-names></name></person-group>. <article-title>Analyses on cleanroom-free performance and transistor manufacturing cycle time of Minimal Fab</article-title>. <source>IEEE Trans Semicond Manuf.</source> (<year>2015</year>) <volume>28</volume>:<fpage>551</fpage>&#x02013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/TSM.2015.2487324</pub-id></citation></ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tanaka</surname> <given-names>H</given-names></name> <name><surname>Ogiso</surname> <given-names>H</given-names></name> <name><surname>Nakano</surname> <given-names>S</given-names></name> <name><surname>Hayami</surname> <given-names>T</given-names></name> <name><surname>Miyazaki</surname> <given-names>T</given-names></name> <name><surname>Khumpuang</surname> <given-names>S</given-names></name> <etal/></person-group>. <article-title>Scallop reduction in Bosch process using a small chamber and rapid gas switching rate</article-title>. <source>IEEJ Trans Sens Micromach.</source> (<year>2016</year>) <volume>136</volume>:<fpage>499</fpage>&#x02013;<lpage>504</lpage>. <pub-id pub-id-type="doi">10.1541/ieejsmas.136.499</pub-id></citation></ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Charles</surname> <given-names>C</given-names></name> <name><surname>Boswell</surname> <given-names>RW</given-names></name> <name><surname>Bish</surname> <given-names>A</given-names></name></person-group>. <article-title>Variable frequency matching to a radiofrequency source immersed in vacuum</article-title>. <source>J Phys D Appl Phys.</source> (<year>2013</year>) <volume>46</volume>:<fpage>365203</fpage>. <pub-id pub-id-type="doi">10.1088/0022-3727/46/36/365203</pub-id></citation></ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takahashi</surname> <given-names>K</given-names></name> <name><surname>Nakano</surname> <given-names>Y</given-names></name> <name><surname>Ando</surname> <given-names>A</given-names></name></person-group>. <article-title>Frequency-tuning radiofrequency plasma source operated in inductively-coupled mode under a low magnetic field</article-title>. <source>J Phys D Appl Phys.</source> (<year>2017</year>) <volume>50</volume>:<fpage>265201</fpage>. <pub-id pub-id-type="doi">10.1088/1361-6463/aa7524</pub-id></citation></ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujita</surname> <given-names>H</given-names></name> <name><surname>Akagi</surname> <given-names>H</given-names></name> <name><surname>Shinohara</surname> <given-names>S</given-names></name></person-group>. <article-title>A 2-MHz 6-kVA voltage-source inverter using low-profile MOSFET modules for low-temperature plasma generators</article-title>. <source>IEEE Trans Power Electronics.</source> (<year>1999</year>) <volume>14</volume>:<fpage>1014</fpage>&#x02013;<lpage>20</lpage>. <pub-id pub-id-type="doi">10.1109/63.803394</pub-id></citation></ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ando</surname> <given-names>A</given-names></name> <name><surname>Komuro</surname> <given-names>A</given-names></name> <name><surname>Matsuno</surname> <given-names>T</given-names></name> <name><surname>Tsumori</surname> <given-names>K</given-names></name> <name><surname>Takeiri</surname> <given-names>Y</given-names></name></person-group>. <article-title>Radio frequency ion source operated with field effect transistor based radio frequency system</article-title>. <source>Rev Sci Instrum.</source> (<year>2010</year>) <volume>81</volume>:<fpage>02B107</fpage>. <pub-id pub-id-type="doi">10.1063/1.3279306</pub-id><pub-id pub-id-type="pmid">20192414</pub-id></citation></ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Takahashi</surname> <given-names>K</given-names></name> <name><surname>Ando</surname> <given-names>A</given-names></name></person-group>. <article-title>Observation of stationary plasma striation and collimated plasma transport in a 100-kHz inductively coupled plasma discharge</article-title>. <source>IEEE Trans Plasma Sci.</source> (<year>2014</year>) <volume>42</volume>:<fpage>2784</fpage>&#x02013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1109/TPS.2014.2308893</pub-id></citation></ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Charles</surname> <given-names>C</given-names></name> <name><surname>Liang</surname> <given-names>W</given-names></name> <name><surname>Raymond</surname> <given-names>L</given-names></name> <name><surname>Rivas-Davila</surname> <given-names>J</given-names></name> <name><surname>Boswell</surname> <given-names>RW</given-names></name></person-group>. <article-title>Vacuum testing of a miniaturized switch mode amplifier powering an electrothermal plasma micro-thruster</article-title>. <source>Front Phys.</source> (<year>2017</year>) <volume>5</volume>:<fpage>36</fpage>. <pub-id pub-id-type="doi">10.3389/fphy.2017.00036</pub-id></citation></ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Lieberman</surname> <given-names>MA</given-names></name> <name><surname>Lichtenberg</surname> <given-names>AJ</given-names></name></person-group>. <source>Principles of plasma discharges and material processing.</source> <edition>2nd edn.</edition> <publisher-loc>Hoboken, NJ</publisher-loc>: <publisher-name>John Wiley &#x00026; Sons, Inc</publisher-name> (<year>2005</year>). <pub-id pub-id-type="doi">10.1002/0471724254</pub-id></citation></ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Godyak</surname> <given-names>VA</given-names></name> <name><surname>Piejak</surname> <given-names>RB</given-names></name> <name><surname>Alexandrovich</surname> <given-names>BM</given-names></name></person-group>. <article-title>Experimental setup and electrical characteristics of an inductively coupled plasma</article-title>. <source>J Appl Phys.</source> (<year>1999</year>) <volume>85</volume>:<fpage>703</fpage>&#x02013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.1063/1.369150</pub-id></citation></ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nagaoka</surname> <given-names>H</given-names></name></person-group>. <article-title>The inductance coefficients of solenoids</article-title>. <source>J Coll Sci Imp Univ Tokyo.</source> (<year>1909</year>) <volume>27</volume>:<fpage>18</fpage>.</citation></ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname> <given-names>C</given-names></name> <name><surname>Wang</surname> <given-names>YF</given-names></name> <name><surname>Kanamori</surname> <given-names>Y</given-names></name> <name><surname>Shih</surname> <given-names>JJ</given-names></name> <name><surname>Kawai</surname> <given-names>Y</given-names></name> <name><surname>Lee</surname> <given-names>CK</given-names></name> <etal/></person-group>. <article-title>Etching submicrometer trenches by using the Bosch process and its application to the fabrication of antireflection structures</article-title>. <source>J Micromech Microeng.</source> (<year>2005</year>) <volume>15</volume>:<fpage>580</fpage>. <pub-id pub-id-type="doi">10.1088/0960-1317/15/3/020</pub-id></citation></ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Perry</surname> <given-names>AJ</given-names></name> <name><surname>Boswell</surname> <given-names>RW</given-names></name></person-group>. <article-title>Fast anisotropic etching of silicon in an inductively coupled plasma reactor</article-title>. <source>Appl Phys Lett.</source> (<year>1989</year>) <volume>55</volume>:<fpage>148</fpage>&#x02013;<lpage>50</lpage>. <pub-id pub-id-type="doi">10.1063/1.102127</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> This work was partially supported by Strategic Core Technology Advancement Program (Supporting Industry Program) by Tohoku Bureau of Economy, Trade and Industry. This is also partially supported by a grant-in-aid for scientific research (18K18907, 18K18746, and 19H00663) from the Japan Society for the Promotion of Science and JAXA.</p></fn>
</fn-group>
</back>
</article>