基于β衰變能譜測量探索確定137Cs衰變分支比的方法

李 冉 張澤濤 秦雨浩 任 贊 張高龍

(1北京航空航天大學物理科學與核能工程學院，北京 100191;2 清華大學工程物理系，北京 100084)

實驗上通過研究原子核的衰變可以建立原子核衰變能級綱圖。衰變綱圖不僅是研究核結構一個十分有用的工具，對它大量的和系統的建立與分析，可以揭示原子核的運動規律；而且對放射性核素的應用也非常重要。它是核物理一個重要的分支-應用核譜學的基礎。β衰變不僅在重核范圍內發生，在全部元素周期表的范圍內都存在β放射性核素。因此，對β衰變的研究非常重要。通過測量原子核的β衰變來確定其能級自旋宇稱和衰變綱圖是近期核物理研究的重要課題之一，通過此方面的研究，也能了解原子核的結構[1,2]。本文基于原子核137Cs的β衰變，從實驗上建立其衰變綱圖。137Cs有兩種β衰變[3]，其中對高能量的β衰變研究較少。本研究中用Si(Li)漂移探測器測量了137Cs的β衰變能譜和內轉換電子譜。利用高能區的低統計計數，通過庫里厄圖和采用更精確的反演數據處理方法能否得到兩β衰變分支的分支比，能否擴展β衰變和內轉換電子譜測量實驗的應用范圍，是本研究的主要目的。進而驗證該技術和方法能否用于原子核結構和衰變性質的研究，也使得學生能夠更加深入地了解原子核有關β衰變和庫里厄圖的有關知識。

1 實驗原理

1.1 β衰變[4]

類比于原子發光的電磁相互作用，β衰變是電子-中微子場與原子核的相互作用，使核子不同狀態之間引起躍遷，發射出電子和中微子，是一種弱相互作用。β躍遷可分為容許躍遷和禁戒躍遷兩種。其中禁戒躍遷又分為一級禁戒躍遷、二級禁戒躍遷等。以ΔI表示衰變前后母核和子核的自旋量子數的變化，即自旋之差：ΔI=Ii-If；以Δπ表示衰變前后母核和子核的宇稱量子數的變化，即宇稱之積：Δπ=Δπi·Δπf，則

在容許躍遷中自旋和宇稱變化滿足

(1)

在一級禁戒躍遷中，

(2)

對于n(n≥2)級禁戒躍遷，有

(3)

利用量子力學的微擾理論，并考慮到庫侖場的影響，可以得到β粒子動量分布的表達式為

(4)

其中，g為弱相互作用常量；Mif為躍遷矩陣元；E為β粒子的動能；Em為β粒子最大動能；p為β粒子的動量；F(Z,E) 為修正因子，在非相對論近似中表示為

在式(6)中，對β衰變取正號，對β+衰變取負號；v為β粒子的速度。本研究中Z為137Cs的核電荷數。137Cs到137Ba 的衰變為β-衰變，故本實驗中x取正號。將式(4)改寫為

(9)

式中，W和W0是以mec2為單位的β粒子總能量及其最大值(量綱為1)，即

(12)

更為簡單判斷躍遷級次的方法是計算logfT1/2的值；其中f為最大衰變能量Em和子核的核電荷數Z的函數，可查閱圖1[4](圖中Z的正負表示β-衰變和β+衰變)；T1/2為母核的半衰期。不同類型和級次的躍遷所對應的 logfT1/2的范圍不同，在實驗上確定Em的值以后，可以計算出logfT1/2的值，進而判斷其躍遷類型和級次。躍遷級次與logfT1/2的對應關系如表1所示。

圖1 函數logf(Z,Em)

躍遷級次logfT1/2超容許躍遷2.9～3.7容許躍遷4.4～6.0一級禁戒(非唯一型)6～9一級禁戒(唯一型)8～10二級禁戒10～13三級禁戒15～18

1.2 衰變綱圖和分支比

原子核的衰變性質一般用衰變綱圖來直觀描述。衰變綱圖是綜合反映某核素放射性衰變的主要特征和數據的示意圖，包含了研究原子核衰變中積累起來的關于核低激發態性質的大量實驗數據。衰變綱圖具有重要的意義和價值。一是在研究原子核運動規律的過程中，需要自旋、宇稱和躍遷級次等數據；二是對于放射性核素應用來說，需要了解某種核素衰變所釋放的α射線、β射線或γ射線能量的具體數值，特別是在核醫學和核輻射防護等實際應用中更為常見，而這些數據都可以在衰變綱圖中找到。

