Gravitational waves (GWs) are subject to gravitational lensing in the same way as electromagnetic radiation. However, to date, no unequivocal observation of a lensed GW transient has been reported. Independently, GW observatories continue to search for the stochastic GW signal which is produced by many transient events at high redshift. We exploit a surprising connection between the lensing of individual transients and limits to the background radiation produced by the unresolved population of binary back hole mergers: we show that it constrains the fraction of individually resolvable lensed binary black holes to less than $\sim 3\times 10^{-5}$ at present sensitivity. We clarify the interpretation of existing, low redshift GW observations (obtained assuming no lensing) in terms of their apparent lensed redshifts and masses and explore constraints from GW observatories at future sensitivity.

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.

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 explore the connection between dust and star formation, in the context of environmental effects on galaxy evolution. In particular, we exploit the susceptibility of dust to external processes to assess the influence of dense environment on star-forming galaxies. We have selected cluster Abell 1758 from the Local Cluster Substructure Survey (LoCuSS). Its complex dynamical state is an ideal test-bench to track dust removal and destruction in galaxies due to merger and accretion shocks. We present a systematic panchromatic study (from 0.15 $\rm \mu$m with GALEX to 500 $\rm \mu$m with Herschel) of spectroscopically confirmed star-forming cluster galaxies at intermediate redshift. We observe that the main subclusters (A1758N and A1758S) belong to two separate large-scale structures, with no overlapping galaxy members. Star-forming cluster members are distributed preferentially outside cluster central regions, and are not grouped in substructures. Rather, these galaxies are being funneled towards the main subclusters along separate accretion filaments. Additionally, we present the first study of dust-to-stellar (DTS) mass ratio used as indicator for local environmental influence on galaxy evolution. Star-forming cluster members show lower mean values (32% at 2.4$\rm \sigma$) of DTS mass ratio and lower levels of infrared emission from birth clouds with respect to coeval star-forming field galaxies. This picture is consistent with the majority of star-forming cluster members infalling in isolation. Upon accretion, star-formation is observed to decrease and warm dust is destroyed due to heating from the intracluster medium radiation, ram-pressure stripping and merger shocks.

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}}$.

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.

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 report direct evidence of pre-processing of the galaxies residing in galaxy groups falling into galaxy clusters drawn from the Local Cluster Substructure Survey (LoCuSS). 34 groups have been identified via their X-ray emission in the infall regions of 23 massive ($\rm \langle M_{200}\rangle = 10^{15}\,M_{\odot}$) clusters at $0.15<z<0.3$. Highly complete spectroscopic coverage combined with 24 $\rm\mu$m imaging from Spitzer allows us to make a consistent and robust selection of cluster and group members including star forming galaxies down to a stellar mass limit of $\rm M_{\star} = 2\times10^{10}\,M_{\odot}$. The fraction $\rm f_{SF}$ of star forming galaxies in infalling groups is lower and with a flatter trend with respect to clustercentric radius when compared to the rest of the cluster galaxy population. At $\rm R\approx1.3\,r_{200}$ the fraction of star forming galaxies in infalling groups is half that in the cluster galaxy population. This is direct evidence that star formation quenching is effective in galaxies already prior to them settling in the cluster potential, and that groups are favourable locations for this process.

Galaxy clusters are expected to form hierarchically in a LCDM universe, growing primarily through mergers with lower mass clusters and the continual accretion of group-mass halos. Galaxy clusters assemble late, doubling their masses since z~0.5, and so the outer regions of clusters should be replete with infalling group-mass systems. We present an XMM-Newton survey to search for X-ray groups in the infall regions of 23 massive galaxy clusters at z~0.2, identifying 39 X-ray groups that have been spectroscopically confirmed to lie at the cluster redshift. These groups have mass estimates in the range 2x10^13-7x10^14Msun, and group-to-cluster mass ratios as low as 0.02. The comoving number density of X-ray groups in the infall regions is ~25x higher than that seen for isolated X-ray groups from the XXL survey. The average mass per cluster contained within these X-ray groups is 2.2x10^14Msun, or 19% of the mass within the primary cluster itself. We estimate that ~10^15Msun clusters increase their masses by 16% between z=0.223 and the present day due to the accretion of groups with M200>10^13.2Msun. This represents about half of the expected mass growth rate of clusters at these late epochs. The other half is likely to come from smooth accretion of matter not bound in halos. The mass function of the infalling X-ray groups appears significantly top-heavy with respect to that of field X-ray systems, consistent with expectations from numerical simulations, and the basic consequences of collapsed massive dark matter halos being biased tracers of the underlying large-scale density distribution.

We investigate the possible presence of active galactic nuclei (AGN) in dwarf galaxies and other nearby galaxies to identify candidates for follow-up confirmation and dynamical mass measurements. We use the Wide-field Infrared Survey Explorer (WISE) All-Sky Release Source Catalog and examine the infrared colours of a sample of dwarf galaxies and other nearby galaxies in order to identify both unobscured and obscured candidate AGN by applying the infrared colour diagnostic. Stellar masses of galaxies are obtained using a combination of three independent methods. Black hole masses are estimated using the bolometric luminosity of the AGN candidates and computed for three cases of the bolometric-to-Eddington luminosity ratio. We identify 303 candidate AGN, of which 276 were subsequently found to have been independently identified as AGN via other methods. The remaining 9% require follow-up observations for confirmation. The activity is detected in galaxies with stellar masses from ~ 10^6 to 10^9 solar masses; assuming the candidates are AGN, the black hole masses are estimated to be ~ 10^3 - 10^6 solar masses, adopting L_bol = 0.1 L_Edd. The black hole masses probed are several orders of magnitude smaller than previously reported for centrally located massive black holes. We examine the stellar mass versus black hole mass relationship in this low galaxy mass regime. We find that it is consistent with the existing relation extending linearly (in log-log space) into the lower mass regime. These findings suggest that CMBH are present in low-mass galaxies and in the Local Universe, and provide new impetus for follow-up dynamical studies of quiescent black holes in local dwarf galaxies.

