塑料閃爍體的閃爍時(shí)間特性

黃展常 楊建倫

（中國工程物理研究院 核物理與化學(xué)研究所 綿陽 621900）

塑料閃爍體的閃爍時(shí)間特性

黃展常 楊建倫

（中國工程物理研究院 核物理與化學(xué)研究所 綿陽 621900）

采用137Cs激發(fā)塑料閃爍體發(fā)光，應(yīng)用單光子計(jì)數(shù)法測(cè)量塑料閃爍體發(fā)光衰減曲線。根據(jù)測(cè)量原理搭建了實(shí)驗(yàn)平臺(tái)，測(cè)量了ST401和EJ232兩種閃爍體的閃爍時(shí)間特性。實(shí)驗(yàn)結(jié)果用雙指數(shù)擬合法進(jìn)行分析可得，ST401發(fā)光衰減常數(shù)為2.9 ns和50 ns，EJ232發(fā)光衰減常數(shù)為1.6 ns和30 ns。實(shí)驗(yàn)平臺(tái)所測(cè)幅度動(dòng)態(tài)范圍大于104，時(shí)間分辨小于0.8 ns。從統(tǒng)計(jì)學(xué)分析單光子計(jì)數(shù)法的物理過程，并采用蒙特卡羅方法模擬整個(gè)探測(cè)過程，探討單光子計(jì)數(shù)法的適用條件。實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果相洽。若終止探測(cè)器平均探測(cè)到的光子數(shù)不大于0.1，則多道測(cè)量結(jié)果可以近似為塑料閃爍體發(fā)光衰減曲線。

塑料閃爍體，閃爍時(shí)間特性，單光子計(jì)數(shù)法，蒙特卡羅方法

與碘化鈉閃爍體相比，塑料閃爍體具有發(fā)光時(shí)間短、發(fā)光強(qiáng)度低的特點(diǎn)。其在強(qiáng)脈沖伽馬場(chǎng)或者強(qiáng)脈沖中子場(chǎng)測(cè)量中有較廣泛的應(yīng)用[1–2]。例如，在氘氘聚變反應(yīng)中，采用飛行時(shí)間法測(cè)量中子能譜時(shí)，有時(shí)候會(huì)關(guān)注初級(jí)中子的散射中子信號(hào)，這就需要了解塑料閃爍體的發(fā)光時(shí)間特性，尤其是慢成分，以便選用合適的閃爍體進(jìn)行脈沖信號(hào)測(cè)量降低初級(jí)中子激發(fā)的慢成分熒光對(duì)散射中子測(cè)量的影響[3]。目前塑料閃爍體的發(fā)光時(shí)間特性測(cè)量方法較多，例如采用激光激發(fā)塑料閃爍體，用高時(shí)間分辨的條紋相機(jī)記錄閃爍體的發(fā)光過程，該方法要求較窄的激光脈寬，測(cè)量結(jié)果受閃爍體的尺寸影響較大[4]。由于單次激發(fā)的熒光數(shù)目有限，該方法可測(cè)動(dòng)態(tài)范圍低于10，主要用于研究快成分。采用激光激發(fā)閃爍體，用高時(shí)間采樣率的示波器記錄探測(cè)器的輸出波形，該方法可測(cè)動(dòng)態(tài)范圍低于10，主要用于研究快成分[5]。采用單光子計(jì)數(shù)法，常用超短脈沖激光或放射源激發(fā)閃爍體，通過統(tǒng)計(jì)大量閃爍體發(fā)光事件實(shí)現(xiàn)時(shí)間的不連續(xù)測(cè)量，由多道統(tǒng)計(jì)譜表征閃爍體發(fā)光衰減曲線[6–9]。該方法應(yīng)用較廣泛，但是較少研究其適用條件，特別是對(duì)減光條件要求的討論較少。文獻(xiàn)[10]采用單光子計(jì)數(shù)法測(cè)量了ST401，主要研究其快成分，圖數(shù)據(jù)顯示動(dòng)態(tài)范圍低于103。本文通過數(shù)值模擬與實(shí)驗(yàn)結(jié)果分析研究塑料閃爍體的閃爍時(shí)間特性，關(guān)注其慢成分的時(shí)間特性，并得到單光子計(jì)數(shù)法應(yīng)用于塑料閃爍體時(shí)間特性研究的適用條件。

1 實(shí)驗(yàn)裝置

