1 Introduction

Helical gyro-TWTs are considered as promising broadband high-power amplifiers for future microwave applications, like advanced radar sensors for space debris observation in low earth orbit (LEO). Due to the enormous progress in the field of high-power microwave technology [1, 2], corresponding radar sensors can also be operated in high frequency bands such as the W-band [3, 4]. Due to the high bandwidth available there, very high resolutions can be achieved [3,4,5]. In near future, even W-band transmission powers in the range of 10 kW to 100 kW may be possible [6]. To realize a W-band radar sensor with such a high transmission power, a suitable high-power amplifier as well as an appropriate transmission line is required. The transmission line connects the output of the high-power amplifier with the antenna feed. Due to the high power, oversized transmission lines are required. A suitable transmission mode is the HE11 hybrid mode, propagating in a corrugated circular waveguide. Low ohmic loss, small mode conversion, high degree of polarization and low-sidelobe radiation patterns can be achieved [7,8,9]. For transmission, the amplified microwave signal has to be converted into the HE11 hybrid mode, after leaving them interaction circuit. In a simple way, this may be realized due to a corrugated horn antenna [10, 11], illuminating the aperture of a highly overmoded circular waveguide [12, 13]. However, since the broadband frequency behaviour of the transmission line components may be sensitive in respect to the mode content [14, 15], a mode purity as best as possible is required, within the considered frequency range. A concept based on a corrugated horn antenna is limited in respect of the mode purity [12, 13]. Therefore, a waveguide mode converter is preferred here, which may attain a significant better mode purity [12]. The interaction mode for the helical gyro-TWT under consideration can be easily transferred to the circular waveguide TE11 mode [16]. Therefore, a TE11 ↦ HE11 mode converter is addressed here. The frequency range from 92 GHz to 100 GHz is considered.

A corresponding mode converter may be realized employing a customized corrugation profile with varying groove depth within a corrugated waveguide [17, 18]. However, due to the small wavelength in the W-band and the required tiny geometry [17, 18], the technical realization of this type of mode transducer may be challenging. Therefore, a smooth-wall mode converter with customized diameter contour is discussed here. Its technical realization will be much simpler.

A similar mode converter, also for a high-power amplifier tube, was considered in [19]. However, there, a horn antenna concept to couple into the TEM00 free-space mode and another approach for optimization was chosen. To the best knowledge of the authors, a similar design procedure as addressed in the present paper has not yet been published elsewhere.

The paper is organized as follows: In Section 2 fundamental aspects are introduced, important for the following sections. Section 3 discusses the proposed design procedure and corresponding simulation results. Section 4 closes with a conclusion and an outlook for future research activities.

2 Fundamentals

To clarify the context of the present paper, Subsection 2.1 shortly discusses fundamental aspects of a helical gyro-TWT. The purpose of use and the location of the addressed mode converter within such a vacuum electron device is emphasized. A comprehensive introduction of helical gyro-TWTs can be found in [20,21,22,23,24,25,26,27]. In Subsection 2.2 calculations on hybrid modes in corrugated circular waveguides are addressed. In Subsection 2.3 equations for the decomposition of an arbitrary electromagnetic wave into the eigenmodes of a corrugated waveguide are given.

2.1 Helical Gyro-TWT

Amplification within a gyro-TWT is based on the electron cyclotron instability [24, 25]. Within the Brillouin diagram [28], those simplified resonance condition represents a straight line [24]. Electron-wave interaction occurs at points of intersection with the interaction mode [21,22,23,24]. A tangential intersection leads to broadband frequency operation [16, 21,22,23,24].

In a helical gyro-TWT, a customized helical interaction waveguide is used to generate an interaction mode which is tangential at close to zero axial wavenumber [1, 2, 20,21,22,23] with a close to constant group velocity [21] within the frequency band of interest. This leads to reduced sensitivity in respect of axial electron velocity spread, broadband frequency behaviour and improved stability [20,21,22,23]. For electron-wave interaction within a three-fold helical waveguide, typically a circularly polarized “quasi” TE21 mode [6] is used. Its transversal field components are similar to the TE21 mode within a circular waveguide, but the dispersion is tailored for broadband electron wave synchronization. After the electron-beam-wave interaction, the “quasi” TE21 mode can be easily transferred to the fundamental circular waveguide TE11 mode [16]. However, for low-loss transmission through to the vacuum window, the HE11 mode is more suitable than the TE11 mode. Therefore, an additional broadband high-power TE11 ↦ HE11 mode converter is required.

