ACCOUNTING FOR DIFFERENT REACTOR DISTANCES Maria Veronica Prado - - PowerPoint PPT Presentation

accounting for different reactor
SMART_READER_LITE
LIVE PREVIEW

ACCOUNTING FOR DIFFERENT REACTOR DISTANCES Maria Veronica Prado - - PowerPoint PPT Presentation

APPEARANCE PROBABILITY MODEL OF ELECTRON ANTI-NEUTRINOS ACCOUNTING FOR DIFFERENT REACTOR DISTANCES Maria Veronica Prado Advisors: Dr. Glenn Horton-Smith Dr. Larry Weaver KamLAND Experiment 56 nuclear reactors and one detector Detector


slide-1
SLIDE 1

APPEARANCE PROBABILITY MODEL OF ELECTRON ANTI-NEUTRINOS ACCOUNTING FOR DIFFERENT REACTOR DISTANCES

Maria Veronica Prado Advisors: Dr. Glenn Horton-Smith

  • Dr. Larry Weaver
slide-2
SLIDE 2

KamLAND Experiment

 56 nuclear reactors and one detector  Detector is located on the island of Honshu, Japan  Each nuclear reactor contains Uranium 235 and 238 &

Plutonium 239 and 241

 Fission occurs:

57.1% from U 235 7.8% from U 238 29.5% from Pu 239 5.6% from Pu 241

slide-3
SLIDE 3

KamLAND Experiment

Source: http://kamland.lbl.gov/Pictures/kamland-ill.html

slide-4
SLIDE 4

KamLAND Experiment

n

moderator reactor with U235 U238 Pu239 Pu 241 two new elements after fission beta decay

n p + e + v

e

n p + e + v

e

The reactors:

slide-5
SLIDE 5

KamLAND Experiment

The detector:

Source: http://kamland.lbl.gov/Pictures/kamland-ill.html

slide-6
SLIDE 6

KamLAND Experiment

 The Liquid Scintillator inside the detector contains C9

H12 (pseudocumene) and C12 H26 (dodecane)

 Some of the anti-neutrinos coming from the reactors

collide with protons found in these molecules

 Inverse beta decay

slide-7
SLIDE 7

KamLAND Experiment

The detector:

Inverse beta decay

v + p n + e

e

+

Positron moves through the LS losing KE as it ionizes atoms Ionization of atoms +

e

UV/Gamma photon unstable atom

Process A

Gamma photons: 10 keV Optical photons: 1-3 eV

γ

slide-8
SLIDE 8

KamLAND Experiment

H H H

Fluorescence The UV photon hits one of the molecules and is absorbed

The detector:

γ γ

A visible photon is emitted from the molecule Optical photon

DETECTED

The visible photon is detected by all the photomultiplier tubes

slide-9
SLIDE 9

KamLAND Experiment

Meanwhile still in the detector…

e + e

+ _

The positron loses energy and comes to a stop where it annihilates with an electron

γ γ

Gamma photons Gamma photons Compton scatter or go through the Photoelectric Effect _

e

_

e

_

e

_

e

_

e

_

e

_

e

_

e

_

e

_

e

_

e

_

e

Go through Process A

DETECTED

Even though there are 3 separate signals in the PMTs, it is detected as only one

slide-10
SLIDE 10

KamLAND Experiment

Simultaneously in the detector…

v + p n + e

e

+

n + p H +

γ

2

The neutron bounces off from the atoms in the LS and moves slower & slower until it is absorbed Recoil (Deuteron)

Go through Process A

Or very rarely

n + p C + γ

13

The neutron interacts with hydrogen (H) from the LS The gamma photon compton scatters or goes through the photoelectric effect with the atoms in the LS. It produces a detected signal called the delayed coincidence.

slide-11
SLIDE 11

KamLAND Experiment

Source: kamland.lbl.gov/Pictures/picgallery.html

The detector:

slide-12
SLIDE 12

Number of Counts

Our main equation:

slide-13
SLIDE 13

A simple derivation

Flux of the anti-neutrino

slide-14
SLIDE 14

A simple derivation

where,

slide-15
SLIDE 15

Terms

The number of counts at each energy prompt The flux of anti-neutrinos expected at the detector The cross section of one proton that could interact with the anti- neutrinos coming into the detector Probability that an electron anti-neutrino will stay an electron anti-neutrino by the time it reaches the detector Probability of detecting a reaction from the reactions that have occurred (due to experimental error) The number of reactions

slide-16
SLIDE 16

Events graph with what was observed in the detector

