ベイズモデルを用いた ブレーザーの時系列偏光データの解析 (Uemura, et al., 2010, PASJ, 62, 69)

Slides:



Advertisements
Similar presentations
かなた望遠鏡・フェルミ衛星によ る活動銀河核の多波長同時連携観 測 伊藤亮介、深沢泰司、かなたチー ム、 フェルミ LAT コラボレーション 多色同時観測 広島大学.
Advertisements

だい六か – クリスマスとお正月 ぶんぽう. て form review ► Group 1 Verbs ► Have two or more ひらがな in the verb stem AND ► The final sound of the verb stem is from the い row.
Report of recent DSSD status R. Kiuchi (SNU) 2012/10/20 E07
小水力班/ Small Hydro Generation Group 研究背景 / Research background
時系列データを手掛かりに2つの偏光成分を分離する ~ブレーザーへの応用~
英語特別講座 疑問文 #1    英語特別講座 2011 疑問文.
日本語の文法 文型(ぶんけい)をおぼえよう!
Chapter 11 Queues 行列.
2010年7月9日 統計数理研究所 オープンハウス 確率モデル推定パラメータ値を用いた市場木材価格の期間構造変化の探求 Searching for Structural Change in Market-Based Log Price with Regard to the Estimated Parameters.
What did you do, mate? Plain-Past
山口大学電波グループ ブレーザー電波データの紹介
X線偏光観測最前線 December 01, 2007 Tsunefumi Mizuno Hiroshima University
日本人の英語文章の中で「ENJOY」はどういうふうに使われているのか
Noun の 間(に) + Adjective Verb てform + いる間(に) during/while.
Let’s discuss in English! What are your opinions? Let’s discuss it!
D. M. Rust and B. J. LaBonte 2005, ApJ, 622, L 年6月6日 太陽雑誌会(速報)
SWIMS Current Status of Development
SP0 check.
Tohoku University Kyo Tsukada
Power Electronics center
GRB 観測 相対論的 Jet の内側を探る 金沢大学 米徳 大輔、村上敏夫 今日のトピックは Inverse Compton
Estimating Position Information by Detecting Network-Connection
The Sacred Deer of 奈良(なら)
Did he/she just say that? Get your head out of the gutter! Oh wait….
On / in / at Honoka Tanno.
know / knows(s) / ___________
VLBI位相遅延を使った飛翔体の位置計測観測 ー位相接続についてー
かなた望遠鏡・Swift/XRT,UVOT・MAXI によるMrk 421の長期観測
Kalman Filter Finite Element Method Applied to Dynamic Motion of Ground Yusuke KATO Department of Civil Engineering, Chuo University.
2018/11/19 The Recent Results of (Pseudo-)Scalar Mesons/Glueballs at BES2 XU Guofa J/ Group IHEP,Beijing 2018/11/19 《全国第七届高能物理年会》 《全国第七届高能物理年会》
P4-21 ネットワーク上の経路に対する 回帰問題について
神戸大大学院集中講義 銀河天文学:講義6 特別編 観測装置の将来計画
Causative Verbs Extensively borrowed from Rubin, J “Gone Fishin’”, Power Japanese (1992: Kodansha:Tokyo) Created by K McMahon.
Photometric properties of Lyα emitters at z = 4
全国粒子物理会 桂林 2019/1/14 Implications of the scalar meson structure from B SP decays within PQCD approach Yuelong Shen IHEP, CAS In collaboration with.
かなた望遠鏡/TRISPECによる変動天体観測
かなた望遠鏡を用いたブレーザーの 可視偏光変動の研究
Optical spectroscopy of flares from the black hole X-ray transient A in quiescence T. Shahbaz, et al., 2004, MNRAS, 354, /10/13(Wed) Wednesday.
情報源:MARA/ARMA 加 工:成田空港検疫所 菊池
Where is Wumpus Propositional logic (cont…) Reasoning where is wumpus
大規模なこと Large scale.
G. Hanson et al. Phys. Rev. Lett. 35 (1975) 1609
Polarization Multiplexing for Bidirectional Imaging
2019年4月8日星期一 I. EPL 84, (2008) 2019年4月8日星期一.
フレアの非熱的成分とサイズ依存性    D1 政田洋平      速報@太陽雑誌会(10/24).
研究会 「LHCが切り拓く新しい素粒子物理学」
FermiによるGRB観測を受けて CTAに期待すること
References and Discussion
著者:六川修一 著者:六川修一 原画像(左画像)は ©METI and JAXA[2007] Distributed by ERSDAC 著者:六川修一.
2019/4/22 Warm-up ※Warm-up 1~3には、小学校外国語活動「アルファベットを探そう」(H26年度、神埼小学校におけるSTの授業実践)で、5年生が撮影した写真を使用しています(授業者より使用許諾済)。
Term paper, report (2nd, final)
Affiliation:基礎工学部 材料工学科 准教授
シミュレーションサマースクール課題 降着円盤とジェット
どのような特徴を見ているのか ― 計算の目的
宇宙粒子線直接観測の新展開 柴田 徹 青学大理工 日本物理学会高知(22/Sep./2013).
ー生命倫理の授業を通して生徒の意識に何が生じたかー
Created by L. Whittingham
東北大 情報科学 田中和之,吉池紀子 山口大 工 庄野逸 理化学研究所 岡田真人
New Sources in the Sgr B & C Regions
Mixed Morphology (MM) SNR が予感するSNR研究の新展開
MO装置開発 Core part of RTR-MOI Photograph of core part.
非等方格子上での クォーク作用の非摂動繰り込み
total flux uncertainty : 11% next : momentum distribution uncertainty
~て しまう.
Measurements of J/ψ with PHENIX Muon Arms in 2003 p+p Collisions
Preflare Features in Radios and in Hard X-Rays
Cluster EG Face To Face meeting 3rd
原始星からのX線発見と課題 (r-Ophの)T-Tauri星からX線放射とフレアーの発見
アノテーションガイドラインの管理を行う アノテーションシステムの提案
Improving Strategic Play in Shogi by Using Move Sequence Trees
Presentation transcript:

