Any successful alternative gravity theory that obviates the need for dark matter must fit our cosmological observations. Measurements of microwave background polarization trace the large-scale baryon velocity field at recombination and show very strong, $O(1)$, baryon acoustic oscillations. Measurements of the large-scale structure of galaxies at low redshift show much weaker features in the spectrum. If the alternative gravity theory's dynamical equations for the growth rate of structure are linear, then the density field growth can be described by a Green's function: $\delta(\vec x,t) = \delta(\vec x,t')G(x,t,t')$. We show that the Green function, $G(x,t,t')$, must have dramatic features that erase the initial baryon oscillations. This implies an acceleration law that changes sign on the $\sim 150$ Mpc scale. On the other hand, if the alternative gravity theory has a large nonlinear term that couples modes on different scales, then the theory would predict large-scale non-Gaussian features in large-scale structure. These are not seen in the distribution of galaxies nor in the distribution of quasars. No proposed alternative gravity theory for dark matter seems to satisfy these 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]

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value.

Supermassive black hole (SMBH) binaries with masses of ~10^8--10^9 Msun are expected to dominate the contribution to the as-yet undetected gravitational wave background (GWB) signal at the nanohertz frequencies accessible to Pulsar Timing Arrays (PTA). We currently lack firm empirical constraints on the amplitude of the GWB due to the dearth of confirmed SMBH binaries in the required mass range. Using HST/WFC3 images, we have discovered a z~0.2 quasar hosted in a merger remnant with two closely separated (0.13'' or ~430pc) continuum cores at the heart of the galaxy SDSSJ1010+1413. The two cores are spatially coincident with two powerful [OIII]-emitting point sources with quasar-like luminosities (L_AGN ~ 5x10^46 erg/s, suggesting the presence of a bound SMBH system, each with M_BH > 4x10^8 Msun. We place an upper limit on the merging timescale of the SMBH pair of 2.5 billion years, roughly the Universe lookback time at z~0.2. There is likely a population of quasar binaries similar to SDSSJ1010+1413 that contribute to a stochastic GWB that should be detected in the next several years. If the GWB is not detected this could indicate that SMBHs merge only over extremely long timescales, remaining as close separation binaries for many Hubble times, the so-called `final-parsec problem'.

We present a sample of accreting supermassive black holes (SMBHs) in dwarf galaxies at $z<1$. We identify dwarf galaxies in the NEWFIRM Medium Band Survey with stellar masses $M_{\star}<3\times 10^{9} M_{\odot}$ that have spectroscopic redshifts from the DEEP2 survey and lie within the region covered by deep (flux limit of $\sim 5\times 10^{-17} - 6\times 10^{-16} \rm{erg \ cm}^{-2} \ \rm{s}^{-1}$) archival Chandra X-ray data. From our sample of $605$ dwarf galaxies, $10$ exhibit X-ray emission consistent with that arising from AGN activity. If black hole mass scales roughly with stellar mass, then we expect that these AGN are powered by SMBHs with masses of $\sim 10^5-10^6 \ M_{\odot}$ and typical Eddington ratios $\sim 5\%$. Furthermore, we find an AGN fraction consistent with extrapolations of other searches of $\sim 0.6-3\%$ for $10^9 \ M_{\odot} \leq M_{\star} \leq 3\times 10^{9} \ M_{\odot}$ and $0.1<z<0.6$. Our AGN fraction is in good agreement with a semi-analytic model, suggesting that as we search larger volumes we may use comparisons between observed AGN fractions and models to understand seeding mechanisms in the early universe.

The observation of GW170817 in both gravitational and electromagnetic waves provides a number of unique tests of general relativity. One question we can answer with this event is: Do large-wavelength gravitational waves and short-frequency photons experience the same number of spacetime dimensions? In models that include additional non-compact spacetime dimensions, as the gravitational waves propagate, they "leak" into the extra dimensions, leading to a reduction in the amplitude of the observed gravitational waves, and a commensurate systematic error in the inferred distance to the gravitational wave source. Electromagnetic waves would remain unaffected. We compare the inferred distance to GW170817 from the observation of gravitational waves, $d_L^\mathrm{GW}$, with the inferred distance to the electromagnetic counterpart NGC 4993, $d_L^\mathrm{EM}$. We constrain $d_L^\mathrm{GW} = (d_L^\mathrm{EM}/\mathrm{Mpc})^\mathrm{\gamma}$ with $\gamma = 1.01^{+0.04}_{-0.05}$ (for the SHoES value of $H_0$) or $\gamma = 0.99^{+0.03}_{-0.05}$ (for the Planck value of $H_0$), where all values are MAP and minimal 68% credible intervals. These constraints imply that gravitational waves propagate in $D=3+1$ spacetime dimensions, as expected in general relativity. In particular, we find that $D = 4.02^{+0.07}_{-0.10}$ (SHoES) and $D = 3.98^{+0.07}_{-0.09}$ (Planck). Furthermore, we place limits on the screening scale for theories with $D>4$ spacetime dimensions, finding that the screening scale must be greater than $\sim 20$ Mpc. We also place a lower limit on the lifetime of the graviton of $t > 4.50 \times 10^8$ yr.

Verlinde (2016) has proposed a new modified theory of gravity, Emergent Gravity (EG), as an alternative to dark matter. EG reproduces the Tully-Fisher relationship with no free parameters, and agrees with the velocity curves of most massive, spiral galaxies well. In its current form, the theory only applies to isolated, spherically symmetric systems in a dark energy-dominated Universe, and thus can only be tested fairly with such systems. Here I develop equations for EG's velocity curve predictions given a realistic, extended mass distribution. I apply this to isolated dwarf galaxies. Then I test the predictions from EG versus the maximum velocity measurements of 81 isolated dwarf galaxies with projected shapes close to circular. I find that EG severely underpredicts the maximum velocities for those galaxies with measured velocities v>165 km/s. Most of these galaxies have greater HI gas masses than stellar masses, and it is seems that EG is unable to describe these systems well. Rotation curves of these isolated, HI gas-rich, nearly spherical dwarf galaxies would provide the definitive test of EG.