Los rayos cósmicos en la región de transición de energías ($10^{15}$–$10^{18}$ eV) siguen siendo un enigma: sus fuentes y composición no se conocen con precisión. Las mediciones directas son imposibles debido al bajo flujo, por lo que se registran cascadas atmosféricas extensas, es decir, cascadas de partículas en la atmósfera. La principal incertidumbre está relacionada con la elección del modelo hadrónico de interacción en las simulaciones. Eventos anómalos raros, como las cascadas dobles, pueden proporcionar información independiente, pero hasta ahora no se habían podido aislar. Este trabajo se sitúa en la periferia: los estudios más cercanos se ocupan de la reconstrucción de energía y de las señales de radio de cascadas ordinarias, así como de ráfagas de radio rápidas, pero nadie había intentado detectar directamente las dobles jorobas. SKA-Low, con su densa red de antenas y su amplio ancho de banda, hace esto posible por primera vez.
Las cascadas se simularon con CORSIKA y la emisión de radio con CoREAS, para protones primarios y helio con energía de $10^{15}$ eV y ángulo cenital de 60°. Para el análisis se utilizó el paquete SMIET, que divide la atmósfera en capas de 5 g/cm² y permite calcular la contribución de cada capa a la señal de radio por separado. Primero se construyó un modelo de dos fuentes puntuales en las posiciones de los máximos de las subcascadas: este relaciona la frecuencia de la muesca de interferencia con la diferencia de tiempos de llegada de las señales. Luego, para la reconstrucción completa, se aplicó un método de fuerza bruta: se probaron los parámetros (profundidades y números de partículas en los máximos), y para cada conjunto se sintetizó rápidamente una cascada doble artificial mediante SMIET, comparando su fluencia con la de referencia mediante la desviación cuadrática media.
El modelo de dos fuentes puntuales reproduce con precisión la posición de las primeras muescas de interferencia en los espectros. Para nueve cascadas dobles simuladas, proporcionó estimaciones de las profundidades de los máximos con errores de 15 a 66 g/cm² en comparación con el ajuste de doble función de Gaisser–Hillas. La reconstrucción completa por fuerza bruta con cuatro parámetros (Xmax1, Xmax2, Nmax1, Nmax2) mostró que el mínimo global de RMS a menudo está desplazado de los valores verdaderos, pero al fijar dos parámetros la precisión mejora: las desviaciones en profundidades van de 6 a 164 g/cm², y las relaciones de números de partículas están en el rango de 0.61 a 1.21. Por lo tanto, incluso un modelo simple permite estimar aproximadamente los parámetros, y el enfoque SMIET proporciona una base para un refinamiento posterior.
Este trabajo demuestra por primera vez que la firma de radio característica de las cascadas dobles puede utilizarse para reconstruir su perfil longitudinal. A diferencia de los estudios vecinos, que se centran en cascadas ordinarias u otras fuentes, aquí se propone un camino concreto para aislar eventos anómalos raros. Esto abre la posibilidad de probar modelos hadrónicos de interacción midiendo la distancia entre los picos, y de obtener una estimación independiente de la composición de masas, especialmente para el helio, cuya partícula líder a menudo se lleva el 25% de la energía de la partícula primaria.
El siguiente paso será añadir ruidos realistas (galáctico, instrumental) y utilizar la configuración completa de la red de antenas de SKA-Low, que es un orden de magnitud más densa que la utilizada en la simulación actual. Esto permitirá comprobar con qué robustez se distinguen las muescas en los espectros en condiciones reales. También se planea ampliar la biblioteca de plantillas SMIET para diferentes atmósferas y energías, con el fin de aplicar el método a los datos de futuras observaciones.
El método impactará en la física de los rayos cósmicos de ultra alta energía y en la radioastronomía: SKA-Low podrá cumplir una doble función, registrando tanto fuentes astrofísicas como rayos cósmicos.
Probar la robustez de la reconstrucción al añadir ruidos y la red completa de antenas; ampliar la biblioteca de plantillas SMIET para diversas condiciones.
En la física de los rayos cósmicos sigue sin resolverse el problema de la elección del modelo hadrónico de interacción a energías ultra altas. Las cascadas dobles ofrecen acceso directo a la longitud de recorrido de la partícula líder, lo que es sensible a los detalles de estos modelos. Además, proponen una nueva forma de determinar la composición de masas que no depende de la calibración con cascadas ordinarias.
La frecuencia de la muesca de interferencia en el espectro de una cascada doble está relacionada con la distancia entre las dos subcascadas: cuanto más separados estén los máximos, menor será la frecuencia de la primera muesca. Para un evento típico con energía de $10^{15}$ eV, la muesca cae en el rango de SKA-Low, alrededor de 100 MHz.
La idea de utilizar la interferencia de señales de fuentes separadas recuerda al principio de la radiointerferometría descrito en la novela de Carl Sagan 'Contacto', donde la protagonista analiza una señal que llega desde diferentes direcciones para aislar una fuente débil sobre el ruido de fondo.