A02 T. Kishimoto Osaka University CANDLES Collaboration Osaka - - PowerPoint PPT Presentation

a02
SMART_READER_LITE
LIVE PREVIEW

A02 T. Kishimoto Osaka University CANDLES Collaboration Osaka - - PowerPoint PPT Presentation

A02 T. Kishimoto Osaka University CANDLES Collaboration Osaka University, Graduate school of science Masoumeh Shokati Temuge Batpurev Ken Lee Keong Bui Tuan Khai


slide-1
SLIDE 1
  • T. Kishimoto

Osaka University

CANDLES Collaboration

Osaka University, Graduate school of science 吉田斉、Masoumeh Shokati、李暁龍、Temuge Batpurev、Ken Lee Keong、芥川一樹、Bui Tuan Khai、 佐藤勇吾、水越彗太、山本康平、宮本幸一郎 Osaka University, RCNP 梅原さおり、能町正治、岸本忠史、竹本康浩、松岡健次、瀧平勇吉、鉄野高之介 Fukui University 玉川洋一、小川泉、中島恭平、戸澤理詞、清水慧悟、清水健生、森勇太、池山佑太、小沢健太、松岡耕平 Tokushima University 伏見賢一 Osaka Sangyo University 硲隆太、中谷伸雄、Noithong Pannipa、田坪博貴 Tsukuba University 飯田崇史 Saga University 大隅秀晃 The Wakasa wan Energy Research Center 鈴木耕拓

Sendai2019/3/7 1

A02

slide-2
SLIDE 2
  • Highest Q value

– 4.27 MeV, (150Nd: 3.3 MeV) – Least BG(g: 2.6 MeV, b: 3.3 MeV) – Large phase space factor

  • Small natural abundance:

– 0.187% – Separated isotope → expensive

  • Next generation

– <mn>~ T-1/2 ~ M-1/2 (no BG) M: mass ~ M-1/4 (BG limited)

– Enrichment: mass+S/N: 500 times – High resolution: bolometer(crystal)

  • Beyond inverted hierarchy

– 48Ca + enrichment + bolometer

2

Why 48Ca

Natural abundance (%) Q value (MeV)

Nuclear matrix element → neutrino mass

__ CANDLES

Sendai2019/3/7

slide-3
SLIDE 3

Super Kamiokande KamLAND

CANDLES

XMASS GDZOOKS!

CANDLES III @ Kamioka

3

Lab D

Kamioka Lab. Map

4m 3m

CANDLES III

CANDLES III

Site: Kamioka U.G.L. ~1000 m Size: 3mΦ × 4mh (water tank) Liquid scintillator

Reservoir tank Purification system(liq-liq)

Sendai2019/3/7

slide-4
SLIDE 4

CANDLES III(UG)

CaF2 scintillator (CaF2(pure)) 305 kg (96 × 3.2kg) t ~ 1msec liquid scintillator (LS) 4π active veto 2m3 t ~ a few 10nsec PMT’s 13inch PMT × 48 20inch PMT × 14 light pipe light collection:energy resolution

CANDLES at Kamioka underground laboratory

CANDLES III

3m 4m

CaF2 Liquid scintillator

water

PMTs light pipe

Veto Pulse shape difference CaF2(pure) : ~1msec Liquid scintillator : a few 10 nsec

Sendai2019/3/7 4

slide-5
SLIDE 5

CANDLES III(UG)

CANDLES at Kamioka underground laboratory

CANDLES III

CaF2 scintillator (CaF2(pure)) 305 kg (96 × 3.2kg) t ~ 1msec liquid scintillator (LS) 4π active veto 2m3 t ~ a few 10nsec PMT’s 13inch PMT × 48 20inch PMT × 14 light pipe light collection:energy resolution

Veto Pulse shape difference CaF2(pure) : ~1msec Liquid scintillator : a few 10 nsec

CaF2 (305kg) Liquid scintillator tank(2m3) PMT Light pipe

Sendai2019/3/7 5

slide-6
SLIDE 6

4p active veto by Liquid scintillator (LS)

  • Rejection of external g–ray background
  • Pulse shape information by 500 MHz Flash ADC.
  • Distinguish event type by offline pulse shape analysis

taking advantage of different decay time.

6

CaF2 (1msec)

