宇宙线高能粒子对测试质量充电机制*

韩瑞龙 蔡明辉† 杨涛 许亮亮 夏清 韩建伟

1) (中国科学院国家空间科学中心,北京 100190)

2) (中国科学院大学天文与空间科学学院,北京 100049)

测试质量是空间引力波测量的核心传感器,宇宙线高能粒子能够穿透航天器屏蔽对其造成电荷注入,进而产生库仑力和洛伦兹力噪声对引力波科学探测造成严重影响.本文采用蒙特卡洛仿真方法,探究了不同宇宙线高能粒子对测试质量的充电过程和机制.研究结果表明,在同一能谱下随着截止能量的降低充电速率逐步增大,充电速率变化约为9%;太阳活动极小年时测试质量的充电速率为39.5 +e/s,其中贡献最大的质子占比约为83.16%,太阳活动极大年时测试质量的充电速率约为12.5 +e/s,1989 年最恶劣的太阳高能粒子事件造成测试质量的充电速率约为120700 +e/s;在太阳活动极小年时,银河宇宙线各成分的充电速率取决于各成分的初级粒子在测试质量中的沉积,其中初级粒子贡献占测试质量总充电速率的73%;太阳活动极小年时,质子的充电贡献主要来自能量为0.1—1 GeV 的区间,占比约为65%.研究结果可用于评估测试质量在轨充电规律,为电荷管理的设计和在轨工作提供依据.

1 引言

引力波是爱因斯坦广义相对论中非常重要的一个预言,探测空间引力波不仅是对爱因斯坦广义相对论直接精确的检验,也对物理学和天文学等基础学科未来的发展有着十分重大的作用和意义[1,2].实现空间引力波探测的核心设备是测试质量,测试质量位于航天器的望远镜筒内.航天器在深空环境中航行时,空间环境中的宇宙线高能粒子能够穿透航天器屏蔽入射到测试质量,导致测试质量带电.带电后的测试质量在电极库仑力或者磁场洛伦兹力的干扰下,产生加速度噪声对引力波科学探测造成严重影响.欧洲和美国的研究人员利用蒙特卡洛模拟仿真工具GEANT (geometry and tracking)对LISA(laser interferometer space antenna)测试质量的充电速率进行了仿真模拟评估[3−8],由于其考虑的几何模型、银河宇宙射线(galactic cosmic rays,GCR)通量和截断能量的不同,计算得到太阳活动极小年时GCR 对测试质量的充电速率为11–52 +e/s.利用另一款蒙特卡洛仿真软件FLUKA(Fluktuierende Kaskade)对LISA 简化几何模型的充电模拟也有很大的进展[9],得到太阳活动极小年能量0.1—1000 GeV 宇宙线质子的充电速率为37 +e/s(截断能量为1 keV).国内的研究者主要是针对单航天器激光天文动力学空间计划(astrodynamical space test of relativity using optical devices,ASTROD)利用GEANT4 开展测试质量的充电机制模拟,仿真模拟了太阳活动极小年时高能宇宙射线入射简化等效几何和相对复杂完整几何模型中的充电速率[10,11],得到太阳活动极小年时充电速率为24.2 +e/s 和33.3 +e/s(截断能量 为250 eV).Grimani 等[12]和Wass 等[13]针 对LISA Pathfinder,利用相同的简化模型仿真模拟的充电速率为38.2 +e/s 和43.7 +e/s,在轨充电速率为TM1:22.9 +e/s;TM2:24.5 +e/s[14].虽然国内外的研究者对测试质量的充电速率已经开展了较多的研究,但是对于不同截断能量阈值对充电速率影响、不同宇宙线的充电能力和初、次级粒子充电速率占比、同一粒子不同能量段的充电贡献等方面的研究尚不完善.

本文依据LISA 航天器建立几何模型,仿真采用CREME96 模型计算得到的GCR 能谱,利用粒子传输仿真工具GEANT4 进行模拟.为了揭示宇宙线不同种类、不同能量粒子对测试质量的充电机制与规律.对不同截断长度下测试质量充电影响进行评估.针对测试质量在太阳活动极小年、太阳活动极大年时能量大于100 MeV/n 的质子、氦核(3He 和4He)和C,N,O 等粒子进行充电模拟,对1989 年太阳高能粒子(solar energetic particle,SEP)事件也进行了充电模拟.太阳极小年时对宇宙线各粒子的初、次级粒子充电率贡献以及质子在不同能量段的充电率贡献进行仿真模拟.

