Amplitude analysis and branching fraction measurement of the decay $D^0 \to K^+K^-\pi^0\pi^0$

BESIII 実験により、$D^0 \to K^+ K^- \pi^0 \pi^0$ 崩壊の絶対分岐比が初めて測定され、その振幅解析から $K^{*}(892)^+ K^{*}(892)^-$ 中間過程が支配的であり、S 波が優位であることが明らかになりました。

要約

この論文は、素粒子物理学の「BESIII 実験」というチームが、「D0 メソン」という小さな粒子がどのように崩壊するかを詳しく調べた報告書です。

専門用語を並べると難しく聞こえますが、実は**「粒子の分解能（分解）と、その中身（成分）の割合」**を調べる面白い実験なのです。

以下に、誰でもわかるような比喩を使って解説します。

---

1. 実験の舞台：「粒子の工場」と「分解機」

まず、実験が行われた場所を想像してください。
北京にある**「BEPCII」**という巨大な加速器は、電子と陽電子（電子の反物質）を光速近くまで加速して衝突させる「粒子の衝突工場」です。

• 衝突の瞬間： 2 つの粒子が激しくぶつかり合うと、エネルギーが熱に変わり、新しい粒子が生まれます。この実験では、**「D0 メソン」**という、少し重い「チャームクォーク」という部品を持った粒子が大量に作られました。
• BESIII 検出器： この工場には、**「BESIII」**という巨大なカメラ（検出器）が設置されています。これは、衝突して生まれた粒子がどう飛び散るかを、360 度すべて撮影する超高精細なカメラです。

2. 調べる対象：「魔法の箱」の分解

今回の実験の目的は、**「D0 メソン」という粒子が、「K+ K- π0 π0」**という 4 つの別の粒子に崩壊する過程を詳しく調べることでした。

これを**「魔法の箱」**に例えてみましょう。

• D0 メソン = 中身が見えない「魔法の箱」。
• 崩壊 = その箱が突然開いて、中から 4 つの異なるおもちゃ（K+, K-, π0, π0）が飛び出すこと。

しかし、問題は**「箱が開く瞬間に、いったいどんな手順で中身が出てきたか」**がわからないことです。
いきなり 4 つが飛び出すのか？
それとも、まず 2 つの「中箱」が出てきて、その中箱がさらに割れて 4 つになったのか？

3. 振幅解析：「中箱」の正体を当てるゲーム

この論文で行われた**「振幅解析（Amplitude Analysis）」とは、まさに「飛び出したおもちゃの動きを分析して、箱が開くプロセス（中箱の正体）を推測する」**というゲームです。

研究者たちは、飛び出した 4 つの粒子の動き（角度やエネルギー）を精密に測定し、以下の仮説を検証しました。

• 仮説 A： 箱から直接 4 つが出る。
• 仮説 B： まず「K*」という中箱と「K*」という中箱が出てきて、それぞれがさらに割れて 4 つになる。
• 仮説 C： 「f1」や「η」という別の中箱が出てくるパターン。

結果：
データを見つめると、**「仮説 B（K* メソンが 2 つ出てくるパターン）」*が圧倒的に多いことがわかりました。
まるで、魔法の箱から「K メソン」という中箱が 2 つ飛び出し、それがさらに割れて最終的な 4 つの粒子になったというストーリーが最も有力だったのです。

4. 重要な発見：「振動する方向」の謎

さらに面白い発見がありました。
「K* メソン」という中箱は、単なる球体ではなく、**「振動している（回転している）」性質を持っています。これを「偏光（ポーラライゼーション）」**と呼びます。

• 縦に振動する（S 波）： 軸に沿って振動する。
• 横に振動する（P 波など）： 軸に対して横に振動する。

これまでの理論（予測）では、「横に振動する方が多いはずだ」と言われていました。
しかし、今回の実験結果は**「縦に振動する（S 波）方が圧倒的に多い」**という、予想とは異なる結果を出しました。

比喩：
「風船が割れるとき、理論的には横に裂けるはずだと言われているのに、実際には縦に裂けることが多い！」という発見です。
これは、粒子の間の「見えない力（強い相互作用）」が、理論家が思っていたよりも複雑に働いていることを示唆しています。

