The Python Sky Model 3 software
read the original abstract
The Python Sky Model (PySM) is a Python package used by Cosmic Microwave Background (CMB) experiments to simulate maps, in HEALPix pixelization, of the various diffuse astrophysical components of Galactic emission relevant at CMB frequencies (i.e. dust, synchrotron, free-free and Anomalous Microwave Emission), as well as the CMB itself. These maps may be integrated over a given instrument bandpass and smoothed with a given instrument beam. PySM 2, released in 2016, has become the de-facto standard for simulating Galactic emission, for example it is used by CMB-S4, Simons Observatory, LiteBird, PICO, CLASS, POLARBEAR and other CMB experiments, as shown by the 80+ citations of the PySM 2 publication. As the resolution of upcoming experiments increases, the PySM 2 software has started to show some limitations, the solution to these issues was to reimplement PySM from scratch focusing on these features: reimplement all the models with the numba Just-In-Time compiler for Python to reduce memory overhead and optimize performance; use MPI through mpi4py to coordinate execution of PySM 3 across multiple nodes and rely on libsharp, for distributed spherical harmonic transforms; employ the data utilities infrastructure provided by astropy to download the input templates and cache them when requested. At this stage we strive to maintain full compatibility with PySM 2, therefore we implement the exact same astrophysical emission models with the same naming scheme. In the extensive test suite we compare the output of each PySM 3 model with the results obtained by PySM 2.
This paper has not been read by Pith yet.
Forward citations
Cited by 6 Pith papers
-
Non-Gaussianity in SMICA
New SMICA formalism and binned bispectrum estimator jointly recover power spectra, spectral parameters, foreground 3-point correlators, and primordial non-Gaussianity constraints from multi-frequency polarization maps...
-
Deep Learning for CMB Foreground Removal and Beam Deconvolution: A U-Net GAN Approach
A U-Net GAN reconstructs CMB T and E maps from Planck-like simulations with foregrounds and systematics, achieving under 1% error outside the Galactic region and demonstrating first-time correction for non-circular be...
-
Forecasts of CMB $E$-mode anomalies for AliCPT-1
A divide-and-conquer framework using QAOA and neural network surrogates accelerates constrained MCMC by factors of 7.6 to 20.3 over classical methods.
-
Synchrotron and free-free mapping with simulated REACH observations between 50-170 MHz
Simulations of REACH observations demonstrate that complex foreground models are needed to recover the 21cm signal but limit accurate free-free mapping while synchrotron is recovered well.
-
Forecasts of CMB $E$-mode anomalies for AliCPT-1
Forecasts indicate AliCPT combined with Simons Observatory can detect injected E-mode dipole modulation at 99% confidence, while AliCPT alone risks biases in alignment and parity tests due to limited sky coverage.
-
BROOM: a python package for model-independent analysis of microwave astronomical data
BROOM is a Python package that applies ILC and GILC techniques for model-independent separation of CMB, SZ, and foreground signals in microwave data along with diagnostic and simulation utilities.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.