Electromagnetic mode profile shaping in waveguides

Electromagnetic mode profile shaping would be a very useful technique, with applications including in accelerator science and data transmission. Two methods are proposed: one varying the permittivity and the other using a wire medium. By the use of Mathieu functions, we demonstrate a flatter, a triangular and a peakier profile. The first two may be used to manipulate longer particle bunches, while the last could be used on shorter bunches and in data transmission. It is also demonstrated how further improvements can be achieved by going beyond Mathieu functions.


Introduction
Electromagnetic modes in a waveguide with uniform cross section have a sinusoidal profile. There are many situations where non-sinusoidal waves may be particularly useful. In accelerator science, for example, one may wish to have a flatter profile. The waveform in Fig. 1 would accelerate a longer bunch, and therefore more electrons, for a given peak power. By contrast a peakier profile, Fig. 2, would accelerate shorter bunches for the same total power. In addition due to the steeper gradient, it would enable one to give a shorter bunch more chirp which is necessary for bunch compression. The peakier profiles would also be useful for data transmission where one can have the higher peak when the signal is being observed. Intermediate between the flatter and peakier profile is a triangular profile as given in Fig 3. This has the advantage of having a significant part of the profile having a constant slope. In terms of interactions with electron bunches, this would give longer bunches a sustained chirp. Figures 1, 2 and 3 are all examples of waveforms based on Mathieu functions, which often have convenient properties. However, we can go beyond Mathieu functions. In Figs. 4 and 5, we demonstrate how an even flatter profile can be achieved.
In this article we modify the mode profile by modifying the constitutive relations, in particular the permittivity. To simplify the analysis, we consider the permittivity to be periodic in z the direction of propagation.
We consider two scenarios. The first is for transverse waves, Sect. 2. In this case, we choose a permittivity which depends on frequency and position [2]. The use of varying the permittivity combined with the use of Mathieu functions was considered by El Haddad [6], where the transmission functions were calculated. However, the author did not consider the use of this to modify the field profile. There are many methods for making such a medium. For example, one can compose the medium from many thin slices of media, each with a different value of , thereby approximating the continuous variation. Another method is to set up a standing wave in a nonlinear material. The linearised permittivity about this solution would have the spatial dependence as desired. However, for the examples in Figs. 1 and 4 given in this article to work would require regions to have negative permittivity, \0 which would require some kind of metamaterial [3,[14][15][16][17]19]. It also does not work for longitudinal waves. An alternative to prescribing the permittivity profile and asking if the resulting electric field profile is desirable is to start with a particular field profile and attempt to see whether it can be constructed with a realistic permittivity. This is investigated in subsequent work [7].
A second scenario, Sect. 3, is to make spatially dispersive as well as inhomogeneous. That is, we set ¼ ðx; k; zÞ where k is the Fourier conjugate variable associated with z. It is natural to ask the meaning of a function depending on both k and z. This is addressed in [2,8,10,11]. The method, as described below, is to replace the permittivity function with a differential equation. A wire medium [1,4,9,22,23] is naturally spatially dispersive, and by varying the radius of the wires periodically in z will naturally give it the correct spatial inhomogeneity (Fig. 12).

Transverse modes periodic inhomogeneous permittivity
Consider the single-frequency transverse mode with E ¼ e ixtẼ ðzÞi, P ¼ e ixtP ðzÞi and H ¼ e ixtH ðzÞj, together with the permittivity ðx; zÞ ¼ 0 r ðx; zÞ and vacuum permeability l 0 . Here fi; j; kg are the unit vectors along fx; y; zg. Then Maxwell's equations givẽ where L is the period of the variation: ðx; zÞ ¼ ðx; z þ LÞ. Then we get the (rescaled) Mathieu's equations It is easy to model spatially varying permittivity numerically using the CST package [5], when [ 0. Since a ¼ A 1 ðqÞ where A n ðqÞ is the Mathieu-A function, then from (2) we require that A 1 ðqÞ [ q j j. This implies that Conveniently, the two limits of this range correspond to the triangular and flatter profiles. To validate these results, numerical simulations of these profiles have been performed using 20 slabs of constant r by the eigen mode solver of CST Microwave Studio [5]. Periodic, PEC and PMC boundary conditions are chosen along the z, x and y directions, respectively. Figures. 6 and 7 show the comparison between the E-field profile for the modelled and simulated modes, where a very good agreement can be observed. The simulated E-field profile on the x À y plane is given in Figs. 8 and 9 for the flatter and triangular modes, respectively. These were implemented directly into CST, by using 20 slabs of constant r . To confirm these modes, we have run a CST simulation as shown in Figs. 6, 7, 8 and 9.
If we are allowed to have negative r , then one can create more extreme shapes. For example, in Figs. 4 and 5 we add a higher harmonic to produce an even flatter profile, whereas in Fig. 2 we obtain a much peakier solution. One must be careful, however, with simulating these. We have used MPB [13] to simulate the fields for q ¼ 0:8 to q ¼ 1, Fig. 10. We see that the region with negative r has numerical instabilities.
The oscillation of the permittivity at the order of the wavelength implies that the structure is similar to a onedimensional photonic crystal or Bragg mirror [12,18,20,21,24]. One can ask what the transmission coefficient is for waves coming into this medium from a vacuum. In Fig. 11 we use CST to calculate the transmission coefficients for the flatter profile as in Fig. 8. We see that at the proposed frequency of 2.9GHz there is a high transmission and we are at the edge of a band gap.  Fig. 6 Comparison of the simulated E x -field (blue) [Fig. 8] for a flatter profile, as in Fig. 1, with the predicted curve (black). However, the two curves are so close that one cannot distinguish them in this figure. In comparison with the cosine (dashed brown) Fig. 7 Comparison of the simulated E x -field (blue) [ Fig. 9] for a triangular profile, as in Fig. 3, with the predicted curve (black). Likewise, the two curves are almost identical. In comparison with the cosine (dashed brown) i.e. D ¼ 0, thus we are looking for epsilon near zero (ENZ) media. When the medium is homogeneous, we will use an empirical model of the permittivity via Combining (5) and (6), we obtain the dispersion relation Taking the Fourier transform of 0 ¼ 0Ê ðkÞ þPðkÞ ¼ ðx; kÞÊðkÞ with respect to k usingPðzÞ ¼ R 1 À1 e ikzP ðkÞdk, one obtains the differential equation   Transmission coefficients for a periodically varying permittivity as in Fig. 8. This is implemented in CST using 10 cells. The transmission is very high and corresponds to the right edge of the band gap at around 2.89 GHz (red) The simplest method to include an inhomogeneity in the permittivity is to let the plasma frequency k p depend on position z, that is k p ¼ k p ðzÞ. Thus (8) becomes Again by choosing the appropriate periodic function for k p ðzÞ, one can replace (9) by the Mathieu equation.
As stated in the introduction, this can be created with a wire medium, as described in more detail in [9].

Conclusion
Two methods of mode profile shaping are suggested. One using transverse waves and a varying permittivity. The second is by using a wire medium with periodic variation in the radius of the wire. In the first method, we observed that for the more extreme profiles such as the peakier profile, Fig. 2, and the flatter profile Fig. 4 require \0. This is challenging for the numerical solvers CST and MPB. In addition, one would suspect it would also be challenging to implement experimentally. Since we have implemented a 1-dimensional photonic crystal, we see from Fig.11 that there are band gaps and the operating frequency is at the edge of this band gap. Fortunately the transmission is high.
We are currently implementing numerically the wire medium to optimise the shape of the wires to form the desired profile shape.