欒天成, 郝景萌

(廣州大學(xué) 物理與材料科學(xué)學(xué)院， 廣東 廣州 510006)

0 引 言

伽馬射線暴(簡稱伽馬暴)是宇宙中最劇烈的能量釋放現(xiàn)象，是一種來自宇宙中的伽馬射線在短時間內(nèi)迅速增強(qiáng)隨后又減弱的現(xiàn)象.關(guān)于伽馬暴的起源以及其形成的物理機(jī)制，從20世紀(jì)60年代被發(fā)現(xiàn)至今，一直都是天體物理學(xué)中的熱門話題.根據(jù)其持續(xù)的時間長短，伽馬暴可以被分為兩種類型：長暴(持續(xù)時間大于2 s)和短暴(持續(xù)時間小于2 s).近年來，有越來越多的觀測證據(jù)表明[1-3]，長暴是來源于大質(zhì)量恒星塌縮的產(chǎn)物[4-5]，而短暴是源自于致密雙星系統(tǒng)(雙中子星或一個中子星一個黑洞)的合并[6].為了能夠更好地理解伽馬暴的前身星系統(tǒng)以及它的中心引擎，對于產(chǎn)生伽馬暴的噴流張角和能量結(jié)構(gòu)的限制是至關(guān)重要的[7-8].在之前的研究中，人們更傾向于假設(shè)伽馬暴噴流能量是存在著特殊的角結(jié)構(gòu)的，而并非是最簡單的各向同性輻射[9-12].特別地，Zhang等[13]考慮了一個具有高斯分布的能量角結(jié)構(gòu)模型并嘗試去解釋伽馬暴的觀測性質(zhì).Pescalli等[14]探討了具有冪律分布結(jié)構(gòu)的伽馬暴噴流對其光度函數(shù)的影響.Salafia等[15]還在其工作中進(jìn)一步探索了這兩種不同角結(jié)構(gòu)模型的性質(zhì).基于觀測值EX/Eγ(各向同性等值的早期X射線余輝能量與伽馬能段能量的比值)，以及長暴的光度函數(shù)，Beniamini等[16]探討了伽馬暴的各種噴流角結(jié)構(gòu)模型.另外，引力波與伽馬暴的聯(lián)合探測也被用來限制不同的噴流角結(jié)構(gòu)模型[17].總的說來，以上提到的高斯、冪律分布的噴流能量結(jié)構(gòu)模型，加上標(biāo)準(zhǔn)均勻分布模型目前正在被廣泛地研究與討論，盡管到目前為止,人們對于最適合的噴流結(jié)構(gòu)模型及其對于伽馬暴各種觀測特征的潛在影響，仍然沒有達(dá)成共識.

本文研究了伽馬暴噴流的各種能量角結(jié)構(gòu)模型對其能量及紅移分布的影響，再利用伽馬暴的可用觀測數(shù)據(jù)并結(jié)合蒙特卡洛模擬的方法，來嘗試限制噴流角結(jié)構(gòu)模型.因為相較于短暴，長暴有著更大并且更完備的樣本，本文的討論僅限定于長暴.本文第一節(jié)中，介紹不同的噴流結(jié)構(gòu)模型和所采用的蒙特卡洛模擬方法；在第二節(jié)中，給出模擬的結(jié)果；在第三節(jié)中，進(jìn)行總結(jié)與討論.本文中所采用的宇宙學(xué)參數(shù)是H0=71 km s-1Mpc-1,ΩM=0.3,ΩΛ=0.7.

1 噴流模型和蒙特卡洛模擬

1.1 噴流模型

在本工作中，將考慮三種之前被廣泛研究和討論的噴流角結(jié)構(gòu)模型.第一種是標(biāo)準(zhǔn)均勻分布噴流模型[9]：該模型假設(shè)單位立體角內(nèi)的能量ε(θ)在噴流錐形張角內(nèi)部是一個常數(shù)，而在張角的外部則為0.其公式如下：

(1)

其中，θj是噴流的半偏角大小.之前的研究一般將該角度的大小限制在5°～15°之間[18-20].本文中，為簡單起見，只選取3個固定的角度值進(jìn)行研究，即θj=5°、15°和30°.另兩種模型就是所謂的有結(jié)構(gòu)的噴流模型，這類模型中的噴流都具有固有的各向異性特征，而伽馬暴的變化則主要取決于不同的觀測角度θv.對于冪律分布和高斯分布噴流模型，ε(θ)作為隨噴流軸夾角變化的函數(shù)，分別為[21]

(2)

和

ε(θ)=ε0e-(θ/θ0)2

