Cosmology

Abstract:

The upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will enable astronomers to discover rare and distant astrophysical transients. Host-galaxy association is crucial for selecting the most scientifically interesting transients for follow-up. LSST deep drilling field (DDF) observations will detect distant transients occurring in galaxies below the detection limits of most all-sky catalogs. Here, we investigate the use of preexisting, field-specific catalogs for host identification in the DDFs and a ranking of their usefulness. We have compiled a database of 70 deep catalogs that overlap with the Rubin DDFs and constructed thin catalogs to be homogenized and combined for transient-host matching. A systematic ranking of their utility is discussed and applied based on the inclusion of information such as spectroscopic redshifts and morphological information. Utilizing this data against a Dark Energy Survey sample of supernovae with pre-identified hosts in the XMM-Large Scale Structure and the Extended Chandra Deep Field-South fields, we evaluate different methods for transient-host association in terms of both accuracy and processing speed. We also apply light data-cleaning techniques to identify and remove contaminants within our associations, such as diffraction spikes and blended galaxies where the correct host cannot be determined with confidence. We use a lightweight machine learning approach in the form of extreme gradient boosting to generate confidence scores in our contaminant selections and associated metrics. Finally, we discuss the computational expense of implementation within the LSST transient alert brokers, which will require efficient, fast-paced processing to handle the large stream of survey data.

Abstract:

Abstract A key question in extragalactic astronomy is how the star-formation rate density (SFRD) evolves over cosmic time. A powerful way of addressing this question is using radio-continuum observations, where the radio waves are unaffected by dust and are able to reach sufficient resolution to resolve individual galaxies. We present an investigation of the 1.4 GHz radio luminosity functions (RLFs) of star-forming galaxies (SFGs) and Active Galactic Nuclei (AGN) using deep radio continuum observations in the COSMOS and XMM–LSS fields, covering a combined area of ∼4 deg2. These data enable the most accurate measurement of the evolution in the SFRD from mid-frequency radio continuum observations. We model the total RLF as the sum of evolving SFG and AGN components, negating the need for individual source classification. We find that the SFGs have systematically higher space densities at fixed luminosity than found in previous radio studies, but consistent with more recent studies with MeerKAT. We attribute this to the excellent low-surface brightness sensitivity of MeerKAT. We then determine the evolution of the SFRD. Adopting the far-infrared – radio correlation results in a significantly higher SFRD at z > 1, compared to combined UV and far-infrared measurements. However, using more recent relations for the correlation between star-formation rate and radio luminosity, based on full spectral energy distribution modelling, can resolve this apparent discrepancy. Thus radio observations provide a powerful method of determining the total SFRD, in the absence of dust-sensitive far-infrared data.

Abstract:

Measuring galaxy rotation curves is critical for inferring the properties of dark-matter haloes in the Lambda Cold Dark Matter (ΛCDM) paradigm. We present H i rotation curves and mass models for 20 galaxies from the MIGHTEE survey. Using extended H i kinematics, we construct resolved mass models that include stellar, gaseous, and dark-matter components. Stellar masses are derived using 3.6 μm imaging under fixed mass-to-light ratio (ϒ* = M/L) assumptions and are complemented, for the first time for a H I-selected sample, by spatially resolved M/L, obtained from multi-wavelength SED fitting. We examine the ratio of baryonic to observed rotation velocity (Vbar/Vobs) at the characteristic radius R2.2. Adopting a fixed ϒ⋆ = 0.5 M⊙/L⊙ yields a clear dependence of V2.2/Vobs on galaxy luminosity, while adopting ϒ⋆ = 0.2 M⊙/L⊙ substantially weakens this trend. In contrast, the resolved M/L analysis preserves the luminosity dependence while modifying the stellar contribution on a galaxy-by-galaxy basis, providing a more accurate representation of the underlying relation. We model the dark-matter haloes using Navarro–Frenk–White profiles and find that the different assumptions for a fixed a M/L systematically shift galaxies relative to the theoretical stellar-to-halo mass and baryonic-to-halo mass relations, while the spatially varying M/L yields the closest agreement with theoretical benchmarks within ΛCDM. We therefore demonstrate that future investigations of the dark matter properties of galaxies using rotation curves need to account for varying M/L across individual galaxy profiles and between galaxies in order to obtain accurate measurements of the dark matter, and therefore test ΛCDM.

Abstract:

If dark matter annihilates into standard model particles with a cross section which is velocity dependent, then Local Group dwarf galaxies will not be the best place to search for the resulting gamma ray emission. A greater flux would be produced by more distant and massive halos, with larger velocity dispersions. We construct full-sky predictions for the gamma ray emission from galaxy- and cluster-mass halos within $\sim 200\mathrm{Mpc}$ using a suite of constrained $N$ -body simulations () based on the Bayesian Origin Reconstruction from Galaxies algorithm. Comparing to observations from the Large Area Telescope and marginalizing over reconstruction uncertainties and other astrophysical contributions to the flux, we obtain constraints on the cross section which are 2 (7) orders of magnitude tighter than those obtained from dwarf spheroidals for $p$ -wave ( $d$ -wave) annihilation. We find no evidence for either type of annihilation from dark matter particles with masses in the range $m_{\chi }=2–500\mathrm{GeV}/c^2$ , for any channel. As an example, for annihilations producing bottom quarks with $m_{\chi }=10\mathrm{GeV}/c^2$ , we find $a_1<2.4\times {10}^{-21}{\mathrm{cm}}^3{\mathrm{s}}^{-1}$ and $a_2<3.0\times {10}^{-18}{\mathrm{cm}}^3{\mathrm{s}}^{-1}$ at 95% confidence, where the product of the cross section, $\sigma$ , and relative particle velocity, $v$ , is given by $\sigma v=a_{\ell }(v/c{)}^{2\ell }$ and $\ell =1$ , 2 for $p$ - and $d$ -wave annihilation, respectively. Our bounds, although failing to exclude the thermal relic cross section for velocity-dependent annihilation channels, are among the tightest to date.

Abstract:

Skew spectra allow us to extract non-Gaussian information by taking the square of a map and finding the power spectrum of this new map with the original map. This allows us to use much of the infrastructure of power spectra and avoid the intricacies of estimating three point statistics. In this paper we present the first extension of skew spectra to arbitrary spin- $s$ fields, as a means to extract non-Gaussian information efficiently from cosmological datasets like cosmic shear or cosmic microwave background polarization. We apply the formalism to weak lensing in the context of large scale structure, and discuss different ways of combining fields to build skew spectra, all while avoiding the problems associated with mass mapping. We provide plots of these new statistics for $\mathrm{\Lambda }$ cold dark matter and vary cosmological parameters.