基于Geant4的質子輸運計算的加速方法研究

姚志明，黑東煒，宋顧周，馬繼明，韓長材，段寶軍，岳志勤

（西北核技術研究所強脈沖輻射環境模擬與效應國家重點實驗室，西安 710024）

基于Geant4的質子輸運計算的加速方法研究

姚志明，黑東煒，宋顧周，馬繼明，韓長材，段寶軍，岳志勤

（西北核技術研究所強脈沖輻射環境模擬與效應國家重點實驗室，西安 710024）

質子在醫學、材料科學、輻射成像等領域有著廣闊的應用前景，質子輸運過程的數值模擬是質子應用理論分析的重要手段。本研究基于Geant4軟件編寫了蒙卡模擬程序，研究了程序的加速方法。結果表明，次級粒子截斷和并行計算可以有效加快程序的運行速度。程序易于擴展，為質子應用的快速模擬奠定了基礎。

質子輸運；蒙特卡羅；加速方法；次級粒子截斷；并行計算

質子在醫學、材料科學、輻射成像等領域有著廣泛的應用。目前，國內有多臺質子加速器已經建成或正在建設當中，它們有著不同的應用前景。例如，中國原子能科學研究院的回旋加速器上開展了很多國防核技術和核物理的研究工作［1］；上海應用物理研究所正在研制用于惡性腫瘤治療的質子治療裝置［2］；在建的中國散裂中子源的質子加速器計劃用于空間輻射效應模擬和質子照相；等等。

數值模擬計算可以為粒子的應用提供理論指導和依據。數值模擬中的蒙特卡羅方法通過隨機抽樣，能夠比較真實地反映粒子的輸運過程［3］。目前，模擬質子輸運的軟件包括SRIM、MCNP、MCNPX、Geant4等，其中Geant4允許用戶二次開發，容易根據實際需要對程序進行改進和擴展。

基于 Geant4的模擬質子輸運的研究比較豐富［4-6］，加速方法的研究卻較少。文獻［7，8］中給出了通用的并行計算加快模擬速度的方法。然而，針對質子輸運過程的特點對Geant4程序進行加速的方法未見報道。

本研究以10 GeV的質子在FTO （French Test Object）中的輸運過程為例，基于Geant4建立蒙卡模擬程序［9，10］，探討加快質子輸運模擬速度的方法。

1　質子與物質的作用規律

質子與物質的作用過程主要有電離能損、多次庫侖散射、與原子核的彈性碰撞和非彈性碰撞［11］。電離能損是質子與很多核外電子碰撞，每次碰撞電子獲得一部分能量，質子的能量隨之減小。多次庫侖散射是當質子與原子最近距離比原子半徑小時，受到原子核電場的作用，在庫侖力作用下運動方向發生改變。與原子核碰撞時質子受到核力的作用，按質子能量是否發生改變分為彈性和非彈性碰撞。彈性碰撞按照原子處于基態和激發態可細分為彈性碰撞和準彈性碰撞，質子的運動方向發生較大偏轉；非彈性碰撞是原子核俘獲質子，發生核反應，核反應過程中釋放次級粒子，導致入射質子的個數損失，并可能產生能量較低的次級質子。本文10 GeV質子在FTO輸運過程中以上四種作用都會發生，導致質子的能量損失、個數衰減和方向偏轉。其中非彈性碰撞產生的中子、γ射線等次級粒子對結果貢獻小，模擬耗時多，截斷后可以提高計算速度。

2　模擬程序的建立

Geant4軟件采用模塊化設計，由不同的模塊分別實現探測器幾何結構定義、粒子與物質相互作用過程描述、跟蹤并記錄粒子的狀態和改變、中止次級粒子的產生等［10］。圖1給出了軟件的一些常用的類、子函數以及它們的功能。

