Raubold-Lynch construction of the phase space in hypertriton three-body mesonic decay

技術要約：ハイパートリトンにおける三体メソン崩壊の位相空間におけるRaubold–Lynch構成

問題提起
ハイパートリトン（$^3_\Lambda$H）は、既知の中で最も軽いハイパー核であり、$\Lambda$－核子相互作用を調査するための重要な実験場である。DalitzやDownsらによる先駆的な研究は、二体メソン崩壊（$^3_\Lambda$H $\to$ $^3$He + $\pi^-$）の理論的枠組みを確立したが、三体メソン崩壊チャネル（$^3_\Lambda$H $\to$ $d + p + \pi^-$ および $^3_\Lambda$H $\to$ $d + n + \pi^0$）はより複雑な運動学的景観を呈している。負のパイオンチャネルは実験的に測定されているが（例：STARコラボレーション）、中性パイオンチャネルについては実験データが存在しない。さらに、連続スペクトラムを計算するための理論的枠組み（Motobaら、Kamadaら）は存在するものの、厳密でモデルに依存しない運動学的ベースラインが必要とされている。このベースラインは、精密な崩壊率や観測量を計算するために必要な、純粋な位相空間効果と、最終状態相互作用や弱崩壊行列要素などの動力学的寄与を分離するために必要である。

手法
本論文では、Raubold–Lynch逐次二体崩壊アルゴリズムを用いた、ハイパートリトンの三体メソン崩壊位相空間のモンテカルロ構成を提示する。核心となる手法は、以下のステップを含む：

1. 再帰的因子分解: 三体崩壊を、仮想的な中間系（$X$）を介した二つの逐次的な二体崩壊へと因子分解する。

   • 荷電チャネルの場合：$^3_\Lambda$H $\to$ $\pi^- + X$、続いて $X \to p + d$。
   • 中性チャネルの場合：$^3_\Lambda$H $\to$ $\pi^0 + X$、続いて $X \to n + d$。
     ここで、$X$はサンプリングされた不変質量（$m_{pd}$ または $m_{nd}$）を持つ、オフシェル粒子としての核子－重陽子系を表す。

2. 運動学的サンプリング:

   • 中間系の不変質量は、親粒子の質量と娘粒子の質量によって決定される、運動学的に許容される全範囲にわたってサンプリングされる。
   • 各サンプリングされた不変質量に対して、アルゴリズムは厳密な相対論的二体運動学を実行し、それぞれの静止系における崩壊運動量を決定する。
     ・等方的な角度分布は、単位球面上での一様分布を保証するため、二つのランダムな回転（z-then-y）の合成を用いて生成される。

3. ローレンツ変換: 中間系 $X$ の静止系において最初に計算された核子と重陽子の四元運動量は、第一段階の崩壊から導かれる $X$ の速度を用いて、ハイパートリトンの静止系へとローレンツブーストされる。これにより、最終的な三体状態に対する厳密な四元運動量保存が保証される。

4. 重み付けと受理: アルゴリズムは、二つの段階の崩壊運動量の積に比例する位相空間の重み（$w \propto p^*_1 p^*_2$）を持つイベントを生成する。受容／拒否の手続きを適用することで、この重み付きサンプルを、ローレンツ不変な三体位相空間測度（$d\Phi_3$）に従う重みなしのアンサンブルへと変換する。

5. 観測量の抽出: 得られた重みなしイベントから、両方のパイオンチャネルについて、運動量スペクトル、対となる運動量相関、開き角分布、およびダリッツ図を構築する。

主な貢献

