Open access peer-reviewed chapter

Concept and Numerical Simulations of Multi-Beam Linear Accelerators EVT with Depressed Collector for Drive Beam

By Vladimir E. Teryaev

Submitted: May 23rd 2017Reviewed: September 13th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.70956

Downloaded: 367

Abstract

The concept of the electron multi-beam linear accelerator electron voltage transformer (EVT) invokes aspects of a two-beam accelerator. Here, effective transformation of energy of several drive beams to energy of a high voltage accelerated beam takes place. It combines an RF generator and essentially an accelerator within the same vacuum envelope. Acceleration occurs in inductively-tuned cavities. The ratios between values of RF beam currents and their voltages are similar to ratios for an electric transformer that is leaded to dub this device an electron voltage transformer (EVT). Properties and high efficiency of the EVT accelerator proves to be true by the example of numerical simulations of three accelerator configurations generating electron beams with energies 1, 2 and 4.3 MeV, respectively. The submitted results of numerical simulations show a resource for the further increase in efficiency of the EVT accelerator up to 70–75% by using a depressed collector for a spent dive beam.

Keywords

  • two-beam accelerator
  • electron gun
  • drive beam
  • accelerated beam
  • RF buncher
  • accelerating cavity
  • depressed collector

1. Introduction

Two-beam linear collider can serve as the most known example use of two-beam concept in the field of high energy physics [1]. Derbenev Ya et al. [2] can serve as an example of consideration of the two-beam concept in the range of low energy (1–10 MeV). Specificity of low energy when the relativistic factor is low is that a change of velocity of particles for conservation of the appropriate phase ratio between drive bunches and accelerated bunches, i.e. for keeping of synchronism in process of beam-beam interactions, shall be taken into account. Derbenev Ya et al. assumed that the use of the H 020 mode in a pill-box cavity where acceleration occurs on an axis of the cavity and the annular drive beam changes its radius during movement [2]. Synchronization is carried out due to change in a focusing magnetic field along the axis of the device and this implies changing of the radius of the annular drive beam. However, this concept has never been developed any further.

Here, a new concept of the electron accelerator briefly named electron voltage transformer (EVT) is considered, where the multi-beam approach is of a crucial importance. The approach to use several beam-lets interacting with the common cavity is well-known and is implemented successfully in multi-beam klystrons (MBK’s) [3, 4]. Such design is also inherent to the greatest degree in the concept of EVT. All beams in EVT are generated in a common electron gun. Then, they are bunched in RF cavities as in a klystron. The energy is transferred from the drive bunches to the accelerated bunch in the accelerating structure. Each cavity of the structure is coupled to all beams, but the cavities are not coupled one to another.

Basic difference from the concept quoted above [2] is that axes of all beams are parallel straight lines and that for beams focusing a homogeneous magnetic field of a solenoid is used. The cavities in EVT are inductively detuned, having the resonant frequencies of their working modes being higher than the operating frequency. As a result, the phase angle between the gap voltage of the cavities and the fundamental harmonic of the drive beam current is close to 90°, as against conditions for the accelerated beam where the phase angle between accelerating gap voltage and a harmonic of a current is close to zero. This operating regime allows one to achieve the desired value of the accelerating field. The most important advantage of using inductively tuned cavities is effective longitudinal bunch focusing, which is exerted by the cavity RF fields of the drive bunches. The drive beam-lets in EVT are weakly relativistic. Therefore, synchronous energy transfer from the drive beams to the accelerated beam is achieved by correct choices of distances ∆z between the adjacent cavity gaps. A similar condition, of two-beam acceleration in inductive tuned cavity, is considered by Kazakov et al. [5], where, in contrast to EVT, the authors have placed the drive beam and the accelerated beam on the same axis. Dolbilov G has discussed a two-beam accelerator of protons, where the use of inductive tuned cavities for an electron drive beam is considered [6, 7]. The accelerated beams and the drive beams move there in opposite directions in order to maintain the condition of synchronization.

Preliminary results of numerical simulations using the known and reliable code MAGIC [8] are presented for an EVT prototype, which operates in S-band, where acceleration up to 1 MeV was predicted at efficiency of 66% [9]. A project version of an EVT accelerator operating at L-band, in which a 6-A, 110-kV beam is accelerated to 2 MeV with an overall efficiency of 60%, is presented below. The further development of the EVT concept is based on a principle of two-stage acceleration, where an ensemble of the drive beam-lets is consisting of two groups. At the first stage, the first group of drive beam-lets transfers their energy to the second group of drive beam-lets. At the second stage, the second group of drive beam-lets transfers their energy to the accelerated beam. Such configuration allows a transformation ratio and energy of the accelerated beam to be increased considerably. Numerical simulations predict acceleration of an initial 1 A, 60 keV beam up to energy of 4.3 MeV at efficiency of 22%. The use of a novel view on a Depressed Collector (DC) for a spent drive beam, where the DC is considered as a separator of a spectrum of particles, increases efficiency in the first two EVT projects up to 74 and 68%, respectively. Numerical results are shown in Figures and Tables.

