# Electromagnetic Waves Dispersion and Interaction of an Annular Beam-Ion Channel System in Plasma Waveguide

Authors: Jixiong Xiao; Zhijiang Wang; Jun Chen; Changhai Liu
Journal: Advances in Mathematical Physics
 (2017)

---

## Abstract

A linear theory for the electromagnetic properties and interactions of an annular beam-ion channel system in plasma waveguide is presented. The dispersion relations for two families of propagating modes, including the electrostatic and transverse magnetic modes, are derived. The dependencies of the dispersion behavior and interaction for different wave modes on the thickness of the annular beam and betatron oscillation frequency are studied in detail by numerical calculations. The results show that the inner and outer radii of the beam have different influences on propagation properties of the electrostatic and electromagnetic modes with different betatron oscillation parameters. In the weak ion channel situation, the two types of electrostatic waves, that is, space charge and betatron modes, have no interaction with the transverse magnetic modes. However, in the strong ion channel situation, the transverse magnetic modes will have two branches and a low frequency mode emerged as the new branch. In this case, compared with the solid beam case, the betatron modes not only can interact with the high frequency branch at small wavenumber but also can interact with the low frequency branch at large wavenumber.


---

## Body

## 1. Introduction

As is well known, plasma filling has a variety of advantages in increasing the space charge limited current, overall energy conversion efficiency, and radiation bandwidth dramatically [1–3] and it has been widely used in many plasma microwave radiation sources, such as traveling-wave tubes (TWT) [4], backward-wave oscillators (BWO) [5], klystrons [6], and gyrotrons [7]. Compared with the vacuum electronic devices, the EM dispersion characteristic and beam-wave energy transfer mechanism will become more complex [8]. As it is the basis of the application of any plasma electronics devices, the EM dispersion behavior and interaction in such plasma apparatus have always been research hot spots in the past two decades. In the past years, most of the published papers concentrated on the case of the relativistic electron beam (REB) with an external magnetic field guiding. In recent years, ion channel guiding has been put forward as an innovative focusing method for guiding the electron bunch transport [9, 10]. To form an ion channel, the electron beam density nb and the plasma density np must satisfy the matching condition nb/γ2≤np<nb [9, 11]. Here, γ is Lorentz factor of the beam. As a REB passes through the preformed plasma, the beam front will push out the plasma electrons continuously by the space charge force induced by beam electron, leaving the almost immobile positive ions to form the so-called ion channel. Then the electron beam can transport reliably with the guidance of the focus force of the ion channel. This new focusing method has been experimentally demonstrated in plasma wave tubes (PWT) [12], free electron lasers (FEL) [13], and charged particle accelerators [14] successfully. In previous works, the high frequency azimuthally symmetrical and nonsymmetrical eigenmodes in a REB with ion-channel guiding have been analyzed by Rouhani et al. [15, 16]. Wang et al. have studied the dispersion relations of EM waves in beam-ion channel system and the mechanism of Cherenkov EM instability [17, 18]. Mirzanejhad et al. have presented the dispersion characteristics of the space charge waves in a uniform and rigid rotation REB in ion channel without betatron oscillations [19].However, these papers mainly focused on the solid beam case. As annular beam has higher space charge limiting current and beam-wave energy conversion efficiency, it is also used widely. Therefore, the aim of the present article is the investigation of the EM wave dispersion and interaction of an annular beam-ion channel system in cylindrical plasma waveguide. There are three types of eigenmodes that will be discussed and the electrostatic (ES) and EM approximation will be used to study the space charge modes, betatron modes, and transverse magnetic (TM) modes, respectively. After dispersion relations are derived, the influences of the beam bunch thickness and betatron oscillation parameter on the dispersion propagation property and interaction of different wave modes are discussed in detail by numerical calculations. The differences between the solid and annular beam are also given by a comparison.The remainder of the article is structured as follows. In Section2, the governing equations and the solutions of the wave equations in different regions for ES and TM waves are derived in the beam coordinate system and the dispersion relations are also presented. In Section 3, a numerical study of azimuthally symmetric eigenmodes is presented. In Section 4, some conclusions are made.


## 2. Basic Equations and Dispersion Relations

