Radiation of a charge moving along a corrugated surface with a relatively small period
Evgeniy S. Simakov, Andrey V. Tyukhtin, Sergey N. Galyamin

TL;DR
This paper investigates electromagnetic radiation emitted by a relativistic charged bunch moving along a corrugated metallic surface, revealing surface wave excitation and absence of volume radiation under certain wavelength conditions.
Contribution
It introduces a method to analyze surface wave excitation and energy losses for a charge moving on a corrugated surface with small period and depth.
Findings
Surface waves are excited by an ultra-relativistic bunch.
Volume radiation is absent for wavelengths larger than the corrugation period.
Spectral density of energy losses is derived and analyzed.
Abstract
We consider electromagnetic radiation of a charged particle bunch moving uniformly along a corrugated planar metallic surface. It is assumed that the wavelengths under consideration are much larger than the period and the depth of corrugation. Using the method of the equivalent boundary conditions we obtain the Fourier-transform of the Hertz vector. It is demonstrated that the ultra-relativistic bunch excites the surface waves, whereas the volume radiation is absent. Fourier-transforms of the surface wave components and spectral density of energy losses are obtained and analyzed.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Radiation of a charge moving along a corrugated surface with a relatively small period
Evgeniy S. Simakov
Andrey V. Tyukhtin
Sergey N. Galyamin
Saint Petersburg State University, 7/9 Universitetskaya nab., St. Petersburg, 199034 Russia
Abstract
We consider electromagnetic radiation of a charged particle bunch moving uniformly along a corrugated planar metallic surface. It is assumed that the wavelengths under consideration are much larger than the period and the depth of corrugation. Using the method of the equivalent boundary conditions we obtain the Fourier-transform of the Hertz vector. It is demonstrated that the ultra-relativistic bunch excites the surface waves, whereas the volume radiation is absent. Fourier-transforms of the surface wave components and spectral density of energy losses are obtained and analyzed.
pacs:
41.60.-m, 41.60.Bq, 41.75.Lx
I Introduction
Electromagnetic interaction of charged particle beams with corrugated structures is of essential interest nowadays for a number of reasons. For example, prospective applications of these structures are connected with microwave and Terahertz (THz) emission. One should mention here a theoretical treatment (based on surface impedance formalism) of bunch wakefields generated in flat and circular metallic waveguides with shallow (fine) corrugation Stupakov and Bane (2012); Bane and Stupakov (2016); Bane et al. (2016); Tyukhtin et al. (2018). Possibilities for microwave and THz radiation in such structures were also studied, including simulations Bane and Stupakov (2012a) and experiments Yalandin et al. (1997); Ginzburg et al. (2016a); Bane et al. (2017); He et al. (2017); Mesyats et al. (2017). Interaction between continuous electron beam and corrugated surface was analyzed in papers Mineo and Paoloni (2010); Ginzburg et al. (2013, 2016b) in view of realization of relativistic surface-wave amplifiers.
Concerning the radiation generated by charged particle bunches moving along periodic structures, one should mention extended researches of various cases where wavelengths are comparable or less than the structure period. This type of radiation is usually referred to as Smith-Purcell radiation, and there are plenty of the literature devoted to this subject (see, for example, books Potylitsyn et al. (2011); Potylitsyn (2011)).
In this paper, we study theoretically the situation, which is principally different from Smith-Purcell case: the wavelengths under consideration are supposed to be much larger than periods of the structure, i.e. we consider so-called “long-wave” (or “low-frequency”) approximation. In this case, the periodic structure can be described by so-called equivalent boundary conditions (EBC), which should be fulfilled on the flat surface (they are also known as Vainstein-Sivov conditions) Nefedov and Sivov (1977).
Contrary to the aforementioned papers dealing with waveguide structures with shallow corrugation Stupakov and Bane (2012); Bane and Stupakov (2016); Bane et al. (2016); Tyukhtin et al. (2018), here we analyze the case of a point charge (or a one-dimensional finite gaussian bunch) moving along infinite flat corrugated surface. This problem has not been solved previously, despite numerous papers dedicated to radiation from charged particle bunches moving in the presence of periodic corrugated structures.
It should be underlined that the paper Tyukhtin et al. (2018) contains the comparison between theoretical results (obtained using the EBC approach) and results of the CST Particle Studio simulations. The authors have demonstrated that the ultra-relativistic Gaussian bunch excites monochromatic wakefield, with the simulated wavelength being equal approximately to mm (for the period of corrugation mm). The difference between theoretical and simulated wavelengths varies from to depending on the width of “hills” of the corrugated structure. The coincidence can be considered as good, and this circumstance justifies the use of the EBC method even in situations where the structure period is only several times shorter than the considered wavelengths.
It is also notable that the problem under consideration is close to similar problems where the periodic structure is not solid. In particular, the “long-wave” radiation from the bunch moving along the planar periodic wire structure has been studied in papers Tyukhtin et al. (2014, 2015). Under this approximation, the structure is described by so-called averaged boundary conditions. Mentioned papers demonstrate that the bunch generates surface waves propagating along the wires, and the surface wave configuration can be used for determination the size of the bunch. However, despite the mentioned similarity, there is a principal difference between EBC and averaged boundary conditions, therefore the problem considered in the present paper is of independent value.
The content of the paper is the following. First, we briefly recall the method of EBC. Next, we obtain the general exact solution of the problem and carry out the asymptotic analysis of this solution. The last section is devoted to the description of the main physical effect, which consists in generation of surface waves.
II Equivalent boundary conditions
We consider a perfectly conductive flat surface having rectangular corrugation (Fig. 1). It is assumed that the period and the depth of corrugation are much less than the wavelength under consideration :
[TABLE]
In this case, we can replace the corrugated surface with a planar surface, on which the so-called equivalent boundary conditions (EBC) are fulfilled Nefedov and Sivov (1977). Note that the EBC are widely applied in the theory of corrugated structures (for instance, in Nefedov and Sivov (1977) one can find a lot of corresponding examples).
The EBC have the following general form for Fourier-transforms of electric and magnetic fields Nefedov and Sivov (1977):
[TABLE]
where and are “impedances”, which are imaginary for perfectly conductive structures. In the case of the structure shown in Fig. 1, we have Nefedov and Sivov (1977):
[TABLE]
where , is a width of groove, and are directional cosines of the incident wave with respect to - and -axis, correspondingly ( is the wave vector of the incident wave, ). The parameter of corrugation is determined by the formula Nefedov and Sivov (1977):
[TABLE]
where parameters and should be found from the following system of transcendent equations:
[TABLE]
Note that the parameter is of the same order as and . In the specific case of a system of thin diaphragms (where ) the parameter is determined by the following simple expression Nefedov and Sivov (1977):
[TABLE]
III General solution of the problem
We assume that the charged particle bunch moves along the corrugated surface with velocity , which is perpendicular to corrugations. The distance between the charge and the surface is (Fig. 1). The charge density is written in the form
[TABLE]
where is the charge distribution along the trajectory of the motion. We will use the Hertz potential and present it as a sum of an “incident” potential and a “reflected” potential :
[TABLE]
where and are corresponding Fourier-transforms. We mean that the incident field is the field in unbounded vacuum, and reflected field is an additional field connected with the influence of the corrugated structure. Relations between the Fourier-transforms of electromagnetic field and Hertz potential are given by the formulae
[TABLE]
It is well known that from Maxwell’s equations one can obtain the Helmholtz equation for the Hertz potential:
[TABLE]
The solution of Eq. (3.4) for the incident field is the well-known Coulomb field of moving charge in unbounded vacuum. The Fourier-transform of the Hertz vector components of this field is
[TABLE]
[TABLE]
where , and . Here is the Fourier-transform of the charge distribution along the trajectory of the motion :
[TABLE]
One can show that the reflected field is described by the two-component Hertz vector
[TABLE]
where
[TABLE]
[TABLE]
Here, and can be called reflection coefficients, which should be found from the boundary conditions.
One can see that, in (III) - (III), the longitudinal component of the wave vector is . Therefore, , and Eq. (2.3) is written in the form
[TABLE]
Satisfying boundary conditions (2.2) one can find the following expressions for the reflection coefficients:
[TABLE]
[TABLE]
where . Thus, we have found the general solution of the problem.
IV Analysis of singularities in the Fourier-integrals
Further, we will perform an asymptotic analysis of the obtained Fourier transforms of the field components. For simplicity, we will consider only the positive frequencies . For negative frequencies, all formulas are easy to obtained by the rule where the asterisk means complex conjugation. This rule is true for Fourier transform of any real function , in particular, for the field components and .
For our consideration, the poles of integrands in (III) and (III) will play an important role. They are found from the dispersion equation
[TABLE]
Taking into account that and are small parameters of the problem, one can find the following approximate solutions of Eq. (4.1):
[TABLE]
As we will see below, if these poles are real then their contributions are surface waves. Since the first term under radical is small then can be real only for ultra-relativistic case (). Approximately, the condition of reality of (4.2) has the form
[TABLE]
The integrands in (III) and (III) have as well other singularities: the function has two branch points . Fig. 2 shows the integration path, the poles , the branch points and the cuts in the complex plane . Taking into account some negligible small losses (i.e., negligible small imaginary part of ) one can show that the poles are located with respect to the integration path as shown in Fig. 2.
To investigate an asymptotic behavior of Fourier-integrals (III) and (III) we can use the saddle point method Felsen and Marcuvitz (2003); Bleistein and Handelsman (1986). First, it is convenient to change the integration variable according to the rule
[TABLE]
Then the phase is written in the following form:
[TABLE]
Saddle points are found from the equation
[TABLE]
It gives two saddle points , which are purely imaginary: . According to (4.4), mapping of the physical sheet of the complex plane is the strip on the complex plane (Fig. 3). There is only one saddle point in this area. The steepest descent path contains the point and can be found from the system
[TABLE]
According to Eqs. (4.7) the steepest descent path is a straight line parallel to the real axis (Fig. 3).
The initial integration path can be transformed to the path with separation of the contribution of the pole for (or for ) that is
[TABLE]
The integrals over the steepest descent path can be asymptotically estimated by the well-known formula (see the books Felsen and Marcuvitz (2003); Bleistein and Handelsman (1986)):
[TABLE]
Formula (IV) demonstrates that the contribution of the saddle point decreases exponentially with increase in and . We can neglect this contribution if the following inequality is satisfied:
[TABLE]
Further, we perform an asymptotic analysis of the field components assuming that and are sufficiently great to satisfy inequality (4.10). In this case, we can neglect the contribution of the saddle point, and Fourier integrals (III) and (III) are determined by the contribution of the pole only (for ; or for ). Underline that exponential decrease of the saddle point contribution means the absence of volume radiation.
The absence of volume radiation can be explained by the following way. The field of the charge on the frequency under consideration must be propagated along the charge trajectory with the phase velocity being equal to the charge velocity (like Cherenkov radiation in the material media). Therefore, at the frequency , we have , and the dispersion equation results in . This equation has no solution if and are real. This predetermines the absence of volume radiation.
It should be underlined that this conclusion refers to the range of relatively low frequencies when the wavelength is much more than the structure period, and we can describe the real structure with help of the EBC. Well-known Smith-Purcell radiation is not generated on these frequencies. Smith-Purcell radiation is an effect of superposition of diffraction radiation from different periods of the structure. The phenomena considered here are not connected with diffraction radiation; they are determined by specific anisotropy of the structure.
V Surface waves
The contribution of the pole for (and for ) is a surface wave under condition that the pole is real, i.e. when inequality (4.3) is satisfied. Accordingly to (4.3), in this case the bunch is ultra-relativistic. Therefore, further we will assume that the parameter is close to unity (). Taking into account this approximation, we can write expressions for poles (4.2) in the following way:
[TABLE]
Calculating the residues we obtain for the Hertz vector of the surface waves
[TABLE]
The corresponding field components are
[TABLE]
[TABLE]
These results are true for . Recall that corresponding formulae for negative frequencies are obtained by the rule .
Note that, in (V), (V), we take into account only the main terms with respect to the small parameter . Note that the components and vanish in the approximation under consideration (we can neglect these components because they are determined by terms of third order of smallness). As we see, all components decrease exponentially with an increase in , i.e. formulae (V) – (V) describe the surface waves propagating in a plane of the structure.
It is also of interest to analyze energy losses of the bunch, which can be obtained by calculation of the energy flow through two parallel half-planes , . This way results in the following expression for the energy losses per unit of the path length:
[TABLE]
where is the -component of the the energy flow density
[TABLE]
The factor “2” in (5.5) is explained by generation of two symmetrical surface waves in positive and negative directions of .
Writing components of electromagnetic field as the Fourier-integrals and substituting these integrals into (5.5), (5.6) we obtain after some transformations the following integral over frequency:
[TABLE]
where is a spectral density of the energy losses
[TABLE]
Substituting formulae (V) and (V) into (5.8) and performing the integration one can obtain the following expression:
[TABLE]
where
[TABLE]
is dimensionless spectral density of the energy losses. Further, we use as well the dimensionless projection of the wave vector:
[TABLE]
We demonstrate below some results of computations performed using the formulas given above. These computations are performed for the conditions providing that the EBC technique is obviously applicable. The small parameter of the problem is from to for the situations considered below. Meanwhile, earlier we have tested the EBC technique for the case of a circular corrugated waveguide with small parameter , and it has been shown that even such moderately small parameter provides the accuracy (compared to CST simulated results) within several percents Tyukhtin et al. (2018).
Fig. 5 shows dependences of the wave vector projection and the spectral density of the energy losses on the depth of corrugations for different periods of the corrugated surface. It is assumed that , , (system of diaphragms). According to the plots, the -component of the wave vector increases with increasing the depth of corrugations , and this dependence is almost linear in the most part of the plot. This is explained by the fact that the parameter only weakly depends on the depth for not small values of (then, accordingly to (5.1), the function is almost linear). Energy losses increase approximately with the square of the depth of corrugations. It should be also emphasized that both the projection of the wave vector and the spectral density of the energy losses decrease with increasing the period .
Fig. 5 shows dependences of the wave vector projection and the spectral density of the energy losses on the depth of corrugations for different values of parameters and (for the same period ). One can see that the -component of the wave vector and the spectral density of the energy losses usually increase with decrease in the width of “hills” , and the system of diaphragms is the most effective for the surface waves generation. However, if the depth of corrugations is small enough then this regularity can be broken.
Concluding this section, we emphasize that the radiation under consideration is fundamentally different from the Smith-Purcell one. Unlike Smith-Purcell radiation, the radiation considered here is surface, has a continuous spectrum and is generated at lower frequencies. However, there are some common features. In particular, the energy of radiation of both types is proportional to the square of the charge and the structure length, and also significantly depends on the geometrical parameters of the cell of the structure.
VI Conclusion
We have investigated radiation of the charged particle bunch moving along the corrugated metallic surface. The principal point here is that we consider radiation at the wavelengths, which are much larger compared to the structure periods (in contrast to traditional consideration for wavelengths which are comparable with the period). The general solution of the problem has been obtained by use of the method of the equivalent boundary conditions. We have performed asymptotic analysis of the electromagnetic field and have shown that the volume radiation is absent in the situation under consideration but the surface waves can be excited. As a rule, these waves are generated if the charge velocity is close to the light velocity. Analysis of the surface waves has demonstrated the main dependences of the energy losses on the problem parameters. In particular, it has been shown that the energy losses increase with increase in the depth of corrugations and, usually, as well with decrease in the width of “hills”.
Acknowledgements.
This work was supported by the Russian Science Foundation (Grant No. 18-72-10137).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Stupakov and Bane (2012) G. Stupakov and K. L. F. Bane, Phys. Rev. ST Accel. Beams 15 , 124401 (2012) . · doi ↗
- 2Bane and Stupakov (2016) K. Bane and G. Stupakov, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 820 , 156 (2016) . · doi ↗
- 3Bane et al. (2016) K. Bane, G. Stupakov, and I. Zagorodnov, Phys. Rev. Accel. Beams 19 , 084401 (2016) . · doi ↗
- 4Tyukhtin et al. (2018) A. V. Tyukhtin, V. V. Vorobev, E. R. Akhmatova, and S. Antipov, Journal of Instrumentation 13 , C 04009 (2018) . · doi ↗
- 5Bane and Stupakov (2012 a) K. L. F. Bane and G. Stupakov, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 677 , 67 (2012 a) . · doi ↗
- 6Yalandin et al. (1997) M. I. Yalandin, S. A. Shunailov, V. G. Shpak, N. S. Ginzburg, I. V. Zotova, A. S. Sergeev, A. D. R. Phelps, A. W. Cross, and P. Aitken, Technical Physics Letters 23 , 948 (1997) . · doi ↗
- 7Ginzburg et al. (2016 a) N. S. Ginzburg, A. M. Malkin, A. S. Sergeev, I. V. Zheleznov, I. V. Zotova, V. Y. Zaslavsky, G. S. Boltachev, K. A. Sharypov, S. A. Shunailov, M. R. Ul’masculov, and M. I. Yalandin, Phys. Rev. Lett. 117 , 204801 (2016 a) . · doi ↗
- 8Bane et al. (2017) K. Bane, G. Stupakov, S. Antipov, M. Fedurin, K. Kusche, C. Swinson, and D. Xiang, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 844 , 121 (2017) . · doi ↗
