PHOTONIC CRYSTAL RIBBON-BEAM TRAVELING WAVE AMPLIFIER
PRIORITY INFORMATION This application claims priority from provisional application Ser. No. 60/483,852 filed June 30, 2003, which is incorporated herein by reference in its entirety.
BACKGROUND OF THE INVENTION The invention relates to the field of optical communication, and in particular to a photonic crystal ribbon-beam traveling wave amplifier. The third-generation (3G) wireless communication standards call for hardware-based upgrade to the second-generation (2G) Global Systems for Mobile Communications using Wideband Code Division Multiple Access (W-CDMA) and Universal Mobile Telephone System (UMTS) as well as software-based upgrade to 2G Code Division Multiple Access (CDMA). The 3G wireless communications require amplifiers operating frequencies that are 1.12 to 3 times that of present frequencies, which are in the range from 900 MHz to 1700 MHz. In general, the bandwidth of a transmitter, which is the most important figure of merit, increases with the central frequency of the amplifier. However, the number of transmitting towers must increase as the square of the central frequency, while keeping the power of the transmitting tower at a constant. This is because the distance between two adjacent transmitting towers is inversely proportional to the frequency. For example, if 1-GHz transmitting towers have a spacing of 10 miles, then 2-GHz transmitting towers must have a spacing of 5 miles. In other words, four 2-GHz transmitting towers are required to cover 100 square miles, whereas only one 1-GHz transmitting tower is needed for the same area. Moreover, the total RF power per unit area increases with increasing data rate. For example, 3G wireless networks are expected to have considerably higher data rate than 2G wireless networks. As a result, the number of power amplifiers increases more dramatically than the square of the carrier frequency. In order for the third-generation (3G) and future wireless communications to be a viable business, it is essential for the telecommunication equipment industry to provide ultra-low-cost amplifiers. At present, most wireless base stations are powered by solid-state power
amplifiers, which operate with efficiencies in the 8-12% range. The cost of solid-state amplifier is about $100 Watt. The cost of power amplifiers per base station at 1.5 kW is $150,000. For the RF power part of a wireless base station, the operating cost is comparable to the capital cost, because of the low operating efficiencies and heat removal. Conventional helix traveling wave tubes (TWTs), which are not employed in any existing wireless base stations, cannot meet the ultra-low-cost requirement set by any potential third-generation wireless infrastructure provider. For example, a 100 W, 2 GHz conventional helix TWT costs $20K a piece or more. There is a need to develop high efficiency, low-cost microwave amplifiers for
3G and future wireless base stations.
SUMMARY OF THE INVENTION According to one aspect of the invention, there is provided a RF amplifier. The RF amplifier includes a RF input section for receiving a RF input signal. At least one slow-wave structure associated with the RF interaction section. An electron ribbon beam interacts with the RF input supported by the at least one slow-wave structure so that the kinetic energy of the electron beam is transferred to the RF fields of the RF input signal, thus amplifying the RF input signal. A RF output section outputs the amplified RF input signal. According to another aspect of the invention, there is provided a method of forming a RF amplifier. The method includes forming RF input section for receiving a RF input signal. Also, the method includes forming al least one photonic crystal for operational control if necessary. An electron ribbon beam is formed that interacts with the RF input supported by the at least one slow-wave structure so that the kinetic energy of the electron beam is transferred to the RF fields of the RF input signal, thus amplifying the RF input signal. Furthermore, the method includes forming a RF output section that outputs the amplified RF input signal. BRIEF DESCRIPTION OF THE DRAWINGS FIGs. 1 A- IB are schematic diagrams of a photonic crystal ribbon-beam traveling wave amplifier (PCRB TWA) in accordance with the invention; FIG. 2 is a schematic diagram illustrating a detailed depiction of a single-sided slow-wave structure;
FIG. 3 is a schematic diagram illustrating a detailed depiction of a double- sided slow-wave structure; FIGs. 4A-4B show plots of the dispersion functions Da {ω, k) for the antisymmetric modes as they vary with frequency at several values of phase shifts; FIGs. 5A-5B show plots of the dispersion functions Ds (ω, k) for the symmetric modes as they vary with frequency at several values of phase shifts; FIG. 6 is a graph illustrating eigenfrequencies for the two lowest bands of antisymmetric modes and the two lowest bands of symmetric modes; FIG. 7 is a graph illustrating the comparison between the eigenfrequencies for the two lowest anti-symmetric modes with the same system parameters as in FIG. 6; FIG. 8 is a schematic diagram illustrating a ribbon electron beam interacting with a single-sided PC slow-wave structure; FIGs. 9A-9B are graphs illustrating the temporal linear growth rate of the lowest anti-symmetric mode as a function of the wavenumber; FIGs. 10A-10B are graphs illustrating the spatial linear growth rate of the lowest anti-symmetric mode as a function of angular frequency; FIG. 11 is a graph demonstrating the maximum temporal growth of the lowest anti-symmetric mode as a function of normalized beam current; FIG. 12 is a graph demonstrating the maximum spatial growth of the lowest anti-symmetric mode as a function of normalized beam current; FIGs. 13A-13B is a graph demonstrating the dispersion characteristics in a two-dimensional dielectric photonic crystal; FIG. 14 is a schematic diagram of the design of the PCRB TWA for 3G wireless base stations; FIG. 15 is a graph of the eigenfrequency versus the phase shift for the two lowest bands in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA structure; FIG. 16 is a graph of the normalized real wavenumbers RekzL versus the normalized frequency for the operating mode in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA structure; FIG. 17 is a graph of the normalized imaginary wavenumbers |lm&zZ,| versus the normalized frequency for the operating mode in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA structure; FIG. 18 is a graph of the real phase shift RekzL versus the frequency / for
the operating mode; FIG. 19 is a graph of the gain G versus the frequency / for the operating mode; and FIG. 20 is a graph demonstration the dispersion characteristics for 4 lowest- order TM modes in a photonic crystal.
DETAILED DESCRIPTION OF THE INVENTION The invention is a novel amplifier that employs two emerging technologies, namely, photonic crystals and low-density ribbon electron beams, in otherwise a conventional vacuum tube millimeter wave amplifier. FIGs. 1A-1B shows a schematic diagram of a photonic crystal ribbon-beam traveling wave amplifier (PCRB TWA) 2 in accordance with the invention. FIG. 1 A shows a double sided PCRB TWA 2 that includes a ribbon electron beam 10 propagating in the z -direction from the emitter 14 and extending out, wiggler magnets 8 for beam focusing, a photonic crystal (PC) slow- wave structure 12 with metallic or dielectric rods and plates, and RF input 4 and output 6 sections. As the electron beam 10 interacts with the RF input 4 supported by slow- wave structure 12, the kinetic energy of the electron beam is transferred to the RF fields, amplifying the RF signal 16. The amplified RF signal 18 exits the amplifier at the RF output 6, and the spent electron beam is collected down stream. FIG. IB shows a single sided PCRB TWA 20 that includes a ribbon electron beam 28 propagating in the z -direction from the emitter -32 and extending out, wiggler magnets 26 for beam focusing, a photonic crystal (PC) slow-wave structure 30 with metallic or dielectric rods and plates, and RF input 22 and output 24 sections. As the electron beam 28 interacts with the RF input signal 34 supported by slow-wave structure 30, the kinetic energy of the electron beam 28 is transferred to the RF fields, amplifying the RF signal 34 to produce the amplified signal 36. Unlike the round beams in conventional helix or .coupled-cavity TWTs, the ribbon electron beam in the PCRB TWA will reduce the magnetic field required for beam focusing, reduce the loading in the amplifier, increase the amplifier efficiency, and improve the amplifier linearity and bandwidth. The ribbon beam will have an aspect ratio of 1 to 10, which effectively lowers the beam perveance (or space-charge) by a factor of 10 in comparison with a round beam. Based on the well-known empirical scaling law, the
' efficiency of the PCRB
TWA is expected to be as high as 80%. Furthermore, because the effective beam perveance is small, the interaction between the electrons and- slow-wave structure is expected to be high. Consequently, a high degree of linearity of the amplifier is expected in the high-efficiency operation. Consider a monochromatic electromagnetic wave propagating in the two- dimensional single-sided slow- wave structure 38 as shown in FIG. 2. The slow- wave structure 38 consists of a metal plate at x = 0 and corrugated vanes 40 located between x = b and d . The period of corrugation is L , and the width of each vane is a . For simplicity, either the periodic metal or dielectric rods that are placed at x > d in the PCRB TWA or the variations in the y -direction is considered. Moreover, the transverse magnetic (TM) modes are of importance with the field distributions E(x,z,t)= e-
iωt [E
x (x,z)e
x + E
2 (x,z)e
2 ], Eq. 1 (x,z,t) = e-
iωlB
y(x,z)e
y . Eq. 2 where Ε(x,z,t) is the electric field, E
x(x,z) is the electric in the x-direction, E
z(x,z) is the electric field in the y-direction, B(x,z,t) is the magnetic field, and treat B
y(x,z) as the generating function. The wave equation can be expressed as
where ω is the frequency of the wave and c is the speed of the light. In cgs units, the electric field can be expressed as
Ex{x,z) = —— y{x,z), Eq. 4 tω dz
E
2(x,z) B (x,z). Eq. 5
Expressing the generating function as a Bloch-wave function of the form B
y(x,z,t) = e^-
a) f „(
M ,
' Eq. 6
one can rewrite the wave equation as
or
where u
n is the Floquet amplitude, and „ is an effective wavenumber.
Note that for
p
n is imaginary. In the corrugated-vane region, the usual approximation is adopted
By{x,z,t)=iEϋei{ sL-ωt)coι ' ~{x-d) Eq.ll
Ex(x,z,t) = 0, Eq.12
E2{x,z,t)≡E0ei{ sL-ωt)s (x-d) Eq.13 c for \z-sL \<al2 and b < Λ: < d. Here, s = 0, ±1, +2,--- Note that Eq. 13 assures
E. x=d = 0. In the drift region with 0 < x < b, Eq.8 implies d2u„ dx +Pn :0, Eq.14
which has a general solution of the form "
» M = Λ„ s (p
nx) + B
n c s{p
nx) , Eq.15 where A and B„ are constants. Because E,\ = 0, one must have A
n =0, Eq.16 or
Therefore, the RF fields in the drift region can be expressed as
By(x,z,t) = e*-— > f c s{pnx)ei2m L; Eq.18
Ex(x,z,t)=e^^∑^ f 2m^ k +- Bn cos(pnx)et2™,L , Eq.19
E2(x,z,i)= e^-^ ∑ —pnBn sin{pnXy Ϊlimz/L Eq. 20
for 0 < x < b . The amplitudes Bn in the drift region are related to the amplitude E0 in the corrugated-vane region by the continuity of the electric field at x = b , i.e., ϊπnzIL E0e''(fei) sin < a/2,
***
n Στr
∞ —ιωPnB P bY
0, otherwise, Εq. 21 where s = 0, ± 1, ±2,— Solving Εq. (3.1.21) for B
n yields
The vacuum dispersion relation can be derived with matching the Poynting flux at x = b . For present purposes, it is useful to introduce the admittance defined by
Substituting Eqs. 11 and 13 to Eq. 22, one obtains
Similarly, substituting Eqs. 18, 20 and 22 into Eq. 23, one obtains
By setting A
+(ω,k
2) = A
~(ω,k
2), Eq. 26 one arrives at the vacuum dispersion relation
for the electromagnetic wave in the single-side slow- wave structure. As pointed out earlier, when inequality of Eq. 10 holds, p
n is imaginary and cot(p
nb)
= coth p
n\b)
In general, Eq. 27 must be solved numerically. As an alternative to the single-side slow-wave structure, a monochromatic wave propagation in a double-side slow-wave structure 42 is considered, as shown in
FIG. 3. The bottom of the double-side symmetric slow- wave structure is a mirror image of the top of the single-side slow-wave structure shown in FIG. 2. In the double-side slow- wave structure, however, there is no metal at x = 0. In the double-side slow-wave structure, the vacuum dispersion relation for the anti-symmetric wave propagation with E(- X, Z, t) = -E(x, Z, t) is the same as the one given in Eq. 27. However, the vacuum dispersion relation for the symmetric wave propagation with E(- x, z, i) = Ε(x, Z, i) has a difference expression. Paralleling the analysis in
Section 3.1, and taking An ≠ 0 and Bn = 0 in Eq. 15, one can show that it is given by
Eq. 29 where p
n can be imaginary in which case,
The mode structures in the double-sided slow-wave structures are qualitatively different from those in the single-sided slow-wave structures, which makes the single- side slow-wave structure suitable for use in a PCRB TWA. To gain some quantitative understanding of the vacuum dispersion characteristics of electromagnetic wave propagation in the .double-side slow-wave
structure shown in FIG. 3, consider the first example with the choice of system parameters corresponding to: L = 0.24 cm, a/L = 0.8, ML -.14 d/L = 6.0.
This choice of system parameters can also represent a single-side slow-wave structure that supports only the anti-symmetric modes. FIGs. 4A-4B show plots of the dispersion functions D
a (ω, k) for the antisymmetric modes as they vary with frequency at several values of phase shifts (or wave numbers). The phase shift is equal to 360° x (k
2L/2π) . For a given phase shift, the zeros of the dispersion function correspond to the eigenfrequencies of the system. In this example, there is one zero below 20 GHz at least, as shown in FIG. 4A. FIGs. 5A-5B show plots of the dispersion functions D
s (ω, k) for the symmetric modes as they vary with frequency at several values of phase shifts (or wave numbers). The phase shift is equal to 360° (k
2L/2π). For a given phase shift, the zeros of the dispersion function correspond to the eigenfrequencies of the system. In this example, there is one zero below 20 GHz at least, as shown in FIG. 5 A. The eigenfrequencies for the two lowest bands of anti-symmetric modes and the two lowest bands of symmetric modes are plotted as a function of phase shift in FIG. 6. h the double-sided slow-wave structure, both anti-symmetric and symmetric modes exist. In the phase shift range from 90° to 270° , the anti-symmetric and symmetric modes are nearly degenerate in the first band as well as in the second band, as shown in FIG. 6. In the single-side slow-wave structure, however, only the antisymmetric modes shown in FIG. 6 exist, and they are well separated in frequency. This sharp contrast in the mode properties makes the single-sided slow-wave structure a better choice for use in the PCRB TWA, because there are no symmetric modes. The eigenfrequencies of the TWA can also be determined using SUPERFISH code. FIG. 7 shows the comparison between the eigenfrequencies calculated with SUPERFISH and Eq. 27 for the two lowest anti-symmetric modes with the same system parameters as in FIG. 6 or in Eq. 31. Only two lowest modes are plotted in FIG. 7. As can be seen from FIG. 7, the theoretical results and SUPERFISH
simulations are in excellent agreement. Therefore, the theory derived for PCRB TWA is applicable for designing an actual device. The ribbon beam, as shown in FIGs. 1 A, is an important component for the efficient operation of the invention. Note the ribbon beam uses a planar wiggler magnetic field. Let us first consider a ribbon beam in the planar wiggler field B
w(x) = -B„ [e
x co k
wx)co$(k
wz)→
2 nh(k
wx)sin(k
wz)] , Eq. 32 where B
w = constant , k
w = 2π/λ
w, and λ
w is the wiggler period. Introducing the vector potential
with the gauge condition B
w(x) = Vx A
w(x), Eq. 34 the Hamiltonian for the single-particle motion in cgs units,
can be expanded for P = mv
y - (e / C)A
W (X, Z) = 0 and \k
wx\ < 1 as H(x,z,P
x,P
y,P
2 )= H
0{P
2)+H
β{x,P
x). Eq. 36 In Eq. 36, the Hamiltonians
and
H
β(x,P
x) = ^-(
x 2 + m
e 2ω
β 2x
2) Eq. 38 2m
e describe the axial motion and the (transverse) betatron oscillations, respectively, and
is the betatron oscillation frequency. For P = 0 and
< 1 , the equations of motion are
Eq. 40 dt m„
Eq. 41 dt me
Eq. 42 dt dPx — - = -mt ωβ 2x . Eq. 43
Since v
2 = β
bc = constant , Eq. 44 the equation of motion for the betatron oscillations can be expressed alternatively as d
2x - + k
βx = 0 , Eq. 45 dz where ω. Eq. 46 ' β
bc is the betatron wavenumber. For simplicity, the coupling between the beam envelopes in the x - and y ■ directions is ignored, and express the beam envelope equation in the x -direction as
where ■ { X Eq. 48 yb ( Eq. 49 κ = 2A Eq. 50 Wl
Here, k
B is the Boltzmarm constant, T is the Kelvin temperature, and (χ) denotes \xfdxdydx' dy' tø ~ Eq. 53 fdxdydx'dy'
with x' = dx/dz , y'
' - dy/dz , and f(x,y,x',y') is the electron distribution function. In the zero-current limit, K = 0 and it follows from Eq. 47 that the equilibrium rms beam envelope is given by
or the wiggler field required for the beam equilibrium amplitude to occur is given by
which allows us to calculate B
w . Note that the wiggler period λ
w does not appear in Eq. 55. However, one must demand 2
wx
b < 1, Eq. 56 or λ
w ≥ 4πx
b Eq. 57 for the approximations in Eq. 51 to be valid. For example, taking x
b = 0.06 cm, β
b = 0.08, V = 1700 eV, Eq. 58 1
7 1 1200
Λ Λ T kj = — x = 0.1 eV, 2 " 40 300 one can obtain from Eq. 55, 5
W = 24.7 G, Eq. 59 which is easily achievable. For a finite beam current, K ≠ 0 and it follows from Eq. 47 that the rms beam envelope is given by
or
where
ξ = ^ Eq. 62 x
b is the aspect ratio of the ribbon beam. In this discussion, the value of ξ is fixed. For example, taking x
b = 0.06 cm, β
b = 0.08,
I
b = 0.076 A, = 10, one obtains from Eq. 63, 5
W = 93.9 G , Eq. 64 which is still easily achievable. Now, a ribbon electron beam 46 interacts with a single-sided PC slow-wave structure 48, similar to structure 38 of FIG. 2, shown in FIG. 8: The equilibrium of the beam 46 is described by = V
b 2 = β
bc
2 , Eq. 65
J = -eσ
bV
bδ(x-h)e
2 = -δ(x-h)e
Υ , Eq. 66 dy where V
b = constant is the equilibrium beam velocity, σ
b ■ is the surface number density of the electrons in the ribbon beam, the sheet x = h specifies the transverse displacement of the beam, and
is the beam current, and y
b is the rms width of the ribbon beam. In this analysis, the variations are ignored in the v -direction, and treat the stability of the ribbon beam 46 as, a two-dimensional problem. The linearized cold- fluid equations are
— + V
b — }δσ + σ
b — δV = 0, Eq. 69 {dt dz)
b dz
τ
70 where the charge and current density perturbations δp(x,z,t) and δJ
2{x,z,t) are defined as δp(x,z,t)=-eδσ(z,t)δ(x-h), Eq.71 δJ
2 (x, z, t) = -eδσ(z, t)δ{x -h)V
b- eσ
bδ(x - h)δV
2 (z, t) . Eq.72 Expressing all perturbations as δψ(x
lz,t)=e*
k*-'
)'∑
lδψ.(x)e
a i2τmzl ' L Eq.73
or
one can rewrite the linearized cold-fluid equations 69-70 as i(ω-k
nV
b)δV
zn=-δE
zn{h), Eq.77 m (<»-KV
b)δσ
n-k
nσ
bδV
zn=0, Eq.78 δ(x — h)
Eq.79 where k
n and pi are defined as 2τm k =k + Eq.80 cr Pn ~ 2
Kn Eq.81
Solving for δVzn and δσn from Eqs.67 and 68 in terms of δEzn (h) yields
le SE2n{h) δV„ = - Eq.82 m(ω-kaVbY ie knσbδE2ll{h) Eq.83 m (ω-knVbf
Substituting Eqs.82 and 83 into Eq.79, one obtains
The solutions to Eq.84 are δE (
x)
= \
A"
sin(P»
x +B"
cos(Pn
x)> 0<x<h, Eq.85 [A sm(p
nx)+B cos(p
nx), h<x<b, where A*, B
* , Aζ and Bζ are constants to be determined by the boundary conditions at x = 0 , x = h , and x — b . The coefficients A
*, Bζ , Aζ and Bζ are determined in the following three steps. First, note the electric field δE
z(x,z ) vanishes at x = 0. This requires SE
z„(0) = 0. Eq.86
Therefore, one must have B = 0. Eq.87 Second, the coefficients Aζ and Bζ are expressed in terms of A
* , using both the continuity of the axial electric field at x = h , for example, Aζs (p
nh)+Bζcos(p
nh)-Aζ n(p
nh) 0 , Eq.88 and the relation
or
A cos{pnh)~ B s (pnh)- Aζ cos{pnh) = πe σ> P" ,2Aζ s (pnh) . m (ω-k„Vb)
Eq.90
Solving Eqs.88 and 90 for Aζ and Bζ gives Aζ ^ = l + an(ω,k)cos{pnh), Eq.91 A
10
-^- = -α
n(fi?, ;)sin ?
nft), Eq.92
An where the function a
n {ω,k) is defined by
Furthermore, the electric and magnetic fields are expressed in the region h<x<b as δE2{x,z,t)= ∑{Aζ sm(p,tx)+Bζ cos(p„x)]e''^-ω() , Eq.94 δEx(x,z,t) = —^-δBy(x,z,t), Eq.95 lω dz ffi,(*, ) = — ∑k oos{pnx)-Bζ c s{pnx)}eiik^!) . Eq.96 C n H=-CO Third, a relationship is derived between the coefficient Aζ and the electric field amplitude E0 in the corrugated-vane region, using the approximations in Εqs. 11-13, ω By(x,z,ή = iE0et{k*sL-ωήcos (x-d) Εq.97 Ex{x,z,t) = 0, Εq.98 E^z^sE.e'^-^sin (x~d) Εq.99
for Iz-sJ < a/2 and b<x<d. Here, s = 0, ±1, ±2,--- . This gives
which relates Aζ and E
0. Using the expressions for the electric and magnetic fields in Εqs.94-99, the continuity of the axial electric field at x = b, and Εq.100, it is readily shown that the loaded admittances A (ω,k
z), defined in the same manner as the unloaded admittances A
i (ω, k
z ) in Εq.23, are given by
and
Eq. 102 Setting At{ω,k
2) = A
L-(ω,k
z), Eq. 103 one obtains the loaded dispersion relation
Eq. 104 where use has been made of Eqs. 91-93, and the functions k
n = k,, (kz)>Pn
= Pn (
ω>kz)
an(l
an ~
a n (
ω>^z)
are defined in Eqs. 80, 81 and 93, respectively. When I
b = 0 (or a
n = 0 ), the dispersion relation (Eq. 104) reduces to the vacuum dispersion relation (Eq. 91) as expected. In the Compton regime, |α « 1 and one may make the approximation
and express the loaded dispersion relation approximately as
where D
a(ω,k
2) is the vacuum dispersion function defined in Eq. (3.1.27).
Furthermore, if only one term on the right-hand side of Eq. (5.2.2) dominates, say n = m, then one can further approximate the loaded dispersion relation as Da{ω,K)(ω-kmVb)2 = εm(ω,k2), Eq.107 where the coupling parameter εm (ω, k2 ) is defined by
To estimate the linear gain and bandwidth in the Compton regime, let ()
c,k
c) denote an intersection point of D
a{ω,k
2) = 0 Eq.109 and ω-{k
2 +2πm/L)V
b=0 Eq.110 in the ω versus k
z diagram. Expanding D
a{ω,k
z) about the point (ω,k
2) = (ω
c,k
c), i.e.,
with υ being the group velocity, and introducing the rescaled coupling parameter
one can express the loaded dispersion relation in the following simplified form [ω-ω
c-υ
g{k
2-k
c)][ω~k
2V
b-(2πm/L)V
b]
2?=ε
m. Eq.113
Because Eq.113 is cubic in either ω or kz , it can be solved analytically. The frequency shift δω and the detuning parameter ΔΩm is defined as δω =ω-kzVb- (2πm I L)Vb , Eq.114 AΩm = ωc +υg{kz -kc)-kzVb-(2mnl L)Vb . Eq.115
The loaded dispersion relation becomes (δω)2(δω - ΔΩ„, ) = εm. Eq.116 Since εm >0 and (δDa /dω)<0 (see FIG.3), where εm <0. Consequently, Eq.116
yields the maximum temporal growth rate / 3 I ιl/3 Im&ø I max = — /} \ Iε m\ I Eq *■. 117
at ΔΩ
m = 0. Since ε
m ∞ I
b , the scaling relation
holds in the Compton regime. In the Raman regime, one must treat the space-charge term in Eq. 104 carefully. In this case, the Eq. 104 is expressed as n
Eq.119 where Da{ω,kz) is the vacuum dispersion function defined in Eq. 27. Substituting Eq.
93 into Eq. 119, and assuming that only one term on the right-hand side of Eq. 106 dominates, say n = m , D
a(ω,k
2)[(ω -k
mV
b)
2
ε
m(ω,k
z), Eq. 120 where the coupling parameter ε
m(ω,k
z) is defined in Eq. 108, and the space-charge parameter {QC)
2 m is defined by
Typically, the space-charge parameter (QC)2 n(ω,kz) is positive in the regime of interest. To estimate the linear gain and bandwidth in the Raman regime, let {ωc,kc) denote an intersection point of Da(ω,k2)= 0 Eq. 122 and ω - (kz + 2πm / ∑)Vb - (QC)m (ω, kz ) = 0 Eq. 123 in the ω versus kz diagram. Making use of the expansion in Eq. I ll, one can express the loaded dispersion relation in the following simplified form
[ω -ω
c -υ
g(k
z
-k
2V
b - {2πm/L)V
b -
Eq. 124
where ε
m is defined in Eq. 112. Following an earlier analysis , the frequency shift δω and the detuning parameter ΔΩ
m is defined as δω = ω -k
2V
b -(2πm/L)V
b -{QC)
m(ω
c,k
c), Eq. 125 ΔΩ
m = ω
c +υ
g(k
2 -k
c)-k
zV
b -{2πm/L)V
b -{QC)
m {ω
c,k
c). Eq. 126
The loaded dispersion relation becomes δω(δω -ΔΩm)= *" . Eq. 127 ΆQCLΨM
If ε
m /{QC)
m (ω
c,k
c)< 0 , the system is unstable, and the maximum temporal growth rate is given by
at ΔΩ
m = 0. Since ε
m /(QC)
m{ω
c,k
c)∞
the scaling relation \ Ilmδω\ I max ∞ I °
1 4 Eq
L. 129 holds in the Raman regime. The dispersion relation in Eq. 104 can be solved numerically using Newton's method to obtain the linear gain. For a real value of the wavenumber k
z , the temporal linear growth rate ω
i = ϊm{ω) > 0 can be obtained from the complex ω that solves
Eq. 104. On the other hand, for a real value of the angular frequency ω , the spatial linear growth rate k. = -Jm(kz) > 0 can be obtained from the complex wavenumber kz that solves Eq. 104. FIGs. 9A-9B show plots of the temporal linear growth rate ωt = ϊm(ω) of the lowest anti-symmetric mode as a function of the wavenumber kz for α /Z, = 0.8 , b/L = 1.0, d/L = 6.0 , h/L = 0.S , ^ = 0.08 and two cases (a) ωp 2 bL2/c2 = 7.0x l0-5 and (b) ωp 2 bL2 /c2 = 7.0x lQ~7 , where ωp 2 bL2 /c2 = 4πe2σbL/c2me = πFbL/βbIcyb wit Ic = 17 kA. The parameters for the RF circuit in FIGs. 9A-9B are identical to those used in FIGs. 6 and 7. As seen from FIGs. 9A and 9B, each curve has a peak near ω » kzVb , which represents a stronger beam-wave interaction. For the purpose of illustration, the beam parameters in these
examples correspond to backward- wave oscillators not traveling-wave amplifiers. For this RF circuit, a more energetic electron beam is required in order to make an amplifier. The high-current case shown in FIG. 9A indicates a temporal growth rate in the Raman regime (i.e., with an upright bell-shaped curve), whereas the low-current case shown in FIG. 9B implies a temporal growth rate in the Compton regime (i.e., with a tilted bell-shaped curve). FIGs. 10A-10B show plots of the spatial linear growth rate kl = -Im(&z) of the lowest anti-symmetric mode as a function of angular frequency ω for a IL = 0.8 , b/L = l.O , d/L = 6.0 , A/J = 0.8 , #, = 0.08 . FIG. 9A shows the case ωp 2 bL2 /c2 = 7.0x l0-5 and FIG. 9B shows the case ωp 2 bL2/c2 = 7.0 x l0"7 . These parameters are identical to those in FIGs. 8A-8B. FIG. 11 shows the maximum temporal growth ωι max = [lm(<y)]max of the lowest anti-symmetric mode as a function of normalized beam current ωp 2 bL2/c2 = 4πe2σbL/c2me = πIbL/βbIcyb for /Z = 0.8 , b/L = l.0, d/L = 6.0 , h IL = 0.8 , and βb = 0.08. FIG. 12 shows the maximum spatial growth kl>nm = -[lm(^z)]max of the lowest anti-symmetric mode as a function of normalized beam current ωp 2 bL2 Jc2 = 4πe2σbL/c2me = πIbLJ βbIcyb for the same parameters as in FIGs. 10A-
10B. FIGs. 10A-10B and 11 show that both peak temporals and spatial growth rates increase with the beam current. It is important to suppress any potential unwanted modes in a microwave amplifier. This is true for the PCRB TWA. In the PCRB- TWA, two techniques are used to suppress unwanted modes. One technique is use of a single-sided slow-wave structure instead of a double-sided slow-wave structure, which eliminates the symmetric modes in the operating band and higher frequency bands. The other technique is use of photonic crystals. Typically, photonic crystals include periodic metallic structures (e.g., periodic metal rods) or periodic dielectric (e.g., periodic dielectric layers, rods or spheres) or a combination of periodic metallic and dielectric structures. They can be one-, two-, or three-dimensional. When designed properly, a photonic crystal acts as a frequency-selective and/or mode-selective filter, which keeps the desired operating mode in the amplifier,
and at same time, allows other modes, especially unwanted modes, to escape from the amplifier. In other words, the photonic crystal effectively damps the unwanted modes. The effectiveness of photonic crystals in both frequency selection and mode selection were demonstrated in an oscillator operating at high-frequencies and using an oversized cavity with its characteristic size greater than the wavelength, but it remains to be seen in amplifier configurations, especially for transverse size less than the wavelength. As an example, the dispersion characteristics of wave propagation in photonic crystals can be calculated using the latest Photonic Band Gap Structure Simulator (PBGSS) code developed at MIT. Shown in FIGs. 13A-13B is an example of the dispersion characteristics in a two-dimensional dielectric photonic crystal, where FIG. 13 A shows the TM modes and FIG. 13B shows the TE modes, with the rod radius r0 and the lattice period p , as calculated using the real-space finite-difference (RSFD) and Fourier transform (FT) methods. In this example, there is only a partial band gap at the X point. Here, r01 p = 0.2 , £, = 3.0 for the rods, ε2 = 8.9 for the background, and a mesh of 21x21 cells and 13x13 plane waves are used in the PBGSS RSFD and FT calculations, respectively. The detailed concept design of the PCRB TWA for 3G wireless base stations. will focus on the frequency range from 1920 to 1980 MHz, which is used the initial rollout of 3G wireless network. The PCRB TWA is a 200W, 1950 MHz, 3% bandwidth structure. The parameters and design results are summarized in Table land a cross-sectional view of the amplifier beam tunnel is shown in FIG. 14.
Table 1 Circuit Parameters: Structure Single-sided slow-wave structure Axial Period L = 0.478 cm Slot Width a = 0.382 cm Slot Depth d = 4.20 cm Tunnel Height b = 0.478 cm Tunnel Width w = 1.6 cm
Beam Parameters: Current = 0.11 A
Energy -m
eβ
2c
2 = 2290 eV (# = 0.0946)
Displacement A = 0.378 cm 2x Full Height 4x
b = 0.12 cm (— -^ = 63% ) d -h Full Width 4v
6 = 1.2 cm Beam Temperature — k
BT = 0.1 eV [Calculated Coupling Constant
T YYi c - - -*- = 17 kA and I
b = 4eσ
iΛ#c ]
Wiggler Field Parameters: £w = 207 G λw = 2L = 0.956 cm Results of Gain Calculation: Center Frequency = 1950 MHz Full Bandwidth Af If = 3.1% Gain G = 1.6-3.1dB/period(3.4- 6.6dB/cm) Efficiency η = 80% Photonic Crystal: Structure Dielectric square lattice Band Gap Width 100 MHz
Shown in FIG. 15 is a plot of the eigenfrequency versus the phase shift for the two lowest bands in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA, as calculated analytically using Eq. 27 with the RF circuit parameters listed in Table 1. The ribbon electron beam is designed to interact with the lowest band at about 120° phase shift to achieve RF signal amplification. Using the loaded dispersion relation in Eq. 104 and the parameters listed in Table 1, the complex wavenumber kz is calculated using the GAIN code. The results are summarized in FIGs. 16-19. Shown in FIGs. 16 and 17 are plots of the normalized real and imaginary wavenumbers RekzL and |lm&_Z,| versus the normalized frequency ωLIc for the operating mode in the 200 W, 1950 MHz, 3% bandwidth PCRB TWA, respectively. In terms of the normalized imaginary wavenumber |lm/crZJ , the small-signal intensity
gain per axial period is expressed as
G = 101og10 e2| lm = 8.69|lm Λ:zZ| Eq. 130 Period Shown in FIGs. 18 and 19 are plots of the real phase shift RekzL and the gain G versus the frequency / for the operating mode. In the 3% bandwidth from 1920 MHz to 1980 MHz, the gain is between 3.4 dB/cm and 6.6dB/cm, which is adequate. It should be point out that the gain curve can be made flat by optimizing the RF circuit, so that the beam interacts it at a smaller value of the phase shift. This has been demonstrated in parametric design studies but will not be further discussed in this report. As the ribbon electron beam interacts with the RF circuit, unwanted modes may be excited. Such unwanted modes could arise from the second or higher bands in the RF circuit. If not suppressed, they could cause the amplifier to self-oscillate. One promising technique to suppress unwanted modes is use of frequency-selective and mode-selective photonic crystals as described herein. There are various photonic crystals, ranging from one- to three-dimensional. For the purpose of illustration, two- dimensional dielectric square lattices are discussed. FIG. 20 show the dispersion characteristics for 4 lowest-order TM modes in a photonic crystal consisting of a dielectric square lattice with rod radius r0 = 0.9 cm, lattice constant a = 4.9 cm, rod dielectric constant εx = 2.3 , and background dielectric constant^ = 8.9. Labels F , X, and M on the horizontal axis follows the conversion in solid state physics. In this design, there is a narrow band gap of 100 MHz at 2 GHz, but no band gaps at lower or higher frequencies. This 100 MHz band gap will confine the operating mode and simultaneously allow all unwanted modes to transmit through, achieving the single-mode operation of the 200 W, 1950 MHz, 3% bandwidth PCRB TWA. The photonic crystal design can still be optimized with larger values of dielectric constants and smaller lattice constants. Because the PCRPB TWA is scalable to higher frequencies, wider bandwidth, and higher power output, the 1950 MHz PCRB TWA can be redesigned as a power amplifier for high-frequency (3-6 GHz) 3G wireless base stations as well as for future wireless base stations. Although the present invention has been shown and described with respect to
several preferred embodiments thereof, various changes, omissions and additions to the form and detail thereof, may be made therein, without departing from the spirit and scope of the invention. What is claimed is: