Statistics of the radiating relativistic electrons

# Statistics of the radiating relativistic electrons

Eugene Bulyak    Nikolay Shul’ga [ National Science Center ‘Kharkov Institute of Physics and Technology’, 1 Academichna str, Kharkiv, Ukraine
July 20, 2019
###### Abstract

We developed a general method for assessing the evolution of the energy spectrum of relativistic charged particles that have undergone small quantum losses, such as the ionization losses when the electrons pass through matter and the radiation losses in the periodic fields. These processes are characterized by a small magnitude of the recoil quantum as compared with the particle’s initial energy. We convey the statistical consideration of the radiating electrons. We demonstrate that in the initial stage, the electron’s spectrum can be described as a composition of consecutive convolutions of the recoil spectrum with itself, taken with the Poisson mass. In this stage, at the relatively small average number of the recoils, the electron spectrum reveals some individual characteristics of the recoil spectrum. Furthermore, the spectrum loses individuality and allows for proximate description in the terms of statistical parameters. The results of this consideration reveal that the width of the electron spectrum is increasing with the number of recoils according to the power law, with the power index being inverse to the stability parameter, which gradually increases with the number of recoils from one to two. Increase of the spectrum width limits the ability of the beam to generate coherent radiation in hard x-ray and gamma-ray region.

###### pacs:
41.60.-m, 41.75.Ht

Also at ]V.N. Karazin National University, 4 Svodody sq., Kharkiv, Ukraine

## Introduction

In a number of processes involving beams of high-energy electrons, such as radiation in the periodic structures, ionization losses in matter, etc., the energy degradation of an incident electron is in the form of the small portions (recoil quanta), which spectrum is almost independent of the electron’s energy. In our previous papers Bulyak and Shul’ga (2017, 2018), we have considered the evolution of the spectrum of such electron beam. It was shown, that the spectrum is determined by the parameters of a single recoil and the average number of the recoils.

We convey the statistical consideration of the radiating electrons. We demonstrate that in the initial stage, the electron’s spectrum can be described as a composition of consecutive convolutions of the recoil spectrum with itself, taken with the Poisson mass. In this stage, at the relatively small average number of the recoils, the electron spectrum reveals some individual characteristics of the recoil spectrum. Furthermore, the spectrum loses individuality and allows for proximate description in the terms of statistical parameters.

This paper is organized as follows. In the first section, we present a method of assessing the evolution of the straggling function that describes the distribution of the energy losses in the interim range of the number of recoils. In the second section, we validate the method by comparing it to the known theories at the limiting cases. The third section presents the results of the study of the kinetics of the radiating electrons in short undulators. The fourth section summarizes the results.

## I Statistics of the Radiating Electrons

### i.1 Preliminaries

Distinguishing features of the considered system are: (i) Big number of the ensemble members (electrons in the radiating bunch) . (ii) Small average number of recoils (defined as ratio of the energy emitted by the electron to the mean energy of the spectrum of the radiation) .

We adopt the assumption that the spectrum of emitting quantum of the radiation inducing energy loss (recoil), is ‘physical’: it has compact support, with being the energy of the recoil. The spectrum is normalized to unity, . (Here and below we drop out the infinite limits in integration.)

In this paper, we use the reduced energy units: for the energy in the straggling spectrum, for the energy of the spectrum of the recoil quantum, all are dimensionless, normalized to the energy unit, e.g., to the charged particle rest energy, Bulyak and Shul’ga (2018). We use a convention for the Fourier transform in the form of:

 (Ff)(s)=^f(s)=∫e−2πωsf(ω)dω

with being the variable in the Fourier transform domain that complements to (or ).

A trajectory of the straggling electron in the plane is presented in Fig. 1. The trajectory is composed of the free paths of random length, with the mean unit value and the (positive) random jumps having the same probability density distribution . Such a process belongs to the subclass of the Subordinate to the Compound Poisson Process, which in turn belongs to the -stable (or Lévy) processes, see, e.g., Applebaum (2004).

As we have shown, the evolution of the bunch’s spectrum is described by the characteristic function (Fourier transform of the distribution density), Bulyak and Shul’ga (2017, 2018):

 ^f=^f0exp[x(ˇw−1)], (1)