(3)

其中,θ0代表著噴流錐的典型半偏角大小，在本文中θ0取值設(shè)為2°.特別地，對于冪率分布噴流模型，考慮a=2和4這兩個不同的值.因此，伽馬暴噴流的總能量Ej可以通過在所有立體角內(nèi)對ε(θ)積分得到：

(4)

伽馬暴的各向同性總能量Eiso則可以簡單地通過下式的計算得到：

Eiso=4πε(θv)

(5)

1.2 蒙特卡洛模擬

為了更好地說明伽馬暴的噴流結(jié)構(gòu)模型是如何影響其能量與紅移分布的，利用蒙特卡洛模擬方法生成了一個長暴的模擬樣本.每一個通過模擬產(chǎn)生的長暴都以一組參數(shù)標(biāo)識，即紅移z、總能量Ej和觀測角度θv.這三個參數(shù)分別根據(jù)以下假設(shè)通過隨機(jī)過程生成：

(1)伽馬暴的紅移由一個概率分布函數(shù)p(z)給出.該函數(shù)可以通過對長暴的固有爆發(fā)率在紅移z=0～10區(qū)間進(jìn)行歸一化得到，其形式為

(6)

其中,Az是歸一化常數(shù)，dV/dz是在單位紅移內(nèi)的共動體積元，而Θ(Zth,z)則表示金屬豐度在閾值Zth以下的恒星質(zhì)量占比[22]，其形式為

(7)

(8)

(2)總能量可以從一個對數(shù)正態(tài)分布的概率分布函數(shù)得到，其形式為

(9)

其中，AE是歸一化常數(shù).上式中正態(tài)分布的中心值Ec和標(biāo)準(zhǔn)差σEj是自由參數(shù)，它們的值可以通過與觀測之間的比較來限制.

(3)分配給每一個事件一個觀測角度，該角度可以通過概率分布函數(shù)p(θv)dθv=sin(θv)dθv在范圍[0～π/2]內(nèi)得到.對于均勻分布模型的情況，如果觀測角度θv>θj，那么這個伽馬暴將無法被探測到.但是對于有結(jié)構(gòu)的噴流模型而言，噴流并不存在一個明確的邊界.

logEpeak=-29.6+0.61logEiso

(10)

(11)

其中,dL是光度距離，k(z)是k改正因子，可以通過下式計算：

(12)

其中,N(E)代表觀測者參考系中的光子數(shù)能譜.

1.3 一致性檢測

1.3.1 觀測樣本的選取

截至2018年7月底，Swift-BAT共計探測到了1 247個伽馬暴(從GRB041217到GRB180728)[32]，從中選取628個通量大于閾值F[15-150]≥10-6erg cm-2的長暴.在這些源中，有271個長暴是有紅移探測的，本文將這271個長暴作為對比與分析的觀測樣本.

1.3.2 模型與觀測的一致性檢驗

模型有兩個自由參數(shù)，即Ec和σEj，用以表征噴流總能量的對數(shù)正態(tài)分布.對于每一對所選取的自由參數(shù)，運(yùn)行一遍模擬程序以生成由104個“被探測到”的伽馬暴所組成的一組模擬樣本.然后，與Kanaan等[33]的方法類似，先計算模擬伽馬暴樣本的通量分布，并通過與Swift-BAT觀測到的真實樣本進(jìn)行比較來限制不同噴流模型下的這兩個自由參數(shù).模擬樣本與觀測樣本的通量分布一致性可以通過2檢驗來判定：

(13)

其中,fi,obs和fi,sim分別代表在第i個區(qū)間塊中的觀測與模擬樣本通量的概率.在不同的自由參數(shù)組合下，通過不斷重復(fù)上述蒙特卡洛模擬步驟，最終找到2值最小的一對自由參數(shù)，而這對參數(shù)值便代表了最優(yōu)的參數(shù)選擇.

2 結(jié) 果

表1給出了與觀測樣本的通量分布符合最好的自由參數(shù)組.對于不同的噴流結(jié)構(gòu)模型，能量中心值Ec的差別并不是很大，基本所有模型都處在50附近，除了a=2的冪律分布模型要偏大約2個數(shù)量級.這些結(jié)果表明,除了a=2的冪律分布模型，伽馬暴的噴流總能量可能要比直接通過觀測推導(dǎo)出的各向同性能量少2個數(shù)量級以上.

表1 卡方最優(yōu)參數(shù)列表

Table 1 Best parameters constrained by 2

表1 卡方最優(yōu)參數(shù)列表

