It has been claimed that the standard model of cosmology (LCDM) cannot easily account for a number of observations on relatively small scales, motivating extensions to the standard model. Here we introduce a new suite of cosmological simulations that systematically explores three plausible extensions: warm dark matter, self-interacting dark matter, and a running of the scalar spectral index of density fluctuations. Current observational constraints are used to specify the additional parameters that come with these extensions. We examine a large range of observable metrics on small scales, including the halo mass function, density and circular velocity profiles, the abundance of satellite subhaloes, and halo concentrations. For any given metric, significant degeneracies can be present between the extensions. In detail, however, the different extensions have quantitatively distinct mass and radial dependencies, suggesting that a multi-probe approach over a range of scales can be used to break the degeneracies. We also demonstrate that the relative effects on the radial density profiles in the different extensions (compared to the standard model) are converged down to significantly smaller radii than are the absolute profiles. We compare the derived cosmological trends with the impact of baryonic physics using the EAGLE and ARTEMIS simulations. Significant degeneracies are also present between baryonic physics and cosmological variations (with both having similar magnitude effects on some observables). Given the inherent uncertainties both in the modelling of galaxy formation physics and extensions to LCDM, a systematic and simultaneous exploration of both is strongly warranted.

We probe the self-interactions of dark matter using observational data of relaxed galaxy groups and clusters. Our analysis uses the Jeans formalism and considers a wider range of systematic effects than in previous work, including adiabatic contraction and stellar anisotropy, to robustly constrain the self-interaction cross section. For both groups and clusters, our results show a mild preference for a nonzero cross section compared with cold collisionless dark matter. Our groups result, $\sigma/m=0.5\pm0.2~\mathrm{cm}^2/\mathrm{g}$, places the first constraint on self-interacting dark matter (SIDM) at an intermediate scale between galaxies and massive clusters. Our clusters result is $\sigma/m=0.19\pm0.09~\mathrm{cm}^2/\mathrm{g}$, with an upper limit of $\sigma / m < 0.35~\mathrm{cm}^2/\mathrm{g}$ (95% CL). Thus, our results disfavor a velocity-independent cross section of order $1~\mathrm{cm}^2/\mathrm{g}$ or larger needed to address small scale structure problems in galaxies, but are consistent with a velocity-dependent cross section that decreases with increasing scattering velocity. Comparing the cross sections with and without the effect of adiabatic contraction, we find that adiabatic contraction produces slightly larger values for our data sample, but they are consistent at the $1\sigma$ level. Finally, to validate our approach, we apply our Jeans analysis to a sample of mock data generated from SIDM-plus-baryons simulations with $\sigma/m = 1~\mathrm{cm}^2/\mathrm{g}$. This is the first test of the Jeans model at the level of stellar and lensing observables directly measured from simulations. We find our analysis gives a robust determination of the cross section, as well as consistently inferring the true baryon and dark matter density profiles.

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

We quantify the performance of mass mapping techniques on mock imaging and gravitational lensing data of galaxy clusters. The optimum method depends upon the scientific goal. We assess measurements of clusters' radial density profiles, departures from sphericity, and their filamentary attachment to the cosmic web. We find that mass maps produced by direct (KS93) inversion of shear measurements are unbiased, and that their noise can be suppressed via filtering with MRLens. Forward-fitting techniques, such as Lenstool, suppress noise further, but at a cost of biased ellipticity in the cluster core and over-estimation of mass at large radii. Interestingly, current searches for filaments are noise-limited by the intrinsic shapes of weakly lensed galaxies, rather than by the projection of line-of-sight structures. Therefore, space-based or balloon-based imaging surveys that resolve a high density of lensed galaxies, could soon detect one or two filaments around most clusters.

In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by simulations of CDM that include baryons, and that a SIDM cross-section of 1 cm$^2/$g is marginally excluded. This constraint is driven mainly by galaxy clusters, with less massive systems being relatively less affected by DM self-interactions with velocity-independent cross-section.

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

