COSMOS: A numerical relativity code specialized for PBH formation

COSMOS est un code de relativité numérique en C++ autonome, parallélisé avec OpenMP, conçu pour simuler la formation de trous noirs primordiaux en résolvant les équations d'Einstein en 3+1 dimensions en utilisant des coordonnées de mise à l'échelle non cartésiennes spécialisées et un raffinement de maillage fixe pour gérer la dynamique gravitationnelle non linéaire.

L'essentiel

La vue d'ensemble : Construire des trous noirs à partir de rien

Imaginez l'univers primitif non pas comme un océan calme et lisse, mais comme une mer agitée avec des vagues géantes. Parfois, ces vagues deviennent si hautes et si lourdes qu'elles s'effondrent sur elles-mêmes, formant des trous noirs primordiaux (PBH). Contrairement aux trous noirs dont on entend habituellement parler — qui sont les cadavres d'étoiles mortes — les PBH sont nés instantanément des « rides » dans le tissu de l'espace et du temps juste après le Big Bang.

Les scientifiques veulent comprendre comment ces trous noirs se forment, mais les mathématiques sont incroyablement difficiles. C'est comme essayer de prédire exactement comment une goutte d'eau spécifique éclaboussera en frappant une flaque, mais la flaque est en train de s'étendre, et l'eau est faite de pure gravité.

La solution : COSMOS (Le laboratoire numérique)

Ce document présente COSMOS, un programme informatique (écrit en C++) conçu spécifiquement pour simuler ces naissances de trous noirs. Considérez COSMOS comme une soufflerie numérique de haute technologie pour la gravité. Tout comme les ingénieurs construisent des modèles réduits de voitures pour tester comment l'air circule autour d'elles, les physiciens utilisent COSMOS pour construire des « modèles réduits » de l'univers primitif afin de voir comment la gravité se comporte quand les choses deviennent désordonnées.

Comment ça marche : L'astuce du « Zoom intelligent »

L'un des plus grands défis de cette simulation est qu'il se passe deux choses de tailles très différentes en même temps :

1. Le point minuscule : La région spécifique où le trou noir s'effondre (très petit).
2. La vue d'ensemble : L'univers entier qui s'étend tout autour (très grand).

Si vous essayez de regarder l'univers entier avec un microscope, vous perdez la vue d'ensemble. Si vous regardez l'univers entier avec un objectif grand angle, vous ne pouvez pas voir les détails minuscules de l'effondrement.

COSMOS résout cela avec une fonction de « Zoom intelligent ».
Imaginez que vous regardez un film d'une étoile en train de s'effondrer. La majeure partie de l'écran montre l'immensité de l'univers vide. Mais alors que l'étoile commence à rétrécir, la caméra zoome automatiquement de très près sur ce point précis, ajoutant plus de « pixels » (résolution) exactement là où ils sont nécessaires. Cela permet à l'ordinateur de gérer l'effondrement minuscule et violent sans avoir besoin d'un supercalculateur de la taille d'une ville juste pour calculer l'espace vide autour.

Les ingrédients : Que contient la simulation ?

Pour faire tourner la simulation, COSMOS mélange deux ingrédients principaux, comme un boulanger préparant un type de pâte spécifique :

1. Le Fluide : Un fluide parfait et lisse qui représente la matière dans l'univers primitif. Il suit des règles simples (comme un ballon qui se gonfle ou se dégonfle).
2. Le Champ Scalaire : Un champ d'énergie fantomatique qui ondule à travers l'espace.

Le programme résout les équations d'Einstein (le livre de règles de la gravité) pour voir comment ces ingrédients interagissent. Il demande : « Si je commence avec une bosse spécifique dans l'univers, va-t-elle s'atténuer ou va-t-elle s'écraser sur elle-même pour former un trou noir ? »

L'avantage du « Sans préparation »

Habituellement, préparer une simulation physique est comme cuire un gâteau où vous devez pré-mélanger la pâte (résoudre des équations mathématiques complexes) avant même d'allumer le four.

COSMOS est différent. Parce qu'il simule un univers qui est déjà rempli de fluide (plutôt que d'un espace vide), il peut commencer avec la pâte déjà mélangée. Les conditions initiales sont si parfaites que l'ordinateur n'a pas besoin de passer du temps à résoudre des équations « elliptiques » difficiles pour démarrer. Il lui suffit d'appuyer sur « exécution » et de regarder l'histoire se dérouler. Cela rend le code plus léger, plus rapide et plus facile à installer et à utiliser pour d'autres scientifiques.

