Astrophysical Foregrounds Research Articles

Gravitational waves from cosmic strings in LISA: reconstruction pipeline and physics interpretation

We initiate the LISA template databank for stochastic gravitational wave backgrounds sourced by cosmic strings. We include two templates, an analytical template, which enables more flexible searches, and a numerical template derived directly from large Nambu-Goto simulations of string networks. Using searches based on these templates, we forecast the parameter space within the reach of the experiment and the precision with which their parameters will be reconstructed, provided a signal is observed. The reconstruction permits probing the Hubble expansion and new relativistic DoF in the early universe. We quantify the impact that astrophysical foregrounds can have on these searches. Finally, we discuss the impact that these observations would have on our understanding of the fundamental models behind the string networks. Overall, we prove that LISA has great potential for probing cosmic string models and may reach tensions as low as Gμ = 10-16 – 10-17, which translates into energy scales of the order 1011 GeV.

Foreground Removal and Angular Power Spectrum Estimation of 21 cm Signal Using Harmonic Space ILC Method

Abstract Mapping the distribution of neutral atomic hydrogen (H i) in the Universe through its 21 cm emission line provides a powerful cosmological probe to map the large-scale structures and shed light on various cosmological phenomena. The baryon acoustic oscillations at low redshifts can potentially be probed by sensitive H i intensity mapping experiments and constrain the properties of dark energy. However, the 21 cm signal detection faces formidable challenges owing to the dominance of various astrophysical foregrounds, which can be several orders of magnitude stronger. Our current work introduces a novel and model-independent internal linear combination (ILC) method in harmonic space using the principal components of the 21 cm signal for accurate foreground removal and power spectrum estimation. We estimate the principal components by incorporating prior knowledge of the theoretical 21 cm covariance matrix. We test our methodology by detailed simulations of radio observations, incorporating synchrotron emission, free–free radiation, extragalactic point sources, and thermal noise. We estimate the full-sky 21 cm angular power spectrum after application of a mask on the full-sky cleaned 21 cm signal by using the mode–mode coupling matrix. These full-sky estimates of angular spectra can be directly used to measure the cosmological parameters. For the first time, we demonstrate the effectiveness of a foreground-model-independent ILC method in harmonic space to reconstruct the 21 cm signal.

A generative modeling approach to reconstructing 21 cm tomographic data

Abstract Analyses of the cosmic 21 cm signal are hampered by astrophysical foregrounds that are far stronger than the signal itself. These foregrounds, typically confined to a wedge-shaped region in Fourier space, often necessitate the removal of a vast majority of modes, thereby degrading the quality of the data anisotropically. To address this challenge, we introduce a novel deep generative model based on stochastic interpolants to reconstruct the 21 cm data lost to wedge filtering. Our method leverages the non-Gaussian nature of the 21 cm signal to effectively map wedge-filtered 3D lightcones to samples from the conditional distribution of wedge-recovered lightcones. We demonstrate how our method is able to restore spatial information effectively, considering both varying cosmological initial conditions and astrophysics. Furthermore, we discuss a number of future avenues where this approach could be applied in analyses of the 21 cm signal, potentially offering new opportunities to improve our understanding of the Universe during the epochs of cosmic dawn and reionization. Code, pre-trained models, and scripts for making plots in this paper can be found here.

Primordial gravitational wave backgrounds from phase transitions with next generation ground based detectors