where is the recoil’s spectrum, and is the initial spectrum. The parameter is the ensemble average number of the recoils undergone by an electron since entering the driving force, Bulyak and Shul’ga (2017, 2018).

The equation (1) may be generalized and simplified due to the model assumption of the independence of recoils on the electron’s energy, as proposed in Landau (1944). Instead of the beam spectrum, we consider the distribution density of losses: so called straggling function Bulyak and Shul’ga (2018). The straggling function, normalized to unity, presents the loss spectrum: only the particles that have undergone at least one recoil contribute to it.

The characteristic function for the straggling function and the Poisson-weighted expansion are

 ^Sx =^wex(^w−1), Sx(ϵ) =∞∑n=0e−xxnΓ(n+1)Fn(ϵ), (2) Fn =Fn−1∗w, F0 =w.

where stands for the convolution operation. The first three moments of the straggling function read:

 ¯ϵ =(1+x)¯¯¯ω; (3a) Var[ϵ] =(1+x)¯¯¯¯¯¯ω2−¯¯¯ω2; (3b) Sk[ϵ] =(1+x)¯¯¯¯¯¯ω3−3¯¯¯¯¯¯ω2¯ω+2¯¯¯ω3. (3c)

Here are the raw moments of the recoil spectrum , .

A universal solution for the straggling function (2) allows for accurate evaluation at the beginning of the process, when the series may be limited to a few self-states ; and in the opposite limit of the large number of recoils, when a few first moments (3) adequately represent the function. (The first moment—mean loss energy—always holds since it presents the energy conservation law.)

### i.2 Statistical properties for finite number of recoils

From the practical point of view, most interesting for the physically realizable systems is a medium number of recoils, when a particle, after entering the system, lost a small fraction of its initial energy in the moderate number of the recoils, . To evaluate the functional dependence of the spectrum width against the average number of recoils, , we compare the distribution (2) with the Lévy -stable distributions, the only ones to which the sum of independent identically distributed variables is attracted, see Nolan (2017).

It should be emphasized that despite the relevance of the trajectory to the -stable class, the bunch of such trajectories do not rigorously match the class since individual trajectories come into stage at different , as is depicted in Fig. 2. Nevertheless, at when almost all of the particles have been recoiled and the width of the distribution exceeds the width of the spectrum of recoils, the bunch of the trajectories is expected to obey a stable law.

The characteristic function of -stable process, see e.g. Nolan (2017), has a general form:

 ^ϕ(s)=exp{−2πisμ−|2πσs|α[1−iβsign(s)Φ]}, (4)

with

 Φ=⎧⎨⎩tan(πα2),α≠1,−2πlog|s|,α=1,

The parameters of the stable distribution are the stability parameter, the scale parameter, the skewness parameter, and the location parameter.

The model under consideration allows for a reduction of range of the parameters: the stability parameter should be in the range due to a finite mean of the recoil spectrum; the location parameter is simply equal to the first moment of the straggling function (3a).

Comparison of the characteristic function of straggling (2) with that of the stable distribution (4) leads to two important consequences: (i) the scaling parameter is determined by the real part of the exponent, and (ii) the Fourier transform of a recoil spectrum in general may not be of the power form, with .

We suggest to evaluate the stability parameter at where the real part of the exponent of the characteristic function equals to unity:

 R[x(1−^w(s∗))]=1=|πσs∗|α. (5)

Here in the square brackets of left-hand side of (5), we intentionally omit the small term .

From this suggestion, an expression for the stability parameter is readily derived:

 α=sDsR[^w]1−R[^w]∣∣∣s=s∗, (6)

where and is the root of (5).

Substituting the explicit expression for ,

 |πσs∗|α =x∫w(ω)[1−cos(−2πs∗ω)]dω =x∫w(ω)[πs∗ω]αdω=xmα[w][πs∗]α,

we get a general dependence of the scale parameter—the width of the straggling distribution—on number of recoils:

 σ(x)=[xmα[w]]1/α, (7)

