We present the first satellite system of the Large Binocular Telescope Satellites Of Nearby Galaxies Survey (LBT-SONG), a survey to characterize the close satellite populations of Large Magellanic Cloud to Milky Way-mass, star-forming galaxies in the Local Volume. In this paper, we describe our satellite finding and completeness measurement methodology and apply this framework to NGC 628, an isolated galaxy with ~1/4 the stellar mass of the Milky Way. We present two new dwarf satellite galaxy candidates: NGC 628 dwA, and dwB with $\text{M}_{\text{V}}$ = $-12.2$ and $-7.7$, respectively. NGC 628 dwA is a classical dwarf while NGC 628 dwB is a low-luminosity galaxy that appears to have been quenched after reionization. Completeness corrections indicate that the presence of these two satellites is consistent with CDM predictions. The satellite colors indicate that the galaxies are neither actively star-forming nor do they have the purely ancient stellar populations characteristic of ultrafaint dwarfs. Instead, and consistent with our previous work on the NGC 4214 system, they show signs of recent quenching, further indicating that environmental quenching can play a role in modifying satellite populations even for hosts smaller than the Milky Way.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the signal of dark matter structure. An alternate approach is to measure quasar nuclear-narrow line emission, which is free from microlensing and present in virtually all quasar lenses. In this paper, we double the number of systems which can be used for gravitational lensing analyses by presenting measurements of narrow-line emission from a sample of 8 quadruply imaged quasar lens systems, WGD J0405-3308, HS 0810+2554, RX J0911+0551, SDSS J1330+1810, PS J1606-2333, WFI 2026-4536, WFI 2033-4723 and WGD J2038-4008. We describe our updated grism spectral modelling pipeline, which we use to measure narrow-line fluxes with uncertainties of 2-10\%, presented here. We fit the lensed image positions with smooth mass models and demonstrate that these models fail to produce the observed distribution of image fluxes over the entire sample of lenses. Furthermore, typical deviations are larger than those expected from macromodel uncertainties. This discrepancy indicates the presence of perturbations caused by small-scale dark matter structure. The interpretation of this result in terms of dark matter models is presented in a companion paper.

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter halos on sub-galactic scales. Using the image positions and flux ratios from eight quadruply-imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhalos and halos along the line of sight, and account for warm dark matter (WDM) free-streaming effects on the mass function and mass-concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated halos, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parameterizing dark matter free-streaming through the half-mode mass $m_{\rm{hm}}$, we constrain the thermal relic particle mass $m_{\rm{DM}}$ corresponding to $m_{\rm{hm}}$. At $95 \%$ CI: $m_{\rm{hm}} < 10^{7.8} M_{\odot}$ ($m_{\rm{DM}} > 5.2 \ \rm{keV}$). We disfavor $m_{\rm{DM}} = 4.0 \rm{keV}$ and $ m_{\rm{DM}} = 3.0 \rm{keV}$ with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrain the projected mass in substructure between $10^6 - 10^{9} M_{\odot}$ near lensed images. At $68 \%$ CI, we infer $2.0 - 6.1 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to mean projected mass fraction $\bar{f}_{\rm{sub}} = 0.035_{-0.017}^{+0.021}$. At $95 \%$ CI, we obtain a lower bound on the projected mass of $0.6 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to $\bar{f}_{\rm{sub}} > 0.005$. These results agree with the predictions of cold dark matter.

The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the image positions and flux ratios from eleven quadruple-image strong gravitational lenses (quads) in the mass range $10^{6} - 10^{10} M_{\odot}$, assuming cold dark matter. Our analysis framework includes both subhalos and line of sight halos, marginalizes over nuisance parameters describing the lens macromodel, accounts for finite source effects on lensing observables, and simultaneously constrains the normalization and logarithmic slope of the mass-concentration relation, and the normalization of the subhalo mass function. At $z=0$, we constrain the concentration of $10^{8} M_{\odot}$ halos $c=12_{-5}^{+6}$ at $68 \%$ CI, and $c=12_{-9}^{+15}$ at $95 \%$ CI. For a $10^{7} M_{\odot}$ halo, we obtain $68 \%$ ($95 \%$) constraints $c=15_{-8}^{+9}$ ($c=15_{-11}^{+18}$), while for $10^{9} M_{\odot}$ halos $c=10_{-4}^{+7}$ ($c=10_{-7}^{+14}$). These results are consistent with the theoretical predictions from mass-concentration relations in the literature, and establish strong lensing by galaxies as a powerful probe of halo concentrations on sub-galactic scales across cosmological distance.

