Hyperfine Structure of $B$ and $D$ Mesons in a QCD-Inspired Potential Model

技術要約：QCDに着想を得たポテンシャルモデルにおけるBおよびD中間子の超微細構造

問題提起
量子色力学（QCD）は、非摂動的な低エネルギー領域において直接解くことが依然として困難である。コーネル・ポテンシャルをはじめとするポテンシャルモデルは、クォークと反クォークのダイナミクスに対して直感的な記述を提供するものの、超微細分裂（ベクトル中間子と擬スカラー中間子の質量差）を正確に再現する上で特定の課題に直面している。文献には顕著な不一致が存在する。例えば、チャームオニウムの$2S$分裂は標準的なポテンシャルモデルの予測よりも著しく小さく、一方で軽クォーク系（$\rho-\pi$など）は標準的な期待値を大幅に上回る分裂を示す。さらに、フェルミ・ブライト・ハミルトニアンにおける従来の接触スピン・スピン相互作用項は、原点（$r=0$）において発散するため、正則化が必要となる。したがって、これらの特異性を解決し、相対論的補正を取り入れ、重い・軽い（$D, B$）系および重い・重い（$c\bar{c}, b\bar{b}, B_c$）系の両方において、基底状態（$1S$）と径方向励起状態（$2S$）の両方に対して安定した予測を提供する一貫した枠組みが求められている。

手法
本研究では、以下のコーネル・アンザッツに基づくQCDに着想を得たポテンシャルモデルを採用している：
$$V(r) = -\frac{4\alpha_s}{3r} + br + c$$
ここで、短距離のクーロン項は一 gluon 交換を表し、長距離の線形項は閉じ込めを表す。

主な手法構成要素は以下の通りである：

1. 解析的波動関数： 著者はダルガノ・ルイス（Dalgarno–Lewis）の摂動法を用いており、クーロン部分を親ハミルトニアンとし、線形閉じ込め項を摂動として扱う。これにより、相対論的補正を考慮した因子 $(r/a_0)^{-\epsilon}$ を含む解析的な波動関数 $\psi(r)$ が得られる。
2. 接触項の正則化： フェルミ・ブライトのスピン・スピン相互作用 $V_{SS}(r) \propto \delta^{(3)}(\vec{r})$ の特異性を解決するため、著者はディラックのデルタ関数をガウス型のスメアリング関数に置き換えている：
   $$\delta^{(3)}(\vec{r}) \rightarrow \frac{\sigma^3}{\pi^{3/2}} e^{-\sigma^2 r^2}$$
   ここで、$\sigma$ は構成クォークの実効的なサイズを表す現象論的なスメアリング・パラメータ（1.8 GeVに設定）である。これにより、有限で物理的に意味のある超微細積分が保証される。
3. 質量スペクトルの計算： 物理的質量は、スピンに依存しない中心質量と、スピン依存項（スピン・軌道、テンソル、および超微細項）の総和として計算される。S波の基底状態において、質量差は主に超微細項 $\Delta E_{HF}$ によって支配される。
4. パラメータ空間： 本モデルは、固定された構成クォーク質量（$m_u=m_d=0.33$ GeV, $m_s=0.50$ GeV, $m_c=1.50$ GeV, $m_b=4.80$ GeV）と、実効的な強い結合定数 $\alpha_s = 0.36$ を使用する。結果の $\alpha_s$（0.34–0.38の範囲内）に対する感度は明示的にテストされている。

主な貢献

• 径方向励起への拡張： 基底状態のみに焦点を当てた従来の多くの研究とは異なり、本研究では $D, D_s, B, B_s$ 中間子の第1径方向励起状態（$2S$）に対する超微細分裂を明示的に計算している。
• 不等質量系に対する統一的枠組み： 本モデルは、同一のパラメータセット（実効 $\alpha_s$、弦張力 $b$、およびスメアリング幅 $\sigma$）を、重い・軽い（$D, B$）系と重い・重い（$c\bar{c}, b\bar{b}$）系の両方に適用することに成功しており、質量の非対称性が異なるシステム間でのパラメータ空間の安定性を示している。
• 正則化を伴う解析的手法： ダルガノ・ルイスの解析的波動関数とガウス型スメアリングの組み合わせは、数値的な発散を回避しつつ、短距離における物理的一貫性を維持する、コンパクトで検証可能な枠組みを提供している。

結果

• 基底状態（$1S$）の精度： 本モデルは、重い・軽い中間子の実験的な基底状態質量を、通常3 MeV未満の偏差で再現している。計算された$1S$状態の超微細分裂（例：$D^*-D \approx 140$ MeV, $B^*-B \approx 45$ MeV, $\Upsilon(1S)-\eta_b(1S) \approx 60$ MeV）は、実験値（PDG 2024）および他の理論モデル（Godfrey-Isgurなど）と密接に一致している。
• 径方向励起（$2S$）の予測： 本モデルは、$1S$状態と比較して$2S$状態の超微細分裂が減少することを予測しており、これは原点における波動関数の密度が低いという期待と一致している。予測値には、$D(2S)$ で $\sim 70$ MeV、$B(2S)$ で $\sim 36$ MeV、チャームオニウム $2S$ で $\sim 57$ MeV が含まれる。
• $B_c$ 系： $B_c$ 系に関する現象論的な推定値も提供されており、$1/(m_q m_{\bar{q}})$ のスケーリング則との整合性を示している。具体的には、$B_c$ の分裂はフェルミ・ブライトのスケーリングを通じて得られ、$B_c(2S)$ の分裂は37 MeVと推定されている。
• 感度分析： 実効結合 $\alpha_s$ を0.34–0.38の範囲内で変化させても、超微細分裂の変化は3%未満であり、結果が小さな結合の変動に対して堅牢であることを示している。

意義と主張
本論文は、基底状態の現象論と励起状態の予測の間の溝を埋める、中間子の超微細構造に関する正確かつ統一的な記述を提供すると主張している。ガウス型スメアリングを組み込むことで、接触相互作用に固有の短距離の特異性を、恣意的なパラメータ化なしに解決している。著者は、本フレームワークがLHCbやBelle IIといった施設での今後の分光測定に向けた「コンパクトで検証可能なベンチマーク」を提供すると断言している。$2S$状態および$B_c$系の包含は、不等質量系および径方向励起に対するポテンシャル・パラメータ空間の安定性を検証する役割を果たしている。本研究は、線形項の摂動的取り扱いがより高い励起状態の精度を制限する可能性があるものの、現在のアプローチは$1S$および$2S$セクターにおける中間子の結合の不可欠なダイナミクスを捉えることに成功していると結論付けている。