ベイズモデルを用いた ブレーザーの時系列偏光データの解析 (Uemura, et al., 2010, PASJ, 62, 69) 植村誠 (広島大学)、他「かなた」望遠鏡チーム ベイズモデルを用いた ブレーザーの時系列偏光データの解析 (Uemura, et al., 2010, PASJ, 62, 69)

Blazars in the optical band Fosati, et al. (1998) blazar Optical-NIR High polarization hope fact A probe of the magnetic field structure But, no universal law has been established

Polarization variations in blazars Erratic variations Random motion in the QU plane? Blinks of a number of polarization components? Systematic variations Increase of the polarization degree (PD) with flares Rotation of the polarization angle (PA) A probe of the magnetic field in jets Are there universal characteristics in blazar polarization? The main theme of this talk We want to find a kind of “rules” in apparently erratic variations. Positive correlation between the flux and PD. (Smith, et al. 1986) Apparently random motion in the QU plane (Moore, et al. 1982)

「かなた/TRISPEC」による ブレーザー多色偏光撮像モニタープロジェクト できればシンプルで普遍的な、”erratic”でない観測的特徴をとらえたい キャンペーン期間:昨年秋~2009年3月 40天体ほど? いくつかの天体は2007年から順次モニター開始していた

「かなた」データの例 光度と偏光度が相関してる short flareもあれば、そうでないのもあり。。。  AO 0235+164:相関してるのもあるが、多少タイムラグがあるものもあるような。一方で全く相関してないのも。(笹田、他、天文学会)  PKS 1749+096:同じような時期、同じような振幅のフレアで相関してるのと、してないのと。

「かなた」データの例(続) BL Lacともあろうものが。。。 やはり、short flareには多少の反応がありそうだが。。。。。 というか、光度曲線に比べて偏光度がバタバタし過ぎ??? 光度変化はさほどでもないのに、突然偏光度があがったり。 OJ49は数か月のトレンドに偏光度が相関してない。。。。。。。

光度と偏光度の相関関係 めちゃくちゃ

“Bluer-when-brighter”: as a universal feature of blazars 32/42 blazars = well-observed objects (> 10 days) No significant correlation can be detected in 10 poorly-observed objects. 23/32 blazars = the bluer-when-brighter objects in their whole data. 2/32 blazars = from redder-when-brighter (faint state) to bluer-when-brighter (bright state) 3C 454.3 and PKS 1510-089 3 blazars = the bluer-when-brighter trends in year-by-year QSO 0454-234, PKS 1502+106, and PKS 0048-097 In total, 28/32 (=88%) blazars showing the bluer-when-brighter trend A sign of the trend can be seen even in the other 4 sources. Color-magnitude diagram of MisV1436

