林 輝，謝 聰，張擁軍,2，熊楨宇，吳東升,*，曹瑞芬，F(xiàn)DS團(tuán)隊(duì)

(1.合肥工業(yè)大學(xué) 電子科學(xué)與應(yīng)用物理學(xué)院，安徽 合肥 230009；2.蚌埠醫(yī)學(xué)院 公共課程部，安徽 蚌埠 233000；3.中國科學(xué)院 核能安全技術(shù)研究所，安徽 合肥 230031)

?

Geant4不同物理模型對放療質(zhì)子束模擬的影響

林 輝1，謝 聰1，張擁軍1,2，熊楨宇1，吳東升1,*，曹瑞芬3，F(xiàn)DS團(tuán)隊(duì)3

(1.合肥工業(yè)大學(xué) 電子科學(xué)與應(yīng)用物理學(xué)院，安徽 合肥 230009；2.蚌埠醫(yī)學(xué)院 公共課程部，安徽 蚌埠 233000；3.中國科學(xué)院 核能安全技術(shù)研究所，安徽 合肥 230031)

Monte Carlo方法常用于質(zhì)子束劑量模擬以優(yōu)化質(zhì)子束臨床治療，合適的物理模型選擇關(guān)系到Monte Carlo模擬結(jié)果的可靠性及模擬效率。本工作選取9個(gè)可用于質(zhì)子模擬的Geant4物理模型，計(jì)算比較了質(zhì)子束在多種介質(zhì)中的劑量及次級粒子產(chǎn)額。結(jié)果表明：標(biāo)準(zhǔn)電磁物理模型、低能電磁Penelope模型和Livermore模型雖在模擬質(zhì)子整體劑量分布上可接受，但在微觀上缺失重粒子生成。參數(shù)化驅(qū)動模型LHEP的模擬時(shí)間最短,但與QGSP相同，均不能產(chǎn)生復(fù)雜重離子。QGSP_BIC_EMY模型較其他模型明顯低估Bragg峰-入口劑量比和絕對劑量。QGSP_BERT模型、QGSP_BIC模型和FTFP_BERT模型較適合放療質(zhì)子束模擬。

物理模型；質(zhì)子放療；微觀產(chǎn)額；劑量特征量

質(zhì)子放療是一新興的腫瘤物理治療技術(shù)。質(zhì)子劑量曲線具有尖銳Bragg峰，在穿越路徑時(shí)僅釋放出少許能量，而達(dá)到治療深度時(shí)釋放出大量能量[1]。放療質(zhì)子的能量范圍通常為60～250 MeV，可治療從淺表層腫瘤(如眼部腫瘤)到深度腫瘤[2]。

蒙特卡羅(MC)方法通過建立求解問題的概率模型，對模型的多次抽樣給出計(jì)算問題的近似統(tǒng)計(jì)解[3]。MC程序可利用機(jī)器隨機(jī)數(shù)模擬粒子輸運(yùn)發(fā)生的復(fù)雜物理過程，給出粒子出射的角分布特征[4]、非均勻介質(zhì)中精確的能量分布[5]以及微觀分析關(guān)鍵物理作用的權(quán)重比例[6]等。但MC方法不同于實(shí)驗(yàn)測量，需基于預(yù)先知道的理論模型對粒子輸運(yùn)作用和方向等進(jìn)行判斷。

目前模擬質(zhì)子輸運(yùn)的程序有MCNPX、Fluka、Geant4[7]等。其中開源程序平臺Geant4包含了描述粒子與物質(zhì)相互作用的多種物理模型，其面向?qū)ο蟮哪K化編程設(shè)計(jì)增強(qiáng)了程序的可移植性。關(guān)鍵的問題是使用者需了解哪些物理模型適合所研究問題[8]。錯(cuò)失粒子類型或物理作用均可造成結(jié)果偏差，但過于細(xì)致的物理模型又會造成模擬速度的降低。文獻(xiàn)[2,5,6,8]比較了不同物理模型對質(zhì)子總沉積劑量的影響，但對不同物理模型對質(zhì)子微觀產(chǎn)額的影響缺乏研究。