A positive feature of the EVT concept is that the accelerator does not require an external high power RF-source, but only a 60–110-kV pulsed source for operation. The duty cycle for the device will be easily adjusted via the pulse duty cycle of the power supply. The absence of an RF-source, utilization of a relatively low-voltage pulsed power supply to drive the machine, and the compact size of the device would ensure a relatively light weight and highly versatile applied accelerators.

2. EVT configuration

The schematic representation of an EVT accelerator is shown in Figure 1 , where the accelerating structure is pop up.

Figure 1.

Schematic representation of an EVT (left); shape of the accelerating cavity and electric field E z along x-axis (right).

An electron gun, a magnetic focusing system and a RF buncher will be provided below in the description of specific projects of EVTs. Figure 1 (right) provides an example of the accelerating cavity. The operating mode is TM 020. An accelerated beam is located in the centre, and a drive beam is located on the periphery, at a peak value of an electric field. The drive beam consists of six beam-lets in the case shown. An additional passive ring-shaped bunching cavity with the fundamental TM 010 mode is located in a long drift tube. Synchronous energy transfer from the drive beams to the accelerated beam can be achieved through correct choices of distances ∆z1 and ∆z2 between gaps of adjacent accelerating cavities.

3. Scientific and technical justification of the EVT concept

It is known that a field Ertin a volume Ωof a cavity generated by a beam with a current density Jrtcan be presented as eigenmode expansion: Ert=λeλtEλr, where eigen vectors Eλrare functions of coordinate rthat are normalized by energy of the cavity as follows:

Wλ=ε02ΩEλr2dΩ.E1

A time-dependent multiplier eλtis the solution of the equation of oscillator with the second member, Eq. (2):

d2eλdt2+ωλQλdeλdt+ωλ2eλ=12WλddtΩJ̇EλE2

Let’s express Eq. (2) in a complex form and consider one of the eigenmodes only, having frequency ω0, replacing: eλtė0expiωt, JrtJrexpiωt. It gives

ė0=iωω2ω02iωω0/Q012W0ΩJ̇rE0rE3

Other terms of expansion generated by the beam current will have small values at high quality Q0of the cavity and at ωωO, and they will not be taken into account. Eqs. (2) and (3) are the most general expressions for determination of fields in cavities [10]. As it takes place ωωO, the equation Eq. (3) can be rewritten as:

ė0=expQ01+ξ0212ωW0ΩJ̇E0E4

where ξ0is a relative detuning:

ξ0=Q0ω2ω02ωω02Q0ω0ωω0=ωω0φ=tan1ξ0E5

Practically, the cavity is described via integral characteristics, such as geometric impedance of drive gap (index 1) and accelerated gap (index 2):

ρ1=R/Q1,ρ2=R/Q2;E6

circuit current İ1of the drive gap and circuit current İ2of the accelerating gap equal to

İ1=exp1TG12FB1I01,İ2=exp2TG22FB2I02;E7

RF beam currents:

IRF1=2FB1I01,IRF2=2FB2I02;E8

voltage V̇1of the drive gap along axis z1and voltage V̇2of the accelerating gap along axis z2:

V̇1=ė0V1=ė0z1E0zdz,V̇2=ė0V2=ė0z2E0zdzE9

where: φ1and φ2are phase angles; TG1, TG2are transit time factors of the drive and accelerating gaps, respectively; FB1, FB2are form-factors of bunches; I01, I02are average drive beam currents of the drive and accelerated beams.

For relativistic beams, when it is possible to neglect relative moving of particles inside a bunch, and taking into account that the integral on the second member of Eq. (2) is not equal to zero in gap 1 and gap 2 only where beam currents exist, the Eq. (2) is rewritten as follows:

ė0=expQ01+ξ0212ωW0İ1V1+İ2V2.E10

Let us add in relationships under consideration

KV=V2/V1=ρ2/ρ1,E11

following from: 2ωW0=V12/ρ1=V22/ρ2.

The description of the equivalent circuit for a case of interaction of relativistic beams with the cavity is

V̇1=expQ0ρ11+ξ02İ1+KVİ2=ZCİE12
V̇2=KVV̇1E13

