APTR质子同步加速器RFQ直线注入器的优化设计

为实现质子治疗装置的国产化和小型化,基于已完成安装调试的上海先进质子治疗装置(APTR),开展质子治疗注入器系统的升级设计研究,利用PARMTEQM设计软件和快聚束策略,针对APTR同步加速器RFQ直线注入器进行动力学设计模拟。RFQ工作频率为325 MHz,流强18 mA,对从离子源引出的低能质子束流进行匹配俘获、横向聚焦、纵向聚束和预加速,引出能量为3.0 MeV。通过优化预注入器RFQ动力学设计方案和极头参数,有效避免参数共振,减小束流损失,使其整体传输效率达到98.0%,在水平和垂直方向上的发射度增长分别为1.2%和3.3%,出口束流满足下一级腔体的注入需求,开展设计模拟验证和相关冗余度分析,为质子同步加速器的治疗设备和直线注入系统提供参照依据。     

小型医用质子同步加速器的慢引出设计

 围绕北京医用质子同步加速器的方案, 研究设计了1/3整数共振的慢引出系统。利用该系统, 可获得较高的引出效率(大于95% )、较小的引出束流发射度(在静电切割器处, 积分发射度为0.15πmm·mrad)、较窄的引出动量分散(0.12% ) 以及较均匀的引出束流。     

LIA注入器发射度增长机理研究及其抑制措施

闪光X光照相是强流直线感应加速器（LIA）目前最重要的应用之一。为了获得更小的焦斑尺寸以达到精密照相的目的，减小LIA的发射度成为主要的技术途径之一，因此首先要求作为LIA束源的注入器能提供低发射度的带电粒子束。本文从非线性作用力导致束剖面电流密度分布变化的角度对注入器发射度的演变进行研究，根据研究揭示的规律提出相应的抑制发射度增长的措施。     

准周期聚焦–加速耦合系统的均温设计方法

强流质子直线加速器要求严格控制束流损失和束流发射度增长.理论已经证明,强的空间电荷作用在不同自由度之间的耦合,会因为其间的温度差异,通过束流的相干不稳定性,使束流发射度增长.因此,有必要按照均温的原则设计强流加速器.但是,由于质子直线加速器的各种加速结构均为准周期耦合系统,使得均温设计十分繁琐而难以达到完全均温.我们利用国际上通用的束流动力学软件TRACE3-D,给它补充了均温设计功能,通过与PARMILA程序的配合使用,可以方便地在加速器设计中实现均温条件.本文将介绍我们对TRACE3-D的修改补充,并以强流质子直线加速器设计实例,说明均温设计的必要性.     

电扫描型单缝单丝强流离子束发射度测量仪

 介绍了一台基于电偏转扫描方法的强流离子束发射度仪，讨论了发射度仪硬件的系统构成和设计要点，该发射度仪采用阶梯型电压对通过前缝的束流进行扫描来获得束流散角信息。并对测量计算软件和数据处理方法进行了介绍，采用了先扣除本底，再设置阈值的方法进行数据处理。处理后的数据利用程序可获得均方根和边界发射度、束斑大小、最大正负散角以及发射度相图。利用该发射度仪对ECR离子源引出的强流质子束的发射度进行了测量，并对测量数据进行了分析。对于引出电压为50 kV、脉冲重复频率为166 Hz、脉冲宽度为1 ms、平均束流强度为4 mA的质子脉冲束，其归一化均方根发射度为0.27 pmm·mrad。     

质子照相中成像磁透镜系统对次级质子的消除

 为了解成像磁透镜系统使用对提高质子照相品质的作用,研究了成像磁透镜系统对核反应产生的次级质子的消除作用。利用MCNPX模拟了次级质子的能谱,结果表明：次级质子的平均能量不足照相质子束的50%。从次级质子的能量特性和成像系统质子运动学方面的分析表明：成像磁透镜系统可消除次级质子的影响。利用程序MCNPX对含成像磁透镜系统的质子照相进行了模拟,结果显示：成像磁透镜系统将次级质子对光程的影响减小到2%以下,显著提高了图像品质。     

高能质子照相的成像方程

与X射线闪光照相相比,GeV级质子照相因质子的穿透能力更强而更能满足流体动力学最佳光程、更小不确定度的要求。质子照相用于提取阴影客体信息的过程比X射线闪光照相复杂。为了便于分析质子照相过程及其特性,以质子照相的输运过程为基础,逐步考虑了质子与物质相互作用的各种过程——核反应、多次库仑散射、准弹性散射和弹性散射,并获得含各种过程的质子照相理论计算表达式。在含准直的系统中,理论计算公式中包含了多次弹性散射和准弹性散射项。通过与QMCPrad计算结果的对比得知弹性散射高斯近似下的理论计算结果与蒙特卡罗方法(MC)计算结果一致,与更准确的"黑体"模型的相对误差为3%~5%。这意味着在某种程序上,当研究质子照相规律或设计系统时,高能质子照相的理论计算公式能够替代MC模拟。     

LIA注入器发射度增长机理研究及其抑制措施