A cylindrical metallic waveguide with radiusR is completely filled with an annular plasma column. This plasma column is with inner radiusa and outer radiusR. A relativistic annular electron beam with inner radiusa and outer radiusb passes through it with an initial velocity vb0z along with the waveguide axis. The annular beam expels the plasma electron continuously and forms an ion channel with radiusc in the background plasma. The schematic diagram has been shown in Figure 1.Figure 1
The cross section of the annular beam-ion channel system. The annular electron beam with inner radiusa and outer radius b, ion channel with radius c, and metallic waveguide with radius R.In the ion channel, the beam is subjected to three forces in the radial direction, which are originated from the actions of the transverse electric fieldEi produced by the positive ion core, the space charge field Ebs of the REB,(1)α=ωb2ω2-γωi2/2-1,Ax,y=J0xN0y-J0yN0x,Bx,y=J0xN1y-J1yN0x,Cx,y=J1xN1y-J1yN1x,Dx,y=I0xK0y-I0yK0x,Ex,y=I0xK1y+I1yK0x,and the azimuthal self-magnetic fieldBbs induced by the beam current. Applying the Gauss’ theorem and Ampere’s theorem, Ei, Ebs, and Bbs can be expressed as follows, respectively:(2)Ei=enir2ε0er,Ebs=-enb2ε0r2-a2rer,Bbs=-evb0znb2c2ε0r2-a2reθ.Here, ni is the ion density and is equal to the plasma density np, nb is the electron beam density, e is the unit charge, r is the radial coordinate, er is the radial unit vector, and eθ is the azimuthal unit vector. The radial force balance equation is(3)γmedr2dt2=-eeni2ε0r-enb2ε0-enbvb0z22ε0c2r2-a2r.If the inner radius a approaches zero, the annular beam will reduce to the solid beam case and the solution has been given in [15], which can be written as r=r0cos⁡ωit/2γ, and r0 denotes the initial radial coordinate. The equation of motion of the electron beam is(4)dpdt=-eE+v×B.Note thatd/dt=∂/∂t+(v·∇). Here, we only discuss linear waves; it means that the amplitude of oscillation of the waves is small and the terms containing higher powers of amplitude factors can be neglected. Therefore, we can assume that all physical quantities can be expressed as f=f0+f1exp⁡[j(ωt-kzz+lθ)], f0 and f1 represent the equilibrium quantity and perturbation quantity, respectively, and f1≪f0. Therefore, we have r=rb0+r1, v=vb0z+vb1, E=Ei+Ebs+E1, and B=Bbs+B1. If introducing cylindrical coordinates (r,θ,z) and z-axis coincides with the waveguide axes, the transverse and longitudinal motion equation of the beam electrons can be, respectively, expressed as(5)∂∂t+vb0z+vb1·∇vb1⊥=-eγmeE1⊥+npe2ε0rb0+rb1-nbe2ε0rb0+rb1-a2errb0+rb1+vb0z+vb1×B1⊥-nbevb0z2ε0c2vb0z+vb1×rb0+rb1-a2errb0+rb1⊥,∂∂t+vb0z+vb1·∇vb1z=-eγmeE1z+vb0z+vb1×B1z-nbevb0z2ε0c2vb0z+vb1×rb0+rb1-a2rb0+rb1z+evb0z2γmec2E1z.In the above derivation we have considered the relativistic effects and must notice that(6)dpdt=γme∂vb1∂t+vb1·∇vb1+mevb0z∂γ∂t,mevb0z∂γ∂t=vb0zc2dWdt.Here, W=γmc2 is the relativistic energy. Using the relation dW/dt=F·V results in (7)mevb0z∂γ∂t=-evboz2c2E1z.In (5), rb0 is used to describe the betatron oscillation in the radial direction and rb1 is used to describe the radial displacement induced by high frequency fields. The subscripts “⊥” and “z” denote the transverse and longitudinal component, respectively. If without disturbances all the beam electrons are nearly on the z-axis at the initial time, the perturbation amplitude caused by betatron oscillation will be nearly identical with the perturbation induced by high frequency fields so that both of them can be considered as a perturbation. Based on this assumption, after linearization, we can obtain the first-order motion equation: (8)∂∂t+vb0z∂∂zvb1=-eγmeE1⊥+vb0z×B1-ωi22γr1-eE1zγ3me.In the comoving coordinate of the beam, although the background plasma will be equivalent to a moving plasma column, the processing of the moving electron beam can be avoided, which is beneficial to simplify the derivation. Transformation to the beam frame (8) can be written as(9)∂∂tvb1=-emeE1-γωi22r1.


