*Article* **Frequency Tuning Characteristics of a High-Power Sub-THz Gyrotron with Quasi-Optical Cavity**

**Xiaotong Guan 1,2 , Jiayi Zhang 2,3, Wenjie Fu 2,3,\* , Dun Lu 2,3 , Tongbin Yang 2,3 , Yang Yan 2,3 and Xuesong Yuan 2,3**


**Abstract:** Motivated by some emerging high-frequency applications, a high-power frequency-tunable sub-THz quasi-optical gyrotron cavity based on a confocal waveguide is designed in this paper. The frequency tuning characteristics of different approaches, including magnetic field tuning, mirror separation adjustment, and hybrid tuning, have been investigated by particle-in-cell (PIC) simulation. Results predict that it is possible to realize a smooth continuous frequency tuning band with an extraordinarily broad bandwidth of 41.55 GHz, corresponding to a relative bandwidth of 18.7% to the center frequency of 0.22 THz. The frequency tunability is provided by varying the separation distance between two mirrors and correspondingly adjusting the external magnetic field. During the frequency tuning, the output power remains higher than 20 kW, which corresponds to an interaction efficiency of 10%. Providing great advantages in terms of broad bandwidth, smooth tuning, and high power, this research may be conducive to the development of high-power frequency-tunable THz gyrotron oscillators.

**Keywords:** gyrotron; quasi-optical cavity; confocal waveguide; frequency tuning; high power; sub-millimeter wave; terahertz

#### **1. Introduction**

A gyrotron is a typical fast-wave vacuum electron device based on the interaction principle between gyrating electrons and the electromagnetic waves propagating in the waveguide [1]. As one of the most powerful radiation sources, a gyrotron performs with the capability of high-power output from the microwave to terahertz (THz) band [2,3]. Up till now, the world power record for a gyrotron is 2.2 MW at 170 GHz applied for electron cyclotron heating and current drive in the International Thermonuclear Experimental Reactor (ITER) [4]. In recent decades, a continuous frequency-tunable gyrotron operating at a single mode has been especially attractive for some modern high-frequency applications [5], such as high-resolution molecular gas spectroscopy [6], nuclear magnetic resonance spectroscopy enhanced by dynamic nuclear polarization (DNP-NMR) [7], and the direct measurement of positronium hyperfine splitting (Ps-HFS) [8], in which radiation sources are required to be high-power and continuously tunable in a wide frequency range.

According to the principles of an electron cyclotron maser (ECM), the gyrating electrons are able to interact with the electromagnetic (EM) waves efficiently only under the cyclotron resonance condition [2].

$$
\omega - k\_z v\_z \approx s \Omega\_\circ, \ \Omega\_\circ = \frac{eB\_0}{\gamma m\_0}, \ \gamma = 1 + \frac{eV\_0}{m\_0 c^2} \tag{1}
$$

**Citation:** Guan, X.; Zhang, J.; Fu, W.; Lu, D.; Yang, T.; Yan, Y.; Yuan, X. Frequency Tuning Characteristics of a High-Power Sub-THz Gyrotron with Quasi-Optical Cavity. *Electronics* **2021**, *10*, 526. https://doi.org/10.3390/ electronics10050526

Academic Editor: Mikhail Glyavin

Received: 30 January 2021 Accepted: 19 February 2021 Published: 24 February 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

where *ω* and *k<sup>z</sup>* are the angular frequency and axial wavenumber of the EM waves in the interaction space (cavity); *v<sup>z</sup>* is the electron axial velocity; Ω<sup>c</sup> is the relativistic cyclotron frequency of the electrons relative to the static magnetic field strength *B*0, the accelerate voltage *V*0, and the cyclotron harmonic number *s*; *γ* is the relativistic factor of the electrons; *m*<sup>0</sup> and *e* are the relativistic electron rest mass and charge; and *c* is the velocity of light in free space.

From Equation (1), there are two possible approaches to controlling the output frequency: changing the cyclotron frequency Ω<sup>c</sup> of the electrons, or changing the EM waves' frequency *ω* in the cavity. The variation of Ω<sup>c</sup> can be easily achieved by altering the magnetic field *B*<sup>0</sup> or the beam voltage *V*0. However, to meet the requirement of high-power output, a conventional gyrotron normally employs a high *Q*-value cavity (in the order of ~1000), which significantly restricts its resonance bandwidth (∆*f* ~ *f* <sup>0</sup>/*Q*). That is why conventional high-power gyrotrons can only achieve a narrowband frequency by tuning at a fixed mode, or discrete broadband tuning for several modes. Therefore, the key issue with frequency-tunable gyrotrons is broadening the cavity resonance bandwidth as much as possible without reducing its *Q*-value.

