How X-rays heat the IGM in different 21-cm simulation codes: a comparison between Licorice and Beorn

Imagina el universo como un gigantesco e invisible océano de gas, compuesto principalmente de hidrógeno, que se extiende en todas las direcciones. Durante los primeros cientos de millones de años después del Big Bang, este océano fue oscuro, frío y silencioso. Luego, las primeras estrellas y galaxias se encendieron como faros en la niebla. Su intensa radiación comenzó a calentar el gas y a arrancar electronos de los átomos, un proceso llamado "reionización". Los astrónomos están desesperados por ver este período desarrollarse porque contiene los secretos de cómo creció nuestro universo. Para lograrlo, están construyendo un enorme radiotelescopio llamado el Square Kilometre Array (SKA), que actúa como un oído gigante escuchando un tenue susurro: la "señal de 21 cm". Esta señal es una onda de radio específica emitida por el hidrógeno neutro, y su intensidad cambia dependiendo de qué tan caliente esté el gas. Si el gas está frío, la señal se ve diferente de si está caliente. Al medir esta señal, los científicos esperan mapear la historia del calentamiento del universo.

Sin embargo, hay un inconveniente. El universo es demasiado vasto y complejo para observarlo directamente en tiempo real; tenemos que construir modelos computacionales para simular lo que sucedió hace miles de millones de años. Estas simulaciones son como máquinas del tiempo digitales. Pero, al igual que cualquier máquina del tiempo, solo son tan buenas como las reglas que siguen. Algunas simulaciones son increíblemente detalladas, rastreando cada partícula de gas y cada rayo de luz, pero requieren semanas de supercomputadoras para ejecutarse. Otras son más rápidas y utilizan atajos, lo que las hace más fáciles de usar pero potencialmente menos precisas. La gran pregunta es: ¿funcionan estos atajos lo suficientemente bien como para confiar en los resultados? Si dos códigos de simulación diferentes nos cuentan historias distintas sobre cuánto se calentó el gas, ¿a cuál deberíamos creerle? Este es el rompecabezas que un nuevo estudio de Romain Meriot y su equipo se propone resolver.

El equipo decidió poner a competir cara a cara a dos "máquinas del tiempo" muy diferentes para ver cómo manejan el calentamiento del gas cósmico mediante rayos X. Por un lado, tenían a Licorice, una simulación 3D de alta definición que actúa como un videojuego de movimiento completo, calculando cómo se mueve el gas y cómo viaja la luz a través de él en tres dimensiones. Por el otro lado, tenían a Beorn, un código más rápido y simple que utiliza un enfoque "unidimensional". Piensa en Licorice como un maestro chef cocinando un estofado complejo, revolviendo cada ingrediente individualmente, mientras que Beorn es como un microondas que calienta toda la comida basándose en una receta promedio. Los investigadores querían ver si el microondas podía producir una comida que supiera tan bien como el estofado del chef, específicamente con respecto a cómo los rayos X de las primeras estrellas calentaron el gas de hidrógeno.

Para que la pelea fuera justa, los científicos no se limitaron a dejar que los dos códigos funcionaran con sus propios ajustes aleatorios. En su lugar, tomaron los "ingredientes" (las estrellas y su producción de luz) de la detallada simulación Licorice y los introdujeron en Beorn. Esto aseguró que ambos códigos estuvieran intentando cocinar exactamente el mismo universo, solo que usando métodos diferentes. Luego compararon los resultados: la temperatura del gas, la intensidad de la señal de 21 cm y los patrones estadísticos de cómo fluctúa esa señal a través del cielo.

Los resultados fueron una mezcla de buenas noticias y una advertencia. Cuando observaron el panorama general —la temperatura promedio del gas y la intensidad general de la señal—, los dos códigos coincidieron muy bien. Era como si tanto el microondas como el chef hubieran logrado que el estofado alcanzara la temperatura promedio correcta. Sin embargo, cuando observaron más de cerca los detalles, las diferencias empezaron a mostrarse. La distribución del calor no era la misma. En la detallada simulación 3D, el calor era irregular, con algunos puntos mucho más calientes y otros más fríos, creando una textura compleja. En la más rápida simulación 1D, el calor era más suave y uniforme.

Estas pequeñas diferencias en la "textura" del calor resultaron ser sorprendentemente importantes. Cuando el equipo calculó el "espectro de potencia" —una forma elegante de medir cuánto fluctúa la señal en diferentes escalas—, encontraron una discrepancia de aproximadamente el 30% entre los dos códigos. Para ponerlo en perspectiva, si estuvieras tratando de adivinar la receta de un pastel basándote en un solo bocado, una diferencia del 30% en el sabor significaría que adivinarías los ingredientes de forma errónea.

Los investigadores luego profundizaron para descubrir por qué los códigos no estaban de acuerdo. Descubrieron que los atajos utilizados en el código más rápido eran los culpables. Específicamente, el código 1D asumía que la densidad del gas era la misma en todas partes (como asumir que el océano es perfectamente plano) y utilizó una fórmula simplificada para la cantidad de calor que los rayos X realmente depositaban. En la realidad, el gas es grumoso, y el calentamiento depende de exactamente dónde se encuentra el gas. Cuando el equipo obligó al código detallado 3D a usar estos mismos atajos simplificados, los resultados se acercaron más al código 1D, confirmando que estas aproximaciones eran la fuente del error.

Finalmente, el equipo se hizo la pregunta definitiva: "Si usamos estas diferentes simulaciones para analizar los datos del futuro telescopio SKA, ¿obtendremos respuestas incorrectas?". Simularon un escenario en el que el telescopio observaría durante 100 horas. Encontraron que las diferencias entre los códigos conducirían a un sesgo significativo en la estimación de las propiedades de las primeras estrellas. En términos estadísticos, el error fue lo suficientemente grande como para ser mayor a 1 sigma (una medida de confianza), e incluso superior para algunos parámetros. Esto significa que, si los astrónomos usaran el código más rápido y simple para interpretar datos reales, podrían concluir con confianza que las primeras estrellas eran diferentes de lo que realmente eran.

El estudio concluye que, si bien los códigos más rápidos son útiles y aciertan la historia del "promedio", aún no son lo suficientemente precisos para la era del SKA. La diferencia del 30% en el espectro de potencia y el sesgo resultante en los datos sugieren que no podemos simplemente confiar en los atajos cuando intentamos medir el universo con extrema precisión. Los autores sugieren que, aunque los códigos son lo suficientemente similares para los telescopios actuales, menos sensibles, se necesita mucho trabajo para refinar los códigos más rápidos o desarrollar nuevas formas de corregir sus errores antes de que el SKA comience a escuchar. Es un recordatorio de que, en la búsqueda por comprender el amanecer del universo, el diablo está realmente en los detalles.