No Arabic abstract
We present a scheme for using stellar catalogues to map the three-dimensional distributions of extinction and dust within our Galaxy. Extinction is modelled as a Gaussian random field, whose covariance function is set by a simple physical model of the ISM that assumes a Kolmogorov-like power spectrum of turbulent fluctuations. As extinction is modelled as a random field, the spatial resolution of the resulting maps is set naturally by the data available; there is no need to impose any spatial binning. We verify the validity of our scheme by testing it on simulated extinction fields and show that its precision is significantly improved over previous dust-mapping efforts. The approach we describe here can make use of any photometric, spectroscopic or astrometric data; it is not limited to any particular survey. Consequently, it can be applied to a wide range of data from both existing and future surveys.
Gaussian processes are the ideal tool for modelling the Galactic ISM, combining statistical flexibility with a good match to the underlying physics. In an earlier paper we outlined how they can be employed to construct three-dimensional maps of dust extinction from stellar surveys. Gaussian processes scale poorly to large datasets though, which put the analysis of realistic catalogues out of reach. Here we show how a novel combination of the Expectation Propagation method and certain sparse matrix approximations can be used to accelerate the dust mapping problem. We demonstrate, using simulated Gaia data, that the resultant algorithm is fast, accurate and precise. Critically, it can be scaled up to map the Gaia catalogue.
Selection effects can bedevil the inference of the properties of a population of astronomical catalogues, unavoidably biasing the observed catalogue. This is particularly true when mapping interstellar extinction in three dimensions: more extinguished stars are fainter and so generally less likely to appear in any magnitude limited catalogue of observations. This paper demonstrates how to account for this selection effect when mapping extinction, so that accurate and unbiased estimates of the true extinction are obtained. We advocate couching the description of the problem explicitly as a Poisson point process, which allows the likelihoods employed to be easily and correctly normalised in such a way that accounts for the selection functions applied to construct the catalogue of observations.
The Galaxy and the stars in it form a hierarchical system, such that the properties of individual stars are influenced by those of the Galaxy. Here, an approach is described which uses hierarchical Bayesian models to simultaneously and empirically determine the mean distance-extinction relationship for a sightline and the properties of stars which populate it. By exploiting the hierarchical nature of the problem, the method described is able to achieve significantly improved precision and accuracy with respect to previous 3D extinction mapping techniques. This method is not tied to any individual survey and could be applied to any observations, or combination of observations available. Furthermore, it is extendible and, in addition, could be employed to study Galactic structure as well as factors such as the initial mass function and star formation history in the Galaxy.
We present a three dimensional (3D) extinction analysis in the region toward the supernova remnant (SNR) S147 (G180.0-1.7) using multi-band photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC), 2MASS and WISE. We isolate a previously unrecognised dust structure likely to be associated with SNR S147. The structure, which we term as S147 dust cloud, is estimated to have a distance $d$ = 1.22 $pm$ 0.21 kpc, consistent with the conjecture that S147 is associated with pulsar PSR J0538 + 2817. The cloud includes several dense clumps of relatively high extinction that locate on the radio shell of S147 and coincide spatially with the CO and gamma-ray emission features. We conclude that the usage of CO measurements to trace the SNR associated MCs is unavoidably limited by the detection threshold, dust depletion, and the difficulty of distance estimates in the outer Galaxy. 3D dust extinction mapping may provide a better way to identify and study SNR-MC interactions.
We present version X of the hammurabi package, the HEALPix-based numeric simulator for Galactic polarized emission. Improving on its earlier design, we have fully renewed the framework with modern C++ standards and features. Multi-threading support has been built in to meet the growing computational workload in future research. For the first time, we present precision profiles of hammurabi line-of-sight integral kernel with multi-layer HEALPix shells. In addition to fundamental improvements, this report focuses on simulating polarized synchrotron emission with Gaussian random magnetic fields. Two fast methods are proposed for realizing divergence-free random magnetic fields either on the Galactic scale where a field alignment and strength modulation are imposed, or on a local scale where more physically motivated models like a parameterized magneto-hydrodynamic (MHD) turbulence can be applied. As an example application, we discuss the phenomenological implications of Gaussian random magnetic fields for high Galactic latitude synchrotron foregrounds. In this, we numerically find B/E polarization mode ratios lower than unity based on Gaussian realizations of either MHD turbulent spectra or in spatially aligned magnetic fields.