Search

Cognition, the process of information processing in form of inference, communication, and memorization, is the central activity of any intelligence. Its physical realization in a brain, computer, or in any other intelligent system requires resources like time, energy, memory, bandwidth, money, and others. Due to limited resources, many real world intelligent systems perform only imperfect cognition. For understanding the trade-off between accuracy and resource investments in existing systems, e.g. in biology, as well as for the resource-aware optimal design of information processing systems, like computer algorithms and artificial neural networks, a quantification of information obtained in an imperfect cognitive operation is desirable. To this end, we propose the concept of achieved information gain (AIG) of a belief update, which is given by the amount of information obtained by updating from the initial knowledge state to the ideal one, minus the amount a change from the imperfect to the ideal state would yield. AIG has many properties desired for quantifying imperfect cognition. The ratio of achieved to ideally obtainable information measures cognitive fidelity and that of AIG to the necessary cognitive effort measures cognitive efficiency. We provide an axiomatic derivation of AIG, illustrate its application at common scenarios of posterior inaccuracies, and discuss the implication of cognitive efficiency for sustainable resource allocation in computational inference., Comment: 14 pages, 5 figures

The EDR and eRASS1 data have already revealed a remarkable number of undiscovered X-ray sources. Using Bayesian inference and generative modeling techniques for X-ray imaging, we aim to increase the sensitivity and scientific value of these observations by denoising, deconvolving, and decomposing the X-ray sky. Leveraging information field theory, we can exploit the spatial and spectral correlation structures of the different physical components of the sky with non-parametric priors to enhance the image reconstruction. By incorporating instrumental effects into the forward model, we develop a comprehensive Bayesian imaging algorithm for eROSITA pointing observations. Finally, we apply the developed algorithm to EDR data of the LMC SN1987A, fusing data sets from observations made by five different telescope modules. The final result is a denoised, deconvolved, and decomposed view of the LMC, which enables the analysis of its fine-scale structures, the creation of point source catalogues of this region, and enhanced calibration for future work.

Understanding the magnetised Universe is a major challenge in modern astrophysics, and cosmic magnetism has been acknowledged as one of the science key drivers of the most ambitious radio instrument ever planned, the SKA telescope. With this work, we aim to investigate the potential of the SKA telescope and its precursors and pathfinders in the study of magnetic fields in galaxy clusters and filaments through diffuse synchrotron radio emission. Galaxy clusters and filaments of the cosmic web are indeed unique laboratories to investigate turbulent fluid motions and large-scale magnetic fields in action and much of what is known about magnetic fields in galaxy clusters comes from sensitive radio observations. Based on cosmological MHD simulations, we predict radio properties (total intensity and polarisation) of a pair of galaxy clusters connected by a cosmic-web filament. We use our theoretical expectations to explore the potential of polarimetric observations to study large-scale structure magnetic fields in the frequency ranges 50-350MHz and 950-1760MHz. We also present predictions for galaxy cluster polarimetric observations with the SKA precursors and pathfinders (LOFAR2.0 and MeerKAT+). Our findings point out that polarisation observations are particularly powerful for the study of large-scale magnetic fields, since they are not significantly affected by confusion noise. The unprecedented sensitivity and spatial resolution of the intermediate frequency radio telescopes make them the favourite instruments for the study of these sources through polarimetric data, potentially allowing us to understand if the energy density of relativistic electrons is in equipartition with the magnetic field or coupled with the thermal gas density. Our results show that low frequency instruments represent as well a precious tool to study diffuse synchrotron emission in total intensity and polarisation., Comment: Accepted for publication on A&A, 12 pages

Reconstructing the electric field from the measured voltages in an antenna, unfolding the antenna response, comes with several problems. Due to the noisiness of the signal it is often necessary to disregard part of the bandwidth of the antenna. It is also not guaranteed, that this system of equations can be inverted at all. In any case, the noise of the measurement will be converted into noise on the electric field. This could be solved by Bayesian inference, however, the electric field is continuous, which would lead to an infinite-dimensional latent space. Information field theory (IFT) has been developed to deal with this problem and allow for Bayesian reasoning on fields. It provides a theoretical backbone and effective tools to approach the inference as a discrete problem in the continuum limit, taking the continuous nature of fields into account. We will present a first working signal model that can be used with IFT-based inference algorithms, which can successfully reconstruct the electric field. The model is based on the current understanding of air shower emission physics, modelling geomagnetic and charge-excess emission and their respective polarisation and spectra separately. Since Bayesian inference provides the posterior distribution, this method also gives an estimate on the uncertainty of the measured field. The performance of this method will be demonstrated with Monte-Carlo simulations of air shower radio signals., Comment: PoS(ARENA2024)