β衰變常常會伴隨γ衰變，這意味著衰變子核一般也具有能級。β衰變的終態可以是滿足某種躍遷級次選擇定則的能級，常常不止一個。這意味著我們測得的β衰變能譜將是兩種衰變路徑能譜的疊加，繪制的庫里厄圖也將不再嚴格滿足直線關系。對于容許躍遷而言，其庫里厄圖將成為一條折線：從高能區向低能區看，每次彎折都相當于追加了一個新的衰變路徑的疊加。利用線性關系，可以將不同路徑分解開，具體方法將在后面的實際數據分析中闡明。利用分解開的庫里厄圖，不僅可以分析各個衰變路徑的特性，還可以測出路徑間總計數的比例，從而得到兩種衰變路徑的比例，標出兩種路徑占全部衰變的比例，也就是分支比。

如果這些路線的躍遷級次不同，那么強度將有幾個數量級的差異：容許躍遷或低級次禁戒躍遷的發生概率比高級次禁戒躍遷高幾個數量級。這樣，在研究低級次躍遷時，完全可以忽略高級次躍遷疊加的效應，在確定了更多關于高級次躍遷的信息后再利用上面的方法研究分支比。后面會看到，我們此次要研究的137Cs 就是這樣具有不同級次禁戒躍遷分支的衰變母核，這意味著其庫里厄圖將不會表現出典型的折線特征。本研究中將對這兩級躍遷分別進行庫里厄圖，利用其數量級差異和γ峰恰好位于兩截止頻率之間的特征，以γ峰和內轉換電子峰為界分離兩個分支，分別完成庫里厄圖，從而實現類似于上述容許躍遷分離的效果。

在此前的工作[5]中，利用內轉換系數方法建立了包括自旋、宇稱和能量信息的衰變綱圖，這一綱圖是下面分析高能β能譜的基礎。基本建立的137Cs的衰變綱圖如圖2所示，其中缺少分支比。

圖2 由實驗數據建立的137Cs衰變綱圖

2 實驗裝置和實驗數據

本實驗在中國原子能科學研究院進行，實驗所用137Cs放射源的活度為3.77×103Bq；Si(Li)漂移探測器為美國CANBERRA ESLB系列，尺寸為300mm2×3mm；主放大器型號為ORTEC 672，前置放大器為主放大器自帶。

圖3 電子學線路圖

Si(Li)漂移探測器產生的信號經過前置放大器和主放大器的成形、放大之后送入數據獲取系統[6]，通過多道脈沖幅度分析器轉化為具體的數字信息，再通過計算機軟件的處理，得到所測量的能譜圖像，儀器的連接示意圖如圖3所示。

3 實驗數據分析和處理

137Cs的β衰變為β-衰變，其衰變子核為137Ba，137Cs經過β衰變到137Ba*的激發態后，會產生γ躍遷或發射內轉換電子。測得衰變能譜如圖4所示，從能譜中可以判斷低能部分最明顯的峰為β衰變電子峰，高能部分(5400～6000道)的3個峰為內轉換電子峰。利用已知的3個內轉換電子峰的能量數據[3]，對所測能譜進行了能量刻度，得到能量E(keV)與道址(Ch)的關系。

根據3個內轉換電子峰的半高寬和所對應的能量，可計算出探測器的能量分辨率，如表2所示。

表2 內轉換電子峰的半高寬和對應的能量分辨率

參考文獻[7]中所制作的Si(Li)電子譜儀對于975. 62keV 能量的內轉換電子譜的半高寬為2.07keV，則能量分辨率為0.21%。本實驗所用Si(Li)漂移探測器對624.22keV、655.67keV和660.36keV內轉換電子的能量分辨率分別為0.28%、0.23%和0.10%，兩者相比，可以看出，本實驗所用探測器的能量分辨率也是較好的。

3.1 高能庫里厄圖

