Neutron Interaction with Germanium Detector
050712 Xinjie Qiu
Updated 050715, 050718, 050719, 050720, 050727*
During the simulation of the charge sharing effect in the CDMS ZIP detector, some processes were not understood. In order to better understand the neutron interaction with Germanium Detector, it was investigated by the Geant4 Mote Carlo Simulation .
Germanium Detector in the simulation:
The actual size of the ZIP detectors in the CDMS experient are 1 cm thick and 3 inches(7.62 cm) in diameter. In order to increase the interaction event number, an artificial huge size detector in the simulation has the folloing dimension: 100m in length, 70m in diameter. Neutron beam from the top incident vertically onto the detector. Around the detector is vaccum, and the temperature is set to 50 mK.
There naturely exists 5 diffirent Germanium isotopes, all of them are found in the simulation.
Name: germanium
Symbol: Ge
Atomic number: 32
Atomic weight: 72.64
| Isotope |
Atomic mass (ma/u) |
Natural abundance (atom %) |
| 70Ge |
69.9242497 (16) |
20.84 (87) |
| 72Ge |
71.9220789 (16) |
27.54 (34) |
| 73Ge |
72.9234626 (16) |
7.73 (5) |
| 74Ge |
73.9211774 (15) |
36.28 (73) |
| 76Ge |
75.9214016 (17) |
7.61 (38) |
Neutron Interaction with nucleus:
Neutrons have no electrical charge and have nearly the same mass as a proton. A neutron is hundreds of times larger than an electron, but one quarter the size of an alpha particle. Because of its size and lack of charge, the neutron is fairly difficult to stop, and has a relatively high penetrating power.
There are five reactions that can occur when a neutron interacts with a nucleus[1,2]. In the first two, known as scattering reactions, a neutron emerges from the reaction. In the remaining reactions, known as absorption reactions, the neutron is absorbed into the nucleus and something different emerges.
| 1. Elastic Scattering -- (n, n) |
2. Inelastic Scattering -- (n, n γ) |
A neutron collides nucleus, transfers some energy to it, and bounces off in a different direction. Target nucleus gains the energy lost by the neutron, and then at an increased speed.
|
A neutron strikes a nucleus and is temporarily absorbed, forming a compound nucleus. This will be in an excited state. It then de-excites by emitting another neutron of lower energy, together with a gamma photon, which takes the remaining energy.
|
 |
 |
3. Transmutation -- (n, p), (n, α)
A nucleus may absorb a neutron forming a compound nucleus, which then de-energizes by emitting a charged particle, either a proton or an alpha particle. This produces a nucleus of a different element. Transmutation is the transformation of one element into another by a nuclear reaction.
| (n, p) | (n, α) |
| |
 |
 |
4. Radiative Capture -- (n, γ)
This is the most common nuclear reaction. The compound nucleus formed emits only a gamma photon. In other words, the product nucleus is an isotope of the same element as the original nucleus. Its mass number increases by one.
5. Fission
While bombarding heavy nuclei with neutrons, a highly excited compund nucleus might be formed. It might split into two nuclei of medium mass. During this process energy is produced and more neutrons are released.
Neutron Interaction with Germanium detector at diffirent energy scale:
The following picutres are view from side. Two red lines on the left and right sides are the top and bottom of the detector. The trajectories of diffirent particles are labeled by color according to their charges. Neutral particles are green, negative red, and positive blue. Most of the green lines in the following plots are nuetrons and gama rays , red lines are electrons, and blue lines are protons or alpha particles(recoil nuclei have very short moving distance before they dessipate all their enery into the crystal, so their trajectories are hardly to see).
It is worth to mention that the neutrons produced by radioactivity in the rock have energyies of a few MeV, while the neutrons from muon interactions in the surrounding rock cavern have energies well above 50 MeV. So We are mainly interested in the discussion of 2MeV and 20 MeV neutrons.
| Incident neutron energy 20keV |
Typical interactions: |
Show plots for:
|
| almost all are | n + Ge --> n + Ge (elastic interaction) |
| multiple scattering: | 2,3,4,5,6,7,8,9 |
|
| Incident neutron energy 200keV |
Typical interactions: |
Show plots for:
|
| most are | n + Ge --> n + Ge (elastic interaction) |
| few | n + Ge --> n + Ge + γ (inelastic interaction) |
| multiple scattering: | 2,5,6,7,9 |
|
| Incident neutron energy 2MeV |
Typical interactions: |
Recoil energy histogram |
Show plots for:
|
| half | n + Ge --> n + Ge (elastic interaction) |
| half | n + Ge --> n + Ge + γ (inelastic interaction) |
| multiple scattering: | 3,6 |
|
 |