なぜ偏光の変動はめちゃくちゃなのか? 3つの仮説 多様性が圧倒的に大きい  フレアが、磁場の揃った(or amplifyされた)ところで起こることもあれば、揃ってないところで起こることもあるし、揃ってるところで起こったとしても途中で磁場の向きが変わることもあれば、あまり変わらないこともある。 色や偏光を観測しても普遍的な描像には届かない。。。   →せっかく観測してきたのに、この結論は避けたい。。。 無数の(ランダムな)細かいフレアの重ね合わせ 1つ1つのフレアは固有の偏光成分をもつ Moore et al. (1982), Impey et al. (1988), Jones et al. (1985, 1988)   →まだマシだけど、「ランダム」だとたいていの観測結果は「たまたま」で説明できてしまうので。。。 フレアに付随する偏光成分とそうでない成分の2成分があって、観測値はその2つが重なって見えている 常に偏光度の大きいような天体だと明らかに重要 Long term成分の推定が問題   →もしこの可能性があるのなら、そのような解の推定には意味があるだろう。

Multiple polarization components? Stokes’ QU parameters of BL Lac for ~20 yr (Hagen-Thorn, et al. 2002) Stokes parameters for linear polarization Systematic variation could be hidden by the presence of another polarization component. Preferential direction of polarization for years Our model assumptions: A flare is always associated with the flare of the polarized flux. There is a long-term variation trend in polarization. (0,0)

Bayesian method to estimate a long-term trend in polarization Q U t flux Prior distribution of the long-term trend (smoother curve is preferred). Posterior distribution of the long-term trend Likelihood function, maximized when the total flux and polarized flux completely correlate. The estimation of the parameters is done with the Markov Chain Monte Carlo (MCMC) method.

Demonstration with artificial data Test 1: Case for low frequency flares pol. flux total flux The long-term trend is reproduced as assumed (0,0) U Q

Demonstration with artificial data Test 1: Case for low frequency flares pol. flux total flux corrected pol. flux The long-term trend is reproduced as assumed long-term trend U Q

Demonstration with artificial data Test 2: Case for High frequency flares pol. flux total flux The long-term trend is successfully estimated. This method can extract a long-term trend even if observed variations are apparently just erratic. U (0,0) Q

Demonstration with artificial data Test 2: Case for High frequency flares pol. flux total flux The long-term trend is successfully estimated. This method can extract a long-term trend even if observed variations are apparently just erratic. U Q

Blazar monitoring with the 1.5-m “Kanata” Telescope Since 2006 TRISPEC developed by Nagoya Univ. Simultaneous optical and near-infrared observation photo-polarimetric mode is available Monitoring of blazars since 2007 42 sources

Case 1: OJ 287 pol. flux Data: Oct. 2008 – May 2009 Observed Q,U distribute NOT around (Q,U)=(0,0) A long-term trend is expected to be present. As expected, our method shows a long-term trend, possibly oscillating in a small range of PA with a time scale of a few tens of days. Good example for our model. total flux (0,0) U Q

Case 1: OJ 287 pol. flux Data: Oct. 2008 – May 2009 Observed Q,U distribute NOT around (Q,U)=(0,0) A long-term trend is expected to be present. As expected, our method shows a long-term trend, possibly oscillating in a small range of PA with a time scale of a few tens of days. Good example for our model. total flux U Q

Case 2: S2 0109+224 total flux Data: Aug.2008 – May 2009 Our model works well, even if the observed Q,U apparently distribute around (Q,U)=(0,0). Growth and decay of long-term trends? Another trend was born in the latter period? pol. flux U (0,0) Q

Case 2: S2 0109+224 total flux Data: Aug.2008 – May 2009 Our model works well, even if the observed Q,U apparently distribute around (Q,U)=(0,0). Growth and decay of long-term trends? Another trend was born in the latter period? pol. flux U Q

Case 3: S5 0716+714 Data: Jul. 2008 – Feb. 2009 The estimated long-term trend leads to no significant improvement of the correlation between the light curve and the polarized flux. No long-term trend (defined in our model) total flux pol. flux U (0,0) Q

Case 3: S5 0716+714 Data: Jul. 2008 – Feb. 2009 The estimated long-term trend leads to no significant improvement of the correlation between the light curve and the polarized flux. No long-term trend (defined in our model) total flux pol. flux U Q

Summary: advantage and disadvantage Advantages of our Bayesian model We can extract a systematic long-term trend from the observed variation, even if it is apparently erratic. The model may be useful for the test of the presence of the long-term trend. The long-term trend is NOT always obtained. Disadvantages of our Bayesian model The model is only valid for the polarized flux, not for the polarization degree in % (, or Q/I, U/I). better definition of the likelihood function or prior distribution? The model do NOT provide evidence for the long-term trend. Confirmation is needed by another kind of observations