COSMOS: A numerical relativity code specialized for PBH formation

COSMOS es un código de relatividad numérica en C++ independiente y paralelizado con OpenMP, diseñado para simular la formación de agujeros negros primordiales mediante la resolución de las ecuaciones de Einstein en 3+1 dimensiones utilizando coordenadas de escalado no cartesianas especializadas y refinamiento de malla fija para manejar la dinámica gravitacional no lineal.

Lo esencial

La visión general: Construyendo agujeros negros desde cero

Imagina el universo temprano no como un océano suave y tranquilo, sino como un mar agitado con olas gigantes. A veces, estas olas se vuelven tan altas y pesadas que colapsan sobre sí mismas, formando Agujeros Negros Primordiales (PBH). A diferencia de los agujeros negros de los que solemos oír hablar —que son los cadáveres de estrellas muertas—, los PBH nacieron instantáneamente de las "arrugas" en el tejido del espacio y el tiempo justo después del Big Bang.

Los científicos quieren entender cómo se forman estos agujeros negros, pero las matemáticas son increíblemente difíciles. Es como intentar predecir exactamente cómo salpicará una gota de agua específica cuando golpee un charco, pero el charco se está expandiendo y el agua está hecha de pura gravedad.

La solución: COSMOS (El laboratorio digital)

Este artículo presenta COSMOS, un programa informático (escrito en C++) diseñado específicamente para simular estos nacimientos de agujeros negros. Piensa en COSMOS como un túnel de viento digital de alta tecnología para la gravedad. Así como los ingenieros construyen modelos a escala de coches para probar cómo fluye el aire a su alrededor, los físicos utilizan COSMOS para construir "modelos a escala" del universo temprano para ver cómo se comporta la gravedad cuando las cosas se ponen complicadas.

Cómo funciona: El truco del "Zoom Inteligente"

Uno de los mayores desafíos en esta simulación es que ocurren dos tamaños muy diferentes al mismo tiempo:

1. El punto diminuto: La región específica donde el agujero negro está colapsando (muy pequeño).
2. El panorama general: Todo el universo expandiéndose a su alrededor (muy grande).

Si intentas mirar todo el universo con un microscopio, pierdes el panorama general. Si miras todo el universo con un lente gran angular, no puedes ver los detalles diminutos del colapso.

COSMOS resuelve esto con una función de "Zoom Inteligente".
Imagina que estás viendo una película de una estrella colapsando. La mayor parte de la pantalla muestra el vasto y vacío universo. Pero a medida que la estrella comienza a encogerse, la cámara hace un zoom automático supercercano en ese punto específico, añadiendo más "píxeles" (resolución) exactamente donde se necesitan. Esto permite que la computadora gestione el colapso diminuto y violento sin necesidad de una supercomputadora del tamaño de una ciudad solo para calcular el espacio vacío a su alrededor.

Los ingredientes: ¿Qué hay en la simulación?

Para ejecutar la simulación, COSMOS mezcla dos ingredientes principales, como un panadero haciendo un tipo específico de masa:

1. El Fluido: Un fluido perfecto y suave que representa la materia en el universo temprano. Sigue reglas simples (como un globo inflándose o desinflándose).
2. El Campo Escalar: Un campo de energía fantasmagórico que ondula a través del espacio.

El programa resuelve las Ecuaciones de Einstein (el libro de reglas de cómo funciona la gravedad) para ver cómo interactúan estos ingredientes. Se pregunta: "Si empiezo con un bulto específico en el universo, ¿se suavizará o se aplastará sobre sí mismo para formar un agujero negro?"

La ventaja de "Sin Preparación"

Normalmente, configurar una simulación de física es como hornear un pastel donde tienes que premezclar la masa (resolver ecuaciones matemáticas complecas) antes de siquiera encender el horno.

COSMOS es diferente. Debido a que simula un universo que ya está lleno de fluido (en lugar de espacio vacío), puede empezar con la masa ya mezclada. Las condiciones iniciales se configuran de tal manera que la computadora no necesita perder tiempo resolviendo ecuaciones "elípticas" difíciles para comenzar. Simplemente presiona "ejecutar" y observa cómo se desarrolla la historia. Esto hace que el código sea más ligero, rápido y fácil de instalar y usar para otros científicos.

