星形成ゼミ 2018/3/2
SFN #301 50-52, #302 1-5
崔仁士

1# 301505152# 302

123452

#51

•
•
•

ALMAでTaurusのprotoplanetary diskをsurvey
disk massを導出し統計的に議論
low mass star (spectarl type M4 – M7) のサンプルを埋めた

観測と結果• ALMA Cycle 1, 885 µm continuum (ALMA Band 7)• samples: Herschelを用いたTBOSS surveyから• M type の Class II のみ24天体• stellar: M4 – M5 の14天体
substellar: M6 – M7の10天体• 空間分解能 ~ 0’’.34x0’’.24• 21天体が8σ以上でdetection，J0414+2811のみ5σでdetectionで，計22天体が受かった• 3σの検出限界は0.27 mJy/beam

Herschel 70 µm flux (mJy)概要

3円盤質量の測定
scaling relation• κ1.3mm = 2.3 cm2 g-1, κν ~ ν0.4
• r = 100 au を仮定→ Mass range~ 0.3−20 Mearth
右上がりの傾向

power-law fitting• Mdust vs M*で正の相関 (∝M*0.94±0.14)• フィッティング結果はAndrews+13 と一致する結果→ low mass でも円盤形成プロセスは同じ• Upper Sco との比較
• slopeはTaurusとほぼ同じ（ ∝M*0.92±0.18，ただし恒星進化モデルに依存）• より進化の進んだUpper Scoではdisk massが減少している4

• MMSN gas massを超えるのは24%• M-star rangeでは17%↔M-dwarf周囲のgas giantの存在頻度は~2-3 %→ gas giantを作る効率は低い• brown dwarf 周りには10Mearthを超えるものも少ない→ brown dwarf周りのplanetary mass companionはdisk起源ではない？
• Kepler で見つかった惑星質量をプロット(Mulders+15，惑星周期は2-50 days)• 見つかった惑星よりもdiskにはもっと質量がある

browndwarf

M-dwarf• Upper ScoではTaurusと比べてdisk massは減少する
が分布は同じ→消えたmassはplanet ?
• 惑星形成は1-2 Myr のTaurusから~10 Myr のUpper
Scoの間で起こっている

5

#1
概要
• 大質量形成過程についてよく知りたい• Infrared dark clouds (IRDCs) の一つである IRDC G034.43+0024 を観測を行った• IRDC G034.43+0024• filament 構造の中に複数の重いcores• dense, quiescentなgas• 速度構造からCloud の形成などを議論

観測••
•
•

IRAM 30 m telescope/EMIR
N2H+ (1-0) と C18O (2-1) 輝線を観測
それぞれdense core，dense coreよりさらに広がった
envelopeをtraceすると期待
空間分解能 ~ 0.5 pc

6

観測結果
surface density

積分強度

速度分布

速度分散

N2H+
N2H+C18O

C18O

コントア: 積分強度

• ピークがいくつかのコアに対応（MM1, MM3 for N2H+，MM1, MM3, MM7, MM8, F1 for C18O）• 東西方向に速度勾配• 速度分散の振る舞いが異なる• 異なる速度成分を分離する• semi-automated gaussian fitting
Mapをいくつかに分けた領域 (SAA) で複数のGaussianをfit
automated hierarchical clusteringalgorithm
Gaussian fitされた各SAAsの分類を行う
1. 各成分のseparationがbeam size以下
2. 速度差が速度分解能の2倍以下
3. 速度分散の差が0.23 km/s 以下を満たすなら同じ成分とみなす•7

Cloud F の速度成分• N2H+は1成分，C18Oは4成分• N2H+，C18Oの間で同じ成分FPPV4に速度差→ <VN2H+-Vc18o> = +0.32 ± 0.06 km/s
• FPPV1とFPPV2の速度差2.9±0.5 km/s

Cloud H (IRDC: G035.39 – 00.33) との比較

• Cloud H にも同じ解析，C18Oで5成分• いくつかの共通点• 同様の2.9±0.3 km/sの速度差をもつ2成分が存在• HPPV4でN2H+とC18Oの間に速度差 +0.26 ± 0.02 km/s
→ 同じ形成シナリオで説明できる？
• Cloud H はfilament merging senario?
• 高速度の成分 (HPPV1,2)と低速度の成分 (HPPV3)がmergingして中間速度の成分 (HPPV4)が形成• N2H+は高速度成分に押されたdenseな場所 → C18Oより高速度になる
1414Barnes, Henshaw, Caselli, Jiménez-Serra, Tan, Fontani, Pon, Ragan

Barnes, Henshaw, Caselli, Jiménez-Serra, Tan, Fontani, Pon, Ragan