Ce que le papier démontre (Les exemples)

Le document inclut trois « essais routiers » pour montrer comment le code fonctionne :

1. L'onde simple : Il simule une petite ride douce dans l'espace pour prouver que le code correspond à des mathématiques simples connues.
2. La sphère parfaite (adiabatique) : Il simule une bulle d'espace parfaitement ronde en train de s'effondrer. Le papier montre une image de la « fonction de lapse » (une mesure de la façon dont le temps s'écoule) et comment l'ordinateur zoome sur le centre au moment où le trou noir se forme.
3. L'onde fantôme (iso-courbure) : Il simule un effondrement causé par ce champ d'énergie « fantomatique » mentionné plus haut.

Dans tous ces cas, le code trouve avec succès le moment où un Trou Noir est né (techniquement appelé « horizon apparent ») et cartographie sa forme.

Pourquoi cela importe

Les auteurs ne fabriquent pas seulement un jouet ; ils construisent un outil spécialisé pour un travail précis. Bien que d'autres outils existent pour les collisions de trous noirs généraux (comme ceux provenant de la fusion d'étoiles), COSMOS est la clé à molette spécialisée conçue spécifiquement pour l'environnement unique, complexe et en expansion de l'univers primitif.

En rendant ce code ouvert et facile d'utilisation, les auteurs espèrent que d'autres scientifiques pourront l'étendre pour poser de nouvelles questions sur l'histoire cachée de l'univers, expliquant potentiellement de quoi la matière noire est faite ou pourquoi nous voyons des ondes gravitationnelles aujourd'hui.

---

En résumé : Ce document présente COSMOS, un programme informatique spécialisé qui agit comme une caméra à « zoom intelligent » pour l'univers primitif. Il permet aux scientifiques de simuler comment de minuscules rides dans l'espace-temps s'effondrent pour former les tout premiers trous noirs, en le faisant efficacement en concentrant leur puissance de calcul exactement là où l'action se déroule.

Résumé technique

Résumé technique de COSMOS : Un code de relativité numérique pour la formation des PBH

Énoncé du problème
Les trous noirs primordiaux (PBH) sont supposés se former dans l'univers primordial à partir de fluctuations primordiales de grande amplitude non linéaire sur l'horizon, contournant ainsi l'évolution stellaire. La simulation de leur formation nécessite la résolution des équations d'Einstein dans des conditions où plusieurs échelles de longueur coexistent : la taille microscopique de la région en effondrement et l'échelle macroscopique de l'expansion cosmologique. Les approches numériques standards peinent souvent à résoudre efficacement ces échelles disparates en raison du coût computationnel. De plus, l'obtention de données initiales pour des espaces-temps remplis de fluides et de champs scalaires, plutôt que des régions de vide asymptotiquement plates, présente des défis spécifiques concernant les équations de contrainte.

Méthodologie
Le document présente COSMOS, un package C++ conçu pour résoudre les équations d'Einstein en 3+1 dimensions spécifiquement pour la formation des PBH. Le code est une traduction et une évolution en C++ du code SACRA, développé à l'origine pour d'autres objectifs.

• Modèle physique : COSMOS modélise la matière comme un fluide parfait avec une équation d'état linéaire ($P = w\rho$) et un champ scalaire sans masse. Le tenseur énergie-impulsion est la somme des contributions du fluide ($T^{FL}_{\mu\nu}$) et du champ scalaire ($T^{SC}_{\mu\nu}$). Le système résout les équations d'Einstein ($G_{\mu\nu} = 8\pi G/c^4 T_{\mu\nu}$), l'équation du champ scalaire ($\nabla_\mu\nabla^\mu\phi = 0$), et les équations de conservation du fluide ($\nabla_\mu T^{FL}_{\mu\nu} = 0$). L'évolution du fluide suit le schéma établi par Kurganov & Tador (2000) et Shibata & Font (2005).
• Techniques numériques : Pour traiter le problème multi-échelle, COSMOS implémente :
  • Des coordonnées de mise à l'échelle non cartésiennes : Pour résoudre efficacement la région en effondrement.
  • Le raffinement de maillage fixe (FMR) : Pour ajuster dynamiquement la résolution là où cela est nécessaire.
  • La parallélisation : OpenMP est utilisé pour atteindre des vitesses de calcul acceptables en pratique.
  • Dépendances : Le code est conçu sans dépendances externes pour faciliter l'installation et l'extension.
