RBS方法中的蒙特卡羅并行計算研究

1、    rbs方法中的蒙特卡羅并行計算研究    付濤+張麗然摘要：本文利用corteo在cluster集群運行全蒙特卡羅計算，其計算速度可以通過增加計算節(jié)點而得到成倍的提升。通過引入并行的corteo程序來模擬計算盧瑟福背散射能譜，可以在較短的時間內(nèi)模擬足夠多的粒子數(shù)來達到要求的精度，很好地解決了rbs方法中由于引入蒙特卡羅模擬而造成的耗時問題。關(guān)鍵詞：氚分析rbs方法；monte carlo； 模擬并行計算1. 引言精確地測量材料中氚的含量及深度分布對于核聚變研究和核技術(shù)應(yīng)用等領(lǐng)域具有重要的意義。目前氚分析的方法有很多，較為常見的如：盧瑟福背散射（r

2、bs）、彈性反沖探測（erda）、核反應(yīng)分析（nra）、增強質(zhì)子背散射分析（epbs）、粒子誘發(fā)x射線熒光分析（pixe）、低能粒子散射（leis）、溝道效應(yīng)（channelling）以及粒子誘發(fā)射線分析（pige）、掃描離子透射顯微鏡（stim）、離子致發(fā)光（il）等。這些技術(shù)的研究與應(yīng)用，極大地促進了交叉學(xué)科的發(fā)展，為新材料、新器件、新能源的開發(fā)做出了巨大的貢獻1-7。盧瑟福背散射（rbs）8作為一種簡便方便的氚分析方法，它通過探測背散射能譜來分析材料中的氚含量及其深度分布。為了精確測量材料表面的氚含量，在測量時將含氚材料保持在氬氣環(huán)境中。背散射能譜的分析通常使用解析軟件simnra9或者

3、利用基于蒙特卡羅方法的模擬軟件corteo10進行分析。2. 蒙特卡羅并行程序corteo是一種高度依賴由不同線程共享的中央數(shù)據(jù)（即散射數(shù)據(jù)和截止能）的多線程程序。所有線程在相同的內(nèi)存中讀取中央數(shù)據(jù)信息。更具體地說corteo支持多線程但不能提供等節(jié)點（node）之間的信息傳遞。所以在集群上運行corteo需要大量的工作。在集群上，需要單獨使用特定于集群的編譯器編譯corteo。我們這里所用的集群的編譯器是由美國intel公司開發(fā)icc編譯器（intel c+ compiler）。這種編譯器計算速度比gcc（gnu compiler collection）快了近510%。此次，我們采用cort

4、eo程序在四川大學(xué)原子核科學(xué)技術(shù)研究所提供的64核（8個節(jié)點，每個節(jié)點8核）計算機進行cluster集群計算。在cluster集群運行全蒙特卡羅corteo首先需要在遠(yuǎn)端服務(wù)器上傳corteo文件夾（包含corteo程序源代碼和運行所需組件），接著上傳輸入文件corteo.in（讀取自corteoui.exe，設(shè)置了發(fā)射源、樣品、探測器等一系列參數(shù)）到cluster所對應(yīng)的corteo子目錄。然后使用特定于集群的編譯器icc編譯corteo。為了估計在集群上計算時間，可以將模擬粒子數(shù)適當(dāng)設(shè)置小一點，例如我們將鎢樣品模擬粒子數(shù)設(shè)置為109計算花了1 h，因此當(dāng)我們將模擬粒子數(shù)設(shè)置為1010時，時

5、間就為10 h（同時也要觀察在測試計算中產(chǎn)生多少數(shù)據(jù)，以免集群在全速計算的時候不堪重負(fù)）。節(jié)點數(shù)設(shè)置為1，因為每個工作是在一個節(jié)點上運行。最后就要用到qsub命令向cluster提交任務(wù)，之前設(shè)置的corteo.pbs就是提交的內(nèi)容。我們向運行界面輸入：qsub corteo.pbs，任務(wù)提交完畢。用于分析的mo基底氚化鈦薄膜的制備方法如下：首先在清潔光滑1 mm厚mo基底表層蒸鍍一層5 ？m鈦膜，置于專用的充氚系統(tǒng)當(dāng)中吸氣到適當(dāng)?shù)碾傲?。鈦膜的吸氚量可以通過充氚系統(tǒng)壓強的變化，可以通過理想氣體狀態(tài)方程pvt計算得到。經(jīng)過pvt計算得到鈦膜樣品中t/ti原子比大約為1.51。在氚樣品中氚的深度分

6、布應(yīng)該均勻否則氚的濃度會隨著深度的增加而減少。在氦氣和氬氣混合氣體中使用磁控濺射方法60在光滑1 mm厚mo（或者si）基體上鍍一層大約1.5 ？m鈦膜。經(jīng)過epbs方法測得鈦薄膜樣品中he/ti原子比大約為0.6。在氦樣品中氦的深度分布應(yīng)該均勻。從圖1中我們可以看出，在pc機上做利用simnra解析計算得到的背射線能譜曲線光滑度很差，尤其是在能譜低能端。而用corteo并行程序在服務(wù)器上模擬計算得到的背散射能譜曲線很光滑。以上研究結(jié)果表明，通過引入corteo并行程序來模擬rbs方法中的背散射，可以在較短的時間內(nèi)模擬足夠多的粒子數(shù)來達到要求的精度，很好地解決了實際應(yīng)用中蒙特卡羅模擬的耗時問題

7、。3. 結(jié)論基于蒙特卡羅串行程序corteo的基礎(chǔ)上發(fā)展了蒙特卡羅并行程序。通過引入蒙特卡羅并行程序來模擬rbs方法中的背線譜，可以很好地解決rbs方法中引入蒙特卡羅模擬而造成的耗時問題，從而使該方法具有了更強的實用性。圖1 2 mev質(zhì)子垂直入射到mo基底表層ti（t）樣品在165°散射角得到的背散射實驗譜與simnra和corteo兩種模擬軟件得到模擬譜的對比。參考文獻1 world survey of accelerator based analytical techniques， iaea.http：//napc/physics/accel

8、surv/index.html.2 malmqvist k g， radiat. phys. and chem.， 2004， 71： 817.3 dran j c， calligaro t， salomon j， et al.， nucl. instr. and meth b， 2004， 7： 219.4 tesmer j r， nastasi m， barbour j c， et al.， handbook of modern ion beam materials analysis， mrs， pittsburg， usa.1995.5 趙國慶，任熾剛，核分析技術(shù)，原子能出版社，北京，1

9、989.6 楊福家，趙國慶，離子束分析，復(fù)旦大學(xué)出版社，上海，1985.7 楊福家，原子物理學(xué)，上?？茖W(xué)技術(shù)出版社，上海，1985.8g. amsel， nucl. instr. meth.b 201 （2003） 325388.9 mayer m. simnra users guide， technical report ipp 9/133， mpi für plasmaphysik， garching， germany， 1997， http：/home.rzg.mpg.de/mam/index.html.10 f. schiettekatte， nucl.instr. meth. b 266 （2008） 1880.endprint科學(xué)與財富2017