燃烧等离子体中子诊断
燃烧等离子体中子诊断:当等离子体自己开始烧
托卡马克的实验做了几十年,绝大部分时间等离子体是被"喂"出来的——中性束注入、射频波加热,外面灌多少功率,里面就烧多少。诊断要回答的问题是:我们喂进去的能量,等离子体把它用得怎么样。但真正意义上的聚变堆不是这样:D+T反应产生的α粒子会留在磁场里,自己加热等离子体,加热功率超过一切外部加热。这个状态叫燃烧等离子体。它一出现,诊断的处境就全变了——你要在反应最猛烈的时刻,回答"α加热到底把它加热成什么样了""燃料离子是不是还老老实实待在热平衡里""氘和氚的比例是不是还在最优点上"。而这些问题,全都要从一样东西里读出来:中子。
中子是聚变唯一能自由逃出磁场的产物。D+T→α+n(14.1 MeV),D+D→³He+n(2.45 MeV)——中子的产额直接就是聚变功率的账本,中子的能谱则记录着反应时燃料离子的运动状态:谱的宽度量离子温度,谱的偏移量等离子体整体流动,谱里那些不属于热分布的鼓包,量的是快离子。所以燃烧等离子体诊断,核心战役就是中子测量。这篇文章把最近几年这条线上的几个标志性进展捋一遍:JET两轮氘氚实验(DTE2、DTE3)上中子相机怎么独立记账、铍钨壁带来的(n,2n)中子倍增怎么搅局、金刚石探测器怎么把2.5和14 MeV两种中子干净地分开、NIF的点火实验又怎么从能谱矩里抓到了超热离子。所有论文、作者、机构都逐条核实过,中间有一处流传很广的数字张冠李戴,我单独点出来。
先交代背景:JET在2021年(DTE2)和2023年(DTE3)跑了两轮氘氚实验,一共产生超过1.5×10^21个14 MeV中子——占JET四十年全部氘氚中子产量的八成。DTE3的纪录脉冲#104522在约6秒里产出2.46×10^19个中子,峰值产额率5.06×10^18每秒,对应约14兆瓦聚变功率,这就是那个著名的69兆焦耳聚变能纪录。这一轮实验同时是一次诊断大考:因为DTE2之前JET把14 MeV中子产额监测器搬进真空室里做了在位绝对标定(用一台经过标定的中子发生器放在73个位置逐点照射),裂变电离室和活化系统的标定精度首次达到了ITER要求的量级。但标定之后,独立测量之间的账却对不齐了——这正是中子相机那篇论文的由来。
JET的中子相机(代号KN3)围着等离子体装了19条准直视线,每条视线配塑料闪烁体探测器,专门盯着14 MeV中子。19条视线合起来,能重建中子发射的二维分布——这是唯一能"看出"聚变反应发生在等离子体哪个部位的诊断。2023年Uppsala大学的团队先做了一半:用第一性原理方法从相机数据独立估算D-D中子产额,从闪烁体光产额标定、堆叠信号处理、中子输运模拟到发射轮廓拟合,一整套不依赖其他任何测量。DD那篇(Rev. Sci. Instrum. 94, 073502, 2023,第一作者Linus Hägg)跟裂变电离室对上了,平均偏差在10%以内,算是开了个好头。2025年他们把同一套方法搬到DT上,就是那篇RSI 96, 063501(2025):对DTE2和DTE3两轮实验的一批炮做独立产额估计,结果出乎意料——中子相机给出的产额比裂变电离室平均高出34%(DTE2)和41%(DTE3),而且偏差在两轮实验里相当稳定。原因目前不明,作者自己把嫌疑对象列了一排:高产额下的信号堆叠修正不完全、光产额标定的增益漂移、钢准直器的透射模型。这篇论文最可贵的地方恰恰是这个"目前不明"——两套各自标定过、号称百分级精度的系统,账差了三分之一还多,谁也说不清是谁的错。对ITER这是实打实的预警:ITER的氘氚中子产额比JET高六个量级,堆叠效应只会更凶,产额测量这条路还没走通。
第二件事,墙的干扰。JET在2011年换了铍钨的"ITER类壁",取代原来的碳壁。对中子物理来说这一换不得了:铍的(n,2n)反应阈值只有1.85 MeV左右,14 MeV中子打在铍上,一个变两个——第二个中子能量不高,恰好落在2.5 MeV附近。也就是说,在氘氚等离子体里,你测到的"DD中子"信号不再只是DD反应的产物,还混着一层14 MeV中子撞墙碎出来的次级中子。这件事在2013年JET聚变技术活动的总结论文(Fusion Engineering and Design 88, 733, 第一作者Paola Batistoni,EFDA-JET/ENEA弗拉斯卡蒂)里写得明明白白:DD运行下(n,2n)倍增可忽略,DT运行下评估为大于25%。同一篇论文还有另一个数:由于铍的弹性散射增强,内辐照端测到的2.5 MeV中子通量被削掉约20%。这里要特别当心一个流传很广的误记——"20%的(n,2n)倍增"是不存在的,20%是2.5 MeV通量因散射下降的数字,(n,2n)倍增是另一个数,大于25%,两件事不能搅在一起。这个倍增对诊断的杀伤在于燃料比测量:ITER要求实时测氘氚密度比,精度20%、时间分辨100毫秒,主力方案就是比较14 MeV与2.5 MeV两个峰的强度。Uppsala团队2010年在RSI上专门论证过(G. Ericsson等,RSI 81, 2010):当DT与DD产额比超过约70倍时,2.5 MeV信号就彻底淹没在14 MeV中子的散射本底里——而倍增效应让这个本底雪上加霜。2022年他们在PPCF上给出解法(B. Eriksson等,PPCF 64, 055008, 2022):氚浓度高时干脆不用DD峰,改用三体TT反应的连续谱定燃料比,绕开了整片雷区。
第三件事,金刚石。探测器家族里最近十年蹿起来一匹黑马:人工合成的单晶金刚石探测器。它本质上是固体电离室,中子打在碳上,沉积能量直接变成电脉冲,能量分辨率比闪烁体高一截,耐高温、抗辐照、计数率上限高。物理上它有个漂亮的巧合:2.5 MeV中子在金刚石里只能走碳核弹性散射这一条路,最大沉积能量约0.71 MeV;14 MeV中子则能走¹²C(n,α)⁹Be反应,在约8.3 MeV处打出一个干净的高斯峰。两个能区之间基本是空的——设一个2 MeV的阈值,DD和DT中子就分家了。米兰CNR等离子体科学与技术研究所的团队(D. Rigamonti、A. Dal Molin、M. Rebai、M. Tardocchi、M. Nocente等,联合Uppsala大学、UKAEA、欧盟联合研究中心)把一套金刚石探测器装上了JET的准直视线,DTE2前完成标定,成果发表在Nuclear Fusion 64, 016016(2024)。这篇论文里最值钱的两个结果:一是低氚浓度炮(氚占比小于3%)里,裂变电离室的"总产额"实际上混进了大量DD中子,而金刚石凭谱分离直接数清了两种中子——裂变室对中子能量"一视同仁",恰恰是它做不了这件事;二是首次在中性束加热的DT等离子体里做出高分辨14 MeV谱学,把热发射和非热发射分开,量出热成分占21%左右——这个数直接关系α自持加热能不能接上班。顺带一提,这套探测器的晶体由Element Six公司供应,表征工作里也有北京大学装置的一份参与。
第四件事,也是物理上最戏剧性的一件:NIF的超热离子。NIF在2022年12月实现点火之后,手里攒下了一批产额极高的DT炮,配套的是五台石英切伦科夫中子飞行时间谱仪——测量精度达到万分之几,中子速度的不确定度只有每秒5公里。热平衡的麦克斯韦分布有一个铁律:中子能谱的均值和方差(对应离子平均动能和离子温度)之间有固定的关系。LLNL的团队用这个关系检查历年数据,发现一旦等离子体开始燃烧、α加热显著起来,DT中子的平均能量就超出这条关系线——离子平均动能比按温度算出来的高。也就是说,真正发生聚变的那部分氘氚离子,不是老老实实的麦克斯韦分布,里面有超热成分。这个结果2016年就露过苗头(M. Gatu Johnson,MIT等离子体科学与聚变中心,PRE 94, 021202:分层DT内爆中DT离子温度系统高于DD离子温度,一维模型解释不了),但真正把案子做实的是2023年1月的Nature Physics(E. P. Hartouni、A. S. Moore等,Nat. Phys. 19, 72-77, 2023,DOI 10.1038/s41567-022-01809-3)。注意这篇论文的第一作者是Hartouni——媒体报道多把Moore放在最前面,正式引用要以Hartouni打头。作者阵容横跨LLNL(二十余人)、帝国理工学院(A. J. Crilly、B. D. Appelbe,理论)、MIT(M. Gatu-Johnson)、通用原子公司(J. Kilkenny)和洛斯阿拉莫斯(K. D. Meaney)。它的含义相当深:α加热不只是把等离子体"烧热",它可能把燃料离子的速度分布本身烧变形了——辐射流体力学程序里那套"分布永远是麦克斯韦"的假设,在燃烧等离子体里开始失效。对聚变能来说这是坏消息也是好消息:坏消息是反应率模型要重算,好消息是这等于免费送了一台"非平衡诊断仪"——中子谱的偏差本身就是α加热强度的直接读数。
把四件事放一起看,燃烧等离子体中子诊断的当前格局很清楚:产额(两套独立系统差34-41%没人能解释)、燃料比(2.5 MeV信号被墙的倍增和散射污染,得换TT连续谱)、成分分离(金刚石提供了干净的DD/DT分家方案)、分布函数(NIF发现麦克斯韦假设失效)。四条线指向同一个结论:燃烧等离子体不是"更热的等离子体",它是一种物理上新的对象——α加热主导、离子分布可以偏离平衡、墙也会参加中子学。JET的DTE2/DTE3和NIF的点火炮是人类仅有的两个拿到过这种对象的实验,它们的诊断都还在跟自己的系统误差搏斗,而且两边的当务之急出奇地一致:把独立测量的账对齐。
中国的位置说实话还在这张地图的边缘。JET那套金刚石探测器里,北京大学参与了晶体表征环节,但探测器的研制主导权在米兰CNR和罗马的组手里;国内EAST是全超导托卡马克,目前只能做氘等离子体,DT级的中子诊断无从谈起;合肥在建的紧凑型装置BEST计划做燃烧等离子体实验,届时中子产额、燃料比、快离子分布这一整套诊断要现学现用——这是未来几年国内聚变诊断界绕不开的一张考卷。好在方向是清楚的:金刚石探测器、切伦科夫nToF、飞行时间谱仪这三样,都是投入可控、可以追赶的硬件;真正难追赶的是JET用四十年、两轮DT实验攒下的那套标定与系统误差控制文化。这可能是比任何单台探测器都贵的资产。
---
核实与纠错记录(本次调研中修正的常见误记)
- "(n,2n)倍增20%"是数字张冠李戴。Batistoni等FED 88, 733 (2013)的原文有两件事:(n,2n)中子倍增在DD运行下可忽略、在DT运行下评估为大于25%;2.5 MeV中子通量因铍弹性散射在内辐照端下降约20%。流传的"20%倍增"把两个数搅在了一起。
- Nature Physics超热离子论文的正式第一作者是E. P. Hartouni(LLNL),不是A. S. Moore(LLNL新闻稿以Moore为叙述主角导致误传)。引用格式:Hartouni, E.P., Moore, A.S., Crilly, A.J. et al., Nat. Phys. 19, 72-77 (2023)。在线发表于2022年11月14日,刊期是2023年1月(vol. 19, issue 1)。
- RSI 96, 063501 (2025)是JET中子相机的DT产额论文(Hägg等);它的DD前作是RSI 94, 073502 (2023),勿混。前作结论是相机与裂变室平均偏差10%以内,DT版结论是偏差34%/41%且原因未明——"原因已查明"的说法是错的,作者明确写了" presently unknown"。
- JET金刚石探测器论文的正式发表信息是Nucl. Fusion 64, 016016 (2024),DOI 10.1088/1741-4326/ad0a49,在线发表于2023年11月23日。第一作者D. Rigamonti,通讯作者A. Dal Molin,单位是意大利国家研究委员会等离子体科学与技术研究所(CNR-ISTP,米兰),不是罗马La Sapienza(La Sapienza/罗马第二大学组是金刚石晶体电子学接触的长期供应方,但本篇署名机构为CNR-ISTP)。
- 裂变电离室对2.5与14 MeV中子响应近乎平坦,因此无法区分DD与DT贡献——这是Rigamonti论文明确给出的低氚炮失效机理,不是探测器的"缺陷"而是原理限制。
- JET纯氘等离子体中氚燃耗贡献的DT中子约占DD中子产额的1%(Villari等FED 217, 115133, 2025),不是百分之几的量级模糊说法;DTE3纪录脉冲为#104522,2.46×10^19中子、峰值5.06×10^18 n/s。
---
参考文献
- L. Hägg, S. Conroy, G. Ericsson, Z. Ghani(UKAEA卡勒姆), L. Giacomelli(欧盟联合研究中心JRC, Ispra), D. Marocco, M. Riva(ENEA弗拉斯卡蒂), A. Milocco(米兰比可卡大学), B. Vingren, G. Robertson, J. Eriksson, E. Andersson Sundén(乌普萨拉大学应用核物理), JET Contributors, *Estimating the neutron yield in a deuterium–tritium plasma with the JET neutron camera*, Rev. Sci. Instrum. 96 (2025) 063501; DOI 10.1063/5.0231639
- L. Hägg, S. Conroy, G. Ericsson等(乌普萨拉大学/UKAEA/ENEA/米兰比可卡大学), *Estimating the neutron yield in a deuterium plasma with the JET neutron camera*, Rev. Sci. Instrum. 94 (2023) 073502; DOI 10.1063/5.0144654
- P. Batistoni(EFDA-JET卡勒姆/ENEA弗拉斯卡蒂), D. Barbier(EFDA CSU卡勒姆), J. Likonen(VTT, 芬兰), JET-EFDA Contributors, *Fusion technology activities at JET in support of the ITER program*, Fusion Eng. Des. 88 (2013) 733-737; DOI 10.1016/j.fusengdes.2012.12.023((n,2n)倍增DD可忽略/DT大于25%、2.5 MeV通量散射削20%的原始出处)
- E.P. Hartouni, A.S. Moore, A.J. Crilly, B.D. Appelbe(帝国理工学院布莱克特实验室惯性聚变研究中心), P.A. Amendt, K.L. Baker, D.T. Casey, D.S. Clark, T. Döppner, M.J. Eckart, J.E. Field, M. Gatu-Johnson(MIT等离子体科学与聚变中心), G.P. Grim, R. Hatarik, J. Jeet, S.M. Kerr, J. Kilkenny(通用原子公司), A.L. Kritcher, K.D. Meaney(洛斯阿拉莫斯国家实验室), J.L. Milovich, D.H. Munro, R.C. Nora, A.E. Pak, J.E. Ralph, H.F. Robey, J.S. Ross, D.J. Schlossberg, S.M. Sepke, B.K. Spears, C.V. Young, A.B. Zylstra(劳伦斯利弗莫尔国家实验室), *Evidence for suprathermal ion distribution in burning plasmas*, Nat. Phys. 19 (2023) 72-77; DOI 10.1038/s41567-022-01809-3
- M. Gatu Johnson(MIT等离子体科学与聚变中心), J.P. Knauer(罗切斯特大学激光能量学实验室), C.J. Cerjan, M.J. Eckart, G.P. Grim, E.P. Hartouni, R. Hatarik, D.T. Casey, T. Döppner, O. Hurricane等(劳伦斯利弗莫尔国家实验室), J.D. Kilkenny(通用原子公司), *Indications of flow near maximum compression in layered deuterium-tritium implosions at the National Ignition Facility*, Phys. Rev. E 94 (2016) 021202(R); DOI 10.1103/PhysRevE.94.021202
- D. Rigamonti, A. Dal Molin(通讯作者), A. Muraro, M. Rebai, G. Gorini, M. Nocente, E. Perelli Cippo, M. Tardocchi(意大利国家研究委员会等离子体科学与技术研究所CNR-ISTP, 米兰;Gorini与Nocente兼米兰比可卡大学), L. Giacomelli(欧盟联合研究中心JRC, Ispra), S. Conroy, G. Ericsson, J. Eriksson(乌普萨拉大学), V. Kiptily, Z. Ghani, Ž. Štancar(英国原子能管理局UKAEA卡勒姆), JET Contributors, *The single crystal diamond-based diagnostic suite of the JET tokamak for 14 MeV neutron counting and spectroscopy measurements in DT plasmas*, Nucl. Fusion 64 (2024) 016016; DOI 10.1088/1741-4326/ad0a49
- G. Ericsson, S. Conroy, M. Gatu Johnson, E. Andersson Sundén, M. Cecconello, J. Eriksson, C. Hellesen, S. Sangaroon, M. Weiszflog(乌普萨拉大学), *Neutron spectroscopy as a fuel ion ratio diagnostic: Lessons from JET and prospects for ITER*, Rev. Sci. Instrum. 81 (2010) 10D324; DOI 10.1063/1.3502323
- B. Eriksson, A. Hjalmarsson, S. Conroy, G. Ericsson, J. Eriksson(乌普萨拉大学), JET Contributors, *Determining the fuel ion ratio for D(T) and T(D) plasmas at JET using neutron time-of-flight spectrometry*, Plasma Phys. Control. Fusion 64 (2022) 055008; DOI 10.1088/1361-6587/ac5a0d
- R. Villari, N. Fonnesu, P. Batistoni等(ENEA弗拉斯卡蒂/UKAEA/JRC Ispra等, 欧盟多机构合作), X. Litaudon, J. Mailloux, Z. Ghani等, *Overview of deuterium-tritium nuclear operations at JET*, Fusion Eng. Des. 217 (2025) 115133; DOI 10.1016/j.fusengdes.2025.115133
- C.F. Maggi(马普等离子体物理研究所/UKAEA)等, *Overview of T and D-T results in JET with ITER-like wall*, Nucl. Fusion 64 (2024) 112012; DOI 10.1088/1741-4326/ad3e16