新星の偏光分光観測新井 彰
京都産業大学 神山天文台
Email: arai6a@cc.Kyoto-su.ac.jp

新星⚫ 激変星の一種⚫ 近接連星系におけるWD表面の熱核暴走反応によって引き起こされる質量放出を伴う増光現象⚫ 銀河系の化学進化に寄与

Ejecta mass: ~10-6 Mo to 10-4 Mo (depends on WD mass)
Ejecta Speed: from ~300 km/sec to ~8000 km/sec

30 Jan. 2021

連星・変光星研究会 20202

新星の光度曲線
可視極大 (T~10,000 K)

・急激な増光（遅いものもある）
・減光速度は、新星ごとに異なる
・減光の途中で、準周期変動やダスト形成よる
減光で可視減光を示す天体もある

火の玉期

静穏期（増光前）

静穏期
新星爆発
時間
Bode & Evans (1989, 2006)

30 Jan. 2021

連星・変光星研究会 20203

新星のスペクトル変化
HαHα[O III]

星雲期(nebular phase)

初期減光期
Hα

V339 Del

Optical Spectrum (FKO)

極大付近

Black bodyFree-Free continuum

時間
30 Jan. 2021

連星・変光星研究会 20204

新星の光度曲線と構造
可視極大 (T~10,000 K)

火の玉期
減光期（星雲期）

静穏期（増光前）

静穏期
新星爆発
時間
Bode & Evans (1989, 2006)
鳴沢真也著「へんな星たち」

30 Jan. 2021

連星・変光星研究会 20205

新星の光度曲線と構造[e.g., McLaughlin 1945, 1947; Payne-Gaposchkin
1957; Friedjung 1987; Williams & Mason 2010;
Arai et al. 2016; Kawakita et al. 2019, Aydi et al. 2019, 2020]v.s. Single flow ejecta (e.g.; Mason et al. 2018)

初期の放出構造が注目されている

火の玉期

？
？
静穏期（増光前）(e.g,; Ribeiro+2009, Pavana+2020)

時間(偏光分光; Kawabata+2006)
30 Jan. 2021

連星・変光星研究会 2020

6

なぜ爆発の初期が面白いのか？
• 爆発直後～極大期のジオメトリはよくわかっていない.
→ いつ、どのようにしてBiploar flow になるのか？
• 初期の増光中の放出構造.
→ Wolf-rayeat 的な外層(Arai+2015)→ 極大前に輝線フェーズあり(e.g.,Tanaka+2011, Imamura+2012)→ イジェクタの構造モデル(e.g.; Hachisu+2016, Li et al. 2017), TNRによる対流モデル(e.g.;
Jose et al. 2020)の検証など• 7Be II, 7Li IIやCN, C2, CO分子の吸収線から同位体元素の組成が測定できる.
ダスト形成を起こす新星もある(全体の~25%).
→ より詳しい放出量を推定するには構造を知りたい

現時点では、球対称を仮定している
新星爆発モデル, TNRモデルの評価
30 Jan. 2021

連星・変光星研究会 2020

7

新星の核合成と銀河の化学組成• 銀河の新星発生: ~ 30 yr-1
• 銀河の年齢: ~1010 yr• 新星の１個あたりのイジェクタ質量：~ 2×10-5 M☉~ 6×106 M☉ ：銀河のガス＆ダストの約1/3000程度

José+(1998)

大部分の元素に関しては、
銀河の化学進化に及ぼす影響は僅か.
しかし、一部の同位体元素(7Li, 13C, 15Nなど) の組成比は太陽組成の 100～1000倍.
→ 銀河の化学進化に寄与
30 Jan. 2021

連星・変光星研究会 20208

新星の化学組成• CNOの過剰• hot-CNO反応によるCNO元素同位体の組成過剰
観測値：

Tarasova et al. (2016)30 Jan. 2021

理論計算：José et al. (2013)

連星・変光星研究会 20209

新星による7Liの合成・放出
理論的な予言:Cameron & Fowler (1971), ApJ 164, 111
Starrfield+(1987), ApJ 222, 600