Lo que muestra el artículo (Los ejemplos)

El artículo incluye tres "pruebas de manejo" para mostrar cómo funciona el código:

1. La onda simple: Simula una pequeña y suave ondulación en el espacio para demostrar que el código coincide con la matemática conocida y sencilla.
2. La esfera perfecta (Adiabática): Simula una burbuja de espacio perfectamente redonda colapsando. El artículo muestra una imagen de la "función de lapso" (una medida de cómo fluye el tiempo) y cómo la computadora hace zoom en el centro a medida que se forma el agujero negro.
3. La onda fantasma (Iso-curvatura): Simula un colapso causado por ese campo de energía "fantasmagórico" mencionado anteriormente.

En todos estos casos, el código encuentra con éxito el momento en que nace un Agujero Negro (técnicamente llamado "horizonte aparente") y mapea su forma.

Por qué esto es importante

Los autores no solo están creando un juguete; están construyendo una herramienta especializada para un trabajo específico. Aunque existen otras herramientas para colisiones generales de agujeros negros (como las de estrellas en fusión), COSMOS es la llave inglesa especializada diseñada específicamente para el entorno único, caótico y en expansión del universo temprano.

Al hacer este código abierto y fácil de usar, los autores esperan que otros científicos puedan extenderlo para plantear nuevas preguntas sobre la historia oculta del universo, explicando potencialmente de qué está hecha la materia oscura o por qué vemos ondas gravitacionales hoy en día.

---

En pocas palabras: Este artículo presenta COSMOS, un programa informático especializado que actúa como una cámara de "zoom inteligente" para el universo temprano. Permite a los científicos simular cómo las diminutas arrugas en el espacio-tiempo colapsan para formar los primeros agujeros negros, haciéndolo de manera eficiente al concentrar su potencia de cálculo exactamente donde ocurre la acción.

Resumen técnico

Resumen Técnico de COSMOS: Un Código de Relatividad Numérica para la Formación de PBH

Planteamiento del Problema
Se hipotetiza que los agujeros negros primordiales (PBH) se forman en el universo temprano a partir de fluctuaciones primordiales de gran amplitud no lineal en escalas sobre-horizonte, evitando la evolución estelar. Simular su formación requiere resolver las ecuaciones de Einstein bajo condiciones donde coexisten múltiples escalas de longitud: el tamaño microscópico de la región en colapso y la escala macroscópica de la expansión cosmológica. Los enfoques numéricos estándar suelen tener dificultades para resolver estas escalas dispares de manera eficiente debido al costo computacional. Además, la obtención de datos iniciales para espaciotiempos llenos de fluido y campos escalares, en lugar de regiones de vacío asintóticamente planas, presenta desafíos específicos respecto a las ecuaciones de restricción.

Metodología
El artículo presenta COSMOS, un paquete en C++ diseñado para resolver las ecuaciones de Einstein en 3+1 dimensiones específicamente para la formación de PBH. El código es una traducción y evolución en C++ del código SACRA, desarrollado originalmente para otros propósitos.

• Modelo Físico: COSMOS modela la materia como un fluido perfecto con una ecuación de estado lineal ($P = w\rho$) y un campo escalar sin masa. El tensor de energía-impulso es la suma de las contribuciones del fluido ($T^{FL}_{\mu\nu}$) y del campo escalar ($T^{SC}_{\mu\nu}$). El sistema resuelve las ecuaciones de Einstein ($G_{\mu\nu} = 8\pi G/c^4 T_{\mu\nu}$), la ecuación del campo escalar ($\nabla_\mu\nabla^\mu\phi = 0$) y las ecuaciones de conservación del fluido ($\nabla_\mu T^{FL}_{\mu\nu} = 0$). La evolución del fluido sigue el esquema establecido por Kurganov & Tadmor (2000) y Shibata & Font (2005).
• Técnicas Numéricas: Para abordar el problema de las multiescalas, COSMOS implementa:
  • Coordenadas de escala no cartesianas: Para resolver eficientamente la región en colapso.
  • Refinamiento de Malla Fija (FMR): Para ajustar dinámicamente la resolución donde sea necesario.
  • Paralelización: Se utiliza OpenMP para lograr velocidades computacionales prácticamente aceptables.
  • Dependencias: El código está diseñado sin dependencias externas para facilitar una instalación y extensión más sencillas.