5. 頻度の測定：「確率」の計算

最後に、この「K* メソン 2 つが出てくるパターン」が、全体の崩壊の中で**何％**を占めるかを計算しました。

• 全体の確率： D0 メソンが K+K-π0π0 に崩壊する確率は、約 0.07%（1000 回に 0.7 回）でした。
• 中身のパターン： その中で、K* メソンが 2 つ出てくるパターンが最も多く、全体の約 40% を占めていました。

これにより、理論モデル（予測）と実験結果を比較することができ、より正確な「粒子の崩壊の地図」が完成しました。

まとめ：なぜこれが重要なのか？

この研究は、単に「粒子がどう崩壊したか」を知るだけでなく、**「自然界の根本的な力（強い力）がどう働いているか」**を理解するための重要なピースです。

• 理論との対決： 既存の理論（予測）と実験結果が一致しない部分（縦振動の多さなど）が見つかりました。これは、物理学者たちが「新しい理論」や「修正」を迫られるきっかけになります。
• CP 対称性の破れ： 将来、物質と反物質の不对称（なぜ宇宙に物質が多いのか）を理解する手がかりになる可能性もあります。

つまり、この論文は**「小さな粒子の箱を開けて、その中身がどう動いたかを詳しく記録し、宇宙の法則を解き明かすための重要なデータ」**を残したという報告なのです。

技術概要

BESIII 実験 collaboration による論文「Amplitude analysis and branching fraction measurement of the decay $D^0 \to K^+K^-\pi^0\pi^0$」の技術的サマリーを以下に日本語で提示します。

1. 研究の背景と課題 (Problem)

チャーム中間子（$D$ メソン）のハドロン崩壊の理論的研究は、チャームクォークの質量が「重いクォーク展開」には重すぎず、「カイラル摂動論」には軽すぎないという中間的な領域にあるため、非常に困難です。この課題に対処するためには、非摂動的な手法を用いたモデルを精密な実験データで制約し、理論予測を検証する必要があります。

特に、$D^0 \to K^+K^-\pi^0\pi^0$ という単一カビボ抑制（SCS）崩壊は、中間過程として $D \to VV$（ベクトルメソン対）や $D \to AP$（軸性ベクトルメソンと擬スカラーメソン）などの複雑な干渉を含むため、その崩壊機構（カラー許容外部 W 放出、カラー抑制内部 W 放出、W 交換など）を解明する上で重要です。しかし、この崩壊の振幅分析はこれまで行われておらず、中間過程の寄与や分岐比、特に $D^0 \to K^*(892)^+K^*(892)^-$ 過程の偏極状態に関する実験的知見は不足していました。

2. 手法 (Methodology)

本研究では、北京電子陽電子衝突型加速器（BEPCII）で生成された $\sqrt{s} = 3.773$ GeV の $e^+e^-$ 衝突データ（積分光度 20.3 fb$^{-1}$）を BESIII 検出器で収集・解析しました。

• イベント選択と信号抽出:

  • ダブルタグ（DT）手法を採用し、片方の $D^0$ をシグナルモード（$K^+K^-\pi^0\pi^0$）、もう片方の $\bar{D}^0$ を既知の単一タグ（ST）モード（$K^+\pi^-$, $K^+\pi^-\pi^0$, $K^+\pi^-\pi^-\pi^+$）で再構成しました。
  • 運動量保存則や既知の質量値（$D^0$, $\pi^0$）を制約条件とした全拘束キネマティックフィットを行い、信号純度を向上させました。
  • 最終的に 791 のイベントがシグナル領域で抽出され、信号純度は 94.2% と推定されました。

• 振幅分析 (Amplitude Analysis):

  • 非束縛最大尤度法を用いた振幅分析を実施しました。
  • 中間共鳴状態を「アイソバールモデル」で記述し、各共鳴の複素係数（大きさ $\rho_n$ と位相 $\phi_n$）を自由パラメータとしてフィットしました。
  • 共鳴状態には $K^*(892)$, $f_1(1420)$, $\eta(1475)$, $\eta(1405)$, $K_1(1400)$, $a_0(980)$ などを考慮し、スピン因子、ブロット・ワイズバックの障壁因子、相対論的 Breit-Wigner 関数（または LASS パラメータ化）を用いて振幅を構築しました。
  • 背景事象の分布は、包括的モンテカルロ（MC）サンプルと XGBoost 機械学習アルゴリズムを用いてモデル化し、シグナルと背景の区別を高精度に行いました。

