BC03 確認のための計算 (3) - 天文学的研究メモ

Calzetti (2000)のダスト吸収がちゃんと計算されているかどうかの確認。
たとえば、5450Åでは

$k'(\lambda) = 2.659 \left( -2.156 + \frac{1.509}{\lambda} - \frac{0.198}{\lambda^2} + \frac{0.011}{\lambda^3} \right) + R_V' \approx 4.0876$

google先生に投げた計算式：
2.659 * (-2.156 + (1.509/0.545) - (0.198/0.545^2) + (0.011/0.545^3) ) + 4.05

ですから、たとえばE(B-V) = 0.5のとき、吸収のファクターは

$-0.4 E_s(B-V) k'(\lambda) \approx - 0.3597$

google先生に投げた計算式：
(- 0.4) * 0.44 * 0.5 * 4.0876

で、実際に5450*4.1Åでの計算の出力を確認したところ、たしかにE(B-V) = 0からE(B-V) = 0.5にすることで上記の吸収が起こってました*1。

*1:この確認をしたおかげで、関数内の条件分岐が間違っていたことを発見できました。