Strong lensing provides a powerful means of investigating the nature of dark matter as it probes dark matter structure on sub-galactic scales. We present an extension of a forward modeling framework that uses flux ratios from quadruply imaged quasars (quads) to measure the shape and amplitude of the halo mass function, including line of sight (LOS) halos and main deflector subhalos. We apply this machinery to 50 mock lenses --- roughly the number of known quads --- with warm dark matter (WDM) mass functions exhibiting free-streaming cutoffs parameterized by the half-mode mass $m_{\rm{hm}}$. Assuming cold dark matter (CDM), we forecast bounds on $m_{\rm{hm}}$ and the corresponding thermal relic particle masses over a range of tidal destruction severity, assuming a particular WDM mass function and mass-concentration relation. With significant tidal destruction, at $2 \sigma$ we constrain $m_{\rm{hm}}<10^{7.9} \left(10^{8.4}\right) M_{\odot}$, or a 4.4 (3.1) keV thermal relic, with image flux uncertainties from measurements and lens modeling of $2\% \left(6\%\right)$. With less severe tidal destruction we constrain $m_{\rm{hm}}<10^{7} \left(10^{7.4}\right) M_{\odot}$, or an 8.2 (6.2) keV thermal relic. If dark matter is warm, with $m_{\rm{hm}} = 10^{7.7} M_{\odot}$ (5.1 keV), we would favor WDM with $m_{\rm{hm}} > 10^{7.7} M_{\odot}$ over CDM with relative likelihoods of 22:1 and 8:1 with flux uncertainties of $2\%$ and $6\%$, respectively. These bounds improve over those obtained by modeling only main deflector subhalos because LOS objects produce additional flux perturbations, especially for high redshift systems. These results indicate that $\sim 50$ quads can conclusively differentiate between warm and cold dark matter.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

We test the ability of strong lensing data to constrain the size of a central core in the dark matter halos of galaxy clusters, using Abell 611 as a prototype. Using simulated data, we show that modeling a cluster halo with ellipticity in the gravitational potential can bias the inferred mass and concentration, which may bias the inferred central density when weak lensing or X-ray data are added. We also the highlight the possibility for spurious constraints on the core size if the radial density profile is different from the assumed model. These systematics can be ameliorated if central images are present in the data. Applying our methodology to Abell 611 and imposing a reasonable prior on the stellar mass-to-light ratio restricts the core size to be less than about 4 kpc, with a minimum reduced $\chi^2$ of 0.28 for 0."2 positional errors. Such small cores imply a constraint on the dark matter self-interaction cross section of the order of $0.1\ \mathrm{cm^2/g}$ at relative velocities of about 1500 km/s.

We consider AGN photoionization of gas in companion galaxies (cross-ionization) as a way to sample the intensity of AGN radiation in both direction and time, independent of the gas properties of the AGN host galaxies. From an initial set of 212 AGN+companion systems, identified with the help of Galaxy Zoo participants, we obtained long-slit optical spectra of 32 pairs which were a priori likely to show cross-ionization based on projected separation or angular extent of the companion. From emission-line ratios, 10 of these systems are candidates for cross-ionization, roughly the fraction expected if most AGN have ionization cones with 70-degree opening angles. Among these, Was 49 remains the strongest nearby candidate. NGC 5278/9 and UGC 6081 are dual-AGN systems with tidal debris, complicating identification of cross-ionization. The two weak AGN in the NGC 5278/9 system ionize gas filaments to a projected radius 14 kpc from each galaxy. In UGC 6081, an irregular high-ionization emission region encompasses both AGN, extending more than 15 kpc from each. The observed AGN companion galaxies with and without signs of external AGN photoionization have similar distributions in estimated incident AGN flux, suggesting that geometry of escaping radiation or long-term variability control this facet of the AGN environment. This parallels conclusions for luminous QSOs based on the proximity effect among Lyman-alpha absorbers. In some galaxies, mismatch between spectroscopic classifications in the common BPT diagram and the intensity of weaker He II and [Ne V] emission lines highlights the limits of common classifications in low-metallicity environments.

