We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction of spatially offset BHs increases with cosmic time, ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near $z=0$, and iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the black hole dynamics of future simulations.

The local distance ladder measurement of the Hubble constant requires a connection between geometric distances at low redshift and Type Ia supernovae in the Hubble flow, which may be achieved through either the Cepheid period--luminosity relation or the luminosity of the Tip of the Red Giant Branch (TRGB) feature of the Hertzsprung--Russell diagram. Any potential solution to the Hubble tension that works by altering the distance ladder must produce consistency of both the Cepheid and TRGB $H_0$ calibrations with the CMB. In this paper we extend our models of screened fifth forces (Desmond et al 2019) to cover the TRGB framework. A fifth force lowers TRGB luminosity, so a reduction in inferred $H_0$ requires that the stars that calibrate the luminosity---currently in the LMC---are on average less screened than those that calibrate the supernova magnitude. We show that even under pessimistic assumptions for the extinction to the LMC, full consistency with Planck can be achieved for a fifth force strength in unscreened RGB stars $\sim$0.2 that of Newtonian gravity. This is allowed by the comparison of Cepheid and TRGB distance measurements to nearby galaxies. Our results indicate that the framework of Desmond et al (2019) is more versatile than initially demonstrated, capable of ameliorating the Hubble tension on a second front.

The discrepancy between the values of the Hubble constant $H_0$ derived from the local distance ladder and the cosmic microwave background provides a tantalising hint of new physics. We explore a potential resolution involving screened fifth forces in the local Universe, which alter the Cepheid calibration of supernova distances. In particular, if the Cepheids with direct distance measurements from parallax or water masers are screened but a significant fraction of those in other galaxies are not, neglecting the difference between their underlying period--luminosity relations biases the local $H_0$ measurement high. This difference derives from a reduction in the Cepheid pulsation period and possible increase in luminosity under a fifth force. We quantify the internal and environmental gravitational properties of the Riess et al. distance ladder galaxies to assess their degrees of screening under a range of phenomenological models, and propagate this information into the $H_0$ posterior as a function of fifth force strength. We consider well-studied screening models in scalar--tensor gravity theories such as chameleon, K-mouflage and Vainshtein, along with a recently-proposed mechanism based on baryon--dark matter interactions in which screening is governed by local dark matter density. We find that a fifth force strength $\sim5-30\%$ that of gravity can alleviate (though not resolve) the $H_0$ tension in some scenarios, around the sensitivity level at which tests based on other distance ladder data can constrain this strength. Although our analysis is exploratory and based on screening models not necessarily realised in full theories, our results demonstrate that new physics-based local resolutions of the $H_0$ tension are possible, supplementing those already known in the pre-recombination era.

We derive the details of a new screening mechanism where the interactions of baryons and dark matter can be screened according to the local dark matter density. In this mechanism, the value of Newton's constant is dark matter density-dependent, allowing for the possibility that astrophysical phenomena are very different in galaxies less dense than the Milky Way. The parameterized post Newtonian parameter $\gamma$, which quantifies the difference between kinematical and lensing probes, also depends on dark matter density. We calculate the effects of varying $G$ on various stages of stellar evolution, focusing on observables that impact cosmology: the Cepheid period--luminosity relation and the supernova Ia magnitude--redshift relation. Other potential tests of the model are also investigated including main-sequence, post-main sequence, and low mass dwarf stars. Finally, we discuss how extragalactic tests of $\gamma$ could provide complementary constraints.

