• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2015-05-04 02:51:28張擁軍熊楨宇吳東升曹瑞芬FDS團隊
    原子能科學技術 2015年7期
    關鍵詞:產(chǎn)額模體強子

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

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

    ?

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

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

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

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

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

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

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

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

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

    1 材料和方法

    1.1 程序發(fā)展

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

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

    1.2 Geant4物理模型

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

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

    Geant4提供以下強子相互作用模型。

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

    2) 級聯(lián)模型(CM,cascade model)。模擬了基本粒子和次級粒子的級聯(lián)輸運過程,包括BIC(binary cascade model(約<10 GeV))和BERT(Bertini cascade model(約<10 GeV))。但級聯(lián)模型不適合低能,因為低能涉及蒸發(fā)和預平衡衰變。僅當粒子的動能高于閾值(75 MeV)時才能使用級聯(lián)模型。

    3) 預復合核模型(PRECO,pre-compound model)。預復合核模型是強子動力學模型的擴展,它將強子動力學模型適用范圍擴展到核-核非彈性碰撞??紤]核反應的預復合核階段直到非平衡態(tài)。從激發(fā)核里釋放核碎片或光子的模擬采用平衡模型。PRECO適用于0~150 MeV范圍內(nèi)的質子、中子模擬。

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

    本工作考察了以下物理模型對放療質子劑量和微觀產(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為電磁模型。為考察強子模型的影響,模型4~9在含有標準電磁模型選項3(即模型1)之外又增添了強子模型處理。模型9為Geant4advancedhadrontherapy例提供的物理模型,在使用QGSP_BIC模型之外添加了G4EmExtraPhysics( )、G4Hadron-ElasticPhysics( )、G4QStoppingPhysics( )、G4IonBinaryCascadePhysics( )、G4Neutron-TrackingCut()處理。除模型1外,其他模型均采用了放射活度衰變離子物理處理。

    1.3 模擬幾何

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

    1.4 模型校驗

    2 結果與討論

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

    為驗證模擬的可靠性,模擬的PDD與文獻[10]數(shù)據(jù)進行了對比(圖2)。從圖2可見,除物理模型9和文獻相差較大(約3.5%)外,其他模型幾乎完全一致(入口處(<5 cm)約1.0%,其他處約0.2%)。表1列出了各模型的次級粒子產(chǎn)額,其中模型1~3僅產(chǎn)生γ和e-,無重子產(chǎ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的一致,這是因為在模擬的能量范圍內(nèi)兩者均使用了BERT和PRECO模型。模型9的H+產(chǎn)額較模型7的增加約38%,但重離子產(chǎn)額增加約79%,反映EMY模型的影響。從模擬時間看,采用簡單參數(shù)化處理的LHEP模型的時間最短,QGSP_BIC_EMY時間最長(LHEP的170%),其他模型時間相仿(約LHEP的120%)。

    圖2 模擬與文獻的PDD對比Fig.2 Comparison of simulated and published PDDs

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

    2.2 質子入射能量影響

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

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

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

    為研究不同能量照射下,不同網(wǎng)格厚度對Bragg峰-入口劑量比的影響。計算了入口和Bragg峰前后1 cm處分別劃分為厚1 mm、0.1 mm或0.01 mm網(wǎng)格的深度劑量曲線和Bragg峰-入口劑量比。圖3示出了60 MeV和高斯能譜對模型4的比較(其他模型與之相似)??梢?,對于高斯能譜,精細的網(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%。因此計算Bragg峰-入口劑量比,對于低能質子(<160 MeV),劑量模擬有必要劃分精細的網(wǎng)格(約0.1 mm);而對于高能質子(≥160 MeV),1 mm的網(wǎng)格可接受。由于一般蒙特卡羅程序算法處理邊界跨越均需耗費大量時間用以判斷,所以過于精細的網(wǎng)格劃分需考慮時間成本。

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

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

    圖5為Bragg峰深度隨入射質子能量和模體材料的變化。由于各模型的Bragg峰位置相同,圖中只代表性地給出了Bragg峰對應的深度。其中水、軟組織和肌肉的Bragg峰深度變化規(guī)律相近,骨的Bragg峰深度變化相對緩慢些。此Bragg峰深度可用1個二次多項式擬合(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峰深度,從圖可見,與本模擬的結果較接近。

    2.5 絕對劑量

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

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

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

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

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

    3 結論

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

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

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

    [2] 劉云鵬,湯曉斌,謝芹,等. 質子放療過程中射束與人眼體作用的蒙特卡羅模擬計算[J]. 原子能科學技術,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] 許淑艷. 蒙特卡羅方法在實驗核物理中的應用[M]. 2 版. 北京:原子能出版社,2006.

    [4] 鞠志萍,曹午飛,劉小偉. 質子散射角分布的蒙特卡羅模擬[J]. 物理學報,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] 謝朝,鄒煉,侯氫,等. 質子束治療中非均勻組織的等效水厚度修正研究[J]. 物理學報,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

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

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

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

    TL84

    A

    1000-6931(2015)07-1290-08

    10.7538/yzk.2015.49.07.1290

    猜你喜歡
    產(chǎn)額模體強子
    一個可靠和準確的光電產(chǎn)額譜模型及應用*
    物理學報(2021年10期)2021-06-01 08:31:40
    基于Matrix Profile的時間序列變長模體挖掘
    基于貝葉斯更新方法的235U熱中子獨立裂變產(chǎn)額協(xié)方差估計
    針對裂變產(chǎn)額和半衰期的燃耗計算靈敏度和不確定度分析方法
    罰你沒商量
    故事會(2020年12期)2020-06-19 08:52:36
    植入(l, d)模體發(fā)現(xiàn)若干算法的實現(xiàn)與比較
    我坐在一塊石頭上(外五首)
    都市(2018年12期)2018-09-10 10:13:17
    基于網(wǎng)絡模體特征攻擊的網(wǎng)絡抗毀性研究
    How to Care about the students’emotions and attitudes and promote their full development
    基于模體演化的時序鏈路預測方法
    自動化學報(2016年5期)2016-04-16 03:38:40
    99久久成人亚洲精品观看| 久久精品影院6| av在线亚洲专区| 久久综合国产亚洲精品| 欧美三级亚洲精品| 亚洲不卡免费看| 欧美性猛交╳xxx乱大交人| 亚洲精品乱码久久久久久按摩| 国产精品永久免费网站| 一级二级三级毛片免费看| 亚洲乱码一区二区免费版| 又黄又爽又刺激的免费视频.| 99国产精品一区二区蜜桃av| 亚洲一区二区三区色噜噜| 免费无遮挡裸体视频| 日本撒尿小便嘘嘘汇集6| 我的老师免费观看完整版| 女的被弄到高潮叫床怎么办| 一区福利在线观看| 国产久久久一区二区三区| 91在线精品国自产拍蜜月| 2021天堂中文幕一二区在线观| 欧美高清性xxxxhd video| 国产精品人妻久久久影院| 婷婷亚洲欧美| 婷婷色av中文字幕| 亚洲一级一片aⅴ在线观看| 国产精品一区二区在线观看99 | 日韩欧美 国产精品| 一级黄片播放器| 91麻豆精品激情在线观看国产| 不卡视频在线观看欧美| 我要看日韩黄色一级片| 不卡一级毛片| 中文在线观看免费www的网站| 深爱激情五月婷婷| 又粗又爽又猛毛片免费看| 中文在线观看免费www的网站| 日本爱情动作片www.在线观看| 亚洲真实伦在线观看| 亚洲在线观看片| 亚洲成人久久爱视频| 久久综合国产亚洲精品| 国产高清激情床上av| 国产免费男女视频| 夜夜夜夜夜久久久久| 搡老妇女老女人老熟妇| 国内揄拍国产精品人妻在线| 亚洲国产高清在线一区二区三| 99精品在免费线老司机午夜| 国产69精品久久久久777片| 免费看美女性在线毛片视频| 国产一区二区三区av在线 | 五月玫瑰六月丁香| 99国产极品粉嫩在线观看| www.av在线官网国产| 最近最新中文字幕大全电影3| 99久久精品国产国产毛片| 久久久久久大精品| 国产一级毛片七仙女欲春2| 又黄又爽又刺激的免费视频.| 99国产精品一区二区蜜桃av| 日韩制服骚丝袜av| 人体艺术视频欧美日本| 成人美女网站在线观看视频| 亚洲欧美清纯卡通| 亚洲熟妇中文字幕五十中出| 午夜福利在线在线| 国内精品一区二区在线观看| 久久人妻av系列| 久久午夜亚洲精品久久| 国产精品精品国产色婷婷| 日本色播在线视频| 中文亚洲av片在线观看爽| 国产日韩欧美在线精品| 亚洲成人久久爱视频| 成人特级av手机在线观看| 色噜噜av男人的天堂激情| 白带黄色成豆腐渣| 国产成人a区在线观看| 中文字幕精品亚洲无线码一区| 听说在线观看完整版免费高清| 3wmmmm亚洲av在线观看| 亚洲国产高清在线一区二区三| 亚洲国产精品国产精品| 好男人在线观看高清免费视频| av免费观看日本| 亚洲中文字幕一区二区三区有码在线看| 波多野结衣巨乳人妻| 美女高潮的动态| 春色校园在线视频观看| 日本爱情动作片www.在线观看| 偷拍熟女少妇极品色| 男女视频在线观看网站免费| 亚洲美女视频黄频| 亚洲成人久久性| 内射极品少妇av片p| 国产精品麻豆人妻色哟哟久久 | 99久久人妻综合| 1024手机看黄色片| 国产黄片视频在线免费观看| 国产一区二区亚洲精品在线观看| 久久亚洲精品不卡| 国产乱人视频| 婷婷亚洲欧美| 波多野结衣巨乳人妻| 一本精品99久久精品77| 1024手机看黄色片| 嫩草影院新地址| 波多野结衣高清作品| 久久久久久久久久久免费av| 99久久精品一区二区三区| 尤物成人国产欧美一区二区三区| 麻豆国产av国片精品| 亚洲一级一片aⅴ在线观看| 男人舔女人下体高潮全视频| 1000部很黄的大片| 色综合色国产| 黄色欧美视频在线观看| 欧美极品一区二区三区四区| 狂野欧美白嫩少妇大欣赏| 亚洲aⅴ乱码一区二区在线播放| 亚洲熟妇中文字幕五十中出| 哪个播放器可以免费观看大片| 一区二区三区高清视频在线| 伦精品一区二区三区| 国产精品蜜桃在线观看 | 成人国产麻豆网| 九草在线视频观看| 亚洲欧美日韩高清在线视频| 日本撒尿小便嘘嘘汇集6| av卡一久久| 国国产精品蜜臀av免费| 中文欧美无线码| 日本黄色视频三级网站网址| 国产精品一区二区在线观看99 | 亚洲中文字幕日韩| 欧美人与善性xxx| 国产精品电影一区二区三区| 99热这里只有是精品50| 99久久人妻综合| 精品99又大又爽又粗少妇毛片| 日韩一区二区三区影片| 毛片女人毛片| 亚洲精品国产av成人精品| 亚洲aⅴ乱码一区二区在线播放| 国国产精品蜜臀av免费| 变态另类丝袜制服| 97超视频在线观看视频| 久久99热这里只有精品18| 成人国产麻豆网| 久久这里只有精品中国| 久久精品国产自在天天线| 亚洲精品日韩av片在线观看| 国产三级中文精品| 国产三级中文精品| 男女啪啪激烈高潮av片| 身体一侧抽搐| 日本黄大片高清| 国产一级毛片七仙女欲春2| 国产精品久久久久久久电影| 亚洲最大成人中文| 国产伦一二天堂av在线观看| 国产伦理片在线播放av一区 | 国产伦精品一区二区三区四那| 联通29元200g的流量卡| 欧美成人免费av一区二区三区| 亚洲中文字幕日韩| av在线蜜桃| 国产伦理片在线播放av一区 | 国产成年人精品一区二区| 深夜a级毛片| 99久久久亚洲精品蜜臀av| av在线蜜桃| 身体一侧抽搐| 久久婷婷人人爽人人干人人爱| 国产成年人精品一区二区| 美女内射精品一级片tv| 黄色视频,在线免费观看| 丰满乱子伦码专区| 亚洲无线在线观看| 99在线人妻在线中文字幕| 久久精品影院6| 国产一区二区激情短视频| 高清日韩中文字幕在线| 人人妻人人澡欧美一区二区| 亚洲av中文av极速乱| 身体一侧抽搐| 精品人妻偷拍中文字幕| 色综合站精品国产| 深夜精品福利| 爱豆传媒免费全集在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成年av动漫网址| 非洲黑人性xxxx精品又粗又长| 午夜老司机福利剧场| 成人午夜高清在线视频| 日韩成人av中文字幕在线观看| 中国国产av一级| 我要搜黄色片| 观看美女的网站| 激情 狠狠 欧美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩亚洲欧美综合| 一区二区三区四区激情视频 | 亚洲国产精品久久男人天堂| 成人二区视频| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕精品亚洲无线码一区| 久久婷婷人人爽人人干人人爱| av天堂在线播放| 91午夜精品亚洲一区二区三区| 亚洲av二区三区四区| 国语自产精品视频在线第100页| 人人妻人人澡人人爽人人夜夜 | 国产精品一区二区性色av| 亚洲一区二区三区色噜噜| 亚洲精品日韩在线中文字幕 | 国产精品久久视频播放| 免费搜索国产男女视频| 国产精品嫩草影院av在线观看| 熟女人妻精品中文字幕| 99热全是精品| 免费看日本二区| 高清毛片免费看| 国产高清三级在线| 国产激情偷乱视频一区二区| 中国美白少妇内射xxxbb| 国产精品一区二区三区四区免费观看| 人体艺术视频欧美日本| 成人毛片a级毛片在线播放| 男人舔奶头视频| 美女大奶头视频| 免费av毛片视频| 黑人高潮一二区| 九九热线精品视视频播放| 亚洲av第一区精品v没综合| 国产精品久久久久久av不卡| 日本一二三区视频观看| 美女 人体艺术 gogo| 九九爱精品视频在线观看| 午夜精品一区二区三区免费看| 国内精品久久久久精免费| 国产成人a区在线观看| 国产免费男女视频| 国产高清不卡午夜福利| 国产成人精品久久久久久| 一进一出抽搐gif免费好疼| 午夜福利在线在线| 成人三级黄色视频| 我的女老师完整版在线观看| 又粗又硬又长又爽又黄的视频 | 麻豆国产97在线/欧美| 91久久精品国产一区二区三区| 国产精品三级大全| 一区二区三区四区激情视频 | 精品久久久久久久久久久久久| 中文字幕久久专区| 人人妻人人看人人澡| 伊人久久精品亚洲午夜| 99视频精品全部免费 在线| 1000部很黄的大片| 3wmmmm亚洲av在线观看| 国产精品一区二区三区四区久久| 少妇裸体淫交视频免费看高清| 久久九九热精品免费| 亚洲国产精品国产精品| 中文资源天堂在线| 亚洲精品乱码久久久久久按摩| 人人妻人人澡欧美一区二区| 热99在线观看视频| 国产成年人精品一区二区| 青青草视频在线视频观看| 男女视频在线观看网站免费| 99在线视频只有这里精品首页| 免费在线观看成人毛片| 国产综合懂色| 99久久精品一区二区三区| 看十八女毛片水多多多| 欧美区成人在线视频| 亚洲精品自拍成人| 日韩制服骚丝袜av| 噜噜噜噜噜久久久久久91| 久久韩国三级中文字幕| 国产精品嫩草影院av在线观看| 91精品国产九色| 成人午夜精彩视频在线观看| 免费无遮挡裸体视频| 中文字幕av成人在线电影| 精品久久久久久久久亚洲| 色哟哟·www| 精品久久久久久久末码| 最近2019中文字幕mv第一页| 精品久久久久久久久久久久久| 久久亚洲国产成人精品v| 精品熟女少妇av免费看| 欧美日韩在线观看h| 亚洲一级一片aⅴ在线观看| 亚洲av第一区精品v没综合| 最近2019中文字幕mv第一页| 国内久久婷婷六月综合欲色啪| 久久久久免费精品人妻一区二区| 国产精品一区二区性色av| 少妇被粗大猛烈的视频| 亚洲最大成人av| 免费观看的影片在线观看| 校园春色视频在线观看| 亚洲性久久影院| 国产熟女欧美一区二区| 国产精品福利在线免费观看| 黑人高潮一二区| 国产高清有码在线观看视频| 亚洲av免费高清在线观看| 久久精品国产99精品国产亚洲性色| 99热全是精品| 国产成人aa在线观看| 特大巨黑吊av在线直播| 精品久久久久久久人妻蜜臀av| 成人永久免费在线观看视频| 1000部很黄的大片| 午夜a级毛片| 丰满乱子伦码专区| 久久久久九九精品影院| 成年免费大片在线观看| 国产乱人偷精品视频| 久久久欧美国产精品| 成人亚洲欧美一区二区av| 日韩欧美精品v在线| 九色成人免费人妻av| 夜夜爽天天搞| 国产人妻一区二区三区在| 天堂影院成人在线观看| 国产精品国产高清国产av| 99久久久亚洲精品蜜臀av| 久久婷婷人人爽人人干人人爱| 欧美性猛交╳xxx乱大交人| 毛片女人毛片| 日本在线视频免费播放| 一个人看视频在线观看www免费| av在线观看视频网站免费| 精品久久国产蜜桃| 毛片女人毛片| 级片在线观看| 国产精品一区二区在线观看99 | 日本与韩国留学比较| 亚洲av.av天堂| 最近中文字幕高清免费大全6| 少妇丰满av| 国语自产精品视频在线第100页| 日本免费一区二区三区高清不卡| 国产精品,欧美在线| 性插视频无遮挡在线免费观看| 欧美3d第一页| 欧美变态另类bdsm刘玥| 免费观看在线日韩| 久久精品夜夜夜夜夜久久蜜豆| av在线播放精品| 国产精品蜜桃在线观看 | 欧美人与善性xxx| 欧美+亚洲+日韩+国产| 久久精品综合一区二区三区| 人妻少妇偷人精品九色| 黄色欧美视频在线观看| 亚洲电影在线观看av| 精品一区二区免费观看| 亚洲国产精品成人综合色| 国产精品蜜桃在线观看 | 亚洲av中文av极速乱| 国产黄a三级三级三级人| 99在线视频只有这里精品首页| 国产av麻豆久久久久久久| 中文字幕av成人在线电影| 少妇的逼好多水| 最后的刺客免费高清国语| 亚洲欧美日韩高清专用| 天堂影院成人在线观看| 插逼视频在线观看| 2022亚洲国产成人精品| 久久久精品94久久精品| 国产视频内射| 国内精品一区二区在线观看| 99视频精品全部免费 在线| 在线观看免费视频日本深夜| 热99在线观看视频| 听说在线观看完整版免费高清| 亚洲欧美中文字幕日韩二区| 久久久久性生活片| 亚洲国产精品成人综合色| 性插视频无遮挡在线免费观看| 久久久色成人| 中出人妻视频一区二区| 久久久欧美国产精品| 亚洲色图av天堂| av国产免费在线观看| 中文精品一卡2卡3卡4更新| 久久久久网色| 九九爱精品视频在线观看| 久久这里有精品视频免费| 中出人妻视频一区二区| 好男人在线观看高清免费视频| 中出人妻视频一区二区| 国产亚洲av片在线观看秒播厂 | 免费看a级黄色片| 蜜桃久久精品国产亚洲av| 我要搜黄色片| 深夜精品福利| 综合色av麻豆| 哪里可以看免费的av片| 久久6这里有精品| 校园春色视频在线观看| 成人午夜精彩视频在线观看| 高清午夜精品一区二区三区 | 乱人视频在线观看| 免费看美女性在线毛片视频| 久久精品夜色国产| 精品久久久久久久久久久久久| 亚洲成人av在线免费| 国产av在哪里看| 国产精品一区二区三区四区久久| 免费观看的影片在线观看| 国产单亲对白刺激| 高清毛片免费看| 欧美激情国产日韩精品一区| 亚洲在久久综合| 少妇人妻一区二区三区视频| 久久人人爽人人爽人人片va| 狂野欧美白嫩少妇大欣赏| 午夜精品国产一区二区电影 | 国产在线男女| 少妇的逼水好多| 国内精品美女久久久久久| avwww免费| 一边摸一边抽搐一进一小说| 婷婷色av中文字幕| 婷婷亚洲欧美| 欧美日本亚洲视频在线播放| 亚洲无线在线观看| 少妇熟女欧美另类| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 男插女下体视频免费在线播放| 成人综合一区亚洲| 一级黄色大片毛片| 国产国拍精品亚洲av在线观看| 国产黄片视频在线免费观看| 给我免费播放毛片高清在线观看| 最近的中文字幕免费完整| 中文字幕av在线有码专区| 一级黄色大片毛片| 3wmmmm亚洲av在线观看| 国产精品国产三级国产av玫瑰| 日本与韩国留学比较| 黄片wwwwww| 夜夜爽天天搞| 亚洲精品乱码久久久v下载方式| 国产伦一二天堂av在线观看| 国产白丝娇喘喷水9色精品| 久久精品久久久久久噜噜老黄 | 国语自产精品视频在线第100页| 久久久午夜欧美精品| 午夜福利在线观看吧| 中文字幕av在线有码专区| 婷婷色综合大香蕉| 亚洲aⅴ乱码一区二区在线播放| 国产 一区 欧美 日韩| 亚洲激情五月婷婷啪啪| 搡老妇女老女人老熟妇| 精品日产1卡2卡| 久久热精品热| 免费观看人在逋| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久免费视频| 观看免费一级毛片| 日韩一本色道免费dvd| 国产精品爽爽va在线观看网站| 欧美成人精品欧美一级黄| 亚洲无线观看免费| 亚洲av中文av极速乱| 嫩草影院入口| 国产又黄又爽又无遮挡在线| 此物有八面人人有两片| 综合色av麻豆| 麻豆av噜噜一区二区三区| 美女被艹到高潮喷水动态| 中文字幕av成人在线电影| 丝袜美腿在线中文| 中文字幕免费在线视频6| 中文在线观看免费www的网站| 成人毛片a级毛片在线播放| 亚洲av免费在线观看| 免费看日本二区| 久久亚洲国产成人精品v| 成人亚洲欧美一区二区av| 欧美一区二区国产精品久久精品| 男插女下体视频免费在线播放| 国产av麻豆久久久久久久| 欧美三级亚洲精品| 亚洲熟妇中文字幕五十中出| 国产日本99.免费观看| 欧美高清成人免费视频www| 婷婷六月久久综合丁香| 天天躁日日操中文字幕| 特大巨黑吊av在线直播| 国产午夜福利久久久久久| 免费观看a级毛片全部| 在线观看美女被高潮喷水网站| 直男gayav资源| 成人特级黄色片久久久久久久| 中文字幕熟女人妻在线| 人妻系列 视频| 久久久欧美国产精品| 91av网一区二区| 男女视频在线观看网站免费| 99久久精品一区二区三区| 2022亚洲国产成人精品| 最近手机中文字幕大全| 国产成人a区在线观看| 中文亚洲av片在线观看爽| 免费大片18禁| 又爽又黄无遮挡网站| 亚洲国产高清在线一区二区三| 久久中文看片网| 免费搜索国产男女视频| 精品久久久久久久人妻蜜臀av| 神马国产精品三级电影在线观看| 午夜福利在线观看吧| 国产老妇伦熟女老妇高清| 欧美日韩国产亚洲二区| 亚洲国产欧美在线一区| 国产精品美女特级片免费视频播放器| 黄色一级大片看看| 亚洲一级一片aⅴ在线观看| 成人鲁丝片一二三区免费| 日韩国内少妇激情av| 热99在线观看视频| 国产精品美女特级片免费视频播放器| 日韩成人伦理影院| av在线观看视频网站免费| 搡老妇女老女人老熟妇| 在线播放国产精品三级| 最好的美女福利视频网| 成人午夜高清在线视频| 一级毛片aaaaaa免费看小| 免费观看a级毛片全部| 精品人妻熟女av久视频| 国产精品电影一区二区三区| 最好的美女福利视频网| 久久久久网色| 亚洲国产日韩欧美精品在线观看| 国产精品免费一区二区三区在线| 性色avwww在线观看| 国产精品蜜桃在线观看 | .国产精品久久| 蜜臀久久99精品久久宅男| 久久午夜亚洲精品久久| av在线亚洲专区| 波多野结衣高清作品| 国产黄色视频一区二区在线观看 | 亚洲国产欧洲综合997久久,| 亚洲图色成人| 久久精品91蜜桃| 日本av手机在线免费观看| 又爽又黄a免费视频| 在线观看66精品国产| 69人妻影院| 别揉我奶头 嗯啊视频| 少妇人妻一区二区三区视频| 在现免费观看毛片| 精品不卡国产一区二区三区| 久久久久久久午夜电影| 久久精品国产99精品国产亚洲性色| 精品人妻偷拍中文字幕| 亚洲国产精品成人综合色| 亚洲,欧美,日韩| 色综合亚洲欧美另类图片| 又爽又黄无遮挡网站| 熟女电影av网| 国产精品日韩av在线免费观看| 最近最新中文字幕大全电影3| 国内精品美女久久久久久| 高清毛片免费观看视频网站| 国产伦在线观看视频一区| 成年女人永久免费观看视频| 国产一级毛片在线| 蜜臀久久99精品久久宅男| 精品久久久久久久末码| 日韩视频在线欧美| 国产精品,欧美在线| 欧美丝袜亚洲另类| 国语自产精品视频在线第100页| 99热网站在线观看| 国产成人aa在线观看| 男的添女的下面高潮视频| 69人妻影院| 亚洲在久久综合| 五月玫瑰六月丁香| 老师上课跳d突然被开到最大视频| 2022亚洲国产成人精品| 少妇熟女欧美另类| av在线蜜桃| 99在线视频只有这里精品首页| 国产成人一区二区在线| 免费av毛片视频| 在线播放国产精品三级| 亚洲成人精品中文字幕电影| 成年免费大片在线观看| 午夜福利在线观看免费完整高清在 | 亚洲av二区三区四区| 丰满人妻一区二区三区视频av| 免费看av在线观看网站| av在线蜜桃| 久久人人爽人人片av| 一个人看的www免费观看视频| 有码 亚洲区|