21cm線の観測理論(島袋先生のXの問題）

島袋隼人先生（天文学者）が今日ゼミでやった導出AIに聞きながら解いてみてと投げかけがあった。いろいろ詰まってるので取り掛かるのは後になるかとおもうけど、問題のみは調べたので書いておく。

「さらに表示」すると「研究しましょう！」となってます。
画像を書き写すと

$$
\frac{dI_\nu}{ds}=-\alpha_\nu I_\nu+j_\nu　\hspace{20mm}(A)
$$

これから

$$
\delta T_b\approx 27x_{HI}(1+\delta)\left(\frac{\Omega_b h^2}{0.023}\right)\left(\frac{0.15}{\Omega_m h^2}\frac{1+z}{10}\right)^{1/2}\left(1-\frac{T_{CMB}}{T_s}\right)\mathrm{mK}　\hspace{20mm}(B)
$$

を導けということです。
(A)は「放射輸送の式」というらしく（B)は21cm線輝度温度の式とのこと。

AIに聞いてみる。

まずは式に現れる項目の説明をお願いします。一見以下のように思えます
(A) s:エントロピー密度？他は不明
（B)$${T_b}$$:バリオン温度？
$${\Omega_b}$$バリオン密度？$${h^2}$$はきっとハッブル定数に関係がありそう
$${\Omega_m}$$物質の密度？バリオンと別にする意味は？
z赤方偏移、ここでは時刻を表している？
$${T_{CMB}}$$CMBの放射の温度
$${T_s}$$何の温度だろう？
単独の$${\delta}$$は何だろう
$${\mathrm{mK}}$$は単位でマイクロケルビン。きっとこれは確定

聞いた結果,色々間違っていた
ます放射輸送の式は

$$
\frac{dI_\nu}{ds}=-\alpha_\nu I_\nu+j_\nu　\hspace{20mm}(A再)
$$

$${I\nu}$$：比強度（specific intensity）単位：エネルギー / (時間・面積・周波数・立体角)
空間の単位面積に届くエネルギー量で、時間当たり周波数当り立体角当たりの単位面積を通るエネルギー量を示す。
エネルギーフラックス$${I_\nu=dE/dt\,dA\,d\nu\, d\Omega}$$
Aは面積、他は慣例的によく使われる奴
$${s}$$経路長　単位：1/長さ：視線方向の距離。なのでdsは光が単位長さ進むときという感じ。
$${\alpha_\nu}$$吸収係数　単位長さ当りの吸収の強さ。つまり強度が落ちる。この逆数が平均自由行程
$${j_\nu}$$放射係数　単位：エネルギー / (時間・体積・周波数・立体角)
単位体積が発する単位時間・単位立体角・単位周波数あたりの自発放射エネルギー。つまり強度が増える。
放射輸送の式は光が$${\Delta s}$$進むとき放射の比強度の変化は$${\Delta I_\nu= j_\nu \Delta s -\alpha_\nu I_\nu \Delta s}$$変化する。この式を$${\Delta s}$$で割って$${\Delta s \longrightarrow 0}$$とすると（A)が得られる
これは放射場の輸送方程式で光子数（生成-消滅）の連続の式を表している。（流体力学の講義大学時代工学部系にはあったような気がするけど理学部系にはなかったか、非常に難解そうな怪しい雰囲気があったので未履修。流体力学関係ありそう）
まあ、光のエネルギーが今どうなっているかの理屈を表している式だと思う。

21cm線輝度温度の式をもう一度書いておくと

$$
\delta T_b\approx 27x_{HI}(1+\delta)\left(\frac{\Omega_b h^2}{0.023}\right)\left(\frac{0.15}{\Omega_m h^2}\frac{1+z}{10}\right)^{1/2}\left(1-\frac{T_{CMB}}{T_s}\right)\mathrm{mK}　\hspace{20mm}(B再)
$$

$${T_b}$$輝度温度、比強度からレイリージーンズの近似によって温度に換算したもの
$${T_{CMB}}$$CMBの温度　電波の強度スペクトルを黒体放射（プランク分布）に当てはめて算出した温度
$${\delta T_b\equiv T_b-T_{CMB}}$$　これは定義でCMBの温度とスピン温度の差を観測するということですね。
$${x_{HI} = \frac{n_{HI}}{n_H} }$$ 水素の中性分率 全水素のうち電離していない割合（0〜1）
$${\delta \equiv \frac{\rho-\bar \rho}{\bar\rho} }$$　密度の平均からのずれの割合
$${\Omega_i}$$　密度パラメータ i=b:バリオン,m:バリオン+ダークマター=物質。
$${h}$$ハッブル定数の無次元化$${H_0=100h \mathrm{km/s/Mpc}}$$$${h \approx 0.67}$$
$${z}$$赤方偏移
$${T_s}$$スピン温度
$${\mathrm{mK}}$$単位、ミリケルビン

式を示して各項を説明して？と教えてもらったAIはChatGPT, gemini, cloude。チャッピーは三人の中では平均的、クローデさんは速いが少しいい加減（次元はしょったのがあって迷路に入った）、ジェミちゃん導出まで説明してくれた。気が利くがもっと掘らないとわからない。全部無課金。
密度パラメータなんて宇宙論パラメータが出てくるところが面白すぎる。
きょうはここまで。