本工作通過開發(fā)Geant4用戶程序，基于Geant4類庫提供的9種物理模型模擬放療質(zhì)子束在多種介質(zhì)中的劑量及次級粒子微觀產(chǎn)額，擬找到平衡精度和速度適合模擬治療質(zhì)子束能量范圍的物理模型。

1 材料和方法

1.1 程序發(fā)展

本用戶程序基于Geant4程序例extendedelectromagneticTestEm3進(jìn)行了二次開發(fā)，借助其內(nèi)外環(huán)套(沿y-z面)的層狀(沿x方向)幾何設(shè)計(jì)，模擬質(zhì)子束照射模體(圖1)。Geant4程序sourcephysics_listslists目錄下提供了常用的模型物理文件，特別是一些字符型模型可進(jìn)行強(qiáng)子輸運(yùn)模擬。若用戶需要還可制作自己的模型物理文件，所基于的初級模塊可在sourcephysics_listsuilderssrc目錄下查找。本工作采用Geant4_9_5_p02版本。

圖1 照射模體Fig.1 Configuration of irradiation

1.2 Geant4物理模型

Geant4程序平臺提供了多種粒子和物理作用模型。Geant4提供了玻色子、輕子、介子、重子、短壽命粒子和離子6類粒子的定義。程序中需考察的初級粒子和次級粒子類型必須在PhysicsList中進(jìn)行ConstructParticle()構(gòu)建。Geant4包含的物理作用有電磁作用、強(qiáng)子相互作用、輸運(yùn)、衰變、光學(xué)作用、光輕子-強(qiáng)子、特定體積參數(shù)化等。程序中需用到的作用類型必須在PhysicsList的ConstructProcess( )中構(gòu)建。

Geant4提供的電磁作用模型有標(biāo)準(zhǔn)電磁模型(standard EM，約1 keV～1 TeV)、低能電磁模型(250 eV～1 PeV)和光量子模型。標(biāo)準(zhǔn)電磁模型提供了G4EmStandardPhysics、G4EmStandardPhysics_option1、G4EmStandard-Physics_option2和G4EmStandardPhysics_option3 4個(gè)物理清單。G4EmStandardPhysics_option3考慮了質(zhì)子的各種電磁相互作用，適合醫(yī)學(xué)物理模擬。低能電磁模型LowE_Penelope和LowE_Livermore對于γ和e-/e+(電子/正電子)輸運(yùn)在0～1 GeV之間使用自己的模型，但在1 GeV～10 TeV之間仍使用standard EM。雖然標(biāo)準(zhǔn)和低能電磁模型未直接模擬強(qiáng)子作用過程，但均包含了一些質(zhì)子電離處理的替代模型，并加入次級粒子產(chǎn)生處理[8]。

Geant4提供以下強(qiáng)子相互作用模型。

1) 參數(shù)化驅(qū)動模型(PDM，parametrization driven model)。模型中僅強(qiáng)子-核子碰撞被詳細(xì)模擬，其他發(fā)生在核子內(nèi)的相互作用是通過產(chǎn)生額外的強(qiáng)子，并當(dāng)作次級粒子進(jìn)行模擬。而額外的強(qiáng)子數(shù)目、類型和分布是通過擬合實(shí)驗(yàn)數(shù)據(jù)獲取的函數(shù)決定的。因此，此模擬需依賴大量的參數(shù)獲取，以保證模型的合理性。目前有兩套PDM可用于模擬高能強(qiáng)子-核子相互作用：低能模型(1～25 GeV)和高能模型(25 GeV～10 TeV)。

