Modélisation de la distribution du redshift des galaxies lentillées d'Euclid pour des analyses au niveau du champ

La mission Euclid fournira des données sur les lentilles gravitationnelles faibles avec une précision sans précédent, ce qui pourrait révolutionner notre compréhension de l'énergie noire et de la croissance des structures cosmiques. Pour en extraire toute la richesse informative, il faut aller au-delà des analyses standard. Afin d'exploiter au mieux ces données, le projet OCAPi analysera les cartes de lentille gravitationnelle d'Euclid directement au niveau des pixels. Cette approche, connue sous le nom d'inférence au niveau du champ, permet de capturer toutes les informations et d'obtenir des contraintes jusqu'à 5 fois plus précises sur les paramètres cosmologiques (Porqueres et al. 2022, 2023).

Cette précision accrue nécessite toutefois une modélisation précise des données. L'un des principaux défis de l'étalonnage dans les relevés de lentilles gravitationnelles faibles est la distribution des galaxies lentillées en fonction du redshift. Les méthodes d'étalonnage actuelles ont été conçues pour les analyses standard et peuvent ne pas être suffisamment précises pour les techniques au niveau du champ. Il est essentiel de quantifier les exigences de précision et de développer des méthodes capables de les atteindre afin de permettre l'analyse au niveau du champ des données d'Euclid et de libérer tout le potentiel scientifique.

L'objectif de ce projet de doctorat est de développer un nouvel échantillonneur de redshifts pour la lentille gravitationnelle faible, conçu pour répondre aux exigences de précision de l'inférence au niveau du champ. Cet échantillonneur combinera des modèles physiques de populations de galaxies avec des techniques d'apprentissage automatique flexibles. La thèse contribuera à maximiser le potentiel des données de lentille gravitationnelle faible d'Euclid et à faire progresser notre compréhension de la formation des structures cosmiques.

RECHERCHES D’EMISSION DIFFUSES EN RAYONS GAMMA DE TRES HAUTE ENERGIE ET PHYSIQUE FONDAMENTALE AVEC H.E.S.S. ET CTAO

Les observations en rayons gamma de très hautes énergies (THE, E>100 GeV) sont cruciales pour la compréhension des phénomènes non-thermiques les plus violents à l’œuvre dans l’Univers. La région centre de la Voie Lactée est une région complexe et active en rayons gamma de THE. Parmi les sources gamma de THE se trouvent le trou noir supermassif Sagittarius A* au coeur de la Galaxie, des vestiges de supernova ou encore des régions de formation d'étoiles. Le centre Galactique (CG) abrite un un accélérateur de rayons cosmiques jusqu’à des énergies du PeV, des émissions diffuses du GeV au TeV dont le « Galactic Center Excess » (GCE) dont l’origine est encore inconnue, de potentielles sources variables au TeV, ainsi que possibles populations de sources non encore résolues (pulsars millisecondes, trous noirs de masses intermédiaires). Le CG devrait être la source la plus brillante d’annihilations de
particules massives de matière noire de type WIMPs. Des candidats matière noire plus légers, les particules de type axions (ALP), pourraient se convertir en photons, et vice versa, dans les champs magnétiques laissant une empreinte d’oscillation dans les spectres gamma de noyaux actifs de galaxies (AGN).
L'observatoire H.E.S.S. situé en Namibie est composé de cinq télescopes imageurs à effet Cherenkov atmosphérique. Il est conçu pour détecter des rayons gamma de quelques dizaines de GeV à plusieurs dizaines de TeV. La région du Centre Galactique est observée par H.E.S.S. depuis vingt ans. Ces observations ont permis de détecter le premier Pevatron Galactique et de poser les contraintes les plus fortes à ce jour sur la section efficace d'annihilation de particules de matière noire dans la plage en masse du TeV. Le futur observatoire CTA sera déployé sur deux sites, l'un à La Palma et l'autre au Chili. Ce dernier composé de plus de 50 télescopes permettra d'obtenir un balayage sans précédent de la région sur Centre Galactique.
Le travail proposé portera sur l'analyse et l'interprétation des observations H.E.S.S. conduites dans la région du Centre Galactique pour la recherche d'émission diffuses (populations de sources non résolues, matière noire massive) ainsi que des observations menées vers un sélection de noyaux actifs de galaxie pour la recherche d'ALPs constituant la matière noire. Ces nouveaux cadres d'analyses seront implémentés pour les analyses CTA à venir. Une implication dans la commissioning des premiers MSTs au Chili ainsi que dans l’analyse des premières données scientifiques sont attendues.