如圖1所示，用137Cs源激發(fā)塑料閃爍體ST401。閃爍體緊貼PMT1光電倍增管的信號(hào)輸入窗。由于PMT1輸出信號(hào)幅度大部分大于1 V，為保護(hù)恒比定時(shí)器CFD935-1，因此把輸出信號(hào)進(jìn)行適當(dāng)衰減后輸入CFD935-1。再把其輸出信號(hào)輸入到時(shí)幅轉(zhuǎn)換器TAC567中，作為計(jì)時(shí)開始信號(hào)。CFD935-1的閾值調(diào)至熱噪聲幅度以上，以降低熱噪聲開啟TAC567的概率。PMT2離開閃爍體一定的距離，滿足單光子收集條件。把其輸出信號(hào)輸入恒比定時(shí)器CFD935-2。為把發(fā)光衰減峰移至?xí)r間幅度轉(zhuǎn)換TAC (Time-to-Amplitude Converter)的最小分辨時(shí)間外，CFD935-2的輸出比CFD935-1多加2 m電纜后輸入到時(shí)幅轉(zhuǎn)換器中，作為計(jì)時(shí)終止信號(hào)。通過大量發(fā)光事件統(tǒng)計(jì)，在多道分析器上得到閃爍體的發(fā)光衰減統(tǒng)計(jì)譜。

實(shí)驗(yàn)前通過多道能譜統(tǒng)計(jì)發(fā)現(xiàn)PMT1和PMT2的熱噪聲幅度在10 mV以下的計(jì)數(shù)率很大，但是大于等于100 mV的平均計(jì)數(shù)率為每秒2個(gè)，熱噪聲幅度大于200 mV的平均計(jì)數(shù)率為5分鐘1個(gè)。無需經(jīng)過快放即可滿足恒比定時(shí)器的輸入幅度要求。實(shí)驗(yàn)時(shí)，衰減器的倍率為10倍，提供起始信號(hào)的恒比定時(shí)器的閾值設(shè)為30 mV，則對(duì)應(yīng)探測(cè)器輸出幅度為300 mV，降低了熱噪聲開門的概率。提供終止信號(hào)的恒比定時(shí)器的閾值設(shè)為100 mV。因此無需制冷裝置。

圖1 實(shí)驗(yàn)框圖Fig.1 Schematic setup of the experiment.

2 測(cè)量原理

假設(shè)伽瑪射線與閃爍體作用后，閃爍體產(chǎn)生N0個(gè)熒光光子。停止探測(cè)器對(duì)熒光光子的總探測(cè)效率為ε，則停止探測(cè)器平均探測(cè)到的熒光數(shù)目為N0ε，記m=N0ε。閃爍體的發(fā)光衰減歸一化曲線為f(t)，f(t)從0到無窮的積分值為1。F(t)等于f(t)從0到t的積分。那么，在[0, t]區(qū)間停止探測(cè)器探測(cè)到0個(gè)光子的概率為：exp(-N0F(t)ε)，在[t, t+dt]區(qū)間探測(cè)到非零個(gè)光子的概率為：N0f(t)εdt。

伽瑪射線與閃爍體作用一次，熒光光子被多道記錄在t-t+dt區(qū)間的概率為：

式中，1-e?m為歸一化系數(shù)；g(t)為多道歸一化譜。

假設(shè)G(t)等于g(t)從0到t的積分，則式(1)兩邊積分整理得：

對(duì)式(2)兩邊求導(dǎo)，可得：

對(duì)式(3)的e?m進(jìn)行泰勒級(jí)數(shù)展開得到：

由式(4)可知，當(dāng)m較小時(shí)，f(t)趨向于g(t)，即由多道統(tǒng)計(jì)的歸一化譜g(t)就可以近似表征為閃爍體的發(fā)光衰減譜f(t)。

伽瑪與閃爍體作用后產(chǎn)生的光子數(shù)N0服從一定的分布。例如本次實(shí)驗(yàn)中，137Cs與塑料閃爍體作用，N0近似服從能量沉積譜分布，因此m的值也服從能量沉積譜分布。但是m具有一定的統(tǒng)計(jì)特征，后面討論與實(shí)驗(yàn)的m是指m所服從的分布平均值。

假設(shè)閃爍的熒光光子服從一個(gè)已知的單指數(shù)分布f(t)，采用蒙特卡羅方法對(duì)探測(cè)過程進(jìn)行模擬，得到不同m值時(shí)，多道統(tǒng)計(jì)譜g(t)與f(t)之間的關(guān)系，如圖2所示。

由圖2可知，當(dāng)m≥0.5時(shí)，歸一化的多道統(tǒng)計(jì)譜g(t)與閃爍體發(fā)光衰減曲線f(t)有較明顯的差別。當(dāng)m≤0.1時(shí)，多道統(tǒng)計(jì)譜g(t)與閃爍體發(fā)光衰減曲線f(t)比較接近，可以采用歸一化后的多道統(tǒng)計(jì)譜g(t)來表征閃爍體發(fā)光衰減曲線f(t)。

