Showing 1–50 of 326 results for author: Enßlin, T A

Spatially Coherent 3D Distributions of HI and CO in the Milky Way

Authors: Laurin Söding, Gordian Edenhofer, Torsten A. Enßlin, Philipp Frank, Ralf Kissmann, Vo Hong Minh Phan, Andrés Ramírez, Hanieh Zhandinejad, Philipp Mertsch

Abstract: The spatial distribution of the gaseous components of the Milky Way is of great importance for a number of different fields, e.g. Galactic structure, star formation and cosmic rays. However, obtaining distance information to gaseous clouds in the interstellar medium from Doppler-shifted line emission is notoriously difficult given our unique vantage point in the Galaxy. It requires precise knowled… ▽ More The spatial distribution of the gaseous components of the Milky Way is of great importance for a number of different fields, e.g. Galactic structure, star formation and cosmic rays. However, obtaining distance information to gaseous clouds in the interstellar medium from Doppler-shifted line emission is notoriously difficult given our unique vantage point in the Galaxy. It requires precise knowledge of gas velocities and generally suffers from distance ambiguities. Previous works often assumed the optically thin limit (no absorption), a fixed velocity field, and lack resolution overall. We aim to overcome these issues and improve previous reconstructions of the gaseous constituents of the interstellar medium of the Galaxy. We use 3D Gaussian processes to model correlations in the interstellar medium, including correlations between different lines of sight, and enforce a spatially coherent structure in the prior. For modelling the transport of radiation from the emitting gas to us as observers, we take absorption effects into account. A special numerical grid ensures high resolution nearby. We infer the spatial distributions of HI, CO, their emission line-widths, and the Galactic velocity field in a joint Bayesian inference. We further constrain these fields with complementary data from Galactic masers and young stellar object clusters. Our main result consists of a set of samples that implicitly contain statistical uncertainties. The resulting maps are spatially coherent and reproduce the data with high fidelity. We confirm previous findings regarding the warping and flaring of the Galactic disc. A comparison with 3D dust maps reveals a good agreement on scales larger than approximately 400 pc. While our results are not free of artefacts, they present a big step forward in obtaining high quality 3D maps of the interstellar medium. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 21 pages, 22 figures

The influence of the 3D Galactic gas structure on cosmic-ray transport and gamma-ray emission

Authors: Andrés Ramírez, Gordian Edenhofer, Torsten A. Enßlin, Philipp Frank, Philipp Mertsch, Vo Hong Minh Phan, Laurin Söding, Hanieh Zhandinejad, Ralf Kissmann

Abstract: Cosmic rays (CRs) play a major role in the dynamics of the interstellar medium (ISM). Their interactions and transport ionize, heat, and push the ISM thereby coupling different regions of it. The spatial distribution of CRs depends on the distribution of their sources as well as the ISM constituents they interact with, such as gas, starlight, and magnetic fields. Particularly, gas interacts closel… ▽ More Cosmic rays (CRs) play a major role in the dynamics of the interstellar medium (ISM). Their interactions and transport ionize, heat, and push the ISM thereby coupling different regions of it. The spatial distribution of CRs depends on the distribution of their sources as well as the ISM constituents they interact with, such as gas, starlight, and magnetic fields. Particularly, gas interacts closely with CRs, influencing CR fluxes and gamma -ray emission. We illustrate the influence of 3D gas structures on CR transport and gamma -ray emission. We use the PICARD code and multiple samples of recent 3D reconstructions of the HI and H$_2$ Galactic gas constituents to investigate the impact on the transport of CRs and emission of gamma -rays. We find the necessary transport parameters to reproduce local measurements of CR fluxes, and see that they depend on the local distribution of gas density and structure. The distribution of CR fluxes exhibits energy-dependent structures that vary for all CR species due to their corresponding loss processes. Regions of enhanced secondary (primary) species are spatially correlated (anti-correlated) with the gas density. We observe a high sensitivity of the gamma -ray emission on the contrast of gas structures, as those determine the 3D spatial distributions of hadronic interactions and bremsstrahlung. We find that corresponding gas-induced structures in the distribution of CR electrons are also visible in Inverse Compton (IC) emission. Due to the aforementioned sensitivity, the analysis of CR data for CR sources and transport parameters requires the usage of accurate 3D gas maps. △ Less

Submitted 2 July, 2024; originally announced July 2024.

fast-resolve: Fast Bayesian Radio Interferometric Imaging

Authors: Jakob Roth, Philipp Frank, Hertzog L. Bester, Oleg M. Smirnov, Rüdiger Westermann, Torsten A. Enßlin

Abstract: Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algo… ▽ More Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algorithm resolve, enabling the application of Bayesian imaging for larger datasets. Methods: By combining computational shortcuts of the CLEAN algorithm with the Bayesian imaging algorithm resolve we developed an accurate and fast imaging algorithm which we name fast-resolve. Results: We validate the accuracy of the presented fast-resolve algorithm by comparing it with results from resolve on VLA Cygnus A data. Furthermore, we demonstrate the computational advantages of fast-resolve on a large MeerKAT ESO 137-006 dataset which is computationally out of reach for resolve. Conclusions: The presented algorithm is significantly faster than previous Bayesian imaging algorithms, broadening the applicability of Bayesian interferometric imaging. Specifically for the single channel VLA Cygnus A datasets fast-resolve is about $144$ times faster than resolve. For the MeerKAT dataset with multiple channels the computational speedup of fast-resolve is even larger. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 24 pages, 14 figures

The "C": The large Chameleon-Musca-Coalsack cloud

Authors: Gordian Edenhofer, João Alves, Catherine Zucker, Laura Posch, Torsten A. Enßlin

Abstract: Recent advancements in 3D dust mapping have transformed our understanding of the Milky Way's local interstellar medium, enabling us to explore its structure in three spatial dimensions for the first time. In this Letter, we use the most recent 3D dust map by Edenhofer et al. to study the well-known Chameleon, Musca, and Coalsack cloud complexes, located about 200 pc from the Sun. We find that thes… ▽ More Recent advancements in 3D dust mapping have transformed our understanding of the Milky Way's local interstellar medium, enabling us to explore its structure in three spatial dimensions for the first time. In this Letter, we use the most recent 3D dust map by Edenhofer et al. to study the well-known Chameleon, Musca, and Coalsack cloud complexes, located about 200 pc from the Sun. We find that these three complexes are not isolated but rather connect to form a surprisingly well-defined half-ring, constituting a single C-shaped cloud with a radius of about 50 pc, a thickness of about 45 pc, and a total mass of about $5 \times 10^{4} \mathrm{M}_{\odot}$, or $9 \times 10^{4} \mathrm{M}_{\odot}$ if including everything in the vicinity of the C-shaped cloud. Despite the absence of an evident feedback source at its center, the dynamics of young stellar clusters associated with the C structure suggest that a single supernova explosion about 4 Myr to 10 Myr ago likely shaped this structure. Our findings support a single origin story for these cloud complexes, suggesting that they were formed by feedback-driven gas compression, and offer new insights into the processes that govern the birth of star-forming clouds in feedback-dominated regions, such as the Scorpius-Centaurus association. △ Less

Submitted 29 June, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: 6 pages, 3 figures, to be published in Astronomy & Astrophysics (A&A)

Re-Envisioning Numerical Information Field Theory (NIFTy.re): A Library for Gaussian Processes and Variational Inference

Authors: Gordian Edenhofer, Philipp Frank, Jakob Roth, Reimar H. Leike, Massin Guerdi, Lukas I. Scheel-Platz, Matteo Guardiani, Vincent Eberle, Margret Westerkamp, Torsten A. Enßlin

Abstract: Imaging is the process of transforming noisy, incomplete data into a space that humans can interpret. NIFTy is a Bayesian framework for imaging and has already successfully been applied to many fields in astrophysics. Previous design decisions held the performance and the development of methods in NIFTy back. We present a rewrite of NIFTy, coined NIFTy.re, which reworks the modeling principle, ext… ▽ More Imaging is the process of transforming noisy, incomplete data into a space that humans can interpret. NIFTy is a Bayesian framework for imaging and has already successfully been applied to many fields in astrophysics. Previous design decisions held the performance and the development of methods in NIFTy back. We present a rewrite of NIFTy, coined NIFTy.re, which reworks the modeling principle, extends the inference strategies, and outsources much of the heavy lifting to JAX. The rewrite dramatically accelerates models written in NIFTy, lays the foundation for new types of inference machineries, improves maintainability, and enables interoperability between NIFTy and the JAX machine learning ecosystem. △ Less

