-
本工作的研究难点在于要满足紧凑型设计的前提下,既要实现粒子能量由入口处的1 MeV/u提高到4 MeV/u;同时又要保证出口束流具有尽可能小的能散以满足同步环的注入要求;还要有较高的束流传输效率,尽可能地减少束流损失来满足辐射防护的要求。因此,关键参数的选择极为重要。
(1) 工作频率的选择
对于常温直线加速器,腔体的品质因数Q反映了在腔体消耗同样功耗条件下,所建立电场的能力,其具体表达式如下:
$$ {{Q}}=\frac{{{\omega}} {{U}}}{{{P}}} \;, $$ (1) 其中:
$\omega =2 {\text{π}} {{f}} $ 为角频率;U为腔体储能;P为功率损耗。腔体单位长度的有效分路阻抗ZT 2有$$ {{ZT}}^2=\frac{\left( {{E}}_0{{T}} \right) ^2{{L}}}{{{P}}} \;, $$ (2) 其中:E0为加速间隙轴向电场;T为渡越时间因子;L为腔体长度。那么(ZT 2)/Q表示了射频结构建立期望电场的能力,且与腔体运行频率成正比。因此,提高频率有利于得到更高的能量转换效率。另外,工作频率提高时,加速单元
$\beta \lambda /2 $ 也会缩短,有助于提高装置紧凑性。但工作频率提高的同时也会带来一系列问题,比如高频率导致腔体体积减小、加工难度提高、RFQ接受度减小[5]等。可见,腔体工作频率是需要综合多种因素后的折中选择。根据已经建成的HI-Linac实际情况,国内现有的实际加工能力和现实中紧凑型要求,并参考国际同类型加速器的工作频率,将工作频率选择为325 MHz。
(2) 腔体
${{K_{\rm{p}}}}$ 值的选择在20世纪50年代,Kilpatrick[6]根据实验数据分析了射频击穿的原因,并定义了高频运行时无击穿现象发生的条件,Kilpatrick的实验结果由Boyd[7]整理成一个公式表示,如下所示:
$$ {{f}}\left( {\rm{MHz}} \right) =1.64{{E}}_{{\rm{k}}}^{2}{\rm{e}}^{-8.5/{{E}}_{{\rm{k}}}} 。$$ (3) 其中:f为高频频率;
$ {{E_{\rm{k}}}} $ 为Kilpatrick Limit,单位是MV/m。根据公式可以得出谐振频率与$ {{E_{\rm{k}}}} $ 的关系曲线,如图1所示。对于我们选定的工作频率325 MHz,计算得到对应的$ {{E_{\rm{k}}}}\approx 17.86 $ MV/m。进一步利用这个方程来估算击穿限值$ {{E_{\rm{s}}}} $ 如下:$$ \begin{array}{l} {{E_{\rm{s}}}}={{K_{\rm{p}}}} {{E_{\rm{k}}}} \end{array} \;,$$ (4) 通常来说,
$ {{K_{\rm{p}}}} $ 值选取在1和2之间,但对于脉冲机器可以适当调高。文献[8]中对腔体$ {{K_{\rm{p}}}} $ 值的极限测试结果,低占空比的脉冲直线加速器的$ {{K_{\rm{p}}}} $ 值选在2.2是较为安全的[8]。 -
在20世纪80年代末,德国法兰克福大学的Ratzinger等[9]提出了一种可同时实现横向与纵向稳定传输的漂移管直线加速器动力学方案,KONUS(KOmbinierte NUll grad Structure)束流动力学。一个完整的KONUS[10]动力学周期主要由三部分组成:(1) 零度同步相位加速段;(2) 横向聚焦段;(3) 纵向聚束匹配段。如图2所示。
图2详细分析了束流在KONUS周期内的完整运动。首先,以0°加速段Section I中的位置a作为起点,束团的相位位于结构定义的同步粒子相位
$ \varphi _{{\rm{s, I}}} $ 的右侧。在Section I内,束团将会不断向负相位滑相,同时由于束团经历的电场小于同步粒子,二者间的能量差距会逐渐减小。随着加速间隙(加速间隙位于两相邻漂移管之间,后文简称gap)数量的增加,束团受到明显的非线性力的作用,导致发射度增大,须在此之前结束Section I。束团在Section I中运动时,高频场对束团存在横向散焦效果,当束流横向尺寸增长到一定程度时,必须要通过b-c段的聚焦元件对包络进行约束。Section II为负同步相位纵向聚束段,该段中束团相位与结构定义的同步粒子相位$ \varphi _{{\rm{s, II}}} $ 重合。经过纵向和横向匹配后,束团进入下一0°加速段Section III,此时结构定义的同步相位为$ \varphi_{{\rm{s, III}}} $ 。束流在一个完整KONUS周期内,按照前后顺序将依次经过主加速段、横向聚焦段、纵向聚束匹配段,便可以实现低发射度增长的高梯度加速。 -
针对KONUS束流动力学理论,法兰克福大学开发了专用的LORASR模拟程序[11]。该程序可以对多种类型结构的漂移管直线加速器进行束流传输模拟,但其主要的模拟对象是基于KONUS的漂移管直线加速器。程序根据同步能量、同步相位、有效电压等一系列输入参数生成相应的加速结构,并基于给定的初始束流分布模拟在加速结构中的传输[12-13]。目前已经成功应用于多个加速器装置的物理设计,如:德国的GSI-HLI[14],欧洲的CERN Linac3[15],中国科学院近代物理研究所的SSC-Linac[16],哈尔滨工业大学的SESRI-Linac[17]和HI-Linac,其可靠性已通过束流实验得到了验证。
-
在本方案中,RFQ出口法兰内表面设定为DTL动力学设计起点。如图3为RFQ出口处的束流相空间分布,其主要参数如表1所列。其中:
$ \varepsilon_{n, {\rm{rms}}} $ 为束流归一化均方根发射度,$ \alpha $ 和$ \beta $ 为束流Twiss参数。由于RFQ出口束流直接注入DTL,该组相空间参数被作为动力学设计的输入条件。束流横向相空间归一化均方根发射度约为0.18 πmm·mrad,X方向为散焦,Y方向为聚焦。表 1 RFQ出口处束流相空间参数
参数 值 $\varepsilon_{x, n, {\rm{rms}}}$ 0.18 πmm·mrad $ \alpha_x $ −1 $ \beta_x $ 0.2 mm/mrad $\varepsilon_{y, n, {\rm{rms}}}$ 0.17 πmm·mrad $ \alpha_y $ 0.98 $ \beta_y $ 0.2 mm/mrad $\varepsilon_{z, n, {\rm{rms}}}$ 31.2 πdeg·keV/u $ \alpha_z $ 0 $ \beta_z $ 1.33 deg/(keV/u) 由于初始束流在X方向为散焦,包络进入DTL即开始增大。为避免束流损失在漂移管内壁,需针对包络变化确定第一段加速段中gap数量。通过模拟第一段加速段gap数量变化对束流包络和传输效率的影响,如图4所示,确定在第一段加速段设置10个gap,其同步相位均为0°。
第一段加速段出口处的X方向相空间分布如图5所示,束团出现了丝化现象,发射度也存在一定程度的增长。
图6给出了gap电场分布。通过分析发现,gap电场的Z向和径向分量在gap中的分布是不均匀的,且越远离轴线其不均匀性越明显,那么束团径向分布就会导致束团内粒子感受到不同强度的Z向和径向电场分量,且横向尺寸越大,电场强度差别也就越大。
根据此前的横向包络分析,束流在X方向为散焦,包络会不断增大,束团内部的粒子也会感受到差别更大的Z向和径向电场分量,从而导致出口处的束团出现丝化。初始束团的非对称分布给束团的横向和纵向运动带来了明显影响,但是束流的传输效率并未受到影响,发射度的增长也还在控制范围,因此可以继续进行后续的物理设计。
在束流横向包络即将增大至接近漂移管内壁时,束流应进入横向聚焦段进行横向匹配。在本文设计中设置了三组四极透镜将束团重新匹配为对称束,三组四极透镜物理参数如表2所列。该段漂移管内半径调整为9 mm,确保束流与内壁间有足够裕量。横向聚焦段结束后束流便将进入负同步相位纵向聚束段,应设置合适的漂移管长度,使束流在进入纵向聚束段第一个gap中心时,高频电场恰好处于设计的聚束相位。该段距离可以通过以下公式计算得到:
表 2 三组内置四极透镜物理参数
磁铁编号 磁场梯度/(T·m−1) 有效长度/mm 孔径/mm Q1 53 80 18 Q2 77 110 18 Q3 80 70 18 $$ {{L}}_{\text {shift }}=\frac{180^{\circ}+\Delta \varphi}{180^{\circ}} \times \frac{\beta \lambda}{2}\;, $$ (5) 其中:
$ \Delta \varphi $ 是漂移段前后gap中心的高频相位差;$ \frac{\beta \lambda }{2} $ 则是加速单元长度;$ \beta $ 为零度同步相位加速段出口处的粒子相对论速度。在横向聚焦段后,通常会设置负同步相位纵向聚束段[参见图2的(c)~(d)]。该段的主要作用是利用负同步相位将纵向发散的束团进行纵向聚束,再注入零度同步相位加速段。聚束单元的个数不宜过多,否则影响加速效率。通过优化,本文确定设置3个负同步相位的聚束单元。在结束了纵向聚束段后,束团将进入一段较长的零度同步相位加速段,可将束团能量加速至4 MeV/u。
各加速间隙的有效加速电压也是优化过程中需重点考虑的[18],应尽可能地使各加速间隙轴上最大电场趋于一致,能量增益达到最大。图7为动力学设计所用到的腔体各gap有效电压及轴上最大电场。
粒子在通过gap时,由于电场随时间发生变化,粒子的能量增益需引入渡越时间因子TTF进行修正。渡越时间因子TTF与粒子能量、电场分布、漂移管结构相关。通过对各加速单元的gap与漂移管长度比以及漂移管内径参数的优化后,全局渡越时间因子TTF如图8所示,整体的TTF维持在0.80~0.84之间,分布较为均匀。
-
基于上述研究,本文完成了CH-DTL的动力学设计。该重离子直线加速器主要参数如表3所列。
表 3 紧凑型重离子直线加速器关键物理参数
装置 参数 数值 工作频率 325 MHz 设计离子 ${} ^{12}{{\rm{C}}} ^{5+}$ 设计流强 200 eμA RFQ 高频结构类型 四翼型(4-vane) 注入能量 8 keV/u 引出能量 1 MeV/u 传输效率 >90% DTL 高频结构类型 交叉杆型(CH) 动力学方案 KONUS 注入能量 1 MeV/u 引出能量 4 MeV/u 引出束流能散 $ \pm $0.3% 传输效率 >90% 束团90%粒子的横向包络演化如图9所示,束流与漂移管内壁间留有充足裕量。DTL段出口处设置了三台四极透镜,可以很容易地将束团匹配至正常尺寸,同时实现了粒子传输效率大于99%的设计目标。
束流在出口处的横向相空间分布如图10所示,束流相空间参数如表4所列。
表 4 动力学设计中出口处束流相空间参数
参数 值 $\varepsilon_{x, n, {\rm{rms}}}$ 0.31 πmm·mrad $ \alpha_x $ 2.02 $ \beta_x $ 0.14 mm/mrad $\varepsilon_{y, n, {\rm{rms}}}$ 0.18 πmm·mrad $ \alpha_y $ 2.42 $ \beta_y $ 0.14 mm/mrad $\varepsilon_{z, n, {\rm{rms}}}$ 68.48 πdeg·keV/u $ \alpha_z $ 1.95 $ \beta_z $ 0.72 deg/(keV/u) 由于gap中非线性高频聚焦和散焦力的存在,在X方向相空间存在畸变,X方向发射度在进入DTL后出现急剧增长,在600 mm后又出现下降;Y方向发射度未发生明显增长;Z方向发射度则呈持续上升趋势。各方向发射度变化如图11所示。
如图12为束流在纵向的相宽、能散包络演化。可以看到,束团的同步相位经过了两次跳相过程,分别发生在横向聚焦段出口和聚束段出口。在加速过程中,束团的相位是一直在由正同步相位滑相至负同步相位。束团出口处相位宽度约为20°,能量分散小于
$\pm 1{\text{%}}$ 。图13为设计中束团中心的演化过程,整体分布较为合理。
根据以上设计,DTL段动力学长度仅为1.3 m,其结构示意图如图14所示,可以将粒子从1 MeV/u加速到4 MeV/u,束流传输效率>99%。相较于同样能量段的HI-Linac中DTL段动力学长度缩短约60%[19],实现了紧凑高效的设计目标。
-
摘要: 基于KONUS动力学设计理论,完成了一台工作在脉冲模式的紧凑型常温交叉杆型漂移管直线加速器(CH-DTL)的束流动力学设计,并对关键问题进行了深入研究。该脉冲直线加速器工作频率选为325 MHz,内含29个加速间隙,动力学长度为1.3 m,可以将
$^{12}{\rm{C}}^{5+}$ 重离子束流能量由1 MeV/u提高至4 MeV/u。对有利于提高装置紧凑性的基本参数,如工作频率、加速电压、Kp值等进行了理论分析和研究。动力学模拟结果显示,粒子在DTL段传输效率大于99%,输出束团相宽约为20°,能量分散小于$\pm 1{\text{%}}$ 。-
关键词:
- DTL /
- 紧凑型强流重离子 /
- KONUS动力学设计
Abstract: In this paper, based on the KONUS dynamics theory, the beam dynamics design of a compact room-temperature cross-bar drift tube linear accelerator(CH-DTL) operating in pulsed mode is completed, and the key issues are carefully investigated. This DTL operates at 325 MHz, contains 29 acceleration gaps, and has a dynamics length of 1.3 m. It can accelerate$^{12}{\rm{C}}^{5+}$ heavy ion beam energy from 1 to 4 MeV/u. The essential parameters that contribute to the compactness of the device, such as operating frequency, voltage, and Kp value, are analyzed and studied theoretically. The simulation results show that the transmission efficiency in the DTL section is greater than 99%, the phase width of the output bunch is about 20°, and the energy spread is less than$\pm 1{\text{%}}$ .-
Key words:
- DTL /
- compact high-intensity heavy ion beam /
- KONUS dynamics design
-
表 1 RFQ出口处束流相空间参数
参数 值 $\varepsilon_{x, n, {\rm{rms}}}$ 0.18 πmm·mrad $ \alpha_x $ −1 $ \beta_x $ 0.2 mm/mrad $\varepsilon_{y, n, {\rm{rms}}}$ 0.17 πmm·mrad $ \alpha_y $ 0.98 $ \beta_y $ 0.2 mm/mrad $\varepsilon_{z, n, {\rm{rms}}}$ 31.2 πdeg·keV/u $ \alpha_z $ 0 $ \beta_z $ 1.33 deg/(keV/u) 表 2 三组内置四极透镜物理参数
磁铁编号 磁场梯度/(T·m−1) 有效长度/mm 孔径/mm Q1 53 80 18 Q2 77 110 18 Q3 80 70 18 表 3 紧凑型重离子直线加速器关键物理参数
装置 参数 数值 工作频率 325 MHz 设计离子 ${} ^{12}{{\rm{C}}} ^{5+}$ 设计流强 200 eμA RFQ 高频结构类型 四翼型(4-vane) 注入能量 8 keV/u 引出能量 1 MeV/u 传输效率 >90% DTL 高频结构类型 交叉杆型(CH) 动力学方案 KONUS 注入能量 1 MeV/u 引出能量 4 MeV/u 引出束流能散 $ \pm $0.3% 传输效率 >90% 表 4 动力学设计中出口处束流相空间参数
参数 值 $\varepsilon_{x, n, {\rm{rms}}}$ 0.31 πmm·mrad $ \alpha_x $ 2.02 $ \beta_x $ 0.14 mm/mrad $\varepsilon_{y, n, {\rm{rms}}}$ 0.18 πmm·mrad $ \alpha_y $ 2.42 $ \beta_y $ 0.14 mm/mrad $\varepsilon_{z, n, {\rm{rms}}}$ 68.48 πdeg·keV/u $ \alpha_z $ 1.95 $ \beta_z $ 0.72 deg/(keV/u) -
[1] KARZMARK C J. Medical Physics, 1984, 11(2): 105. doi: 10.1118/1.595617 [2] LEE Y, HAHN G, KIM E S, et al. Design Study and Multi-particle Tracking Simulation of the IH-DTL with KONUS Beam Dynamics for KHIMA Project[C]//International Particle Accelerator Conference: Number 7 Proc. of International Particle Accelerator Conference (IPAC’16), May 8-13, 2016, Busan, Korea. Switzerland: JACoW Publishing, 2016: 3299. [3] DU H, KONG Q, LI Z, et al. RF Design and Cold Model Measurement of an IH-DTL for HIMM Injector[C]//Linear Accelerator Conference: Number 29 Proc. 29th Linear Accelerator Conference (LINAC’18), September 16-21, 2018, Beijing, China. Switzerland: JACoW Publishing, 2019: 329. [4] BATYGIN Y K. Nucl Instr and Meth A, 2005, 539(3): 455. doi: 10.1016/j.nima.2004.10.029 [5] 王科栋. 紧凑型重离子射频四极场直线加速器的理论设计与研 究[D]. 兰州: 中国科学院近代物理研究所, 2019. WANG Kedong. Theoretical and Experimental Study on Compact Heavyion RFQ Linac, IMP, CAS, 2019. (in Chinese) [6] KILPATRICK W D. Review of Scientific Instruments, 1957, 28(10): 824. doi: 10.1063/1.1715731 [7] BOYDT J. Kilpatrick’s Criterion[Z]. 1982. [8] 李钟汕. 重离子射频四极场直线加速器的物理设计与实验研究 [D]. 兰州: 中国科学院近代物理研究所, 2017. LI Zhongshan. Theoretical Design and Experimental Studies of Heavy Ion Radiofrequency Quadrupole Linacs[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2017. (in Chinese) [9] RATZINGER U. A Low Beta rf Linac-structure of the ih-type with Improved Radial Acceptance[C]//Proceedings of the 1988 Linear Accelerator Conference, October 3-7, 1988, Williamsburg, Virginia, USA. Switzerland: JACoW Publishing, 1989. [10] 乔舰, 殷学军, 原有进, 等. 原子核物理评论, 2017, 34(2): 170. doi: 10.11804/NuclPhysRev.34.02.170 QIAO Jian, YIN Xuejun, YUAN Youjin, et al. Nuclear Physics Review, 2017, 34(2): 170. (in Chinese) doi: 10.11804/NuclPhysRev.34.02.170 [11] TIEDE R, CLEMENTE G, MINAEV S, et al. LORASR Code Development[EB/OL].[2022-03-16].https://cds.cern.ch/record/1078929. [12] TANG R, XING Q, ZHENG S, et al. Nucl Instr and Meth A, 2019, 920: 50. doi: 10.1016/j.nima.2018.10.125 [13] ZHAO L, PANG J, HE X, et al. Nucl Instr and Meth A, 2016, 806: 75. doi: 10.1016/j.nima.2015.10.007 [14] RATZINGER U, HÄHNEL H, TIEDE R, et al. Phys Rev Accel Beams, 2019, 22: 114801. doi: 10.1103/PhysRevAccelBeams.22.114801 [15] ANGERT N, BLEUEL W, GAISER H, et al. The ih Linac of the Cern Lead Injector[Z]. 1994. [16] 杜衡. IH-DTL 直线加速器的设计及实验研究[D]. 兰州: 中国科学院近代物理研究所, 2017. DU Heng. Design and Experiment Study of IH-DTL Linac[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2017. (in Chinese) [17] JIANG H P, LIU J L, HAO H F, et al. Journal of Physics: Conference Series, 2019, 1350(1): 012081. doi: 10.1088/1742-6596/1350/1/012081 [18] LIU J, LI H, LU Y, et al. International Journal of Modern Physics A, 2018, 33(11): 1850062. doi: 10.1142/S0217751X18500628 [19] 杜衡, 原有进, 殷学军, 等. 原子核物理评论, 2018, 35(1): 034. doi: 10.11804/NuclPhysRev.35.01.034 DU Heng, YUAN Youjin, YIN Xuejun, et al. Nuclear Physics Review, 2018, 35(1): 034. doi: 10.11804/NuclPhysRev.35.01.034