Baba et al. (2017) Data

2019/05/02 作成

2020/03/12 追記

はじめに

  • 具体的な研究課題の共同研究のために限定公開します。→ 馬場までお問い合わせください

  • データの無断二次配布・転用はご遠慮ください (ご相談いただければ対応します)。

  • このデータを用いた論文を執筆する際には、以下の論文を引用して頂くようお願い致します。

  • 本研究は「若手研究 B (代表・馬場淳一; 研究課題番号: 26800099)」の支援のもと実施したものです。

シミュレーションの概要 (Baba et al. 2017)

  • 動的渦状腕モデル (e.g. Baba et al. 2013) における多相星間ガスの数値シミュレーション(ただし、棒状構造あり)

  • Baba et al. 2017で解析している時刻の周辺 200 Myr程度 (0.25 Myr刻み) のデータあり

  • コード:ASURA-2 (Saitoh et al., 2008; Saitoh & Makino 2009, 2010)

  • 計算機:Cray XC30 (アテルイ; 国立天文台 CfCA)

  • 計算時間:約 3ヶ月程度 (512コア)

  • 粒子数:~80 millions

  • 粒子質量:~650 Msun

  • 重力ソフトニング長:3 pc

  • 物理過程:恒星系力学、星間ガスの自己重力、星形成、放射冷却 (20 K < T < 10^8 K)、FUV加熱、Hii領域加熱、超新星爆発加熱、H2生成

  • 映像 (4D2U):「天の川銀河紀行」[Youtube]

データ形式1(2019/05/02)

  • ガス粒子データ (2.5 GB)

    • x, y, z [kpc], vx, vx, vz [km/s], Metallicity, m [Msun], ax, ay, az [km/s^2], rho [H/cc], T [K]

    • http://****/Baba2017/MilkyWayGasOffset090Ro08.0PA23.0FileID10501.0000.0001.dat

  • 星粒子データ (6.2 GB)

    • x, y, z [kpc], vx, vx, vz [km/s], Metallicity, age [Gyr], m [Msun], ax, ay, az [km/s^2]

    • http://****/Baba2017/MilkyWayStarsOffset090Ro08.0PA23.0FileID10501.0000.0001.dat

  • 太陽系を (x,y) = (0,-8) kpc付近におくと、バーの角度が25度となるように回転変換済み

データ形式2 (2020/03/12):Synthetic observation of the Galactic-scale multiphase ISM

  • ガス粒子データ (~6 GB)

    • http://****/PostProcessRe_rot_light.dat (FileID = 10501)

    • http://****/read.c (データ読み込みプログラム)

      • % gcc -lm read.c -o read

      • % ./read

  • 太陽系を (x,y) = (0,-8) kpc付近におくと、バーの角度が25度となるように回転変換済み

  • Baba et al. (2017) データの星間環境 (密度、温度、FUV強度など) に対して、星間化学反応 (Wolfire et al. 1995; Koyama & Inutsuka 2000) を計算し、H2量、CO量 (化学平衡) などをポストプロセス計算

    • すべてのSPH粒子に対して solar metallicity = 0.019 を仮定 (重元素量勾配を考慮していない)

    • 各SPH粒子の面密度 Nt = nH × 100 pc を仮定

    • G0:各SPH粒子のG0値 (Baba et al. 2017の計算法で評価; 最小値 0.001を設定)

    • X0 = 1.0

  • 各SPH粒子に対して観測量を簡易計算(注:詳細な輻射輸送計算ではない)

    • 微細構造線輝線強度:[CII]158um, [OI]63um (optically thin、統計平衡を仮定)

    • 中性水素21cm輝線:4πJν [erg/s/cc]、α [1/cm]

    • CO輝線:RADEX (LVG近似) を利用

      • H2 orth/para = 1.0

      • Tbg = 0.0 (i.e. mean Galactic background radiation fieldに対応)

      • deltav = 1.0 km/s

  • l-v図の例