Submitted 15 June, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: 10 pages, 2 figures, published in Journal of Open Source Software (JOSS)

Journal ref: Journal of Open Source Software, volume 9(98), year 2024, page 6593

Sardinia Radio Telescope observations of the Coma Cluster

Authors: M. Murgia, F. Govoni, V. Vacca, F. Loi, L. Feretti, G. Giovannini, A. Melis, R. Concu, E. Carretti, S. Poppi, G. Valente, A. Bonafede, G. Bernardi, W. Boschin, M. Brienza, T. E. Clarke, F. de Gasperin, T. A. Ensslin, C. Ferrari, F. Gastaldello, M. Girardi, L. Gregorini, M. Johnston-Hollitt, E. Orru', P. Parma , et al. (3 additional authors not shown)

Abstract: We present deep total intensity and polarization observations of the Coma cluster at 1.4 and 6.6 GHz performed with the Sardinia Radio Telescope. By combining the single-dish 1.4 GHz data with archival Very Large Array observations we obtain new images of the central radio halo and of the peripheral radio relic where we properly recover the brightness from the large scale structures. At 6.6 GHz we… ▽ More We present deep total intensity and polarization observations of the Coma cluster at 1.4 and 6.6 GHz performed with the Sardinia Radio Telescope. By combining the single-dish 1.4 GHz data with archival Very Large Array observations we obtain new images of the central radio halo and of the peripheral radio relic where we properly recover the brightness from the large scale structures. At 6.6 GHz we detect both the relic and the central part of the halo in total intensity and polarization. These are the highest frequency images available to date for these radio sources in this galaxy cluster. In the halo, we find a localized spot of polarized signal, with fractional polarization of about 45%. The polarized emission possibly extends along the north-east side of the diffuse emission. The relic is highly polarized, up to 55%, as usually found for these sources. We confirm the halo spectrum is curved, in agreement with previous single-dish results. The spectral index is alpha=1.48 +/- 0.07 at a reference frequency of 1 GHz and varies from alpha ~1.1, at 0.1 GHz, up to alpha ~ 1.8, at 10 GHz. We compare the Coma radio halo surface brightness profile at 1.4 GHz (central brightness and e-folding radius) with the same properties of the other halos, and we find that it has one of the lowest emissivities observed so far. Reanalyzing the relic's spectrum in the light of the new data, we obtain a refined radio Mach number of M=2.9 +/- 0.1. △ Less

Submitted 11 February, 2024; originally announced February 2024.

Comments: 27 pages, 18 figures, this article has been accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society

RMTable2023 and PolSpectra2023: standards for reporting polarization and Faraday rotation measurements of radio sources

Authors: C. L. Van Eck, B. M. Gaensler, S. Hutschenreuter, J. Livingston, Y. K. Ma, C. J. Riseley, A. J. M. Thomson, B. Adebahr, A. Basu, M. Birkinshaw, T. A. Ensslin, G. Heald, S. A. Mao, N. M. McClure-Griffiths

Abstract: Faraday rotation measures (RMs) have been used for many studies of cosmic magnetism, and in most cases having more RMs is beneficial for those studies. This has lead to development of RM surveys that have produced large catalogs, as well as meta-catalogs collecting RMs from many different publications. However, it has been difficult to take full advantage of all these RMs as the individual catalog… ▽ More Faraday rotation measures (RMs) have been used for many studies of cosmic magnetism, and in most cases having more RMs is beneficial for those studies. This has lead to development of RM surveys that have produced large catalogs, as well as meta-catalogs collecting RMs from many different publications. However, it has been difficult to take full advantage of all these RMs as the individual catalogs have been published in many different places, and in many different formats. In addition, the polarization spectra used to determine these RMs are rarely published, limiting the ability to re-analyze data as new methods or additional observations become available. We propose a standard convention for RM catalogs, RMTable2023, and a standard for source-integrated polarized spectra of radio sources, PolSpectra2023. These standards are intended to maximize the value and utility of these data for researchers and to make them easier to access. To demonstrate the use of the RMTable2023 standard, we have produced a consolidated catalog of 55 819 RMs collected from 42 published catalogs. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 44 pages, 5 figures, accepted by ApJS

Bayesian radio interferometric imaging with direction-dependent calibration

Authors: Jakob Roth, Philipp Arras, Martin Reinecke, Richard A. Perley, Rüdiger Westermann, Torsten A. Enßlin

Abstract: Context: Radio interferometers measure frequency components of the sky brightness, modulated by the gains of the individual radio antennas. Due to atmospheric turbulence and variations in the operational conditions of the antennas these gains fluctuate. Thereby the gains do not only depend on time but also on the spatial direction on the sky. To recover high quality radio maps an accurate reconstr… ▽ More Context: Radio interferometers measure frequency components of the sky brightness, modulated by the gains of the individual radio antennas. Due to atmospheric turbulence and variations in the operational conditions of the antennas these gains fluctuate. Thereby the gains do not only depend on time but also on the spatial direction on the sky. To recover high quality radio maps an accurate reconstruction of the direction and time-dependent individual antenna gains is required. Aims: This paper aims to improve the reconstruction of radio images, by introducing a novel joint imaging and calibration algorithm including direction-dependent antenna gains. Methods: Building on the \texttt{resolve} framework, we designed a Bayesian imaging and calibration algorithm utilizing the image domain gridding method for numerically efficient application of direction-dependent antenna gains. Furthermore by approximating the posterior probability distribution with variational inference, our algorithm can provide reliable uncertainty maps. Results: We demonstrate the ability of the algorithm to recover high resolution high dynamic range radio maps from VLA data of the radio galaxy Cygnus A. We compare the quality of the recovered images with previous work relying on classically calibrated data. Furthermore we compare with a compressed sensing algorithm also incorporating direction-dependent gains. Conclusions: Including direction-dependent effects in the calibration model significantly improves the dynamic range of the reconstructed images compared to reconstructions from classically calibrated data. Compared to the compressed sensing reconstruction, the resulting sky images have a higher resolution and show fewer artifacts. For utilizing the full potential of radio interferometric data, it is essential to consider the direction dependence of the antenna gains. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures

Journal ref: A&A 678, A177 (2023)

Disentangling the Faraday rotation sky

Authors: Sebastian Hutschenreuter, Marijke Haverkorn, Philipp Frank, Nergis C. Raycheva, Torsten A. Enßlin

Abstract: Magnetic fields permeate the diffuse interstellar medium (ISM) of the Milky Way, and are essential to explain the dynamical evolution and current shape of the Galaxy. Magnetic fields reveal themselves via their influence on the surrounding matter, and as such are notoriously hard to measure independently of other tracers. In this work, we attempt to disentangle an all sky map of the line-of-sight… ▽ More Magnetic fields permeate the diffuse interstellar medium (ISM) of the Milky Way, and are essential to explain the dynamical evolution and current shape of the Galaxy. Magnetic fields reveal themselves via their influence on the surrounding matter, and as such are notoriously hard to measure independently of other tracers. In this work, we attempt to disentangle an all sky map of the line-of-sight parallel component of the Galactic magnetic field from the Faraday effect, utilizing several tracers of the Galactic thermal electron density. Additionally, we aim to produce a Galactic electron dispersion measure map and quantify several tracers of the structure of the ionized medium of the Milky Way. We rely on compiled catalogs of extragalactic Faraday rotation measures and Galactic pulsar dispersion measures, a well as data on bremsstrahlung and the hydrogen $α$ spectral line to trace the ionized medium of the Milky Way. We present the first full sky map of the line-of-sight averaged Galactic magnetic field. Within this map, we find LoS parallel and LoS-averaged magnetic field strengths of up to 4 $μ$G, with an all-sky root-mean-square of 1.1 $μ$G, which is consistent with previous local measurements and global magnetic field models. Additionally, we produce a detailed electron dispersion measure map, which agrees with already existing parametric models at high latitudes, but suffers from systematic effects in the disk. Further analysis of our results with regard to the 3D structure of $n_{th}$ reveals that it follows a Kolmogorov-type turbulence for most of the sky. From the reconstructed dispersion measure and emission measure maps we construct several tracers of variability of $n_{th}$ along the LoS. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 24 pages, 16 Figures, submitted to A&A

