論文「Light-cone sum rules with B-meson distribution amplitudes for the B → p form factors in B-mesogenesis models」の技術的サマリー

1. 研究の背景と問題設定

**B メソゲネシス（B-mesogenesis）**シナリオは、宇宙のバリオン非対称性とダークマターの存在量を同時に説明しようとする理論モデルです。このシナリオでは、B メソンが軽いバリオン（陽子など）と「見えないダーク反バリオン（$\Psi$）」に崩壊する過程（$B \to p\Psi$）が特徴的なシグネチャとして予言されています。

本研究の主要な課題は、この崩壊の分岐比を理論的に見積もるために必要な**$B \to p$ 遷移形状因子（form factors）**を、より高精度かつ独立した手法で計算することです。
既存の研究（例：[10, 12]）では、**核子（陽子）の分布振幅（DAs）**を用いた QCD ライトコーン和則（LCSRs）が適用されていましたが、以下の問題点が指摘されていました：

1. 高次ツイスト（twist-6 まで）の寄与が非常に大きく、展開の収束性が遅い。
2. 核子 DAs のパラメータに大きな不確実性が残っている。
3. 結果として、形状因子の予測値に大きなばらつきが生じている。

これらの課題を解決し、実験的な検証可能性を高めるため、本研究では**B メソンの分布振幅（B-meson DAs）**を用いた新しい LCSRs 手法を適用し、$B \to p$ 形状因子を再計算しました。

2. 手法と理論的枠組み

2.1. 有効相互作用とモデル

B メソゲネシスモデルの 2 つのバージョン、モデル (d) とモデル (b) を対象とします。

• モデル (d): $d$ 夸子がダーク反バリオン $\Psi$ と重媒介粒子に直接結合。
• モデル (b): $b$ 夸子が $\Psi$ と重媒介粒子に直接結合。
  これらは有効 4 フェルミ相互作用として記述され、$B^+ \to p\Psi$ 崩壊振幅は $B \to p$ 遷移行列要素に比例します。

2.2. 相関関数の構成と OPE

本研究の核心は、相関関数の構成を逆転させる点にあります。

• 従来の手法 (LCSR-N): 真空から陽子への相関関数（陽子 DAs を使用）。
• 本研究の手法 (LCSR-B): 真空から B メソンへの相関関数を構成します。
  • 相関関数: $F(P, q) = i \int d^4x e^{iP\cdot x} \langle 0 | T \{ \eta_p(x), O(0) \} | B^+(P+q) \rangle$
  • ここで $\eta_p$ は陽子の補間電流（Ioffe 電流）、$O$ は有効相互作用演算子です。
  • この相関関数の Operator Product Expansion (OPE) を、**B メソンの分布振幅（HQET 枠組みで定義）**を用いて展開します。

2.3. 計算の詳細

• ツイスト展開: 精度をツイスト 5まで高めて計算を行いました。
• 主要な結果:
  • 2 粒子 B メソン DAs（twist-2, 3, 4, 5）の寄与は計算されました。
  • 驚くべき結果: 3 粒子（クォーク・反クォーク・グルーオン）DAs の寄与は、選択した陽子補間電流と有効相互作用のスカラー構造、および質量less な $u, d$ 夸子の近似により、完全にゼロになることが示されました。
  • 非摂動入力として、B メソン崩壊定数 $f_B$、逆モーメント $\lambda_B$、およびパラメータ比 $R$ を使用しました。

2.4. 形状因子の抽出と外挿

• 計算された OPE 側とハドロン側（陽子極と連続状態）を比較し、8 つの独立した LCSRs を導出しました。
• 得られた形状因子は、空間的領域（$q^2 < 0$）で計算され、その後 **z 展開（z-expansion）**を用いて物理的領域（$q^2 = m_\Psi^2 > 0$）へ外挿されました。
• 外挿には $\Lambda_b$ ボーの極を除去する変換を施し、解析性を保ちながら物理的な質量領域へ接続しました。