Self-interacting dark matter (SIDM) is an able alternative to collisionless dark matter. If dark matter does have self-interactions, we would expect this to cause a separation between the collisionless stars and the dark matter halo of a galaxy as it falls through a dark matter medium. For stars arranged in a disk, this would generate a U-shaped warp. The magnitude of this warping depends on the SIDM cross section, type of self-interaction, relative velocity of galaxy and background, halo structure, and density of the dark matter medium. In this paper, we set constraints on long-range (light mediator) dark matter self-interaction by means of this signal. We begin by measuring U-shaped warps in $3,213$ edge-on disk galaxies within the Sloan Digital Sky Survey. We then forward-model the expected warp from SIDM on a galaxy-by-galaxy basis by combining models of halo structure, density and velocity field reconstructions, and models for the dark matter interactions. We find no evidence for a contribution to the warps from SIDM. Our constraints are highly dependent on the uncertain velocities of our galaxies: for a normalized Rutherford-like cross section we find $\tilde{\sigma}/m_{\rm{DM}} \lesssim 3\times 10^{-13}~\rm{cm}^2/\rm{g}$ at fixed velocity $v = 300~\rm{km/s}$ -- a bound that scales roughly linearly with increasing $v$. In the appendix we translate these bounds into limits on the momentum transfer cross section, finding $\sigma_T(300~\rm{km/s})/m_{\rm{DM}} \lesssim 0.1~\rm{cm}^2/\rm{g}$. [abridged]

Cross-correlating two different types of galaxy gives rise to parity breaking in the correlation function that derives from differences in the galaxies' properties and environments. This is typically associated with a difference in galaxy bias, describing the relation between galaxy number density and dark matter density, although observational effects such as magnification bias also play a role. In this paper we show that the presence of a screened fifth force adds additional degrees of freedom to the correlation function, describing the effective coupling of the force to the two galaxy populations. These are also properties of the galaxies' environments, but with different dependence in general to galaxy bias. We derive the parity-breaking correlation function analytically as a function of fifth-force strength and the two populations' fifth-force charges, and explore the result numerically using Hu-Sawicki $f(R)$ as a toy model of chameleon screening. We find that screening gives rise to an octopole, which, in the absence of magnification bias, is not present in any gravity theory without screening and thus is a qualitatively distinct signature. The modification to the dipole and octopole can be $\mathcal{O}(10\%)$ and $\mathcal{O}(100\%)$ respectively at redshift $z \gtrsim 0.5$ due to screening, but decreases towards lower redshift. The change in the background power spectrum in $f(R)$ theories induces a change in the dipole of roughly the same size, but dominant to the effect of screening at low $z$. While current data is insufficient to measure the parity-breaking dipole or octopole to the precision required to test these models, future surveys such as DESI, Euclid and SKA have the potential to probe screened fifth forces through the dipole.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