Fast-Cadence High-Contrast Imaging with Information Field Theory

Authors: Jakob Roth, Gianluca Li Causi, Vincenzo Testa, Philipp Arras, Torsten A. Enßlin

Abstract: Although many exoplanets have been indirectly detected over the last years, direct imaging of them with ground-based telescopes remains challenging. In the presence of atmospheric fluctuations, it is ambitious to resolve the high brightness contrasts at the small angular separation between the star and its potential partners. Post-processing of telescope images has become an essential tool to impr… ▽ More Although many exoplanets have been indirectly detected over the last years, direct imaging of them with ground-based telescopes remains challenging. In the presence of atmospheric fluctuations, it is ambitious to resolve the high brightness contrasts at the small angular separation between the star and its potential partners. Post-processing of telescope images has become an essential tool to improve the resolvable contrast ratios. This paper contributes a post-processing algorithm for fast-cadence imaging, which deconvolves sequences of telescope images. The algorithm infers a Bayesian estimate of the astronomical object as well as the atmospheric optical path length, including its spatial and temporal structures. For this, we utilize physics-inspired models for the object, the atmosphere, and the telescope. The algorithm is computationally expensive but allows to resolve high contrast ratios despite short observation times and no field rotation. We test the performance of the algorithm with point-like companions synthetically injected into a real data set acquired with the SHARK-VIS pathfinder instrument at the LBT telescope. Sources with brightness ratios down to $6\cdot10^{-4}$ to the star are detected at $185$ mas separation with a short observation time of $0.6\,\text{s}$. △ Less

Submitted 7 February, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 12 pages, 6 figures

Journal ref: AJ 165 86 (2023)

Puzzling large-scale polarization in the galaxy cluster Abell 523

Authors: Valentina Vacca, Federica Govoni, Matteo Murgia, Richard A. Perley, Luigina Feretti, Gabriele Giovannini, Ettore Carretti, Fabio Gastaldello, Filippo Cova, Paolo Marchegiani, Elia Battistelli, Walter Boschin, Torsten A. Ensslin, Marisa Girardi, Francesca Loi, Federico Radiconi

Abstract: Large-scale magnetic fields reveal themselves through diffuse synchrotron sources observed in galaxy clusters such as radio halos. Total intensity filaments of these sources have been observed in polarization as well, but only in three radio halos out of about one hundred currently known. In this paper we analyze new polarimetric Very Large Array data of the diffuse emission in the galaxy cluster… ▽ More Large-scale magnetic fields reveal themselves through diffuse synchrotron sources observed in galaxy clusters such as radio halos. Total intensity filaments of these sources have been observed in polarization as well, but only in three radio halos out of about one hundred currently known. In this paper we analyze new polarimetric Very Large Array data of the diffuse emission in the galaxy cluster Abell 523 in the frequency range 1-2 GHz. We find for the first time evidence of polarized emission on scales of ~ 2.5 Mpc. Total intensity emission is observed only in the central part of the source, likely due to observational limitations. To look for total intensity emission beyond the central region, we combine these data with single-dish observations from the Sardinia Radio Telescope and we compare them with multi-frequency total intensity observations obtained with different instruments, including the LOw Frequency ARray and the Murchison Widefield Array. By analysing the rotation measure properties of the system and utilizing numerical simulations, we infer that this polarized emission is associated with filaments of the radio halo located in the outskirts of the system, in the peripheral region closest to the observer. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 14 pages, 12 figures, accepted for publication on MNRAS

The Galactic 3D large-scale dust distribution via Gaussian process regression on spherical coordinates

Authors: R. H. Leike, G. Edenhofer, J. Knollmüller, C. Alig, P. Frank, T. A. Enßlin

Abstract: Knowing the Galactic 3D dust distribution is relevant for understanding many processes in the interstellar medium and for correcting many astronomical observations for dust absorption and emission. Here, we aim for a 3D reconstruction of the Galactic dust distribution with an increase in the number of meaningful resolution elements by orders of magnitude with respect to previous reconstructions, w… ▽ More Knowing the Galactic 3D dust distribution is relevant for understanding many processes in the interstellar medium and for correcting many astronomical observations for dust absorption and emission. Here, we aim for a 3D reconstruction of the Galactic dust distribution with an increase in the number of meaningful resolution elements by orders of magnitude with respect to previous reconstructions, while taking advantage of the dust's spatial correlations to inform the dust map. We use iterative grid refinement to define a log-normal process in spherical coordinates. This log-normal process assumes a fixed correlation structure, which was inferred in an earlier reconstruction of Galactic dust. Our map is informed through 111 Million data points, combining data of PANSTARRS, 2MASS, Gaia DR2 and ALLWISE. The log-normal process is discretized to 122 Billion degrees of freedom, a factor of 400 more than our previous map. We derive the most probable posterior map and an uncertainty estimate using natural gradient descent and the Fisher-Laplace approximation. The dust reconstruction covers a quarter of the volume of our Galaxy, with a maximum coordinate distance of $16\,\text{kpc}$, and meaningful information can be found up to at distances of $4\,$kpc, still improving upon our earlier map by a factor of 5 in maximal distance, of $900$ in volume, and of about eighteen in angular grid resolution. Unfortunately, the maximum posterior approach chosen to make the reconstruction computational affordable introduces artifacts and reduces the accuracy of our uncertainty estimate. Despite of the apparent limitations of the presented 3D dust map, a good part of the reconstructed structures are confirmed by independent maser observations. Thus, the map is a step towards reliable 3D Galactic cartography and already can serve for a number of tasks, if used with care. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Multi-Component Imaging of the Fermi Gamma-ray Sky in the Spatio-spectral Domain

Authors: Lukas I. Platz, Jakob Knollmüller, Philipp Arras, Philipp Frank, Martin Reinecke, Dominik Jüstel, Torsten A. Enßlin

Abstract: We perform two distinct spatio-spectral reconstructions of the gamma-ray sky in the range of 0.56-316 GeV based on Fermi Large Area Telescope (LAT) data. Both describe the sky brightness to be composed of a diffuse-emission and a point-source component. The first model requires minimal assumptions and provides a template-free reconstruction as a reference. It makes use of spatial and spectral corr… ▽ More We perform two distinct spatio-spectral reconstructions of the gamma-ray sky in the range of 0.56-316 GeV based on Fermi Large Area Telescope (LAT) data. Both describe the sky brightness to be composed of a diffuse-emission and a point-source component. The first model requires minimal assumptions and provides a template-free reconstruction as a reference. It makes use of spatial and spectral correlations to distinguish between the different components. The second model is physics-informed and further differentiates between diffuse emission of hadronic and leptonic origin. For this, we assume parametric, but spatially varying energy spectra to distinguish between the processes and use thermal Galactic dust observations to indicate the preferred sites of hadronic interactions. To account for instrumental effects we model the point-spread, the energy dispersion, and the exposure of the telescope throughout the observation. The reconstruction problem is formulated as a Bayesian inference task, that is solved by variational inference. We show decompositions of the Gamma-ray flux into diffuse and point-like emissions, and of the diffuse emissions into multiple physically motivated components. The diffuse decomposition provides an unprecedented view of the Galactic leptonic diffuse emission. It shows the Fermi bubbles and their spectral variations in high fidelity and other areas exhibiting strong cosmic ray electron contents, such as a thick disk in the inner Galaxy and outflow regions. Furthermore, we report a hard spectrum gamma ray arc in the northern outer bubble co-spatial with the reported X-ray arc by the eROSITA collaboration. All our spatio-spectral sky reconstructions and their uncertainty quantification are publicly available. △ Less

Submitted 6 December, 2023; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: Please see the journal version for a language-edited release of this manuscript (open-access): https://doi.org/10.1051/0004-6361/202243819

Journal ref: A&A 680, A2 (2023)

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for th… ▽ More In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

Geometric variational inference

Authors: Philipp Frank, Reimar Leike, Torsten A. Enßlin