As already mentioned, a smooth-wall mode converter with customized diameter contour is addressed here. Besides a simpler technical realization, this also allows to integrate an electron collector within the mode converter. This further simplifies the tube design. However, the corresponding collector design is out of scope of the present paper and will be addressed in more details in an upcoming paper.

At the mode converter output, a corrugated circular waveguide, able to propagate the HE11 hybrid mode under the balanced hybrid mode condition [10], is flanged. No critical reflection or mode conversion has to be expected thereby. The vacuum window will be integrated within the corrugated waveguide.

2.2 Hybrid Modes in Corrugated Circular Waveguides

Figure 1 shows a simplified sketch of a corrugated circular waveguide with the radius a, the corrugation period c, the slot width b and the corrugation depth d.

Fig. 1
figure 1

Corrugated waveguide with the waveguide radius a, the corrugation period c, the slot width b and the corrugation depth d

The electromagnetic field components for r ≤ a may be written as:

$$ {\displaystyle \begin{array}{c}\overrightarrow{E}={\overrightarrow{E}}^{\mathrm{TE}}+\psi \cdotp {\overrightarrow{E}}^{\mathrm{TM}}\\ {}\overrightarrow{H}={\overrightarrow{H}}^{\mathrm{TE}}+\psi \cdotp {\overrightarrow{H}}^{\mathrm{TM}}\end{array}} $$
(1)

with:

$$ {\displaystyle \begin{array}{c}{H}_z^{\mathrm{TE}}={A}^{\mathrm{TE}}\cdotp {J}_m\left(\frac{X}{a}\cdotp r\right)\cdotp \exp \left[- j\omega t\right]\cdotp \left(\genfrac{}{}{0pt}{}{\sin \left( m\varphi \right)}{\cos \left( m\varphi \right)}\right)\\ {}{E}_z^{\mathrm{TM}}={A}^{\mathrm{TM}}\cdotp {J}_m\left(\frac{X}{a}\cdotp r\right)\cdotp \exp \left[- j\omega t\right]\cdotp \left(\genfrac{}{}{0pt}{}{\cos \left( m\varphi \right)}{\sin \left( m\varphi \right)}\right)\end{array}} $$
(2)

and ψ = ATM/ATE. The function Jm(x) is the Bessel function of the first kind and the order m. The function \( {J}_m^{\prime }(x) \) is its first derivation in respect of x. Using the Maxwell equations [28], all remaining field components in Eq. (1) can be calculated based on Eq. (2) [28]. To avoid Bragg reflection, the angular field component at r = a has to vanish [10, 11]. This leads to:

$$ \psi =\frac{m\overset{\sim }{\beta }}{X{Z}_0}\cdotp \frac{J_m(X)}{J_m^{\prime }(X)} $$
(3)

with m as the azimuthal mode index, \( \overset{\sim }{\beta } \) as the normalized phase constant with \( \overset{\sim }{\beta }=\beta /k \), k = 2π/λ as the wavenumber and Z0 = 120π Ω as the free space wave impedance. Ohmic losses are neglected here.

In the following, the eigenmodes are always normalized to Pz = 1 W, with Pz as the mode power in the direction of propagation. For this purpose, the pointing vector [28] is integrated over the waveguide cross-section:

$$ {P}_z=\frac{1}{2}\cdotp {\int}_A\left(\overrightarrow{E}\times {\overrightarrow{H}}^{\ast}\right)\times \overrightarrow{z}\kern0.5em \mathrm{d}\overrightarrow{A} $$
(4)

with \( \overrightarrow{z}={\left(0,0,1\right)}^{\mathrm{T}} \) as the normal vector in the direction of propagation. Then, the field components in Eq. (1) are normalized in respect to \( \sqrt{P_z} \).