Tajitsu+2015, 2016; Izzo+2015, 2018,
Molaro+2016,2020; Selveli+2019; Arai+ in prep.
爆発+7～80d という比較的初期に検出される.
V339 Del (+47 d)

Credit : NAOJ
30 Jan. 2021

Tajitsu+ (2015), Nature, 518, 381
連星・変光星研究会 202010

新星と銀河のリチウム量の進化
観測量

理論予測

Arai et al.
submitted to ApJ.

リチウムの組成比は理論予想よりもかなり多い
30 Jan. 2021

(クレジット：国立天文台)
https://www.subarutelescope.org/jp/results/2015/02/18/936.html

小質量星起源のリチウムは、新星が主要な起源だろう[e.g.; Cescutti & Molaro 2019, MNRAS, 482, 4372]
連星・変光星研究会 202011

炭素・窒素の同位体：プレソーラー粒子との関係• 爆発初期に現れる分子スペクトルから、
14N/15Nや12C/13Cの同位体組成比が得られる
AGB• プレソーラー粒子と近いものもある→ 太陽系の起源に関係？

CarbonStar?

Novae?

(e.g., Kawakita+2015, Kawakita & Arai 2017)

SNe II
Nova V2676 Oph
12C/13C〜4
14N/15N〜2

The figure of Isotopic ratios of pre-solar
grains are quoted from Jose +(2004).
30 Jan. 2021

Nova

12C/13C

14N/15N

Molecule

References

DQ Her

≥ 1.5

≥2

CN

Sneden & Lambert (1975)

NQ Vul> 3COFerland et al. (1979)

V842 Cen

2.9±0.4

COWichmann et al. (1990, 1991)

V705 Cas≥5

CO

Evans et al. (1996)

V2274 Cyg

1.2±0.3

CO

Rudy et al. (2003)

V2615 Oph>2

CO

Das et al. (2009)

V496 Sct

≥ 1.5

CO

Raj et al. (2010)

1.3

CO

Rudy et al. (2009, 2012)

C2, CN

Kawakita et al. (2015)

CO

Banerjee et al. (2016)〜

V2676 Oph

連星・変光星研究会
2020
V5668
Sgr

4

〜

1.5

〜

2

〜

12

新星の爆発放出物⚫ 多様性に富む爆発放出物の幾何（爆発後の直接撮像や分光観測より)⚫ 極大付近でのガンマ線の検出 → 爆発から衝撃波までの時間差を示唆

Gill & O’Brien 1998, 2000; O’Brien & Dunlop 1995
30 Jan. 2021

Credit: NASA/DOE/Fermi LAT Collaboration
連星・変光星研究会
2020
Credit: NASA/DOE/Fermi
LAT Collaboration

13新星の物質放出モデル（２段階放出シナリオ）
Phase I (爆発直後)

Phase II (可視極大付近)

L2 (L3) ラグランジュ点から放出される
高密度, 低速度な公転面にスパイラル構造の放出物

低密度,高速度の等方的な放出物→ 低速度放出物との衝突でガンマ線輝線の放出+ 温度が冷えてくるとダスト生成

Top view

Side view
Side view
14 /26Li et al. 2017, Nat. Astron., 1, 697
30 Jan. 2021

連星・変光星研究会 2020

新星の物質放出モデル（２段階放出シナリオ）
Phase I (爆発直後)

Phase II (可視極大付近)

L2 (L3) ラグランジュ点から放出される
高密度, 低速度な公転面にスパイラル構造の放出物

低密度,高速度の等方的な放出物→ 低速度放出物との衝突でガンマ線輝線の放出+ 温度が冷えてくるとダスト生成⚫ このシナリオは以下を説明可能(1) ガンマ線輝線(2) ダスト形成(3) 新星爆発放出物の非対称構造⚫ しかしながら、新星爆発直後の爆発放出物の幾何について観測的にほとんど明らかとなっていない– 大部分の新星は大望遠鏡で撮像観測をしても点源である
Top view

Side view
Side view

Li et al. 2017, Nat. Astron., 1, 697
30 Jan. 2021

15 /26連星・変光星研究会 2020