For a conventional gyrotron with a cylindrical cavity, several approaches have been proposed and demonstrated for the problem of frequency tuning. One of the most important mechanisms relies on exciting a series of high-order axial modes (HOAMs) in a long gyrotron cavity. With a well-elaborated selection of cavity length and beam current, the frequency region obtained by operating in one axial mode is able to overlap with the frequency region in another axial mode [9]. In principle, continuous broadband frequency tuning can be accessible for gyrotron operation by increasing the axial mode indices. Nowadays, lots of continuous frequency-tunable THz gyrotrons have been successfully developed and applied to DNP-NMR applications at the Massachusetts Institute of Technology (MIT, Cambridge, MA, USA) [10], the Bruker Biospin company in collaboration with the Communications & Power Industries Company (CPI, Palo Alto, CA, USA) [11], and the Research Center for the Development of the Far-Infrared Region of University of Fukui (FIR-UF, Fukui, Japan) [12]. However, to maintain a wide bandwidth, the operating beam currents for these HOAMs are limited at several hundred milliamperes, resulting in a medium power level (less than 100 W). There is the same weakness in other frequency tuning methods, including using an improved multi-section cavity [13], cathode-end power output [14], and using backward-wave components [15]. Recently, the Terahertz Research Center of University of Electronic Science and Technology of China (TRC-UESTC) reported a frequency-tunable HOAMs gyrotron operating at a high beam current and with experimental output higher than 0.45 kW over a 0.79 GHz frequency range [16], which was high-power but not continuously tunable.

As for the gyrotron cavity based on a multi-conductor waveguide, there is another approach for acquiring smooth frequency tunability by adjusting its structural parameters. For example, in a coaxial gyrotron cavity, the cavity eigenfrequency depends on the ratio of the radii of the external and internal conductors. It is possible to realize a continuous frequency tuning by moving the tapered inner conductor longitudinally. Researchers at the Institute of Applied Physics of the Russian Academy of Sciences (IAP, Moscow, Russia) have numerically investigated this tuning mechanism. One result presented the possibility of frequency tuning by 8 GHz at around 394.6 GHz (within a frequency band of about 2%), with an output power of about several hundred watts [17], while another result demonstrated the smooth frequency tuning at one mode by no less than 3.5% around 330 GHz with about 10 kW output power [18]. Lately, a similar mechanical tuning method was applied to a high-power large orbit gyrotron (LOG) with a variable slit-cavity [19]. By changing the slit width transversely, theoretical calculation results suggested a relative frequency bandwidth of 8.5% around 316 GHz with an output power of over 10 kW.

Recently, high-power continuous frequency-tunable radiation over a 3.1 GHz bandwidth around 0.2 THz has been experimentally observed in a quasi-optical gyrotron with a straight confocal waveguide [20], which does not contain a resonance structure. Although

lacking a comprehensive theoretical model to explain it, the experimental result still points out the possibility of generating high-power frequency-tunable THz radiation from a quasi-optical gyrotron with a confocal waveguide.

A quasi-optical cavity provides many good characteristics, such as a high power capacity and a low mode density, which bring positive effects to high-frequency and highpower gyrotron design. In the 1990s, a quasi-optical gyrotron based on a Fabry–Perot cavity generated an output power of 90 kW at 100 GHz, operated at the fundamental cyclotron [21]. The cylindrical confocal waveguide is another quasi-optical structure. Utilizing a confocal waveguide as the interaction structure for a gyrotron was first proposed at MIT, and experimentally demonstrated later by a 140 GHz fundamental gyrotron oscillator [22] and a 140 GHz gyrotron traveling wave amplifier (gyro-TWA) [23]. As for the harmonic gyrotron, a 0.4 THz gyrotron with a confocal cavity was developed at TRC-UESTC and, experimentally, achieved an output power of 6.44 kW operating at the second cyclotron harmonic [24].

Furthermore, as a type of open structure, a quasi-optical waveguide presents a natural frequency tunability related to the separation distance between two mirrors, which introduces another method for frequency tuning. In this paper, we propose a high-power broadband continuous frequency-tunable gyrotron cavity based on a confocal waveguide. Its frequency tuning characteristics fall under three different strategies, namely magnetic field tuning, mirror separation adjustment, and the hybrid tuning of the above two parameters, which are both investigated by particle-in-cell (PIC) simulation. Results suggest that the proposed quasi-optical cavity is able to generate high power of no less than 20 kW over a smoothly continuous frequency tuning band with an extraordinarily broad bandwidth of 41.55 GHz around 0.22 THz. Compared with other frequency tuning approaches, this method provides advantages in terms of high power, broad band, and smooth continuity.

