Approximate Bayesian Neural Doppler Imaging

08/20/2021
by   A. Asensio Ramos, et al.
6

The non-uniform surface temperature distribution of rotating active stars is routinely mapped with the Doppler Imaging technique. Inhomogeneities in the surface produce features in high-resolution spectroscopic observations that shift in wavelength depending on their position on the visible hemisphere. The inversion problem has been systematically solved using maximum a-posteriori regularized methods assuming smoothness or maximum entropy. Our aim in this work is to solve the full Bayesian inference problem, by providing access to the posterior distribution of the surface temperature in the star. We use amortized neural posterior estimation to produce a model that approximates the high-dimensional posterior distribution for spectroscopic observations of selected spectral ranges sampled at arbitrary rotation phases. The posterior distribution is approximated with conditional normalizing flows, which are flexible, tractable and easy to sample approximations to arbitrary distributions. When conditioned on the spectroscopic observations, they provide a very efficient way of obtaining samples from the posterior distribution. The conditioning on observations is obtained through the use of Transformer encoders, which can deal with arbitrary wavelength sampling and rotation phases. Our model can produce thousands of posterior samples per second. Our validation of the model for very high signal-to-noise observations shows that it correctly approximates the posterior, although with some overestimation of the broadening. We apply the model to the moderately fast rotator II Peg, producing the first Bayesian map of its temperature inhomogenities. We conclude that conditional normalizing flows are a very promising tool to carry out approximate Bayesian inference in more complex problems in stellar physics, like constraining the magnetic properties.

READ FULL TEXT

page 3

page 7

page 9

page 11

page 12

research
04/21/2023

Machine Learning and the Future of Bayesian Computation

Bayesian models are a powerful tool for studying complex data, allowing ...
research
09/13/2019

Adversarial α-divergence Minimization for Bayesian Approximate Inference

Neural networks are popular models for regression. They are often traine...
research
02/12/2021

Stochastic Convergence Rates and Applications of Adaptive Quadrature in Bayesian Inference

We provide the first stochastic convergence rates for adaptive Gauss–Her...
research
02/19/2019

Scalable Thompson Sampling via Optimal Transport

Thompson sampling (TS) is a class of algorithms for sequential decision-...
research
10/04/2022

Amortized Bayesian Inference of GISAXS Data with Normalizing Flows

Grazing-Incidence Small-Angle X-ray Scattering (GISAXS) is a modern imag...
research
06/22/2020

Constraining subglacial processes from surface velocity observations using surrogate-based Bayesian inference

Basal motion is the primary mechanism for ice flux outside Antarctica, y...
research
11/06/2017

Fast amortized inference of neural activity from calcium imaging data with variational autoencoders

Calcium imaging permits optical measurement of neural activity. Since in...

Please sign up or login with your details

Forgot password? Click here to reset