基于復合結構的氣體電子倍?增器增益模擬和實驗研究

張余煉1)2)3) 祁輝榮2)3)? 胡碧濤1)2)? 溫志文1)2)3) 王海云2)3)4)歐陽群2)3) 陳元柏2)3) 張建2)3)

1)(蘭州大學核科學與技術學院,蘭州 730000)

2)(核探測與核電子學國家重點實驗室,北京 100049)

3)(中國科學院高能物理研究所,北京 100049)

4)(中國科學院大學,北京 100049)

基于復合結構的氣體電子倍?增器增益模擬和實驗研究

張余煉1)2)3) 祁輝榮2)3)? 胡碧濤1)2)? 溫志文1)2)3) 王海云2)3)4)歐陽群2)3) 陳元柏2)3) 張建2)3)

1)(蘭州大學核科學與技術學院,蘭州 730000)

2)(核探測與核電子學國家重點實驗室,北京 100049)

3)(中國科學院高能物理研究所,北京 100049)

4)(中國科學院大學,北京 100049)

(2017年3月16日收到;2017年4月24日收到修改稿)

氣體電子倍增器(GEM)作為高性能的微結構氣體探測器在高能物理相關領域內得到了廣泛的研究和應用.其中增益是GEM探測器基本性能研究中的一個重要參數,該值的精確測量至關重要.增益的測量一般采用電流測量或者能譜測量方法,但均存在精度較低或者過程繁瑣的問題,且無法精確測量低增益值.針對GEM探測器增益的精確測量,本文提出了一種由GEM探測器與微網結構氣體探測器(MM)級聯構成的復合結構探測器(GEM-MM).利用GEM-MM結構以相對方法實現GEM增益的精確測量.該方法既可以省去傳統方法中復雜的電子學標定過程,同時不需要進行原初電離電子數的估算,保證了增益的精確測量,并且可以實現GEM低增益的測量.基于GEM-MM測量GEM增益的原理,本文首先對GEM-MM電荷輸運過程進行了模擬研究,優化了合適的工作電壓.比較了三種不同類型和配比工作氣體下GEM增益模擬結果,并在Ar/iC4H10(95/5)氣體中測量了單層GEM在3—24范圍內的有效增益.不同Penning系數下GEM增益的模擬結果表明,Penning系數為0.32時GEM增益的模擬結果與實驗測量結果符合得很好.由此可以確定一個大氣壓下的Ar/iC4H10(95/5)氣體中,Penning系數為0.32±0.01.

氣體電子倍增器,微網結構氣體探測器,增益

1 引 言

氣體電子倍增器(gas electron multiplier, GEM)[1]是一種高精度分辨的微結構氣體探測器,由Sauli于1997年提出.此后,國內外對GEM開展了大量研究,促進了其在不同領域的應用[2].由于它所具有的高位置分辨率(＜70μm)[3]、優秀的能量分辨(10%)[4,5]、可以在高計數率下長時間穩定工作[6]以及易于制作安裝的特性.GEM技術正在研究用于高能物理實驗中,例如直線對撞機時間投影室(LCTPC)、緊湊型μ介子螺線管探測器(CMS)以及大型離子對撞機實驗等[7-9].GEM膜是雙面覆銅的聚酰亞胺(Kapton)膜,膜上刻蝕出貫穿的微孔結構(孔徑70μm,孔間距140μm).在GEM膜兩側銅電極上加上合適的電壓差,微孔內會形成高強的電場(30—50 kV/cm).一般以單層或者多層級聯的GEM膜作為倍增電極,在膜的兩側設置漂移極和讀出極,密閉在工作氣體(一般為惰性以氣體為主的混合氣體)環境中構成探測器.單層GEM探測器結構如圖1(a)所示.漂移極與GEM膜之間為漂移區,GEM膜與讀出極之間為感應區.保持漂移區以及感應區電場方向與孔內電場方向相同.帶電粒子在漂移區與氣體介質電離產生的原初電離電子,在電場作用下向GEM微孔內運動,在孔內的強場中發生雪崩倍增.倍增后大量的電子離開GEM膜到達下一級的倍增電極或者被收集到讀出電極,最終由讀出電子學讀出.為保證GEM的穩定性,單層GEM膜的增益通常在幾十到幾百左右,在實際應用中需要級聯多層的GEM膜以獲取高的增益.微網結構氣體探測器(micro-mesh gaseous structure,MM)[10]是另外一種被廣泛研究的微結構氣體探測器,由Saclay實驗室于1996年提出.MM探測器的結構較GEM更為簡單,包括漂移極、讀出極和金屬絲網,如圖1(b)所示.金屬絲網與讀出極之間的區域為雪崩區, MM工作時雪崩區電場強度約為30 kV/cm.原初電離電子產生于漂移極與金屬絲網之間的漂移區,在電場作用下穿過金屬絲網到達雪崩區.為保證雪崩區距離的均勻,金屬微網與讀出PCB之間為支撐柱結構.倍增后的電子被讀出極收集讀出.

圖1 探測器結構示意圖 (a)MM;(b)單層GEMFig.1.Diagram of the Micromegas and GEM detectors：(a)Micromegas detector;(b)single GEM detector.