GMC scales for Cloud F• large scale の観測と今回の観測がよく一致• Cloud F はより大きなGiant
Molecular Filamentの一部の可能性がある

8

#2概要

• 分子雲から円盤形成の数値シミュレーション• 初の原始星円盤のpopulation synthesis studyを行なった

計算について• 3D SPHコードでradiative hydrodynamicsの計算• 初期条件• 一様密度の分子雲• 質量: 500 Msun
• スケール: 0.404 pc• 密度: 1.2x10-19 g cm-3
• tff = 6.0x1012 s ~ 1.90x105 yr• T = 10.3 K• 乱流あり（rms: M=13.7）••••• 円盤の条件• r < 2000 au かつ軌道の離心率e < 0.3のSPH粒子• 次のものは除外- Mdisk < 0.003 Msun 以下，全円盤質量の63%を含む半径が300 au以上，63%を含む半径が50%を含む半径の3倍以上• 円盤半径- 全円盤質量の63.2%を含む半径

計算は1.20 tff ~ 0.23 Myr まで
最も進化した protostar は 0.9x105 yr
分解能は円盤サイズ 10 au
最終的に 84のsingle protostarと40のmultiple systems（28 binary, 5 triple, 7 quadruple）が形成

9

結果
多様な円盤
• misaligned systems, spiral, warped disks, etc.

Figure 12. The evolution of the disc radii of isolated protostars (sink particles) as a function of age. Each line represents the evolution of the discaround a particular protostar. Lines may stop and start, for example, if theprotostar becomes part of a multiple system, or is expelled from a multiple
system, respectively.

• 典型的に• M* ~ 0.02 - 0.4 Msun
• Mdisk ~ 0.003 - 0.5 Msun
• Rdisk ~ 20 – 100 au
Figure 13. The evolution of the ratio of the disc mass to protostellar (sinkparticle) mass as a function of age for isolated protostars. Each line represents the evolution of the disc around a particular protostar. Lines may stopand start, for example, if the protostar becomes part of a multiple system,
or is expelled from a multiple system, respectively.

Figure 11. The time evolution of the disc masses of isolated protostars (sinkparticles) during the calculation. The top panel gives the disc mass versus
the (linear) time in the calculation, while the bottom panel gives the discmass versus the age of the protostar (using a logarithmic time scale). For
comparison with the top panel, the middle panel provides the time evolution of the protostellar mass (i.e. sink particle mass) of the isolated protostars whose disc masses are given in the top panel. Each line represents theevolution of the disc around a particular isolated protostar. Lines may stopand start, for example, if the protostar becomes part of a multiple system,
or is expelled from a multiple system, respectively.
c 0000 RAS, MNRAS 000, 000–000

disk mass• 時間とともに増加• 典型的に103 yr – 104 yrで3倍程度まで成長• 104 yrより先はdisk massが減少• Md ∝ M*~0.85 (M* < 0.5 Msun)

Note that there are few isolated protostars older than 30,000 yrs
when the calculation is stopped. Also, after ⇡ 104 yrs, it is clearthat the disc masses around some protostars rapidly decline (dueto accretion, encounters with other protostars, and/or ram-pressurestripping). Some protostars also suddenly become isolated as they
are ejected from multiple systems and these usually have low discmasses.
In Fig. 12 we plot the disc radii of isolated protostars versus
age. As with disc mass, there is a general trend for disc radii to
get larger with time. They range from radii of 10–50 au at 103 yrs
old to 20–100 au at 104 yrs old. Note that even though the calculation does not treat magnetic fields (which could provide angularmomentum transport by magnetic braking), the disc radii are notunusually large. We will discuss this further in Section 6.
In Fig. 13 we plot the ratio of the disc mass to the protostellar

Diversity and properties of protostellar discs

21

Figure 19. Distributions of the relative orientation angle between the two circumstellar discs in bound protostellar pairs. Pairs include both binaries and bound11

観測との比較

dust massvs Class II disks• disk mass: ~ 104 yrの若いdiskの方がClass IIのdiskに比べて30 – 300倍重い• disk radius: 分布は観測とシミュレーションでよく一致

vs Class 0 disks

gas mass• disk mass: 観測から0.05 – 0.2 Msun (Jorgensen+9,Enoch+11), 0.09 - 0.36 Msun (Segura-Cox+16)• シミュレーション（0.02 - 0.4 Msun）とコンシステント• disk radius も比較的コンシステント• disk mass は disk radiusにはあまり依らない• Class 0からClass IIへはdisk radiusを保ったまま質量のみを減らして進化する

12

#301 50

13

52

14

#302 315

416

517