2 仿真方法

2.1 几何建模

针对引力波探测中测试质量的充电机制与规律进行仿真模拟,对LISA 的复杂航天器结构模型进行简化[15].LISA 航天器整体模型和惯性传感器模型如图1 所示.

图1 LISA 航天器模型 (a)整体模型;(b)惯性传感器模型Fig.1.LISA spacecraft model:(a) Overall model;(b) inertial sensor model.

本文简化的等效航天器模型如图2 所示,模型包括:测试质量(test mass,TM),为边长46 mm的立方体;钼电极(Mo),为内边长74 mm、外边长为86 mm 的立方体壳层;钛室(Ti),为内半径长75 mm、外半径长80 mm 的球壳层;碳外壳(C)(航天器其它结构如电池、望远镜筒和支架等),为内半径长80 mm、外半径长100 mm 的球壳层;测试质量(TM)与钼电极(Mo)之间和钼电极(Mo)与钛室(Ti)之间的区域均为真空,具体航天器几何尺寸和材料如表1 所列.在本文的模拟中,真空的密度为1.0×10–25g/cm3,设置宇宙线粒子是从一个半径为120 mm 的球面上各项同性均匀地入射到整个航天器模型.

图2 本文航天器模型图 (a) geant4 模型图;(b)模型平面示意图Fig.2.The spacecraft model diagram in this paper:(a) Geant4 model diagram;(b) schematic diagram of the model plane.

表1 航天器几何尺寸和材料构成Table 1.Spacecraft geometric dimensions and material composition.

2.2 辐射环境

引力波探测航天器在类地球轨道绕太阳飞行[16],将受到GCR 和SEP 的轰击,导致测试质量带电.GCR 能谱随太阳周期的变化强度会有所不同,在太阳极小年GCR 的通量最高.在太阳极小年,GCR 能谱组成为:大约90%的质子、8%的氦核(3He 和4He),1%的重核和1%的电子.本文主要考虑太阳极小年通量相对较大的质子、氦核(3He和4He)和C,N,O 等粒子和太阳极大年质子、氦核(3He 和4He)以及对于典型的1989 年SEP 峰值能谱进行充电机制仿真分析.

利用CREME96 模型计算得到太阳活动极大、太阳活动极小年距太阳1AU 处各粒子的能谱如图3 所示,1989 年9 月29 日的SEP 峰值能谱通量如图4 所示.由于只有大于100 MeV/n 的粒子才能穿透航天器测试质量外围的屏蔽,故只考虑能量大于100 MeV/n 的粒子.CREME96 模型计算得到的宇宙线能谱数据没有区分3He 和4He,需要对宇宙线He 能谱进行参数化.针对3He/4He 比例C 的测量试验已经获得比较准确细致的结果[17−20],表2 和表3 中展示了太阳活动极大年参数C(M)、太阳活动极小年参数C(m)的值[21],根据公式由3+4He 的通量可以分别得到3He 和4He 的通量[22,23]:

图3 1 AU 处GCR 各粒子微分能谱Fig.3.Differential energy spectra of each element of GCR at 1 AU.

图4 1AU 处1989 年SEP 微分能谱Fig.4.1989 SEP differential energy spectra at 1AU.

其中 C 为3He 与4He 的比例;F 的单位为m–2·s–1·Sr–1·(GeV/n)–1.

根据表2 和表3 太阳活动极大、极小年的参数化3He/4He 比例C,代入3He 和4He 的通量计算公式(1a)式和(1b)式中,结合由CREME96 模型计算得到的3+4He 的通量,计算得到太阳活动极大、太阳活动极小年3He 和4He 的通量,如图5 所示.

图5 宇宙线Proton 和He(3He 和4He)微分能谱Fig.5.Differential energy spectra of cosmic rays Proton and He (3He and 4He).

表2 太阳极小年3He/4He 的参数化比例C(m)Table 2.The parameterized ratio C(m) of 3He/4He in solar minimum.

表3 太阳极大年3He/4He 的参数化比例C(M)Table 3.The parameterized ratio C(M) of 3He/4He during solar maximum.

2.3 物理模型

