Explore the vast range of titles available.

Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of Δβ d ≲ 0.01 and thus test models of dust composition that predict that β d in polarization differs from that measured in total intensity; measure the correlation coefficient between polarized dust and synchrotron emission with a factor of two greater precision than current constraints; exclude the nonexistence of exo-Oort clouds at roughly 2.9σ if the true fraction is similar to the detection rate of giant planets; map more than 850 molecular clouds with at least 50 independent polarization measurements at 1 pc resolution; detect or place upper limits on the polarization fractions of CO(2–1) emission and anomalous microwave emission at the 0.1% level in select regions; and measure the correlation coefficient between optical starlight polarization and microwave polarized dust emission in 1° patches for all lines of sight with N H ≳ 2 × 1020 cm−2. The goals and forecasts outlined here provide a roadmap for other microwave polarization experiments to expand their scientific scope via Milky Way astrophysics. 37 37 A supplement describing author contributions to this paper can be found at https://simonsobservatory.org/wp-content/uploads/2022/02/SO_GS_Contributions.pdf.

Polarized foreground emission from the Galaxy is one of the biggest challenges facing current and upcoming cosmic microwave background (CMB) polarization experiments. We develop new models of polarized Galactic dust and synchrotron emission at CMB frequencies that draw on the latest observational constraints; that employ the “polarization fraction tensor” framework to couple intensity and polarization in a physically motivated way; and that allow for stochastic realizations of small-scale structure at subarcminute angular scales currently unconstrained by full-sky data. We implement these models into the publicly available Python Sky Model (PySM) software and additionally provide PySM interfaces to select models of dust and CO emission from the literature. We characterize the behavior of each model by quantitatively comparing it to observational constraints in both maps and power spectra, demonstrating an overall improvement over previous PySM models. Finally, we synthesize models of the various Galactic foreground components into a coherent suite of three plausible microwave skies that span a range of astrophysical complexity allowed by current data. Author contributions to this paper can be found at the end of this work.

We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a mapmaker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the mapmaker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions, and input beams. We additionally show the reconstruction quality as a function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within an ∼1.5% error for a multipole range ℓ = 30–700 and an ∼0.5% error for a multipole range ℓ = 50–200. We conclude by using a harmonic-domain component separation algorithm to verify that the beam reconstruction errors and biases observed in our analysis do not significantly bias the Simons Observatory r-measurement

We derived astroparticle constraints in different dark matter scenarios that are alternatives to cold dark matter (CDM): thermal relic warm dark matter, WDM; fuzzy dark matter, ψDM; self-interacting dark matter, SIDM; sterile neutrino dark matter, νDM. Our framework is based on updated determinations of the high-redshift UV luminosity functions for primordial galaxies to redshift z∼10, on redshift-dependent halo mass functions in the above DM scenarios from numerical simulations, and on robust constraints on the reionization history of the Universe from recent astrophysical and cosmological datasets. First, we built an empirical model of cosmic reionization characterized by two parameters, namely the escape fraction fesc of ionizing photons from primordial galaxies, and the limiting UV magnitude MUVlim down to which the extrapolated UV luminosity functions steeply increased. Second, we performed standard abundance matching of the UV luminosity function and the halo mass function, obtaining a relationship between UV luminosity and the halo mass, whose shape depends on an astroparticle quantity X specific to each DM scenario (e.g., WDM particle mass); we exploited such a relationship to introduce (in the analysis) a constraint from primordial galaxy formation, in terms of the threshold halo mass above which primordial galaxies can efficiently form stars. Third, we performed Bayesian inference on the three parameters fesc, MUVlim, and X via a standard MCMC technique, and compared the outcomes of different DM scenarios on the reionization history. We also investigated the robustness of our findings against educated variations of still uncertain astrophysical quantities. Finally, we highlight the relevance of our astroparticle estimates in predicting the behavior of the high-redshift UV luminosity function at faint, yet unexplored magnitudes, which may be tested with the advent of the James Webb Space Telescope.

We present ForSE (Foreground Scale Extender), a novel Python package which aims at overcoming the current limitations in the simulation of diffuse Galactic radiation, in the context of Cosmic Microwave Background experiments (CMB). ForSE exploits the ability of generative adversarial neural networks (GANs) to learn and reproduce complex features present in a set of images, with the goal of simulating realistic and non-Gaussian foreground radiation at sub-degree angular scales. This is of great importance in order to estimate the foreground contamination to lensing reconstruction, de-lensing and primordial B-modes, for future CMB experiments. We applied this algorithm to Galactic thermal dust emission in both total intensity and polarization. Our results show how ForSE is able to generate small scale features (at 12 arc-minutes) having as input the large scale ones (80 arc-minutes). The injected structures have statistical properties, evaluated by means of the Minkowski functionals, in good agreement with those of the real sky and which show the correct amplitude scaling as a function of the angular dimension. The obtained thermal dust Stokes Q and U full sky maps as well as the ForSE package are publicly available for download.

