CC BY
11
ТЕОРЕТИЧЕСКАЯ ФИЗИКА
DOI: 10.5862/JPM.225.1 8 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
ELECTRON AND POSITRON PROPAGATION IN STRAIGHT AND PERIODICALLY BENT AXIAL AND PLANAR SILICON CHANNELS
In this paper the results of simulations of axial and planar channeling of electrons and positrons in straight and periodically bent Si crystals are presented. Simulations with direct calculation of trajectories of projectiles accounting for all-atom interactions were carried out using the MBN Explorer software package. The full atomistic approach for particle trajectories simulation allows to quantitatively compare axial and planar channeling processes. The results of the simulations show significantly lower dechanneling length and number of channeling projectiles in the axial channeling case. For this case the dependence of characteristics of the channeling process on the choice of an axis direction and on a direction of the crystal bending have been investigated.
electron, positron, channeling, silicon crystal, periodically bent channel, trajectory.
1. Introduction
Channeling is an effect of propagation of charged relativistic projectiles in a crystalline medium along crystal planes and axes [1]. Motion of projectiles in the channels can be characterized by the average length of staying in a channel, the fraction of particles that are captured in a channel and a spectrum of radiation emitted by the particles.
The concept of a crystalline undulator (CU) as a source of undulator-like electromagnetic radiation in the high energy range up to the MeV region was formulated in Ref. [2] and further studied in Ref. [3, 4]. In crystalline undulators the projectile particles follow periodically bent channels and emit undulator radiation in addition to channeling radiation characteristic for the case of channeling in linearly oriented crystals.
In recent years several experiments were performed [5 — 7] to detect the radiation from electron- or positron-based CU units.
The most recent works in this field are now in progress at the Mainz Microtron (Germany) facility for the 195—855 MeV electrons, and at the SLAC facility (Stanford Linear Accelerator Center, USA) with 10-20 GeV electron beams. In these experiments planar channeling of electrons is studied in bent [8] and periodically bent crystals.
In order to simulate processes of planar and axial channeling a new module was developed for MBN Explorer [9, 10]. MBN Explorer is a versatile software package for simulating molecular systems of various degrees of complexity. MBN Explorer utilizers a broad variety of interatomic potentials to describe different molecular systems, such as atomic clusters, fullerenes, nanotubes, polypeptides, proteins, DNA, composite systems, nanofractals and many more.
Computer simulations of planar channeling using the MBN Explorer software were performed in the previous works [11-15].
The effect of axial channeling in straight crystals was studied by different groups both experimentally [16] and numerically using the averaged potential method [17—19].
This work studies axial channeling and compares it with the planar case. The results of numerical simulation of axial channeling of electrons and positrons in straight and periodically bent silicon (Si) crystals are presented. The beam energy range of 10— 20 GeV corresponds to the beam energies available at the SLAC facility.
The modeling was carried out using the full-atom simulation of three-dimensional motion of projectiles in the crystalline medium. The interaction of a projectile with atoms is modeled using the Molière approximation. Such model itself is not limited to a specific structure of channels and direction. The use of this approach allows simulations of both axial and planar channeling using the same methodology and the direct comparison of the results. The efficient implementation of simulation algorithms allows modeling trajectories in thick crystals up to macroscopic sizes (of ~1 cm).
It is shown that the characteristics of the channeling effect depend strongly on the selection of the beam direction in a crystal. Channeling parameters for three axes (<100>, <110> and <111>) are presented and compared with the parameters for planar channeling in a (110) plane. Simple analysis of the average potential of channeling for three axes is also presented. The optimal directions of bending are chosen for the case of periodically bent crystals. It is shown that the effects of both dechanneling and rechanneling depend strongly on local curvature of the channel for periodically bent channels. This leads to the periodicity in the derivative of the number of channeling projectiles on the penetration distance.
2. Physical model
The motion of an ultra-relativistic projectile of the charge q and the mass m in an external electrostatic field with potential energy U (r) can be described with relativistic equations of motion written in the following form:
dr / dt = v, dp / dt = -qdU / dr.
The momentum p, written in terms of velocity, reads p = myv, where y is the Lorentz factor
Y
= V 1 - V2/ с2 =8 /
mc
with s being the projectile energy.
The differential equations (1) are to be integrated for t > 0 using the initial values of the coordinates (x0, y0, z0) and the velocity components (vx0, vy0, vz0) of the particle. To ensure an accurate numerical integration the fourth-order Runge-Kutta scheme and a time step variation algorithm were implemented.
The description of the electrostatic field is done using the Molière approximation [20] of the electrostatic potential of a neutral atom:
UM(p) = — S(p);
S (p) = X3=1 ae
(2)
-Pip/OTF
(1)
where Z is the atomic number, the coefficients in the screening function £(p) are as follows: a12,3 = (0.35; 0.55; 0.10) (so that £^a, = 1) and P123 = (0.3; 1.2; 6.0). The Thomas-Fermi radius is related to the Bohr radius an via
TF 0
aTF = 0.8853Z-1/3a0.
In order to simulate the motion of a particle in the medium we used a dynamic simulation box technique (see Ref. [11]). With this approach the crystalline medium is generated in the close vicinity of the particle taking the predefined unit cell of the crystal and the set of transformations: rotation, displacement and bending.
In order to take into account thermal vibration of atoms in the grid structure, atoms are randomly shifted from their nodal positions. Each component of the displacement vector is normally distributed with the root-mean-square amplitude uT = 0.075 A which corresponds to the room temperature [21].
Through such an approach, the motion in any crystalline system can be described by setting the appropriate crystal unit cell, rotation transformation of the system and a root-mean-square amplitude of thermal vibrations. Characteristics of the channeling process strongly depend on the properties of the medium, the amplitude of the vibrations
Fig. 1. Axial channels structure for three directions: <100> (a), <110> (b), <111> (c). Closed circles correspond to the lines of Si atoms which are the centers of the channels for electrons. Dots correspond to the centers of channels for positrons. Channels themselves for positrons are highlighted with light grey,
channels for electrons — with dark grey
and the direction of the beam in the medium.
In this article, the following directions were chosen for comparison: axes <100>, <110>, <111>, and a plane (110). The structure of channels for axial directions is illustrated in Fig. 1. This structure is a result of the Voronoi decomposition of the plane being perpendicular to the beam direction. Coordinates of the projections of atoms are used as seeds of this decomposition for electrons. As for positrons, those with the minimum of the potential energy surface are taken as seeds of the decomposition. This description of channels can be applied to any axial direction and any type of lattice.
The axial channels' structure depends on the projectile charge. For electrons the channel center coincides with the line of Si atoms. For positrons the channel center is situated between the lines of atoms. The projectile was considered to be in a channel if the distance to this channel center is the shortest among the ones to all channels for the given type of projectiles. If the projectile stays in one channel and changes the velocity sign in the X or Y directions four times, then it is considered as captured in a channeling mode.
The concept of the crystalline undulator implies a motion of a projectile in a periodically bent crystal. This bending leads to a change of the radiation spectrum and produces undulator radiation in addition to channeling radiation. In simulations, the periodic crystal bending is implemented using a simple transformation of
the coordinates of the particle:
r ^ r + d sin(2np • r + s),
where d is a bending amplitude, p is an inverse period vector modulus and s is a phase shift.
Such a transformation can be applied to either axial or planar channels in order to simulate the crystalline undulator.
3. Numerical results
In order to analyze the effect of axial channeling, numerical calculations oftrajectories were performed using random sampling with different initial positions of projectiles in the channel and with random positioning of Si atoms in the medium. The average number of trajectories for each case was at a level of 104, which allowed reducing the statistical error of obtained quantitative characteristics to a level of few percent.
In all calculations, the Z axis was aligned with the main direction of channeling (the beam direction). The Z coordinate of the particle is a penetration distance. The used combinations of the Z—Y axes were as follows:
<100> - <011>; <110> -<111> - <011>.
<110>;
For planar channeling the taken axis directions were as follows:
Z(<10,10, 1>), Y(<110>).
The beam direction along the Z axis was taken in the (110) plane.
Fig. 2. Examples of simulated trajectories of an electron (a) and a positron (b) in the <110> axial channel in the Si single crystal with Face Centered Cubic (FCC) structure. The trajectories illustrate the channeling regime of both particles (a, b) and a motion outside the channel and the electron rechanneling (a). Dots mark the channel centers
3.1. Axial channeling in straight crystals
A trajectory of the projectile is the direct result of the channeling simulation. Examples of obtained trajectories for electrons and positrons are shown in Fig. 2.
Several characteristic regimes of the motion can be seen: the channeling motion in the vicinity of the channel center, the rechanneling events, i.e. a capture into the channeling mode, and, finally, the dechanneling of a projectile from a channel and the subsequent motion above the potential barrier through the crystal. The channeling segment is the part of the trajectory between the rechanneling and the dechanneling events. Primary channeling segment is the one that starts at the crystal entrance.
The following quantities were used in order to characterize the channeling effect: N0 is the total number of simulated projectiles, Nc (z) is
the number of projectiles in a channeling mode at a depth z, Nc1(z) is the number of projectiles in a primary channeling mode, n(z) = Nc (z) / N0 is a fraction of particles in a channeling mode in any channel, n(z) = Nc1(z) / N0 is a fraction of particles in a channeling mode in the primary channel, i. e. the number of channeled particles excluding those that experience rechanneling.
In order to quantitatively describe a channeling we used the following characteristics: A is an acceptance (a fraction of particles Nc (0) / N0 that start channeling right from the crystal entrance); L 1, Lp2 are mean values of the primary and of all ^channeling segments, respectively; Ld is a characteristic dechanneling length in exponential approximation
n(z) = A exp(-z / Ld),
where n (z) is the fraction of projectiles in a primary channel.
Table 1
The obtained channeling characteristic values in different axial directions for 10 GeV electrons and positrons in the straight 320 ^m thick Si single crystal
Projectile, axis A,% Lp p Ит lpV ит LcP ит
e-, <100> 49 ± 2 21 ± 1 23 ± 1 15 ± 1
e-, <110> 36 ± 2 19 ± 1 21 ± 1 13 ± 1
e-, <111> 48 ± 2 21 ± 1 23 ± 1 14 ± 1
e+, <100> 65 ± 2 246 ± 5 111 ± 3 2760+50
e+, <110> 77 ± 2 252 ± 5 121 ± 3 3670+70
e+, <111> 40 ± 2 212 ± 8 63 ± 2 2660+50
Fig. 3. Channeling fractions for 10 GeV electrons (a) and positrons (b) in a straight 320-^m-thick silicon crystal versus the penetration distance for the indicated axial directions. For each direction we calculated
two fractions: with and without rechanneling
Let us examine the axial channeling of 10 GeV electrons and positrons in a Si single crystal in three axial directions <100>, <110>, <111> (see Table 1). Three axial cases exhibit very different acceptance values, and less different dechanneling lengths. These values also depend strongly on the charge sign of the projectile. The <110> axis is the most favorable for positrons but the least favorable for electrons. The dechanneling length for electrons is at least an order of magnitude lower than that for positrons due to a closer passage of projectiles about the axis during channeling. The exponential approximation is quite good for electrons, but it does not describe channeling fraction for positrons (see Fig. 3).
For positrons the values of L and Lp2 are more limited by the total crystal thickness (of 320 ^m) than by the real dechanneling
length. Fig. 3 shows that the channeling fractions exhibit a very small decrease even at the crystal exit. The difference between
L , and L., and low values of L . indicate a
pi p2' p2
great number of rechanneling events and short channeling segments.
In order to describe the difference between channeling in various channels let us discuss the picture of average potential energy of interaction between positrons and atomic strings along different axial directions (see Fig. 4). For electrons the interaction energy values have an opposite sign.
For a positron channel, the most important parameter is the difference between the potential energies calculated at two points: (i) the center of the channel (which is the equilibrium position, i.e. the potential minimum) and (ii) the boundary of two adjacent channels. This
г з 4 s X, А
Fig. 4. 3D picture of the average potential energy (in eV scale with tones) of interaction between positrons and atomic strings along the <100> (a), <110> (b), <111> (c) axes, and the (110) atomic plane (d) as well. Dashed lines mark the channels for positrons. Arrows point the characteristic directions of dechanneling from the channel centers
parameter, U0, is a potential well depth and its growth lowers the probability of the projectile dechanneling (see Table 2).
For the axes <100> and <111> the barriers are the same in any direction, for <110> the barriers are very different for different directions and much bigger than those for other axes. This leads to better channeling in the <110> axis and better channeling in the case of bent crystals for bending in the Y direction.
3.2. Axial channeling in periodically bent crystals
For periodically bent crystals, channeling in each direction was simulated with a different amplitude of bending a, the bending period was 1 = 39.6 ^m. These results were compared to those in straight crystals (see Table 3). For electrons the channeling length was found to
be lower than the bending period, so this case was omitted.
For periodically bent crystals the difference between the channels is much more than that for straight crystals. In all channels the growth of the bending amplitude leads to a decrease in both acceptance and mean dechanneling lengths values. For axis <111> this decrease is the greatest and can be attributed to the least value of the mean potential well depth. For the rest of the axes in question the drop in dechanneling length is less and the length value is much more.
Due to the dependence of the potential barrier height on direction in case of axis <110> the effect of crystal bending also depends on the bending direction. Two variants of crystal bending are given in Table 3. Bending in the Y direction leads to lesser losses in both
Table 2
The values of the potential well depth calculated in different directions for 10 GeV positrons channeling in a straight 320-^m-thick Si crystal
Direction <100> <110> <111> Y
-U , eV 2.89 5.04 0.43 23.28
Table 3
The channeling characteristics in different axial directions for 10 GeV positrons in a periodically bent
320-^m-thick Si single crystal
Axis a, À ( Y direction) A, % L^ цт
0 65 ± 2 246 ± 5 111 ± 3
<100> 1 54 ± 3 168 ± 12 65 ± 4
2 34 ± 3 113 ± 11 43 ± 2
0 77 ± 2 252 ± 5 121 ± 3
1 69 ± 3 177 ± 10 84 ± 4
<110> 2 57 ± 1 171 ± 4 72 ± 1
1 (X direction) 51 ± 4 152 ± 13 64 ± 4
2 (X direction) 33 ± 3 75 ± 8 36 ± 1
0 40 ± 2 212 ± 8 63 ± 2
<111> 1 27 ± 3 35 ± 3 25 ± 1
2 12 ± 2 19 ± 1 19 ± 1
Notes: the zero bending amplitude a corresponds to the straight crystal, corresponding values are the same as those in Table 1 and are presented for comparison; the bending period X = 39.6 ^m.
acceptance and channeling length values than those for the X direction.
3.3. Rechanneling in periodically bent crystals
The finite crystal curvature increases the importance of both dechanneling and rechanneling effects for two reasons. First, the curvature of the crystal plane causes a centrifugal force and results in the dechanneling. Second, the bend of the channel axis leads to the change in the angle between the axis tangent and the velocity of over-barrier projectiles which, in turn, causes the volume capture and the
rechanneling. Both effects depend periodically on the projectile coordinate.
Let us discuss the dependence of channeling fractions on penetration distance with and without rechanneling for the periodically bent crystal (Fig. 5, a). The bending is done along the <110> axis with an amplitude of 2 A. The curve of channeling fraction with rechanneling consists of a normal decreasing part and includes an oscillatory part with the bending period of the crystal.
To analyze this effect quantitatively let us consider Fig. 5, b where the derivative of
u
ГО
Ш
с
с m -С
с;
z, Mm Z, Mm
Fig. 5. The channeling fraction (a) and its derivative (b) versus the penetration distance for 10 GeV positrons in a silicon crystal (320 ^m thick) periodically bent along the <110> axis; a = 2 À, X = 39.6 ^m
the channeling fraction (with rechanneling) is compared with the corresponding approximation that follows the fitting formula
dn / dz = n(z) sin(2nz / X) + const,
where n(z) is the channeling fraction, X is the period of crystal bending.
This fitting indicates both an exponential decay of the number of particles in the channel and a dependence of the bending profile on the value of the first derivative. It indicates that the probability of dechanneling is proportional to the angle between the beam direction and the tangent of the channel centerline.
3.4. A comparison of the axial and planar channeling
In order to compare the axial and the planar channeling, some simulations were performed. In these simulations channeling in straight and periodically bent crystals was compared for (110)
planar and <110> axial channels of a silicon crystal. In the planar case the channels are formed by equidistant planes spaced at 1.92 A intervals, in the axial case the channels are two-dimensional (see Fig. 2). The comparison of the channeling parameters calculated for the two cases is presented in Table 4. Both acceptance and length characteristics are higher in the planar case. Lower relative values of Lp2 in case of axial channeling also indicate a high number of short channeling segments.
The results of the simulations show that the channeling of positrons in plane (110) is characterized by a much longer channeling length and a higher acceptance. The only drawback of the channeling in this case is the angular distribution of the outgoing particles in the direction perpendicular to both the normal of the plane and the direction of the beam. This distribution is limited for axial channeling and is much broader for the planar case.
Table 4
Comparison between the calculated channeling characteristics of two channels for 10 GeV positrons in a periodically bent 320-^m-thick Si crystal
Direction a, Â A, % lpi, Lp2, Hm
<110> 0 ll ± 2 252 ± 5 121 ± З
2 51 ± 1 111 ± 4 12 ± 1
0 91 ± 1 З02 ± 4 212 ± б
(110) 2 90 ± 2 З01 ± 5 25б ± 8
4 82 ± З 281 ± 1 209 ± 9
Table 5
Comparison of acceptance and channeling lengths for axial and planar channeling of electrons with different
energies in a straight 320-^m-thick Si crystal
Direction e, GeV A, % Lpl, Hm Lp2, Hm Ld, nm
0.855 З2 ± 4 5 ± 1 5 ± 1 З ± 1
Axial 1.6o З4 ± З б ± 1 1 ± 1 4 ± 1
<110> б.1 Зб ± З 15 ± 1 11 ± 1 10 ± З
10 Зб ± 2 19 ± 1 21 ± 1 1З ± 1
20 З1 ± 4 З1 ± З З1 ± 2 21 ± З
0.855 59 ± 4 1З ± 1 14 ± 1 9 ± З
Planar 1.6o б1 ± 4 19 ± 1 21 ± 1 14 ± 1
(110) б.1 б4 ± 4 б2 ± 5 5б ± З 49 ± 4
10 б5 ± 4 81 ± 1 1З ± 4 б8 ± 10
20 б8 ± З 129 ± 10 112 ± б 110 ± 10
For electrons the difference between the axial and the planar channelings turned out to be more pronounced (see Table 5). The acceptances in the axial and planar cases differ by a factor of 2-3, whereas the dechanneling lengths differ by a factor of 2-5 depending on the calculation procedure. In this case the dechanneling length is shorter than the bending period of the crystal and therefore the results on the periodically bent crystal are omitted.
The effect of the shorter dechanneling length in the axial case is caused by a higher probability of collision with the atoms of the medium. Trajectories of channeling electrons are located in the close vicinity of atom strings (marked by dark and light tones in Fig. 4). The local density of atoms in this region is several times higher than that in the case of planar channeling, and this difference results in a higher dechanneling rate.
To quantitatively estimate this effect, let us consider a cylinder disposed along a row of atoms along the <110> axis (Fig. 6, a). The local density of atoms in such a cylinder depends on its radius and for high radius values approaches the average density of the crystal. The relationship between the local density of the medium and the cylinder radius for planar and axial channels is shown in Fig. 6, b. This density value is the same for two cases at the large radius
limit (the average crystal density is the same for both cases), but becomes several times greater in axial case for the radius values less than 0.5 A. For small radius values (e.g. of 0.3 A), such a cylinder contains 3.5 times more atoms than the same one disposed in the (110) plane. This dependence results in a higher probability of hard collision of projectiles leading to very low impact parameter and dechanneling.
4. Summary
In this paper, the general simulation approach of the channeling process using MBN Explorer has been described. The propagation of projectiles through the medium was simulated by solving a three-dimensional equation of motion of a relativistic projectile interacting with individual atoms of the medium. This approach allows high-precision simulation of any channeling types in various media including periodically bent crystals.
The results of the numerical simulation of the axial channeling of positrons and electrons in straight and periodically bent crystals are presented. It is shown that the main characteristics of this effect (an acceptance and a channeling length) depend strongly on the axis direction and the sign of the projectile charge. The results of simulation of position channeling in periodically bent crystals are also
Fig. 6. Schematic illustration of the local density calculation for planar (top) and axial (bottom) channels (a), and comparison plots of the local density of atoms versus inverse cylinder radius for planar (1) and axial (2) channels (b)
analyzed and presented. The optimal beam and bending directions obtained on the basis of numerical simulations and an analysis of the averaged potential energy are given as well. Channeling in this mode is characterized by the highest values of the acceptance and the dechanneling length.
The universal approach to simulation of medium allows to directly compare the planar and the axial channeling types. It is shown that planar channeling is characterized by values of acceptance and dechanneling length higher than
those for the axial type due to the difference in the local medium density.
Acknowledgments
The work was supported by the European Commission (the CUTE-IRSES project, grant GA-2010-269131). The opportunity 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] J. Lindhard, K. Dan, Stopping power of electron gas and equipartition rule, Vidensk. Selsk. Mat. Fys. Medd. 34 (1965) 1-64.
[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] M. Tabrizi, A.V. Korol, A.V. Solov'yov, W. Greiner, Feasibility of an electron-based crystalline undulator, Phys. Rev. Lett. 98 (2007) 164801.
[4] A.V. Korol, A.V. Solov'yov, W. Greiner, Channeling and Radiation in Periodically Bent Crystals, Springer-Verlag, Berlin, Heidelberg, 2013.
[5] V.T. Baranov, S. Bellucci, V.M. Biryukov, et al., First results of investigation of radiation from positrons in a crystalline undulator, JETP Letters. 82 (2005) 562-564.
[6] V.T. Baranov, S. Bellucci, V.M. Biryukov, et al, Preliminary results on the study of radiation from positrons in a periodically deformed crystal, Nucl. Instrum. Methods B. 252 (2006) 32-35.
[7] H. Backe, D. Krambrich, W. Lauth, et al., Channeling and Radiation of Electrons in Silicon Single Crystals and Si1—^Ge^ Crystalline Undula-tors, Journal of Physics: Conference Series. 438(1) (2013) 012017.
[8] u. Wienands, T. Markiewicz, J. Nelson, et al., Observation of a Remarkable Deflection of Multi-GeV Electron, SLAC Scientific Publications, SLAC-PUB-15952 (2014).
[9] I.A. Solov'yov, A.V. Yakubovich, P.V. Niko-laev, et al., MesoBioNano Explorer — a universal program for multiscale computer simulations of complex molecular structure and dynamics, J. Comp. Chem. 33 (2012) 2412—2439.
[10] http://www.mbnexplorer.com
[11] G.B. Sushko, V.G. Bezchastnov, I.A. So-lovyov, 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 (1) (2013) 012018.
[13] G. Sushko, V. Bezchastnov, A. Korol, et al., Simulation of electron channeling in bent silicon crystal, Journal of Physics: Conference Series. 438 (1) (2013) 012019.
[14] R.G. Polozkov, V.K. Ivanov, G.B. Sushko, et al., Radiation emission by electrons channeling in bent silicon crystals, Eur. Phys. J. D 68(9) (2014) 268
[15] V.G. Bezchastnov, A.V. Korol, A.V. Solovy-
ov, Radiation from multi-GeV electrons and positrons in periodically bent silicon crystal, Journal of Physics B: Atomic, Molecular and Optical Physics. 47(19) (2014) 195401.
[16] J. Bak, J.A. Ellison, B. Marsh, et al., Channeling radiation from 2-55 GeV electrons and positrons: (I) planar case, Nucl. Phys. B. 302 (1988) 525-558.
[17] O. Bogdanov, E. Fiks, K. Korotchenko, et
al., Basic channeling with Mathematica©: A new computer code, in: Journal of Physics: Conference Series, 236 (2010) 012029.
[18] B. Azadegan, S. Dabagov, W. Wagner, Computer simulation of the radiation of electrons axially channeled in a thin Ge single crystal. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms. 269(19) (2011) 2098-2106.
[19] S. Abdrashitov, O. Bogdanov, S. Dabagov, et al., Total yield of channeling radiation from rela-tivistic electrons in thin Si and W crystals, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms. 309 (2013) 59-62.
[20] G. Molière, Theorie der Streuung schnel-
ler geladener Teilchen I: Einzelstreuung am abge- [21] D.S. Gemmel, Channeling and related ef-
schirmten Coulomb-Feld. Z. f. Naturforsch. A 2 fects in the motion of charged particles through (1947) 133-145. crystals, Rev. Mod. Phys. 46 (1947) 129-227.
THE AUTHORS
SUSHKO Gennady B.
MBN Research Center at Frankfurt Innovation Center of Biotechnology Altenhoferallee 3, 60438 Frankfurt am Main, Germany sushko@mbnexplorer.com
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
korol@fias.uni-frankfurt.de
SOLOV'YOV Andrei V.
MBN Research Center at Frankfurt Innovation Center of Biotechnology Altenhoferallee 3, 60438 Frankfurt am Main, Germany solovyov@mbnresearch.com
Сушко Г.Б., Король A.B., Соловьев A.B. распространение электронов и
позитронов в прямых и периодически изогнутых плоскостных и аксиальных каналах в монокристаллах кремния.
В статье представлены результаты численных расчетов плоскостного и аксиального типов кана-лирования электронов и позитронов в прямых и периодически изогнутых каналах в монокристаллах кремния. Расчет траекторий частиц учитывал их взаимодействия с соседними атомами кристалла и выполнялся с помощью пакета программ MBN Explorer. Использованный подход позволил сравнивать характеристики плоскостного и аксиального типов каналирования. Результаты расчетов показали, что длина деканалирования и число каналирующих частиц заметно ниже в аксиальном случае. Для последнего изучена зависимость характеристик процесса от выбора направлений оси и изгиба кристалла.
электрон, позитрон, каналирование, кристалл кремния, периодически изогнутый канал, траектория.
СПИСОК ЛИТЕРАТУРЫ
[1] Lindhard J., Dan K. Stopping power of electron gas and equipartition rule // Vidensk. Selsk. Mat. Fys. Medd. 1965. Vol. 34. Pp. 1-64.
[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] Tabrizi M., Korol A.V., Solov'yov A.V., Greiner W. Feasibility of an electron-based crystalline undulator // Phys. Rev. Lett. 2007. Vol. 98. Pp. 164801 (1-4).
[4] Korol A.V., Solov'yov A.V., Greiner W. Channeling and Radiation in Periodically Bent Crystals. Springer-Verlag, Berlin Heidelberg, 2013.
[5] Baranov V.T., Bellucci S., Biryukov V.M., et al. First results of investigation of radiation from
positrons in a crystalline undulator // JETP Letters. 2005. Vol. 82. Pp. 562-564.
[6] Baranov V.T., Bellucci S., Biryukov V.M., et al. Preliminary results on the study of radiation from positrons in a periodically deformed crystal // Nucl. Instrum. Methods B. 2006. Vol. 252. Pp. 32-35.
[7] Backe H., Krambrich D., Lauth W., et al.
Channeling and radiation of electrons in silicon single crystals and Si^Ge^ crystalline undulators // Journal of Physics: Conference Series. 2013. Vol. 438. Pp. 012017.
[8] Wienands U., Markiewicz T., Nelson J., et al. Observation of a Remarkable Deflection of Multi-GeV Electron, SLAC Scientific Publications, SLAC-PUB-15952 (2014).
[9] Solov'yov I.A., Yakubovich A.V., Nikolaev
P.V., et al. MesoBioNano Explorer — a universal program for multiscale computer simulations of complex molecular structure and dynamics // J. Comp. Chem. 2012. Vol. 33. Pp. 2412—2439.
[10] http://www.mbnexplorer.com
[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., et al. Sub-GeV electron and positron channeling in straight, bent and periodically bent silicon crystals // Journal of Physics: Conference Series. 2013. Vol. 438. Pp. 012018.
[13] Sushko G., Bezchastnov V., Korol A., et al. Simulation of electron channeling in bent silicon crystal // Journal of Physics: Conference Series. 2013. Vol. 438. Pp. 012019.
[14] Polozkov R.G., Ivanov V.K., Sushko G.B., et al. Radiation emission by electrons channeling in bent silicon crystals // Eur. Phys. J. D. 2014. Vol. 68. Pp. 268.
[15] 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 Phys-
ics. 2014. Vol. 47. Pp. 195401.
[16] Bak J., Ellison J.A., Marsh B., et al.
Channeling radiation from 2-55 GeV electrons and positrons: (I) planar case // Nucl. Phys. B. 1988. Vol. 302. Pp. 525—558.
[17] Bogdanov O., Fiks E., Korotchenko K., et al. Basic channeling with Mathematica©: A new computer code // Journal of Physics: Conference Series. 2010. Vol. 236. Pp. 012029.
[18] Azadegan B., Dabagov S., Wagner W. Computer simulation of the radiation of electrons axially channeled in a thin Ge single crystal // Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms. 2011. Vol. 269. Pp. 2098—2106.
[19] Abdrashitov S., Bogdanov O., Dabagov S., et al. Total yield of channeling radiation from rela-tivistic electrons in thin Si and W crystals // Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms. 2013. Vol. 309. Pp. 59—62.
[20] G. Molière. Theorie der Streuung schneller geladener Teilchen I: Einzelstreuung am abgeschirmten Coulomb-Feld // Z. f. Naturforsch. A. 1947. Vol. 2. Pp. 133—145.
[21] Gemmel D.S. Channeling and related effects in the motion of charged particles through crystals // Rev. Mod. Phys. 1947. Vol. 46. Pp. 129—227.
СВЕДЕНИЯ ОБ АВТОРАХ
СУШКО Геннадий Борисович — научный сотрудник Научно-исследовательского центра мезобиона-носистем (MBN); г. Франкфурт-на-Майне, Германия. Altenhoferallee 3, 60438 Frankfurt am Main, Germany sushko@mbnexplorer.com
КОРОЛь Андрей Владимирович — кандидат физико-математических наук, доцент кафедры физики Санкт-Петербургского государственного морского технического университета, научный сотрудник Научно-исследовательского центра мезобионаносистем (MBN); г. Франкфурт-на-Майне, Германия. 198262, Российская Федерация, г. Санкт-Петербург, Ленинский пр., 101. korol@fias.uni-frankfurt.de
СОЛОВьЕВ Андрей Владимирович — доктор физико-математических наук, профессор, директор Научно-исследовательского центра мезобионаносистем (MBN); г. Франкфурт-на-Майне, Германия. Altenhoferallee 3, 60438 Frankfurt am Main, Germany solovyov@mbnresearch.comая ул., 7 larionov ivanvl@mail.ru
© Санкт-Петербургский политехнический университет Петра Великого, 2015