GCR 和SEP 的高能性和强子性意味着入射到航天器中会与航天器材料发生复杂的物理过程,并产生大量次级带电粒子.将这些次级带电粒子追踪到尽可能低的能量,就能更加准确的模拟测试质量的充电过程.本文利用GEANT4 模拟包含从数百eV 到TeV 能量范围的物理模型,并且对每个粒子和它产生的所有次级粒子进行全程跟踪.不同种类和能量的粒子与相应物理过程紧密联系,考虑的物理过程为:低能电磁过程、强子过程、衰变过程和光核过程,随着粒子能量的增高与其对应的物理过程也不断变化.对于低能电磁过程,GEANT4默认电磁过程次级粒子(e–,e+,γ)产生最低阈值为250 eV(对应的截断长度约为50 nm),低于这个能量的粒子将不再进行跟踪.

3 仿真结果

3.1 GEANT4 不同截断长度的充电模拟

在GEANT4 中,电磁过程产生次级粒子的最低阈值为250 eV.本文利用太阳活动极小年质子能谱对不同截断长度下的测试质量充电率进行模拟计算,当50 万质子入射航天器模型,分别得到截断长度为1 mm,100 µm,10 µm,1 µm,100 nm和50 nm 时的充电速率,如图6 所示.

图6 不同截断长度的充电速率Fig.6.Charging rate with different cut-off lengths.

如图6 所示,不同截断长度对太阳活动极小年质子的充电速率会产生影响,总体来说,截断长度越低(次级粒子产生的能量阈值越低)模拟得到的充电速率会越高.在截断长度为50 nm 时充电速率最大,为32.821 +e/s.随着截断长度的减小,截断长度从1 mm 减小到50 nm 充电率增长了8.74%.表明截断长度的变化不会显著的影响测试质量的充电速率.

3.2 太阳活动极大、极小年GCR 和SEP的充电模拟

在太阳活动极小年下,仿真模拟的宇宙线Proton,3He,4He,C,N,O,Ne,Mg,Si 和Fe 粒子的参数如表4 所列,获得这些粒子入射到航天器后造成测试质量的充电速率如图7 所示.因为C,N 和O 等重离子在银河宇宙线总的通量占比仅为1%,与Proton,3He 和4He 相比可以忽略,导致各重离子入射航天器得到测试充电率非常小,所以在太阳极大年将不考虑重离子的贡献,仅考虑Proton,3He和4He 对充电率的贡献.在太阳活动极大年下,仿真模拟的宇宙线Proton,3He 和4He 粒子的参数如表5 所列,获得这些粒子入射到航天器后造成测试质量的充电率分布如图7 所示.

表4 太阳极小年宇宙线主要粒子仿真参数Table 4.The main particle simulation parameters of cosmic rays during solar minimum.

表5 太阳极大年宇宙线主要粒子仿真参数Table 5.The main particle simulation parameters of cosmic rays during solar maximum.

如图7 所示,在太阳活动极小年银河宇宙线各粒子入射的总充电速率为39.469 +e/s:入射Proton的充电速率为32.821 +e/s;入射3He 的充电速率为1.166 +e/s;入射4He 的充电速率为5.024 +e/s;入射其它重离子的充电速率为0.458 +e/s.其中质子充电率占比83.16%;3He 的充电率占比为2.95%;4He 的充电率占比为12.73%;其它重离子的充电率占比为1.16%.在太阳活动极小年,Proton的充电速率贡献占主导地位,质子和4He 的充电速率贡献占总充电速率的95.89%,其他重离子的贡献仅占充电率的1.16%,可以忽略.

图7 宇宙线各粒子的充电率 (a) 太阳活动极小年;(b) 太阳活动极大年Fig.7.The charging rate of each particle of cosmic rays:(a) Solar minimum year;(b) solar maximum year.

在太阳活动极大年银河宇宙线各主要贡献粒子入射的总充电速率为12.531 +e/s:入射Proton的充电速率为9.143 +e/s;入射3He 的充电速率为0.435 +e/s;入射4He 的充电速率为2.953 +e/s.太阳活动极大年的总充电率为太阳极小年(GCR最恶劣情况)的三分之一.其中质子充电率占比72.96%;3He 的充电率占比为3.47%;4He 的充电率占比为23.57%.可以看到,太阳活动极大年质子的充电率贡献占比最大,质子和4He 的充电速率贡献占总充电速率的96.53%.