Many advances in astronomy and astrophysics originate from accurate images of the sky emission across multiple wavelengths. This often requires reconstructing spatially and spectrally correlated signals detected from multiple instruments. To facilitate the high-fidelity imaging of these signals, we introduce the universal Bayesian imaging kit (UBIK). Specifically, we present J-UBIK, a flexible and modular implementation leveraging the JAX-accelerated NIFTy.re software as its backend. J-UBIK streamlines the implementation of the key Bayesian inference components, providing for all the necessary steps of Bayesian imaging pipelines. First, it provides adaptable prior models for different sky realizations. Second, it includes likelihood models tailored to specific instruments. So far, the package includes three instruments: Chandra and eROSITA for X-ray observations, and the James Webb Space Telescope (JWST) for the near- and mid-infrared. The aim is to expand this set in the future. Third, these models can be integrated with various inference and optimization schemes, such as maximum a posteriori estimation and variational inference. Explicit demos show how to integrate the individual modules into a full analysis pipeline. Overall, J-UBIK enables efficient generation of high-fidelity images via Bayesian pipelines that can be tailored to specific research objectives.

The galaxy M87 is one of the prime targets for high resolution radio imaging pursuing the ringlike shadow of its supermassive black hole, the innermost regions of accretion flow, and the formation of the relativistic jet. However, it remains challenging to observe both jointly. Only recently, global mm-VLBI array (GMVA)+ALMA observations at 86 GHz in 2018 were able to reconstruct the M87 black hole shadow and the extended jet emission simultaneously. In order to analyze the ring and jet of M87, conventional CLEAN algorithms were mainly employed alongside the RML method SMILI in the previous work. To test the robustness of the reconstructed structures of M87 GMVA+ALMA observations at 86GHz, we estimate the ring diameter, width, and the extended jet emission with the possible central spine by two different novel imaging algorithms: resolve and DoG-HiT. Overall reconstructions are consistent with the results reported in the previous paper. The ring structure of the M87 is resolved at higher resolution and the posterior distribution of M87 ring features is explored. The resolve images show that the ring diameter is 60.9 +- 2.2 muas and width is 16.0 +- 0.9 muas. The ring diameter is 61.0 muas and width is 20.6 muas by DoG-HiT. The ring diameter is therefore in agreement with the estimation (64+4-8 muas) by SMILI and the geometrical modeling. Two bright spots in the ring are reconstructed by four independent imaging methods, the substructure in the ring is therefore most likely originated from the data. A consistent limb-brightened jet structure is reconstructed by resolve and DoG-HiT, albeit with a less pronounced central spine. Modern data-driven imaging methods confirm the ring and jet structure in M87, complementing traditional VLBI methods with novel perspectives on the significance of recovered features. They confirm the result of the previous report., Comment: submitted to A&A

Nested sampling (NS) is a stochastic method for computing the log-evidence of a Bayesian problem. It relies on stochastic estimates of prior volumes enclosed by likelihood contours, which limits the accuracy of the log-evidence calculation. We propose to transform the prior volume estimation into a Bayesian inference problem, which allows us to incorporate a smoothness assumption for likelihood-prior volume relations. As a result, we aim to increase the accuracy of the volume estimates and thus improve the overall log-evidence calculation using NS. The method presented works as a post-processing step for NS and provides posterior samples of the likelihood-prior-volume relation, from which the log-evidence can be calculated. We demonstrate an implementation of the algorithm and compare its results with plain NS on two synthetic datasets for which the underlying evidence is known. We find a significant improvement in accuracy for runs with less than one hundred active samples in NS, but are prone to numerical problems beyond this point.

Faraday rotation contains information about the magnetic field structure along the line of sight and is an important instrument in the study of cosmic magnetism. Traditional Faraday spectrum deconvolution methods such as RMCLEAN face challenges in resolving complex Faraday dispersion functions and handling large datasets. We develop a deep learning deconvolution model to enhance the accuracy and efficiency of extracting Faraday rotation measures from radio astronomical data, specifically targeting data from the MeerKAT Galaxy Cluster Legacy Survey (MGCLS). We use semi-supervised learning, where the model simultaneously recreates the data and minimizes the difference between the output and the true signal of synthetic data. Performance comparisons with RMCLEAN were conducted on simulated as well as real data for the galaxy cluster Abell 3376. Our semi-supervised model is able to recover the Faraday dispersion with great accuracy, particularly for complex or high-RM signals, maintaining sensitivity across a broad RM range. The computational efficiency of this method is significantly improved over traditional methods. Applied to observations of Abell 3376, we find detailed magnetic field structures in the radio relics, and several AGN. We also apply our model to MeerKAT data of Abell 85, Abell 168, Abell 194, Abell 3186 and Abell 3667., Comment: 20 pages, 20 figures, submitted to Astronomy & Astrophysics