Abstract Third generation ground-based gravitational wave (GW) detectors, such as Einstein Telescope and Cosmic Explorer, will operate in the (few-104) Hz frequency band, with a boost in sensitivity providing an unprecedented reach into primordial cosmology. 
Working concurrently with pulsar timing arrays in the nHz band, and LISA in the mHz band, these 3G detectors will be powerful probes of beyond the standard model particle physics on scales T ∼106 GeV and above. Here we focus on their ability to probe phase transitions (PTs) in the early universe. We first overview the landscape of detectors across frequencies, discuss the relevance of astrophysical foregrounds, and provide convenient and up-to-date power-law integrated sensitivity curves for these detectors. We then present the constraints expected from GW observations on both first order PTs, and on topological defects that may be formed when a symmetry is broken irrespective of the order of the phase transition (strings and domain walls). These constraints can then be applied to specific models leading to first order PTs and/or topological defects. In particular we discuss the implications for axion models, which solve the strong CP problem by introducing a spontaneously broken Peccei-Quinn (PQ) symmetry. For post-inflationary breaking, the PQ scale must lie in the 108-1011 GeV range, and so the signal from a first order PQ PT falls within reach of ground based 3G detectors. A scan in parameter space of signal-to-noise ratio in a representative model reveals their large potential to probe the nature of the PQ transition. Additionally, in heavy axion type models domain walls form, which can lead to a detectable GW background. We discuss their spectrum and summarise the expected constraints on these models from 3G detectors, together with SKA and LISA.

Improving the detection sensitivity to primordial stochastic gravitational waves with reduced astrophysical foregrounds. II. Subthreshold binary neutron stars

Improving the detection sensitivity to primordial stochastic gravitational waves with reduced astrophysical foregrounds. II. Subthreshold binary neutron stars

Forward modeling fluctuations in the DESI LRGs target sample using image simulations

We usethe forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.

First Constraints on the Epoch of Reionization Using the Non-Gaussianity of the Kinematic Sunyaev-Zel'dovich Effect from the South Pole Telescope and Herschel-SPIRE Observations.

We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel'dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and Herschel-SPIRE experiments over a 100 deg^{2} field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 μK-arcmin in bands centered at 95, 150, and 220GHz. For SPIRE, we include data from the 600 and 857GHz bands. We reconstruct the velocity-induced large-scale correlation of the small-scale kSZ signal with a quadratic estimator that uses two cosmic microwave background (CMB) temperature maps, constructed by optimally combining data from all the frequency bands. We reject the null hypothesis of a zero trispectrum at 10.3σ level. However, the measured trispectrum contains contributions from both the kSZ and other undesired components, such as CMB lensing and astrophysical foregrounds, with kSZ being sub-dominant. We use the agora simulations to estimate the expected signal from CMB lensing and astrophysical foregrounds. After accounting for the contributions from CMB lensing and foreground signals, we do not detect an excess kSZ-only trispectrum and use this nondetection to set constraints on reionization. By applying a prior based on observations of the Gunn-Peterson trough, we obtain an upper limit on the duration of reionization of Δz_{re,50}<4.5 (95% confidence level). We find these constraints are fairly robust to foregrounds assumptions. This trispectrum measurement is independent of, but consistent with, Planck's optical depth measurement. This result is the first constraint on the epoch of reionization using the non-Gaussian nature of the kSZ signal.

Impact of extragalactic point sources on the low-frequency sky spectrum and cosmic dawn global 21-cm measurements

Abstract Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements of the luminosity function and the angular correlation function of extragalactic point sources, we develop a model for the contribution of these sources to the low-frequency sky spectrum. Using this model, we find that the contribution of sources with flux density &amp;gt;10−6 Jy to the sky-averaged spectrum is smooth and of the order of a few kelvins at 50–200 MHz. We combine this model with measurements of the galactic foreground spectrum and weigh the resultant sky by the beam directivity of the conical log-spiral antenna planned as part of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) project. We find that the contribution of point sources to the resultant spectrum is $\sim 0.4\%$ of the total foregrounds, but still larger by at least an order of magnitude than the standard predictions for the cosmological 21-cm signal. As a result, not accounting for the point-source contribution leads to a systematic bias in 21-cm signal recovery. We show, however, that in the REACH case, this reconstruction bias can be removed by modelling the point-source contribution as a power law with a running spectral index. We make our code publicly available as a Python package labelled epspy.