### 2.1. Electrostatic Modes

The Poisson equation is(10)∇2φ~=enb1ε0.The linearized continuity equation is(11)∂∂tnb1+nb0b∇·vb1=0,where nb0b=nb/γ is the unperturbed density of the REB in the beam frame.Noticing thatv⊥=dr/dt=-jωr, the transverse component of (9) can be written as(12)∂∂t-γωi22jωvb1⊥=-eme∇⊥φ~.By using (10), (11), and (12), the ES waves equation in the beam region can be obtained as (13)∇⊥2φ~+T2φ~=0,where(14)T2=kz2ω2-γωi2/2ω2-ωb2ω2-ω2+ωb2+γωi2/2and ωb is the beam frequency in beam frame. For the vacuum and ion channel regions, we can get(15)∇⊥2φ~-kz2φ~=0.Therefore, the suitable solutions in the different regions in the waveguide can be expressed as, respectively,(16)φ~=c1Ilkzr0<r≤ac2JlTr+c3NlTra<r≤bc4Ilkzr+c5Klkzrb<r≤R,where Jl is the Bessel function of the first kind, Nl is the Bessel function of the second kind, and Il and Kl are modified Bessel functions. Applying the appropriate boundary conditions [20], we can get a set of equations for five unknown coefficients. Eliminating the unknown coefficients, the ES wave dispersion relation can be obtained as(17)I1kzaI0kza=αTkzkzEkzR,kzbBTb,Ta+TαDkzR,kzbCTb,TakzEkzR,kzbATb,Ta-TαDkzR,kzbBTa,Tb.Here, only symmetric modes (l=0) are considered and, in the above equation, we have introduced(18)α=ωb2ω2-γωi2/2-1,Ax,y=J0xN0y-J0yN0x,Bx,y=J0xN1y-J1yN0x,Cx,y=J1xN1y-J1yN1x,Dx,y=I0xK0y-I0yK0x,Ex,y=I0xK1y+I1yK0x.Equation (17) is the dispersion relation for space charge and betatron waves with ES approximation. If the annular beam fills the waveguide completely, the dispersion relation will reduce to (19)I1kzaI0kza=TαkzBTR,TaATR,Ta.If the inner radius of the annular beam is close to zero, the annular beam will reduce to the solid beam case (20)J1TbJ0Tb=kzTαEkzR,kzbDkzR,kzb.


### 2.2. TM Modes

