申請書
Memo
ストーリー
高密度トレーサーを見る動機：
複数の速度成分を持つ（＝衝突の兆候がある）分子雲について、より高密度（＝星形成に密接にかかわる）領域について観測し解析したい
みたいな感じ？
複数の速度成分を持つ（＝衝突の兆候がある）の部分を、
- 分子雲衝突が階層的に起こる（今の自分の解析結果・解釈）
- 既存の分子雲衝突の研究
の二つを合わせて深堀したい
既存の分子雲衝突の解析のダメなところを挙げる （星形成領域近傍で衝突の傾向がみられるのに、運動学的には衝突頻度はそこまで多くない）
→ 階層構造だったら（dendro解析だったら）それが解決できる （高密度ガスの周囲に低密度ガスの階層構造、低密度ガスの収縮により内部高密度ガスの衝突しやすい環境が作られる）
→ COのdendro解析でもまだわからないことを挙げる （階層的衝突の兆候が見えるエリアで本当にmassive star形成につながるような高密度領域が形成されているか？ or 複数の速度成分があり衝突が示唆されているが、COのdendro解析ではそれが見えない）
→ 高密度トレーサー観測のdendro解析が必要
HCNとHCO+の有意性（ChatGPT）
① 密度感度のちょうどよい重なりと差異（臨界密度の特性）
HCNとHCO⁺は高密度トレーサーの中でも「密度感度が似ていながら、明確に異なる」珍しい組み合わせです。
- HCN: 臨界密度約
- HCO⁺: 臨界密度約
この微妙な臨界密度の差異がポイントです。
両者は密度に敏感な範囲が部分的に重なりますが、特に高密度（星形成コアやクランプの内部）ではHCNが卓越し、それに比べやや低密度（クランプの周辺部やフィラメントとの境界層）ではHCO⁺が支配的になります。
他の高密度トレーサー（例：N₂H⁺, CS, SiO）は、 - N₂H⁺はCO凍結の影響が大きすぎる（低温ガスに特化しすぎ）。
- CSは広すぎる密度範囲をトレースし、星形成活性度を直接区別しにくい。
- SiOは主に衝撃波のみに限定され、化学進化を全般的に捉えにくい。
一方、HCNとHCO⁺は密度トレーサーとしての微妙な違いに加え、「星形成領域における化学的進化と密度構造の同時解析」に非常にバランスが良いトレーサーなのです。
② 化学的安定性と感度のバランスの良さ（化学進化を診断できる）
HCNとHCO⁺は、
- 安定した豊富な分子で、広範な天体・領域で検出可能。
- 星形成の進行による化学進化に対し感度が高く、若い星形成コアからやや進化した活発な星形成領域まで一貫して観測できる。
- 両者は星形成に伴うガス加熱、紫外線照射、乱流運動の影響を明確に受けるため、輝線比が化学的進化やダイナミクスを強く反映します。
他の分子輝線（例えばN₂H⁺/CSなど）では、低温環境や凍結環境に過敏だったり、あるいは星形成の特定ステージに偏りが強すぎるため、幅広い環境や進化段階を横断的に比較するのに難があります。
一方、HCNとHCO⁺なら、「星形成初期段階〜活発に進行中の領域」までを一貫して、比較的安定した感度で観測可能です。
③ 星形成活動を診断するための実証的裏付けの豊富さ
HCN/HCO⁺輝線比は星形成天文学において広く研究され、観測的および理論的研究によって次のことが示されています。
- 活発な星形成領域ほどHCNが顕著に強まり、HCN/HCO⁺比が高くなる（）。
- 比較的若く活動が弱い領域では、HCO⁺が優勢でHCN/HCO⁺比が低く（<1）なる。
- この傾向が銀河系内から銀河系外の星形成領域に至るまで多くの観測で確認されており、信頼できる指標
これに対し、他の輝線比（例えばCS/N₂H⁺やCS/HCO⁺）は、天体間でのばらつきが大きく、輝線比の意味を一義的に解釈するのが困難になるケースがあります。
HCNとHCO⁺は輝線強度と星形成指標（赤外線放射強度、YSOの存在、アウトフロー）の間で明確な経験則的相関が確立されているため、観測データから星形成状況を具体的に診断する強力な根拠を提供します。
④ 実際のデータ解析での具体的優位性
FUGIN 観測データから導出した密度構造（フィラメント、クランプ候補）をさらに詳しく分析する場合、HCNとHCO⁺の観測が特に相性が良い理由は以下の通りです。
- データで得られる密度（）よりも一段階高密度をトレースするトレーサーとして、HCN/HCO⁺が密度構造の階層的な関係（dense core, clump, filament間の密度差）を捉えやすい。
- HCNとHCO⁺の輝線強度や速度構造は比較的強く明瞭で、観測効率が良い（短時間観測で良質なデータを取得可能）。
- 輝線比解析を行う際に必要なデータ精度や感度の要求が、他の高分子トレーサーよりも緩やかであり、観測の現実性（プロポーザル成功率）も高まります。
観測対象選定
G18
HCNもHCO+も未観測、近傍のSNRはHCNでもHCO+でも調べられている（Shen+2025、Paron+2012）
複数の速度成分に乗るleafがtrunkに囲まれる理想的な解析結果
これを基準に以下のように選定（ChatGPT）
【本命候補】W51
- 複数の速度成分それぞれに高密度構造
- 成分間の“またがり構造”が存在 → CCCの階層性が明瞭
- 全体を囲うtrunk構造がある → 全球収縮との整合
📌 階層的衝突＋massive core形成の描像にもっとも合致
➡︎ 本命ターゲットとして堂々と据えられる
【優秀対抗】Sh2-48
- G18と同様の構造 → すでにプロトタイプとして分析済
- 構造の入り組みが適度で、視認性が高い
📌 G18に準じた観測設計が転用可能
➡︎ 第二ターゲットまたは検証対象群として非常に良い
【潜在本命】W33, W39
- 複数速度成分 × 各に階層構造あり
- 視線方向の重なり or 広がり次第で階層性の理解が深まる
- W33は過去観測も多く補助データあり（※今回は除外）
📌 動的階層進化の別タイプとしての比較対象に最適
【参考対照】G45, N14, W43, W49
- 「10 km/s またがる構造」「細長い trunk」など
- 速度幅と密度構造が合致 → 収縮系構造の典型
📌 “衝突ではなく収縮由来のdense core”という対照として使える
【やや弱い】M16, M17, N4, G47
- 単一成分上に構造 → CCC的な解釈にはやや弱い
- 階層構造の“衝突的含意”を出しにくい
📌 重力収縮のみの階層進化例としては使えるが、仮説検証力は落ちる
観測パラメータ設定
観測パラメータの決定ガイド
- Receiver
– FOREST（4ビーム群）を選ぶのが一般的。単一ビーム → 感度・速度範囲の選定で不利。 - Polarization
– “Double” を選ぶと両偏波の合成で S/N が改善されます（S/N ≈ √2UP）。 - Frequency (GHz)
– 狙う分子線の周波数（例えば HCN/J=1–0 → ≈ 88.63 GHz, HCO⁺ → ≈ 89.19 GHz）。ドップラー補正で回転系速度を考慮。 - Maximum Elevation (deg)
– 観測対象の天体の最大地平高度。天体の赤緯・赤経から観測地（野辺山）での最大高度を計算。80°超の場合は80°と入力。 - Length of mapping area (scan)
– 観測範囲を arcsec 単位で入力。FORESTの場合はプロポーショナルに＋50″。たとえば、中心領域10′＝600″を狙うなら、
→ スキャン方向：650″、垂直方向：650″ - Length of mapping area (par.)
上参照 - Time for one scan (sec)
– 通常20 sec/sweep を推奨。総面を網羅しながら、時間も延びすぎないバランス。 - Dump Time (sec)
– “0.04 s” を選べばサンプリングが密になりノイズも改善される。過剰すぎる場合は“0.1 s”でもOK。 - Number of ONs per OFF
– 典型的には”2”。ベースライン安定性とスキャン速度のバランスを見て選択。 - Separation between the scans (arcsec)
– ビームの HPBW(≈17″) の 1/3：6″程度を推奨。Nyquist + オーバーサンプリング。 - Map grid (arcsec)
– スキャン間隔以上、HPBW/2 以下。例えば 6″（スキャン間隔と合致） - OFF-point separation (arcmin)
– 対象マップから OFF 位置までの距離。例えば、8′ 以上離すと対象ガスライン干渉回避しやすい。 - Tsys (K)
– ステータスレポート（冬期）から OBS の周波数帯・偏波状態で目安値を読みとる。
– たとえば90 GHz 帯では Tsys ≈ 150–250 K 程度が一般的（EL, 天候次第で変動）。 - Weather/external parameters（自動？）
– オフラインでは考慮不要。サイト側で自動取込されることが多い。 - Overhead factors
– Tsys 補正・校正・位置切替にかかる時間を含め、≒1.2〜1.3×を乗じるのが実際的。
Time estimates
G18 HCN
* Parameters :
1.Receiver = FOREST
2.Pol = Double
3.Frequency = 88.63 GHz
4.Max Elevation [deg] = 60
5.Length along the scans [arcsec] = 2880
6.Length perp. to the scans [arcsec] = 3000
7.Time for scan [sec] = 20
8.Dump time [sec] = 0.04
9.Num. of ONs per OFF = 2
10.Separation between scans [arcsec] = 6
11.Map grid [arcsec] = 6
12.OFF-point separation [arcmin] = 8
13.Tsys_inp [K] = 200.0
14.Resolution
- Frequency resolution (kHz) = 244.14
- Velocity resolution (km/s) =
15.Convolution function = Bessel*Gauss* Results :
1.theta [deg] = 0.0
2.Beam overlapped area = 2830.0 ["] * 2950.0 ["]
3.vscan is 144.0 ["/sec] (5.8 per 0.04s sample)
4.Nrow = 501
5.tapp [sec] = 10, ttran [sec] = 10, ttranoff [sec] = 8
6.tOFF = 3.9 [sec] -> 4 [sec]
7.tOH = 23.0 [sec]
8.ttot(ON) [min] = 167.0
9.ttot(OBS) [min]/[hour] = 400.8 / 6.68
10.eta(ON/OBS) = 0.42
11.tcell(ON) [sec] = 1.43
12.tcell(OFF) [sec] = 32.00
13.Tsys_inp [K] = 200.0
14.Tsys_obs_sm[K] = 223.7
15.dTa*_sm [K] = 0.4393
16.Version = 2.190718G18 HCO+
* Parameters :
1.Receiver = FOREST
2.Pol = Double
3.Frequency = 89.19 GHz
4.Max Elevation [deg] = 60
5.Length along the scans [arcsec] = 2880
6.Length perp. to the scans [arcsec] = 3000
7.Time for scan [sec] = 20
8.Dump time [sec] = 0.04
9.Num. of ONs per OFF = 2
10.Separation between scans [arcsec] = 6
11.Map grid [arcsec] = 6
12.OFF-point separation [arcmin] = 8
13.Tsys_inp [K] = 200.0
14.Resolution
- Frequency resolution (kHz) = 244.14
- Velocity resolution (km/s) =
15.Convolution function = Bessel*Gauss* Results :
1.theta [deg] = 0.0
2.Beam overlapped area = 2830.0 ["] * 2950.0 ["]
3.vscan is 144.0 ["/sec] (5.8 per 0.04s sample)
4.Nrow = 501
5.tapp [sec] = 10, ttran [sec] = 10, ttranoff [sec] = 8
6.tOFF = 3.9 [sec] -> 4 [sec]
7.tOH = 23.0 [sec]
8.ttot(ON) [min] = 167.0
9.ttot(OBS) [min]/[hour] = 400.8 / 6.68
10.eta(ON/OBS) = 0.42
11.tcell(ON) [sec] = 1.43
12.tcell(OFF) [sec] = 32.00
13.Tsys_inp [K] = 200.0
14.Tsys_obs_sm[K] = 223.7
15.dTa*_sm [K] = 0.4393
16.Version = 2.190718W51 HCN
* Parameters :
1.Receiver = FOREST
2.Pol = Double
3.Frequency = 88.63 GHz
4.Max Elevation [deg] = 50
5.Length along the scans [arcsec] = 3600
6.Length perp. to the scans [arcsec] = 3650
7.Time for scan [sec] = 20
8.Dump time [sec] = 0.04
9.Num. of ONs per OFF = 2
10.Separation between scans [arcsec] = 6
11.Map grid [arcsec] = 6
12.OFF-point separation [arcmin] = 8
13.Tsys_inp [K] = 200.0
14.Resolution
- Frequency resolution (kHz) = 244.14
- Velocity resolution (km/s) =
15.Convolution function = Bessel*Gauss* Results :
1.theta [deg] = 0.0
2.Beam overlapped area = 3550.0 ["] * 3600.0 ["]
3.vscan is 180.0 ["/sec] (7.2 per 0.04s sample)
4.Nrow = 610
5.tapp [sec] = 10, ttran [sec] = 8, ttranoff [sec] = 8
6.tOFF = 3.5 [sec] -> 3 [sec]
7.tOH = 22.0 [sec]
8.ttot(ON) [min] = 203.3
9.ttot(OBS) [min]/[hour] = 471.7 / 7.86
10.eta(ON/OBS) = 0.43
11.tcell(ON) [sec] = 1.15
12.tcell(OFF) [sec] = 24.00
13.Tsys_inp [K] = 200.0
14.Tsys_obs_sm[K] = 229.9
15.dTa*_sm [K] = 0.5054
16.Version = 2.190718W51 HCO+
* Parameters :
1.Receiver = FOREST
2.Pol = Double
3.Frequency = 89.19 GHz
4.Max Elevation [deg] = 50
5.Length along the scans [arcsec] = 3600
6.Length perp. to the scans [arcsec] = 3650
7.Time for scan [sec] = 20
8.Dump time [sec] = 0.04
9.Num. of ONs per OFF = 2
10.Separation between scans [arcsec] = 6
11.Map grid [arcsec] = 6
12.OFF-point separation [arcmin] = 8
13.Tsys_inp [K] = 200.0
14.Resolution
- Frequency resolution (kHz) = 244.14
- Velocity resolution (km/s) =
15.Convolution function = Bessel*Gauss* Results :
1.theta [deg] = 0.0
2.Beam overlapped area = 3550.0 ["] * 3600.0 ["]
3.vscan is 180.0 ["/sec] (7.2 per 0.04s sample)
4.Nrow = 610
5.tapp [sec] = 10, ttran [sec] = 8, ttranoff [sec] = 8
6.tOFF = 3.5 [sec] -> 3 [sec]
7.tOH = 22.0 [sec]
8.ttot(ON) [min] = 203.3
9.ttot(OBS) [min]/[hour] = 471.7 / 7.86
10.eta(ON/OBS) = 0.43
11.tcell(ON) [sec] = 1.15
12.tcell(OFF) [sec] = 24.00
13.Tsys_inp [K] = 200.0
14.Tsys_obs_sm[K] = 229.9
15.dTa*_sm [K] = 0.5054
16.Version = 2.190718Sh2-48 HCN
* Parameters :
1.Receiver = FOREST
2.Pol = Double
3.Frequency = 89.19 GHz
4.Max Elevation [deg] = 50
5.Length along the scans [arcsec] = 3600
6.Length perp. to the scans [arcsec] = 3650
7.Time for scan [sec] = 20
8.Dump time [sec] = 0.04
9.Num. of ONs per OFF = 2
10.Separation between scans [arcsec] = 6
11.Map grid [arcsec] = 6
12.OFF-point separation [arcmin] = 8
13.Tsys_inp [K] = 200.0
14.Resolution
- Frequency resolution (kHz) = 244.14
- Velocity resolution (km/s) =
15.Convolution function = Bessel*Gauss* Results :
1.theta [deg] = 0.0
2.Beam overlapped area = 2110.0 ["] * 1080.0 ["]
3.vscan is 108.0 ["/sec] (4.3 per 0.04s sample)
4.Nrow = 190
5.tapp [sec] = 9, ttran [sec] = 8, ttranoff [sec] = 8
6.tOFF = 4.4 [sec] -> 4 [sec]
7.tOH = 21.0 [sec]
8.ttot(ON) [min] = 63.3
9.ttot(OBS) [min]/[hour] = 145.2 / 2.42
10.eta(ON/OBS) = 0.44
11.tcell(ON) [sec] = 1.91
12.tcell(OFF) [sec] = 32.00
13.Tsys_inp [K] = 200.0
14.Tsys_obs_sm[K] = 223.7
15.dTa*_sm [K] = 0.3831
16.Version = 2.190718Sh2-48 HCO+
* Parameters :
1.Receiver = FOREST
2.Pol = Double
3.Frequency = 89.19 GHz
4.Max Elevation [deg] = 60
5.Length along the scans [arcsec] = 2160
6.Length perp. to the scans [arcsec] = 1130
7.Time for scan [sec] = 20
8.Dump time [sec] = 0.04
9.Num. of ONs per OFF = 2
10.Separation between scans [arcsec] = 6
11.Map grid [arcsec] = 6
12.OFF-point separation [arcmin] = 8
13.Tsys_inp [K] = 200.0
14.Resolution
- Frequency resolution (kHz) = 244.14
- Velocity resolution (km/s) =
15.Convolution function = Bessel*Gauss* Results :
1.theta [deg] = 0.0
2.Beam overlapped area = 2110.0 ["] * 1080.0 ["]
3.vscan is 108.0 ["/sec] (4.3 per 0.04s sample)
4.Nrow = 190
5.tapp [sec] = 9, ttran [sec] = 8, ttranoff [sec] = 8
6.tOFF = 4.4 [sec] -> 4 [sec]
7.tOH = 21.0 [sec]
8.ttot(ON) [min] = 63.3
9.ttot(OBS) [min]/[hour] = 145.2 / 2.42
10.eta(ON/OBS) = 0.44
11.tcell(ON) [sec] = 1.91
12.tcell(OFF) [sec] = 32.00
13.Tsys_inp [K] = 200.0
14.Tsys_obs_sm[K] = 223.7
15.dTa*_sm [K] = 0.3831
16.Version = 2.190718