Natural Saturation Of The Sterile Neutrino Dark Matter Resonant Production By a High Lepton Flavor Asymmetry In Primordial Plasma

Resumen Técnico: Saturación Natural de la Producción Resonante de Materia Oscura de Neutrinos Estériles

Planteamiento del Problema
Los neutrinos estériles en el rango de masa de los keV son candidatos viables para la materia oscura, producibles en el Universo temprano mediante oscilaciones activo-estériles. Si bien la producción térmica no resonante (mecanismo de Dodelson-Widrow) está ampliamente excluida por las observaciones de rayos X, la producción resonante (mecanismo de Shi-Fuller), impulsada por una asimetría de leptones primordial, sigue siendo una vía viable. Los tratamientos estándar sugieren que aumentar la asimetría de leptones inicial ($L$) potencia el potencial de materia, promoviendo así la conversión resonante y aumentando la abundancia final de neutrinos estériles.

Sin embargo, esta relación lineal se rompe ante grandes asimetrías debido a dos efectos físicos omitidos:

1. Resonancia Retrasada: Una mayor asimetría inicial desplaza la época de producción resonante hacia tiempos cosmológicos posteriores (temperaturas más bajas).
2. Depleción de la Oscilación Activa: A estas temperaturas más bajas ($T \lesssim 20$ MeV), las oscilaciones de neutrinos activos se vuelven eficientes, redistribuyendo las asimetrías de sabor entre los tres sabores activos. Esto puede reducir sustancialmente o casi borrar la asimetría de sabor específica necesaria para impulsar la resonancia antes o durante la época de producción.

Además, en el caso de asimetrías muy grandes, el potencial químico de la carga eléctrica puede aproximarse a la masa del pion cargado, lo que podría inducir la condensación de piones, alterando la condición de neutralidad de carga y el trasfondo termodinámico. Estudios previos (p. ej., Ref. [18]) estuvieron limitados a asimetrías $L < 0.1$ y no consideraron plenamente la depleción dinámica de las asimetrías ni la condensación de piones.

Metodología
Los autores desarrollaron un marco numérico acoplado para evolucionar de forma autoconsistente el trasfondo cosmológico y la producción de neutrinos estériles:

• Trasfondo Termodinámico: El sistema resuelve para los potenciales químicos ($\mu_{\nu_e}, \mu_{\nu_\mu}, \mu_{\nu_\tau}, \mu_{\text{Q}}$) bajo las restricciones de neutralidad de carga y asimetría bariónica. El sector hadrónico se modela utilizando nucleones y mesones ligeros.
• Condensación de Piones: Para asimetrías grandes donde $|\mu_{\text{Q}}| \approx m_\pi$, los autores incorporan una fase de condensación de piones utilizando un modelo de quarks-mesones (siguiendo la Ref. [19]). Emplean una prescripción mesónica simplificada donde se desprecia la contribución fermiónica al potencial termodinámico, lo que permite una solución analítica para la amplitud del condensado. Esto modifica la densidad de carga y los potenciales químicos que entran en el potencial de materia.
• Oscilaciones de Neutrinos Activos: Los autores utilizan una versión modificada del código COFLASY para resolver las ecuaciones cinéticas cuánticas de las matrices de densidad de neutrinos activos. Este código está acoplado con un resolvedor de trasfondo en Python.
  • Esquema de Acoplamiento: La evolución progresa en una rejilla de temperatura decreciente. En cada paso $T_i$, las asimetrías de sabor se pasan al resolvedor de trasfondo para actualizar los potenciales químicos y las cantidades termodinámicas. Estos se pasan luego a COFLASY para evolucionar el sector de neutrinos activos hacia $T_{i+1}$. El cambio resultante en las asimetrías de sabor ($\delta L_\alpha$) actualiza el estado para el siguiente paso.
  • Inicialización: Se prueban dos regímenes: "continuación" (preservando el vector de estado dinámico completo) y "paso-fresco" (re-inicializando la matriz de densidad basándose en las asimetrías actualizadas). El enfoque de paso-fresco resulta ser más estable para asimetrías grandes donde surgen inestabilidades numéricas en el régimen de continuación.
• Producción de Neutrinos Estériles: La abundancia de neutrinos estériles se calcula mediante una prescripción de probabilidad efectiva generalizada (Ref. [18]) que reproduce las ecuaciones cinéticas cuánticas en el régimen de grandes asimetrías. La producción se trata analíticamente mediante la aproximación de ancho estrecho de la resonancia, integrando la ecuación de Boltzmann sobre los cruces de la resonancia.

Contribuciones Clave

