⾃⼰回帰モデルを⽤いた重⼒波データ解析：
ブラックホール合体のリングダウン波形の抽出
真貝寿明，山本峻 元M2（大阪工大）
H. Shinkai, S. Yamamoto (OIT)
MergerRingdown

InspiralBH 準固有モード(quasi-normal modes)⇦ BH 摂動理論⇨ (M, a)
強い重⼒場の表れ⇨ ⼀般相対論の検証ができる

新しい方法を提案． テンプレートを使わず，データから波形を再構築．

2018/9/14 物理学会 ＠ 信州大学

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Ideal vs Reality(Theory vs Data Analysis)
GW150914(S/N=23.7)

h(t)

16.25

16.45

challenging for data analysis
GW data is with noise
signal quickly decays(M=60Msun, a=0.75 ̶> 300Hz, tau = 3 ms)
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

2

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Mock data example (0) : QNM extraction contest

modiﬁed ringdown signals from GRwith LIGO detectorʼs noise

真⾙ 寿明 (⼤阪⼯⼤)

1. Standard Matched-ﬁltering method
2. Improved Matched-ﬁltering method
3. Hilbert-Huang transformation method
4. Auto-Regressive method
5. Neural network method

2018/09/14 物理学会 @ 信州⼤学

3

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

1. Auto-Regressive model (Method, general) I
Fitting data with linear func.

e.g. xn = Ae rn t cos(!n t)
Z1 = e (r j!) t
Z2 = e (r+j!) t
1.00.8

0.6

0.4

0.5

0.2

5101520253051015

-0.2

-0.4
-0.5

-0.6

can be applied also to noisy data by adjusting M
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

4

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

1. Auto-Regressive model (Method, general) II
Fitting data with linear func.

• ﬁnd aj (Burg method)• ﬁnd M (FPE ﬁnal prediction error method)• re-construct wave signal from ﬁtted function• apply FFT with arbitrary precision.
power spectrum

・・・・・・・・・・

真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

5

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Auto-Regressive modelvs Short FFTsegment = 1/128 sec = 32 points

h(t)shift = 1/1024 sec = 4 points
The order M can be ﬁxed at 2〜8.
Even for short segment,
AR model shows precise
power-spectrum.

P(f)FFTARfreq.

[mock data, SNR=40, inspiral part]
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

6

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

1. Auto-Regressive model (Method, general) IIIFitting data with linear func.

• ﬁnd aj (Burg method)• ﬁnd M (FPE ﬁnal prediction error method)• re-construct wave signal from ﬁtted function• apply FFT with arbitrary precision.
power spectrum

characteristic eq.

真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

7

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Mock data example (1) ﬁtting well
Mock Data (Nakano02 p)

freq.

h(t)

spectrogramx original data, ̶ ﬁtted

8"mockN02p_segmentH.txt""mockN02p_segmentR.txt"

6420-2-4-6
1.593 1.594 1.595 1.596 1.597 1.598 1.599

1.6

1.601 1.602

Fitting data with linear func.

a.1
a.2a.3
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

= -2.235e+00= 1.869e+00= -5.545e-01
8⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Mock data example (2) spectrogram
Mock Data (Nakano02 p)

freq.

spectrogram

power spectrum

真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

9

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Mock data example (3) characteristic eq.
Mock Data (Nakano02 p)

Fitting data with linear func.

a.1
a.2
a.3

characteristic eq.

= -2.235e+00= 1.869e+00= -5.545e-01