Modélisation de la réponse instrumentale des télescopes spatiaux avec un modèle optique différentiable

Contexte

L'effet de lentille gravitationnelle faible [1] est une sonde puissante de la structure à grande échelle de notre univers. Les cosmologistes utilisent l'effet de lentille faible pour étudier la nature de la matière noire et sa distribution spatiale. Les missions d'observation de l'effet de lentille faible nécessitent des mesures très précises de la forme des images de galaxies. La réponse instrumentale du télescope, appelée fonction d'étalement du point (PSF), produit une déformation des images observées. Cette déformation peut être confondue avec les effets d'un faible effet de lentille sur les images de galaxies, ce qui constitue l'une des principales sources d'erreur systématique lors de la recherche sur les faibles effets de lentille. Par conséquent, l'estimation d'un modèle de PSF fiable et précis est cruciale pour le succès de toute mission de faible lentille [2]. Le champ de la PSF peut être interprété comme un noyau convolutionnel qui affecte chacune de nos observations d'intérêt, qui varie spatialement, spectralement et temporellement. Le modèle de la PSF doit être capable de gérer chacune de ces variations. Nous utilisons des étoiles spécifiques considérées comme des sources ponctuelles dans le champ de vision pour contraindre notre modèle PSF. Ces étoiles, qui sont des objets non résolus, nous fournissent des échantillons dégradés du champ de la PSF. Les observations subissent différentes dégradations en fonction des propriétés du télescope. Ces dégradations comprennent le sous-échantillonnage, l'intégration sur la bande passante de l'instrument et le bruit additif. Nous construisons finalement le modèle de la PSF en utilisant ces observations dégradées et utilisons ensuite le modèle pour déduire la PSF à la position des galaxies. Cette procédure constitue le problème inverse mal posé de la modélisation de la PSF. Voir [3] pour un article récent sur la modélisation de la PSF.

La mission Euclid récemment lancée représente l'un des défis les plus complexes pour la modélisation de la PSF. En raison de la très large bande passante de l'imageur visible (VIS) d'Euclid, allant de 550 nm à 900 nm, les modèles de PSF doivent capturer non seulement les variations spatiales du champ de PSF, mais aussi ses variations chromatiques. Chaque observation d'étoile est intégrée avec la distribution d'énergie spectrale (SED) de l'objet sur l'ensemble de la bande passante du VIS. Comme les observations sont sous-échantillonnées, une étape de super-résolution est également nécessaire. Un modèle récent appelé WaveDiff [4] a été proposé pour résoudre le problème de modélisation de la PSF pour Euclid et est basé sur un modèle optique différentiable. WaveDiff a atteint des performances de pointe et est en train d'être testé avec des observations récentes de la mission Euclid.

Le télescope spatial James Webb (JWST) a été lancé récemment et produit des observations exceptionnelles. La collaboration COSMOS-Web [5] est un programme à grand champ du JWST qui cartographie un champ contigu de 0,6 deg2. Les observations de COSMOS-Web sont disponibles et offrent une occasion unique de tester et de développer un modèle précis de PSF pour le JWST. Dans ce contexte, plusieurs cas scientifiques, en plus des études de lentille gravitationnelle faible, peuvent grandement bénéficier d'un modèle PSF précis. Par exemple, l'effet de lentille gravitationnel fort [6], où la PSF joue un rôle crucial dans la reconstruction, et l'imagerie des exoplanètes [7], où les speckles de la PSF peuvent imiter l'apparence des exoplanètes, donc la soustraction d'un modèle de PSF exact et précis est essentielle pour améliorer l'imagerie et la détection des exoplanètes.

Projet de doctorat

Le candidat visera à développer des modèles PSF plus précis et plus performants pour les télescopes spatiaux en exploitant un cadre optique différentiable et concentrera ses efforts sur Euclid et le JWST.

Le modèle WaveDiff est basé sur l'espace du front d'onde et ne prend pas en compte les effets au niveau du pixel ou du détecteur. Ces erreurs au niveau des pixels ne peuvent pas être modélisées avec précision dans le front d'onde car elles se produisent naturellement directement sur les détecteurs et ne sont pas liées aux aberrations optiques du télescope. Par conséquent, dans un premier temps, nous étendrons l'approche de modélisation de la PSF en tenant compte de l'effet au niveau du détecteur en combinant une approche paramétrique et une approche basée sur les données (apprises). Nous exploiterons les capacités de différenciation automatique des cadres d'apprentissage automatique (par exemple TensorFlow, Pytorch, JAX) du modèle WaveDiff PSF pour atteindre l'objectif.

