This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We present an analysis of 12 optically selected dual AGN candidates at $z < 0.34$. Each candidate was originally identified via double-peaked [O III] $\lambda$5007 emission lines, and have received follow-up $Chandra$ and $HST$ observations. Because the X-ray data are low-count ($<100$ counts) with small separations ($<1$"), a robust analysis is necessary for classifying each source. Pairing long-slit [O III] observations with existing $Chandra$ observations, we re-analyze the X-ray observations with ${\tt BAYMAX}$ to determine whether the X-ray emission from each system is more likely a single or dual point source. We find that 4 of the 12 sources are likely dual X-ray point source systems. We examine each point source's spectra via a Monte Carlo method that probabilistically identifies the likely origin of each photon. When doing so, we find that (i) the secondary X-ray point sources in 2 of the systems have $L_{\mathrm{X}}<10^{40}$ erg s$^{-1}$, such that we cannot rule out a non-AGN origin, (ii) one source has a secondary with $L_{\mathrm{X}}>10^{40}$ erg s$^{-1}$ but a spectrum that is too soft to definitively preclude being X-ray emitting diffuse gas that was photoionized by the primary AGN, and (iii) one system (SDSS J1126+2944) is a dual AGN. Additionally, using complementary $HST$ observations, we analyze a sub-sample of systems that are visually identified as merging. Our results suggest that dual AGNs may preferentially reside in mergers with small separations, consistent with both simulations and observations.

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We present the Advanced Camera for Surveys Active Galactic Nuclei (ACS-AGN) Catalog, a catalog of 2585 active galactic nucleus (AGN) host galaxies that are at redshifts 0.2<z<2.5 and that were imaged with the Hubble Space Telescope's Advanced Camera for Surveys (ACS). Using the ACS General Catalog (ACS-GC) as our initial sample, we select an AGN subsample using Spitzer and Chandra data along with their respective established AGN selection criteria. We then gather further multi-wavelength photometric data in order to construct spectral energy distributions (SEDs). Using these SEDs we are able to derive multiple AGN and host galaxy properties, such as star formation rate, AGN luminosity, stellar mass, and nuclear column density. From these data, we show that AGN host galaxies tend to lie below the star-forming main sequence, with X-ray-selected AGN host galaxies being more offset than IR-selected AGN host galaxies. This suggests that there is some process, possibly negative feedback, in AGN host galaxies causing decreased star formation. We also demonstrate that there is a positive trend between star formation rate and AGN luminosity in AGN host galaxies, in individual redshift bins and across all redshift bins, and that both are correlated with the stellar mass of their galaxies. This points towards an underlying link between the stellar mass, stellar growth, and SMBH growth in a galaxy.

Dynamical black hole mass measurements in some gas-rich galaxy mergers indicate that they are overmassive relative to their host galaxy properties. Overmassive black holes in these systems present a conflict with the standard progression of galaxy merger - quasar evolution; an alternative explanation is that a nuclear concentration of molecular gas driven inward by the merger is affecting these dynamical black hole mass estimates. We test for the presence of such gas near the two black holes in NGC 6240 using long-baseline ALMA Band 6 observations (beam size 0"06 $\times$ 0"03 or 30 pc$\times$15 pc). We find (4.2-9.8) $\times10^{7}$ M$_{\odot}$ and (1.2-7.7) $\times10^{8}$ M$_{\odot}$ of molecular gas within the resolution limit of the original black hole mass measurements for the north and south black holes, respectively. In the south nucleus, this measurement implies that 6-89% of the original black hole mass measurement actually comes from molecular gas, resolving the tension in the original black hole scaling relations. For the north, only 5% to 11% is coming from molecular gas, suggesting the north black hole is actually overmassive. Our analysis provides the first measurement of significant molecular gas masses contaminating dynamical black hole mass measurements. These high central molecular gas densities further present a challenge to theoretical black hole accretion prescriptions, which often assume accretion proceeds rapidly through the central 10 pc.

Hyper-luminous X-ray sources (HLXs; L_X>10^41 erg s^-1) are off-nuclear X-ray sources in galaxies and strong candidates for intermediate-mass black holes (IMBHs). We have constructed a sample of 169 HLX candidates by combining X-ray detections from the Chandra Source Catalog (Version 2) with galaxies from the Sloan Digital Sky Survey and registering individual images for improved relative astrometric accuracy. The spatial resolution of Chandra allows for the sample to extend out to z~0.9. Optical counterparts are detected among one-fourth of the sample, one-third of which are consistent with dwarf galaxy stellar masses. The average intrinsic X-ray spectral slope indicates efficient accretion, potentially driven by galaxy mergers, and the column densities suggest one-third of the sample has significant X-ray absorption. We find that 144 of the HLX candidates have X-ray emission that is significantly in excess of the expected contribution from star formation and hot gas, strongly suggesting that they are produced by accretion onto black holes more massive than stars. After correcting for an average background or foreground contamination rate of 8%, we estimate that at least ~20 of the HLX candidates are consistent with IMBH masses, and this estimate is potentially several times higher assuming more efficient accretion. This catalog currently represents the largest sample of uniformly-selected, off-nuclear IMBH candidates. These sources may represent scenarios in which a low-mass galaxy hosting an IMBH has merged with a more massive galaxy and provide an excellent sample for testing models of low-mass BH formation and merger-driven growth.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

We present the first results from BAYMAX (Bayesian AnalYsis of Multiple AGN in X-rays), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual AGNs is to use X-ray observations, only sources that are widely separated relative to the instrument PSF are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandra's angular resolution (<1''). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ~0.3'' from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find the probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0.3'' to be very low. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, and studying the dual AGN population where past spatial resolution limits have prevented systematic analyses.

Gravitational wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. Multi-wavelength observations with next-generation telescopes, including systematic large-area surveys, will provide unprecedented opportunities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales. We highlight ways in which these observations will open a new avenue for multi-messenger studies in the dawn of low-frequency (~ nHz - mHz) GW astronomy.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black hole growth lies in a thoroughly-studied census of dual AGN across cosmic history, which will be enabled by next-generation observational capabilities, theoretical advances, and simulations. This white paper outlines the key questions in galaxy mergers, dual and offset AGN, and proposes multiwavelength solutions using future high-resolution observatories in the X-rays (AXIS, Lynx), near and mid-infrared (30 meter class telescopes, JWST), and submillimeter (ALMA).

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Galaxy pairs with separations of only a few kpc represent important stages in the merger-driven growth of supermassive black holes (SMBHs). However, such mergers are difficult to identify observationally due to the correspondingly small angular scales. In Paper I we presented a method of finding candidate kpc-scale galaxy mergers that is leveraged on the selection of X-ray sources spatially offset from the centers of host galaxies. In this paper we analyze new Hubble Space Telescope (HST) WFC3 imaging for six of these sources to search for signatures of galaxy mergers. The HST imaging reveals that four of the six systems are on-going galaxy mergers with separations of 1.2-6.6 kpc (offset AGN). The nature of the remaining two spatially offset X-ray sources is ambiguous and may be associated with super-Eddington accretion in X-ray binaries. The ability of this sample to probe small galaxy separations and minor mergers makes it uniquely suited for testing the role of galaxy mergers for AGN triggering. We find that galaxy mergers with only one AGN are predominantly minor mergers with mass ratios similar to the overall population of galaxy mergers. By comparison, galaxy mergers with two AGN are biased toward major mergers and larger nuclear gas masses. Finally, we find that the level of SMBH accretion increases toward smaller mass ratios (major mergers). This result suggests the mass ratio effects not only the frequency of AGN triggering but also the rate of SMBH growth in mergers.

Double-peaked narrow emission lines in active galactic nucleus (AGN) spectra can be produced by AGN outflows, rotation, or dual AGNs, which are AGN pairs in ongoing galaxy mergers. Consequently, double-peaked narrow AGN emission lines are useful tracers of the coevolution of galaxies and their supermassive black holes, as driven by AGN feedback and AGN fueling. We investigate this concept further with follow-up optical longslit observations of a sample of 95 Sloan Digital Sky Survey (SDSS) galaxies that have double-peaked narrow AGN emission lines. Based on a kinematic analysis of the longslit spectra, we confirm previous work that finds that the majority of double-peaked narrow AGN emission lines are associated with outflows. We also find that eight of the galaxies have companion galaxies with line-of-sight velocity separations < 500 km/s and physical separations <30 kpc. Since we find evidence of AGNs in both galaxies, all eight of these systems are compelling dual AGN candidates. Galaxies with double-peaked narrow AGN emission lines occur in such galaxy mergers at least twice as often as typical active galaxies. Finally, we conclude that at least 3% of SDSS galaxies with double-peaked narrow AGN emission lines are found in galaxy mergers where both galaxies are resolved in SDSS imaging.

Gravitational-wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. The ngVLA offers unparalleled capabilities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales, opening a new avenue for multi-messenger studies in the dawn of low-frequency GW astronomy.

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Theoretical models and numerical simulations have established a framework of galaxy evolution in which galaxies merge and create dual supermassive black holes (with separations of one to ten kiloparsecs), which eventually sink into the centre of the merger remnant, emit gravitational waves and coalesce. The merger also triggers star formation and supermassive black hole growth, and gas outflows regulate the stellar content. Although this theoretical picture is supported by recent observations of starburst-driven and supermassive black hole-driven outflows, it remains unclear how these outflows interact with the interstellar medium. Furthermore, the relative contributions of star formation and black hole activity to galactic feedback remain unknown. Here we report observations of dual outflows in the central region of the prototypical merger NGC 6240. We find a black-hole-driven outflow of [O III] to the northeast and a starburst-driven outflow of H{\alpha} to the northwest. The orientations and positions of the outflows allow us to isolate them spatially and study their properties independently. We estimate mass outflow rates of 10 and 75 solar masses per year for the H{\alpha} bubble and the [O III] cone, respectively. Their combined mass outflow is comparable to the star formation rate, suggesting that negative feedback on star formation is occurring.

We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.

Galaxy mergers are capable of triggering both star formation and active galactic nuclei (AGN) and therefore may represent an important pathway in the co-evolution of galaxies and supermassive black holes (SMBHs). However, correlated enhancements of merger-induced star formation and AGN triggering may be hidden by the variable conditions and timescales during which they occur. In Paper I, we presented evidence of merger-triggered AGN in a sample of six late-stage galaxy mergers (2-8 kpc nuclear separations). In this follow-up work, we use multi-wavelength Hubble Space Telescope imaging and additional archival data to examine their star-forming properties to test for merger-triggered star formation, and if it is correlated with SMBH growth. We find that the morphological asymmetries are correlated with enhanced specific star formation rates, indicating the presence of merger-triggered star formation. Additionally, the stellar populations become younger with increasing radius from the nucleus, indicating that the merger-induced star formation primarily occurs on global scales. However, we also find that the star formation rate enhancements are consistent with or lower than those of larger separation galaxy pair samples. This result is consistent with simulations predicting a decline of the global star formation rates in late-stage galaxy mergers with <10 kpc nuclear separations. Finally, we find that enhancements in specific star formation rate and AGN luminosity are positively correlated, but that an average temporal delay of >10^8 yrs likely exists between the peak of global star formation and the onset of AGN triggering in 80% of the systems.

We analyze the evolution of the dual AGN in ESO 509-IG066, a galaxy pair located at z=0.034 whose nuclei are separated by 11 kpc. Previous observations with XMM-Newton on this dual AGN found evidence for two moderately obscured ($N_H\sim10^{22}$ cm$^{-2}$) X-ray luminous ($L_X\sim10^{43}$ erg/s) nuclear sources. We present an analysis of subsequent Chandra, NuSTAR and Swift/XRT observations that show one source has dropped in flux by a factor of 10 between 2004 and 2011, which could be explained by either an increase in the absorbing column or an intrinsic fading of the central engine possibly due to a decrease in mass accretion. Both of these scenarios are predicted by galaxy merger simulations. The source which has dropped in flux is not detected by NuSTAR, which argues against absorption, unless it is extreme. However, new Keck/LRIS optical spectroscopy reveals a previously unreported broad H-alpha line which is highly unlikely to be visible under the extreme absorption scenario. We therefore conclude that the black hole in this nucleus has undergone a dramatic drop in accretion rate. From AO-assisted near-infrared integral-field spectroscopy of the other nucleus, we find evidence that the galaxy merger is having a direct effect on the kinematics of the gas close to the nucleus of the galaxy, providing a direct observational link between the galaxy merger and the mass accretion rate on to the black hole.

We present the discovery of an active galactic nucleus (AGN) that is turning off and then on again in the z=0.06 galaxy SDSS J1354+1327. This episodic nuclear activity is the result of discrete accretion events, which could have been triggered by a past interaction with the companion galaxy that is currently located 12.5 kpc away. We originally targeted SDSS J1354+1327 because its Sloan Digital Sky Survey spectrum has narrow AGN emission lines that exhibit a velocity offset of 69 km s$^{-1}$ relative to systemic. To determine the nature of the galaxy and its velocity-offset emission lines, we observed SDSS J1354+1327 with Chandra/ACIS, Hubble Space Telescope/Wide Field Camera 3, Apache Point Observatory optical longslit spectroscopy, and Keck/OSIRIS integral-field spectroscopy. We find a ~10 kpc cone of photoionized gas south of the galaxy center and a ~1 kpc semi-spherical front of shocked gas, which is responsible for the velocity offset in the emission lines, north of the galaxy center. We interpret these two outflows as the result of two separate AGN accretion events; the first AGN outburst created the southern outflow, and then $<10^5$ yrs later the second AGN outburst launched the northern shock front. SDSS J1354+1327 is the galaxy with the strongest evidence for an AGN that has turned off and then on again, and it fits into the broader context of AGN flickering that includes observations of AGN light echoes.

We apply an analytic Markov Chain Monte Carlo model to a sample of 18 AGN-driven biconical outflows that we identified from a sample of active galaxies with double-peaked narrow emission lines at z < 0.1 in the Sloan Digital Sky Survey. We find that 8/18 are best described as asymmetric bicones, 8/18 are nested bicones, and 2/18 are symmetric bicones. From the geometry and kinematics of the models, we find that these moderate-luminosity AGN outflows are large and energetic. The biconical outflows axes are randomly oriented with respect to the photometric major axis of the galaxy, implying a randomly oriented and clumpier torus to collimate the outflow, but the torus also allows some radiation to escape equatorially. We find that 16/18 (89%) outflows are energetic enough to drive a two-staged feedback process in their host galaxies. All of these outflows geometrically intersect the photometric major axis of the galaxy, and 23% of outflow host galaxies are significantly redder or have significantly lower specific star formation rates when compared to a matched sample of active galaxies.

While 2% of active galactic nuclei (AGNs) exhibit narrow emission lines with line-of-sight velocities that are significantly offset from the velocity of the host galaxy's stars, the nature of these velocity offsets is unknown. We investigate this question with Chandra/ACIS and Hubble Space Telescope/Wide Field Camera 3 observations of seven velocity-offset AGNs at z<0.12, and all seven galaxies have a central AGN but a peak in emission that is spatially offset by < kpc from the host galaxy's stellar centroid. These spatial offsets are responsible for the observed velocity offsets and are due to shocks, either from AGN outflows (in four galaxies) or gas inflowing along a bar (in three galaxies). We compare our results to a velocity-offset AGN whose velocity offset originates from a spatially offset AGN in a galaxy merger. The optical line flux ratios of the offset AGN are consistent with pure photoionization, while the optical line flux ratios of our sample are consistent with contributions from photoionization and shocks. We conclude that these optical line flux ratios could be efficient for separating velocity-offset AGNs into subgroups of offset AGNs -- which are important for studies of AGN fueling in galaxy mergers -- and central AGNs with shocks -- where the outflows are biased towards the most energetic outflows that are the strongest drivers of feedback.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

Galaxy mergers are likely to play a role in triggering active galactic nuclei (AGN), but the conditions under which this process occurs are poorly understood. In Paper I, we constructed a sample of spatially offset X-ray AGN that represent galaxy mergers hosting a single AGN. In this paper, we use our offset AGN sample to constrain the parameters that affect AGN observability in galaxy mergers. We also construct dual AGN samples with similar selection properties for comparison. We find that the offset AGN fraction shows no evidence for a dependence on AGN luminosity, while the dual AGN fractions show stronger evidence for a positive dependence, suggesting that the merger events forming dual AGN are more efficient at instigating accretion onto supermassive black holes than those forming offset AGN. We also find that the offset and dual AGN fractions both have a negative dependence on nuclear separation and are similar in value at small physical scales. This dependence may become stronger when restricted to high AGN luminosities, though a larger sample is needed for confirmation. These results indicate that the probability of AGN triggering increases at later merger stages. This study is the first to systematically probe down to nuclear separations of <1 kpc (~0.8 kpc) and is consistent with predictions from simulations that AGN observability peaks in this regime. We also find that the offset AGN are not preferentially obscured compared to the parent AGN sample, suggesting that our selection may be targeting galaxy mergers with relatively dust-free nuclear regions.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We present optical longslit observations of the complete sample of 71 Type 2 active galactic nuclei (AGNs) with double-peaked narrow emission lines at $z < 0.1$ in the Sloan Digital Sky Survey. Double-peaked emission lines are produced by a variety of mechanisms including disk rotation, kpc-scale dual AGNs, and NLR kinematics (outflows or inflows). We develop a novel kinematic classification technique to determine the nature of these objects using longslit spectroscopy alone. We determine that 86% of the double-peaked profiles are produced by moderate luminosity AGN outflows, 6% are produced by rotation, and 8% are ambiguous. While we are unable to directly identify dual AGNs with longslit data alone, we explore their potential kinematic classifications with this method. We also find a positive correlation between the narrow-line region (NLR) size and luminosity of the AGN NLRs (R$_{\mathrm{NLR}}\propto \; {\mathrm{L}_{\mathrm{[OIII]}}}^{0.21 \pm 0.05}$), indicating a clumpy two-zone ionization model for the NLR.

We obtained Keck/OSIRIS near-IR adaptive optics-assisted integral-field spectroscopy to probe the morphology and kinematics of the ionized gas in four velocity-offset active galactic nuclei (AGNs) from the Sloan Digital Sky Survey. These objects possess optical emission lines that are offset in velocity from systemic as measured from stellar absorption features. At a resolution of ~0.18", OSIRIS allows us to distinguish which velocity offset emission lines are produced by the motion of an AGN in a dual supermassive black hole system, and which are produced by outflows or other kinematic structures. In three galaxies, J1018+2941, J1055+1520 and J1346+5228, the spectral offset of the emission lines is caused by AGN-driven outflows. In the remaining galaxy, J1117+6140, a counterrotating nuclear disk is observed that contains the peak of Pa$\alpha$ emission 0.2" from the center of the galaxy. The most plausible explanation for the origin of this spatially and kinematically offset peak is that it is a region of enhanced Pa$\alpha$ emission located at the intersection zone between the nuclear disk and the bar of the galaxy. In all four objects, the peak of ionized gas emission is not spatially coincident with the center of the galaxy as traced by the peak of the near-IR continuum emission. The peaks of ionized gas emission are spatially offset from the galaxy centers by 0.1"-0.4" (0.1-0.7 kpc). We find that the velocity offset originates at the location of this peak of emission, and the value of the offset can be directly measured in the velocity maps. The emission-line ratios of these four velocity-offset AGNs can be reproduced only with a mixture of shocks and AGN photoionization. Shocks provide a natural explanation for the origin of the spatially and spectrally offset peaks of ionized gas emission in these galaxies.

We present a sample of 18 optically-selected and X-ray detected spatially offset active galactic nuclei (AGN) from the Sloan Digital Sky Survey (SDSS). In 9 systems, the X-ray AGN is spatially offset from the galactic stellar core that is located within the 3'' diameter SDSS spectroscopic fiber. In 11 systems, the X-ray AGN is spatially offset from a stellar core that is located outside the fiber, with an overlap of 2. To build the sample, we cross-matched Type II AGN selected from the SDSS galaxy catalogue with archival Chandra imaging and employed our custom astrometric and registration procedure. The projected angular (physical) offsets span a range of 0."6 (0.8 kpc) to 17."4 (19.4 kpc), with a median value of 2."7 (4.6 kpc). The offset nature of an AGN is an unambiguous signature of a galaxy merger, and these systems can be used to study the properties of AGN in galaxy mergers without the biases introduced by morphological merger selection techniques. In this paper (Paper I), we use our sample to assess the kinematics of AGN photoionized gas in galaxy mergers. We find that spectroscopic offset AGN selection may be up to 89% incomplete due to small projected velocity offsets. We also find that the magnitude of the velocity offsets are generally larger than expected if our spatial selection introduces a bias toward face-on orbits, suggesting the presence of complex kinematics in the emission line gas of AGN in galaxy mergers.

In the last few years, it became possible to observationally resolve galaxies with two distinct nuclei in their centre. For separations smaller than 10kpc, dual and offset active galactic nuclei (AGN) are distinguished: in dual AGN, both nuclei are active, whereas in offset AGN only one nucleus is active. To study the origin of such AGN pairs, we employ a cosmological, hydrodynamic simulation with a large volume of (182Mpc)^3 from the set of Magneticum Pathfinder Simulations. The simulation self-consistently produces 35 resolved black hole (BH) pairs at redshift z=2, with a comoving distance smaller than 10kpc. 14 of them are offset AGN and nine are dual AGN, resulting in a fraction of (1.2 \pm 0.3)% AGN pairs with respect to the total number of AGN. In this paper, we discuss fundamental differences between the BH and galaxy properties of dual AGN, offset AGN and inactive BH pairs and investigate their different triggering mechanisms. We find that in dual AGN the BHs have similar masses and the corresponding BH from the less massive progenitor galaxy always accretes with a higher Eddington ratio. In contrast, in offset AGN the active BH is typically more massive than its non-active counterpart. Furthermore, dual AGN in general accrete more gas from the intergalactic medium than offset AGN and non-active BH pairs. This highlights that merger events, particularly minor mergers, do not necessarily lead to strong gas inflows and thus, do not always drive strong nuclear activity.

We have examined a subset of 18 active galactic nuclei (AGNs) drawn from a sample of 81 galaxies that possess double-peaked narrow optical emission line spectra in the Sloan Digital Sky Survey, have two optical AGN emission components separated by >0.2", and are detected in the Faint Images of the Radio Sky at Twenty-centimeters survey. Without follow-up observations, the sources of the double-peaked narrow emission lines are uncertain, and may be produced by kpc-scale separation dual active supermassive black holes, AGN outflows, or disk rotation. In this work, we propose a new methodology to characterize double-peaked narrow emission-line galaxies based on optical long-slit spectroscopy and high resolution multi-band Very Large Array observations. The nature of the radio emission in the sample galaxies is varied. Of the 18 galaxies, we detect two compact flat-spectrum radio cores with projected spatial separations on the sky between 0.6-1.6 kpc in three galaxies: J1023+3243, J1158+3231, and J1623+0808. The two radio sources are spatially coincident with the two optical components of ionized gas with AGN-like line ratios, which confirms the presence of dual AGNs in these three galaxies. Dual AGNs account for only ~15% (3/18) of the double-peaked AGNs in our sample. Gas kinematics produce ~75% (13/18) of the double-peaked narrow emission lines, distributed in the following way: 7 AGN wind-driven outflows, 5 radio-jet driven outflows, and one rotating narrow-line region. The remaining ~10% (2/18) are ambiguous cases. Our method demonstrates the power of spatially resolved spectroscopy and high resolution radio observations for the identification of AGN outflows and AGN pairs with angular separations as small as 0.18".

Dual active galactic nuclei (AGNs) and offset AGNs are kpc-scale separation supermassive black holes pairs created during galaxy mergers, where both or one of the black holes are AGNs, respectively. These dual and offset AGNs are valuable probes of the link between mergers and AGNs but are challenging to identify. Here we present Chandra/ACIS observations of 12 optically-selected dual AGN candidates at z < 0.34, where we use the X-rays to identify AGNs. We also present HST/WFC3 observations of 10 of these candidates, which reveal any stellar bulges accompanying the AGNs. We discover a dual AGN system with separation of 2.2 kpc, where the two stellar bulges have coincident [O III] and X-ray sources. This system is an extremely minor merger (460:1) that may include a dwarf galaxy hosting an intermediate mass black hole. We also find six single AGNs, and five systems that are either dual or offset AGNs with separations < 10 kpc. Four of the six dual AGNs and dual/offset AGNs are in ongoing major mergers, and these AGNs are 10 times more luminous, on average, than the single AGNs in our sample. This hints that major mergers may preferentially trigger higher luminosity AGNs. Further, we find that confirmed dual AGNs have hard X-ray luminosities that are half of those of single AGNs at fixed [O III] luminosity, on average. This could be explained by high densities of gas funneled to galaxy centers during mergers, and emphasizes the need for deeper X-ray observations of dual AGN candidates.

Offset active galactic nuclei (AGNs) are AGNs that are in ongoing galaxy mergers, which produce kinematic offsets in the AGNs relative to their host galaxies. Offset AGNs are also close relatives of dual AGNs. We conduct a systematic search for offset AGNs in the Sloan Digital Sky Survey, by selecting AGN emission lines that exhibit statistically significant line-of-sight velocity offsets relative to systemic. From a parent sample of 18314 Type 2 AGNs at z<0.21, we identify 351 offset AGN candidates with velocity offsets of 50 km/s < |v| < 410 km/s. When we account for projection effects in the observed velocities, we estimate that 4% - 8% of AGNs are offset AGNs. We designed our selection criteria to bypass velocity offsets produced by rotating gas disks, AGN outflows, and gravitational recoil of supermassive black holes, but follow-up observations are still required to confirm our candidates as offset AGNs. We find that the fraction of AGNs that are offset candidates increases with AGN bolometric luminosity, from 0.7% to 6% over the luminosity range 43 < log(L_bol) [erg/s] < 46. If these candidates are shown to be bona fide offset AGNs, then this would be direct observational evidence that galaxy mergers preferentially trigger high-luminosity AGNs. Finally, we find that the fraction of AGNs that are offset AGN candidates increases from 1.9% at z=0.1 to 32% at z=0.7, in step with the growth in the galaxy merger fraction over the same redshift range.

We present observations of extended, 20-kpc scale soft X-ray gas around a luminous obscured quasar hosted by an ultra-luminous infrared galaxy caught in the midst of a major merger. The extended X-ray emission is well fit as a thermal gas with a temperature of kT ~ 280 eV and a luminosity of L_X ~ 10^42 erg/s and is spatially coincident with a known ionized gas outflow. Based on the X-ray luminosity, a factor of ~10 fainter than the [OIII] emission, we conclude that the X-ray emission is either dominated by photoionization, or by shocked emission from cloud surfaces in a hot quasar-driven wind.

This whitepaper describes how the VLASS could be designed in a manner to allow the identification of candidate dual active galactic nuclei (AGN) at separations <7 kpc. Dual AGN represent a clear marker of two supermassive black holes within an ongoing merger. A dual AGN survey will provide a wealth of studies in structure growth and gravitational-wave science. Radio wavelengths are ideal for identifying close pairs, as disturbed stellar and gaseous material can obscure their presence in optical and shorter wavelengths. With sufficiently high resolution and sensitivity, a large-scale radio imaging survey like the VLASS will uncover many of these systems and provide the means to broadly study the radio properties of candidate dual systems revealed at other wavelengths. We determine that the ideal survey for our purposes will be at as high a resolution as possible, with significantly more science return in A array at L-band or higher, or B array at C-band or higher. We describe a range of potential survey parameters within this document. Based on the analysis outlined in this whitepaper, our ideal survey would create a catalogue of $\gtrsim$100 dual AGN in either: 1) a medium-sensitivity (~1 mJy detection threshold), wide-field (few thousand square degree) survey, or 2) a high-sensitivity (~10 uJy threshold) survey of several hundred square degrees.