Abstract: Efficiently accessing the information contained in non-linear and high dimensional probability distributions remains a core challenge in modern statistics. Traditionally, estimators that go beyond point estimates are either categorized as Variational Inference (VI) or Markov-Chain Monte-Carlo (MCMC) techniques. While MCMC methods that utilize the geometric properties of continuous probability dist… ▽ More Efficiently accessing the information contained in non-linear and high dimensional probability distributions remains a core challenge in modern statistics. Traditionally, estimators that go beyond point estimates are either categorized as Variational Inference (VI) or Markov-Chain Monte-Carlo (MCMC) techniques. While MCMC methods that utilize the geometric properties of continuous probability distributions to increase their efficiency have been proposed, VI methods rarely use the geometry. This work aims to fill this gap and proposes geometric Variational Inference (geoVI), a method based on Riemannian geometry and the Fisher information metric. It is used to construct a coordinate transformation that relates the Riemannian manifold associated with the metric to Euclidean space. The distribution, expressed in the coordinate system induced by the transformation, takes a particularly simple form that allows for an accurate variational approximation by a normal distribution. Furthermore, the algorithmic structure allows for an efficient implementation of geoVI which is demonstrated on multiple examples, ranging from low-dimensional illustrative ones to non-linear, hierarchical Bayesian inverse problems in thousands of dimensions. △ Less

Submitted 2 July, 2021; v1 submitted 21 May, 2021; originally announced May 2021.

Comments: 42 pages, 18 figures, accepted by Entropy

The Galactic Faraday rotation sky 2020

Authors: Sebastian Hutschenreuter, Craig S. Anderson, Sarah Betti, Geoffrey C. Bower, Jo-Anne Brown, Marcus Brüggen, Ettore Carretti, Tracy Clarke, Andrew Clegg, Allison Costa, Steve Croft, Cameron Van Eck, B. M. Gaensler, Francesco de Gasperin, Marijke Haverkorn, George Heald, Charles L. H. Hull, Makoto Inoue, Melanie Johnston-Hollitt, Jane Kaczmarek, Casey Law, Yik Ki Ma, David MacMahon, Sui Ann Mao, Christopher Riseley , et al. (12 additional authors not shown)

Abstract: This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This h… ▽ More This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This has culminated in an all-sky data set of 55,190 data points, which is a significant expansion on the 41,330 used in previous works, hence making an updated separation of the Galactic component a promising venture. The increased source density allows us to present our results in a resolution of about $1.3\cdot 10^{-2}\, \mathrm{deg}^2$ ($46.8\,\mathrm{arcmin}^2$), which is a twofold increase compared to previous works. As for previous Faraday rotation sky reconstructions, this work is based on information field theory, a Bayesian inference scheme for field-like quantities which handles noisy and incomplete data. In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand $\mathrm{rad}\,\mathrm{m}^{-2}$. The improvements can mainly be attributed to the new catalog of Faraday data, but are also supported by advances in correlation structure modeling within numerical information field theory. We furthermore give a detailed discussion on statistical properties of the Faraday rotation sky and investigate correlations to other data sets. △ Less

Submitted 8 October, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: accepted in A&A; 15 pages, 12 Figures; results at https://wwwmpa.mpa-garching.mpg.de/~ensslin/research/data/faraday2020.html and http://cutouts.cirada.ca/rmcutout

Journal ref: A&A 657, A43 (2022)

Reconstructing non-repeating radio pulses with Information Field Theory

Authors: Christoph Welling, Philipp Frank, Torsten A. Enßlin, Anna Nelles

Abstract: Particle showers in dielectric media produce radio signals which are used for the detection of both ultra-high energy cosmic rays and neutrinos with energies above a few PeV. The amplitude, polarization, and spectrum of these short, broadband radio pulses allow us to draw conclusions about the primary particles that caused them, as well as the mechanics of shower development and radio emission. Ho… ▽ More Particle showers in dielectric media produce radio signals which are used for the detection of both ultra-high energy cosmic rays and neutrinos with energies above a few PeV. The amplitude, polarization, and spectrum of these short, broadband radio pulses allow us to draw conclusions about the primary particles that caused them, as well as the mechanics of shower development and radio emission. However, confidently reconstructing the radio signals can pose a challenge, as they are often obscured by background noise. Information Field Theory offers a robust approach to this challenge by using Bayesian inference to calculate the most likely radio signal, given the recorded data. In this paper, we describe the application of Information Field Theory to radio signals from particle showers in both air and ice and demonstrate how accurately pulse parameters can be obtained from noisy data. △ Less

Submitted 17 March, 2021; v1 submitted 30 January, 2021; originally announced February 2021.

Fast magnetic field amplification in distant galaxyclusters

Authors: Gabriella Di Gennaro, Reinout J. van Weeren, Gianfranco Brunetti, Rossella Cassano, Marcus Brüggen, Matthias Hoeft, Timothy W. Shimwell, Huub J. A. Röttgering, Annalisa Bonafede, Andrea Botteon, Virginia Cuciti, Daniele Dallacasa, Francesco de Gasperin, Paola Domínguez-Fernández, Torsten A. Ensslin, Fabio Gastaldello, Soumyajit Mandal, Mariachiara Rossetti, Aurora Simionescu

Abstract: In the present-day Universe, magnetic fields pervade galaxy clusters, with strengths of a few microGauss obtained from Faraday Rotation. Evidence for cluster magnetic fields is also provided by Megaparsec-scale radio emission, namely radio halos and relics. These are commonly found in merging systems and are characterized by a steep radio spectrum. It is widely believed that magneto-hydrodynamical… ▽ More In the present-day Universe, magnetic fields pervade galaxy clusters, with strengths of a few microGauss obtained from Faraday Rotation. Evidence for cluster magnetic fields is also provided by Megaparsec-scale radio emission, namely radio halos and relics. These are commonly found in merging systems and are characterized by a steep radio spectrum. It is widely believed that magneto-hydrodynamical turbulence and shock-waves (re-)accelerate cosmic rays, producing halos and relics. The origin and the amplification of magnetic fields in clusters is not well understood. It has been proposed that turbulence drives a small-scaledynamo that amplifies seed magnetic fields (primordial and/or injected by galactic outflows, as active galactic nuclei, starbursts, or winds). At high redshift, radio halos are expected to be faint, due to the Inverse Compton losses and dimming effect with distance. Moreover, Faraday Rotation measurements are difficult to obtain. If detected, distant radio halosprovide an alternative tool to investigate magnetic field amplification. Here, we report LOFAR observations which reveal diffuse radio emission in massive clusters when the Universe was only half of its present age, with a sample occurrence fraction of about 50%. The high radio luminosities indicate that these clusters have similar magnetic field strengths to those in nearby clusters, and suggest that magnetic field amplification is fast during the first phases ofcluster formation. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: Published in Nature Astronomy on 2 November 2020. The published version is available at this URL https://www.nature.com/articles/s41550-020-01244-5#citeas

Towards Bayesian Data Compression

Authors: Johannes Harth-Kitzerow, Reimar Leike, Philipp Arras, Torsten A. Enßlin

Abstract: In order to handle large data sets omnipresent in modern science, efficient compression algorithms are necessary. Here, a Bayesian data compression (BDC) algorithm that adapts to the specific measurement situation is derived in the context of signal reconstruction. BDC compresses a data set under conservation of its posterior structure with minimal information loss given the prior knowledge on the… ▽ More In order to handle large data sets omnipresent in modern science, efficient compression algorithms are necessary. Here, a Bayesian data compression (BDC) algorithm that adapts to the specific measurement situation is derived in the context of signal reconstruction. BDC compresses a data set under conservation of its posterior structure with minimal information loss given the prior knowledge on the signal, the quantity of interest. Its basic form is valid for Gaussian priors and likelihoods. For constant noise standard deviation, basic BDC becomes equivalent to a Bayesian analog of principal component analysis. Using Metric Gaussian Variational Inference, BDC generalizes to non-linear settings. In its current form, BDC requires the storage of effective instrument response functions for the compressed data and corresponding noise encoding the posterior covariance structure. Their memory demand counteract the compression gain. In order to improve this, sparsity of the compressed responses can be obtained by separating the data into patches and compressing them separately. The applicability of BDC is demonstrated by applying it to synthetic data and radio astronomical data. Still the algorithm needs further improvement as the computation time of the compression and subsequent inference exceeds the time of the inference with the original data. △ Less

Submitted 29 December, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 39 pages, 15 figures, 1 table, for code, see https://gitlab.mpcdf.mpg.de/jharthki/bdc

Journal ref: ANNALEN DER PHYSIK 2021, 2000508

Efficient wide-field radio interferometry response

