液体アルゴン

液体アルゴンTPC検出器の開発
早稲田大学,KEKA,岩手大学B,ETHZC
岡本飛鳥,岡本迅人,長坂優志,永野間淳二,三谷貴志,寄田浩平,
笠見勝祐A,木村誠宏A,小林隆A,田中雅士A,西川公一郎A,
長谷川琢哉A,牧宗慶A,丸山和純A,吉岡正和A,
内藤祐貴B,成田晋也B,
A.BadertscherC,A.CurioniC,S.DiLuiseC,U.DegundaC,L.EpprechtC,L.EspositoC,
A.GendottiC,S.HorikawaC,L.KnechtC,C.LazzaroC,D.LussiC,A.MarchionniC,
A.MeregagliaC,G.NattererC,F.PetroloC,F.ResnatiC,A.RubbiaC,C.StrabelC,T.ViantC
2011年5月13日(金)高エネルギー物理 春の学校@彦根ビューホテル
Introduction
TPC(Time Projection Chamber)とは
• 荷電粒子による電離電子の読み出しによ
り、その飛跡を3次元で高精度に再構成で
きる検出器
• 電離信号の大きさから、液体アルゴン中で
のdE/dxを測定することも可能
Ionization
Electron
E
Cherenkov Light
液体アルゴンの特徴・性質
•
•
•
•
電離電子(~5×104個/cm MIP)
シンチレーション光(~4×104photon/cm MIP 128nm)
電子捕獲をしない(O2の除去、純度が重要)
密度 1.39g/cm3,沸点-186℃
Physics Motivation
 νμ→νe振動におけるCP位相の測定
 陽子崩壊 p→ νK+ (250LArTPC,K1.1BR実験)
 暗黒物質探索(10LArTPC,早稲田開発研究)