The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parameterization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free streaming scale $m_{\rm{hm}}$. We forecast bounds on dark matter warmth for 120-180 lenses, attainable with future surveys, at typical lens (source) redshifts of 0.5 (1.5) in early-type galaxies with velocity dispersions of 220-270 km/sec. We demonstrate that limits on $m_{\rm{hm}}$ deteriorate rapidly with increasing uncertainty in image fluxes, underscoring the importance of precise measurements and accurate lens models. For our forecasts, we assume the deflectors in the lens sample do not exhibit complex morphologies, so we neglect systematic errors in their modeling. Omitting the additional signal from line of sight halos, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding the destruction of all subhalos within the host scale radius, we forecast $2\sigma$ bounds on $m_{\rm{hm}}$ (thermal relic mass) of $10^{7.5} \ (5.0)$, $10^{8} \ (3.6)$, and $10^{8.5} \ (2.7) \ M_{\odot} \left(\rm{keV}\right)$ for flux errors of $2\%$, $4\%$, and $8\%$. With a higher normalization, these constraints improve to $10^{7.2} \ (6.6)$, $10^{7.5} \ (5.3) $, and $10^{7.8} \ (4.3) \ M_{\odot} \left(\rm{keV}\right)$ with 120 systems. We are also able to measure the normalization of the mass function, which has implications for baryonic feedback models and tidal stripping.

The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years and constrain the systematic uncertainties in non-varying parameters such as the black hole mass. We find that the BLR geometry of a thick disk viewed close to face-on is stable over this time, although the size of the BLR grows by a factor of $\sim 2$. The dynamics of the BLR are dominated by inflow and the inferred black hole mass is consistent for the three datasets, despite the increase in BLR size. Combining the inference for the three datasets yields a black hole mass and statistical uncertainty of $\log_{10}($M$_{\rm BH}/\rm{M}_{\odot})=6.82^{+0.09}_{-0.09}$ with a standard deviation in individual measurements of 0.13 dex.

A fundamental prediction of the cold dark matter cosmology is the existence of a large number of dark subhalos around galaxies, most of which should be entirely devoid of stars. Confirming the existence of dark substructures stands among the most important empirical challenges in modern cosmology: if they are found and quantified with the mass spectrum expected, then this would close the door on a vast array of competing theories. But in order for observational programs of this kind to reach fruition, we need robust predictions. Here we explore substructure predictions for lensing using galaxy lens-like hosts at z=0.2 from the Illustris simulations both in full hydrodynamics and dark matter only. We quantify substructures more massive than ~ 10^9 M_sun, comparable to current lensing detections derived from HST, Keck, and ALMA. The addition of full hydrodynamics reduces the overall subhalo mass function by about a factor of two. Even for the dark matter only runs, most (~ 85%) lines of sight through projected cylinders of size close to an Einstein radius contain no substructures larger than 10^9 M_sun. The fraction of empty sight lines rises to ~ 95% in full physics simulations. This suggests we will likely need hundreds of strong lensing systems suitable for substructure studies, as well as predictions that include the effects of baryon physics on substructure, to properly constrain cosmological models. Fortunately, the field is poised to fulfill these requirements.

Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant radio emission which is extended enough avoid significant lensing by stars in the plane of the lens galaxy, and red enough to be minimally affected by differential dust extinction. Narrow-line emission is a promising alternative as it is also extended and, unlike radio, detectable in virtually all optically selected AGN lenses. We present first results from a WFC3 grism narrow-line survey of lensed quasars, for the quadruply lensed AGN HE0435-1223. Using a forward modelling pipeline which enables us to robustly account for spatial blending, we measure the [OIII] 5007 \AA~ flux ratios of the four images. We find that the [OIII] fluxes and positions are well fit by a simple smooth mass model for the main lens. Our data rule out a $M_{600}>10^{8} (10^{7.2}) M_\odot$ NFW perturber projected within $\sim$1\farcs0 (0\farcs1) arcseconds of each of the lensed images, where $M_{600}$ is the perturber mass within its central 600 pc. The non-detection is broadly consistent with the expectations of $\Lambda$CDM for a single system. The sensitivity achieved demonstrates that powerful limits on the nature of dark matter can be obtained with the analysis of $\sim20$ narrow-line lenses.