Self-calibration methods with the CLEAN algorithm have been widely employed in Very Long Baseline Interferometry (VLBI) data processing in order to correct antenna-based amplitude and phase corruptions present in the data. However, human interaction during the conventional CLEAN self-calibration process can impose a strong effective prior, which in turn may produce artifacts within the final image and hinder the reproducibility of final results. In this work, we aim to demonstrate a combined self-calibration and imaging method for VLBI data in a Bayesian inference framework. The method corrects for amplitude and phase gains for each antenna and polarization mode by inferring the temporal correlation of the gain solutions. We use Stokes I data of M87 taken with the Very Long Baseline Array (VLBA) at 43GHz, pre-calibrated using the rPICARD CASA-based pipeline. For antenna-based gain calibration and imaging, we use the Bayesian imaging software resolve. To estimate gain and image uncertainties, we use a Variational Inference method. We obtain a high-resolution M87 Stokes I image at 43GHz in conjunction with antenna-based gain solutions using our Bayesian self-calibration and imaging method. The core with counter-jet structure is better resolved, and extended jet emission is better described compared to the CLEAN reconstruction. Furthermore, uncertainty estimation of the image and antenna-based gains allows us to quantify the reliability of the result. Our Bayesian self-calibration and imaging method is able to reconstruct robust and reproducible Stokes I images and gain solutions with uncertainty estimation by taking into account the uncertainty information in the data., Comment: accepted for publication in A&A

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., Comment: 21 pages, 22 figures

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.

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., Comment: 24 pages, 14 figures

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., Comment: 6 pages, 3 figures, to be published in Astronomy & Astrophysics (A&A)

(abridged) Ultra-high energy cosmic rays (UHECRs) are extremely energetic charged particles with energies surpassing $10^{18}$ eV. Their sources remain elusive, obscured by deflections caused by the Galactic magnetic field (GMF). This challenge is further complicated by our limited understanding of the three-dimensional structure of the GMF, as current GMF observations consist primarily of quantities integrated along the line-of-sight (LOS). Nevertheless, data from upcoming stellar polarisation surveys along with Gaia's stellar parallax data are expected to yield local GMF measurements.. In this work, we employ methods of Bayesian statistical inference in order to sample the posterior distribution of the GMF within part of the Galaxy. By assuming a known rigidity and arrival direction of an UHECR, we backtrack its trajectory through various GMF configurations drawn from the posterior distribution. Our objective is to rigorously evaluate our algorithm's performance in scenarios that closely mirror the setting of expected future applications. In pursuit of this, we condition the posterior to synthetic integrated LOS measurements of the GMF, in addition to synthetic local POS-component measurements. In this proof of concept work, we assume the ground truth to be a magnetic field produced by a dynamo simulation of the Galactic ISM. Our results demonstrate that for all locations of the observed arrival direction on the POS, our algorithm is able to substantially update our knowledge on the original arrival direction of UHECRs with rigidity $E/Z = 5 \times 10^{19}$ eV, even in the case of complete absence of LOS information. If integrated data is included in the inference, then the regions of the celestial sphere where the maximum error occurs diminishes greatly. Even in those regions the maximum error is diminished by a factor of about $3$ in the specific setting studied., Comment: 11 pages, 7 figures

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., Comment: 10 pages, 2 figures, published in Journal of Open Source Software (JOSS)

Nested sampling provides an estimate of the evidence of a Bayesian inference problem via probing the likelihood as a function of the enclosed prior volume. However, the lack of precise values of the enclosed prior mass of the samples introduces probing noise, which can hamper high-accuracy determinations of the evidence values as estimated from the likelihood-prior-volume function. We introduce an approach based on information field theory, a framework for non-parametric function reconstruction from data, that infers the likelihood-prior-volume function by exploiting its smoothness and thereby aims to improve the evidence calculation. Our method provides posterior samples of the likelihood-prior-volume function that translate into a quantification of the remaining sampling noise for the evidence estimate, or for any other quantity derived from the likelihood-prior-volume function.

An ESO internal ALMA development study, BRAIN, is addressing the ill-posed inverse problem of synthesis image analysis employing astrostatistics and astroinformatics. These emerging fields of research offer interdisciplinary approaches at the intersection of observational astronomy, statistics, algorithm development, and data science. In this study, we provide evidence of the benefits of employing these approaches to ALMA imaging for operational and scientific purposes. We show the potential of two techniques, RESOLVE and DeepFocus, applied to ALMA calibrated science data. Significant advantages are provided with the prospect to improve the quality and completeness of the data products stored in the science archive and overall processing time for operations. Both approaches evidence the logical pathway to address the incoming revolution in data rates dictated by the planned electronic upgrades. Moreover, we bring to the community additional products through a new package, ALMASim, to promote advancements in these fields, providing a refined ALMA simulator usable by a large community for training and/or testing new algorithms., Comment: 9 pages, 5 figures, MaxEnt2023 conference

