Astrophysics

Up-to-date cosmological data analyses have shown that \textit{(a)} a closed universe is preferred by the Planck data at more than 99% CL, and \textit{(b)} interacting scenarios offer a very compelling solution to the Hubble constant tension. In light of these two recent appealing scenarios, we consider here an interacting dark matter-dark energy model with a non-zero spatial curvature component and a freely varying dark energy equation of state in both the quintessential and phantom regimes. When considering Cosmic Microwave Background data only, a phantom and closed universe can perfectly alleviate the Hubble tension, without the necessity of a coupling among the dark sectors. Accounting for other possible cosmological observations compromises the viability of this very attractive scenario as a global solution to current cosmological tensions, either by spoiling its effectiveness concerning the H 0 problem, as in the case of Supernovae Ia data, or by introducing a strong disagreement in the preferred value of the spatial curvature, as in the case of Baryon Acoustic Oscillations.

We observed 51 sources in the Q-U-I JOint TEnerife (QUIJOTE) cosmological fields which were brighter than 1 Jy at 30 GHz in the Planck Point Source Catalogue (version 1), with the Very Large Array at 28 -- 40 GHz, in order to characterise their high-radio-frequency variability and polarization properties. We find a roughly log-normal distribution of polarization fractions with a median of 2%, in agreement with previous studies, and a median rotation measure (RM) of ??1110 rad m ?? with one outlier up to ??64000 rad m ?? which is among the highest RMs measured in quasar cores. We find hints of a correlation between the total intensity flux density and median polarization fraction. We find 59% of sources are variable in total intensity, and 100% in polarization at 3? level, with no apparent correlation between total intensity variability and polarization variability. This indicates that it will be difficult to model these sources without simultaneous polarimetric monitoring observations and they will need to be masked for cosmological analysis.

The Large-aperture Experiment to detect the Dark Age (LEDA) was designed to measure the 21-cm signal from neutral hydrogen at Cosmic Dawn, z??15-30. Using observations made with the ??200 m diameter core of the Owens Valley Long Wavelength Array (OVRO-LWA), we present a 2-D cylindrical spatial power spectrum for data at 43.1-53.5 MHz ( z median ??8 ) incoherently integrated for 4 hours, and an analysis of the array sensitivity. Power from foregrounds is localized to a "wedge" within k ??, k ??space. After calibration of visibilities using 5 bright compact sources including VirA, we measure ? 2 (k)??? 10 12 mK 2 outside the foreground wedge, where an uncontaminated cosmological signal would lie, in principle. The measured ? 2 (k) is an upper limit that reflects a combination of thermal instrumental and sky noise, and unmodelled systematics that scatter power from the wedge, as will be discussed. By differencing calibrated visibilities for close pairs of frequency channels, we suppress foreground sky structure and systematics, extract thermal noise, and use a mix of coherent and incoherent integration to simulate a noise-dominated power spectrum for a 3000 h observation and z= 16-37. For suitable calibration quality, the resulting noise level, ? 2 (k)??00 mK 2 (k = 0.3 Mpc ?? ), would be sufficient to detect peaks in the 21-cm spatial power spectrum due to early Ly- α and X-ray sources, as predicted for a range of theoretical model parameters.

The Cosmic Microwave Background (CMB) radiation B mode polarization signal contains the unique signature of primordial metric perturbations produced during the inflation. The separation of the weak CMB B-mode signal from strong foreground contamination in observed maps is a complex task, and proposed new generation low noise satellite missions compete with the weak signal level of this gravitational background. In this article, for the first time, we employ a foreground model-independent internal linear combination (ILC) method to reconstruct the CMB B mode signal using simulated observations over large angular scales of the sky of 6 frequency bands of future generation CMB mission Probe of Inflation and Cosmic Origins (PICO). We estimate the joint CMB B mode posterior density following the interleaving Gibbs steps of B mode angular power spectrum and cleaned map samples using the ILC method. We extend and improve the earlier reported Bayesian ILC method to analyze weak CMB B mode reconstruction by introducing noise bias corrections at two stages during the ILC weight estimation. By performing 200 Monte Carlo simulations of the Bayesian ILC method, we find that our method can reconstruct the CMB signals and the joint posterior density accurately over large angular scales of the sky. We estimate Blackwell-Rao statistics of the marginal density of CMB B mode angular power spectrum and use them to estimate the joint density of scalar to tensor ratio r and a lensing power spectrum amplitude A lens . Using 200 Monte Carlo simulations of the delensing approach, we find that our method can achieve an unbiased detection of the primordial gravitational wave signal r with more than 8 σ significance for levels of r⩾0.01 .

The estimation of the polarization P of extragalactic compact sources in Cosmic Microwave Background images is a very important task in order to clean these images for cosmological purposes -- as, for example, to constrain the tensor-to-scalar ratio of primordial fluctuations during inflation -- and also to obtain relevant astrophysical information about the compact sources themselves in a frequency range, ν??0 -- 200 GHz, where observations have only very recently started to be available. In this paper we propose a Bayesian maximum a posteriori (MAP) approach estimation scheme which incorporates prior information about the distribution of the polarization fraction of extragalactic compact sources between 1 and 100 GHz. We apply this Bayesian scheme to white noise simulations and to more realistic simulations that include CMB intensity, Galactic foregrounds and instrumental noise with the characteristics of the QUIJOTE experiment Wide Survey at 11 GHz. Using these simulations, we also compare our Bayesian method with the frequentist Filtered Fusion method that has been already used in WMAP data and in the \emph{Planck} mission. We find that the Bayesian method allows us to decrease the threshold for a feasible estimation of P to levels below ??00 mJy (as compared to ??00 mJy that was the equivalent threshold for the frequentist Filtered Fusion). We compare the bias introduced by the Bayesian method and find it to be small in absolute terms. Finally, we test the robustness of the Bayesian estimator against uncertainties in the prior and in the flux density of the sources. We find that the Bayesian estimator is robust against moderate changes in the parameters of the prior and almost insensitive to realistic errors in the estimated photometry of the sources.

