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.
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.
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.
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.
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.
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.
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.
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.
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
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.
Travaux voisins dans notre archive
L'archive est dense autour de ce travail : 3 analyses voisines, la plus proche à 0.72. Le terrain est fréquenté, et il y a de quoi se comparer.
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.