The paper is organized as follows: the design principles of the quasi-optical cavity, including the introduction of a quasi-optical waveguide, the cavity design, and the cold cavity frequency tunability, are described in Section 2. Section 3 presents the detailed PIC simulation results for three different frequency tuning approaches, covering magnetic tuning, mirror separation adjustment, and hybrid tuning. Finally, a summary for this work is reported in Section 4.

## **2. Cavity Design Principles**

#### *2.1. Quasi-Optical Waveguide*

As shown in Figure 1, the open quasi-optical waveguide is composed of two identical cylindrical mirrors with a finite aperture of 2*a* and a curvature radius of *Rc*. When the separation distance between the two mirrors *L*<sup>⊥</sup> is equal to *Rc*, the two mirrors form a confocal system, which is called a confocal waveguide.

**Figure 1.** Cross-section scheme of the quasi-optical waveguide and the electric field distribution for the TE<sup>06</sup> mode.

Under the geometrical optical approximation for a high frequency wave, the membrane function for the transverse electric (TE) mode in the open waveguide can be obtained by solving the Helmholtz equation in the elliptic coordinate system [24]. The numerically calculated results of the electron field distribution for the TE<sup>06</sup> mode are shown in Figure 1. The eigen-wavenumber *kmn* and the cut-off frequency *fcmn* for the TE*mn* mode can be deduced as

$$k\_{mn} = \frac{\pi}{L\_\perp} \left( n + \frac{2m+1}{\pi} \arcsin\sqrt{\frac{L\_\perp}{2R\_c}} \right) \tag{2}$$

$$f\_{cmn} = \frac{k\_{mn}c}{2\pi} = \frac{c}{2L\_{\perp}} \left( n + \frac{2m+1}{\pi} \arcsin\sqrt{\frac{L\_{\perp}}{2R\_{c}}} \right) \tag{3}$$

According to Equations (2) and (3), the eigenfrequency for the EM mode in a quasioptical waveguide is almost linear with the mirror separation distance *L*⊥. It appears that a smoothly continuous variation of the operating frequency can be easily realized in a quasi-optical gyrotron cavity. By moving the mirrors smoothly, the cut-off frequency for the operating mode will change, leading to the possibility of the continuous tuning of the oscillation frequency in the quasi-optical cavity. This brings a new approach to frequency tunability that cannot be accomplished in traditional gyrotrons based on closed waveguides.

On the other hand, the quasi-optical waveguide provides an impressive mode selection feature. Since lacking sidewalls, as shown in Figure 1, some EM modes will be diffracted out and undergo a large diffraction loss. Previous researchers have demonstrated that only the TE0*<sup>n</sup>* mode could be stably propagated by selecting a small mirror aperture. Thus, the frequency separation between neighboring modes in a quasi-optical waveguide is about ∆*f* = *c*/(2*L*⊥), which is a much greater isolation than that of a cylindrical or coaxial waveguide. The lower mode density provides a distinct advantage for quasi-optical gyrotrons in realizing broadband frequency tuning.

#### *2.2. Frequency Tuning Characteristics in a Cold Cavity*

As shown in Figure 2, a TE<sup>06</sup> mode frequency-tunable quasi-optical cavity for a highpower sub-THz gyrotron is designed and studied in this paper. This cavity is similar to the 0.4 THz second harmonic confocal cavity reported previously [24]. The mirror radius in the straight section is set to 4.20 mm, corresponding to a cut-off frequency of 223.06 GHz for a TE<sup>06</sup> mode under a rigorously confocal situation. The detailed structural parameters of the designed quasi-optical cavity are listed in Table 1.

**Figure 2.** Structural configuration of the quasi-optical gyrotron cavity.

**Table 1.** Structural parameters of the proposed quasi-optical gyrotron cavity.


<sup>1</sup> Under a rigorously confocal condition.

For a cold cavity with the absence of an electron beam, the oscillation frequency *fosc* can be obtained by solving the following differential equation

$$\frac{d^2f(z)}{dz^2} + k\_z^2(z)f(z) = 0\tag{4}$$

