すばる中間赤外線装置COMICSでさぐる大質量星形成
岡本 美子 （北里大学）
１．イントロダクション
２．超コンパクト(UC)HII領域の電離星
３．超コンパクト(UC)HII領域の電離ガス運動

イントロ：
大質量星形成研究の意義• 星形成としての興味– 小中質量星とはどのように異なるのか– 星形成の各段階の対応天体とタイムスケール– 初期状態と途中のプロセス• スタートの質量、途中の質量降着・放出（アウトフロー）、合体の可能性、周辺物質との相互作用 など• クラスター形成と孤立星形成との違い– 環境による大質量星形成の違い• Metallicity、分子雲質量、外圧などの影響• 系外銀河における星形成の理解の基礎– 銀河の活動性を担う1

• 大質量星形成に関係する天体– Cold cores– Warm cores– Hot cores• 0.1pc, 100K, 105cm-3, NH3etc– Hypercompact HII regions– Ultracompact HII regions←Protostarの形成？
←主系列星の形成？Ｓ5• 0.1pc, ne>~ 105cm-3, EM>~ 106cm-6 pc, 10 yr– OB stars/associations• 広いmass range:– ３桁のluminosity range, 6桁の電離光子数range– 質量と形成ステージがカップル• 込合った領域→ＳＥＤを天体ごとに分離できていない

Cloudの性質
ｖｓ形成される星• 割とよく性質がわかっているのはUCHIIsとhot cores– わかっているのは全体量• 温度、密度、電離フラックス、cloud質量など• いったいそこから・そこにどのくらいの質量の星が生まれる（生まれている）のかはよくわかっていない• Cloudの進化を考える上で重要• 若い段階から追うのは困難

Kim & Koo 2001?

– →UCHII領域段階で見てみる
Hot cores/UCHIIRs2

すばる／COMICS• Nバンド(8-13µm),Qバンド(16-25µm)• 観測モード：– 撮像(30”x40”),– ロングスリット分光(40”; R250/2500/10000)• 高分散分光は特定輝線波長のみ• 感度（1σ1s)：– 撮像 30mJy (N), 400mJy (Q)– N分光 340mJy(R250), 2.5Jy(R2500), 7Jy(R10000)– Q分光 3Jy(R2500), 7.5Jy(R5000)• 空間分解能– 8.2mの回折限界 (0.3”@10µm, 0.7”@20µm)

１．超コンパクト質量HII領域の電離星

3UCHII領域の大質量星の質量分布• 電離星の質量と数の決定• ２つの情報から– トータル電離フラックス←電波
連続波– 電離放射場の硬さ←金属イオ
ンMIR禁制線• [NeII]12.8µm 22eV
• [ArIII]8.99µm 28eV
• [SIV]10.51µm 35eV

– モデルとかなり異なる– モデルによってもかなり異なる– クラスター効果– Metallicityの影響• 系外銀河の埋もれたＨＩＩ領域にも応用可能

Okamoto et al. 2001

すばるによるMultiplicityの直接検証
J,H,K(Subaru/CISCO)

K3-50A
8.8,9.7,11.7(Subaru/COMICS)

５”~0.2pc4

Multiplicityの直接検証• Subaru+COMICS 0.4” resolution• K3-50A （Okamoto et al. 2003)– O5.5→O6-O8+O9.5[ArIII] 8.99μm[SIV] 10.51μm[NeII] 12.81μm

２．超コンパクト質量HII領域の電離ガスの運動＊＊一部、北里学術奨励研究助成（第16期）の助成による

5

形成シナリオ：降着 ｖｓ 合体
 形成過程について、小質量星のような確立したパラダイ

ムがない• 放射圧が降着を妨げる• 降着率が高ければＯＫ？
– >10-3Msun/yr, Dust opacityにも依存• 連星率・クラスタ密度との関係– オリオンでは大質量星の方が小質量星よりも連星率が大• 降着円盤に関しては最近Ｂ型星のまわりに発見したとする例がいくつか出てきた（IRAS20126+4104；Zhang et al. 1998など）– サイズ～10000AU• 形成中・直後の星周構造を観測的に明らかにする必要がある– 大質量星形成に関連したどのような星周構造が形成されるのか？
（Disk? Outflow?)

