质子碰撞硼原子非辐射的电荷转移过程.pdf

质子碰撞硼原子非辐射的电荷转移过程*朱宇豪1)袁翔2)吴勇2)王建国2)1)(西安建筑科技大学理学院,西安710055)2)(北京应用物理与计算数学研究所,计算物理重点实验室,北京100084)(2023年 3月 27 日收到;2023年 6月 3 日收到修改稿)103103eV/u重粒子碰撞中的电子转移涉及复杂的电子关联机制,极大地影响等离子体中电荷态平衡,也是 X 射线的辐射的重要来源之一.电子转移截面与速率系数是国防工业发展核聚变等离子体所需要的重要原子参数.基于全量子的非辐射分子轨道密耦合方法,系统研究了质子碰撞硼原子在 能量区间内的硼原子电子转移过程.计算采用多参考组态方法得到总共 15 个电子转移、激发以及弹性散射的通道,每个通道对应的分子态能量得到了与实验符合较好的结果.分子态的绝热势能曲线间的避免交叉现象明显,构成了电子转移的主要途径.计算发现,质子碰撞硼原子过程中 2s 轨道的电子转移是占主导地位,2p 轨道的电子转移贡献较小.在低能区,电子转移截面出现了明显的量子共振现象,这些共振主要来源于不同能量通道的耦合.此外,还计算了不同温度下的质子碰撞硼原子的电子转移速率,该速率可为复杂等离子体环境的模拟诊断提供重要的原子参数支持.关键词：电子转移,重粒子碰撞,势能曲线PACS：34.10.+x,34.20.Cf,31.10.+zDOI:10.7498/aps.72.202304701引言重粒子碰撞中的电子转移涉及多电荷动力学以及多体相互作用等复杂物理过程,因为其在天体和等离子体物理中有非常重要的应用,吸引了广泛的理论和实验研究1.在磁性约束核聚变中,质子对电子的俘获极大影响了整个系统的电荷态平衡以及等离子体冷却过程.电荷态平衡是诊断复杂等离子体系统基础的指标,从中可以获得元素丰度、电荷密度等物理信息2.在天体环境中,X 射线谱是识别天体可观测气体元素重要的手段之一,并且各类元素的 X 射线辐射都与电子转移有密切的关系.特别对于内壳层电子的转移过程,是原子系统X 射线辐射的重要来源3.所以,为了模拟上述X 射线谱和等离子体中电荷态平衡等过程,获得重粒子碰撞中的电子转移态选择截面和总截面是很有意义的.那么,重粒子碰撞的电子转移截面等相关信息已经成为不可或缺的原子参数,为国防工业发展核聚变技术,以及为等离子体物理等学科中演化模拟复杂等离子体系统提供强有力的原子参数支持4,5.近年来,对重粒子碰撞的理论研究也迅猛发展,在不同的入射能区,理论计算可以做到与实验测量符合较为精确的结果.如在中高能区的原子轨道密耦合方法69,能够精确计算重粒子碰撞的电子激发和转移的态选择截面以及总截面.因为其只对电子做全量子化的处理,而原子核依然采用经典力学的运动方式,所以该方法为半经典的理论方法.在低能入射区,量子效应更加明显,需要采用全量子的碰撞方法才能考虑足够的物理信息得到精确的电子转移截面.目前,较为常见的全量子的理*国家自然科学基金(批准号:11934004)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-1论方法是非辐射分子轨道密耦合方法,Liu 等6,1014基于此方法成功计算出了 He,Li 以及 Be 等原子系统的电子转移截面,并与相关精密实验测量结果符合得较好.全量子的分子轨道方法,因为分子轨道模型可以准确描述较低能区的散射态波函数,并能够有效处理各个分子轨道之间的耦合,其在计算低能区对电子转移的截面计算属于国内外公认的最高精度方法.当然,分子轨道方法的精确度主要取决于分子态能量,即分子电子态结构的精度.这要求在计算分子结构时,需考虑足够大的组态空间和足够多的电子组态,充分考虑电子关联对分子结构的影响14.110 keV/u0.11.0 keV/u在实验研究方面,重粒子碰撞是阐明量子多体关联,揭示靶粒子内部结构信息的重要科学手段.国内外已经开展了数十年的重粒子碰撞实验,诸多质子或氢原子与重元素原子不同价态系统的碰撞实验都相继展开.比如,Bruhns 等15,16使用合并束技术研究了氢原子与 Si3+离子在地能区碰撞的电子转移过程,并精密测量了电子转移截面.Gieler等17,18研究了质子碰撞 K 原子的低能区电子转移总截面,其中 Gieler 等17的碰撞能区在 ,Ebel 等18的研究能区在 .虽然诸多原子系统与质子碰撞中电子转移过程被实验陆续研究1922,但是据目前调研所知,还从未开展对 B原子的实验研究,然而 B 原子在天体中也具有较高丰度,因此迫切需要 B 原子参数.此外,B 原子与质子的碰撞过程也是研究多体相互作用系统的良好平台,其中涉及多体相互作用的电子转移机制还没有被发掘.因此,质子碰撞 B 原子的理论研究,不仅可以更深一步了解其电子转移机制,还能认识碰撞过程中多电子关联效应.本文基于全量子非辐射的分子轨道密耦合方法,系统研究大范围碰撞能区内质子碰撞 B 原子电子转移过程.考虑足够多的电子态,以确保整个计算中通道的收敛性.第 2 节详细阐述计算的理论方法;第 3 节展示计算结果并讨论相关电子转移的物理机制;第 4 节对研究工作做简短的总结.2理论方法 2.1 分子电子态结构计算分子电子态结构的计算,采用 Buenker 等2325提出的 Table-directCI 结合多参考单双激发组态相互作用方法(MRDCI).该方法在自洽场计算(HF-SCF)之后,能基于单双激发的参考组态,可以最优化选取占比重要的组态,缩小了组态空间,且可以保证精度.则说明 MRDCI 可以快速优化组态构成,得到精确的电子态波函数,其具体计算公式为(PJM)=nr=1cr(r)r(PJM),(1)(PJM)r(PJM)cr(r)其中,是电子态波函数,为参考组态考虑单双激发后的优化的组态波函数,是混合系数.实际计算中,通常选取对内、外壳层和极化作用考虑较好的芯-价关联极化的多Zeta 基组,加入弥散基来弥补由于高激发态离核较远导致电子云分布较大的影响.做完自洽场计算后,构建组态空间,对于选定的参考组态进行单双电子激发,但是这样经常会导致组态空间偏大,因此可以通过设定阈值、结合“有效势场”控制活动电子数量来缩小组态空间,选取贡献较大的组态.2.2 重粒子碰撞动力学计算非辐射分子轨道密耦合方法(QMOCC)对于重离子散射系统,采用质心坐标和原子单位,两原子体系碰撞系统的总哈密顿量表示为H(R,r)=122R12mNi=12ri12Mi=ji j+V,(2)RVM=M1+M2=M1M2/Mm=M/(M+1)式中,为入射离子相对靶粒子核间距,为所有电子以及原子核的库仑相互作用之和,质量项,注意(2)式忽略了极化相互作用项.iR 令 为电子总波函数,在 时候,电子相对于原子核运动,以核为中心,电子总波函数能在绝热表象下有如下形式:i(R,r)=i(R,r)expi(R,r),(3)expi(R,r)式中,称电子转移因子.在微扰的稳态近似下,可以忽略电子转移因子,把系统的总波函数近似用绝热电子波函数代替,利用 Born-Oppenheimer 近似,则总波函数可以用电子绝热波函数展开:(R,r)=iFi(R)i(R,r),(4)物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-2iii其中对 的求和限于完备的绝热矢量集合 中的子集.而且,每一个 在渐进区对于这不同的原子态,这些原子态也称为通道.将(4)式代入能量本征方程,可得到微扰的稳态方程:2RI+2 u EIF(R)=M(R)+P(R)RF(R),(5)Iuuij=jijjM(R)P(R)式中,为单位矩阵,为对角的绝热势矩阵,是绝热态的本征能量,和 是耦合矩阵分别满足:Mij(R)=dri(R,r)2Rj(R,r),(6)Pij(R)=2dri(R,r)Rj(R,r).(7)在球坐标系里,(5)式可以写为d2dR2+J2 2R2+2(u E)bRF(R)=VR(R)+VC(R)bRF(R),(8)bRRJi|L2z|j=2jijVR(R)VC(R)其中,为 方向的单位矢量,为系统总角动量,是对角矩阵,且满足 .和 分别为径向耦合矩阵和转动耦合矩阵,具体为VRij(R)=i|2R2+L2x+L2yR2|j2i|R|jRij,(9)VCij(R)=1R22i|iLy|j 2i|iLx|j1sin 2ji|iLx|jcot.(10)对波函数作分波展开:Fj(R)=R1J=0fJj(R)DJMj(,),(11)DJMj(,)其中 为角向波函数,代入(8)式,则得到径向耦合方程:d2dR2+J(J+1)2R2+2 u(R)EfJ(R)=VR(R)+VC(R)fJ(R).(12)glfl=CglCfg对于求解(12)式,定义一个新的函数 ,满足,其中 是正交变换矩阵.若 是(12)式的一个解,那么 同样是另一个解,因此可以这样C构建 :dCdR+AC=0.(13)limR C(R)I取边界条件 ,(12)式可以变为d2dR2+J(J+1)2R2 2EgJ+2U,(R)gJ=0,(14)U,(R)C(u P)C1,(15)uU,P式中 是对角的绝热势矩阵,是非绝热势矩阵,是转动耦合矩阵.(14)式满足以下边界条件:limRgJ(R)1k,jJ(kR)+KJ,J(kR),(16)k=2E (),(17)jJJKl其中 和 分别为通道的规则和非规则的库仑函数,为实对称的 K 矩阵.从 K 矩阵可以得到S 矩阵:Sl=I+iKlI iKl.(18)i因此,从初态 到末态 j 的散射截面为ij=k2il(2l+1)?Sl?2i,j,(19)Sll其中 表示第 个分波对应的 S 矩阵.对于速率系数的计算,采用 Maxwellian-Bolt-zmann 分布对散射截面求平均,那么反应速率为(T)=(8)1/2(1kT)3/20E(E)exp(EkT)dE,(20)kT其中,为玻尔兹曼常数,为温度,是约化质量.3计算结果和讨论2+2对于 BH+分子的电子态,总共计算了 9 个 态和 6 个 态,这 15 个电子态的绝热势能曲线如图 1 所示,横坐标为原子核间距,纵坐标是电子态能量.从图 1(a)可看到较低的 5 个分子态之间有较大的能量差,较高的激发态能级比较密集,处于近简并的状态.电子态之间的避免交叉近区相护作物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-3R=3.8 3242R=2.3,3.2,4.6 4252用很强烈,这些避免交叉现象在耦合矩阵元中也有明显的表现.图1(b)中,在 处,和 态有一个显著的避免交叉,而分别在 处,和 也有比较明显的避免交叉现象.对于 QMOCC 的计算,这些避免交叉能量相近的点起到关键的作用,因为绝热分子态的耦合构成了非辐射电子转移的主要途径.为了保证后续动力学计算的准确性,表 1 列出了这 15 个电子态的理论结果和实验结果,以及各个不同分子态解离极限所对应的原子态.在QMOCC 的理论框架下,散射末态对应的原子能级与分子态解离极限对应的原子态一一对应.上述计算构建的散射矩阵,可以清晰说明不同散射通道之间跃迁的物理信息.从表 1 还可以看到,分子态32+22的解离极限涉及的原子态能级最大到 2p 电子轨道,电子从 B 原子的转移主要来源于 2s 和 2p 电子轨道.除了非弹性的散射通道之外,计算还包括一些弹性散射通道,如分子态 和 .这些弹性散射通道与非弹性散射通道的耦合,对最终散射态之间的跃迁有不可忽略的影响.103103eV/u3P1P1P表 1 中所有的散射通道都被考虑在 QMOCC计算内,计算能区从 .图 2 和图 3分别为 B 原子 2s 和 2p 轨道的电子转移截面随入射能量变化的函数图像.在图 2 中,对分子的态的精细结构也进行区分,分别为 和 态.因为 态的散射通道能量较高,所以需要较大的入射能量,该通道才被打开.在图 3 中,将两个 2p 轨道的电子转移截面加和,得到 2p 轨道的电子转移总截-25.0-24.8-24.6-24.4-24.2-24.0-23.8-23.6(a)2468101214161820AEnergy/a.u.12S+22S+32S+42S+52S+62S+72S+82S+92S+B(1s22s22p)+H+B+(1s22s2p)+H(1s)B+(1s22s2)+H(1s)-25.0-24.8-24.6-24.4-24.2-24.0-23.8-23.6(b)2468101214161820AEnergy/a.u.12P22P32P42P52P62PB(1s22s22p)+H+B+(1s22s2p)+H(1s)图1(a)BH+分子 态势能曲线;(b)BH+分子 态势能曲线Fig.1.(a)EnergycurveofBH+moleculewith states;(b)energycurveofBH+moleculewith states.表1不同分子态对应的解离极限下的原子态以及态能量的理论结果和实验结果Table1.Atomicstatescorrespondingtodifferentmolecularstatesunderdissociationlimit,andthetheoreticalandexperi-mentalresultsofstateenergy.原子态分子态能量/eV理论实验误差B+(1s22s2)+H(1s)12+000B+(1s22s2p)+H(1s)22+12 4.62264.62300.0004B(1s22s22p)+H+32+22 5.30425.30030.0039B+(1s22s2p)+H(1s)42+32 9.06639.10010.0338B+(1s22s2)+H(2p)52+42 10.176910.19870.0218B+(1s22s2)+H(2s)62+10.192910.19870.0058B(1s22s23s)+H+72+10.261510.25970.0018B(1s22s2p2)+H+82+52 11.249111.23370.0154B(1s22s23s)+H+92+62 11.362311.32730.0350物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-4103eV/u面.然而,从图 3 可以得到,2s 和 2p 轨道的电子转移截面在计算的能区范围内,有至少两个数量级的差异.所以内壳层的 2s 轨道的电子转移在整个电子转移过程中占主导地位.因为 QMOCC 考虑各种电子态的散射通道之间的耦合,从图 3 可以看到,在低于 1eV/u 能量,电子转移截面具有较多的共振现象.尤其在靠近 能量附近,2s 和 2p 轨道的电子转移截面都有强烈的共振,这是低能区碰撞过程量子效应比较强的一种表现.重粒子碰撞过程中,往往内壳层的电子转移都会伴随退激发光辐射现象.12+22+12+42+32+42+52+62+82+52+32+22对于质子碰撞 B 原子,如表 1 所列,电子转移通道包括 ,和 这 10 个绝热分子态,其中涉及内壳层 2s 电子转移的通道有 6 个,导致外壳层 2p 电子转移的通道有 4 个.弹性通道有,两个绝热分子态,其开通道能量为 5.3eV72+92+62+52+62+22+32+12+22+42+52+左右,所有其在低能区并未与电子转移通道存在竞争关系.此外,还有 3 个电子从 2p 轨道跃迁到 3s轨道的激发通道,分别是 ,和 绝热分子态,激发通道都需要较大的能量才能打开,主要贡献都在于与其他电子转移通道耦合.图 4 和图 5为部分耦合较强的绝热分子态之间的径向耦合矩阵元和转动耦合矩阵元随核间距的变化.从图 4 可以看到,径向耦合中,分子态 之间的耦合最强,这两个通道都对 2p 轨道的电子转移有贡献.而且,除了分子态 的耦合为 2s轨道的电子转移有影响,剩下几个耦合比较强的径向耦合矩阵元,比如 和 也都对 2p 轨道的电子转移有影响.然而,2s 轨道的电子转移涉及较多的绝热分子态通道,所以其总体截面大于 2p 轨道电子转移.虽然,影响 2p 轨道电子转移的径向耦合强度较大,但是图 3 的计算说明,这些耦合更多的造成 2p 轨道电子转移截面的减小.5/ARadial matrix element/a.u.010152040353025201510512S+-22S+22S+-32S+32S+-42S+42S+-52S+52S+-62S+62S+-72S+72S+-82S+82S+-92S+1.00.5Radial matrixelement/a.u.0-0.5-1.05/A图4径向耦合矩阵元随核间距的变化Fig.4.Radial coupling matrix element changes with thenucleardistance.42525262425252622s42526242同样图 5 中最大的两个转动耦合矩阵元分别是分子态 和 .其中 之间的转动耦合也是影响 2p 轨道电子转移,但 之间的耦合却对 轨道的电子转移有贡献.,三个态是较高的激发态,因为要打开最低能量的态 也需要一定入射能量,所以转动耦合虽然大,但是对电子转移过程在低能区的影响极其有限.涉及 2s 轨道的电子转移通道,虽然耦合贡献不及 2p 轨道电子,但是这些通道对应的分子态大多处于较低能级,且数量较多,所以最终计算的 2s 轨道的电子转移截面远大于 2p 轨道电子.10310-310-2100101/(eVSu-1)102103Crossing section/(10-16 cm2)10-610210110010-110-210-310-410-510-1B+3P(2s2p)B+1P(2s2p)图2质子碰撞 B 原子过程中 2p 轨道电子转移截面Fig.2.Electrontransfercrosssectionof2porbitalsintheprocessofprotonscollidingwithBatom.10-310-2100101/(eVSu-1)102103Crossing section/(10-16 cm2)10-610210110010-110-210-310-410-510-1B+(1s22s2)B+(1s22s12p1)图3质子碰撞 B 原子过程中 2s 和 2p 轨道电子转移截面Fig.3.Electrontransfercrosssectionof2sand2porbitalsintheprocessofprotonscollidingwithBatom.物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-5(T)在表 2 中,考虑整个系统服从麦克斯韦-玻尔兹曼能量分布.计算的电子总的转移速率系数,计算的温度点量级从 101107,这个温度跨度范围比较广泛,可以反映宏观系统电子转移的物理性质,注意计算的电子转移的速率将包含 2s 和2p 轨道.从表 2 还可以发现,因为相对高能态,低能态的布局对速率系数的影响较大,随着温度的升高,电子转移总速率系数增长幅度越来越平缓.当温度高于 105K,系统的电子转移总速率变化幅度明显小于低温区.虽然我们考虑总的电子转移速率,但是 2s 轨道电子转移占绝对主导地位,所以总的电子转移速率,也可以在一定程度上反映内壳层电子转移后导致的 X 射线辐射的强度.4结论2+2采用 QMCOCC 全量子的理论方法,研究了质子碰撞 B 原子电荷转移过程.通过 MRDCI 的第一性原理计算,得到了耦合分子 9 个 态和6 个 态,取得了与实验测量较好的结果.得到的103103eV/u分子态势能曲线可避免交叉,这些能量点的避免交叉对绝热分子态耦合具有重要的意义.本文所计算的 15 个散射通道中,不仅包含弹性散射通道,还有电子激发和转移的非弹性散射通道.其中,电子转移通道都集中在中低能区,而电子激发处于较高能区.经过计算,得到了从 的质子入射能区的电子转移截面,并发现在低能区出现较强的共振现象,这主要是因为低能区通道耦合造成的散射量子效应.电子转移截面分为 2s 和 2p 轨道的电子转移截面,其中 2s 轨道的电子转移截面占据绝对主导地位.这些内壳层的电子转移过程,往往伴随这高能量的辐射退激发发光想象.从 15 个通道分析,发现包含 2s 轨道的电子转移的绝热分子态数量最多,且这些通道的所需要的开通道能量也比较小.此外,还计算了大能量跨度的电子转移速率系数,为等离子体模拟提供重要的原子参数.对于宏观大尺度系统的电子转移过程,低能区的散射截面占比较大.通过计算也验证了这点,研究发现电子转移速率随着系统温度的升高,增长的幅度越来越小.众所周知,B 元素在天体中的丰度也较高,所以第一性原理计算不仅可以提供精确的原子参数,也可以帮助人们深入认识质子碰撞 B 原子的物理机制.感谢西安建筑科技大学理学院和北京应用物理与计算数学研究所对本工作的支持.参考文献 JanevRK,PresnyakovLP,ShevelkoVP1985Physics ofHighly Charged Ions(Berlin-Heidelberg:Springer)pp161LoarteA,LipschultzB,KukushkinAS2007Nucl.Fusion47S2032KallmanTR,PalmeriP2007Rev.Mod.Phys.79793LiuCH2009Ph.D.Dissertation(Beijing:Universityofthe4表2电子转移总的速率系数,其中包含 2s 和 2p 轨道的电子转移Table2.Totalratecoefficientofelectrontransfer,includingelectrontransferof2sand2porbitals.T/K 2 5 10 40/(1012cm3s1)1.72 108 5.09 104 1.40 102 1.24 101 T/K 100 500 1000 6000/(1014cm3s1)14.7 10.3 8.62 46.9 T/(105K)0.1 0.5 1.0 5.0/(1012cm3s1)1.72 8.91 101 3.86 102 4.74 103 T/(106K)1.0 5.0 10 50/(109cm3s1)9.38 25.8 26.8 7.19 515/A20Radial matrix element/a.u.-3.03.02.52.01.51.00.50-0.5-1.0-1.5-2.0-2.51012P-22P22P-32P32P-42P42P-52P52P-62P图5转动耦合矩阵元随核间距的变化Fig.5.Rotation coupling matrix element varies with thenucleardistance.物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-6ChineseAcademyofSciences)(inChinese)刘春华2009博士学位论文(北京:中国科学院大学)LinXH2019Ph.D.Dissertation(Beijing:BeijingInstituteofTechnology)(inChinese)林晓贺2019博士学位论文(北京:北京理工大学)5LiuL,LiuCH,WangJG,JanevRK2011Phys.Rev.A840327106Gao J W,Wu Y,Sisourat N,Wang J G,Dubois A 2017Phys.Rev.A960527037Gao J W,Wu Y,Wang J G,Sisourat N,Dubois A 2018Phys.Rev.A970527098ZhangYW,GaoJG,WuY,ZhouFY,WangJG,SisouratN,DuboisA2020Phys.Rev.A1020228149LiuCH,WangJG,JanevRK2012Phys.Rev.A8504271910YanLL,WuY,QuYZ,WangJG,BuenkerRJ2013Phys.Rev.A8802270611LiuCH,LiuL,WangJG2014Phys.Rev.A9001270812LiuCH,WangJG2013Phys.Rev.A8704270913LiuCH,WangJG,JanevRK2014Phys.Rev.A8906271914BruhnsH,KreckelH,SavinDW,seelyDG,HavenerCC152008Phys.Rev.A77064702KimHJ,PhaneufRA,MeyerFW,StelsonPH1987Phys.Rev.A1785416GielerM,ZiegelwangerP,AumayrF,WinterH,FritschW1991J.Phys.B2464717EbelF,SalzbornE1987J.Phys.B20453118Hayakawa S,Kadomura K,Kimura M,Dutta C M 2004Phys.Rev.A7002270819Hayakawa H,Helm H,Briggs J S,Salzborn E 2007 Phys.Rev.Lett9917320220KimHK,SchofflerMS,HouamerS,ChuluunbaatarO,TitzeJN,SchmidtLPH,JahnkeT,Schmidt-BockingH,GalstyanA,PopovYV,DornerR2012Phys.Rev.A8502270721BayfieldJE,KhayrallahGA1975Phys.Rev.A1192022Buenker R J,Liebermann H P,Whitten J L 2001 Chem.Phys.265123BuenkerRJ,LiebermannHP,IzgorodinaEI2003Chem.Phys.29111524AlcheevPG,BuenkerRJ,ChernovV,ZonBA2003J.Mol.Spectrosc.21819025物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-7Non-radiative charge transfer process ofproton impcating B atom*ZhuYu-Hao1)YuanXiang2)WuYong2)WangJian-Guo2)1)(School of Science,Xian University of Architecture and Technology,Xian 710055,China)2)(Key Laboratory of Computational Physics,Institute of Applied Physics and Computational Mathematics,Beijing 100084,China)(Received27March2023;revisedmanuscriptreceived3June2023)Abstract103 103eV/uElectrontransferinheavyparticlecollisions,involvingcomplicatedelectroncorrelationmechanisms,greatlyaffectsthechargestatebalanceinplasma,andisalsooneofimportantsourcesforX-rayradiation.Electrontransfercrosssectionandratecoefficientareimportantatomicparametersrequiredforthedevelopmentofnuclear fusion plasma in the national defense industry.We systematically investigate the electron transferprocessofprotonimpactingboronatominanenergyrangeof basedonafullyquantumnon-radiativemolecularorbitalclose-couplingmethod.Atotalof15channelsofelectrontransfer,excitation,andelastic scattering are obtained by using the multi-reference configuration interaction method,and thecorresponding molecular energy of each channel is in good agreement with the experimental results.Thephenomenonofavoidingcrossingbetweentheadiabaticpotentialenergycurvesofmolecularstatesisobvious,whichconstituentsthemainpathwayofelectrontransfer.Thecalculationshowsthattheelectrontransferofthe2sorbitalisdominatedintheprocessofprotonimpactingboronatom,whiletheelectrontransferofthe2porbitalcontributesalittle.Inthelowenergyregion,thereareobviousquantumresonancesintheelectrontransfercrosssection,whichoriginatefromthecouplingbetweenhighenergychannelandlowenergychannel.Inaddition,wecalculatetheelectrontransferratesofproton-impactboronatomatdifferenttemperatures,whichcanprovideimportantatomicparameterwhichsupportforsimulatinganddiagnosingcomplexplasmaenvironments.Keywords:electrontransfer,heavyioncollision,potentialcurvePACS:34.10.+x,34.20.Cf,31.10.+zDOI:10.7498/aps.72.20230470*ProjectsupportedbytheNationalNaturalScienceFoundationofChina(GrantNo.11934004).Correspondingauthor.E-mail:物理学报ActaPhys.Sin.Vol.72,No.16(2023)163401163401-8