The Simons Observatory: impact of bandpass, polarization angle and calibration uncertainties on small-scale power spectrum analysis

We study the effects due to mismatches in passbands, polarization angles, and temperature and polarization calibrations in the context of the upcoming cosmic microwave background experiment Simons Observatory (SO). Using the SO multi-frequency likelihood, we estimate the bias and the degradation of constraining power in cosmological and astrophysical foreground parameters assuming different levels of knowledge of the instrumental effects. We find that incorrect but reasonable assumptions about the values ofall the systematics examined here can have significant effects on cosmologicalanalyses, hence requiring marginalization approaches at the likelihood level.When doing so, we find that the most relevant effect is due to bandpass shifts. When marginalizing over them, the posteriors of parameters describing astrophysical microwave foregrounds (such as radio point sources or dust) get degraded, while cosmological parameters constraints are not significantly affected.Marginalization over polarization angles with up to 0.25° uncertainty causes an irrelevant bias ≲ 0.05 σ in all parameters.Marginalization over calibration factors in polarization broadens the constraints on the effective number of relativistic degrees of freedom Neff by a factor 1.2, interpreted here as a proxy parameter for non standard model physics targeted by high-resolution CMB measurements.

Primordial black holes from conformal Higgs

Scale-invariant extensions of the electroweak theory are not only attractive because they can dynamically generate the weak scale, but also due to their role in facilitating supercooled first-order phase transitions. We study the minimal scale-invariant U(1)D extension of the standard model and show that Primordial Black Holes (PBHs) can be abundantly produced. The mass of these PBHs is bounded from above by that of the moon due to QCD catalysis limiting the amount of supercooling. Lunar-mass PBHs, which are produced for dark Higgs vev vϕ≃20TeV, correspond to the best likelihood to explain the HSC lensing anomaly. For vϕ≳400TeV, the model can explain hundred per cent of dark matter. At even larger hierarchy of scales, it can contribute to the 511keV line. While the gravitational wave (GW) signal produced by the HSC anomaly interpretation is large and detectable by LISA above astrophysical foreground, the dark matter interpretation in terms of PBHs can not be entirely probed by future GW detection. This is due to the dilution of the signal by the entropy injected during the decay of the long-lived U(1)D scalar. This extended lifetime is a natural consequence of the large hierarchy of scales.

Identifying frequency de-correlated dust residuals in B-mode maps by exploiting the spectral capability of bolometric interferometry

Context. Astrophysical polarized foregrounds represent the most critical challenge in cosmic microwave background (CMB) B-mode experiments, requiring multifrequency observations to constrain astrophysical foregrounds and isolate the CMB signal. However, recent observations indicate that foreground emission may be more complex than anticipated. Not properly accounting for these complexities during component separation can lead to a bias in the recovered tensor-to-scalar ratio. Aims. In this paper we investigate how the increased spectral resolution provided by band-splitting in bolometric interferometry (BI) through a technique called spectral imaging can help control the foreground contamination in the case of an unaccounted-for Galactic dust frequency de-correlation along the line of sight (LOS). Methods. We focused on the next-generation ground-based CMB experiment CMB-S4 and compared its anticipated sensitivity, frequency, and sky coverage with a hypothetical version of the same experiment based on BI (CMB-S4/BI). We performed a Monte Carlo analysis based on parametric component separation methods (FGBuster and Commander) and computed the likelihood of the recovered tensor-to-scalar ratio, r. Results. The main result is that spectral imaging allows us to detect systematic uncertainties on r from frequency de-correlation when this effect is not accounted for in the component separation. Conversely, an imager such as CMB-S4 would detect a biased value of r and would be unable to spot the presence of a systematic effect. We find a similar result in the reconstruction of the dust spectral index, and we show that with BI we can more precisely measure the dust spectral index when frequency de-correlation is present and not accounted for in the component separation. Conclusions. The in-band frequency resolution provided by BI allows us to identify dust LOS frequency de-correlation residuals where an imager with a similar level of performance would fail. This creates the possibility of exploiting this potential in the context of future CMB polarization experiments that will be challenged by complex foregrounds in their quest for B-mode detection.