Dans une deuxième direction, nous envisagerons l'estimation conjointe du champ de la PSF et des densités d'énergie spectrale (SED) stellaires en exploitant des expositions répétées ou des dithers. L'objectif est d'améliorer et de calibrer l'estimation originale de la SED en exploitant les informations de modélisation de la PSF. Nous nous appuierons sur notre modèle PSF, et les observations répétées du même objet changeront l'image de l'étoile (puisqu'elle est imagée sur différentes positions du plan focal) mais partageront les mêmes SED.

Une autre direction sera d'étendre WaveDiff à des observatoires astronomiques plus généraux comme le JWST avec des champs de vision plus petits. Nous devrons contraindre le modèle de PSF avec des observations de plusieurs bandes pour construire un modèle de PSF unique contraint par plus d'informations. L'objectif est de développer le prochain modèle de PSF pour le JWST qui soit disponible pour une utilisation généralisée, que nous validerons avec les données réelles disponibles du programme COSMOS-Web JWST.

La direction suivante sera d'étendre les performances de WaveDiff en incluant un champ continu sous la forme d'une représentation neuronale implicite [8], ou de champs neuronaux (NeRF) [9], pour traiter les variations spatiales de la PSF dans l'espace du front d'onde avec un modèle plus puissant et plus flexible.

Enfin, tout au long de son doctorat, le candidat collaborera à l'effort de modélisation de la PSF par les données d'Euclid, qui consiste à appliquer WaveDiff aux données réelles d'Euclid, et à la collaboration COSMOS-Web pour exploiter les observations du JWST.

Références

[1] R. Mandelbaum. “Weak Lensing for Precision Cosmology”. In: Annual Review of Astronomy and Astro- physics 56 (2018), pp. 393–433. doi: 10.1146/annurev-astro-081817-051928. arXiv: 1710.03235.
[2] T. I. Liaudat et al. “Multi-CCD modelling of the point spread function”. In: A&A 646 (2021), A27. doi:10.1051/0004-6361/202039584.
[3] T. I. Liaudat, J.-L. Starck, and M. Kilbinger. “Point spread function modelling for astronomical telescopes: a review focused on weak gravitational lensing studies”. In: Frontiers in Astronomy and Space Sciences 10 (2023). doi: 10.3389/fspas.2023.1158213.
[4] T. I. Liaudat, J.-L. Starck, M. Kilbinger, and P.-A. Frugier. “Rethinking data-driven point spread function modeling with a differentiable optical model”. In: Inverse Problems 39.3 (Feb. 2023), p. 035008. doi:10.1088/1361-6420/acb664.
[5] C. M. Casey et al. “COSMOS-Web: An Overview of the JWST Cosmic Origins Survey”. In: The Astrophysical Journal 954.1 (Aug. 2023), p. 31. doi: 10.3847/1538-4357/acc2bc.
[6] A. Acebron et al. “The Next Step in Galaxy Cluster Strong Lensing: Modeling the Surface Brightness of Multiply Imaged Sources”. In: ApJ 976.1, 110 (Nov. 2024), p. 110. doi: 10.3847/1538-4357/ad8343. arXiv: 2410.01883 [astro-ph.GA].
[7] B. Y. Feng et al. “Exoplanet Imaging via Differentiable Rendering”. In: IEEE Transactions on Computational Imaging 11 (2025), pp. 36–51. doi: 10.1109/TCI.2025.3525971.
[8] Y. Xie et al. “Neural Fields in Visual Computing and Beyond”. In: arXiv e-prints, arXiv:2111.11426 (Nov.2021), arXiv:2111.11426. doi: 10.48550/arXiv.2111.11426. arXiv: 2111.11426 [cs.CV].
[9] B. Mildenhall et al. “NeRF: Representing Scenes as Neural Radiance Fields for View Synthesis”. In: arXiv e-prints, arXiv:2003.08934 (Mar. 2020), arXiv:2003.08934. doi: 10.48550/arXiv.2003.08934. arXiv:2003.08934 [cs.CV].

À la recherche des exo-aurores

Les aurores sont des phénomènes optiques bien connus dans les planètes du système solaire. Elles ont une grande valeur diagnostique, car leurs émissions révèlent la composition atmosphérique des planètes, la présence de champs magnétiques et les conditions du vent solaire à l'orbite de la planète. La recherche d'aurores sur les exoplanètes et les naines brunes constitue la prochaine étape. Une première avancée dans cette direction a récemment eu lieu, avec la détection d'une émission de CH4 attribuée à l'excitation aurorale sur la naine brune W1935. Cette détection, ainsi que la perspective d'observer d'autres caractéristiques aurorales avec les télescopes existants et à venir, sont à l'origine de ce projet. Nous construirons notamment le premier modèle dédié à l'étude des émissions aurorales de CH4 et H3+ sur les exoplanètes et les naines brunes. Ce modèle sera utilisé pour étudier les conditions sur W1935 et pour prédire la détectabilité des aurores sur d'autres objets substellaires.