2) 級聯(lián)模型(CM，cascade model)。模擬了基本粒子和次級粒子的級聯(lián)輸運(yùn)過程，包括BIC(binary cascade model(約<10 GeV))和BERT(Bertini cascade model(約<10 GeV))。但級聯(lián)模型不適合低能，因?yàn)榈湍苌婕罢舭l(fā)和預(yù)平衡衰變。僅當(dāng)粒子的動能高于閾值(75 MeV)時(shí)才能使用級聯(lián)模型。

3) 預(yù)復(fù)合核模型(PRECO，pre-compound model)。預(yù)復(fù)合核模型是強(qiáng)子動力學(xué)模型的擴(kuò)展，它將強(qiáng)子動力學(xué)模型適用范圍擴(kuò)展到核-核非彈性碰撞。考慮核反應(yīng)的預(yù)復(fù)合核階段直到非平衡態(tài)。從激發(fā)核里釋放核碎片或光子的模擬采用平衡模型。PRECO適用于0～150 MeV范圍內(nèi)的質(zhì)子、中子模擬。

此外，適合高能強(qiáng)子模擬的模型還有QGS(quark gluon string model(約>20 GeV))和FTF(Fritiof model(約>10 GeV))，它們與低能模型PRECO結(jié)合成為QGSP和FTFP才適合治療質(zhì)子束模擬。

本工作考察了以下物理模型對放療質(zhì)子劑量和微觀產(chǎn)額的影響。

模型1，G4EmStandardPhysics_option3；模型2，LowE_Penelope；模型3，LowE_Livermore；模型4，LHEP(低能與高能PDM)+模型1；模型5，QGSP(QGS+PRECO)+模型1；模型6，QGSP_BERT(QGS+PRECO+BERT)+模型1；模型7，QGSP_BIC(QGS+PRECO+BIC)+模型1；模型8，F(xiàn)TFP_BERT (FTF+PRECO+BERT)+模型1；模型9，QGSP_BIC_EMY(QGSP_BIC模型+EMY處理)+模型1。

其中，模型1～3為電磁模型。為考察強(qiáng)子模型的影響，模型4～9在含有標(biāo)準(zhǔn)電磁模型選項(xiàng)3(即模型1)之外又增添了強(qiáng)子模型處理。模型9為Geant4advancedhadrontherapy例提供的物理模型，在使用QGSP_BIC模型之外添加了G4EmExtraPhysics( )、G4Hadron-ElasticPhysics( )、G4QStoppingPhysics( )、G4IonBinaryCascadePhysics( )、G4Neutron-TrackingCut()處理。除模型1外，其他模型均采用了放射活度衰變離子物理處理。

1.3 模擬幾何

由于質(zhì)子側(cè)向散射較小，一般更關(guān)注其劑量曲線深度方向的分布。本模型采用一維長方體模型，模體大小為xcm×10 cm×10 cm，x為照射方向，可根據(jù)質(zhì)子束能量調(diào)節(jié)(圖1)。計(jì)數(shù)探測器沿x方向(即照射中心軸)平行排列，在垂直于照射軸面上的大小為5 cm×5 cm，沿照射方向根據(jù)部位不同，分別取非關(guān)鍵深度1 mm以及入口和Bragg深度處1 mm、0.1 mm或0.01 mm。照射野為5 cm×5 cm，垂直入射。各粒子的截?cái)嗑嚯x取1 mm。Geant4中截?cái)嗑嚯x對應(yīng)了產(chǎn)生δ電子和軔致輻射光子的射程閾值。Geant4中過小的射程截?cái)鄷@著增加模擬時(shí)間和計(jì)算結(jié)果的漲落[9]。Geant4advancedhadrontherapy例中相對于0.2 mm的網(wǎng)格厚度也采用了1 mm的粒子射程截?cái)?。為平衡精度和速度的矛盾，模擬前期進(jìn)行了截?cái)嘀颠x擇的驗(yàn)證工作，發(fā)現(xiàn)1 mm差異可接受。模體中分別充滿均勻介質(zhì)水、脂肪、肺、肌肉和骨(密度分別為1、0.99、0.30、1.04、1.49 g/cm3，數(shù)據(jù)來自于advancedhuman_phantom例)。每種情況跟蹤2×105個(gè)粒子，數(shù)據(jù)不確定度<1%。整個(gè)計(jì)算是在2.67 GHz CPU/2G內(nèi)存的計(jì)算機(jī)上進(jìn)行。