The eigenvalues within a corrugated circular waveguide can be calculated as discussed in detail in [10, 11, 29, 30]. For the HE11 hybrid mode, fulfilling the balanced hybrid mode condition, follows: X = 2.404826.

2.3 Mode decomposition

An arbitrary electromagnetic field shall be represented as:

$$ {\displaystyle \begin{array}{c}\overrightarrow{E}=\sum \limits_{\infty}\left(\ {\vartheta}_n^{+}+{\vartheta}_n^{-}\ \right)\cdotp {\overrightarrow{e}}_n\\ {}\overrightarrow{H}=\sum \limits_{\infty}\left(\ {\vartheta}_n^{+}-{\vartheta}_n^{-}\ \right)\cdotp {\overrightarrow{h}}_n\end{array}} $$
(5)

with \( {\vartheta}_n^{+/-} \) as the nth complex mode amplitude in forward/backward direction and \( {\overrightarrow{e}}_n/{\overrightarrow{h}}_n \) as the electromagnetic field components of the nth corresponding eigenmode. As introduced in section 2.2, all eigenmodes are normalized to Pz = 1 W, here. With fundamental manipulations and the power orthogonality condition, \( {\int}_A\left({\overrightarrow{e}}_{\perp n}\cdotp {\overrightarrow{e}}_{\perp m}^{\ast}\right)\ \mathrm{d}\overrightarrow{A}=0 \) for n ≠ m [31, 32], follows:

$$ {\vartheta}_m^{+}=\frac{1}{4}\cdotp \left(\ \frac{1}{Z_m\ }\cdotp {\int}_A\left(\ {\overrightarrow{E}}_{\perp}\cdotp {\overrightarrow{e}}_{\perp m}^{\ast }\ \right)\ \mathrm{d}\overrightarrow{A}+{\mathrm{Z}}_m\cdotp {\int}_A\left(\ {\overrightarrow{H}}_{\perp}\cdotp {\overrightarrow{h}}_{\perp m}^{\ast }\ \right)\ \mathrm{d}\overrightarrow{A}\ \right) $$
(6)
$$ {\vartheta}_m^{-}=\frac{1}{4}\cdotp \left(\ \frac{1}{Z_m\ }\cdotp {\int}_A\left(\ {\overrightarrow{E}}_{\perp}\cdotp {\overrightarrow{e}}_{\perp m}^{\ast }\ \right)\ \mathrm{d}\overrightarrow{A}-{\mathrm{Z}}_m\cdotp {\int}_A\left(\ {\overrightarrow{H}}_{\perp}\cdotp {\overrightarrow{h}}_{\perp m}^{\ast }\ \right)\ \mathrm{d}\overrightarrow{A}\kern0.5em \right) $$
(7)

with:

$$ {\mathrm{Z}}_m=\frac{1}{2}\cdotp {\int}_A{\left|{\overrightarrow{e}}_{\perp m}\right|}^2\ \mathrm{d}\overrightarrow{A}={\left(\frac{1}{2}\cdotp {\int}_A{\left|{\overrightarrow{h}}_{\perp m}\right|}^2\ \mathrm{d}\overrightarrow{A}\right)}^{-1} $$
(8)

as the normalized mode impedance of the corresponding eigenmodes. The integrals are executed in respect of the waveguide cross-section (for r ≤ a).

Based on Eqs. (6) and (7), an arbitrary electromagnetic field can be decomposed in corresponding forward/backward travelling eigenmodes. Just the transversal field components are required, thereby.

3 Design procedure

In this section the design procedure for the broadband high-power TE11 ↦ HE11 mode converter is addressed. The Subsection 3.1 introduces the principle design approach. Corresponding simulation results are addressed in the Subsection 3.2.

3.1 Principle approach