In chameleon $f(R)$-gravity, the fifth force will lead to `upturns' in galaxy rotation curves near the screening radius. The location of the upturn depends on the cosmic background value of the scalar field $f_{R0}$, as well as the mass, size and environment of the galaxy. We search for this signature of modified gravity in the SPARC sample of measured rotation curves, using an MCMC technique to derive constraints on $f_{R0}$. Assuming NFW dark matter haloes and with $f_{R0}$ freely varying for each galaxy, most galaxies prefer $f(R)$ gravity to $\Lambda$CDM, but there is a large spread of inferred $f_{R0}$ values, inconsistent with a single global value. Requiring instead a consistent $f_{R0}$ value for the whole sample, models with $\log_{10}|f_{R0}| > -6.1$ are excluded. On the other hand, models in the range $-7.5<\log_{10}|f_{R0}|<-6.5$ seem to be favoured with respect to $\Lambda$CDM, with a significant peak at -7. However, this signal is largely a result of galaxies for which the $f(R)$ signal is degenerate with the core/cusp problem, and when the NFW profile is replaced with a cored halo profile, $\Lambda$CDM gives better fits than any given $f(R)$ model. Thus, we find no convincing evidence of $f(R)$ gravity down to the level of $|f_{R0}| \sim 6 \times 10^{-8}$, with the caveat that if cored halo density profiles cannot ultimately be explained within $\Lambda$CDM, a screened modified gravity theory could possibly provide an alternative solution for the core/cusp problem. However, the $f(R)$ models studied here fall short of achieving this.

We study the baryonic Tully-Fisher relation (BTFR) at z=0 using 153 galaxies from the SPARC sample. We consider different definitions of the characteristic velocity from HI and H-alpha rotation curves, as well as HI line-widths from single-dish observations. We reach the following results: (1) The tightest BTFR is given by the mean velocity along the flat part of the rotation curve. The orthogonal intrinsic scatter is extremely small (6%) and the best-fit slope is 3.85+/-0.09, but systematic uncertainties may drive the slope from 3.5 to 4.0. Other velocity definitions lead to BTFRs with systematically higher scatters and shallower slopes. (2) We provide statistical relations to infer the flat rotation velocity from HI line-widths or less extended rotation curves (like H-alpha and CO data). These can be useful to study the BTFR from large HI surveys or the BTFR at high redshifts. (3) The BTFR is more fundamental than the relation between angular momentum and galaxy mass (the Fall relation). The Fall relation has about 7 times more scatter than the BTFR, which is merely driven by the scatter in the mass-size relation of galaxies. The BTFR is already the "fundamental plane" of galaxy discs: no value is added with a radial variable as a third parameter.

The mass--velocity--size relation of late-type galaxies decouples into independent correlations between mass and velocity (the Tully-Fisher relation), and between mass and size. This behaviour is different to early-type galaxies which lie on a Fundamental Plane. We study the coupling of the Tully-Fisher and mass-size relations in observations (the SPARC sample) and in empirical galaxy formation models based on halo abundance matching, and rotation curve fits with a hydrodynamically motivated halo profile. We systematically investigate the correlation coefficient between the Tully-Fisher residuals $\Delta V_r$ and mass-size residuals $\Delta R$ as a function of the radius $r$ at which the velocity is measured, and thus present the $\Delta V_r-\Delta R$ relation across rotation curves. We find no significant correlation in either the data or models for any $r$, aside from $r \ll R_\text{eff}$ where baryonic mass dominates. We show that this implies an anticorrelation between galaxy size and halo concentration (or halo mass) at fixed baryonic mass, and provides evidence against the hypothesis that galaxy and halo specific angular momentum are proportional. Finally, we study the $\Delta V_r-\Delta R$ relations produced by the baryons and dark matter separately by fitting halo profiles to the rotation curves. The balance between these components illustrates the "disk-halo conspiracy" required for no overall correlation.

Extensions of the standard models of particle physics and cosmology often lead to long-range fifth forces with properties dependent on gravitational environment. Fifth forces on astrophysical scales are best studied in the cosmic web where perturbation theory breaks down. We present constraints on chameleon- and symmetron-screened fifth forces with Yukawa coupling and megaparsec range -- as well as unscreened fifth forces with differential coupling to galactic mass components -- by searching for the displacements they predict between galaxies' stars and gas. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we set upper bounds on fifth-force strength relative to Newtonian gravity from $\Delta G/G_N < \text{few} \: \times 10^{-4}$ for range $\lambda_C = 50$ Mpc, to $\Delta G/G_N \lesssim 0.1$ for $\lambda_C = 500$ kpc. In $f(R)$ gravity this requires $f_{R0} < \text{few} \: \times \: 10^{-8}$. The analogous bounds without screening are $\Delta G/G_N < \text{few} \: \times 10^{-4}$ and $\Delta G/G_N < \text{few} \times 10^{-3}$. These are the tightest and among the only fifth-force constraints on galaxy scales. We show how our results may be strengthened with future survey data and identify the key features of an observational programme for furthering fifth-force tests beyond the Solar System.

We present a new empirical relation between galaxy dark matter halo mass (${\rm M_{halo}}$) and the velocity along the flat portion of the rotation curve (${\rm V_{flat}}$), derived from 120 late-type galaxies from the SPARC database. The orthogonal scatter in this relation is comparable to the observed scatter in the baryonic Tully-Fisher relation (BTFR), indicating a tight coupling between total halo mass and galaxy kinematics at $r\ll R_{\rm vir}$. The small vertical scatter in the relation makes it an extremely competitive estimator of total halo mass. We demonstrate that this conclusion holds true for different priors on $M_*/L_{[3.6\mu]}$ that give a tight BTFR, but requires that the halo density profile follows DC14 rather than NFW. We provide additional relations between ${\rm M_{halo}}$ and other velocity definitions at smaller galactic radii (i.e. ${\rm V_{2.2}}$, ${\rm V_{eff}}$, and ${\rm V_{max}}$) which can be useful for estimating halo masses from kinematic surveys, providing an alternative to abundance matching. Furthermore, we constrain the dark matter analog of the Radial Acceleration Relation and also find its scatter to be small, demonstrating the fine balance between baryons and dark matter in their contribution to galaxy kinematics.

Intra-galaxy signals contain a wealth of information on fundamental physics, both the dark sector and the nature of gravity. While so far largely unexplored, such probes are set to rise dramatically in importance as upcoming surveys provide data of unprecedented quantity and quality on galaxy structure and dynamics. In this paper, we use warping of stellar disks to test the chameleon- or symmetron-screened fifth forces which generically arise when new fields couple to matter. We take r-band images of mostly late-type galaxies from the Nasa Sloan Atlas and develop an automated algorithm to quantify the degree of U-shaped warping they exhibit. We then forward-model the warp signal as a function of fifth-force strength $\Delta G/G_N$ and range $\lambda_C$, and the gravitational environments and internal properties of the galaxies, including full propagation of the non-Gaussian uncertainties. Convolving this fifth-force likelihood function with a Gaussian describing astrophysical and observational noise and then constraining $\Delta G/G_N$ and $\lambda_C$ by Markov Chain Monte Carlo, we find the overall likelihood to be significant increased ($\Delta\log(\mathcal{L}) \simeq 20$) by adding a screened fifth force with $\lambda_C \simeq 2$ Mpc, $\Delta G/G_N \simeq 0.01$. The variation of $\Delta\log(\mathcal{L})$ with $\lambda_C$ is quantitatively as expected from the correlation of the magnitude of the fifth-force field with the force's range, and a similar model without screening achieves no increase in likelihood over the General Relativistic case $\Delta G=0$. Although these results are in good agreement with a previous analysis of the same model using offsets between galaxies' stellar and gas mass centroids (Desmond et al. 2018), we caution that the effects of confounding baryonic and dark matter physics must be thoroughly investigated for the results of the inference to be unambiguous.

One of the most common consequences of extensions to the standard models of particle physics or cosmology is the emergence of a fifth force. While generic fifth forces are tightly constrained at Solar System scales and below, they may escape detection by means of a screening mechanism which effectively removes them in dense environments. We constrain the strength $\Delta G/G_N$ and range $\lambda_C$ of a chameleon- or symmetron-screened fifth force with Yukawa coupling -- as well as an unscreened fifth force with differential coupling to galactic mass components -- by searching for the displacement it predicts between galaxies' stellar and gas mass centroids. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we find $6.6\sigma$ evidence for $\Delta G>0$ at $\lambda_C \simeq 2$ Mpc, with $\Delta G/G_N = 0.025$ at maximum-likelihood. A similar fifth-force model without screening gives no increase in likelihood over the case $\Delta G = 0$ for any $\lambda_C$. Although we validate these results by several methods, we do not claim screened modified gravity to provide the only possible explanation for the signal: this would require knowing that "galaxy formation" physics could not be responsible. We show also the results of a more conservative -- though less well motivated -- noise model which yields only upper limits on $\Delta G/G_N$, ranging from $\sim10^{-1}$ for $\lambda_C \simeq 0.5$ Mpc to $\sim \: \text{few} \times 10^{-4}$ at $\lambda_C \simeq 50$ Mpc. We show how these constraints may be improved by future galaxy surveys and identify the key features of an observational programme for directly constraining fifth forces on galactic scales. This paper provides a complete description of the analysis summarised in Desmond et al. (2018b).

In a $\Lambda$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems. Assuming that haloes initially accrete the cosmological baryon fraction, we show that the maximum energy available from supernovae is typically not enough to completely eject these "missing baryons" from a halo, but it is often sufficient to heat them to the virial temperature. The energy available from supernovae has the same scaling with galaxy mass as the energy needed to heat or eject the "missing baryons", indicating that the coupling efficiency of the feedback to the ISM may be constant with galaxy virial mass. We further find that the energy available from supernova feedback is always enough to convert a primordial cusp into a core and has magnitude consistent with what is required to heat the "missing baryons" to the virial temperature. Taking a census of the baryon content of galaxies with ${\rm 10^9<M_{vir}/M_{\odot}<10^{12}}$ reveals that $\sim86\%$ of baryons are likely to be in a hot phase surrounding the galaxies and possibly observable in the X-ray, $\sim7\%$ are in the form of cold gas, and $\sim7\%$ are in stars.

Tests of gravity at the galaxy scale are in their infancy. As a first step to systematically uncovering the gravitational significance of galaxies, we map three fundamental gravitational variables -- the Newtonian potential, acceleration and curvature -- over the galaxy environments of the local universe to a distance of approximately 200 Mpc. Our method combines the contributions from galaxies in an all-sky redshift survey, halos from an N-body simulation hosting low-luminosity objects, and linear and quasi-linear modes of the density field. We use the ranges of these variables to determine the extent to which galaxies expand the scope of generic tests of gravity and are capable of constraining specific classes of model for which they have special significance. Finally, we investigate the improvements afforded by upcoming galaxy surveys.

In recent work, Lelli et al. (2016) argue that the tightness of the baryonic Tully-Fisher relation (BTFR) of the SPARC galaxy sample, and the weakness of the correlation of its residuals with effective radius, pose challenges to LCDM cosmology. In this Letter we calculate the statistical significance of these results in the framework of halo abundance matching, which imposes a canonical galaxy-halo connection. Taking full account of sample variance among SPARC-like realisations of the parent halo population, we find the scatter in the predicted BTFR to be 3.6 sigma too high, but the correlation of its residuals with galaxy size to be naturally weak. Further, we find abundance matching to generate BTFR curvature in 3.0 sigma disagreement with the data, and a fraction of galaxies with non-flat rotation curves somewhat larger than observed.

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.

We use the mass discrepancy-acceleration relation (the correlation between the ratio of total-to-visible mass and acceleration in galaxies; MDAR) to test the galaxy-halo connection. We analyse the MDAR using a set of 16 statistics that quantify its four most important features: its shape, its scatter, the presence of a "characteristic acceleration scale," and the correlation of its residuals with other galaxy properties. We construct an empirical framework for the galaxy-halo connection in LCDM to generate predictions for these statistics, starting with conventional correlations (halo abundance matching; AM) and introducing more where required. Comparing to the SPARC data, we find that: 1) the approximate shape of the MDAR is readily reproduced by AM, and there is no evidence that the acceleration at which dark matter becomes negligible has less spread in the data than in AM mocks; 2) even under conservative assumptions, AM significantly overpredicts the scatter in the relation and its normalisation at low acceleration, and furthermore positions dark matter too close to galaxies' centres on average; 3) the MDAR affords 2 sigma evidence for an anticorrelation of galaxy size and Hubble type with halo mass or concentration at fixed stellar mass. Our analysis lays the groundwork for a bottom-up determination of the galaxy-halo connection from relations such as the MDAR, provides concrete statistical tests for specific galaxy formation models, and brings into sharper focus the relative evidence accorded by galaxy kinematics to LCDM and modified gravity alternatives.

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 onto {mass, velocity} space. In this work we redeploy and expand the framework of Desmond & Wechsler (2015) to ask whether abundance matching-based LCDM models that 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 normalisation 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 (IMF), 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 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.

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.

Recent work by Efstathiou (2014) highlighted the importance of outliers in the period-luminosity (PL) relation of Cepheid data on the distance ladder. We present a statistical framework designed to address this difficulty, and apply it to the Cepheid data from the Milky Way (MW), the Large Magellanic Cloud (LMC), and the Riess et al. (2011) (hereafter R11) dataset. We consider two possible models of the outlier population in the R11 Cepheid dataset. One of these models exhibits tension between the PL relation of the R11 cepheids and the MW+LMC cepheids, while the other does not. We extend our models to adequately account for tension between the cepheid data sets when appropriate. Our outlier treatment has a significant impact on the distance scales to Supernovae hosts with Cepheid distances, increasing the uncertainty in these distances by a median factor of ~30%. We further find that our Cepheid outlier treatment translates into a modest, but non-negligible increase in the statistical uncertainty of H0, adding in quadrature 1.2 km/s/Mpc. Combined with the increased scatter in the Hubble diagram reported by Jones et al. (2015), we find H0=72.6+/-2.8 km/s/Mpc, corresponding to a 3.8% uncertainty in H0 from local measurements. This value is fully consistent with both the Planck and inverse-distance ladder H0 constraints.

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.