1.4 模型校驗(yàn)

2 結(jié)果與討論

2.1 PDD與次級粒子產(chǎn)額

為驗(yàn)證模擬的可靠性，模擬的PDD與文獻(xiàn)[10]數(shù)據(jù)進(jìn)行了對比(圖2)。從圖2可見，除物理模型9和文獻(xiàn)相差較大(約3.5%)外，其他模型幾乎完全一致(入口處(<5 cm)約1.0%，其他處約0.2%)。表1列出了各模型的次級粒子產(chǎn)額，其中模型1～3僅產(chǎn)生γ和e-，無重子產(chǎn)生，表明標(biāo)準(zhǔn)電磁模型和低能電磁模型中植入的替代強(qiáng)子輸運(yùn)模型較簡單，未考慮重子生成。雖然三者的γ和e-產(chǎn)額有少許差異，但較相近。模型4和5可產(chǎn)生重子n、H+、2H、3H、α和輕子e+，且產(chǎn)額較接近，但不能產(chǎn)生一些重離子(如Li+、Be2+、B3+、C4+、N3-、O2-等及其同位素)。模型6～9可產(chǎn)生重離子，說明級聯(lián)模型BIC和BERT加入了重離子產(chǎn)生模擬。其中模型6的n、H+產(chǎn)額較模型7高20%～30%，但2H、3H、α產(chǎn)額較模型7降低了30%～50%。模型6和模型8的一致，這是因?yàn)樵谀M的能量范圍內(nèi)兩者均使用了BERT和PRECO模型。模型9的H+產(chǎn)額較模型7的增加約38%，但重離子產(chǎn)額增加約79%，反映EMY模型的影響。從模擬時(shí)間看，采用簡單參數(shù)化處理的LHEP模型的時(shí)間最短，QGSP_BIC_EMY時(shí)間最長(LHEP的170%)，其他模型時(shí)間相仿(約LHEP的120%)。

圖2 模擬與文獻(xiàn)的PDD對比Fig.2 Comparison of simulated and published PDDs

表1 高斯能譜次級粒子產(chǎn)額比較Table 1 Comparison of secondary yield for Gaussian spectrum

2.2 質(zhì)子入射能量影響

表2列出了60 MeV和200 MeV的次級粒子產(chǎn)額，其總體特點(diǎn)與高斯能譜類似，但具體產(chǎn)額比有些差異，如n、H+產(chǎn)額方面，對于60 MeV，模型6較模型7高30%～50%，對于200 MeV僅高約20%；2H和3H產(chǎn)額方面，對于60 MeV，模型6較模型7降低了80%和5%，對于200 MeV降低了30%和27%，而α產(chǎn)額降低差不多(約45%),因此BERT模型(模型6)較BIC模型(模型7)對于低能傾向于產(chǎn)生n、H+，而對于高能更傾向于產(chǎn)生3H、α(但仍較BIC模型少)。對于60 MeV和200 MeV，模型9的H+產(chǎn)額較模型7增加了約49%和35%；而重離子產(chǎn)額增加了96%和69%，因此對比高斯能譜可見，EMY模型的影響隨入射質(zhì)子能量的減小而增大，尤其是對于治療眼部腫瘤的60 MeV的影響較顯著。

表2 60 MeV和200 MeV次級粒子產(chǎn)額比較Table 2 Comparison of secondary yields for 60 MeV and 200 MeV

2.3 材料影響

