CC BY
33
UDC 539.1.076, 539.122.2, 538.97
Laser Acceleration of Channeled Light Particles in Periodically
Modulated Crystal
V. V. Serov, T. A. Sergeeva
Chair of Theoretical Physics Saratov State University 83 Astrakhanskaya str, Saratov 410026, Russia
An approach for the laser acceleration of the charged particles channeled in a periodically bent crystal have been investigated under the accounting for the transverse energy leveles quantization.
Key words and phrases: channeling, laser acceleration, metacrystallic structures.
1. Introduction
At the present time, particle accelerators are widely used in a number of researching fields both of the fundamental and applied science. However, the modern acceleration mechanisms are limited in terms of the practically attainable acceleration gradients. This limitation causes an enormous increase in accelerators length that is appearing to be a significant challenge. It gives a reason to investigate and develop some possible aproaches for channeled particles acceleration, alternative to conventional ones.
As early as at the end of the previous century a number of ideas for such alternative acceleration technologies has emerged [1]. Among the basic concepts that should be emphazised are the possibility to utilize a wave in plasma [2,3], and diverse ways to make use of superpower laser radiation [4,5]. The latter ones are namely the direct laser acceleration (DLA) in a lens focus [6], and the acceleration in photonic crystals (both two-dimensional [7] and three-dimensional [8]), and also the inverse free electron laser (IFEL) acceleration [9-11]. A number of these ideas has been successfully implemented in an experiment [12-14]. The approaches have in common the concept of creation of an extremely high power accelerating field by some means or other. On the other hand, these fields already exist within the crystal structures. It is concerned with the channeling phemonena that hence can be used in a particle accelerator.
Channeling is a phenomena of guiding of a charged particle along the major crys-tallographic directions called channels. Provided that the particle velocity is large enough and an angle between the motion direction and the crystallographic planes, or the atomic strings, is small enough (smaller than so called critical Lindhard angle [15,16]), it experiences the multiple scattering from the crystallographic planes or strings (in the cases of the planar or axial channeling, respectively). The particle interaction with the planes or strings of atoms in a crystal lattice nodes is described by the continuous averaged Lindhard potential [15,16]. A positively charged particle, being channeled in a crystal, moves in a region of the minimal electron density, consequently it penetrates at much larger distance than in the absence of the channeling.
An important feature of the channeling phenomena is the preserving of this motion mode in a slightly bent crystal. Therefore one of the channeling applications is the ability of the deviation of charged high energy particles with the help of the bent crystals. Notably, the deflection angle achieved by such a way could be much larger than its maximal possible magnitude enabled by the use of a magnetic field. The concept was first proposed in [17]. Subsequent works have shown the capability of this effect to be more widely used for the steering of the high energy particles beams, namely the beams focusing, splitting, collimation, extraction from the accelerators and even for the measurements of short-lived particle magnetic moments [18]. In the
Received 29th June, 2012.
This work was supported by the President of Russian Federation Grant No. MK-2344.2010.2 and Grant No. RFBR 11-01-00523-a.
number of cases it is the only possibility to carry out the experiments in the collider mode and with a fixed target as well.
The first concepts of the possibility to achieve any advantages in accelerators design by means of the channeling features usage were also suggested as far back as in 80th years of the last century [5,19,20]. Since a bent crystal may be used instead of a magnetic field for the charged particle deflection, an idea to utilize a crystal with the periodically bent channels as an undulator in a FEL have emerged rather long time ago. The next step in this technique development appeared to be the acceleration on the principle of IFEL with the aid of a crystalline undulator. That is exactly what has been put forward in the Bogacz work [21]. He suggested to use an effect of particles channeling in the periodically bent channels under the laser field having a polarization transverse to a channel direction for the purpose of particles acceleration. This method is based upon the laser pumpling of the particle transverse degree of freedom and its subsequent transfer to the longitudinal one due to a crystal periodical bending.
Since the Bogacz paper [21] publication techniques of metacrystallic structures production and utilizing have been developed extensively. In particular, a paradigm of the crystalline undulators construction in the form of the crystals with the channels periodically bent by means of the mechanical stress has already been realized experimentally [22] (it implies the deformation creating through the making use of the regular tranches on the surface of a thin crystal silicium plate). Other schemes for the obtaining of such structures are widely investigated [23]. Moreover, practically feasible magnitudes of the laser radiation power have increased significantly. In total, it gives a convincing reason to continue the theoretical research of the subject in this field.
In the Bogacz work [21] completely classical approach was used which is rather applicable in the context of the heavy particle channeling. But for the light ones the quantum effects, namely the transverse energy levels quantization, cannot be negliged in the course of the transverse degree of freedom considering. Meanwhile, the energy levels quantization may lead to the acceleration regimes discretization. The advantage can be taken from this in order to produce the coherent particle beams. In the current work a method of applying the quantum treatment for the transverse degree of freedom and the classical approach for the longitudinal one is explicated. This method will be used to examine a feasibility of the laser acceleration of a light charged particle channeled in a crystal with the periodically bent channels, and also a possibility of the particle filtering for the aid of the monochromatic beams extraction.
The work structure is as follows. In the section 2 the quantum-classical equations describing a particle motion in a crystal with the periodically bent channels under the complementary external laser field are derived from the Schrodinger equation. Next, in the section 3 we obtain in terms of the two level approximation the solution to these equations yielding a dependence of a particle acceleration on the parameters of the laser field and a channel bending, and also on the initial conditions. This is followed by the numerical calculation results illustrating the obtained theoretical estimations in the section 4. The section 5 contains the numerical evaluations of the above-stated acceleration parameters for the different kinds of particles and crystals. Also, the relevant possibilities of the proposed method application are discussed. Finally, in the conclusion the general summary is presented.
Atomic units will be used throughout the paper.
2. The quantum-classical approximation
The motion of a quantum particle is described by the Schrodinger equation
- 2ro V2 + U ^
(1)
where ^(r, t) is the particle wavefunction and U(r, t) is a potential affecting it. The latter consists of a crystal lattice potential and the external laser field summed
U(r, i) = Uc[p - Po(z)] + V(p, t), (2)
where p denotes a transversal coordinate, z being a longitudinal one, and U0(p) is an effective potential of a crystal channel, V(p, t) is a potential of the laser field, where Po(z) is a position of the center of the channel bent for instance by a shear sound wave. The potential might be rewritten as
U(r, t) = Uo(p) + V(p, t) + W(p, z), (3)
where Uo(p - p0(z)) - U0(p) ~ -(p)po(z) ^ -^(0)(po(z) • p) is a potential approximately linear over p.
Further, one can obtain the solutions to the stationary Schrodinger equation for the transversal coordinates having the form
h ±<pn(p) = £n<pn(p), where h± = V^ + U0(p)
2m
is a transversal Hamiltonian, en are the eigenenergies of the stationary transversal states, and <pn(p) are the wavefunctions of these states that may be used subsequently for the deriving of the required quantities.
Next, a problem appraisal could be made for the purpose of achieving any simplifications to the equation being solved. First, for the fast particles the quantum effects are likely to appear only for the transverse degree of freedom. In contrast, for the longitudinal one the wavelength magnitude is much less than the typical system size gauge, namely the channel width. As a consequence, it makes a sense to handle the longitudinal degree of freedom approximately. Specifically, it can be described in terms of the Gaussian wavepacket. In order to obtain the motion equations in the classical approximation for the particle's longitudal degree of freedom from the general Schrodinger equation, one may use the Petrov-Galerkin approach.
First, let us assume a trial function for the solution of the Eq.(1) to be expressed
as
^(r, t)
J
ax/K
(t) exp
ikn(t)[z - zn(t)\ - [z Of - i J En(t')df
<-Pn(P).
Here amplitudes of transversal levels population denoted by cn(t), longitudinal positions of the wave packet center labelled as zn(t), average longitudinal packet momenta written as kn(t) and effective energies En(t) are the unknown parameters. In order to evaluate them, the Petrov-Galerkin method will be applied. For this purpose, let us introduce the test functions in the following way:
$n0(r, t)
1
exp
ikn(t)[z - zn(t)] - [z 2a^2(t)]2 - i j E»(Odi'
pn (p);
t
t
$nl(r, t) =
{¿WK1
[ z - zn(t)]x
x exp
ikn(t)[z - zn(t)] - [z 2a"2(t)]2 - i J En(t')dt'
The Petrov-Galerkin method is confined in the orthogonalization of the discrepancy to the test function, namely
d
1 di -
i d2
- 2 ^ + h^(p) + V(P, f) + W(P, z)
= 0; n - 1, 2; a = 0,1.
(4)
Further, by collecting the terms before in equations from (4) obtained through
V2 a
the projection on $n1, one arrives to
¿n - kn .
In other words, as expected, the position of the wave packet "'center of mass"' and its average momentum are related to each other in the same manner as a position and a momentum of a classical particle. Besides, the term containing cn which is given by projection on $n0 might be removed from equations on the assumption that
En = - Y + 4a +en'
Next, as mentioned above, the method is applicable only on the condition that the positions of packets in the different transversal states are much less then the packet width, namely (zn< — zn ^ a. Furthermore, the difference of momenta Akn<n = kn — kn should also be small, viz A kn'n/kn < 1. Consequently, the approximated expressions for the potentials matrix elements can be obtained in this way:
Vnn' (z, k, t) - ($no|V |$n0> - (n|V |n' ) e-i (k»'z»'
Wnn(z, k, t) - ($n0|W|$no) - (n|Wln')(znn>)e- «v»^'-^»'*»'
where the denotations un/n(t) = En(t) — En(t) and znn' = (zn< + zn)/2 are used. In turn, the expanding of W(p, z) and exp(«Akn/nz) in Taylor series by z up to the first order yields the following:
<$ni|V|^no> ^ iAknn^(nlvK>e-ifiu»'»(t')dt'-i-knZn)
($nl|W|$n0> -
V2
d(n\W\nf> , ,Ai . /w_ ,
^ ( %nn' ) + 1 ^ kn'n {n W ^ >(Znn' )
az
o-i lo ^n'n (t')dt '-i (kn'zn>-knzn)
Having the terms of the same order of a collected, one arrives to the system of equations
c - -i[V(z, k, t) + W(z, k, t)]c,
Zn — kn,
(5)
(6)
| Cn|2k n c*nFnn' (z, k, t) Cn c*n[Vnn' (z, k, t) + Wnn' (z, k, t)]i(kn -kn) cn, (7)
where
Fnn' (z, k, t) = - (W) e-i 10 ^
oz
may be interpreted as a force affecting the particle due to the presence of the spatial field W.
The further problem simplification can then be achieved through the introducing of the average longitudinal position and the momentum of the particle
{^ ^cnzncn, {k) ^ ^cnkncn.
n n
Next, let us multiply the equation (6) by |cn\2 and subsequently perform the summation over n in the Eq.(6) and the Eq.(7). Thereupon, considered that zn & {z) and kn & {k) in the right hand sides of equations (6) and (5), the final result is
C = -i[V(i) + W({ z), t)] • c;
{*) = {k); (8)
{ k) = Y, cnFnW ({z), t)cn>,
n,n!
which appear to be much easier to solve than the common Schrodinger equation. Here the coefficients of the wavefunction expansion in the terms of the transversal states <pn(p) having a physical meaning of population amplitudes and denoted by cn, the wavepacket center position along the z axis labelled as { z), and { k) standing for the averaged longitudinal momentum of the particle are the system variables. Further, the system contains the following values:
Vnn,(t) = {n\V(p, t)\n')eiu™'(9)
Wnn,(z, t) = {n|W(p, z)\n')e™™'' (10)
being the matrix elements of the summand potentials from (3) and
Fnn' (,, t) = - (11)
o
is the one of a longitudinal force F defined as the force corresponding to the spatial potential affecting the particle. Besides, unn< = en - en' are the transition frequencies, n and n' stand for the stationary states wavefunctions <pn(p) referred to above.
3. Analitycal solution by means of the two-level
approximation
Provided that the potentials (9,10,11) explicit form is specified, the solution of the system (8) might be obtained. So, let us assume that
V(p, t) = Sop sin(ut); (12)
W(p, z) = -Vop sin(Kz), (13)
and also the fields amplitudes £0 fl Here u denotes the external field frequency and k means the periodical spatial potential wavenumber. In this case
Vnn' (t) = —Sodnn' sin(wt)1;
Wnn (z, t) = Vodnn sin(KZ)e^™'l; (14)
Fnn' (z, t) = —vVodnn* cos(kz)\
where dnn = (pn|x|pn> is the dipole transition matrix element, at that Ox ft £0. On condition that the channel potential U0(p) is even the diagonal element dnn = 0.
Next, consider the particular case of w = w2i. Granted that (k) is large and changes rather slightly, true is the expression z — ( k)t + z0. Thereupon, if ( k> = kres, where
kres = w2i/k, (15)
then the z-dependent potential affects the particle with an effective frequency xkres = w2i. In such instance the transitions between the levels 1 and 2 are the most probable, hence it is possible to use the two-level approximation. When the quickly oscillating parts of the matrix elements (14) are omitted one arrives to
Vn(t) - 21 ^21;
W2i(krest, t) - — 1 P0d2ie-iSa; F2i(krest, t) - — 1 Kp0d2ie-iSo,
where 50 = xz0 is the spatial potential initial phase. In that event, namely the two-level approximation, the solution for the population amplitudes is obtained in the form
c1(i) = cos Qt cos ^ — sin Qt sin ^ ez( S2-">); c2(t) = — cosQtsin ^ e%&2 — sinQtcos ^ e%1,
where
1
Q = 21(£0 —v0 e-s° )d2i|
is the Rabi frequency,
£0 — V0 e-i d2i
e"'
o%l -
|£o -T>oe-is° | |d2i|
is the effective potential phase factor, ¿i stands for an initial phase of Rabi oscillations, and ô2 denotes an initial phase shift between the population amplitudes. In particular, on the assumption that
£o - T>o, (16)
and ôo - 0, the Rabi frequency fi - 0. Consequently, the coefficients ci,2(i) are constant; therefore the longitudinal force
( F> - ^ C*nFnn' (krest, t)Cn' - ksin(2<5i) cos 02,
n, n
is constant too. It is easy to see that the maximal acceleration is achieved on the condition that ^ = k/4 and 62 = 0, i.e., when ci = — c2.
Now turn back to the more general case of Q = 0. The force might be averaged by time over a single Rabi period in the following manner:
2ir/n
Q
{F) — 2 J ^ ^ CnFnni (krest, i)cn> —
= KVld21 sin 281 [cOs(82 + ¿o) - COs(82 - ¿0 - 27)] . As long as the condition (16) is satisfied, 7 = (i - 80)/2 and thus
(p/ = K^d21 sin 2¿1 [cos(¿0 + ¿2) + cos ¿2] .
It is clear from the last expression that the average force is maximal in the event of ^ = i/4 and 62 = 50 = 0. Besides, its magnitude is equal to
—— = KPod21 \ / max o .
In such a manner, the conditions of resonance (15) and amplitudes equality (16) have been met, the equation for the particle momentum k from the system under consideration (8) can be derived in the form
k = X sin(2<51) [cos(82 + ¿0) + cos 82] .
Its solution is _
k(t) = + at;
22
Z(f) = 3a( k2 + at)3/2 - + Z0,
where denotation a = w12£0d21 sin(2<51) [cos(82 + 80) + cos 82] /4 is used. From these expressions one can obtain the momentum dependence on z as follows:
k( ) —
, 3 3 a, ,
k3 + y(z — Zo)
1/3
In turn, it may be substituted to the resonance condition (15). Therefore, in order to supply the resonance of the transversal transition with z-dependent potential its "wave number" should be equal to
/ n / 3 3a \ 1/3
k( z) = W21 I K + —z I ,
where k = k0 is the particle initial momentum. Thereat, instead of (13), the spatial potential should have the form
W(p, z) — -V0p sin I J k(z')dz— -V0p sin —21 ^k3 + ~az^ — k2
, (17)
where now k and a denote the metacrystallic structure parameters. Besides, it is easy to see that the parameter a should meet the condition 0 < a < u12D0d21/2 to provide the acceleration of particles with the initial momentum k0 = k. Thuswise we have obtained the potential yielding the resonance and, as a consequence, enabling the acceleration.
Now let us turn to the exploring of the resonance mode stability. To this end consider the phase difference between z-dependent potential and transversal oscillations; it might be expressed in this way:
ic(i) = ^
a
+ ^¿(i)) / -
W21Z0 K - U2it ~ —^ = ¿0(0). \fat y at
So, as long as 50(0) is small enough to keep the acceleration during rather long time period, S0 converges to zero and the acceleration becomes stable. Thus, the acceleration mode is stable while S0 detunes slightly from zero.
Finally, the energy of the particle having the maximal acceleration is equal to
E = £0 + i.
Thuswise, the maximal energy increase per a distance unit, i.e. the accelerating gradient, can be estimated as
d^=¿o. (18)
dz v/8£
4. Numerical calculations
In order to demonstrate the particle acceleration process, the system (8) was solved, the potentials have been assumed to have the forms (12), (17) proved above to provide the resonance mode. In the Fig.1 one can see the results of this numerical solution.
13,2-,
12,5-1-,-1-,-1-,-1-,-1
0 5000 10000 15000 20000
t, a.u.
Figure 1. The particle momentum dependence on time for the different initial conditions: = ^/4, S0 = S2 = 0 (solid curve), S0 = ^/2, = ^"/4, S2 = 0 (dashed curve), S0 = 0, = ^"/4, S2 = w/2 (dotted curve), S0 = ^"/2, = ^"/4, S2 = ^/3 (dash-dotted curve), S0 = ^"/2, = ^/12, S0 = S2 = 0 (short dotted curve)
The system parameters used in the calculations are as follows: d 12 = 1, w12 = 2^/50, S0 = V0 = 0.01, k = 4^, a = u12£0d21/2 at k(t = 0) = k and different initial conditions specified by the phases ¿0, ¿1, S2. The choice of the constant a standing for the nonuniformity of the spatial potential period provides the permanent resonance with the accelerating force magnitude being 2 times less than its maximal value. Therefore in the event of the initial conditions yielding the maximal force, namely S1 = ^/4, S0 = 62 = 0, the particle momentum periodically becomes larger
than its resonance value, as a consequence the resonance vanishes and the particle decelerates, whereat it enters back into the resonance mode. In sum, the acceleration 2 times less than the maximum possible is produced. The population density evolution for this case is displayed in the Fig.2.
0,60-
0,35-1-,-1-,-1-,-1-,-1
0 5000 10000 15000 20000
t, a.u.
Figure 2. The lower transversal level population dependence on time for the case of the maximal accelerating force
It is apparent that in proximity of the resonance point the population density curve nearly reaches the plateau, while at the momentum magnitude strongly different from its resonance value it oscillates. The Fig.1 also presents the curves for several sets of initial conditions yielding the average force magnitude equal to the half of its maximal value. In other words, our choice for a should enable the permanent resonance observed, just as it is clear from the linearity of the momentum increase over time. The plot shows as well the curve corresponding to the initial conditions producing the average force magnitude equal to zero, i.e., as a consequence, the absence of the acceleration. Hereby the acceleration mode stability has been ensured, besides it is provided by the fact that the force magnitude is greater than some value corresponding to the chosen a at the initial acceleration stage.
The Fig.3 presents the dependence of the probability of the affecting at the initial acceleration stage of the force {F), such that {F)/{F)max > f. At that the distribution over either of the phases 50, Si, 62 is assumed to be uniform. It is clear from the figure that the particle entering the crystal is rather likely to get into the regime such that the accelerating force magnitude would be close to the maximal possible value.
In the Bogacz work [21] the particle acceleration was investigated in the case that it is allowable to use the harmonic approximation for the channel transverse potential as well as the classical approach for the particle transverse motion. In terms of the quantum mechanics it implies a cascade of transitions between a number of equidistant transverse motion energy levels. In contrast, we examine the situation of the two-level approximation applicability. Therefore in this case one should take into consideration not only the resonance between the laser radiation and the periodically spatial potential which always occurs at the certain particle momentum. Also the resonance of the laser radiation with the frequency of the transition between the transverse levels requires taking into account.
The fig.4 demonstrates the evolution of the momentum of the particle channeled in a medium with a constant "wavelength" of the periodically spatial potential. At that, the particle is assumed to have an initial momentum
kres = w/k,
(19)
f
Figure 3. The probability for the particle to be affected by the initial conditions such that the normalized average accelerating force (F}/ (F}max > f
providing the resonance between the laser radiation and the periodically spatial potential. Presented are curves corresponding to the various frequency w detuning from the frequency of the tranverse transition The acceleration range width in the frequency space is apparent to have the same order as the width of the resonance broadening due to the acceleration itself.
(0
12,85 12,80 12,75 12,7012,65 12,6012,55 12,5012,45 12,40 -12,3512,30
/X
—I-1-1-
5000 10000
t.a.u.
15000
20000
0,005w12 -0,005co1;
-0,01to12 0,1co12 0,02œ12 -0,02co
Figure 4. The accelerated particle longitudinal momentum dependence on time in the presence of the resonance with the periodically spatial potential for the different laser radiation frequency detuning from the frequency of the tranverse transition: w = w 12 (thick solid curve), w = 0.005^12 (dashed curve), w = —0.005^12 (dotted curve), w = — 0.01w12 (dash-dotted curve), w = 0.1w12 (dash-dot-dotted curve), w = 0.02w12 (short dashed curve), w = —0.02w12 (short dotted curve)
All the previous cases have been considered under the neglecting of both ionizing and radiation energy losses of the particle travelling through the crystal. Meanwhile, these losses are necessary to be regarded in the common crystals. The paper [24] gives a half-emphirical dependence expression for the energy loss per the length unity for a positron passing through the matter. The channeling losses are some less than in the course of the propagation in a random direction. The value of the coefficient of the stopping power reducing in the channeling mode has been shown to be roughly equal to 0.7 [25]. So, in order to take the energy loss into account we have added in the last equation of the system (8) the stopping power calculated by the formula from [24] and additionally multiplied by 0.7.
The fig. 5 presents the accelerated particle longitudinal momentum dependence on the longitudinal coordinate meaning with provision for the stopping power for the constant "wavelength" of the periodically spatial potential and the different values of the external laser field intensity.
z, a.u.
Figure 5. The accelerated particle longitudinal momentum dependence on the longitudinal coordinate meaning with provision for the stopping power for the different values of the external laser field intensity: £0 = 1 (solid curve), £0 = 1/2 (dashed curve), £0 = 1/4 (dotted curve), £0 = 1/8 (dash-dotted curve)
At that, the phases initial meanings have been set to grant the maximal laser acceleration. It is clear that the laser acceleration can be observed only at very large S0 magnitudes yielding its efficiency greater than the stopping power. Besides, after a number of momentum oscillations with respect to the resonance value at a certain moment system inevitably transgress to the clean deceleration mode in cosequence of the stopping power presence. For the case of the laser acceleration force less than the stopping power this transgression occurs immediately. One may conclude that for the small acceleration force, contrary to the apparent expectations, the laser acceleration efficiency is not "substracted" from the stopping power, but simply disappears. This effect arises as a consequence of the rapid departing out of the resonance with the periodically spatial potential.
Further, in order to compensate the resonance vanishing by virtue of the stoping power presence, the variable "wavelength" of the periodically spatial potential have been introduced. This have been done in such a manner that the resonance would nevertheless remain in the presence of the stopping power and the infinitesimal laser field as well. That is, the periodically spatial potential "wavelength" decreases with the z increasing in such a way that if the particle entering the crystal had the resonance
momentum meaning then the momentum decreasing due to the stopping power would not result in the disappearance of the resonance with the periodically spatial potential. The relevant accelerated particle momentum dependencies on the distance from the crystal entering point for the different values of the external laser field intensity are displayed in the fig. 6.
z,a.u.
Figure 6. The accelerated particle longitudinal momentum dependence on the longitudinal coordinate meaning on the condition of the periodically spatial potential "wavelength" providing the continuous resonance in the absence of the laser field for the different values of the external laser field intensity: £o = 0.125 (dotted curve), £0 = 0.01 (dashed curve), £0 = 0 (solid curve)
The plots have been obtained on the conditions that the periodically spatial potential has a variable "wavelength" supporting the continuous resonance under the absence of the laser field and also the stopping power is present. It is apparent that for the S0 values not leading to the crystal destroy (see the sec. 5) the laser acceleration provides some increase in the particle penetration depth. However, since this increase is extremely small, the idea of the monochromatic positrons producing by means of the laser acceleration appears to be practically unrealizeable. Hereafter, for rather large S0 values (but still less than the meanings yielding the laser acceleration efficiency exceeding the stopping power) the momentum oscillations with respect to the resonance values arise, just as well as in the case of the constant periodically spatial potential "wavelength" in the absence of the stopping power.
5. Estimations
Now it is necessary to discuss a matter of the practical aspects of the suggested acceleration method implementation. It is known that the optical damage threshold for a Si crystal for 1550 nm laser pulse lasting less than 1 ps is Ith ^ 1011 W/cm2 [26], that corresponds to S0 = 1.7 x 10-3[a.u.]= 0.9 GeV/m. Next, let us assume the transversal dipole transition matrix element to be |d12| ~ 1. For a nonrelativistic proton in Si w0 ~ 10-2 a.u. In the case of the proton energy value being equal to
d^
E = 10 MeV the proton velocity is v = 20 a.u., that yields — = 0.1 MeV/m. It is negligible with regard to the stopping power of 10 MeV-protons that is on order
d^
of — 10 GeV/m. Moreover, this quantity of —— is small in comparison with the one
dz
yielded by an usual linear accelerator, which acceleration rate may be estimated as — 10 MeV/m.
For the channeled nonrelativistic positron u0 ~ 1 a.u., and the velocity value equal to v = 20 a.u. (and, as a conqequence, the energy value of E = 5.4 KeV) implies d^
-j— = 10 MeV/m. Though the latter energy increment is much greater than that for
the proton case, it appears to be still unsufficient in terms of the direct application of the considered effect for an accelerator designing on this framework.
Furthermore, the usage of carbon nanotubes instead of the Si crystals strongly reduces the stopping power (due to much larger channel radius values). The effective transversal potential of the tube can be approximated by a circular well potential. If one considers the two lowest transversal states in the tube, its radius R increase would lead to the transition frequency wi2 decreasing in accordance with wi2 — 1/R2. On the other hand, the dipole matrix element of the transition between thses states would increase as di2 — R. Thus, in sum the acceleration gradient (18) depends on the tube radius as
dE So
- r--j -.
dz mvR
However, since the stopping power is inversely related to the material density, true is the following:
dE ~dz
stop
1
hence there is no possibility to make the accelerating efficiency larger than the stopping power through the only increasing of the tube radius. Therefore the utilizing of the periodically bent carbon nanotubes instead of the crystals has the only advantage of the production easiness in contrast to a crystal with the periodically bent channels.
r^j
6. Conclusion
In this work we have studied an approach of the charged particles acceleration in a crystal with the periodically bent channels through the use of the laser radiation field. The concept has been originally suggested by Bogacz [21], and developed here for the case of the transverse energy levels quantization being significant. The technique as it stands have been demonstrated to be inefficient for the purpose of the particles acceleration upon using either the Si crystals or the carbon nanotubes. This appeared to be true at least for the laser radiation intensities meanings not causing a medium destroy. The considered scheme turned out to be nonapplicable in terms of its usage with the aim of the monochromatic particle beams obtaining as well.
Nevertheless, the scientific field of the particle acceleration, particularly using the channeling special features, is widely researched at the time. It comprises a number of perspectives and possibilities of its implementation. The development of particle acceleration approaches, especially non-conventional innovative ones, appears to be of crucial importance. So, further thorough investigations to the branch are greatly desireable.
References
1. Carrigan R. A. J. Studying and Applying Channeling at Extremely High Bunch Charges // Nuclear Instruments and Methods in Physics Research, Section B. — 2005. — Vol. 234. — Pp. 116-121.
2. Tajima T., Dawson J. M. Laser Electron Accelerator // Physical Review Letters. — 1979. — Vol. 43, No 4. — Pp. 267-270.
3. Esarey E., Sprangle P. Overview of Plasma-Based Accelerator Concepts // IEEE Transactions on Plasma Science. — 1996. — Vol. 24, No 2. — Pp. 252-288.
4. Shimoda K. Proposal for an Electron Accelerator using an Optical Maser // Applied Optics. — 1962. — Vol. 1, No 1. — Pp. 33-35.
5. Palmer R. A Laser-Driven Grating Linac // Particle Accelerators. — 1980. — Vol. 11. — Pp. 81-90.
6. Esarey E., Sprangle P., Krall J. Laser Acceleration of Electrons in Vacuum // Physical Review E. — 1995. — Vol. 52, No 5. — Pp. 5443-5453.
7. Lin X. E. Photonic Band Gap Fiber Accelerator // Physical Review Special Topics - Accelerators and Beams. — 2001. — Vol. 4. — P. 051301.
8. Cowan B. M. Three-Dimensional Dielectric Photonic Crystal Structures for Laser-Driven Acceleration // Physical Review Special Topics - Accelerators and Beams. — 2008. — Vol. 11. — P. 011301.
9. Palmer R. Interaction of Relativistic Particles and Free Electromagnetic Waves in the Presence of a Static Helical Magnet // Journal of Applied Physics. — 1972. — Vol. 43, No 7. — Pp. 3014-3024.
10. Courant E. D., Pellegrini C., Zakowicz W. High-Energy Inverse Free-Electron-Laser Accelerator // Physical Review A. — 1985. — Vol. 32, No 5. — Pp. 28132823.
11. Observation of Energy Gain at the BNL Inverse Free-Electron-Laser Accelerator / A. Van Steenbergen, J. Gallardo, J. Sandweiss, J.-M. Fang // Physical Review Letters. — 1996. — Vol. 77, No 13. — Pp. 2690-2623.
12. Esarey E., Schroeder C. B., Leemans W. P. Physics of Laser-Driven Plasma-Based Electron Accelerators // Reviews of Modern Physics. — 2009. — Vol. 81, No 3. — Pp. 1229—1285.
13. Demonstration of High-Trapping Efficiency and Narrow Energy Spread in a Laser-Driven Accelerator / W. D. Kimura, M. Babzien, I. Ben-Zvi et al. // Physical Review Letters. — 2004. — Vol. 92, No 5. — P. 054801.
14. Geddes C. G. R. et al. Laser Plasma Particle Accelerators: Large Fields for Smaller Facility Sources // SciDAC Review. — 2009. — Vol. 13. — Pp. 13—21.
15. Линдхард Й. Влияние кристаллической решетки на движение быстрых заряженных частиц // Успехи Физических Наук. — 1969. — Т. 99, № 2. — С. 249296. [Lindkhard Yj. Vliyanie kristallicheskoyj reshetki na dvizhenie bihstrihkh zaryazhennihkh chastic // Uspekhi Fizicheskikh Nauk. — 1969. — T. 99, No 2. — S. 249-296. ]
16. Кумахов М. А., Ширмер Г. Атомные столкновения в кристаллах. — М.: Атом-издат, 1980. — 192 с. [Kumakhov M. A., Shirmer G. Atomnihe stolknoveniya v kristallakh. — M.: Atomizdat, 1980. — 192 s. ]
17. Таратин А. М, Цыганов Е. Н., Воробьев С. А. Поворот пучков заряженных частиц изогнутым монокристаллом. Численный эксперимент // Письма в ЖТФ. — 1978. — Т. 4. — С. 947-950. [Taratin A. M, Cihganov E. N., Vorobjev S. A. Povorot puchkov zaryazhennihkh chastic izognutihm monokristallom. Chislennihyj ehksperiment // Pisjma v ZhTF. — 1978. — T. 4. — S. 947-950. ]
18. Biryukov V. M., Chesnokov Y. A., Kotov V. I. Crystal Channeling and its Application at High Energy Accelerators. — Berlin: Springer, 1997. — 219 p.
19. Carrigan R. A. J., Ellison J. A. NATO Science Series B: Physics. V. 165. Relativistic Channeling. — NY: Plenum Press, 1987. — 527 p.
20. Carrigan R. A. J. Developments in Relativistic Channeling // Proceedings of the 7th Workshop Advanced Accelerator Concepts / Ed. by S. Chattopadhyay. — 1996.
21. Bogacz S. Inverse FEL Proton Accelerator via Periodically Modulated Crystal Structure // Particle Accelerators. — 1993. — Vol. 42, No 3-4. — Pp. 181-197.
22. S. B. et al. Experimental Study for the Feasibility of a Crystalline Undulator // Physical Review Letters. — 2003. — Vol. 90, No 3. — Pp. 034801-034804.
23. Korol A. V., Solovyov A. V., Greiner W. Channeling of Positrons through Periodically Bent Crystals: on Feasibility of Crystalline Undulator and Gamma-Laser // International Journal of Modern Physics E. — 2004. — Vol. 13, No 5. — Pp. 867-916.
24. Pal P. B., Varshney V. P., Gupta D. K. Semiempirical Stopping Power Equation for Positrons // Journal of Applied Physics. — 1986. — Vol. 60, No 1. — Pp. 461— -463.
25. Tielens T. M. H. E., Gerrit B. E. W., Stoof T. H. Physics of Laser-Driven Plasma-Based Electron Accelerators // Physical Review B. — 1994. — Vol. 49, No 8. — Pp. 5741—5744.
26. Cowan B. M. Optical Damage Threshold of Silicon for Ultrafast Infrared Pulses // SLAC-PUB. — 2007. — Vol. 12960. — Pp. AARD-493.
УДК 539.1.076, 539.122.2, 538.97
Лазерное ускорение легких частиц, каналированных в периодически модулированных кристаллах
В. В. Серов, Т. А. Сергеева
Кафедра теоретической физики Саратовский государственный университет им. Н. Г. Чернышевского 410026, Россия, Саратов, Астраханская, 83
Рассмотрен метод лазерного ускорения заряженных частиц, каналированных в периодически искривлённых кристаллах, с учётом дискретности уровней поперечных энергий.
Ключевые слова: каналирование, лазерное ускорение, метакристаллические структуры.
