Cosmology

Abstract:

We demonstrate the efficacy of symbolic regression (SR) to probe models of particle physics Beyond the Standard Model (BSM), by considering the so-called Constrained Minimal Supersymmetric Standard Model (CMSSM). Like many incarnations of BSM physics this model has a number (four) of arbitrary parameters, which determine the experimental signals, and cosmological observables such as the dark matter relic density. We show that analysis of the phenomenology can be greatly accelerated by using symbolic expressions derived for the observables in terms of the input parameters. Here we focus on the Higgs mass, the cold dark matter relic density and the contribution to the anomalous magnetic moment of the muon. We find that SR can produce remarkably accurate expressions. Using them we make global fits to derive the posterior probability densities of the CMSSM input parameters which are in good agreement with those performed using conventional methods. Moreover, we demonstrate a major advantage of SR, which is the ability to make fits using differentiable methods rather than sampling methods. We also compare the method with neural network (NN) regression. SR produces more globally robust results, while NNs require data that is focused on the promising regions in order to be equally performant. This article is part of the discussion meeting issue 'Symbolic regression in the physical sciences'.

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:

We present JWST observations of the radio galaxy TGSSJ1530+1049, spectroscopically confirmed at $z=4.0$ . NIRCam images and NIRSpec/IFU spectroscopy ( $R=2700$ ) show that TGSSJ1530+1049 is part of one of the densest-known structures of continuum and line-emitting objects found at these redshifts. NIRCam images show a number of distinct continuum objects and evidence of interactions traced by diffuse emission, and the NIRSpec IFU cube reveals further strong line emitting regions. We identify six continuum and four additional strong Halpha emitting sources with weaker or no underlying continuum within the 3’‘x3’’ IFU field. From spatial alignment with high-resolution radio data and emission line profiles, the radio AGN host galaxy is clearly identified. The bright Halpha emission (but not the optical components) is distributed remarkably linearly along the radio axis, suggestive of a biconical illumination zone by a central obscured AGN. The emission line kinematics indicate jet-gas interactions on scales of a few kpc. However, due to large relative velocities and presence of underlying continuum, the alignment with the radio structure appears to be, at least partly, caused by a particular configuration of interacting galaxies. At least four objects within a 10x10 (projected) kpc area which includes the radio source have high stellar masses (log( $M_{\star }/M_{\odot })>10.3$ ) and star formation rates in the range 70-163 $M_{\odot }$ yr. Using a stellar mass-based analysis, we predict a total dark matter halo mass of $\approx {10}^{13}{M}_{\odot }$ . Based on the physical separations and velocity differences between the galaxies, it is expected that these galaxies will merge to form a massive galaxy within a few Gyr. The system qualitatively resembles the forming brightest cluster galaxies in cosmological simulations that form early through a rapid succession of mergers.

Abstract:

We present a quadratic estimator that detects and reconstructs spatially varying multiplicative (m-) bias in weak lensing shear measurements, by exploiting the mode coupling that it generates. The method combines E and B modes with inverse-variance weights, to yield an unbiased reconstruction of to first order. We study the ability of future Stage IV surveys to obtain an unbiased reconstruction of the m-bias in differing scenarios, considering differing bias morphologies, and characteristic scales, as well as differing metrics to quantify the signal-to-noise ratio of the reconstructed map. We consider an m pattern repeating on sky patches, as might be the case for an m field caused by focal-plane systematics. With a Euclid-like redshift distribution, we find that root mean square (rms) variations in m-bias may be detected at the 20 level, after stacking between and patches (rising to between and for 1 per cent rms variations, data volumes that are becoming available with upcoming surveys), depending on the morphology of the m pattern. We show that these results are robust against the cosmological model assumed in the reconstruction, as well as the presence of intrinsic alignments or baryonic effects, and that the method shows no spurious response to additive (c-) bias. These results demonstrate that percent-level, spatially varying m-bias can be detected at high significance, enabling diagnosis and mitigation in the Stage IV weak lensing era.

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 () assumptions and are complemented, for the first time for a H i-selected sample, by spatially resolved , obtained from multiwavelength spectral energy distribution fitting. We examine the ratio of baryonic to observed rotation velocity () at the characteristic radius . Adopting a fixed yields a clear dependence of on galaxy luminosity, while adopting substantially weakens this trend. In contrast, the resolved 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 systematically shift galaxies relative to the theoretical stellar-to-halo mass and baryonic-to-halo mass relations, while the spatially varying 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 across individual galaxy profiles and between galaxies in order to obtain accurate measurements of the dark matter, and therefore test CDM.