留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

能散对直线感应加速器焦斑分布的影响

王科 代志勇 夏连胜 张篁 李劲 范培亮 杨治勇

王科, 代志勇, 夏连胜, 等. 能散对直线感应加速器焦斑分布的影响[J]. 强激光与粒子束, 2022, 34: 094002. doi: 10.11884/HPLPB202234.210450
引用本文: 王科, 代志勇, 夏连胜, 等. 能散对直线感应加速器焦斑分布的影响[J]. 强激光与粒子束, 2022, 34: 094002. doi: 10.11884/HPLPB202234.210450
Wang Ke, Dai Zhiyong, Xia Liansheng, et al. Influence of energy spread on the transverse profile of the focused electron beam[J]. High Power Laser and Particle Beams, 2022, 34: 094002. doi: 10.11884/HPLPB202234.210450
Citation: Wang Ke, Dai Zhiyong, Xia Liansheng, et al. Influence of energy spread on the transverse profile of the focused electron beam[J]. High Power Laser and Particle Beams, 2022, 34: 094002. doi: 10.11884/HPLPB202234.210450

能散对直线感应加速器焦斑分布的影响

doi: 10.11884/HPLPB202234.210450
基金项目: 国家自然科学基金项目(11875242, 11805191)
详细信息
    作者简介:

    王 科,caepacc@sina.com

    通讯作者:

    杨治勇,fantasy13089@163.com

  • 中图分类号: TL501