xn 1 = zxn<latexit sha1_base64="N6OULxLyV4n+YDLj4gzs3pFIM1g=">AAACcHichVHLSsNAFD2N7/qquhFcWC2KCJYbERRBEN249NUq1FKSOK3BvEzSYlv8AX/AhRsVRMTPcOMPuPATxJ0V3LjwNg2IFvUOM3PmzD13zsyojqF7PtFTRGppbWvv6OyKdvf09vXHBgbTnl10NZHSbMN2d1XFE4ZuiZSv+4bYdVyhmKohdtTD1fr+Tkm4nm5b237ZEVlTKVh6XtcUn6nsca5qzcgnS5X4cc7KxRKUpCDizUAOQQJhrNuxG+xhHzY0FGFCwILP2IACj1sGMggOc1lUmXMZ6cG+wAmirC1yluAMhdlDHgu8yoSsxet6TS9Qa3yKwd1lZRwT9Ei3VKMHuqNn+vi1VjWoUfdS5lltaIWT6z8d3nr/V2Xy7OPgS/WnZx95LARedfbuBEz9FlpDX6qc1bYWNyeqk3RFL+z/kp7onm9gld606w2xeY4of4D887mbQXo2KVNS3phLLK+EX9GJEYxjit97HstYwzpSfO4RznCBy8irNCyNSmONVCkSaobwLaTpT1LLjq8=</latexit>z = exp[ 2⇡if t]<latexit sha1_base64="qLlUL+TuUzXdGJgBPG/gNGBc+b0=">AAACgnichVFNSxtBGH5ctWpsa9SLIMJgUKTQ8K4KiiiIevDoV1RIQtjdTuzgZnfZnQQ1ePLmH/DgSUFEvNmf0It/wIM/QTwq9NJD32wWxErtO8zMM8+8zzvPzNiBqyJNdN9itLa1f+jo7Ep1f/z0uSfd27cZ+dXQkTnHd/1w27Yi6SpP5rTSrtwOQmlVbFdu2buLjf2tmgwj5Xsbej+QxYq146mycizNVCk9dDBXkHuByH8dF4VACSXKorAkXW0JLYqldIayFId4C8wEZJDEip++RAHf4MNBFRVIeNCMXViIuOVhghAwV0SduZCRivclDpFibZWzJGdYzO7yuMOrfMJ6vG7UjGK1w6e43ENWCozQHV3RE93SNT3Q73/Wqsc1Gl72ebabWhmUeo4H1n/9V1XhWeP7i+pdzxplTMdeFXsPYqZxC6eprx2cPK3PrI3UR+mcHtn/Gd3TT76BV3t2Llbl2ilS/AHm38/9FmyOZ03KmquTmfmF5Cs6MYhhjPF7T2Eey1hBjs89whVu8MNoM74YpjHRTDVaEk0/XoUx+wct+pLf</latexit>

12
3x.r
0.962
0.962
1.447

x.i
0.566
-0.5660.000

f_R[Hz]346.800-346.800
0.000
真⾙ 寿明 (⼤阪⼯⼤)|x|
f_I[Hz]8.025e-01 71.721
8.025e-01 71.721
4.775e-01 240.931
2018/09/14 物理学会 @ 信州⼤学

10

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Mock data example (4) identify ring-down freq.
Mock Data (Nakano02 p)

characterestic eq.

xn 1 = zxn<latexit sha1_base64="N6OULxLyV4n+YDLj4gzs3pFIM1g=">AAACcHichVHLSsNAFD2N7/qquhFcWC2KCJYbERRBEN249NUq1FKSOK3BvEzSYlv8AX/AhRsVRMTPcOMPuPATxJ0V3LjwNg2IFvUOM3PmzD13zsyojqF7PtFTRGppbWvv6OyKdvf09vXHBgbTnl10NZHSbMN2d1XFE4ZuiZSv+4bYdVyhmKohdtTD1fr+Tkm4nm5b237ZEVlTKVh6XtcUn6nsca5qzcgnS5X4cc7KxRKUpCDizUAOQQJhrNuxG+xhHzY0FGFCwILP2IACj1sGMggOc1lUmXMZ6cG+wAmirC1yluAMhdlDHgu8yoSsxet6TS9Qa3yKwd1lZRwT9Ei3VKMHuqNn+vi1VjWoUfdS5lltaIWT6z8d3nr/V2Xy7OPgS/WnZx95LARedfbuBEz9FlpDX6qc1bYWNyeqk3RFL+z/kp7onm9gld606w2xeY4of4D887mbQXo2KVNS3phLLK+EX9GJEYxjit97HstYwzpSfO4RznCBy8irNCyNSmONVCkSaobwLaTpT1LLjq8=</latexit>z = exp[ 2⇡if t]<latexit sha1_base64="qLlUL+TuUzXdGJgBPG/gNGBc+b0=">AAACgnichVFNSxtBGH5ctWpsa9SLIMJgUKTQ8K4KiiiIevDoV1RIQtjdTuzgZnfZnQQ1ePLmH/DgSUFEvNmf0It/wIM/QTwq9NJD32wWxErtO8zMM8+8zzvPzNiBqyJNdN9itLa1f+jo7Ep1f/z0uSfd27cZ+dXQkTnHd/1w27Yi6SpP5rTSrtwOQmlVbFdu2buLjf2tmgwj5Xsbej+QxYq146mycizNVCk9dDBXkHuByH8dF4VACSXKorAkXW0JLYqldIayFId4C8wEZJDEip++RAHf4MNBFRVIeNCMXViIuOVhghAwV0SduZCRivclDpFibZWzJGdYzO7yuMOrfMJ6vG7UjGK1w6e43ENWCozQHV3RE93SNT3Q73/Wqsc1Gl72ebabWhmUeo4H1n/9V1XhWeP7i+pdzxplTMdeFXsPYqZxC6eprx2cPK3PrI3UR+mcHtn/Gd3TT76BV3t2Llbl2ilS/AHm38/9FmyOZ03KmquTmfmF5Cs6MYhhjPF7T2Eey1hBjs89whVu8MNoM74YpjHRTDVaEk0/XoUx+wct+pLf</latexit>