where is the raw generalized -moment of the recoil spectrum:

 mα[w]≡∫ωαw(ω)dω.

Thus, the width of the spectrum increases with the average number of the recoils as . The stability parameter , in turn, increases with from unity to two.

It should be noted, that the scale parameter is equal to half-width of the distribution at of the maximum. At when the distribution approaches the normal (Gaussian) distribution, the scale parameter approaches the square root of Gaussian variance divided by two, .

## Ii Verification of the Method

The two known functional limits of the considered process, with being the Landau distribution and the Fokker-Plank (diffusive limit), may be considered the benchmarks of the method, see, e.g., Ashrafi and Golmankhameh (2017).

### ii.1 The diffusion limit

As stated in the Central Limit Theorem, the sum of independent identically distributed variables with the finite variance should approach the normal (Gaussian) distribution, which is a limiting case of the stable distributions with . Directly following from (2), at , the real part of the exponent approaches Gaussian:

 R[^w−1]≈2π2s2¯¯¯¯¯¯ω2.

The same result, , directly stems from (6) since when .

### ii.2 Landau’s function

A particular case of the stable distributions, the Landau function Landau (1944) (, ), is of special importance since it has undergone extensive study and experimental validation, see Vavilov (1957); Bichsel (1988). The process of ionization losses described by the Landau distribution, agrees with the assumption of the small recoils, whose spectrum is independent of the energy of the particles. The problem is the employment of the idealized unbound recoil spectrum of dependence on energy. This spectrum—the Rutherford cross section—can not be normalized (it has infinite moments).

To avoid the divergence, we consider a truncated recoil spectrum, , then take the limits , , and keeping the total energy losses finite. A “physical” Rutherford cross section is sketched in Fig. 3.

The normalized truncated spectrum reads:

 SL(ω)=sign(ω−a)−sign(ω−b)2ω2ab(b−a). (8)

Fourier transform for this cross-section is:

 F(Rcs) =1(b−a)√2π× (9) {bcos(as)−acos(bs)+sab[Si(as)−Si(bs)]− i[sab(Ci(as)−Ci(bs))−bsin(as)+asin(bs)]}

where is the integral sinus, and is the integral cosine.

Explicitly, we have for the model:

 αL(s;a,b)=1+cos(2πbs)b−cos(2πas)a+2πsab[Si(2bπs)−Si(2aπs)]b−a (10)

The stability parameter (10) has two limits:

 lims→0αL(s;a,b) =2,0

For the physically grounded cases of ionization losses, both the Landau and the Vavilov formulas are valid. The Landau distribution evolves into the Gaussian well beyond the physical region, see Fig. 4.

## Iii Radiation in Periodic Structure

As an example of application of the method to a practical case, we consider evolution of the straggling function due to emission of the undulator radiation, see, e.g., Howells and B.M.Kincaid (1992). Evolution of the straggling function for a long undulator and approximates the diffusion process Agapov and Geloni (2014). On the other hand, for the dipole radiation , and a short undulator (or the entrance section of a long undulator), this function is asymmetric and non-Gaussian Bulyak and Shul’ga (2015, 2017).

Figure 5 represents the straggling function profiles computed in accordance to (2) for a small average number of recoils. It shows that the straggling function resembles the spectrum of recoil at , then it gradually spreads out and smoothes, approximating to some degree the Landau distribution.

The stability parameter against the number of recoils, computed based on (6) for different undulator parameters is presented in Fig. 6. As it can be seen from the figure, the wider the recoil spectrum the later the stability parameter approaches the diffusion limit of .

When the stability parameter approaches the ‘diffusion’ value of (still remaining below it), the skewness , the third centered moment (3c) stays positive and increases with . We can derive practical information about the mode of the distribution. Making use of Pearson’s skewness for a distribution close to the normal, see, e.g. Kenney and Keeping (1962),

 ¯ϵ−ϵmodeσ=μ32σ3,

where is the maximum of the distribution density, and is the centered third moment, we get:

 ϵmode =¯ϵ−μ32σ2 =(1+x)¯¯¯ω−(1+x)¯¯¯¯¯¯ω3−3¯¯¯¯¯¯ω2¯¯¯ω+2¯¯¯ω32[(1+x)¯¯¯¯¯¯ω2−¯¯¯ω2]. (12)