where ZCis impedance of the drive gap, φis a phase angle between gap voltage V̇1and equivalent circuit current İ. Eq. (13) shows properties of the chosen operating eigenmode of a cavity when phases of fields in accelerating and drive gaps are shifted by angle π.

At interaction of a cavity with low relativistic beam it is necessary to consider also beam-loading effect or a beam admittance YB=GB+iBB. Therefore, Eq. (12) is revised as

V̇1=11/ZC+YBİ1+KVİ2E14

Let’s set out well-known analytical formulas for a rectangular bunch with angular length ψBand for the flat cavity gap having a transit time angle ψC:

TC=sinψC/2ψC/2;FB=sinψB/2ψB/2;E15
GB=I0VB12TG21ψC2ctanψC2;BB=I0VBsinψC1+cosψCψC/2ψC2,E16

where I0is an average beam current, VBis an average beam voltage in the corresponding gap (drive or accelerating), TGis the corresponding transit time factor and FBis a form-factor of the corresponding beam.

3.1. Conditions for preset drive bunching maintenance

As it is known, a bunching cavity has inductive detuning, i.e. ω<ω0, at which the phase angle between a gap voltage and an equivalent current of a cavity circuit becomes close to π/2. Figure 2 shows mutual position of bunches, currents and gap voltages on a phase plane in an accelerating cavity. The phase position of a head Ḣheadand tail Ṫtailof a drive bunch is shown here, too. Let’s notice, that complex RF beam current İRF1is displaced at angle πrelative to the circuit current İ1.

Figure 2.

Electric field in cavity gaps (left) and position of complex values on the phase plane (right), providing phase focusing of a drive bunch and acceleration of the beam.

3.2. An impedance of the drive gap loaded with the accelerated beam

As it is evidently from Figure 2 , suitable conditions for an accelerating cavity shall be as follows: φπ/2, expi. It is the chosen phase of drive current: φ1=0. The phase of the accelerated beam current equal to φ2π/2provides a bunch being in the electric field close to the peak value, i.e. İ1=I1and İ2=iI2. Thus, using Eqs. (12) and (13), the required impedance of the loaded drive gap is

Z1=V̇1İ1=Q0ρ11+ξ02KVKTr+iE17

where KTr=I1/I2is a transformation ratio. An impedance of the accelerating gap is similarly found:

Z2=V̇2İ2=Q0ρ11+ξ02KVKTriKV2E18

On the other hand for the equivalent circuit of the cavity the impedance Z1is given as:

Z1=r+ix=Qloadρ11+ξload21iξloadE19

where Qloadand

ξload2Qloadωload/ωE20

are parameters of the equivalent loaded cavity circuit. Equating the real and imaginary parts of Eqs. (18) and (19), the loading quality Qloadand relative detuning ξloadcan be represented as follows:

Qload=Q01+ξ02KTrKV1+KV2KTr2=Q01+ξ02KTrKVξload=KTrKVE21

Considering Eq. (21) together with Eq. (20) leads to:

ωloadω=1+ξ022Q011+KV2/KTr2ω0ωE22

where a high value of detuning ω0is suggested, which results in 1+ξ02ξ0. In addition, it is assumed that the transformation ratio is large enough, i.e., KV2/KTr21. For example, KV2, KTr20, ω0/ω103, Q08000, ξ016.

Eq. (22) shows that loading of the accelerated beam practically does not change the imaginary part of the impedance of the cavity, and, thus, it is possible to model loading as a certain equivalent external load quality. This feature is used below in the description of numerical model, see Section 4.

3.3. Voltage transformation ratio

Let’s find the voltage VB1lost by the drive beam and the voltage VB2added by the accelerated beam, which are connected to a voltage V1of the drive gap and a voltage V2of an accelerating gap as assigned at the corresponding moment of time (see Figure 2 ):

VB1TG1V1,VB2TG2V2E23

where transit-time factors TG1, TG2are determined the same as in Eq. (7). Let us find voltages V1and V2using Eqs. (12) and (13) as follows:

V1=ReZ1I1=Q0ρ1ρ21+ξ02I2V2=ReZ2I2=Q0ρ1ρ21+ξ02I1E24

where the transformation ratio KTris:

KTr=I1I2=V2V1E25

In view of Eqs. (7), (8), and (23), the transformation ratio can be written through the RF beam currents as:

KTr=IRF1IRF2VB2VB1,E26

or, at a high value of transit time factors TG1, Eq. (26) can be written as

KTr=I01I02,E27