表3列出了模型4、6、7 3個(gè)具有代表性模型在不同材料模體中的次級粒子產(chǎn)額。由于肺密度較小，射程較長(Bragg深度約60 cm)，為對比需要，表3中給出了肺模體厚30 cm和90 cm兩種情況，其他材料模體均厚30 cm?？梢?，除γ和e-外，其他與水的類似，特別是重子產(chǎn)額表現(xiàn)出參數(shù)化LHEP模型與所研究材料的依賴性不大。從模擬時(shí)間看，在模體厚度一定的情況下，材料的密度越大，模擬時(shí)間越短。部分原因是：密度越大，粒子射程越短(如肺-0.1 cm/30 cm(網(wǎng)格厚0.1 cm，總厚30 cm)粒子輸運(yùn)達(dá)30 cm，但水、軟組織等僅約17 cm)，輸運(yùn)時(shí)間也就降低了。此外，肺材料中易發(fā)生散射的e-產(chǎn)額特別高(是其他材料的116～260倍)。若增大肺模體的厚度和網(wǎng)格厚度(0.3 cm)，追蹤粒子模擬直至其“死亡”或離開模體，其模擬時(shí)間不及精細(xì)網(wǎng)格(厚度為0.1 cm) 30 cm厚的模體的長，這說明Geant4中粒子模擬在邊界跨越上花費(fèi)的時(shí)間更長，較模體厚度對粒子輸運(yùn)的影響大。QGSP_BIC模型和QGSP_BERT模型與材料的關(guān)聯(lián)性增加，其中骨中各重子產(chǎn)額相對較高。整體上看，在模體厚度一定的情況下，材料的密度越大，模擬時(shí)間越短。

2.4 網(wǎng)格厚度與Bragg峰-入口劑量比

第二，資金匱乏。我國的養(yǎng)老服務(wù)項(xiàng)目在目前主要依靠政府購買的方式來維持，由于政府購買服務(wù)的資金比重較為偏低，導(dǎo)致服務(wù)受眾群體范圍變窄，很多需要幫助的老人不能享受到相應(yīng)的服務(wù)，從而使養(yǎng)老的福利性嚴(yán)重不足[9]。另外，由于智慧養(yǎng)老直接服務(wù)群體和間接服務(wù)群體的支付能力有限等因素的制約，再加上智能設(shè)備的購買和維護(hù)成本高等原因，限制了養(yǎng)老服務(wù)對象的積極參與；專業(yè)服務(wù)人才隊(duì)伍建設(shè)成本高；政府資金投入力度不夠，資金獲取渠道單一；系統(tǒng)編程制作與系統(tǒng)維護(hù)成本高，支撐系統(tǒng)軟件的硬件設(shè)施成本高等，都是當(dāng)前智慧養(yǎng)老服務(wù)模式建設(shè)所面臨的挑戰(zhàn)。

為研究不同能量照射下，不同網(wǎng)格厚度對Bragg峰-入口劑量比的影響。計(jì)算了入口和Bragg峰前后1 cm處分別劃分為厚1 mm、0.1 mm或0.01 mm網(wǎng)格的深度劑量曲線和Bragg峰-入口劑量比。圖3示出了60 MeV和高斯能譜對模型4的比較(其他模型與之相似)。可見，對于高斯能譜，精細(xì)的網(wǎng)格(0.1 mm或0.01 mm)和1 mm厚網(wǎng)格對劑量曲線無影響。但對于60 MeV，0.1 mm與0.01 mm的較一致(差異<1%)，但1 mm的與它們的有偏移。>160 MeV的情況的與高斯能譜的相同，本文未給出。本工作也研究了其他能量1 mm和0.1 mm厚網(wǎng)格差異，發(fā)現(xiàn)直到160 MeV，兩者的差異才小于2%。因此計(jì)算Bragg峰-入口劑量比，對于低能質(zhì)子(<160 MeV)，劑量模擬有必要?jiǎng)澐志?xì)的網(wǎng)格(約0.1 mm)；而對于高能質(zhì)子(≥160 MeV)，1 mm的網(wǎng)格可接受。由于一般蒙特卡羅程序算法處理邊界跨越均需耗費(fèi)大量時(shí)間用以判斷，所以過于精細(xì)的網(wǎng)格劃分需考慮時(shí)間成本。