After linearization, Maxwell’s equations are(21)∇×E1=jωB1,(22)∇×B1=-μ0enb0bvb1-jωμ0ε0E1,(23)∇·E1=-enb1ε0,(24)∇·B1=0.By using (21) and (22), we can find the transverse and longitudinal equation for E1 as follows:(25)-∂∂z∂E1r∂z-∂E1z∂r-ω2c2E1r=-jωμ0enb0bvb1r,1r∂∂rr∂E1r∂z-∂E1z∂r-ω2c2E1z=-jωμ0enb0bvb1z.From (9) and (25), the TM wave equation in the beam region can be obtained as(26)∇⊥2E1z+T2E1z=0,(27)T2=ω2-ωb2kz2γωi2/2-ω2-ω2/c2ωb2-ω2+γωi2/2ω2ω2-ωb2-γωi2/2.In the beam coordinate system, the quasi-static background plasma can be viewed as a moving medium with a velocity vb0z along the z-axis negative direction. Using linear theory, we can get the expressions of the transverse and longitudinal perturbation current after a simple derivation, respectively,(28)Jp1⊥=-jωε0-ωp2ω2Ep1⊥-jωp2ω+kzvb0zvb0zω2∇⊥Ep1z,Jp1z=-jωε0-ωp2ω+kzvb0zγ2+ωp2vb0z2ω+kzvb0z2ω2∇⊥2Ep1z-jωp2vb0zω+kzvb0zω2∇⊥·Ep1⊥.Ampere’s circuital law is (29)∇×B1=μ0J1-jωε0μ0E1=-jωμ0ε→→·E1.Therefore, the equivalent dielectric tensor of the plasma column can be obtained as (30)ε→→=ε0ε10jε2∇r0ε1jε2∇θjε2∇rjε2∇θε3,where ε1=1-ωp2/ω2, ε2=-ωp2vb0z/ω+kzvb0zω2, and ε3=1-ωp2/ω+kzvb0z2γ2+ωp2vboz2/ω+kzvb0z2ω2∇⊥2.And the wave equation in the plasma region can be obtained as(31)∇⊥2E1z+p2E1z=0;here, p2=ω2-ωp2/c2-kz2.For high frequency EM waves, the ion channel region can be regarded as vacuum and the wave equation in this region is(32)∇⊥2E1z+τ2E1z=0and τ2=ω2/c2-kz2.The appropriate solutions for (26), (31), and (32) in the different regions of the waveguide are of the form(33)E1z=c1Jlτr0<r≤a,E2z=c2JlTr+c2NlTra<r≤b,E3z=c4Jlτr+c5Nlτrb<r≤c,E4z=c6Jlpr+c7Nlprc<r≤R.Applying the appropriate boundary conditions [20], the TM modes dispersion equation can be obtained as follows: (34)J0′τaJ0τa=A0A1BpR,pcX1+ApR,pcX2ApR,pcX3-A1BpR,pcX4.Here, only symmetric modes (l=0) are considered. In the above equation, we have introduced(35)A0=τTω2-ωb2ω2,A1=ε1-ε2kzτp,X1=A0CTb,TaAτb,τc+BTb,TaBτc,τb,X2=BTb,TaCτb,τc-A0CTb,TaBτb,τc,X3=A0BTa,TbBτb,τc-ATa,TbCτb,τc,X4=ATa,TbBτc,τb+A0BTa,TbAτb,τc.In laboratory coordinate system the expression of the ion channel radius is of the form r=ΔRnb/np1/2 [21] and ΔR is the thickness of the beam. It will increase with the thickness of the beam ΔR or the reciprocal of the plasma neutralization factor f=np/nb. If it fills the waveguide completely, in this case, the dispersion relation will be rewritten as(36)J0′τaJ0τa=τTω2-ωb2ω2τ/Tω2-ωb2/ω2CTb,TaAτR,τb-BτR,τbBTb,TaBτR,τbATa,Tb-τ/Tω2-ωb2/ω2BTa,TbAτR,τb.When the waveguide filled with the annular beam completely, (36) will reduce to (37)J0′τaJ0τa=τTω2-ωb2ω2BTR,TaATR,Ta.If the inner radius of the annular beam is close to zero, the annular beam will reduce to the solid beam case (38)J0′TbJ0Tb=1A0A1BpR,pcBτc,τb+ApR,pcCτb,τcBτb,τcApR,pc-A1BpR,pcAτb,τc.


## 2.1. Electrostatic Modes