氣體探測器的增益即氣體的放大倍數是探測器的一個重要參數.其中包括以各種方式(氣體對電子的吸附以及GEM膜對電子的吸收等)損失從而未能被讀出電極收集的電子.以GEM探測器為例,根據氣體放大機制理論[11],在一定范圍內, GEM的增益隨極電壓差的增大呈指數上升.通過增益與電壓關系的測量可以確定優化的工作電壓范圍.傳統方法通過對讀出電極上平均電流的測量確定GEM的增益[12].例如在X射線入射時,產生的原初電離電子數、X射線計數率以及增益之間的乘積與讀出電極上的電流大小成正比.該方案中需要精確測量電流大小,而通常電流是很小的.例如在55Fe X射線源照射下,電流通常約為nA量級.因此實現高精度的電流測量比較困難.實際測量中,此方法適用于在高通量的X光照下(例如X光機)的增益測量.還可以采用前置放大器加主放大器的方案實現電荷測量,但需要進行放大器的標定,確定其輸出電壓與輸入電荷的關系.這兩種方案都需要計算原初電離電子數.氣體產生一對電子離子對所需的平均能量為氣體的平均電離能(如氬氣中約為26 eV),根據帶電粒子在氣體中的能量損失由實驗測量得到.混合氣體中的平均電離能為按氣體體積百分比計算的平均值,因此根據平均電離能計算原初電離電子數會帶來較大誤差(氬氣中為10%)[13].增益為電子經過GEM放大前后的電子數的比值,無論以何種方式測量增益,估算原初電離電子數帶來的誤差都是不可避免的問題.此外,由于測量電子學對信噪比的要求,GEM需要有較高的增益.這就導致在低增益下GEM增益的測量十分困難.

另外一種GEM增益的測量方案是利用GEM與其他電子放大結構級聯實現的相對測量方法.對GEM與微條氣體室(MSGC)級聯[14],及GEM與平面微間隙探測器(PMGC)級聯[15]的研究表明,級聯后探測器總的增益顯著低于GEM增益與MSGC或者PMGC增益的乘積.這是因為MSGC以及PMGC的陰極和陽極為同一平面上的左右結構,其陰極和陽極上不同的電壓導致GEM與它們之間傳輸區電場不均勻,整體場強降低.導致GEM膜的電子出孔率降低,從而使級聯探測器的增益顯著下降.此外MSGC與PMGC的微電極附近電子倍增累積的空間電荷效應也是增益下降的一個重要原因.而且MSGC以及PMGC極易受到打火影響而損傷.為避免傳統測量方法和相對測量方法的問題,本文介紹了一種GEM結合MM的復合結構(GEM-MM)測量GEM增益的方案.這種結構利用了GEM膜電子倍增放大與讀出電極相互獨立這一特性,將GEM膜級聯于MM的金屬微網與漂移電極之間(參見圖1).因為位于讀出極與GEM膜之間的MM微網為平面結構,且微網與讀出極為上下結構.只需在微網上加合適的電壓就可以實現對其兩側區域電場的控制.改變MM工作電壓時,MM與GEM之間的傳輸區電場不受影響.通過優化GEM以及MM的工作電壓,在GEM增益很低時(＜10)也可以容易地區分出MM的增益以及GEM的增益.而在MM中電子倍增發生在雪崩區的均勻空間內,空間電荷效應影響很小[16].因此GEM-MM結構較之前提出的GEM與其他探測器級聯的結構有明顯優勢.尤其是在GEM增益較低情況下,通過增加MM的增益,采用相對方法能夠有效地測量GEM增益并保證一定的測量精度. GEM-MM結構還可以降低打火率的同時提高MM的增益[17],研究表明這種復合結構有效地提高了整體探測器的增益,同等增益下顯著降低了打火率,且沒有出現GEM-MSGC以及GEM-PMGC結構中增益降低的現象.

下文中首先介紹利用有限元分析軟件以及Gar field++程序對GEM-MM探測器的模擬研究結果;通過比較探測器中各部分區域場強及工作電壓對電子在各個電極上的收集比例的影響,優化選擇了GEM-MM合適的工作電壓;同時模擬了GEM探測器在不同氣體中的增益,并進行了比較;實驗部分介紹了該方法測量GEM增益的原理以及實驗測量裝置,隨后給出了GEM增益的實驗測量結果以及誤差分析,最后對GEM增益測量的實驗結果與模擬結果進行了對比,并確定了Ar/iC4H10(95/5)氣體的Penning系數.

2 模擬研究

2.1 建模及模擬方法

