Radiative Corrections to the Direct Detection of Inelastic Scattering of Higgsino-like Neutralino Dark Matter

技術要約：ヒッグシーノ的なニュートラリーノの非弾性散乱の直接検出に対する放射補正

問題提起
Rパリティが保存される最小超対称標準模型（MSSM）は、最も軽いニュートラリーノ（$\tilde{\chi}^0_1$）を生存可能なダークマター（DM）候補として提供する。ほぼ純粋なヒッグシーノ的な $\tilde{\chi}^0_1$ は、自然性の問題の解決および効率的な対消滅による観測されたDM存在量（$\Omega_{DM}h^2 \simeq 0.12$）の適合において特に有望である。しかし、このようなヒッグシーノ的なDMの直接検出（DD）は困難である。なぜなら、ヒッグス粒子および$Z$ボソンへのツリーレベルの結合が強く抑制されており、スピン非依存（SI）の散乱断面積が極めて小さくなるためである。

弾性散乱（$\tilde{\chi}^0_1 q \to \tilde{\chi}^0_1 q$）については次世代（NLO）電弱（EW）補正が計算されているが、非弾性散乱プロセス（$\tilde{\chi}^0_1 q \to \tilde{\chi}^0_2 q$）は、しばしば見落とされる重要な側面である。圧縮されたスペクトルのシナリオでは、最も軽いニュートラリーノ（LSP）と次に軽いニュートラリーノ（NLSP）の間の質量分裂（$\delta = m_{\tilde{\chi}^0_2} - m_{\tilde{\chi}^0_1}$）は、極めて小さく（$\sim \mathcal{O}(10\text{--}100)$ keV）なることがある。典型的な運動エネルギー $\sim \mathcal{O}(100)$ keV を持つ銀河系DMにとって、これらのほぼ縮退した重い状態への遷移は、動力学的にアクセス可能となる。LUX-ZEPLIN（LZ）実験による248 keV付近の核反跳イベントのような最近の観測は、非弾性チャネルが支配的になるか、あるいは総イベント率に大きく寄与する可能性があるTeVスケールのヒッグシーノDMの再評価を促している。

手法
著者らは、完全な1ループ電弱放射補正を含む、弾性および非弾性チャネルの両方を組み込んだSI直接検出断面積の包括的な計算を行う。その手法は以下の通りである：

1. 繰り込みスキーム： 著者らは、シャージノ・ニュートラリーノ・セクターに対してオンシェル繰り込みスキームを採用している。これは、非弾性散乱率が質量分裂 $\delta$ の精密な値に強く依存するため、極めて重要である。このスキームは、物理的なニュートラリーノおよびシャージノの質量が1ループレベルで正しく定義されることを保証する。
2. 頂点補正： 著者らは、3点頂点 $\tilde{\chi}^0_1 \tilde{\chi}^0_i h_j$ （ここで $i, j \in \{1, 2\}$ および $h_j$ は軽いおよび重いヒッグス粒子）に対する1ループ補正を計算している。これには以下が含まれる：
   • フェルミオン（$F$）、スカラー（$S$）、およびベクトルボソン（$V$）を含むループ図。具体的には、クォーク・スクワークのループ（トップ・ストップ・ループが支配的）や、シャージノ-$W$ループ、およびニュートラリーノ-$Z$ループに焦点を当てている。
   • 波動関数繰り込みおよび質量繰り込みから導出されるカウンタータームの寄与。
3. 計算フレームワーク：
   • MSSMスペクトルを生成し、2ループ・ヒッグス質量および$\overline{DR}$スキームにおける1ループ質量を計算するために SPheno と SARAH を使用する。
   • 頂点の振幅を生成し、それらをパスリーノ・ヴェルトマン関数に簡約するために FeynArts と FormCalc を用いる。
   • 数値評価のために LoopTools を使用する。
   • 補正された頂点と質量を micrOMEGAs (v6.3.0) に実装し、弾性散乱断面積を計算する。
   • 非弾性散乱振幅は、ヒッグス媒介のスカラー演算子のみを考慮した（twist-2やその他の副次的な寄与は無視した）有効ラグランジアン・アプローチに基づくカスタムルーチンを用いて計算される。
4. イベント率の計算： LZ実験（露光量 4.2 tonne-yrs）に対する微分イベント率を、弾性チャネルと非弾性チャネルの両方の寄与を合算して計算する。解析では、$M_1$-$M_2$ 平面を様々なヒッグシーノ質量パラメータ（$\mu = \pm 0.5, \pm 1$ TeV）に対してスキャンし、$\delta \lesssim 100$ keV となる領域を特定する。

