偏極と音響振動

次に強結合近似の1次の方程式を考える．上で得られた0次の解を式 (12.6.285)-(12.6.294)の右辺に代入することにより，

$${\delta^{\rm (GI)}_\gamma}' - \frac43 k^2 v^{\rm (GI)}_\gamma + 4 {\mit\Psi}' = 0,$$ (L.6.298)

$$v^{\rm (GI)}_{\rm b} = v^{\rm (GI)}_\gamma + \tau_{\rm c} \left( {v^{\rm (GI)}_\gamma}' + \frac14 \delta^{\rm (GI)}_\gamma + {\mit\Phi} \right),$$ (L.6.299)

$${\mit\Theta}_2 = - \frac{8}{45} \tau_{\rm c} k^2 v^{\rm (GI)}_\gamma,$$ (L.6.300)

$${\mit\Theta}_l = 0, \quad (l \geq 3),$$ (L.6.301)

$${\mit\Theta}_{{\rm P} 0} = - \frac29 \tau_{\rm c} k^2 v^{\rm (GI)}_\gamma,$$ (L.6.302)

$${\mit\Theta}_{{\rm P} 1} = 0,$$ (L.6.303)

$${\mit\Theta}_{{\rm P} 2} = - \frac{2}{45} \tau_{\rm c} k^2 v^{\rm (GI)}_\gamma,$$ (L.6.304)

$${\mit\Theta}_{{\rm P} l} = 0, \quad (l \geq 3)$$ (L.6.305)

$${\delta_{\rm b}^{\rm (GI)}}' - k^2 v_{\rm b}^{\rm (GI)} + 3 {\mit\Psi}' = 0,$$ (L.6.306)

$$v^{\rm (GI)}_{\rm b} = v^{\rm (GI)}_\gamma - \tau_{\rm c} R \left( {v^{\rm (GI)}_\gamma}' + {\cal H}v^{\rm (GI)}_\gamma + {\mit\Phi} \right)$$ (L.6.307)

を得る．式(12.6.305)を見ると，強結合近似の1次の効果として４重極 成分 ${\mit\Theta}_2$ が発生していることがわかる．これは散乱と散乱の間に光子 が自由運動することで分布の双極子成分 $v_\gamma^{\rm (GI)} \propto {\mit\Theta}_1$ が４重極子成分を生み出すためである．そうして生み出された４重 極成分は散乱によりさらに偏極成分を生みだす．式(12.6.305), (12.6.307), (12.6.309)から

$${\mit\Theta}_{{\rm P} 2} = \frac15 {\mit\Theta}_{{\rm P} 0} = \frac14 {\mit\Theta}_2$$ (L.6.308)

という関係があるので，偏極は分布の４重極成分に比例し，

$${\mit\Theta}_{\rm P} = \frac54 \left[1 - P_2(\mu)\right] {\mit\Theta}_2$$ (L.6.309)

で与えられる．

こうして，光子の双極子成分が４重極成分と偏極成分を決定づけていることが わかる．ここでポテンシャルを外場とみなせば，双極子成分は式 (12.6.303) により単極子成分 $\delta^{\rm (GI)}_\gamma \propto {\mit\Theta}_0$ から求められる．そこで，単極子成分の発展方程式を求めると都合 がよいことがわかるだろう．他に単極子成分を含む方程式は式 (12.6.304)であるが，これはバリオンの速度場 $v^{\rm (GI)}_{\rm b}$ を含んでいるので，さらに式(12.6.312)を用いて $v^{\rm (GI)}_{\rm b}$ を消去すると，

$${v^{\rm (GI)}_\gamma}' + \frac{R'}{1 + R} v^{\rm (GI)}_\gamma + \frac14 \frac{1}{1 + R} \delta^{\rm (GI)}_\gamma = - {\mit\Phi}$$ (L.6.310)

が得られる．ただしここで，式(12.4.238)により示される式 $R' = {\cal H}R$ を用いた．この式と式(12.6.303)から $v^{\rm (GI)}_\gamma$ を消去すれば，単極子成分に関する発展方程式

$${\delta^{\rm (GI)}_\gamma}'' + \frac{R'}{1+R} {\delta^{\rm (GI)}_... ...amma = - \frac{4k^2}{3} {\mit\Phi}- \frac{4R'}{1+R} {\mit\Psi}' - 4{\mit\Psi}''$$ (L.6.311)

を得る． この方程式は左辺を強制項とし，まさつのある振動子を表す線形２階微分方程 式となっている．つまり，左辺で表される変動重力ポテンシャ ルの影響を受けつつ，光子のエネルギー密度が自己重力と圧力により振動する 音響振動の様子を記述していて、その角振動数は $c_{\rm s}k$ で与えられるこ とが読み取れる．

一般に，1変数の線形微分方程式では，強制項のない斉次方程式の解をすべて求 めることができれば，グリーン関数が構成でき，もとの方程式の一般解を求め ることができる．一般に強制項$F(\tau)$ を持つ２階微分方程式について，そ の斉次方程式の２つの独立解を$X_1(s)$ , $X_2(s)$ とすれば，一般解は斉次方 程式の一般解に非斉次方程式の特解を加えた形

$$X(\tau) = C_1 X_1(\tau) + C_2 X_2(\tau) + \int_{\tau_0}^\tau F(s)... ...c{X_1(s) X_2(\tau) - X_1(\tau) X_2(s)} {X_1(s) {X_2}'(s) - {X_1}'(s) X_2(s)} ds$$ (L.6.312)

で与えられる．ここで$\tau_0$ は初期値を与える点である．そこで，式 (12.6.316)の解を求めるにはその斉次方程式の一般解を求めればよい ことになる．

だが正確な一般解を求めることは難しいので，$R$ の時間変化よりも十分速く 振動する解を考えてWKB近似によって求めることにする．節 6.3.4と同様に計算することにより，斉次方程式の独立な2解と して，

$$\delta^{\rm (GI)}_\gamma = \left\{ \begin{array}{c} (1 + R)^{-1/4} \cos(k r_{\rm s}) (1 + R)^{-1/4} \sin(k r_{\rm s}) \end{array} \right.$$ (L.6.313)

が得られる．ただし，

$$r_{\rm s}(\tau) = \int_0^\tau c_{\rm s} \frac{d\tau}{c}$$ (L.6.314)

は音響ホライズン (sound horizon) ，す なわち音速が宇宙年齢の間に到達できる距離である．この量は式 (4.2.29)および式(12.6.296)の近似のもとでこれらの式，および $R = R_{\rm eq} a/a_{\rm eq}$ となることを使えば積分でき，

$$r_{\rm s}(\tau) = \frac{1}{\sqrt{3}} \int_0^\tau \frac{d\tau'}{\s... ...\left( \frac{\sqrt{1+R} + \sqrt{R + R_{\rm eq}}}{1 + \sqrt{R_{\rm eq}}} \right)$$ (L.6.315)

となる．ただしここで， $R_{\rm eq} \equiv R(\tau_{\rm eq})$ は等密度時の $R$ の値， $k_{\rm eq} \equiv a(\tau_{\rm eq}) H(\tau_{\rm eq})$ は等密 度時のホライズンスケールの波数である．WKB解(12.6.318)を式 (12.6.317)に用いれば式(12.6.316)の近似解が求まり，

$$\delta^{\rm (GI)}_\gamma(\tau) = \frac{4}{\left[1 + R(\tau)\right... ...C_1 \cos\left[k r_{\rm s}(\tau)\right] + C_2 \sin\left[k r_{\rm s}(\tau)\right]$$

$$\qquad\qquad\qquad\qquad\qquad - \frac{\sqrt{3}}{k} \int_0^\tau... ...+ \frac{R'(\tau')}{1+R(\tau')} {\mit\Psi}'(\tau') + {\mit\Psi}''(\tau') \right]$$

$$\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times \left[1 + ... ...\right]^{3/4} \sin\left[ k r_{\rm s}(\tau) - k r_{\rm s}(\tau')\right] \Biggr\}$$

となり，この式は近似を用いずに数値計算をして得られる振動の様子を非常に よく表すことが知られている^L4． ここで$C_1$ , $C_2$ は積分定数であり，上に求めた超ホライズンスケールの振 る舞いから決まる初期条件から定まる． このために，上の式を微分してみると，いま，WKB近似を用いているので，振動 部以外の因子として出てくる$R$ の微分は無視して，

$${\delta^{\rm (GI)}_\gamma}'(\tau) = - \frac{4k}{3\left[1 + R(\tau... ...C_1 \sin\left[k r_{\rm s}(\tau)\right] - C_2 \cos\left[k r_{\rm s}(\tau)\right]$$

$$\qquad\qquad\qquad\qquad\qquad + \frac{\sqrt{3}}{k} \int_0^\tau... ...+ \frac{R'(\tau')}{1+R(\tau')} {\mit\Psi}'(\tau') + {\mit\Psi}''(\tau') \right]$$

$$\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times \left[1 + ... ...\right]^{3/4} \cos\left[ k r_{\rm s}(\tau) - k r_{\rm s}(\tau')\right] \Biggr\}$$

となる．したがって，式(12.5.279), (12.5.283)により， 断熱ゆらぎに対し，

$$C_1 = -\frac12 {\mit\Phi}_0,\qquad C_2 = 0$$ (L.6.316)

また等曲率ゆらぎに対し，

$$C_1 = 0, \qquad C_2 = \frac{3\sqrt{2}k_{\rm H,eq}}{16k} \left(1 + \frac{2}{15}R_\nu\right)S_0$$ (L.6.317)

と係数が求まる．

この近似解(12.6.321)から、斉次解の効かない領域では，断熱ゆらぎ のとき $\cos[kr_{\rm s}(\tau)]$ ，等曲率ゆらぎのとき $\sin[kr_{\rm s}(\tau)]$ で表される振動をし，ゆらぎのタイプによって振動の位相が異な ることがわかる． 断熱ゆらぎの場合には光子-バリオン系のゆらぎの波数が $k = n\pi/r_{\rm s}$ , ( $n = 1,2,\ldots$ )を満たす位置にピークがあり，等 曲率ゆらぎの場合には $k = (n-1/2)\pi/r_{\rm s}$ の位置にピークがある．こ の違いは宇宙の晴れ上がりのときの背景放射のゆらぎに反映される．宇宙の初 期条件に対して観測可能な予言を与える極めて重要な事実である．また，物質 優勢期を考えると重力ポテンシャルはダークマターによってほとんど決定づけ られるのでポテンシャルは振動しない．すると，式(12.6.303)におい て，双極子成分 $v^{\rm (GI)}_\gamma$ も波数の関数として振動していて，そ の位相は ${\delta^{\rm (GI)}_\gamma}'$ の振動の位相と一致する． これは単 極子成分 $\delta^{\rm (GI)}_\gamma$ と位相が$\pi/2$ ずれている．さらに式 (12.6.305)より双極子成分の位相は4重極子成分と位相が一致し，また， 式(12.6.314)より，偏極成分とも位相が一致する．つまり，光子の密 度ゆらぎと偏極のゆらぎはスケールの関数として位相が$\pi/2$ ずれることが わかる．実際には非斉次解があるのでこれよりも多少複雑であるが，定性的 な結果はあまり変わらない．これはゆらぎの初期条件によらない重要な帰結で あり，やはり宇宙の背景放射の観測によって検証可能なのである．

バリオンのゆらぎも光子のゆらぎと全く同様に振る舞う．実際，ここで考えて いる強結合の1次近似では式(12.6.304), (12.6.312)の右辺に おいて0次近似の解 $v_\gamma^{\rm (GI)}=v_{\rm b}^{\rm (GI)}$ , $\delta_\gamma^{\rm (GI)}=3\delta_{\rm b}^{\rm (GI)}/4$ を代入してもよ い．したがって，光子の場合の式(12.6.315)に対応して，

$${v^{\rm (GI)}_{\rm b}}' + \frac{R'}{1 + R} v^{\rm (GI)}_{\rm b} + \frac13 \frac{1}{1 + R} \delta^{\rm (GI)}_{\rm b} = - {\mit\Phi}$$ (L.6.318)

が得られる．この式と式(12.6.311)により $v_{\rm b}^{\rm (GI)}$ を 消去すれば，バリオンゆらぎの発展方程式

$${\delta^{\rm (GI)}_{\rm b}}'' + \frac{R'}{1+R} {\delta^{\rm (GI)}... ...m (GI)}_{\rm b} = - k^2 {\mit\Phi}- \frac{3R'}{1+R} {\mit\Psi}' - 3{\mit\Psi}''$$ (L.6.319)

が得られる．この式は光子に対する発展方程式(12.6.316)において $\delta_\gamma^{\rm (GI)}=4\delta_{\rm b}^{\rm (GI)}/3$ を代入したもの となっている．したがって，初期条件が同じであれば，ゆらぎの振動の様子は 全く同じである．一般には，その差 $\delta_\gamma^{\rm (GI)}-4\delta_{\rm b}^{\rm (GI)}/3$ は上式の斉次方程式の解である．だが，式 (12.5.273)で見たように初期条件においても近似的に $\delta_\gamma^{\rm (GI)}=4\delta_{\rm b}^{\rm (GI)}/3$ の関係があるか ら、強結合の1次近似でもバリオンと光子はほぼ一体となって振る舞う．　

Footnotes

... よく表すことが知られている^L4

W. Hu and N. Sugiyama, Phys. Rev. D 51, 2599 (1995); Astrophys. J. 444, 489 (1995)