Fully non-Gaussian Scalar-Induced Gravitational Waves

Scalar-Induced Gravitational Waves (SIGWs) represent a particular class of primordial signals which are sourced at second-order in perturbation theory whenever a scalar fluctuation of the metric is present. They form a guaranteed Stochastic Gravitational Wave Background (SGWB) that, depending on the amplification of primordial scalar fluctuations, can be detected by GW detectors. The amplitude and the frequency shape of the scalar-induced SGWB can be influenced by the statistical properties of the scalar density perturbations. In this work we study the intuitive physics behind SIGWs and we analyze the imprints of local non-Gaussianity of the primordial curvature perturbation on the GW spectrum. We consider all the relevant non-Gaussian contributions up to fifth-order in the scalar seeds without any hierarchy, and we derive the related GW energy density ΩGW(f). We perform a Fisher matrix analysis to understand to which accuracy non-Gaussianity can be constrained with the LISA detector, which will be sensitive in the milli-Hertz frequency band. We find that LISA, neglecting the impact of astrophysical foregrounds, will be able to measure the amplitude, the width and the peak of the spectrum with an accuracy up to 𝒪(10-4), while non-Gaussianity can be measured up to 𝒪(10-3). Finally, we discuss the implications of our non-Gaussianity expansion on the fraction of Primordial Black Holes.

Strategy for mitigation of systematics for EoR experiments with the Murchison Widefield Array

Abstract Observations of the 21 cm signal face significant challenges due to bright astrophysical foregrounds that are several orders of magnitude higher than the brightness of the hydrogen line, along with various systematics. Successful 21 cm experiments require accurate calibration and foreground mitigation. Errors introduced during the calibration process such as systematics can disrupt the intrinsic frequency smoothness of the foregrounds, leading to power leakage into the Epoch of Reionisation window. Therefore, it is essential to develop strategies to effectively address these challenges. In this work, we adopt a stringent approach to identify and address suspected systematics, including malfunctioning antennas, frequency channels corrupted by radio frequency interference, and other dominant effects. We implement a statistical framework that utilises various data products from the data processing pipeline to derive specific criteria and filters. These criteria and filters are applied at intermediate stages to mitigate systematic propagation from the early stages of data processing. Our analysis focuses on observations from the Murchison Widefield Array Phase I configuration. Out of the observations processed by the pipeline, our approach selects 18%, totalling 58 h, that exhibit fewer systematic effects. The successful selection of observations with reduced systematic dominance enhances our confidence in achieving 21 cm measurements.

Retrieving the 21-cm signal from the Epoch of Reionization with learnt Gaussian process kernels

ABSTRACT Direct detection of the Cosmic Dawn and Epoch of Reionization via the redshifted 21-cm line of neutral Hydrogen will have unprecedented implications for studying structure formation in the early Universe. This exciting goal is challenged by the difficulty of extracting the faint 21-cm signal buried beneath bright astrophysical foregrounds and contaminated by numerous systematics. Here, we focus on improving the Gaussian Process Regression (GPR) signal separation method originally developed for LOFAR observations. We address a key limitation of the current approach by incorporating covariance prior models learnt from 21-cm signal simulations using variational autoencoder (VAE) and interpolatory autoencoder (IAE). Extensive tests are conducted to evaluate GPR, VAE–GPR, and IAE–GPR in different scenarios. Our findings reveal that the new method outperforms standard GPR in component separation tasks. Moreover, the improved method demonstrates robustness when applied to signals not represented in the training set. It also presents a certain degree of resilience to data systematics, highlighting its ability to effectively mitigate their impact on the signal recovery process. However, our findings also underscore the importance of accurately characterizing and understanding these systematics to achieve successful detection. Our generative approaches provide good results even with limited training data, offering a valuable advantage when a large training set is not feasible. Comparing the two algorithms, IAE–GPR shows slightly higher fidelity in recovering power spectra compared to VAE–GPR. These advancements highlight the strength of generative approaches and optimize the analysis techniques for future 21-cm signal detection at high redshifts.

