the simulation of gas avalanche in a micro pixel chamber...
TRANSCRIPT
Va=460Vr=0.31
The Simulation of Gas Avalanche in a Micro Pixel Chamber using Garfield++
A. Takada, T. Tanimori, H. Kubo, J. D. Parker, T. Mizumoto, Y. Mizumura, S. Iwaki, T. Sawano, K. Nakamura, K. Taniue, N. Higashi, Y. Matsuoka, S. Komura, Y. Sato, S. Namamura, M. Oda,
S. Sonoda, D. Tomono, K. Miuchi1, S. Kabuki2, Y. Kishimoto3, S. Kurosawa4
Kyoto University, 1Kobe University, 2Tokai University, 3KEK, 4Tohoku University
1. Micro Pixel Chamber (μ-PIC)Characteristics of μ-PIC[1]
A gaseous 2D imaging detector with strip read outCu electrodes and polyimide substrateEach pixel is placed with a pitch of 400 μmLarge detection area 10×10 cm2, and 30×30 cm2
Based on the Print Circuit Board technologyMaximum gas gain of ~15,000Fine position resolution (RMS ~120µm)Good gas gain uniformity (RMS ~5%) at the whole areaStable operation for >1000 hours with gas gain ~6,000
Applications of μ-PIC are … Electron-Tracking Compton Camera
in astronomy[2], Medical imaging[3]Neutron imaging[4]Direct dark matter search[5]X-ray crystallography[6]Gas photomultiplier[7]
Time Projection Chamber with TOT
4. Signals of μ-PIC
Recently, more precise tracks are required by some applications, we therefore developed a new data acquisition system for a Time Projection Chamber (TPC) using Time-Over-Threshold (TOT). By the TOT-TPC, we can already obtain the high resolution tracks, like the following figures.
A simulator, which is necessary for the accurate study of application detectors, requires the wave form, the uncertainty of gas gain, and the dependence on the anode voltage. For the purpose, we should simulate the avalanche in microscopic scale.
Needs from Applications
2. Geometry & Electric fieldLet’s simulate with Garfield++[8] !!
Definition of GeometryWe defined the geometry of a pixel electrode as a unit cell in the area of 400×400 μm2, and generated the three dimensional finite element mesh using Gmsh[9]. The gas area had the depth of 1.5 mm above the electrodes, and was filled by Ar 90% + Ethane 10% with the pressure of 1.0 atm. In the Garfield++, we placed this unit cell periodically. The initial electrons generated at random position in the area of 400×400 μm2 at 1.0 mm above the electrodes.
Calculation of Electric Field
Avalanche Size
Single Electron Spectrum
Map of Termination Points
Wave Form of Single Seed Electron
For the calculation of the maps of electric field, voltage, and weighting field, we adopted Elmer[10]. The drift electric field was fixed at 1.0 kV/cm. The dielectric constants of gas, polyimide substrate, and copper electrode are 1, 3.5, and 1010, respectively. To define these maps in Garfield++, we used ComponentElmer class[11].
[V] [V/cm]
VoltageVA=560 V
Electric fieldVA=560 V
[ 1] A. Ochi+, NIM A, 478 (2002), 196; T. Nagayoshi+, NIM A, 525 (2004), 20; A. Takada+, NIM A, 573 (2007), 195.[ 2] T. Tanimori+, New Astron. Rev., 48 (2004), 263; A. Takada+, ApJ, 733 (2011), L13.[ 3] S. Kabuki+, NIM A, 623 (2010), 606.[ 4] J. D. Parker+, NIM A, 697 (2013), 23.[ 5] K. Miuchi+, Phys. Let. B, 686 (2010), 11.[ 6] K. Hattori+, J. Syncrotron. Rad., 16 (2009), 231.[ 7] H. Sekiya+, JINST, 4 (2009), P11006.[ 8] http://garfieldpp.web.cern.ch/garfieldpp[ 9] http://geuz.org/gmsh
5. Summary
References
First, we simulated the avalanche size by counting the number of generated electrons in an avalanche using microscopic tracking. The obtained avalanche size was shown in the Figure a) with various anode voltage and Penning effect rate 𝑟𝑟. For comparison, the measured effective gas gain was also plotted in this figure (open squares and filled triangles). We fitted these plots with
𝐹𝐹𝑟𝑟(𝑥𝑥) = exp(𝛼𝛼 + 𝛽𝛽𝑥𝑥)at each 𝑟𝑟, and the dependences of the obtained 𝛼𝛼 and 𝛽𝛽 on 𝑟𝑟 are shown in the Figures b) and c). In the case of Ar 90% + Ethane 10%, the Penning effect rate is 0.31 [12], we hence interpolated 𝛼𝛼 and 𝛽𝛽under the assumption that the dependences of 𝛼𝛼 and 𝛽𝛽 are the linear functions of 𝑟𝑟. The calculated avalanche size using the interpolated 𝛼𝛼 and 𝛽𝛽 is represent the blue solid line in Figure a), and it is well consistent with the measured value.
Parameter Value
Anode
Pillar diameter 70 μm (upper)50 μm (lower)
cap Diameter: 60 μmHeight: 15 μm
Strip Width: 285 μmThickness: 15 μm
CathodeOpening Diameter: 256 μm
Strip Width: 340 μmThickness: 15μm
SubstrateThickness 75 μmmaterial polyimide
gasElectric Field 1 kV/cm
material Ar 90% + C2H6 10%, 1 atm
3. Avalanche in μ-PIC
[10] http://www.csc.fi/english/pages/elmer[11] Private support by J. Renner.[12] Ö. Sahin+, JINST, 5 (2010), P05002.[13] T. Nagayoshi+, NIM A, 546 (2005), 457.[14] O. Sasaki+, IEEE TNS, 46 (1999), 1871.[15] R. Orito+, IEEE TNS, 51 (2004), 1337.[16] S. Agostinelli+, NIM A, 506 (2003), 250.
In the case of a proportional counter, the single electron spectrum, which represents the uncertainty of gas gain, is described by Polyadistribution
𝑔𝑔 𝑥𝑥 = 𝑥𝑥(1+𝜃𝜃)�𝐴𝐴
𝜃𝜃exp −𝑥𝑥(1+𝜃𝜃)�𝐴𝐴 ,
where �̅�𝐴 is the average of gas gain. The single electron spectrum of μ-PIC is shown in right figures, and it is fitted by 𝑔𝑔(𝑥𝑥). From these figures, it can be explained with Polya distribution, although the electric field map of μ-PIC is complicated in comparison with that of a proportional counter. The theoretical limit of energy resolution is obtained using Fano factor 𝑓𝑓 and 𝜃𝜃 as follows:
𝜎𝜎𝐸𝐸𝐸𝐸 = 𝑓𝑓+
11+𝜃𝜃�𝑁𝑁 ,
where �𝑁𝑁 is the average number of seed electrons. 𝜃𝜃 of μ-PIC is approximately 0.65, which is independent on the anode voltage.
The left figures are maps of generated point of the secondary electrons at the anode voltage of 460 V and 560 V, respectively. The color indicates means the number density of secondary electrons, and the contour means the number density of the gravity center of each avalanche. The majority of avalanche is made in the limited area at the distance of 100 μm from anode electrode, and the center of gravity of avalanche are concentrated above anode
Emulation of X-ray Irradiation
Time over Threshold
The left figure is the termination points map of all electron. The initial electrons were made at random position in the red box above the substrate at the distance of 1 mm. 2% of the electrons in avalanche are terminated on the substrate, and the electron collection efficiency of μ-PIC is 98%. This result is roughly consistent with the previous simulation result[13] which said that most of electrons reach to the anode electrode and the rest drift onto the substrate.
electrode at the distance of approximately 5 μm. Because the avalanche points are almost not changed by the anode voltage, the signal waveform is insensitive to the anode voltage.
We defined the three dimensional finite element mesh using Gmsh, and calculated the map of voltage, electric field, and weighting field with Elmer.Using Garfield++, we simulated the avalanche in μ-PIC.The simulated avalanche size is well consistent with the measured gas gain.The single electron spectrum of μ-PIC is described by a Polya distribution.The electron collection efficiency is roughly consistent with the previous simulation result.We also obtain the wave form of single electron incident.The pulse height dependence on the time constant of preamplifier is well explained.With Geant4 and ASD emulator, we can roughly explain the TOT hits distribution.By these results, the simulation for the applications will be improved. Moreover, the Garfield++ simulation will make the μ-PIC development more easily.
We simulated the induced current using the weighting field of anode electrode. The left figure is a sample of single electron signal. The signal waveform consists of two components. One is sharp spike induced by electrons, and the other induced by ions has slow decay time.
Ion componentElectron component Electron component
0 100 200 300 400Time [ns]
0
-0.02
-0.04
Curr
ent
[μA
]
1
10
102
103
0 0.4 X [mm]-0.4-0.8-0.8
-0.4
0
0.4
Y[m
m]
0-0.01
Anode: 460Vr = 0.31
-0.02Charge induced by electron [fC]
-0.03
Ava
lanc
he S
ize
0
1000
2000
3000
4000
5000
• The charge induced by electron is proportional to the avalanche size.
• The pulse width of electron component is approximately 1 ns.
0-0.1-0.2Charge induced by ion [fC]-0.3
Ava
lanc
he S
ize
0
1000
2000
3000
4000
5000
Ion component• The charge induced by ion is proportional to the
avalanche size, and it is approximately 9 times larger than the charge of electron component.
• As like the left figure, the ion component takes various pulse shape, which depends on the ion transportation paths.0 100 200 300
Time [ns]
0
-2
-4
Curr
ent/
Ava
lanc
he S
ize
[pA
]
Only ion component
Anode: 460Vr = 0.31
Anode: 460Vr = 0.31
□: SN040426-1▲: SN050921-1
Ava
lanc
he s
ize
102
103
104
105
460 480 500 520 540 560Anode Voltage [V]
α β
0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6
0.028
0.020-3.4
-3.0
0 40 80 120Radius [μm]
0
40
80
Hei
ght
[μm
]
CathodeAnode
1
10
0 40 80 120Radius [μm]
0
40
80
Hei
ght
[μm
]
CathodeAnode
1
10
102
103
[#][#] Avalanche AreaVa=560Vr=0.31
In order to emulate the signal output of u-PIC, we simulated 2000 single electron signals as the wave form templates. Using this template, we calculated with the processes as follows:
1. Create seed electron with a Gaussian having mean of 𝐸𝐸𝑤𝑤 and sigma of ⁄𝑓𝑓𝑓𝑓 𝑤𝑤, where 𝐸𝐸, 𝑓𝑓, and 𝑤𝑤 are the deposit energy, Fano factor, and w value of gas, respectively.
2. Define the uncertainty of arrival time using the longitudinal diffusion for each seed electron.3. Because a μ-PIC has the fluctuation in the gain uniformity, define the average avalanche size �̅�𝐴 with a
Gaussian having the sigma of 𝜎𝜎�̅�𝐴, for each pixel.4. Obtain the avalanche size 𝐴𝐴 using a Polya distribution with mean of �̅�𝐴 and 𝜃𝜃 of 0.65 for each seed electron.5. Select a wave form from 2000 signal templates, and multiplied by 𝐴𝐴.6. Obtain an output signal of u-PIC with summing the single electron signal for all seed electron.7. Emulate the readout circuit, and compare with experimental data.
r r
0.024-3.2
-3.6
-2.8
0.018
a)
b) c)
0 0.1 0.2X [mm]
0.3-0.1-0.2
0
-0.1
-0.2-0.3
0.1
0.2
Z[c
m]
1
10
102
103
104
105
0 0.1 0.2X [mm]
0.3-0.1-0.2-0.3
0
-0.1
-0.2
0.1
0.2
Z[c
m]
0
100
200
300
400
500
20 mV
10 mV
400 ns
55Fe irradiation (5.89 keV)Ar 80% + C2H6 20%
1 atm, Gain ~3000
16 nsec ASD
80 nsec ASD
Anode: 460Vr = 0.31
𝐸𝐸 = 5890 eV𝑓𝑓 = 0.17, 𝑤𝑤 = 26 eV�̅�𝐴 = 3000, 𝜎𝜎�̅�𝐴 = 7%𝜃𝜃 = 0.65, 𝜏𝜏 = 16 ns
𝐸𝐸 = 5890 eV𝑓𝑓 = 0.17, 𝑤𝑤 = 26 eV�̅�𝐴 = 3000, 𝜎𝜎�̅�𝐴 = 7%𝜃𝜃 = 0.65, 𝜏𝜏 = 80 ns
0 200 400 600 800Time [ns]
0
-10
-20
-30
Emul
ated
pre
ampl
ifie
r ou
tput
[mV]
200 400 600 800Time [ns]
0
-10
-20
-30Emul
ated
pre
ampl
ifie
r ou
tput
[mV]
0
100 pF
1 pF
16 kΩ
38 pF
VA
1 GΩ
As a readout, we use an amplifier-shaper-discriminator (ASD) chip having a time constant of 16 ns[14], which was designed for the Thin Gap Chamber in ATLAS, or a redesigned ASD having a time constant 80 ns[15]. By the measurement with X-ray irradiation of 55Fe, the pulse height obtained by 80 ns ASD is about two times higher than that of 16 ns ASD. We tried to confirm this pulse height difference. The emulated signal pulses are shown in below. At the gas gain of 3000, the emulated pulse heights of 16 ns preamplifier and 80 ns preamplifier are 11 mV and 23 mV, respectively. The ratio of pulse height is approximately 2.2, which is well consistent with the measurement.
100 pF
1 pF
80 kΩ
38 pF
VA
1 GΩ
Cosmic muon
104Va = 460 Vr = 0.31
recoil electronin Compton scattering
scattering
stoppedδ-ray
Emulated pulse Emulated pulse
Finally, we simulated the TOT distribution depending on the deposit energy using the hit pixel number of the TOT-track images.The TOT distribution of electrons in the experimental data has two components. One is the component of fully contained events, and another is the component of escape events. At the energy of 40 keV and 100 keV, the typical TOT hits for fully contained electrons are about 500 pixels and 1500 pixels, respectively. For the purpose of emulating the TOT-tracks, we made a simulator of μ-PIC TPC using Geant4[16]. By this simulator, we obtained the energy deposit of electrons, calculated signals at each pixel, emulated digital pulses, and compared with the experimental TOT tracks. The emulated TOT distribution also has two components, and the number of emulated TOT hits is nearly equal to that of experimental data at each energy. Therefore, we can say that this emulator well describes the TOT distribution, and we will be able to realize a precise simulator forμ-PIC application.
103
102
TPC Size:Gas:Gas gain:Readout:Threshold:Source:Drift field:Induction:
TPC Size:Gas:w value:Fano factor:Gas gain:Readout:Threshold:Electron energy:Drift velocity:Diffusion:
7.5×7.5×14 cm3
Ar 90% + C2H6 10%, 1.5 atm24000 (μ-PIC + GEM)ASD (τ = 80 ns)-15 mV137Cs irradiation170 V/cm1.1 kV/cm
10×10×14 cm3
Ar 90% + C2H6 10%, 1.5 atm26 eV0.1724000 (θ = 0.65)ASD (τ = 80 ns)-15 mV< 200 keV3.2 cm/μs450 μm/√cm for transverse350 μm/√cm for longitudinal
0
1000
2000
3000
0
1000
2000
3000
0 40 80 120Deposit energy [keV]
0 40 80 120Deposit energy [keV]
1
10
102
1
10
102
Deposit energy [keV]
0
1000
2000
3000
0 40 80 120Deposit energy [keV]
Experiment setup Simulation parameters
All Events (experiment)
Fully Contained Events(experiment)
Fully Contained Events (simulation)
1
10
ToT
0
1000
2000
3000
0 40 80 120
All Events (simulation)
1
10
ToT
ToT
ToT
102
102
Experiment
Experiment
CathodeAnode (cap)
Anode (strip)
Samples of various waveforminduced by ions
Electron component
Ion compoment
Enlarged imagearound anode
Enlarged imagearound anode