Measurement of ambient neutrons in an underground laboratory at - - PowerPoint PPT Presentation
Measurement of ambient neutrons in an underground laboratory at - - PowerPoint PPT Presentation
Measurement of ambient neutrons in an underground laboratory at Kamioka Observatory Keita Mizukoshi Kobe University TAUP2019 at Toyama International Conference Center 9 Sep. 2019 Introduction Neutron is the one of the most serious
Keita Mizukoshi Kobe Univ.
Introduction
- Neutron is the one of the most serious backgrounds
(BG) for experiments in underground.
- Direct Dark matter search
- Neutrino-less double beta decay search
- To evaluate and shield neutron BG,
it is very important to evaluate ambient neutron flux
- Such neutron BG has not measured systematically.
- Low rate in underground i→ Required high efficiency
- Unknown generated points → Energy unknown
- Our goal is quantitative neutron flux in the underground
2
Keita Mizukoshi Kobe Univ.
Detector (He-3 proportional counter)
- We used a 3He
proportional counter.
- The energy of the
exothermal reaction in the neutron capture can be obtained.
- This detector is
sensitive to thermal neutrons (~0.025 eV), and cannot measure an initial neutron energy.
3
3He proportional counter
DAQ PC
3He + n → 3H + p + 0.76 MeV
3He - 10 atm
SUS
thermal fast
“Setup A”
380mm φ52mm
- K. Mizukoshi et al., PTEP 123C01
~5 k barn at thermal
Keita Mizukoshi Kobe Univ.
Setup for fast neutron
thermal fast “Setup A”
Polyethylene
3He
Boron sheet
thermal fast “Setup B”
510mm
t5mm
3He
Efficiency estimated by Geant4
9 −
10
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 Generated neutron energy (MeV) 20 40 60 80 100 120 140 160 )
2
Counts/(neutron/cm
Setup A Setup B Setup B w/o B-sheet
Thermal neutron ~0.025eV
- To measure high energy neutron,
we used a moderator (polyethylene).
- Boron sheet captured thermal
neutrons and reduce its effect.
4
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
Results
- Full energy peak is 0.76 MeV.
If 3H or p escapes, continuum region will be made in a low energy (Wall effect).
- Low energy region below 0.3 MeV is
dominated by electric noise for ambient neutron measurement.
- We counted events up to 0.85 MeV
and down to 0.5 MeV, then the number
- f total events was estimated by a
clear spectrum observed using 252Cf.
5
Setup
A B
Count rate (×10-3cps)
1.295 ± 0.034 0.446 ± 0.018
Live time (day)
14.03 19.27
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Measured energy (MeV) 200 400 600 800 1000 1200 Counts/bins 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Measured energy (MeV) 20 40 60 80 100 Counts/bins
Spectrum of Source (252Cf) Measured spectrum in setup B
3He + n → 3H + p + 0.76 MeV
Count rate in each setup
- The count rate of Setup A
(RA) and B (RB) involves a detection of thermal and fast neutron, respectively.
Electric noise peak Wall effect
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
Simulation
- To convert from the count rates (RA, RB) to
ambient neutron flux, the spectral shape was
- required. The shape cannot measured by He-3
detector thus estimated by simulation.
- We considered the source of the neutrons
made from (α,n) reaction of U/Th series decay.
- Neutron induced by cosmic muon is negligible.
- We picked three types of rocks as samples,
they had much different abundance of chemical compositions.
- The difference affects much the yield of
neutrons.
6
1 2 3 4 5 6 7 Neutron energy (MeV) 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2
6 −
10 × Number of nentrons /chain decay/0.1 MeV
, n) α JR-1 ( , n) α JA-3 ( , n) α Sample 1 ( Sample 1 + 3% of hydrogen , n) α Sample 2 ( , n) α Sample 3 ( U fission)
238
Sample 1 ( Sample 1 (Th series)
Bin width 0.1 MeV
Generated neutron in vary rocks
(wt. %) sample1 sample2 sample3
O 40.5 37.9 35.6 Ca 28.0 24.3 29.7 Si 16.6 15.6 12.0 Fe 7.6 16.6 13.5 Al 5.2 0.3 0.1 Mn 0.8 3.5 2.9
Main components in each samples
JR-1 and JA-3 are geometrical reference database
Much difference!
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
Data driven analysis
- We cannot investigate the all wall rock
components in details.
- Especially amount of water contents in the
rock and chemical composition including Hydrogen will much affect thermalization of fast neutrons.
- Thus, thermalization in the rock was
unknown.
- We regarded the percentage of hydrogen
(%of h. e.) in simulation as a thermalization parameter.
- %of h.e. was derived by the experimental
result (the ratio between setups A and B) in each rock component.
- The most likely spectra (made from
experimental data) in each sample are almost same.
- This is not affected by uncertainty of
Simulation. 7
1 2 3 4 5 6 7 % of h. e. 0.5 1 1.5 2 2.5 3 3.5 4
B
/R
A
Ratio of count rates R
Experimental ratio Error band Sample 1 Sample 2 Sample 3 JR-1 JA-3 9 −
10
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 Neutron energy (MeV)
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 /s
2
Counts /MeV /cm
Sample 1 with 3% of h.e. JR-1 with 1% of h.e. JA-3 with 1% of h.e. 1/E
The most likely spectrum
Experimental ratio v.s. parameter
Thermalization parameter
- btained by measurement
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
Obtained spectrum
- We obtained the most likely
spectrum of the ambient neutron.
- We compared the fluxes
(the previous study fluxes in other underground laboratories).
- They are the same order of
magnitude.
- It is difficult to compare
the result simply because there are many difference in these measurement (e.g., detector, assumption
- f spectral shape, and
definition of flux) 8
9 −
10
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 Neutron energy (MeV)
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 /s
2
Counts /MeV /cm
The most likely spectrum 2% of hydrogen 4% of hydrogen 1/E
Flux (×10-6 cm-2 s-1) Thermal Non-thermal
Kamioka (This result, Mizukoshi) 7.9 ± 0.23 +0.7
- 0.7
15.6 ± 0.5 +1.2
- 1.4
Kamioka (Minamino 2004) 8.26 ± 0.58 11.5 ± 1.2 Gran Sasso (A. Lindi 1988)※ 13.3 ± 1.5 10.2 ± 1.1 LSM (K. Eitel 2012)※ 14.3 ± 1.3 4.2 ± 2.8
The most likely spectrum Neutron fluxes in previous researches
※They used the different definition of flux. We adjusted the same definition of us.
Thermal neutron ~0.025eV
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
New interest
- In the previous research, rough
spectral shape was assumed (e.g., Boltzmann distribution and 1/E).
- The most likely spectrum
suggests the excess in a few MeV.
- The excess is interesting
for direct dark matter search.
- The excess should be
confirmed by a liquid scintillator which has a sensitivity for the neutron.
- Even such basic information
has not confirmed…
9
9 −
10
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 Neutron energy (MeV)
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 /s
2
Counts /MeV /cm
The most likely spectrum 2% of hydrogen 4% of hydrogen 1/E
The most likely spectrum
Excess
- Since the cross section of
high energy neutrons is small, it continues to be a high energy neutron.
- Once it lose energy, the
cross section increases. it continues to lose energy.
- Therefore, the excess will
remain at several MeV.
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
Summary
- We evaluated an ambient neutron spectrum and
- btained the flux (23.5 ± 0.7stat. sys.×10-6 cm-2 s-1)
at the Kamioka Observatory.
- using 3He proportional counter and moderator effectively
- with data-driven analysis and simulation
- considering systematic errors
- Spectral excess around a few MeV was
- suggested. It should be confirmed by a sensitive
detector for non-thermal neutron.
- We are preparing a low BG liquid scintillator.
10
+1.9
- 2.1
Backup slides
Keita Mizukoshi Kobe Univ.
Main backgrounds of our experiment
- Rate events search experiments are
placed in underground to reduce BGs ~ Cosmic muon
- For Others
- Remaining high-energy muon
← active veto by scintillator
- Ambient gamma
← shield / self-shielding / PSD
- Alpha from U/Th chain
← very careful washing
- Neutron
← It’s difficult to reduce
- Shield with materials which
have large cross-section for neutron
- Sometimes neutron makes
- ther BGs (gamma in detector)
12
Fiducial Volume (FV) for events
cosmic muon External gamma
General ways to reduce BGs in underground experiment
Muon Veto counter Shield (Pb)
(Conceptional fig)
Detector neutron
Keita Mizukoshi Kobe Univ.
How much neutron is in underground
- I worked for a neutrino-less
double beta decay experiment (CANDLES 3+ Experiment) at Kamioka Observatory.
- This experiment reduces
ambient neutron with Boron shield (~5000 barn for thermal neutron).
- Ambient neutron (flux,
spectrum) was not understood well.
- I would like to show how to
measure ambient neutron to demonstrate how difficult to handle neutron as a rare BG.
13
CANDLES Experiment
for neutrinoless double beta decay using 48CaF2 Scintillator
Photo CANDLES Collaboration
Keita Mizukoshi Kobe Univ.
Main neutron detectors for underground
14
He-3 Counter Liquid Scintillator Gd-Water Čerenkov
Čerenko
Emulsion Thermal Neutron Very good No Very good No Fast Neutron No sensitivity† Good Good Good Energy Sensitivity No Possible No Possible n/γ rejection Good Sometime Good No Good Internal background Low Sometime bad Low Low Handling Analysis Easy Difficult Difficult Very difficult Detection Efficiency Very good so-so so-so Very good
†Generally.
Keita Mizukoshi Kobe Univ.
Neutron source (Not U/Th)
- We can consider the ambient
neutron made from cosmic muon.
- It can make high energy
neutron (>10 MeV), the number of neutrons by muon is 100 times less than the
- nes by U/Th series.
- In this research, we ignored
the contribution of muon.
15
Generated neutrons from cosmic muon Spectrum of each source for sample1
100 200 300 400 500 600 700 800 900 1000 Neutron energy(MeV)
16 −
10
15 −
10
14 −
10
13 −
10
12 −
10
11 −
10 /s
3
Number of produced neutrons /MeV /cm
Sample 1 Sample 2 Sample 3
Bin width 10 MeV 9 −
10
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10
2
10 Neutron energy (MeV)
14 −
10
13 −
10
12 −
10
11 −
10
10 −
10
9 −
10
8 −
10
7 −
10
6 −
10
5 −
10
4 −
10
3 −
10
2 −
10
1 −
10 1 10 /s
2
Counts /MeV /cm
Summed spectrum w/o muon , n) α U ( , n) α Th ( U fission
238Muon
Keita Mizukoshi Kobe Univ.
Simulation uncertainty for neutron
- Geant4 - Monte-Carlo toolkit
widely used in Particle Physics
- PHITs - Boltzmann equation
solver used in Nuclear Physics
- These simulation tools produce
huge uncertainty.
- These spectrum simulate
neutrons made in rocks around underground
- laboratory. These neutron
source is Uranium and Thorium chain nuclei.
16
0.05 0.1 0.15 0.2 0.25 0.3
6 −
10 × Energy(MeV) 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 current /cm^2 /MeV /s 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 Neutron energy (MeV) 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
6 −
10 × Counts /s /cm^2 /MeV
Example: Differences of simulation tools for transported neutron energy spectrum PHITs Geant4 Thermal neutron Fast neutron PHITs Geant4
- K. Mizukoshi et al., PTEP 123C01
Keita Mizukoshi Kobe Univ.
All components of the rocks
17
(wt. %) sample1 sample2 sample3
sampleA sampleC
KamRock Si
16.6 15.6 12.0 29.1 27.8 18.5
Ti
0.2 0.0 0.0 0.5 0.5 0.1
Al
5.2 0.3 0.1 7.1 8.3 10.6
Fe
7.6 16.6 13.5 4.6 4.5 1
Mn
0.8 3.5 2.9 0.1 0.1
Mg
0.6 1.1 0.7 2.2 1.3 0.3
Ca
28.0 24.3 29.7 4.5 5.2 1.8
Na
0.0 0.2 0.0 2.4 2.6 3.9
K
0.0 0.1 0.0 1.2 1.5 2.1
P
0.2 0.0 0.0 0.1 0.1 0.1
S
0.0 0.1 1.2 0.0 0.0
Zn
0.0 0.1 4.3 0.0 0.0
Sr
0.1 0.0 0.0 0.0 0.0
Nb
0.0 0.0 0.0 0.0 0.0
Sn
0.1 0.0 0.0 0.0 0.0
Pb
0.0 0.0 0.0 0.0 0.0
O
40.5 37.9 35.6 48.3 48.1 60.7
Keita Mizukoshi Kobe Univ.
He-3 Cross section for neutron
- He cross section is much
large for thermal neutrons.
- Cross sections of the rock
components have the same trend.
- Since the cross section of
high energy neutrons is small, it continues to be a high energy neutron.
- Once it lose energy, the
cross section increases. it continues to lose energy.
- Therefore, the dip will
remain at several MeV.
18
10-7 10-6 10-5 10-4 10-3 10-2 10-1 100 101 102 103 104 105 10-2 10-1 100 101 102 103 104 105 106 107
Cross Section (barns) Neutron Energy (eV)
He-3
total elastic capture (n,p)
Keita Mizukoshi Kobe Univ.
Definition of Flux
- Two types of definition are
used.
- (1)Number of particle through
the sphere (radius r)/ the area of grate circle(πr2)
- Widely used in Nuclear physics
- We use that.
- (2)Number of particle through
the circle (radius r)/ the area (πr2)
- Widely used in Particle physics
- LSM and Gran Sasso would
use this definition.
19
R r
- Def. 1
- Def. 2