05/13/2011
Charged
Particles
Scintillation Light
液相
2
Outline
• J-PARC K1.1BR実験
– 研究動機
– 250LAr検出器とデータ、ビームライン
– 取得信号例、シミュレーションデータ
– 純度測定
• 10LArTPC開発研究@早稲田大学
– 研究動機
– シンチレーション光観測
– 電離電子増幅
05/13/2011
3
J-PARC K1.1BR 実験
05/13/2011
4
研究動機
• 核子崩壊
– 核子の寿命は新しい物理に対して重要なパラメーター
LAr検出器はp→K+νのK+の信号を直接検出可能(SUSY GUTで大きな分岐比)
シミュレーション
K+340MeV/c
Ar中の飛程~ 20 cm
停止直前のK+のdE/dxは大きい!
Simulation Study(JHEP 0704:041,2007)
Background reduction ~105
Kaon efficiency ~97%
K1.1BR実験では
LAr検出器のdE/dxによるK/π分離能力を実験的に検証する
05/13/2011
5
250LAr検出器とデータ
時間
ADC count
Anode
入射粒子の飛跡
75ch
0ch
Particle
Drift
e e
z
Real Data
Scintillation light
e
Ionization
e
Electric field
停止点
(崩壊点)
LAr
Cathode
z
y
x
崩壊粒子の
飛跡
PMT
x
• dE/dxによって、K/π分離能力を測定するために…
– Tracking Algorithmの確立(荷電粒子の飛跡、停止点の決定)
– 減衰による電荷補正←宇宙線μによる純度測定
– チャンネル間補正(電場補正、アンプ特性)
…それぞれシミュレーションデータによる検証が必要。
05/13/2011
6
取得信号例
05/13/2011
7
Beam Line Setup
Fitch Cherenkov
(π-ring,K-ring)
250LArTPC
LArBDC
BDC=Beam Defining Counter
…Beam’s Final Focus (FF)
Gas Cherenkov
TOF2
particle
TOF1
Degrader
K,π,p,e
空気
• K1.1BRのビーム
– 運動量:800MeV/c (+degrader)
– K/π~1/4、3個/1 beam bunch(6s)
• 粒子同定
– TOF1、TOF2
– Fitch Cherenkov(K/π)
– Gas Cherenkov
• トリガー条件
BDC & TOF1 & TOF2 & LArBDC
8
Simulation Setup (Beam Line)
~5.3m
Degrader
BDC
K Beam
TOF1
Fitch Cherenkov(FC)
物質
厚さ
[mm]
Fitch
Cherenkov
アクリル
40
BDC
プラスチック
5
TOF
プラスチック
20,25
Degrader
鉛ガラス
125
鉛ブロック
25
LArBDC
プラスチック
5
250LAr
SUS304(蓋*)
50
液体アルゴン
836.4
TOF2
250LAr
• 250LAr検出器とBeam Line上のそれぞれの
検出器をGeant4を用いてSimulationに実装。
黒い部分は空気で満たした。
05/13/2011
検出器
検出器の物質と厚さ
*ただし、中心のビーム窓(210×210mm2)はハニカム構造
になっている…放射長が短い(radiation length ~ 0.16X0)
9
Simulation (K+→μ++νμ)
μ+
250LAr
K+
νμ
05/13/2011
Geant4によるシミュレーションイベント例
10
Simulation Data
Simulation
Real
TPC Ch10
TPC Ch10
シミュレーションデータでTracking Algorithm等の評価を行っている。
05/13/2011
11
純度解析
• LAr中の不純物によりドリフト電子が捕獲される
– N (t Dri )  N 0 exp  t Dri  e 
– 40cmのドリフトには1ppb以下の純度が必要!
(1ppb ⇔ e=300s)
ある1事象の電荷とドリフト時間の関係
• 解析方法
– 容器内PMTコインシデンスによりトリガー
– 電荷とドリフト時間の関係をプロットし、
τeを見積もる。
• 純度解析結果
実験を通して(7日間)、τe>300μsを達成!
→1ppb以下の純度を常に保持!
実験開始時 : τ≒670μs( 約 0.45ppb )
実験終了時 : τ≒385μs( 約 0.78ppb )
毎時 0.002ppb 程度の純度の悪化
05/13/2011
純度と液体Ar注入からの時間の関係
Preliminary Plot
τ = 671 – 1.90 ( t - t0 )
ここまでのまとめ(まだ終わりません)
• LArTPC検出器での核子崩壊探索感度の実験的検証のために
J-PARCにてビームテストを行った。
• Simulation data(K1.1BR 、宇宙線信号)を作成
• 宇宙線信号から液体アルゴンの純度を測定
– 全実験期間で必要な純度(<1ppb)を達成
• Simulation dataを用いて、Tracking Algorithm等の評価
–
–
–
–
ノイズの影響
崩壊点の決定
dE/dxの見積もり
入射粒子と崩壊粒子の決定
• 実データに、Tracking Algorithmを適用し、dE/dxを算出し、
K/π分離能力を評価する。
05/13/2011
13
10LArTPC開発研究
@早稲田大学
05/13/2011
14
暗黒物質探索
• 暗黒物質探索
2相型ArTPCで2種類のシンチレーション光
を捉えることにより、直接暗黒物質探索を
行うことが出来る。
2種類のシンチレーション光
 入射粒子との反応で発生
1次シンチレーション光 ・・・
S1
 電離電子が高電場中で加速されたときに発生
2次シンチレーション光 ・・・ S2
e,γ
WIMP
anode
-
e-e e- S2
grid
particle
drift time
( S2 / S1 ) WIMP
<<
e-
drift time
( S2 / S1 ) e,γ
cathode
E
PMT1
S1・S2の比で背景事象を取り除くことができる。
05/13/2011
ee- S1
PMT2
15
シンチレーション光観測のセットアップと信号例
 ガスAr中でのS1・S2信号の確認