satisfying certain boundary conditions at the left end (*z* = *zin*) and the right end (*z* = *zout*).

$$\left. \frac{df(z)}{dz} \right|\_{z=z\_{in}} - jk\_z(z\_{in})f(z\_{in}) = 0 \tag{5}$$

$$\left. \frac{df(z)}{dz} \right|\_{z=z\_{\text{out}}} + jk\_z(z\_{\text{out}})f(z\_{\text{out}}) = 0 \tag{6}$$

Based on the complex oscillation frequency *fosc* and the axial field distribution *f* (*z*), the quality factor *Q* for the metallic quasi-optical cavity can be calculated by

$$\frac{1}{Q} = \frac{1}{Q\_{diff}} + \frac{1}{Q\_{diff\perp}} + \frac{1}{Q\_{ohm}}\tag{7}$$

$$\|Q\_{diff}\| = \frac{\text{Re}(f\_{\text{osc}})}{2\text{Im}(f\_{\text{osc}})},\\\|Q\_{diff}\| = \frac{k\_{\text{min}}L\_{\perp}}{\Lambda},\\\|Q\_{\text{ohm}} = \frac{2}{\delta} \frac{\iiint\_{V} |H|^{2}dV}{\iint\_{S\_{\text{mirror}}} |H\_{\perp}|^{2}dS} \tag{8}$$

where *δ* is the skin depth on the metallic mirror surface *Smirror*, and the diffraction loss parameter Λ can be approximated as

$$\log\_{10} \Lambda = -0.0069 \mathcal{C}\_F^2 - 0.7088 \mathcal{C}\_F + 0.5443, \text{ for TE}\_{0t} \text{ mode, } \mathcal{C}\_F = k\_{mn} a^2 / L\_\perp \tag{9}$$

For a rigorously confocal cavity where *L*<sup>⊥</sup> = *R<sup>c</sup>* = 4.20 mm, the cold cavity characteristics for axial modes TE06*<sup>q</sup>* (*q* = 1, 2, 3, 4) are numerically calculated. The cold cavity oscillation frequency *f*, the *Q*-value, and the frequency band for each TE06*<sup>q</sup>* mode are listed in Table 2, while Figure 3 shows their normalized axial field profiles. Theoretical results suggest that the output frequency of the confocal cavity may be tuned in a frequency range of 4.33 GHz, from 223.20 to 227.53 GHz, if the first four HOAMs can be excited. It is a wider band compared to the previously reported HOAMs gyrotron with a cylindrical cavity, in Reference [16].

**Table 2.** Cold cavity characteristics for the TE06*<sup>q</sup>* modes.


**Figure 3.** Calculated axial field profiles of the TE06*<sup>q</sup>* (*q* = 1, 2, 3, 4) modes in the cold cavity.

Besides, the dependencies of the cold cavity oscillation frequency and the *Q*-value for the TE<sup>061</sup> mode on the mirror separation distance *<sup>L</sup>*<sup>⊥</sup> are calculated and plotted in Figure 4. It should be pointed out that the adjusting of mirror separation *<sup>L</sup>*<sup>⊥</sup> in this work is improved by moving the total upper and lower parts of the cavity with geometrical symmetry in the *y*-direction, rather than just the straight sections reported in References [25,26].

**Figure 4.** Calculated results of the cold cavity oscillation frequency and the *Q*-value for the operating TE<sup>061</sup> mode depending on the mirror separation distance *L*⊥. The adjustment method of *L*<sup>⊥</sup> is schematically shown in the illustration.

As illustrated in Figure 4, the calculation results predict that the cold cavity oscillation frequency can be smoothly varied by adjusting the separation distance *L*⊥. If *L*<sup>⊥</sup> changes from 3.80 to 4.60 mm, within a ±0.40 mm displacement compared to the confocal condition *L*<sup>⊥</sup> = 4.20 mm, the oscillation frequency for the TE<sup>061</sup> mode will be tuned from 246.14 to 204.38 GHz continuously, corresponding to a significant wideband of 41.76 GHz (about 19% around 0.22 THz). At the same time, the cavity's *Q*-value will gradually reduce from 4016 to 2520, but still be retained at a high level. The great frequency characteristics of the cold cavity make it possible to achieve broadband frequency tuning in gyrotron operations.

#### **3. PIC Simulation**

To investigate the output performance of the proposed quasi-optical cavity driven by a gyrating electron beam, especially for the frequency tuning characteristics, the designed model has been built and simulated with the help of a 3D particle-in-cell (PIC) code, CHIPIC [27]. The initial operating parameters of the electron beam are listed in Table 3, and are based on the linear gyrotron theory [24]. During the simulation procedure, the beam velocity spread and the cavity ohmic loss were not taken into consideration.