COSMOS J100043.15+020637.2 is a merger remnant at z = 0.36 with two optical nuclei, NW and SE, offset by 500 mas (2.5 kpc). Prior studies suggest two competing scenarios for these nuclei: (1) SE is an active galactic nucleus (AGN) lost from NW due to a gravitational-wave recoil. (2) NW and SE each contain an AGN, signaling a gravitational-slingshot recoil or inspiralling AGNs. We present new images from the Very Large Array (VLA) at a frequency nu = 9.0 GHz and a FWHM resolution theta = 320 mas (1.6 kpc), and the Very Long Baseline Array (VLBA) at nu = 1.52 GHz and theta = 15 mas (75 pc). The VLA imaging is sensitive to emission driven by AGNs and/or star formation, while the VLBA imaging is sensitive only to AGN-driven emission. No radio emission is detected at these frequencies. Folding in prior results, we find: (a) The properties of SE and its adjacent X-ray feature resemble those of the unobscured AGN in NGC 4151, albeit with a much higher narrow emission-line luminosity. (b) The properties of NW are consistent with it hosting a Compton-thick AGN that warms ambient dust, photoionizes narrow emission-line gas and is free-free absorbed by that gas. Finding (a) is consistent with scenarios (1) and (2). Finding (b) weakens the case for scenario (1) and strengthens the case for scenario (2). Follow-up observations are suggested.