For a big number of the recoils, , the shift of the mode from the mean is almost independent of the number of recoils:

 ϵmode−¯ϵ≈−¯¯¯¯¯¯ω32¯¯¯¯¯¯ω2.

The mode—position of the maximum—is shifted from the mean to smaller energy losses by the constant value, which is determined by the raw moments of the recoil spectrum.

## Iv Summary

A general dependence of the distribution of energy losses by the relativistic electrons due to radiation in periodic structures or ionization losses in matter was analyzed. The straggling function—distribution density of fluctuations—is exclusively determined by the ensemble-average number of recoils having undergone by the particle since entering the field (or medium in the case of ionization losses), and the spectrum of the recoil.

The straggling function was compared to the Lévy stable process as the only attractor of such processes according to the Generalized Central Limit Theorem. The results of this consideration reveal that the width of the electron spectrum is increasing with the number of recoils according to the power law, with the power index being inverse to the stability parameter, i.e., linearly with the number of recoils at the beginning of the process, and in proportion to the square root from the number of recoils at the diffusion limit.

Increase of the spectrum width limits the ability of the beam to generate coherent radiation in hard x-ray and gamma-ray region.

Despite of the assumed independence of the recoil spectrum on the electron’s energy, a ‘negligible’ (from the electron’s point of view) change in this spectrum may play an important role in the reduction of the brightness of sources of hard x-rays and gamma-rays, which employed relativistic electrons. It occurs due to the fact that only a small fraction of the spectrum is used: the pin-hole fraction of the spectrum has strong dependency upon the energy spread of the electrons. The attainable width of the pin-hole collimated radiation—upper limit of it—will exceed doubled electron’s energy spread, Bulyak and Urakawa (2014).

## References

• Bulyak and Shul’ga (2017) E. Bulyak and N. Shul’ga, Nucl Instrum Methods in Physics Research, Sect B 402, 121 (2017).
• Bulyak and Shul’ga (2018) E. Bulyak and N. Shul’ga, Journal of Instrumentation 13, C02051 (2018).
• Applebaum (2004) D. Applebaum, Notices of the American Mathematical Society. Providence, RI: American Mathematical Society 51, 1336 (2004).
• Landau (1944) L. Landau, J. Phys. USSR 8, 201 (1944).
• Nolan (2017) J. P. Nolan, Stable Distributions - Models for Heavy Tailed Data (Birkhauser, Boston, 2017) in progress, Chapter 1 online at fs2.american.edu/jpnolan/www/stable/stable.html.
• Ashrafi and Golmankhameh (2017) S. Ashrafi and A. K. Golmankhameh, Turkish Journal of Physics 1705, 1 (2017).
• Vavilov (1957) P. Vavilov, Soviet Physics JETP 5, 749 (1957).
• Bichsel (1988) H. Bichsel, Rev. Mod. Phys. 60, 663 (1988).
• Howells and B.M.Kincaid (1992) M. Howells and B.M.Kincaid, The properties of undulator radiation, Tech. Rep. (LBL-34751, UC-406, 1992).
• Agapov and Geloni (2014) I. Agapov and G. Geloni, Phys. Rev. ST Accel. Beams 17, 110704 (2014).
• Bulyak and Shul’ga (2015) E. Bulyak and N. Shul’ga, “Electron spectra and coherence of radiation in undulators,” http://arXiv.org/1506.03255v2 (2015).
• Kenney and Keeping (1962) J. F. Kenney and E. S. Keeping, Mathematics of Statistics (Pt. 1, 3rd ed. Van Nostrand, Princeton, NJ, 1962).
• Bulyak and Urakawa (2014) E. Bulyak and J. Urakawa, Journal of Physics: Conference Series 517, 012001 (2014).
Comments 0
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters

Loading ...
326823

You are asking your first question!
How to quickly get a good answer:
• Keep your question short and to the point
• Check for grammar or spelling errors.
• Phrase it like a question
Test
Test description