コンテンツにスキップ

2024年度

25/01/07 進捗報告

さて、今日は昨日発覚したプログラムのミス、つまり適切に(MWの)サンプルセレクションをすること。

本当に、Massive mass halo はAGNによってDenistyを上げているのか?(MWの場合と逆になっているのか?)を調べる。

ToDo;

  1. Massive mass haloのAGN vs SNのエネルギーを調べる ✅
  2. AccretionRateを調べる ✅
  3. ~~適切なサンプルセレクションに変更する~~ ❌

目次

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 で割る。

結果

image-20250107195114909

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

SFR

image-20250107200335648

Accretion Rate

image-20250107200709384

BH mass

image-20250107200917995

Massive mass galaxiesのAGN feedbackについて

AGN feedbackは明らかにMWよりも大きくなっているといえる。

で、あればDensityは下げられるはず?温度分布も見てみる。

Radial Profile Temperature

image-20250107202654127

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

image-20250107203222446

ほとんどなさそう。