Une vision unifiée de la matière noire et de la cosmologie : les modèles de diffusion dans le lentillage gravitationnel ✛ ✦

Original: "Joint inference of weak lensing convergence map and cosmology with diffusion models"
arXiv:2606.31988 · 2026-06-30 · CC BY 4.0 · 4 min
La nouvelle méthode JADE combine la reconstruction des cartes de lentillage gravitationnel et l'inférence bayésienne des paramètres cosmologiques sans simulateur différentiable.

Contexte

Les relevés célestes modernes, comme LSST, collectent d'énormes quantités de données sur les formes de milliards de galaxies. L'effet de lentillage gravitationnel faible, prédit par Einstein dans le cadre des équations d'Einstein, se manifeste par des distorsions cohérentes de ces formes sous l'effet de la gravité de la matière noire. L'extraction de l'information cosmologique de ces données repose traditionnellement sur le calcul des fonctions de corrélation, mais cela perd une part significative de l'information codée dans les caractéristiques non gaussiennes de la distribution de matière. Le passage à une analyse au niveau des champs promet une amélioration radicale de la précision, mais nécessite de résoudre simultanément deux problèmes : la reconstruction des cartes de masse cachées et l'estimation des paramètres cosmologiques, ce qui est réalisé dans ce travail.

Méthodes

Les auteurs ont proposé la méthode JADE, basée sur des modèles de diffusion — des réseaux de neurones génératifs qui apprennent à transformer du bruit en échantillons de données réalistes. Contrairement aux approches traditionnelles, le modèle est entraîné directement sur des simulations de champs de convergence log-normaux, reproduisant l'agglutination non linéaire de la matière et dépendant de six paramètres cosmologiques du modèle wCDM : la densité de matière noire, celle de énergie noire, l'amplitude des fluctuations et d'autres. L'architecture utilise un transformeur qui traite des sections de carte aplaties en jetons, ainsi que les paramètres cosmologiques, comme une séquence unique, ce qui permet de tenir compte de leur influence mutuelle. À l'étape d'inférence, le modèle entraîné, recevant le champ de cisaillement observé, génère des échantillons a posteriori conjoints des cartes de convergence et des paramètres, en résolvant une équation différentielle stochastique.

Résultats

Le modèle JADE, testé sur des données simulées du relevé LSST Y10 (champ de 5×5 degrés, 128×128 pixels, 5 bins tomographiques en décalage vers le rouge), a montré un excellent accord avec l'échantillonnage classique par MCMC. Les distributions a posteriori cosmologiques pour les paramètres Ωc, σ8, w0 et d'autres ont coïncidé avec les chaînes NUTS, et la métrique de calibration MIRA a atteint 0.635±0.017 (idéal 2/3). Les cartes de convergence reconstruites reproduisent non seulement les corrélations à deux points (spectre de puissance avec une erreur relative de 6.7%), mais aussi correctement la statistique non gaussienne, y compris les distributions à un point avec une erreur relative de 5.2%. Il est important de noter que toutes les vérifications ont été effectuées image par image : pour chaque échantillon, une simulation indépendante a été générée avec les paramètres correspondants, confirmant la cohérence interne de la distribution a posteriori conjointe.

Implications

Ce travail démontre que l'apprentissage profond peut unifier l'inférence de champs et de paramètres dans un schéma bayésien unique, éliminant la nécessité d'avoir une réponse directe différentiable — longtemps un goulot d'étranglement pour de telles analyses. Cela signifie que l'on peut désormais utiliser des simulateurs aussi complexes que l'on veut, y compris des codes complets à N corps et hydrodynamiques, sans perdre la possibilité d'une inférence bayésienne honnête. La méthode jette les bases d'une analyse des données préservant toute l'information non gaussienne, ce qui est particulièrement précieux aux petites échelles, où l'évolution non linéaire des structures génère une riche structure de corrélation d'ordre supérieur.

Développements futurs

Le développement ultérieur de la méthode prévoit le passage à des simulations plus réalistes à N corps et l'extension aux sphères célestes courbes au lieu de l'approximation plane. Un axe important sera l'inclusion dans le processus d'apprentissage des gradients de modélisation directe, si disponibles, pour améliorer l'efficacité d'utilisation des simulations coûteuses. De plus, l'architecture pourrait être adaptée pour reconstruire conjointement les conditions initiales de l'Univers, ce qui permettrait de tester les modèles cosmologiques en prédisant des champs non observables.

Impact

L'approche proposée aura un impact sur l'astronomie extragalactique et la cosmologie, en fournissant un outil pour extraire le maximum d'information des relevés à venir. Elle pourrait également être utilisée dans d'autres problèmes nécessitant une reconstruction conjointe de champs et de paramètres globaux, par exemple en astrophysique des hautes énergies ou en physique de l'expansion de l'univers.

Prochaines étapes

Les prochaines étapes incluent l'application de JADE à des simulations plus réalistes, la prise en compte du masquage et des effets systématiques, ainsi que l'exploration de la possibilité d'une analyse multi-sondes combinant le lentillage avec d'autres observables cosmologiques.

Problèmes ouverts clés

