Spectrum of radiation from global strings and the relic axion density

Resumen Técnico: Espectro de Radiación de Cuerdas Globales y la Densidad de Axiones Reliquia

Planteamiento del Problema
La producción de axiones mediante defectos topológicos (cuerdas cósmicas globales) es un candidato principal para explicar la densidad observada de Materia Oscura Fría (CDM); sin embargo, la masa del axión predicha ($m_a$) y la constante de decaimiento ($f_a$) siguen siendo altamente inciertas debido a escenarios teóricos contradictorios respecto al espectro de radiación emitido por estas cuerdas.

• Escenario A: Basado en la acción de Kalb-Ramond (KR), que trata las cuerdas como objetos tipo Nambu acoplados a un campo tensor antisimétrico, este escenario predice un espectro "suave" donde la radiación tiene su pico cerca del modo fundamental de las perturbaciones de la cuerda ($q > 1$).
• Escenario B: Respaldado por diversas simulaciones numéricas de teoría de campos, este escenario sugiere un espectro "duro" con una divergencia logarítmica ($q \approx 1$), lo que implica una emisión significativa en altas frecuencias hasta el ancho del núcleo de la cuerda.

La discrepancia entre estos escenarios conduce a predicciones muy distintas para la densidad de axiones reliquia. Además, las simulaciones numéricas existentes suelen carecer del rango dinámico ($\log(\Delta/\delta) \approx 70$) necesario para resolver completamente el núcleo de la cuerda frente a la escala de perturbación, y pueden estar midiendo inadvertidamente el campo propio de la cuerda en lugar de la radiación de axiones emitida.

Metodología
Los autores emplean un enfoque de dos vertientes que combina el modelado analítico y la simulación numérica:

1. Marco Analítico:

   • Los autores derivan expresiones para la densidad de axiones reliquia ($\Omega_a$) considerando dos canales de producción: axiones emitidos directamente por la red de cuerdas largas y axiones emitidos por el decaimiento de bucles de cuerda.
   • Introducen una función espectral $g(z)$ para modelar el espectro de radiación, parametrizada por un índice de ley de potencia ($q$ o $p$) y un corte ($n_*$ o $m_*$).
   • Definen un factor espectral $G_2(q, n_*)$ que encapsula la dependencia de la densidad de reliquia sobre el espectro de radiación. Este factor se calcula tanto para modelos de ley de potencia como para modelos de caída exponencial.
   • El modelo incorpora la "eficiencia de corte" ($c$) y la fracción de energía ($F_\ell$) emitida en bucles frente a la emitida directamente en axiones, reconociendo la incertidumbre en estos parámetros entre las simulaciones de cuerdas de Nambu y las de teoría de campos.

2. Simulaciones Numéricas de Cuerdas Rectas Perturbadas:

   • Los autores realizan simulaciones numéricas de alta resolución del modelo de Goldstone $U(1)$ utilizando una rejilla de hasta $801^3$ puntos para resolver el núcleo de la cuerda y la radiación.
   • Descomposición de Campo: Una innovación metodológica crítica es la separación explícita del campo de fase total ($\alpha$) en el campo propio de la cuerda ($\alpha_{str}$) y el campo de axión radiado ($\Delta\alpha$). Los autores demuestran que las salidas estándar de las simulaciones suelen confundir estos dos, los cuales están dominados por el campo propio cerca del núcleo de la cuerda.
   • Técnicas de Sustracción: Comparan dos métodos para aislar la radiación de axiones:
     1. Sustracción directa de un ansatz de cuerda estática ($\alpha \approx \theta + \Delta\alpha$).
     2. Cálculo de la proyección del flujo de momento ($P_\alpha$) del tensor de energía-momento.
   • Enmascaramiento: Para asegurar que el espectro medido represente axiones libres y no el campo propio de la cuerda, los autores aplican una máscara cilíndrica (excluyendo una región $r < r_{crit} \approx 10\Delta x$) alrededor de la cuerda antes de realizar el análisis espectral.

