Orientation-dependent sensitivity of the all-sky all-frequency gravitational-wave radiometer search to continuous-wave sources

技術要約：ASAF重力波ラディオメーター探索における方位依存の感度

問題提起
全天・全周波数（ASAF）重力波ラディオメーター探索は、検出器のペアから得られるデータを相互相関させることで、変形した中性子星などの持続的な狭帯域信号（連続波：CW）を検出するように設計されている。標準的なASAF統計量は、$+$ および $\times$ のテンソル偏極に対して等しい重みを割り当てる。この重み付けは、非偏極の確率的重力波背景放射（SGWB）に対しては最適であるが、CW源に対しては最適ではない。CW源は、その傾斜角 $\iota$ と偏極角 $\psi$ によって決定される偏極振幅で放射される。そのため、光源がフェイスオン（$\iota=0$）でない限り、等重みの統計量は、同等の固有振幅を持つフェイスオンの光源が生成するよりも少ないパワーしか測定できない。これにより、「方位バイアス」が生じ、探索は光源の固有の強さを過小評価し、信号対雑音比（SNR）を低下させる可能性がある。これは、実際の信号が検出閾値を下回る原因となり得る。

手法
著者らは、このバイアスの光源の向きへの正確な依存性を導出し、それをASAF解析パイプラインを通じて伝播させている。

1. 投影係数 ($\mu^2$) の導出:
   論文では、投影係数 $\mu^2(\iota, \psi)$ を、ASAF統計量によって測定される実効歪みパワーと、CW源の固有の歪みパワー ($h_0^2$) との比として定義している。

   • 著者らは、$\mu^2(\iota, \psi)$ の閉形式の表現（式33）を導出し、それが傾斜角のみに依存する項と、偏極角 $\psi$ で単一調和関数として変調する項から構成されることを示している。
   • 変調度 $\epsilon$ は、ベースラインの幾何学的配置と光源の赤緯によって決定される。
   • 著者らは、この係数の全方位に対する等方的な平均が、空の位置、ベースラインの幾何学的配置、ノイズの重み付け、または周波数に依存せず、正確に $2/5$ (0.4) であることを証明している。

2. 探索量への伝播:
   この定式化は、ASAF探索の主要な出力に適用される。

   • 歪み振幅: 固有の振幅 $h_0$ は、報告される実効歪み $h_{\text{eff}}$ を用いて $h_0 = h_{\text{eff}} / \sqrt{\mu^2}$ と関係付けられる。
   • 信号対雑音比 (SNR): 特定の向きにある光源に対する回収されたSNRは、$\rho_{\text{rec}} = \mu^2 \rho_{\text{ref}}$ である。ここで $\rho_{\text{ref}}$ は、同じ固有強度を持つフェイスオンの光源のSNRである。
   • 上限値 (ULs): ベイズ的な上限値 $h_0$ は、未知の向きのパラメータ（$\iota, \psi$）について尤度を周辺化することによって導出される。その際、標準的な円偏極の仮定（$\mu^2=1$）を $\mu^2$ の分布に置き換える。

3. 検証と適用:

   • O2 ハードウェア注入: 本定式化は、既知のパラメータと多様な向き（ほぼフェイスオンからほぼエッジオンまで）を持つ、第2観測ラン（O2）からの4つのシミュレーション・パルサー信号（ハードウェア注入）に対してテストされている。
   • O4a 候補: 本手法は、フォローアップ調査のために選択された、第4観測ランの第1部分（O4a）からの505個のサブスレッショルド（閾値下）候補に適用されている。

主な結果

• 解析的証明: 論文は、完全なASAF統計量において平均投影係数が任意の空の位置で正確に $2/5$ であることを確認している。また、$\psi$ に対する残差依存性の正確な閉形式を提供し、最も不利な向きと、与えられた空の位置に対する最大検出可能傾斜角を特定している。
• ハードウェア注入解析:
  • ほぼエッジオンの注入（$\iota \approx 94.6^\circ$）の場合、測定されたSNRは5.20であり、検出閾値の6.19を下回っていた。しかし、投影係数（$\mu^2 \approx 0.167$）で補正すると、同じ信号がフェイスオンで観測された場合、SNRは $31.1 \pm 6.0$ になることが明らかになった。
  • 向きに関して周辺化された、固有の歪み振幅 $h_0$ のベイズ上限値は、標準的なフェイスオン（円偏極）の上限値よりも2.7倍から2.9倍大きかった。
  • 解析は、$\mu^2$ による補正および既知のスペクトル損失（ドップラー効果およびピクセル位相効果）を考慮することで、注入された振幅を正常に回収した。
• O4a 候補解析:
  • 505個のサブスレッショルド候補について、周辺化された上限値は、公表されている円偏極の上限値と比較して2.5倍から2.9倍増加した。
  • もし候補が実際のCW源であるならば、その傾斜角がフェイスオンから $34.7^\circ$ から $50.3^\circ$ の範囲内（特定の候補のSNRに依存）にある場合にのみ、検出閾値に達することになる。これは、等方的な空のわずか18%から36%の向きの範囲に対応する。
  • 「ワーストケース」の上限値（最も不利な向きを想定した場合）は、円偏極の上限値の最大で $\sim 3.6$ 倍に達する。

意義と主張
本論文は、向きのバイアス $\mu^2$ を考慮することが、候補の統計的有意性を変えることなく、CW源に関するASAFの結果の解釈を根本的に変えるものであると主張している。

1. 非検出の再解釈: 低い測定SNRは、必ずしも固有の光源が弱いことを意味するわけではない。それは、不利な向きにある強い光源であることを示している可能性がある。あるSNRの背後にある固有の強さは、フェイスオンの推定値よりも $1/\mu^2_{\text{min}} \approx 8$ から $13$ 倍大きくなる可能性がある。
2. 報告基準: 著者らは、将来のASAF解析において、上限値を単一のフェイスオン値としてではなく、向きの範囲（または周辺化された値）として報告することを推奨している。また、測定されたSNRおよび実効歪みとともに、最大検出可能傾斜角（$\iota_{\text{max}}$）および向きの範囲（$C_{\text{orient}}$）を併せて報告することも提案している。
3. 選択バイアス: 測定されたSNRのみに基づいて候補を選択することは、フェイスオンの光源へのバイアスを生じさせる。フォローアップ探索においては、サブスレッショルドの候補が高傾斜角の強い光源である可能性があることを考慮すべきである。

論文は、候補の測定されたSNRおよびp値自体は変化しないものの、歪み振幅および検出確率の物理的な解釈には、向き依存の投影係数の包含が必要であると結論付けている。