Comprendre l’origine de la remarquable efficacité de formation des galaxies lointaines

Le télescope spatial James Webb révolutionne notre compréhension de l’univers lointain. Un résultat s’impose qui questionne nos modèles : la très grande efficacité de formation d’étoiles des galaxies lointaines. Mais ce constat est dérivé de manière indirecte : on mesure la masse d’étoiles dans les galaxies, pas leur taux de formation d’étoiles. C’est la principale faiblesse du James Webb. Le but de cette thèse est de remédier à cette faiblesse du James Webb en utilisant sa capacité de résolution angulaire, qui n’a jusqu’ici pas été prise en compte afin d’obtenir une mesure plus robuste du SFR des galaxies distantes. On en déduira une loi qui permettra d’améliorer la robustesse de la détermination du SFR grâce aux propriétés morphologiques et en combinant les données du James Webb avec celles d’ALMA (z=1-3). Puis on l’appliquera à l’univers lointain (z=3-6, 2e partie) et on l’utilisera comme benchmark pour les simulations numériques (3e partie).

Exploration des exoplanètes rocheuses avec le JWST

L’un des principaux objectifs du télescope spatial JWST est de caractériser, pour la première fois, les atmosphères des exoplanètes rocheuses et tempérées, une étape clé dans la recherche de mondes potentiellement habitables. Les exoplanètes rocheuses tempérées accessibles au JWST sont principalement celles en orbite autour d’étoiles de type M. Cependant, une question majeure demeure quant à la capacité des planètes orbitant autour de ces étoiles froides à conserver leur atmosphère. En 2024, un programme exceptionnel de 500 heures de temps discrétionnaire du directeur (Director’s Discretionary Time, DDT), intitulé Rocky Worlds, a été consacré à cette thématique, soulignant son importance stratégique au plus haut niveau (NASA, STScI).
L’objectif principal de ce projet de doctorat est : 1) Analyser l’ensemble des données d’éclipses JWST/MIRI disponibles pour les exoplanètes rocheuses issues du programme Rocky Worlds et d’autres programmes publics, en utilisant un cadre d’analyse cohérent et homogène ; 2) Rechercher des tendances à l’échelle des populations dans les observations et de les interpréter à l’aide de simulations atmosphériques tridimensionnelles.
À travers ce travail, nous visons à identifier les processus physiques qui régissent la présence et la composition des atmosphères des exoplanètes rocheuses tempérées.

Etude numérique de la turbulence interstellaire à l'heure de l'exascale

Ce projet de thèse vise à mieux comprendre la turbulence du milieu interstellaire, un phénomène clé pour la formation des étoiles et des structures galactiques. Cette turbulence, à la fois magnétisée, supersonique et multiphasique, influence la manière dont l’énergie se propage et se dissipe, régulant ainsi l’efficacité de la formation stellaire à travers l’histoire de l’Univers. Son étude est complexe car elle implique une vaste gamme d’échelles spatiales et temporelles, difficile à reproduire numériquement. Les progrès du calcul haute performance, notamment l’arrivée des supercalculateurs exascale à GPU permet désormais d’envisager des simulations beaucoup plus fines.

Le code Dyablo, développé à l’IRFU, sera utilisé pour réaliser des simulations tridimensionnelles de très grande taille, avec un maillage adaptatif pour affiner les zones de dissipation d’énergie. L’étude progressera par étapes : d’abord des écoulements simples et isothermes, puis des modèles incluant chauffage, refroidissement, champ magnétique et gravité. Les propriétés turbulentes seront analysées via spectres de puissance, fonctions de structure et distributions de densité, afin de mieux comprendre la formation des zones denses propices à la naissance des étoiles. Enfin, une extension du travail à l’échelle galactique, en collaboration avec d’autres instituts français, visera à explorer la cascade d’énergie turbulente à grande échelle dans les galaxies entières.

Cosmologie avec la forêt Lyman-alpha du grand relevé cosmologique DESI.