Liquid Scintillator (~10nsec) CaF2+ Liquid Scintillator Ca CaF2 Liquid Scintillator CaF CaF2 CaF CaF2 Liquid Scintillator

b-ray g-ray

bb signal !? External g BG External g BG

Sendai2019/3/7

slide-7
SLIDE 7

Internal backgrounds and reduction

  • External BGs were reduced by

LS active shild.

  • Remaining BGs originate from

internal radioactivity of Th chain (208Tl and 212Bi-212Po).

  • 2nbb is not serious BG in

current sensitivity. (it will be major BG after 48Ca enrichment)

  • We reject remaining BGs by

analysis.

b+g decay Q=5MeV b+a decay Qvis=5.3MeV T1/2=300ns

Sendai2019/3/7 7

Double pulse rejection Preceding a rejection

slide-8
SLIDE 8

8

Reduction

100MHz FADC (old)

DT > 30ns(3ch) ; ~5%

500MHz FADC . . . DT > 10ns ; ~2%

Rejection of Double Pulse

Prompt Delayed

Typical Pulse Shape

900ns 50ns

212Bi

212Po

T1/2 = 0.299msec

64%

Qa = 8.95MeV

Qb = 3.27MeV Qb = 2.25MeV

Qa = 7.83MeV

b a

208Pb

BG in Qbb region: Sum E Eα(1/3 Quench) +Eβ≃5.3 MeV

Sendai2019/3/7

slide-9
SLIDE 9

208Tl event cut

  • 1. Find parent 212Bi a-decay candidate by pulse shape analysis.
  • 2. Apply 12min veto from 212Bi candidate in the same crystal.

208Tl

T1/2 =3.05min

Qb = 5.00MeV

Qa = 6.21MeV

208Pb

stable

a b a-b崩壊

212Bi

T1/2 =60.6m

232Th

Th系列 Eneryg spectrum of prompt events

212Bi candidate Accidental event 212Bi candidate-Accidental

Delayed : 3.1-5MeV Time gate : ~180sec T½ = 178±55sec

Sendai2019/3/7 9

slide-10
SLIDE 10

External backgrounds

  • - Neutron source run --
  • To confirm our assumption that high E gamma ray BG’s

are from (n, g) reactions, 252Cf neutron source was set on the detector and data were taken.

2m

  • Spectra for neutron source run and physics run are consistent.
  • MC simulation of (n,g) can well reproduce the BG spectrum.

We identified main BG as (n,g) !!

Sendai2019/3/7 10

Neutron flux@Kamioka see K. Mizukoshi’s Poster

slide-11
SLIDE 11

Shield for (n,g) background reduction

CANDLES tank CANDLES shield overview n g Pb shield (7-12cm)

Reduce g-ray from surrounding rock Effect of Pb (n,g) is one order smaller than that of stainless tank

Boron sheet (4-5mm)

Reduce n captured by stainless tank

  • (n,g) BGs in CANDLES is expected

to become 1/80 by MC.

  • Expected number of backgrounds

after shield installation:

Rock : 0.34±0.14 e .14 even vent/year ear Tank : 0.4±0.2 .2 even vent/year ear

Sendai2019/3/7 11

slide-12
SLIDE 12

Pb shield construction

12

  • Pb shield construction was started from

March 2015.

  • All the collaborators worked very hard!

Bottom Pb shield Side Pb shield Top Pb shield

Sendai2019/3/7

slide-13
SLIDE 13

Position reconstruction and crystal selection

  • Position of each event is reconstructed by weighted mean of
  • bserved charge in each PMT.
  • Crystal separation is ~7s peak to peak.
  • Crystal selection criteria is within 3s from the peak.
  • 27 clean crystals (Th contamination < 10 mBq/kg) out of 96

crystals are selected and the results are compared to all crystals.

X axis [mm]

±3s Th activity [mBq/kg]

Sendai2019/3/7 13

slide-14
SLIDE 14

Energy Spectra & Event Selection

14

LiveTime : 131 days 95 crystals 27 crystals (232Th <10uBq/kg)

  • Exp. Data

212BiPo Cut

LS Cut

208Tl Cut

# e eve vent 95 95 crys ystal als 27 c crys ystal als

β

4-5MeV 5.5-6.5MeV Qβ

β

4-5MeV 5.5-6.5MeV

LS Cut 115 257 8 12 23 1

