Ion beam forming dynamics in an injector taking into account a plasma boundary
The developed program code taking into account a shape and position of an emitting plasma boundary had been used to examine a deuterium beam forming in an injector. It has been shown that the beams formed by a three-electrode ion optical system have a complex density distribution in a beam phase spa...
Gespeichert in:
Datum: | 2001 |
---|---|
Hauptverfasser: | , |
Format: | Artikel |
Sprache: | English |
Veröffentlicht: |
Національний науковий центр «Харківський фізико-технічний інститут» НАН України
2001
|
Schriftenreihe: | Вопросы атомной науки и техники |
Online Zugang: | http://dspace.nbuv.gov.ua/handle/123456789/79260 |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Zitieren: | Ion beam forming dynamics in an injector taking into account a plasma boundary / P.O. Demchenko, M.G. Shulika // Вопросы атомной науки и техники. — 2001. — № 3. — С. 144-146. — Бібліогр.: 7 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-79260 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-792602015-03-31T03:02:13Z Ion beam forming dynamics in an injector taking into account a plasma boundary Demchenko, P.O. Shulika, M.G. The developed program code taking into account a shape and position of an emitting plasma boundary had been used to examine a deuterium beam forming in an injector. It has been shown that the beams formed by a three-electrode ion optical system have a complex density distribution in a beam phase space. For the narrow current range it is possible to have a beam crossover at the injector output with a radius of 1mm and divergence less then 20 mrad. 2001 Article Ion beam forming dynamics in an injector taking into account a plasma boundary / P.O. Demchenko, M.G. Shulika // Вопросы атомной науки и техники. — 2001. — № 3. — С. 144-146. — Бібліогр.: 7 назв. — англ. 1562-6016 PACS numbers: 29.20.Bd http://dspace.nbuv.gov.ua/handle/123456789/79260 en Вопросы атомной науки и техники Національний науковий центр «Харківський фізико-технічний інститут» НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
description |
The developed program code taking into account a shape and position of an emitting plasma boundary had been used to examine a deuterium beam forming in an injector. It has been shown that the beams formed by a three-electrode ion optical system have a complex density distribution in a beam phase space. For the narrow current range it is possible to have a beam crossover at the injector output with a radius of 1mm and divergence less then 20 mrad. |
format |
Article |
author |
Demchenko, P.O. Shulika, M.G. |
spellingShingle |
Demchenko, P.O. Shulika, M.G. Ion beam forming dynamics in an injector taking into account a plasma boundary Вопросы атомной науки и техники |
author_facet |
Demchenko, P.O. Shulika, M.G. |
author_sort |
Demchenko, P.O. |
title |
Ion beam forming dynamics in an injector taking into account a plasma boundary |
title_short |
Ion beam forming dynamics in an injector taking into account a plasma boundary |
title_full |
Ion beam forming dynamics in an injector taking into account a plasma boundary |
title_fullStr |
Ion beam forming dynamics in an injector taking into account a plasma boundary |
title_full_unstemmed |
Ion beam forming dynamics in an injector taking into account a plasma boundary |
title_sort |
ion beam forming dynamics in an injector taking into account a plasma boundary |
publisher |
Національний науковий центр «Харківський фізико-технічний інститут» НАН України |
publishDate |
2001 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/79260 |
citation_txt |
Ion beam forming dynamics in an injector taking into account a plasma boundary / P.O. Demchenko, M.G. Shulika // Вопросы атомной науки и техники. — 2001. — № 3. — С. 144-146. — Бібліогр.: 7 назв. — англ. |
series |
Вопросы атомной науки и техники |
work_keys_str_mv |
AT demchenkopo ionbeamformingdynamicsinaninjectortakingintoaccountaplasmaboundary AT shulikamg ionbeamformingdynamicsinaninjectortakingintoaccountaplasmaboundary |
first_indexed |
2025-07-06T03:18:06Z |
last_indexed |
2025-07-06T03:18:06Z |
_version_ |
1836865958471270400 |
fulltext |
ION BEAM FORMING DYNAMICS IN AN INJECTOR TAKING INTO
ACCOUNT A PLASMA BOUNDARY
P.O. Demchenko, M.G. Shulika
Institute of Plasma Electronics and New Methods of Acceleration,
National Scientific Center "Kharkov Institute of Physics and Technology"
Kharkov, Ukraine, demchenko@kipt.kharkov.ua
The developed program code taking into account a shape and position of an emitting plasma boundary had been
used to examine a deuterium beam forming in an injector. It has been shown that the beams formed by a three-elec-
trode ion optical system have a complex density distribution in a beam phase space. For the narrow current range it
is possible to have a beam crossover at the injector output with a radius of 1mm and divergence less then 20 mrad.
PACS numbers: 29.20.Bd
1 INTRODUCTION
Now in the world some proton and deuteron linacs
up to energy 1-1.5 GeV and an average beam current of
5÷100 mA are being developed [1]. These accelerators
are intended for production of high neutron flows or for
tritium production. For realization of these accelerators
the high current injectors are needed with a high beam
luminosity. In all designs it is supposed the injection of
ion beams with the energy 50÷100 keV in an initial part
of the accelerator on the basis of a radio-frequency
quadrupole focusing (RFQ). That imposes some re-
quirements on a value and shape of the beam phase vol-
ume at an RFQ input.
The requirements to the ion optical performances of
the injected beams have stimulated the mathematical
modeling of ion beam forming from the plasma of gas
discharge sources. One of the modeling problems is the
definition of the shape and location of a self-consistent
plasma boundary, from which ions are extracted, and
initial conditions for starting particles. In the present
work for examination of ion beam forming dynamics
from plasma with the help of three-electrode ion optical
system the code INJECTOR was used [2]. The code
takes into account a location and shape of an emitting
plasma boundary. The results of the numerical modeling
are given for a deuterium beam with the output energy
100 keV. The beam is intended for injection to a
deuteron linac with the output energy 10÷12 MeV for
production of a medicine isotope-generator 99Mo as a re-
sult of a nuclear reaction 98Mo(d,p)99Mo [3].
2 THE PROGRAM CODE FOR SIMULA-
TION OF ION BEAM FORMING
The description of the INJECTOR code is given in
[2]. The code is based on the solution of a motion prob-
lem for ions, emitted from a plasma boundary, in a self-
consistent electric field produced by forming system
electrodes with given potentials and a beam space
charge. The emitting plasma boundary is supposed
smooth and crossing an aperture in a plasma electrode 1,
Fig. 1, which is contacting with plasma.
The shape of the boundary was supposed to be a sur-
face of revolution of the second degree curve (spherical,
elliptic or parabolic segment). In zero approach a plas-
ma boundary position (a zb point of crossing with a lon-
gitudinal z-axis) was defined from the balance of a ki-
netic plasma momentum and momentum of a vacuum
electric field produced by injector electrodes: nkTe=ε
oE2/2, where n is the plasma density, Te is the electron
temperature, E is the electric field strength, k is the
Boltzmann constant, ε0 is the dielectric constant of vac-
uum.
Fig. 1. Deuteron trajectories for the beam current
I=40mA.
Since the initial parameters of the problem are sup-
posed: a beam current I, electron Te and ion Ti tempera-
tures, the plasma density was determined from the ex-
pression for an ion saturation current: I=0.4Zen
(2kTe/M) 1/2S [4], where Ze and M are the charge and the
mass of an ion, respectively, S is the area of an emitting
plasma surface.
The plasma boundary potential was guessed equal to
the plasma electrode potential. The potential jump in a
double layer was taken into account in the start ion ve-
locity. I.e. it was supposed, that ions are starting along a
normal line to the plasma surface with a velocity
vo=(kTe/M)1/2, defined by the Bohm criterion [4].
The ion velocity component, tangential to plasma
surface, was defined by the ion temperature Ti and was
(kTi/M)1/2. Knowledge of the plasma boundary and the
initial velocities of ions allows to calculate their trajec-
tory in an injector ion optical system and respectively a
space charge distribution ρ(r,z). For the calculation of
ВОПРОСЫ АТОМНОЙ НАУКИ И ТЕХНИКИ. 2001. №3.
Серия: Ядерно-физические исследования (38), с. 144-146.
144
mailto:demchenko@kipt.kharkov.ua
the space charge the method of macroparticles was used
[5].
Fig. 2. Deuteron trajectories for the beam current
I=90 mA.
For successive iterations of the electric field deter-
mination the Poisson equation was solved using a new
space charge distribution. Every time the position of the
plasma boundary was defined more exactly using the re-
quirement that the electric field strength would be ap-
proximately zero in a point zb where the boundary was
crossed by z-axis. During the simulation the minimum
boundary field strength was taken kTe/ehz, where hz is a
step along z-direction of a computing grid, e is elemen-
tary charge.
This process had been repeated until the converging
solution had been succeeded.
3 RESULTS OF DEUTERON BEAM FORM-
ING SIMULATION
The INJECTOR program code was used to simulate
the forming process for a deuterium beam with the ener-
gy 100 keV. The ion optical system with an accel-deccel
potential difference between electrodes was considered.
The respective ion injector is supposed to use for beam
injection in RFQ of a deuteron linac of energy 10÷
12MeV for production of a medicine isotope-generator
99Mo. This radionuclide may be produced due to a nu-
clear reaction 98Mo(d,p)99Mo under bombardment by
deuterons of a target from a natural molybdenum or one
enriched by a natural isotope 98Mo [3].
The electron Te and ion Ti plasma temperatures were
chosen Te=5eV and Ti=1eV, respectively, as the most
probable for many gas discharge plasma sources [6].
The geometry of the injector forming system is shown
in Fig. 1.
The plasma electrode 1, Fig. 1, directly contacting
with the gas discharge plasma, had potential U1=100 kV
relative to a grounded electrode 3. Radius of the input
aperture in the plasma electrode 1 was ro=0.4 cm. A po-
tential of an extraction electrode 2, Fig. 1, was
U2 = -10 kV. A length of an extraction gap between
electrodes 1 and 2, Fig. 1, was d=2 cm. i. e. the length d
was above the minimum value causing a vacuum break
down: d≥1.4⋅10-3U3/2, where d is in cm, U=U1-U2 is in
kV [6].
In Fig. 1 and Fig. 2, the trajectories of deuterons are
shown for two values of a total beam current I or ac-
cordingly for different plasma densities n. In Fig. 3 and
Fig. 4 the respective distributions of particle density in a
transverse phase plane of the formed beams are given.
Fig. 3. Density distribution in a transverse plane of
the beam phase space for I=40 mA.
Here in Fig. 3 and Fig. 4 the root mean square values
(rms) of a beam width Xrms and its divergence X'rms at
an injector output are indicated and also scales Xm and
X'm of the respective coordinate axes.
Fig. 4. Density distribution in a transverse plane of
the beam phase space for I=90 mA.
As well as in case of a proton beam in ref. [2], for a
deuteron beam there are three modes of its forming, de-
pending on a total beam current I or according to plasma
density n. At small currents I with regard to the Child-
Langmuir value Ich=(4/9) (Ze/M)1/2εoU3/2S/d2, the diver-
gent beams are formed with a crossover in the extrac-
tion gap and a concave emitting boundary of plasma,
Fig. 5, I=10 mA.
With increase of the beam current the crossover dis-
places to the injector output. At some current value
I=Iopt, in this case Iopt≈40÷50 mA, at the injector output,
Fig. 1, the minimum beam width, is produced, Fig. 3.
The beam radius does not exceed 1mm. Thus the nor-
malized emittance is εn=2⋅10-2 π⋅cm⋅mrad. Simultane-
ously the emitting plasma boundary approaches to flat
in limits of the plasma electrode aperture, Fig. 5,
145
I=40 mA.
At I>Iopt the plasma boundary becomes convex to
the extraction gap, Fig. 5, I=90 mA, and the formed
beams are strongly divergent and have no crossover,
Fig. 2 and Fig. 4.
Fig. 5. Shape of the plasma boundary versus a beam
current.
Fig. 6. Spherical (1) and elliptical (2) shapes of the
plasma boundary for the beam current I=90 mA.
It is necessary to note that with the help of the
smooth second-order curves it is impossible to obtain
the solution with the constant electric field strength
along the plasma boundary.
The results of modeling have given an oscillating
character of the field strength. At currents of a formed
beam of about I≈Ich the minimum field oscillations and
also the higher convergence rate of the solution are ob-
served when the plasma boundary has an ellipsoid form
with the relation of a vertical semi-axis to the longitudi-
nal one of order 4, curve 2, Fig. 6.
4 CONCLUSION
The results of the mathematical simulation indicate
that the chosen ion optical system allows to produce a
deuterium beam with a radius no more than 1 mm and
divergence less then 20 mrad at an injector output, but
in a very narrow interval of a beam current of 40÷
50 mA at a deuteron energy 100 keV. It is obviously pos-
sible to calculate a RFQ section, which will transport
such a beam without a matching unit. At other currents
of a beam or its ion optical performances, a matching
unit, for example similar LEBT (Low Energy Beam
Transport) used in [7], is needed. The nonlinear effects
caused by a vacuum electric field, and a space charge
field of a high current beam, are resulting in a complex
particle density distribution in a beam phase space and
are increasing a beam effective emittance.
REFERENCES
1. I.D.Schneider. Overview of High Power CW Pro-
ton Accelerators // Proceedings of EPAC 2000, Vi-
enna, Austria, 2000, p. 118-121.
2. P.O.Demchenko, M.G.Shulika. A Simulation Code
to Model an Ion Beam Forming in an Injector //
Proc. of the Seventh European Particle Accelerator
Conference, Vienna, Austria, 2000, p. 1363-1365.
3. P.O.Demchenko, V.A.Voronko, V.Ya.Migalenja
and et al. Application of compact linacs for the
medical purposes // Problems of Atomic Science
and Technology, Issue: Nuclear Physics Research-
es (31-32). 1997, v. 4-5, p. 168-170 (in Russian).
4. A.T.Forrester. Large Ion Beams. New York: John
Wiley and Sons, 1988.
5. D.Potter. Computational Physics. New York: John
Wiley and Sons, 1973.
6. G.Bolme, L.Hansborough, T.Hardek, et al. Half-
Power Test of a CW Proton Injector with a
1.25 MeV RFQ // Proceedings of the 1999 Particle
Accelerator Conference, New York, 1999,
p. 1946-1948.
7. H.V.Smith, T.Figueroa, L.D.Hansborough, et al.
Comparison of Beam Simulations with Measure-
ments for the LEDA LEBT H + Beam // Proceed-
ings of the 1999 Particle Accelerator Conference,
New York, 1999, p. 1929-1931.
146
|