Incorporating curved geometry in cosmological simulations

Este artículo presenta un marco relativista completo para ejecutar simulaciones cosmológicas con geometría espacial curva, resolviendo el problema de las condiciones de contorno mediante la inmersión de una tapa esférica en un exterior plano para permitir una modelización coherente de observables a grandes distancias.

Lo esencial

Imagina que el universo es una inmensa fiesta cósmica. Durante mucho tiempo, los físicos han asumido que el "suelo" de esta fiesta es perfectamente plano, como una mesa de billar infinita. Pero, ¿y si en realidad el suelo tiene una ligera curvatura, como si fuera la superficie de una pelota gigante (cerrada) o de una silla de montar (abierta)?

Este es el gran misterio que Julian Adamek y Renan Boschetti intentan resolver en su nuevo trabajo. Quieren crear simulaciones por computadora que no asuman que el universo es plano, sino que incluyan esa posible curvatura.

Aquí te explico cómo lo hacen, usando analogías sencillas:

1. El problema: La caja de juguetes cuadrada

Los científicos usan superordenadores para simular cómo se forman las galaxias. Para que los cálculos sean posibles, estos ordenadores suelen trabajar dentro de "cajas" imaginarias con paredes invisibles. La forma más fácil de manejar estas paredes es hacerlas periódicas: si un juguete sale por la derecha, entra por la izquierda, como en el videojuego Pac-Man.

El problema es que estas cajas son cuadradas y planas. Si el universo real es curvo (como una esfera), no puedes meter una esfera perfecta dentro de un cubo cuadrado sin deformarla o romper las reglas de la física. Hasta ahora, los científicos usaban un "truco": simulaban un pedacito pequeño de universo curvo, pero ignoraban la curvatura geométrica real, asumiendo que la distancia era recta. Esto funciona bien para cosas pequeñas, pero falla cuando miramos muy lejos, donde la curvatura del espacio se nota.

2. La solución: El "pastel de cumpleaños" cósmico

Los autores proponen una idea genial: incrustar una esfera curvada dentro de un espacio plano.

Imagina que tienes un pastel de cumpleaños (el espacio plano de la simulación). En lugar de ponerle una fruta encima, cortas un agujero en el centro y colocas dentro una bola de gelatina curvada (el universo curvo que queremos estudiar).

• La gelatina (el parche curvo): Representa el universo con curvatura. Aquí es donde viven los observadores y donde ocurren las cosas interesantes.
• La masa del pastel (el exterior plano): Rodea a la gelatina. Esta parte es plana y permite que el simulador use sus reglas habituales (las paredes cuadradas) sin romperse.

Gracias a las ecuaciones de Einstein, saben exactamente cómo "pegar" la gelatina al pastel sin que se vea la costura. Es como si la gravedad hiciera que la transición fuera suave y natural.

3. Los observadores: Mirando desde el centro o desde el borde

En su simulación, pueden poner a un "observador" (nosotros) en dos lugares:

• Observador A (en el centro): Mira hacia afuera en todas direcciones. Ve un universo curvo perfecto, como si estuviera dentro de una esfera.
• Observador B (cerca del borde): Está casi tocando la "costura" donde la gelatina se une al pastel. Si mira hacia el centro, ve la curvatura. Si mira hacia afuera, ve el espacio plano.

Lo increíble es que, aunque el espacio es una mezcla de curvo y plano, la física funciona perfectamente. La luz viaja, las galaxias se mueven y las distancias se calculan correctamente, respetando la geometría real.

4. ¿Por qué es importante?

Antes, si queríamos predecir qué verían los telescopios modernos (como el Euclid o el Rubin) en un universo curvo, teníamos que usar aproximaciones que fallaban a grandes distancias. Era como intentar medir la circunferencia de la Tierra usando una regla recta: a corta distancia funciona, pero si intentas medir el globo entero, te equivocas.

Con este nuevo método:

• Precisión total: Pueden simular la luz viajando a través de un universo curvo y calcular exactamente cómo se vería desde la Tierra.
• Prueba de la inflación: Si el universo es curvo, significa que la teoría de la "inflación" (que dice que el universo se expandió rapidísimo y se volvió plano) podría estar incompleta o necesitar ajustes.
• El futuro: Esto ayuda a entender si el universo es finito (como una esfera) o infinito, y cómo eso afecta la masa de los neutrinos y la energía oscura.