模擬采用的主要工具是Gar field++程序[18],它是對氣體探測器精確模擬的專用程序.利用Gar field++可以計算并繪制靜電場、電子和離子在氣體中的運動以及計算讀出電極上的感應信號.氣體中電荷的運動依賴于探測器的電場配置,因此電場計算是Gar field++模擬研究的基礎.但是Gar field++本身只能處理簡單的二維電場的計算問題,針對微結構的GEM電場計算需要采用專用的電場計算軟件.模擬中采用了開源的有限元分析和計算程序Gmsh[19]和Elmer[20].Gmsh可以進行探測器的幾何建模和網格劃分.GEM-MM結構的模型建立如圖2(a)所示.漂移區距離為4 mm, GEM膜與金屬網之間的傳輸區距離為1.4 mm.圖2(b)所示為GEM探測器的模型,模擬中漂移區間距及感應區間距均為1 mm.Elmer通過讀取Gmsh的網格劃分進行電場的計算.由于模擬中探測器建模為選取的可重復小單元,電場計算時,邊界條件的設置十分重要.合適的邊界條件設定可以確保對電場計算結果在x和y方向重復時的準確性.在Gmsh對GEM-MM建模中,設置各個電極的表面為不同電壓的邊界條件.長方體單元的相對側面設置為對稱邊界條件.

將電場計算的結果導入到Gar field++中,由其所提供的Magboltz接口定義探測器的工作氣體.由于工作氣體參數,如漂移速度、擴散系數、吸附系數以及倍增系數等對電子輸運過程有重要的影響,利用Magboltz程序計算了氣體的相關性質.由于這些氣體參數依賴于電場大小,最后由AvalancheMicroscopic類提供的方法計算電子在氣體中的輸運過程時根據電子所處位置的電場大小考慮了氣體參數的影響.Gar field++提供了ROOT接口,計算的結果保存為ROOT文件,進行相應的數據處理和分析.

模擬研究主要包括兩方面的內容.首先是對電子在GEM-MM探測器中電荷輸運過程的模擬,目的在于確定電子在GEM-MM中各個電極上的收集情況及其與GEM-MM探測器電場配置的關系.模擬中,電子隨機放置在圖2(a)所示的漂移區內,距離GEM膜上表面為0.5 mm.計算電子在不同的探測器電場配置下的漂移以及放大過程,統計5000次電子事例及其對應的所有倍增電子在探測器中的位置信息并保存為ROOT文件.對GEM增益的計算是模擬研究的第二個重要內容.在圖2(b)所示模型的漂移區內距離GEM膜上表面為0.5 mm的平面上隨機放置原初電子,在不同的GEM電壓下模擬電子在經過GEM膜后的倍增數目.統計5000個原初電子分別對應的電子倍增數目分布,分析得到GEM的增益.模擬中工作氣體設置為Ar/iC4H10(95/5)混合氣體,溫度為20°C,氣壓為1 atm.計算了GEM在三種不同類型和配比氣體下的增益.

圖2 (網刊彩色)使用Gmsh建立的探測器幾何構型 (a)GEM-MM結構,自上至下為漂移極、GEM膜、金屬絲網和讀出陽極;(b)GEM探測器結構,包括漂移極、GEM膜和讀出陽極.圖中單位為cmFig.2.(color online)The structure of detector module by Gmsh software tool(a)GEM-MM detector：drift cathode,GEM,stainless steel mesh and readout anode from top to bottom(b)GEM detector：drift anode, GEM and readout anode(in cm).

2.2 模擬結果

2.2.1 電場計算

探測器中電子的運動與其所處的電場環境密切相關,因此對探測器的電場計算是后續模擬研究的基礎,電場計算結果如圖3所示.圖3(a)為探測器在x-z平面的電場大小分布,為了更好地展示出GEM膜以及MM金屬微網附近的電場結構,對傳輸區的距離做了調整.從圖中可以看出,僅在GEM以及MM金屬微網附近存在不均勻的電場.且這種不均勻的電場向GEM膜以及金屬網的微孔聚攏,這一特性有利于提高電子的透過率.圖3(b)為圖3(a)所示平面內沿x=0的直線上電場大小的分布,橫軸表示直線上距離最下方讀出陽極的相對距離.從圖中可見MM雪崩區、GEM微孔的高強電場以及感應區和漂移區的電場.GEM-MM結構中電場配置分別是漂移區場強(Ed),GEM工作電壓(VGEM)、傳輸區場強(Et)以及MM工作電壓(Vmesh).圖3所示的電場計算結果對應的電場配置為：Ed=1 kV/cm,VGEM=300 V, Et=1.5 kV/cm,Vmesh=440 V.

圖3(網刊彩色)由ANSYS計算的GEM-MM探測器電場強度 (a)在x-z平面的分布;(b)沿穿過GEM孔中心的直線路徑上的分布Fig.3.(color online)The electric field of GEM-MM detector calculated by ANSYS：(a)Distribution in the x-z plane;(b)distribution along the center of the GEM detector’s hole.

2.2.2 電子透過率

模擬中電子由漂移區向GEM膜運動,經過GEM的雪崩放大后繼續向MM傳輸,經過第二次倍增后被讀出電極收集.該過程中電子有可能終止在其他電極或者氣體介質內,而這一過程與探測器的電場配置相關.因此需要優化確定探測器電場配置以保證電子在讀出電極上的收集效率.模擬中根據電子終止位置,確定出其最后是被探測器的哪一個電極所收集.探測器中電子在各個電極的收集情況與GEM-MM各部分電場(或工作電壓)關系的模擬結果如圖4所示,圖中縱軸表示探測器中在各個電極上收集的電子數與總電子數(包括初始電子事例以及倍增產生的電子)的比值.圖4(a)—圖4(d)所示為分別改變四個電場變量中的Ed,VGEM,Et以及Vmesh,保持其他三個變量不變時電子在各個電極的收集情況隨電場及電壓的變化.如改變MM的工作電壓Vmesh時,GEM膜上下表面以及漂移極的電壓均隨之改變,但保持傳輸區、漂移區兩端以及GEM膜兩側電極壓差不變,從而保持VGEM,Et以及Vmesh不變.四部分場強及電壓分別設置為Ed=100 V/cm,VGEM=340 V,Et=100 V/cm及Vmesh=420 V.圖中未給出電極上電子收集比例為零或者十分接近于零.

