maxsmooth: Rapid Maximally Smooth Function Fitting With Applications in Global 21-cm Cosmology

Abstract

Maximally Smooth Functions (MSFs) are a form of constrained functions in which there are no inflection points or zero crossings in high order derivatives. Consequently, they have applications to signal recovery in experiments where signals of interest are expected to be non-smooth features masked by larger smooth signals or foregrounds. They can also act as a powerful tool for diagnosing the presence of systematics. The constrained nature of MSFs makes fitting these functions a non-trivial task. We introduce maxsmooth, an open source package that uses quadratic programming to rapidly fit MSFs. We demonstrate the efficiency and reliability of maxsmooth by comparison to commonly used fitting routines and show that we can reduce the fitting time by approximately two orders of magnitude. We introduce and implement with maxsmooth Partially Smooth Functions, which are useful for describing elements of non-smooth structure in foregrounds. This work has been motivated by the problem of foreground modelling in 21-cm cosmology. We discuss applications of maxsmooth to 21-cm cosmology and highlight this with examples using data from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) and the Large-aperture Experiment to Detect the Dark Ages (LEDA) experiments. We demonstrate the presence of a sinusoidal systematic in the EDGES data with a log-evidence difference of $86.19\pm 0.12$ when compared to a pure foreground fit. MSFs are applied to data from LEDA for the first time in this paper and we identify the presence of sinusoidal systematics. maxsmooth is pip installable and available for download at: https://github.com/htjb/maxsmooth

arXiv comment21 pages, 16 figures

References in corpus (43)

• Wilkinson Microwave Anisotropy Probe (WMAP) Three Year Results: Implications for Cosmology
• Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• The Formation of the First Star in the Universe
• An absorption profile centred at 78 megahertz in the sky-averaged spectrum
• Importance Nested Sampling and the MultiNest Algorithm
• 21-cm cosmology
• Evidence for Reionization at z ~ 6: Detection of a Gunn-Peterson Trough in a z=6.28 Quasar
• ARCADE 2 Measurement of the Extra-Galactic Sky Temperature at 3-90 GHz
• Insights on Dark Matter from Hydrogen during Cosmic Dawn
• Enhanced global signal of neutral hydrogen due to excess radiation at cosmic dawn
• Severely Constraining Dark Matter Interpretations of the 21-cm Anomaly
• Concerns about Modelling of the EDGES Data
• Signs of Dark Matter at 21-cm?
• Energy Feedback from X-ray Binaries in the Early Universe
• Modeling the Radio Background from the First Black Holes at Cosmic Dawn: Implications for the 21 cm Absorption Amplitude
• Early-Universe Constraints on Dark Matter-Baryon Scattering and their Implications for a Global 21cm Signal
• Charting the Parameter Space of the Global 21-cm Signal
• Signature of Excess Radio Background in the 21-cm Global Signal and Power Spectrum
• Constraining the unexplored period between reionization and the dark ages with observations of the global 21 cm signal
• What does the first highly-redshifted 21-cm detection tell us about early galaxies?
• Design and characterization of the Large-Aperture Experiment to Detect the Dark Age (LEDA) radiometer systems
• Tighter Limits on Dark Matter Explanations of the Anomalous EDGES 21cm Signal
• Probing the Dark Ages at Z~20: The SCI-HI 21 cm All-Sky Spectrum Experiment
• Testing for calibration systematics in the EDGES low-band data using Bayesian model selection
• BIGHORNS - Broadband Instrument for Global HydrOgen ReioNisation Signal
• The Radio Background Below 100 MHz
• The redshifted 21-cm signal in the EDGES low-band spectrum
• The Rise of the First Stars: Supersonic Streaming, Radiative Feedback, and 21-cm Cosmology
• SARAS 2 constraints on global 21-cm signals from the Epoch of Reionization
• Decoding the X-ray Properties of Pre-Reionization Era Sources
• SARAS 2: A Spectral Radiometer for probing Cosmic Dawn and the Epoch of Reionization through detection of the global 21 cm signal
• Foreground Model and Antenna Calibration Errors in the Measurement of the Sky-Averaged λ21 cm Signal at z~20
• A General Bayesian Framework for Foreground Modelling and Chromaticity Correction for Global 21cm Experiments
• Emulating the Global 21-cm Signal from Cosmic Dawn and Reionization
• An MCMC approach to extracting the global 21-cm signal during the cosmic dawn from sky-averaged radio observations
• Results from EDGES High-Band: III. New Constraints on Parameters of the Early Universe
• On the detection of spectral ripples from the Recombination Epoch
• Radio background and IGM heating due to Pop III supernovae explosions
• HI Epoch of Reionization Arrays
• Modeling the Radio Foreground for detection of CMB spectral distortions from Cosmic Dawn and Epoch of Reionization