La méthode est directement liée au problème de la matière noire et de l'énergie noire, permettant de mesurer plus précisément les paramètres qui gouvernent l'évolution de l'Univers. Elle nous rapproche également de l'utilisation complète de l'information contenue dans la structure à grande échelle, ce qui est l'un des défis clés de la cosmologie moderne.

Le regard de la machine du savoir

Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.

Ce que fait ce travail

Les auteurs présentent JADE, un modèle de diffusion pour l'inférence conjointe des cartes de convergence et des paramètres cosmologiques à partir des données de lentillage faible. Le modèle est entraîné sur des simulations de champs log-normaux et montre un accord avec MCMC sur les distributions a posteriori des paramètres, tout en reconstruisant la statistique non gaussienne des champs.

Où est sa force

JADE combine l'inférence de champ et l'inférence paramétrique dans un schéma bayésien unifié, sans nécessiter de réponse directe différentiable, ce qui lève une limitation clé des approches précédentes.

Vers où avancer

À essayer

  • Tester la robustesse de JADE face à des bruits et des masques réalistes, en utilisant des approches issues de l'article 2607.09867.
  • Explorer la possibilité d'utiliser JADE pour générer rapidement des cartes de convergence synthétiques avec des paramètres cosmologiques donnés.
  • Évaluer comment l'inclusion d'informations provenant des ondes gravitationnelles pourrait améliorer les contraintes sur les paramètres de l'énergie noire.

Pourquoi cela peut compter

La combinaison de l'inférence de champ et de l'inférence paramétrique ouvre la voie à une utilisation plus complète de l'information dans les données de lentillage faible, ce qui est crucial pour les futurs relevés comme LSST, et pourrait conduire à des contraintes plus précises sur la nature de la matière noire et de l'énergie noire.

Résumé original

tel que rédigé par les auteurs, en anglais arXiv:2606.31988

We present a method for joint inference of cosmological parameters and convergence maps from weak lensing observations, targeting the full posterior conditioned on the observed shear field. Our approach uses implicit inference with diffusion models, learning the joint distribution from simulations, without the need to have an explicit and differentiable forward model for gradient-based MCMC sampling. We introduce a transformer-based architecture that operates in pixel space and treats cosmological parameters as additional tokens in a unified sequence, enabling efficient multimodal processing within a single network. At inference time, the trained model generates posterior samples of joint convergence maps and cosmological parameters conditioned on observed noisy shear fields. We demonstrate the method on simulated weak lensing data generated from log-normal fields in a wcdm cosmology. The model accurately reconstructs convergence maps and recovers cosmological posteriors that agree with traditional MCMC, while remaining well calibrated across the prior, with a MIRA calibration score of $$0.635 \pm 0.017$$ on the joint posterior (where $$0.667$$ is optimal). The inferred fields reproduce the correct two-point statistics as well as non-Gaussian statistics such as the one-point distribution. This work establishes diffusion-based implicit inference as a viable route toward full field-level cosmological analyses, paving the way for applications to more realistic, non-differentiable simulators.

Pour entraîner le modèle, il n'a fallu que 100 000 simulations, alors qu'une seule exécution de MCMC pour une observation nécessite environ un million d'appels au simulateur. De plus, la génération d'un échantillon a posteriori ne prend que 0.2 seconde.

En science-fiction, l'idée de cartographier des structures invisibles par leur influence gravitationnelle a été exploitée dans le roman « Blindness » de Peter Watts, où des vaisseaux extraterrestres se camouflent mais se trahissent par la distorsion de l'arrière-plan.

\kappa(\boldsymbol{\vartheta}) = \frac{3 H_0^2 \Omega_m}{2c^2} \int_0^{\chi_{\lim}} d\chi \, \frac{q(\chi)}{a(\chi)} f_K(\chi) \delta(f_K(\chi) \boldsymbol{\vartheta}, \chi)
Équation (1) de l'article : la convergence comme intégrale pondérée du contraste de densité le long de la ligne de visée.
\tilde{\kappa}_E + i\tilde{\kappa}_B = \left( \frac{k_1^2 - k_2^2}{k^2} + i \frac{2k_1 k_2}{k^2} \right) (\tilde{\gamma}_1 + i\tilde{\gamma}_2)
Équation (3) : relation entre le cisaillement et la convergence dans l'espace de Fourier en approximation plane.

Chiffres clés

  • Nombre de simulations pour l'entraînement: 100 000
  • Taille du champ: 5×5 degrés
  • Temps de génération d'un échantillon: 0.2 seconde
  • Métrique MIRA: 0.635±0.017
  • Erreur relative du spectre de puissance: 6.7%
Scientifiques
Fritz ZwickyVera RubinRoberto D. PecceiHelen QuinnR. Brent TullyBohdan Paczyński
Notions
cisaillement gravitationnel faible distribution a posteriori apprentissage profond en imagerie médicale apprentissage faiblement supervisé modélisation de lentille gravitationnelle cartographie de l'image médicale au génome moindres carrés imagerie à pixel unique matrice de confusion modèles de diffusion spectromètre intégré régression symbolique reconstruction de microlentille par réseau de neurones convergence composite analysis échantillonnage imbriqué irm de diffusion rééchantillonnage jackknife échantillonnage par rejet autoregressive model