1989 年9 月29 日SEP 峰值能谱模拟仿真参数如表6 所列,模拟了能量范围为0.1—30 GeV 的太阳质子.仿真得到在1989 年9 月29 日的太阳爆发极端环境下的充电速率为120700 +e/s,是太阳活动极小年时充电率的3000 多倍.对于零星发生的SEP 事件,对航天器造成极大的充电风险,极短的时间内就会达到引力波探测能忍受的最大加速度噪声.

表6 1989 年9 月29 日太阳高能粒子事件仿真参数Table 6.Simulation parameters of the SEP event on September 29,1989.

3.3 太阳活动极小年GCR 初、次级粒子贡献和充电能力模拟

在太阳活动极小年,计算了银河宇宙线各粒子入射航天器产生的初、次级粒子造成的测试质量充电速率,如表7 所列.宇宙线各粒子的初级粒子充电速率比例都占总充电率的60%以上(除通量过小的Fe).对于充电率贡献最大的质子和4He,初级粒子充电率占比分别为67.07%和100%.所有宇宙线的初级粒子充电率占比为72.86%,表明测试质量充电速率的主要贡献来源于初级粒子直接沉积在测试质量中.

表7 宇宙线初、次级粒子造成的充电率Table 7.Charge rate caused by primary and secondary particles of cosmic rays.

如图8 所示,在太阳活动极小年,模拟计算得到宇宙线各粒子的充电能力.值得注意的是,在初级宇宙线中占比为8%的He(3He 和4He),贡献了总充电率为15.68%;在初级宇宙线中占比为1%的重粒子,贡献了总充电率的1.16%.针对ASTROD I也有类似的结果[11].这是由于银河宇宙线入射航天器测试质量的充电能力不同(各粒子的电离能力和射程的不同),4He 充电能力约为质子的2 倍,重粒子的充电能力基本都大于质子的充电能力.

图8 宇宙线各粒子的充电能力Fig.8.The charging ability of each particle of the cosmic ray.

3.4 太阳活动极小年GCR 质子不同能量段的充电模拟

在太阳活动极小年,仿真模拟从0.1—100 GeV范围内不同能量的银河宇宙线质子各向同性入射到航天器模型,对于每个单能点仿真模拟10000 个事例(其它条件与能谱仿真时保持相同),得到了不同能量质子的充电能力(单个粒子入射造成的测试质量净充电个数).如图9 所示,能量为200 MeV时充电能力有峰值0.0498 +e/primary,在高能部分(大于10 GeV)充电能力迅速增大为另一个峰值0.2077 +e/primary.

图9 太阳活动极小充电能力 (a) 不同能量充电能力和质子通量;(b) 充电能力和结合质子通量的总充电比例Fig.9.Solar minimum charging capacity:(a) Different energy charging capacity and proton flux;(b) charging capacity and total charge ratio of combined proton flux.

结合太阳活动极小年时入射的质子能谱,可以得到对于太阳极小年不同能量的初级宇宙线质子对测试质量充电的贡献.设仿真模拟的粒子数为N0,测试质量的净充电个数为N,充电能力为CA,充电率为CR,各能量点质子的微分通量为F,各向同性粒子出射的半径为R.

由(2a)式和(2b)式可到:

其中E 为质子的能量,单位为GeV.

由(2a)式和(2b)式得到(3)式,可以反应不同能量下的充电率贡献.结合质子能谱的初级宇宙线质子对测试质量充电的贡献如图9 所示,充电能力峰值的高能部分由于能谱通量远低于低能部分,导致最终的充电率贡献很低.计算可以得到,在能量区间0.1—10 GeV 内,对测试质量充电的占比约为94%.在能量区间0.1—1 GeV,对测试质量充电的占比约为65%.仿真结果可为在轨粒子探测器的设计提供依据.

4 讨论

