熱的残存粒子 (ねつてきざんぞんりゅうし, thermal relic) とは、現代宇宙論 において、初期宇宙 においては熱平衡 にあったものの、宇宙の歴史 のある時点において相互作用が実効的に無視できるようになり、それ以降自由粒子として振る舞うような粒子の総称である[ 1] [ 2] 。「残存物」ともいう[ 3] 。相互作用が無視できるようになった時点でその粒子が相対論的であるか非相対論的であるかに応じて「熱い」残存粒子または「冷たい」残存粒子と呼ばれる。例えば宇宙ニュートリノ背景 は熱い残存粒子であり、暗黒物質 の候補であるWIMPは冷たい残存粒子モデルに基づくものである。
概要 ある粒子 χ {\displaystyle \chi } が自身の反粒子 χ ¯ {\displaystyle {\bar {\chi }}} との対消滅 チャンネルを持つ場合を考える。この対消滅反応の反応率は反粒子の数密度 n χ ¯ {\displaystyle n_{\bar {\chi }}} 、対消滅断面積 σ {\displaystyle \sigma } 、粒子速度 v {\displaystyle v} の積
P = n χ ¯ σ v {\displaystyle P=n_{\bar {\chi }}\sigma v} である[ 4] 。高温高密度状態にある初期宇宙 においてはこの反応が頻繁に起き、粒子 χ {\displaystyle \chi } はこの反応を通じて熱平衡にある。やがて宇宙が膨張し温度と数密度が下がると、反応率 P {\displaystyle P} が減少し宇宙の膨張率(ハッブルパラメータ ) H {\displaystyle H} を下回る。これにより対生成・対消滅反応が実効的に停止し(これを「脱結合 (decoupling)」と呼ぶ)、粒子 χ {\displaystyle \chi } の存在量が変化しなくなる[ 2] 。その後この粒子が別の過程により崩壊しない限り、このようにして初期に熱平衡にあった粒子は現在の宇宙にも一定量が存在することになる。これが熱的残存粒子である。
例えばニュートリノ は弱い相互作用 により初期宇宙においてはレプトン(特に電子 ・陽電子 )と相互作用し熱平衡にあるが、およそ温度 1.5 MeV で脱結合し、その後は自由粒子として振る舞う[ 5] 。このように初期宇宙において熱的に生成されたニュートリノは現在の宇宙でも宇宙ニュートリノ背景 として存在し続けていると考えられており、間接的な証拠はあるものの、その数密度は一世代あたり112個/cm3 と小さく[ 5] 、またエネルギーが小さいため、2016年現在それを実験的に検出することはできていない[ 6] 。
あるいは、熱的残存粒子は暗黒物質 を生成するシナリオのひとつとして有力視されている。特に相互作用断面積や質量として電弱相互作用 から予想される値を仮定するとき、冷たい残存粒子シナリオはその粒子の現在の量として暗黒物質の存在量を再現するため、この種の粒子(WIMPとして知られる)は特に有力な暗黒物質の候補となっている[ 7] 。
熱的残存粒子の存在量 本節では脱結合後の熱的残存粒子の存在量がどのように決定されるのかを述べる。なお通常の熱的残存粒子シナリオでは問題の粒子 χ {\displaystyle \chi } とその反粒子 χ ¯ {\displaystyle {\bar {\chi }}} の性質は同じであると仮定する。本記事もこれを仮定するが、バリオン数生成 の理論に動機づけられて粒子・反粒子間に非対称性を導入するモデルもあり、これはAsymmetric dark matterとして知られている[ 8] 。
ボルツマン方程式 問題の粒子 χ {\displaystyle \chi } はその反粒子 χ ¯ {\displaystyle {\bar {\chi }}} との対生成 ・対消滅 チャンネルを持つものと仮定する[ 9] 。
χ + χ ¯ ⇌ X + X ¯ {\displaystyle \chi +{\bar {\chi }}\rightleftharpoons X+{\bar {X}}} ただしここではすべてのチャンネルで反応先の粒子 X {\displaystyle X} , X ¯ {\displaystyle {\bar {X}}} は輻射温度 T {\displaystyle T} に等しい熱平衡にあるものとする[ 1] 。
この系は一般相対論的なボルツマン方程式 によって記述されるが、それを積分することにより粒子数密度 n χ {\displaystyle n_{\chi }} の時間変化を記述する次の方程式が導かれる[ 10] [ 11] [ 4] 。
d n χ d t + 3 H n χ = − ⟨ σ v ⟩ [ n χ 2 − ( n χ E Q ) 2 ] {\displaystyle {\frac {dn_{\chi }}{dt}}+3Hn_{\chi }=-\langle \sigma v\rangle \left[n_{\chi }^{2}-\left(n_{\chi }^{\mathrm {EQ} }\right)^{2}\right]} 左辺第2項は宇宙膨張による希釈を、右辺は対生成・対消滅による粒子数密度の変化を表す。なお右辺の ⟨ σ v ⟩ {\displaystyle \langle \sigma v\rangle } は対消滅断面積と速度の積を熱平均したものであり[ 4] 、 n χ E Q {\displaystyle n_{\chi }^{\mathrm {EQ} }} は熱平衡時の数密度
n χ E Q = g χ ∫ 1 e β ε p ∓ 1 d 3 p ( 2 π ℏ ) 3 = g χ ζ ( 3 ) k B 3 π 2 ℏ 3 c 3 T 3 × 1 2 ζ ( 3 ) ∫ 0 ∞ y 2 exp x 2 + y 2 ∓ 1 d y {\displaystyle n_{\chi }^{\mathrm {EQ} }=g_{\chi }\int {\frac {1}{e^{\beta \varepsilon _{\mathbf {p} }}\mp 1}}{\frac {d^{3}p}{(2\pi \hbar )^{3}}}=g_{\chi }{\frac {\zeta (3)k_{\mathrm {B} }^{3}}{\pi ^{2}\hbar ^{3}c^{3}}}T^{3}\times {\frac {1}{2\zeta (3)}}\int _{0}^{\infty }{\frac {y^{2}}{\exp {\sqrt {x^{2}+y^{2}}}\mp 1}}dy} である。後者は粒子 χ {\displaystyle \chi } が相対論的であるときには n χ = g χ ′ ζ ( 3 ) k B 3 π 2 ℏ 3 c 3 T 3 {\displaystyle n_{\chi }=g'_{\chi }{\frac {\zeta (3)k_{\mathrm {B} }^{3}}{\pi ^{2}\hbar ^{3}c^{3}}}T^{3}} , 非相対論的であるときには n χ = g χ ( m χ k B T 2 π ℏ 2 ) 3 2 e − β m χ c 2 {\displaystyle n_{\chi }=g_{\chi }\left({\frac {m_{\chi }k_{\mathrm {B} }T}{2\pi \hbar ^{2}}}\right)^{\frac {3}{2}}e^{-\beta m_{\chi }c^{2}}} と書ける。ここに ζ ( s ) {\displaystyle \zeta (s)} はリーマンゼータ関数 ( ζ ( 3 ) {\displaystyle \zeta (3)} はアペリーの定数 )、 g χ ′ {\displaystyle g'_{\chi }} は粒子 χ がボース粒子 であるときその内部自由度 g χ {\displaystyle g_{\chi }} 、フェルミ粒子 であるとき 3 4 g χ {\displaystyle {\frac {3}{4}}g_{\chi }} である。
独立変数として時刻 t {\displaystyle t} の代わりに温度 T {\displaystyle T} を用いる。ただし温度は宇宙膨張とともに減少するため、粒子 χ の質量 m χ {\displaystyle m_{\chi }} で無次元化した変数
x := m χ c 2 k B T {\displaystyle x:={\frac {m_{\chi }c^{2}}{k_{\mathrm {B} }T}}} を導入する[ 12] [ 13] 。このとき、 x ≲ 1 {\displaystyle x\lesssim 1} のとき粒子 χ {\displaystyle \chi } は相対論的であり、 x ≫ 1 {\displaystyle x\gg 1} のとき非相対論的である[ 3] [ 13] 。さらに、宇宙膨張による希釈の効果には興味がないため、(あらゆる粒子の寄与を考慮した)エントロピー密度 s {\displaystyle s} で規格化した、単位エントロピーあたりの粒子数
Y := n χ s {\displaystyle Y:={\frac {n_{\chi }}{s}}} を用いる[ 12] [ 14] 。このとき、上の数密度に関する方程式は次のように書き直せる[ 15] [ 16] 。
d Y d x = − π 45 g ∗ S g ∗ ⟨ σ v ⟩ x 2 [ Y 2 − Y E Q 2 ] {\displaystyle {\frac {dY}{dx}}=-{\sqrt {\frac {\pi }{45}}}{\frac {g_{*S}}{\sqrt {g_{*}}}}{\frac {\langle \sigma v\rangle }{x^{2}}}\left[Y^{2}-Y_{\mathrm {EQ} }^{2}\right]} ここに g ∗ {\displaystyle g_{*}} , g ∗ S {\displaystyle g_{*S}} は相対論的な粒子の有効自由度である。
熱い残存粒子 時刻 x f ≲ 1 {\displaystyle x_{f}\lesssim 1} にこの粒子が脱結合したとすると、それ以降 Y ( x ) = n χ ( x ) / s ( x ) {\displaystyle Y(x)=n_{\chi }(x)/s(x)} は近似的に保存する(これを「凍結 (freeze out)」と表現する[ 15] )。そこで現在の宇宙での値と脱結合の瞬間での値を等置する。
n χ ( x f ) s ( x f ) = n χ 0 s 0 {\displaystyle {\frac {n_{\chi }(x_{f})}{s(x_{f})}}={\frac {n_{\chi 0}}{s_{0}}}} 相対論的な粒子の数密度の表式 n χ = g χ ′ ζ ( 3 ) k B 3 π 2 ℏ 3 c 3 T 3 {\displaystyle n_{\chi }=g'_{\chi }{\frac {\zeta (3)k_{\mathrm {B} }^{3}}{\pi ^{2}\hbar ^{3}c^{3}}}T^{3}} およびエントロピー密度の表式 s ( x ) = g ∗ S ( x ) 2 π 2 k B 4 45 ℏ 3 c 3 T 3 {\displaystyle s(x)=g_{*S}(x){\frac {2\pi ^{2}k_{\mathrm {B} }^{4}}{45\hbar ^{3}c^{3}}}T^{3}} を代入すると、左辺の値は Y f := 45 ζ ( 3 ) 2 π 4 k B g χ ′ g ∗ S ( x f ) {\displaystyle Y_{f}:={\frac {45\zeta (3)}{2\pi ^{4}k_{\mathrm {B} }}}{\frac {g'_{\chi }}{g_{*S}(x_{f})}}} と評価される[ 17] 。
現在の宇宙では粒子 χ {\displaystyle \chi } は非相対論的であると仮定すると、その現在のエネルギー密度は
ρ χ 0 c 2 = m χ c 2 n χ 0 = m χ c 2 s 0 Y f {\displaystyle \rho _{\chi 0}c^{2}=m_{\chi }c^{2}n_{\chi 0}=m_{\chi }c^{2}s_{0}Y_{f}} と書ける。従って、粒子 χ {\displaystyle \chi } の密度パラメータ Ω χ 0 {\displaystyle \Omega _{\chi 0}} は
Ω χ 0 = 2 × 8 π G ρ χ 0 3 H 0 2 = 16 ζ ( 3 ) 3 π G m χ ℏ 3 c 3 g ∗ S 0 ( k B T 0 ) 3 H 0 2 g ′ g ∗ S ( x f ) {\displaystyle \Omega _{\chi 0}=2\times {\frac {8\pi G\rho _{\chi 0}}{3H_{0}^{2}}}={\frac {16\zeta (3)}{3\pi }}{\frac {Gm_{\chi }}{\hbar ^{3}c^{3}}}g_{*S0}{\frac {(k_{\mathrm {B} }T_{0})^{3}}{H_{0}^{2}}}{\frac {g'}{g_{*S}(x_{f})}}} により与えられる[ 17] 。なお係数 2 は反粒子 χ ¯ {\displaystyle {\bar {\chi }}} の寄与を考慮することを表す。CMB 温度 T 0 {\displaystyle T_{0}} およびハッブル定数 H 0 {\displaystyle H_{0}} の観測値を代入すると
Ω χ 0 = 0.15 h − 2 g ′ g ∗ S ( x f ) ( m χ c 2 1 e V ) {\displaystyle \Omega _{\chi 0}=0.15h^{-2}{\frac {g'}{g_{*S}(x_{f})}}\left({\frac {m_{\chi }c^{2}}{1\,\mathrm {eV} }}\right)} という結論を得る[ 17] 。
冷たい残存粒子 冷たい残存粒子の存在量 Y {\displaystyle Y} の温度変化の数値解[ 18] 。断面積としては温度に依存しないs波対消滅を仮定した。実線、破線、一点短鎖線が異なる λ {\displaystyle \lambda } に対する数値解であり、点線は熱平衡が保たれると仮定した場合の値。熱平衡の値 Y E Q {\displaystyle Y_{\mathrm {EQ} }} は指数関数的に減少する。しかし数値解は脱結合により断面積に依存したある時点で熱平衡からの値から逸脱し、ある値 Y ( ∞ ) {\displaystyle Y(\infty )} で「凍結」する。 粒子 χ {\displaystyle \chi } が温度 T {\displaystyle T} の熱平衡にあるとき、それが非相対論的である ( k B T ≪ m χ c 2 {\displaystyle k_{\mathrm {B} }T\ll m_{\chi }c^{2}} ) ならば、その数密度 n χ {\displaystyle n_{\chi }} は
n χ = g χ ( m χ k B T 2 π ℏ 2 ) 3 2 e − β m χ c 2 {\displaystyle n_{\chi }=g_{\chi }\left({\frac {m_{\chi }k_{\mathrm {B} }T}{2\pi \hbar ^{2}}}\right)^{\frac {3}{2}}e^{-\beta m_{\chi }c^{2}}} と書ける。従って熱平衡にあるときの Y {\displaystyle Y} の値は
Y E Q ( x ) = 45 4 2 π 7 / 2 g χ g ∗ S ( x ) x 3 2 e − x {\displaystyle Y_{\mathrm {EQ} }(x)={\frac {45}{4{\sqrt {2}}\pi ^{7/2}}}{\frac {g_{\chi }}{g_{*S}(x)}}x^{\frac {3}{2}}e^{-x}} と表示できる[ 19] 。この場合の最終的な粒子 χ {\displaystyle \chi } の存在量は Y {\displaystyle Y} に関する方程式を数値的に解くことによって求められる。
しばしば ⟨ σ v ⟩ {\displaystyle \langle \sigma v\rangle } が温度のべき乗の依存性を持つと仮定される[ 19] [ 20] : ⟨ σ v ⟩ ∝ T n {\displaystyle \langle \sigma v\rangle \propto T^{n}} 。この場合、定数 λ {\displaystyle \lambda } を
λ := s H ⟨ σ v ⟩ | x = 1 {\displaystyle \lambda :=\left.{\frac {s}{H}}\langle \sigma v\rangle \right|_{x=1}} により定義すると、 Y {\displaystyle Y} に関する方程式は
d Y d x = − λ x n + 2 [ Y 2 − Y E Q 2 ] {\displaystyle {\frac {dY}{dx}}=-{\frac {\lambda }{x^{n+2}}}\left[Y^{2}-Y_{\mathrm {EQ} }^{2}\right]} と書き直せる[ 19] 。いくつかの λ {\displaystyle \lambda } に対する数値解を図に示す。ここからわかるように、初期には解 Y ( x ) {\displaystyle Y(x)} は熱平衡の場合の値 Y E Q {\displaystyle Y_{\mathrm {EQ} }} に一致するが、ある時刻 x f {\displaystyle x_{f}} でそこから逸脱する。この時刻はガモフの基準により方程式
x f n − 1 2 e x f = 45 4 2 π 7 / 2 g χ g ∗ S λ {\displaystyle x_{f}^{n-{\frac {1}{2}}}e^{x_{f}}={\frac {45}{4{\sqrt {2}}\pi ^{7/2}}}{\frac {g_{\chi }}{g_{*S}}}\lambda } を満足する値 x f {\displaystyle x_{f}} として概算できる[ 21] 。その後、対生成・対消滅反応が停止し、 Y {\displaystyle Y} は最終的な値 Y ( ∞ ) {\displaystyle Y(\infty )} に「凍結」する。その値は
Y ( ∞ ) ∼ Y E Q ( x f ) = x f n + 1 λ {\displaystyle Y(\infty )\sim Y_{\mathrm {EQ} }(x_{f})={\frac {x_{f}^{n+1}}{\lambda }}} である[ 22] 。
現在の宇宙における冷たい残存粒子の密度パラメータ Ω χ 0 {\displaystyle \Omega _{\chi 0}} は次のように求まる[ 22] 。
Ω χ 0 = 2 × 8 π G 3 H 0 2 m χ s 0 Y ( ∞ ) = 32 π 5 / 2 9 5 m P l 3 g ∗ S 0 T 0 3 H 0 2 g ∗ g ∗ S ⟨ σ v ⟩ 1 x f n + 1 {\displaystyle \Omega _{\chi 0}=2\times {\frac {8\pi G}{3H_{0}^{2}}}m_{\chi }s_{0}Y(\infty )={\frac {32\pi ^{5/2}}{9{\sqrt {5}}m_{\mathrm {Pl} }^{3}}}{\frac {g_{*S0}T_{0}^{3}}{H_{0}^{2}}}{\frac {\sqrt {g_{*}}}{g_{*S}\langle \sigma v\rangle _{1}}}x_{f}^{n+1}}
観測的制限
熱い残存粒子 熱い残存粒子が現在の宇宙に非相対論的な粒子として存在するならば、それはホットダークマター として知られるタイプの暗黒物質 として振る舞うが、その存在は宇宙論的観測によって棄却されている[ 23] 。そのため、熱い残存粒子の密度パラメータは現在の暗黒物質の量 Ω c 0 = 0.2582 {\displaystyle \Omega _{c0}=0.2582} に比べて十分小さい必要があり、このことから熱い残存粒子の質量 m χ {\displaystyle m_{\chi }} に上限が与えられる。例えばニュートリノ の場合、この考察から m ν ≪ 10 e V {\displaystyle m_{\nu }\ll 10\,\mathrm {eV} } という制限が得られる[ 17] 。ただしこの種の制限は g ∗ ( x f ) {\displaystyle g_{*}(x_{f})} の値を通じて断面積の大きさに依存し、断面積が小さく g ∗ ( x f ) ≥ 106.75 {\displaystyle g_{*}(x_{f})\geq 106.75} で凍結する場合には熱い残存粒子の質量の上限は 10 e V {\displaystyle 10\,\mathrm {eV} } 程度まで緩和される[ 19] 。
冷たい残存粒子 冷たい残存粒子はコールドダークマター と呼ばれるタイプの暗黒物質となる[ 19] 。冷たい残存粒子の存在量は強く対消滅断面積に依存するため、質量に関する制限もまたその相互作用モデルに強く依存する[ 22] が、宇宙論的に意味のある量の残存粒子が生成されるためには断面積は極めて小さい必要がある[ 24] 。一般論としては、理論のユニタリー性により相互作用断面積は粒子質量と
σ ≲ 4 π m χ 2 {\displaystyle \sigma \lesssim {\frac {4\pi }{m_{\chi }^{2}}}} という関係にあるため、 Ω χ ≲ 0.2 {\displaystyle \Omega _{\chi }\lesssim 0.2} は m χ ≲ 120 T e V {\displaystyle m_{\chi }\lesssim 120\,\mathrm {TeV} } を要求する[ 25] 。
具体的に相互作用断面積およびその質量として電弱相互作用 から示唆される値を用いるとき、冷たい残存粒子シナリオが予測する暗黒物質量は現在の観測値を「奇跡的」に再現する[ 7] 。この事実は「WIMPの奇跡[ 26] (WIMP miracle)」として知られており[ 7] 、暗黒物質が weakly interacting massive particle (WIMP) と呼ばれる種類の素粒子であると考える根拠のひとつとなっている。
脚注 ^ a b 松原, p. 39. ^ a b Profumo, p. 34-35. ^ a b 現代の天文学2, pp. 91-92. ^ a b c Profumo, pp. 47-54. ^ a b 松原隆彦『宇宙論の物理 上』東京大学出版会、2014年、92-93頁。ISBN 978-4130626156。 ^ 中家剛『ニュートリノ物理 ―ニュートリノで探る素粒子と宇宙―(基本法則から読み解く物理学最前線 9)』共立出版、2016年、41頁。ISBN 9784320035294。 ^ a b c Profumo, pp. 39-41. ^ Profumo, pp. 206-210. ^ Profumo, p. 36. ^ 松原, pp. 27-29. ^ 現代の天文学2, pp. 86-90 ^ a b 松原, p. 40. ^ a b Profumo, p. 39. ^ 現代の天文学2, p. 90. ^ a b 松原, p. 41. ^ Profumo, p. 54. ^ a b c d 松原, p. 42. ^ 松原, p. 44, 図6.5. ^ a b c d e 松原, p. 43. ^ Profumo, p. 55. ^ 松原, p. 44. ^ a b c 松原, p. 45. ^ 松原, p. 41-42. ^ 現代の天文学2, p. 92. ^ Profumo, p. 42. ^ カイ・マルテンス. “なぜWIMPを探すのか、どうやって捕らえるのか”. 2020年8月1日 閲覧。
参考文献 松原隆彦『宇宙論の物理 下』東京大学出版会、2014年。ISBN 978-4130626163。 佐藤勝彦, 二間瀬敏史 編『シリーズ現代の天文学2 宇宙論1』(第2版)日本評論社、2012年。ISBN 978-4535607392。 Profumo, Stefano (2017). An Introduction To Particle Dark Matter . World Scientific Publishing Europe. ISBN 978-1786340016
関連項目