Influence of energy spread on the transverse profile of the focused electron beam

  • 摘要: 介绍了均方根(RMS)、半高全宽(FWHM)和50%调制传递函数(MTF)等效均匀分布等焦斑尺寸评价方法。针对一台直线感应加速器,通过建立理论模型和开展逐束片PIC模拟研究了长脉冲(约100 ns)电子束的聚焦过程,分析了能量、流强与平顶区差异较大的脉冲上升/下降沿对FWHM和MTF焦斑尺寸的影响。研究结果表明:能散度较大的束流上升/下降沿会导致束流焦斑,尤其是MTF焦斑显著增长,在模拟结果中FWHM尺寸增长约9%,而MTF尺寸增长达到约24%,是目前导致束流底宽偏大的主要因素之一。
  • 在X光闪光照相实验中[1-6],光源的尺寸及分布是影响成像质量的一个主要因素。光源点扩展函数的半高全宽(FWHM)作为一种描述焦斑尺寸的常用参数,多年来一直是表征光源性能的关键指标[7]。而近些年来,随着闪光照相技术的持续发展和对图像分辨率要求的不断提升,仅靠FWHM来描述光源的空间分辨率有其固有的局限,相关成像实验发现,如果光源点扩展函数不同,即使其FWHM相同,图像的分辨情况也会不同,因为束流分布的底宽也会对成像质量造成影响,为了充分考虑这一因素,美国LANL提出利用空间频域50%的调制传递函数(MTF)来评价光源尺寸[8-9],该参数也将成为下一步光源底宽优化的重要评价指标。

    高剂量X光源通常来自高能电子束轰击高Z靶材产生的轫致辐射,光源分布很大程度上取决于击靶电子束的焦斑分布[10],优化电子束的焦斑分布是改善光源分布的重要手段。本文通过研究一台直线感应加速器中电子束的聚焦过程,着重分析束流上升/下降沿能量、流强不一致性等对束流焦斑FWHM与MTF尺寸的影响,为束流焦斑的进一步优化提供理论基础。考虑到直线感应加速器中电子束长度达到约30 m(持续时间约100 ns),全束流模拟需要大量计算资源,往往难以实现,本文将高能(约18 MeV)电子束流分割成约1.1 ns的束片,通过动态网格的PIC代码[11](含束流自身电磁力)模拟逐束片的聚焦过程,然后叠加形成最终束流分布。

    这里以常用的光斑尺寸描述方法来描述电子束打靶焦斑,设电子束的横向分布函数为f(x,y),则其对应的线扩展函数fL(x)=f(x,y)dy,点扩展函数为fs(x)=f(x,0)。FWHM尺寸指其点扩展函数fs的半高全宽,记fs(x)的最大值为fs.max,且有fs(x1)=fs(x2)=0.5fs.max,则束流的FWHM尺寸为|x2x1|,该参数主要反映有效应用区域的尺寸,是实验中常用的评价指标。调制传递函数(MTF)等效均匀分布尺寸描述线扩展函数fL(x)的空间频域特性,设线扩展函数fL(x)的傅里叶变换为F(k),且有F(k0)=0.5max(F(k)),则束流的MTF尺寸为0.705/k0(0.705为均匀圆盘分布情况下k0与半高全宽的乘积),该参数对于分布的底宽较为敏感。

    直线感应加速器中常用的束流分布表征函数有均匀分布(KV)、高斯分布(Gaussian)和本涅特(Bennett)分布,在相同的FWHM时,三种分布的底宽和MTF焦斑依次增大,其分布曲线(线扩展函数)及对应的空间频域曲线如图1所示。在三种分布函数下,FWHM尺寸σFWHM和MTF尺寸σMTF有线性对应关系,如表1所示。均匀布时σMTF/σFWHM=1,高斯分布和本涅特分布的底宽依次增大,σMTF/σFWHM分别为1.6和2.7,可见σMTF/σFWHM值与分布底宽呈正相关,后文的分析中将通过σMTF/σFWHM来判断束流的分布及其底宽特征。

    图  1  均匀分布、高斯分布、本涅特分布曲线(a),及对应的空间频域曲线F(k)(b)
    Figure  1.  The curvatures of Kapchiskij-Vladimirskij (KV) distribution, Gaussian distribution, Bennett distribution (a), and their spatial Fourier expansion (b)
    表  1  不同分布函数下,FWHM尺寸与MTF尺寸的对应关系
    Table  1.  Relationship between the FWHM beam size and the MTF beam size in various distributions
    f(x)σFWHM/mmk0/mm−1σMTF/mm
    KV1.000.7051.00
    Gaussian1.000.4411.60
    Bennett1.000.2582.70
    下载: 导出CSV 
    | 显示表格

    忽略空间电荷效应以及透镜像差,单能电子束在聚焦平面处的均方根尺寸主要受限于束流发射度,近似表示为σf.ϵ=fϵx/σi,其中fσiϵx分别为焦距,束流初始均方根半径及发射度[12]。焦距与电子束的能量正相关即fp0p0为参考电子的能量),一个能量为p的电子与参考电子的相对能量差异为δ=(pp0)/p0,经过聚焦线圈后焦距差异为Δf=fδ,如图2所示。一般而言,电磁透镜的焦距远大于焦斑尺寸,则根据图2易得初始位置为xi、相对能量差异为δ的电子在焦点处的横向位置为ΔxΔfθ=fxiδ,则受能散限制的焦斑尺寸可近似表示为σf.δ=σiσδ,其中σiσδ分别为束流的均方根半径和均方根能散,即xiδ的均方根值。由于能散与发射度的影响相互独立,最终焦斑尺寸可以表示为

    图  2  焦距差异导致的束流横向位置差异示意图
    Figure  2.  The position deviation of electrons induced by focusing length difference
    σf=f2ϵ2x/σ2i+σ2iσ2δ
    (1)

    公式(1)的推导中忽略了束流聚焦前后发射度变化,这在横向非线性力较弱时对于薄束片而言是适用的。以一组实验参数作为计算输入条件:束流的能散分布由加速腔腔压叠加并考虑束流负载效应后获得,中心束片能量17 MeV,如图3虚线所示,流强波形通过B-dot测量获得如图3实线所示,切片发射度ϵx.slice=200 µm,束流聚焦前尺寸σi=10mm,根据螺线管磁场求得中心束片的焦距f=1/e2B2z/(2mcβγ)2dz=209mm,其中me分别为电子的质量和电荷量,Bz为聚焦磁场强度,β为电子经光速归一化的速度,γ为洛伦兹系数,z为轴向坐标变量。将图3所示束流分为91个1.1 ns的切片(从左到右依次编号,中心束片为#46),将其能量差分别代入公式(1)可估算逐束片聚焦后的均方根尺寸。

    图  3  根据实验数据获得的束流的能量和流强分布
    Figure  3.  Current and energy distribution of the electron beam measured in experiment

    假设聚焦后单束片呈高斯分布,将束片分布按权重叠加即可得到聚焦后整个束流的横向分布,束片权重为其归一化电荷量。根据公式(1)计算得到的中心束片(#46)及整个束流(#1-#91束片叠加)的焦斑分布如图4所示,中心束片以实线表示,整个束流以虚线表示,可见整个束流的底宽明显大于中心束片,但是半高宽差异并不明显。

    图  4  根据公式(1)计算得到的中心束片和整个束流聚焦后的横向分布
    Figure  4.  Transverse profiles of the center slice and of the whole beam calculated with formula (1)

    将中心束片聚焦后的MTF和FWHM尺寸分别记为σ0.MTFσ0.FWHM,将#(46-n)号束片至#(46+n)号束片聚焦后叠加分布的MTF和FWHM尺寸分别记为σn.MTFσn.FWHM。根据计算结果可得在焦点位置处σ0.MTF=0.46mmσ0.FWHM=0.29mm,为了分析上升/下降沿叠加后对束流尺寸、分布的影响,进一步计算了n分别取值1~45时叠加分布的MTF/FWHM尺寸相对于中心束片的增长量,将相对增长量表示为Δσn.MTFΔσn.FWHM,其中Δσn.MTF=(σn.MTFσ0.MTF)/σ0.MTFΔσn.FWHM=(σn.FWHMσ0.FWHM)/σ0.FWHM,计算结果如图5所示。可见上升/下降沿束片的叠加使得束流绝对尺寸逐渐增大,但MTF尺寸的增长幅度明显高于FWHM尺寸,即Δσn.MTF/Δσn.FWHM比值随n增大,束流底宽随n的增大逐渐增长。全束流的MTF/FWHM尺寸与中心束片相比分别增长了约18.6%/10.1%,单束片呈高斯分布有σ0.MTF/σ0.FWHM=1.60,而全束片叠加后的MTF尺寸与FWHM尺寸比值增长至1.73,对比表1数据可知,全束流分布介于高斯分布和本涅特分布之间。

    图  5  根据公式(1)计算得到第#(46−n) 至 #(46+n)号束片叠加后的MTF尺寸(右),其中#46为中心束片
    Figure  5.  MTF size of the electron bunch composed with slice #(46−n) to #(46+n) (right), the results are calculated with formula (1)

    上述估算中束片间的唯一差异为能量,可见减小束流能散可以有效的控制束流底宽,减小束流的MTF尺寸。但是由于中间束片的电荷量较大,而上升/下降沿束片电荷量小,束片内的空间电荷力并不一致,因此对束片内分布的影响也不一致,定量评估能散/流强对束流聚焦后尺寸影响需要考虑空间电荷力。

    为了定量评估能量/流强差异较大的上升/下降沿对束流焦斑尺寸的影响,利用PIC代码ASTRA对逐束片的束流动力学行为进行了模拟仿真(含空间电荷力)。初始束片分布设为高斯分布,每个束片含50000个宏粒子,但在叠加获得聚焦后全束流分布时各束片粒子数按其电荷量加入权重因子。

    根据模拟结果,焦斑位置处,中心束片的尺寸为σ0.MTF=0.67mmσ0.FWHM=0.34mm,发射度从200 µm增长至612 µm,尺寸的绝对值与理论分析值存在较大的差异,且σ0.MTF/σ0.FWHM=1.98,中心束片底宽特征介于高斯分布和本涅特分布之间。可见约2 kA的流强下,空间电荷力会改变单束片聚焦后的分布特征,其非线性力已不可忽略。

    为了评估上升/下降沿束片叠加导致的束流尺寸增长,与图4图5的解析结果进行比较,统计了中心束片与全束流聚焦后的分布(如图6所示),以及#(46−n)至#(46+n)号束片叠加分布的尺寸相对于中心束片尺寸的增长(如图7所示)。增长趋势与图5理论计算结果基本一致,全束流的MTF/FWHM尺寸与中心束片相比分别增长了约28.0%、9.8%,全束流的MTF尺寸与FWHM尺寸比值为2.3,更接近于本涅特分布的底宽特征。

    图  6  PIC模拟得到的中心束片和整个束流聚焦后的横向分布
    Figure  6.  Transverse profiles of the center slice and of the whole beam simulated with ASTRA
    图  7  PIC模拟得到的第#(46−n) 至 #(46+n)号束片叠加后的MTF尺寸,其中#46为中心束片
    Figure  7.  MTF size of the electron bunch composed with slice #(46−n) to #(46+n) , the results are calculated with the PIC code ASTRA

    在理论估算(图4)中忽略了空间电荷力,并且假设焦点处束片呈高斯分布,是导致与模拟结果(图5)存在差异的主要原因。对比理论估算结果和模拟结果,虽然聚焦后束流尺寸绝对值差异较大,但是上升/下降沿叠加导致的束流相对尺寸增长趋势却比较一致,即理论模型可用于快速估算上升/下降沿能量差异导致的束流尺寸、底宽相对增长。

    理论模型和模拟结果均表明,中心38 ns电子(n=17)的Δσn.MTFΔσn.FWHM增长量相等,对应的底宽特征与中心束片一致。中心38 ns区域束流强度基本一致,均方根能散为1.7%,则束流平顶区能散小于1.7%时,能散对底宽的影响可以忽略,这需要通过后续的注入器及腔压波形优化来实现。同时注意到从n=37至n=45之间,Δσn.FWHM变化缓慢只增长约1.7%,而Δσn.MTF增长8.5%,对束流底宽的增长贡献较大,该区间对应于上升/下降沿各10 ns的电子,流强小,对最终X射线照射量贡献较低,在后续研究中可考虑使用快速踢束器踢除、利用偏转段的能量—位置差异准直、或者利用不同能量束团的感应加速器振荡差异进行准直。

    本文根据理论计算和数值模拟结果分析了一台直线感应加速器中上升/下降沿电子能量/流强差异对束流焦斑尺寸、底宽的影响。理论模型忽略了空间电荷力,与模拟结果存在一定差异,但能有效的反应上升/下降沿(能量差异较大)叠加导致的束流尺寸增长趋势,可以作为快速估算的手段。通过对模拟结果的分析,从提升注入器/腔压稳定性和踢除上升/下降沿角度,提出了减小束流底宽的具体指标,为下一步实验设计提供了理论指导。而聚焦前束片大小、位置、会聚角等不一致性对焦斑分布的影响与全束线磁场配置、准直误差等有关,需要从头至尾的PIC模拟来分析,将在后续工作中开展。

  • 图  1  均匀分布、高斯分布、本涅特分布曲线(a),及对应的空间频域曲线F(k)(b)

    Figure  1.  The curvatures of Kapchiskij-Vladimirskij (KV) distribution, Gaussian distribution, Bennett distribution (a), and their spatial Fourier expansion (b)

    图  2  焦距差异导致的束流横向位置差异示意图

    Figure  2.  The position deviation of electrons induced by focusing length difference

    图  3  根据实验数据获得的束流的能量和流强分布

    Figure  3.  Current and energy distribution of the electron beam measured in experiment

    图  4  根据公式(1)计算得到的中心束片和整个束流聚焦后的横向分布

    Figure  4.  Transverse profiles of the center slice and of the whole beam calculated with formula (1)

    图  5  根据公式(1)计算得到第#(46−n) 至 #(46+n)号束片叠加后的MTF尺寸(右),其中#46为中心束片

    Figure  5.  MTF size of the electron bunch composed with slice #(46−n) to #(46+n) (right), the results are calculated with formula (1)

    图  6  PIC模拟得到的中心束片和整个束流聚焦后的横向分布

    Figure  6.  Transverse profiles of the center slice and of the whole beam simulated with ASTRA

    图  7  PIC模拟得到的第#(46−n) 至 #(46+n)号束片叠加后的MTF尺寸,其中#46为中心束片

    Figure  7.  MTF size of the electron bunch composed with slice #(46−n) to #(46+n) , the results are calculated with the PIC code ASTRA

    表  1  不同分布函数下,FWHM尺寸与MTF尺寸的对应关系

    Table  1.   Relationship between the FWHM beam size and the MTF beam size in various distributions

    f(x)σFWHM/mmk0/mm−1σMTF/mm
    KV1.000.7051.00
    Gaussian1.000.4411.60
    Bennett1.000.2582.70
    下载: 导出CSV
  • [1] 邓建军, 丁伯南, 王华岑, 等. “神龙一号”直线感应加速器物理设计[J]. 强激光与粒子束, 2003, 15(5):502-504. (Deng Jianjun, Ding Bo’nan, Wang Huacen, et al. Physical design of the Dragon-I linear induction accelerator[J]. High Power Laser and Particle Beams, 2003, 15(5): 502-504

    Deng Jianjun, Ding Bonan, Wang Huacen, et al. Physical design of the Dragon-I Linear induction accelerator[J]. High Power Laser and Particle Beams, 2003, 15(5): 502-504
    [2] Klasky M L. Physics based radiography[R]. LA-CP-18-00890, 2018.
    [3] 石金水, 邓建军, 章林文, 等. 神龙二号加速器及其关键技术[J]. 强激光与粒子束, 2016, 28:010201. (Shi Jinshui, Deng Jianjun, Zhang Linwen, et al. Dragon-II accelerator and its key technology[J]. High Power Laser and Particle Beams, 2016, 28: 010201 doi: 10.11884/HPLPB201628.010201

    Shi Jinshui, Deng Jianjun, Zhang Linwen, et al. Dragon-II accelerator and its key technology[J]. High Power Laser and Particle Beams, 2016, 28: 010201 doi: 10.11884/HPLPB201628.010201
    [4] Ekdahl C. Electron-beam dynamics for an advanced flash-radiography accelerator[J]. IEEE Transactions on Plasma Science, 2015, 43(12): 4123-4129. doi: 10.1109/TPS.2015.2496499
    [5] 施将君. 高能闪光照相引论[M]. 绵阳: 中物院流体物理研究所科协, 1998

    Shi Jiangjun. Introduction of high-energy flash radiography[M]. Mianyang: Scientific Book Series of IFP, 1998
    [6] Klasky M L. A correct flat field model for DARHT[R]. LA-UR-19-29164, 2019.
    [7] 王毅, 李勤, 刘云龙, 等. 用于直线感应加速器光源焦斑测量的双锥厚针孔结构设计[J]. 强激光与粒子束, 2019, 31:065102. (Wang Yi, Li Qin, Liu Yunlong, et al. Structure design of thick pinhole with double cones for spot size measurement of the linear induction accelerator light source[J]. High Power Laser and Particle Beams, 2019, 31: 065102 doi: 10.11884/HPLPB201931.180291

    Wang Yi, Li Qin, Liu Yunlong, et al. Structure design of thick pinhole with double cones for spot size measurement of the linear induction accelerator light source[J]. High Power Laser and Particle Beams, 2019, 31: 065102 doi: 10.11884/HPLPB201931.180291
    [8] Wu Yuanhui, Chen Y J. ENSOLVE: A simulation code for FXR LIA downstream section[C]//Proceedings of the 9th International Particle Accelerator Conference. 2018: 2271-2273.
    [9] Ekdahl C. Characterizing flash-radiography source spots[J]. Journal of the Optical Society of America A, 2011, 28(12): 2501-2509. doi: 10.1364/JOSAA.28.002501
    [10] 王毅, 李勤, 代志勇. 蒙特卡罗模拟分析电子束发射度对照射量空间分布影响[J]. 强激光与粒子束, 2017, 29:065006. (Wang Yi, Li Qin, Dai Zhiyong. Analysis on influence of beam emittance on spatial distribution of exposure using Monte Carlo simulation[J]. High Power Laser and Particle Beams, 2017, 29: 065006 doi: 10.11884/HPLPB201729.170029

    Wang Yi, Li Qin, Dai Zhiyong. Analysis on influence of beam emittance on spatial distribution of exposure using Monte Carlo simulation[J]. High Power Laser and Particle Beams, 2017, 29: 065006 doi: 10.11884/HPLPB201729.170029
    [11] Pöplau G, van Rienen U, Flöettmann K. 3D space charge calculations for bunches in the tracking code Astra[C]//Proceedings of EPAC. 2006: 2203-2205.
    [12] Zhu Jun. Design study for generating sub-femtosecond to femtosecond electron bunches for advanced accelerator development at SINBAD[D]. Hamburg: Universität Hamburg, 2017.
  • 期刊类型引用(1)

    1. 顾余辉. 直线加速器维修维护及质控措施研究分析. 中国设备工程. 2024(07): 48-50 . 百度学术

    其他类型引用(0)

  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  645
  • HTML全文浏览量:  235
  • PDF下载量:  45
  • 被引次数: 1
出版历程
  • 收稿日期:  2021-10-24
  • 修回日期:  2022-04-02
  • 录用日期:  2022-04-19
  • 网络出版日期:  2022-04-21
  • 刊出日期:  2022-06-17

目录

/

返回文章
返回