12月24日
現状報告
POとMCDの2つのモデルが逆転しているのではないか
POのノルムを1に固定するか
初期値に澤口さんの値を入れてみる
パワースペクトルについてはQPOの論文
「Detection of Strong Short-term Variability in NGC6946 X-1」(ADS)
と同じ観測データを使っているので
問題はやはり設定か？
ちなみに上記の論文では
こんな形
パワースペクトル
12月25日 (log)
Raoさんの論文の図を目指す
> The power spectral densities (PSDs) for
> observations 0500730201 and 0500730101 from photons in the
> energy range of 0.3–10 keV are shown in Figures 2(a) and (b),
> respectively, and are normalized following Leahy et al.
(1983).
とあったのでその論文ADS
0.3〜10.0keVでフィルターをかけたPNデータのパワースペクトル(PNを使っているのは時間分解能が良いため)
これは101
似ているように見るが位置が違う(右がRaoさんの論文のパワースペクトル)
やはりbinの数が違う？
澤口さんの論文に
> フーリエ変換に使うライトカーブを
> 作る際にも、ビンまとめによってノイズが入ることがあるため、パワースペクトル作成の
> ためのライトカーブは、呼びだし時間 (=時間分解能) そのものを 1 ビンの大きさとした
とあったのでライトカーブ作成からbin sizeは0.0734 sでやる必要がある
(これまでの200という数字はバックグラウンドノイズを調べるためのライトカーブ作成時のbin sizeだった)
(一度何もない領域のパワースペクトルを調べる必要があるかも？)
0.0734でやった場合
バックグラウンドのパワースペクトルを調べる必要がある？
1月4日 (log)
あけましておめでとうございます
初期値に澤口さんの結果の値を入れてフィッティングした
101(2番目に長い) ログ
初期値は
phabs:nH>0.28
diskbb:Tin>0.21
norm>1786 (内縁半径から逆算)
powerlaw:PhoIndex>2.03
norm>
↓結果
========================================================================
Model phabs<1>(diskbb<2> + powerlaw<3>) Source
No.: 1 Active/On
Model Model Component Parameter
Unit Value
par comp
Data group: 1
1 1
phabs
nH
10^22 0.387381 +/-
4.29401E-02
2 2
diskbb
Tin
keV 0.164112
+/- 1.38771E-02
3 2
diskbb
norm
258.903 +/-
191.272
4 3 powerlaw
PhoIndex
2.18818 +/- 9.32003E-02
5 3 powerlaw
norm
1.89249E-04 +/- 2.22158E-05
Data group: 2
6 1
phabs
nH
10^22 0.387381 = 1
7 2
diskbb
Tin
keV 0.164112 =
2
8 2
diskbb
norm
258.903 = 3
9 3 powerlaw
PhoIndex
2.18818 = 4
10 3 powerlaw
norm
1.89249E-04 = 5
Data group: 3
11 1
phabs
nH
10^22 0.387381 = 1
12 2
diskbb
Tin
keV 0.164112 =
2
13 2
diskbb
norm
258.903 = 3
14 3 powerlaw
PhoIndex
2.18818 = 4
15 3 powerlaw
norm
1.89249E-04 = 5
________________________________________________________________________
Fit statistic : Chi-Squared
= 303.21 using 277
PHA bins.
Test statistic : Chi-Squared
= 303.21 using 277
PHA bins.
Reduced chi-squared
= 1.1147 for 272 degrees of
freedom
Null hypothesis probability = 9.368278e-02
内縁半径は
R=1365.3[km]
光度は
Model Luminosity 2.8328e+37 ergs/s
201(一番長い、長谷川と同じデータ) ログ
初期値は
phabs:nH>0.33
diskbb:Tin>0.18
norm>5818.6
powerlaw:PhoIndex>2.34
norm>
↓結果
========================================================================
Model phabs<1>(diskbb<2> + powerlaw<3>) Source
No.: 1 Active/On
Model Model Component Parameter
Unit Value
par comp
Data group: 1
1 1
phabs
nH
10^22 0.462084 +/-
5.03741E-02
2 2
diskbb
Tin
keV 0.137430
+/- 1.08148E-02
3 2
diskbb
norm
1115.53 +/-
919.507
4 3 powerlaw
PhoIndex
2.52976 +/- 8.08012E-02
5 3 powerlaw
norm
2.85265E-04 +/- 2.76502E-05
Data group: 2
6 1
phabs
nH
10^22 0.462084 = 1
7 2
diskbb
Tin
keV 0.137430 =
2
8 2
diskbb
norm
1115.53 = 3
9 3 powerlaw
PhoIndex
2.52976 = 4
10 3 powerlaw
norm
2.85265E-04 = 5
Data group: 3
11 1
phabs
nH
10^22 0.462084 = 1
12 2
diskbb
Tin
keV 0.137430 =
2
13 2
diskbb
norm
1115.53 = 3
14 3 powerlaw
PhoIndex
2.52976 = 4
15 3 powerlaw
norm
2.85265E-04 = 5
________________________________________________________________________
Fit statistic : Chi-Squared
= 277.33 using 322
PHA bins.
Test statistic : Chi-Squared
= 277.33 using 322
PHA bins.
Reduced chi-squared
= 0.87485
for 317 degrees of freedom
Null hypothesis probability = 9.473627e-01
内縁半径はR=2834.0 [km]
Model Luminosity 2.7112e+37 ergs/s
光度は相変わらずe+37のまま
他にも試してみた時のログ
101と201で内縁半径が倍ほども違うが
先行研究でもそうなのでキニシナイ
1月6日 (log)
中間報告レポートを書いた
1月7日 (log)
初期値をいじってみる