La distribution de matière à grande échelle dans l'univers est utilisée pour tester nos modèles cosmologiques. On utilise pour cela avant tout les oscillations acoustiques de baryons (BAO) mesurée dans la fonction de corrélation à deux points de cette distribution. Cependant l'ensemble du champ de matière contient des informations à diverses échelles, qui permettent de mieux contraindre nos modèles que le BAO seul. A redshift z > 2, la meilleure sonde de cette distribution de matière est la forêt Lyman-alpha, un ensemble de raies d'absorption mesurées dans les spectres de sources lointaines. Le grand relevé spectroscopique DESI a collecté environ un million de ces spectres. Avec un lot partiel de données "DR2", nous avons déjà mesuré le BAO avec une précision de 0.7%, contraignant ainsi fortement le taux d'expansion de l'univers au cours des premiers milliards d'années de son évolution.

Le but de cette thèse est d'exploiter l'ensemble complet des données Lyman-alpha à grande échelle de DESI pour obtenir les meilleures contraintes possibles sur les modèles cosmologiques. Pour cela, dans une première étape l'étudiant appliquera pour la première fois une méthode dite de reconstruction, qui permet d'améliorer la précision de la mesure du BAO, en exploitant l'information du champ de densité de matière. Dans la suite de sa thèse, en lien avec des efforts similaires menés dans notre groupe avec les galaxies de DESI, il implémentera une nouvelle méthode dite de simulation-based inference: dans cette approche l'ensemble du champ de matière est utilisé directement pour estimer les paramètres cosmologiques en particulier l'énergie noire. L'étudiant apportera ainsi une pierre importante aux mesures cosmologiques finales de DESI avec le Lyman-alpha.

Cette thèse sera de préférence précédée d'un stage.

Estimation impartiale du cisaillement pour Euclid avec modélisation automatiquement différentiable et accélérée par GPU

Le projet de thèse porte sur la mesure sans biais du cisaillement gravitationnel faible, un effet dû à la déviation de la lumière des galaxies par la matière présente sur la ligne de visée. Cette technique est essentielle pour étudier la matière noire, l’énergie noire et la gravité, et constitue un pilier de la mission spatiale Euclid, lancée en 2023. Les méthodes classiques de mesure des formes de galaxies introduisent des biais systématiques dans l’estimation du cisaillement. L’objectif de cette thèse est de développer une approche innovante de modélisation directe permettant d’inférer le cisaillement sans passer par la mesure de forme, en simulant des images de galaxies réalistes à l’aide d’architectures d’apprentissage profond. Le doctorant participera à l’adaptation de cette méthode aux données réelles d’Euclid, en intégrant la complexité du système de traitement des données (SGS) et en optimisant le calcul sur GPU et supercalculateurs. Ce travail s’inscrit dans un contexte très dynamique, coïncidant avec la première diffusion publique des données d’Euclid prévue pour 2026. Les résultats attendus sont une estimation du cisaillement plus précise et robuste, ouvrant la voie à des analyses cosmologiques de nouvelle génération.

Analyse spectro-temporelle de l'émission rémanente des sursaut gamma cosmiques détectés avec SVOM

Les sursauts gamma cosmiques (Gamma-Ray Bursts, GRBs) sont les explosions les plus puissantes de l'Univers. Ils durent quelques dizaines de secondes et émettent autant d'énergie que le Soleil pendant toute sa durée de vie. Leur émission en rayons gamma est suivie d'une émission de longue durée (de quelques heures à plusieurs jours) allant des rayons X à la bande radio. Cette émission « rémanente » est riche en informations sur l'environnement proche du GRB et sur la galaxie hôte. SVOM (Space based astronomical Variable Object Monitor) est une mission sino-française dédiée à l'étude des GRB, qui a été lancée avec succès en juin 2024. Elle embarque une charge utile multi-longueurs d'onde couvrant les rayons gamma/rayons X/optiques et comprend deux télescopes robotiques terrestres dédiés au Mexique et en Chine.
Le projet de doctorat est axé sur l'exploitation des données SVOM pour les GRB. Le candidat retenu rejoindra l'équipe scientifique MXT au DAp. MXT est un nouveau type de télescope à rayons X, dont le DAp est responsable et dont le centre d'instrumentation est également hébergé au DAp.
Le doctorant participera activement à l'analyse spectrale et temporelle des données MXT. Ces données seront comparées
aux autres données acquises par la collaboration SVOM, en particulier dans les domaines optique et infrarouge.
Cet ensemble de données sera utilisé comme support à l'interprétation physique des sursauts gamma. Plus précisément, les aspects liés à la modélisation de l'injection d'énergie dans les premières phases de la rémanence seront utilisés pour déterminer la nature de l'objet compact à l'origine du flux relativiste, générant l'émission électromagnétique observée.

Top