Inclusive $\bar B_s\mapsto X_{\bar sc} \ell \bar \nu$ decays from lattice QCD: computational strategy and a first physical result

技術要約：格子QCDによる包括的 $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ 崩壊

問題の定義
CKM（Cabibbo–Kobayashi–Maskawa）行列要素 $|V_{cb}|$ の決定は、フレーバー物理学における中心的な課題であり、現在は、排他的（exclusive）な決定と包括的（inclusive）な決定との間の顕著な不一致によって特徴付けられている。包括的な $B$ メソンの半レプトン崩壊は、演算子積展開（OPE）を介して理論的に堅牢な経路を提供するが、$\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ の包括的崩壊率 $\Gamma$ を第一原理に基づく格子QCDで計算することは、歴史的に「二つのスケールの問題」（$a m_b \ll 1$ と $\Lambda_{QCD} L \gg 1$ を同時に満たすこと）および、重いクォーク質量が増大するにつれて四点相関関数の信号対雑音比（S/N比）が急速に悪化することから、困難であると考えられてきた。

手法および計算戦略
著者らは、物理的な軽、ストレンジ、およびチャームクォーク質量を持つ $N_f = 2+1+1$ 拡張ツイスト質量（ETMC）ゲージ・アンサンブルを用いて、この包括的崩壊率を第一原理から計算するための新しい戦略を提示している。その手法は、主に以下の3つの柱に基づいている。

1. OPEによる質量補間： 物理的な $b$ クォーク質量での直接的なシミュレーションは（カットオフ効果のため）計算コストが極めて高いため、著者らは物理的な $b$ クォーク質量よりも軽い重いクォーク質量（$M_I < M_{\bar{B}_s}$）でシミュレーションを行う。この際、位相空間の体積を一定に保つために、初期メソンと最終メソンの質量の比（$r_{IF} = M_F/M_I = M_{D_s}/M_{\bar{B}_s}$）を固定する。物理的な結果は、非摂動的な格子データ点を、静的極限（$m_b \to \infty$）において厳密となるOPE予測を用いて、物理点へと補間することによって得られる。これにより、外挿の問題が補間の問題へと変換される。

2. 四点相関関数に対する新しい数値戦略： アルゴリズムにおける重要な革新は、必要な四点ユークリッド相関関数を計算するための新しい戦略の導入である。著者らは、「旧来の戦略」（ソースとシンクの演算子間の時間間隔を固定する方法）が、重いメソンの信号対雑音比の壊滅的な喪失を招くことを指摘している。新しい戦略では、弱電流の挿入時刻とソース／シンク演算子の間の時間間隔（$t_{snk} - t = 0 - t_{src}$）を固定しつつ、電流挿入時刻 $t$ を変化させる。このアプローチは、計算量は増えるものの（より多くの逆行列計算を必要とする）、信号対雑音比の問題を克服することに成功しており、旧来の手法では失敗した大きな重いメソン質量における相関関数の解像を可能にしている。

3. スペクトル再構成と解析： 崩壊率は、切断された（amputated）ユークリッド相関関数からHansen–Lupo–Tantalo (HLT) アルゴリズムを用いて抽出される。著者らは、極限の順序（連続体極限 $a \to 0$、スミアリング $\sigma \to 0$、および数値積分）を変化させ、また物理点への2つの異なる補間アプローチを採用することで、系統誤差を厳密に分析している：

   • フィッティング関数アプローチ： 非摂動的な寄与に対する有理関数と、摂動的なOPE項を組み合わせる手法。
   • 関数空間アプローチ： 摂動的なOPEの結果によってアンカーされた事前分布を用いたガウス過程を用いる手法。

主要な結果
計算は、3つのETMCアンサンブル（B64, C80, D96）を用いて行われた。これらは格子間隔 $a \approx 0.08, 0.07, 0.6$ fm および空間体積 $L \approx 5.5$ fm を持つ。著者らは、包括的崩壊率の予備的な決定値を得た：
$$\Gamma[\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}] = 2.47(13)(10)[17] \times 10^{-11} \text{ GeV}$$
ここで、誤差はそれぞれ統計誤差、系統誤差（解析ステップに由来）、および合計誤差を表す。合計相対誤差は約7%である。

この結果を用い、$X_{\bar{s}u}$ の寄与を無視すると（$|V_{ub}|^2/|V_{cb}|^2$ による抑制により正当化される）、著者らは $|V_{cb}|$ の値を導出する：
$$|V_{cb}|_{\text{latt-incl}} = 40.9(1.2)_{\text{latt}}(1.7)_{\text{exp}}[2.1]_{\text{tot}} \times 10^{-3}$$
この値は、排他的な決定（$39.62(47) \times 10^{-3}$）およびOPEに基づく包括的な決定（$41.64(47) \times 10^{-3}$）の両方と整合しているが、現在の統計精度は、既存のこれら手法間の緊張関係を解決するには不十分である。

意義と主張
本論文は、以下の分野における重要な貢献を主張している：

• アルゴリズムの突破口： 四点相関関数に対する新しい数値戦略は、重いメソンに関する長年の信号対雑音比の問題を解決し、格子上での包括的 $B$ 崩壊の研究を可能にした。
• 手法の妥当性確認： 本研究は、非摂動的データとOPEの制約を組み合わせることにより、包括的な半レプトン崩壊率が、現象論的に関連のある精度で、第一原理に基づく格子QCDから計算可能であることを示している。
• 独立した検証： 異なる格子離散化（ツイスト質量）および解析手法（HLTアルゴリズム）を用いることで、本研究は、同時並行的に進められている他の取り組み（例：文献[15]）とは完全に独立した $\bar{B}_s \to X_{\bar{s}c}\ell\bar{\nu}$ の第一原理に基づく結果を提供している。
• 今後の展望： 著者らは、現在の7%の誤差が、限定的な格子アンサンブルの数と、データ駆動型の有限体積研究の欠如に起因することを強調している。彼らは、より細かい格子間隔（$a \approx 0.05$ fm）やより大きな体積を追加することで、合計誤差を大幅に削減でき、結果として $|V_{cb}|$ の謎を決定的に解決できる可能性があると主張している。

著者らは、現在の結果については「最初の物理的な結果」であるとし、有限体積効果および連続体極限に関連する系統誤差には、今後の出版物においてさらなる洗練が必要であると述べており、控えめなトーンを維持している。