Many distant objects can only be detected, or become more scientifically valuable, if they have been highly magnified by strong gravitational lensing. We use EAGLE and BAHAMAS, two recent cosmological hydrodynamical simulations, to predict the probability distribution for both the lens mass and lens redshift when point sources are highly magnified by gravitational lensing. For sources at a redshift of two, we find the distribution of lens redshifts to be broad, peaking at z=0.6. The contribution of different lens masses is also fairly broad, with most high-magnification lensing due to lenses with halo masses between 10^12 and 10^14 solar masses. Lower mass haloes are inefficient lenses, while more massive haloes are rare. We find that a simple model in which all haloes have singular isothermal sphere density profiles can approximately reproduce the simulation predictions, although such a model over-predicts the importance of haloes with mass <10^12 solar masses for lensing. We also calculate the probability that point sources at different redshifts are strongly lensed. At low redshift, high magnifications are extremely unlikely. Each z=0.5 source produces, on average, 5x10^-7 images with magnification greater than ten; for z =2 this increases to about 2x10^-5. Our results imply that searches for strongly lensed optical transients, including the optical counterparts to strongly lensed gravitational waves, can be optimized by monitoring massive galaxies, groups and clusters rather than concentrating on an individual population of lenses.

We present a cosmological hydrodynamical simulation of a 10^13 Msun galaxy group and its environment (out to 10 times the virial radius) carried out using the EAGLE model of galaxy formation. Exploiting a novel technique to increase the resolution of the dark matter calculation independently of that of the gas, the simulation resolves dark matter haloes and subhaloes of mass 5x10^6 Msun . It is therefore useful for studying the abundance and properties of the haloes and subhaloes targeted in strong lensing tests of the cold dark matter model. We estimate the halo and subhalo mass functions and discuss how they are affected both by the inclusion of baryons in the simulation and by the environment. We find that the halo and subhalo mass functions have lower amplitude in the hydrodynamical simulation than in its dark matter only counterpart. This reflects the reduced growth of haloes in the hydrodynamical simulation due to the early loss of gas by reionisation and galactic winds and, additionally, in the case of subhaloes, disruption by enhanced tidal effects within the host halo due to the presence of a massive central galaxy. The distribution of haloes is highly anisotropic reflecting the filamentary character of mass accretion onto the cluster. As a result, there is significant variation in the number of structures with viewing direction. The median number of structures near the centre of the halo, when viewed in projection, is reduced by a factor of two when baryons are included.

We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effects responsible for the previously claimed inconsistency. First, galaxy clusters containing baryonic matter have higher central densities than their counterparts simulated with only dark matter. Second, a sample of clusters selected because of the presence of pronounced gravitational lensing arcs preferentially finds centrally concentrated clusters with large Einstein radii. Third, lensed arcs are usually straighter than critical curves, and the chosen image analysis method (fitting circles through the arcs) overestimates the Einstein radii. After accounting for these three effects, LCDM predicts that galaxy clusters should produce giant lensing arcs that match those in the observed Universe.

Motivated by discovering strongly-lensed supernovae, gravitational waves, and kilonovae in the 2020s, we investigate whether to build a watch-list of clusters based on observed cluster properties (i.e. lens-plane selection) or on the detectability of strongly-lensed background galaxies (i.e. source-plane selection). First, we estimate the fraction of high-redshift transient progenitors that reside in galaxies that are themselves too faint to be detected as being strongly-lensed. We find $\sim15-50$ per cent of transient progenitors reside in $z = 1-2$ galaxies too faint to be detected in surveys that reach ${\rm AB}\simeq23$, such as the Dark Energy Survey. This falls to $\ls10$ per cent at depths that will be probed by early data releases of LSST (${\rm AB}\simeq25$). Second, we estimate a conservative lower limit on the fraction of strong lensing clusters that will be missed by magnitude limited searches for multiply-imaged galaxies and giant arcs due to the faintness of such images. We find that DES-like surveys will miss $\sim75$ per cent of $10^{15}$M$_\odot$ strong lensing clusters, rising to $\sim100$ per cent of $10^{14}$M$_\odot$ clusters. Deeper surveys, such as LSST, will miss $\sim40$ per cent at $10^{15}$M$_\odot$, and $\sim95$ per cent at $10^{14}$M$_\odot$. Our results motivate building a cluster watch-list for strongly-lensed transients that includes those found by lens-plane selection.

We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.

The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We show that the stellar mass varies systematically with resolution by over a factor of two, in a manner that depends on the final stellar mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. 70~per~cent of WDM haloes of mass $>10^{8}M_{\odot}$ collapse early enough to form stars before $z=6$, compared to 90~per~cent of CDM and SIDM galaxies. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization, in order to use these observables to distinguish between dark matter models.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

