CC BY
27
UDC 539.1.03
V Y. Epp, M. A. Sosedova
RADIATION OF THE ATOMIC CHAINS EXITED BY A CHANNELING PARTICLE
Basic properties of radiation of the atomic chains excited by a channeling particle are considered. Using a very simple two-dimensional model of a crystal lattice we have shown that the main part of this radiation is generated on the frequency of oscillations of a channeling particle between the crystal planes, shifted by the Doppler effect. Angular distribution of the radiation of the chain of oscillating atoms is sharply peaked in the direction of the velocity of channeling particle because of coherence of the fields, produced by individual atoms.
Key words: radiation, channeling particle, crystal lattice, angular distribution, coherence.
1. Introduction
Physics of channeling of the accelerated particles in crystal, is a quickly developing field of science [1]. The effect of channeling has served as the base of development of new experimental methods of research on the crystal structure. By use of the channeling effect it is possible to study thermal oscillations and displacement of atoms in a lattice, distribution of electronic density in interatomic space of crystals and to measure exact orientation of the crystal planes. Channeling of light particles is used for generation of intensive monochromatic X-ray radiation. Using experiments and the theories which describe the orientation effects in crystals, a new sources of X-ray and the gamma radiations consisting of accelerators and precisely oriented crystals are created.
In the course of interaction between the channeling particle and a crystal the electromagnetic radiation is generated. In the framework of classical electrodynamics one can consider this radiation as radiation from different sources - radiation of the channeling particle oscillating between the atomic planes or around an crystal axis, radiation of the electronic gas excited by a channeling particle (weik fields of electrons) and the radiation which is produced by the excited atoms of a crystal lattice. Radiation of a channeling particle is investigated in details (see for example [2, 3]), the experimental research of radiation of wakefield trace can be found in [4]. Radiation of the atoms excited by a channeling particle has not been studied so far. In this article we consider part of this radiation caused by oscillations of the atoms embedded at specific sites in the lattice of a crystal, leaving aside the radiation caused by quantum transitions of electrons from excited to principal states. We suppose that the conditions of implementation of classical electrodynamics are fulfilled. Namely, the energy of excitation of atom is much more than distance between the energy levels; we consider the part of a radiation spectrum of single atom where the energy of photons is much less than energy of excitation of atom. Specific feature of considered radiation is that the phases of oscillations of atoms correlate between themselves as oscillations are excited by the same
channeling particle. Therefore radiation should be considered as coherent.
2. Dynamics of atoms in a crystal lattice
Let us consider very simplified model of a crystal lattice which allows, nevertheless, to find out the basic properties of radiation of the excited atomic chains. We will consider a crystal as a two-dimensional lattice of right-angled structure. Let the relativistic, positively charged particle having a charge q1 enter in a crystal between two atomic chains under small angle to the crystal axes. We consider the energy of the particle much above the energy loss on radiation and interaction with the crystal. Then the velocity of the charged particle in the crystal remains to be constant and equal to the initial velocity V .
Now we are interested in oscillations of the atoms at the lattice points which are caused by collision with the channeling particle, and in the radiation related to these oscillations. In spite of the fact that the atom, as a whole, is neutral, only part of its electrons is bound to the atom. The peripheral electrons belong to the whole crystal. Colliding with an atom, a channeling particle interacts with the screened potential of the atom nucleus, disturbing it from the equilibrium. Oscillations of a nucleus with interior electrons cause time-dependent polarization of atom and hence electromagnetic radiation. At calculation of interaction of a particle with individual atom we assume that the channeling particle moves parallel to atomic chains with a constant velocity. Rather slow oscillations of a particle between atomic planes we will take into account by the distance between a particle and nucleus of atoms, which is slowly changing along a trajectory of channeling particle. The model of a crystal lattice considered here has been used in [5] for calculation of polarization of atomic chains induced by the channeling particle.
Let the axis OX of a coordinate system be parallel to the atomic chains and lying in the middle between the neighboring atomic chains of a crystal (Fig. 1). We denote the distance between the next atoms lying in a chain by b and the distance between the atomic chains by 2D. ”
Fig. 1. Model of a two-dimensional crystal lattice
Then the atom with number i of the first chain has coordinates (xt, D), and the opposite atom of the second chain has coordinates (xi, -D). We choose the origin of coordinates system so that x0 = 0 .
The field of a relativistic channeling particle in the laboratory reference system is concentrated basically in a plane which is perpendicular to a direction of motion [6].
Hence, the particle interacts effectively with the atom during a short interval of time when the atom is in the vicinity of this plane. One can expect that this interaction will be reduced to the momentum transfer from the particle to atom in a direction orthogonal to the velocity of channeling particle. In order to prove it. we find x and y components of the momentum transmitted to the atom of a lattice. The electric field of a charged particle moving with constant velocity is set by the formula [6]:
E =
(i -P2)
(1)
smf :
R
yla2 + V 2t2
(2)
Time t is equal to zero when the particle is at the minimal distance from a nucleus. The atom is acted on by the field of moving positively charged particle with the Coulomb’s force
F = Eq2, (3)
where q2 is an effective charge of a nucleus of atom which is screened by the interior electrons. It is obvious that atoms of a crystal lattice are neutral, but a channeling particle approaching the lattice point of a crystal appears to be inside the atom and interacts therefore with the nucleus partially screened by the in-
terior electrons. Consequently, an effective charge q2 is less than the charge of a nucleus.
One can see from the formula (1) that the field of an ultrarelativistic particle (P ~ 1) has important quantity in a vicinity sin y ~ 1 or more exactly cos y ~ y-1,
where y-1 = yjl - p2 is the relativistic factor. Then, it follows from (2) that interaction of a particle with a nucleus occurs during time:
a a ...
t ^ «_. (4)
V y cy
Integrating coordinate components of formula (3) over time t e [-<», ®], we find components of momentum Px and Py which is transmitted to atom of chain along the axes OX and OY respectively:
P =
Mi
aV
Vi-p7,
P = q1q2
y
aV
It follows from the ratio
(5)
(6)
P.
1
:Y>
R (-p2sin2^)
where R is a vector pointing from the charged particle to the atom, P = V / c, c is the velocity of light, y is an angle between the volocity of the moving particle and vector R . Vector n is a unit vector of direction defined by n = R / R . The angle y depends on parameter a which is the shortest distance between the trajectory of channeling particle and a nucleus of atom.
Px V1 -P2
that the momentum transferred to the atom in Y direction is much greater then that in X direction. Hence, in case of an ultrarelativistic particle motion of the atom along the axis OX can be neglected and only oscillations along the axis OY can be considered. The passing particle excites harmonic oscillations of the atom at frequency © which is defined by crystal properties. For the moment we consider oscillations of atom as free ones, attenuation of these oscillations we shall take into account later. The law of motion of the i-th atom of the first chain can be written as:
/
= A1
sin ro
v
ib
t------
V
\
t > *. V
(7)
Accordingly, the law of motion for the atom of the second chain will be written as:
A ib V
yi
' sin ©
t----I-D
ib t >—, V
(8)
where A/(j = 1,2) are constants defined by initial conditions. The origin of coordinates is accepted so that t = 0 when the channeling particle hits the atom number i = 0. The amplitude of oscillations is defined by an impulse transmitted to atom. It is obvious that the laws of motion (7) and (8) should be completed by equation
ib
= 0,
t <—.
V
Now, in order to calculate the amplitude of oscillations we suppose that the particle transmits part of its momentum to the atom in no time. Such assumption is valid, if effective time of interaction of a particle with
an atom is much less than period of its own oscillations. From the formula (4) follows that this condition takes form
©a — < Y, c
which obviously is fulfilled for a nonrelativistic atom as ©a is of the same order of magnitude as the velocity of atom oscillating with amplitude a. The amplitude of oscillations we find from initial condition that the amplitude of momentum of the atom m&Aj is equal to momentum (6) accepted from the particle. As a result we have:
A - , (9)
mra Vai
where aj is the shortest distance between the particle and i-th atom of chain number j. As the velocity of atom is the nonrelativistic, one can use dipole approximation at calculation of radiation of an oscillating charge. We will combine two atoms which are in different chains opposite each other and are simultaneously excited by a passing particle into one dipole. The dipole moment of obtained system is equal to
pi = (2 Pyi,0), where
Py(t ) = <h(4- 42)sin ®
/
V
ib
t------
V
(10)
yc = K cos
Qx
___c_
V
Accordingly, parameters a/ for the atoms having coordinates xi = ib, are equal to
, Qx o Qx.
a. = D - K cos—, a. = D + Kcos—. (11)
' V ' V
Substituting these parameters in formulas (9) and
(10), we obtain amplitude of oscillations of i-th dipole
(10)
A, = q( 4-A;) = cos
mVD © V
(12)
Let us find dependence of parameter a/ on coordinate x. Motion of a channeling particle in case of small amplitude of oscillations (it means the amplitude of oscillations in Y axis direction is much less than interatomic distance 2D) to a good approximation can be described by the equations:
= K cos(Qt + ), xc = Vt,
where yc and xc are particle coordinates, K and ^0 is amplitude and initial phase. The initial phase can be putted equal to zero by corresponding choice of the origin of X axis1. Then the equation of trajectory of a channeling particle takes the form:
And the net result for the dipole moment generated by two atoms of the neighbor chains will be written as:
Pyi (x;, t) = Ai sin t - V j, pxi (X;) = 2xtqv (13)
In principle, the form of expression (13) was obvious from the beginning. All previous reasoning has been directed to calculation of amplitude of oscillations (12). At this point we can correct the primitive model of a crystal lattice to some extend setting the amplitude Ai from some empirical data. For example, we have not set the quantity of screened charge of nucleus q2, have not considered thermal oscillations of atom etc. Some of such corrections can be made by inserting of corresponding coefficients in formula (12).
The most essential improvement we are going to make at this step - is to take into account the attenuation of atom oscillations. The main cause of attenuations is the transmission of the energy of oscillations to surrounding atoms and to a lattice as a whole. In terms of quantum mechanics this process is described by radiation of phonons. This process is studied, for example, in [7]. By assuming that the energy of phon-on is much less than the energy of oscillations we can consider attenuation of oscillations by putting an exponential factor into amplitude of oscillations. With account of attenuation the equation for the dipole moment (13) becomes:
Pyi (Xi , t) =
-a(t)
Ae
sin ©
V
x. ) x.
t —11, t VJ V
(14)
where a is the attenuation coefficient. We will ignore the constant component of the dipole moment pxi as the radiation field is proportional to the second derivative on time from the dipole moment.
3. Radiation of the excited chains of atoms
Electric field of radiation of the dipole moment is defined by the formula [6]:
E'•(rj ) = ~cT ^ [n,p'(t ,
(15)
where the unit vector n is defined by equality n = r / r , r is the vector pointing from the dipole to the point for which the vector Et is written down, r =| r |, pt is to be evaluated at the retarded time t'
t = t —.
1 We have already defined the coordinate origin just after Eq. (8), bounding it to some specific atom. Now we can shift the origin by integer number of distance b. Note that b << V/ Q..
Further we assume that the distance between the crystal and the observer is much greater than sizes of the crystal. Then the vector n and distance r in Eq. (15) are constant and do not depend on x. Dependence on x should be considered only in the phase of an electromagnetic wave, namely in t'. We denote the vector from^ the coordinates origin to the point of observation by R, and radius-vector of the dipole moment by rp as shown in Fig. 2.
Fig. 2. System of coordinates
As R » rp, it is possible to write approximately:
Vc
r = \/(R - rp)
R - (nrp X
(16)
where n = R / R . Accordingly, for the retarded time t' we have:
R - xn
t' = t---------------------------------------------------x . (17)
c
Let us find the components of vector E(t) in a spherical system of coordinates shown in Fig. 3. Direction of the unit vector n is defined by the polar angle 0 and azimuthal angle 9. The Cartesian coordinates of the unit vectors of spherical coordinate system can be written as:
n = (sin 9 sin 9, cos 9, sin 9 cos 9), e9 = (cos 9 sin 9, - sin 9, cos 9 cos 9), h =(cos 9,0, - sin 9).
(18)
(19)
(20)
Then it follows from Eq. (15) that:
e
Ei (R, t)
Rc‘
-p yi (t >m 9.
(21)
Taking twice derivative from Eq. (14) with respect to time, we obtain
P yi(xi,t') =
Aie aT|^(a2-©2)sin©x-2a©cos©x], t'
>-
x
V‘
(22)
0, t'<-
V
where
T = f'-V = f - - - V (1 -$nx), P = L
V c V c
Let the atomic chain consist of large, but finite number of atoms N. The net field of radiation is the discrete sum of fields of individual atoms of the crystal
E (R, t ) = YEt (R, t).
(23)
As the distance between the atoms is much less than the wavelength of radiation, it is possible to turn to continuous distribution of the dipole moment along the X-axis, and replace the sum by integral. There is one dipole per interval of length b in atomic chain, hence, the quantity s(x, t) = Et6 (R, t) / b has the meaning of the field generated by unit length of the atomic chain. Substituting into this expression the field of Eq. (21) and taking into account Eq. (22), we obtain s (x, t) for a time interval where it is not equal to zero:
z(x,t) = A0e~az cos^|^(a2 -m2)smmx-2arncosmxj, (24) here
A _ 2q1q^K sin 9 _ R
c
mVD &bRc
V
(1 -p nx). (25)
Trigonometric functions in the square brackets can be combined, using the phase y:
2a©
sin y© = —2-------t-
a + ©
cos y©:
2 2
a -©
2 2 *
a + ©
As a result we get:
Qx
(26)
(27)
s(x,t) = A0(a + © )e aT cos —sin©(T - V).
Integrating this expression over x, we find the field generated by pair of atomic chains
2
E (t )e = Js(x, t)dx.
(28)
The lower limit x1 of integral is coordinate of the crystal origin. The upper limit x2 is a time function.
i=0
x
While the channeling the particle moves in a crystal, the upper limit coincides with coordinate of the particle according to Eq. (22). After the particle leaves the crystal, radiation is generated by all the chain of atoms, and x2 coincides with coordinate of the end of
chain x20:
x2{t ) =
(ct-, t<- + x(1 -Pn),
(1 -pnx) c V
R x _ .
X20> t + — (1 -Pnx ).
cV
(29)
Integration in Eq. (28) gives
E,(t) = -LA0(a2 + ®Va(t°-kx) x 2k
/
a sin(®tj + kxA) - A_ cos(ratj + kxA_) a2 + A2
-
a sin(ro?j - kxA+) + A+ cos(ra^ - kxA+)
+
a2 + A+
(30)
where
o
1 -$nx
words, we consider only the first term in Eq. (32) where x2 = t0 / k . This is the main term in case of a “long crystal”. Then the formula for a field of radiation becomes
E(t) = ^ 4>(a 2 +®2)x 2k
x\ sin O't
a cos ray - A_ sin ray
a2 + A2
a cos ray + A+ sin ray
a2 + A+
+
+ cos O't
-a sin ray _ A_ cos ray
a2 + A2
a sin ray _ A+ cos ray
a2 + A+
(33)
R R
t0 t , ti t ^,
c c
k = -1(1 -$nx), A± = n' + ffl. (31)
We see from this expression that radiation is generated on two frequencies, namely on frequency of atoms oscillations © and the frequency of oscillations of channeling particle which is shifted by the Doppler effect:
Eq. (30) can be written as:
Ee(t) = e~at° [F(x2,-F(xj,t)eahj ], (32)
where F (x, t) is a periodic function of time. Considering exponential factor in last formula, we see that while
R x „ „ „ t ^~ + ^7(1 -$nx). c V
the amplitude of the electric field of radiation is growing accordingly to
Ee 0) = F(x2, t)eafa2 - F(Xj, t)ea(kXj-tfl),
and approaching a constant level if the crystal is long enough. It means that the time it takes for the particle to cross the crystal is much greater then the lifetime of oscillating atoms 1/a. We shall refer such a crystal as a “long” one. After a particle leaves the crystal, the amplitude of field exponentially decreases. Let us next neglect the boundary effects and consider only radiation with constant amplitude of the field. In other
As we see in this case the radiation is generated only on frequency Q'. Let us find angular distribution of intensity of radiation. Intensity of radiation dl in an element of space angle do is defined by the formula [6]:
*=cREL. (34)
do 4n
Averaging expression (33) over the period of radiation and substituting in (34) we obtain
dl _ q1q2K2 sin2 0
do 2 nm2 b2 C D4(1 — p nx )2 4 O'2 a2 + (a2 + ra2)2 X ((a2 — ra2 +O'2)2 + 4a2 ra2)'
From last formula follows that dependence of in tensity of radiation on frequencies Q' and © has reso nant nature. At a reasonably small a the resonance oc curs at O' = ra or
(35)
O
= ra.
- (36)
1 -K
If the channeling particle is ultrarelativistic one, rather wide range of values of the relation ra / O' satisfies the resonance condition (36):
1 © ~ 2 -< — <2y2.
2 O'
4. Conclusion
In this work the basic properties of radiation of the atomic chains excited by a channeling particle are considered. On a simplified model of a crystal lattice it is shown that this radiation is generated on two frequencies, namely on the frequency of atom oscillations, and on the frequency of oscillations of a channeling particle, shifted by Doppler effect. In case of a reasonably long crystal the radiation with frequency of oscillations
x
of a channeling particle prevails. Occurrence of the factor (1 - Pnx) in denominator of the formula for angular distribution of intensity of radiation (35) results in relativistic effect - the basic part of radiation is generated into a cone around the direction of velocity of a relativistic channeling particle. Hence, the frequency of radiation and angular distribution are similar to corresponding characteristics of radiation of a channeling particle. This can rise a problem of experimental separation of one form of radiation from another one and a problem of identification of considered radiation.
The interesting effect showing that a set of not rela-tivistic oscillators at rest or oscillating atoms produces
radiation with the properties which are specific for radiation of the relativistic charged particle, is obviously a result of the fact that oscillations of atoms are coherent. Because phases of their oscillation are tuned by the channeling particle. The interference of fields of radiation of individual atoms gives the properties of the net radiation, described above.
Acknowledgment
The authors would like to thank Professor V. G. Ty-uterev for valuable discussions. This work has been supported by the grant for LRSS, project No 3558.2010.2.
References
1. Baryshevsky V. G. Channeling, Radiation and Reactions in Crystals at High Energy (in Russian). M.: Moscow University Press, 1982. P. 256.
2. Lindhard J. // Phys. Let. 1964. Vol. 12. P. 126.
3. Ahiezer A. I., Shulga N. F. Electrodynamics of High-Energy Particles in Matter (in Russian). M.: Nauka, 1993. P. 344.
4. Shchelkunov S. V., Marshall T. C., Hirshfield J. L. et al. Experimental observation of constructive superposition of wake fields generated by electron bunches in a dielectric-lined waveguide // Phys. Rev. Sc. and Technol., Accelerators and Beams 9, 011301 (2006).
5. Gogolev S. Yu. XII All-Russian conference of students, postgraduate students and young scientists “Science and education” (In Russian). Vol. 1. Science. P. 1. Physics and mathematics. Tomsk: TSPU, 2008. P. 243.
6. Landau L. D., Lifshitz E. M. The Classical Theory of Fields. Butterworth Heinemann. Vol. 2. 4th ed., 1994. P. 193.
7. Tyuterev V. G. Quasi-stationary regime of electron-lattice relaxation of highly excited electrons in dielectrics // Tomsk State Pedagogical University Bulletin. 2006. Issue 6. P. 7-10 (in Russian).
Epp V. Y.
Tomsk State Pedagogical University.
Ul. Kievskaya, 60, Tomsk, Russia, 634061.
Sosedova M. A.
Tomsk State Pedagogical University.
Ul. Kievskaya, 60, Tomsk, Russia, 634061.
Received 14.03.2011.
В. Я. Эпп, М. А. Соседова ИЗЛУЧЕНИЕ ЦЕПОЧКИ АТОМОВ, ВОЗБУЖДЕННЫх КАНАЛИРУЮщЕЙ ЧАСТИЦЕЙ
Рассмотрены основные свойства излучения атомных цепочек, возбужденных каналирующей частицей. Используя упрощенную модель двумерной кристаллической решетки, мы показали, что основная часть этого излучения генерируется на частоте колебаний каналирующей частицы, смещенной эффектом Доплера. Угловое распределение излучения цепочки колеблющихся атомов усиливается в направлении скорости каналирующей частицы из-за когерентности полей, создаваемых отдельными атомами.
Ключевые слова: излучение, каналирование частиц, кристаллическая решетка, угловое распределение, когерентность.
Эпп В. Я., доктор физико-математических наук, профессор.
Томский государственный педагогический университет.
Ул. Киевская, 60, Томск, Россия, 634061.
Соседова М. А., аспирант.
Томский государственный педагогический университет.
Ул. Киевская, 60, Томск, Россия, 634061.