Figure 2 shows the principle flow diagram of the proposed design procedure. The kernel is a genetic algorithm [33] used to solve a global optimization problem. The optimization problem is defined due to a cost-function with the scalar output parameter C. The principle concept of a genetic algorithm is inspired by a natural selection process [33,34,35]: Based on an initial population of individuals, an evolution is set in motion. This is an iterative process, where each iteration leads to a new generation. Within each iteration the cost function for each individual is carried out and evaluated. Due to a stochastic process, the most suitable individuals are selected and used as the starting point for the upcoming generation. For the next generation, the selected individuals are reproduced and mutated (e.g. randomly). Then, again, for each individual, the cost function is carried out, and the most suitable individuals are selected stochastically for the next generation. The optimization terminates, if convergence is detected. Besides a prober initial population, the genetic algorithm is initialized with prober boundary conditions. Here, these are the waveguide diameters at the input/output and a permitted range of the mode converter length L. A more detailed consideration of the genetic algorithm is out of scope of the present paper. A comprehensive consideration is found in [33,34,35]. Here, the MATLAB implementation discussed in detail in [36] is used.

Fig. 2
figure 2

Principle flow diagram of the design procedure.

The parameters, optimized due to a genetic algorithms, are named chromosomes [33]. Here, the chromosomes represent the geometry of the considered mode converter. However, to reduce the required amount of calculation, the chromosomes are limited to a few net points, here. The complete geometry is calculated based on an interpolation within the cost-function.

The cost-function parameter C is calculated as:

$$ 2C=1-\left({\left(\ {\left|{\vartheta}^{+}\right|}^2\ \right)}_{z=L}-{\left(\ \sum \limits_m\ {\left|{\Gamma}_m\right|}^2\ \right)}_{z=0}\ \right) $$
(9)

with ϑ+ as the complex HE11 mode amplitude in forward direction at the mode converter output (z = L) and Γm as the complex reflection coefficient of the mth TE/TM eigenmode at the input of the mode converter (z = 0).

Within the optimization, Eq. (9) is carried out at a single frequency. In principle, also a multifrequency cost-function could be considered. However, the design procedure based on Eq. (9) showed inherent broadband frequency behaviour for corresponding mode converters. Therefore, to reduce the required amount of calculation, a multifrequency cost function was waived. Equation (9) is carried out at the centre frequency of 96 GHz.

Figure 3 shows the principle calculation flow of the cost-function parameter C. The flow chart represents the grey circle in Fig. 2.

Fig. 3
figure 3

Principle flow chart of the cost function calculations

The input parameters are the chromosomes, used for optimization. As mentioned, those represent single net points of the mode converter geometry. The complete geometry is calculated based on an interpolation, which is carried out in the first step in Fig. 3. Based on the mode converter geometry, the corresponding scattering matrix [28] of the structure is calculated. The scattering matrix links the incoming and outcoming electromagnetic waves of the structure. To consider multimode devices, the scattering matrix can be extended to the generalized scattering matrix [37,38,39]. For corresponding calculations, a KIT implementation [39] of the scattering matrix formalism, based on [37, 38], was extended and modified for the present application. As output, the reflection and transmission coefficients of the considered eigenmodes as well as their coupling within the structure are calculated. A more detailed consideration of such a scattering matrix formalism is out of scope of the present paper. A comprehensive consideration is found in [37,38,39]. The utilized scattering matrix code was verified in [39] with multiple simulation structures and different full wave simulation tools (e.g. CST Microwave Studio [40]).

To calculate the cost-function parameter C, the single reflection coefficients at the mode converter input as well as the HE11 mode content at the mode converter output are required. The reflection coefficients are directly provided by the scattering matrix formalism [37, 39]. All of the considered reflection coefficients at the input (z = 0) are summed up to get the overall power reflection at the mode converter input. In ideal applies: ∑|Γm|2 = 0 at z = 0, which implies no reflected power at the mode converter input.

To calculate the HE11 mode content, first the electromagnetic field components at the mode converter output are calculated as superposition of the calculated TE/TM mode content, provided by the scattering matrix formalism [37, 39]. Using the equations from subsection 2.3 and the derived HE11 eigenmodes from subsection 2.2, the required HE11 mode content can be calculated. In ideal applies: |ϑ+|2 = 1 at z = L, which implies a pure HE11 mode at the mode converter output.

Based on the calculations on the reflection coefficients and the HE11 mode content, Eq. (9) can be carried out and used for optimization. In ideal applies: C = 0, which implies that the whole TE11 mode power at the mode converter input is coupled into the HE11 mode at the mode converter output.