• 分岐比測定:

  • DT 手法に基づき、シグナル事象数、タグ事象数、検出効率、および量子相関補正パラメータを用いて絶対分岐比を算出しました。

3. 主要な貢献と結果 (Key Contributions and Results)

A. 分岐比の精密測定

• $D^0 \to K^+K^-\pi^0\pi^0$ の絶対分岐比を初めて高精度で測定しました。
  • 結果: $(0.73 \pm 0.03_{\text{stat}} \pm 0.01_{\text{syst}}) \times 10^{-3}$
  • これは以前の BESIII の測定値（$0.69 \pm 0.07 \times 10^{-3}$）と 1$\sigma$ 以内で一致しつつ、精度を 2.5 倍向上させたものです。

B. 支配的な中間過程の同定と分岐比

• 振幅分析の結果、支配的な中間過程は $D^0 \to K^*(892)^+K^*(892)^-$ であることが確認されました。
• この過程の分岐比は $(2.79 \pm 0.13_{\text{stat}} \pm 0.11_{\text{syst}}) \times 10^{-3}$ と測定されました。
• 他の主要な寄与として、$D^0 \to \eta(1475)\pi^0$ や $D^0 \to \eta(1405)\pi^0$ などが観測されました。

C. 偏極状態の測定（重要な発見）

• $D^0 \to K^*(892)^+K^*(892)^-$ 過程において、縦偏極分率（Longitudinal Polarization Fraction, $F_L$）を初めて測定しました。
  • 結果: $F_L = 0.468 \pm 0.046_{\text{stat}} \pm 0.011_{\text{syst}}$
• この値は、既存の非相対論的構成クォークモデル（FSI を含む・含まず）が予測する「横偏極優位（$F_L \approx 0.3$）」とは異なり、縦偏極が支配的であることを示しています。
• 測定値は理論予測（Ref. [4]）よりも 3$\sigma$ 以上高い値であり、現在の理論モデルの修正や、最終状態相互作用（FSI）のダイナミクスに関する新たな洞察を必要としています。

D. 振幅の相対的な大きさと位相

• 各中間過程の相対的な大きさ（Magnitude）と位相（Phase）を初めて決定しました。
• $D^0 \to K^*(892)^+K^*(892)^-$ 過程は S 波が支配的であることが明らかになりました。
• 干渉項を含むフィット分数（Fit Fractions）を算出し、全過程の寄与を定量化しました。

4. 意義 (Significance)

本研究は、$D^0 \to K^+K^-\pi^0\pi^0$ 崩壊に関する最初の包括的な振幅分析であり、以下の点で極めて重要です。

1. 理論モデルの検証: 測定された分岐比と偏極状態は、既存のハドロン崩壊モデル（NRCQM, SU(3) 対称性モデルなど）に対して厳しい制約を与えます。特に、$F_L$ の測定値と理論予測の大きな乖離は、チャームセクターにおけるハドロンダイナミクスや最終状態相互作用（FSI）の理解を深めるための重要な手がかりとなります。
2. CP 対称性の破れ研究への寄与: $D$ メソンの CP 対称性の破れを研究する際、崩壊振幅の位相と大きさを正確に知ることは不可欠です。本分析で得られた振幅パラメータは、将来の CP 対称性の破れ探索の基礎データとして機能します。
3. 中間共鳴状態の特性解明: SCS 崩壊における軸性ベクトルメソン（$K_1$, $f_1$ など）やスカラーメソンの寄与を定量的に評価し、これらの粒子の性質や混合角に関する情報を提供しました。

総じて、この研究はチャーム物理の分野において、実験データに基づく精密な振幅分析の新たな基準を示し、理論と実験の間のギャップを埋める重要な一歩となりました。