Primordial deuterium abundance from calculations of $p(n,γ)$ and $d(p,γ)$ reactions within potential-model approach

本論文は、マルフリート・ジョン相互作用に基づくポテンシャルモデルを用いて一貫した枠組みで$p(n,\gamma)$および$d(p,\gamma)$反応を解析し、得られた核反応断面積からビッグバン元素合成における原始デューテリウム存在量$\mathrm{D/H} = 2.479^{+0.350}_{-0.177}\times 10^{-5}$を導出したところ、金属量の少ないダプド・ライマン-$\alpha$系から推定される値と良好に一致したことを報告しています。

要約

1. 宇宙の料理：重水素という「焦げ」

想像してください。ビッグバン直後の宇宙は、超高温の巨大なオーブン（鍋）のようなものでした。そこには、水素（最も単純なブロック）が大量に溢れていました。

• 水素（プロトン） ＝ 単なるレゴブロック。
• 中性子 ＝ もう一つのレゴブロック。
• 重水素 ＝ 水素と中性子がくっついた「2 個セットのブロック」。

このオーブンの中で、ブロック同士が激しくぶつかり合い、くっついて大きなもの（ヘリウムなど）を作ろうとしていました。しかし、**「重水素」**は非常に壊れやすく、すぐにバラバラになってしまいます。

ここで重要なのが、「重水素の瓶詰め」（核反応）です。

1. p(n, γ) 反応：水素と中性子がくっついて重水素を作る「瓶詰め」作業。これが始まると、宇宙の料理は本格的に始まります（これを「デューテリウムのボトルネックの解消」と呼びます）。
2. d(p, γ) 反応：できた重水素が、さらに水素とくっついてヘリウムに変わってしまう「消費」作業。

この研究は、**「この 2 つの作業が、宇宙の温度や圧力の中で、どれくらい速く、正確に行われるか」**を計算し直しました。

2. 研究の手法：魔法の「拡大縮小レンズ」

これまでの研究では、この反応の速さを正確に測るのが難しかったです。実験室で低温の宇宙を再現するのは、氷を溶かさずに測るようなものだからです。

そこで、この論文の著者たちは**「ポテンシャルモデル（可能性のモデル）」**という、少し違うアプローチを取りました。

• 従来の方法：一つ一つの粒子の動きを、複雑なコンピュータで細かくシミュレーションする（高価で時間がかかる）。
• この論文の方法：**「魔法の拡大縮小レンズ（スケーリングファクター λ）」**を使う。

彼らは、**「p(n, γ) 反応（瓶詰め）」**の実験データに合わせて、この「レンズ」の倍率を微調整しました。

• 「あ、この倍率なら実験結果とぴったり合うな」と倍率を決めます。
• その同じ倍率を、**「d(p, γ) 反応（消費）」**にもそのまま適用します。

まるで、**「料理の味付け（塩分）を、スープ（p(n, γ)）で調整したら、同じレシピでメインディッシュ（d(p, γ)）も作れる」**という考え方です。これにより、実験データが少ない部分でも、信頼できる予測が可能になりました。

3. 発見：宇宙のレシピは「繊細」だった

この研究でわかった最も重要なことは、**「この魔法のレンズの倍率（λ）を少し変えるだけで、宇宙の出来上がりが大きく変わる」**ということです。

• 倍率を少し変える ＝ 重水素が作られる速さや消える速さが変わる。
• 結果：宇宙に残る重水素の量（D/H 比）が、劇的に変わります。

彼らが計算した結果、「重水素と水素の比率は、約 10 万分の 2.5」という値になりました。
これは、天文学者が遠くの古いガス雲（金属が少ないダップド・ライマン・アルファ・システム）を望遠鏡で観測して推測した値と驚くほど一致しました。

4. なぜこれが重要なのか？

この研究は、単に数字を合わせただけではありません。

• 宇宙の「バネ」の強さを測る：重水素の量は、宇宙にどれだけの物質（バリオン）があったかを示す「物差し」です。この研究は、その物差しが正しいことを裏付けました。
• 理論のシンプルさ：複雑な計算をしなくても、シンプルなモデル（拡大縮小レンズ）を使えば、宇宙の初期の化学反応を正しく説明できることを示しました。

まとめ

この論文は、**「宇宙という巨大な料理鍋の中で、重水素という『焦げ』がどれくらい残ったか」を、「実験結果に合わせて調整した魔法のレンズ」**を使って計算し直した物語です。

その結果、計算された重水素の量は、実際の宇宙観測と完璧に一致しました。これは、私たちが宇宙の誕生について理解している「レシピ」が、間違いなく正しい方向にあることを強く示唆しています。

一言で言えば：
「宇宙の赤ちゃんの頃、重水素がどれだけ生き残ったかを、シンプルな『拡大鏡』を使って計算し直したら、実際の宇宙の姿とぴったり合っていたよ！」という発見です。

技術概要

この論文「Potential-model approach における p(n, γ) および d(p, γ) 反応の計算から導かれる原始重水素存在量」の技術的な要約を以下に示します。

1. 研究の背景と課題

ビッグバン核合成（BBN）は、宇宙初期の物理と観測される軽元素の存在量を結びつける重要なプロセスであり、その中で重水素（D）の存在量は宇宙論的プローブとして極めて高精度です。しかし、BBN における重水素の生成と破壊を支配する核反応率、特に以下の 2 つの放射性捕獲反応の低エネルギー領域における断面積の精度が、宇宙バリオン密度の制約を改善する上でボトルネックとなっています。