Authors: Philipp Arras, Martin Reinecke, Rüdiger Westermann, Torsten A. Enßlin

Abstract: Radio interferometers do not measure the sky brightness distribution directly but rather a modified Fourier transform of it. Imaging algorithms, thus, need a computational representation of the linear measurement operator and its adjoint, irrespective of the specific chosen imaging algorithm. In this paper, we present a C++ implementation of the radio interferometric measurement operator for wide-… ▽ More Radio interferometers do not measure the sky brightness distribution directly but rather a modified Fourier transform of it. Imaging algorithms, thus, need a computational representation of the linear measurement operator and its adjoint, irrespective of the specific chosen imaging algorithm. In this paper, we present a C++ implementation of the radio interferometric measurement operator for wide-field measurements which is based on "improved $w$-stacking". It can provide high accuracy (down to $\approx 10^{-12}$), is based on a new gridding kernel which allows smaller kernel support for given accuracy, dynamically chooses kernel, kernel support and oversampling factor for maximum performance, uses piece-wise polynomial approximation for cheap evaluations of the gridding kernel, treats the visibilities in cache-friendly order, uses explicit vectorisation if available and comes with a parallelisation scheme which scales well also in the adjoint direction (which is a problem for many previous implementations). The implementation has a small memory footprint in the sense that temporary internal data structures are much smaller than the respective input and output data, allowing in-memory processing of data sets which needed to be read from disk or distributed across several compute nodes before. △ Less

Submitted 28 January, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 13 pages, 8 figures

Journal ref: A&A 646, A58 (2021)

Probabilistic simulation of partial differential equations

Authors: Philipp Frank, Torsten A. Enßlin

Abstract: Computer simulations of differential equations require a time discretization, which inhibits to identify the exact solution with certainty. Probabilistic simulations take this into account via uncertainty quantification. The construction of a probabilistic simulation scheme can be regarded as Bayesian filtering by means of probabilistic numerics. Gaussian prior based filters, specifically Gauss-Ma… ▽ More Computer simulations of differential equations require a time discretization, which inhibits to identify the exact solution with certainty. Probabilistic simulations take this into account via uncertainty quantification. The construction of a probabilistic simulation scheme can be regarded as Bayesian filtering by means of probabilistic numerics. Gaussian prior based filters, specifically Gauss-Markov priors, have successfully been applied to simulation of ordinary differential equations (ODEs) and give rise to filtering problems that can be solved efficiently. This work extends this approach to partial differential equations (PDEs) subject to periodic boundary conditions and utilizes continuous Gaussian processes in space and time to arrive at a Bayesian filtering problem structurally similar to the ODE setting. The usage of a process that is Markov in time and statistically homogeneous in space leads to a probabilistic spectral simulation method that allows for an efficient realization. Furthermore, the Bayesian perspective allows the incorporation of methods developed within the context of information field theory such as the estimation of the power spectrum associated with the prior distribution, to be jointly estimated along with the solution of the PDE. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 13 pages, 4 figures, submitted to Physical Review E

Bayesian decomposition of the Galactic multi-frequency sky using probabilistic autoencoders

Authors: Sara Milosevic, Philipp Frank, Reimar H. Leike, Ancla Müller, Torsten A. Enßlin

Abstract: All-sky observations of the Milky Way show both Galactic and non-Galactic diffuse emission, for example from interstellar matter or the cosmic microwave background (CMB). The different emitters are partly superimposed in the measurements, partly they obscure each other, and sometimes they dominate within a certain spectral range. The decomposition of the underlying radiative components from spectr… ▽ More All-sky observations of the Milky Way show both Galactic and non-Galactic diffuse emission, for example from interstellar matter or the cosmic microwave background (CMB). The different emitters are partly superimposed in the measurements, partly they obscure each other, and sometimes they dominate within a certain spectral range. The decomposition of the underlying radiative components from spectral data is a signal reconstruction problem and often associated with detailed physical modeling and substantial computational effort. We aim to build an effective and self-instructing algorithm detecting the essential spectral information contained Galactic all-sky data covering spectral bands from $γ$-ray to radio waves. Utilizing principles from information theory, we develop a state-of-the-art variational autoencoder specialized on the adaption to Gaussian noise statistics. We first derive a generic generative process that leads from a low-dimensional set of emission features to the observed high-dimensional data. We formulate a posterior distribution of these features using Bayesian methods and approximate this posterior with variational inference. The algorithm efficiently encodes the information of 35 Galactic emission data sets in ten latent feature maps. These contain the essential information required to reconstruct the initial data with high fidelity and are ranked by the algorithm according to their significance for data regeneration. The three most significant feature maps encode astrophysical components: (1) The dense interstellar medium (ISM), (2) the hot and dilute regions of the ISM and (3) the CMB. The machine-assisted and data-driven dimensionality reduction of spectral data is able to uncover the physical features encoding the input data. Our algorithm is able to extract the dense and dilute Galactic regions, as well as the CMB, from the sky brightness values only. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 25 pages, 8 figures, 3 tables. Submitted to Astronomy & Astrophysics

Journal ref: A&A 650, A100 (2021)

Planck intermediate results. LV. Reliability and thermal properties of high-frequency sources in the Second Planck Catalogue of Compact Sources

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, C. Burigana, E. Calabrese, P. Carvalho, H. C. Chiang, B. P. Crill, F. Cuttaia, A. de Rosa, G. de Zotti , et al. (95 additional authors not shown)