Obtaining accurately calibrated redshift distributions of photometric samples is one of the great challenges in photometric surveys like LSST, Euclid, HSC, KiDS, and DES. We combine the redshift information from the galaxy photometry with constraints from two-point functions, utilizing cross-correlations with spatially overlapping spectroscopic samples. Our likelihood framework is designed to integrate directly into a typical large-scale structure and weak lensing analysis based on two-point functions. We discuss efficient and accurate inference techniques that allow us to scale the method to the large samples of galaxies to be expected in LSST. We consider statistical challenges like the parametrization of redshift systematics, discuss and evaluate techniques to regularize the sample redshift distributions, and investigate techniques that can help to detect and calibrate sources of systematic error using posterior predictive checks. We evaluate and forecast photometric redshift performance using data from the CosmoDC2 simulations, within which we mimic a DESI-like spectroscopic calibration sample for cross-correlations. Using a combination of spatial cross-correlations and photometry, we show that we can provide calibration of the mean of the sample redshift distribution to an accuracy of at least 0.002(1+z) , consistent with the LSST-Y1 science requirements for weak lensing and large-scale structure probes.

Emission from the interstellar medium can be a significant contaminant of measurements of the intensity and polarization of the cosmic microwave background (CMB). For planning CMB observations, and for optimizing foreground-cleaning algorithms, a description of the statistical properties of such emission can be helpful. Here we examine a machine learning approach to inferring the statistical properties of dust from either observational data or physics-based simulations. In particular, we apply a type of neural network called a Variational Auto Encoder (VAE) to maps of the intensity of emission from interstellar dust as inferred from Planck sky maps and demonstrate its ability to a) simulate new samples with similar summary statistics as the training set, b) provide fits to emission maps withheld from the training set, and c) produce constrained realizations. We find VAEs are easier to train than another popular architecture: that of Generative Adversarial Networks (GANs), and are better-suited for use in Bayesian inference.

We report on a search for dark matter axion conversion photons from the magnetosphere of the Galactic Center magnetar PSR J1745-2900 using spectra obtained from the Karl G. Jansky Very Large Array. No significant spectral features are detected. Using a hybrid model for PSR J1745-2900 and canonical assumptions about the dark matter density profile, we exclude axion models with axion-photon coupling g aγγ > 6-34 x 10 −12 GeV −1 with 95% confidence over the mass ranges 4.2-8.4, 18.6-26.9, 33.0-41.4, 53.7-62.1, and 126.0-159.3 μ eV. If there is a dark matter cusp, the limits reduce to g aγγ > 6-34 x 10 −14 GeV −1 , which overlap some axion models for the observed mass ranges > 33 μ eV. These limits may be improved by modeling the stimulated emission that can boost the axion-photon conversion process.

The identification of the electromagnetic counterpart candidate ZTF19abanrhr to the binary black hole merger GW190521 opens the possibility to infer cosmological parameters from this standard siren with a uniquely identified host galaxy. The distant merger allows for cosmological inference beyond the Hubble constant. Here we show that the three-dimensional spatial location of ZTF19abanrhr calculated from the electromagnetic data remains consistent with the updated sky localization of GW190521 provided by the LIGO-Virgo Collaboration. If ZTF19abanrhr is associated with the GW190521 merger and assuming a flat wCDM model we find that H 0 = 48 +24 −10 km/s/Mpc, Ω m = 0.39 +0.38 −0.29 , and w 0 =− 1.29 +0.63 −0.50 (median and 68% credible interval). If we use the Hubble constant value inferred from another gravitational-wave event, GW170817, as a prior for our analysis, together with assumption of a flat Λ CDM and the model-independent constraint on the physical matter density ω m from Planck, we find H 0 = 69.1 8.7 −6.0 km/s/Mpc.

We present a new theory for the hierarchical clustering of dark matter (DM) halos based on stochastic differential equations, that constitutes a change of perspective with respect to existing frameworks (e.g., the excursion set approach); this work is specifically focused on the halo mass function. First, we present a stochastic differential equation that describes fluctuations in the mass growth of DM halos, as driven by a multiplicative white (Gaussian) noise dependent on the spherical collapse threshold and on the power spectrum of DM perturbations. We demonstrate that such a noise yields an average drift of the halo population toward larger masses, that quantitatively renders the standard hierarchical clustering. Then, we solve the Fokker-Planck equation associated to the stochastic dynamics, and obtain the Press & Schechter mass function as a (stationary) solution. Moreover, generalizing our treatment to a mass-dependent collapse threshold, we obtain an exact analytic solution capable of fitting remarkably well the N-body mass function over a wide range in mass and redshift. All in all, the new perspective offered by the theory presented here can contribute to better understand the gravitational dynamics leading to the formation, evolution and statistics of DM halos across cosmic times.