Given the urgency to reduce fossil fuel energy production to make climate tipping points less likely, we call for resource-aware knowledge gain in the research areas on Universe and Matter with emphasis on the digital transformation. A portfolio of measures is described in detail and then summarized according to the timescales required for their implementation. The measures will both contribute to sustainable research and accelerate scientific progress through increased awareness of resource usage. This work is based on a three-days workshop on sustainability in digital transformation held in May 2023., Comment: 20 pages, 2 figures, publication following workshop 'Sustainability in the Digital Transformation of Basic Research on Universe & Matter', 30 May to 2 June 2023, Meinerzhagen, Germany, https://indico.desy.de/event/37480

Sagittarius A* (Sgr A*), the supermassive black hole at the heart of our galaxy, provides unique opportunities to study black hole accretion, jet formation, and gravitational physics. The rapid structural changes in Sgr A*'s emission pose a significant challenge for traditional imaging techniques. We present dynamic reconstructions of Sgr A* using Event Horizon Telescope (EHT) data from April 6th and 7th, 2017, analyzed with a one-minute temporal resolution with the Resolve framework. This Bayesian approach employs adaptive Gaussian Processes and Variational Inference for data-driven self-regularization. Our results not only fully confirm the initial findings by the EHT Collaboration for a time-averaged source but also reveal intricate details about the temporal dynamics within the black hole environment. We find an intriguing dynamic feature on April 6th that propagates in a clock-wise direction. Geometric modelling with ray-tracing, although not fully conclusive, indicates compatibility with high-inclination configurations of about $\theta_o = 160^\circ$, as seen in other studies.

We explore quantum measurement in the context of Everettian unitary quantum mechanics and construct an explicit unitary measurement procedure. We propose the existence of prior correlated states that enable this procedure to work and therefore argue that correlation is a resource that is consumed when measurements take place. It is also argued that a network of such measurements establishes a stable objective classical reality., Comment: 12 pages, 6 figures. Presented at the MaxEnt 23 conference

Supernovae are an important source of energy in the interstellar medium. Young remnants of supernovae have a peak emission in the X-ray region, making them interesting objects for X-ray observations. In particular, the supernova remnant SN1006 is of great interest due to its historical record, proximity and brightness. It has therefore been studied by several X-ray telescopes. Improving the X-ray imaging of this and other remnants is important but challenging as it requires to address a spatially varying instrument response in order to achieve a high signal-to-noise ratio. Here, we use Chandra observations to demonstrate the capabilities of Bayesian image reconstruction using information field theory. Our objective is to reconstruct denoised, deconvolved and spatio-spectral resolved images from X-ray observations and to decompose the emission into different morphologies, namely diffuse and point-like. Further, we aim to fuse data from different detectors and pointings into a mosaic and quantify the uncertainty of our result. Utilizing prior knowledge on the spatial and spectral correlation structure of the two components, diffuse emission and point sources, the presented method allows the effective decomposition of the signal into these. In order to accelerate the imaging process, we introduce a multi-step approach, in which the spatial reconstruction obtained for a single energy range is used to derive an informed starting point for the full spatio-spectral reconstruction. The method is applied to 11 Chandra observations of SN1006 from 2008 and 2012, providing a detailed, denoised and decomposed view of the remnant. In particular, the separated view of the diffuse emission should provide new insights into its complex small-scale structures in the center of the remnant and at the shock front profiles.

We present results of a wideband high-resolution polarization study of Hydra A, one of the most luminous FR I radio galaxies known and amongst the most well-studied. The radio emission from this source displays extremely large Faraday rotation measures (RM), ranging from -12300 rad m$^{-2}$ to 5000 rad m$^{-2}$, the majority of which are believed to originate from magnetized thermal gas external to the radio tails. The radio emission from both tails strongly depolarizes with decreasing frequency. The depolarization, as a function of wavelength, is commonly non-monotonic, often showing oscillatory behavior, with strongly non-linear rotation of the polarization position angle with $\lambda^2$. A simple model, based on the RM screen derived from the high frequency, high resolution data, predicts the lower frequency depolarization remarkably well. The success of this model indicates the majority of the depolarization can be attributed to fluctuations in the magnetic field on scales $< 1500$ pc, suggesting the presence of turbulent magnetic field/electron density structures on sub-kpc scales within a Faraday rotating (FR) medium., Comment: 21 pages, 19 figures, Accepted for publication in The Astrophysical Journal