Abstract: We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution, which leads to emission at both lower and higher frequencies, and adjusts the parameters of the sour… ▽ More We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution, which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1 per cent of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. The results of the BeeP extension of PCCS2, which are made publicly available via the PLA, will enable the study of the thermal properties of well-defined samples of compact Galactic and extra-galactic dusty sources. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 55 pages. Accepted for publication in A&A. The BeeP catalogue will be published in the Planck Legacy Archive (https://pla.esac.esa.int/pla)

Journal ref: A&A 644, A99 (2020)

Comparison of classical and Bayesian imaging in radio interferometry

Authors: Philipp Arras, Hertzog L. Bester, Richard A. Perley, Reimar Leike, Oleg Smirnov, Rüdiger Westermann, Torsten A. Enßlin

Abstract: CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions… ▽ More CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions; and its results are highly dependent on the so-called weighting scheme as well as on any human choice of CLEAN masks to guiding the imaging. Here, we present the Bayesian imaging algorithm resolve which solves the above problems and naturally leads to super-resolution. We take a VLA observation of Cygnus~A at four different frequencies and image it with single-scale CLEAN, multi-scale CLEAN and resolve. Alongside the sky brightness distribution resolve estimates a baseline-dependent correction function for the noise budget, the Bayesian equivalent of weighting schemes. We report noise correction factors between 0.4 and 429. The enhancements achieved by resolve come at the cost of higher computational effort. △ Less

Submitted 25 January, 2021; v1 submitted 26 August, 2020; originally announced August 2020.

Comments: 23 pages, 16 figures, 4 tables, data published at https://doi.org/10.5281/zenodo.4267057

Journal ref: A&A 646, A84 (2021)

Resolving nearby dust clouds

Authors: R. H. Leike, M. Glatzle, T. A. Enßlin

Abstract: Aims: Mapping the interstellar medium in 3D provides a wealth of insights into its inner working. The Milky Way is the only galaxy for which detailed 3D mapping can be achieved in principle. In this paper, we reconstruct the dust density in and around the local super-bubble. Methods: The combined data from surveys such as Gaia, 2MASS, PANSTARRS, and ALLWISE provide the necessary information to m… ▽ More Aims: Mapping the interstellar medium in 3D provides a wealth of insights into its inner working. The Milky Way is the only galaxy for which detailed 3D mapping can be achieved in principle. In this paper, we reconstruct the dust density in and around the local super-bubble. Methods: The combined data from surveys such as Gaia, 2MASS, PANSTARRS, and ALLWISE provide the necessary information to make detailed maps of the interstellar medium in our surrounding. To this end, we used variational inference and Gaussian processes to model the dust extinction density, exploiting its intrinsic correlations. Results: We reconstructed a highly resolved dust map, showing the nearest dust clouds at a distance of up to 400pc with a resolution of 1pc. Conclusions: Our reconstruction provides insights into the structure of the interstellar medium. We compute summary statistics of the spectral index and the 1-point function of the logarithmic dust extinction density, which may constrain simulations of the interstellar medium that achieve a similar resolution. △ Less

Submitted 5 August, 2020; v1 submitted 14 April, 2020; originally announced April 2020.

Journal ref: A&A 639, A138 (2020)

Planck intermediate results. LVI. Detection of the CMB dipole through modulation of the thermal Sunyaev-Zeldovich effect: Eppur si muove II

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, C. Burigana, E. Calabrese, J. -F. Cardoso, B. Casaponsa, H. C. Chiang, C. Combet, D. Contreras, B. P. Crill , et al. (104 additional authors not shown)

Abstract: The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular… ▽ More The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first, and significant, detection of this signal in the tSZ maps and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, which is modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales. △ Less

Submitted 7 September, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 15 pages, 8 figures. Added references, small clarifying and language edits. All results remain the same

Journal ref: A&A 644, A100 (2020)

Probing Cosmic Ray Transport with Radio Synchrotron Harps in the Galactic Center

Authors: Timon Thomas, Christoph Pfrommer, Torsten A. Enßlin

Abstract: Recent observations with the MeerKAT radio telescope reveal a unique population of faint non-thermal filaments pervading the central molecular zone (CMZ). Some of those filaments are organized into groups of almost parallel filaments, seemingly sorted by their length, so that their morphology resembles a harp with radio emitting "strings". We argue that the synchrotron emitting GeV electrons of th… ▽ More Recent observations with the MeerKAT radio telescope reveal a unique population of faint non-thermal filaments pervading the central molecular zone (CMZ). Some of those filaments are organized into groups of almost parallel filaments, seemingly sorted by their length, so that their morphology resembles a harp with radio emitting "strings". We argue that the synchrotron emitting GeV electrons of these radio harps have been consecutively injected by the same source (a massive star or pulsar) into spatially intermittent magnetic fiber bundles within a magnetic flux tube or via time-dependent injection events. After escaping from this source, the propagation of cosmic ray (CR) electrons inside a flux tube is governed by the theory of CR transport. We propose to use observations of radio harp filaments to gain insight into the specifics of CR propagation along magnetic fields of which there are two principle modes: CRs could either stream with self-excited magneto-hydrodynamical waves or diffuse along the magnetic field. To disentangle these possibilities, we conduct hydrodynamical simulations of either purely diffusing or streaming CR electrons and compare the resulting brightness distributions to the observed synchrotron profiles of the radio harps. We find compelling evidence that CR streaming is the dominant propagation mode for GeV CRs in one of the radio harps. Observations at higher angular resolution should detect more radio harps and may help to disentangle projection effects of the possibly three-dimensional flux-tube structure of the other radio harps. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 7 pages, 5 figures, submitted to ApJL. Comments welcome!

Jets, Bubbles, and Heat Pumps in Galaxy Clusters

Authors: Yi-Hao Chen, Sebastian Heinz, Torsten A. Enßlin

Abstract: Feedback from AGN jets has been proposed to counteract the catastrophic cooling in many galaxy clusters. However, it is still unclear which physical processes are acting to couple the energy from the bi-directional jets to the ICM. We study the long-term evolution of rising bubbles that were inflated by AGN jets using MHD simulations. In the wake of the rising bubbles, a significant amount of low-… ▽ More Feedback from AGN jets has been proposed to counteract the catastrophic cooling in many galaxy clusters. However, it is still unclear which physical processes are acting to couple the energy from the bi-directional jets to the ICM. We study the long-term evolution of rising bubbles that were inflated by AGN jets using MHD simulations. In the wake of the rising bubbles, a significant amount of low-entropy gas is brought into contact with the hot cluster gas. We assess the energy budget of the uplifted gas and find it comparable to the total energy injected by the jets. Although our simulation does not include explicit thermal conduction, we find that, for reasonable assumptions about the conduction coefficient, the rate is fast enough that much of the uplifted gas may be thermalized before it sinks back to the core. Thus, we propose that the AGN can act like a heat pump to move low-entropy gas from the cluster core to the heat reservoir and will be able to heat the inner cluster more efficiently than would be possible by direct energy transfer from jets alone. We show that the maximum efficiency of this mechanism, i.e. the ratio between the conductive thermal energy and the work needed to lift the gas, $ξ_{\mathrm{max}}$ can exceed 100 per cent. While $ξ$ < $ξ_{\mathrm{max}}$ in realistic scenarios, AGN-induced thermal conduction has the potential to significantly increase the efficiency with which AGN can heat cool-core clusters and transform the bursty AGN activities into a smoother and enduring heating process. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: 11 pages, 10 figures. Accepted for publication in MNRAS

Planck 2018 results. V. CMB power spectra and likelihoods

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (143 additional authors not shown)

Abstract: This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low… ▽ More This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $τ$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $τ$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $Λ$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$σ$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged) △ Less

Submitted 15 September, 2020; v1 submitted 30 July, 2019; originally announced July 2019.

Comments: Revised to match version published in Astronomy & Astrophysics

Journal ref: A&A 641, A5 (2020)

hammurabi X: Simulating Galactic Synchrotron Emission with Random Magnetic Fields

Authors: Jiaxin Wang, Tess R. Jaffe, Torsten A. Enßlin, Piero Ullio, Shamik Ghosh, Larissa Santos

Abstract: 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… ▽ More 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. △ Less

Submitted 10 February, 2020; v1 submitted 29 June, 2019; originally announced July 2019.

Comments: 25 pages, 18 figures, Accepted by ApJS with joint JOSS publication at https://doi.org/10.21105/joss.01889, the numerical package can be found at https://bitbucket.org/hammurabicode/hamx or https://github.com/hammurabi-dev/hammurabiX

Journal ref: ApJS (2020) 247 18

A radio ridge connecting two galaxy clusters in a filament of the cosmic web

Authors: F. Govoni, E. Orrù, A. Bonafede, M. Iacobelli, R. Paladino, F. Vazza, M. Murgia, V. Vacca, G. Giovannini, L. Feretti, F. Loi, G. Bernardi, C. Ferrari, R. F. Pizzo, C. Gheller, S. Manti, M. Brüggen, G. Brunetti, R. Cassano, F. de Gasperin, T. A. Enßlin, M. Hoeft, C. Horellou, H. Junklewitz, H. J. A. Röttgering , et al. (4 additional authors not shown)

Abstract: Galaxy clusters are the most massive gravitationally bound structures in the Universe. They grow by accreting smaller structures in a merging process that produces shocks and turbulence in the intra-cluster gas. We observed a ridge of radio emission connecting the merging galaxy clusters Abell 0399 and Abell 0401 with the Low Frequency Array (LOFAR) at 140 MHz. This emission requires a population… ▽ More Galaxy clusters are the most massive gravitationally bound structures in the Universe. They grow by accreting smaller structures in a merging process that produces shocks and turbulence in the intra-cluster gas. We observed a ridge of radio emission connecting the merging galaxy clusters Abell 0399 and Abell 0401 with the Low Frequency Array (LOFAR) at 140 MHz. This emission requires a population of relativistic electrons and a magnetic field located in a filament between the two galaxy clusters. We performed simulations to show that a volume-filling distribution of weak shocks may re-accelerate a pre-existing population of relativistic particles, producing emission at radio wavelengths that illuminates the magnetic ridge. △ Less

Submitted 18 June, 2019; originally announced June 2019.

Comments: Published in Science

Journal ref: Science 2019, Volume 364, Issue 6444, pp. 981-984

Planck 2018 results. VII. Isotropy and Statistics of the CMB

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso, B. Casaponsa, H. C. Chiang , et al. (125 additional authors not shown)

Abstract: Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $Λ$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on la… ▽ More Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $Λ$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $Λ$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales. △ Less

Submitted 14 September, 2020; v1 submitted 6 June, 2019; originally announced June 2019.

Comments: Paper VII of the Planck 2018 release, revised to closely match version published in Astronomy and Astrophysics

Journal ref: A&A 641, A7 (2020)

Planck 2018 results. IX. Constraints on primordial non-Gaussianity

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso, B. Casaponsa, A. Challinor , et al. (135 additional authors not shown)

Abstract: We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polariz… ▽ More We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged). △ Less

Submitted 14 May, 2019; originally announced May 2019.

Comments: 50 pages, 20 figures

The Galactic Faraday depth sky revisited

Authors: Sebastian Hutschenreuter, Torsten A. Enßlin