新星の放出物の構造を探る方法
爆発初期
・干渉計：全体の構造
・吸収線：視線方向の構造
・偏光分光：全体構造・視線方向の構造
爆発中〜後期・偏光分光：全体像・非対称放出物の構造
・星雲輝線の解析([O III]など)：全体的な電離ガスの構造
・直接撮像（> 10年後）：全体像
偏光分光
干渉計
吸収線
輝線

直接撮像

30 Jan. 2021

連星・変光星研究会 2020

偏光分光観測は情報量は多い.
初期の構造を探る上で有用
16

新星の偏光分光観測
これまでの新星の偏光分光観測
Nova Cyg 1986 Whitney & Clayton (1989), AJ, 98, 297
BY CirJohnson et al. (1997), AJ, 113, 2200
Nova Cyg 1992 Bjorkman et al. (1994), ApJ, 425, 247 ●V4444 Sge
Kawabata et al. (2000), ApJ, 540, 429
U Sco, NSgr1998 Ikeda et al. (2000), A&A, 335, 256
V1494 Aql
Kawabata et al. (2001), ApJ, 552, 782 ●V475 Sct
Kawabata et al. (2006), AJ, 132, 433 ●V339 Del

30 Jan. 2021

Kawakita, Ahinnnaka, Arai, Arasaki, Ikeda (2019), ApJ, 872, 120 ●連星・変光星研究会 202017

偏光の機構• 双極子散乱
入射光により電気振動子が励起されて放射– 電子散乱– レイリー散乱– ミー散乱– 反射・屈折

川端弘治 氏 “偏光観測チュートリアル
一般的事項と可視近赤外観測”(公開資料)より借用
30 Jan. 2021

連星・変光星研究会 202018

偏光スペクトル：q-u 平面
規格化強度, I

North

ストークスQ: 0°/90° 成分
Q

U

Q =

45°

–
0 deg

規格化したストークスq = Q/I

ストークスU : 45°/135° 成分

規格化したストークスu = U/IU =–45 deg
偏光度, P

90 degu

135 degq-u 平面図

q = Q/I
u = U/I𝒒, 𝒖
偏光方位角, θ

偏光フラックス, I xP𝑃 =𝑞 2 + 𝑢22𝜃 = arctan𝑢𝑞q North

Kawakita et al. 2019, ApJ, 872, 120
30 Jan. 2021

連星・変光星研究会 202019

偏光スペクトル：q-u 平面
North

規格化強度, I

ストークスQ: 0°/90° 成分
Q =

–
0 deg

規格化したストークスq = Q/I

ストークスU : 45°/135° 成分

規格化したストークスu = U/IU =45 deg

連続光偏光 (星間 + 新星固有)偏光度, P

90 deg–135 degq-u 平面図

u𝒒, 𝒖
偏光方位角, θ

線偏光(新星固有)
連続光偏光(星間 + 新星固有)q North

偏光フラックス, I xP

Kawakita et al. 2019, ApJ, 872, 120
30 Jan. 2021

連星・変光星研究会 202020

連続光の偏光と幾何光球面付近で生じる電子散乱の模式図
偏光した光

電子

光球面

光球面

観測者

無偏光の光

光球

恒星

球対称な形状の場合→ 無偏光

非対称 → 偏光あり

⚫ 新星固有の連続光偏光は有効光球面の非対構造を反映# 星間偏光の推定が必要⚫ 連続光の偏光から爆発直後の新星の光球面の形状を推定する→ 幾何は偏光方位角と垂直方向に偏る
30 Jan. 2021

連星・変光星研究会 2020

21

線偏光と幾何のモデル計算
NA) 双極ジェット

Hoflich (1991), A&A, 246, 481
Recent studies: Tanaka et al. 2012, ApJ, 754, 63; 2017, ApJ, 837, 105B) トーラスC) クランプ

E

PA = θ= 90 deg
線偏光を用いる利点⚫ 光源周囲で卓越する幾何の非対構造と速度場が得られる– 偏光方位角は有効幾何の方位角を反映– 撮像観測で空間分解できない天体についても可能⚫ 星間偏光の影響を無視できる– 一般に星間偏光の高精度の推定は難しい⚫ 異なる成分から異なる領域の幾何を推定できる– 分子・原子により、イオン化エネルギーなどが異なるため、成分ごとに異なる領域が見える