COMICSによるUCHII領域
電離ガス運動の観測– MIRイオン輝線の高分散分光
観測R10000• [NeII]12.8µm (22eV), [ArIII]8.99µm (28eV), [SIV] 10.5µm(35eV)• RRLよりも高い電離エネルギー

Thermal broadning + broadning
due to electron impacts– →より星近傍をトレース可能• 30km/sの分解能だがプロファイル
フィットにより中心速度はa few
km/sで決定可能– UCHIIRｓの線幅は~20-60km/s• ロングスリット分光→速度プロファ
イルの空間分布
Garay & Lizano 1999

6

G45.12+0.13の場合

Fig4 13CO J=1→0 outflow (Hunteret al. 1997)▲: Radio peak of UCHII region.
×: Suggested origin of the 2nd flow

• Cometary UCHII (WC89a)

Outflow2– ただし、multipleに近い• 9.5kpc → 1”~0.046pc• かなりmassiveな星の生成– Log(Nc') 49.53 s-1--> O5– IRAS FIR →O5 (Cluster仮定)

赤青• ne: 1.13x105cm-3 (peak)• EM: 1.5x109 cm-6pc (peak)• Ltotal=1.7 x106Lsun• 2’scaleのCOアウトフロー
Outflow1

G45.12+0.13のslit scanによる[NeII]輝線チャンネルマップ
VLSR=-42km/s

59km/s-15km/s

79km/s-12km/s

106km/s

39km/s

これまでに観測した４UCHIIｓのうちで最も線幅大7速度幅のアウトフロー方向へのenhancement
１成分gaussianでフィットしたときの[NeII]輝線幅(FWHM)
CO outflow (赤)CO outflow (青)物理量：アウトフローｖｓ電離ガス• アウトフロー駆動とUCHII現象は段階が異なる現象– エネルギー規模,年齢に大きな隔たり• UCHII形成時点でoutflowで生じた構造が残存– UCHIIのradio/MIR/NIR放射が、outflow軸に関連したようなアーク構造を持っている ＆ 電離ガス運動のoutflowとの関連

Energy [erg]Outflow (HunterUCHII region

et al.1997)(from ∆v of [NeII])

5.6x1048

1.9x1045(1/2 Mionized gas ∆v2)

Moment[Msun km/s]

31000

τdynamical [yr]2x10550(Mionized gas ∆v)

6x103(UCHII size / ∆v)

電離ガス質量～1.3Msun：電子密度とサイズからすべてH+として求めた

8

複数のアーク構造• アウトフローに沿ったような複数のアーク構造が、Nband域
各波長、[NeII]放射、NIR（J,H,K)放射で観測される
７“～0.3pc=60000AU

MIR/NIR/電波観測 1.過去のアウトフロー活動で、から考えられる中心 アウトフローローブの方向に周辺よりも密度の低い部の構造領域を形成
アウトフロー赤ローブ
MIRアークMIRアーク (warm dustemission）
emission）
低密度領域
へ膨張する
電離ガス
過去のアウトフロー活動に関連した高密度領域
NIR アーク（hotアーク（hot dustemission, reflection,or ionized gas)2.大質量星からのUV放射開始後、電離ガスが低密度領域へ膨張。運動はシャンパン流的になっており、この低密度領域で[NeII]速度幅が増加。はき寄せられた物質がMIR/NIRアークとして見えている大質量星+星周高密度物質？

→O型星でも小質量星的な降着による形成？

アウトフロー青ローブ

9

Simpleな形状を持つUCHII：• G35.20-1.74 (3.3kpc(Araya et al. 2002))– 典型的なcometary UCHIIR– 2slit位置 (cometaryの軸に並行と垂直)– 分子線LSR速度：中心は42-46km/s、ウイング込 37-51km/s– ＲＲＬLSR速度47-48km/s
G35 2cm 10”x10”– Radio continuum → O7.5相当• G33.92+0.11 (7.1kpc (Fish et al 2003))– Radioではcore-haloだがNバンドではbipolar– 1slit位置 (bipolarの軸)– 分子線LSR速度
106-108km/s– ＲＲＬはLSR速度
98.7-103.3km/s– Radio continuum → O7.5相当• これらの[NeII]ガス運動をCOMICS高分散
分光で観測

G33 2cm 10”x10”
Wood & Churchwell 1989