| Incident neutron energy 20MeV |
Typical interactions: |
Show plots for:
|
| half | n + Ge --> n + Ge (elastic interaction) |
| half | n + Ge --> n + Ge + γ (inelastic interaction) |
| few | n+ Ge --> n + p + Ga + γ (transmutation) |
| multiple scattering: | 8 |
|
 |
| Incident neutron energy 200MeV |
Typical ineractions: |
Show plots for:
|
| some | n + Ge --> n + Ge (elastic interaction) |
| some | n + Ge --> n + Ge + γ (inelastic interaction) |
| some | n+ Ge --> n + p + Ga + γ (transmutation) |
| multiple scattering: | 4,7,8 |
|
| Incident neutron energy 2000MeV |
Typical ineractions: |
Show plots for:
|
| some | n + Ge --> n + Ge (elastic interaction) |
| some | n + Ge --> n + Ge + γ (inelastic interaction) |
| lots of | n+ Ge --> n + p + XYZ + γ (transmutation) |
|
General structure of neutron cross sections[3]:
1. Most neutron absorption occurs in low energy (< 1 eV) range. The absorption cross section in this range tends to be inversely proportional to neutron velocity (i.e., 1/v). Germanium provides extremely high thermal neutron absorption cross sections.
2. The energy range from about 1 eV to the first inelastic scattering threshold (which is isotope-dependent) is dominated by elastic scattering , except in the vicinity of resonances.
3. Above the energy threshold for the lowest level of inelastic scattering, this mechanism tends to dominate the slowing-down processes, especially for heavier nuclides. This is true even if the cross section value for inelastic scattering is lower than elastic, because of the higher fractional energy loss associated with inelastic scatter. The secondary neutrons that emerge from inelastic scattering reactions tend to be isotropic in the COM system due to the formation of a compound nucleus -- and the subsequent "forgetting" of the original direction of the incoming neutron.
4. For details about the neutron interaction Q-values and threshold of each isotope of Germanium, see appendix.
Conclusions:
At low energy(10 keV), the interaction are more likely to be elastic interaction. The multiple neutron interaction with Germaninum occurs at such energy range. As energy becomes higher and higher(~MeV), more and more interactions are inelastic. Above 10MeV energy level, nuclear interactions happen, mostly transmutation, produce protons, which will kick off the electron from the Germanum atoms to produce free electrons. No fission is oberserved.
Reference
[1] http://canteach.candu.org/library/20040706.pdf Neutrons and Neutron Interactions
[2] http://canteach.candu.org/library/20040707.pdf Fission
[3] http://web.utk.edu/~rpevey/NE406 lesson 8 and lesson 10
[4] http://t2.lanl.gov/data/qtool.html Q-values and Thresholds
Update changes
050715 correct the wrong energy level, more simulation performed at diffirent energy levels.
050718 general neutron nucleus interaction explaination and pictures added, interaction type classified for each enery level, conclusion and reference added
050719 neutron scattering and absorption properties table added
050720 Q values and threshold energy added.
050727 e-notes structure changed, some minor errors corrected, linked to charge sharing notes.
Appendix
Neutron interaction Q-values and Thresholds[4]:
For neutrons, the interaction is with the target nucleus itself. It is possible, therefore that some of the kinetic energy of the particles will be used to excite the nucleus. This excitation, however, must correspond to particular "states" that are unique for each nuclide; of particular importance to us is the energy of the first excited state -- if there is not enough energy to excite this state, then the scattering collision must be elastic. The amount of energy that a particle must have to excite a given excited is always slightly greater than the energy of the state itself, because of the energy that must go into the recoil of the target. This minimum energy is called the threshold energy. Allowing for the negative Q value, the threshold energy is given by: Et = -(A+1)/A*Q(note: Et is possitive).
| n + 70Ge |
| Reaction Products | Q-Value (MeV) | Threshold (MeV) |
| 71Ge + γ | 7.41596 | 0.00000 |
| 67Zn + α | 2.96327 | 0.00000 |
| 70Ge + n | 0.00000 | 0.00000 |
| 70Ga + p | -0.87327 | 0.88587 |
| 66Zn + n + α | -4.08896 | 4.14795 |
| 69Ga + n + p | -8.52843 | 8.65145 |
| 69Ge + 2 n | -11.53809 | 11.70453 |
| 68Ge + 3 n | -19.72617 | 20.01073 |
|
| n + 72Ge |
| Reaction Products | Q-Value (MeV) | Threshold (MeV) |
| 73Ge + γ | 6.78296 | 0.00000 |
| 69Zn + α | 1.47580 | 0.00000 |
| 72Ge + n | 0.00000 | 0.00000 |
| 72Ga + p | -3.21672 | 3.26184 |
| 68Zn + n + α | -5.00648 | 5.07669 |
| 71Ga + n + p | -9.73778 | 9.87435 |
| 71Ge + 2 n | -10.75207 | 10.90286 |
| 70Ge + 3 n | -18.16803 | 18.42282 |
|
| n + 73Ge |
| Reaction Products | Q-Value (MeV) | Threshold (MeV) |
| 74Ge
+ γ
| 10.19627 | 0.00000 |
| 70Zn
+ α
| 3.90873 | 0.00000 |
| 73Ge
+ n
| 0.00000 | 0.00000 |
| 73Ga
+ p
| -0.81095 | 0.82216 |
| 69Zn
+ n
+ α
| -5.30716 | 5.38056 |
| 72Ge
+ 2 n
| -6.78296 | 6.87678 |
| 72Ga
+ n
+ p
| -9.99968 | 10.13800 |
| 71Ge
+ 3 n
| -17.53503 | 17.77757 |
|
| n + 74Ge |
| Reaction Products | Q-Value (MeV) | Threshold (MeV) |
| 75Ge +γ | 6.50527 | 0.00000 |
| 74Ge + n | 0.00000 | 0.00000 |
| 71Zn +α | -0.45393 | 0.46012 |
| 74Ga + p | -4.58568 | 4.64825 |
| 70Zn + n +α | -6.28754 | 6.37333 |
| 73Ge
+ 2 n
| -10.19627 | 10.33540 |
| 73Ga
+ n
+ p
| -11.00722 | 11.15742 |
| 72Ge
+ 3 n
| -16.97923 | 17.21092 |
|
| n + 76Ge |
| Reaction Products | Q-Value (MeV) | Threshold (MeV) |
| 77Ge
+ γ
| 6.07262 | 0.00000 |
| 76Ge
+ n
| 0.00000 | 0.00000 |
| 73Zn
+ α
| -2.15651 | 2.18516 |
| 76Ga
+ p
| -6.22770 | 6.31044 |
| 72Zn
+ n
+ α
| -7.50946 | 7.60922 |
| 75Ge
+ 2 n
| -9.42839 | 9.55365 |
| 75Ga
+ n
+ p
| -12.03776 | 12.19769 |
| 74Ge
+ 3 n
| -15.93366 | 16.14535 |
|