Dr Harry Desmond

Abstract:

Empirical models of galaxy formation require assumptions about the correlations between galaxy and halo properties. These may be calibrated against observations or inferred from physical models such as hydrodynamical simulations. In this Letter, we use the EAGLE simulation to investigate the correlation of galaxy size with halo properties. We motivate this analysis by noting that the common assumption of angular momentum partition between baryons and dark matter in rotationally supported galaxies overpredicts both the spread in the stellar mass–size relation and the anticorrelation of size and velocity residuals, indicating a problem with the galaxy–halo connection it implies. We find the EAGLE galaxy population to perform significantly better on both statistics, and trace this success to the weakness of the correlations of galaxy size with halo mass, concentration and spin at fixed stellar mass. Using these correlations in empirical models will enable fine-grained aspects of galaxy scalings to be matched.

Abstract:

The Fundamental Plane (FP) describes the relation between the stellar mass, size, and velocity dispersion of elliptical galaxies; the Faber–Jackson relation (FJR) is its projection on to {mass, velocity} space. In this work, we re-deploy and expand the framework of Desmond & Wechsler to ask whether abundance matching-based Λ-cold dark matter models which have shown success in matching the spatial distribution of galaxies are also capable of explaining key properties of the FJR and FP, including their scatter. Within our framework, agreement with the normalization of the FJR requires haloes to expand in response to disc formation.We find that the tilt of the FP may be explained by a combination of the observed non-homology in galaxy structure and the variation in mass-to-light ratio produced by abundance matching with a universal initial mass function, provided that the anisotropy of stellar motions is taken into account. However, the predicted scatter around the FP is considerably increased by situating galaxies in cosmologically motivated haloes due to the variations in halo properties at fixed stellar mass and appears to exceed that of the data. This implies that additional correlations between galaxy and halo variables may be required to fully reconcile these models with elliptical galaxy scaling relations.

Abstract:

The Tully–Fisher relation (TFR) expresses the connection between rotating galaxies and the dark matter haloes they inhabit, and therefore contains a wealth of information about galaxy formation. We construct a general framework to investigate whether models based on halo abundance matching are able to reproduce the observed stellar mass TFR and mass–size relation (MSR), and use the data to constrain galaxy formation parameters. Our model tests a range of plausible scenarios, differing in the response of haloes to disc formation, the relative angular momentum of baryons and dark matter, the impact of selection effects, and the abundance matching parameters. We show that agreement with the observed TFR puts an upper limit on the scatter between galaxy and halo properties, requires weak or reversed halo contraction, and favours selection effects that preferentially eliminate fast-rotating galaxies. The MSR constrains the ratio of the disc to halo specific angular momentum to be approximately in the range 0.6–1.2. We identify and quantify two problems that models of this nature face. (1) They predict too large an intrinsic scatter for the MSR, and (2) they predict too strong an anticorrelation between the TFR and MSR residuals. We argue that resolving these problems requires introducing a correlation between stellar surface density and enclosed dark matter mass. Finally, we explore the expected difference between the TFRs of central and satellite galaxies, finding that in the favoured models this difference should be detectable in a sample of ∼700 galaxies.

Abstract:

Providing a theoretical basis for the baryonic Tully-Fisher Relation (BTFR; baryonic mass vs rotational velocity in spiral galaxies) in the LCDM paradigm has proved problematic. Simple calculations suggest too low a slope and too high a scatter, and recent semi-analytic models and numerical galaxy simulations typically fail to reproduce some aspects of the relation. Furthermore, the assumptions underlying one model are often inconsistent with those behind another. This paper aims to develop a rigorous prediction for the BTFR in the context of LCDM, using only a priori expected effects and relations, a minimum of theoretical assumptions, and no free parameters. The robustness of the relation to changes in key galactic parameters will be explored. I adopt a modular approach, taking each of the stand alone galaxy relations necessary for constructing the BTFR from up-to-date numerical simulations of dark halos. These relations -- and their expected scatter -- are used to describe model spirals with a range of masses, resulting in a band in the space of the BTFR that represents the current best guess for the LCDM prediction. Consistent treatment of expected LCDM effects goes a large way towards reconciling the naive slope-3 LCDM prediction with the data, especially in the range 10^9 M_sun < M_bar < 10^11 M_sun. The theoretical BTFR becomes significantly curved at M_bar > 10^11 M_sun, but this is difficult to test observationally due to the scarcity of extremely high mass spirals. Low mass gas-rich galaxies have systematically lower rotational velocity than the LCDM prediction, although the relation used to describe baryon mass fractions must be extrapolated in this regime. The fact that the BTFR slope derived here is significantly greater than in early predictions is a direct consequence of a corresponding increase in the expected sensitivity of baryon mass fraction to total halo mass.

Abstract:

Galactic scale tests have proven to be powerful tools in constraining fundamental physics in previously under-explored regions of parameter space. The astrophysical regime which they probe is inherently complicated, and the inference methods used to make these constraints should be robust to baryonic effects. Previous analyses have assumed simple empirical models for astrophysical noise without detailed calibration or justification. We outline a framework for assessing the reliability of such methods by constructing and testing more advanced baryonic models using cosmological hydrodynamical simulations. As a case study, we use the Horizon-AGN simulation to investigate warping of stellar disks and offsets between gas and stars within galaxies, which are powerful probes of screened fifth forces. We show that the degree of `U'-shaped warping of galaxies is well modelled by Gaussian random noise, but that the magnitude of the gas-star offset is correlated with the virial radius of the host halo. By incorporating this correlation we confirm recent results ruling out astrophysically relevant Hu-Sawicki $f(R)$ gravity, and identify a $\sim 30\%$ systematic uncertainty due to baryonic physics. Such an analysis must be performed case-by-case for future galactic tests of fundamental physics.