A determination of the backscattering probability of low-energy antiprotons

技術要約：低エネルギー反陽子の後方散乱確率の決定

問題提起
GBAR（静止状態における反物質の重力挙動）実験において、反水素（$\bar{H}$）の生成断面積は、反陽子（$\bar{p}$）がポジトニウム雲との電荷交換反応にどれだけ関与するかを測定することによって決定される。この測定における極めて重要な入力値は、反応キャビティに入射する反陽子の精密なフラックスである。このフラックスは通常、「パススルー（通過）」ランを用いて定量化される。このランでは、反陽子がキャビティを通過し、下流のマイクロチャネルプレート（MCP）検出器（MCPH）上で消滅する。発生した消滅生成物（荷電パイ中間子）は補完型金属酸化物半導体（CMOS）センサによって検出され、較正を通じて入射反陽子数が推論される。

重大な系統誤差として、低エネルギーの反陽子（4–6 keV）がMCP表面に衝突した際、直ちに消滅せずに、材料から大きな角度でクーロン散乱を受け、後方散乱（バックスキャッタリング）される可能性が挙げられる。もし後方散乱が発生した場合、反陽子はMCP上でではなく、上流（例：ビームパイプ上）で消滅する可能性がある。CMOSセンサはMCP上での消滅による生成物を検出するように配置されているため、これらの後方散乱されたイベントに対して「盲目」になる可能性がある。電子、陽電子、陽子の後方散乱確率はよく研究されているが、低陽子の後方散乱確率は大部分が未知であり、クーロン散乱、エネルギー損失、および反陽子特有の消滅の間の特有の競合状態があるため、他の粒子から直接推論することもできない。過去の報告では、大きな後方散乱分率が存在する可能性が示唆されていたが、もしそれが事実であれば、反陽子フラックスの較正および導出される$\bar{H}$生成断面積に対して補正が必要となる。

手法
この不確実性に対処するため、著者らはGBARビームライン内の専用の実験セットアップを用いた、インサイチュ（in-situ）による後方散乱分率の測定を実施した。

1. 実験セットアップ： 実験では、ELENA減速器の下流に配置された可動式MCP（MCPup）を利用した。通常は反応チャンバーのMCPHの上にあるCMOSセンサを、一時的にMCPupの上方に移動させた。センサはレールに取り付けられており、ビーム軸からの距離（$R$）を9.3 cmから94.3 cmの範囲で変化させることが可能である。データ収集は、2025年5月および6月のランにおいて、ビームエネルギー4 keV、6 keV、および100 keVで行われた。
2. 測定原理： 本手法の核心は、CMOSセンサにおける信号の幾何学的依存性に依拠している。
   • MCP上で直接消滅する反陽子は、幾何学的受容性（geometric acceptance）に従い、$1/R^2$に比例してスケールする確率 $P_{annihil}(R)$ でCMOSセンサに信号を与える。
   • MCPで後方散乱された後に、MCPの上流にあるビームパイプ上で消滅する反陽子は、信号を与える確率 $P_{back}(R)$ を持つ。シミュレーションによれば、消滅点は上流に偏位している（平均して約9 cm）ため、直接消滅の場合と比較して、特に小さな$R$において幾何学的受容性が異なる。
   • $R$の関数としてのCMOSセンサ内のクラスター数の変化を測定し、それをシミュレーションと比較することで、後方散乱された反陽子の分率（$f_{back}$）を抽出することができる。
3. シミュレーション（GEANT4）： ビームラインの幾何学、MCPupのアセンブリ（支持リングやコーティングを含む）、およびCMOSセンサをモデル化するために、詳細なGEOT4（バージョン11.3）シミュレーションが行われた。
   • シミュレーションでは、反陽子の消滅のためにINCLXX物理リストを、電磁多重散乱のためにWentzelVIモデルを使用した。
   • 絶対的な後方散乱分率の予測自体はモデル依存性が高いことが判明しているが、後方散乱粒子の「角度分布」（これにより粒子がどこにビームパイプを叩くかが決まる）のシミュレーションを堅牢にするために、散乱角の制限（$\theta_{lim}$）やレンジカットなどのパラメータを調整した。
   • シミュレーションは、様々なセンサオフセットおよびビームエネルギーに対する確率 $P_{annihil}(R)$ および $P_{back}(R)$ を予測した。

主な貢献と結果

1. 後方散乱の上限値： 実験データ（クラスター数 vs $1/R^2$）は、4 keV/100 keVおよび6 keV/100 keVのデータセットに対する結合フィットを用いて、GEANT4の予測値にフィットされた。フィットパラメータには、後方散乱分率、正規化係数、および幾何学的アライメントパラメータ（縦方向のシフトおよびレールの角度）が含まれた。
   • 結果として、4 keVまたは6 keVにおける後方散乱の証拠は見られなかった。最適フィットによる後方散乱分率は、両方のケースでゼロであった。
   • 4 keVにおいて14%、6 keVにおいて11.5%の後方散乱確率の上限値が、68%の信頼区間で設定された。
2. CMOS検出器の較正： 本研究は、シミュレーションが信号の$R$依存性を正しく記述していることを確認した。測定されたクラスター数とELENAビーム強度モニターを比較することで、較正係数が導出された。
   • 解析の結果、シミュレーションは絶対的なクラスター数を約7.4%過小評価していることが判明した（これはモンテカルロの正規化または検出器効率の不確実性に起因する）。
   • 未知の分率が0–20%の範囲であることを想定し、保守的な範囲として、測定場所（MCPup）と生成場所（MCPH）の間での後方散乱分率の違いを考慮し、CMOS較正係数に5%の系統誤差を割り当てた。
3. 系統誤差の評価： 本論文は、反陽子フラックスの決定に対する後方散乱の影響を定量化している。仮に20%の後方散達分率があったとしても、較正係数に生じる変動は5%に限定される。入射反陽子数の全系統誤差は約8.4%と推定されており、これは後方散乱よりもむしろビーム位置の不確かさや強度モニターの誤差に支配されている。

意義
本論文は、低エネルギーの反陽子がMCP表面で即時消滅するという仮定が、GBAR実験で要求される精度内で妥当であることを確立した。後方散乱確率に対して厳格な上限値を設定することで、本研究は、反陽子生成断面積の決定における主要な潜在的系統誤差の源を除去した。本研究は、後方散乱分率が無視できるほど小さく、したがって反陽子フラックス測定の精度を著しく低下させることはなく、それによって反物質の重力挙動に関するGBARの結果の信頼性を裏付けている。距離依存的な信号変化を用いて後方散乱効果を分離するという手法は、将来の低エネルギー反陽子実験のための堅牢なテクニックとして提示されている。