25/01/07 進捗報告
さて、今日は昨日発覚したプログラムのミス、つまり適切に(MWの)サンプルセレクションをすること。
本当に、Massive mass halo はAGNによってDenistyを上げているのか?(MWの場合と逆になっているのか?)を調べる。
ToDo;
- Massive mass haloのAGN vs SNのエネルギーを調べる
- AccretionRateを調べる
- ~~適切なサンプルセレクションに変更する~~
目次
AGN vs SN Energy
\[
\dot{E}_{\text{AGN}} = \varepsilon_{\text{FB}} \varepsilon_r \dot{M} \cdot c^2 = \left(\frac{\varepsilon_{\text{FB}}}{0.15}\right) \left(\frac{\varepsilon_r}{0.1}\right) \dot{M} c^2
\]
\[
\dot{E}_{\text{SN}} = \eta_{\text{SN}} \cdot \text{SFR} = \left(\frac{\eta_{\text{SN}}}{7e48 \, \text{erg}/M_{\odot}}\right) \cdot \text{SFR}
\]
を用いて計算する。SFRに関しては、CROCODILEでは Msun/yr で出てくる。
accretion rateである M dot は
UnitLength_in_cm = 3.08568e24
UnitMass_in_g = 1.989e43
UnitVelocity_in_cm_per_s = 100000
UnitTime_in_s = UnitLength_in_cm / UnitVelocity_in_cm_per_s
kei = 1.989e33 / 3.15576e7
aa = UnitMass_in_g / UnitTime_in_s / kei
aa
aa をかけてあげれば Msun/yr になる。c^2に関しては
8.987551787e+20 # erg/g
で良い。単位合わせで \(\eta_{\text{SN}}\) を1.989e+33 g/Msun で割る。
結果

redshift 4 から急激に上昇してSuperNovaのエネルギーの10倍近いエネルギーをfeedbackしていることがわかる。
SFR

Accretion Rate

BH mass

Massive mass galaxiesのAGN feedbackについて
AGN feedbackは明らかにMWよりも大きくなっているといえる。
で、あればDensityは下げられるはず?温度分布も見てみる。
Radial Profile Temperature

ビリアル半径/Halo Massとかに影響があったのかな…?

ほとんどなさそう。