High-resolution 3D maps of interstellar dust are critical for probing the underlying physics shaping the structure of the interstellar medium, and for foreground correction of astrophysical observations affected by dust. We aim to construct a new 3D map of the spatial distribution of interstellar dust extinction out to a distance of 1.25 kpc from the Sun. We leveraged distance and extinction estimates to 54 million nearby stars derived from the Gaia BP/RP spectra. Using the stellar distance and extinction information, we inferred the spatial distribution of dust extinction. We modeled the logarithmic dust extinction with a Gaussian process in a spherical coordinate system via iterative charted refinement and a correlation kernel inferred in previous work. In total, our posterior has over 661 million degrees of freedom. We probed the posterior distribution using the variational inference method MGVI. Our 3D dust map has an angular resolution of up to 14' (Nside = 256), and we achieve parsec-scale distance resolution, sampling the dust in 516 logarithmically spaced distance bins spanning 69 pc to 1250 pc. We generated 12 samples from the variational posterior of the 3D dust distribution and release the samples alongside the mean 3D dust map and its corresponding uncertainty. Our map resolves the internal structure of hundreds of molecular clouds in the solar neighborhood and will be broadly useful for studies of star formation, Galactic structure, and young stellar populations. It is available for download in a variety of coordinate systems online and can also be queried via the publicly available dustmaps Python package., Comment: 27 pages, 24 figures, 1 table, to be published in Astronomy & Astrophysics (A&A)

The concept of attention, numerical weights that emphasize the importance of particular data, has proven to be very relevant in artificial intelligence. Relative entropy (RE, aka Kullback-Leibler divergence) plays a central role in communication theory. Here we combine these concepts, attention and RE. RE guides optimal encoding of messages in bandwidth-limited communication as well as optimal message decoding via the maximum entropy principle (MEP). In the coding scenario, RE can be derived from four requirements, namely being analytical, local, proper, and calibrated. Weighted RE, used for attention steering in communications, turns out to be improper. To see how proper attention communication can emerge, we analyze a scenario of a message sender who wants to ensure that the receiver of the message can perform well-informed actions. If the receiver decodes the message using the MEP, the sender only needs to know the receiver's utility function to inform optimally, but not the receiver's initial knowledge state. In case only the curvature of the utility function maxima are known, it becomes desirable to accurately communicate an attention function, in this case a by this curvature weighted and re-normalized probability function. Entropic attention communication is here proposed as the desired generalization of entropic communication that permits weighting while being proper, thereby aiding the design of optimal communication protocols in technical applications and helping to understand human communication. For example, our analysis shows how to derive the level of cooperation expected under misaligned interests of otherwise honest communication partners., Comment: 25 pages, 4 figures, re-submitted

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., Comment: 13 pages, 9 figures

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 $\alpha$ 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 $\mu$G, with an all-sky root-mean-square of 1.1 $\mu$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., Comment: 27 pages, 18 Figures, accepted in A&A The results are accessible via https://zenodo.org/records/10523170 . Major differences to the previous versions are a new appendix on test with synthetic data and a more detailed discussion of the method. Version 1 comments: 24 pages, 16 Figures, submitted to A&A

(abridged) Ultra-high energy cosmic rays (UHECRs) are highly energetic charged particles with energies exceeding $10^{18}$ eV. Identifying their sources and production mechanism can provide insight into many open questions in astrophysics and high energy physics. However, the Galactic magnetic field (GMF) deflects UHECRs, and the high uncertainties in our current understanding of the $3$-dimensional structure of the GMF does not permit us to accurately determine their true arrival direction on the plane of the sky (PoS). This difficulty arises from the fact that currently all GMF observations are integrated along the line-of-sight (LoS). Upcoming stellar optopolarimetric surveys as well as Gaia data on stellar parallaxes, are expected to provide local measurements of the GMF in the near future. In this paper, we evaluate the reconstruction of the GMF in a limited region of the Galaxy given sparse and local GMF measurements within that region, through Bayesian inference using principles of Information Field Theory. We backtrack UHECRs through GMF configurations drawn from the posterior to improve our knowledge of their true arrival directions. We show that, for a weakly turbulent GMF, it is possible to correct for its effect on the observed arrival direction of UHECRs to within $\sim 3^\circ$. For completely turbulent fields, we show that our procedure can still be used to significantly improve our knowledge on the true arrival direction of UHECRs., Comment: 12 pages, 10 figures

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}$., Comment: 12 pages, 6 figures

This article explores how measurement can be understood in the context of a universe evolving according to unitary (reversible) quantum dynamics. A unitary measurement procedure is developed consistent with the non-measurement axioms of quantum mechanics wherein the system being measured and the observer become correlated. It is argued that for this to work the correlation necessarily has to be transferred from somewhere else. Thus, correlation is a resource that is consumed when measurements take place. It is also argued that a network of such measurements establishes a stable objective classical reality, especially in the context of repeatability of experiments., Comment: 13 pages, 5 figures