The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data.

The empirical correlation between the mass of a super-massive black hole (MBH) and its host galaxy properties is widely considered to be evidence of their co-evolution. A powerful way to test the co-evolution scenario and learn about the feedback processes linking galaxies and nuclear activity is to measure these correlations as a function of redshift. Unfortunately, currently MBH can only be estimated in active galaxies at cosmological distances. At these distances, bright active galactic nuclei (AGN) can outshine the host galaxy, making it extremely difficult to measure the host's luminosity. Strongly lensed AGNs provide in principle a great opportunity to improve the sensitivity and accuracy of the host galaxy luminosity measurements as the host galaxy is magnified and more easily separated from the point source, provided the lens model is sufficiently accurate. In order to measure the MBH-L correlation with strong lensing, it is necessary to ensure that the lens modelling is accurate, and that the host galaxy luminosity can be recovered to at least a precision and accuracy better than that of the typical MBH measurement. We carry out extensive and realistic simulations of deep Hubble Space Telescope observations of lensed AGNs obtained by our collaboration. We show that the host galaxy luminosity can be recovered with better accuracy and precision than the typical uncertainty on MBH(~ 0.5 dex) for hosts as faint as 2-4 magnitudes dimmer than the AGN itself. Our simulations will be used to estimate bias and uncertainties on the actual measurements to be presented in a future paper.

It is widely believed that the large discrepancy between the observed number of satellite galaxies and the predicted number of dark subhalos can be resolved via a variety of baryonic effects which suppress star formation in low mass halos.Supporting this hypothesis, numerous high resolution simulations with star formation, and associated feedback have been shown to reproduce the satellite luminosity function around Milky Way-mass simulated galaxies at redshift zero. However, a more stringent test of these models is their ability to simultaneously match the satellite luminosity functions of a range of host halo masses and redshifts. In this work we measure the luminosity function of faint (sub-Small Magellanic Cloud luminosity) satellites around hosts with stellar masses 10.5$<\log_{10}$M$_*$/M$_\odot<11.5$ to an unprecedented redshift of 1.5. This new measurement of the satellite luminosity function provides powerful new constraining power; we compare these results with predictions from four different simulations and show that although the models perform similarly over-all, no one model reproduces the satellite luminosity function reliably at all redshifts and host stellar masses. This result highlights the continued need for improvement in understanding the fundamental physics that governs satellite galaxy evolution.

In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the AGN continuum, and measurements of the broad H-beta line widths in mean and root-mean square (rms) spectra. For the most highly variable AGNs we also measured broad H-beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H-beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H-beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593 the broad H-beta velocity shifted by ~250 km/s over a one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise.

Strong gravitational lenses can be used to detect low mass subhalos, based on deviations in image fluxes and positions from what can be achieved with a smooth mass distribution. So far, this method has been limited by the small number of (radio-loud, microlensing free) systems which can be analysed for the presence of substructure. Using the gravitational lens B1422+231, we demonstrate that adaptive optics integral field spectroscopy can also be used to detect dark substructures. We analyse data obtained with OSIRIS on the Keck I Telescope, using a Bayesian method that accounts for uncertainties relating to the point spread function and image positions in the separate exposures. The narrow-line [OIII] fluxes measured for the lensed images are consistent with those measured in the radio, and show a significant deviation from what would be expected in a smooth mass distribution, consistent with the presence of a perturbing low mass halo. Detailed lens modelling shows that image fluxes and positions are fit significantly better when the lens is modelled as a system containing a single perturbing subhalo in addition to the main halo, rather than by the main halo on its own, indicating the significant detection of substructure.The inferred mass of the subhalo depends on the subhalo mass density profile: the 68% confidence interval for the perturber mass within 600 pc are: 8.2$^{+0.6}_{-0.8}$, 8.2$^{+0.6}_{-1}$ and 7.6$\pm0.3$ $\log_{10}$[M$_{\rm{sub}}$/M$_\odot$] respectively for a singular isothermal sphere, a pseudo-Jaffe, and an NFW mass profile. This method can extend the study of flux ratio anomalies to virtually all quadruply imaged quasars, and therefore offers great potential to improve the determination of the subhalo mass function in the near future.

