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:
20keV: 1 2 3 4 5 6 7 8 9
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:
200keV: 1 2 3 4 5 6 7 8 9
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:
2MeV: 1 2 3 4 5 6 7 8 9
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:
20MeV: 1 2 3 4 5 6 7 8 9
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:
200MeV: 1 2 3 4 5 6 7 8 9
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:
2000MeV: 1 2 3 4 5 6 7 8 9
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