Barnard's Loop is a famous arc of H$\alpha$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Brice\~{n}o-1 (OBP-B1) cluster in the middle of a large dust cavity. The large-scale gas in the region also appears to be expanding from a central point, originally proposed to be Orion X. OBP-B1 appears to serve as another possible center, and we evaluate whether Orion X or OBP-B1 is more likely to be the cause of the expansion. We find that neither cluster served as the single expansion center, but rather a combination of feedback from both likely propelled the expansion. Recent 3D dust maps are used to characterize the 3D topology of the entire region, which shows Barnard's Loop's correspondence with a large dust cavity around the OPB-B1 cluster. The molecular clouds Orion A, Orion B, and Orion $\lambda$ reside on the shell of this cavity. Simple estimates of gravitational effects from both stars and gas indicate that the expansion of this asymmetric cavity likely induced anisotropy in the kinematics of OBP-B1. We conclude that feedback from OBP-B1 has affected the structure of the Orion A, Orion B, and Orion $\lambda$ molecular clouds and may have played a major role in the formation of Barnard's Loop., Comment: 25 pages, 10 figures. Submitted to ApJ

The Atacama Large Millimeter/submillimeter Array (ALMA) is currently revolutionizing observational astrophysics. The aperture synthesis technique provides angular resolution otherwise unachievable with the conventional single-aperture telescope. However, recovering the image from the inherently undersampled data is a challenging task. The CLEAN algorithm has proven successful and reliable and is commonly used in imaging the interferometric observations. It is not, however, free of limitations. Point-source assumption, central to the CLEAN is not optimal for the extended structures of molecular gas recovered by ALMA. Additionally, negative fluxes recovered with CLEAN are not physical. This begs to search for alternatives that would be better suited for specific science cases. We present the recent developments in imaging ALMA data using Bayesian inference techniques, namely the RESOLVE algorithm This algorithm, based on information field theory \cite{Ensslin2013}, has been already successfully applied to image the Very Large Array data. We compare the capability of both CLEAN and RESOLVE to recover known sky signal, convoluted with the simulator of ALMA observation data and we investigate the problem with a set of actual ALMA observations., Comment: Proceedings of International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, IHP, Paris, July 18-22, 2022

The Atacama Large Millimeter/submillimeter Array with the planned electronic upgrades will deliver an unprecedented amount of deep and high resolution observations. Wider fields of view are possible with the consequential cost of image reconstruction. Alternatives to commonly used applications in image processing have to be sought and tested. Advanced image reconstruction methods are critical to meet the data requirements needed for operational purposes. Astrostatistics and astroinformatics techniques are employed. Evidence is given that these interdisciplinary fields of study applied to synthesis imaging meet the Big Data challenges and have the potentials to enable new scientific discoveries in radio astronomy and astrophysics., Comment: 8 pages, 5 figures, proceedings International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, IHP, Paris, July 18-22, 2022

Gaussian Processes (GPs) are highly expressive, probabilistic models. A major limitation is their computational complexity. Naively, exact GP inference requires $\mathcal{O}(N^3)$ computations with $N$ denoting the number of modeled points. Current approaches to overcome this limitation either rely on sparse, structured or stochastic representations of data or kernel respectively and usually involve nested optimizations to evaluate a GP. We present a new, generative method named Iterative Charted Refinement (ICR) to model GPs on nearly arbitrarily spaced points in $\mathcal{O}(N)$ time for decaying kernels without nested optimizations. ICR represents long- as well as short-range correlations by combining views of the modeled locations at varying resolutions with a user-provided coordinate chart. In our experiment with points whose spacings vary over two orders of magnitude, ICR's accuracy is comparable to state-of-the-art GP methods. ICR outperforms existing methods in terms of computational speed by one order of magnitude on the CPU and GPU and has already been successfully applied to model a GP with $122$ billion parameters.

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., Comment: 14 pages, 12 figures, accepted for publication on MNRAS

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., Comment: Please see the journal version for a language-edited release of this manuscript (open-access): https://doi.org/10.1051/0004-6361/202243819