Strong gravitational lensing can be used to directly measure the mass function of their satellites, thus testing one of the fundamental predictions of cold dark matter cosmological models. Given the importance of this test it is essential to ensure that galaxies acting as strong lenses have dark and luminous satellites which are representative of the overall galaxy population. We address this issue by measuring the number and spatial distribution of luminous satellites in ACS imaging around lens galaxies from the Sloan Lens Advanced Camera for Surveys (SLACS) lenses, and comparing them with the satellite population in ACS imaging of non lens galaxies selected from COSMOS, which has similar depth and resolution to the ACS images of SLACS lenses. In order to compare the samples of lens and non lens galaxies, which have intrinsically different stellar mass distributions, we measure, for the first time, the number of satellites per host as a continuous function of host stellar mass for both populations. We find that the number of satellites as a function of host stellar mass, as well as the spatial distribution are consistent between the samples. Using these results, we predict the number of satellites we would expect to find around a subset of the Cosmic Lens All Sky Survey (CLASS) lenses, and find a result consistent with the the number observed by Jackson et al. 2010. Thus we conclude that within our measurement uncertainties there is no significant difference in the satellite populations of lens and non lens galaxies.

The standard cosmological model based on cold dark matter (CDM) predicts a large number of subhalos for each galaxy-size halo. It is well known that matching the subhalos to the observed properties of luminous satellites of galaxies in the local universe poses a significant challenge to our understanding of the astrophysics of galaxy formation. We show that the cosmic evolution and host mass dependence of the luminosity function of satellites provides a powerful new diagnostic to disentangle astrophysical effects from variations in the underlying dark matter mass function. We illustrate this by comparing the results of recent observations of satellites out to $z=0.8$ based on Hubble Space Telescope images with the predictions of three different sets of state-of-the art semi-analytic models with underlying CDM power spectra and one semi-analytic model with an underlying Warm Dark Matter (WDM) power spectrum. We find that even though CDM models provide a reasonable fit to the local luminosity function of satellites around galaxies comparable or slightly larger than the Milky Way, they do not reproduce the data as well for different redshift and host galaxy stellar mass. This tension indicates that further improvements are likely to be needed in the description of star formation if the models are to be reconciled with the data. The WDM model matches the observed mass dependence and redshift evolution of satellite galaxies more closely than any of the CDM models, indicating that a modification of the underlying power spectrum may offer an alternative solution to this tension. We conclude by presenting predictions for the color magnitude relation of satellite galaxies to demonstrate how future observations will be able to further distinguish between these models and help constrain baryonic and non-baryonic physics.

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region.

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

We infer the normalization and the radial and angular distributions of the number density of satellites of massive galaxies ($\log_{10}[M_{h}^*/M\odot]>10.5$) between redshifts 0.1 and 0.8 as a function of host stellar mass, redshift, morphology and satellite luminosity. Exploiting the depth and resolution of the COSMOS HST images, we detect satellites up to eight magnitudes fainter than the host galaxies and as close as 0.3 (1.4) arcseconds (kpc). Describing the number density profile of satellite galaxies to be a projected power law such that $P(R)\propto R^{\rpower}$, we find $\rpower=-1.1\pm 0.3$. We find no dependency of $\rpower$ on host stellar mass, redshift, morphology or satellite luminosity. Satellites of early-type hosts have angular distributions that are more flattened than the host light profile and are aligned with its major axis. No significant average alignment is detected for satellites of late-type hosts. The number of satellites within a fixed magnitude contrast from a host galaxy is dependent on its stellar mass, with more massive galaxies hosting significantly more satellites. Furthermore, high-mass late-type hosts have significantly fewer satellites than early-type galaxies of the same stellar mass, likely a result of environmental differences. No significant evolution in the number of satellites per host is detected. The cumulative luminosity function of satellites is qualitatively in good agreement with that predicted using subhalo abundance matching techniques. However, there are significant residual discrepancies in the absolute normalization, suggesting that properties other than the host galaxy luminosity or stellar mass determine the number of satellites.

