Approximate Bayesian neural Doppler imaging

Asensio Ramos, A.; Díaz Baso, C. J.; Kochukhov, O.
Bibliographical reference

Astronomy and Astrophysics

Advertised on:
2
2022
Number of authors
3
IAC number of authors
1
Citations
9
Refereed citations
7
Description

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.
Related projects
Project Image
Solar and Stellar Magnetism
Magnetic fields are at the base of star formation and stellar structure and evolution. When stars are born, magnetic fields brake the rotation during the collapse of the mollecular cloud. In the end of the life of a star, magnetic fields can play a key role in the form of the strong winds that lead to the last stages of stellar evolution. During
Tobías
Felipe García
Project Image
Magnetism, Polarization and Radiative Transfer in Astrophysics
Magnetic fields pervade all astrophysical plasmas and govern most of the variability in the Universe at intermediate time scales. They are present in stars across the whole Hertzsprung-Russell diagram, in galaxies, and even perhaps in the intergalactic medium. Polarized light provides the most reliable source of information at our disposal for the
Tanausú del
Pino Alemán