模擬程序的幾何布局如圖2所示。Construct函數中定義物體的材料和幾何結構。物體選為FTO［11］，它是高能閃光照相中使用的一種通用測試對象，是由鎢、銅和塑料組成的多層同心球，內半徑分別為1 cm、4.5 cm、6.5 cm，塑料的外半徑為22.5 cm。物體后1 mm處放置1 mm厚的足夠大長方體，定義為敏感探測器。Geant4中的物理過程比較豐富，采用的物理模型也有多種選擇，為方便用戶的使用，Geant4開發人員將物理模型打包，定義了多種physic list，其中包含了質子與物質作用的全部物理過程和產生的中子、伽瑪等次級粒子的作用過程。本文選取了physic list中的QGSP_BERT物理過程列表。GeneratePrimaries函數中設置發射粒子的參數信息。粒子類型為質子，初始能量為 10 GeV，發射方向沿 z軸正方向。ProcessHits、Print、EndOfEvent函數中定義敏感探測器的功能，用于粒子信息的記錄。當粒子打到敏感探測器上，在ProcessHits函數中將該粒子的出射位置、運動方向、類型、動能和輸運編號存儲在一個Hit中。在Print函數中定義Hit中各信息的輸出格式。當發射的一個質子輸運完成后，EndOfEvent函數調用Print函數將各信息存儲在文本文件中，供分析使用。

3　加速方法與結果分析

3.1次級粒子截斷

質子輸運過程中發生核反應會產生大量中子、γ射線和低能量的質子。除低能量的質子外，還會產生一定數量的電子、μ介子、π介子等其他帶電粒子。這些次級粒子的輸運過程模擬會耗費大量時間。Geant4中的G4VUserStackingAction類可以對粒子的信息進行判斷，實現次級粒子的截斷，達到節省計算時間的目的。然而，次級粒子的截斷可能會影響質子輸運的模擬，本文選取了敏感探測器記錄到的質子總個數和直穿質子個數兩個參數來評估截斷對模擬準確性的影響。

Geant4軟件運行在聯想ThinkCentre M8380t型號的計算機上，安裝了windows7 32位操作系統，CPU頻率為3.4 GHz。

首先，不截斷次級粒子，程序將模擬所有直穿質子和次級粒子的輸運過程。發射105個質子，時間由G4Timer類記錄，模擬用時2 538 s。表1給出了敏感探測器記錄的粒子種類和個數。直穿質子和次級質子，可以由輸運編號區分，輸運編號為1的是直穿質子，大于1的是次級質子。

表1　不截斷產生的粒子個數Tab.1 Particle generation without cut

圖1　Geant4常用類和函數結構圖Fig.1 Common class and function of Geant4

圖2　模擬程序的幾何布局示意圖Fig.2 Geometry of the simulation program

然后，在ClassifyNewTrack函數中分別實現γ射線、中子和其他帶電粒子的截斷，并考慮了只模擬質子和只模擬直穿質子的情況。以截斷γ射線為例，ClassifyNewTrack函數中的程序代碼如下：

表2給出了各情況下模擬所需的時間，通過與表1比較，給出了直穿質子個數和總質子個數減少的百分比，并計算了截斷前后時間的比值，即加速比。

由表2可以得出以下結論：第一，截斷次級粒子可以有效地減少計算時間，截斷不同的粒子時加速比均達到了2倍以上。第二，截斷不同的粒子后對直穿質子的個數的影響均小于1%，可以認為不會影響直穿質子的模擬結果。第三，截斷γ射線對總質子個數的影響小于1%，截斷中子和其他帶電粒子后總質子個數減少了10%以上，分析原因是由于γ射線與物質作用后不會產生質子，截斷不會影響質子的模擬結果，而中子和一些帶電粒子與物質作用會進一步產生次級質子，截斷這些粒子會使得出射的次級質子有所減少，在需要考慮次級質子輸運過程的應用中不能采用。

表2　截斷后的計算速度測試結果Tab.2 Results of calculation speed testing with cut

3.2并行化加速方法

有關Geant4程序并行化的研究已較為豐富［7，8］。本研究將并行化方法應用于質子輸運模擬，采用文獻［8］中的方法對程序進行了并行化改造。并行計算在銀河計算機上進行，該計算機安裝了紅帽公司的Linux服務器操作系統，CPU頻率為1.6 GHz。表3給出了并行化后的測試結果。需要說明的是，當模擬的質子個數較多時，如果將所有質子的信息存儲在一個文件中，得到的數據文件將會很大，不便于后續數據的統計分析。在print函數中可以首先讀取CPU的編號，用編號為輸出文件命名，將不同CPU的計算結果存儲在不同的文件當中，解決數據文件過大的問題。實現語句如下：