The Poisson equation is(10)∇2φ~=enb1ε0.The linearized continuity equation is(11)∂∂tnb1+nb0b∇·vb1=0,where nb0b=nb/γ is the unperturbed density of the REB in the beam frame.Noticing thatv⊥=dr/dt=-jωr, the transverse component of (9) can be written as(12)∂∂t-γωi22jωvb1⊥=-eme∇⊥φ~.By using (10), (11), and (12), the ES waves equation in the beam region can be obtained as (13)∇⊥2φ~+T2φ~=0,where(14)T2=kz2ω2-γωi2/2ω2-ωb2ω2-ω2+ωb2+γωi2/2and ωb is the beam frequency in beam frame. For the vacuum and ion channel regions, we can get(15)∇⊥2φ~-kz2φ~=0.Therefore, the suitable solutions in the different regions in the waveguide can be expressed as, respectively,(16)φ~=c1Ilkzr0<r≤ac2JlTr+c3NlTra<r≤bc4Ilkzr+c5Klkzrb<r≤R,where Jl is the Bessel function of the first kind, Nl is the Bessel function of the second kind, and Il and Kl are modified Bessel functions. Applying the appropriate boundary conditions [20], we can get a set of equations for five unknown coefficients. Eliminating the unknown coefficients, the ES wave dispersion relation can be obtained as(17)I1kzaI0kza=αTkzkzEkzR,kzbBTb,Ta+TαDkzR,kzbCTb,TakzEkzR,kzbATb,Ta-TαDkzR,kzbBTa,Tb.Here, only symmetric modes (l=0) are considered and, in the above equation, we have introduced(18)α=ωb2ω2-γωi2/2-1,Ax,y=J0xN0y-J0yN0x,Bx,y=J0xN1y-J1yN0x,Cx,y=J1xN1y-J1yN1x,Dx,y=I0xK0y-I0yK0x,Ex,y=I0xK1y+I1yK0x.Equation (17) is the dispersion relation for space charge and betatron waves with ES approximation. If the annular beam fills the waveguide completely, the dispersion relation will reduce to (19)I1kzaI0kza=TαkzBTR,TaATR,Ta.If the inner radius of the annular beam is close to zero, the annular beam will reduce to the solid beam case (20)J1TbJ0Tb=kzTαEkzR,kzbDkzR,kzb.


## 2.2. TM Modes

After linearization, Maxwell’s equations are(21)∇×E1=jωB1,(22)∇×B1=-μ0enb0bvb1-jωμ0ε0E1,(23)∇·E1=-enb1ε0,(24)∇·B1=0.By using (21) and (22), we can find the transverse and longitudinal equation for E1 as follows:(25)-∂∂z∂E1r∂z-∂E1z∂r-ω2c2E1r=-jωμ0enb0bvb1r,1r∂∂rr∂E1r∂z-∂E1z∂r-ω2c2E1z=-jωμ0enb0bvb1z.From (9) and (25), the TM wave equation in the beam region can be obtained as(26)∇⊥2E1z+T2E1z=0,(27)T2=ω2-ωb2kz2γωi2/2-ω2-ω2/c2ωb2-ω2+γωi2/2ω2ω2-ωb2-γωi2/2.In the beam coordinate system, the quasi-static background plasma can be viewed as a moving medium with a velocity vb0z along the z-axis negative direction. Using linear theory, we can get the expressions of the transverse and longitudinal perturbation current after a simple derivation, respectively,(28)Jp1⊥=-jωε0-ωp2ω2Ep1⊥-jωp2ω+kzvb0zvb0zω2∇⊥Ep1z,Jp1z=-jωε0-ωp2ω+kzvb0zγ2+ωp2vb0z2ω+kzvb0z2ω2∇⊥2Ep1z-jωp2vb0zω+kzvb0zω2∇⊥·Ep1⊥.Ampere’s circuital law is (29)∇×B1=μ0J1-jωε0μ0E1=-jωμ0ε→→·E1.Therefore, the equivalent dielectric tensor of the plasma column can be obtained as (30)ε→→=ε0ε10jε2∇r0ε1jε2∇θjε2∇rjε2∇θε3,where ε1=1-ωp2/ω2, ε2=-ωp2vb0z/ω+kzvb0zω2, and ε3=1-ωp2/ω+kzvb0z2γ2+ωp2vboz2/ω+kzvb0z2ω2∇⊥2.And the wave equation in the plasma region can be obtained as(31)∇⊥2E1z+p2E1z=0;here, p2=ω2-ωp2/c2-kz2.For high frequency EM waves, the ion channel region can be regarded as vacuum and the wave equation in this region is(32)∇⊥2E1z+τ2E1z=0and τ2=ω2/c2-kz2.The appropriate solutions for (26), (31), and (32) in the different regions of the waveguide are of the form(33)E1z=c1Jlτr0<r≤a,E2z=c2JlTr+c2NlTra<r≤b,E3z=c4Jlτr+c5Nlτrb<r≤c,E4z=c6Jlpr+c7Nlprc<r≤R.Applying the appropriate boundary conditions [20], the TM modes dispersion equation can be obtained as follows: (34)J0′τaJ0τa=A0A1BpR,pcX1+ApR,pcX2ApR,pcX3-A1BpR,pcX4.Here, only symmetric modes (l=0) are considered. In the above equation, we have introduced(35)A0=τTω2-ωb2ω2,A1=ε1-ε2kzτp,X1=A0CTb,TaAτb,τc+BTb,TaBτc,τb,X2=BTb,TaCτb,τc-A0CTb,TaBτb,τc,X3=A0BTa,TbBτb,τc-ATa,TbCτb,τc,X4=ATa,TbBτc,τb+A0BTa,TbAτb,τc.In laboratory coordinate system the expression of the ion channel radius is of the form r=ΔRnb/np1/2 [21] and ΔR is the thickness of the beam. It will increase with the thickness of the beam ΔR or the reciprocal of the plasma neutralization factor f=np/nb. If it fills the waveguide completely, in this case, the dispersion relation will be rewritten as(36)J0′τaJ0τa=τTω2-ωb2ω2τ/Tω2-ωb2/ω2CTb,TaAτR,τb-BτR,τbBTb,TaBτR,τbATa,Tb-τ/Tω2-ωb2/ω2BTa,TbAτR,τb.When the waveguide filled with the annular beam completely, (36) will reduce to (37)J0′τaJ0τa=τTω2-ωb2ω2BTR,TaATR,Ta.If the inner radius of the annular beam is close to zero, the annular beam will reduce to the solid beam case (38)J0′TbJ0Tb=1A0A1BpR,pcBτc,τb+ApR,pcCτb,τcBτb,τcApR,pc-A1BpR,pcAτb,τc.


