Estimating amplitude of matter density fluctuations in solar and supernova models using neutrino flavor evolution

Resumen Técnico: Estimación de la Amplitud de las Fluctuaciones de la Densidad de la Materia en Modelos Solares y de Supernovas mediante la Evolución del Sabor de los Neutrinos

Planteamiento del Problema
Los neutrinos producidos en fuentes astrofísicas, específicamente el Sol y las supernovas de colapso de núcleo (CCSN), experimentan una evolución de sabor influenciada por los perfiles de densidad de materia que atraviesan. Mientras que los modelos estándar suelen asumir distribuciones de materia suaves, los entornos astrofísicos reales probablemente contienen irregularidades debidas a la convección, la turbulencia y los campos magnéticos. La literatura previa indica que tales fluctuaciones en la densidad de la materia pueden alterar significativamente la evolución del sabor de los neutrinos, induciendo potencialmente la decoherencia, suprimiendo la conversión de sabor o creando transiciones no adiabáticas. Sin embargo, determinar hasta qué punto estas irregularidades pueden ser identificadas puramente a partir de las observaciones del sabor de los neutrinos sigue siendo un desafío abierto, particularmente dada la escasez de datos de supernovas (limitados a SN1987A) y la compleja dinámica no lineal de los entornos de CCSN.

Metodología
Los autores emplean la Asimilación Estadística de Datos (SDA), una técnica de inferencia desarrollada originalmente para el pronóstico meteorológico, para determinar si las señales de neutrinos pueden revelar la amplitud de las fluctuaciones de la densidad de la materia. El enfoque trata la evolución del sabor de los neutrinos como un sistema dinámico gobernado por ecuaciones diferenciales, donde el objetivo es encontrar el conjunto de variables de estado y parámetros que minimicen una función de costo (acción) que equilibre la dinámica del modelo con las restricciones observacionales.

• Modelos Físicos: El estudio utiliza una aproximación de dos sabores para ambos escenarios, solar y CCSN.
  • Modelo Solar: Utiliza cuatro modos de energía (7, 9, 11, 13 MeV) y parámetros correspondientes a la mezcla solar ($\theta_{12}, \delta m^2_{12}$). Se desprecian las interacciones neutrino-neutrino. La densidad electrónica base se deriva del modelo solar estándar BS05(OP).
  • Modelo CCSN: Utiliza un sistema interactivo mínimo de dos modos (14 y 11 MeV) con parámetros correspondientes a la mezcla atmosférica ($\theta_{13}, \delta m^2_{13}$). Este modelo incluye interacciones no lineales $\nu$–$\nu$ mediante un potencial de "bulbo de neutrinos", basado en una simulación 1D de un progenitor de 9.6 $M_\odot$.
• Implementación de las Fluctuaciones: Las fluctuaciones de la densidad de la materia se introducen como un ruido multiplicativo $\Sigma(r)$ aplicado al potencial de materia $V(r)$. El ruido se genera a partir de una distribución gaussiana con una desviación estándar $\sigma$, filtrada para eliminar fluctuaciones por debajo de una escala de longitud mínima (100 pasos de rejilla) para asegurar la resolución, permaneciendo más rápido que la escala de evolución MSW.
• Proceso de Inferencia: El procedimiento SDA minimiza una función de acción $A_0$, compuesta por un término de modelo ($A_{model}$) que penaliza la desviación de las ecuaciones de movimiento y un término de medición ($A_{meas}$) que penaliza la desviación de los datos simulados.
  • Mediciones Simuladas: Los datos "reales" se generan integrando hacia adelante las ecuaciones de evolución con amplitudes de ruido específicas. Las mediciones consisten en la probabilidad de supervivencia del neutrino electrónico en la Tierra ($P_{\nu e, \oplus}$) y la condición inicial en la fuente ($P_z=1$).
  • Optimización: Se utiliza un procedimiento de recocido simulado (simulated annealing) para navegar por el espacio de estados, aumentando gradualmente el peso del término del modelo para guiar las trayectorias desde el ajuste de la medición hacia la consistencia dinámica completa.
  • Escenarios: El estudio prueba dos configuraciones: (1) la SDA recibe la forma exacta del perfil de ruido utilizada en la simulación directa, y (2) la SDA recibe una forma de perfil de ruido diferente y aleatorizada.

