Issue |
A&A
Volume 658, February 2022
|
|
---|---|---|
Article Number | A162 | |
Number of page(s) | 14 | |
Section | Stellar atmospheres | |
DOI | https://doi.org/10.1051/0004-6361/202142027 | |
Published online | 16 February 2022 |
Approximate Bayesian neural Doppler imaging
1
Instituto de Astrofísica de Canarias,
C/Vía Láctea s/n,
38205
La Laguna,
Tenerife,
Spain
e-mail: andres.asensio@iac.es
2
Departamento de Astrofísica, Universidad de La Laguna,
38206
La Laguna,
Tenerife,
Spain
3
Institute for Solar Physics, Dept. of Astronomy, Stockholm University, AlbaNova University Centre,
10691
Stockholm,
Sweden
4
Department of Physics and Astronomy, Uppsala University,
Box 516,
Uppsala
75120,
Sweden
Received:
16
August
2021
Accepted:
9
December
2021
Aims. 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 because of the Doppler effect, 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 compatible with the observations.
Methods. 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, these normalizing flows provide a very efficient way of obtaining samples from the posterior distribution. The conditioning on observations is achieved through the use of Transformer encoders, which can deal with arbitrary wavelength sampling and rotation phases.
Results. Our model can produce thousands of posterior samples per second, each one accompanied by an estimation of the log-probability. Our exhaustive validation of the model for very high-signal-to-noise observations shows that it correctly approximates the posterior, albeit 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 for carrying out approximate Bayesian inference in more complex problems in stellar physics, such as constraining the magnetic properties using polarimetry.
Key words: stars: atmospheres / stars: activity / line: profiles / methods: data analysis / stars: individual: II Peg
© ESO 2022
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.