Dual supermassive black holes (SMBHs) with kiloparsec scale separations in merger-remnant galaxies are informative tracers of galaxy evolution, but the avenue for identifying them in large numbers for such studies is not yet clear. One promising approach is to target spectroscopic signatures of systems where both SMBHs are fueled as dual active galactic nuclei (AGNs), or where one SMBH is fueled as an offset AGN. Dual AGNs may produce double-peaked narrow AGN emission lines, while offset AGNs may produce single-peaked narrow AGN emission lines with line-of-sight velocity offsets relative to the host galaxy. We search for such dual and offset systems among 173 Type 2 AGNs at z<0.37 in the AGN and Galaxy Evolution Survey (AGES), and we find two double-peaked AGNs and five offset AGN candidates. When we compare these results to a similar search of the DEEP2 Galaxy Redshift Survey and match the two samples in color, absolute magnitude, and minimum velocity offset, we find that the fraction of AGNs that are dual SMBH candidates increases from z=0.25 to z=0.7 by a factor of ~6 (from 2/70 to 16/91, or 2.9% to 18%). This may be associated with the rise in the galaxy merger fraction over the same cosmic time. As further evidence for a link with galaxy mergers, the AGES offset and dual AGN candidates are tentatively ~3 times more likely than the overall AGN population to reside in a host galaxy that has a companion galaxy (from 16/173 to 2/7, or 9% to 29%). Follow-up observations of the seven offset and dual AGN candidates in AGES will definitively distinguish velocity offsets produced by dual SMBHs from those produced by narrow-line region kinematics, and will help sharpen our observational approach to detecting dual SMBHs.