A non-parametric higher-order Jeans analysis method, GravSphere, was recently used to constrain the density profiles of Local Group dwarf galaxies, classifying them into those that are more likely to have an inner dark matter cusp and those that are likely to have a core (Read et al.). In this work we test this method using 31 simulated galaxies, comparable to Fornax, selected from the APOSTLE suite of cosmological hydrodynamics simulations, which include CDM and Self-Interacting Dark Matter (SIDM) cosmologies. We find that the mass profiles of the simulated dwarfs are often, but not always, well recovered by GravSphere. The less successful cases may be identified using a chi-squared diagnostic. Although the uncertainties are large in the inner regions, the inferred mass profiles are unbiased and exhibit smaller scatter than comparable Jeans methods. We find that GravSphere recovers the density profiles of simulated dwarfs below the half-light radius and down to the resolution limit of our simulations with better than 10+-30 per cent accuracy, making it a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.

Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed been observed - but only interpreted via a simplified, analytic model of oscillations. We compare these observations to the BAHAMAS-SIDM suite of cosmological simulations, which include SIDM and a fully hydrodynamical treatment of star formation and feedback. We predict that the median offset of BCGs increases with the SIDM cross-section, cluster mass and the amount of stellar mass within 10kpc, while CDM exhibits no trend in mass. Interpolating between the simulated cross-sections, we find that the observations (of 10 clusters) is consistent with CDM at the ~1.5$\sigma$ level, and prefer cross-section $\sigma$/m < 0.12(0.39)cm$^2$/g at 68% (95%) confidence level. This is on the verge of ruling out velocity-independent dark matter self-interactions as the solution to discrepancies between the predicted and observed behaviour of dwarf galaxies, and will be improved by larger surveys by Euclid or SuperBIT.

N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos.

We present observations of two strong-lensing galaxy clusters located within the $90$ per cent credible sky localization maps released following LIGO-Virgo's discovery of the binary black hole (BH-BH) gravitational wave (GW) source GW170814. Our objectives were (1) to search for candidate electromagnetic (EM) counterparts to GW170814 under the hypothesis that it was strongly-lensed, and thus more distant and less massive than inferred by LIGO-Virgo, and (2) to demonstrate the feasibility of rapid target of opportunity observations to search for faint lensed transient point sources in crowded cluster cores located within GW sky localizations. Commencing $20$ hours after discovery, and continuing over $12$ nights, we observed Abell 3084 ($z=0.22$) and SMACSJ0304.3$-$4401 ($z=0.46$) with GMOS on the Gemini-South telescope, and Abell 3084 with MUSE on ESO's Very Large Telescope. We detect no candidate EM counterparts in these data. Calibration of our photometric analysis methods using simulations yield $5\sigma$ detection limits for transients in difference images of the cores of these clusters of $i=25$. This is the most sensitive photometric search to date for counterparts to GW sources, and rules out the possibility that GW170814 was lensed by these clusters with a kilonova-like EM counterpart. Based on the detector frame masses of the compact objects, and assuming that at least one Neutron Star (NS) is required in the merging system to produce a kilonova-like counterpart, implies that GW170814 was neither a NS-NS nor NS-BH merger at $z>8$ lensed by either of these clusters. Also, in the first ever emission line search for counterparts to GW sources, we detected no lines down to a $5\sigma$ detection limit of $5\times10^{-17}{\rm erg\,s^{-1}\,cm^{-2}}$.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

LSST's wide-field of view and sensitivity will revolutionize studies of the transient sky by finding extraordinary numbers of new transients every night. The recent discovery of a kilonova counterpart to LIGO/Virgo's first detection of gravitational waves (GWs) from a double neutron star (NS-NS) merger also creates an exciting opportunity for LSST to offer a Target of Opportunity (ToO) mode of observing. We have been exploring the possibility of detecting strongly lensed GWs, that would enable new tests of GR, extend multi-messenger astronomy out to $z\gtrsim1$, and deliver a new class of sub-millisecond precision time-delay constraints on lens mass distributions. We forecast that the rate of detection of lensed NS-NS mergers in the 2020s will be $\sim0.1$ per Earth year, that the typical source will be at $z\simeq2$, and that the multiply-imaged kilonova counterpart will have a magnitude of ${\rm AB}\simeq25.4$ in $g/r/i$-band filters - i.e. fainter than the sensitivity of a single LSST WFD visit. We therefore advocate (1) creating a flexible and efficient Target of Opportunity programme within the LSST observing strategy that is capable of discovering sources fainter than single-visit depth, and (2) surveying the entire observable extragalactic sky as rapidly as possible in the WFD survey. The latter will enable a very broad range of early science that relies on wide survey area for detection of large samples of objects and/or maximizing the fraction of sky over which reference imaging is available. For example, it will enable prompt discovery of a uniform and all-sky sample of galaxy/group/cluster-scale lenses that will underpin LSST strong-lensing science. This white paper complements submissions from DESC, SLSC, and TVSSC, that discuss kilonova, GW, and strong lensing.