根据直线感应加速器(LIA)的特点采用了轴对称直流模型,结合E.P.Lee的理论导出了发射度增长的微分方程,提出了一种直线感应加速器发射度增长的机理,根据该机理计算了从均匀分布到高斯分布的演变过程中发射度增长的大小,指出直线感应加速器(LIA)注入器阶段束流发射度增长主要是由于非线性效应导致的束刨面电流密度分布变化引起的.根据研究结果,提出了相应的抑制发射度增长的措施.     

Weibel不稳定性自生电磁场对探针质子束的偏转作用研究

Weibel不稳定性的自生电磁场对于等离子体能量输运、无碰撞冲击波形成等物理过程具有关键的影响.实验上往往采用质子束照相来诊断其电磁场结构.一般认为,探针质子束的轨迹偏转主要来自于磁场,而自生电场的作用被认为可忽略不计.本文利用三维粒子模拟程序研究了典型参数下的Weibel不稳定性发展过程,并使用径迹追踪法评估了Weibel不稳定性的质子束照相过程中电场和磁场对探针质子束的偏转作用.对比分析发现,引起探针质子束偏转的主要因素并不是磁场,而是过去研究中常被忽略的电场.主要原因为:Weibel不稳定性的自生磁场往往成管状结构,在使用探针质子束对其进行侧向照相时,磁场的作用会被自身中和并抵消.该认识将有助于深入理解Weibel不稳定性质子照相的实验结果.     

ADS强流质子加速器低能传输段发射度匹配研究北大核心CSCD

研制了一套强流质子源及低能传输线(LEBT)注入器用于ADS质子直线加速器。质子源产生35 ke V强流束经过低能传输段聚焦进入射频四极(RFQ)入口。低能传输段不匹配是强流RFQ中引起束流丢失的主要原因。不同加速段的束流匹配是减少束流损失与抑制发射度增长的重要手段。束流损失导致RFQ电极表面受热变形进而引起高频打火,降低RFQ长期运行的稳定性。针对以上问题,研究LEBT发射度在不同的实验条件下如何实现加速器更好的匹配。研究结果表明,LEBT出口束流在35 ke V,10 m A下,束流发射度小于0.2πmm·mrad,当LEBT螺线管电流为210和270 A时,束流在RFQ入口满足匹配条件。     

粒子束团状态对测量束流发射度及能量的影响

为了验证国产质子注入器的参数是否满足需求,注入器团队设计了束流测量系统用于测量国产质子直线注入器束流的流强、发射度、能量以及能散等关键指标。此测量系统包含了采用变聚焦法测量发散度、采用分析磁铁测量束流能量和能散的主要功能。利用束流输运线设计软件Tracewin(版本2.11.4.1)进行了系统束线的物理设计,对束测系统测量质子束流的发散度和能量的精度进行了模拟计算。由于经过RFQ-(APF)DTL加速后的粒子束团为“拖尾”的非理想粒子束团,需要针对非理想束团对束测系统测量发射度和能量产生的影响进行分析。通过对模拟计算结果的分析,发现相对于测量理想粒子束团的结果非理想粒子束团对束测系统测量发射度精度影响较大;非理想粒子束团对束测系统测量能量精度影响较小。     

ATPF–a dedicated proton therapy facility

A proton therapy facility based on a linac injector and a slow-cycling synchrotron is proposed. To obtain good treatments for different cancer types, both the spot scanning method and the double-scattering method are adopted in the facility, whereas the nozzles include both gantry and fixed beam types. The proton accelerator chain includes a synchrotron of 250 MeV in maximum energy, an injector of 7 MeV consisting of an RFQ and a DTL linac, with a repetition rate of 0.5 Hz. The slow extraction using the third-order resonance and together with the RFKO method is considered to be a good method to obtain a stable and more-or-less homogenous beam spill. To benefit the spot scanning method, the extraction energy can be as many as about 200 between 60 MeV and 230 MeV. A new method – the emittance balancing technique of using a solenoid or a quadrupole rotator is proposed to solve the problem of unequal emittance in the two transverse planes with a beam slowly extracted from a synchrotron. The facility has been designed to keep the potential to be upgraded to include the carbon therapy in the future.     

强流ECR离子源引出系统研究

为了提高强流ECR 离子源的引出束流品质,分别设计了1# 和2# 引出系统,利用束流引出模拟软件PBGUNS 对1# 和2# 引出系统进行了质子束流引出与传输的模拟计算,结合实际测得的发射度数据分析引出系统,发现2# 引出系统比1# 引出系统引出束流品质高。对ECR 离子源引出系统的电势等位线分布等参数引起的球差进行了简单数学推导及MATLAB 绘图,并结合1# 和2# 引出系统束流相图模拟结果证明了球差会使引出束流品质有效发射度增长,通过适当加大电极孔径可改善束流聚焦情况,得到了束流光学聚焦较好的束流引出系统设计。To improve the quality of extracted ion beam from a high current ECR ion source, 1# and 2# extraction systems were designed and tested. The PBGUNS code was used to simulate the 1# and 2# extraction systems of proton ion beam. The emittance measurement results with the two different extraction systems were compared and analyzed with the simulation, the conclusion that more high quality beam extracted from 2# system than 1# system was got. The formula derivation of ECR ion source extraction system spherical aberration and MATLAB drawing was done by the analyzing on the distribution of extraction field equipotentials, effective emittance increasing caused by spherical berration was proved by 1# and 2# extraction systems beam phase space simulation result, beam focusing would be improved if electrode hole size increasing appropriately and a general concept on good optics focusing of ion beam extraction system was proposed finally.     