We study the radial dependence in stellar populations of 33 nearby early-type galaxies with central stellar velocity dispersions sigma* > 150 km/s. We measure stellar population properties in composite spectra, and use ratios of these composites to highlight the largest spectral changes as a function of radius. Based on stellar population modeling, the typical star at 2 R_e is old (~10 Gyr), relatively metal poor ([Fe/H] -0.5), and alpha-enhanced ([Mg/Fe]~0.3). The stars were made rapidly at z~1.5-2 in shallow potential wells. Declining radial gradients in [C/Fe], which follow [Fe/H], also arise from rapid star formation timescales due to declining carbon yields from low-metallicity massive stars. In contrast, [N/Fe] remains high at large radius. Stars at large radius have different abundance ratio patterns from stars in the center of any present-day galaxy, but are similar to Milky Way thick disk stars. Our observations are thus consistent with a picture in which the stellar outskirts are built up through minor mergers with disky galaxies whose star formation is truncated early (z~1.5-2).

We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.

We present a sample of 131 quasars from the Sloan Digital Sky Survey at redshifts 0.8<z<1.6 with double peaks in either of the high-ionization narrow emission lines [NeV]3426 or [NeIII]3869. These sources were selected with the intention of identifying high-redshift analogs of the z<0.8 active galactic nuclei (AGN) with double-peaked [OIII]5007 lines, which might represent AGN outflows or dual AGN. Lines of high-ionization potential are believed to originate in the inner, highly photoionized portion of the narrow line region (NLR), and we exploit this assumption to investigate the possible kinematic origins of the double-peaked lines. For comparison, we measure the [NeV]3426 and [NeIII]3869 double peaks in low-redshift (z<0.8) [OIII]-selected sources. We find that [NeV]3426 and [NeIII]3869 show a correlation between line-splitting and line-width similar to that of [OIII]5007 in other studies; and the velocity-splittings are correlated with the quasar Eddington ratio. These results suggest an outflow origin for at least a subset of the double-peaks, allowing us to study the high-ionization gas kinematics around quasars. However, we find that a non-neligible fraction of our sample show no evidence for an ionization stratification. For these sources, the outflow scenario is less compelling, leaving the dual AGN scenario as a viable possibility. Finally, we find that our sample shows an anti-correlation between the velocity-offset ratio and luminosity ratio of the components, which is a potential dynamical argument for the presence of dual AGN. Therefore, this study serves as a first attempt at extending the selection of candidate dual AGN to higher redshifts.