G35.20-1.74[NeII]PV図OFFSET [arcsec]

10G35.20-1.74
シャンパン流的cometary HII領域としてよく説明できる• 「NeII]輝線– V(LSR)が40-50km/s– cometaryの軸上ではtail側にかけて赤くなる– PV図でどのCutでも赤い側にウイング• 分子線 (Hofner et al. 2000等)– V(LSR)はおよそ42-45km/s– 青い側にウイング

350µmdust/分子UCHII
Tail方向に
ガス加速• 350µmダスト放射 (Hunter et al. 2000)– 北側と東側にバー状に分布– Cometaryのheadが北東• 速度範囲は20km/s程度→視線向こう側に開いたシャンパン流膨張

G33.92+0.11の場合：bipolar軸にそった[NeII]ガス加速• 分解能に近い線幅– プロファイルに1成分gaussian fit• 南北極方向への加速
輝度分布(+95)と中心速度11G33.92+0.11• [NeII]は分子線よりも青い– RRLも青いがRRLより分子との差が大(数km/sと15km/s)• 中心で一定に近く、その外で距離に対して線形に加速、さらに外では再び一定速度– ただし、全体として9km/s/pc程度の加速にのっているよう– 加速量は15km/s、率 125km/s/pc• ほぼedge-onで見たbipolar ionized outflowか？
– 北側が赤いので向こう側– 11.7µm像で北側が弱いのと矛盾しない• 11.7µmより8.8µmの方が中心の2ピークの距離が小– おそらくbipolarの中心部に電離星• 電離フローの開口角77度• NIRでは中心星は見えていない。
– 350mm peakにUCHIIがあるので矛盾はしない

構造減光体？

• 分子ガスよりも赤いほうが見えていないのは不思議– 青と赤の2成分系に単純化して考えた場合、赤側が見えなくなるのにおよそ10%程度までそちらが減光しているとするとこれはτにすると、2程度以上に相当– Avにすると～70mag、– N(H)=1.1e23 [cm-2]相当。少し大きすぎる？

観測• 向こう側への運動をさえぎるものがある？
• 輝線のτの問題？
[NeII]gas12

加速メカニズムは何か？
• 加速領域の大きさ・加速が線形であること・加速度・速度分布(中心速度、速度幅)が重要だろう• 考えられるメカニズム• 降着円盤による磁気圧：開口角が大きすぎでは• 放射圧– この場合、点源放射源による放射圧加速は、点源から遠くなるほど光子との衝突頻度が落ちるため、加速度が小さくなりそう• 星風：もとをたどれば放射圧考えられる加速メカニズム• 圧力勾配（シャンパン流）– 数十km/sまでの速度を実現可能。
– 速度プロファイルも高密度物質塊の切れ目付近で距離によっておよそ線形の速度勾配を実現可能。その外では加速がゆるくなり、やがて一定速度になる。G33で観測された速度勾配によく似ている。
– 速度分布から考えて、中心に東西に伸びたような高密度塊が必要• 光蒸発ディスク– disk状の減光は検出できていないので、あったとしても非常に小さなディスク– 速度と中心星質量から考えて～100AU半径のdiskからのフローというのは考え得る– ただし、Hollenbach et al. 1994のモデル論文には速度場の情報がない。

→disk的な平坦高密度構造がbipolarに垂直な方向

13シャンパン流の場合
Tenorio-Tagle 1979

30km/s初期密度勾配の境界

Bodenheimer et al. 1979

速度距離まとめ• UCHII領域の電離ガスの運動をCOMICSのNバンド高分散分光モードを用いて観測– [NeII]ガス運動の様子が見えてきた– G45.12+0.13：過去のアウトフロー活動に関係した構造– G35.20-1.70: わかりやすいシャンパン流UCHII領域。分子よりも赤い方(視線奥側）へ物質が流出。
– G33.92+0.11: MIRではbipolarなUCHIIに見え、特徴的な速度分布。
• 中心片側～0.1pcの範囲に渡って距離に対して線形の加速。
• 板状高密度ガス塊の存在可能性• 今後、電離エネルギーの違う輝線を使うと奥行き方向などをさらに分解して運動を調べられる– [ArIII]8.99µm (28eV), [SIV]10.51µm (35eV)

14