Cross Correlation of Pencil-beam Galaxy Surveys and Line-intensity Maps: An Application of the James Webb Space Telescope

Line-intensity mapping (IM) experiments seek to perform statistical measurements of large-scale structure with spectral lines such as 21 cm, CO, and Lyα. A challenge in these observations is to ensure that astrophysical foregrounds, such as galactic synchrotron emission in 21 cm measurements, are properly removed. One method that has the potential to reduce foreground contamination is to cross correlate with a galaxy survey that overlaps with the IM volume. However, telescopes sensitive to high-redshift galaxies typically have small field of views compared to IM surveys. Thus, a galaxy survey for cross correlation would necessarily consist of pencil beams that sparsely fill the IM volume. In this paper, we develop the formalism to forecast the sensitivity of cross correlations between IM experiments and pencil-beam galaxy surveys. We find that a random distribution of pencil beams leads to very similar overall sensitivity as a lattice spaced across the IM survey and derive a simple formula for random configurations that agrees with the Fisher matrix formalism. We explore examples of combining high-redshift James Webb Space Telescope (JWST) observations with both an SPHEREx-like Lyα IM survey and a 21 cm experiment based on the Hydrogen Epoch of Reionization Array (HERA). We find that the JWST-SPHEREx case is promising, leading to a total signal-to-noise ratio of ∼5 after 100 total hours of JWST (at z = 7). We find that HERA is not well-suited for this approach owing to its drift-scan strategy, but that a similar experiment that can integrate down on one field could be.

Improving the detection sensitivity to primordial stochastic gravitational waves with reduced astrophysical foregrounds

One of the primary targets of third-generation (3G) ground-based gravitational wave (GW) detectors is detecting the stochastic GW background (SGWB) from early universe processes. The astrophysical foreground from compact binary mergers will be a major contamination to the background, which must be reduced to high precision to enable the detection of primordial background. In this work, we revisit the limit of foreground reduction computed in previous studies, point out potential problems in previous foreground cleaning methods and propose a novel cleaning method subtracting the approximate signal strain and removing the average residual power. With this method, the binary black hole foreground is reduced with fractional residual energy density below 10 −4 for frequency f ∈ (10, 10 2 ) Hz, below 10 −3 for frequency f ∈ (10 2 , 10 3 ) Hz and below the detector sensitivity limit for all relevant frequencies in our simulations. Similar precision is achieved to clean the foreground from binary neutron stars (BNSs) that are above the detection threshold, so that the residual foreground is dominated by sub-threshold BNSs, which will be the next critical problem to solve for detecting the primordial SGWB in the 3G era.

Precision cosmology with primordial GW backgrounds in presence of astrophysical foregrounds

The era of Gravitational-Wave (GW) astronomy will grant the detection of the astrophysical GW background from unresolved mergers of binary black holes, and the prospect of probing the presence of primordial GW backgrounds. In particular, the low-frequency tail of the GW spectrum for causally-generated primordial signals (like a phase transition) offers an excellent opportunity to measure unambiguously cosmological parameters as the equation of state of the universe, or free-streaming particles at epochs well before recombination. We discuss whether this programme is jeopardised by the uncertainties on the astrophysical GW foregrounds that coexist with a primordial background. We detail the motivated assumptions under which the astrophysical foregrounds can be assumed to be known in shape, and only uncertain in their normalisation. In this case, the sensitivity to a primordial signal can be computed by a simple and numerically agile procedure, where the optimal filter function subtracts the components of the astrophysical foreground that are close in spectral shape to the signal. We show that the degradation of the sensitivity to the signal in presence of astrophysical foregrounds is limited to a factor of a few, and only around the frequencies where the signal is closer to the foregrounds. Our results highlight the importance of modelling the contributions of eccentric or intermediate-mass black hole binaries to the GW background, to consolidate the prospects to perform precision cosmology with primordial GW backgrounds.