30 Jan. 2021

直線

反転

2θ= 180 deg

2θ= 180 deg

連星・変光星研究会 2020

ループ

22 /26

偏光分光観測
線偏光
光源周囲の物質の幾何と運動を反映

連続光の偏光
光球面の幾何を反映

光球面

球対称 (→ 無偏光)

非対称 (→ 偏光あり)

30 Jan. 2021

連星・変光星研究会 202023

V339 Delの線偏光分光観測
VESPolA (ベスポラ)
Very precious Echelle SpectroPolarimeter on Araki-telescope(Arasaki, Ikeda, Shinnaka et al. 2015, PASJ, 67, 35)

特徴⚫ 高波長分解能 (R~8000)かつ高偏光決定精度(δP<0.1%)で複数成分を一度に取得できる

VESPolA
荒木望遠鏡に搭載したVESPolA
京都産業大学 神山天文台
Longitude 35° 04’ 13’’N
Latitude 135° 45’ 29’’E
Altitude 162 m⚫ 小口径望遠鏡 (神山天文台 1.3m-荒木望遠鏡)を
主力望遠鏡とすることで、変光星や突発天体等の時間変動する天体のモニター観測が可能
項目

波長分解能
波長範囲
偏光決定精度
限界等級 (𝛿P < 0.1%, 4 h)
30 Jan. 2021

連星・変光星研究会 2020

仕様

R ~ 8,000(~1 Å or ~37 km/s at H𝛼)𝜆 = 550 – 830 nm
𝛿P < 0.1%

Mv = 7.8 mag
24

V339 Delの線偏光分光観測

発見

0.9 d

→
→

2.2

2.9 4.1

4.9

→
→

→

→
→

発見からの日数 (発見日：2013 Aug 14.56)
0
1.4
3.4
9.4
5.4
7.4

5.9

7.0

観測開始時刻(世界時)

発見からの日数 (日)

積分時間(s)2013 Aug 15.49

0.91 d6,000

2013 Aug 16.73

2.15 d1,200

2013 Aug 17.51

2.93 d3,600

2013 Aug 18.69

4.11 d1,200

2013 Aug 19.49

4.91 d2,800

2013 Aug 20.48

5.90 d1,600

2013 Aug 21.53

6.95 d2,000# Discovery: 2013 Aug 14.56 (by K. Itagaki)
From AAVSO database

Date in 2013

可視極大前の偏光スペクトルは、新星爆発におけるWD表面でのTNRの様子を詳細に研究するにあたり本質的に重要– 可視極大前は新星の光球面と爆発放出物の最前面とが一致する (火の玉期)
30 Jan. 2021

連星・変光星研究会 202025

古典新星V339 Del (= いるか座新星2013)
発見からの日数 (発見時刻：2013 Aug 14.56)

0

10

24

38

52

66

80

© The Virtual Telescope⚫ 板垣公一氏により発見
発見(2013 Aug 14.58)⚫ 近年で非常に明るかった新星⚫ 爆発直後から様々な時期で様々な観測手法による多くの観測

30 Jan. 2021

連星・変光星研究会 202026

古典新星V339 Del (いるか座新星2013)
発見からの日数 (発見時刻：2013 Aug 14.56)

0

10

24

38

52

66

80

可視光極大 (mV = 4.46 mag at t = 1.9 d)
Munari et al. 2013, IBVS, 6080, 1

γ-線輝線の検出 at t ~4–13 dAckermann et al. 2014, Science, 345,554

ダスト形成 from t ~30 d
Burlak et al. 2015, BaltA,
24, 109© The
Virtual Telescope⚫ 板垣公一氏により発見
発見(2013 Aug 14.58)

7Beの検出 at t =48, 49 d

Tajitsu et al. 2015, Nature, 518, 381⚫ 過去10年で最も明るい新星⚫ 爆発直後から様々な時期で様々な観測手法による多くの観測