where I0is the average current of the corresponding beam. It reflects the law of conservation of energy for the idealized conditions applied in the beginning of section 3.3 and corresponding to small resistive loss in a cavity. In view of losses it is resulted:

KTr=I1I2=V2V11ηBE28

where ηBis efficiency of power transfer, defined in the Section 3.4.

3.4. Efficiency of power transfer between beams and other parameters of the EVT

From conservation of energy, the efficiency of power transfer from one beam to another can be estimated as:

ηB=V2I2V2I2+V22/2Rs2E29

where Rs2=Q0ρ2is the shunt impedance of the accelerating gap. For achieving a high efficiency of energy transfer, it is important to have a large enough current I2of the accelerated beam, i.e. I2V2/2Rs2.

Based on the above consideration and Figure 2 , it can see that in order to reach the effective bunching of drive beam, with an angle between the current İ1and the voltage V̇1close to π/2, two conditions shall be satisfied. The first is to have a large enough value of the cavity detuning ω0; the second is I1/KVI2. The angle between İ1and İ2should be close to π/2for effective acceleration since the accelerated bunch should pass through an accelerating gap at the time when the voltage on the gap is near its maximum. Thus, conditions

I1/KVI2V2/2RsE30

are the necessary factors of effective performance of the EVT, which are readily achieved when multiple drive beam-lets are used. For example, I1=60A, KV=2, I2=3A, V2=80kV, V2=100kV, Q08000, ρ2=60ohm, ηB=0.95.

We define the efficiency ηof the accelerator as a ratio of the power of the accelerated beam at the EVT output to the total beam power in the gun, namely:

η=PAccelIGunVGunE31

where IGunis the total current of the gun, and VGunis the gun voltage. The high EVT efficiency is confirmed by the results of numerical simulations, as shown in Table 1 .

Multi-beam linear acceleratorS-bandL-band
Operating frequency, F28561300MHz
Voltage of common electron gun60110kV
Output energy of accelerated beam1.12MeV
Total gun current of drive beam72180A
Gun current of accelerated beam36A
Focusing magnetic field in solenoid8001000G
Energy spread of accelerated beam±10±10%
Micro-perveance of one drive beam-let0.8160.82μA/V3/2
Transmission of drive beam0.990.99
Transmission of accelerated beam0.900.87
Transformation ratio2430
Voltage gain factor≈18≈18
Typical R/Q factor of accelerating gap5854Ω
Typical cavity detuning, ∆F = F0 - F8 ÷ 127MHz
Typical KV factor of gaps of cavity (Eq. 11)4.55
Efficiency of transfer of power (Eq. 29)9497%
Number of beam-lets of drive beam66
Number of cavities of RF bunchers44
Number of accelerating cavities1412
Number of passive bunching cavities17
Length of RF part of accelerator0.792.4m
Efficiency of accelerator EVT≈66≈60%
Efficiency of EVT with DC for drive beam7266%

Table 1.

The simulated parameters of the S-band and L-band linear accelerators EVT.

We define a voltage gain factor to be a ratio of the voltage of the accelerated beam at the EVT output (the final beam energy) to the gun voltage. It is:

KEVT=VAccel/VGun1+KTrηE32

where ηand KEVTare shown in the Table 1 , as a result of numerical simulations.

3.5. Condition of synchronization and phasing of bunches

The drive beam-lets in the EVT are essentially low relativistic. Therefore, a synchronous energy transfer from the drive beams to the accelerated beam is achieved through correct choices of distances ∆z between the adjacent accelerating cavities gaps. In case of equality of distances between gaps shown in the Figure 1 , a condition of synchronism is as follows:

z=z1=z2=2πcωβ1β2β1β2E33

where β1cand β2care velocities of the drive and accelerated bunches in space zbetween cavities gaps. Phasing capabilities considerably extend if distances z1and z2between a gaps of adjacent cavities are not equal to each other. More general Eq. (34) is corresponding to this case:

z1β1z2β2=Nλ;N=0;1;2;;λ=2πc/ωE34

It is supposed to use RF bunchers of klystron type for generation of drive bunches. Further bunches should be decelerated in cavities, gradually giving the energy to the accelerated bunches. There is a reasonable question for how long the drive bunches can keep their form and phase length. Results of 2D simulations show long enough conservation of the generated bunch (see below projects of EVT). The similar research confirming these results can be found in [6, 7].

4. Representation of 3D objects via equivalent 2D objects

3D codes have been used for simulation of elements of the EVT (cavities, magnetic system), but for the analysis and optimization of dynamics of ensemble of the interacting beams moving in numerous drift tubes the direct and slow 3D analysis is too intricate problem. Difficulty is aggravated also with necessity of carrying out numerous optimizations.