**Table 3.** Initial operating parameters of the electron beam.


Under the initial beam parameters, the simulation results of the rigorously confocal cavity, where *L*<sup>⊥</sup> = *R<sup>c</sup>* = 4.20 mm, are illustrated in Figure 5. As can be seen, the confocal cavity could generate a stable output power of 21.2 kW at a single frequency of 222.7 GHz, corresponding to an interaction efficiency of 10.6%. The field distributions shown in Figure 5c,d predict that the operating mode is a TE<sup>061</sup> mode, as expected.

**Figure 5.** Typical simulation results of a rigorously confocal cavity: (**a**) instantaneous and average output power; (**b**) spectrum of the output magnetic field component *Bz*; (**c**) transverse distribution of the electric field component *Ex*; (**d**) axial distribution of the magnetic field component *Bz*.

#### *3.1. Magnetic Field Tuning for Confocal Cavity*

The frequency tuning characteristic of the proposed quasi-optical gyrotron cavity is first studied by changing the operating magnetic field, as the tuning method for traditional gyrotrons does. For the mirror separation *<sup>L</sup>*<sup>⊥</sup> <sup>=</sup> *<sup>R</sup><sup>c</sup>* = 4.20 mm in a rigorously confocal cavity, the external magnetic field *B*<sup>0</sup> varying from 8.26 to 9.02 T is simulated when the other beam

parameters are the same as the initial setting listed in Table 3. The oscillation frequency and the output power are as shown in Figure 6. As *B*<sup>0</sup> increases from 8.30 to 8.98 T, the oscillation frequency increases from 222.55 to 226.15 GHz with a frequency tuning range of 3.6 GHz, while the output power changes between 37.2 kW (at 8.32 T) and 1.1 kW (at 8.58 T). On the other hand, the excited transverse mode is the TE<sup>15</sup> mode for *B*<sup>0</sup> at less than 8.30 T, and the TE<sup>16</sup> mode for *B*<sup>0</sup> at larger than 8.98 T, rather than the expected TE<sup>06</sup> mode.

**Figure 6.** Simulation results of the oscillation frequency and the output power depending on the magnetic field *B*<sup>0</sup> for a confocal cavity.

The variation trends of both frequency and power in a confocal cavity are similar to those of the HOAMs excited in a long cylindrical cavity as reported in Reference [28]. In fact, the observation results of axial field distribution suggest that the first four axial modes of TE06*<sup>q</sup>* (*q* = 1, 2, 3, 4) are excited in the confocal cavity during the whole simulation. There are two obvious frequency jump points around 8.73 and 8.85 T. The reason for this is that the operating axial mode varies from one to another in these regions, somewhat affecting the continuity of frequency tuning. However, compared with the performances of conventional gyrotrons operating in this frequency band or output power level, a 3.6 GHz frequency tuning range has provided a great advantage, which should be attributed to the good mode-selective characteristics of a quasi-optical waveguide.

#### *3.2. Mirror Separation Adjustment*

Another kind of frequency tunability for the open quasi-optical cavity related to the mirror separation distance, as introduced in Section 2, is also investigated by PIC simulation.

By parallel shifting the upper and lower parts of the cavity in the *y*-direction, the separation distance between the two mirrors *L*<sup>⊥</sup> is adjusted from 4.10 to 4.60 mm. With the fixed operating parameters listed in Table 3, the simulation results for the oscillation frequency and output power are shown in Figure 7. At a fixed magnetic field *B*<sup>0</sup> = 8.4 T, with an increase in mirror separation distance *<sup>L</sup>*<sup>⊥</sup> from 4.13 mm to 4.56 mm, the oscillation frequency continuously decreases from 227.15 GHz to 210.25 GHz. The output power fluctuates between 29.5 kW (at 4.16 mm) and 1.3 kW (at 4.32 mm). Like in magnetic tuning, the spurious TE<sup>15</sup> and TE<sup>16</sup> modes are also excited at the lower and higher frequency edge, respectively.

**Figure 7.** Simulation results of the oscillation frequency and the output power depending on the mirror separation *<sup>L</sup>*<sup>⊥</sup> at a fixed magnetic field 8.4 T.