漂移區電場對電子輸運影響較小,改變不同漂移場強,電子在各電極收集的比例基本保持不變,如圖4(a)所示.GEM與MM的工作電壓影響則非常明顯,如圖4(b)及圖4(d)所示.當GEM工作電壓升高保持傳輸區場強不變時,GEM膜內倍增電子的出孔率會相應降低,導致更多的電子被收集于GEM膜的下表面電極.在GEM-MM探測器中,為保持GEM膜對電子的出孔率不變,需要保持GEM的工作電壓與傳輸區兩側的電壓差的比例不變.即在改變GEM工作電壓時,傳輸區電場也需要隨之改變.MM工作電壓越高,金屬網對電子的透過率越高.由此可知,要提高電子的收集效率,MM需要工作在較高的工作電壓下,而GEM工作電壓不宜過高.圖4(c)描述了傳輸區電場對電子收集比例的影響.在GEM工作電壓為340 V時,傳輸區場強高于500 V/cm后,GEM的電子出孔率基本達到最大并保持恒定.

2.2.3 增 益

對GEM增益的模擬是研究的一個主要內容.由于Penning效應對增益有著一定的貢獻,在Gar field++的模擬中以系數r描述了Penning效應帶來的氣體電離.r表示具有激發能高于氣體電離能的激發態所引發的氣體電離概率.在設定的Ar/iC4H10(95/5)混合氣體中,r參考值選取為0.40[21].模擬中漂移區和感應區的場強分別是Ed=250 V/cm和Ei=1 kV/cm.GEM電子倍增對應的增益可以用Polya分布描述[22]

式中C0為常數;G0為單電子倍增分布的平均增益; θ是決定Polya分布方差的參數,其大小與能量超過電離閾值的電子比例有關.圖5所示為GEM電壓為300 V時,5000個電子經過GEM膜倍增后電子數量分布的模擬結果.紅色線為按照(1)式擬合的結果,可以看出Polya分布可以很好地描述GEM的增益.

圖5 模擬計算Ar/iC4H10(95/5)混合氣體中GEM單電子倍增數分布,VGEM=300 V,Ed=250 V/cm, Ei=1 kV/cm,r=0.40Fig.5.Simulated single electron avalanche size distribution of GEM in Ar/iC4H10(95/5)with VGEM= 300 V,Ed=250 V/cm,Ei=1 kV/cm,Penning transfer rate r=0.40.

但是當GEM電壓較低時,模擬中發現增益并不能用Polya很好地描述.圖6是GEM電壓為150 V時,5000個電子經過GEM膜倍增后總的電子數量分布模擬結果.此時倍增電子數分布沒有Polya分布中明顯的最大值,這時用指數分布描述更為合適[23],

式中C1為常數.圖6所示的擬合結果中Slope= -1/G0,其中G0為單電子倍增后的平均增益.

圖6 模擬計算Ar/iC4H10(95/5)混合氣體中GEM單電子倍增數分布,VGEM=150 V,Ed=250 V/cm, Ei=1 kV/cm,r=0.40Fig.6. Simulated single electron avalanche size distribution ofGEM in the operation gas of Ar/iC4H10(95/5)at VGEM=150 V,Ed=250 V/cm, Ei=1 kV/cm,r=0.40.

模擬中對GEM的總增益和有效增益做了區分.總增益(G)即每個電子發生雪崩后的總電子數,其中包括以各種方式(氣體對電子的吸附以及GEM膜結構對電子的吸收等)損失從而未能被讀出電極收集的電子.有效增益(Geff)則指原初電子經GEM膜倍增后到達讀出電極的電子數,不包括在到達讀出電極之前損失掉的電子.圖5及圖6所示為總增益,本文中實驗部分增益均為有效增益,模擬中提到的增益除特別指出外,也指有效增益.模擬的GEM單電子增益分布在不同的GEM電壓下采取了不同的擬合方式(指數擬合以及Polya擬合).從單電子倍增數分布中直接獲取平均值也可以得到增益.Ar/iC4H10(95/5)氣體中GEM增益的模擬結果列于表1.從結果來看,總增益與有效增益區別較大,這主要是由GEM膜(包括兩側金屬電極以及中間Kapton膜)對電子的吸收造成的.模擬給出了Ar/iC4H10(95/5), Ar/CO2(90/10)(r=0.452)以及Ar/CO2(70/30) (r=0.57)[24]混合氣體中的增益特性,結果如圖7所示.從圖中可以看出,在GEM相同工作電壓下,Ar/iC4H10(95/5)較其他兩種氣體具有最高增益,且與GEM工作電壓之間表現出很好的指數相關性.GEM增益的實驗測量中,選取了Ar/iC4H10(95/5)工作氣體.