3. 主要な結果

3.1. 形状因子の計算値

モデル (d) と (b) に対する $B \to p$ 形状因子（$F_{B\to p}^R$ と $\tilde{F}_{B\to p}^L$）を $q^2=0$ において評価しました。

• モデル (d) と (b) の関係: 本研究の手法では、両モデルの形状因子が単純な係数 2 の関係（$F^{(b)} = 2 F^{(d)}$）で結ばれることが示されました。これは従来の核子 DAs を用いた手法では見られなかった明確な関係です。
• 高次ツイストの挙動: 従来の手法（LCSR-N）では高次ツイストの寄与が支配的でしたが、本研究（LCSR-B）では低次ツイスト（twist-2, 3）が支配的であり、高次項（twist-4, 5）は抑制されています。これは OPE の収束性が良好であることを示唆しています。
• 不確実性: 主要な不確実性は、B メソン DAs の逆モーメント $\lambda_B$ と核子崩壊定数 $\lambda_p$ に由来しますが、従来の手法と比較して予測範囲が狭まり、系統的に小さな値を示す傾向があります。

3.2. 分岐比の予測

計算された形状因子を用いて、$B^+ \to p\Psi$ の分岐比をダーク反バリオン質量 $m_\Psi$ の関数として予測しました。

• 有効結合定数の上限（LHC 検索結果に基づく）を仮定した場合、分岐比は $m_\Psi$ に依存して変化します。
• 従来の手法と比較して、本研究の予測値はより小さく、かつ不確実性が小さい範囲に収束しています。

3.3. 排他的崩壊と包括的崩壊の比率

重要なステップとして、排他的崩壊幅（$B \to p\Psi$）と包括的崩壊幅（$B \to X_N \Psi$）の比率を計算しました。

• 包括的崩壊幅は、重クォーク展開（HQE）を用いて [20, 21] で計算された結果を使用し、$m_\Psi \approx 2.0$ GeV (モデル d) および $3.0$ GeV (モデル b) 以下の領域で信頼性が高いと判断しました。
• この比率を用いることで、未知の有効結合定数を相殺し、結合定数に依存しない下限を導出できます。
• 得られた下限は、現在の BaBar や Belle/Belle II による実験的上限と比較されました。

4. 結論と意義

4.1. 科学的意義

1. 手法の独立性と検証: B メソン DAs を用いた LCSRs は、核子 DAs を用いた従来の手法とは独立したアプローチであり、両者の結果が大きな不確実性の中で一致することは、B メソゲネシスモデルの理論的基盤の堅牢性を裏付けています。
2. OPE の収束性: 本研究の手法では高次ツイスト項が抑制されており、OPE 展開の信頼性が従来の手法よりも高いことを示しました。
3. 実験への指針: 現在の実験的上限（$O(10^{-6})$）と、本研究で導出された理論的下限（$O(10^{-8} \sim 10^{-7})$）の間にまだギャップがあります。B メソゲネシスの決定的な検証のためには、実験的な感度を10 桁以上（$10^{-8} \sim 10^{-7}$ レベル）向上させる必要があることが示されました。

4.2. 今後の展望

• 本研究で確立された「B メソン DAs を用いた LCSRs」は、陽子以外の重バリオン（$\Lambda, \Lambda_c$ など）への遷移形状因子の計算にも容易に拡張可能です。これにより、B メソゲネシスモデルの他の崩壊モード（$B \to \Lambda \Psi$ など）の予測も可能になります。
• 格子 QCD による $B \to p$ 形状因子の計算がまだ存在しない現状において、本研究の手法はハドロン入力を提供する唯一の QCD 基礎的手法として極めて重要です。

総じて、本論文は B メソゲネシスという新物理シナリオの検証可能性を、より厳密な QCD 計算によって評価し、将来の実験計画に重要な指針を提供するものです。