24/12/09 進捗報告
ブラックホールの成長・時間進化を見る。
まず、調査対象のsubhaloの100 kpc内に存在するBHのの一覧がこちらである：BHParticleIDData
（※ちなみにunitはMsunではなく、10^10 Msunである。）
注目してほしいのが、MWID:4と5、MWID:16と17である。
Center Positionが同じ値を出している。それによって、メインのブラックホールも同じものを捉えてしまっているのである。そこで、追跡するブラックホールのIDは以下のものとする。合計で61個のBHを追跡する。
249242294
169680463
265333883
172423807
246163043
170226224
155281906
208612453
223522582
157207430
156680727
146628891
213210494
159826993
266035094
169407686
160050077
214741514
245545475
188488601
135789381
219420739
150041765
199195569
168125978
163211497
254732697
233275827
136073089
199852754
184939576
212631350
220545004
138560963
257525875
161365927
162704975
171705571
183881157
154430399
149768689
220215486
188932950
156212985
143077627
221603749
173849085
208078675
176628935
212817526
164657014
193837684
260688852
241173556
253733277
260836379
157809534
220506187
146332175
234944912
257318483
BH accretion Rate evolution
BH mass growth
Re-check from \(\delta M / \delta t\)
ちなみにtime intervalの計算では次のコードを使用した。
import numpy as np
from scipy.integrate import quad
# 宇宙定数の設定
Om0 = 0.3 # 現在の物質密度パラメータ
Ok0 = 0.0 # 現在の曲率密度パラメータ
Ol0 = 0.7 # 現在の暗黒エネルギー密度パラメータ
g = 9.779 / 0.67742 # to Gyr
def E(z):
return np.sqrt(Om0 * (1 + z)**3 + Ok0 * (1 + z)**2 + Ol0)
def integrand(z):
return 1.0 / ((1 + z) * E(z))
def time_from_z(z):
integral, _ = quad(integrand, 0, z)
return integral * g
# 各 z での宇宙年齢を計算
ages = np.array([time_from_z(z) for z in redshift])
# 時間間隔を計算
time_intervals = np.diff(ages)
# 結果を表示
print("Time intervals between z values (Gyr):", time_intervals)