Difference between revisions of "Alpha Ionization"

From New IAC Wiki
Jump to navigation Jump to search
Line 61: Line 61:
 
=Number of Primary  ionization events per cm=
 
=Number of Primary  ionization events per cm=
  
,Inside the chamber
+
;Inside the chamber
  
 
[[File:PhotonEnergy_numberof_primary_electrons_Drift.png | 200px]]
 
[[File:PhotonEnergy_numberof_primary_electrons_Drift.png | 200px]]
Line 69: Line 69:
 
The figures show the threshold energy for photons travel in the detector and the number of the absorbed photons within the detector's drift area. Photons travel through the the drift area and mainly interact by photoelectric effect with the gas's atoms and molecules; a photon's energy may increase up to 1.1 MeV and still the most probable interaction with the medium is the photoelectric effect. Figure xx1 shows the same numbers of primary electrons produced by an incident photon simulation by GEANT4; once when all photon's interactions are considered (red), and when only the photoelectric effect only considered (blue). It also shows the all the photons are absorbed by Ar/CO2 gaseous medium. Since the incident photon energy varies from 10 keV up to 1.1 MeV, only the absorbed photons within a 1 cm will affect the detector signal. GEANT4 simulation shows the dominant interaction is the photoelectric effect even if the incident photon energy reaches to 1.1 MeV, and 25 percent of the  photons will be absorbed when their energy is 25 keV, as the energy increases the number of photons decreases to reach 0.03 when the photon energy is 96.1 keV. As a result, gamma energy, emitted by U-233, or bremsrahlung photons, appear when the lineac is in operation, should have a maximum energy of 96.1 keV as it is in the detector's drift area to leave a clear effect on the detector's signal.  
 
The figures show the threshold energy for photons travel in the detector and the number of the absorbed photons within the detector's drift area. Photons travel through the the drift area and mainly interact by photoelectric effect with the gas's atoms and molecules; a photon's energy may increase up to 1.1 MeV and still the most probable interaction with the medium is the photoelectric effect. Figure xx1 shows the same numbers of primary electrons produced by an incident photon simulation by GEANT4; once when all photon's interactions are considered (red), and when only the photoelectric effect only considered (blue). It also shows the all the photons are absorbed by Ar/CO2 gaseous medium. Since the incident photon energy varies from 10 keV up to 1.1 MeV, only the absorbed photons within a 1 cm will affect the detector signal. GEANT4 simulation shows the dominant interaction is the photoelectric effect even if the incident photon energy reaches to 1.1 MeV, and 25 percent of the  photons will be absorbed when their energy is 25 keV, as the energy increases the number of photons decreases to reach 0.03 when the photon energy is 96.1 keV. As a result, gamma energy, emitted by U-233, or bremsrahlung photons, appear when the lineac is in operation, should have a maximum energy of 96.1 keV as it is in the detector's drift area to leave a clear effect on the detector's signal.  
  
,7cm Outside the chamber
+
;7cm Outside the chamber
  
  

Revision as of 16:16, 9 March 2014

Energy threshold for electron and gamma entering detector

Detemine the energy threshold that an electron and a photon need to surpass in order to pass through the kapton, gas, the copper cathode, and the FR4.

Draw a picture with dimensions of all the layers


Energy -vs- counts per incident particle


Photon and electron Energy -vs- distance through the detector


energy straggling plot at the point where particle just makes it through the detector

Alpha Particles ionization simulation using GEANT4

GEANT4 simulates the ionization of alpha particles in Ar/CO2 90/10 gas. Geant4 can simulate the ionization process for alpha particles. Unfortunately the value of the step function underestimates the number of delta electrons even after decreasing the step cut to 1 nm. Also,using GEANT4 overestimates the range of alpha particles in Ar/CO2 gas when compared to those that srim calculates File:Alpha range ArCo2.txt, the following table shows the maximum range of alpha particles that are emitted from the U-233, and the ranges calculated by srim.

Alpha Energy (MeV) G4 Range (cm) Srim Range (um)
1.0 0.56599 129.49
2.0 1.1467 255.91
3.0 1.9024 417.27
4.0 2.8012 612.45
5.0 3.8425 839.91

Based on the previous table, GEANT4 failed to calculate the expected alpha range for most alpha energies, and underestimated the number of alpha's delta electrons emitted through that range, but it is still useful tool to simulate the primary delta electrons when negative beta particles pass through a defined medium.

Calculating the number of the delta electrons without using GEANT4

There is another way to calculate the number of delta electrons without using GEANT4. It starts by calculating the average energy loss [math] \Delta E_a[/math] by the alpha particles and the average energy loss per unit length [math] {\zeta} [/math] in Ar/CO2 gas using the Bethe-Block equation. It then uses the the following equation:

[math] \lambda = \frac{ \Delta E - \Delta E_a}{\zeta} [/math]

to calculate the actual energy loss by ionization [math] \Delta E [/math], where [math] \lambda [/math] represents random landau number.

By dividing the energy loss by the minimum energy for producing a pair of ion/electron pair W, this equation yields the number of electrons emitted by ionization.

Srim can simulate the motion of an alpha particle in Ar/CO2 gas by allowing for the change in the stopping power per unit length. It can also show the ionization energy loss as shown in the following figure:

4.8MeV ionization ArCo2.jpeg

This figure shows the ionization energy loss of 999 alpha particles as they pass through Ar/CO2 gas. Integrating the area under the curve and dividing by (W*1000) will estimate the number of delta electrons produced by an alpha particle.