由表3可以看出，加速比隨CPU個數的增加近似成線性增長，與文獻［8］中的結果一致。與CPU為3.4 GHz的聯想計算機的結果比較，在只模擬直穿質子的條件下，聯想計算機用時65 s，銀河計算機用時295 s。銀河計算機單個CPU運行的速率不及CPU頻率較高的聯想計算機，然而，銀河計算機上最多可提供的CPU數達1 024個，可以根據需要進一步增加CPU的使用個數來提高質子輸運的模擬速度。

表3　并行程序的計算速度測試結果Tab.3 Results of calculation speed testing of parallel program

4　結論

本研究基于Geant4軟件建立了模擬10 GeV能量的質子穿透FTO過程的蒙卡程序。重點探討了程序的加速方法，在程序中分別實現了次級粒子截斷和并行計算。截斷γ、中子和其他帶電粒子后，加速比均達到了2倍以上，截斷中子和一些帶電粒子會使透射的次級質子個數有所減少，當不需要考慮次級質子輸運過程時可以采用。質子輸運模擬程序并行化后的加速比隨CPU使用個數的增加呈線性增長，與文獻結果一致。與單個CPU不采用粒子截斷的情況相比，只記錄初級質子的16個CPU并性程序加速比可達610倍。

［1]張天爵，樊明武，李振國，等.強流回旋加速器技術的研究進展［J］.科學通報，2009，54（15）：2169-2176.

［2]楊朝霞，李德明，張滿洲.基于Geant4模擬的質子治療束配系統的束流光學設計［J］.核技術，2013，36（7）：070201.

［3]許淑艷.蒙特卡羅方法在實驗核物理中的應用［M］.北京：原子能出版社，2006.

［4]張慧，王振超，關世榮，等.質子能量歧離的Geant4模擬［J］.黑龍江科學，2014，1（5）：14-16.

［5]路偉.Geant4的開發及其在空間輻射效應分析中的應用［D］.長沙:國防科學技術大學研究生院，2007：13-15.

［6]聶鵬煊，汪一夫，鄭濤，等.CsI（TI）對高能質子能量響應的蒙特卡羅研究［J］.武漢科技學院學報.2010，23（2）:27-32.

［7]王力，李煜，岳騫，等.基于Geant4的射線探測仿真并行化研究［J］.核電子學與探測技術，2011，11（11）：1237-1240.

［8]謝紅剛，朱金輝，黃流興.Geant4在高性能計算機上的并行［C］//第十屆全國蒙特卡羅方法及其應用學術會論文：148-151.

［9]張慧，梅雪松，關世榮，等.4He離子盧瑟福背散射的Geant4模擬［J］.黑龍江科學，2013，1（4）：22-24.

［10]Geant4 user's guide for application developers［EB10L］.（292-11-30）［2015-03-06］http://geant4.web.cern.ch/geant4/UserDocumentation/ UserGuides/ForApplication Developer/html/index.html.

［11]ZIOCK H J，ADAMS K J，ALRICK K R，et al.The Proton Radiography Concept［R］.Los Alamos National Report，LA-UR-98-1368，1998：1-17.

Research on Accelerating Method of Proton Transportation Simulation Program Based on Geant4

YAOZhi-ming，HEI Dong-wei，SONGGu-zhou，MAJi-ming，HANChang-cai，DUANBao-jun，YUE Zhi-qin
（State Key Laboratory of Intense Pulsed Radiation Simulation and Effect，Northwest Institute ofNuclear Technology，Xi'an 710024，China）

Proton has bright future in the field of medical，material analysis and radiography.Numerical simulation provides a good way to analyze the process protons transmitting through materials.This article establishes a Monte Carlo program based on Geant4.The accelerating method is introduced.The result implies that secondary particles cut and program paralleled can shorten the time cost effectively.The program can be easily extended which provides a good foundation for fast simulation ofproton transportation.

Proton transportation；Monte Carlo；Acceleration method；Secondaryparticle cut；Parallel program

TL99

A

1674-8646（2015）07-0007-03

2015-04-09

超快物理過程射線診斷及應用研究（51311060502）

姚志明（1989-），男，吉林輝南人，碩士研究生，助理工程師，主要從事質子掃描成像的原理探究。