anode
0.8cm
grid1
波長変換剤(TPB)を用い、Arのシンチレーション光
grid2
(128nm)を検出する。
S2
induction
field
低温で動作可能なPMTを使用。
MEG実験で使用されたものの改良版
使用したガス:純Ar(G1,純度99.9999%以上)
真空引きを行った後、ガスを流入。(常温、1気圧)
3.2cm
S1
drift field
0.7cm
e-
cathode
α線源
E
2.0cm
PMT self triggerで信号取得
PMT
信号例
induction
field
S2
grid1
S1
~8μs
grid2
α線源
Cathode(grid)
16
S2信号の電場特性
Drift E = 500 V/cm(fixed)
S2
それぞれの電場で500事象ずつデータを取得。
Averageを取りS2信号の電場特性を確認した。
S1
Induction電場とS2の関係
Induction電場 大
↓
Ar分子を励起する電離電子の数 増加
↓
S2が増加する!
Induction 電場とS2の関係
Induction E = 2.6 kV/cm(fixed)
Drift電場とS2の関係
Drift 電場 大
↓
電離電子の再結合 低下
電離電子のdrift速度 大
↓
S2が増加しピーク時間も早くなる!
………定量的な解析は現在遂行中である。
05/13/2011
S2
S1
Drift 電場とS2の関係
GEM + Wire読み出しセットアップ
• 暗黒物質探索を行うために、
微小な電離電子信号を取得したい。
→GEM+Wire読み出しのゲイン測定
• 使用したガス:純Ar(G1,純度99.9999%以上)
真空引きを行った後、
ガスを流入させる。(1気圧、常温、封じ切り)
• GEM間電圧 ΔV を変化させ、
それぞれ200事象ずつ信号を取得した。
grid1
Induction Field
GEM
ΔV
grid2
Drift Field
e-
e- eCondenser
(2.2nF)
Wire
Xe Flash Lamp
Cathode
Resistance
(100MΩ)
-HV
時定数τ=220msのハイパスフィルタ
• ガス増幅によるGainを次のように定義。
anode信号量
Gain 
cathode信号量
05/13/2011
18
Gain測定結果
(2)
• ワイヤー電圧は2kV/cm、2.5kV/cmとした。
GEM+ワイヤー読み出しによるGainの
ΔVGEM特性は右のようになる。
(1)
放電電圧
2kV/cm (1)
2.5kV/cm (2)
1010V
900V
(1) GEMによる放電
(2) 10L容器での放電
→クエンチャーが入れられないため 放電 しやすい。
放電電圧の直前までで、
Wire 電圧
2.0kV/cm : Gain~508 (ΔVGEM = 1000V)
2.5kV/cm : Gain~568 (ΔVGEM = 900V)
ΔVGEM
参照
現在までに、ガスAr中にて
GEM単体
Gain~89 (ΔVGEM=970V)
Wire単体
Gain~38 (2.5kV/cm)
を達成している。
05/13/2011
19
Summary & Future Plan
シンチレーション光
電離電子増幅
10L容器、ガス中にてS2を確認!
10L容器、GAr 常温・1気圧中にて、
GEM+Wire読み出しを用いて
Gain ~ 500を達成!
 Induction電場の増加に対する、
S2増幅を確認。
 Drift電場の増加に対して、
S2のピークが早くなり、
信号量が増幅することを確認した。
次の目標
気液2相型でS2を確認する!
次の目標
より高い増幅率に向けて、
 多段GEMなどを試す
 放電対策
※低温でのGain特性も調べる。
低温での高Gain・安定動作を目指す!
これらを組み合わせることで・・・
2相型TPCとして確立させ、暗黒物質探索を行いたい!
05/13/2011
20
Backup
05/13/2011
21
探索感度
22
Geant4 Simulation(Beam)
sx = 0.8 cm sy = 0.5 cm
@ BDC(FF,by simulation)
from TREK experiment
x_mean=5.2cm,σx=8.2cm
y_mean=0.2cm,σy=8.3cm
@the front of 250LAr Detector
(measurement by scintillation counter)
Beam
Beam Line
~5m
• BDCでのBeam分布から1点、250LAr検出器でのBeam分布から1点をランダムに
選び、直線を結ぶことによってビーム方向を決定した。
→Beamの拡がりをSimulationに実装
• この条件の上で、LArBDCを通過した事象をSimulation Dataとして選択する。 23
擬似データ作成
NQ
1.
2.
3.
粒子のエネルギー損失をストリップ毎に計算。
drift
W値から、生じる電離電子数を見積もる。
N0
ドリフト速度から、ドリフト時間(tDri)を計算。
4.
ドリフトし減少した後に各ストリップで読み出
される電子量(NQ)を見積もる。
5.
6.
7.
NQ と tDriから信号を正規分布で生成。
電子量から、ADCカウントに変換する。
#electron→charge→voltage→ ADC data
デッドチャンネルを考慮した。
ノイズを付加する。