The drive beam-lets are located on a circle, see Figure 1 (right), and it can be quite replaced with single equivalent 2D annular beam. It is possible to choose the corresponding width of the ring filled with a beam, providing close values of plasma wave-lengths of 2D and 3D models of the beam. In the same way it is replaced real 3D cavities with equivalent 2D cavities. There are chose parameters (gap, R/Q) of interaction of a 2D beam with equivalent 2D cavities similar to interaction of a 3D beam with 3D cavities. Certainly, such equivalent model does not describe exact dynamics of the drive beam, but reflects quite well the properties of the EVT accelerator and can serve as initial approach. The results of parameter optimization for the EVT accelerators using the 2D code MAGIC [8], are provided in Sections 5–7.

It is obvious that, on the one hand, implementation of a pure annular concept of the drive beam is embarrassing technically, and, on the other hand, the annular model is rather rough to describe dynamics of a set of cylindrical drive beam-lets. This simulation is provided here to illustrate a method of acceleration used in the EVT as it correctly represents interaction of two electron beams without any additional assumptions.

The following, more adequate model is appropriate for description of a 3D problem, which can be named 2 × 2D model. This model is based on thesis about an impedance of the loaded drive gap considered in Section 3.3 and on the following circumstances:

  1. Each of beam-lets and electric RF fields of cavities near to a beam-let can be considered as axially symmetric relative to a local axis of the given beam-let;

  2. Parameters of drive beam-lets are identical;

  3. Accelerating cavities are weakly loaded, and, therefore, RF fields of cavities are rather close to a field of an operating RF eigenmode. I.e. the loaded quality of a cavity remains high enough. In other words, the RF power circulating in a cavity PRF=V122R/Q1is much higher than power P2=Poutput2Pinput2transferred to the accelerated beam, for example, Qload=PRF/P2100.

According to these points, it is considered a 3D problem as two connected 2D problems, see Figure 3 . The drive beam goes through the first equivalent 2D cavity, having the R/Q of drive gap equal to R/Q of corresponding 3D cavity. This cavity has an RF load in which RF power Pload=Pdriveis absorbed.

Figure 3.

Two 2D cavities, the parameters of which are related via the shown relationships, model processes in a 3D cavity.

The accelerated beam goes through the second cavity having a gap voltage Vgap2and a phase angle φgap2, which are calculated based on a voltage Vgap1on the gap of the cavity 1 for the drive beam and where power of the accelerated beam increases at the value P2=Poutput2Pinput2. Positions of gaps of these two equivalent cavities correspond to the position of gaps of the 3D cavity. Combination of these two processes will consist in maintenance of equations:

P2=Paccel=6PdrivePlossVgap2=KVVgap1φgap2=φgap1+πE35

where N=6is the number of drive beam-lets, Plossis resistive loss in the 3D cavity. Eq. (35) is carried out by choosing the corresponding conductivity of an external RF load of the first (drive) cavity; the value KV=Vgap2Vgap13Dis defined (see Eq. 11) upon the results of the 3D cavity eigenmode analysis. Numerical simulations using the code MAGIC show fast convergence of the Eq. (35) process. 2–3 iterations are practically enough.

5. Design of S-band 1 MeV linear accelerator EVT

The design of a linear accelerator EVT operating in S-band (F = 2856 MHz) and generating 1 MeV electron beam is presented in this section. One accelerated beam and 6 drive beam-lets are generated by a common electron gun fed from a 60 kV power supply. The magnetic lenses around the gun and two matching lenses ensure formation of beams without ripples in the same way, as it is described in Ref. [3]. Bunches are formed in the bunching section energized by two coherent microwaves via RF couplers. The generated drive beam-lets interact with the accelerated beam in the cavities of the energy transfer section (accelerating structure), with the cavities working at the TM020 mode. For a layout of the S-band EVT, please see Figure 4 . The results of the initial numerical analysis using code MAGIC are given below in Table 1 . Two models described above in Section 4 have been applied for the analysis. The models have shown close results. Beam dynamics of the annular model is more evident and results in Figure 5 . Figure 6 shows energy distribution of particles of beams in the 1st accelerating cavity (z > 314 mm) and at the output of the accelerating structure (z > 787 mm). Energy of 1.1 MeV, averaged on the size of a bunch, is achieved. You can see disruption of the spent drive bunches at the accelerator output. The achieved 66% efficiency is twice the best efficiency accessible for the facility consisting of a high power klystron and a high efficient linear accelerator [11].

Figure 4.

Lay-out of the accelerator EVT.

