CC BY
9
DOI: 10.5862/JPM.225.1 9 UDC 538.911; 539.1.03
G.B. Sushko', A.V. Korol'-2, A.V. Solov'yov'
' MBN Research Center at Frankfurt Innovation Center of Biotechnology 2 State Marine Technical University of St. Petersburg
A SMALL-AMPLITUDE CRYSTALLINE UNDULATOR BASED ON 20 GEV ELECTRONS AND POSITRONS: SIMuLATIONS
This paper presents the results of numerical simulations of a crystalline undulator based on the channeling of 20 GeV electrons and positrons. The device considered is characterized by a small amplitude and a short period of periodic bending. Calculations have been performed accounting for all-atom interactions using the MBN Explorer software package. The effect of low crystal thickness (less than a channeling oscillations period) on radiation spectrum was studied. A new scheme to product a high-energy radiation was proposed. It is based on the short-period small-amplitude crystalline undulator and allows decreasing the intensity of the non-undulator part of the spectrum.
small-amplitude crystalline undulator, electron, positron, channeling, simulation, crystalline target.
1. Introduction
A crystalline undulator (CU) [1] is a device producing high-energy radiation by propagating charged relativistic particles through a periodically bent oriented crystalline medium. The feasibility of creating such a device has been predicted theoretically [2 — 4] and now is being studied experimentally [5]. This concept is based on the channeling effect that involves the propagation of projectiles in oriented crystals along crystalline planes or axes. Projectiles in channels interact with crystalline medium and oscillate around the center of the channel. The shape of the channel is responsible for the amplitude and the period of channeling oscillations and the corresponding radiation. Channeling in bent crystals leads to a special type of radiation which depends on the parameters of crystal bending. In recent years two main types of crystalline undulators have been discussed: small-amplitude short-period (SASP) undulators [5 — 7], and large-amplitude large-period (LALP) undulators [5, 8 — 10]. The bending amplitude is considered high if its value is more than the interplanar distance in the crystal. The bending period is considered large if its value is larger than that of the period of channeling oscillations in the channel. Short-period undulators can be used
to generate the radiation with photon energies higher than the characteristic energy of the channeling radiation.
In recent series of experiments at Mainzer Microtron [5] with 600 and 855 MeV electrons the effect of a small-amplitude short-period undulator was observed. Another set of experiments with diamond crystalline undulators is planned within the E-212 collaboration at the SLAC facility (Stanford Linear Accelerator Center, USA) with 10 - 20 GeV electron beam. The current experiments with small-amplitude short-period undulators are based on thin silicon (or diamond) crystals doped with small amount of germanium atoms. These crystals are produced using the Molecular Beam Epitaxy (MBE) technology.
The cases of LALP and SASP undulators can be distinguished by the value of the parameter C = z/(RU' ). In the case of an LALP undulator,
^ max' '
the value of this parameter is considered small, C << 1, because the projectile has to follow bends of crystalline planes or axes to produce appropriate radiation. In this case an increase in the bending amplitude leads to an increase in the radiation intensity, but also causes an increase in C being a limiting factor. In SASP undulators the C parameter can also be formally calculated and its value is more than 1. In this
case projectiles are unable to follow bends of the crystalline medium. What produces the undulator radiation here is the periodic force acting on the projectile due to the rapid change of the interplanar potential. Then the growth of the bending amplitude also boosts the intensity of radiation but reduces the number of channeling projectiles. The amplitude in this instance is limited to a half of the interplanar distance.
In previous papers [11 — 13] the methodology of simulation of channeling with MBN Explorer was presented and applied for straight and bent crystals for sub-GeV energies. The simulations for multi-GeV energies were also described for bent [14] and periodically bent [6] crystals.
MBN Explorer implements a full-atom model of the three-dimensional motion of projectiles in the crystalline medium. In this work, the crystal parameters were taken from experiment [5] as a starting point for the simulations. With these parameters the simulation of propagation of projectiles through the crystals was performed. Using a quasi-classical approach to the calculation of a radiation spectrum the radiation of projectiles was calculated and compared for different crystals. It is shown that for crystals with a thickness less than channeling oscillations period the channeling radiation effect can be suppressed, while undulator radiation can still be generated. The parameters of this radiation are studied numerically.
2. Physical model
The simulation of propagation of a relativistic projectile in a crystalline medium was performed by solving the classical relativistic equations of motion. The force acting on a projectile was calculated as a sum of its interactions with neighboring atoms of the medium. Each atom of the medium was considered as fixed in space as its velocity is much lower than the speed of projectile. Interaction of the projectile with atoms was modeled using the classical Molière [15] interaction potential for screened charges. The atom positions are conditioned by the crystal grid and random displacement due to thermal vibrations. We took the crystalline grid parameters and thermal vibrations from literature [16].
In order to simulate the propagation of a whole beam, the parameters of individual projectiles are randomly sampled. These parameters include an entry point of the projectile into the crystal and an angle between a particle initial velocity and a beam direction. The width of the angular distribution is determined by beam emittance. A random shift of the atoms of the medium due to thermal vibration also causes the sampling of the crystal parameters.
An accelerated motion of charged projectiles in crystal channels produces radiation which can be characterized by a radiation spectrum differential with respect to the photon energy ho> and integrated over the given angular aperture 04.
The simulated trajectories were used to compute the spectral distribution of the emitted radiation. For each set of simulated trajectories of the total number N0 the spectral distribution emitted within the cone 0 < 04 with respect to the incident beam was calculated as follows:
dE (e<e- > = -L £ f2n 4 4,f we-^ ,(1> h4a N0 j~iio Jo k4a4n
here 43Ej / k4a4n is the spectral-angular distribution emitted by a particle moving along the j-th trajectory.
The sum is carried out over all simulated trajectories, i.e. it takes into account the contribution of the channeling segments of the trajectories as well as of those corresponding to the non-channeling regime.
In order to calculate 43EJ / k4a4n a general quasi-classical method developed by Baier and Katkov [17] was used. The quasi-classical approach explicitly takes into account the quantum corrections due to the radiative recoil. The method is applicable in the whole range of the emitted photon energies, except for the extreme high-energy tail of the spectrum where (1 - ka / e) << 1.
Within the framework of the quasi-classical approach, the spectral distribution in energy radiated in the given direction n by an ultra-relativistic particle is given by the following expression (see Ref. [18] for the details):
Jl^ = = r 4^ r 4t eawo-vft» x k4 a 4 n 8n2 J-» 1J 2
xi (1 + (1 + M)2) f ^G^nM - l1 + 41, (2)
I I C J Y J
where
h® < n' r(t) . _ n j. \
u =-, y(t) = t--—, a — a(1 + u).
s - hra c
In the classical limit u ^ 0, a' = a. These equations allow us to compute the emission spectrum for each simulated trajectory. Averaging the spectrum over all trajectories allows calculating the spectrum of the beam, and determines the statistical error.
3. Numerical results
Let us consider planar channeling in straight and periodically bent crystals of silicon and diamond. The thickness of the crystal is 4 ^m, the bending period rku is 0.4 ^m and the bending amplitude a is 0.4 A which is lower than a half of the interplanar distance in both cases. The motion of charged projectiles in the straight crystal corresponds to a classical channeling regime. The second case corresponds to a
a)
SASP crystalline undulator regime. The beam is oriented in a (110) plane, avoiding axial channeling directions.
The results of the simulation of radiation of 20 GeV electrons and positrons are compared for the straight and periodically bent diamond and silicon crystals (Fig. 1). At this energy, the natural emission angle is 0c = 1/y = 25.6 ^rad. The taken value of beam emittance was y = 5 ^rad. The taken angular aperture of the detector 0max = 150 ^rad was 5.8 times higher than the natural emission angle, and covering most of the radiation of projectiles.
In both pictures (see Fig. 1) the rise of each curve (below 1 GeV) corresponds to the channeling radiation, while the peak at ~6 GeV as well as other ones corresponds to the undulator radiation. Bending of a crystal leads to a significant suppression of the channeling peak. This effect can be described in terms of an average potential model of channeling. With increase in the bending amplitude the potential well depth decreases but its width grows. This leads to a decrease in the number of channeling projectiles
0.4
0
0
0 0.3
X
"з 0.2
■«
Ti 0.1
Ш
T3
0.0
. up to 2.0 \up to 3,1 C(110)
A ^ ( \ \ M \ \ A i ff 1 \ x \ У —■- - electrons, straight + electrons, bent -A positrons, straight " , -y- positrons, bent
\ ч ^ ' ■ j
Fig. 1. Radiation spectra for 20 GeV electrons and positrons in the straight and the periodically bent diamond (a) and silicon (b) crystals; a = 0.4 À, I = 0.4 ^m
2.0-
о о о 1.5-
'з ïS 1.0-
LU ТЗ 0.5-
0.0
\ 1 1 1 ' 1 1 1 -
2 _
// -
0.1
10
Photon energy, GeV
Fig. 2. A comparison between radiation spectra (20 GeV electrons) for two bent diamond crystals with
different thickness values: 24 ^m (1) and 4 ^m (2).
The characteristic channeling oscillations period is equal to 9 ^m and falls between the thickness values. For the sake
of comparison, the curve 2 is multiplied by the factor of 6.
and in frequencies of channeling oscillations.
Another factor of suppression of the channeling radiation arises from the high parameter values being unfit for thin crystals; for example, the characteristic periods of channeling oscillations in simulated trajectories for positrons are equal to 9 ^m in C(110) channels and 11.5 ^m in Si(110) channels at the energy of projectiles of 20 GeV; these periods are longer than some crystal thickness values. This difference causes the suppression of the channeling effect in thin crystals. The comparison between the radiation spectra for a thick (24 ^m) crystal and a thin (4 ^m) one is shown in Fig. 2. The values for the thinner crystal are multiplied by a factor of 6 for comparison. As can be seen from the plots, the undulator peaks of both curves nearly coincide, while the curve shapes in the region of lower
photon energies are significantly different. The channeling produces a sharp peak of radiation around the energy of 150 MeV for a thick crystal, while the corresponding peak is absent for a less thickness, and the radiation is produced by a synchrotron effect.
The radiation emission by projectiles in this simulation is directed strongly along the beam. For the small-aperture undulator the radiation is emitted within a narrow energy range. The dependence of the emitted radiation spectrum on the detector aperture 0 is shown in
* max
Fig. 3. In the case of the narrow aperture 0max = 10 ^rad the peak in the spectrum is rather sharp. With an increase in the aperture the absolute value of the number of high-energy photons grows, but the peak width grows as well. At 0max = 150 ^rad the radiation spectrum saturates and does not change at higher aperture values.
Fig. 3. A comparison between the radiation spectra (20 GeV electrons) in CU for different apertures 0m
1 ^rad (1), 30 ^rad (2), 150 ^rad (3) "
3.1. Crystalline undulator stack
The effect of suppression of channeling radiation in short-periodically bent crystals can be used in order to produce undulator radiation with a higher efficiency. To increase the energy of undulator radiation without increasing channeling radiation, it is possible to use a stack of short crystalline undulators instead of one long undulator. The scheme of such a crystalline system is shown in Fig. 4. In this system, the projectile passes several layers of periodically bent crystalline medium, the radiation produced at all layers adds to the total radiation produced by projectile. For an SASP undulator the thickness of the layers could be taken in the interval between the bending period 1u and the characteristic channeling period of projectile. Such choice of the parameters leads to the absence of full channeling oscillation periods in each channeling segment of trajectory of projectile which results in the suppression of channeling radiation. The effect of undulator radiation in the system remains and grows with an increase in the number of layers.
To simulate the effect of channeling in CU stack the following system was modeled. A set of layers (l — 4 ^m) of periodically bent crystals was put into the simulation box with gaps (l' — 4 ^m) between layers, the period of crystal bending was set to 1 — 0.4 ^m. Such a
Photon energy, GeV
Fig. 5. Radiation spectra for the small (a) and of layers (l = 4 ^m) in a stack: 1 (curve 1), 4 (2),
a = 0.4 Â.
1 2 ... n
l Г Л.
Fig. 4. The scheme of the CU stack made of 1, 2, ..., n periodically bent crystal layers, each of the thickness l. The quantity l stands for the gap between the layers, Xu is the bending period.
Arrowed lines illustrate the trajectories of projectiles
system corresponds to a case of the crystalline undulator stack. In this instance the projectiles are captured in the first crystal, escape from the crystal at some point, and have to be captured again in the next crystal layer. The process of the recapturing of the projectiles in the channeling mode leads to the increase in angular dispersion in velocity of projectiles and the decrease in the number of channeling particles as the number of layers increases.
The thickness of layers in this simulation is set below the channeling oscillations period (11.5 ^m) but above the period of crystal
Photon energy, GeV
large (b) apertures calculated for different number 8 (4), 12 (5); 9max — 250 ^rad (a) and 15.6 ^rad (b);
, 1u — 0.4 ^m max
Photon energy, GeV
Fig. 6. A comparison between radiation spectra formed by positrons in a 24-^m-thick crystal (1, 3) and in a stack of 6 crystals of l = 4 ^m (2, 4). The spectra are calculated for two different apertures: 15.6 ^rad (3, 4) and 250 ^rad (1, 2).
The beam energy is 20 GeV; a = 0.4 Â, Xu = 0.4 ^m
bending (0.4 ^m). Such a choice of the crystal thickness leads to the suppression of channeling radiation due to the lack of full oscillation periods, but the undulator radiation remains.
The comparison between radiation spectra for the different number of layers in a stack is given in Fig. 5. The energy of both synchrotron and undulator radiation grows linearly with an increase in the number of stack layers up to 4 (the case for a small aperture). For the higher number of stack layers the angular distribution of projectile velocities becomes wider and the radiation in a narrow cone saturates. For a higher aperture the growth of the peak intensity continues.
The radiation from a single crystal undulator (the crystal is 24 ^m thick) is compared to a stack of them (6 * 4 ^m) in Fig. 6. It can be seen, that in the stack case the channeling peak in the radiation spectrum is suppressed, while the undulator peak is nearly the same as the one for the case of a single crystal.
4. Summary
This paper deals with the simulation of
photon emission by 20 GeV electrons and positrons in a small-amplitude crystalline undulator. The obtained results have been presented and thoroughly discussed. The simulations were performed using the MBN Explorer software package with accounting for all-atom interactions of projectiles. It is shown that for high-energy projectiles in thin periodically bent silicon and diamond crystals the effect of channeling is highly suppressed, but the motion in a periodically bent medium generates intensive radiation. This regime is favorable for the construction of a crystalline undulator based on light sources as the energy of a projectile is transferred mostly to the undulator radiation.
Acknowledgments
The work was supported by the European Commission (the CUTE-IRSES project, grant GA-2010-269131). The possibility to perform computer simulations at the Frankfurt Center for Scientific Computing is gratefully acknowledged. AVK acknowledges the support from the Alexander von Humboldt Foundation.
REFERENCES
[1] A.V. Korol, A.V. Solov'yov, W. Greiner,
Channeling and Radiation in Periodically Bent Crystals, Springer-Verlag Berlin Heidelberg, second edition, 2014.
[2] A.V. Korol, A.V. Solov'yov, W. Greiner, Coherent radiation of an ultrarelativistic charged particle channeled in a periodically bent crystal, J. Phys. B. 24 (1998) L45-L53.
[3] A.V. Korol, A.V. Solov'yov, W. Greiner, Photon emission by an ultra-relativistic particle channeling in a periodically bent crystal, Int. J.
Mod. Phys. 8 (1999) 49-100.
[4] A.V. Korol, A.V. Solov'yov, W. Greiner, Channeling and Radiation in Periodically Bent Crystals, Springer-Verlag Berlin Heidelberg, 2013.
[5] T.N. Wistisen, K.K. Andersen, S. Yilmaz, et al., Experimental realization of a new type of crystalline undulator, Physical review letters. 112 (2014) 254801.
[6] V.G. Bezchastnov, A.V. Korol, A.V. Solovyov, Radiation from multi-GeV electrons and positrons in periodically bent silicon crystal, Journal of
Physics B: Atomic, Molecular and Optical Physics. 47 (2014) 195401.
[7] V.V. Tikhomirov, A benchmark construction of positron crystal undulator, arXiv preprint arXiv:1502.06588, 2015.
[8] V.G. Bezchastnov, A.V. Korol, A.V. Solovyov, Radiation from multi-GeV electrons and positrons in periodically bent silicon crystal, Journal of Physics B: Atomic, Molecular and Optical Physics. 47 (2014) 195401.
[9] A. Mazzolari, E. Bagli, L. Bandiera, et al., Steering of a sub-Gev electron beam through planar channeling enhanced by rechanneling, Physical Review Letters. 112 (2014) 135503.
[10] Kostyuk, Crystalline undulator with a small amplitude and a short period, Physical Review Letters. 110 (2013) 115503.
[11] G.B. Sushko, V.G. Bezchastnov, I.A. Solovyov, et al., Simulation of ultra-relativistic electrons and positrons channeling in crystals with MBN Explorer, Journal of Computational Physics. 252 (2013) 404-418.
[12] G. Sushko, A. Korol, W. Greiner, A. Solov'yov, Sub-GeV electron and positron channeling in straight, bent and periodically bent
silicon crystals, Journal of Physics: Conference Series. 438 (2013) 012018.
[13] G. Sushko, V. Bezchastnov, A. Korol, et al., Simulations of electron channeling in bent silicon crystal, Journal of Physics: Conference Series. 438 (2013) 012019.
[14] G.B. Sushko, A.V. Korol, A.V. Solovyov, Multi-GeV electron and positron channeling in bent silicon crystals, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 2015.
[15] G. Molière, Theorie der Streuung schneller geladener Teilchen I: Einzelstreuung am abgeschirmten Coulomb-Feld, Z. f. Naturforsch. A 2 (1947) 133-145.
[16] D.S. Gemmel, Channeling and related effects in the motion of charged particles through crystals, Rev. Mod. Phys. 46 (1947) 129-227.
[17] V.N. Baier, V.M. Katkov, Processes involved in the motion of high energy particles in magnetic field, Zh. Eksp. Teor. Fiz. 53 (1967) 1478-1491.
[18] V.N. Baier, V.M. Katkov, V.M. Strakhovenko, Electromagnetic Processes at High Energies in Oriented Single Crystals, World Scientific, Singapore, 1998.
THE AUTHORS
SuSHKO Gennady B.
MBN Research Center at Frankfurt Innovation Center of Biotechnology Altenhoferallee 3, 60438 Frankfurt am Main, Germany [email protected]
KOROL Andrei V.
State Marine Technical University of St. Petersburg
3, Lotsmanskaya St., St. Petersburg 190008, Russian Federation,
and MBN Research Center at Frankfurt Innovation Center of Biotechnology
Altenhoferallee 3, 60438 Frankfurt am Main, Germany
SOLOV'YEV Andrei V.
MBN Research Center at Frankfurt Innovation Center of Biotechnology Altenhoferallee 3, 60438 Frankfurt am Main, Germany [email protected]
Сушко Г.Б., Король A.B., Соловьев A.B. мОДЕЛИРОВАНИЕ мАЛОАмПЛИТУДНОГО КРИСТАЛЛИЧЕСКОГО ОНДУЛЯТОРА НА ОСНОВЕ ЭЛЕКТРОНОВ И ПОЗИТРОНОВ С ЭНЕРГИЕЙ 20 ГэВ.
В статье приводятся результаты численного моделирования малоамплитудного кристаллического ондулятора на основе электронов и позитронов с энергией 20 ГэВ. Расчеты были проведены с использованием пакета программ MBN Explorer, учитывающего взаимодействие налетаюших частиц с кристаллом на атомарном уровне. Рассмотрен случай тонких кристаллических мишеней, толщина которых не превышает характерного периода осцилляции частиц при каналировании. Предложена схема получения интенсивного излучения высоких энергий, на основе малоамплитудного ондулятора, которая позволяет снизить уровень фонового излучения.
малоамплитудный кристаллический ондулятор, электрон, позитрон, каналирование, численное моделирование, кристаллическая мишень.
СПИСОК ЛИТЕРАТУРЫ
[1] Korol A.V., Solov'yov A.V., Greiner W.
Channeling and Radiation in Periodically Bent Crystals, Springer-Verlag, Berlin, Heidelberg, second edition, 2014.
[2] Korol A.V., Solov'yov A.V., Greiner W. Coherent radiation of an ultrarelativistic charged particle channeled in a periodically bent crystal // J. Phys. B. 1998. Vol. 24. Pp. L45-L53.
[3] Korol A.V., Solov'yov A.V., Greiner W. Photon emission by an ultra-relativistic particle channeling in a periodically bent crystal // Int. J. Mod. Phys. 1998. Vol. 8. Pp. 49-100.
[4] Korol A.V., Solov'yov A.V., Greiner W. Channeling and Radiation in Periodically Bent Crystals, Springer-Verlag, Berlin, Heidelberg, 2013.
[5] Wistisen T.N., Andersen K.K., Yilmaz S., et al. Experimental realization of a new type of crystalline undulator // Physical Review Letters. 2014. Vol. 112. P. 254801.
[6] Bezchastnov V.G., Korol A.V., Solovyov A.V. Radiation from multi-GeV electrons and positrons in periodically bent silicon crystal // Journal of Physics B: Atomic, Molecular and Optical Physics. 2014. Vol. 47. P. 195401.
[7] Tikhomirov V.V. A benchmark construction of positron crystal undulator // arXiv preprint arXiv:1502.06588. 2015.
[8] Bezchastnov V.G., Korol A.V., Solovyov A.V. Radiation from multi-GeV electrons and positrons in periodically bent silicon crystal // Journal of Physics B: Atomic, Molecular and Optical Physics. 2014. Vol. 47. P. 195401.
[9] Mazzolari A., Bagli E., Bandiera L., et al. Steering of a sub-GeV electron beam through planar channeling enhanced by rechanneling // Physical Review Letters. 2014. Vol. 112. P. 135503.
[10] Kostyuk A., Crystalline undula-
tor with a small amplitude and a short period // Physical Review Letters. 2013. Vol. 110. P. 115503.
[11] Sushko G.B., Bezchastnov V.G., Solovyov
I.A., et al. Simulation of ultra-relativistic electrons and positrons channeling in crystals with MBN Explorer // Journal of Computational Physics. 2013. Vol. 252. Pp. 404-418.
[12] Sushko G., Korol A., Greiner W., Solov'yov A., Sub-GeV electron and positron channeling in straight, bent and periodically bent silicon crystals // Journal of Physics: Conference Series. 2013. Vol. 438. P. 012018.
[13] Sushko G., Bezchastnov V., Korol A., et al. Simulations of electron channeling in bent silicon crystal // Journal of Physics: Conference Series. 2013. Vol. 438. Pp. 012019.
[14] Sushko G.B., Korol A.V., Solovyov A.V. Multi-GeV electron and positron channeling in bent silicon crystals // Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms. 2015 (in print).
[15] Molière G. Theorie der Streuung schneller geladener Teilchen I: Einzelstreuung am abgeschirmten Coulomb-Feld // Z. f. Naturforsch. A. 1947. Vol. 2. Pp. 133-145.
[16] Gemmel D.S. Channeling and related effects in the motion of charged particles through crystals // Rev. Mod. Phys. 1947. Vol. 46. Pp. 129-227.
[17] Baier V.N., Katkov V.M. Processes involved in the motion of high energy particles in magnetic field // Zh. Eksp. Teor. Fiz. 1967. Vol. 53. Pp. 1478-1491.
[18] Baier V.N., Katkov V.M., Strakhovenko V.M. Electromagnetic Processes at High Energies in Oriented Single Crystals, World Scientific, Singapore, 1998.
СВЕДЕНИЯ ОБ АВТОРАХ
СУШКО Геннадий Борисович — научный сотрудник Научно-исследовательского центра мезобиона-носистем (MBN); г. Франкфурт-на-Майне, Германия. Altenhoferallee 3, 60438 Frankfurt am Main, Germany [email protected]
КОРОЛь Андрей Владимирович — кандидат физико-математических наук, доцент кафедры физики Санкт-Петербургского государственного морского технического университета, научный сотрудник Научно-исследовательского центра мезобионаносистем (MBN); г. Франкфурт-на-Майне, Германия. 198262, Российская Федерация, г. Санкт-Петербург, Ленинский пр., 101. [email protected]
СОЛОВьЕВ Андрей Владимирович — доктор физико-математических наук, профессор, директор Научно-исследовательского центра мезобионаносистем (MBN); г. Франкфурт-на-Майне, Германия. Altenhoferallee 3, 60438 Frankfurt am Main, Germany [email protected]
© Санкт-Петербургский политехнический университет Петра Великого, 2015
