Spectrum of radiation from global strings and the relic axion density

技術要約：グローバル・ストリングからの放射スペクトルと遺物アクシオン密度

問題提起
トポロジカル欠陥（グローバル宇宙ひも）を介したアクシオン生成は、観測された冷暗黒物質（CDM）密度を説明するための主要な候補である。しかし、これらひもから放出される放射のスペクトルに関する理論的シナリオの相違により、予測されるアクシオン質量（$m_a$）および崩壊定数（$f_a$）には依然として高い不確実性が存在する。

• シナリオA： ひもを、反対称テンソル場に結合したナブー様（Nambu-like）の物体として扱うカルブ・ラムンド（KR）作用に基づくこのシナリオは、放射がひもの摂動の基本モード（$q > 1$）付近でピークを持つ「ソフト」なスペクトルを予測する。
• シナリオB： 様々な数値場理論シミュレーションによって支持されているこのシナリオは、対数的な発散（$q \approx 1$）を伴う「ハード」なスペクトルを示唆しており、ひものコア幅に至る高周波数領域での大幅な放射を意味している。

これらのシナリオ間の相違は、遺物アクシオン密度に対して極めて異なる予測をもたらす。さらに、既存の数値シミュレーションは、ひものコアと摂動スケールを完全に解像するために必要なダイナミックレンジ（$\log(\Delta/\delta) \approx 70$）を欠いていることが多く、また、放出されたアクシオン放射ではなく、ひもの自己場を誤って測定している可能性がある。

手法
著者らは、解析的モデリングと数値シミュレーションを組み合わせた二段構えのアプローチを採用している。

1. 解析的枠組み：

   • 著者らは、2つの生成チャネル（長大なひもネットワークから直接放出されるアクシオンと、ひもループの崩壊によって放出されるアクシオン）を考慮して、遺物アクシオン密度（$\Omega_a$）の式を導出している。
   • 放射スペクトルをモデル化するために、べき乗則指数（$q$ または $p$）とカットオフ（$n_*$ または $m_*$）でパラメータ化されたスペクトル関数 $g(z)$ を導入している。
   • 放射スペクトルの依存性を集約するスペクトル因子 $G_2(q, n_*)$ を定義し、これはべき乗則モデルと指数関数的減衰モデルの両方について計算されている。
   • モデルは、「チョッピング効率」（$c$）および、ナブー・ストリングと場理論シミュレーション間の不確実性を認識した上での、ループへのエネルギー放出分（$F_\ell$）を取り入れている。

2. 摂動を受けた直線ひもの数値シミュレーション：

   • 著者らは、ひものコアと放射を解像するために、最大 $801^3$ 個の格子点を持つ $U(1)$ ゴールドストーン・モデルの高解像度数値シミュレーションを実行している。
   • 場の分解： 重要な手法上の革新は、全場の位相（$\alpha$）を、ひもの自己場（$\alpha_{str}$）と放射されたアクシオン場（$\Delta\alpha$）に明示的に分離することである。著者らは、標準的なシミュレーションの出力が、これら2つ（特にひものコア付近で支配的な自己場）を混同していることを示している。
   • 減算手法： 著者らは、放出されたアクシオン放射を孤立させるための2つの方法を比較している。
     1. 静的なひも・アンザッツ（$\alpha \approx \theta + \Delta\alpha$）による直接減算。
     2. エネルギー・運動量テンソルの運動量フラックス射影（$P_\alpha$）の計算。
   • マスキング： 測定されるスペクトルが、ひもの自己場ではなく自由なアクシオンを表すことを確実にするため、スペクトル解析を行う前に、ひもの周囲の円筒状領域（$r < r_{crit} \approx 10\Delta x$）を排除する円筒マスクを適用している。

主な結果

1. スペクトルの特性評価：

   • シミュレーションの結果、マスキングを行わない場合、測定されるスペクトルはひもの自己場によって支配されており、$k^{-2.25}$ のべき乗則減衰を示すことが明らかになった。
   • 自己場減算とマスキングを適用した後、放出されたアクシオン放射（$\Delta\alpha$）のスペクトルは、指数関数的減衰と整合しており、これはカルブ・ラムンド作用の予測と一致している。
   • この結果は、アンマスキングされたネットワーク・シミュレーションでしばしば報告される「ハード」なスペクトル（$q \approx 1$）を明確に否定している。著者らは、以前の「ハード」なスペクトルは、放射されたアクシオンではなく自己場を測定したことによるアーティファクトであると主張している。

2. スペクトル・パラメータへの依存性：

   • 遺物密度は、スペクトル指数 $q$ およびカットオフ $n_*$ に対して非常に敏感であり、特に $1 < q < 2$ の範囲において顕著である。
   • $q \le 1$ の場合、密度は高周波モード（コア幅）によって支配されるが、$q > 1$ の場合は低周波モード（摂動波長）によって支配される。

3. 修正されたアクシオン質量の予測：

   • 著者らは、遺物密度の計算における不確実性を定量化している。予測されるアクシオン質量範囲は、初期ミスマ alignment メカニズム（IMM）から導かれるものよりも広いことが判明した。
   • シナリオ1（ソフトなスペクトル）： 主な崩壊モードが、バックリアクション・スケール付近で生成されたループからのソフトなスペクトルを介する場合、著者らは $m_a \approx 160 \, \mu\text{eV}$ （周波数 $f \approx 38 \, \text{GHz}$）と推定している。
   • シナリオ2（シミュレーションによるハードなスペクトル）： 通常のアンマスキングされた場理論シミュレーションのスペクトルを用いて、ネットワークから直接アクシオンが放出される場合、$m_a \approx 4 \, \mu\text{eV}$ （$f \approx 1 \, \text{GHz}$）となる。
   • シナリオ3（修正されたスペクトルを持つ振動するひも）： 著者らの修正されたスペクトル（指数関数的／ソフト）を用いた場合、推定値は $m_a \approx 125 \, \mu\text{eV}$ （$f \approx 30 \, \text{GHz}$）へとシフトする。

意義と主張
本論文は、遺物アクシオン密度の不確実性の大部分は、現在の数値シミュレーションが、ひもの自己場と放出されたアクシオン放射を区別できないことに起因すると主張している。

• 手法の修正： 著者らは、アクシオンを正しく解析するためには、「大幅な修正」（相当の大きさの領域の排除）が必要であると断言している。彼らは、標準的な場理論シミュレーションはしばしばこれを怠り、ハードなスペクトルという誤った結論を導いていると論じている。
• 不一致の解消： 放射スペクトルが（ハードではなく）指数関数的（ソフト）であることを示すことで、本論文は、カルブ・ラムンド作用の予測と場理論シミュレーションとの間の緊張関係が、ナブー近似の根本的な失敗ではなく、より優れたデータ解析手法によって解決される可能性があることを示唆している。
• 検出への影響： もしソフトなスペクトルが正しいのであれば、アクシオン質量は、修正されていないシミュレーションから導かれる値よりも高くなる（$\sim 100\text{--}160 \, \mu\text{eV}$）可能性があり、これにより、アクシオン検出実験のターゲット周波数が $\sim 30\text{--}40 \, \text{GHz}$ の範囲へとシフトすることを意味している。

著者らは、ループへのエネルギー放出分（$F_\ell$）や放射スペクトルの正確なパラメータに関して依然として大きな不確実性が残っていることを認めつつも、スペクトルの形状が実行可能なアクシオン質量窓を決定する最も支配的な要因であることを強調し、謙虚に結論づけている。