High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro...

50
High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke Sakurai Supervisor: Daniel Mazin, Masaaki Hayashida

Transcript of High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro...

Page 1: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

High energy gamma-raysand

Lorentz invariance violation

Gamma-ray team A – data analysisTakahiro Sudo,Makoto Suganuma,

Kazushi Irikura,Naoya Tokiwa,Shunsuke Sakurai

Supervisor: Daniel Mazin, Masaaki Hayashida

Page 2: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Introduction

Page 3: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

What can we learn from γ-rays ?

• Motivation: to see whether the special relativity holds at high energy scale.

• Is there Quantum Gravitational effect, which modifies space-time structure and cause Lorentz invariance violation?

3

Page 4: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

How we measured:

• If QG makes space not flat, γ-rays of shorter wavelength are more affected, so higher energy γ-rays travel slower.

• Then, the speed of light is not constant!• So the arrival times of γ-rays emitted simultaneously

depend on their energies.

4

Page 5: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

What we measured:

• We measured arrival times of γ-rays of higher energies and lower energies.

• We determined ΔE, got Δt from data, and calculated “quantum gravity energy scale”

• We compared EQG of n=1 and 2 with Planck Energy scale.

5

Page 6: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

What we learn in this research:

• The meaning of EQG is the energy scale at which QG effects begin to appear.

• So if EQG is less than Planck energy scale, it means QG effect is detected

• The birth of a new physics!

6

Page 7: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Fermi Analysis

7

Page 8: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

About Fermi

launched from Cape Canaveral 11 June 2008

The Fermi satellite is in orbit around the earth today.

8

Page 9: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

About Fermi-Two Gamma-Ray        detectorsLAT(Large Area Telescope)->High energy rangeDetects Gamma-Rays of 20MeV-300GeV

GBM(Gamma-Ray Burst Monitor)->Low energy rangeDetects Gamma-Rays of 8keV-40MeV

http://fermi.gsfc.nasa.gov

9

Page 10: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Gamma-ray Burst Monitor(GBM)・ Detects Gamma-Rays of 8keV-40MeV  ( Low energy range )・ Views entire unoccupied sky

Instrument

GBM

Scintillator

10

Page 11: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Instrument

Large Area Telescope(LAT)

LAT

Detects Gamma-Rays of 20MeV-300GeV( High energy range )

Gamma-Ray converts in LAT to an electron and a positron.->1. Direction of the photon 2. energy of the photon 3.arrival time of the photon

11

Page 12: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Target Object(GRB)• GRB080916C(z=4.35±0.15

)– Hyper nova (Long

Burst≃a few 10 s) – (119.847,-56.638)

• GRB090510B(z=0.903±0.001)– The Neutron star

merging (Short Burst≃1s)

– (333.553,-26.5975)• Gamma-ray emission

mechanism not well understood

12

Page 13: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

GRB080916C Skymap

“Relative time” = Relative time to the onboard event trigger time. 13

Page 14: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Method

• Low energy range(GBM data)– How to decide to arrival time (tlow)

• High energy range(LAT data) – How to select photon– (Check a direction of photon’s source)– Decide to arrival time(thigh)

dt = thigh - tlow

14

Page 15: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

How to decide to arrival time(tlow)

σ=21 count

Here is tlow

Probability that count of noise is more than 5σ ~ 0.000001

15

Page 16: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

How to select highest energy photon

Use this photon

Here is tHigh

16

Page 17: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Result(Fermi)

GRB080916C(long-burst)Red Shift: z = 4.350.15Photon’s high energy : Ehigh = 13.31.07 [GeV]

Time lag: thigh – tlow ≈ 16.7[sec]

n=1:• EQG1(Lower limit) =[GeV]

n=2• EQG2(Lower limit) =[GeV]

17

Page 18: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Result(Fermi)

GRB090510B(short-burst)Red Shift: z = 0.9030.001Photon’s high energy : Ehigh = 31.1 2.5[GeV]

Time lag: thigh – tlow ≈ 0.84 [sec]

n=1• EQG1(Lower limit) = [GeV]

n=2• EQG2(Lower limit) = [GeV]

18

Page 19: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

MAGIC analysisGamma-rays from Blazars

Page 20: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

What’s MAGIC?NAME:

Major Atmospheric Gamma-ray Imaging Cherenkov(=MAGIC) Telescope

SYSTEM:Two 17 m diameter Imaging Atmospheric Cherenkov Telescope

ENERGY THRESHOLD:50 GeV

Page 21: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Atmospheric Cherenkov• Gamma-ray shower:

spreading narrow• Hadron shower:

spreading wide, background• Measuring Cherenkov Light:

both of showers make CL

21

Page 22: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Difference of image• Gamma-ray shower:

an ELLIPSE image, main axis points toward to the arrival direction

• Hadron shower:captured as somehow RANDOM image, using to reduce background

22

Page 23: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Stereo telescope• Ellipse image:

detectable direction• Stereo system:

compare MASIC1 with MASIC2 to detecting point

23

Page 24: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Targets• Mrk421:

An AGN, blazar, high peaked BL Lac, 11h04m27.3s +38d12m32s, z=0.030,Data got 2013/04/13

• S30218:An AGN, blazar, high peaked BL Lac, 02h21m05.5s +35d56m14s, z=0.944,Data got 2014/07/23-31

