Optimization of the S-band side-coupled cavities for proton acceleration

2020-03-19 03:26HaoYunLiXinMiaoWanWeiChenChenHuiShiZhiHuiLi
Nuclear Science and Techniques 2020年3期

Hao-Yun Li· Xin-Miao Wan· Wei Chen· Chen-Hui Shi·Zhi-Hui Li

Abstract The proton beam with energy around 100 MeV has seen wide applications in modern scientific research and in various fields.However,proton sources in China fall short for meeting experimental needs owing to the vast size and expensive traditional proton accelerators.The Institute of Nuclear Science and Technology of Sichuan University proposed to build a 3 GHz side-coupled cavity linac(SCL)for re-accelerating a 26 MeV proton beam extracted from a CS-30 cyclotron to 120 MeV. We carried out investigations into several vital factors of S-band SCL for proton acceleration,such as optimization of SCL cavity geometry,end cell tuning, and bridge coupler design. Results demonstrated that the effective shunt impedance per unit length ranged from 22.5 to 59.8 MΩ/m throughout the acceleration process,and the acceleration gradient changed from 11.5 to 15.7 MV/m when the maximum surface electric field was equivalent to Kilpatrick electric field.We obtained equivalent circuit parameters of the biperiodic structures and applied them to the end cell tuning;results of the theoretical analysis agreed well with the 3D simulation.We designed and optimized a bridge coupler based on the previously obtained biperiodic structure parameters, and the field distribution un-uniformness was <1.5% for a two-tank module. The radio frequency power distribution system of the linac was obtained based on the preliminary beam dynamics design.

Keywords Proton beam · Side-coupled cavity linac ·Accelerating cavity · Biperiodic structure · Bridge coupler

1 Introduction

Since the successful first operation of a particle accelerator at the beginning of the last century, particle accelerator physics and technology have made remarkable advances [1]. Energetic particle beams extracted from accelerators have become an indispensable tool for exploring the basic structure and interactions of matter.Applications of particle accelerators have penetrated into almost all fields of the national economy and contributed significantly to the improvement of people’s living standards [2]. However, there are far fewer traditional proton and ion accelerators than there are electron accelerators[3].Further, there are only a few laboratories able to operate them owing to their vast size and high cost. Therefore,research is ongoing to reduce costs and sizes of proton and ion accelerators.

At the beginning of the 1990s, a research group, motivated by the need of a small hadron therapy device suitable for hospitals under the TERA Foundation,proposed to use S-band acceleration structures to accelerate protons[4].A scheme with the cyclotron as an injector (extraction energy of 30 MeV and 60 MeV) and SCL as the booster was proposed [4–8]. The feasibility of an S-band SCL structure for a 60 MeV proton beam accelerator was confirmed experimentally with an effective acceleration gradient of 28.5 MV/m, which was considerably higher than that obtained in conventional proton acceleration structures[9–11]. This structure paved the way for compact and cheaper proton facilities.

The Institute of Nuclear Science and Technology of Sichuan University is operating a CS-30 cyclotron, which can accelerate protons (extraction energy of 26 MeV),deuterons, and alpha particles. The cyclotron is currently used for studies on radioisotope production and material irradiation [12]. To achieve better utilization of the CS-30 and satisfy the demands of fundamental experimental research, such as the single-event effect and the biological effect of radiation, we proposed to build an S-band SCL booster with a total length of about 10 m to accelerate proton beams to 120 MeV.This paper mainly discussed the optimization design of SCL cavity structures and related work of the RF design.

2 Design and optimization of a single cavity

The SCL belongs to the biperiodic structures composed of the accelerating cavities (ACs) and coupling cavities(CCs)[3].ACs are located on the beam axis,while CCs are located between ACs and off-axis up and down alternately,as shown in Fig. 1. When operating at π/2-mode, the electric field distribution in ACs was just as in the π-mode,while the electromagnetic field was approximately zero in CCs for the π/2 phase shifts compared with that in the adjacent ACs [13]. The high-frequency characteristics of the SCL structures were mainly defined by ACs and had both a higher acceleration gradient as π-mode structures and good stability as π/2-mode structures, which afforded the overall structure the advantages of a high acceleration gradient, good stability, insensitivity to machining errors,superior temperature control, and assembly diversity[14, 15]. For easy fabrication, the SCL usually adopts a graded β structure,i.e.,AC length is kept constant as βgλ/2 in each tank, where βgis called geometry β. For different tanks, βgincreases gradually as the particle energy increases for high acceleration efficiency [3], and it is roughly equal to the average speed of protons through the tank.However,the shape of cavities does not change much with the smoothness of the electric field distribution and convenience for cavity fabrication and post-processes.

