Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology
Accurate detection of the cosmological 21 cm global signal requires galactic foreground models that can remove power over 10 ^6 . Although foreground and global signal models unavoidably exhibit overlap in their vector spaces inducing bias error in the extracted signal, a second source of bias and e...
Main Authors: | , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
IOP Publishing
2023-01-01
|
Series: | The Astrophysical Journal |
Subjects: | |
Online Access: | https://doi.org/10.3847/1538-4357/ad06b3 |
_version_ | 1797393602056290304 |
---|---|
author | Joshua J. Hibbard David Rapetti Jack O. Burns Nivedita Mahesh Neil Bassett |
author_facet | Joshua J. Hibbard David Rapetti Jack O. Burns Nivedita Mahesh Neil Bassett |
author_sort | Joshua J. Hibbard |
collection | DOAJ |
description | Accurate detection of the cosmological 21 cm global signal requires galactic foreground models that can remove power over 10 ^6 . Although foreground and global signal models unavoidably exhibit overlap in their vector spaces inducing bias error in the extracted signal, a second source of bias and error arises from inadequate foreground models, i.e., models that cannot fit spectra down to the noise level of the signal. We therefore test the level to which seven commonly employed foreground models—including nonlinear and linear forward models, polynomials, and maximally smooth polynomials—fit realistic simulated mock foreground spectra, as well as their dependence upon model inputs. The mock spectra are synthesized for an EDGES-like experiment and we compare all models’ goodness of fit and preference using a Kolmogorov–Smirnov (K-S) test of the noise-normalized residuals in order to compare models with differing, and sometimes indeterminable, degrees of freedom. For a single local sidereal time (LST) bin spectrum and p -value threshold of p = 0.05, the nonlinear forward model with four parameters is preferred ( p = 0.99), while the linear forward model fits well with six to seven parameters ( p = 0.94, 0.97, respectively). The polynomials and maximally smooth polynomials, like those employed by the EDGES and SARAS3 experiments, cannot produce good fits with five parameters for the experimental simulations in this work ( p < 10 ^−6 ). However, we find that polynomials with six parameters pass the K-S test ( p = 0.4), although a nine-parameter fit produces the highest p -value ( p ∼ 0.67). When fitting multiple LST bins simultaneously, we find that the linear forward model outperforms (a higher p -value) the nonlinear model for 2, 5, and 10 LST bins. Importantly, the K-S test consistently identifies best-fit and preferred models. |
first_indexed | 2024-03-09T00:04:58Z |
format | Article |
id | doaj.art-5e7b8bf96dd5495db6b3a00887509fe8 |
institution | Directory Open Access Journal |
issn | 1538-4357 |
language | English |
last_indexed | 2024-03-09T00:04:58Z |
publishDate | 2023-01-01 |
publisher | IOP Publishing |
record_format | Article |
series | The Astrophysical Journal |
spelling | doaj.art-5e7b8bf96dd5495db6b3a00887509fe82023-12-12T13:26:30ZengIOP PublishingThe Astrophysical Journal1538-43572023-01-01959210310.3847/1538-4357/ad06b3Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm CosmologyJoshua J. Hibbard0https://orcid.org/0000-0002-9377-5133David Rapetti1https://orcid.org/0000-0003-2196-6675Jack O. Burns2https://orcid.org/0000-0002-4468-2117Nivedita Mahesh3https://orcid.org/0000-0003-2560-8023Neil Bassett4https://orcid.org/0000-0001-7051-6385Center for Astrophysics and Space Astronomy, Department of Astrophysical and Planetary Science, University of Colorado Boulder , CO 80309, USA ; Joshua.Hibbard@colorado.eduCenter for Astrophysics and Space Astronomy, Department of Astrophysical and Planetary Science, University of Colorado Boulder , CO 80309, USA ; Joshua.Hibbard@colorado.edu; NASA Ames Research Center , Moffett Field, CA 94035, USA; Research Institute for Advanced Computer Science, Universities Space Research Association , Washington, DC 20024, USACenter for Astrophysics and Space Astronomy, Department of Astrophysical and Planetary Science, University of Colorado Boulder , CO 80309, USA ; Joshua.Hibbard@colorado.eduCalifornia Institute of Technology , Pasadena, CA 91125, USACenter for Astrophysics and Space Astronomy, Department of Astrophysical and Planetary Science, University of Colorado Boulder , CO 80309, USA ; Joshua.Hibbard@colorado.eduAccurate detection of the cosmological 21 cm global signal requires galactic foreground models that can remove power over 10 ^6 . Although foreground and global signal models unavoidably exhibit overlap in their vector spaces inducing bias error in the extracted signal, a second source of bias and error arises from inadequate foreground models, i.e., models that cannot fit spectra down to the noise level of the signal. We therefore test the level to which seven commonly employed foreground models—including nonlinear and linear forward models, polynomials, and maximally smooth polynomials—fit realistic simulated mock foreground spectra, as well as their dependence upon model inputs. The mock spectra are synthesized for an EDGES-like experiment and we compare all models’ goodness of fit and preference using a Kolmogorov–Smirnov (K-S) test of the noise-normalized residuals in order to compare models with differing, and sometimes indeterminable, degrees of freedom. For a single local sidereal time (LST) bin spectrum and p -value threshold of p = 0.05, the nonlinear forward model with four parameters is preferred ( p = 0.99), while the linear forward model fits well with six to seven parameters ( p = 0.94, 0.97, respectively). The polynomials and maximally smooth polynomials, like those employed by the EDGES and SARAS3 experiments, cannot produce good fits with five parameters for the experimental simulations in this work ( p < 10 ^−6 ). However, we find that polynomials with six parameters pass the K-S test ( p = 0.4), although a nine-parameter fit produces the highest p -value ( p ∼ 0.67). When fitting multiple LST bins simultaneously, we find that the linear forward model outperforms (a higher p -value) the nonlinear model for 2, 5, and 10 LST bins. Importantly, the K-S test consistently identifies best-fit and preferred models.https://doi.org/10.3847/1538-4357/ad06b3CosmologyObservational cosmologyH I line emissionReionizationInterstellar synchrotron emissionGalactic radio sources |
spellingShingle | Joshua J. Hibbard David Rapetti Jack O. Burns Nivedita Mahesh Neil Bassett Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology The Astrophysical Journal Cosmology Observational cosmology H I line emission Reionization Interstellar synchrotron emission Galactic radio sources |
title | Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology |
title_full | Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology |
title_fullStr | Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology |
title_full_unstemmed | Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology |
title_short | Fitting and Comparing Galactic Foreground Models for Unbiased 21 cm Cosmology |
title_sort | fitting and comparing galactic foreground models for unbiased 21 cm cosmology |
topic | Cosmology Observational cosmology H I line emission Reionization Interstellar synchrotron emission Galactic radio sources |
url | https://doi.org/10.3847/1538-4357/ad06b3 |
work_keys_str_mv | AT joshuajhibbard fittingandcomparinggalacticforegroundmodelsforunbiased21cmcosmology AT davidrapetti fittingandcomparinggalacticforegroundmodelsforunbiased21cmcosmology AT jackoburns fittingandcomparinggalacticforegroundmodelsforunbiased21cmcosmology AT niveditamahesh fittingandcomparinggalacticforegroundmodelsforunbiased21cmcosmology AT neilbassett fittingandcomparinggalacticforegroundmodelsforunbiased21cmcosmology |