綜上可得，當(dāng)m值較大時(shí)，低道計(jì)數(shù)升高，高道計(jì)數(shù)降低，表現(xiàn)為衰減快。當(dāng)m值較小時(shí)，需要測(cè)量的時(shí)間長，穩(wěn)定性要求高。由m值即可預(yù)估在單光子測(cè)量實(shí)驗(yàn)中需要進(jìn)行減光的倍率。假設(shè)取m=0.1，若一次伽瑪事件平均產(chǎn)生103個(gè)熒光，那么需要減光倍率應(yīng)該大于104倍（包含熒光從閃爍體的透出率、幾何減光和濾片減光的總和）。在這種減光條件下，多道統(tǒng)計(jì)譜可近似為發(fā)光衰減時(shí)間譜。m值的提出為減光程度提供了一個(gè)方向，而不是盲目地減光，在保證結(jié)果可靠的前提下又能保證一定的符合計(jì)數(shù)率。

圖2 不同m值多道統(tǒng)計(jì)譜與閃爍體發(fā)光衰減曲線Fig.2 Comparison of scintillator decay curve with simulated multi-channel analyzer outputs for various values of parameter m.

3 實(shí)驗(yàn)結(jié)果及討論

據(jù)圖1布局，把尺寸為?40 mm×20 mm的ST401緊貼PMT1進(jìn)行實(shí)驗(yàn)。通過調(diào)節(jié)PMT2與塑料閃爍體之間的距離以及通過在PMT2輸入窗面緊貼不同直徑的小孔的光闌來改變PMT2接收光子的立體角，進(jìn)而改變其收集光子個(gè)數(shù)。圖3(a)是較典型的多道統(tǒng)計(jì)波形，由圖3(a)可知，第0–20道無計(jì)數(shù)，這是由于TAC自身固有時(shí)間分辨造成。圖3(a)的總計(jì)數(shù)為650 110，TAC開啟次數(shù)估計(jì)為650110/0.05。0.05為符合率m。由前面已知PMT2噪聲計(jì)數(shù)率為每秒2個(gè)，因此噪聲的符合計(jì)數(shù)為：2s-1×100 ns×650110/0.05=3。由于高道址的統(tǒng)計(jì)漲落較大，主要記錄慢成分的時(shí)間行為，可以適當(dāng)?shù)貭奚鼤r(shí)間分辨，因此對(duì)多道統(tǒng)計(jì)譜400道以后的數(shù)據(jù)進(jìn)行合并，每10道取均值并把該均值作為這10道中的第5道的計(jì)數(shù)值，得到如圖3(b)所示。

圖3 ST401塑料閃爍體發(fā)光的多道統(tǒng)計(jì)譜(a) 原始數(shù)據(jù)，(b) 處理數(shù)據(jù)Fig.3 ASPEC927 statistical spectra for ST401 scintillator. (a) Original data, (b) Re-binned data

實(shí)驗(yàn)過程中，m的估計(jì)值等于start信號(hào)與stop信號(hào)符合計(jì)數(shù)率除以start信號(hào)的計(jì)數(shù)率。不同m值時(shí)的多道歸一化統(tǒng)計(jì)譜如圖4所示。實(shí)驗(yàn)過程中，每道的時(shí)間間隔為0.1 ns。

由圖4可知，當(dāng)m=0.5時(shí)，其發(fā)光衰減峰比m=0.05時(shí)的發(fā)光衰減峰尖銳。對(duì)于下降沿，與m=0.05相比，m=0.5的統(tǒng)計(jì)譜表現(xiàn)為低道址的計(jì)數(shù)值升高，高道址的計(jì)數(shù)值降低。這與§2的數(shù)值模擬結(jié)果相洽。由數(shù)值模擬結(jié)果可知，當(dāng)m=0.05時(shí)，采用單光子計(jì)數(shù)法得到的多道統(tǒng)計(jì)譜可以近似為閃爍體的發(fā)光衰減時(shí)間譜。針對(duì)m=0.05的多道統(tǒng)計(jì)譜，把峰位設(shè)為0時(shí)刻并對(duì)峰值歸一化后，進(jìn)行y–t雙指數(shù)擬合[11]，得到y(tǒng)=0.999× exp(-t/2.9)+0.001×exp(-t/50.0)。即ST401的快成分發(fā)光衰減常數(shù)為2.9 ns，與文獻(xiàn)[10]結(jié)果相洽。雙指數(shù)擬合過程：對(duì)于如圖3(b)譜的下降沿，先采用Origin進(jìn)行單指數(shù)擬合得到快成分衰減常數(shù)的范圍，再用Origin進(jìn)行雙指數(shù)擬合并輸入快成分的取值范圍得到快慢成分衰減常數(shù)的取值。通過MATLAB對(duì)快慢成分的衰減常數(shù)做變動(dòng)調(diào)節(jié)，進(jìn)行試解，當(dāng)擬合結(jié)果與譜線的相關(guān)系數(shù)為0.99時(shí)停止擬合。