## 3. Numerical Results

In this part, a numerical calculation is used to analyze the effects of the inner and outer radii on the dispersion curves of the azimuthally symmetric(l=0) ES and TM modes and the interaction between them under different EM characteristics of the ion channel. For the convenience of the numerical calculation, The dispersion frequency ω and the wavenumber kz are normalized by the beam frequency ωb and the speed of the light in vacuum; that is, ω~=ω/ωb and k~z=kzc/ωb. The normalized waveguide radius has the value Rωb/c=1.2, which corresponds toR= 2 mm and ωb = 1.8 × 1011 rad/s with the beam density nb0b = 1.033 × 1013 cm−3 in the beam coordinate system and these values are consistent with [15].From (14) and (27), we can find that the transverse oscillation characteristic frequency ωβ=γ/2ωi has important influence on the solution of the ES and TM mode dispersion equation. For example, ωβ is larger or smaller than the beam frequency ωb; the dispersion equations of the ES and TM modes will have different solutions and it leads to the wave modes in both cases having different dispersion characteristics. Based on this point it can be classified as strong and weak ion channel and in the cases of the weak and strong ion channel, which meet 1/γ<ωp/ωb≤3/γ and 3/γ<ωp/ωb≤γ, respectively.The asymptotic behaviors and the influences of the inner and outer radii on the space charge and betatron modes for strong ion channel are presented in Figures2–5. In this case the dispersion curves of the SC01 waves will shift to the beam frequency ωb as the wavenumber is close to infinity and the Be01 waves have upper hybrid frequency ωH=(ωb2+γωi2/2)1/2 [15]. It can also be found that, as the inner radius decreased or the outer radius increased, the SC01 modes experience obviously shifts close to higher frequencies but the case of the Be01 modes is just to the opposite and its dispersion curves shift down clearly. Besides, compared with the solid beam case, the SC01 modes have higher frequencies and the Be01 modes have lower frequencies in the annular beam case.Figure 2
The influence of the ratioa/b on the dispersion curves of the SC01 waves for strong ion channel with ωβ=2.3195. The outer radius of the annular beam is taken as 0.5R and a/b=0 represents the solid beam case.Figure 3
The influence of the ratioa/b on the dispersion curves of the Be01 waves for strong ion channel with ωβ=2.3195. The outer radius of the annular beam is taken as 0.5R and a/b=0 represents the solid beam case.Figure 4
The influence of the ratiob/R on the dispersion curves of the SC01 waves for strong ion channel with ωβ=2.3195. The inner radius of the annular beam is taken as 0.25R. b/R=1 represents the waveguide completely filled with the annular beam.Figure 5
The influence of the ratiob/R on the dispersion curves of the Be01 waves for strong ion channel with ωβ=2.3195. The inner radius of the annular beam is taken as 0.25R. b/R=1 represents the waveguide completely filled with the annular beam.Figures6–9 have shown the asymptotic behaviors and the effects of the inner and outer radii on space charge and betatron modes for weak ion channel. In this case SC01 waves have upper hybrid frequency ωH=(ωb2+γωi2/2)1/2 and will approach the beam frequency ωb with the wavenumber being close to infinity, and Be01 waves have the cutoff frequency ωβ. We can also see that from the figures, in contrast to the strong channel case, the roles of the space charge and betatron modes will be exchanged in weak ion channel so that the effects of the inner and outer radii on the space charge and betatron modes are also exchanged.Figure 6
The effect of the ratioa/b on the dispersion curves of the SC01 waves for weak ion channel with ωβ=0.8761. The outer radius of the annular beam is taken as 0.5R and a/b=0 represents the solid beam case.Figure 7
The effect of the ratioa/b on the dispersion curves of the Be01 waves for weak ion channel with ωβ=0.8761. The outer radius of the annular beam is taken as 0.5R and a/b=0 represents the solid beam case.Figure 8
The effect of the ratiob/R on the dispersion curves of the SC01 modes for weak ion channel with ωβ=0.8761. The inner radius of the annular beam is taken as 0.25R. b/R=1 represents the waveguide completely filled with the annular beam.Figure 9
The effect of the ratiob/R on the dispersion curves of the Be01 modes for weak ion channel with ωβ=0.8761. The inner radius of the annular beam is taken as 0.25R. b/R=1 represents the waveguide completely filled with the annular beam.Figures2–9 also show that, in both of the strong and weak ion channel case, the asymptotic behaviors of the SC01 and Be01 modes for the case of the solid beam in [15] coincide with the annular beam case in this paper. Besides, the beam radius of the solid beam and the outer radius of the annular beam have a similar effect on the space and betatron modes. As the solid beam can be regarded as a special case of the inner radius of the annular beam close to zero, we can theoretically expect such results.The influences of the inner radius on the TM01 waves in the cases of strong and weak ion channel are demonstrated in Figures 10 and 11. With the increasing of the inner radius (the ratio a/b), the dispersion curves of the TM01 modes will undergo slightly frequency shifts toward higher frequencies for both of the strong and weak ion channel and the upward shifts are more obscurely in strong ion channel case. In general, the EM eigenmode in an annular beam system has a higher phase velocity than the solid beam case with the same condition [22]. Therefore, the TM modes have higher frequencies for the solid case.Figure 10
The effect of the ratioa/b on the dispersion curves of the TM01 modes for strong ion channel with ωβ=1.6202. The outer radius of the annular beam is taken as b=0.4R and a/b=0 represents the solid beam case.Figure 11
The effect of the ratioa/b on the dispersion curves of the TM01 modes for weak ion channel with ωβ=0.8832. The outer radius of the annular beam is taken as b=0.4R and a/b=0 represents the solid beam case.The TM modes normalized dispersion frequencies as a function of the ratio of the outer radiusb to the waveguide radius R are illustrated in Figure 12. In this figure, ωβ/ωb=1.6205 and ωβ/ωb=0.8832 correspond to ωi/ωb=1.1456 and ωi/ωb=0.6245, respectively. As the radius of the ion channel is proportional to the thickness, that is, r=ΔRγωb/ωp in beam frame, the ion channel will fill the waveguide gradually with the increasing of the out radius. In this case, we can find that if b/R>0.46 or b/R>0.6625, the radius of the ion channel will be larger than the waveguide radiusR, which means that it will fill the waveguide completely. As this figure shows, while the waveguide is only filled with it partially, the TM modes frequencies will shift down gradually as the outer radius increased. However, if the waveguide is filled with it completely, this case is just to the opposite of the partially filled case and the frequencies will shift up.Figure 12
The effect of the ratiob/R on the dispersion curves of the TM modes for strong and weak ion channel with γ=4. The inner radius of the annular beam is taken as a=0.1R.Figure13 shows the influences of the betatron oscillation characteristic parameter ωβ on the TM dispersion properties for strong ion channel. With the increasing of ωβ, there will be two branches for TM modes and a low frequency (LF) mode with a lower cutoff frequency emerged as the new branch. We can also find that the dispersion behavior of this new branch is similar to the high frequency branch at small wavenumber, but it is similar to betatron modes at large wavenumber. It is supposed that the reason of this phenomenon is that the dispersion equation (34) will have a new solution with the increase of the betatron oscillation frequency. Besides, as the longitudinal wavenumber kz approaches infinity, the new LF branch will have longitudinal wave-like character, which may result in the new LF branch having betatron-like character. This is very different from the solid beam case and we did not find similar phenomenon in that case.Figure 13
The effect of betatron oscillation parameterωβ on the dispersion characteristics of the TM01 modes for strong ion channel. The inner and outer radius of the annular beam are taken as a=0.25R and b=0.5R. The Lorentz factor is γ=6.The interactions between the TM and ES modes for the cases of the annular beam and solid are presented in Figures14–17. As can be seen from the figures, the cutoff frequencies for the SC01 and Be01 waves are lower than the TM wave in the weak ion channel case so that there is no coupling between the three types propagating modes. However, the strong ion channel case is completely different. Although the asymptotic frequency ωb of the SC01 modes is far smaller than the cutoff frequency of the TM modes and it cannot interact with the TM modes, the Be01 modes can interact with the TM modes for both the annular and solid beam cases. Compared with the solid beam, from Figure 15, it can also be seen that the dispersion curve of the LF branch for the TM modes is close to the betatron modes gradually and the two dispersion curves almost coincident at large number wave. This is an interesting phenomenon. As the betatron waves are backward waves and have negative group velocity, the coupling between betatron and TM wave may produce absolute instability and this interaction can be used for microwave generation [23].Figure 14
In annular beam case, there is no coupling between the TM modes and ES waves in weak ion channel case. The inner and outer radius of the annular beam are taken asa=0.25R and b=0.5R; the Lorentz factor is γ=6.Figure 15
In annular beam case, the TM wave of the two branches can interact with the betatron wave in strong ion channel case. The inner and outer radius of the annular beam are taken asa=0.25R and b=0.5R; the Lorentz factor is γ=6.Figure 16
In solid beam case, there is no coupling between the TM and ES waves in weak ion channel case. The beam radius is taken asa=0.25R; the Lorentz factor is γ=6.Figure 17
For solid beam, the TM wave only has one branch and can interact with the betatron wave in strong ion channel case. The beam radius is taken asa=0.25R; the Lorentz factor is γ=6.