Some active galactic nuclei (AGN) are surrounded by extended emission-line regions (EELRs), which trace both the illumination pattern of escaping radiation and its history over the light-travel time from the AGN to the gas. From a new set of such EELRs, we present evidence that the AGN in many Seyfert galaxies undergo luminous episodes 20,000-200,000 years in duration. Motivated by the discovery of the spectacular nebula known as Hanny's Voorwerp, ionized by a powerful AGN which has apparently faded dramatically within ~ 100,000 years, Galaxy Zoo volunteers have carried out both targeted and serendipitous searches for similar emission-line clouds around low-redshift galaxies.We present the resulting list of candidates and describe spectroscopy identifying 19 galaxies with AGN-ionized regions at projected radii > 10 kpc. This search recovered known EELRs and identified additional previously unknown cases, one with detected emission to r = 37 kpc. At least 14/19 are in interacting or merging systems; tidal tails are a prime source of extraplanar ionized gas. We see a mix of one- and two-sided structures, with observed cone angles from 23-112 degrees. We consider the energy balance in the ionized clouds, with lower and upper bounds on ionizing luminosity from recombination and ionization-parameter arguments, and estimate the luminosity of the core from the far-infrared data. The implied ratio of ionizing radiation seen by the clouds to that emitted by the nucleus, for a constant nuclear source, ranges from 0.02 to > 12; 7/19 exceed unity. Small values imply heavily obscured AGN. However, large values may require that the AGN has faded over tens of thousands of years, giving us several examples of systems in which such dramatic long-period variation has occurred; this is the only current technique for addressing these timescales in AGN history. (Abridged)

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler mission. In order to determine the black hole mass in Zw 229-015 from H-beta reverberation mapping, we have carried out nightly observations with the Kast Spectrograph at the Lick 3m telescope during the dark runs from June through December 2010, obtaining 54 spectroscopic observations in total. We have also obtained nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 0.9m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect strong variability in the source, which exhibited more than a factor of 2 change in broad H-beta flux. From cross-correlation measurements, we find that the H-beta light curve has a rest-frame lag of 3.86(+0.69,-0.90) days with respect to the V-band continuum variations. We also measure reverberation lags for H-alpha and H-gamma and find an upper limit to the H-delta lag. Combining the H-beta lag measurement with a broad H-beta width of sigma = 1590+/-47 km/s measured from the root-mean-square variability spectrum, we obtain a virial estimate of M_BH = 1.00(-0.24,+0.19)*10^7 solar masses for the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability characteristics in active galactic nuclei.

We study the spatial distribution of faint satellites of intermediate redshift (0.1<z<0.8), early-type galaxies, selected from the GOODS fields. We combine high resolution HST images and state-of-the-art host subtraction techniques to detect satellites of unprecedented faintness and proximity to intermediate redshift host galaxies (up to 5.5 magnitudes fainter and as close as 0."5/2.5 kpc to the host centers). We model the spatial distribution of objects near the hosts as a combination of an isotropic, homogenous background/foreground population and a satellite population with a power law radial profile and an elliptical angular distribution. We detect a significant population of satellites, Ns =1.7 (+0.9,-0.8) that is comparable to the number of Milky Way satellites with similar host-satellite contrast.The average projected radial profile of the satellite distribution is isothermal, gamma_p= -1.0(+0.3,-0.4), which is consistent with the observed central mass density profile of massive early-type galaxies. Furthermore, the satellite distribution is highly anisotropic (isotropy is ruled out at a >99.99% confidence level). Defining phi to be the offset between the major axis of the satellite spatial distribution and the major axis of the host light profile, we find a maximum posterior probability of phi = 0 and |phi| less than 42 degrees at the 68% confidence level. The alignment of the satellite distribution with the light of the host is consistent with simulations, assuming that light traces mass for the host galaxy as observed for lens galaxies. The anisotropy of the satellite population enhances its ability to produce the flux ratio anomalies observed in gravitationally lensed quasars.