表3 高斯能譜模型4、6、7不同材料模體中的次級粒子產(chǎn)額Table 3 Secondary yields of model 4, 6 and 7 with Gaussian in different material phantoms

圖3 網(wǎng)格厚度對模型4 PDD的影響Fig.3 Influence of grid thickness on PDD of model 4

圖4為Bragg峰-入口劑量比隨入射質(zhì)子能量和模體材料的變化。由于模型4～8均較相近，僅給出了模型5作代表。從圖可見，隨入射質(zhì)子能量的增加，Bragg峰-入口劑量比單調(diào)降低，說明隨著質(zhì)子射程的增加，更多的能量損失在傳播路徑上，使得Bragg峰處的治療優(yōu)勢大為降低。其中，模體材料的變化對Bragg峰-入口劑量比的影響很小，遠(yuǎn)小于模型9的EMY處理的影響(5%～13%)。圖4中也給出了Wroe等[6]模擬的結(jié)果，Wroe等使用了適用范圍0～150 MeV的PRECO模型。其中60 MeV時(shí)，Wroe等和模型5的結(jié)果較接近，但200 MeV時(shí)，Wroe等模擬的結(jié)果較本工作采用模型的高，這可能是Wroe等超限使用了PRECO模型。

圖4 入射質(zhì)子能量和模體材料對Bragg峰-入口劑量比的影響Fig.4 Influence of H+ energy and phantom material on dose ratio of Bragg peak to entrance

圖5為Bragg峰深度隨入射質(zhì)子能量和模體材料的變化。由于各模型的Bragg峰位置相同，圖中只代表性地給出了Bragg峰對應(yīng)的深度。其中水、軟組織和肌肉的Bragg峰深度變化規(guī)律相近，骨的Bragg峰深度變化相對緩慢些。此Bragg峰深度可用1個(gè)二次多項(xiàng)式擬合(R2>0.99)，如對軟組織，Bragg峰深度=-17.475 88+0.536 97E+0.004 15E2；對骨，Bragg峰深度=-11.681 63+0.369 25E+0.002 98E2。其中，能量E的單位為MeV，深度的單位為mm。圖5中也給出了Wroe等[6]的Bragg峰深度，從圖可見，與本模擬的結(jié)果較接近。

2.5 絕對劑量

表4給出了Bragg峰處60 MeV和250 MeV單個(gè)質(zhì)子沉積的絕對劑量。模型6和8在研究的能量范圍內(nèi)均使用了PRECO+BERT模型，數(shù)據(jù)相同，在表中未進(jìn)行區(qū)分?？梢娰|(zhì)子能量越小，Bragg峰處沉積的劑量越大。60 MeV的Bragg峰絕對劑量約為250 MeV的5倍，而60 MeV的Bragg峰-入口劑量比為250 MeV的2倍，因此60 MeV的入口劑量僅為250 MeV的2.5倍。靶區(qū)獲取的絕對劑量決定著腫瘤療效，在靶區(qū)劑量閾值一定的情況下，60 MeV質(zhì)子束前端對正常組織的輻射損傷相對較好。

圖5 入射質(zhì)子能量和模體材料對Bragg峰深度的影響Fig.5 Influence of H+ energy and phantom material on Bragg depth

表4 Bragg峰處單個(gè)入射質(zhì)子沉積劑量Table 4 Deposit dose for single H+ at Bragg peak

但從射程角度看，60 MeV的Bragg峰深度在水及等效介質(zhì)中僅3 cm，而250 MeV的高達(dá)37.5 cm，因此在治療深度腫瘤方面，高能量的質(zhì)子束具有優(yōu)勢。此外，通過旋轉(zhuǎn)分割照射，也可減少高能質(zhì)子對相同部位正常組織的輻射損傷。