Bethe block formula calculates the average ionization energy loss as the following: [math]- \frac{dE}{dx} = \frac{4 \pi}{m_e c^2} \cdot \frac{nz^2}{\beta^2} \cdot \left(\frac{e^2}{4\pi\varepsilon_0}\right)^2 \cdot \left[\ln \left(\frac{2m_e c^2 \beta^2}{I \cdot (1-\beta^2)}\right) - \beta^2\right][/math]

for an electron of rest mass [math] mc^2 = 0.511 MeV [/math],

Number of Primary ionization events per cm

Inside the chamber

PhotonEnergy numberof primary electrons Drift.png PhotonEnergy numberofphotonsabsorbedwithin1cm.png


The figures show the threshold energy for photons travel in the detector and the number of the absorbed photons within the detector's drift area. Photons travel through the the drift area and mainly interact by photoelectric effect with the gas's atoms and molecules; a photon's energy may increase up to 1.1 MeV and still the most probable interaction with the medium is the photoelectric effect. Figure xx1 shows the same numbers of primary electrons produced by an incident photon simulation by GEANT4; once when all photon's interactions are considered (red), and when only the photoelectric effect only considered (blue). It also shows the all the photons are absorbed by Ar/CO2 gaseous medium. Since the incident photon energy varies from 10 keV up to 1.1 MeV, only the absorbed photons within a 1 cm will affect the detector signal. GEANT4 simulation shows the dominant interaction is the photoelectric effect even if the incident photon energy reaches to 1.1 MeV, and 25 percent of the photons will be absorbed when their energy is 25 keV, as the energy increases the number of photons decreases to reach 0.03 when the photon energy is 96.1 keV. As a result, gamma energy, emitted by U-233, or bremsrahlung photons, appear when the lineac is in operation, should have a maximum energy of 96.1 keV as it is in the detector's drift area to leave a clear effect on the detector's signal.

7cm Outside the chamber


PhotonEnergy numberofphotonsabsorbedwithin7cm out.png


from beta particle as a function of beta particle energy using Garfield and GEANT4

GEANT4 can simulate the primary number of delta electrons for a negative beta particle penetrating ArCO2. Using TestEm10 example and choosing appropriate default cut, GEANT4 counts the same number of primary delta electrons for a 1.1 MeV negative beta particles, and determines the energy of the delta electrons energy and momentum depending on Moller or Bhabha's scattering depending on the value of the kinetic energy cut used.<ref name = "Urban">Urbán, L. (1998, 10 09). Geant4 physics reference manual. Retrieved from http: //geant4.web.cern.ch/geant4/G4UsersDocuments/UsersGuides/PhysicsReferenceManual/html/node41.html </ref>

The following figure shows the number of the primary electrons that are emitted from negative beta with energy E, passing through a 90/10 Ar/CO2 gas:


Neg beta vs numberof primary electrons Garf.png G4 electron primaryeN.png Garf G4 Pelec ion.png

Overlay with dE/dx formula and fill in gaps

something is wrong with G4 result for electron energies below 1 keV.

from photon as a function of photon particle energy using Garfield and GEANT4

Photon eprimary.png

Photoionization Ar.png


Chamber Geometry in GEANT4

The chamber geometry for the GEM chamber is set for Calculating the created charge in the detector's drift area is shown below:

G4 chamber design.png

Charge particles and photons primary ionization

The physical interactions that take place inside the detector is simulated by GEANT4.9.6 to estimation the primary charge collected by the drift electric in GEM detector. GEM detector depends mainly on the particles interaction with Ar/CO2 atoms and molecules, especially the ones that directly (or indirectly) produce free charge that can be collected by the drift electric field. Specifically for the GEM detector described in section XX, the detector signal can be produced by gamma rays, electrons, or neutrons. GEANT4.9.6 has been used to estimate the primary charge caused by those particles in case of electrons and neutrons, but for photons, GENAT4 is used to produce the photoionization energy spectrum for Ar/CO2 gas.

The simulation started by determining the number of primary electrons produced by charged particles in Ar/CO2. The number of the primary electrons depends on the particle's charge and mass and kinetic energy. As a result alpha particles have the highest number of primary electrons, electrons, then photons which can only produce an electron per photon assuming the photoabsorption is dominant process in Ar/CO2 gas for an energy range of 1 to 300eV.

By modifying EmTest10 example in GEANT4.9.6, the number of primary electrons produced by an electron travels 1 cm is shown in the figure below:

Garf G4 Pelec ion.png

The figure also shows Garfield simulation for an electron's primary ionization, it is noticed that GENAT4.9.6 results meets with those of Garfield when the minimum incident electron energy is [math] 2x10^4[/math] eV. Such an energy is much lower than the threshold kinetic energy which an electron needs to make an ionization in the detector's drift region (the threshold kinetic energy for an electron is more than 2 MeV 5cm away from the detector's window).

GEANT4.9.6 predicts the number of the primary electrons for a photon. Photon ionization has the least number of electrons, even producing an free electron is constrained to the electronic structure of Ar/CO2 gas.

Photon eprimary.png Photoionization Ar.png

GENAT4.9.6 does not predict the Argon photoabsorption spectrum specifically for each energy as shown above, the photon energy may have [math] \pm 1[/math]eV shift. It is worth to mention if the photon is not in the drift region, then the detector is completely blind for any photon of energy of 20-300 eV (their range in fm). In case of a 50 keV photon or more, the photon will easily pass the detector chamber without affecting its medium if only photoionization is considered, within that energy range Compton scattering cross section increases, it will produce an electron and a photon in the drift region that can easily be detected. Being in am accelerator environment, photon interaction becomes dominant by Compton Scattering interaction in a way that makes the detector in efficient to detect any other particle.

Total charge on anode

Now use the primary above and include detector gain to calculate total charge hitting the anode and getting readout.


<References/>