Referring to the SCL cavities model developed in the Superfish program [16], we designed the initial model of ACs and CCs, as shown in Fig. 2. Figure 2a shows a simplified diagram of the AC structure projected in the (r,z) plane, with 14 control parameters provided for cavity optimization. Among them, three parameters a, b, and c were used to identify the shape and position of the tuning ring in the tuning procession and would not apply in the initial cavity optimization. The cavity length, L, was determined by the average particle velocity crossing the tank βg.S is the wall thickness and is mainly defined by the requirement of cooling and mechanical stability. Dacis the diameter of the AC and is determined by the resonate frequency. In our optimization, parameters S and Dacwere kept constant in all ACs for the convenience of cavity processing and welding. The maximum surface electric field was another crucial factor as it determined the field level at which the cavity could be operated safely. The maximum surface electric field usually locates at the chamfer of the nose cone region, so R1and R0were the essential parameters in our optimizations. We also had to consider that the linac could accept the beam extracted from the cyclotron,which had a relatively larger emittance than that from linac, so the aperture radius Rbrequired optimization to the largest possible size. Figure 2b shows the schematic diagram of the CC structure projected in the(r, z) plane. Since in the ideal biperiodic structure, the power level in CC was almost zero,even in real operation,it was also much lower than that in ACs. The only necessary procedure was to keep it resonant at 3 GHz by adjusting the gap width gc. Therefore, our main task concentrated on the optimization of ACs, and all work was based on numerical simulation with CST Studio Suite.

The maximum electric field is a function of RF frequency [3, 5]:

When the working frequency was 3 GHz, we obtained the maximum electric field Ekof 46.8 MV/m from Eq (1).The experimental results demonstrated that in the real case,maximum electric field maintainable was also determined by the pulse length,and for short pulse machine,it could be several times larger than Kilpatrick limits Ek. The acceleration gradient was obtained as a function of work conditions and was proportional to the maximum surface electric field.We defined Eacc,Kas the acceleration gradient when the maximum electric field was Ekas in Eq (2):

The Epeakis the maximum surface electric field, E(s, β,t) is the electric field along the cavity axis, and Ekis the Kilpatrick limit. Obviously, the maximum surface electric field is most severe for the smallest cavity. Therefore, our optimization focused on the shortest cavity with length L at 11.5 mm, which corresponded to the CS-30 extraction energy of 26 MeV (β = 0.23). The effective shunt impedance was the ratio of the voltage squared to the power loss in the cavity that described the cavity acceleration efficiency:

From 3D electromagnetic simulation results, it was evident that the maximum electric field was located at R1.Figure 3 shows the relationship of the effective shunt impedance Reffand the acceleration gradient Eacc,Kas a function of parameter R1, while keeping the resonant frequency at 3 GHz by adjusting the cavity diameter Dac.It is shown that when R1changed from 0.1 to 1.9 mm, the Reffdecreased from 23.9 to 21.1 MΩ/m, but the Eacc,Kincreased from 5 to 11.6 MV/m. In the R1adjustment process,the Dacvaried from 74.5 to 73.4 mm.The value of R1= 1.1 mm was finally chosen for a small maximum surface electric field with good acceleration efficiency.

Figure 4 shows the relationship between Rb/L, Reff, and Eacc,K.When Rb/L changed from 0 to 0.35(corresponds Rbfrom 0 to 4 mm),the Reffand Eacc,Kwere both decreasing.When Rb/L was <0.13 (Rb= 1.5 mm), the decreasing trend was relatively flat; however, when Rb/L was >0.13,both Reffand Eacc,Kdecreased more rapidly as Rb/L increased. As a trade-off between large acceptance and good RF properties, Rb= 3.5 mm was chosen as our research value. The corresponding Reffand Eacc,Kwere 22 MΩ/m and 10.8 MV/m, respectively. Compared with those when Rb= 2.5 mm, they decreased by 22.4% and 15.2%,respectively.

As energy increased, the cavity length should also be increased. We adjusted the gap width g rather than cavity diameter Dacto keep the eigenfrequency at 3 GHz. Figure 5 shows the trend of Reffand Eacc,Kfor the entire accelerating range (β value ranged from 0.23 to 0.46,corresponding energy changed from 26 to 120 MeV). Reffranged from 22.5 to 59.8 MΩ/m and Eacc,Kincreased from 11 MV/m at 26 MeV to 14.7 MV/m at 120 MeV.The final cavity geometry parameters are listed in Table 1.

Table 1 Some critical parameters of the ACs

3 Equivalent circuit parameters and end cell tuning

In SCL, ACs and CCs were coupled through coupling slots and formed a biperiodic structure with the equivalent circuit shown in Fig. 6.Subscript A(C)represents AC(CC).The circuit was characterized by the followed circuit parameters:LAand LCwere the circuit inductance,CAand CCwere circuit capacitances, and RAand RCrepresented circuit resistances. The nearest-neighbor coupling was k,while the second nearest-neighbor couplings were kAand kCfor AC–AC and CC–CC,respectively.According to the cavity equivalent circuit model, we wrote the resonance equations of the infinite biperiodic structure with non-lossy[17–19]:

We solved the homogeneous equation by assuming that:

where φ = πq/2N is the phase shift between adjacent cavities; N is the number of ACs, q = 0, 1, 2, …; and 2 N corresponds to different eigenfrequencies of different eigenmodes. Substituting Eq (6) into Eq (5), the following relationship is obtained:

Eliminating A and B, the dispersion relation is found[18–20]:

Five periodic parameters k,ωA,ωC,kA,and kCin Eq (8)were critical for studying the properties of the biperiodic structure. Therefore, according to Eq (8), five different modes corresponding to five different phase shifts(φ1~φ5)were needed to calculate these five parameters.The shortest period structures in our research model of CST simulations were two complete ACs, one complete CC,and two half CCs.Four eigenmodes were obtained for CST calculation when the periodic boundary conditions were applied in the longitudinal direction with phase shift at γ for each per period:

where m is the number of cavities in the model and n ranges from 0 to m - 1. In our case, m = 4. When n equals to 0,1,2,and 3,the corresponding φ1~φ4are γ/4,π/2 + γ/4,π + γ/4,and 3π/2 + γ/4,respectively.When we limited the phase shift within 0 to π,the phase shifts of the four modes were γ/4,π/2 + γ/4,π-γ/4,and π/2-γ/4.We could easily to get parts 1 and 2 in Fig. 7 corresponding to phases φ1and φ3as the 0- and π-mode electric field distributions were calculated when setting γ = 0.Furthermore,according to the symmetry of the dispersion curve, when setting the ranges of γ/4 from 0 to π/2, we could achieve a full dispersion curve as seen in Fig. 7. It was easy to understand parts 3 and 4 corresponding to φ2and φ4,respectively.

Based on the phase relation deduced above, five different phase shifts between adjacent cavities could be calculated by setting two different phases in the shortest periodic structure. First, we calculated φ1= 0, φ2= π/2(the electric field distribution existed in CC or AC),φ3= π/2 (the electric field distribution existed in AC or CC), and φ4= π with corresponding frequencies ω1, ω2,ω3, and ω4by setting the structure period phase to 0.Second, we set the periodic structure phase shift to π and calculated φ5= π/4, φ6= 3π/4 with corresponding frequencies of ω5and ω6. Finally, we applied phases φ1= 0,φ2= φ3= π/2, φ4= π, φ5= π/4, φ6= 3π/4, and ω1,ω2= ω3, ω4, ω5, and ω6to Eq. (8), and it was easy to obtain five circuit parameters.

As shown in Fig. 7, the ‘‘stopband’’ would appear for the untuned biperiodic structure, and this situation broke the stability at π/2-mode operation. By adjusting the gap capacitance of both AC and CC (g and gc), the stopband was eliminated. Nine situations would appear during the stopband adjustment: ω2could have been >, <, or = 3 GHz,and so could ω3.ω2and ω3were independent of each other during the adjustment process. We only needed to find the frequency corresponding to cavity(AC or CC)and adjusted the cavity gap width(g or gc)to meet a frequency at 3 GHz. Usually, the difference between ω2and ω3would be <1 MHz to maintain stability under π/2-mode operation [18].

The illustration in Fig. 8 is the shortest period structure in our research. In order to more intuitively show the change of the coupling coefficient k with the size of the coupling slot and ensure the Dacand the relative position between ACs and CCs kept constant for all tanks, we applied the cube coupling slot for simulation. As shown in Fig. 8, coupling slots are established between ACs and CCs. The coupling coefficient could be changed by adjusting the slot parameter L.Further,the graph shows the relationship among parameters L,k,ωa,and ωc.The range of k was between 0.0034 and 0.097 when L changed from 7.5 mm to 48 mm, and the effective shunt impedance Reffdecreased from 31 to 22 MΩ/m as k increased. The coupling coefficient k was proportional to the slope of the dispersion curve at the π/2-mode. The larger the slope, the larger the frequency difference was between adjacent modes and the more stable was the systems. In the SCL structure,when the ACs length was the same as our model,the Reffwas usually >25 MΩ/m. L = 32 mm was chosen with a corresponding k at 0.04.At this time,Reffwas about 29.6 MΩ/m; five parameters are shown in Table 2.

In the practical application, two end cells broke the original biperiodic properties[17],as shown by the dashed line in Fig. 10. This situation led to the electric field distributes un-uniformness. We needed to choose the right frequency for the end cell to compensate for field distribution caused by the loss of periodic properties.