影响X光源特性的参数研究

针对高能电子束撞击重金属靶产生轫致辐射光子, 利用蒙特卡罗方法详细研究了电子束半径与发射度以及靶的厚度对X光源特性的影响. 结果表明：电子束半径与发射度不仅是影响照射量的主要因素, 也是造成照射量分布不均匀的主要原因; 给定电子束, 存在一个靶厚度使得靶前1 m处的照射量最大. 因此, 在X射线成像系统的设计和模拟过程中, 应综合考虑电子束半径与发射度和靶厚的影响.     

质子照相中基于能量损失的密度重建

提出了一种质子能量在中高能时利用能量损失进行密度重建的方法,并利用Bethe-Bolch公式给出了利用能量损失进行密度重建的方程及条件.针对1.6 GeV的质子能量,通过定量计算常见材料的阻止本领,得出质子能量在1.45–1.6 GeV范围内时,材料的阻止本领的变化率小于1%,可近似为常数.最后,通过理论计算和Geant 4模拟,得出质子能量在1.6 GeV时,可以对面密度为113 g/cm²的缩比法国实验客体进行密度重建.     

Stripping extraction calculation and simulation for CYCIAE-100

A 100 MeV H- compact cyclotron is under construction at China Institute of Atomic Energy (CYCIAE-100). The proton beams of 75 MeV-100 MeV at an intensity of 200 μA will be extracted in dual opposite directions by charge exchange stripping devices. The crossing point at the switching magnet center is fixed inside the magnet yoke and the stripping points for various extraction energies are calculated by the code CYCTRS. With the code GOBLIN, we can calculate the transfer matrix including the dispersion effects from the stripping point to the switch magnet. The beam distribution just after stripping foil can be obtained from the multi-particle tracking code COMA and the extracted beam parameters after the switch magnet such as emittance, envelope, dispersion, energy spread, bunch length, etc. are given by the extraction orbit simulations.     

Transmission calculation by empirical numerical model and Monte Carlo simulation in high energy proton radiography of thick objects

An empirical numerical model that includes nuclear absorption, multiple Coulomb scattering and energy loss is presented for the calculation of transmission through thick objects in high energy proton radiography. In this numerical model the angular distributions are treated as Gaussians in the laboratory frame. A Monte Carlo program based on the Geant4 toolkit was developed and used for high energy proton radiography experiment simulations and verification of the empirical numerical model. The two models are used to calculate the transmission fraction of carbon and lead step-wedges in proton radiography at 24 Ge V/c, and to calculate radial transmission of the French Test Object in proton radiography at 24 Ge V/c with different angular cuts. It is shown that the results of the two models agree with each other, and an analysis of the slight differences is given.     

Analytical and simulation studies for diode and triode ion beam extraction systems

This work is concerned with ion beam dynamics and compares the emittance to aberration ratios of two-and three-electrode extraction systems. The study is conducted with the aid of Version 7 of SIMION 3D ray-tracing software. The beam dependence on various parameters of the extraction systems is studied and the numerical results lead to qualitative conclusions.     

Monoenergetic proton beams accelerated by a radiation pressure driven shock

We report on the acceleration of impurity-free quasimononenergetic proton beams from an initially gaseous hydrogen target driven by an intense infrared (λ=10 μm) laser. The front surface of the target was observed by optical probing to be driven forward by the radiation pressure of the laser. A proton beam of ～MeV energy was simultaneously recorded with narrow energy spread (σ～4%), low normalized emittance (～8 nm), and negligible background. The scaling of proton energy with the ratio of intensity over density (I/n) confirms that the acceleration is due to the radiation pressure driven shock.     

Laser Accelerated, High Quality Ion Beams

Intense beams of protons and heavy ions have been observed in ultra-intense laser-solid interaction experiments. Thereby, a considerable fraction of the laser energy is transferred to collimated beams of energetic ions (e.g. up to 50 MeV protons; 100 MeV fluorine), which makes these beams highly interesting for various applications. Experimental results indicate very short pulse duration and an excellent beam quality, leading to beam intensities in the TW range. To characterize the beam quality and its dependence on laser parameters and target conditions, we performed experiments at several high-power laser systems. We found a strong dependence on the target rear surface conditions allowing to tailor the ion beam by an appropriate target design. We also succeeded in the generation of heavy ion beams by suppressing the proton amount at the target surface. We will present recent experimental results demonstrating a transverse beam emittance far superior to accelerator-based ion beams. Finally, we will discuss the prospect of laser-accelerated ion beams as new diagnostics in laser-solid interaction experiements. Special fields of interest are proton radiography, electric field imaging, and relativistic electron transport inside the target.