Abstract: The Galactic Faraday depth sky is a tracer for both the Galactic magnetic field and the thermal electron distribution. It has been previously reconstructed from polarimetric measurements of extra-galactic point sources. Here, we improve on these works by using an updated inference algorithm as well as by taking into account the free-free emission measure map from the Planck survey. In the future,… ▽ More The Galactic Faraday depth sky is a tracer for both the Galactic magnetic field and the thermal electron distribution. It has been previously reconstructed from polarimetric measurements of extra-galactic point sources. Here, we improve on these works by using an updated inference algorithm as well as by taking into account the free-free emission measure map from the Planck survey. In the future, the data situation will improve drastically with the next generation Faraday rotation measurements from SKA and its pathfinders. Anticipating this, the aim of this paper is to update the map reconstruction method with the latest development in imaging based on information field theory. We demonstrate the validity of the new algorithm by applying it to the Oppermann et al. (2012) data compilation and compare our results to the previous map.\\ Despite using exactly the previous data set, a number of novel findings are made: A non-parametric reconstruction of an overall amplitude field resembles the free-free emission measure map of the Galaxy. Folding this free-free map into the analysis allows for more detailed predictions. The joint inference enables us to identify regions with deviations from the assumed correlations between the free-free and Faraday data, thereby pointing us to Galactic structures with distinguishably different physics. We e.g. find evidence for an alignment of the magnetic field within the line of sights along both directions of the Orion arm. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 16 pages, 15 figures

Journal ref: A&A 633, A150 (2020)

Metric Gaussian Variational Inference

Authors: Jakob Knollmüller, Torsten A. Enßlin

Abstract: Solving Bayesian inference problems approximately with variational approaches can provide fast and accurate results. Capturing correlation within the approximation requires an explicit parametrization. This intrinsically limits this approach to either moderately dimensional problems, or requiring the strongly simplifying mean-field approach. We propose Metric Gaussian Variational Inference (MGVI)… ▽ More Solving Bayesian inference problems approximately with variational approaches can provide fast and accurate results. Capturing correlation within the approximation requires an explicit parametrization. This intrinsically limits this approach to either moderately dimensional problems, or requiring the strongly simplifying mean-field approach. We propose Metric Gaussian Variational Inference (MGVI) as a method that goes beyond mean-field. Here correlations between all model parameters are taken into account, while still scaling linearly in computational time and memory. With this method we achieve higher accuracy and in many cases a significant speedup compared to traditional methods. MGVI is an iterative method that performs a series of Gaussian approximations to the posterior. We alternate between approximating the covariance with the inverse Fisher information metric evaluated at an intermediate mean estimate and optimizing the KL-divergence for the given covariance with respect to the mean. This procedure is iterated until the uncertainty estimate is self-consistent with the mean parameter. We achieve linear scaling by avoiding to store the covariance explicitly at any time. Instead we draw samples from the approximating distribution relying on an implicit representation and numerical schemes to approximately solve linear equations. Those samples are used to approximate the KL-divergence and its gradient. The usage of natural gradient descent allows for rapid convergence. Formulating the Bayesian model in standardized coordinates makes MGVI applicable to any inference problem with continuous parameters. We demonstrate the high accuracy of MGVI by comparing it to HMC and its fast convergence relative to other established methods in several examples. We investigate real-data applications, as well as synthetic examples of varying size and complexity and up to a million model parameters. △ Less

Submitted 30 January, 2020; v1 submitted 30 January, 2019; originally announced January 2019.

Comments: Code is part of NIFTy5 release at https://gitlab.mpcdf.mpg.de/ift/NIFTy

Charting nearby dust clouds using Gaia data only

Authors: R. H. Leike, T. A. Enßlin

Abstract: Aims: Highly resolved maps of the local Galactic dust are an important ingredient for sky emission models. In nearly the whole electromagnetic spectrum one can see imprints of dust, many of which originate from dust clouds within 300pc. Having a detailed 3D reconstruction of these local dust clouds enables detailed studies, helps to quantify the impact on other observables and is a milestone neces… ▽ More Aims: Highly resolved maps of the local Galactic dust are an important ingredient for sky emission models. In nearly the whole electromagnetic spectrum one can see imprints of dust, many of which originate from dust clouds within 300pc. Having a detailed 3D reconstruction of these local dust clouds enables detailed studies, helps to quantify the impact on other observables and is a milestone necessary to enable larger reconstructions, as every sightline for more distant objects will pass through the local dust. Methods: To infer the dust density we use parallax and absorption estimates published by the Gaia collaboration in their second data release. We model the dust as a log-normal process using a hierarchical Bayesian model. We also infer non-parametrically the kernel of the log-normal process, which corresponds to the physical spatial correlation power spectrum of the log-density. Results: Using only Gaia data of the second Gaia data release, we reconstruct the 3D dust density and its spatial correlation spectrum in a 600pc cube centered on the Sun. We report a spectral index of the logarithmic dust density of $3.1$ on Fourier scales with wavelengths between 2pc and 125pc. The resulting 3D dust map as well as the power spectrum and posterior samples are publicly available for download. △ Less

Submitted 17 July, 2019; v1 submitted 17 January, 2019; originally announced January 2019.

Journal ref: A&A 631, A32 (2019)

Bayesian parameter estimation of miss-specified models

Authors: Johannes Oberpriller, T. A. Enßlin

Abstract: Fitting a simplifying model with several parameters to real data of complex objects is a highly nontrivial task, but enables the possibility to get insights into the objects physics. Here, we present a method to infer the parameters of the model, the model error as well as the statistics of the model error. This method relies on the usage of many data sets in a simultaneous analysis in order to ov… ▽ More Fitting a simplifying model with several parameters to real data of complex objects is a highly nontrivial task, but enables the possibility to get insights into the objects physics. Here, we present a method to infer the parameters of the model, the model error as well as the statistics of the model error. This method relies on the usage of many data sets in a simultaneous analysis in order to overcome the problems caused by the degeneracy between model parameters and model error. Errors in the modeling of the measurement instrument can be absorbed in the model error allowing for applications with complex instruments. △ Less

Submitted 19 December, 2018; originally announced December 2018.

The LOFAR Two-metre Sky Survey - II. First data release

Authors: T. W. Shimwell, C. Tasse, M. J. Hardcastle, A. P. Mechev, W. L. Williams, P. N. Best, H. J. A. Röttgering, J. R. Callingham, T. J. Dijkema, F. de Gasperin, D. N. Hoang, B. Hugo, M. Mirmont, J. B. R. Oonk, I. Prandoni, D. Rafferty, J. Sabater, O. Smirnov, R. J. van Weeren, G. J. White, M. Atemkeng, L. Bester, E. Bonnassieux, M. Brüggen, G. Brunetti , et al. (82 additional authors not shown)

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,μ$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: 16 figures, 1 table and 22 pages. This paper is part of the LOFAR surveys data release 1 and has been accepted for publication in a special edition of A&A that will appear in Feb 2019, volume 622. The catalogues and images from the data release will be publicly available on lofar-surveys.org upon publication of the journal

Determining the composition of radio plasma via circular polarization: the prospects of the Cygnus A hot spots

Authors: Torsten A. Enßlin, Sebastian Hutschenreuter, Gopal-Krishna

Abstract: The composition of the relativistic plasma produced in active galactic nuclei and ejected via powerful jets into the interstellar/intergalactic medium is still a major unsettled issue. It might be a positron-electron plasma in case the plasma was created by pair production in the intense photon fields near accreting super-massive black holes. Alternatively, it might be an electron-proton plasma in… ▽ More The composition of the relativistic plasma produced in active galactic nuclei and ejected via powerful jets into the interstellar/intergalactic medium is still a major unsettled issue. It might be a positron-electron plasma in case the plasma was created by pair production in the intense photon fields near accreting super-massive black holes. Alternatively, it might be an electron-proton plasma in case magnetic fields lift and accelerate the thermal gas of accretion discs into relativistic jets as the recent detection of $γ$-rays from blazars indicates. Despite various attempts to unambiguously establish the composition of the relativistic jets, this remains a major unknown. Here, we propose a way to settle the question via sensitive measurements of circular polarization (CP) in the radio emission of the hot spots of bright radio galaxies like Cygnus A. The CP of synchrotron emission is determined by the circular motions of the radiating relativistic leptons. In case of charge symmetric energy spectra of a electron-positron plasma, it should be exactly zero. In case of an electron-proton plasma the electrons imprint their gyration onto the CP and we expect the hot spots of Cygnus A to exhibit a fractional CP at a level of $10^{-3}\,(ν/\mbox{GHz})^{-{1}/{2}}$, which is challenging to measure, but not completely unfeasible. △ Less