本文使用Geant4 全面仿真了引力波探测航天器在距太阳1 AU 处受到宇宙线高能粒子的轰击导致测试质量充电的规律.本文的仿真计算考虑了更为全面的初级和次级粒子、更低的能量阈值,与LISA 探路者飞行试验结果相吻合,可为后续引力波探测任务中在轨电荷实时监测、预报和管理提供重要依据.Grimani 等[12]利用Fluka 仿真模拟了参数化的太阳极大、极小年质子和氦核能谱注入LISA 探路者相对复杂几何模型中的充电率,得到太阳极大年充电率为15.13 +e/s,太阳极小年充电率为38.2 +e/s.该结果与本文的研究结果高度一致,但是由于Fluka 软件的最低截断能量为1 keV 和Grimani 等[12]使用的能谱为参数化能谱,进而会造成一定的误差.Araújo 等[15]利用Geant4仿真模拟了LISA 模型在太阳极大、极小年测试质量的充电率,但此工作中没有考虑到e+在测试质量中的沉积,最终得到太阳极大年充电率为23.7 +e/s,太阳极小年充电率为49.9 +e/s,而本文全面地考虑3—20 号元素、proton,3He,4He,e–,e+、氘、氚和正负π 介子,更综合地分析了测试质量的充电情况.空间验证试验结果也证明了本文仿真模拟的正确性,LISA 探路者作为LISA 的验证计划于2015 年发射,已经得到其上搭载的2 个测试质量的充电率,TM1 和TM2 的充电率分别为22.9 +e/s 和24.5 +e/s[14],此次的测量是在太阳极小前3—4 年,测试质量的充电率在本文仿真的预期范围内12.531—39.469 +e/s.在此基础上,本文进一步研究分析了典型SEP 充电率、太阳活动极小年银河宇宙线的初级粒子充电率占比以及太阳活动极小年银河宇宙线质子在不同能量区间的充电率贡献.研究表明典型SEP 峰值能谱通量下的充电速率为120700 +e/s,认为零星发生的SEP事件对航天器造成极大的充电风险,太阳活动极小年银河宇宙线的初级粒子充电率占比最大,针对太阳活动极小年测试质量总充电贡献占比最大的银河宇宙线质子,其低能部分充电率贡献较大.

本文仿真方法更切合实际地模拟分析了宇宙线高能粒子对测试质量的充电率,更进一步研究了太阳极小年各银河宇宙线粒子的充电率贡献以及对银河宇宙线质子在不同能量区间的贡献.此项工作可以讨论由于测试质量充电而引起的加速度噪声和它的频谱分析,以及对未来引力波探测计划中粒子探测器的粒子探测种类和能量区间提供一些启示,最终设置合理的放电方案对引力波探测的航天器进行电荷控制,以降低由于测试质量充电所引起的噪声.本文的研究结果可为未来中国引力波探测计划测试质量的加速度噪声预测、粒子探测器的设计和放电方案提供理论依据.

5 结论

本文研究了高能宇宙线粒子对测试质量的充电机制,建立了等效航天器模型来仿真空间引力波探测中的航天器的结构与材料.通过GEANT4 对航天器在太阳活动极大、太阳活动极小年银河宇宙线各粒子和SEP 的轰击下测试质量充电过程的仿真模拟,可以得出以下结论:

1) 在同一能谱不同截断能量下,随着截断长度的降低充电速率逐步增大,截断长度从1 mm 减小到50 nm 充电率增长了8.74%.在本文几何模型、物理过程和能谱设置条件下,截断长度的变化不会显著的影响测试质量的充电速率;

2) 太阳活动极小年时的总充电速率为39.469+e/s,除Proton,3He 和4He 之外,其他重粒子对充电率的贡献仅为1.16%;太阳活动极大年的总充电速率为12.531 +e/s,为太阳活动极小年的三分之一;1989 年9 月29 日SEP 峰值能谱通量下的充电速率为120700 +e/s,零星发生的SEP 事件对航天器造成极大的充电风险;

3) 太阳活动极小年时,银河宇宙线的初级粒子充电率占比为72.86%,其中充电率贡献最大的是质子和4He,初级粒子充电率占比分别为67.07%和100%,表明测试质量充电速率的主要贡献来源于初级粒子直接沉积在测试质量中;

4) 太阳活动极小年时Proton 在总充电贡献中占比最大,为83.16%,He 的总充电率贡献为15.68%,是初级宇宙线中He 通量占比的2 倍,其他重离子对测试质量的充电率贡献不大;

5) 太阳活动极小年时,宇宙线质子低能部分充电率贡献最大,在能量区间0.1—1 GeV 的充电率贡献约占总充电率贡献的65%,在能量区间0.1—10 GeV 内,对测试质量充电率的占比约为94%.

感谢国家科技基础条件平台-国家空间科学数据中心(https://www.nssdc.ac.cn)对本文空间环境数据的支持.