We present ForSE+, a Python package that produces non-Gaussian diffuse Galactic thermal dust emission maps at arcminute angular scales and that has the capacity to generate random realizations of small scales. This represents an extension of the ForSE (Foreground Scale Extender) package, which was recently proposed to simulate non-Gaussian small scales of thermal dust emission using generative adversarial networks (GANs). With the input of the large-scale polarization maps from observations, ForSE+ has been trained to produce realistic polarized small scales at 3' following the statistical properties, mainly the non-Gaussianity, of observed intensity small scales, which are evaluated through Minkowski functionals. Furthermore, by adding different realizations of random components to the large-scale foregrounds, we show that ForSE+ is able to generate small scales in a stochastic way. In both cases, the output small scales have a similar level of non-Gaussianity compared with real observations and correct amplitude scaling as a power law. These realistic new maps will be useful, in the future, to understand the impact of non-Gaussian foregrounds on the measurements of the cosmic microwave background (CMB) signal, particularly on the lensing reconstruction, de-lensing, and the detection of cosmological gravitational waves in CMB polarization B-modes.

We describe a novel method for the application of Convolutional Neural Networks (CNNs) to fields defined on the sphere, using the HEALPix tessellation scheme. Specifically, We have developed a pixel-based approach to implement convolutional layers on the spherical surface, similarly to what is commonly done for CNNs in Euclidian space. The algorithm is fully integrable with existing libraries for NNs (e.g., PyTorch or TensorFlow). We present two applications: (i) recognition of handwritten digits projected on the sphere; (ii) estimation of cosmological parameter from Cosmic Microwave Background (CMB) simulated maps. We have built a simple NN architecture, consisting in four convolutional+pooling layers, and have used it for all the applications explored herein. For what concerns the handwritten digits, our CNN reaches an accuracy of about 95%, comparable with other existing spherical CNNs. For CMB applications, we have tested the CNN on the estimation of a \"mock\" parameter, defining the angular scale at which the power spectrum of a Gaussian field projected on the sphere peaks. We have estimated this parameter directly from maps, in several cases: temperature and polarization, presence of noise and partial sky coverage. In all the cases, the NN performances are comparable with those from standard spectrum-based bayesian methods. We demonstrate, for the first time, the capability of CNNs to extract information from polarization fields and to distinguish between E and B-modes. Lastly, we have applied our CNN to the estimation of the Thomson scattering optical depth at reionization (tau) from simulated CMB maps. Even without any specific optimization of the NN architecture, we reach an accuracy comparable with standard bayesian methods. This work represents a first step towards the exploitation of NNs in CMB parameter estimation and demonstrates the feasibility of our approach.

We present BROOM, a new python package for the application of blind, minimum-variance component-separation techniques to microwave observations. The package enables the reconstruction of signals with known spectral energy distributions, such as the Cosmic Microwave Background (CMB), Sunyaev--Zeldovich distortions, or foreground moments, in both temperature and polarization through a suite of Internal Linear Combination (ILC) implementations, in the presence of astrophysical and instrumental contaminants. In addition, BROOM supports the blind reconstruction of coherent emission components with unknown covariance properties via a Generalized ILC (GILC) framework. Beyond component separation, the package provides tools to diagnose foreground complexity and to estimate residual contamination leaking into reconstructed maps across angular scales and sky regions. It also includes utilities to generate realistic microwave simulations for arbitrary CMB experiments and to compute angular power spectra of the resulting products. We present a comprehensive description and validation of the implemented pipelines in two representative experimental configurations: a full-sky satellite mission and a ground-based experiment. BROOM is publicly available, fully documented, and easily installable at https://github.com/alecarones/broom

Observations of the Cosmic Microwave Background (CMB) B-mode polarisation provide a unique probe of inflationary physics. Extracting a reliable constraint on the tensor-to-scalar ratio \\(r\\) nonetheless demands stringent suppression of diffuse Galactic foregrounds, whose residuals can bias the inferred signal. This work introduces and evaluates two extensions of the Needlet Internal Linear Combination (NILC) framework aimed at reducing foreground-induced biases on \\(r\\). The first extension implements the deprojection of selected foreground moments directly within the component-separation step. The second performs a likelihood-level marginalisation over residual foreground power using a data-driven template. Using Simons Observatory Small Aperture Telescope (SO-SAT) - like simulations, we show that both methods effectively control residual contamination, yielding unbiased estimates of \\(r\\) and a consistent reconstruction of the lensing B-mode amplitude. These results indicate that enhanced foreground-mitigation strategies will be useful for next-generation CMB polarisation analyses seeking a robust detection of primordial B-modes.