• Construction des données initiales : Contrairement aux méthodes standards pour les espaces-temps de vide asymptotiquement plats qui nécessitent la résolution d'équations différentielles elliptiques, COSMOS génère des données initiales en adoptant des solutions de modes croissants à grande longueur d'onde (jusqu'à l'ordre suivant de l'étalon de l'expansion $\epsilon = k/(aH) \ll 1$).
  • Pour les fluctuations adiabatiques, les données initiales sont caractérisées par la perturbation de courbure $\zeta(\vec{x})$.
  • Pour les perturbations d'iso-courbure, elles sont caractérisées par $\Upsilon(\vec{x})$ associé au champ scalaire sans masse.
  • La distribution du fluide est générée pour satisfaire les équations de contrainte d'Einstein à la précision de la machine près, éliminant ainsi le besoin de solveurs elliptiques.

Contributions clés

1. Publication d'un logiciel spécialisé : Le document présente la version publique de COSMOS, un outil spécifiquement adapté à la simulation de la formation des PBH, se distinguant des codes de relativité numérique à usage général.
2. Génération efficace des données initiales : Le code implémente une méthode pour satisfaire les équations de contrainte de manière analytique via des approximations de grande longueur d'onde, évitant la surcharge computationnelle des solveurs elliptiques pour les données initiales dans les fonds de fluides cosmologiques.
3. Documentation et exemples : Le package inclut trois exemples de démonstration destinés à des fins pédagogiques :
   • Évolution d'une perturbation sinusoïdale à mode unique (validation par rapport à la théorie des perturbations linéaires).
   • Fluctuations initiales adiabatiques à symétrie sphérique (démontrant la formation d'un horizon apparent).
   • Perturbations d'iso-courbure à symétrie sphérique pilotées par un champ scalaire sans masse.
   • Ces exemples utilisent une faible résolution pour démontrer la fonctionnalité plutôt que des résultats de haute précision.

Résultats
Le document fournit des résultats visuels et qualitatifs issus des exemples inclus :

• Validation linéaire : L'évolution temporelle de la trace de la courbure extrinsèque ($trK$) pour une petite fluctuation sinusoïdale montre qu'elle correspond à la solution de l'équation de perturbation linéaire.
• Formation d'horizon : Les simulations de perturbations adiabatiques et d'iso-courbure démontrent avec succès la formation d'horizons apparents. Le code génère des fichiers de données (par exemple, ini_all.dat) permettant aux utilisateurs de reconstruire la géométrie et la distribution de la matière au moment où un horizon apparent est détecté.
• Raffinement de maillage : Les visualisations confirment le fonctionnement des couches de raffinement de maillage fixe, montrant des régions de résolution distinctes (couches la plus basse, 1ère et 2ème) autour de la région en effondrement.

Signification et revendications
Les auteurs positionnent COSMOS comme un outil spécialisé pour soutenir l'intérêt croissant pour les PBH, particulièrement dans le contexte des observations d'ondes gravitationnelles. Le document affirme que le code offre une approche rationalisée pour simuler la dynamique gravitationnelle non linéaire dans l'univers primordial en :

• Fournissant un environnement C++ autonome, facile à installer et à étendre pour divers contextes scientifiques une fois le code source compris.
• Offrant une méthode robuste pour gérer le problème spécifique des valeurs initiales des fluides cosmologiques sans solveurs elliptiques.
• Servant de fondation pour des recherches antérieures et futures, comme en témoigne son utilisation dans de nombreuses publications des auteurs et de leurs collaborateurs (ex. Yoo et al., 2013 ; Okawa et al., 2014 ; Escrivà & Yoo, 2025).

Le document note modestement que les exemples inclus sont à des fins de démonstration et d'instruction, la résolution étant intentionnellement maintenue basse. Il précise également que certaines fonctions avancées utilisées dans les publications antérieures des auteurs peuvent ne pas être présentes dans la version publique, impliquant que la reproduction de tous les résultats passés peut nécessiter des implémentations supplémentaires non publiques.