模型4～8的結(jié)果較接近，模型9的絕對劑量明顯偏低，且此差異隨質(zhì)子能量的增加而增大(如60 MeV降低4%，250 MeV降低16%)。模型9使用QGSP_BIC模型外，添加了G4EmExtraPhysics( )、G4HadronElasticPhysics( )、G4QStoppingPhysics( )、G4IonBinaryCascade-Physics( )、G4NeutronTrackingCut( )處理，致使其模擬結(jié)果偏離其他物理模型及文獻(xiàn)發(fā)表數(shù)據(jù)。盡管Geant4hadrontherapy例物理清單中指出此為簡潔版物理模型，但I(xiàn)NFN使用的復(fù)雜物理模型不易獲取。此例題與Geant4新版本反復(fù)發(fā)布的簡潔物理模型，易使模擬者錯(cuò)誤采納，應(yīng)引起模擬者注意。

3 結(jié)論

本工作選取9個(gè)Geant4常用于治療質(zhì)子束能量范圍(60～250 MeV)的物理模型進(jìn)行比較研究。結(jié)果發(fā)現(xiàn)：若僅使用電磁物理模型，雖整體劑量分布上可接受，但微觀重子產(chǎn)額有缺失；參數(shù)化驅(qū)動模型LHEP的模擬時(shí)間最短，但不能產(chǎn)生復(fù)雜重子；級聯(lián)模型BERT和BIC的生物損傷輻射防護(hù)的預(yù)警性更強(qiáng)，尤以BIC模型更甚，因此QGSP_BERT模型、QGSP_BIC模型和FTFP_BERT模型較適合放療質(zhì)子束模擬。從平衡網(wǎng)格劃分與模擬效率的角度看，<160 MeV的質(zhì)子需0.1 mm厚網(wǎng)格，≥160 MeV的質(zhì)子1 mm厚網(wǎng)格即可。低能量(如60 MeV)質(zhì)子束的Bragg峰-入口劑量較高，前端對正常組織的輻射損傷相對較好。然而，由于能量與射程的二次多項(xiàng)式增長關(guān)系，在治療深度腫瘤方面，高能量的質(zhì)子束具有優(yōu)勢。通過旋轉(zhuǎn)照射也可減少高能量的質(zhì)子束對正常組織的輻射損傷。

Geant4中hadrontherapy例提供的物理模型QGSP_BIC_EMY，與Geant4新版本多次發(fā)布，易被錯(cuò)誤采納。本工作驗(yàn)證了hadrontherapy例物理模型與其他物理模型的差別，以提醒使用者。

[1] 唐勁天. 腫瘤質(zhì)子放射治療學(xué)[M]. 北京：中國醫(yī)藥科技出版社，2004.

[2] 劉云鵬，湯曉斌，謝芹，等. 質(zhì)子放療過程中射束與人眼體作用的蒙特卡羅模擬計(jì)算[J]. 原子能科學(xué)技術(shù)，2010，44(增刊)：602-607.

LIU Yunpeng, TANG Xiaobin, XIE Qin, et al. Monte-Carlo simulation of proton radiotherapy for human eye[J]. Atomic Energy Science and Technology, 2010, 44(Suppl.): 602-607(in Chinese).

[3] 許淑艷. 蒙特卡羅方法在實(shí)驗(yàn)核物理中的應(yīng)用[M]. 2 版. 北京：原子能出版社，2006.

[4] 鞠志萍，曹午飛，劉小偉. 質(zhì)子散射角分布的蒙特卡羅模擬[J]. 物理學(xué)報(bào)，2009,58(1):174-177.

JU Zhiping, CAO Wufei, LIU Xiaowei. Study of scattering angular distribution of proton using Monte-Carlo method[J]. Acta Physics Sinica, 2009, 58(1): 174-177(in Chinese).

