Heavy neutral leptons beyond the BBN bound: probing the lepton asymmetry of the Universe

Resumen Técnico: Leptones neutros pesados más allá del límite de la BBN: sondeando la asimetría de leptones del Universo

Planteamiento del Problema
La cosmología estándar impone una restricción estricta sobre los Leptones Neutros Pesados (HNL, por sus siglas en inglés) con canales de decaimiento hadrónicos y tiempos de vida $\tau_N \gtrsim 0.02$ s. En el Universo temprano, los HNL con masas $m_N > m_l + m_\pi$ decaen en piones cargados ($N \to \pi^+ l^-$ y $\bar{N} \to \pi^- l^+$). En una población simétrica ($n_N = n_{\bar{N}}$), los mesones inyectados impulsan interacciones fuertes ($\pi^- + p \to n + \pi^0/\gamma$) que convierten protones en neutrones mucho más rápido que las tasas débiles. Esto impulsa la relación neutrón-protón ($n/p$) hacia un valor isospin-simétrico de $\sim 1$, provocando una sobreproducción catastrófica de helio ($Y_p \approx 0.6$) y deuterio, violando las observaciones de la Nucleosíntesis del Big Bang (BBN). En consecuencia, este "límite de la BBN" excluye una parte significativa del espacio de parámetros al que apuntarán los próximos experimentos de aceleradores como SHiP, DUNE e Hyper-K, particularmente para $m_N \lesssim 1$ GeV y $\tau_N \lesssim 1$ s. Los intentos previos para evadir este límite se basaron en suprimir toda la abundancia de los HNL mediante un recalentamiento de baja escala, lo cual no deja firmas cosmológicas correlacionadas.

Metodología
Los autores proponen un mecanismo donde una gran asimetría de sabor leptónico primordial ($L_\alpha$) permite que los HNL evadan el límite de la BBN sin suprimir su población. El estudio emplea un flujo numérico acoplado de tres etapas para rastrear la evolución del Universo desde la producción de los HNL hasta la BBN:

1. Etapa de Producción ($T \sim \text{GeV} \to 20$ MeV): Una evolución de Boltzmann resuelta en el momento calcula la producción de HNL de Dirac ($N$) y sus antipartículas ($\bar{N}$) mediante la mezcla con neutrinos activos. Crucialmente, el cálculo tiene en cuenta la redistribución de las cargas conservadas (bariónica, eléctrica y de sabor leptónico) entre todas las especies del plasma. La presencia de grandes asimetrías de sabor crea una diferencia de potencial químico, lo que resulta en una relación de población desigual $r \equiv n_N / n_{\bar{N}} \neq 1$.
2. Etapa de Oscilación y Dinámica de Sabor ($20$ MeV $\to 5$ MeV): Utilizando un marco de cinética cuántica (COFLASY), los autores evolucionan el sistema de sabor de los neutrinos en presencia de los HNL en decaimiento. Esta etapa rastrea cómo los decaimientos asimétricos inyectan asimetrías de sabor leptónico en el sector activo y cómo las oscilaciones de neutrinos redistribuyen estas asimetrías.
3. Etapa de Transporte y BBN ($T < 5$ MeV): Se utiliza un enfoque híbrido. Un código de transporte de Boltzmann resuelto en el momento ($\nu$DSMC) sigue la evolución espectral de neutrinos, antineutrinos y productos de decaimiento (mesones, muones) para calcular con precisión las tasas débiles de $p \leftrightarrow n$ y el número efectivo de especies de neutrinos ($N_{\text{eff}}$). Esto se acopla con un código de red nuclear para calcular las abundancias primordiales ($Y_p$, D/H). El código rastrea explícitamente la inyección de neutrinos no térmicos y las tasas de conversión impulsadas por mesones.

Contribuciones Clave y Mecanismo
La contribución central es la demostración de que una asimetría de carga en la población de HNL puede cancelar los efectos destructivos de la inyección de mesones en la BBN:

• Desequilibrio de Carga de Mesones: En un escenario de Dirac con $r \neq 1$, los decaimientos inyectan números desiguales de $\pi^+$ y $\pi^-$. Dado que el $\pi^-$ impulsa la conversión $p \to n$ (sobreproduciendo neutrones) y el $\pi^+$ impulsa la conversión $n \to p$, un exceso de $\pi^+$ puede revertir el enriquecimiento de neutrones.
• Restauración de las Abundancias Estándar: Al ajustar la relación de asimetría $r$ a un valor específico ($r_{\text{req}} \approx 4\text{--}10$), el efecto neto de la inyección de mesones puede ajustarse para dejar la relación $n/p$ en su valor estándar de la BBN cuando los mesones desaparecen.
• Cancelación de $N_{\text{eff}}$: Las mismas asimetrías de sabor leptónico que crean el desequilibrio de los HNL también aumentan la densidad de energía de los neutrinos (elevando $N_{\text{eff}}$). Por el contrario, el decaimiento de los HNL en neutrinos no térmicos y el subsiguiente calentamiento del plasma electromagnético tienden a reducir $N_{\text{eff}}$. Los autores muestran que estos dos efectos pueden cancelarse, manteniendo el cambio $\Delta N_{\text{eff}}$ dentro de los límites observacionales.
• Ventana Cinemática: El mecanismo es efectivo para $m_l + m_\pi \lesssim m_N \lesssim 1.5$ GeV. Por encima de $\sim 1.5$ GeV, los estados finales de múltiples hadrones y los kaones neutros ($K_L$) diluyen la correlación de carga, impidiendo el desequilibrio necesario.

Resultos
El estudio identifica cuatro regiones distintas en el espacio de parámetros de los HNL ($m_N, \tau_N$) que son consistentes con las restricciones de la BBN y el Fondo Cósmico de Microondas (CMB):

1. Región 1 (Invisible): Parámetros donde $Y_p$, D/H y $N_{\text{eff}}$ permanecen indistinguibles de la cosmología estándar.
2. Región 2 (Déficit de $N_{\text{eff}}$): Un cambio negativo detectable en $N_{\text{eff}}$ (por ejemplo, $\Delta N_{\text{eff}} \approx -0.08$), que potencialmente explica las pistas observacionales actuales de $N_{\text{eff}} < 3.044$.
3. Región 3 (Fondo de Neutrinos Relictos Modificado): Un escenario donde el mar de neutrinos relictos se ve incrementado en más de un $30\%$ en densidad numérica debido a la asimetría superviviente y las colas no térmicas, incluso si $N_{\text{eff}}$ parece estándar. Esto ofrece un objetivo para experimentos como PTOLEMY.
4. Región 4 (Ondas Gravitacionales): Un subconjunto del espacio de parámetros donde las grandes asimetrías de sabor leptónico impulsan la transición de fase de QCD a ser de primer orden, pudiendo originar un fondo estocástico de ondas gravitacionales.

El artículo proporciona puntos de referencia (Tabla I) que demuestran que estos escenarios son compatibles con los límites actuales de SHiP, DUNE e Hyper-K, reabriendo efectivamente el dominio previamente excluido para las búsquedas de HNL.

Significado y Afirmaciones
El artículo afirma identificar el primer mecanismo por el cual una población producida por mezcla de partículas con decaimiento hadrónico puede evadir las restricciones de la BBN sin suprimir la abundancia de la partícula.

• Implicaciones del Descubrimiento: Un descubrimiento de un HNL en el dominio previamente excluido no sería meramente la detección de una nueva partícula; serviría como una sonda de las asimetrías de sabor leptónico primordiales del Universo.
• Firmas Correlacionadas: La detección de un HNL constreñiría las asimetrías primordiales a una familia de un solo parámetro. La distribución específica de la asimetría entre los sabores electrón y tau determinaría si la transición de QCD fue de primer orden (señal de GW) y la magnitud del incremento del fondo de neutrinos relictos (señal de PTOLEMY).
• Consistencia Teórica: El escenario requiere que el HNL sea de Dirac (para mantener poblaciones distintas de $N$ y $\bar{N}$) y que el número leptónico se conserve a un alto grado ($\delta m \lesssim 10^{-15}$ eV). Esto vincula las búsquedas en aceleradores con preguntas fundamentales sobre la naturaleza de la masa de los neutrinos y la historia térmica del Universo.

Los autores enfatizan que su cálculo es el primero en rastrear simultáneamente las grandes asimetrías de sabor leptónico a través de la dinámica de sabor, la evolución de neutrinos resuelta en el momento y la nucleosíntesis impulsada por mesones, estableciendo la existencia de estas soluciones y sus improntas cosmológicas correlacionadas.