1. Identificación de un Mecanismo de Saturación: El artículo demuestra que aumentar la asimetría de leptones primordial no conduce a un aumento indefinido en la abundancia de neutrinos estériles. Más allá de un cierto umbral, la resonancia se desplaza a temperaturas donde las oscilaciones activas deprecian la asimetría de sabor impulsora, causando que la tasa de producción caiga o se sature.
2. Inclusión de la Condensación de Piones: El estudio incorpora los efectos de la condensación de piones en asimetrías grandes, mostrando cómo esto modifica la condición de neutralidad de carga y el potencial de materia, lo cual es crucial para el seguimiento preciso de la resonancia para $L \gtrsim 0.4$.
3. Marco de Evolución Acoplado: Los autores presentan una implementación numérica autoconsistente que evoluciona simultáneamente el trasfondo cosmológico (incluyendo condensados de piones) y las oscilaciones de sabor de los neutrinos activos, un paso necesario para modelar la producción resonante a altas asimetrías.
4. Restricciones sobre los Ángulos de Mezcla: Al mapear la abundancia máxima alcanzable frente a la asimetría inicial, los autores derivan restricciones actualizadas sobre el ángulo de mezcla activo-estéril mínimo ($\sin^2 2\theta$) requerido para explicar la densidad de materia oscura observada.

Resultos

• Saturación de la Abundancia: Para configuraciones de sabor representativas (escenarios de Jerarquía Invertida y Jerarquía Normal), la abundancia de neutrinos estériles $\Omega_{\nu_s}$ crece inicialmente con la asimetría inicial $L$. Sin embargo, a medida que $L$ aumenta (típicamente $L \gtrsim 0.5-1.0$), la resonancia se desplaza hacia temperaturas más bajas donde las oscilaciones activas reducen rápidamente la asimetría de sabor relevante. Esto provoca una caída brusca en la eficiencia de producción. Incluso con asimetrías mayores, la resonancia puede entrar en un plateau de asimetría reducida o en la época de las oscilaciones del sabor electrón (donde todas las asimetrías son llevadas a cero), impidiendo un mayor crecimiento en la abundancia.
• Ángulo de Mezcla Mínimo: El efecto de saturación implica un límite inferior para el ángulo de mezcla requerido. Incluso con la asimetría máxima, el ángulo de mezcla no puede ser arbitrariamente pequeño. Los autores identifican un espacio de parámetros de $\sin^2 2\theta$ que abarca aproximadamente dos órdenes de magnitud (p. ej., $10^{-18}$ a $10^{-16}$ para $m_s \sim 40$ keV) que es consistente con el mecanismo resonante y los límites observacionales actuales.
• Espectros de Momento: Las distribuciones de momento resultantes son no térmicas. A medida que la asimetría aumenta, el momento promedio aumenta (materia oscura "caliente"), aproximándose a $\langle x \rangle \approx 3.2$. En los casos donde la producción se solapa con las oscilaciones activas, emergen características espectrales complejas (p. ej., picos dobles). No obstante, los autores señalan que las incertidumbres numéricas en estos regímenes de solapamiento impiden una determinación totalmente robusta de la forma espectral detallada, particularmente para el Caso I (mezcla de sabor tau).
• Estabilidad Numérica: El método de inicialización de "paso-fresco" resulta más robusto que el de "continuación" para asimetrías grandes, evitando las inestabilidades numéricas observadas en el segundo cuando las oscilaciones activas son rápidas.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona un "mecanismo de saturación natural" para la producción resonante de neutrinos estériles, actuando como una restricción adicional para los modelos que dependen de grandes asimetrías de leptones. La interacción entre la resonancia retrasada y la depleción por oscilación activa limita la densidad máxima de materia oscura alcanzable mediante este mecanismo.

El artículo concluye que:

• Las sensibilidades actuales de los telescopios de rayos X son insuficientes para explorar plenamente este mecanismo; se requiere una mejora de un orden de magnitud en la sensibilidad al ángulo de mezcla para detectar los picos espectrales predichos.
• La determinación precisa de la abundancia máxima y los espectros de momento está actualmente limitada por la complejidad numérica de la evolución acoplada. El trabajo futuro requiere un marco numérico unificado que integre la evolución del trasfondo y las ecuaciones cinéticas cuánticas para resolver las épocas de solapamiento de la oscilación activa y la producción estéril.
• Las grandes asimetrías primordiales pueden inducir una transición de fase de QCD de primer orden mediante la condensación de piones, lo que potencialmente genera un fondo estocástico de ondas gravitacionales. Detectar tal señal podría proporcionar evidencia independiente de las grandes asimetrías requeridas para la materia oscura de neutrinos estériles por resonancia.

Los autores mantienen la modestia respecto a los espectros de momento detallados en los regímenes más complejos (Caso I), reconociendo que las irregularidades de pequeña escala en sus resultados numéricos actuales no deben asignárseles una interpretación física sin realizar comprobaciones de convergencia adicionales.