06/30/10
Teppei Katori Massachusetts Institute of Technology CPT and Lorentz symmetry 2010 Bloomington, Indiana, USA, June 30, 10
Test for Lorentz violation with MiniBooNE low energy excess
1 Teppei Katori, MIT
Test for Lorentz violation with MiniBooNE low energy excess Teppei - - PowerPoint PPT Presentation
Test for Lorentz violation with MiniBooNE low energy excess Teppei Katori Massachusetts Institute of Technology CPT and Lorentz symmetry 2010 Bloomington, Indiana, USA, June 30, 10 06/30/10 Teppei Katori, MIT 1 Test for Lorentz violation
06/30/10
Teppei Katori Massachusetts Institute of Technology CPT and Lorentz symmetry 2010 Bloomington, Indiana, USA, June 30, 10
1 Teppei Katori, MIT
06/30/10
Teppei Katori Massachusetts Institute of Technology CPT and Lorentz symmetry 2010 Bloomington, Indiana, USA, June 30, 10
2 Teppei Katori, MIT
06/30/10
3 Teppei Katori, MIT
06/30/10
Neutrino oscillation is an interference experiment (cf. double slit experiment)
νµ
If 2 neutrino Hamiltonian eigenstates, ν1 and ν2, have different phase rotation, they cause quantum interference.
νµ
4 Teppei Katori, MIT
06/30/10
Neutrino oscillation is an interference experiment (cf. double slit experiment) If 2 neutrino Hamiltonian eigenstates, ν1 and ν2, have different phase rotation, they cause quantum interference. If ν1 and ν2, have different coupling with Lorentz violating field, interference fringe (oscillation pattern) depend on the sidereal motion.
νµ ν1 ν2 Uµ1 Ue1
*
ν2 ν1 νµ
5 Teppei Katori, MIT
06/30/10
Neutrino oscillation is an interference experiment (cf. double slit experiment) If 2 neutrino Hamiltonian eigenstates, ν1 and ν2, have different phase rotation, they cause quantum interference. If ν1 and ν2, have different coupling with Lorentz violating field, interference fringe (oscillation pattern) depend on the sidereal motion. The measured scale of neutrino eigenvalue difference is comparable the target scale of Lorentz violation (<10-19GeV).
νµ νe ν1 ν2 Uµ1 Ue1
*
νe
6 Teppei Katori, MIT
06/30/10 Booster
K+
target and horn detector dirt absorber
primary beam tertiary beam secondary beam
(protons) (mesons) (neutrinos)
π+
decay region FNAL Booster
MiniBooNE neutrino oscillation experiment at Fermilab is looking for νµ to νe oscillation Signature of νe event is the single isolated electron like events
νµ
→ νe
e (electron-like Cherenkov) p n W
νµ
7 Teppei Katori, MIT
06/30/10
Booster Target Hall
MiniBooNE extracts beam from the 8 GeV Booster
Booster
K+
target and horn detector dirt absorber
primary beam tertiary beam secondary beam
(protons) (mesons) (neutrinos)
π+
decay region FNAL Booster
MiniBooNE collaboration, PRD79(2009)072002
νµ
8 Teppei Katori, MIT
06/30/10
within a magnetic horn (2.5 kV, 174 kA) that increases the flux by × 6 8GeV protons are delivered to a 1.7 λ Be target Magnetic focusing horn
Booster primary beam tertiary beam secondary beam
(protons) (mesons) (neutrinos)
K+
π+
target and horn dirt absorber detector decay region FNAL Booster π+ π+ π- π-
MiniBooNE collaboration, PRD79(2009)072002
νµ
9 Teppei Katori, MIT
νµ
The decay of mesons make the neutrino beam. The neutrino beam is dominated by νµ (93.6%), of this, 96.7% is made by π+-decay
Booster primary beam tertiary beam secondary beam
(protons) (mesons) (neutrinos)
K+
π+
target and horn dirt absorber detector decay region FNAL Booster π+ π+ π- π- Predicted νµ-flux in MiniBooNE
06/30/10
MiniBooNE collaboration, PRD79(2009)072002
π+ →µ+ + νµ
νe?
10 Teppei Katori, MIT
06/30/10 Booster
K+
target and horn detector dirt absorber
primary beam tertiary beam secondary beam
(protons) (mesons) (neutrinos)
π+
decay region FNAL Booster
MiniBooNE detector is the spherical Cherenkov detector
MiniBooNE collaboration, NIM.A599(2009)28
νµ νe?
11 Teppei Katori, MIT
– Sharp, clear rings
– Scattered rings
MiniBooNE collaboration, NIM.A599(2009)28
µ ν
– Sharp, clear rings
– Scattered rings
MiniBooNE collaboration, NIM.A599(2009)28
ν e
06/30/10
14 Teppei Katori, MIT
06/30/10
Blind analysis MiniBooNE perform ~5 years blind analysis. νe candidate data is not used to tune MC. Background errors are constraint from measurements by MiniBooNE detector. e.g.
NCπo measurement of MiniBooNE provide precise estimation for this type of background.
πο
MiniBooNE collaboration PLB664(2008)41 PRL100(2008)032301
πο NCπo production
NCπo production with asymmetric decay
15 Teppei Katori, MIT
06/30/10
Blind analysis MiniBooNE perform ~5 years blind analysis. νe candidate data is not used to tune MC. Background errors are constraint from measurements by MiniBooNE detector. e.g.
νµ measurement of MiniBooNE provide precise prediction of νe from µ-decay
MiniBooNE collaboration PLB664(2008)41 PRL100(2008)032301
µ → e νµ νe π → µ νµ
Predicted neutrino flux
νµ νe muon neutrino
νe from µ-decay
16 Teppei Katori, MIT
06/30/10
MiniBooNE first oscillation result There is no νe candidate excess in the analysis region (where the LSND signal is expected from 1 sterile neutrino interpretation). However there is visible excess at low energy region, which is inconsistent with two massive neutrino oscillation hypothesis. Is this excess real physics?
MiniBooNE collaboration, PRL98(2007)231801
17 Teppei Katori, MIT
06/30/10
MiniBooNE collaboration, PRL102(2009)101802
New MiniBooNE oscillation result After ~1 year careful reanalysis, again no excess in oscillation candidate region, but low energy excess is confirmed. The energy dependence of νe low energy excess is not consistent with ~1/E, hence it is not consistent with two massive neutrino
Lorentz violating neutrino oscillation has various energy dependences, therefore it is interesting to test Lorentz violation!
18 Teppei Katori, MIT
06/30/10
19 Teppei Katori, MIT
06/30/10
The examples of model independent features that represent characteristic signals of Lorentz violation for neutrino oscillation (1) Spectral anomalies (2) L-E conflict (3) Sidereal variation Any signals cannot be mapped on Δm2- sin22θ plane (MS-diagram) could be Lorentz violation, since under the Lorentz violation, MS diagram is no longer useful way to classify neutrino oscillations LSND is the example of this class of signal.
Kostelecky and Mewes PRD69(2004)016005
20 Teppei Katori, MIT
06/30/10
The examples of model independent features that represent characteristic signals of Lorentz violation for neutrino oscillation (1) Spectral anomalies (2) L-E conflict (3) Sidereal variation Any signals do not have 1/E oscillatory dependence could be Lorentz violation. Lorentz violating neutrino oscillation can have various type of energy dependences. MiniBooNE signal falls into this class. usual term (3X3) additional terms (3X3) effective Hamiltonian of neutrino oscillation (direction averaged)
(heff)ab ~ 1 2E(m2)ab +(a)ab+(c)abE+
Kostelecky and Mewes PRD69(2004)016005
21 Teppei Katori, MIT
06/30/10
The examples of model independent features that represent characteristic signals of Lorentz violation for neutrino oscillation (1) Spectral anomalies (2) L-E conflict (3) Periodic variation sidereal variation of the neutrino
Lorentz violation This signal is the exclusive smoking gun of Lorentz violation. Test for Lorentz violation in MiniBooNE is to find sidereal time variation from low energy excess νe candidate data example of sidereal variation for LSND signal
Kostelecky and Mewes PRD69(2004)016005
22 Teppei Katori, MIT
06/30/10
Test for Lorentz violation in MiniBooNE follows LSND sidereal time analysis (1) fix the coordinate system (2) write down Lagrangian including Lorentz violating terms under the formalism (3) write down the observables using this Lagrangian
LANSCE to LSND detector is almost east to west χ = 54.1o, colatitude of detector θ = 99.0o, zenith angle of beam φ = 82.6o, azimuthal angle of beam
beam dump detector 30.8m
LSND collaboration, PRD72(2005)076004
23 Teppei Katori, MIT
06/30/10
Test for Lorentz violation in MiniBooNE follows LSND sidereal time analysis (1) fix the coordinate system (2) write down Lagrangian including Lorentz violating terms under the formalism (3) write down the observables using this Lagrangian
Booster neutrino beamline is almost south to north χ = 48.2o, colatitude of detector θ = 90.0o, zenith angle of beam φ = 180o, azimuthal angle of beam
Be target detector 541m
24 Teppei Katori, MIT
06/30/10
Test for Lorentz violation in MiniBooNE follows LSND sidereal time analysis (1) fix the coordinate system (2) write down Lagrangian including Lorentz violating terms under the formalism (3) write down the observables using this Lagrangian
Kostelecky and Mewes PRD69(2004)016005
SME parameters Modified Dirac Equation (MDE)
i(ΓAB
ν ∂ν − MAB)νB = 0
ΓAB
ν = γνδAB + cAB µν γµ + dAB µν γµγ5 + eAB ν
+ ifAB
ν γ5 + 1
2 gAB
λµνσλµ
MAB = mAB + im5ABγ5 + aAB
µ γµ + bAB µ γ5γµ + 1
2HAB
µν σµν
CPT odd CPT even
25 Teppei Katori, MIT
06/30/10
Test for Lorentz violation in MiniBooNE follows LSND sidereal time analysis (1) fix the coordinate system (2) write down Lagrangian including Lorentz violating terms under the formalism (3) write down the observables using this Lagrangian Pνe→νµ ~ | (heff)eµ |2 L2 (c)2 = L c
2
| (C)eµ +(As)eµ sinw⊕T
⊕ +(Ac)eµ cos w⊕T ⊕ +(Bs)eµ sin2w⊕T ⊕ +(Bc)eµ cos2w⊕T ⊕ |2
sidereal frequency sidereal time
w⊕ = 2π 23h56m4.1s T⊕
Kostelecky and Mewes PRD70(2004)076002
Sidereal variation of neutrino oscillation probability for MiniBooNE
26 Teppei Katori, MIT
06/30/10
Expression of 5 observables (14 SME parameters)
NX NY NZ = cos χsinθcosφ − sinχcosθ sinθsinφ −sinχsinθcosφ − cos χcosθ (C)eµ = (aL)eµ
T −NZ(aL)eµ Z +E − 1
2 (3−NZNZ)(cL)eµ
TT + 2NZ(cL)eµ TZ + 1
2 (1− 3NZNZ)(cL)eµ
ZZ
(As)eµ = NY(aL)eµ
X −NX(aL)eµ Y +E −2NY(cL)eµ TX + 2NX(cL)eµ TY + 2NYNZ(cL)eµ XZ − 2NXNZ(cL)eµ YZ
(Ac)eµ = −NX(aL)eµ
X −NY(aL)eµ Y +E 2NX(cL)eµ TX + 2NY(cL)eµ TY − 2NXNZ(cL)eµ XZ − 2NYNZ(cL)eµ YZ
(Bs)eµ = E NXNY (cL)eµ
XX − (cL)eµ YY
( ) − (NXNX −NYNY)(cL)eµ
XY
(Bc)eµ = E − 1 2 (NXNX −NYNY) (cL)eµ
XX − (cL)eµ YY
( ) − 2NXNY(cL)eµ
XY
Kostelecky and Mewes PRD70(2004)076002
Coordinate vector (aL)µ = aµ + bµ (cL)µν = cµν + dµν Neutrino oscillation depends on coordinate
(χ = 48.2o, θ = 90.0o, φ = 180o)
27 Teppei Katori, MIT
06/30/10
28 Teppei Katori, MIT
06/30/10
Data set
lowE event and oscE event distribute equally in all run period. Unbinned K-S test: P(lowE,oscE)=0.73 (compatible)
29 Teppei Katori, MIT
Event distribution over runs
Preliminary
06/30/10
Data set
lowE event and oscE event distribute equally in GMT time. Unbinned K-S test: P(lowE,oscE)=0.76 (compatible)
30 Teppei Katori, MIT
Event distribution over GMT time
Preliminary
06/30/10
31 Teppei Katori, MIT
Event distribution over sidereal time
Preliminary
Data set
lowE event and oscE event distribute equally in sidereal time. Unbinned K-S test: P(lowE,oscE)=0.59 (compatible)
06/30/10
Systematic error study Time varying backgrounds are potentially dangerous…
Detector response may have day-night effect (temperature, etc), however effect is expected to be small
Proton beam follows Fermilab accelerator complex run plan, namely MiniBooNE receive high intensity beam
Maximum variation is ±6%, but this effect is washed out in sidereal time distribution. Full systematic error study is ongoing
proton on target day-night distribution
±6%
before POT correction after POT correction GMT lowE excess distribution Sidereal lowE excess distribution
32 Teppei Katori, MIT
06/30/10
Before SME parameter fit… Statistical test
33 Teppei Katori, MIT
Pearson’s χ2 test
06/30/10
Before SME parameter fit… Statistical test
All samples are consistent with flat However, data is not inconsistent with small variation, so we proceed to fit Pearson’s χ2 test d.o.f χ2 P(χ2) unbinned K-S test P(K-S) lowE GMT 107 107.6 0.47 0.42
83 69.6 0.85 0.81 lowE sidereal 107 106.0 0.51 0.13
82 76.2 0.66 0.64
34 Teppei Katori, MIT
Preliminary
06/30/10
LogL = −(µsig[SMEs]+ µbkgd)+ [µsig[SMEs]⋅ ωsig[SMEs,Eν
QE,T⊕ ]+ µbkgd ⋅ ωbkgd] i=1 n
Unbinned likelihood method
parameters (SME parameters) from unbinned likelihood fit. µsig[SMEs] : predicted number of signal event, function of SME parameters (=14). µbkgd: predicted number of background event, from MiniBooNE public data ωsig[SMEs,Eν
QE,T ]: probability density function (PDF) of signal distribution with
sidereal time, function of SMEs and reconstructed neutrino energy. ωbkgd: PDF of background (=const, background is assumed time independent).
35 Teppei Katori, MIT
06/30/10
parameter space data and best fit curve parameter space parameter space
Unbinned likelihood method, simultaneous fit result
Since data has good χ2 with flat hypothesis, fit improve goodness-of-fit
After fit χ2/dof=5.7/9, P(χ2)=0.77 Flat hypothesis χ2/dof=9.3/11, P(χ2)=0.60
low energy excess (200MeV<EνQE<475MeV
36 Teppei Katori, MIT
06/30/10
Only (C)eµ (sidereal time independent term) is statistically significantly non-zero extraction for (aL)µ and (CL)µν is rather simple due to special coordinate of MiniBooNE Unbinned likelihood method, simultaneous fit result
low energy excess (200MeV<EνQE<475MeV
parameter (C)eµ (As)eµ (Ac)eµ (Bs)eµ (Bc)eµ SME parameters (aL)T+0.75(aL)Z+0.35[-1.22(CL)TT-1.49(CL)TZ-0.33(CL)ZZ] 0.67(aL)Y+0.35[-1.33(CL)TY-0.99(CL)YZ] 0.67(aL)X+0.35[-1.33(CL)TX-0.99(CL)XZ] 0.35[-0.45(CL)XY] 0.35[-0.22((CL)XX-(CL)YY)] fit value
0.4
0.9 1.2 1.3
37 Teppei Katori, MIT
Preliminary
06/30/10
(aL)T (aL)X (aL)Y (aL)Z fit value (GeV)
0.6×10-20
statistic error (GeV) 0.9×10-20 1.9×10-20 1.8×10-20 1.2×10-20 (cL)TT (cL)TX (cL)TY (cL)TZ (cL)XX (cL)XY (cL)XZ (cL)YY (cL)YZ (cL)ZZ fit value 7.2×10-20
1.3×10-20 5.9×10-20
2.6×10-19 statistic error 2.1×10-20 2.8×10-20 2.6×10-20 1.7×10-20
0.8×10-19
38 Teppei Katori, MIT
low energy excess (200MeV<EνQE<475MeV
Unbinned likelihood method, simultaneous fit result
zero.
parameters correspond to sidereal time independent solution.
06/30/10
parameter space data and best fit curve parameter space parameter space
Unbinned likelihood method, simultaneous fit result
Since data has very little excess, best fit solution is consistent with zero. After fit χ2/dof=9.3/9, P(χ2)=0.41 Flat hypothesis χ2/dof=11.1/11, P(χ2)=0.44
(475MeV<EνQE<1300MeV
39 Teppei Katori, MIT
06/30/10
Since excess is small, all parameters are consistent with zero. Unbinned likelihood method, simultaneous fit result
(475MeV<EνQE<1300MeV
parameter (C)eµ (As)eµ (Ac)eµ (Bs)eµ (Bc)eµ SME parameters (aL)T+0.75(aL)Z+0.35[-1.22(CL)TT-1.49(CL)TZ-0.33(CL)ZZ] 0.67(aL)Y+0.35[-1.33(CL)TY-0.99(CL)YZ] 0.67(aL)X+0.35[-1.33(CL)TX-0.99(CL)XZ] 0.35[-0.45(CL)XY] 0.35[-0.22((CL)XX-(CL)YY)] fit value
0.4
1.0 1.7 1.4
40 Teppei Katori, MIT
Preliminary
06/30/10
LSND result
LSND collaboration, PRD72(2005)076004
parameter (C)eµ (As)eµ (Ac)eµ (Bs)eµ (Bc)eµ SME parameters (aL)T+0.19(aL)Z+0.04[-1.48(CL)TT-0.39(CL)TZ+0.44(CL)ZZ] 0.98(aL)X+0.053(aL)Y+0.04[-1.96(CL)TX-0.11(CL)TY-0.38(CL)XZ-0.021(CL)YZ] 0.053(aL)X-0.98(aL)Y+0.04[-0.11(CL)TX+1.96(CL)TY-0.021(CL)XZ+0.38(CL)YZ] 0.04[-0.052((CL)XX-(CL)YY)+0.96(CL)XY] 0.04[0.48((CL)XX-(CL)YY)+0.10(CL)XY] fit value
4.0 1.9
1.0 1.4 1.8
LSND data sidereal 3 parameter fit distribution 5 parameter fit
parameter spaces show 2 solutions in 1-sigma of the fit
41 Teppei Katori, MIT
LSND solution (0MeV<Ee+<60MeV)
06/30/10
LSND result
from published result (3 parameter fit)
SME parameters correspond to sidereal time dependent solution. (aL)T (aL)X (aL)Y (aL)Z fit value (GeV) 0.2×10-19 4.2×10-19
1.0×10-19 statistic error (GeV) 1.0×10-19 1.5×10-19 1.8×10-19 5.4×10-19 (cL)TT (cL)TX (cL)TY (cL)TZ (cL)XX (cL)XY (cL)XZ (cL)YY (cL)YZ (cL)ZZ fit value 0.3×10-18
2.1×10-18 1.3×10-18
statistic error 1.8×10-18 1.9×10-18 2.2×10-18 6.7×10-18
5.9×10-18
LSND collaboration, PRD72(2005)076004
LSND data sidereal 3 parameter fit distribution 5 parameter fit
42 Teppei Katori, MIT
LSND solution (0MeV<Ee+<60MeV)
06/30/10
Anti-neutrino oscillation data
MiniBooNE collaboration, Neutrino 2010, Athens
tomorrow (July 1) morning session
43 Teppei Katori, MIT
06/30/10
Lorentz and CPT violation has been shown to occur in Planck scale physics. MiniBooNE low energy excess νe data suggest Lorentz violation is an interesting solution of neutrino oscillations. Low energy excess events are statistically consistent with no sidereal variation. SME parameters are extracted under short baseline approximation. The flat solution is favoured. Full systematic study is ongoing. Recently antineutrino result will be analysed under SME formalism, too.
44 Teppei Katori, MIT
Teppei Katori, MIT 06/30/10
University of Alabama Bucknell University University of Cincinnati University of Colorado Columbia University Embry Riddle Aeronautical University Fermi National Accelerator Laboratory Indiana University University of Florida Los Alamos National Laboratory Louisiana State University Massachusetts Institute of Technology University of Michigan Princeton University Saint Mary's University of Minnesota Virginia Polytechnic Institute Yale University
45
06/30/10
46 Teppei Katori, MIT