• p(n, γ) 反応: 中性子と陽子の放射性捕獲。重水素の「ボトルネック」を解消し、核合成を開始させる。
• d(p, γ) 反応: 重水素と陽子の放射性捕獲。重水素の燃焼（破壊）を制御する。

実験的には、特に d(p, γ) 反応はクーロン障壁により低エネルギーでの断面積が極めて小さく、BBN 領域（約 100 keV 付近）でのデータは依然として限られています。また、p(n, γ) 反応についても低エネルギーでの測定データは不足しています。したがって、理論的な計算による高精度な反応率の提供が不可欠です。

2. 手法と理論的枠組み

本研究は、一貫した**2 体ポテンシャルモデル（Potential Model, PM）**の枠組みを用いて、p(n, γ) および d(p, γ) 反応を解析しました。

• 相互作用ポテンシャル:
  • Malfliet-Tjon (MT) 相互作用に基づき、スピン軌道相互作用を無視した中心力ポテンシャルを使用しました。
  • ポテンシャルは、束縛状態（bound state）と散乱状態（scattering state）で定義されます。
  • スケーリング因子 $\lambda$: 散乱ポテンシャルは、束縛ポテンシャルをスケーリング因子 $\lambda$ 倍することで得られます（$V^s_{MT} = \lambda V^b_{MT}$）。
• 遷移の計算:
  • 電磁双極子遷移（E1）と磁気双極子遷移（M1）の両方を考慮しました。
  • 波動関数は、中心-of-mass 系における 2 体シュレーディンガー方程式を解いて得られました。
  • 行列要素は、束縛状態と散乱状態の波動関数を用いて計算され、E1 および M1 遷移の寄与を統合しました。
• パラメータの制約:
  • p(n, γ) 反応: 熱中性子捕獲断面積（$E_n = 25.3$ meV）の実験値（$332.6 \pm 0.7$mb）を用いて、スケーリング因子$\lambda_{p(n,\gamma)}$を制約しました。これにより$\lambda \approx 0.734$ が得られました。
  • d(p, γ) 反応: クラスターモデルの考え方に基づき、pd 相互作用は np 相互作用と pp 相互作用の和（$V_{pd} \approx 2V_{np}$）と近似しました。これにより、d(p, γ) のスケーリング因子は $\lambda_{d(p,\gamma)} \approx 2\lambda_{p(n,\gamma)} \approx 1.470$ と設定されました。
• 反応率の導出:
  • 計算された断面積から天体物理学的 S 因子を導出し、Maxwell 平均反応率を計算しました。
  • 得られた反応率は解析的なパラメータ化を行い、BBN コード「PArthENoPE」に実装しました。

3. 主要な結果

• p(n, γ) 反応の断面積:
  • 低エネルギー領域（100 keV 以下）では M1 遷移が支配的であり、計算値は実験データ（Suzuki et al., Nagai et al.）とよく一致しました。
  • 高エネルギーでは E1 遷移の寄与が増加しますが、モデル全体として実験データと整合性がありました。
  • スケーリング因子 $\lambda$ の不確かさ（$\delta\lambda/\lambda \approx 3.1\%$）が断面積の不確かさに直接反映され、特に低エネルギー領域で感度が高いことが示されました。
• d(p, γ) 反応の天体物理学的 S 因子:
  • 低エネルギー領域で S 因子は $S(0) = 0.228 \pm 0.015$ eV b となりました。これは最近の決定値より約 4% 高い値ですが、誤差範囲内で一致しています。
  • 計算された S 因子は、LUNA 実験などの高精度データおよび ab initio 計算（Ref. [31]）と概ね一致しました（ab initio 計算とは約 5% の差があり、これは多体電流や核相関の扱いの違いに起因すると考えられます）。
• 反応率とパラメータ化:
  • 0.001 GK から 10 GK までの温度範囲で反応率を計算し、解析的なパラメータ化（多項式展開）を行いました。
  • 高温（8 GK 以上）では反応率の温度依存性が弱くなる傾向が確認されました。

4. 原始重水素存在量（D/H）への影響

• PArthENoPE を用いた BBN 計算により、制約された $\lambda$ の変動を伝播させて原始重水素存在量を評価しました。
• 得られた結果は D/H = $2.479^{+0.350}_{-0.177} \times 10^{-5}$ となりました。
• この値は、金属の少ないダンプド・ライマン$\alpha$系（Damped Lyman-$\alpha$ systems）から推定される観測値と良好に一致しています。
• 感度分析: $\lambda$ のわずかな変動（実験誤差範囲内）が、予測される D/H 比および軽元素存在量に有意な変化をもたらすことが示されました。これは、低エネルギー散乱ダイナミクスが BBN 予測において極めて重要であることを裏付けています。

5. 結論と意義

• 一貫性の確立: 本研究は、p(n, γ) と d(p, γ) という 2 つの重要な反応を、単一のスケーリング因子 $\lambda$ で制御された一貫したポテンシャルモデル枠組みで記述することに成功しました。
• 実用性: 複雑な第一原理計算（ab initio）と比較して計算コストは低いものの、低エネルギーデータに制約された簡素な枠組みでありながら、BBN に必要な高精度な核入力値を提供できることを示しました。
• 宇宙論への貢献: 得られた反応率と D/H 存在量は、宇宙バリオン密度の制約を強化し、BBN 理論の信頼性を高める上で重要な役割を果たします。特に、核モデルの不確かさが宇宙論的観測量にどのように伝播するかを定量的に評価した点が重要です。

この論文は、理論的な核反応率の計算が、観測的な宇宙論パラメータの精度向上にどのように寄与できるかを明確に示す重要な研究です。