3.2 Results

Figure 4 shows the geometry of the optimized TE11 ↦ HE11 mode converter. The red circles represent the chromosomes, used for optimization. Note that the geometry is not described due to an analytical function but corresponds to numerical optimized data points.

Fig. 4
figure 4

Radius profile along the optimized TE11 ↦ HE11 mode converter. The red circles represent the single chromosomes, used for optimization

The specified input diameter is D = 2.749 mm which corresponds to a standard waveguide flange [41] which fits to an optimized helical interaction waveguide. The specified output diameter is D = 31.75 mm which corresponds to a standard corrugated waveguide diameter [42]. The position and the value of the first net point were fixed (D1 = 2.749 mm, z1 = 0) as well as the value of the last net point (D6 = 31.75 mm). Therefore, five net point positions and four net point values had to be optimized. The mode converter length was optimized to L = 39.81 cm with z6 = L. Within the cost function, the geometry is considered with 1000 steps, which ensures convergence in respect of a smooth waveguide profile [39].

Figure 5 shows the corresponding simulations results in respect of the HE11 mode content at the mode converter output, within the frequency range from 90 GHz to 100 GHz. The squared absolute value (|ϑ+|2, colored in brown) as well as the phase history (angle(ϑ+), colored in blue) is shown. Due to the applied normalization, the squared absolute value of the complex mode amplitude corresponds directly to the relative mode content.

Fig. 5
figure 5

Complex HE11 mode amplitude at the mode converter output. Phase response (blue circles). Squared absolute value (brown curve)

The calculated mode content shows the broadband frequency behaviour of the design. Within the frequency range from 92 GHz to 100 GHz, the HE11 mode content is always higher than 98.6 %. Inside this frequency band, the mode content is even better and reaches ≈ 99.7 % at ≈95 GHz.

For broadband frequency operation, a high mode content as well as a linear phase response is essential. A non-linear phase response would lead to a varying group velocity [28] which may cause critical distortions. In respect of a broadband radar sensor, this may impair the geometrical radar resolution [43, 44]. However, Fig. 5 shows a close to ideal linear phase response of the designed mode converter and therefore, again, broadband frequency behaviour. No critical distortion for broadband frequency operation has to be expected.

For the shown mode converter geometry, the optimization procedure leads to a cost-function parameter C= 4.2041e-3. The overall reflection coefficient at the mode converter input, calculated in respect of all considered eigenmodes, is far below −30 dB and can be neglected. In total, 122 TE/TM modes were considered within the scattering matrix formalism. Due to the symmetry within the circular structure and the TE11 input mode, the calculations could be limited to the TE1n/ TM1n modes [37,38,39] with n = 1…61.

Note that the present paper focused on the elaborated design procedure for broadband high-power TE11 ↦ HE11 mode converters. The presented mode converter design will be experimentally proven and discussed in more details later, when the mode converter is manufactured as well as the corresponding test setup is built up. Due to the direct coupling into a corrugated waveguide, simple far-field radiation pattern measurements will not lead to adequate results. An advanced test setup has to be build up, which will be addressed in an upcoming paper as well as corresponding measurements components like advanced broadband mode coupler for precise mode content determination within a corrugated waveguide.

4 Conclusion

The present paper addresses the design procedure for a broadband high-power TE11 ↦ HE11 mode converter. The mode converter is realized within a smooth circular waveguide which simplifies the manufacturing and allows the integration of an electron collector within the mode converter. The mode converter is designed for a W-band (92 GHz to 100 GHz) helical gyro-TWT which is intended to be used for an advanced broadband high-power radar sensor for space debris observation in the low earth orbit (LEO). The design procedure is based on numerical optimization using a genetic algorithm and the scattering matrix formalism.

Within the considered frequency band, a high HE11 mode content as well as a linear phase response could be achieved. This is essential for broadband frequency operation.

Future research activities will address the manufacturing of the mode converter as well as the design and the construction of a prober test system at KIT for such mode converters. Furthermore, more detailed simulation in respect of the integration of an electron collector will be carried out.