We modify the chemo-dynamical code GEAR to simulate the impact of self-interacting dark matter on the observable quantities of 19 low mass dwarf galaxies with a variety star forming properties. We employ a relatively high, velocity independent cross-section of $\sigma/m = 10$cm$^2$/g and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-sight velocities, the chemical abundance and the star formation history. We find that despite the creation of large cores at the centre of the dark matter haloes, the impact of SIDM on the observable quantities of quenched galaxies is indiscernible, dominated mostly by the stochastic build up of the stellar matter. As such we conclude that it is impossible to make global statements on the density profile of dwarf galaxies from single or small samples. Although based mostly on quenched galaxies, this finding supports other recent work putting into question the reliability of inferred cored density profiles that are derived from observed line-of-sight velocities.

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

The Dark Energy Spectroscopic Instrument (DESI) is under construction to measure the expansion history of the universe using the baryon acoustic oscillations technique. The spectra of 35 million galaxies and quasars over 14,000 square degrees will be measured during a 5-year survey. A new prime focus corrector for the Mayall telescope at Kitt Peak National Observatory will deliver light to 5,000 individually targeted fiber-fed robotic positioners. The fibers in turn feed ten broadband multi-object spectrographs. We describe the ProtoDESI experiment, that was installed and commissioned on the 4-m Mayall telescope from August 14 to September 30, 2016. ProtoDESI was an on-sky technology demonstration with the goal to reduce technical risks associated with aligning optical fibers with targets using robotic fiber positioners and maintaining the stability required to operate DESI. The ProtoDESI prime focus instrument, consisting of three fiber positioners, illuminated fiducials, and a guide camera, was installed behind the existing Mosaic corrector on the Mayall telescope. A Fiber View Camera was mounted in the Cassegrain cage of the telescope and provided feedback metrology for positioning the fibers. ProtoDESI also provided a platform for early integration of hardware with the DESI Instrument Control System that controls the subsystems, provides communication with the Telescope Control System, and collects instrument telemetry data. Lacking a spectrograph, ProtoDESI monitored the output of the fibers using a Fiber Photometry Camera mounted on the prime focus instrument. ProtoDESI was successful in acquiring targets with the robotically positioned fibers and demonstrated that the DESI guiding requirements can be met.

Discovery of strongly-lensed gravitational wave (GW) sources will unveil binary compact objects at higher redshifts and lower intrinsic luminosities than is possible without lensing. Such systems will yield unprecedented constraints on the mass distribution in galaxy clusters, measurements of the polarization of GWs, tests of General Relativity, and constraints on the Hubble parameter. Excited by these prospects, and intrigued by the presence of so-called "heavy black holes" in the early detections by LIGO-Virgo, we commenced a search for strongly-lensed GWs and possible electromagnetic counterparts in the latter stages of the second LIGO observing run (O2). Here, we summarise our calculation of the detection rate of strongly-lensed GWs, describe our review of BBH detections from O1, outline our observing strategy in O2, summarize our follow-up observations of GW170814, and discuss the future prospects of detection.

We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

Simulations of self-interacting dark matter (SIDM) predict that dark matter should lag behind galaxies during a collision. If the interaction is mediated by a high-mass force carrier, the distribution of dark matter can also develop asymmetric dark matter tails. To search for this asymmetry, we compute the gravitational lensing properties of a mass distribution with a free {\em skewness} parameter. We apply this to the dark matter around the four central galaxies in cluster Abell~3827. In the galaxy whose dark matter peak has previously been found to be offset, we tentatively measure a skewness $s=0.23^{+0.05}_{-0.22}$ in the same direction as the peak offset. Our method may be useful in future gravitational lensing analyses of colliding galaxy clusters and merging galaxies.