Submitted 10 July, 2019; v1 submitted 21 August, 2018; originally announced August 2018.

Comments: 6 pages, no figures, accepted by JCAP, erratum with respect to accepted version applied, numerical corrections leading to a stronger effect

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. I. R. Alves, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher, C. Burigana, E. Calabrese , et al. (138 additional authors not shown)

Abstract: We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on… ▽ More We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $ψ$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas. △ Less

Submitted 7 March, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 641, A12 (2020)

Planck 2018 results. X. Constraints on inflation

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (151 additional authors not shown)

Abstract: We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index… ▽ More We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (φ) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation. △ Less

Submitted 2 August, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: References added and minor improvements. BICEP2/Keck Array BK15 is used in the place of BICEP2/Keck Array BK14

Journal ref: A&A 641, A10 (2020)

Planck 2018 results. VIII. Gravitational lensing

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron , et al. (133 additional authors not shown)

Abstract: We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,σ$ to $9\,σ$. Combined with temperature, lensing is detected at $40\,σ$. We present an extensive set of tests of the robustness of the lensing-potential… ▽ More We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,σ$ to $9\,σ$. Combined with temperature, lensing is detected at $40\,σ$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{ΛCDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $σ_8 Ω_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,σ$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $σ_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $Ω_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $σ_8$ to better than $1\,\%$ precision, finding $σ_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged). △ Less

Submitted 29 July, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: Abstract abridged for arxiv submission. Lensing data products available at https://wiki.cosmos.esa.int/planck-legacy-archive/index.php/Lensing. Matches version accepted by A&A, with minor updates from v1

Journal ref: A&A 641, A8 (2020)

Planck 2018 results. VI. Cosmological parameters

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, R. Battye, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese , et al. (157 additional authors not shown)

Abstract: We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $Λ$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $Λ$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis g… ▽ More We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $Λ$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $Λ$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $Ω_c h^2 = 0.120\pm 0.001$, baryon density $Ω_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $τ= 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100θ_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$Λ$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $Ω_m = 0.315\pm 0.007$; and matter fluctuation amplitude $σ_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$Λ$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_ν< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$Λ$CDM at over $2\,σ$, which pulls some parameters that affect the lensing amplitude away from the base-$Λ$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged) △ Less

Submitted 9 August, 2021; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 73 pages; Updated with published reionization result corrigendum on p59. Parameter tables and chains available at https://wiki.cosmos.esa.int/planck-legacy-archive/index.php/Cosmological_Parameters

Journal ref: A&A 641, A6 (2020)

Planck 2018 results. IV. Diffuse component separation

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron, B. Casaponsa, A. Challinor, L. P. L. Colombo , et al. (128 additional authors not shown)

Abstract: We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow those described in earlier papers, adopting four methods (Comm… ▽ More We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects. △ Less

Submitted 26 September, 2020; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 74 pages, A&A, 641, A4

Planck 2018 results. III. High Frequency Instrument data processing and frequency maps

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron, A. Challinor , et al. (130 additional authors not shown)

Abstract: This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations… ▽ More This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35μ$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: Accepted for publication on A&A (AA/2018/32909)

Planck 2018 results. II. Low Frequency Instrument data processing

Authors: Planck Collaboration, Y. Akrami, F. Argüeso, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (126 additional authors not shown)

Abstract: We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters.… ▽ More We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters. △ Less

Submitted 11 September, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Journal ref: A&A 641, A2 (2020)

Planck 2018 results. I. Overview and the cosmological legacy of Planck

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, R. Battye, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese , et al. (166 additional authors not shown)

Abstract: The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Plan… ▽ More The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances. △ Less

Submitted 3 December, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 61 pages, 40 figures, matches version accepted by A&A

Journal ref: A&A 641, A1 (2020)

HMCF - Hamiltonian Monte Carlo Sampling for Fields - A Python framework for HMC sampling with NIFTy

Authors: Christoph Lienhard, Torsten A. Enßlin

Abstract: HMCF "Hamiltonian Monte Carlo for Fields" is a software add-on for the NIFTy "Numerical Information Field Theory" framework implementing Hamiltonian Monte Carlo (HMC) sampling in Python. HMCF as well as NIFTy are designed to address inference problems in high-dimensional spatially correlated setups such as image reconstruction. HMCF adds an HMC sampler to NIFTy that automatically adjusts the many… ▽ More HMCF "Hamiltonian Monte Carlo for Fields" is a software add-on for the NIFTy "Numerical Information Field Theory" framework implementing Hamiltonian Monte Carlo (HMC) sampling in Python. HMCF as well as NIFTy are designed to address inference problems in high-dimensional spatially correlated setups such as image reconstruction. HMCF adds an HMC sampler to NIFTy that automatically adjusts the many free parameters steering the HMC sampling machinery. A wide variety of features ensure efficient full-posterior sampling for high-dimensional inference problems. These features include integration step size adjustment, evaluation of the mass matrix, convergence diagnostics, higher order symplectic integration and simultaneous sampling of parameters and hyperparameters in Bayesian hierarchical models. △ Less

Submitted 2 December, 2018; v1 submitted 7 July, 2018; originally announced July 2018.

Comments: 36 pages, 4 figures, available at https://gitlab.mpcdf.mpg.de/ift/HMCF, see also arXiv:1708.01073

Separating diffuse from point-like sources - a Bayesian approach

Authors: Jakob Knollmüller, Philipp Frank, Torsten A. Enßlin

Abstract: We present the starblade algorithm, a method to separate superimposed point sources from auto-correlated, diffuse flux using a Bayesian model. Point sources are assumed to be independent from each other and to follow a power-law brightness distribution. The diffuse emission is described as a non-parametric log-normal model with a priori unknown correlation structure. This model enforces positivity… ▽ More We present the starblade algorithm, a method to separate superimposed point sources from auto-correlated, diffuse flux using a Bayesian model. Point sources are assumed to be independent from each other and to follow a power-law brightness distribution. The diffuse emission is described as a non-parametric log-normal model with a priori unknown correlation structure. This model enforces positivity of the underlying emission and allows for variation in the order of magnitudes. The correlation structure is recovered non-parametrically in addition to the diffuse flux and is used for the separation of the point sources. Additionally many measurement artifacts appear as point-like or quasi-point-like effects, not compatible with superimposed diffuse emission. An estimate of the separation uncertainty can be provided as well. We demonstrate the capabilities of the derived method on synthetic data and data obtained by the Hubble Space Telescope, emphasizing its effect on instrumental artifacts as well as physical sources. The performance of this method is compared to the background estimation of the SExtractor method, as well as to a denoising auto-encoder. △ Less

Submitted 6 August, 2018; v1 submitted 16 April, 2018; originally announced April 2018.

Information theory for fields

Authors: Torsten A. Enßlin

Abstract: A physical field has an infinite number of degrees of freedom since it has a field value at each location of a continuous space. Therefore, it is impossible to know a field from finite measurements alone and prior information on the field is essential for field inference. An information theory for fields is needed to join the measurement and prior information into probabilistic statements on field… ▽ More A physical field has an infinite number of degrees of freedom since it has a field value at each location of a continuous space. Therefore, it is impossible to know a field from finite measurements alone and prior information on the field is essential for field inference. An information theory for fields is needed to join the measurement and prior information into probabilistic statements on field configurations. Such an information field theory (IFT) is built upon the language of mathematical physics, in particular on field theory and statistical mechanics. IFT permits the mathematical derivation of optimal imaging algorithms, data analysis methods, and even computer simulation schemes. The application of IFT algorithms to astronomical datasets provides high fidelity images of the Universe and facilitates the search for subtle statistical signals from the Big Bang. The concepts of IFT might even pave the road to novel computer simulations that are aware of their own uncertainties. △ Less

Submitted 1 September, 2018; v1 submitted 10 April, 2018; originally announced April 2018.

Comments: 21 pages, 7 figures, resubmitted

Journal ref: Annalen der Physik 2019, vol. 531, issue 3, p. 1800127