表1 Ar/iC4H10(95/5)(r=0.40)混合氣體中GEM增益模擬結果Table 1.Total gain and e ff ective gain of GEM simulation in the Ar/iC4H10(95/5)with r=0.40.

圖7 (網刊彩色)不同氣體中GEM增益與GEM工作電壓關系的模擬結果,Ed=250 V/cm,Ei=1 kV/cmFig.7.(color online)Simulated gain if GEM in the different gas mixtures at Ed=250 V/cm,Ei=1kV/cm.

3 測量原理與實驗裝置

3.1GEM-MM結構和增益測量原理

GEM-MM的結構如圖8所示,自上而下為漂移極、GEM膜、金屬絲網和讀出陽極.GEM膜為歐洲核子中心提供的標準結構GEM.MM為與Saclay實驗室合作的基于Bulk工藝的探測器,其雪崩區間距為128μm.金屬絲網絲徑22μm,絲間距40μm.漂移區和傳輸區間距分別為4和1.4 mm.

圖8GEM-MM測器結構示意圖Fig.8.The schematic diagram of GEM-MM detector.

實驗采用的工作氣體為氬基混合氣體(Ar/iC4H10(95/5)). 光電效應是能量EX為5.9 keV的55Fe X射線與氬原子相互作用的主要物理過程.發生光電效應時,光子與氬原子內層電子或束縛電子發生相互作用,其能量被完全吸收,用于電離或者激發氬原子.電離過程中,獲得能量的電子發射出來成為自由電子(光電子).發射特征X射線能量為

式中EK和EL分別為K和L殼層電子的結合能.另外一種方式為發射俄歇電子,這兩種過程同時存在.光電效應過程中產生的所有電子(包括光電子和俄歇電子)由于能量較低,其能量全部沉積在探測器的氣體介質中.因此在實驗測量的能譜中存在一個與俄歇電子發射相對應的X光子能量全部沉積形成的全能峰,以及特征X射線過程相對應的逃逸峰.

X射線由探測器上方的入射窗進入工作氣體內.5.9 keV X射線在厚度為4 mm的氬氣中的轉換效率約13%[25],部分光子在漂移區沉積能量,還有部分光子在傳輸區與氣體發生相互作用.因此在漂移區和傳輸區都會產生光電子.漂移區產生的光電子在電場的作用下經過GEM膜和MM兩級的倍增放大,最終在讀出極收集為Nt1.而在傳輸區產生的光電子只經過MM一級的放大,最后被收集在讀出極上為Nt2.因此,實驗測量得到的55Fe能譜包含有兩組全能峰和逃逸峰.設GEM膜和MM的增益分別為Geff和Gmm,GEM-MM的增益為Ggem-mm,則有

由此,GEM的增益可表示為

而讀出極收集的電荷量正比于多道的道數,因此

式中CH1和CH2分別為漂移區和傳輸區光子能量沉積對應的多道的道數,GEM膜的增益測量就轉換為多道道數比值的測量,無需電子學的標定以及對原初電離電子數的估算.

3.2 實驗裝置

增益測量實驗裝置如圖9所示.整體的GEMMM探測器有效面積為50 mm×50 mm,讀出電極為面積為50 mm×50 mm的單個Pad.MM雪崩區的倍增電子在金屬網和讀出陽極感應出大小相等極性相反的信號,因而從兩個電極都可以進行測量.但是實驗測量中發現,若從陽極讀出,由于金屬絲網電極高壓端需要加保護電阻,而保護電阻會使得陽極讀出的信號遠小于不接保護電阻時的信號.此外由于所使用的電荷靈敏前置放大器輸入端存在隔直電容,因此陽極讀出時還需要加接地電阻,以利于電荷泄放,避免探測器工作不穩定.鑒于此,在只有單個讀出單元的GEM-MM探測器測試中,采取了Pad接地,從金屬絲網讀出信號的方式.其余電極采用了CAEN N471型高壓插件提供高壓,信號讀出前端采用ORTEC 142IH電荷靈敏型前置放大器,主放為ORTEC 572 A,成形時間設置為1μs.數據最后由ORTEC ASPEC927多道分析器收集.射線源活度為5mCi的55Fe X放射源.采用流氣式工作模式,氣流量為30mL/min.

圖9 復合結構探測器實物圖Fig.9.Photos of the GEM-MM detector in lab.

4 測量結果分析

4.1 能譜測量

讀出Pad上的感應電荷,經過前置放大器和主放大器放大成形后,由多道分析器采集.圖10為GEM-MM探測器的能譜測量結果.GEM和MM的工作電壓分別為240和370 V,隨漂移區和傳輸區場強分別為300和625 V/cm.從圖中可以明顯地看到四個峰位,前兩個分別為傳輸區光子事例的逃逸峰和全能峰,后兩個則分別對應漂移區光子事例的逃逸峰和全能峰.圖中同時給出了高斯擬合的結果,擬合后可確定出全能峰對應的多道的道值及其誤差.利用傳統的增益測量方法,可以確定出MM和GEM-MM探測器的增益.GEM的增益根據(7)式由能譜直接測量得到.即分別擬合出兩個全能峰的道值,后面道數與前面道數的比值就是GEM的增益.