Cited by in corpus (39)

• The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28
• A General Bayesian Framework for Foreground Modelling and Chromaticity Correction for Global 21cm Experiments
• GLOBALEMU: A novel and robust approach for emulating the sky-averaged 21-cm signal from the cosmic dawn and epoch of reionisation
• Quantifying Ionospheric Effects on Global 21-cm Observations
• Elevated UV luminosity density at Cosmic Dawn explained by non-evolving, weakly mass-dependent star formation efficiency
• Antenna beam characterisation for the global 21cm experiment LEDA and its impact on signal model parameter reconstruction
• Informing antenna design for sky-averaged 21-cm experiments using a simulated Bayesian data analysis pipeline
• Use of Time Dependent Data in Bayesian Global 21cm Foreground and Signal Modelling
• Array Element Coupling in Radio Interferometry I: A Semi-Analytic Approach
• A Bayesian approach to modelling spectrometer data chromaticity corrected using beam factors -- I. Mathematical formalism
• Exploring the cosmic dawn and epoch of reionization with 21cm line
• A General Bayesian Framework to Account for Foreground Map Errors in Global 21-cm Experiments
• Sky-averaged 21-cm signal extraction using multiple antennas with an SVD framework: the REACH case
• Bayesian Data Analysis for Sky-averaged 21-cm Experiments in the Presence of Ionospheric Effects
• Extracting the Global 21-cm signal from Cosmic Dawn and Epoch of Reionization in the presence of Foreground and Ionosphere
• The effects of the antenna power pattern uncertainty within a global 21 cm experiment
• Bayesian evidence-driven diagnosis of instrumental systematics for sky-averaged 21-cm cosmology experiments
• Joint constraints on reionization: a framework for combining the global 21cm signal and the kinetic Sunyaev-Zel'dovich effect
• Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
• Bayesian evidence-driven likelihood selection for sky-averaged 21-cm signal extraction
• Direction-dependent effects on global 21-cm detection
• Rapid and Late Cosmic Reionization Driven by Massive Galaxies: a Joint Analysis of Constraints from 21-cm, Lyman Line & CMB Data Sets
• Comparing sampling techniques to chart parameter space of 21 cm Global signal with Artificial Neural Networks
• Hunting for bumps in the margins
• Improvements to the Search for Cosmic Dawn Using the Long Wavelength Array
• Global signal in the redshifted hydrogen 21-cm line from the Dark Ages and Cosmic Dawn: Dependence on the nature of dark matter and modeling of first light
• A general Bayesian model-validation framework based on null-test evidence ratios, with an example application to global 21-cm cosmology
• A Bayesian Method to Mitigate the Effects of Unmodelled Time-Varying Systematics for 21-cm Cosmology Experiments
• Probing the global 21-cm background by velocity-induced dipole and quadrupole anisotropies
• A Bayesian approach to modelling spectrometer data chromaticity corrected using beam factors -- II. Model priors and posterior odds
• Forecast measurement of the 21 cm global spectrum from Lunar orbit with the Vari-Zeroth-Order Polynomial (VZOP) method
• Exploring Machine Learning Regression Models for Advancing Foreground Mitigation and Global 21cm Signal Parameter Extraction
• Beamforming approaches toward detecting the 21-cm global signal from Cosmic Dawn with radio array telescopes
• The magnitude of the dark ages 21-cm signal in the context of existing early and late time constraints on $\Lambda$CDM
• Recovering 21cm Monopole Signals Without Smoothness
• Calibrating Bayesian Tension Statistics using Neural Ratio Estimation
• On the accuracy of posterior recovery with neural network emulators
• A Comprehensive Bayesian re-analysis of the SARAS2 data from the Epoch of Reionization
• Optimizing Foreground Modelling for Global 21cm Cosmology with GPU-Accelerated Nested Sampling