The galaxy cluster Abell 523 (A523) hosts an extended diffuse synchrotron source historically classified as a radio halo. Its radio power at 1.4 GHz makes it one of the most significant outliers in the scaling relations between observables derived from multi-wavelength observations of galaxy clusters: it has a morphology that is different and offset from the thermal gas, and it has polarized emission at 1.4 GHz typically difficult to observe for this class of sources. A magnetic field fluctuating on large spatial scales (~ 1 Mpc) can explain these peculiarities but the formation mechanism for this source is not yet completely clear. To investigate its formation mechanism, we present new observations obtained with the LOw Frequency ARray at 120-168 MHz and the Jansky Very Large Array at 1-2 GHz, which allow us to study the spectral index distribution of this source. According to our data the source is observed to be more extended at 144 MHz than previously inferred at 1.4 GHz, with a total size of about 1.8 Mpc and a flux density S_144MHz = (1.52 +- 0.31) Jy. The spectral index distribution of the source is patchy with an average spectral index alpha ~ 1.2 between 144 MHz and 1.410 GHz, while an integrated spectral index alpha ~ 2.1 has been obtained between 1.410 GHz and 1.782 GHz. A previously unseen patch of steep spectrum emission is clearly detected at 144 MHz in the south of the cluster. Overall, our findings suggest that we are observing an overlapping of different structures, powered by the turbulence associated with the primary and a possible secondary merger., Comment: 19 pages, 11 figures, accepted for publication on MNRAS

Information field theory (IFT), the information theory for fields, is a mathematical framework for signal reconstruction and non-parametric inverse problems. Artificial intelligence (AI) and machine learning (ML) aim at generating intelligent systems including such for perception, cognition, and learning. This overlaps with IFT, which is designed to address perception, reasoning, and inference tasks. Here, the relation between concepts and tools in IFT and those in AI and ML research are discussed. In the context of IFT, fields denote physical quantities that change continuously as a function of space (and time) and information theory refers to Bayesian probabilistic logic equipped with the associated entropic information measures. Reconstructing a signal with IFT is a computational problem similar to training a generative neural network (GNN) in ML. In this paper, the process of inference in IFT is reformulated in terms of GNN training. In contrast to classical neural networks, IFT based GNNs can operate without pre-training thanks to incorporating expert knowledge into their architecture. Furthermore, the cross-fertilization of variational inference methods used in IFT and ML are discussed. These discussions suggests that IFT is well suited to address many problems in AI and ML research and application., Comment: 12 pages, three figures, invited talk at MaxEnt2020/2021, reviewed and published by Entropy

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds on and improves past and currently proposed missions at such energies. We outline the motivations for combining a coded mask and a Compton telescope and we define the scientific goals of such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which enable identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena., Comment: 16 pages. Contact authors: Elena Orlando, Eugenio Bottacini

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles., Comment: arXiv admin note: substantial text overlap with arXiv:1910.08028, arXiv:1707.01348

Neural Networks play a growing role in many science disciplines, including physics. Variational Autoencoders (VAEs) are neural networks that are able to represent the essential information of a high dimensional data set in a low dimensional latent space, which have a probabilistic interpretation. In particular the so-called encoder network, the first part of the VAE, which maps its input onto a position in latent space, additionally provides uncertainty information in terms of a variance around this position. In this work, an extension to the Autoencoder architecture is introduced, the FisherNet. In this architecture, the latent space uncertainty is not generated using an additional information channel in the encoder, but derived from the decoder, by means of the Fisher information metric. This architecture has advantages from a theoretical point of view as it provides a direct uncertainty quantification derived from the model, and also accounts for uncertainty cross-correlations. We can show experimentally that the FisherNet produces more accurate data reconstructions than a comparable VAE and its learning performance also apparently scales better with the number of latent space dimensions., Comment: 25 pages, 11 figures

We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well-described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is $\approx 3:1$. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, $n$, falling off like $n^{-2}$ at large radii. Using Bayesian model selection, we find that parameterizing the clouds' profiles using a single Gaussian is disfavored. We compare our results with 2D dust extinction maps, finding that the 3D dust recovers the total cloud mass from integrated approaches with fidelity, deviating only at higher levels of extinction ($A_V \gtrsim 2 - 3$ mag). The 3D cloud structure described here will enable comparisons with synthetic clouds generated in simulations, offering unprecedented insight into the origins and fates of molecular clouds in the interstellar medium., Comment: Published in ApJ. For data behind the figures see https://dataverse.harvard.edu/dataverse/3D_cloud_structure . For interactive gallery of local clouds see https://faun.rc.fas.harvard.edu/czucker/Paper_Figures/3D_Cloud_Topologies/gallery.html . Companion paper to Bialy et al. 2021 (arXiv:2109.09763). For press see https://www.cfa.harvard.edu/news/gigantic-cavity-space-sheds-new-light-how-stars-form