ペデスタルの揺らぎ(σ=2~5ADC count)
or
実データを用いたノイズ
を導入。
40cm
ee e e e e eee
e ee ee e
e
e
e
e
K
e e e
ee
μ
Electric field
Cathode
電子量
N Q  N 0 exp  t Dri τe 
Anode
σ= 3.6μs
NQ
t
tDri
W value[eV]
Drift Velocity[mm/μs]
23.3
0.8
24
宇宙線信号例
シャワーが
トリガされた例
宇宙線がほぼ垂直
に抜けた例
(赤線:トリガタイミング)
ほぼ全てのチャンネル
にわたり、まっすぐな
飛跡が見えた例
多重散乱で飛跡
がまがっている例
δ線も見える
05/13/2011
25
ノイズ除去について
(1)生データ例
(2)FFTによる周波数強度 @ ch8
信号
ノイズ
200kHz以上に
大きなノイズ成分
(3)FFTフィルター前後 @ ch8
(4)FFTフィルタ後
- フィルタ前
- フィルタ後
高周波ノイズが大幅に改善
05/13/2011
飛跡がクリアに
26
ヒット・クラスタ探索手順
ヒット1
• 閾値を超える連続した
ACDカウントの数が
一定数を超えた場合
にヒットとする。
(下図の各BOX)
ヒット2
ヒット3
閾値
ヒット3
• 隣接するヒットをつなげ
クラスタを構成
(下図の各色)
05/13/2011
ヒット2
ヒット1
27
Comparison with data
Simulation (300事象の重ね合わせ)
0
0
Time[μs]
Real data (300事象の重ね合わせ)
300
200
100
0
20
40
60
TPC Channel
20
40
60
TPC Channel
• Kのデータとの比較
– LAr中のビームの拡がりや,ビームの到達距離などが、
Simulationによって、おおよそ再現できている。
• 今後、このシミュレーションを用いてTracking Algorithmの評価
を行っている(現在進行中)。
28
気液2相型のセットアップと結果 anode
 気・液2相のAr中でのS1・S2信号の確認
grid1
extraction
4.4×10-4Paまで真空を引いた後、容器を冷やし
field grid2
ガスArを流入、10L容器を1気圧に保ったまま液化。
• 液面をextraction領域に保持。(誤差1mm以下)
(液面計・白金抵抗)
drift field
cathode
S1
extraction E
6kV/cm
S2信号が見えていない。
05/13/2011
S1
extraction E
0 V/cm
白金抵抗1
0.8cm
0.7cm
白金抵抗2
3.2cm
α線源
•extraction領域に高電場をかける。
(~5kV/cm : 電離電子を気相に取り出すため )
液面 : grid2 から 0.3cm
drift電場 : 900V/cm 固定
液面
2.0cm
PMT
PMT self triggerでそれぞれ
1000イベント データを取得し、
Averageを取った。
純度が悪いため電離電子が液中で
減衰しS2が見えていないと思われる。
今後純度を改善し、2相型で
電子を気相に取り出してS2
S1信号の遅い成分が見えていない。
信号を確認する。
(純度が悪いと考えられる。)
JINST 5:P05003,2010
E
29
10LArTPC中のPMT
Hamamatsu R6041-06MOD
変換剤:TPB(テトラフェニルブタジエン)
05/13/2011
30
ドリフト電場とドリフト速度の関係
31
シンチレーション光(S1・S2に関して)
S2
S1
S1
anode
S2
0.8cm
grid1
S2
grid2
e光電効果
3.2cm
S1
ecathode
0.7cm
α線源
E
2.0cm
PMT
ある頻度でS1がgridに当たり光電効果を起こすと考えられる。
05/13/2011
32
電離電子増幅
暗黒物質探索を行うために、微小な電離信号を取得したい。
5kV/cm
反跳エネルギーを50keVと仮定すると、
生成される電離電子は、~10個ほど (by ArDM group)
→ 電荷量1.6e-3[fc]
PreAmp(Amptek A250)
input charge 1fC → output height 1mV
1kV/cm
Recoil Energy (MeV)
暗黒物質探索のためには、Gainが~1000必要!
増幅のためのツール
厚型GEM (T-GEM-100-400/700)
穴径 : 300 μm
ピッチ : 700 μm
厚み : 400 μm
自作の読み出し基板
材質 : 金メッキタングステン(金3%)
直径 : 30 μm
ピッチ : 5 mm
張力 : 0.45 N
GEM + Wire読み出しセットアップ
Wire – grid1
6mm
2 or 2.5kV/cm
GEM(top) – Wire
6mm
2 or 2.5kV/cm
GEM(bottom – top)
ΔV = variable
grid2-
GEM(bottom)
7.5mm
500V/cm
Cathode – grid2
10mm
200V/cm
05/13/2011
34