MonteCarlo方法及相關(guān)軟件在實(shí)驗(yàn)核物理中的應(yīng)用

1、Monte Carlo方法及應(yīng)用軟件在實(shí)驗(yàn)核物理中的應(yīng)用,內(nèi)容,Monte Carlo簡(jiǎn)介 核技術(shù)特點(diǎn)及Monte Carlo方法的應(yīng)用 Monte Carlo應(yīng)用軟件 Geant4具體應(yīng)用實(shí)例,將一根長(zhǎng)度為l的針，隨機(jī)的投放在兩條間距為d(dl)的兩條平行線中間，計(jì)算針與兩條平行線相交的概率。 通過解析法，容易得到 針與平行線相交的概率 由此可以得到圓周率,Buffon投針問題,Buffon投針問題,Buffon投針問題,Buffon投針問題,Buffon投針問題,Buffon投針問題,Monte Carlo,蒙特卡羅方法又稱統(tǒng)計(jì)模擬方法，是利用隨機(jī)數(shù)進(jìn)行數(shù)值模擬的方法 可以將蒙特卡羅法看作

2、利用隨機(jī)試驗(yàn)的方法計(jì)算積分，所計(jì)算的積分可看做服從某種分布的密度函數(shù)為f(r)的隨機(jī)變量g(r)的數(shù)學(xué)期望 而通過某種試驗(yàn)得到N個(gè)觀察值r1，r2.rN等，將N個(gè)隨機(jī)變量的值g(r1), g(r2), g(rN)的算術(shù)平均值作為積分估計(jì)值,Monte Carlo,Monte Carlo,根據(jù)概率論知識(shí)，可知,優(yōu)點(diǎn) 對(duì)于具有隨機(jī)性質(zhì)的事件或物理過程其物理意義逼真 受幾何限制小 收斂速度與問題的維數(shù)無關(guān) 誤差容易確定 程序結(jié)構(gòu)簡(jiǎn)單，易于實(shí)現(xiàn) 缺點(diǎn) 收斂速度慢 誤差具有概率性 在粒子輸運(yùn)問題中，計(jì)算結(jié)果與系統(tǒng)大小有關(guān),Monte Carlo,蒙特卡羅算法主要組成部分 概率密度函數(shù)(pdf) 必須給出

3、描述一個(gè)物理系統(tǒng)的一組概率密度函數(shù); 隨機(jī)數(shù)產(chǎn)生器能夠產(chǎn)生在區(qū)間0,1上均勻分布的隨機(jī)數(shù) 抽樣規(guī)則如何從在區(qū)間0,1上均勻分布的隨機(jī)數(shù)出發(fā),隨機(jī)抽取服從給定的pdf的隨機(jī)變量; 模擬結(jié)果記錄記錄一些感興趣的量的模擬結(jié)果 誤差估計(jì)必須確定統(tǒng)計(jì)誤差（或方差）隨模擬次數(shù)以及其它一些量的變化； 減少方差的技術(shù)利用該技術(shù)可減少模擬過程中計(jì)算的次數(shù)； 并行和矢量化可以在先進(jìn)的并行計(jì)算機(jī)上運(yùn)行的有效算法,Monte Carlo,蒙特卡羅算法模擬的主要步驟,Monte Carlo,Monte Carlo,數(shù)值 解法,Monte Carlo,Monte Carlo,Monte Carlo模擬,+,統(tǒng)計(jì)特性 高樣

4、本量,非常有利于Monte Carlo方法的應(yīng)用,在核物理中的應(yīng)用 而其中粒子輸運(yùn)問題是蒙特卡羅模擬的基本核物理過程,Monte Carlo模擬,蒙卡模擬粒子輸運(yùn)問題的主要步驟,Monte Carlo模擬,確定粒子輸運(yùn)過程中涉及的物理過程,蒙卡模擬粒子輸運(yùn)問題的主要步驟 多粒子情況下常用的有字典編輯分支法,Monte Carlo模擬,確定所用的蒙卡技巧（單粒子）,蒙卡模擬粒子輸運(yùn)問題的主要步驟,Monte Carlo模擬,確定粒子狀態(tài)參數(shù)及狀態(tài)序列,蒙卡模擬粒子輸運(yùn)問題的主要步驟 以容易實(shí)現(xiàn)，所用計(jì)算時(shí)間少為標(biāo)準(zhǔn),Monte Carlo模擬,確定粒子輸運(yùn)過程中狀態(tài)分布的抽樣方法,在通過對(duì)粒子輸

