Offre en lien avec l’Action/le Réseau : – — –/– — –

Laboratoire/Entreprise : Observatoire de Lyon
Durée : 36 mois
Contact : ferreol.soulez@univ-lyon1.fr
Date limite de publication : 2021-05-01

Contexte :
Natural sciences such as astrophysics, geophysics and nuclear physics often use numerical simulations to model highly complex physical systems. These simulations are now more and more accurate thanks to the computational power available. For example, 3D convection models can simulate the thermochemical evolution and structure of stars and planets.
However, to disentangle different models and to estimate physical parameters (e.g., initial conditions), the outputs of these simulations have to be compared to observations. This confrontation of simulations to observations is a major challenge in natural sciences. Indeed, numerical simulations are now able to model quite accurately objects that are impossible to observe directly (e.g., interior of stars and planets, stars and black holes environment …). As for the observations, although their quality and quantity are rapidly increasing, they are often only indirectly related to the actual parameters of interest (e.g., seismic waves observations are used to construct images of the earth mantle, measured interferometric visibilities are used to characterize planet forming disk …).
To infer simulation parameters from observations is very challenging. When a single simulation is computationally intensive, it is impossible to use either stochastic or continuous optimization methods to infer parameters. In most cases, one can only rely on finding the best fits on a low dimensional pre-computed grid of model parameters.

Sujet :
The ultimate goal of the proposed thesis is to build a fast interpolation method on a grid of computational physics simulated images (in a broad sense as it can also be 3D volumes or spectra). With such a method, we will quickly have an approximation of a simulated image from any possible set of parameters, without having to run the expensive simulation. It then will be possible to use any method (optimization, Bayesian inference) to derive the so sought-after distribution of parameters.
The main idea is to use a deep learning framework to build the interpolator. Indeed, all possible modeled images are concentrated on a lower-dimensional subspace or manifold. Deep neural networks such as Generative Adversarial Networks (GAN) [1] appear to be very efficient to model manifolds and could be much more efficient interpolators than classical polynomial interpolators. Trained on a grid on simulated images, these deep neural networks will produce continuous approximations of the images. As a toy example, in a properly defined manifold, the images of a single circle vary continuously with the circle radius. Interpolation between two images of circles with different radius must follow this manifold whereas any polynomial interpolation will produce an image with two circles rather than an image of a single circle with intermediate radius.
Grids of models are quite ubiquitous in physics, and hence such a project can have important impact. To ensure that it will be both robust and useful in practice, the deep learning based interpolator will be developed for two different applications: (i) planet forming disk characterization using VLTI in collaboration with J. Kluska (KU Leuven) and (ii) reconstruction of mantle structure based on geophysical surface observations.

Profil du candidat :
The candidates must hold a national master degree or equivalent in Applied Math, Signal/Image Processing or Machine Learning.
He/She should have a good coding proficiency ( Julia / Python) and a real interest in the underlying astrophysical or geophysical questions.

Formation et compétences requises :
The candidates must hold a national master degree or equivalent in Applied Math, Signal/Image Processing or Machine Learning.

Adresse d’emploi :
Observatoire de Lyon
9 avenue Charles André
Saint Genis Laval