Through adjusting mirror separation, the quasi-optical cavity covers a 17.1 GHz smooth frequency tuning band (about a relative bandwidth of 7.8% to the center frequency of 220 GHz), which is much broader than the 3.6 GHz acquired using a magnetic tuning mechanism. Although the continuous transitions of the operating axial mode from TE<sup>061</sup> to TE<sup>064</sup> have still been observed, the frequency tunability in mirror separation adjustment exhibits a much smoother tendency than that of the magnetic tuning shown in Figure 6.

## *3.3. Hybrid Tuning*

As mentioned in Section 1, for the two items of the electron cyclotron frequency Ω<sup>c</sup> and the EM wave frequency *ω*, there is only one variable that is alternated during the above simulations. Due to the mismatch of the cyclotron resonance conditions, the output powers are not stable and present a wide fluctuation in both frequency tuning manners. In this case, to maintain the output power at a relatively stable level, a hybrid tuning approach adjusting both the mirror separation *L*<sup>⊥</sup> and the magnetic field *B*<sup>0</sup> is simulated. During this simulation, under each separation distance *L*⊥, the magnetic field value *B*<sup>0</sup> is scanned to find the peak output power while the other beam parameters are still as the values listed in Table 3.

The simulated dependencies of the peak output power and its corresponding oscillation frequency on the mirror separation distance *L*<sup>⊥</sup> are plotted in Figure 8, which also illustrates the corresponding values of the magnetic field *B*0. It can be found that the output power is able to be kept at a high level, over 20 kW (an interaction efficiency of 10%), by adjusting the magnetic field correspondingly, while the maximum power is about 34.4 kW at *L*<sup>⊥</sup> = 4.22 mm. The oscillation frequency covers an extraordinarily broad bandwidth of 41.55 GHz, from 245.50 to 203.95 GHz, when *L*<sup>⊥</sup> varies within ±0.40 mm around the confocal distance *L*<sup>⊥</sup> = *R<sup>c</sup>* = 4.20 mm. By this tuning method, the frequency tuning bandwidth achieves 18.9% relative to the center frequency of 220 GHz. The frequency tuning characteristics predicted by PIC simulation are consistent with those of cold cavity analysis in Figure 4. Besides this, the axial field distribution results predict that the TE<sup>061</sup> mode has always been the operating mode of every optimized parameter for hybrid tuning. That is why the oscillation frequency in Figure 8 varies almost linearly with the mirror separation distance. Therefore, the frequency tunability of hybrid tuning possesses great characteristics of smooth continuity, wide bandwidth, and high power.

**Figure 8.** Simulation results of the peak output power and the corresponding oscillation frequency depending on the mirror separation *<sup>L</sup>*<sup>⊥</sup> during hybrid tuning. The green labels are the values of the corresponding magnetic field *B*<sup>0</sup> in Tesla.

The frequency tuning range is much wider than the requirements of high-power THz applications at present. In principle, it is possible to achieve a lower output frequency by further increasing the separation distance between two mirrors. PIC simulation result predicts that the quasi-optical cavity could still be operated in the TE<sup>061</sup> mode if *L*<sup>⊥</sup> = 5.20 mm, and generate an output power of 21 kW at 181.05 GHz.

#### **4. Conclusions**

In this paper, the frequency tuning characteristics of a high-power sub-THz quasioptical gyrotron cavity based on a confocal waveguide have been investigated by PIC simulation. The possibility of smooth continuous frequency tuning over a broader band in a high-power sub-THz gyrotron is demonstrated by the simulation results. By continuously varying the mirror separation distance of a quasi-optical cavity and adjusting the operating magnetic field correspondingly, an extraordinarily wide frequency tuning band of 18.9% around 0.22 THz can be realized in a single mode of gyrotron operation, while the output power is kept at more than 20 kW. With the advantages of broad bandwidth, smooth tuning, and high power, this frequency tuning approach can be promoted to other bands and may be beneficial to the development of high-power frequency-tunable THz gyrotrons for some modern emerging applications.

**Author Contributions:** X.G. and W.F. contributed to the overall study design, analysis, computer simulation, and writing of the manuscript. J.Z., D.L., T.Y., Y.Y. and X.Y. provided technical support and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported in part by the National Key Research and Development Program of China No. 2019YFA0210202, in part by the National Natural Science Foundation of China under Grant 61971097 and Grant 6201101342, in part by the Sichuan Science and Technology Program No. 2018HH0136, and in part by the Terahertz Science and Technology Key Laboratory of Sichuan Province Foundation No. THZSC201801.

**Data Availability Statement:** Data sharing not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