圖4 不同m值時(shí)ST401塑料閃爍體發(fā)光的多道歸一化統(tǒng)計(jì)譜Fig.4 ASPEC-927 normalized statistical spectra of different m for ST401 scintillator.

在m=0.05的條件下，把ST401塑料閃爍體換成尺寸為?50 mm×20 mm的EJ232塑料閃爍體進(jìn)行實(shí)驗(yàn)，結(jié)果如圖5所示。把峰位設(shè)為0時(shí)刻并對(duì)峰值歸一化后，進(jìn)行y–t雙指數(shù)擬合，y=0.98× exp(-t/1.6)+0.02×exp(-t/30.0)。即EJ232的快成分發(fā)光衰減常數(shù)為1.6 ns，與文獻(xiàn)[12]結(jié)果相洽。

圖5 EJ232塑料閃爍體發(fā)光的多道統(tǒng)計(jì)譜(a) 原始數(shù)據(jù)，(b) 處理數(shù)據(jù)Fig.5 ASPEC-927 statistical spectra for EJ232 scintillator. (a) Original data, (b) Re-binned data

4 系統(tǒng)時(shí)間分辨標(biāo)定

保持圖1的連線不變，用同步機(jī)DG645觸發(fā)皮秒激光器輸出激光信號(hào)。用光纖分路盒把光信號(hào)一分為二，一路輸出靠近PMT1信號(hào)窗，模擬閃爍體發(fā)光光子輸入PMT1；另一路輸出適當(dāng)減光靠近PMT2信號(hào)窗，模擬閃爍體發(fā)光光子輸入PMT2。把PMT1的輸出信號(hào)輸入恒比定時(shí)器為TAC提供開啟信號(hào)。把PMT2的輸出信號(hào)采用衰減器進(jìn)行適當(dāng)衰減使之與實(shí)驗(yàn)中采用放射源時(shí)輸出信號(hào)的平均幅度相近，再輸入恒比定時(shí)器為TAC提供終止信號(hào)。其他實(shí)驗(yàn)條件不變。多道統(tǒng)計(jì)的結(jié)果如圖6所示。由圖6可知，時(shí)間分辨統(tǒng)計(jì)峰的半高寬為4道。在此次標(biāo)定實(shí)驗(yàn)中，標(biāo)定值為每道0.2ns，因此時(shí)間分辨統(tǒng)計(jì)峰的半高寬為0.8ns。

圖6 系統(tǒng)時(shí)間分辨統(tǒng)計(jì)譜Fig.6 ASPEC927 statistical spectrum of system time resolution.

5 結(jié)語

本文通過數(shù)值模擬分析了單光子計(jì)數(shù)法的可行性，提出了用多道統(tǒng)計(jì)譜表征閃爍體發(fā)光衰減曲線的條件。在分析過程中引入了m值，解決了實(shí)驗(yàn)中減光倍率和計(jì)數(shù)率的協(xié)調(diào)問題。實(shí)驗(yàn)研究了ST401和EJ232塑料閃爍體的發(fā)光時(shí)間特性。實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果相洽，實(shí)驗(yàn)平臺(tái)運(yùn)行穩(wěn)定可靠，為研究塑料閃爍體快、慢成分熒光時(shí)間特性提供了技術(shù)參考。

致謝 感謝朱學(xué)斌、李忠寶、鄭普等同事提供的硬件支持。

1 Lindl J D, Amendt P, Berger R L, et al. The physics basis for ignition using indirect-drive targets on the national ignition facility[J]. Physics of Plasmas, 2004, 11(2): 339-491

2 胡孟春, 葉文英, 周殿忠, 等. 微秒級(jí)寬脈沖信號(hào)下能輸出大于 1.5 A 線性電流的光電倍增管[J]. 核技術(shù), 2004, 27(4): 313-316

HU Mengchun, YE Wenying, ZHOU Dianzhong, et al. A kind of photomultiplier tube with the highest linear current larger than 1.5 A for measuring microsecond pulse signals[J]. Nuclear Techniques, 2004, 27(4): 313-316