Figure 5.

The geometry of a compound RF buncher, accelerating structure and beams dynamics are shown for the S-band EVT.

Figure 6.

Particle energy distribution in the 1st accelerating cavity (left) and at the output of the accelerating structure (right).

5.1. Features of the multi-beam electron gun

The radius of a ring on which drive beam-lets are located is defined by a frequency and configuration of the cavities (see Figure 1 (right)). Six cathodes with a current of 12 A are freely located on this radius at the gun voltage 60 kV, Figure 7 . Clearly, a drive beam-let perveance should not be too high to not interfere with phase stability and synchronization of bunches.

Figure 7.

The multi-beam cathode of the electron gun. The peripheral cathodes generate the drive beam-lets, and the central cathode generates the accelerated beam.

Distance between the drive beam-let cathodes is big enough in comparison with distance between the cathode and anode therefore transverse electrostatic fields on an axis of drive beam-lets are insignificant. Thus the electron optics of each of beam-lets can be quite considered as axial symmetric. Figure 8 shows the 2D geometry and an example of simulation of a drive beam-let. The double lens around the gun allows obtain the desired profile of the magnetic field to focus the beam-lets and minimize their ripples.

Figure 8.

The geometry and an example of 2D simulation of the drive beam-let gun.

Cells of magnetic lenses are formed by flat iron pole-pieces and system of apertures (see Figure 4 ). Properties of cells for drive beam-lets are identical. The sizes of apertures and distances between them are such that there are no appreciable transverse fields generated on an axis of beam-lets. Thus, the magnetic field quite possesses axial symmetry near the beam-let axis.

6. Specification of L-band 2 MeV linear accelerator EVT

The S-band EVT project has shown achievement of a high efficiency, but has simultaneously shown limited capabilities for an acceleration rate and for a factor of transformation as decay of bunches after Ncav=15cavities has been observed. Numerical results show difficulty of bunch confinement at a voltage of a drive gap Vgap1>12Vgun. It results in the following limiting estimated condition:

VAccel<12KVNcavVgun,E36

where KVis from Eq. (11). This estimation is close to observable result VAccel=1.1MeVachieved in the S-band EVT.

Thus, an advance in a path of higher voltage of the gun and a low frequency is meant to increase a voltage of an accelerated beam in the described L-band EVT. Its layout does not differ from Figure 4 . Numerical simulations have shown feasibility of reaching the energy of the accelerated beam about 2 MeV with efficiency about 60% at a gun voltage equal to 110 kV. More detailed data of the project are submitted in Table 1 . The accelerating structure includes 12 accelerating cavities and 7 passive cavities (see Figure 9 ). The code MAGIC 2D and both models described in Section 4 have been applied for the analysis. For the results of beam dynamics of the annular model, please see Figure 9 , too. Figure 10 shows particle energy distribution.

Figure 9.

Geometry of the RF structure and the result of numerical simulations are shown for L-band EVT. For more details on the dynamics of the beams, please see below.

Figure 10.

Particle energy distribution at the output of the 1st accelerating cavity (left) and at the output of the accelerator section (right).

7. Two-stage multi-beam linear accelerator EVT

The further increase in energy of the accelerated beam will consist in increase in a voltage of a drive beam. It is suggested to consider two-stage acceleration. Let’s assume that the set of the drive beam-lets consists of two groups. At the first stage, the first group of drive beam-lets that form drive beam 1 transfer their energy to the second group of beam-lets that form drive beam 2. At the second stage, the second group of drive beam-lets transfers their energy to the accelerated beam. Such configuration is allowed to increase a transformation ratio and energy of the accelerated beam considerably. Preliminary results of numerical simulations of two-stage EVT operating in S-band with a 60 kV gun and generating a 1 A, 4.3 MV beam at its output, with an efficiency of 22%, are presented. See Figure 12 and Table 2 .

Operating frequency2856MHz
Voltage of common electron gun60kV
Output energy of accelerated bunch4.3MeV
Average energy of accelerated beam3.6MeV
Gun current of accelerated beam1A
Total gun current of drive beams252A
Solenoid magnetic field1000G
Length of RF part of accelerator1.9m
Transmission of accelerated beam90%
Total efficiency of two-stage EVT≈22%
Stage1st2nd
Gun current of drive beam21636A
Perveance of drive beam-let0.410.41μA/V3/2
Transformation ratio636
Voltage gain factor of stage≈3.5≈21
Number of drive beam-lets366
Transmission of drive beam9890%
Efficiency of stage≈55≈45%

Table 2.

Simulated parameters of the two-stage EVT.