• ハイパートリトンのためのRaubold–Lynchの実装: 本論文は、荷電（$\pi^-$）および中性（$\pi^0$）の両チャネルをカバーする、ハイパートリトンの三体メソン崩壊に特化したRaubold–Lynchアルゴリズムの完全な実装を提供する。
• 運動学的ベースラインの生成: ローレンツ不変な位相空間測度を厳密に再現する、大規模な重みなしモンテカルロイベントのアンサンブルを生成する。これは「純粋な位相空間」のベースラインとして機能し、動力学的仮定から自由である。
• 中性チャネルの特性評価: 中性パイオンの崩壊に関する実験データが存在しないため、本論文は相対論的運動学に基づいた、$^3_\Lambda$H $\to$ $d + n + \pi^0$ 崩壊の詳細な運動学的特性（運動量限界、相関、ダリッツ図）を初めて提示する。
• モジュール型フレームワーク: 本研究は、運動学と動力学の分離を確立している。生成された位相空間サンプルには、将来的に任意の弱崩壊行列要素を用いてイベントごとに重み付けを行うことができ、運動学を再生成することなく崩壊率を計算することが可能である。

結果
本研究は、両方の崩壊モードについて詳細な運動学的分布を示す：

• 不変質量分布: 中間系不変質量（$m_{pd}$ および $m_{nd}$）の分布は、位相空間密度が消失する運動学的端点において特徴的な抑制を示す。中性チャネルは、荷電チャネル（$Q \approx 37.3$ MeV）と比較して広い利用可能エネルギー範囲（$Q \approx 40.6$ MeV）を有しており、その結果、より広い運動学的限界を示す。
• 運動量スペクトル: 両チャネルのすべての三体崩壊生成物は、長い低運動量テールと鋭い高運動量カットオフを持つ非対称な運動量スペクトルを示す。
  • 荷電チャネルにおいて、パイオンの平均運動量は $\approx 70$ MeV/cであり、陽子と重陽子はそれぞれ $\approx 148$ および $\approx 154$ MeV/cである。
  • 三体連続体の端点は、対応する二体崩壊の単色運動量から明確に分離されており（例：$\pi^-$ の端点は二体ピークの $\approx 8.4$ MeV/c 下方にある）、運動学的な識別可能性を確保している。
• 相関と角度:
  • 運動量相関: パイオンと核子の間には強い反相関（ピアソンの $r \approx -0.9$）が観察される一方、核子－核子ペアの間には中程度の正の相関（$r \approx +0.6$）が見られる。
  • 開き角: 核子－核子の開き角分布は $180^\circ$（背中合わせ）付近に鋭くピークを持ち、これは軽いパイオンに対する重い核子ペアの反跳を反映している。パイオン－核子の角度は広く分布している。
• ダリッツ図: 運動エネルギーおよび相対論的不変質量二乗の慣習の両方において、運動学的に許容された領域全体にわたって一様な分布を示し、重みなしの位相空間を正しく再現していることを確認している。

意義と主張
本論文の主要な意義は、そのモジュール性と計算効率にあると主張している。ローレンツ不変な位相空間の生成を動力学的行列要素からデカップリング（分離）することにより、本手法は将来の理論計算のための再利用可能な、モデルに依存しない基礎を提供する。

著者は、生成された重みなしイベントが以下のための必須のリファレンスとなることを強調している：

1. 動力学的効果の定量化: 最終状態相互作用や特定の波動関数を含む将来の計算を、この純粋な位相空間ベースラインと比較することで、動力学的な修正を孤立させることが可能となる。
2. 中性チャネルへの予測: 中性パイオン崩壊に関する実験データが存在しないため、ローレンツ不変な位相空間は理論的予測のための有用なベースラインとなる。
3. 崩壊率の計算: 重みなしサンプルを用いることで、スピン平均された二乗行列要素を適用するだけで、モンテカルロ積分を介して部分崩壊幅を直接計算できる。これにより、異なる動力学モデルに対して運動学を解き直す必要がなくなる。

本研究は、ハイパートリトンの崩壊自体の動力学的問題を解決することを目的としているのではなく、むしろそれを解くために必要な、厳密な運動学的インフラストラクチャを提供することを目的としている。