En resumen

Los autores han creado un "puente" matemático que les permite meter un universo curvo dentro de una caja de simulación cuadrada, sin romper las reglas de la física. Es como si pudieran simular el interior de una pelota de fútbol dentro de una caja de zapatos, y que todo funcionara como si la pelota fuera el universo real.

Esto nos permite hacer "modelos a futuro" (forward modelling) mucho más precisos, preparándonos para cuando los telescopios del futuro nos envíen datos que podrían cambiar nuestra comprensión de cómo nació y cómo es nuestro universo.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "Incorporating curved geometry in cosmological simulations" de Julian Adamek y Renan Boschetti, estructurado según los puntos solicitados.

1. El Problema

La curvatura espacial es un parámetro cosmológico fundamental que debe ser medido y no asumido como cero (plano) a priori. Aunque las observaciones actuales (como las de DESI, Euclid y el CMB) imponen restricciones estrictas, la detección de una curvatura no nula tendría implicaciones profundas para la teoría de la inflación y el origen del universo.

El desafío principal radica en la modelización hacia adelante (forward modelling) de observaciones de gran escala. Las simulaciones cosmológicas estándar (N-cuerpos) suelen utilizar condiciones de contorno periódicas en cajas cúbicas, lo que implica una geometría espacial euclidiana (plana) más allá de la escala de periodicidad. Esto presenta dos problemas:

1. Geometría no euclidiana: A grandes distancias (cuando los fotones viajan desde alto corrimiento al rojo), la relación entre distancia y volumen en un universo curvo difiere significativamente de la euclidiana. El enfoque de "universo separado" (separate universe), que ajusta la tasa de expansión pero mantiene la geometría espacial plana, es insuficiente para calcular observables en el cono de luz pasado a grandes escalas.
2. Condiciones de contorno: Es difícil implementar condiciones de contorno no periódicas en códigos de simulación existentes sin reescribir toda la infraestructura.

2. Metodología

Los autores proponen un marco relativista completo para ejecutar simulaciones cosmológicas con geometría espacial curva, superando las limitaciones de los códigos Newtonianos estándar.

• Estrategia de Empotramiento (Embedding):

  • Se embebe un "parche" esférico de un espacio-tiempo FLRW curvo (que contiene al observador y las observaciones) dentro de una solución exterior espacialmente plana.
  • Esta configuración permite utilizar condiciones de contorno periódicas en la caja de simulación exterior, manteniendo la compatibilidad con códigos existentes.
  • La unión entre el parche curvo y el exterior plano se realiza mediante una región de vacío descrita por la solución de Schwarzschild, generalizando el modelo de Einstein-Straus.

• Transformación de Coordenadas y Gauge:

  • La solución exacta para la materia sin presión (polvo) en el parche curvo se describe mediante la métrica Lemaître-Tolman-Bondi (LTB) en coordenadas sincrónicas.
  • Los códigos de simulación (como gevolution) operan en el gauge de Poisson. Los autores derivan la transformación de coordenadas entre el gauge sincrónico y el de Poisson hasta segundo orden en teoría de perturbaciones. Esto es crucial para mantener la precisión cuando la curvatura es grande y para capturar el acoplamiento entre perturbaciones y curvatura.
  • Se definen dos observadores: uno en el centro del parche (Observador A) y otro cerca del borde (Observador B), permitiendo estudiar la isotropía y los efectos de perspectiva.

• Condiciones Iniciales y Defecto de Masa Relativista:

  • Para establecer las condiciones iniciales, se debe tener en cuenta el defecto de masa relativista. En un universo cerrado, el volumen propio es mayor que en uno plano para las mismas dimensiones de contorno. Si se usa una distribución de partículas homogénea estándar, la densidad efectiva dentro del parche curvo sería incorrecta.
  • Los autores corrigen esto ajustando la masa de las partículas (o distribuyendo la masa faltante) para satisfacer la restricción hamiltoniana, asegurando que el campo de desplazamiento sea compatible con las condiciones de contorno periódicas.

• Perturbaciones de Materia:

  • Se utiliza un método aproximado para inicializar perturbaciones, asumiendo que los modos en el parche curvo pueden aproximarse por modos de Fourier del exterior plano (similar al enfoque de universo separado), pero con una conversión cuidadosa de las unidades de los números de onda y los parámetros cosmológicos.

