Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators
The results of analytical studies and numerical simulations of wakefield excitation in a plasma-dielectric structure are presented. In the linear theory approximation (overdense plasma) it is shown that at a certain plasma density the superposition of the plasma wave and the dielectric waves allows...
Збережено в:
Дата: | 2013 |
---|---|
Автори: | , , , , , |
Формат: | Стаття |
Мова: | English |
Опубліковано: |
Національний науковий центр «Харківський фізико-технічний інститут» НАН України
2013
|
Назва видання: | Вопросы атомной науки и техники |
Теми: | |
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/111929 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators / R.R. Kniaziev, O.V. Manuilenko, P.I. Markov,T.C. Marshall, I.N. Onishchenko, G.V. Sotnikov // Вопросы атомной науки и техники. — 2013. — № 4. — С. 84-89. — Бібліогр.: 15 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-111929 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1119292017-01-16T03:04:00Z Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators Kniaziev, R.R. Manuilenko, O.V. Markov, P.I. Marshal, T.C. Onishchenko, I.N. Sotnikov, G.V. Новые методы ускорения заряженных частиц The results of analytical studies and numerical simulations of wakefield excitation in a plasma-dielectric structure are presented. In the linear theory approximation (overdense plasma) it is shown that at a certain plasma density the superposition of the plasma wave and the dielectric waves allows the acceleration of the witness bunch with simultaneous focusing. Also, we carried out a PIC simulation of the underdense (“blowout”) and overdense regime of wakefield excitation in the unit. Представлені результати аналітичних досліджень і чисельного моделювання збудження кільватерних полів у плазмово-діелектричної структурі. У лінійному наближенні (надщільна плазма) показано, що при певній щільності плазми суперпозиція плазмової і діелектричної хвиль дозволяє прискорювати тестовий згусток з його одночасної фокусуванням. Також ми виконали моделювання методом "частинка в комірці" збудження кільватерних полів для випадків розрідженій і надщільної плазми. Представлены результаты аналитических и численных исследований возбуждения кильватерных полей в плазменно-диэлектрической структуре. В линейном приближении (сверхплотная плазма) показано, что при определенной плотности плазмы суперпозиция плазменной и диэлектрической волн позволяет ускорять тестовый сгусток с его одновременной фокусировкой. Также мы выполнили моделирование методом “частица в ячейке” возбуждения кильватерных полей для случаев разреженной и сверхплотной плазмы. 2013 Article Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators / R.R. Kniaziev, O.V. Manuilenko, P.I. Markov,T.C. Marshall, I.N. Onishchenko, G.V. Sotnikov // Вопросы атомной науки и техники. — 2013. — № 4. — С. 84-89. — Бібліогр.: 15 назв. — англ. 1562-6016 PACS: 41.75.Ht, 41.75.Lx, 41.75.Jv, 96.50.Pw, 533.9. http://dspace.nbuv.gov.ua/handle/123456789/111929 en Вопросы атомной науки и техники Національний науковий центр «Харківський фізико-технічний інститут» НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
topic |
Новые методы ускорения заряженных частиц Новые методы ускорения заряженных частиц |
spellingShingle |
Новые методы ускорения заряженных частиц Новые методы ускорения заряженных частиц Kniaziev, R.R. Manuilenko, O.V. Markov, P.I. Marshal, T.C. Onishchenko, I.N. Sotnikov, G.V. Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators Вопросы атомной науки и техники |
description |
The results of analytical studies and numerical simulations of wakefield excitation in a plasma-dielectric structure are presented. In the linear theory approximation (overdense plasma) it is shown that at a certain plasma density the superposition of the plasma wave and the dielectric waves allows the acceleration of the witness bunch with simultaneous focusing. Also, we carried out a PIC simulation of the underdense (“blowout”) and overdense regime of wakefield excitation in the unit. |
format |
Article |
author |
Kniaziev, R.R. Manuilenko, O.V. Markov, P.I. Marshal, T.C. Onishchenko, I.N. Sotnikov, G.V. |
author_facet |
Kniaziev, R.R. Manuilenko, O.V. Markov, P.I. Marshal, T.C. Onishchenko, I.N. Sotnikov, G.V. |
author_sort |
Kniaziev, R.R. |
title |
Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators |
title_short |
Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators |
title_full |
Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators |
title_fullStr |
Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators |
title_full_unstemmed |
Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators |
title_sort |
focusing of electron and positron bunches in plasma-dielectric wakefield accelerators |
publisher |
Національний науковий центр «Харківський фізико-технічний інститут» НАН України |
publishDate |
2013 |
topic_facet |
Новые методы ускорения заряженных частиц |
url |
http://dspace.nbuv.gov.ua/handle/123456789/111929 |
citation_txt |
Focusing of electron and positron bunches in plasma-dielectric wakefield accelerators / R.R. Kniaziev, O.V. Manuilenko, P.I. Markov,T.C. Marshall, I.N. Onishchenko, G.V. Sotnikov // Вопросы атомной науки и техники. — 2013. — № 4. — С. 84-89. — Бібліогр.: 15 назв. — англ. |
series |
Вопросы атомной науки и техники |
work_keys_str_mv |
AT kniazievrr focusingofelectronandpositronbunchesinplasmadielectricwakefieldaccelerators AT manuilenkoov focusingofelectronandpositronbunchesinplasmadielectricwakefieldaccelerators AT markovpi focusingofelectronandpositronbunchesinplasmadielectricwakefieldaccelerators AT marshaltc focusingofelectronandpositronbunchesinplasmadielectricwakefieldaccelerators AT onishchenkoin focusingofelectronandpositronbunchesinplasmadielectricwakefieldaccelerators AT sotnikovgv focusingofelectronandpositronbunchesinplasmadielectricwakefieldaccelerators |
first_indexed |
2025-07-08T02:58:55Z |
last_indexed |
2025-07-08T02:58:55Z |
_version_ |
1837045945201590272 |
fulltext |
ISSN 1562-6016. ВАНТ. 2013. №4(86) 84
FOCUSING OF ELECTRON AND POSITRON BUNCHES
IN PLASMA-DIELECTRIC WAKEFIELD ACCELERATORS
R.R. Kniaziev1, O.V. Manuilenko2, P.I. Markov2,
T.C. Marshall3, I.N. Onishchenko2, G.V. Sotnikov2
1V.N. Karazin Kharkov National University, Kharkov, Ukraine;
2National Science Center “Kharkov Institute of Physics and Technology”, Kharkov, Ukraine;
3Columbia University, New York; USA
E-mail: sotnikov@kipt.kharkov.ua
The results of analytical studies and numerical simulations of wakefield excitation in a plasma-dielectric struc-
ture are presented. In the linear theory approximation (overdense plasma) it is shown that at a certain plasma density
the superposition of the plasma wave and the dielectric waves allows the acceleration of the witness bunch with si-
multaneous focusing. Also, we carried out a PIC simulation of the underdense (“blowout”) and overdense regime of
wakefield excitation in the unit.
PACS: 41.75.Ht, 41.75.Lx, 41.75.Jv, 96.50.Pw, 533.9.
INTRODUCTION
Using of dielectric structures for acceleration of
charged particles by the wakefields excited in them by
relativistic electron bunches, is actively developing di-
rection of new methods of acceleration [1, 2]. The ex-
perimental researches carried out, in ANL and SLAC
confirmed [3, 4] performability of this method of accel-
eration of charged particles. Now the dielectric wake-
field accelerator is considered as the promising candi-
date for future electron-positron colliders of TeV range
of energy [5].
Despite shown both theoretically, and experimen-
tally possibilities of obtaining a high acceleration rates,
one problem, that is not solved completely, remains -
difficulties with stabilization of transverse motion of the
drive and accelerated bunches and, thereof, obtaining
the accelerated bunches of particles with a small emit-
tance. In this article possibility of using for this purpose
of plasma filling the drift channel of dielectric structure
is considered. Such plasma can be created as a result of
the capillary discharge in a dielectric tube [6]. Plasma
use for focusing of an accelerated bunch is not the new
proposition. Focusing properties of plasma were inves-
tigated in PWFA both in the linear condition [7, 8] and
in a non-linear regime [9, 10]. But in the linear condi-
tion the peak of an accelerating field corresponds the
zero focusing field, and in a non-linear regime the re-
gion of acceleration is localized only near a drive bunch
because of a destruction of a non-linear plasma wave.
As we will demonstrate below, using of plasma-
dielectric structure allows avoiding these restrictions.
1. STATEMENT OF THE PROBLEM
Let's consider the dielectric waveguide of an annular
cross-section surrounded with a metal sheath of radius
b , the inner radius of the dielectric bushing is equal to
a . A dispersion of the dielectric tube permittivity is
absent and is equal to dε . The interior of such dielectric
waveguide r a< (the accelerating channel) is com-
pletely filled with isotropic plasma of density pn . The
solid monoenergetic electron bunch with initial velocity
0v and charge 0Q starts being injected into the acceler-
ating channel, parallel to a waveguide axis, at some
time. The electron bunch builds up electromagnetic
fields which can be used for increasing of energy of
accelerated bunches in so called wake method of accel-
eration of the charged particles. The object of the pre-
sent research is finding of longitudinal (accelerat-
ing/braking) forces, and also transverse forces which are
responsible for transverse dynamics of drive and driven
bunches. These forces will be used further at an simula-
tion of a motion of test charged particles.
2. LINEAR APPROXIMATION
In the linear approximation we will assume that ions
of plasma are immobile, and perturbations of plasma
electrons, caused by motion of electrons of plasma, are
linear. We will neglect the change of speed of bunch
electrons, also as thermal motions of electrons of plas-
ma and a bunch. In view of a linearity of plasma its dis-
persion properties can be described an dielectric permit-
tivity 2 21p pε ω ω= − , 2 24p pe n mω π= , where e− and
m are the charge and the mass of electron, pn is a non-
perturbed density of plasma, ω is an angular frequency.
The set of Maxwell equations, describing excitation of
the axially-symmetrical electromagnetic fields by elec-
tron bunches in a plasma-dielectric waveguide of a cy-
lindrical configuration, looks like:
( ) , ,
1 1, ,
1 1 4 , ( ) ,
r z r
z
z r z r z
H HE E D
z r c t z c t
DrH j D E d
r r c t c
ϕ ϕ
ω
ϕ
π ε ω ω
∞
−∞
∂ ∂∂ ∂ ∂
− = − − =
∂ ∂ ∂ ∂ ∂
∂∂
= + =
∂ ∂ ∫
(1)
where
( ), if ,
( )
, otherwise
p
d
r aε ω
ε ω
ε
<⎧
= ⎨
⎩
, ,r zEω is Fourier trans-
form from ,r zE , zj is a current density of bunches.
Let's find at first the electromagnetic field excited by
annular infinitely thin in axial direction bunch. For such
bunch the current density zj is described by expression:
0
0 0( ) ( )
2z
Qj r r t
r
δ δ τ
π
= − − , (2)
where 0t z vτ = − , 0r is a bunch radius, 0t is a arrival
time of bunch into waveguide ( 0z = ),δ is the Dirac
delta function.
ISSN 1562-6016. ВАНТ. 2013. №4(86) 85
For finding of the solution of the set of equations
(1), (2) we will use well-known procedure (see, for ex-
ample, [8, 11]). It consists in a decomposition of the
equations (1), (2) in a Fourier integrals over time t and
the longitudinal coordinate z , solution of obtained or-
dinary differential equations set relatively a variable r
by method of partial regions and a performing of an
inverse Fourier transform with use of the calculus of
residues. As a result we will obtain:
L d C
z z z zE E E E= + + , (3)
L d C
r r r rE E E E= + + , (4)
d CH H Hϕ ϕ ϕ= + . (5)
Components of fields with the "d" index correspond
to a wakefield of a dielectric waveguide, they are exist
in case of the vacuum drift channel ( )0pn = also, being
modified when filling with its plasma. They look like:
0 0 0 02 ( ) ( ) ( ) cos ( )d s s
z z z s
s
E E t e r e r tτ ω τ= − Θ − −∑ , (6)
0 0 0 02 ( ) ( ) ( )sin ( )d s s
r z r s
s
E E t e r e r tτ ω τ= Θ − −∑ , (7)
0 0 0 02 ( ) ( ) ( )sin ( )d s s
z s
s
H E t e r h r tϕ ϕτ ω τ= Θ − −∑ , (8)
where eigenfunctions s
ze , s
re , shϕ are defined as follows:
0
1/2
0
0
0
( )
, if
( )
( )
'( ) ( , ) , otherwise
( , )
s
p
s
ps
z s s
s s d d
s s
d d
I r
r a
I aae r
D F r b
F a b
κ
κ
ω ω κ κ
κ κ
⎧
≤⎪
⎛ ⎞ ⎪= ⎨⎜ ⎟
⎝ ⎠ ⎪
⎪⎩
, (9)
1
1/2 2
00
1
2
00
( )1 ,
( )1 ( )
( )
'( ) ( , )1 ,
( , )1
s
p
s
pp ss
r s s
s s d d
s s
d dd
I r
r a
I aae r
D F r b r a
F a b
κ
κβ ε ω
ω ω κ κ
κ κβ ε
⎧ −
≤⎪
−⎪⎛ ⎞
= − ⎨⎜ ⎟
⎝ ⎠ ⎪ >⎪ −⎩
,
(10)
0( ) ( ) ( ).s s
s rh r e rϕ β ε ω= (11)
In expressions (6)–(11) the notations are used:
1/22
0 01 ( )s
p p s s vκ β ε ω ω⎡ ⎤= −⎣ ⎦ , ( )1/22
0 01s
d d s vκ β ε ω= −
[ ]0 0( , ) ( 1) ( ) ( ) ( ) ( )n
n n nF x y J x Y y Y x J y= − − , nJ and nY
are the Bessel and the Weber functions of n -th order,
nI are the modified Bessel function of n -th order, Θ is
the Heaviside function; 2
0 02E Q a= , 0 0v cβ = ; sω
are eigenfrequencies of a dielectric waveguide which
are defined from the dispersion equation:
( )
1 1
0 0
( ) ( , )( ) 0
( ) ,
p p d d d
p p d d d
I a F a bD
I a F a b
ε κ ε κ κω
κ κ κ κ κ
≡ + = . (12)
'( )sD ω in expressions (9) - (10) is a derivative of dis-
persion function ( )D ω with respect to frequency ω ,
calculated at sω ω= , a root of equation (12).
Components of fields with the "C" index in expres-
sions (3) - (5) describe a quasistatic field of a bunch.
Components of quasistatic fields are equal:
0 0 0 0sgn( ) ( ) ( ) exp[ | |]C s s
z z z s
s
E E t e r e r tτ ω τ= − − − −∑ , (13)
0 0 0( ) ( ) exp[ | |]C s s
r z r s
s
E E e r e r tω τ= − −∑ , (14)
0 0 0( ) ( ) exp[ | |]C s s
z s
s
H E e r h r tϕ ϕ ω τ= − −∑ , (15)
where:
0
1/2
0
0
0
( )
, if
( )
( )
'( ) ( , ) , else
( , )
s
p
s
ps
z s s
s s d d
s s
d d
J r
r a
J aae r
D r b
a b
κ
κ
ω ω κ κ
κ κ
⎧
≤⎪
⎛ ⎞ ⎪= ⎨⎜ ⎟
Δ⎝ ⎠ ⎪
⎪Δ⎩
, (16)
1
1/2 2
00
1
2
00
( )1 , if
( )1
( )
'( ) ( , )1 , else
( , )1
s
p
s
pps
r s s
s s d d
s s
d dd
J r
r a
J aae r
D r b
a b
κ
κβ ε
ω ω κ κ
κ κβ ε
⎧
≤⎪
−⎪⎛ ⎞
= ⎨⎜ ⎟
Δ⎝ ⎠ ⎪
⎪ Δ−⎩
,
(17)
0( ) ( ) ( )s s
s rh r e rϕ β ε ω= (18)
and
1/22
0 01s
p p s vκ β ε ω⎡ ⎤= −⎣ ⎦ ,
1/22
0 01s
d d s vκ β ε ω⎡ ⎤= −⎣ ⎦ ,
0 0( , ) ( ) ( ) ( 1) ( ) ( )n
n n nx y I x K y K x I yΔ = − − , nK are
Macdonald functions; sgn(x) is a function, equal to 1, 0
and –1, if x>0, x=0 and x<0, correspondingly;
2 21p p sε ω ω= + , sω are eigenvalues which are defined
from the equation:
( )
1 1
0 0
( ) ( , )( ) 0
( ) ,
p p d d d
p p d d d
J a a bD
J a a b
ε κ ε κ κω
κ κ κ κ κ
Δ
≡ − =
Δ
. (19)
'( )sD ω in expressions (16)–(17) is a derivative of func-
tion ( )D ω with respect to ω , calculated at sω ω= .
It should be noted that the equation (19) can be ob-
tained from the equation (12) if in the last to make re-
placement iω ω= , and to impose a condition
2 1/2
0 0 0 0, (1 )s pω ω β γ γ β −> = − . (20)
In consequence of this condition, for ultrarelativistic
electron bunch and high densities of plasma (namely
such cases are of interest for perspective wake accelera-
tors), the quasistatic fields (13) - (15) very quickly fall
down from a bunch and are small inside a bunch. The
numerical analysis of these fields is given in other pa-
per [12].
And at last, components with the "L" index in ex-
pressions (3) - (4) describe plasma (Langmuir) wake-
field with a frequency pω ω= which is localized in the
plasma channel. The electric field of a plasma wakefield
is described by expressions:
2
0 0 0
0 0 0 0 0
0 0 0 0 0
2 cos[ ( )] ( )
( , ) ( ) ( ) , if
( ) ( , ) ( ) , if
L
z p p
p p p p
p p p p
E Q k t t
k a k r I k r I k a r r
I k r k a k r I k a r r a
ω τ τ= − − Θ −
Δ ≤⎧
×⎨ Δ < ≤⎩
, (21)
2
0 0 0
0 0 1 0 0
0 0 1 0 0
2 sin[ ( )] ( )
( , ) ( ) ( ) , if
( ) ( , ) ( ) , if
L
r p p
p p p p
p p p p
E Q k t t
k a k r I k r I k a r r
I k r k a k r I k a r r a
ω τ τ= − Θ −
Δ <⎧
×⎨− Δ < <⎩
. (22)
The field of Langmuir wave has no component of a
magnetic field, unlike the wakefield of dielectric wave-
guide (6) - (8) where the magnetic field decreases trans-
ISSN 1562-6016. ВАНТ. 2013. №4(86) 86
verse defocusing force by 2 1
0[1 ( )]p sβ ε ω −− times. As the
plasma wave period in the general case does not coin-
cide with the period of a dielectric wakefield, the plas-
ma wave gives possibility for independent focusing of
accelerated particles. At that we assume that accelera-
tion is provided by the wakefield of a dielectric wave.
Expressions (6) - (22) describe the electromagnetic
field excited by a thin annular bunch. To obtain the
fields excited by a bunch of the finite sizes, it is neces-
sary to integrate these expressions over ring’s arrival
time into the waveguide 0t and by ring’s positions with
the corresponding charge-density distribution function.
Let's consider the simple, square profile of a charge
density of a bunch in the longitudinal and transverse
directions. Then resultant expressions look like:
0 ||2 ( ) ( ) ( )d s s
z s b z
s
E E R r e r τ= − Ψ∑ , (23)
02 ( ) ( ) ( )d s s
r s b r
s
E E R r e r τ⊥= Ψ∑ , (24)
02 ( ) ( ) ( )d s s
s b
s
H E R r h rϕ ϕ τ⊥= Ψ∑ , (25)
0 ||( ) ( ) ( )C s C
z s b z
s
E E R r e r τ= Ψ∑ , (26)
0 ( ) ( ) ( )C s C
r s b r
s
E E R r e r τ⊥= Ψ∑ , (27)
0 ( ) ( ) ( )C s C
s b
s
H E R r h rϕ ϕ τ⊥= Ψ∑ , (28)
0
1
00
||
0
0
0
( )1 ( , ),
( )4 ( )
( )
( , ),
( )
p
p b p b
p b pL p
z
p bb b
p p b
p
I k r
k r k a r r
k r I k aQE
I k rr L
k a k r r r a
I k a
τ
⎧
− Δ <⎪
⎪= − Ψ ⎨
⎪ Δ < <⎪⎩
(29)
1
1
00
1
1
0
( )
( , ),
( a)4 ( )
( )
( r, a),
( a)
p
p b p b
pL p
r
p bb b
p p b
p
I k r
k r k a r r
I kQE
I k rr L
k k r r a
I k
τ⊥
⎧
Δ <⎪
⎪= − Ψ ⎨
⎪ Δ < <⎪⎩
(30)
In expressions (23)-(30) functions ||Ψ and ⊥Ψ de-
scribe longitudinal profile of excited fields:
, 0
|| ,
,
,
( ) sin( ) ( )
sin ( ) ( ) ,
p s
p s
p s b
p s b b
v
L
τ ω τ τ
ω
ω τ τ τ τ
⎡Ψ = Θ⎣
⎤− − Θ − ⎦
(31)
[ ]0
|| ( ) exp( | |) exp( | |)C
s s b
s b
v
L
τ ω τ ω τ τ
ω
Ψ = − − − − , (32)
, 0
,
,
,
( ) (1 cos ) ( )
(1 cos ( )) ( ) ,
p s
p s
p s b
p s b b
v
L
τ ω τ τ
ω
ω τ τ τ τ
⊥ ⎡Ψ = − Θ⎣
⎤− − − Θ − ⎦
(33)
[
]
0( ) sgn( )(1 exp( | |))
sgn( )(1 exp( | |) ,
C
s
s b
b s b
v
L
τ τ ω τ
ω
τ τ ω τ τ
⊥Ψ = − −
− − − − −
(34)
0b bL vτ = , bL is a bunch length, and functions sR and
sR define a transverse form-factor of a solid bunch of
radius br :
1/2
1
0
( )2( ) ,
'( ) ( )
s
p b
s b s s
p b s s p
I raR r
r D I a
κ
κ ω ω κ
⎛ ⎞
= ⎜ ⎟
⎝ ⎠
(35)
1/2
1
0
( )2( )
'( ) ( )
s
p b
s b s s
p b s s p
J raR r
r D J a
κ
κ ω ω κ
⎛ ⎞
= ⎜ ⎟
⎝ ⎠
. (36)
Let's now proceed to numerical investigations of the
obtained expressions (23) - (30). First of all the possible
rate of acceleration (the value of longitudinal force be-
hind a leading bunch) and possibility of focusing (value
of transverse force) of accelerated bunch is in the interest.
The numerical calculations in the gigahertz fre-
quency range of excited wakefield waves have shown
that it is rather easy to place an accelerated bunch in a
focusing phase. In this case a focusing is provided by a
field of the plasma wave (30) in which there is no com-
ponent of a magnetic field, and acceleration is provided
by a field of eigenmodes of a dielectric waveguide (23).
One of such variants is given in the paper [11]:
1.1cma = , b = 4.3 cm, 2,1dε = , 10 -3m10 cpn = , en-
ergy of bunch electrons is 5 MeV, 0 0.32 nCQ = ,
1.7cmbL = , bunch cross-section area is 20.91cm− . The
vacuum wavelength resonant with a bunch of the fun-
damental mode of dielectric structure is equal to
~ 11cm, and plasma wavelength is ~ 33 cm. Just essen-
tial difference in lengths of two types of waves provides
possibility of focusing of an accelerated bunch at its
corresponding placing behind a drive bunch.
The case investigated in paper [11] corresponds to
an approximation of the linear plasma, as 0.04b pn n = .
As the field frequency increases, it becomes more and
more difficult to fulfill the strong inequality 1b pn n .
The reason consists in the unequal growth rate of den-
sity of electrons in a bunch bn and plasma density pn .
Really, the frequency of dielectric mode grows in in-
verse proportion to thickness of a dielectric tube
( )~ 1s b aω − . (37)
The electron bunch sizes must be less than a half of
wavelength therefore the required bunch density at de-
creasing of the sizes of structure grows by the cubic law
( )3~ 1bn b a− . (38)
At the same time the density of electrons of plasma
changes in inverse proportion to a square of the size of
the accelerating channel
2~ 1pn a . (39)
The condition (39) follows from necessity to work
near a maximum of amplitude of a transverse electric
field of a plasma wave that is carried out under condi-
tion ~ 1paω [8, 13].
If to refuse the requirement of fulfillment the strong
inequality 1b pn n , it is possible to provide an focus-
ing of an accelerated bunch also in more high-frequency
range, than it is investigated in paper [12]. Longitudinal
and transverse profiles of forces for one of such possible
cases are shown in Fig. 1 and Fig. 2.
For the results of calculations presented in these fig-
ures the quartz tube ( )3.75dε = with an outer radius
b =0.6 mm and inner radius b =0.5 mm was used. The
energy of bunch electrons was 5 GeV, the bunch charge
was 3 nC, its length 0.2bL = mm, the bunch radius
0.45br = mm. The sizes of a bunch and its charge give
ISSN 1562-6016. ВАНТ. 2013. №4(86) 87
density of electrons in it 141.47 10bn = ⋅ cm-3. Plasma
density used in calculations 14.414 10pn ⋅= cm-3. In that
way, 1 3b pn n = . In the following section we will show
that this dimensions is sufficient for validity of the re-
sults of the linear calculations given in Figs. 1-2. For the
presented calculations the wavelength of the fundamen-
tal dielectric mode is equal to ∼1 mm, and a plasma
wavelength is 1.6 mm.
Fig. 1. Axial profile of the axial force (black line) and
axial profile of transverse force (red line) at the
distance r=0.45 mm from waveguide axis. Drive bunch
(yellow rectangle) moves from right to left. Cyan
rectangle shows possible location of electron witness
bunch and green rectangles show possible location
of positron witness bunch. Plasma density p bn 3n=
Fig. 1 shows axial dependence of the longitudinal
force, acting on the test particle. It follows from the de-
pendence, given in Fig. 1, that we can ensure accelera-
tion of charged particles with their simultaneous radial
focusing by placing the testing bunch at some distance
from the drive bunch head. As it can be seen in the Fig-
ure, the radial force almost harmoniously depends on
the axial coordinate with the period of approximately
0.16 cm, i.e. the Langmuir wave makes the greatest con-
tribution into the radial force. At the same time, its con-
tribution into the axial force, accelerating test particles,
is predominantly small. The axial force is predomi-
nantly determined by the eigen modes of the dielectric
waveguide; its complex behavior from the axial coordi-
nate is caused by excitation several radial modes of the
dielectric waveguide.
Fig. 2 shows the radial dependence of axial and
transverse forces, acting on the test particle, placed in
the first of the maximums of the accelerating field, at
the distance of 2.3 mm behind the drive bunch head.
The axial force changes insufficiently in the transport
channel cross section, while the radial force remains
focusing along all the channel section.
For the used in analytical calculations parameters of
the dielectric waveguide, bunch and plasma, the focus-
ing force amplitude is approx. 300 MeV/m, which
equals the focusing magnetic field induction ~1T. The
acceleration gradient, G > 200 MeV/m, is attractive, and
also is the prospect of possibly stabilizing the motion of
both drive and witness bunches, hitherto a difficulty for
the single-channel DWA [14].
Fig.2. Transverse profile of the longitudinal force
(black line) and transverse one (red line), acting
on a witness electron, located at a distance of 2.3 mm
( 0 0v v t zξ τ≡ = − ) from the head of the drive bunch
It should be noted the undoubted advantage of the
considered plasma-dielectric wakefield accelerator: it
can be used both for acceleration of electron bunches,
and for acceleration of positron bunches. From Fig. 1 it
is seen that for this purpose accelerated electron and
positron bunches must have different delay relative to
drive bunch. Focusing will be provided for both bunch-
es. This property of charge symmetry is topical one for
perspective electron-positron colliders.
3. NUMERICAL SIMULATIONS
The analytical expressions and the calculations
which have been carried out with using them and have
been presented in previous section are true for the ap-
proximation of linear plasma b pn n (overdense re-
gime). As pointed above, this strong inequality difficult
to carry out in high frequency range of the excited wake
fields. Therefore the terahertz example given of previ-
ous section which used for the calculation the simple
inequality / /b pn n 1 3 1= < requires further verification
of the validity of the numerical results. With this object
in view we have made full particle-in-cell (PIC) nu-
meric simulation of wakefields excitation in the plasma
dielectric structure under investigation. For calculations
we used both our own PIC code and XOOPIC code real-
ized for Linux [15]. The results of modeling, made with
both codes, coincide well.
In Fig. 3 the longitudinal distributions of transverse
and longitudinal forces acting on a test electron, calcu-
lated by XOOPIC code are presented. For numerical
simulation the same parameters of structure and a
bunch, as in the previous section were used. The input
and output ends of the waveguide of length 8l = mm
was short-circuited by the conductive planes, i.e. the
boundary condition 0,| 0r z lE = = was used. Electrons and
the ions (hydrogen) of plasma striking the dielectric
surface were moved away from calculation domain.
Comparison of the curves in Fig. 3 with correspond-
ing curves in Fig. 1 confirms acceptable coincidence of
the results of PIC modeling at self-consistent account of
the plasma dynamics and analytical results.
ISSN 1562-6016. ВАНТ. 2013. №4(86) 88
Fig. 3. Axial profile of the axial force (black line) and
axial profile of transverse force (red line)
at the distance r=0.45 mm from waveguide axis at time
t=26.688 psec. Drive bunch moves from left to right,
its head is located at z=8 mm. Distance (x-axis) is
measured in mm, forces (y-axis) are in MeV/m
The accelerating and focusing fields coincide quite
well. The greatest difference is observed in defocusing
areas of the wake field. This difference can be explained
for by the pushing out of plasma electrons to the dielec-
tric surface. The remaining plasma ions reduce the de-
focusing field. This is confirmed by configuration space
of plasma electrons (see Fig. 4) and configuration space
of ions. Outer boundary of plasma electrons is deeply
modulated with space period equal to plasma wave-
length. At that time outer boundary of ions remains
smooth. Outer boundary of electrons is deeply modu-
lated with space period equal to plasma wavelength. At
the same time outer boundary of ions remains smooth.
Thus near the dielectric surface the periodic cavities,
having total positive charge, arise. These positive
charged cavities decrease defocusing force acting on
test electrons. It should be noted that in this simulation
the loss of plasma electrons was not very high, about
~8%.
Fig. 4. Configuration space (z,r), of plasma electrons
(yellow dots) at time t=26.688 psec. Dielectric tube
is shown by orange color
Fig. 5 shows axial profiles of the accelerating and
focusing forces created with the drive bunch of 1nC
charge (ratio b pn n 1 9 1= ). As it followed from
Fig. 4 we observe a more exact coincidence with the
analytical calculations (taking into account the normali-
zation for the bunch charge), given in Fig. 1.
Another, in comparison with Figs. 3, 4, extreme case
b pn n (underdense or blowout regime [9, 10]) with a
focusing provided by plasma ions, remaining in the
transport channel after plasma electrons have been
pushed out of it by the intense drive bunch.
Fig. 5. The same as in Fig. 3 only for the drive bunch
charge equal to 1nC
In Fig. 6 longitudinal profiles of longitudinal and
transverse forces with a case 30b pn n = are shown.
For these calculations the quartz tube with an outer ra-
dius 0.5b = mm and inner radius 0.2a = mm was used.
Radius of an electron bunch was equal to 0.2 mm, rest
parameters correspond to Figs. 1 - 4. Density of plasma
used in calculations, given in Fig. 6, 1410pn = cm-3.
Fig. 6. Axial profile of the axial force (black line) and
axial profile of transverse force (red line) in the case of
blowout excitation regime of plasma-dielectric wave-
guide: b pn n =30. Forces are recorded at the distance
r=0.1 mm from waveguide axis at time t=44.54 psec.
Drive bunch moves from left to right, its location is
shown green rectangle
As follows from Fig. 6 the focusing force is almost
axially homogeneous behind the drive bunch and is
equal ~ 120 MeV/m. Longitudinal force weakly
changed in comparison with a vacuum case, i.e. it, as
expected, is formed by eigen waves of a dielectric tube.
If to increase plasma density to the value corresponding
to Figs. 1 - 4, it is possible to expect amplitude of the
focusing force ~ 530 MeV/m, that corresponds to an
induction of a focusing magnetic field ~1.8 T.
At last, we will consider transverse motion of the
test electrons accelerated by a wakefield wave. In Fig. 7
are shown the positions of trailing electrons during
moving under action of longitudinal and transverse
forces presented in Figs. 1, 2. Initial longitudinal posi-
tions of test electrons are within the electron bunch
shown in Fig. 1 ( )2.2 0.1mmz = ± . In a transverse di-
rection the test electrons are uniformly distributed from
0 to br .
ISSN 1562-6016. ВАНТ. 2013. №4(86) 89
Fig. 7. Radii of witness bunch electrons during
the acceleration by wakefield of drive bunch. Initially
(at s=0) center of witness bunch is located at the
distance of 2.2 cm behind of drive bunch head
From Fig. 7 follows that for the distance ~10 cm the
accelerated electron bunch is focused almost twice. At
that a spread of the transverse coordinates of particles in
the head and in the tail of bunch is less 12%.
ACKNOWLEDGEMENTS
This work was supported by the us department of
energy/nnsa through the global initiatives for prolifera-
tion prevention (gipp) program in partnership with the
science and technology center in ukraine (project anl-t2-
247-ua and stcu agreement p522).
REFERENCES
1. Wei Gai. Advanced accelerating structures and their
interaction with electron beams // AIP Conf. Proc.
№ 1086 (AIP, New York, 2009), p. 3-11.
2. Eric R. Colby. Present limits and future prospect for
dielectric acceleration // Proc. of 35th Int. Conf. on
high energy physics (Paris, France, 2010).
3. W. Gai et al. Experimental demonstration of dielec-
tric structure based two beam acceleration // AIP Conf.
Proc. № 569 (AIP, New York, 2001), p. 287-293.
4. M.C. Thompson et al. Breakdown limits on giga-
volt-per-meter electron-beam-driven wakefields in
dielectric structures // Phys. Rev. Lett. 2008, v. 100,
p. 214801.
5. W. Gai, J.G. Power, and C. Jing. Short-pulse dielec-
tric two-beam acceleration // J. Plasma Physics.
2012, v. 78, № 4, p. 339-345.
6. L.C. Steinhauer and W.D. Kimura. Quasistatic capil-
lary discharge plasma model // Phys. Rev. ST Accel-
erator and Beams 9, 2006, p. 081301.
7. R.D. Ruth, A.W. Chao, P.L. Morton, P.B. Wilson. A
plasma wake field accelerator // Particle Accelera-
tors. 1985, v. 17, p. 171-189.
8. V.A. Balakirev, N.I. Karbushev, A.O. Ostrovsky,
Yu.V. Tkach. Theory of Cherenkov amplifiers and
generators based on relativistic beams. Kiev: «Nau-
kova dumka», 1993, p. 161-165 (in Ukrainian).
9. J.B. Rosenzweig, B. Breizman, T. Katsouleas, and
J.J. Su. Acceleration and focusing of electrons in
two-dimensional nonlinear plasma wake fields //
Phys. Rev. A. 1991, v. 44, № 10, p. R6189-R6192.
10. N. Barov and J. B. Rosenzweig. Propagation of short
electron pulses in underdense plasmas // Phys. Rev.
E. 1994, v. 49, № 5, p. 4407-4416.
11. R.R. Kniazev, G.V. Sotnikov. Focusing of electron
bunches wake fields in a plasma-dielectric wave-
guide // The Journal of Kharkov National Univer-
sity. Physical series «Nuclei, Particles, Fields».
2012, v. 54, № 2, p. 64-68.
12. R.R. Kniaziev, G.V. Sotnikov. Quasistatic field in-
fluence on bunches focusing by wakefields in the
plasma-dielectric waveguide // Proc. IPAC2013,
Shanghai, China, 12-17 May 2013.
13. V.A. Balakirev et al. Excitation of wake fields by
relativistic electron bunches in a plasma with radial
density variation // Plasma Physics Reports. 1997,
v. 23, № 4, p. 290-298.
14. The Argone wakefield accelerator facility.
http://www.hep.anl.gov/awa/index.html
15. J.P. Verboncoeur, A.B. Langdon, N.T. Gladd. An
object-oriented electromagnetic PIC code // Com-
puter Physics Communications. 1995, v. 87, № 1-2,
p. 199-211.
Article received 13.05.2013.
ФОКУСИРОВКА ЭЛЕКТРОННЫХ И ПОЗИТРОННЫХ СГУСТКОВ
В ПЛАЗМЕННО-ДИЭЛЕКТРИЧЕСКОМ КИЛЬВАТЕРНОМ УСКОРИТЕЛЕ
Р.Р. Князев, О.В. Мануйленко, П.И. Марков, Т.К. Маршалл, И.Н. Онищенко, Г.В. Сотников
Представлены результаты аналитических и численных исследований возбуждения кильватерных полей в
плазменно-диэлектрической структуре. В линейном приближении (сверхплотная плазма) показано, что при
определенной плотности плазмы суперпозиция плазменной и диэлектрической волн позволяет ускорять тес-
товый сгусток с его одновременной фокусировкой. Также мы выполнили моделирование методом “частица
в ячейке” возбуждения кильватерных полей для случаев разреженной и сверхплотной плазмы.
ФОКУСУВАННЯ ЕЛЕКТРОННИХ І ПОЗИТРОННИХ ЗГУСТКІВ
В ПЛАЗМОВО-ДІЕЛЕКТРИЧНОМУ КІЛЬВАТЕРНОМУ ПРИСКОРЮВАЧІ
Р.Р. Князєв, О.В. Мануйленко, П.І. Марков, Т.К. Маршалл, І.М. Онiщенко, Г.В. Сотнiков
Представлені результати аналітичних досліджень і чисельного моделювання збудження кільватерних по-
лів у плазмово-діелектричної структурі. У лінійному наближенні (надщільна плазма) показано, що при пев-
ній щільності плазми суперпозиція плазмової і діелектричної хвиль дозволяє прискорювати тестовий згус-
ток з його одночасної фокусуванням. Також ми виконали моделювання методом "частинка в комірці" збу-
дження кільватерних полів для випадків розрідженій і надщільної плазми.
|