Figure 11 (left) shows a multi-beam cathode system. 36 peripheral cathodes generate the drive 1 beam-lets, and 6 cathodes generate the drive 2 beam-lets. The central cathode generates the accelerated beam. A rather low perveance (0.41 μA/V3/2) of drive beam-lets is chosen. It is allowed to simplify the magnetic system and electron optics of the gun. The drive beam is frozen in a homogeneous magnetic field of the focusing solenoid. You can see the presence of ripples having an allowable value in Figure 11 (right).

Figure 11.

Multi-beam cathode system of the two-stage EVT (left) and an example of 2D simulation of the drive beam-let gun (right).

MAGIC 2D simulations of beam dynamics were carried out using an annular model of beams and cavities. Please see the results of numerical simulations in Figures 12 14 and in Table 2 .

Figure 12.

MAGIC 2D results of simulation of beam dynamics in the bunching part and at the 1st stage.

Figure 13.

Accelerating structure of 2nd stage. Dynamics of drive 2 and accelerated beams.

Figure 14.

Energy distribution at the output of the two-stage EVT.

Figure 14 shows particle energy distribution at the accelerator output. For the change in beams energy vs. distance, please see Figure 15 .

Figure 15.

Change of beams voltage along the accelerator axis.

8. Depressed collector (DC) for spent drive beam

The efficiency of 66% achieved in preliminary numerical simulations of the S-band EVT is a rather high value. The further increase in efficiency of the EVT accelerators is possible at the use of a depressed collector. It is suggested considering a novel concept of a DC, which allows theoretical estimation of DC properties. The given concept is based on the thesis of existence of an extremum of the DC efficiency versus DC voltage. The concept does not struggle with the energy distribution of electrons having a place after the device output, but uses properties of energy distribution. Let us name this concept a depressed collector with a spectra separator, DC&SS. The schematic circuit of the tube with a one-stage DC is shown in Figure 16 (left). Table 3 shows the basic relationships describing the circuit of the EVT with a DC.

Figure 16.

Schematic circuit of a tube with depressed collector. One-stage DC (left); two-stage DC (right).

Beam voltageVB
Voltage of stageVC , V1 , V 2
Current of potential part of stageIC , I1 , I2
Beam powerPB = IBVB
RF power, transfer to accelerated beamPRF
Loss power in ground part of collectorPG
Loss power in potential partPC , PC1 , PC2 ,
EVT electron efficiencyη e = PRF/PB
Beam power at output of tubePOut = (1- η e)PB
Recovery powerPR = VC IC
Total power generating by high-voltage sourcePt = VB IB - VC IC
Efficiency of depressed collectorη c = PR/POut = VC IC /POut
Total tube efficiency with DCη t = PRF /Pt  = PRF /(PB - VC IC )

Table 3.

Notations and relationships describing a DC.

Substituting the above resulted relationships, it is derived the total tube efficiency with a DC:

ηt=ηe11ηeηcE37

Eq. (36) is well known (see, for example [12]), but we shall regard it now jointly with the presence of an energy spectrum of particles leaving a tube. Variable Ic=IcVcis remarkable since it is a function of a collector voltage as a consequence of the presence of energy distribution. Our subsequent problem is to study the function ηcVcin order to find an extremum at some assumptions concerning the energy distribution jV, which is normalized on the beam current IB:

0jVdV=IBE38

In the gap of the separator, the beam is divided into low-voltage and high-voltage fractions. The current

IcVc=VcjVdVE39

that has reached the potential part of the collector is the fraction of the current at which V>VC. In the same way, it is written down the equation for the beam power leaving the accelerating structure:

Pout=0VjVdVE40

Thus, it is evaluated the fundamental formula Eq. (41) for efficiency ηcof a depressed collector included in the DC&SS concept:

ηcVc=VcVcjVdV0VjVdVE41

Figure 6 (right) above shows a numerical example of the energy distribution for the spent drive beam of the S-band EVT versus distance along beam axes. Figure 17 shows the energy distribution of the same beam. The gamma distribution xex(curved line) can serve as a close approximation for numerical distribution.

Figure 17.

A numerical simulated histogram of particle distribution vs. voltage and an approximating gamma distribution are shown for the S-band EVT.

Thus, let us approximate jVby the gamma distribution:

jV=IBV0xexE42

where x=V/V0, V0=VB1ηe/2. Further, it is defined Icaccording to Eq. (39):

Icxc=IBexCxC+1, where xc=Vc/V0, and it is defined Poutaccording to Eq. (40):

Pout=2IBV0E43

Substituting the found expressions in Eq. (41), it is obtained:

ηcxc=12xcexcxc+1E44

The function ηcxchas the maximum equal to 0.42 at xc=1+52=1.618.

Thus, it is derived the collector efficiency and other parameters at the optimum point:

ηc=0.42;Ic=0,519IB;IG=0,481IB;Vc=1.618V0.E45

As an example of the S-band EVT ( Table 1 ), one have:

VB=6104;ηe=0.66;V0=VB1ηe/2=1.02104;Vc=1.65104.E46

The total efficiency (according to Eq. 37) at the optimum point is ηt=0.77, which is the upper theoretical limit!

For a DC having m stages i.e. m gaps and consequently m gap voltages (see Figure 16 ), efficiency is as follows:

ηt=ηe11ηeηcηc=n=1mηnE47

where ηcis DC efficiency in multi-stage case. The collector efficiency ηnof the nth stage is calculated according to the Eq. (41), where the lower limit of integration is equal to potential Un=i=nmViof the given stage concerning the ground. This Eq. (48) defines the concept of a DC&SS in multi-stage case:

ηnV1V2Vn=VnUnjVdV0VjVdVE48

Using the approximation (see Eq. 42) for beam distribution jV, leads to:

ηn=12xn1+i=nmxiei=nmxiηc=n=1mηnE49

where: xi=Vi/V0, xn=Vn/V0, V0=VB1ηe/2, Viis gap voltage of the i th stage, VBis beam voltage. For a schematic circuit of a tube with a two-stage DC, Figure 17 (right), the total collector efficiency will be as follows:

ηcx1x2=η1+η2=12x11+x1ex1+x21+x1+x2ex1+x2E50

The extremum of this function is equal to ηc=0.586at: x1=1.118; x2=1.398.

At ηe=0.66, according to Eq. (37), it will correspond to the total tube efficiency ηt=0.824.

Table 4 shows the results of efficiency ηcoptimization for a collector that contains one to six stages.

mx1x2x3x4x5x6ηc
11.6180.420
21.1181.3980.586
30.8840.9411.3190.676
40.7430.7320.8761.2760.734
50.6480.6080.6750.8411.2490.773
60.5780.5260.5580.6450.8191.2290.803

Table 4.

The optimization result for the function ηcx1x2via its parameters x1,x2,.

The theoretical premise of an increase in efficiency of a device with a multi-stage DC are rather optimistic. However, simulations of beam dynamics in conditions close to real do not show any noticeable increase in the device’s efficiency while engineering complexity of such multi-stage DC tangibly increases. Thus, use of a one-stage DC is, perhaps, more preferable in practice.

The analytical treatment described above is rather idealized as it does not include the beam particle movement dynamics; however, it shows parameters to which there can try to approach. The result of initial simulation of a one-stage DC for a spent drive beam of the S-band EVT is shown in Figure 18 . No special optimization of the DC geometry and the magnetic field shape was carried out. Secondary emission was not taken into consideration either. The annular 2D model of the drive beam has been used.

Figure 18.

Geometry of a one-stage DC for the S-band EVT. A picture of particles energy is added at the bottom of the figure.

The results of initial simulations have shown a DC efficiency about 25%, which increases the total efficiency of the S-band EVT up to 72%. It certainly is substantially lower than the idealized theoretical results of the previous section and allows further optimization.

9. Conclusion

Theoretical and numerical research of properties of the linear accelerator EVT has shown that there is sufficient comprehension about it properties. There is no visible vagueness in technical problems. Certainly, the EVT accelerator does not show a record rate of acceleration; however, it shows that to reach a high efficiency of 70% or more is possible. It is much higher than the total efficiency of a klystron loaded with an ordinary linear accelerator. The two-stage configuration of the accelerator considerably increases the transformation ratio, but the accelerator’s efficiency decreases; therefore, further optimization of this configuration is preferable. Thus, there are all preconditions in place to proceed with the engineering design.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Vladimir E. Teryaev (December 20th 2017). Concept and Numerical Simulations of Multi-Beam Linear Accelerators EVT with Depressed Collector for Drive Beam, Accelerator Physics - Radiation Safety and Applications, Ishaq Ahmad and Maaza Malek, IntechOpen, DOI: 10.5772/intechopen.70956. Available from:

chapter statistics

367total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Optically Controlled Laser-Plasma Electron Acceleration for Compact γ-Ray Sources

By Serge Y. Kalmykov, Xavier Davoine, Isaac Ghebregziabher and Bradley A. Shadwick

Related Book

First chapter

Introductory Chapter: Introduction to Ion Implantation

By Ishaq Ahmad and Waheed Akram

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us