We investigate how self-interacting dark matter (SIDM) with anisotropic scattering affects the evolution of isolated dark matter haloes as well as systems with two colliding haloes. For isolated haloes, we find that the evolution can be adequately captured by treating the scattering as isotropic, as long as the isotropic cross-section is appropriately matched to the underlying anisotropic model. We find that this matching should not be done using the momentum transfer cross-section, as has been done previously. Matching should instead be performed via a modified momentum transfer cross-section that takes into account that dark matter particles can be relabelled after they scatter, without altering the dynamics. However, using cross-sections that are matched to give the same behaviour in isolated haloes, we find that treating dark matter scattering as isotropic under-predicts the effects of anisotropic dark matter scattering when haloes collide. In particular, the DM-galaxy offset induced by SIDM in colliding galaxy clusters is larger when we simulate the underlying particle model, than if we use a matched isotropic model. On the other hand, well motivated particle models with anisotropic scattering typically have cross-sections with a strong velocity dependence, and we discover a previously unrecognised effect that suppresses DM-galaxy offsets in colliding clusters making it hard for these systems to provide competitive constraints on such particle models.

We perform numerical simulations of the merging galaxy cluster 1E 0657-56 (the Bullet Cluster), including the effects of elastic dark matter scattering. In a similar manner to the stripping of gas by ram pressure, dark matter self-interactions would transfer momentum between the two galaxy cluster dark matter haloes, causing them to lag behind the collisionless galaxies. The absence of an observed separation between the dark matter and stellar components in the Bullet Cluster has been used to place upper limits on the cross-section for dark matter scattering. We emphasise the importance of analysing simulations in an observationally-motivated manner, finding that the way in which the positions of the various components are measured can have a larger impact on derived constraints on dark matter's self-interaction cross-section than reasonable changes to the initial conditions for the merger. In particular, we find that the methods used in previous studies to place some of the tightest constraints on this cross-section do not reflect what is done observationally, and overstate the Bullet Cluster's ability to constrain the particle properties of dark matter. We introduce the first simulations of the Bullet Cluster including both self-interacting dark matter and gas. We find that as the gas is stripped it introduces radially-dependent asymmetries into the stellar and dark matter distributions. As the techniques used to determine the positions of the dark matter and galaxies are sensitive to different radial scales, these asymmetries can lead to erroneously measured offsets between dark matter and galaxies even when they are spatially coincident.

We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution.

If dark matter interacts, even weakly, via non-gravitational forces, simulations predict that it will be preferentially scattered towards the trailing edge of the halo during collisions between galaxy clusters. This will temporarily create a non-symmetric mass profile, with a trailing over-density along the direction of motion. To test this hypothesis, we fit (and subtract) symmetric halos to the weak gravitational data of 72 merging galaxy clusters observed with the Hubble Space Telescope. We convert the shear directly into excess {\kappa} and project in to a one dimensional profile. We generate numerical simulations and find that the one dimensional profile is well described with simple Gaussian approximations. We detect the weak lensing signal of trailing gas at a 4{\sigma} confidence, finding a mean gas fraction of Mgas/Mdm = 0.13 +/- 0.035. We find no evidence for scattered dark matter particles with a estimated scattering fraction of f = 0.03 +/- 0.05. Finally we find that if we can reduce the statistical error on the positional estimate of a single dark matter halo to <2.5", then we will be able to detect a scattering fraction of 10% at the 3{\sigma} level with current surveys. This poten- tially interesting new method can provide an important independent test for other complimentary studies of the self-interaction cross-section of dark matter.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

Transits of exoplanets observed in the near-UV have been used to study the scattering properties of their atmospheres and possible star-planet interactions. We observed the primary transits of 15 exoplanets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-16b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-12b, WASP-33b, WASP-36b, WASP-44b, WASP-48b, and WASP-77Ab) in the near-UV and several optical photometric bands to update their planetary parameters, ephemerides, search for a wavelength dependence in their transit depths to constrain their atmospheres, and determine if asymmetries are visible in their light curves. Here we present the first ground-based near-UV light curves for 12 of the targets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-33b, WASP-36b, WASP-48b, and WASP-77Ab). We find that none of the near-UV transits exhibit any non-spherical asymmetries, this result is consistent with recent theoretical predictions by Ben-Jaffel et al. and Turner et al. The multi-wavelength photometry indicates a constant transit depth from near-UV to optical wavelengths in 10 targets (suggestive of clouds), and a varying transit depth with wavelength in 5 targets (hinting at Rayleigh or aerosol scattering in their atmospheres). We also present the first detection of a smaller near-UV transit depth than that measured in the optical in WASP-1b and a possible opacity source that can cause such radius variations is currently unknown. WASP-36b also exhibits a smaller near-UV transit depth at 2.6$\sigma$. Further observations are encouraged to confirm the transit depth variations seen in this study.