• Construcción de Datos Iniciales: A diferencia de los métodos estándar para espaciotiempos de vacío asintóticamente planos que requieren resolver ecuaciones diferenciales elípticas, COSMOS genera datos iniciales adoptando soluciones de modo creciente de longitud de onda larga (hasta el siguiente orden de magnitud en el parámetro de expansión $\epsilon = k/(aH) \ll 1$).
  • Para fluctuaciones adiabáticas, los datos iniciales se caracterizan por la perturbación de curvatura $\zeta(\vec{x})$.
  • Para perturbaciones de iso-curvatura, se caracterizan por $\Upsilon(\vec{x})$ asociado al campo escalar sin masa.
  • La distribución del fluido se genera para satisfacer las ecuaciones de restricción de Einstein con precisión de máquina, eliminando la necesidad de resolvedores elípticos.

Contribuciones Clave

1. Lanzamiento de Software Especializado: El artículo presenta el lanzamiento público de COSMOS, una herramienta diseñada específicamente para simular la formación de PBH, distinguiéndola de los códigos de relatividad numérica de propósito general.
2. Generación Eficiente de Datos Iniciales: El código implementa un método para satisfacer las ecuaciones de restricción analíticamente mediante aproximaciones de longitud de onda larga, evitando la sobrecarga computacional de los resolvedores elípticos para datos iniciales en fondos de fluidos cosmológicos.
3. Documentación y Ejemplos: El paquete incluye tres ejemplos de demostración destinados con fines instructivos:
   • Evolución de una perturbación sinusoidal de modo único (validación contra la teoría de perturbación lineal).
   • Fluctuaciones iniciales adiabáticas con simetría esférica (demostración de la formación de un horizonte aparente).
   • Perturbaciones de iso-curvatura con simetría esférica impulsadas por un campo escalar sin masa.
   • Estos ejemplos utilizan una resolución baja para demostrar la funcionalidad en lugar de resultados de alta precisión.

Resultados
El artículo proporciona resultados visuales y cualitativos de los ejemplos incluidos:

• Validación Lineal: Se muestra que la evolución temporal de la traza de la curvatura extrínseca ($trK$) para una pequeña fluctuación sinusoidal coincide con la solución de la ecuación de perturbación lineal.
• Formación de Horizonte: Las simulaciones de perturbaciones adiabáticas e iso-curvatura demuestran con éxito la formación de horizontes aparentes. El código genera archivos de datos (por ejemplo, ini_all.dat) que permiten a los usuarios reconstruir la geometría y la distribución de la materia en el momento en que se detecta un horizonte aparente.
• Refinamiento de Malla: Las visualizaciones confirman la operación de las capas de refinamiento de malla fija, mostrando regiones distintas de resolución (capas 0, 1 y 2) alrededor de la región en colapso.

Significancia y Reivindicaciones
Los autores posicionan a COSMOS como una herramienta especializada para apoyar el creciente interés en los PBH, particularmente en el contexto de las observaciones de ondas gravitacionales. El artículo afirma que el código ofrece un enfoque simplificado para simular la dinámica gravitacional no lineal en el universo temprano al:

• Proporcionar un entorno C++ autónomo que es fácil de instalar y extender para diversos entornos científicos una vez comprendido el código fuente.
• Ofrecer un método robusto para manejar el problema específico de los valores iniciales de los fluidos cosmológicos sin resolvedores elípticos.
• Servir como base para investigaciones previas y futuras, como lo evidencia su uso en numerosas publicaciones de los autores y colaboradores (por ejemplo, Yoo et al., 2013; Okawa et al., 2014; Escrivà & Yoo, 2025).

El artículo señala modestamente que los ejemplos incluidos son para demostración e instrucción, con la resolución mantenida intencionalmente baja. También aclara que algunas funciones avanzadas utilizadas en las publicaciones previas de los autores pueden no estar presentes en la versión pública, lo que implica que reproducir todos los resultados pasados podría requerir implementaciones adicionales no públicas.