Resultos Clave

1. Caracterización Espectral:

   • Las simulaciones revelan que, sin el enmascaramiento, el espectro medido está dominado por el campo propio de la cuerda, el cual exhibe una caída de ley de potencia de $k^{-2.25}$.
   • Tras aplicar la sustracción del campo propio y el enmascaramiento, se encuentra que el espectro de la radiación de axiones emitida ($\Delta\alpha$) es compatible con una caída exponencial, consistente con las predicciones de la acción de Kalb-Ramond tipo Nambu.
   • Los resultados contradicen explícitamente el "espectro duro" ($q \approx 1$) reportado frecuentemente en simulaciones de redes sin enmascaramiento. Los autores argumentan que los espectros "duros" previos fueron artefactos de la medición del campo propio en lugar de los axiones radiados.

2. Dependencia de los Parámetros Espectrales:

   • Se muestra que la densidad de reliquia es altamente sensible al índice espectral $q$ y al corte $n_*$, particularmente en el rango $1 < q < 2$.
   • Para $q \le 1$, la densidad está dominada por modos de alta frecuencia (ancho del núcleo), mientras que para $q > 1$, está dominada por modos de baja frecuencia (longitud de onda de perturbación).

3. Revisiones de las Predicciones de Masa del Axión:

   • Los autores cuantifican las incertidumbres en el cálculo de la densidad de reliquia. Encuentran que el rango de masa del axión predicho es más amplio que el derivado del Mecanismo de Desalineación Inicial (IMM).
   • Escenario 1 (Espectro Suave): Si el modo de decaimiento primario es a través de un espectro suave de bucles producidos cerca de la escala de retroacción, los autores estiman $m_a \approx 160 \, \mu\text{eV}$ (frecuencia $f \approx 38 \, \text{GHz}$).
   • Escenario 2 (Espectro Duro de Simulaciones): Si los axiones son emitidos directamente por la red usando espectros de simulaciones típicas de teoría de campos sin enmascaramiento, $m_a \approx 4 \, \mu\text{eV}$ ($f \approx 1 \, \text{GHz}$).
   • Escenario 3 (Cuerda Oscilante con Espectro Corregido): Usando los espectros corregidos de los autores para una cuerda oscilante (exponencial/suave), la estimación se desplaza a $m_a \approx 125 \, \mu\text{eV}$ ($f \approx 30 \, \text{GHz}$).

Significado y Reivindicaciones
El artículo sostiene que la incertidumbre en la densidad de reliquia del axión es impulsada en gran medida por la incapacidad de las simulaciones numéricas actuales para distinguir entre el campo propio de la cuerda y la radiación de axiones emitida.

• Corrección de la Metodología: Los autores afirman que se requiere una "corrección significativa" (excluir una región considerable alrededor de la cuerda) para analizar los axiones correctamente. Argumentan que las simulaciones estándar de teoría de campos a menudo fallan en esto, llevando a la conclusión errónea de un espectro duro.
• Resolución de la Discrepancia: Al demostrar que el espectro de radiación es exponencial (suave) en lugar de duro, el artículo sugiere que la tensión entre las predicciones de la acción de Kalb-Ramond y las simulaciones de teoría de campos puede resolverse mediante mejores técnicas de análisis de datos en lugar de un fallo fundamental de la aproximación tipo Nambu.
• Impacto en la Detección: Los hallazgos implican que, si el espectro suave es el correcto, la masa del axión es probablemente mayor ($\sim 100\text{--}160 \, \mu\text{eV}$) que los valores derivados de simulaciones no corregidas, desplazando la frecuencia objetivo para los experimentos de detección de axiones al rango de $\sim 30\text{--}40 \, \text{GHz}$.

Los autores concluyen con modestia, reconociendo que persisten incertidumbres significativas respecto a la fracción de energía emitida en bucles ($F_\ell$) y los parámetros exactos del espectro de radiación, pero enfatizan que la forma espectral es el factor dominante para determinar la ventana viable de la masa del axión.