Resultados Clave
La eficacia del método SDA fue evaluada por su capacidad para inferir la amplitud de ruido verdadera ($\tilde{\sigma}$) a partir de las mediciones simuladas.

• Neutrinos Solares: El método SDA demostró una alta fiabilidad.
  • Cuando se le proporcionó la misma forma del perfil de ruido, las 10 trayectorias iniciales aleatorizadas convergieron a la amplitud de ruido verdadera para niveles de ruido bajos (0.283%). Para amplitudes más altas (1.414% y 4.243%), un subconjunto de las trayectorias convergió al valor verdadero, mientras que otras quedaron atrapadas en mínimos locales con niveles de acción significativamente más altos.
  • Crucialmente, el nivel de acción sirvió como un diagnóstico robusto: las trayectorias que convergieron a la solución verdadera exhibieron consistentemente los niveles de acción asintótica más bajos.
  • Incluso cuando se le proporcionó una forma de perfil de ruido diferente y aleatorizada, el método identificó con éxito la amplitud correcta para niveles de ruido más altos, siendo los caminos de menor acción los que convergieron al valor verdadero.
• Neutrinos de CCSN: El método mostró una sensibilidad limitada a las fluctuaciones de baja amplitud, pero se mantuvo efectivo para amplitudes altas.
  • Para amplitudes de ruido bajas (0.283% y 1.414%), el procedimiento SDA falló consistentemente en inferir la amplitud no nula verdadera, convergiendo a menudo a valores cercanos a cero.
  • Sin embargo, para la amplitud de ruido más alta probada (4.243%), el método identificó con éxito el valor correcto en una fracción de las trayectorias, distinguidas nuevamente por los niveles de acción más bajos.
  • De manera similar al caso solar, el nivel de acción filtró eficazmente las trayectorias que estimaban amplitudes de ruido significativamente mayores que la verdad.
  • El estudio señala que la sensibilidad del modelo CCSN probablemente está restringida por las longitudes de oscilación más cortas y la complejidad de las interacciones no lineales en comparación con el caso solar.

Significancia y Reivindicaciones
El artículo presenta una exploración de prueba de concepto demostrando que la Asimilación Estadística de Datos puede utilizarse para inferir información sobre las fluctuaciones de la densidad de la materia en entornos astrofísicos basándose en las observaciones del sabor de los neutrinos.

• Límites Superiores: Aunque el método no siempre pudo determinar la amplitud exacta de las fluctuaciones (particularmente en el caso de CCSN con amplitudes bajas), demostró ser capaz de excluir consistentemente las trayectorias que predecían amplitudes significativamente mayores que la verdad. Esto sugiere que el método puede establecer eficazmente límites superiores en la amplitud de tales fluctuaciones.
• Utilidad Diagnóstica: El nivel de acción se identifica como una métrica crítica para distinguir las soluciones verdaderas de los mínimos locales, una característica esencial para aplicaciones en el mundo real donde la verdad de base es desconocida.
• Potencial Futuro: Los autores postulan que este enfoque podría servir como un paso inicial para analizar modelos más realistas y de gran escala de los entornos de CCSN. Sugieren que futuros eventos de supernovas galácticas, que proporcionarían órdenes de magnitud más datos que la SN1987A, podrían permitir un sondeo más profundo de los entornos dinámicos y las irregularidades de la materia dentro de estas explosiones cataclísmicas.

El estudio concluye que, si bien los modelos simplificados actuales muestran limitaciones en los regímenes de baja amplitud, el marco de la SDA ofrece una vía viable para extraer restricciones físicas sobre los perfiles de materia astrofísica a partir de futuros datos de neutrinos.