A major question in the field of star formation is how molecular clouds form out of the diffuse Interstellar Medium (ISM). Recent advances in 3D dust mapping are revolutionizing our view of the structure of the ISM. Using the highest-resolution 3D dust map to date, we explore the structure of a nearby star-forming region, which includes the well-known Perseus and Taurus molecular clouds. We reveal an extended near-spherical shell, 156 pc in diameter, hereafter the "Per-Tau Shell", in which the Perseus and Taurus clouds are embedded. We also find a large ring structure at the location of Taurus, hereafter, the "Tau Ring". We discuss a formation scenario for the Per-Tau Shell, in which previous stellar and supernova (SN) feedback events formed a large expanding shell, where the swept-up ISM has condensed to form both the shell and the Perseus and Taurus molecular clouds within it. We present auxiliary observations of HI, H$\alpha$, $^{26}$Al, and X-rays that further support this scenario, and estimate Per-Tau Shell's age to be $\approx 6-22$ Myrs. The Per-Tau shell offers the first three-dimensional observational view of a phenomenon long-hypothesized theoretically, molecular cloud formation and star formation triggered by previous stellar and SN feedback., Comment: ApJ Letters, accepted. Main text: 8 pages, 4 figures. An Interactive + Augmented Reality Figure at https://faun.rc.fas.harvard.edu/czucker/Paper_Figures/sbialy/pertau_superbubble.html

Reputation is a central element of social communications, be it with human or artificial intelligence (AI), and as such can be the primary target of malicious communication strategies. There is already a vast amount of literature on trust networks addressing this issue and proposing ways to simulate these networks dynamics using Bayesian principles and involving Theory of Mind models. The main issue for these simulations is usually the amount of information that can be stored and is usually solved by discretising variables and using hard thresholds. Here we propose a novel approach to the way information is updated that accounts for knowledge uncertainty and is closer to reality. In our game, agents use information compression techniques to capture their complex environment and store it in their finite memories. The loss of information that results from this leads to emergent phenomena, such as echo chambers, self-deception, deception symbiosis, and freezing of group opinions. Various malicious strategies of agents are studied for their impact on group sociology, like sycophancy, egocentricity, pathological lying, and aggressiveness. Even though our modeling could be made more complex, our set-up can already provide insights into social interactions and can be used to investigate the effects of various communication strategies and find ways to counteract malicious ones. Eventually this work should help to safeguard the design of non-abusive AI systems., Comment: 36 pages main part + 17 pages appendix, 28 figures, 3 tables, latex glitches fixed, revised and resubmitted

The viral load of patients infected with SARS-CoV-2 varies on logarithmic scales and possibly with age. Controversial claims have been made in the literature regarding whether the viral load distribution actually depends on the age of the patients. Such a dependence would have implications for the COVID-19 spreading mechanism, the age-dependent immune system reaction, and thus for policymaking. We hereby develop a method to analyze viral-load distribution data as a function of the patients' age within a flexible, non-parametric, hierarchical, Bayesian, and causal model. The causal nature of the developed reconstruction additionally allows to test for bias in the data. This could be due to, e.g., bias in patient-testing and data collection or systematic errors in the measurement of the viral load. We perform these tests by calculating the Bayesian evidence for each implied possible causal direction. The possibility of testing for bias in data collection and identifying causal directions can be very useful in other contexts as well. For this reason we make our model freely available. When applied to publicly available age and SARS-CoV-2 viral load data, we find a statistically significant increase in the viral load with age, but only for one of the two analyzed datasets. If we consider this dataset, and based on the current understanding of viral load's impact on patients' infectivity, we expect a non-negligible difference in the infectivity of different age groups. This difference is nonetheless too small to justify considering any age group as noninfectious.

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., Comment: 42 pages, 18 figures, accepted by Entropy

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., Comment: 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

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.

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., Comment: 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

Knowledge on evolving physical fields is of paramount importance in science, technology, and economics. Dynamical field inference (DFI) addresses the problem of reconstructing a stochastically driven, dynamically evolving field from finite data. It relies on information field theory (IFT), the information theory for fields. Here, the relations of DFI, IFT, and the recently developed supersymmetric theory of stochastics (STS) are established in a pedagogical discussion. In IFT, field expectation values can be calculated from the partition function of the full space-time inference problem. The partition function of the inference problem invokes a functional Dirac function to guarantee the dynamics, as well as a field-dependent functional determinant, to establish proper normalization, both impeding the necessary evaluation of the path integral over all field configurations. STS replaces these problematic expressions via the introduction of fermionic ghost and bosonic Lagrange fields, respectively. The action of these fields has a supersymmetry, which means there exists an exchange operation between bosons and fermions that leaves the system invariant. In contrast to this, measurements of the dynamical fields do not adhere to this supersymmetry. The supersymmetry can also be broken spontaneously, in which case the system evolves chaotically. This affects the predictability of the system and thereby make DFI more challenging. We investigate the interplay of measurement constraints with the non-linear chaotic dynamics of a simplified, illustrative system with the help of Feynman diagrams and show that the Fermionic corrections are essential to obtain the correct posterior statistics over system trajectories., Comment: 20 pages, 2 figures, 6 Feynman diagrams, 162 numbered equations

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., Comment: 39 pages, 15 figures, 1 table, for code, see https://gitlab.mpcdf.mpg.de/jharthki/bdc