We present PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.

We have used the Caltech Submillimeter Observatory (CSO) to follow-up a sample of WISE-selected, hyperluminous galaxies, so called W1W2-dropout galaxies. This is a rare (~ 1000 all-sky) population of galaxies at high redshift (peaks at z=2-3), that are faint or undetected by WISE at 3.4 and 4.6 um, yet are clearly detected at 12 and 22 um. The optical spectra of most of these galaxies show significant AGN activity. We observed 14 high-redshift (z > 1.7) W1W2-dropout galaxies with SHARC-II at 350 to 850 um, with 9 detections; and observed 18 with Bolocam at 1.1 mm, with five detections. Warm Spitzer follow-up of 25 targets at 3.6 and 4.5 um, as well as optical spectra of 12 targets are also presented in the paper. Combining WISE data with observations from warm Spitzer and CSO, we constructed their mid-IR to millimeter spectral energy distributions (SEDs). These SEDs have a consistent shape, showing significantly higher mid-IR to submm ratios than other galaxy templates, suggesting a hotter dust temperature. We estimate their dust temperatures to be 60-120 K using a single-temperature model. Their infrared luminosities are well over 10^{13} Lsun. These SEDs are not well fitted with existing galaxy templates, suggesting they are a new population with very high luminosity and hot dust. They are likely among the most luminous galaxies in the Universe. We argue that they are extreme cases of luminous, hot dust-obscured galaxies (DOGs), possibly representing a short evolutionary phase during galaxy merging and evolution. A better understanding of their long-wavelength properties needs ALMA as well as Herschel data.