5、運(yùn)問題的進(jìn)行模擬的基礎(chǔ)上，考慮實(shí)際情況進(jìn)一步對(duì)其結(jié)果進(jìn)行分析，可以針對(duì)不同應(yīng)用得出不同的物理結(jié)果，如屏蔽問題中的光通量，核輻射探測(cè)器中的探測(cè)效率、能量沉積譜、能量分辨率等。,Monte Carlo模擬,Monte Carlo應(yīng)用軟件,全名MultigroupOak Ridge Stochastic Experiment，是美國橡樹嶺國家實(shí)驗(yàn)室從60年代開始研制的大型、多功能、多群中子-光子耦合輸運(yùn)程序 特點(diǎn) 具有三維幾何能力 具有多種功能，可解決中子、光子、中子-光子耦合輸運(yùn)問題 使用多群截面數(shù)據(jù) 模塊結(jié)構(gòu) 包括幾種有效地減小方差技巧和計(jì)算技巧 程序具有很大的靈活性,Monte Carlo應(yīng)用

6、軟件MORSE,全名Monte Carlo Neutron andPhoton Transport Code，由美國洛斯阿拉莫斯國家實(shí)驗(yàn)室編制的具有目前最高水平的中子-光子輸運(yùn)程序，被稱為“超級(jí)蒙特卡羅程序”。 特點(diǎn) 程序中的幾何可以是三維任意態(tài) 使用精細(xì)的點(diǎn)截面數(shù)據(jù) 可用于中子、光子、中子-光子耦合輸運(yùn)問題，及臨界系統(tǒng)特征值問題 減小方差技巧比較全面 程序的通用性很強(qiáng) 使用Fortran語言,Monte Carlo應(yīng)用軟件MCNP,全名為electron-gamma shower，由美國Stanford Linear Accelerator Center提供，用以模擬在任意幾何中能量從幾ke

7、V到幾TeV的電子-光子簇射過程，目前最新版本為EGS-5 特點(diǎn) 元素和材料介質(zhì)齊全 帶電粒子和光子的輸運(yùn)均采用隨機(jī)游動(dòng)模式進(jìn)行 光子的能量可以從1keV到幾千GeV 反應(yīng)類型非常齊全 使用Mortran語言,Monte Carlo應(yīng)用軟件EGS,全名FLUktuierende Kaskade，是INFN（意大利核物理研究所）和CERN（歐洲粒子物理研究中心）聯(lián)合開發(fā)的一種蒙特卡羅模擬計(jì)算程序，質(zhì)子、電子加速器屏蔽設(shè)計(jì)，量熱計(jì)，活化，計(jì)量學(xué)，探測(cè)器設(shè)計(jì)，宇宙射線，中微子物理及放射治療等 特點(diǎn) 采用微觀模型，物理過程嚴(yán)謹(jǐn) 約60種不同的粒子可以使用FLUKA進(jìn)行模擬，能量范圍寬 使用改進(jìn)的CG（Combinatorial Geometry）軟件包，可處理復(fù)雜幾何結(jié)構(gòu) 針對(duì)于大多數(shù)的應(yīng)用，F(xiàn)LUKA使用者不需要編寫程序 完全免費(fèi),Monte Carlo應(yīng)用軟件FLUKA,全名GEometry ANd Tracking，由上世紀(jì)90年代CERN（歐洲核子研究組織）和KEK（日本高能加速器研究中心）組織了來自歐洲、俄羅斯、加拿大、日本和美國的幾十個(gè)實(shí)驗(yàn)室、高校和研究機(jī)構(gòu)的超過100名科學(xué)家和工程師進(jìn)行合作編譯出來的具有開放結(jié)構(gòu)的蒙卡程序 特點(diǎn)： 采用面向