[5] 謝朝，鄒煉，侯氫，等. 質(zhì)子束治療中非均勻組織的等效水厚度修正研究[J]. 物理學(xué)報(bào)，2013,62(6):068701-5.

XIE Chao, ZOU Lian, HOU Qing, et al. Study of inhomogeneous tissue equivalent water thickness correction method in proton therapy[J]. Acta Physics Sinica, 2013, 62(6): 068701-5(in Chinese).

[6] WROE A J, CORNELIUS I M, ROSENFELD A B, et al. The role of nonelastic reactions in absorbed dose distributions from therapeutic proton beams in different medium[J]. Med Phys, 2005, 32(1): 37-41.

[7] Geant4 is a toolkit for the simulation of the passage of particles through matter[R/OL]. http:∥geant4.cern.ch/.

[8] CIRRONE G A P, CUTTONE G, ROSA F D, et al. Validation of Geant4 physics models for the simulation of the proton Bragg peak[C]∥Conf Rec 2006 IEEE Nucl Sci Symp. [S. l.]: [s. n.], 2006: 788-792.

[9] TOSHITO T, BAGULYA A L, IVANCHENKO A V, et al. Validation of new Geant4 electromagnetic physics models for ion therapy applications[J]. Progress in Nuclear Science and Technology, 2011, 2: 918-922.

[10]PAGANETTI H. Nuclear interactions in proton therapy: Dose and relative biological effect distributions originating from primary and secondary particles[J]. Phys Med Biol, 2002, 47: 747-764.

Effect of Different Geant4 Physical Models on Simulation of Radiotherapy Proton Beam

LIN Hui1, XIE Cong1, ZHANG Yong-jun1,2, XIONG Zhen-yu1,WU Dong-sheng1,*, CAO Rui-fen3, FDS Team3

(1.SchoolofElectronicScience&ApplicationPhysics,HefeiUniversityofTechnology,Hefei230009,China; 2.DepartmentofBasicCourses,BengbuMedicalCollege,Bengbu233000,China;3.InstituteofNuclearEnergySafetyTechnology,ChineseAcademyofSciences,Hefei230031,China)

The Monte Carlo method is usually used to simulate the feature of the external radiotherapy of proton beam to optimize the clinical therapy scenario. The selection of a suitable physical model is very critical to ensure the correction and the efficiency of the simulation. This work studied nine Geant4 physical models, which are often used in the simulation of radiotherapy of proton beam with energy of 60-250 MeV. The dosage features and the microcosmic secondary yields in different materials were compared. The result shows that, although the standard EM, the low energy Penelope model and the low energy Livermore model can output the right dose distribution, they can’t simulate the baryon yields. The parametrization driven model LHEP and the Quark gluon string-pre-compound model both can’t simulate the complex ion transportation. Thus the above five models are not enough for the research calculation in the radiation protection and the radiation damage. The QGSP_BIC_EMY model in Geant4 hadron therapy example is very bad for underestimating the dose ratio of Bragg peak to entrance and absolute dose. The QGSP_BERT model, QGSP_BIC model and FTFP_BERT model are suitable for the radiotherapy proton beam with energy of 60-250 MeV.

physical model; proton radiotherapy; microscopic yield; dosage characteristic

2014-03-10;

2014-06-25

中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)資助項(xiàng)目(XDA03040000)；合肥工業(yè)大學(xué)科學(xué)發(fā)展研究基金資助項(xiàng)目(2013HGXJ0193，J2014HGXJ0093)；合肥工業(yè)大學(xué)自主創(chuàng)新項(xiàng)目資助(2012HGZY0007)

林 輝(1973—)，女，安徽阜陽人，副教授，博士，核技術(shù)醫(yī)學(xué)應(yīng)用專業(yè)

*通信作者：吳東升，E-mail: wudongsheng@hfut.edu.cn

TL84

A

1000-6931(2015)07-1290-08

10.7538/yzk.2015.49.07.1290