3. Contribuciones Clave

• Marco Relativista para Geometría Curva: Es la primera vez que se presenta un marco para simular la formación de estructuras en un universo con curvatura global, considerando explícitamente los efectos de distancia y volumen no euclidianos en el cono de luz.
• Solución Exacta de Empalme: Proporcionan una solución exacta para la métrica de Schwarzschild en el gauge de Poisson dentro de un espacio-tiempo FLRW en expansión, derivando expresiones analíticas exactas para los potenciales de Bardeen ($\phi$ y $\psi$) en la región de vacío.
• Correcciones de Segundo Orden: La derivación de la transformación de gauge hasta segundo orden permite una precisión superior a la de las aproximaciones lineales, esencial para curvaturas significativas.
• Compatibilidad con Códigos Existentes: Demuestran que su metodología puede integrarse en códigos N-cuerpos relativistas (como gevolution) y, con aproximaciones, en códigos Newtonianos, requiriendo cambios principalmente en las condiciones iniciales y el post-procesamiento.

4. Resultados Numéricos

Los autores validaron su marco mediante tres experimentos numéricos utilizando el código gevolution:

1. Simulación Einstein-de Sitter (Curvatura Alta):

   • Se ejecutó un modelo con $\tilde{\Omega}_K = -0.25$ (curvatura cerrada fuerte).
   • Resultado: Las relaciones distancia-corrimiento al rojo ($d_L-z$) inferidas por observadores centrales y periféricos coinciden con la predicción exacta del modelo curvo con una precisión de una parte por mil hasta $z \approx 1.7$.
   • Comparación: El enfoque de "universo separado" (que ignora la geometría espacial) muestra errores del ~1% a $z=1$, demostrando la necesidad de la geometría curvada.

2. Simulación $\Lambda$CDM Curvo (Sin Perturbaciones):

   • Modelo con $\tilde{\Omega}_K = -0.1$, $\tilde{\Omega}_m = 0.4$, $\tilde{\Omega}_\Lambda \approx 0.7$.
   • Resultado: Se confirma la concordancia con el modelo de entrada curvo para ambos observadores. Se observa que el observador periférico puede ver a corrimientos al rojo más altos antes de que el cono de luz cruce a la región plana. El error residual es de ~0.05%, atribuible a la precisión del parámetro de Hubble.

3. Simulación $\Lambda$CDM Curvo (Con Perturbaciones):

   • Se incluyeron perturbaciones de materia iniciales.
   • Resultado: Se midió el espectro de potencia angular ($C_\ell$) de la agrupación de materia. Los resultados coinciden con las predicciones lineales de CLASS.
   • Principio Cosmológico: Se verificó que las estadísticas de agrupación son independientes de la posición del observador (central vs. periférico) y de la dirección, confirmando que el parche curvo satisface el principio cosmológico.

5. Significado e Impacto

Este trabajo representa un avance significativo en la cosmología teórica y observacional:

• Precisión en Futuras Encuestas: Con la llegada de datos de alta precisión de misiones como Euclid, DESI y el Observatorio Vera C. Rubin, los errores sistemáticos debidos a la suposición de geometría plana podrían dominar sobre el ruido estadístico. Este método permite modelar estos efectos geométricos correctamente.
• Prueba de la Inflación: Permite realizar pruebas rigurosas de la predicción de inflación de un universo plano. Si se detecta curvatura, este marco es esencial para interpretar correctamente los datos y distinguir entre modelos de inflación y física pre-inflacionaria.
• Más Allá de la Aproximación de Universo Separado: Demuestra que la modificación de la historia de expansión (enfoque de universo separado) no es suficiente para predecir observables en el cono de luz a grandes escalas; la geometría espacial no euclidiana tiene un impacto directo y medible (del orden del 1-10% en ciertos estadísticos) que debe ser incluido en las simulaciones.
• Herramienta para la Comunidad: Al implementar esto en una versión actualizada de gevolution y mostrar que es compatible con códigos Newtonianos mediante correcciones de condiciones iniciales, hacen accesible esta capacidad a una amplia gama de investigadores.

En resumen, el artículo proporciona la primera herramienta robusta y relativista para simular la formación de estructuras en universos con curvatura espacial, cerrando la brecha entre la teoría de perturbaciones de alto orden y las simulaciones numéricas de gran escala.