KamLAND events graph

http://www.awa.tohoku.ac.jp/KamLAND/4th_result_data_release/4th_result_data_release.html

slide-17
SLIDE 17

KamLAND appearance probability graph

http://arxiv.org/pdf/1009.4771v2.pdf

Appearance probability from an average reactor length

slide-18
SLIDE 18

Our theoretical research

Make a change of variables:

Q Q P P N N

slide-19
SLIDE 19

Our theoretical research

For simplicity, we shall call this:

We need to minimize for N(Ep): where,

N(Ep) has a Poisson distribution because of the rare amount of interactions at the detector

what we want to find empirically

slide-20
SLIDE 20

Our theoretical research

By taking the derivative of and setting it equal to zero, we get

where,

transformation of the Q matrix Our observed values

slide-21
SLIDE 21

Small proof

slide-22
SLIDE 22

Why we want to do this

 Prove neutrino oscillations and KamLAND’s conclusions

empirically

 Gain knowledge about how neutrinos behave, which

could lead to a better understanding of dark matter

 Gain knowledge about neutrinos to be able to control

nuclear reactors efficiently by monitoring neutrinos that leave

slide-23
SLIDE 23

Forming the Q matrix

For ex:

  • Test for as many l’s as possible, binning them
  • If the lies between 1.8 MeV-10 MeV, then plug the values into

the Q equation

  • If the lies outside of that range, it does not contribute to the

detector, so we input zero for that matrix element

  • Obtain a different Q matrix for each reactor
  • Superpose all the Q matrices
slide-24
SLIDE 24

Forming the Q matrix

slide-25
SLIDE 25

Our ‘no oscillations’ graph

  • ur ‘no oscillations’ graph (without

taking into account certain small factors)

KamLAND’s ‘no oscillations’ graph

slide-26
SLIDE 26

Setting up the test

where,

appearance probability if there were no oscillations

C C Y Y

where,

slide-27
SLIDE 27

Setting up the test

Where R is a matrix containing the

  • rthonormal eigenvectors for each

eigenvalue, and D is a diagonal matrix containing all the eigenvalues

  • f C

Since,

Has the smallest eigenvalue element equal to zero

slide-28
SLIDE 28

The Binning

 The greater the counts per bin, the smaller

the relative error

 As a result, approaches a Gaussian

Why bin the Ep’s? Why bin the l’s?

 More functions than unknowns  A higher sum in each l column will provide

for a smaller error

Why find the eigenvalues of C?

 If product of eigenvalues is big, error is

small when inverting C

 If difference is big, magnifies error

slide-29
SLIDE 29

Why test it this way

 Accounting for bias by using N0 prime to

calculate V inverse instead of N1 prime:

  • This method gives each element in N1

prime their corresponding importance according to how many number of counts they each contribute and, therefore, how much data they contain

For example: while

slide-30
SLIDE 30

Why omit the smallest eigenvalue?

Contains background noise Contains inverse eigenvalues The smaller the eigenvalue, the more noise error it contributes

slide-31
SLIDE 31

Error

 Background noise in the data N1 from the experiment  Approximation of l values due to the l binning in Q0

Biggest error contributors:

 Create N’1true with a specific P(l)  Add randomized background noise to N’1true  Create 1000 different P(l)s, each using a different

randomized N’1observed

 Find the average P(l) and its standard deviation to obtain

different error bars for each P(l) entry

Producing reasonable error bars for our test of specific P(l)s :

slide-32
SLIDE 32

Omitting the smallest eigenvalue

16 Ep bins,11 l bins Smallest error bars so far, but not as small as were expected

slide-33
SLIDE 33

16 Ep bins,11 l bins Smaller error bars

Omitting vs not omitting smallest eigenvalue

slide-34
SLIDE 34

Testing KamLAND’s N

slide-35
SLIDE 35

Comparing Chi Squares

Chi square of the N between our estimate and the N observed:

6.68

Chi square of the P between our estimate and the closest straight line

  • f 0.44 without taking into account covariance:

9.65

The above chi square with covariance:

67.95

slide-36
SLIDE 36

The Covariance Matrix

slide-37
SLIDE 37

Conclusions

 Obtained appearance probabilities for

11 values of L/Eν without assuming an average L

 Appearance probability cannot be

constant

 Predicted N matched KamLAND’s

  • bserved N
slide-38
SLIDE 38

Acknowledgements

 Dr. Horton-Smith & Dr. Weaver  Dr. Corwin  KSU HEP Department  Kansas State University  NSF