208Tl Cut

19 49 6 3 6 1 10 34 6 2 1

☑ No event in high purity crystals is confirmed.

β

β

Sendai2019/3/7

slide-15
SLIDE 15

Results

15

PANIC2017

95 C CaF aF2 27 C CaF aF2 Livetime 131 0ν β βeff. 0.39 ± 0.06 Event in ROI 10 Expected BG ~11 ~1.2 >3.8x1022 > 6.2x1022

22

Sensitivity (yr) 6.2x1022 3.6x1022

22

  • Exp. Data and BG MC

In 27 CaF2

β

χ2β<1.5, -3σ<SI<1σ

  • 2σ<position cut<2σ

Pileup cut > 20ns

208Tl cut

  • 1σ<0ν

β βwindow<2σ

CANDLES is now giving the best lifetime limit! ・further measurement ・developments for future

Sendai2019/3/7

slide-16
SLIDE 16

Statistics : 504 days

The obtained spectra as expected from BG estimation We have ~ 300 days more statistics (not yet finished analysis) BG from (n,γ) is reduced by ~ 100 with shield installation.

0νββ analysis

CaF2 Crystal x 21

228Th contents within crystal < 10 μBq/kg

All BG cuts are applied, but cut condition is not optimized yet.

LS veto & β-events cut

212Bi-Po sequential decay, 208Tl veto after 212Po-decay (18 min.)

Very Preliminary 21 crystals (228Th < 10 μBq/kg)

Livetime 504 days

Very Preliminary 93 crystals

Livetime 504 days

BG from (n,γ) Sendai2019/3/7

16

Replace crystals

slide-17
SLIDE 17

48Ca enrichment

  • Natural abundance of 48Ca is 0.187%.
  • 48Ca has a room of 500 times improvement (S & S/N)

by enrichment

  • Commercial 48Ca  too expensive (M$/10g but kg-

ton)

  • Enrichment is crucial for large volume 48Ca DBD

search.

  • Challenges in CANDLES:

– Crown ether resin + chromatography

  • 1.3 times

– Crown ether + micro reactor – Laser separation – Multi-channel counter current electrophoresis (MCCCE)

Sendai2019/3/7 17

slide-18
SLIDE 18

Laser Isotope Separation of 48Ca

Principle

Absorption spectrum of Ca isotope shift:a few hundred MHz laser line width:<1 MHz

Experiment Ionization Method

Optimization of various parameters

  • Excitation laser

power density

  • Ionization laser

wavelength power

Deflection Method

We confirmed the enrichment of 48Ca in the deflected atomic beam Optimization of laser power density

hν1 e.s. g.s hν2 i.p.

photon momentum

48Ca

emission (random) photon

48Ca

absorption TOF meas. of ion spatial distribution

  • meas. by TOF

ionization deflection

Two lasers ・selective excitation (CW, diode laser) ・ionization (pulse, dye laser) Single laser ・deflection (CW, diode laser) For TOF measurement ・ionization (pulse, YAG laser)

  • 2.5
  • 2
  • 1.5
  • 1
  • 0.5

0.5 Caイオン信号量[mV] 質量数 diode+dye dye

achieved 90% (nat. ×480)

40Ca 48Ca

40Ca 44Ca

Time[μs] ion signal[V]

40Ca 44Ca

Time[μs] ion signal[V]

48Ca 48Ca

same as 44Ca (nat.×11) TOF spectra at 4.5mm from the center of original atomic beam

ionization deflection

High concentration Small duty factor

pulse(10nsec, 10Hz)

  • Mid. Concentration

Continuous

  • peration

For the future mass production…

Development of collection system of 48Ca We continue R&D of Increase atomic beam/laser power

Deflection Method

Optimization of various parameters

Position of the collection plate will be adjusted to the

  • ptimal position

More details... see K. Matsuoka’s Poster

Sendai2019/3/7 18

slide-19
SLIDE 19

Multi-channel counter current electrophoresis

1 2 3 4 145 155 165 175 185

R(MCCCE)