From AAVSO
database(LCG: https://www.aavso.org/lcg)
30 Jan.
2021

連星・変光星研究会 202027

V339 Delの線偏光分光観測
新星 V339 Del

Si II

Hα

Kawakita, Shinnaka et al. 2019, ApJ, 872, 120

30 Jan. 2021

連星・変光星研究会 202028

Hα線の線偏光 (t = 0.91 d)
P Cyg profile
連続光偏光を差し引き後B) トーラス2θ
2θ= 180 deg

2θ=119.8 ± 2.2 deg→ PA ~60 deg

PA = θ= 90 deg
Tanaka et al. 2017, ApJ, 837, 105⚫ “光球面周囲に膨張するトーラス構造”
が極大前の新星の爆発放出物の幾何として最もありうる幾何構造⚫ 推定されるPAtorus = 59.9 ± 1.1 deg
Shinnaka et al. 2019, ApJ, 872, 120
30Kawakita,
Jan. 2021

連星・変光星研究会 202029

連続光偏光 (Hα波長)⚫ V339 Delの光球面の有効幾何が日ごとに変化⚫ 火の玉期 (t <~3 d)では、連続光成分の偏光は爆発放出面の幾何を反映(t = 2.2 d)– この期間は光球面と新星爆発放出物の先端とがほぼ一致

ISP⚫ Aug 15 → Aug 16で有効幾何が90度変化 (q-u平面で原点に対して180度移動)(t = 2.9 d)(t = 0.91 d)– Aug 15 (t = 0.9 d): θ~150 deg– Aug 16 (t = 2.2 d): θ~60 deg– Aug 17 (t = 2.9 d): ~ISP
# V339 DelのRバンド偏光撮像観測から同新星の星間偏光を推定 (Shakhovskoy et al. 2017)(q, u)ISM = (–0.07%, 0.35%)

Kawakita, Shinnaka et al. 2019, ApJ, 872, 120
30 Jan. 2021

連星・変光星研究会 202030

Si II線の線偏光 (t = 0.91 d)C) Radii of the clumps

Loop⚫ Si IIは光球面上に非均質に存在–Tphotosphere = 9,000– 10,000 K at Aug 15(Skopal et al. 2014, A&A, 569, A112)– Si II (I.P.= 16.3 eV) was ionized to Si III⚫ 新星の光球面は非一様な温度分布を持つ→ WD表面での非一様なTNR (熱核暴走反応)
Kawakita, Shinnaka et al. 2019, ApJ, 872, 120
30 Jan. 2021

連星・変光星研究会 202031 /26

偏光分光観測から推測したV339 Delの幾何(A) t = 0.91 d (2013 Aug 15)(B) t = 2.15 d (2013 Aug 16)
Bipolarshaped ejecta

Fast expanding
photosphere
PAtorus~60 deg

有効幾何が90度変化

Slowly expanding torus

PAbipolar~150 deg

Kawakita, Shinnaka et al. 2019, ApJ, 872, 120⚫ 光球面の幾何が爆発直後は1日以下の時間スケールで変化している– 可視光極大前はトーラス構造が存在する (Hαの線偏光より)– t = 0.91 dとt = 2.15 dで有効光球面が90度回転 (連続光偏光より)⚫ V339 Delの爆発メカニズムとして2段階放出モデル (Li e al. 2017など) を支持
30 Jan. 2021

連星・変光星研究会 202032V339 Delに見られた爆発初期の非対称構造
CHARAアレイによる赤外線干渉計 (H-band)

PA~115 deg

PA~149 deg

PA~157 deg

PA = 128°±30°(爆発後2–9日の平均値)

Schaefer et al. 2014, Nature, 515, 243
30 Jan. 2021

連星・変光星研究会 202033

議論： Hバンド干渉計による直接撮像発見からの日数 (発見日：2013 Aug 14.58)
0

1.4

3.4

5.4

7.4

9.4

t = 2.2 d(PA~150 deg)→NE

t = 9.0 d
PA = 139 ± 16 deg→

t = 5.0 d
PA = 149 ± 2 deg

t = 2.0 d

発見

PA = 82 ± 32 deg
t = 0.9 d(PA~60 deg)t = 3.0 d
PA = 115 ± 11 deg