3 王寶林, 劉鳳玲. 塑料閃爍體在中子輻照下發(fā)光特性的研究[J]. 核技術(shù), 1991, 14(9): 535-540

WANG Baolin, LIU Fengling. Plastic scintillators irradiation pulse emissivity brightness temperature[J]. Nuclear Techniques, 1991, 14(9): 535-540

4 彭曉世, 王峰, 徐濤, 等. 快響應(yīng)塑料閃爍體EJ232的時(shí)間特性研究[J]. 光學(xué)學(xué)報(bào), 2011, 12(31): 1216001

PENG Xiaoshi, WANG Feng, XU Tao, et al. Study on time characteristics of fast time response plastics scintillator EJ232[J]. Acta Optica Sinica, 2011, 12(31): 1216001

5 郭存, 徐榮昆, 李正宏, 等. 一種新型快塑料閃爍體的性能研究[J]. 物理學(xué)報(bào), 2004, 53(5): 1331-1334

GUO Cun, XU Rongkun, LI Zhenghong, et al. The capacity study of fast plastics scintillator for X-ray energy[J]. Acta Physica Sinica, 2004, 53(5): 1331-1334

6 Bollinger L M, Thomas G E. Measurement of the time dependence of scintillation intensity by a delayedcoincidence method[J]. Review of Scientific Instruments, 1961, 32(9): 1044-1050

7 Moszyński M, Bengtson M B. Status of timing with plastic scintillation detectors[J]. Nuclear Instruments and Methods, 1979, 158: 1-31

8 Zhong W L, Li Z H, Yang C G, et al. Measurement of decay time of liquid scintillator[J]. Nuclear Instruments and Methods in Physics Research, 2008, A587(2): 300-303

9 Hatarik R, Bernstein L A, Caggiano J A, et al. Characterizing time decay of bibenzyl scintillator using time correlated single photon countinga[J]. Review of Scientific Instruments, 2012, 83(10): 10D911

10 汲長松, 賀宣慶, 王婷婷, 等. ST401型塑料閃爍體閃爍時(shí)間測(cè)量[J]. 核電子學(xué)與探測(cè)技術(shù), 2013, 33(4): 462-466

JI Changsong, HE Xuanqing, WANG Tingting, et al. The decay time measurement of ST401 type plastic scintillator[J]. Nuclear Electronics & Detection Technology, 2013, 33(4): 462-466

11 林久令, 董志山, 張秀峰, 等. 熒光衰減曲線的一種多指數(shù)擬合方法[J]. 發(fā)光學(xué)報(bào), 2002, 23(1): 93-96

LIN Jiuling, DONG Zhishan, ZHANG Xiufeng, et al. A method for multi-exponential fitting to the fluorescence decay[J]. Chinese Journal of Luminescence, 2002, 23(1): 93-96

12 EJ-232 plastic scintillator[EB/OL]. http://www. southernscientific.co.uk, 2015

Decay time characteristics of plastic scintillator

HUANG Zhanchang YANG Jianlun
(Institute of Nuclear Physics and Chemistry, Chinese Academy of Engineering Physics, Mianyang 621900, China)

Background: The plastic scintillator is mostly used for the inertia confinement fusion determination systems. It is necessary to know the scintillator decay time. Purpose: The aim is to measure the decay time of the plastic scintillator and analyze the applicability of single photon method. Method: The plastic scintillator is excited by137Cs, and the decay time of the plastic scintillator is measured by single photon method. The ST401 and EJ232 plastic scintillators are measured by this device. The theory of single photon method is analyzed via statistics and the measurement process is simulated by Monte Carlo method. Results: The results are 2.9 ns with 50 ns for ST401 and 1.6 ns with 30 ns for EJ232 via double exponential fitting. The device has a dynamic range of over 104and an intrinsic time resolution under 0.8 ns. The experiment result agrees to the simulated one. Conclusion: The comparison between the experiment result and the simulation result indicates that with the mean counts of stop detector no more than 0.1, the result measured by single photon method can be approximated to the decay time of plastic scintillator.

Plastic scintillator, Flashing time characteristics, Single photon counting method, Monte Carlo method CLC TL812

TL812

10.11889/j.0253-3219.2015.hjs.38.060203

中國工程物理研究院科學(xué)技術(shù)發(fā)展基金(No.2013B0102008)資助

黃展常，男，1989年出生，2012年畢業(yè)于清華大學(xué)，現(xiàn)為碩士研究生，主要從事慣性約束聚變物理診斷等方面的研究

2015-01-28，

2015-03-29