Beecroft Institute for Particle Astrophysics and Cosmology

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 (ϒ* = 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:

We present results of Population III (Pop III) formation in the megatron suite of simulations, which self-consistently follows radiation and non-equilibrium chemistry, and resolves gas at near-pc resolution in a Milky Way-mass progenitor at Cosmic Dawn. While the very first Pop III stars form in haloes with masses well below the atomic cooling limit, the majority of Pop III stars form in more massive systems above the K atomic cooling threshold as a Lyman–Werner (LW) background of is rapidly established. We find that the global Pop III star formation rate stabilizes to a value of at . Among the three processes that quench Pop III star formation in minihaloes, the LW background, gas starvation, and external chemical enrichment, the LW background is most important. A small fraction of haloes undergo multiple episodes of Pop III star formation when the earlier forming stars all directly collapse to black holes. If the haloes become massive enough, they can form up to Pop III stars in a single burst, which may be observable by James Webb Space Telescope with moderate gravitational lensing. Pop III stars form at a wide range of distances from UV-bright galaxies, with only per cent of Pop III stars forming within the virial radius of galaxies with . Finally, by tracking Pop III star remnants down to , we find that per cent reside in the stellar halo of our simulated Milky Way analogue, while the remainder are gravitationally bound to lower mass systems, including satellite haloes.

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.