Merger-remnant galaxies with kpc-scale separation dual active galactic nuclei (AGNs) should be widespread as a consequence of galaxy mergers and triggered gas accretion onto supermassive black holes, yet very few dual AGNs have been observed. Galaxies with double-peaked narrow AGN emission lines in the Sloan Digital Sky Survey are plausible dual AGN candidates, but their double-peaked profiles could also be the result of gas kinematics or AGN-driven outflows and jets on small or large scales. To help distinguish between these scenarios, we have obtained spatial profiles of the AGN emission via follow-up long-slit spectroscopy of 81 double-peaked narrow-line AGNs in SDSS at 0.03 < z < 0.36 using Lick, Palomar, and MMT Observatories. We find that all 81 systems exhibit double AGN emission components with ~kpc projected spatial separations on the sky, which suggests that they are produced by kpc-scale dual AGNs or kpc-scale outflows, jets, or rotating gaseous disks. In addition, we find that the subsample (58%) of the objects with spatially compact emission components may be preferentially produced by dual AGNs, while the subsample (42%) with spatially extended emission components may be preferentially produced by AGN outflows. We also find that for 32% of the sample the two AGN emission components are preferentially aligned with the host galaxy major axis, as expected for dual AGNs orbiting in the host galaxy potential. Our results both narrow the list of possible physical mechanisms producing the double AGN components, and suggest several observational criteria for selecting the most promising dual AGN candidates from the full sample of double-peaked narrow-line AGNs. Using these criteria, we determine the 17 most compelling dual AGN candidates in our sample.

We present the Advanced Camera for Surveys General Catalog (ACS-GC), a photometric and morphological database using publicly available data obtained with the Advanced Camera for Surveys (ACS) instrument on the Hubble Space Telescope. The goal of the ACS-GC database is to provide a large statistical sample of galaxies with reliable structural and distance measurements to probe the evolution of galaxies over a wide range of look-back times. The ACS-GC includes approximately 470,000 astronomical sources (stars + galaxies) derived from the AEGIS, COSMOS, GEMS, and GOODS surveys. Galapagos was used to construct photometric (SExtractor) and morphological (Galfit) catalogs. The analysis assumes a single S\'ersic model for each object to derive quantitative structural parameters. We include publicly available redshifts from the DEEP2, COMBO-17, TKRS, PEARS, ACES, CFHTLS,and zCOSMOS surveys to supply redshifts (spectroscopic and photometric) for a considerable fraction (~74%) of the imaging sample. The ACS-GC includes color postage stamps, Galfit residual images, and photometry, structural parameters, and redshifts combined into a single catalog.