V (Applied voltage)

) ( 48 / 43 ) ( 48 / 43 ) ( natural Ca Ca MCCCE Ca Ca MCCCE R 

Enrichment (48/43): 3.08 (48/40): 6

  • Separation using difference of migration

speed between 40Ca / 48Ca.

  • High power + effective heat removal

– Migration path: thermal conductor and insulator (BN)

  • Pulsed flow to get uniform flow speed

19 Sendai2019/3/7

Counter current Electric field

40Ca 48Ca

BN plate 10 mm thick 0.8mmΦ, every 4 mm

Ion exchange membrane

BN; Insulator but high thermal conductivity

inlet

  • utlet
slide-20
SLIDE 20

Sendai2019/3/7 20

Tabletop instrument

slide-21
SLIDE 21

history

  • 2012: got MCCCE idea
  • 2015: 10mm BN ~3 48Ca/43Ca, (6 48Ca/40Ca) PTEP

– then faced difficulty

  • 2017 year end

– After 2 years struggle, results become reproducible

  • 2018 February: ~10 times
  • 2018 April: modification to give uniform T and E

– ~a few 10’s times – May: ~ 100 times

  • Condition

– BN 20 mm

Sendai2019/3/7 21

slide-22
SLIDE 22

Highly enriched samples

Sendai2019/3/7 22

48Ca/43Ca (Enrichment/natural)

Obtained samples

TK, T. Ohata, K. Matsuoka

slide-23
SLIDE 23

48Ca/40Ca ratio

  • 48Ca/43Ca is so high then 48Ca/40Ca?

– We usually measure 48Ca/43Ca, since no interference – Similar nat. ab. 48Ca: 0.187%, 43Ca: 0.135%

  • 40Ar forbids 40Ca measurement in ICP-MS

– Reaction(collision)-cell ICP-MS + reaction-gas (H2, He, NH3)

23

Ar+ Ca+ H2 H+

40Ar+ 40Ar

Ar+ + H H+ + Ar

Reduce Ar+40Ca

40Ca

+

Sendai2019/3/7

slide-24
SLIDE 24

Enrichment

  • Migration distance

– μ:mobility difference – Separation

  • Diffusion: deteriorate separation

– Diffusion constant:D

  • Enrichment
  • Yield ~5% (concentration)

– Migration speed difference ~ 5% – Long Migration distance ~ 20/0.05 ~ 400 mm – Enrichment and yields are consistent

Sendai2019/3/7 24

Dt 2  s

) ( ) (

48 40

Ca Ca m m m   D

Et l m 

  E t E 1 1   D s

increase of E t (ℓ)

m m D D  ~ % 5 ~ v v Y

Et m D  D

slide-25
SLIDE 25

Enrichment

  • Enrichment → Reduction of 40Ca

– 10mm BN 1/6 → σ~10mm – 20mm BN 1/50 → σ~10mm – Width

  • Hagen-Poiseuille flow → (send pause, I sec )
  • Every 2mm times 200 → s small (1/√200)

– 40mm BN 1/300 → σ~14mm

  • 99.7% enrichment is possible
  • Practical goal is set to 80% enough for DBD

Sendai2019/3/7 25

Next step

v0 2v0

v0 2v0

thickness

BN  s

mm 115 ~ 12 400  s mm 8 ~ 200 12 2  s

10mm 20mm 40mm

80%

slide-26
SLIDE 26

Production of enriched 48Ca

  • Current system

– 16(10)% (48Ca/40Ca) – 12 cm2, 0.01 N → 0.1mg/day

  • Next system

– 80% or more – 1.2m2, 0.03 N → 0.3g/day → 100g/year – Tons; require plant → further needs

  • Our field
  • Other fields (beam, medical use,.. )

– CANDLES works for 80%

Sendai2019/3/7 26

plant

slide-27
SLIDE 27

本研究会では、不安定核物理、加速 器、二重ベータ崩壊、放射化学、医 学利用、異分野連携、新同位体濃縮 法といった視点からの講演を通して、 濃縮同位体を用いた研究と今後の発 展を概観することを目指しています。

“Workshop for Isotope Enrichment and Basic Science”

  • Mar. 21 Osaka

Sendai2019/3/7 27

slide-28
SLIDE 28
  • Sei Yoshida, ….
  • Collaborative research with Korean colleague

Yong-Hamb Kim (IBS & KRISS) Minkyu Lee (KRISS) Inwook Kim Do-Hyoung Kwon Hyejin Lee Hye-Lim Kim

Sendai2019/3/7

28

slide-29
SLIDE 29

Tail of 2νββ spectrum Improving energy resolution

48CaXX internal radioactivities

Th-chain(β-α sequential decays)  Bolometer Th-chain(208Tl)  Segmentation, Multi-crystal Environmental neutrons Improving resolution +Multi-crystal

But... new BG candidate

Q value of 48Ca : 4267.98(32) keV @ arXiv:1308.3815 Q-value of 238U (α-decay) : 4270 keV

Impossible to avoid  required particle ID

Scintillator  Bolometer Possible to further reduce the BG by developing Bolometer  Developing CaF2 Scintillating Bolometer

Sendai2019/3/7

29

slide-30
SLIDE 30

First Challenge using CaF2(pure) and MMC

CaF2 crystal Light Detector Crystal: CaF2(pure)

Volume: 300g (5cmφ×5cm) Emission peak : 280nm Light output: 25,000 photons/MeV β/γ α’s (226Ra, 222Rn, 218Po) β-α (214Bi-Po) μ

Problem

UV scintillation of CaF2 is absorbed on Au-deposit for heat signal. There is position dependence of scintillation absorption.  make worse E-resolution.

Sendai2019/3/7

30

slide-31
SLIDE 31

New trial to overcome UV absorption

CaF2(Eu) +Ag-deposit instead of CaF2(pure) + Ag-deposit

Au absorption spectrum Ag absorption spectrum CaF2(pure) emission CaF2(Eu) emission

222Rn 226Ra 218Po

β/γ and μ β-α (214Bi-Po)  Improved light signal properties.  In the heat channel, peaks of α’s are widely spread. (due to position dependence)  Due to doping Eu ? see Xiaolong Li’s Poster

Sendai2019/3/7

31

slide-32
SLIDE 32

Evaluate energy resolution w/o position dependence

We use contaminated CaF2 crystal for R&D, 226Ra ; ~ 30 mBq within crystal Delayed coincidence (222Rn  218Po  214Pb) Apply energy , PID parameter, Δtime cuts

3.10 min.

0 < Time difference < 3min Correlated events (α-α events ) from same position

Evaluated energy resolution without position dependence, < 0.2% (σ) @ Qββ

0.18 % (σ) @ 5.85 MeV

using maximum pulse height

Sendai2019/3/7

32

slide-33
SLIDE 33

Improving E-resolution of CaF2(pure) scintillating bolometer

Radio-pure CaF2(pure) crystal had been developed. Doping Eu may affect phonon propagation in CaF2 crystal.

New trial in the next step

CaF2(pure) crystal with smaller but thicker Au-deposit phonon collector.

Smaller  reducing scintillation absorption effect Thicker  increasing the strong electron-phonon interaction.

Sendai2019/3/7

33

slide-34
SLIDE 34

CANDLES project

  • CANDLES:

– CANDLES III(UG): current detector

  • First result gives the best limit
  • Measurement and analysis underway
  • CaF2 crystals with low BG (will be replaced)
  • Future prospect

– Enrichment 48Ca

  • MCCCE works for future tons

– Bolometer

  • CaF2 Scintillating Bolometer

34

CANDLES to normal hierarchy

Still lot to do Promising

Sendai2019/3/7

slide-35
SLIDE 35

35 Sendai2019/3/7

slide-36
SLIDE 36
  • T. Kishimoto

Osaka University

CANDLES Collaboration

Osaka University, Graduate school of science 岸本忠史、吉田斉、鈴木耕拓、角畑秀一、Wang Wei、Chan Wei Min、 Van Trang、 石川貴志、田中大樹、田中美穂、土井原正明、前田剛、太畑貴綺、鉄野高之介 Osaka University, RCNP 能町正治、味村周平、梅原さおり、中島恭平、飯田崇史、松岡健次 Fukui University 玉川洋一、小川泉、中島恭平、川村篤史、富田翔悟、藤田剛志、原田知優、坂本康介、 吉澤真敦、犬飼祐司 Tokushima University 伏見賢一 Osaka Sangyo University 硲隆太、中谷伸雄 Tsukuba University 飯田崇史

Sendai2019/3/7 38

slide-37
SLIDE 37

Sendai2019/3/7 39

0.5 1 1.5 2 2.5 3 3.5 4 4.5 0.100 1.000 10.000 100.000 1000.000

48/40 Ca 48Caの濃度[ppm]

5% 10% 20% 35% 70%