我們所畫出的衰變綱圖(圖2)并不完整，該綱圖涉及兩條衰變分支，必須補全選擇各個分支的概率。對于一般的容許躍遷，只利用折線形的庫里厄圖就可以拆分各個衰變階段。然而，本實驗中，由于受到實驗中計數時間不能太長，以及兩種躍遷是不同級的禁戒躍遷等問題的制約，不能用傳統的方法進行分離。這里，利用之前工作中得到的衰變綱圖，可以將數據拆分成兩張庫里厄圖，從而將兩者分開、并能更加精確地分析低計數率的高能電子譜。

利用圖2中繪制的衰變綱圖[4]，可以看出直接到達基態的衰變自旋改變ΔI=-2，宇稱不改變，這意味著該衰變為非唯一型二級禁戒躍遷，強度會比另一條β衰變低一到兩個數量級；對于衰變到基態的分支，僅僅利用已經得到的實驗數據和上面分析得到的能級參量，就可以比較容易地從能譜中找到這種β衰變的證據，并計算出其分支比。

圖4 137Cs 的β衰變能譜

圖5 二級禁戒躍遷庫里厄圖

為了避免γ躍遷和內轉換電子的干擾，選取大于662.5keV的數據，這部分扣除本底后僅包括綱圖中直接躍遷到基態這一路徑的β粒子，這一躍遷是二級禁戒躍遷，因此可以利用二級禁戒修正的庫里厄圖。由于這一區域計數率相當低，我們采取將100道并為1道的方法降低統計誤差。

可以看出這一段的庫里厄圖接近線性關系，得到其中的線性系數k=(1.02±0.17)×10-5，這樣分布N(p)=Fp2S2k2(E0-E)2中的所有參量全部確定，可以數值積分得到衰變路徑的總計數

(13)

式中同樣只計入了統計誤差。這樣就得到了高能β射線的總計數。在工作[5]中，已經對低能β躍遷進行了一級禁戒居里描繪，得到低能衰變路徑的總計數為

(14)

這一結果只計入了計數誤差。

3.2 分支比

根據得到的兩條β衰變路徑的事件數，可以得到高能β衰變(137Cs基態-137Ba基態)粒子所占分支比為(4.4±0.9)%，與參考值(5.3±0.2)%[8]相比偏小，這里可能存在高能粒子相對比低能粒子更容易直接穿過探測器而不被吸收，探測器在高能區域的吸收率低于低能區域導致的系統誤差。

完整的衰變綱圖如圖6所示，其中標明了各分支的分支比。

圖6 完整的137Cs衰變綱圖

根據圖6，直線與橫軸的交點值為(1022±274)keV，這里均只考慮統計誤差，文獻[8]中給出的參考值為(1175.63±0.17)keV。(1022±274)keV在置信區間內，因此高能β能譜的庫里厄圖可以視為統計誤差允許范圍內的合理結果，得出的分支比在統計誤差范圍內應當是可信的。

另一方面，已知137Cs的半衰期T1/2為30.08a[3]，其衰變子核Ba的原子序數為56，查閱圖1可知1.6

4 結論

在此研究的實驗和數據分析中，我們成功地得到了137Cs的衰變綱圖，利用137Cs和137Ba對應能級的能量和自旋宇稱研究了低計數率的高能區域，得到了兩個衰變分支的分支比。得出結論：137Cs可以發生兩種β衰變，分別到達137Ba的基態3/2+和激發態11/2-，分支比分別約為(95.6±0.9) %和(4.4±0.9) %，參考值基本在這一結論統計誤差的置信區間，因此該結果基本具備在統計誤差內的可信度。

依托二級禁戒躍遷庫里厄圖還可以推出的另一個結論：高能β射線最大能量在(1022±274)keV，同樣在置信區間內，驗證了我們采取的這種利用參數和低計數率實驗數據進行二級禁戒躍遷庫里厄圖的方法仍然是適用的，因此，這種方法是合理的，可以用來繼續分析其他衰變問題的數據。

在原子核結構中，通常用高純鍺探測器只測量γ射線，研究γ躍遷級次。但是，只測量γ衰變的能譜不能準確確定γ躍遷的級次和推斷原子核的自旋宇稱。我們同時測量了β和γ躍遷級次，聯合確定躍遷級次和自旋宇稱，確定β衰變類型，從而能夠準確構造修正函數，分析超低計數率的高能β射線譜數據。這樣，配合內轉換系數分析方法，我們僅需要較短時間的測量數據即可完全確定衰變綱圖，無需長時間測量高能β譜，更適合開發為教學實驗。