Introduction
Electron linear accelerators are commonly used in both industry and medicine, i.e., for radiotherapy, nondestructive testing, industrial irradiation, and sterilizing [1]. Currently, most electron linear accelerators operate in the S band. Their large weight and volume cause various difficulties during testing and installation, such as a lack of space in the operating room for intraoperative irradiation. Linear accelerators operating at C- and X bands have been proposed to overcome these problems [2-4].
In recent years, X-band electron linear accelerators have attracted considerable attention worldwide. Compared with S- and C-band linear accelerators, X-band linear accelerators have advantages such as a compact structure, light weight, and high accelerating gradient. In particular, the Accelerator Laboratory at Tsinghua University has conducted research in this field [5]. Over the past ten years, S- and C-band standing-wave (SW) linear accelerators driven by domestic magnetrons have been successfully studied [6-9].
For radiotherapy applications, a narrow energy spread is important because it can significantly affect the dosage rate. The consumption of radiofrequency (RF) power by the linear accelerator is expressed as [10]:
With a 2.4 MW input RF power source, we designed a compact X-band SW linear accelerator at Tsinghua University for medical applications, including X-ray and intraoperative electron beam radiotherapies. The design goal was to achieve a 100 mA pulse current with the kinetic energy of 9.0 MeV. The length of the accelerating structure is less than 600 mm. In addition, the SW biperiodic on-axis coupled structure was selected and operated at π/2 mode and 9.3 GHz. RF focusing was introduced using noses in the accelerating cells, and an external solenoid was unnecessary. A narrow energy spread was achieved using a thermal DC gun and specifically optimized bunching cells. The basic components of a linear accelerator are illustrated in Fig. 1, including an RF power source, circulator, accelerating structure, thermal DC gun, and titanium window. In addition, Table 1 presents some critical parameters of the accelerator.
Parameter | Value |
---|---|
Number of accelerating cavities | 29 |
Length (m) | 0.59 |
Input RF power (MW) | 2.4 |
Pulse beam current (mA) | 100 |
Capture ratio | 32% |
RMS beam spot radius (mm) | 0.5 |
Energy dispersion (keV) | 250 |
Average accelerating gradient (MeV/m) | 15 |
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F001.jpg)
This low- and high-power accelerator was designed, fabricated, and tested in our laboratory. The frequency of each cell, distribution of the on-axis electric field, and coupling coefficient were measured during the cold test. The corresponding results are in good agreement with those of our previous simulation. A high-power RF test was conducted to measure the performance of the accelerator, including the output beam current, capture ratio, energy spectrum, size of the beam spot, and breakdown rate (BDR), among other useful data.
This paper presents the RF design, simulation analysis, fabrication, and test results of the proposed SW linear accelerator. The remainder of the paper is organized as follows. In section 2, the RF design and analysis of the X-band accelerating structure are described. The results of the calculations for the electronic gun simulation and beam dynamics are also presented. Section 3 introduces the fabrication, tuning, and low-power RF testing of the accelerating structure. Section 4 details the setup and results of the high-power RF experiment. Finally, the conclusions and outlook are provided in section 5.
RF design and analysis
A high accelerating gradient is necessary to improve the accelerating efficiency and increase the output beam energy without excessive length. The longitudinal shunt impedance of the accelerator is expected to be as large as possible within a reasonable range of the RF BDR. Accordingly, in this study, we aimed to determine a group of appropriate structural parameters to satisfy these requirements simultaneously. For the design process, the cavity geometry was optimized using SUPERFISH, the beam dynamics were investigated using ASTRA, and full-scale RF simulations were performed using CST codes [11-12].
Cavity geometry optimization
Geometric optimization of normal cells is essential for achieving a high accelerating efficiency. In a linear accelerator, the kinetic energy of an electron gain can be expressed as [13]:
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F002.jpg)
The π/2 mode was selected as the working mode of the electron linear accelerator for maximal adjacent mode separation and frequency robustness during operation. According to the biperiodic structure theory, to maximize the longitudinal shunt impedance, the coupling cell should have a small volume and a very weak internal electromagnetic field. In our model, the coupling cells are only 1.5 mm in length; therefore, their effective shunt impedance and quality factor can be ignored.
The RF characteristics of the normal cells are presented in Table 2: fa indicates the frequency of the accelerating cells, fc represents the frequency of the coupling cells, Zeff is the effective longitudinal shunt impedance, and Q0 is the unloaded quality factor.
Parameter | fa(GHz) | fc (GHz) | Zeff (MΩ/m) | Q0 |
---|---|---|---|---|
Value | 9.3 | 9.3 | 165 | 8900 |
Design and simulation of cavity chain
As the nose weakens the electric coupling through the beam hole, magnetic coupling was adopted instead. We implemented magnetic coupling slots between the accelerating and coupling cells as portrayed in Fig. 3.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F003.jpg)
First, the coupling coefficient between cavities was determined. To avoid the subadjacent coupling effect, the adjacent magnetic coupling slots were arranged orthogonally, as shown in Fig. 3. The dispersion equation for a coupled cavity chain can be expressed as [13]:
Design of external coupler and waveguide window
To feed RF power into the accelerator, an external coupler comprising a unique coupler cell and an external coupling hole was connected to the rectangular waveguide. A vacuum model of the accelerator with a coupler was built into the CST, as shown in Fig. 4. The results, including the on-axis electric field distribution and S11 parameter, are displayed in Fig. 5.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F004.jpg)
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F005.jpg)
The external coupling coefficient is given by [13]:
Waveguide windows such as klystrons, magnetrons, and accelerators are widely incorporated in electric vacuum devices. The functions of the waveguide window are to translate RF power with low loss, and provide isolation between the vacuum within the accelerator and the SF6 inside the waveguide. For this accelerator, we optimized and fabricated a waveguide window based on a 2.4 mm-thick quartz dielectric slice, as shown in Fig. 6.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F006.jpg)
According to classical transmission line theory, when RF breakdown occurs in the accelerator, the terminal is mismatched, and the incident wave is fully reflected. Currently, transmission waveguides contain SW components. The RF window should be positioned far from the antinode of the SW and near the wave node, as indicated in Fig. 7, to reduce the probability of dielectric breakdown on the surface of the quartz window.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F007.jpg)
Accelerating electric field buildup time
The field buildup time of an RF structure is closely related to its loaded quality factor. Unloaded quality factor Q0 of the cavity chain was obtained via CST simulation. The loaded quality factor can be directly calculated using Q0 and external coupling coefficient β.
Based on the results of the RF analysis, QL was calculated as approximately 3055. In the field construction process inside the accelerator, basic parameters, such as the electromagnetic field and input impedance within the structure, change exponentially. In particular, the amplitude of the electric field in the cavity increases exponentially with time, as:
Further, at the end of the transient time
We can obtain that
Simulation of thermal DC gun and beam dynamics
Thermal DC gun
We adopted a thermal cathode DC gun as the electron source, constituting a cathode, heating filament, focusing electrode, anode, insulating ceramic shell, and connector to the accelerator. The emission model was simulated using the CST particle module. The emission current at 12.5 kV high voltage is 320 mA, and the perveance is 0.23 μP. The simulation results and basic parameters are shown in Fig. 8, while the basic parameters of the DC gun are listed in Table 3.
Parameter | Ve (kV) | Ie (mA) | Pe (μP) | Rw (mm) | Lw (mm) | Rh (mm) |
---|---|---|---|---|---|---|
Value | 12.5 | 320 | 0.23 | 0.316 | 13.05 | 1.5 |
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F008.jpg)
Beam dynamics study
Our design process involved a beam dynamics study using the ASTRA codes. The apertures were carefully set based on the nose position and the radius. We matched the actual accelerating structure and achieved an accurate beam loss calculation. The pulse beam current of 100 mA and capture ratio of 32% were achieved. The results of the beam dynamics calculations are displayed in Fig. 9.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F009.jpg)
The accelerating structure generates weak RF focusing/defocusing on the electron beams during their transit. Depending on the electric field distribution of the symmetric accelerating cell, the electrons were alternately focused and defocused. However, their momentum is greater when the electrons are focused than when they are defocused, leading to an overall focusing effect.
Along the light-speed section, the energy of the accelerating electrons gradually increases, and the space charge force decreases. Benefiting from the RF focusing effect, application of an additional solenoid is unnecessary in our design for reducing the additional space cost and inconvenience. Fig. 10 shows the beam spot and the transverse and longitudinal phase–space distributions at the exit of the accelerator.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F010.jpg)
Low-power RF test
Thirty-three pairs of copper disks were manufactured in total, each with four cooling water channels. Owing to their thinness, all copper disks were fabricated without tuning of the holes or pins. Each disk was processed using a smaller cavity radius. After each cold test, unsatisfactory copper disks were reprocessed using a high-precision milling machine. In addition, the structure was reprocessed for cavities with a large deviation away from 9.3 GHz to achieve finetuning. All RF cavities were tuned into an appropriate range around the operating frequency of 9300±5 MHz. An RF cold test was performed at each fabrication step. The test results for the last machining round agreed well with the simulation results.
A few disks and the low-power RF test setup are shown in Fig. 11. The cavity chain was assembled using a fixture, the on-axis field distribution was measured via the bead-pulling method, and the S11 data were sampled [15,16,17,18]. The results were compared with our previous simulation results, as can be seen in Fig. 12. The measured on-axis field distributions are compared in Fig. 5 (a). The cold test results show a good agreement with the simulation and actual measurement results.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F011.jpg)
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F012.jpg)
High-power RF experiment
High-power experiment platform and conditioning
The assembled cavity chain, waveguide coupler, and quartz waveguide window were brazed with a silver-based alloy. Subsequently, a cooling water pipe, thermal DC gun, and titanium window were welded to the accelerator structure using an argon arc. After vacuum-leak detection and cathode activation of the thermal DC gun, the accelerator structure was complete, as shown in Fig. 13.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F013.jpg)
Using the waveguide and titanium windows, the vacuum inside the accelerator tube was sealed and isolated from the SF6 gas in the high-power waveguide transmission system. A CF63 flange was installed to connect to the magnetic spectrometer and measure the energy spectrum of the output beam. A titanium getter pump was welded onto the waveguide to maintain low vacuum in the accelerating structure. The vacuum degree of the accelerating structure should remain lower than 10-6 Pa. With fixed performance of the thermal DC gun, the corresponding applied high voltage can be adjusted to change the emission current: The emission and pulse beam currents increased as the voltage increased. The mismatch between the initial kinetic energy and accelerating field caused a drop in the capture ratio at excessive voltages (>13 kV). Therefore, we adjusted the high voltage in the range of 12–13 kV for the high-power experiment.
Subsequently, the RF characteristics of the completed accelerator were determined. For this, the reflection coefficient, SW ratio, quality factor, and external coupling coefficient were calculated. An R&S ZVA60 Vector Network Analyzer was employed as the measuring instrument. No significant frequency shift was observed after brazing. The simulation results were slightly higher than the external coupling coefficient of 1.8. The reason for this may be an error in the waveguide port calibration before the cold test, as well as the unsatisfactory electrical contact between the waveguide window and measurement port.
Based on a study on the RF breakdown phenomenon of a high-gradient accelerating structure, we know that the RF BDR is closely related to the accelerating gradient and RF pulse width. The relationship between the BDR, accelerating gradient, and pulse width can be expressed as [19-21]:
In our high-power experiment, the klystron pulse was long, up to 10 μs, and the ignition rate was high. Moreover, the pulse must be conditioned in advance to reduce the BDR and improve the stability during high-power experiments. The RF power fed into the accelerator was gradually increased. To avoid irreversible damage to the cavity wall, the accelerating structure was fully conditioned to control the BDR below 10-4 /pulse, with the input RF power of 2.4 MW [22-23].
We built a high-power experimental platform to obtain such results for the accelerator. A compact 2.4 MW X-band high-power multibeam klystron was installed as the RF power source to provide input RF power for the accelerator. The circulator separates the klystron and accelerator, and these components are connected by a waveguide filled with high-pressure SF6 gas. We installed several directional couplers to monitor the incident wave in the accelerator and the wave reflected from the accelerator. To measure the output beam current, a Faraday cup with an integrated circuit was installed at the exit of the accelerator.
Additionally, a Tektronix DPO4034 oscilloscope was employed to sample electrical signals, including the klystron's emitting current and the collector current of the thermal DC gun, the emission current and high voltage, and the output beam current collected by the Faraday cup. Moreover, the electrical signals were denoised to eliminate high-voltage interference from the high-voltage modulator. A Rode & Schwarz NRP-Z81 power meter was incorporated to measure the incident and reflected wave signals of the accelerator. The beam current waveform, energy spectrum, capture coefficient, and output beam spot size were measured and evaluated during the high-power experiments.
Results of high-power experiment
The signals of the output beam current and high voltage applied to the thermal DC gun were sampled using a Tektronix DPO4034 oscilloscope. Accordingly, the best data are displayed in Fig. 14, for which the pulse width is 9.8 μs, the pulse beam current is 100 mA, and the repetition rate is 50 Hz (i.e., the operating duty ratio is 0.5‰). Because of some mismatch between the klystron and modulator, the input RF power is not strictly flat and decreases slowly with pulse time. This imperfect input RF pulse causes a slight decline in the waveform of the output beam current. The ratio between the output beam current and emission current from the thermal DC gun was regarded as an approximate estimate of the capture ratio. Under the conditions in Fig. 14, the capture ratio of (31.2±0.5)% was determined, which is just below the predicted value of 32%, possibly owing to processing and measurement errors.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F014.jpg)
A magnetic spectrometer was installed at the exit of the linear accelerator, as depicted in Fig. 15, to measure the beam energy spread. Further, to reduce the beam spot size and improve the energy resolution, a slit with the length of 38 mm and width of 0.2 mm was inserted at the entry of the magnetic spectrometer during the test. The size of the slit significantly affected the relative brightness of the screen, thereby affecting the energy resolution. With varying magnetic field, electrons with different kinetic energies exhibit unique drifting trajectories. The beam energy spread was evaluated based on the position and brightness of the YAG screen. The energy spectra of the output beams corresponding to different input powers are presented in Fig. 16. For the 2.4 MW input power, the kinetic energy of the output beam was 9.04 MeV and the FWHM of the energy spectrum is approximately 3.5%.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F015.jpg)
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F016.jpg)
A mirror was placed at the angle of 45° to the beamline, and the light was reflected from the YAG screen. A CCD camera with a focusing lens was used to capture images. The CCD camera has a resolution of approximately 50 m/px. The spot size of the output beam was estimated using an image on a YAG screen. The brightness distribution of the images can be used to determine the RMS diameter of the beam because the brightness of the optical transition radiation is proportional to the number of incoming particles. For the input RF power of 2.4 MW and pulse beam current of 100 mA, the beam spot image was captured by the CCD camera, and the calculated RMS radius of the RMS beam spot was approximately 0.5 mm, as depicted in Fig. 17.
-202307/1001-8042-34-07-012/alternativeImage/1001-8042-34-07-012-F017.jpg)
Conclusion
An X-band 9 MeV SW linear accelerator powered with a 2.4 MW klystron was developed at Tsinghua University. Its advantages include RF focusing without a solenoid, low-energy spread with a thermal DC gun, and no tuning pins.
SUPERFISH, CST, and ASTRA codes can be employed to facilitate predesign of the accelerator RF parameters, cavity geometry optimization, cavity chain RF simulation, and beam dynamics study. In addition, we conducted a low-power RF test after fabricating and tuning the copper disks. The results of the cold test and simulation are in good agreement.
A high-power RF experiment was further performed using an accelerator. During the high-power test, the output beam current, capture ratio, beam energy spectrum, and beam spot size were measured at the accelerator exit. With 2.4 MW input RF power, we obtained the pulse beam current of 100 mA and pulse width of approximately 9.8 μs (under 0.1‰ duty factor). The output RMS spot radius was 0.5 mm and the output kinetic energy was 9.04 MeV, with the spectral FWHM of 3.5%. The results of the high-power RF experiment demonstrate the good performance of this accelerator. In future, we will examine its performance with a high duty factor of 1. Further explorations will also be conducted regarding the proposed linear accelerator for higher-order modes.
Present status of the accelerator industry in Asia
, inMedical application of C-band accelerator technologies
, in Proceedings of the 19th International Linear Accelerator Conference,Low-energy linacs and their applications in Tsinghua university
, in Proceedings of the 23rd International Linear Accelerator Conference,Development of a new concept automatic frequency controller for an ultra-small C-band linear accelerator guide
. Med. Phys. 34, 1797 (2007). doi: 10.1118/1.2752581Fabrication, tuning, and high-gradient testing of an X-band traveling-wave accelerating structure for VIGAS
. Nucl. Sci. Tech. 33, 102 (2022). doi: 10.1007/s41365-022-01086-yDesign of a C-band 6 MeV standing-wave linear accelerator structure
, in Proceedings of the 2nd International Particle Accelerator Conference,Fabrication and high power RF test of a C-band 6 MeV standing-wave linear accelerating structure
, in Proceedings of the 3rd International Particle Accelerator Conference,High power test of a C-band 6 MeV standing-wave linear accelerator
, in Proceedings of the 4th International Particle Accelerator Conference,Development of a C-band 6 MeV standing-wave linear accelerator
, Phys. Rev. Accel. Beams. 16, 090102 (2013). doi: 10.1103/PhysRevSTAB.16.090102Monte Carlo modelling of electron beams from medical accelerators
. Phys. Med. Biol. 44, 157 (1999). doi: 10.1088/0031-9155/44/12/201Development of a 7-cell S-band standing-wave RF-deflecting cavity for Tsinghua Thomson scattering X-ray source
. Nucl. Sci. Tech. 32, 36 (2021). doi: 10.1007/s41365-021-00871-5Design, fabrication, and testing of low-group-velocity S-band traveling-wave accelerating structure
. Nucl. Sci. Tech. 33, 147 (2022). doi: 10.1007/s41365-022-01124-9Tuning of CLIC accelerating structure prototypes at CERN
, in Proceedings of the 28rd International Linear Accelerator Conference,Tuning of X-band traveling-wave accelerating structures
. Nucl. Instrum. Meth. A. 704, 14 (2013). doi: 10.1016/j.nima.2012.11.182Development and high-gradient test of a two-half accelerator structure
. Nucl. Sci. Tech. 32, 60 (2021). doi: 10.1007/s41365-021-00895-x.Study of basic breakdown phenomena in high gradient vacuum structures
. in Proceedings of Linear Accelerator Conference,Geometric dependence of radio-frequency breakdown in normal conducting accelerating structures, Appl
. Phys. Lett. 97, 171501 (2010). doi: 10.1063/1.3505339New local field quantity describing the high gradient limit of accelerating structures
. Phys. Rev. Spec. Top-Ac. 12, 102001 (2009). doi: 10.1103/PhysRevSTAB.12.102001Comparison of the conditioning of High Gradient Accelerating Structures
. Phys. Rev. Accel. Beams. 19, 032001 (2016). doi: 10.1103/PhysRevAccelBeams.19.032001Jia-Ru Shi is an editorial board member for Nuclear Science and Techniques and was not involved in the editorial review, or the decision to publish this article. All authors declare that there are no competing interests.