24

Page 25: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

MAGIC Data analysis

1. Distinguishing -ray shower from other showers.– Using the shape of shower and Montecarlo simulation.

• Use MAGIC standard software to verify the gamma-ray signal and the source position

• Reconstruct spectrum and light curve to select significant energy ranges and look for features in the light curve.

• S30218• Mrk421

25

Page 26: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

plot.

• Mrk421 • S30218We can use this energy range.

26

Page 27: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

MAGIC Data analysis(2)

2. Selecting energy bin for Flux vs. Time plot(Light curve).– Energy bins should be good-detection energy range.

500-2000 GeV and 2000-10000GeV (Mrk421) [GeV]

70-130 GeV and 130-200 GeV (S30218) [GeV]

3. Reconstruct light curve in determined energy bins.

27

Page 28: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

MAGIC Data analysis (3)

4. Normalize the light curve to the mean flux in the corresponding energy bin

5. Fitting the Light curve.– Using Gaussian and Linear function.

We allow these functions only to slide (strictly same shape) If these bins have the same origin, light curve must be the same.

– Calculate the delay of time• Simply we calculate the difference of Gaussian peak or point the linear

function crosses the time-axis(:crossing point).

28

Page 29: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Result (Mrk421)

[sec]

Actual Flux

Normalized Flux

Actual Flux

Normalized Flux

29

Page 30: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Result (S30218)

Actual Flux

Normalized Flux

[sec]

Actual Flux

Normalized Flux

30

Page 31: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Result of calculation and Estimate of

• The error of is too large… and in S30218 … So we estimate and discuss LOWER LIMIT of (LL:lower limit, UL:upper

limit)

• [GeV], [GeV] (Mrk421)• [GeV], [GeV] (S30218)

n=1 n=2

31

Page 32: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Discussion

Page 33: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Combined Result

– LL E_QG = Lower Limit of E_QG E_pl = Planck Energy scale = 2.435 e+18 GeV

33

Page 34: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Discussion

• In this research, we could not determine the value of E_QG.• We set lower limit for E_QG for n=1,2.• It’s possible quantum gravitational effect appears at energy scale

higher than 1.4 e+18 GeV

We can almost reach Planck Energy scale in gamma-ray astronomy!

34

Page 35: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Discussion

• Fermi data is the best for linear term(n=1).

Fermi

MAGIC

35

Page 36: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Discussion

• MAGIC data is the best for quadratic term(n=2).

MAGIC

36

Page 37: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Summary

• We analysed data from Fermi and MAGIC to calculate quantum gravitational energy scale.

• We set lower limits for E_QG and E_QG for n=1 and 2.Our limit for n=1 is close to Planck Energy Scale!!

• Fermi is the best for linear term while MAGIC is the best for quadratic term.

• We still have room for improvement especially for n=2.

More data from CTA will help!!

37

Page 38: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Back up

38

Page 39: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Odie

• To get the plot.• is proportional

to width of the shower.

To separate ONevent from OFF event.

39

Page 40: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

FLUTE(2)

• To determine energy bins using plot.

• And also to get light curve.(as I said in my presentation)

40

Page 41: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

plot

• What is chosen as -ray shower is “ON event”• What is thrown away as other shower is “OFF

event”Then simply calculate

When (), we detected the shower in a Energy range.

• Need to show image.41

Page 42: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Uncertainty of

• Flux follows the power law • So represent value of Energy bin() should be

logarithmic mean.

• Then, • And we assume that the uncertainty is 18%

Uncertainty of photons energy (Using energy bins) is 15%

Uncertainty by using logarithmic mean is 10%42

Page 43: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Uncertainty of

• We fit light curves using the Gaussian and linear function. Get the value of Gaussian peak and crossing point AND their error.

• , • is the mean value of and

• The Uncertainty of is calculated by using error propagation.

Each value has each error.

43

Page 44: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Determination of upper/lower limit

• Generally, • However, can be zero…• So we can determine only

• What is upper/lower limit?We use 95% and ,so we can’t simply use

44

Page 45: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

plot

• What is chosen as -ray shower is “ON event”• What is thrown away as other shower is “OFF

event”Then simply calculate

When (), we detected the shower in a Energy range.

• Need to show image.45

Page 46: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Error of Physics quantity

• Arrival time • Photon’s energy– dE/E = 0.08 in this energy range

46

Page 47: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Where is photon’s truly source?

• The probability of ratio that photon came from GRB source is 0.9999971

• There are background sources like as galactic gamma-rays and isotropic gamma-rays

47

Page 48: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Check a direction of photon’s source(1)

• A ・ B= |A||B|(sinθAcosφAsinθBcosφB

+ sinθAsinφAsinθBsinφB + cosθAcosθB)=|A||B|cosθ

B=(|B|sinθBcosφB, |B|sinθBsinφB, |B|cosθB)

A=(|A|sinθAcosφA, |A|sinθAsinφA, |A|cosθA)

θ

Difference of degree is 0.1degree!48

Page 49: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

Check a direction of photon’s source(2)

49

Page 50: High energy gamma-rays and Lorentz invariance violation Gamma-ray team A – data analysis Takahiro Sudo,Makoto Suganuma, Kazushi Irikura,Naoya Tokiwa, Shunsuke.

What kappa

50