主な貢献

• 非弾性チャネルの包含： 弾性散乱に焦点を当てるか、あるいは $\tilde{\chi}^0_1 \tilde{\chi}^0_2 h_j$ 頂点補正を無視していた先行研究とは異なり、本研究では放射補正された非弾性散乱断面積を明示的に計算している。
• オンシェル質量分裂： 著者らは、ニュートラリーノとシャージノの1ループ補正質量を計算するためにオンシェル繰り込みスキームを一貫して適用しており、これにより、非弾性動力学を支配するパラメータである質量分裂 $\delta$ の精密な決定を保証している。
• 補正の相互作用： 論文は、NLO補正が弾性と非弾性の相対的な断面積の大きさを劇的に変え得ることを示している。ヒッグシーノ的なDMに対するツリーレベルの結合は抑制されているが、ループ誘起の補正はそれらを数桁オーダーで増大させ得る。特定のシナリオでは、非弾性断面積が弾性断面積を上回ることもあるが、最終的なイベント率は動力学的抑制に依存する。
• ベンチマーク解析： 質量分裂と放射補正の異なるレジームを例示するために、5つのベンチマークポイント（BP）が選択され、LO（Leading Order）およびNLO予測がLZの制約と比較されている。

結果

• 補正の大きさ： 頂点因子（$\tilde{\chi}^0_1 \tilde{\chi}^0_i h_j$）に対するNLO補正は、特にLOの値が小さい場合、$\sim 10\%$ から特定のベンチマークポイントでは $500\%$ 以上に達することもあり、非常に大きなものとなる。
• 断面積 vs イベント率： 特定のパラメータ領域（例：ベンチマークポイント2）では、NLOにおける非弾性散乱断面積が弾性散歴断面積を上回ることがある。しかし、イベント率は、重い状態を励起するために必要な最小速度（$v_{min}^{inel}$）によるボルツマン抑制の影響を受ける。BP2においては、NLOでの非弾性散乱断面積が弾性NLO断面積と同等か、あるいはわずかに大きいものの、動力学的抑制のため、総イベント率への非弾性チャネルの寄与は弾性チャネルよりもわずかに低くなる。他のポイント（例：BP1）では、非弾性率は弾性率に対して強く抑制されている。
• ボルツマン抑制： 断面積が大きくなる可能性があるにもかかわらず、非弾性イベント率は、質量分裂 $\delta$ を克服するために高い最小DM速度を必要とするため、弾性率に対して減少することが多い。
• LZ制約： 検討されたすべてのベンチマークポイントは、90%信頼区間において現在のLZ実験のヌル結果と一致している。しかし、非弾性チャネルとNLO補正を含めることで、予測されるイベント率が変化し、排除限界や将来の検出に必要な露光量に影響を与える。
• ループ寄与： 放射補正の主要な寄与はトップ・スクワーク・ループ（$\tilde{t}q\tilde{q}$ 型）であり、次いでシャージノ-$W$ループである。重いヒッグス（$h_2$）の寄与は、その大きな質量（$> 6$ TeV）により抑制されている。

意義と主張
本論文は、特に質量分裂が小さい場合、ヒッグシーノ的なDMの直接検出に関する精密な理論的取り扱いは、弾性と非弾性の両方の散乱チャネルを考慮しない限り不完全であると主張している。著者らは以下を強調している：

1. 放射補正は無視できない： これらは散乱率を大幅に増大させ、LZ、XENONnT、PandaXのような現在および将来の実験の射程内にヒッグシーノ的なDMを捉える可能性がある。
2. 非弾性断面積が弾性断面積を上回り得る： 圧縮されたスペクトルの特定のシナリオでは、非弾性散乱成分が弾性成分を上回ることがある。しかし、ボルツマン抑制のため、総イベント率は依然として弾性チャネルによって支配される可能性があり、断面積と動力学的因子の両方を正確に計算することの重要性が浮き彫りになる。
3. 最近のデータへの関連性： 本フレームワークは、最近のLZ核反跳イベント（例：248 keVイベント）の解釈、およびTeVスケールのヒッグシーノDM解釈の妥当性を評価するために直接的に関連している。
4. 自然なSUSYのプローブ： これらの研究は、圧縮されたスペクトルと抑制されたツリーレベル相互作用を特徴とし、他の手段ではアクセスが困難な「自然な」超対称的シナリオを探索するための補完的な経路を提供する。

著者らは、将来の直接検出実験は、放射補正を含む精密な理論予測と組み合わされることで、ヒッグシーノ的なDMの発見の可能性を堅牢に評価するために不可欠であると結論付けている。