We present SCUBA-2 follow-up of 61 candidate high-redshift Planck sources. Of these, 10 are confirmed strong gravitational lenses and comprise some of the brightest such submm sources on the observed sky, while 51 are candidate proto-cluster fields undergoing massive starburst events. With the accompanying Herschel-SPIRE observations and assuming an empirical dust temperature prior of $34^{+13}_{-9}$ K, we provide photometric redshift and far-IR luminosity estimates for 172 SCUBA-2-selected sources within these Planck overdensity fields. The redshift distribution of the sources peak between a redshift of 2 and 4, with one third of the sources having $S_{500}$/$S_{350} > 1$. For the majority of the sources, we find far-IR luminosities of approximately $10^{13}\,\mathrm{L}_\odot$, corresponding to star-formation rates of around $1000$ M$_\odot \mathrm{yr}^{-1}$. For $S_{850}>8$ mJy sources, we show that there is up to an order of magnitude increase in star-formation rate density and an increase in uncorrected number counts of $6$ for $S_{850}>8$ mJy when compared to typical cosmological survey fields. The sources detected with SCUBA-2 account for only approximately $5$ per cent of the Planck flux at 353 GHz, and thus many more fainter sources are expected in these fields.

We present a study of the star formation and central black hole accretion activity of the galaxies hosted in the two nearby (z$\sim$0.2) rich galaxy clusters Abell 983 and 1731. Aims: We are able to quantify both the obscured and unobscured star formation rates, as well as the presence of active galactic nuclei (AGN) as a function of the environment in which the galaxy is located. Methods: We targeted the clusters with unprecedented deep infrared Spitzer observations (0.2 mJy @ 24 micron), near-IR Palomar imaging and optical WIYN spectroscopy. The extent of our observations ($\sim$ 3 virial radii) covers the vast range of possible environments, from the very dense cluster centre to the very rarefied cluster outskirts and accretion regions. Results: The star forming members of the two clusters present star formation rates comparable with those measured in coeval field galaxies. The analysis of the spatial arrangement of the spectroscopically confirmed members reveals an elongated distribution for A1731 with respect to the more uniform distribution of A983. The emerging picture is compatible with A983 being a fully evolved cluster, in contrast with the still actively accreting A1731. Conclusions: The analysis of the specific star formation rate reveals evidence of on-going galaxy pre-processing along A1731's filament-like structure. Furthermore, the decrease in the number of star forming galaxies and AGN towards the cluster cores suggests that the cluster environment is accelerating the ageing process of galaxies and blocking further accretion of the cold gas that fuels both star formation and black hole accretion activity.

Clusters of galaxies are embedded in halos of optically thin, gravitationally stratified, weakly magnetized plasma at the system's virial temperature. Due to radiative cooling and anisotropic heat conduction, such intracluster medium (ICM) is subject to local instabilities, which are combinations of the thermal, magnetothermal and heat-flux-driven buoyancy instabilities. If the ICM rotates significantly, its stability properties are substantially modified and, in particular, also the magnetorotational instability (MRI) can play an important role. We study simple models of rotating cool-core clusters and we demonstrate that the MRI can be the dominant instability over significant portions of the clusters, with possible implications for the dynamics and evolution of the cool cores. Our results give further motivation for measuring the rotation of the ICM with future X-ray missions such as ASTRO-H and ATHENA.

We present the results of a study of a statistically significant sample of galaxies which clearly demonstrate that supermassive black holes are generically present in all morphological types. Our analysis is based on the quantitative morphological classification of 1.12 million galaxies in the SDSS DR7 and on the detection of black hole activity via two different methods, the first one based on their X-ray/radio emission and the second one based on their mid-infrared colors. The results of the first analysis confirm the correlation between black hole and total stellar mass for 8 galaxies and includes one galaxy classified as bulgeless. The results of our second analysis, consisting of 15,991 galaxies, show that galaxies hosting a supermassive black hole follow the same morphological distribution as the general population of galaxies in the same redshift range. In particular, the fraction of bulgeless galaxies, 1,450 galaxies or 9 percent, is found to be the same as in the general population. We also present the correlation between black hole and total stellar mass for 6,247 of these galaxies. Importantly, whereas previous studies were limited to primarily bulge-dominated systems, our study confirms this relationship to all morphological types, in particular, to 530 bulgeless galaxies. Our results indicate that the true correlation that exists for supermassive black holes and their host galaxies is between the black hole mass and the total stellar mass of the galaxy and hence, we conclude that the previous assumption that the black hole mass is correlated with the bulge mass is only approximately correct.

We study simple models of massive galaxy clusters in which the intracluster medium (ICM) rotates differentially in equilibrium in the cluster gravitational potential. We obtain the X-ray surface brightness maps, evaluating the isophote flattening due to the gas rotation. Using a set of different rotation laws, we put constraint on the amplitude of the rotation velocity, finding that rotation curves with peak velocity up to \sim 600 km s^-1 are consistent with the ellipticity profiles of observed clusters. We convolve each of our models with the instrument response of the X-ray Calorimeter Spectrometer on board the ASTRO-H to calculate the simulated X-ray spectra at different distance from the X-ray centre. We demonstrate that such an instrument will allow us to measure rotation of the ICM in massive clusters, even with rotation velocities as low as \sim 100 km s^-1