ModellogEcσEj 2Const_5dgre49.791.790.129Const_15dgre50.791.930.127Const_30dgre50.532.540.132Gaussian49.611.930.118PL_2a52.711.550.131PL_4a49.861.870.113

在確定最優(yōu)自由參數(shù)組之后，可以進(jìn)一步得到各種噴流模型下模擬伽馬暴樣本的紅移和各向同性的總能量分布，并與觀測進(jìn)行比較.圖1給出了均勻分布模型下伽馬暴各向同性能量Eiso的分布與觀測的比較.其中噴流張角分別取值5°、15°和30°，并以實線、虛線和點(diǎn)虛線表示.然后利用Kolmogorov-Smirnov(K-S)測試來檢驗?zāi)M與觀測樣本之間的一致性.結(jié)果發(fā)現(xiàn)，在3個角度中，只有當(dāng)θj的取值為30°的時候，K-S 檢驗的結(jié)果給出P>0.05，這意味著Eiso的分布與觀測比較一致，而對于5°和15°的情況，K-S 檢驗只給出P<0.05.高斯分布模型和冪律分布模型與觀測的比較分別見圖2～圖3.對所有模型K-S檢驗的結(jié)果都小于0.05.

圖3 冪律噴流分布模型模擬的各向同性能量概率密度分布圖

圖2 高斯噴流分布模型模擬的各向同性能量概率密度分布圖

圖1 均勻噴流分布模型模擬的各向同性能量概率密度分布圖

對于紅移而言，三種噴流張角結(jié)構(gòu)模型(常數(shù)、高斯和冪率)的紅移與觀測紅移的累積分布分別見圖4～圖6，其中對應(yīng)的與觀測紅移分布之間的K-S檢驗P值已在各自圖中右下角標(biāo)出.無論是對于均勻分布模型，還是對于高斯分布和冪率分布模型，給出的結(jié)果區(qū)別都不大，K-S檢驗都得出了一個相當(dāng)?shù)偷母怕?P<0.05).所有模型都低估了低紅移伽馬暴的數(shù)量.

圖4 均勻噴流分布模型模擬的紅移累積分布圖

圖5 高斯噴流分布模型模擬的紅移累積分布圖

圖6 冪律噴流分布模型模擬的紅移累積分布圖

3 總結(jié)與討論

了解伽馬暴的噴流張角和能量結(jié)構(gòu)對于進(jìn)一步理解伽馬暴的前身星系統(tǒng)以及它的中心引擎是至關(guān)重要的.本文利用蒙特卡洛模擬的方法，研究了均勻分布、高斯分布及冪律分布三種伽馬暴噴流的能量角結(jié)構(gòu)模型對其各向同性能量及紅移分布的影響.首先，通過模擬樣本與觀測數(shù)據(jù)的通量分布的比較，限制了不同模型下假設(shè)為對數(shù)正態(tài)分布的伽馬暴噴流總能量分布的特征參數(shù).發(fā)現(xiàn)能量中心值分布在49.790.05).然而對于紅移分布，所有模型的結(jié)果都比較相似，并且與觀測符合得不是很好.模擬與觀測的區(qū)別主要來自于低紅移端，說明所有模型都低估了低紅移伽馬暴的數(shù)量.總之，目前的觀測數(shù)據(jù)可能還不足以區(qū)分三種噴流模型.

模擬結(jié)果與觀測的差異可能來自于以下兩個方面：①模擬所采用的模型比較簡單.一方面，模型中很多參數(shù)都被設(shè)為了定值，例如，金屬豐度閾值等；另一方面，為了簡化計算，忽略了噴流的相對論效應(yīng)、伽馬暴光度函數(shù)和噴流張角等可能存在的紅移演化；②簡化了Swift-BAT探測器的探測標(biāo)準(zhǔn).在本文中，采用了一個單一的通量閾值作為伽馬暴能否被探測到的標(biāo)準(zhǔn).而事實上，Swift-BAT的伽馬暴觸發(fā)機(jī)制非常復(fù)雜.因此，低通量伽馬暴的截斷有可能會導(dǎo)致低紅移伽馬暴的低估.

在后續(xù)工作中，筆者會持續(xù)改進(jìn)工作，例如，考慮一個更復(fù)雜的模型，采用更好的閾值標(biāo)準(zhǔn)，使用馬爾科夫鏈?zhǔn)厦商乜?MCMC)方法替代蒙特卡洛方法直接從離散觀測樣本分布中采樣，還有增加觀測數(shù)據(jù)上的限制，進(jìn)一步探究各種因素對噴流模型的影響，從而能對模型進(jìn)行更好的限制.