Detecting non-Gaussian gravitational wave backgrounds: A unified framework

We describe a novel approach to the detection and parameter estimation of a non-Gaussian stochastic background of gravitational waves. The method is based on the determination of relevant statistical parameters using importance sampling. We show that it is possible to improve the Gaussian detection statistics by simulating realizations of the expected signal for a given model. While computationally expensive, our method improves the detection performance, leveraging the prior knowledge on the expected signal, and can be used in a natural way to extract physical information about the background. We present the basic principles of our approach, characterize the detection statistic performances in a simplified context, and discuss possible applications to the detection of some astrophysical foregrounds. We argue that the proposed approach, complementarily to the ones available in literature might be used to detect suitable astrophysical foregrounds by currently operating and future gravitational wave detectors.

Mock data study for next-generation ground-based detectors: The performance loss of matched filtering due to correlated confusion noise

The next-generation (3G/XG) ground-based gravitational-wave (GW) detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE) will begin observing in the next decade. Due to the extremely high sensitivity of these detectors, the majority of stellar-mass compact-binary mergers in the entire Universe will be observed. It is also expected that 3G detectors will have significant sensitivity down to 2-7 Hz; the observed duration of binary neutron star signals could increase to several hours or days. The abundance and duration of signals will cause them to overlap in time, which may form a confusion noise that could affect the detection of individual GW sources when using naive matched filtering; matched filtering is only optimal for stationary Gaussian noise. We create mock data for CE and ET using the latest population models informed by the GWTC-3 catalog and investigate the performance loss of matched filtering due to overlapping signals. We find the performance loss mainly comes from a deviation in the noise's measured amplitude spectral density. The redshift reach of CE (ET) can be reduced by 15%-38% (8%-21%) depending on the merger rate estimate. The direct contribution of confusion noise to the total signal-to-noise ratio (SNR) is generally negligible compared to the contribution from instrumental noise. We also find that correlated confusion noise has a negligible effect on the quadrature summation rule of network SNR for ET, but might reduce the network SNR of high detector-frame mass signals for detector networks including CE if no mitigation is applied. For ET, the null stream can mitigate the astrophysical foreground. For CE, we demonstrate that a computationally efficient, straightforward single-detector signal subtraction method suppresses the total noise to almost the instrument noise level; this will allow for near-optimal searches.

A Measurement of the Cosmic Optical Background and Diffuse Galactic Light Scaling from the R < 50 au New Horizons-LORRI Data

Direct photometric measurements of the cosmic optical background (COB) provide an important point of comparison to both other measurement methodologies and models of cosmic structure formation, and permit a cosmic consistency test with the potential to reveal additional diffuse sources of emission. The COB has been challenging to measure from Earth due to the difficulty of isolating it from the diffuse light scattered from interplanetary dust in our solar system. We present a measurement of the COB using data taken by the Long-Range Reconnaissance Imager on NASA's New Horizons mission, considering all data acquired to 47 au. We employ a blind methodology where our analysis choices are developed against a subset of the full data set, which is then unblinded. Dark current and other instrumental systematics are accounted for, including a number of sources of scattered light. We fully characterize and remove structured and diffuse astrophysical foregrounds including bright stars, the integrated starlight from faint unresolved sources, and diffuse galactic light. For the full data set, we find the surface brightness of the COB to be nW m−2 sr−1. This result supports recent determinations that find a factor of 2–3× more light than expected from the integrated light from galaxies and motivate new diffuse intensity measurements with more capable instruments that can support spectral measurements over the optical and near-IR.