We use the Mitchell Spectrograph (formerly VIRUS-P) on the McDonald Observatory 2.7m Harlan J. Smith Telescope to search for the chemical signatures of massive elliptical galaxy assembly. The Mitchell Spectrograph is an integral-field spectrograph with a uniquely wide field of view (107x107 sq arcsec), allowing us to achieve remarkably high signal-to-noise ratios of ~20-70 per pixel in radial bins of 2-2.5 times the effective radii of the eight galaxies in our sample. Focusing on a sample of massive elliptical galaxies with stellar velocity dispersions sigma* > 150 km/s, we study the radial dependence in the equivalent widths (EWs) of key metal absorption lines. By twice the effective radius, the Mgb EWs have dropped by ~50%, and only a weak correlation between sigma* and Mgb EW remains. The Mgb EWs at large radii are comparable to those seen in the centers of elliptical galaxies that are approximately an order of magnitude less massive. We find that the well-known metallicity gradients often observed within an effective radius continue smoothly to 2.5R_e, while the abundance ratio gradients remain flat. Much like the halo of the Milky Way, the stellar halos of our galaxies have low metallicities and high alpha-abundance ratios, as expected for very old stars formed in small stellar systems. Our observations support a picture in which the outer parts of massive elliptical galaxies are built by the accretion of much smaller systems whose star formation history was truncated at early times.

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor {\alpha}CO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log({\alpha}CO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

The X-ray source CXOXBJ142607.6+353351 (CXOJ1426+35), which was identified in a 172 ks Chandra image in the Bootes field, shows double-peaked rest-frame optical/UV emission lines, separated by 0.69" (5.5 kpc) in the spatial dimension and by 690 km s^-1 in the velocity dimension. The high excitation lines and emission line ratios indicate both systems are ionized by an AGN continuum, and the double-peaked profile resembles that of candidate dual AGN. At a redshift of z=1.175, this source is the highest redshift candidate dual AGN yet identified. However, many sources have similar emission line profiles for which other interpretations are favored. We have analyzed the substantial archival data available in this field, as well as acquired near-infrared (NIR) adaptive optics (AO) imaging and NIR slit spectroscopy. The X-ray spectrum is hard, implying a column density of several 10^23 cm^-2. Though heavily obscured, the source is also one of the brightest in the field, with an absorption-corrected 2-10 keV luminosity of ~10^45 erg s^-1. Outflows driven by an accretion disk may produce the double-peaked lines if the central engine accretes near the Eddington limit. However, we may be seeing the narrow line regions of two AGN following a galactic merger. While the AO image reveals only a single source, a second AGN would easily be obscured by the significant extinction inferred from the X-ray data. Understanding the physical processes producing the complex emission line profiles seen in CXOJ1426+35 and related sources is important for interpreting the growing population of dual AGN candidates.

Type IIn SNe show spectral evidence for strong interaction between their blast wave and dense circumstellar material (CSM) around the progenitor star. SN2010jl was the brightest core-collapse SN in 2010, and it was a Type IIn explosion with strong CSM interaction. Andrews et al. recently reported evidence for an IR excess in SN2010jl, indicating either new dust formation or the heating of CSM dust in an IR echo. Here we report multi-epoch spectra of SN2010jl that reveal the tell-tale signature of new dust formation: emission-line profiles becoming systematically more blueshifted as the red side of the line is blocked by increasing extinction. The effect is seen clearly in the intermediate-width (400--4000 km/s) component of H$\alpha$ beginning roughly 30d after explosion. Moreover, we present near-IR spectra demonstrating that the asymmetry in the hydrogen-line profiles is wavelength dependent, appearing more pronounced at shorter wavelengths. This evidence suggests that new dust grains had formed quickly in the post-shock shell of SN 2010jl arising from CSM interaction. Since the observed dust temperature has been attributed to an IR echo and not to new dust, either (1) IR excess emission at $\lambda < 5 \mu$m is not a particularly sensitive tracer of new dust formation in SNe, or (2) some assumptions about expected dust temperatures might require further study. Lastly, we discuss one possible mechanism other than dust that might lead to increasingly blueshifted line profiles in SNeIIn, although the wavelength dependence of the asymmetry argues against this hypothesis in the case of SN2010jl.

We report Chandra observations of a double X-ray source in the z=0.1569 galaxy SDSS J171544.05+600835.7. The galaxy was initially identified as a dual AGN candidate based on the double-peaked [O III] emission lines, with a line-of-sight velocity separation of 350 km/s, in its Sloan Digital Sky Survey spectrum. We used the Kast Spectrograph at Lick Observatory to obtain two longslit spectra of the galaxy at two different position angles, which reveal that the two AGN emission components have not only a velocity offset, but also a projected spatial offset of 1.9 kpc/h70 on the sky. Chandra/ACIS observations of two X-ray sources with the same spatial offset and orientation as the optical emission suggest the galaxy most likely contains Compton-thick dual AGN, although the observations could also be explained by AGN jets. Deeper X-ray observations that reveal Fe K lines, if present, would distinguish between the two scenarios. The observations of a double X-ray source in SDSS J171544.05+600835.7 are a proof of concept for a new, systematic detection method that selects promising dual AGN candidates from ground-based spectroscopy that exhibits both velocity and spatial offsets in the AGN emission features.

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 use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

Scaling relations of observed galaxy cluster properties are useful tools for constraining cosmological parameters as well as cluster formation histories. One of the key cosmological parameters, sigma8, is constrained using observed clusters of galaxies, although current estimates of sigma8 from the scaling relations of dynamically relaxed galaxy clusters are limited by the large scatter in the observed cluster mass-temperature (M-T) relation. With a sample of eight strong lensing clusters at 0.3 < z <0.8, we find that the observed cluster concentration-mass relation can be used to reduce the M-T scatter by a factor of 6. Typically only relaxed clusters are used to estimate sigma8, but combining the cluster concentration-mass relation with the M-T relation enables the inclusion of unrelaxed clusters as well. Thus, the resultant gains in the accuracy of sigma8 measurements from clusters are twofold: the errors on sigma8 are reduced and the cluster sample size is increased. Therefore, the statistics on sigma8 determination from clusters are greatly improved by the inclusion of unrelaxed clusters. Exploring cluster scaling relations further, we find that the correlation between brightest cluster galaxy (BCG) luminosity and cluster mass offers insight into the assembly histories of clusters. We find preliminary evidence for a steeper BCG luminosity - cluster mass relation for strong lensing clusters than the general cluster population, hinting that strong lensing clusters may have had more active merging histories.

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

We present strong evidence for dual active galactic nuclei (AGN) in the z=0.36 galaxy COSMOS J100043.15+020637.2. COSMOS Hubble Space Telescope (HST) imaging of the galaxy shows a tidal tail, indicating that the galaxy recently underwent a merger, as well as two bright point sources near the galaxy's center. Both the luminosities of these sources (derived from the HST image) and their emission line flux ratios (derived from Keck/DEIMOS slit spectroscopy) suggest that both are AGN and not star-forming regions or supernovae. Observations from zCOSMOS, Sloan Digital Sky Survey, XMM-Newton, Very Large Array, and Spitzer fortify the evidence for AGN activity. With HST imaging we measure a projected spatial offset between the two AGN of 1.75 +- 0.03 kpc/h, and with DEIMOS we measure a 150 +- 40 km/s line-of-sight velocity offset between the two AGN. Combined, these observations provide substantial evidence that COSMOS J100043.15+020637.2 is a dual AGN in a merger-remnant galaxy.

We present a new technique for observationally identifying galaxy mergers spectroscopically rather than through host galaxy imaging. Our technique exploits the dynamics of supermassive black holes (SMBHs) powering active galactic nuclei (AGNs) in merger-remnant galaxies. Because structure in the universe is built up through galaxy mergers and nearly all galaxies host a central SMBH, some galaxies should possess two SMBHs near their centers as the result of a recent merger. These SMBHs spiral to the center of the resultant merger-remnant galaxy, and one or both of the SMBHs may power AGNs. Using the DEEP2 Galaxy Redshift Survey, we have examined 1881 red galaxies, of which 91 exhibit [O III] and Hbeta emission lines indicative of Seyfert 2 activity. Of these, 32 AGNs have [O III] emission-line redshifts significantly different from the redshifts of the host galaxies' stars, corresponding to velocity offsets of ~50 km/s to ~300 km/s. Two of these AGNs exhibit double-peaked [O III] emission lines, while the remaining 30 AGNs each exhibit a single set of velocity-offset [O III] emission lines. After exploring a variety of physical models for these velocity offsets, we argue that the most likely explanation is inspiralling SMBHs in merger-remnant galaxies. Based on this interpretation, we find that roughly half of the red galaxies hosting AGNs are also merger remnants, which implies that mergers may trigger AGN activity in red galaxies. The AGN velocity offsets we find imply a merger fraction of ~30% and a merger rate of ~3 mergers/Gyr for red galaxies at redshifts 0.34 < z < 0.82.

While both isothermal and NFW-based mass models for galaxy clusters are widely adopted in strong-lensing studies, they cannot easily be distinguished based solely on observed positions of arcs and arclets. We compare the magnifications predicted for giant arcs obtained from isothermal and NFW profiles, taking axially-symmetric and asymmetric mass distributions into account. We find that arc magnifications can differ strongly between the two types of density profiles even if the image morphology is well reproduced. Magnifications by lenses with NFW density profiles are usually larger than those for lenses with singular or nearly singular isothermal density profiles, unless the latter have large cores. Asymmetries play an important role. We illustrate our results with the two well-studied clusters MS 2137 and A~370. We confirm earlier results showing that both isothermal and NFW mass models can very well reproduce the observed arcs, radial arcs and other arclets. While the mass model for MS 2137 is not very well constrained, the two types of mass models produce strongly differing critical curves and caustics for A~370. We find that the NFW mass profile is preferred for A~370. We identify new arclet candidates in the field of A~370. Redshift estimates allowed by the lens model are consistent with results in the literature, where available. Three newly found counter-images are suggested to arise from an active, dust-enshrouded star-forming galaxy at $z\approx1.1$.

The properties of clusters of galaxies offer key insights into the assembly process of structure in the universe. Numerical simulations of cosmic structure formation in a hierarchical, dark matter dominated universe suggest that galaxy cluster concentrations, which are a measure of a halo's central density, decrease gradually with virial mass. However, cluster observations have yet to confirm this correlation. The slopes of the run of measured concentrations with virial mass are often either steeper or flatter than predicted by simulations. In this work, we present the most complete sample of observed cluster concentrations and masses yet assembled, including new measurements for 10 strong lensing clusters, thereby more than doubling the existing number of strong lensing concentration estimates. We fit a power law to the observed concentrations as a function of virial mass, and find that the slope is consistent with the slopes found in simulations, though our normalization factor is higher. Observed lensing concentrations appear to be systematically larger than X-ray concentrations, a more pronounced effect than found in simulations. We also find that at fixed mass, the bulk of observed cluster concentrations are distributed log-normally, with the exception of a few anomalously high concentration clusters. We examine the physical processes likely responsible for the discrepancy between lensing and X-ray concentrations, and for the anomalously high concentrations in particular. The forthcoming Millennium simulation results will offer the most comprehensive comparison set to our findings of an observed concentration-mass power law relation.

Although N-body simulations of cosmic structure formation suggest that dark matter halos have density profiles shallower than isothermal at small radii and steeper at large radii, whether observed galaxy clusters follow this profile is still ambiguous. We use one such density profile, the asymmetric NFW profile, to model the mass distributions of 11 galaxy clusters with gravitational arcs observed by HST. We characterize the galaxy lenses in each cluster as NFW ellipsoids, each defined by an unknown scale convergence, scale radius, ellipticity, and position angle. For a given set of values of these parameters, we compute the arcs that would be produced by such a lens system. To define the goodness of fit to the observed arc system, we define a chi^2 function encompassing the overlap between the observed and reproduced arcs as well as the agreement between the predicted arc sources and the observational constraints on the source system. We minimize this chi^2 to find the values of the lens parameters that best reproduce the observed arc system in a given cluster. Here we report our best-fit lens parameters and corresponding mass estimates for each of the 11 lensing clusters. We find that cluster mass models based on lensing galaxies defined as NFW ellipsoids can accurately reproduce the observed arcs, and that the best-fit parameters to such a model fall within the reasonable ranges defined by simulations. These results assert NFW profiles as an effective model for the mass distributions of observed clusters.

Recent discoveries by the Sloan Digital Sky Survey (SDSS) of four bright redshift z=6 quasars could constrain the mechanism by which the supermassive black holes powering these sources are assembled. Here we compute the probability that the fluxes of the quasars are strongly amplified by gravitational lensing, and therefore the likelihood that the black hole masses are overestimated when they are inferred assuming Eddington luminosities. The poorly-constrained shape of the intrinsic quasar luminosity function (LF) at redshift z=6 results in a large range of possible lensing probabilities. If the LF is either steep, or extends to faint magnitudes, the probability for amplification by a factor of ten ore more (and with only one image detectable by SDSS) can reach essentially 100 percent. We show that future observations, in particular, either of the current four quasars at the high angular resolution provided by the Hubble Space Telescope, or of an increased sample of about twenty redshift six quasars at the current angular resolution, should either discover several gravitational lenses, or else provide interesting new constraints on the shape of the z=6 quasar LF.

We show that the observed mass-to-light ratio of galaxy clusters increases with cluster temperature as expected from cosmological simulations. Contrary to previous observational suggestions, we find a mild but robust increase of M/L from poor (T = 1 - 2 keV) to rich (T = 12 keV) clusters; over this range, the mean M/L_V increases by a factor of about 2. The best fit relation satisfies M/L_V = (170 +- 30) T(keV)^(0.3 +- 0.1) h at z=0, with a large scatter. This trend confirms predictions from cosmological simulations which show that the richest clusters are antibiased, with a higher ratio of mass per unit light than average. The antibias increases with cluster temperature. The effect is caused by the relatively older age of the high-density clusters, where light has declined more significantly than average since their earlier formation time. Combining the current observations with simulations, we find a global value of M/L_V = 240 +- 50 h, and a corresponding mass density of the universe of Omega_m = 0.17 +- 0.05.