## 4. Conclusion

In present analysis, the EM dispersion properties and interactions of an annular beam-ion channel system in plasma waveguide are studied by linear perturbation theory. The influences of the inner and outer radii of the beam and betatron oscillation parameter on the dispersion characteristics of space charge modes, betatron modes, and TM modes are revealed by numerical calculations. It is shown that the inner and outer radii have different influences on the three types of waves under different betatron characteristic frequencies and the TM modes will have two branches in strong ion channel case with appropriate betatron frequency. The investigation of the interaction between the ES waves and TM waves for strong and weak ion channel indicates that there is no coupling between the two families propagating modes in weak ion channel case but the Be01 modes can interact with the both of the two branches of the TM01 waves in the strong ion channel case. The coupling between the low frequency branch and the betatron wave especially is a new phenomenon and the details still need further study. Obviously, it can be found from the conclusions of this paper that the ion channel on the one hand plays an important role of guiding the beam to ensure the transport effectively. On the other hand, different betatron oscillation characteristic frequency of the ion channel can make the propagating wave modes in the system show different EM characteristics and change the interaction between them. In addition, the physical model used in this paper also has some limitations; for example, the influence of the electron duration is not taken into account. On the one hand, for a short pulse electron beam, when the beam is injected into the plasma, a return current will be excited in it and the direction of the return current is opposite to the beam current so that it may have a significant impact on the balance of the electron beam in the radial direction and also can affect the eigenstructure of the system. On the other hand, for a long pulse electron, there may be oscillating charge on the boundary of the ion channel and then the oscillating current may be formed [24]. It may also affect the dispersion characteristics of the system. Therefore, these works still need further analyses in the future.



---

*Source: 8141549-2017-11-20.xml*