As shown in Fig. 9, with a complete AC (No. 0 cavity)as the end cell,we filled the missing cavity(No.-1 CC and No.-2 AC)by the dotted line in the biperiodic structure as there was only the secondary adjacent coupling between the AC and the CC in our calculation. The boundarycondition was antisymmetric as X-1-n= -X-1+nfor the electromagnetic level at zero in No. -1 CC. The end cell circuit equation is as follows [21–23]:

Table 2 Five basic parameters of the biperiodic structure

Taking the boundary condition X0= -X-2into Eq. (10),

The end cell should be satisfied with

We then obtained ωe= 3006.3 MHz by substituting the circuit parameters in Table 2 into Eq. (13).That is,the end cell needed to satisfy their frequency at ω′A= 3006.3 MHz in the biperiodic structure with other conditions, such as coupling slot and CC cavity, remaining constant. By optimizing AC control parameters a,b,and c link to the tuning ring, we achieved the final result of ω′A= ωe=-3006.3 MHz, which corresponded to a single cavity frequency equal to 3051.5 MHz. We employed six ACs models with end cells at both ends in our investigation to verify the final tuning result. The electric field distribution shows a solid line in Fig. 10, and compared with the field distribution without tuning (dashed line in Fig. 10), the field un-uniformness was <1%.

4 Design of the bridge coupler

A bridge coupler(BC)was applied to connect two tanks in a module and maintain the possibility of having enough space to install a quadrupole magnet; these are usually comprised of three parts: two half cavities (BCC) and a complete central cavity (BCA) [9, 24–26]. Each BCC combines with a half CC from a single side of the tank to form a new complete cavity; meanwhile, the BCA couples with the waveguide to provide an RF source for the accelerator. Figure 11a is a cross section of the BC,wherein each BCC is connected to the BCA through two coupling slots. Usually, two pairs of coupling slots are adjusted 90° to each other to reduce the second nearestneighbor coupling. As shown in Fig. 11b, the coupling coefficient between BCA and BCC can be adjusted by changing the angle or width of the slot.

In order to keep the uniformness distribution of the electric field in the two tanks,BC periodic parameters must satisfy the periodic conditions when only AC and CC are existing. In this paper, the BCA played the role of AC in the periodic structure but did not provide the acceleration electric field for the particles; meanwhile, the BCC played a role as CC. The steps of the BC design could be divided into the following steps: (1) design of the BCA and BCC,(2)connection of three elements by coupling slots, and (3)adjustment of the coupling coefficient and calculate biperiodic circuit parameters. At the same time, the following two conditions were required: (1) selection of the suitable coupling slot to reduce the energy reflection into the waveguide and (2) ensuring that the resonance frequencies of BCA and BCC were equal to ωAand ωC,respectively [17].

We used the same biperiodic structure model and parameter calculation method as AC-CC. Five periodic parameters were ωBCA, ωBCC, kBCA–BCC, kBCA–BCA, and kBCC–BCC, where ωBCAand ωBCCwere the frequencies of BCA and BCC in the periodic structure, respectively;kBCA–BCCwas the coupling coefficient between BCA and BCC, and the second nearest-neighbor coupling was kBCA–BCAand kBCC–BCC.By adjusting the coupling slot and cavity capacitance, we obtained the basic parameters of BC, as shown in Table 3.

Note that we only discuss the first three parameters in Table 3 as the second nearest-neighbor coupling for BCC–BCC and BCA–BCA would not exist in one BC. The kBCA–BCCvalue usually chosen for reducing store energy in BAC is 0.05[17].However,in cavity simulation processes,it is impossible to tune(ωBCA,ωBCC)to be precisely equalto (ωA, ωC) for the different shapes and the coupling method between BCA-BCC and AC-CC. As shown in Fig. 12,a BC was applied to connect the two tanks,where each tank contained five ACs with end cells at both ends to verify the final electric field distribution.

Table 3 Five basic parameters of bridge coupler

The peak electric field unevenness was controlled below 1.5% after end cell and BC were added in the biperiodic structure.

5 Conclusion

Single-cell cavity optimization was performed for the βgrange of 0.238–0.467. The effective shunt impedance ranged from 22.5 to 59.8 MΩ/m; the corresponding acceleration gradient was 11 MV/m at 26 MeV and increased to 14.7 MV/m at 120 MeV when the maximum surface electric field was equal to the Kilpatrick field.Based on the beam dynamic simulation, it is feasible that by applying a permanent magnetic quadrupole with focusing and defocusing system focusing structure and acceleration gradients of 15 MV/m in all cavities, a total number of 24 17-cell tanks would be needed to accelerate a proton beam to 120 MeV. Four tanks were grouped into a module and connected by BCs and powered by a klystron to improve the klystron utilization. The peak power of these six modules ranged from 5.05 to 6.36 MW.The total length of the SCL structure was approximately 8.7 m. The overall layout is shown in Fig. 13.