圖10(網刊彩色)GEM-MM55Fe放射源X射線能譜測量圖 VGEM=220 V,VMM=370 V,Et=625 V/cm, Ed=300 V/cmFig.10. (color online)Energy spectrum of GEMMM detector with55Fe source@VGEM=220 V, VMM=370 V,Et=625 V/cm,Ed=300 V/cm.

采用相對方法測量GEM增益誤差來源于對兩個全能峰峰位的測量.如圖10所示,對GEM-MM對應的全能峰和MM對應的全能峰分別做高斯擬合后,平均值為所得峰位取高斯擬合的標準差為峰位測量的誤差(分別為σCH1.在圖10所示的能譜測量中,利用擬合結果可得GEMMM相對方法測量GEM增益的相對誤差為

另一方面,傳統增益測量方法中,由原處電離電子計算帶來的誤差主要來源于平均電離能的誤差σw、刻度電子學σe的誤差以及能譜測量的誤差σch.總的誤差可表示為

由文獻給出的平均電離能誤差,在混合氣體中,平均電離能的計算誤差可由單一氣體的計算誤差得到,在Ar-iC4H10混合氣體中為14%;σe小于1%, σch即為(8)式中σCH1.由此可見,采用相對方法測量更加精確.

4.2 增益測量

GEM-MM探測器的增益與GEM和MM的工作電壓直接相關.此外漂移區和傳輸區的場強同樣會對增益有影響.實驗中首先測量了增益隨漂移區和傳輸區場強的變化關系,結果如圖11所示.圖中增益為GEM-MM的增益,通過傳統方法測量得到.

優化測量中,MM和GEM的工作電壓分別保持在VMM=370 V和VGEM=240 V.首先傳輸區的場強固定在Et=625 V/cm,測量了增益隨Ed的變化,如圖11(a)所示.從圖中可以看出,隨著漂移區場強的增加,增益先增大后減小.這是因為在漂移區場強較低時,電離電子不能完全進入到GEM中得到放大,有相當一部分的電子在漂移區被吸收.隨著電場強度的增加,電離電子到達GEM膜的比例達到最大值.但是當漂移區電場強度進一步增大時,電子被收集到GEM膜上表面電極的比例增加,從而導致增益下降.優化選取Ed=300 V/cm.隨后Ed保持300 V/cm不變,測量增益隨傳輸區場強變化關系,如圖11(b)所示.隨著Et的增加,增益呈現出先增大后減小的變化.原因是MM金屬絲網對電子的透過率隨著傳輸區場強的增加是先升高后降低的.優化選取最大增益處Et=625 V/cm,據此確定GEM電壓與漂移區兩側電壓差的比值以及傳輸區兩側壓差分別與VGEM及VMM的比值.優化工作電壓后,在后續增益測量中,改變GEM或MM工作電壓時,保持Et和Ed同比例變化,以保證GEM-MM最大的電子透過率.

圖11 GEM-MM測量增益 (a)隨漂移區場強變化,VGEM=240 V,VMM=370 V,Et=625 V/cm;(b)隨傳輸區場強變化,VGEM=240 V,VMM=370 V,Ed=300 V/cmFig.11.Measurement of GEM-MM detector’s gain：(a)Ed10–1500 V/cm@VGEM=240 V,VMM=370 V, Et=625 V/cm;(b)Et100–3000 V/cm@VGEM=240 V,VMM=370 V,Ed=300 V/cm.

圖12 探測器測量增益 (a)隨GEM工作電壓的變化,VMM=370 V;(b)隨MM工作電壓的變化,VGEM= 240 V;Ed及Et隨GEM或MM工作電壓按比例變化Fig.12.Measurement of the detector’s gain：(a)VGEM190–270 V@VMM=370 V;(b)VMM300–390 V@ VGEM=240 V and Ed/Etaccording to the optimized voltage ratio.

根據能譜測量結果,進行了GEM-MM增益與GEM和MM工作電壓關系的測量,結果圖12所示.其中GEM增益根據能譜中全能峰道值直接得到, GEM-MM探測器的增益和MM的增益采用傳統方法經過電子學的標定測量得到.

5 測量結果與模擬結果對比

利用Gar field++計算的Ar/iC4H10(95/5)氣體中GEM增益與利用GEM-MM復合結構采用相對方法測量的GEM增益進行了對比,結果如圖13所示.圖中Geff( fi tted mean)為采用擬合方法獲取的增益,Geff(mean)表示增益分布直方圖的平均值.

圖13 (網刊彩色)GEM增益隨工作電壓變化的模擬與實驗測量結果比較Fig.13.(color online)Comparison of GEM gain simulation and measurement results.

模擬的有效增益結果與GEM電壓的關系如圖13所示,可以看出增益越低,模擬中采用取平均值和擬合得到的增益二者結果差別越大.尤其是在增益較低時(小于5),增益與GEM電壓的關系偏離指數關系.這是因為在GEM工作電壓較低時,有很大一部分的電子經過GEM膜后沒有發生倍增,且會大概率地被GEM膜所收集.從而導致了有效增益小于1,擬合方法得到的結果會有較大偏差.因此在增益的模擬中,在GEM增益較低時取單電子增益分布的平均值作為平均增益更為合適.圖中紅色三角符號標識為GEM增益實驗測量結果,藍色三角符號標識為r=0.40時增益的模擬結果,可以看到模擬增益較測量結果偏高.這是因為文獻[21]中只給出了Ar/iC4H10(90/10)氣體中Penning系數的取值(r=0.40),而實驗中工作氣體為Ar/iC4H10(95/5).于是進行了不同r取值的增益計算,圖中黑色標識的實心圓點為r=0.32時的增益模擬結果,可以看到模擬的增益結果與實驗測量結果基本一致.通過對Penning系數取值的進一步微調,當r=0.32±0.01時模擬的增益值在實驗測量1σ的誤差范圍之內.由此得到Ar/iC4H10(95/5)混合氣體中(一個大氣壓),Penning系數取值為r=0.32±0.01.

6 結 論

本文提出了一種由GEM與MM級聯構成的復合結構探測器GEM-MM.GEM與MM級聯的方案較之前的與PMGC以及MSGC級聯方案有明顯的優勢.在該結構探測器的能譜測量中可以清楚地區分GEM和MM對整體增益的貢獻.在55Fe X射線入射時,測量的能譜中出現兩個全能峰及逃逸峰,兩個全能峰的比值即為GEM的有效增益.適當增加MM的增益后,以這種相對方法可以實現GEM低增益的測量.此外這種方法測量GEM的增益避免了對電子學繁雜的標定過程,且無需進行原初電離電子數的估算.模擬研究中確定了其合適的工作電壓以及工作氣體.實驗中研制了GEMMM探測器,并在Ar/iC4H10(95/5)混合氣體中進行了GEM增益的相對方法測量.得到了GEM在3—24范圍內的有效增益.實際上,單層GEM增益為20時,三層GEM級聯探測器的總增益為6000.這正是三層GEM級聯探測器的通常工作增益,本文中測量的單層GEM增益范圍包含了多層GEM級聯探測器的通常工作增益范圍.高的探測器增益是為了在電子學噪聲水平較高時提高信噪比,而對于低噪聲水平的電子學,較低的探測器增益就能滿足信噪比的要求.因此低GEM增益的測量對于低噪聲水平的電子學設計具有重要的參考意義.在GEM和MM的工作電壓分別為220和370 V時,這種方法測量GEM增益的相對誤差為14.96%.這一誤差小于傳統增益測量方法的測量誤差,因而更加精確.結合實驗測量結果與模擬研究,通過調整Penning系數獲得與實驗測量GEM增益相符的結果,確定了Ar/iC4H10(95/5)氣體中的Penning系數為0.32±0.01.基于本文的研究工作,下一步會開展對Ar-CF4-iC4H10(95：3：2)(所謂的T2K氣體)以及其他工作氣體Penning系數以及多層級聯GEM增益精細測量的研究.

感謝清華大學工物系高原寧教授對復合結構氣體探測器研制中的支持以及對于本文寫作的中肯建議;感謝法國原子能研究院CEA-Saclay氣體探測器組負責人Roy教授對微網結構氣體探測器研制和研究工作的幫助及支持;感謝中國科學院高能物理研究所苑長征研究員和唐光毅博士對本文寫作的建議.

[1]Sauli F 1997 Nucl.Instrum.Methods A 386 531

[2]Sauli F 2016 Nucl.Instrum.Methods A 805 2

[3]Ketzer B,Weitzel Q,Paul S,Sauli F,Ropelewski L 2004 Nucl.Instrum.Methods A 535 314

[4]Bressan A,de Oliveira R,Gandi A,Labbé J C 1999 Nucl.Instrum.Methods A 425 254

[5]Ketzer B 2013 Nucl.Instrum.Methods A 732 237

[6]Benlloch J M,Dokoutchaeva V,Malakhov N,Menzione A,Munar A 1998 Nucl.Instrum.Methods A 419 410

[7]Tsionou D 2017 Nucl.Instrum.Methods A 845 309

[8]Abbaneo D,Abbas M,Abbrescia M,et al.2017 Nucl. Instrum.Methods A 845 298

[9]Lippmann C 2016 Nucl.Instrum.Methods A 824 543

[10]Giomataris Y,Rebourgeard Ph,Robert J,Charpak G 1996 Nucl.Instrum.Methods A 376 29

[11]Blum W,Rolandi L 1993 Particle Detection with Drift Chambers(Berlin：Springer)p125

[12]Sn?ll J 2016 M.S.Thesis(Lund：Lund University)

[13]Binks W 1954 Acta Radiologica 41 85

[14]Benlloch J,Bressan A,Buttner C,Capeans M,Gruwe M,Hoch M,Labbe J C,Placci A,Ropelewski L,Sauli F,Sharma A,Veenhof R 1998 IEEE Trans.Nucl.Sci. 45 234

[15]Bellazzini R,Brez A,Gariano G,Latronico L,Lumb N, Spandre G,Massai M M,Ra ff o R,Spezziga M A 1998 Nucl.Instrum.Methods A 419 429

[16]Charpak G,Derré J,Giomataris Y,Rebourgeard P 2002 Nucl.Instrum.Methods A 478 26

[17]Kane S,May J,Miyamoto J,Shipsey I 2003 Nucl.Instrum.Methods A 515 261

[18]Schindler H,Veenhof R Gar field++,https：//gar fieldppwebcernch/gar fieldpp/[2016-10-04]

[19]Geuzaine C,Remacle J F 2009 Int.J.Numer.Methods Engineer.79 1309

[20]CSC-IT Center for Science LTD,Elmer,https：//wwwcsc fi/web/elmer/elmer[2016-10-04]

[21]Sahin ?,Tapan I,?zmutlu E N,Veenhof R 2010 JINST 5 05002

[22]Zerguerras T,Genolini B,Lepeltier V,Peyré J,Pouthas J,Rosier P 2009 Nucl.Instrum.Methods A 608 397

[23]Mir J A,Maia J M,Concei??o A S,et al.2008 IEEE Trans.Nucl.Sci.55 2334

[24]Sahin ?,Kowalski T Z,Veenhof R 2014 Nucl.Instrum. Methods A 768 104

[25]Xie Y G,Chen C,Wang M,Lü J G,Meng X C,Wang F,Gu S D,Guo Y N 2003 Nuclear Detector and Data Acquisition(Beijing：Science Press)p628(in Chinese) [謝一岡,陳昌,王曼,呂軍光,孟祥承,王峰,顧樹棣,過雅南2003 粒子探測器與數據獲取 (北京：科學出版社)第628頁]

PACS:29.40.Cs DOI:10.7498/aps.66.142901

Measurement and simulation of the hybrid structure gaseous detector gain?

Zhang Yu-Lian1)2)3)Qi Hui-Rong2)3)?Hu Bi-Tao1)2)?Wen Zhi-Wen1)2)3)Wang Hai-Yun2)3)4)Ouyang Qun2)3)Chen Yuan-Bo2)3)Zhang Jian2)3)
1)(School of Nuclear Science and Technology,Lanzhou University,Lanzhou 730000,China)
2)(State Key Laboratory of Particle Detection and Electronics,Beijing 100049,China)
3)(Institute of High Energy Physics,Chinese Academy of Sciences,Beijing 100049,China)
4)(University of Chinese Academy of Sciences,Beijing 100049,China)