12
3x.r
0.962
0.962
1.447

x.i
0.566
-0.5660.000

f_R[Hz]346.800-346.800
0.000|x|f_I[Hz]8.025e-01 71.721
8.025e-01 71.721
4.775e-01 240.931

t
f_R (z_plane) f_I (z_plane) f_Rh(spectr) f_Rmax(spectr)
0.159375E+010.363837E+030.280414E+020.340000E+030.363000E+030.159668E+010.344258E+03
0.166608E+020.331000E+030.344000E+03
0.159766E+010.346800E+030.717212E+020.240000E+030.329000E+030.161230E+010.357677E+030.122067E+030.213000E+030.338000E+030.161328E+010.361098E+030.948919E+020.261000E+030.350000E+030.161523E+010.379918E+03
0.772796E+02
0.304000E+03
0.373000E+03
average & variance zfr =0.359E+03
0.118E+02 fr(sp) =0.350E+03
average & variance zfi =0.684E+02
0.365E+02
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

f_Rh(spectr)
0.384000E+03
0.357000E+03
0.382000E+03
0.431000E+03
0.422000E+03
0.432000E+03
0.148E+02

11⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Application to the LIGO/Virgo data

List of Detected GW eventsref.
GW150914
LVT151012
GW151226
GW170104
GW170608
GW170814
GW170817

M1+M2=Mf,spin

Mpc

Mdiﬀ/Mtotala̲ﬁnalz

PRL116, 061102

36.2+29.1=62.3+3.0(2016/2/11)

4.59%(2016/2/11)

23+13=35+1.5
2.78%

PRL116, 241103

14.2+7.5=20.8+0.9(2016/6/15)

4.15%

PRL118, 221101

31.2+19.4=48.7+1.9(2017/6/1)

3.75%

ApJ 851, L35

12+7=18.0+1.0(2017/12/18)

5.2%

PRL119,141101

30.5+25.3=53.2+2.6(2017/10/6)

4.66%

PRL119, 161101

1.36〜1.60 + 1.17〜1.36(2017/10/16)= 2.74 + ?

真⾙ 寿明 (⼤阪⼯⼤)

0.68

0.66

0.74

0.64

0.690.70?

2018/09/14 物理学会 @ 信州⼤学

410Mpc0.09
1000Mpc
0.20
440Mpc0.09880Mpc
0.18340Mpc
0.07
540Mpc0.1140Mpc

SNRdeg^2

23.7

600

9.7

13.0

850

13

1300

135201860

32.4

2812⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Ringdown wave of GW150914

LIGO paper

4096 sampling rate
100-400 Hz ﬁlter
1 segment = 1/64 sec= 64 points
1 shift = 1/512 sec = 8 points

AR modelHanfordP(f)max M = 3
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

freq.
13

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出Ringdown wave of GW150914

16384 sampling rate

AR modelHanford4096 sampling rate
100-400 Hz ﬁlter
1 segment = 1/64 sec= 64 points
1 shift = 1/512 sec = 8 points

AR modelHanfordP(f)max M = 3
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

freq.
14

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Ringdown wave of GW150914

16384 sampling rate

AR modelHanford真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

15

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Ringdown wave of GW150914
tf realf imag

We see QNM at 300Hz, 0.003s after the merger.

Berti, Cardoso & Will PRD 73, 064030 (2006).

真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

16

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Ringdown wave of GW150914, GW170814, GW170104

真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

17

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Ringdown wave of GW150914, GW170814, GW170104

Ｎ/A
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

18

⾃⼰回帰モデルを⽤いた重⼒波データ解析：ブラックホール合体のリングダウン波形の抽出

Summary & Outlook
⾃⼰回帰モデル x(t)短いデータ (〜 30 pts) に対しても精度よく周波数・減衰率を特定できる．
シグナルを⾒つけるのにテンプレートは不要．
LIGO/Virgo の 3 実イベントデータに適⽤，リングダウン部分の抽出を試みた．
LIGO/Virgoのは発表している (M, a) とGW150914, GW170814よく⼀致
S/Nの低い GW170104 ではいまいち(S/N=23.7)(S/N=18)(S/N=13)★リングダウン部分だけを検出できれば，強い重⼒場の重⼒理論検証ができる．
他の⽅法と組み合わせ，検出データ数が稼げれば，相対論検証ができるだろう．
★テンプレートを使わない⽅法は，今後，未知の重⼒波シグナルの候補検出に役⽴つかも．
真⾙ 寿明 (⼤阪⼯⼤)

2018/09/14 物理学会 @ 信州⼤学

19