t = 6.9 d
PA = 157 ± 4 deg

干渉計観測と偏光分光観測から推定された有効幾何のPAの変化はほぼ一致(PA~60 deg → ~150 deg)
From AAVSO database(LCG: https://www.aavso.org/lcg)
30 Jan. 2021

Date in 2013
連星・変光星研究会 2020

Schaefer et al. 2014, Nature, 515, 24334

議論：ダスト生成領域
R-バンド偏光撮像観測

ダスト形成前 (t =7–25 d)
ダスト形成後 (t = 25–65 d)
ダスト形成後 (t >65 d)(Shakhovskoy et al. 2017)ダスト形成

2θ= 330 ± 10 deg→ θ= 165 ± 5 deg# V339 DelのRバンド偏光撮像観測から同新星の星間偏光を推定 (Shakhovskoy et al. 2017)(q, u)ISM = (–0.07%, 0.35%)
30 Jan. 2021

連星・変光星研究会 202035

議論：ダスト生成領域

N
PA~60度

双極放出物
ダスト形成前 (t =7–25 d)
ダスト形成後 (t = 25–65 d)
ダスト形成後 (t >65 d)

E

トーラス(スパイラル)2θ= 330 ± 10 deg→ θ= 165 ± 5 deg# V339 DelのRバンド偏光撮像観測から同新星の星間偏光を推定 (Shakhovskoy et al. 2017)(q, u)ISM = (–0.07%, 0.35%)
30 Jan. 2021

⚫ ダストはθ=165 ± 5 deg (PA~75deg) 方向で生成したようだ– 偏光観測から推定されるトーラス構造(PA~60 deg)と似た方向⚫ スパイラル放出物中での生成を示唆– 低速度成分と高速度成分との衝突でダストが生成する (Derdzinski et al. 2017)
連星・変光星研究会 2020

36

V339 Delの高分散偏光分光のまとめ
まとめ (詳しくは Kawakita et al. 2019, ApJ, 872, 120)
⚫ “偏光分光観測” が新星の爆発直後の爆発放出物の空間分布を明らかにする非常に強力な手法であることを示した– 爆発直後は視直径が非常に小さく撮像観測による空間分解が困難な天体でも可能– ISPを正確に推定できなくても推定可能⚫ 得られた新星V339 Delの幾何構造– 可視光極大前：トーラス形状 (PA~60 deg、Hα線の線偏光より)– 可視光極大後：有効幾何のPAが90度変化 (連続光偏光より)→ 2段階放出モデル (低速スパイラル放出 + 高速等方放出) を支持– WD表面における非均一なTNRを示唆 (Si II線の線偏光より)

30 Jan. 2021

連星・変光星研究会 202037

新星の光度曲線と構造~2dKawakita+2019は構造変化の様子を観測的に追加

火の玉期

減光期（星雲期）

静穏期（増光前）(e.g,; Ribeiro+2009, Pavana+2020)

30 Jan. 2021

連星・変光星研究会 2020

時間38

先行研究との比較
Bjorkman et al. (1994) for nova Cyg 1992

30 Jan. 2021

連星・変光星研究会 202039

先行研究との比較
Kawabata et al. (2001) for nova V1494 Aql

連続光の偏光の急激なPAの回転→ 初期の構造変化を捉えている
V339 Delの結果と似ている

30 Jan. 2021

連星・変光星研究会 202040

先行研究との比較
Kawabata et al. (2006) for nova V475 Sct

赤方偏移した偏光成分を検出.
双極上のwindの存在を示唆
30 Jan. 2021

連星・変光星研究会 202041

新星の偏光分光観測：まとめ• 偏光分光観測は、干渉計観測と並び、初期の新星のイジェクタの非対称形状を探る強力な手段の１つ• 観測例が極めて少ない• 偏光度が小さいので、高精度で観測可能な、明るい新星にしか実施できない• 装置が非常に少ない（ただし、日本にはある）• 偏光分光観測の展望• V339 Delなどで見られた形状進化は、普遍的なのか？
• 新星放出物の構造モデルの検証• ダスト形成領域の制限にも有効である

30 Jan. 2021

連星・変光星研究会 202042