16 March 2017;revised manuscript

24 April 2017)

As one of the most popular micro pattern gaseous detectors,gas electron multiplier(GEM)has been extensively studied and applied in recent years.The studies of the detector gain measurement and simulation are important, especially on a low gain scale.Traditionally,the gain measurement is realized by measuring the current or the pulse height spectrum.The former needs complicated electronic chain calibration and the latter needs necessarily to calculate the primary electron number.In this paper,an alternative method to determine the e ff ective gain of GEM is introduced. The GEM gain can be precisely achieved through a gaseous detector of hybrid structure which combines GEM with micromesh gaseous structure(MM).The hybrid structure is called GEM-MM for short.The GEM-MM detector consists of drift cathode,standard GEM foil,stainless steel micro mesh,and readout anode.In this detector,the space between the cathode and the GEM foil is called drift gap and the other space between the GEM foil and the mesh is named transfer gap.When the X-rays irradiate into the gas volume of GEM-MM,the primary ionization occurs in both regions. Photoelectrons in the drift gap transfer from the drift region to ampli fi cation sensitive areas of the GEM and the MM detector while those in the transfer region are only ampli fi ed by the MM detector.In the energy spectrum of55Fe,there is a clear energy pro fi le including two sets of peaks.The gain of GEM can be easily obtained from the energy spectrum. Meanwhile,detailed simulations are carried out with Gar field++software package.Simulation of the electron transport parameters has been optimized.and the gains of GEM detector are also calculated for three di ff erent gas mixtures. Experimental results of the gains ranging from 3 to 24 are obtained.The gains of GEM under di ff erent working voltages are studied precisely from the spectrum measurements.The Penning transfer rate could reach 0.32±0.01 when the simulated value matches the measurement within 1σ error.

gas electron multiplier,micromegas,gain

:29.40.Cs

10.7498/aps.66.142901

?國家重點研發計劃“大科學裝置前沿研究”重點專項(批準號：2016YFA0400400)、國家自然科學基金(批準號：11675197)和中國科學院高能物理研究所創新基金資助的課題.

?通信作者.E-mail：qihr@ihep.ac.cn

?通信作者.E-mail：hubt@lzu.edu.cn

?2017中國物理學會Chinese Physical Society

http://wulixb.iphy.ac.cn

*ProjectsupportedbytheNationalKeyProgrammeforS&TResearchandDevelopment,China(GrantNo. 2016YFA0400400),the National Natural Science Foundation of China(Grant No.11675197),and the Innovation Fund of Institute of High Energy Physics,Chinese Academy of Sciences.

?Corresponding author.E-mail：qihr@ihep.ac.cn

?Corresponding author.E-mail：hubt@lzu.edu.cn