We use the "Evolution and Assembly of GaLaxies and their Environments" ( EAGLE ) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent $\epsilon = 700$ pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxy's direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with $\sigma \leq 196$ pc. Since astrophysical effects produce no feasible analogues for the $1.62^{+0.47}_{-0.49}$ kpc offset recently observed in Abell 3827, the observational result is in tension with the collisionless cold dark matter model assumed in our simulations.

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

Recently, Lattelais et al. (2009) have interpreted aggregated observations of molecular isomers to suggest that there exists a "minimum energy principle'', such that molecular formation will favor more stable molecular isomers for thermodynamic reasons. To test the predictive power of this principle, we have fully characterized the spectra of the three isomers of C$_{3}$H$_{2}$O toward the well known molecular region Sgr B2(N). Evidence for the detection of the isomers cyclopropenone (c-C$_{3}$H$_{2}$O) and propynal (HCCCHO) is presented, along with evidence for the non-detection of the lowest zero-point energy isomer, propadienone (CH$_2$CCO). We interpret these observations as evidence that chemical formation pathways, which may be under kinetic control, have a more pronounced effect on final isomer abundances than thermodynamic effects such as the minimum energy principle.

We observed nine primary transits of the hot Jupiter TrES-3b in several optical and near-UV photometric bands from 2009 June to 2012 April in an attempt to detect its magnetic field. Vidotto, Jardine and Helling suggest that the magnetic field of TrES-3b can be constrained if its near-UV light curve shows an early ingress compared to its optical light curve, while its egress remains unaffected. Predicted magnetic field strengths of Jupiter-like planets should range between 8 G and 30 G. Using these magnetic field values and an assumed B_star of 100 G, the Vidotto et al. method predicts a timing difference of 5-11 min. We did not detect an early ingress in our three nights of near-UV observations, despite an average cadence of 68 s and an average photometric precision of 3.7 mmag. However, we determined an upper limit of TrES-3b's magnetic field strength to range between 0.013 and 1.3 G (for a 1-100 G magnetic field strength range for the host star, TrES-3) using a timing difference of 138 s derived from the Nyquist-Shannon sampling theorem. To verify our results of an abnormally small magnetic field strength for TrES-3b and to further constrain the techniques of Vidotto et al., we propose future observations of TrES-3b with other platforms capable of achieving a shorter near-UV cadence. We also present a refinement of the physical parameters of TrES-3b, an updated ephemeris and its first published near-UV light curve. We find that the near-UV planetary radius of Rp = 1.386+0.248-0.144 RJup is consistent with the planet's optical radius.

Using high spatial resolution Hubble Space Telescope/Advanced Camera for Surveys archival imaging observations of Arp 116, centred on the elliptical galaxy NGC 4649, we explore the novel technique of pixel-by-pixel analysis of the galaxy's colour-magnitude diagramme to search for any evidence of recent enhanced star formation due to the apparent tidal interaction with its spiral companion, NGC 4647. From a detailed analysis of the system's geometry, and based on additional circumstantial evidence from extant multi-wavelength observations, we conclude that, while there may be grounds for the tidal-interaction assumption for this system, any interaction has thus far been of insufficient strength to trigger an enhanced level of recent star formation in the elliptical component, although close inspection of our colour images shows a faint excess of bluer pixels (a ~0.20 mag bluer "loop") in the elliptical galaxy on the side of the spiral companion. Given that there appears to be a moderate reservoir of available gas for ongoing star formation (although at low column density), this suggests that we are currently witnessing the onset of the tidal interaction between NGC 4647 and NGC 4649. In addition, the triggering of new star formation in NGC 4649 may be significantly impeded due to the much lower mass of the spiral component.