skip to main content


Title: The evolution of galaxy intrinsic alignments in the MassiveBlackII universe
ABSTRACT We investigate the redshift evolution of the intrinsic alignments (IAs) of galaxies in the MassiveBlackII (MBII) simulation. We select galaxy samples above fixed subhalo mass cuts ($M_h\gt 10^{11,12,13}\,\mathrm{M}_{\odot }\, h^{-1}$) at z = 0.6 and trace their progenitors to z = 3 along their merger trees. Dark matter components of z = 0.6 galaxies are more spherical than their progenitors while stellar matter components tend to be less spherical than their progenitors. The distribution of the galaxy–subhalo misalignment angle peaks at ∼10 deg with a mild increase with time. The evolution of the ellipticity–direction (ED) correlation amplitude ω(r) of galaxies (which quantifies the tendency of galaxies to preferentially point towards surrounding matter overdensities) is governed by the evolution in the alignment of underlying dark matter (DM) subhaloes to the matter density of field, as well as the alignment between galaxies and their DM subhaloes. At scales $\sim 1~\mathrm{Mpc}\, h^{-1}$, the alignment between DM subhaloes and matter overdensity gets suppressed with time, whereas the alignment between galaxies and DM subhaloes is enhanced. These competing tendencies lead to a complex redshift evolution of ω(r) for galaxies at $\sim 1~\mathrm{Mpc}\, h^{-1}$. At scales $\gt 1~\mathrm{Mpc}\, h^{-1}$, alignment between DM subhaloes and matter overdensity does not evolve significantly; the evolution of the galaxy–subhalo misalignment therefore leads to an increase in ω(r) for galaxies by a factor of ∼4 from z = 3 to 0.6 at scales $\gt 1~\mathrm{Mpc}\, h^{-1}$. The balance between competing physical effects is scale dependent, leading to different conclusions at much smaller scales ($\sim 0.1~\mathrm{Mpc}\, h^{-1}$).  more » « less
Award ID(s):
1716131
NSF-PAR ID:
10286739
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
491
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
4116 to 4130
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We introduce a novel technique for constraining cosmological parameters and galaxy assembly bias using non-linear redshift-space clustering of galaxies. We scale cosmological N-body simulations and insert galaxies with the SubHalo Abundance Matching extended (SHAMe) empirical model to generate over 175 000 clustering measurements spanning all relevant cosmological and SHAMe parameter values. We then build an emulator capable of reproducing the projected galaxy correlation function at the monopole, quadrupole, and hexadecapole level for separations between $0.1\, h^{-1}\, {\rm Mpc}$ and $25\, h^{-1}\, {\rm Mpc}$. We test this approach by using the emulator and Monte Carlo Markov Chain (MCMC) inference to jointly estimate cosmology and assembly bias parameters both for the MTNG740 hydrodynamic simulation and for a semi-analytical model (SAM) galaxy formation built on the MTNG740-DM dark matter-only simulation, obtaining unbiased results for all cosmological parameters. For instance, for MTNG740 and a galaxy number density of $n\sim 0.01 h^{3}\, {\rm Mpc}^{-3}$, we obtain $\sigma _{8}=0.799^{+0.039}_{-0.044}$ and $\Omega _\mathrm{M}h^2= 0.138^{+ 0.025}_{- 0.018}$ (which are within 0.4 and 0.2σ of the MTNG cosmology). For fixed Hubble parameter (h), the constraint becomes $\Omega _\mathrm{M}h^2= 0.137^{+ 0.011}_{- 0.012}$. Our method performs similarly well for the SAM and for other tested sample densities. We almost always recover the true amount of galaxy assembly bias within 1σ. The best constraints are obtained when scales smaller than $2\, h^{-1}\, {\rm Mpc}$ are included, as well as when at least the projected correlation function and the monopole are incorporated. These methods offer a powerful way to constrain cosmological parameters using galaxy surveys.

     
    more » « less
  2. ABSTRACT

    Our understanding of reionization has advanced considerably over the past decade, with several results now demonstrating that the intergalactic medium transitioned from substantially neutral at z = 7 to largely reionized at z = 6. However, little remains known about the sizes of ionized bubbles at z ≳ 7 as well as the galaxy overdensities which drive their growth. Fortunately, rest-ultraviolet (UV) spectroscopic observations offer a pathway towards characterizing these ionized bubbles thanks to the resonant nature of Lyman-alpha photons. In a previous work, we presented Ly α detections from three closely separated Lyman-break galaxies at z ≃ 6.8, suggesting the presence of a large (R > 1 physical Mpc) ionized bubble in the 1.5 deg2 COSMOS field. Here, we present new deep Ly α spectra of 10 UV-bright ($\mathrm{\mathit{ M}}_{\mathrm{UV}}^{} \le -20.4$) z ≃ 6.6–6.9 galaxies in the surrounding area, enabling us to better characterize this potential ionized bubble. We confidently detect (S/N > 7) Ly α emission at z = 6.701–6.882 in nine of ten observed galaxies, revealing that the large-scale volume spanned by these sources (characteristic radius R = 3.2 physical Mpc) traces a strong galaxy overdensity (N/〈N〉 ≳ 3). Our data additionally confirm that the Ly α emission of UV-bright galaxies in this volume is significantly enhanced, with 40 per cent (4/10) showing strong Ly α emission (equivalent width >25 Å) compared to the 8–9 per cent found on average at z ∼ 7. The median Ly α equivalent width of our observed galaxies is also ≈2 times that typical at z ∼ 7, consistent with expectations if a very large (R ∼ 3 physical Mpc) ionized bubble is allowing the Ly α photons to cosmologically redshift far into the damping wing before encountering H i.

     
    more » « less
  3. ABSTRACT

    We perform high-resolution simulations of an MW-like galaxy in a self-interacting cold dark matter model with elastic cross-section over mass of $1~\rm cm^2\, g^{-1}$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disc and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances, and radial distributions within the virial radius. The SIDM host halo is denser in the centre than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size ${\sim}7 \, \rm kpc$). The most massive subhaloes (with $V_{\mathrm{peak}}\gt 20 \, \rm km\, s^{-1}$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogues at radii less than 500 pc but the deviation diminishes for less massive subhaloes. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhaloes do not display the too-big-to-fail problem. However, the least dense among the massive subhaloes in both these simulations tend to have the smallest pericenter values, a trend that is not apparent among the bright MW satellite galaxies.

     
    more » « less
  4. null (Ed.)
    ABSTRACT Galaxy intrinsic alignments (IAs) have long been recognized as a significant contaminant to weak lensing-based cosmological inference. In this paper we seek to quantify the impact of a common modelling assumption in analytic descriptions of IAs: that of spherically symmetric dark matter haloes. Understanding such effects is important as the current generation of IA models are known to be limited, particularly on small scales, and building an accurate theoretical description will be essential for fully exploiting the information in future lensing data. Our analysis is based on a catalogue of 113 560 galaxies between z = 0.06 and 1.00 from massiveblack-ii, a hydrodynamical simulation of box length $100 \, h^{-1}$ Mpc. We find satellite anisotropy contributes at the level of $\ge 30\!-\!40{{\ \rm per\ cent}}$ to the small-scale alignment correlation functions. At separations larger than $1 \, h^{-1}$ Mpc the impact is roughly scale independent, inducing a shift in the amplitude of the IA power spectra of $\sim 20{{\ \rm per\ cent}}$. These conclusions are consistent across the redshift range and between the massiveblack-ii and the illustris simulations. The cosmological implications of these results are tested using a simulated likelihood analysis. Synthetic cosmic shear data are constructed with the expected characteristics (depth, area, and number density) of a future LSST-like survey. Our results suggest that modelling alignments using a halo model based upon spherical symmetry could potentially induce cosmological parameter biases at the ∼1.5σ level for S8 and w. 
    more » « less
  5. ABSTRACT Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass haloes, accelerated if haloes are subhaloes of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of ${\lesssim}10^8 \, \mathrm{M}_\odot$ subhaloes in lens-mass systems (${\sim}10^{13}\, \mathrm{M}_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhaloes with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhaloes (≲1/1000 host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhaloes to core-collapse, finding that it is nearly impossible to drive core-collapse in subhaloes in SIDM models with constant cross-sections. Any discovery of ultracompact dark substructures with future substructure lensing observations favours additional degrees of freedom, such as velocity-dependence, in the cross-section. 
    more » « less