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).

Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsecs. Here we report optical integral field observations of the massive but compact galaxy SDSS J211824.06+001729.4. The oxygen [O II] lines at wavelengths of 3726 and 3729 angstroms reveal an ionized outflow spanning 80 by 100 square kiloparsecs, depositing metal-enriched gas at 10,000 kelvin through an hourglass-shaped nebula that resembles an evacuated and limb-brightened bipolar bubble. We also observe neutral gas phases at temperatures of less than 10,000 kelvin reaching distances of 20 kiloparsecs and velocities of around 1,500 kilometres per second. This multi-phase outflow is probably driven by bursts of star formation, consistent with theory.

We present a sample of 48 nearby galaxies with central, biconical outflows identified by the Mapping Nearby Galaxies at APO (MaNGA) survey. All considered galaxies have star formation driven bi-conical central outflows (SFB), with no signs of AGN. We find that the SFB outflows require high central concentration of the star formation rate. This increases the gas velocity dispersion over the equilibrium limit and helps maintain the gas outflows. The central starbursts increase the metallicity, extinction, and the alpha/Fe ratio in the gas. Significant amount of young stellar population at the centers suggests that the SFBs are associated with the formation of young bulges in galaxies. More than 70% of SFB galaxies are group members or have companions with no prominent interaction, or show asymmetry of external isophotes. In 15% SFB cases stars and gas rotate in the opposite directions, which points at the gas infall from satellites as the primary reason for triggering the SFB phenomena.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

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.

We present Atacama Large Millimeter/submillimeter Array observations of a massive (M_stars~10^11 M_Sun) compact (r_e,UV~100 pc) merger remnant at z=0.66 that is driving a 1000 km/s outflow of cool gas, with no observational trace of an active galactic nucleus (AGN). We resolve molecular gas on scales of approximately 1-2 kpc, and our main finding is the discovery of a wing of blueshifted CO(2-1) emission out to -1000 km/s relative to the stars. We argue that this is the molecular component of a multiphase outflow, expelled from the central starburst within the past 5 Myr through stellar feedback, although we cannot rule out previous AGN activity as a launching mechanism. If the latter is true, then this is an example of a relic multiphase AGN outflow. We estimate a molecular mass outflow rate of approximately 300 M_Sun/yr, or about one third of the 10 Myr-averaged star formation rate. This system epitomizes the multiphase 'blowout' episode following a dissipational major merger - a process that has violently quenched central star formation and supermassive black hole growth.

The circumgalactic medium (CGM) close to ~L* star-forming galaxies hosts strong MgII 2796 absorption (with equivalent width W_2796>0.1 Ang) with a near-unity covering fraction. To characterize the spatial coherence of this absorption, we analyze the W_2796 distribution in the CGM of 27 star-forming galaxies detected in deep spectroscopy of bright background (b/g) galaxies first presented in Rubin et al. (2018). The sample foreground (f/g) systems have redshifts 0.35<z<0.8 and stellar masses 9.1<log M_*/M_sun<11.1, and the b/g galaxies provide spatially-extended probes with half-light radii 1.0 kpc<R_eff<7.9 kpc at projected distances R_perp<50 kpc. Our analysis also draws on literature W_2796 values measured in b/g QSO spectroscopy probing the halos of f/g galaxies with a similar range in M_* at z ~ 0.25. By making the assumptions that (1) samples of like galaxies exhibit similar circumgalactic W_2796 distributions; and that (2) the quantity log W_2796 has a Gaussian distribution with a dispersion that is constant with M_* and R_perp, we use this QSO-galaxy pair sample to construct a model for the log W_2796 distribution in the CGM. We then demonstrate the dependence of this distribution on the ratio of the surface area of the b/g probe to the projected absorber surface area (x_A=A_G/A_A), finding that distributions which assume x_A>=15 are statistically inconsistent with that observed toward our b/g galaxies at a 95% confidence level. This limit, in combination with the b/g galaxy sizes, requires that the length scale over which W_2796 does not vary (the "coherence scale" of MgII absorption) is l_A>1.9 kpc. This novel constraint on the morphology of cool, photoionized structures in the inner CGM suggests that either these structures each extend over kiloparsec scales, or that the numbers and velocity dispersion of these structures are spatially correlated over the same scales.

We present XMM-Newton imaging spectroscopy of ten weak emission-line quasars (WLQs) at $0.928\leq z \leq 3.767$, six of which are radio quiet and four which are radio intermediate. The new X-ray data enabled us to measure the power-law photon index, at rest-frame energies $>2$ keV, in each source with relatively high accuracy. These measurements allowed us to confirm previous reports that WLQs have steeper X-ray spectra, suggesting higher accretion rates with respect to 'typical' quasars. A comparison between the photon indices of our radio-quiet WLQs and those of a control sample of 85 sources shows that the first are significantly higher, at the >~$3\sigma$ level. Collectively, the four radio-intermediate WLQs have lower photon indices with respect to the six radio-quiet WLQs, as may be expected if the spectra of the first group are contaminated by X-ray emission from a jet. Therefore, in the absence of significant jet emission along our line of sight, these results are in agreement with the idea that WLQs constitute the extreme high end of the accretion rate distribution in quasars. We detect soft excess emission in our lowest-redshift radio-quiet WLQ, in agreement with previous findings suggesting that the prominence of this feature is associated with a high accretion rate. We have not detected signatures of Compton reflection, Fe K$\alpha$ lines, or strong variability between two X-ray epochs in any of our WLQs, which can be attributed to their relatively high luminosity.

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.

Spectroscopy of background QSO sightlines passing close to foreground galaxies is a potent technique for studying the circumgalactic medium (CGM). QSOs are effectively point sources, however, limiting their potential to constrain the size of circumgalactic gaseous structures. Here we present the first large Keck/LRIS and VLT/FORS2 spectroscopic survey of bright (B_AB < 22.3) background galaxies whose lines of sight probe MgII 2796, 2803 absorption from the CGM around close projected foreground galaxies at transverse distances 10 kpc < R_perp < 150 kpc. Our sample of 72 projected pairs, drawn from the PRIsm MUlti-object Survey (PRIMUS), includes 48 background galaxies which do not host bright AGN, and both star-forming and quiescent foreground galaxies with stellar masses 9.0 < log M_*/M_sun < 11.2 at redshifts 0.35 < z_f/g < 0.8. We detect MgII absorption associated with these foreground galaxies with equivalent widths 0.25 Ang < W_2796 < 2.6 Ang at >2sigma significance in 20 individual background sightlines passing within R_perp < 50 kpc, and place 2sigma upper limits on W_2796 of <0.5 Ang in an additional 11 close sightlines. Within R_perp < 50 kpc, W_2796 is anticorrelated with R_perp, consistent with analyses of MgII absorption detected along background QSO sightlines. Subsamples of these foreground hosts divided at log M_*/M_sun = 9.9 exhibit statistically inconsistent W_2796 distributions at 30 kpc < R_perp < 50 kpc, with the higher-M_* galaxies yielding a larger median W_2796 by 0.9 Ang. Finally, we demonstrate that foreground galaxies with similar stellar masses exhibit the same median W_2796 at a given R_perp to within <0.2 Ang toward both background galaxies and toward QSO sightlines drawn from the literature. Analysis of these datasets constraining the spatial coherence scale of circumgalactic MgII absorption is presented in a companion paper.

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.

We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We present a study of the kinematics of the extraplanar ionized gas around several dozen galaxies observed by the Mapping of Nearby Galaxies at the Apache Point Observatory (MaNGA) survey. We considered a sample of 67 edge-on galaxies out of more than 1400 extragalactic targets observed by MaNGA, in which we found 25 galaxies (or 37%) with regular lagging of the rotation curve at large distances from the galactic midplane. We model the observed $H\alpha$ emission velocity fields in the galaxies, taking projection effects and a simple model for the dust extinction into the account. We show that the vertical lag of the rotation curve is necessary in the modeling, and estimate the lag amplitude in the galaxies. We find no correlation between the lag and the star formation rate in the galaxies. At the same time, we report a correlation between the lag and the galactic stellar mass, central stellar velocity dispersion, and axial ratio of the light distribution. These correlations suggest a possible higher ratio of infalling-to-local gas in early-type disk galaxies or a connection between lags and the possible presence of hot gaseous halos, which may be more prevalent in more massive galaxies. These results again demonstrate that observations of extraplanar gas can serve as a potential probe for accretion of gas.

Diffuse Ionized Gas (DIG) is prevalent in star-forming galaxies. Using a sample of 365 nearly face-on star-forming galaxies observed by MaNGA, we demonstrate how DIG in star-forming galaxies impacts the measurements of emission line ratios, hence the interpretation of diagnostic diagrams and gas-phase metallicity measurements. At fixed metallicity, DIG-dominated low H\alpha\ surface brightness regions display enhanced [SII]/H\alpha, [NII]/H\alpha, [OII]/H\beta, and [OI]/H\alpha. The gradients in these line ratios are determined by metallicity gradients and H\alpha\ surface brightness. In line ratio diagnostic diagrams, contamination by DIG moves HII regions towards composite or LI(N)ER-like regions. A harder ionizing spectrum is needed to explain DIG line ratios. Leaky HII region models can only shift line ratios slightly relative to HII region models, and thus fail to explain the composite/LI(N)ER line ratios displayed by DIG. Our result favors ionization by evolved stars as a major ionization source for DIG with LI(N)ER-like emission. DIG can significantly bias the measurement of gas metallicity and metallicity gradients derived using strong-line methods. Metallicities derived using N2O2 are optimal because they exhibit the smallest bias and error. Using O3N2, R23, N2=[NII]/H\alpha, and N2S2H\alpha\ (Dopita et al. 2016) to derive metallicities introduces bias in the derived metallicity gradients as large as the gradient itself. The strong-line method of Blanc et al. (2015; IZI hereafter) cannot be applied to DIG to get an accurate metallicity because it currently contains only HII region models which fail to describe the DIG.

Using spatially resolved spectroscopy from SDSS-IV MaNGA we have demonstrated that low ionisation emission line regions (LIERs) in local galaxies result from photoionisation by hot evolved stars, not active galactic nuclei. LIERs are ubiquitous in both quiescent galaxies and in the central regions of galaxies where star formation takes place at larger radii. We refer to these two classes of galaxies as extended LIER (eLIER) and central LIER (cLIER) galaxies respectively. cLIERs are late type galaxies located around the green valley, in the transition region between the star formation main sequence and quiescent galaxies. These galaxies display regular disc rotation in both stars and gas, although featuring a higher central stellar velocity dispersion than star forming galaxies of the same mass. cLIERs are consistent with being slowly quenched inside-out; the transformation is associated with massive bulges, pointing towards the importance of bulge growth via secular evolution. eLIERs are morphologically early types and are indistinguishable from passive galaxies devoid of line emission in terms of their stellar populations, morphology and central stellar velocity dispersion. Ionised gas in eLIERs shows both disturbed and disc-like kinematics. When a large-scale flow/rotation is observed in the gas, it is often misaligned relative to the stellar component. These features indicate that eLIERs are passive galaxies harbouring a residual cold gas component, acquired mostly via external accretion. Importantly, quiescent galaxies devoid of line emission reside in denser environments and have significantly higher satellite fraction than eLIERs. Environmental effects thus represent the likely cause for the existence of line-less galaxies on the red sequence.

We have conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above the midplane. With 16 galaxies we can extend this analysis out to about 9 kpc, i.e. a distance of ~2R_e, vertically from the midplane. In the halo, the surface brightnesses of the [OII] and Halpha emission lines are comparable, unlike in the disk where Halpha dominates. When we split the sample by specific star formation rate, concentration index, and stellar mass, each subsample's emission line surface brightness profiles and ratios differ, indicating that extra-planar gas properties can vary. The emission line surface brightnesses of the gas around high specific star formation rate galaxies are higher at all distances, and the line ratios are closer to ratios characteristic of HII regions compared with low specific star formation rate galaxies. The less concentrated and lower stellar mass samples exhibit line ratios that are more like HII regions at larger distances than their more concentrated and higher stellar mass counterparts. The largest difference between different subsamples occurs when the galaxies are split by stellar mass. We additionally infer that gas far from the midplane in more massive galaxies has the highest temperatures and steepest radial temperature gradients based on their [NII]/Halpha and [OII]/Halpha ratios between the disk and the halo.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

We measure the clustering of X-ray, radio, and mid-IR-selected active galactic nuclei (AGN) at 0.2 < z < 1.2 using multi-wavelength imaging and spectroscopic redshifts from the PRIMUS and DEEP2 redshift surveys, covering 7 separate fields spanning ~10 square degrees. Using the cross-correlation of AGN with dense galaxy samples, we measure the clustering scale length and slope, as well as the bias, of AGN selected at different wavelengths. Similar to previous studies, we find that X-ray and radio AGN are more clustered than mid-IR-selected AGN. We further compare the clustering of each AGN sample with matched galaxy samples designed to have the same stellar mass, star formation rate, and redshift distributions as the AGN host galaxies and find no significant differences between their clustering properties. The observed differences in the clustering of AGN selected at different wavelengths can therefore be explained by the clustering differences of their host populations, which have different distributions in both stellar mass and star formation rate. Selection biases inherent in AGN selection, therefore, determine the clustering of observed AGN samples. We further find no significant difference between the clustering of obscured and unobscured AGN, using IRAC or WISE colors or X-ray hardness ratio.

We present results on gas flows in the halo of a Milky Way-like galaxy at z=0.413 based on high-resolution spectroscopy of a background galaxy. This is the first study of circumgalactic gas at high spectral resolution towards an extended background source (i.e., a galaxy rather than a quasar). Using longslit spectroscopy of the foreground galaxy, we observe spatially extended H alpha emission with circular rotation velocity v=270 km/s. Using echelle spectroscopy of the background galaxy, we detect Mg II and Fe II absorption lines at impact parameter rho=27 kpc that are blueshifted from systemic in the sense of the foreground galaxy's rotation. The strongest absorber EW(2796) = 0.90 A has an estimated column density (N_H>10^19 cm-2) and line-of-sight velocity dispersion (sigma=17 km/s) that are consistent with the observed properties of extended H I disks in the local universe. Our analysis of the rotation curve also suggests that this r=30 kpc gaseous disk is warped with respect to the stellar disk. In addition, we detect two weak Mg II absorbers in the halo with small velocity dispersions (sigma<10 km/s). While the exact geometry is unclear, one component is consistent with an extraplanar gas cloud near the disk-halo interface that is co-rotating with the disk, and the other is consistent with a tidal feature similar to the Magellanic Stream. We can place lower limits on the cloud sizes (l>0.4 kpc) for these absorbers given the extended nature of the background source. We discuss the implications of these results for models of the geometry and kinematics of gas in the circumgalactic medium.

Many astronomers now participate in large international collaborations, and it is important to examine whether these structures foster a scientific climate that is inclusive and diverse. The Committee on the Participation of Women in the Sloan Digital Sky Survey (CPWS) was formed to evaluate the demographics and gender climate within SDSS-IV, one of the largest and most geographically distributed astronomical collaborations. In April 2014, the CPWS administered a demographic survey to establish a baseline for the incipient SDSS-IV. We received responses from 250 participants (46% of the active membership). Half of the survey respondents were located in the US or Canada and 30% were based in Europe. Eleven percent of survey respondents considered themselves to be an ethnic minority at their current institution. Twenty-five percent of the SDSS-IV collaboration members are women, a fraction that is consistent with the US astronomical community, but substantially higher than the fraction of women in the IAU (16%). Approximately equal fractions of men and women report holding positions of leadership. When binned by academic age and career level, men and women also assume leadership roles at approximately equal rates, in a way that increases steadily for both genders with increasing seniority. In this sense, SDSS-IV has been successful in recruiting leaders that are representative of the collaboration. Yet, more progress needs to be made towards achieving gender balance and increasing diversity in the field of astronomy, and there is still room for improvement in the membership and leadership of SDSS-IV. For example, at the highest level of SDSS-IV leadership, women disproportionately assume roles related to education and public outreach. The CPWS plans to use these initial data to establish a baseline for tracking demographics over time as we work to assess and improve the climate of SDSS-IV.

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

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

Recent observations have revealed that starburst galaxies can drive molecular gas outflows through stellar radiation pressure. Molecular gas is the phase of the interstellar medium from which stars form, so these outflows curtail stellar mass growth in galaxies. Previously known outflows, however, involve small fractions of the total molecular gas content and are restricted to sub-kiloparsec scales. It is also apparent that input from active galactic nuclei is in at least some cases dynamically important, so pure stellar feedback has been considered incapable of aggressively terminating star formation on galactic scales. Extraplanar molecular gas has been detected in the archetype starburst galaxy M82, but so far there has been no evidence that starbursts can propel significant quantities of cold molecular gas to the same galactocentric radius (~10 kpc) as the warmer gas traced by metal absorbers. Here we report observations of molecular gas in a compact (effective radius 100 pc) massive starburst galaxy at z=0.7, which is known to drive a fast outflow of ionized gas. We find that 35 per cent of the total molecular gas is spatially extended on a scale of approximately 10 kpc, and one third of this has a velocity of up to 1000 km/s. The kinetic energy associated with this high-velocity component is consistent with the momentum flux available from stellar radiation pressure. This result demonstrates that nuclear bursts of star formation are capable of ejecting large amounts of cold gas from the central regions of galaxies, thereby strongly affecting their evolution.

We investigate the process of rapid star formation quenching in a sample of 12 massive galaxies at intermediate redshift (z~0.6) that host high-velocity ionized gas outflows (v>1000 km/s). We conclude that these fast outflows are most likely driven by feedback from star formation rather than active galactic nuclei (AGN). We use multiwavelength survey and targeted observations of the galaxies to assess their star formation, AGN activity, and morphology. Common attributes include diffuse tidal features indicative of recent mergers accompanied by bright, unresolved cores with effective radii less than a few hundred parsecs. The galaxies are extraordinarily compact for their stellar mass, even when compared with galaxies at z~2-3. For 9/12 galaxies, we rule out an AGN contribution to the nuclear light and hypothesize that the unresolved core comes from a compact central starburst triggered by the dissipative collapse of very gas-rich progenitor merging disks. We find evidence of AGN activity in half the sample but we argue that it accounts for only a small fraction (<10%) of the total bolometric luminosity. We find no correlation between AGN activity and outflow velocity and we conclude that the fast outflows in our galaxies are not powered by on-going AGN activity, but rather by recent, extremely compact starbursts.

We present a study of Spitzer/IRAC and X-ray active galactic nuclei (AGNs) selection techniques in order to quantify the overlap, uniqueness, contamination, and completeness of each. We investigate how the overlap and possible contamination of the samples depends on the IR and X-ray depths. We use Spitzer/IRAC imaging, Chandra and XMM X-ray imaging, and PRism MUlti-object Survey (PRIMUS) spectroscopic redshifts to construct galaxy and AGN samples at 0.2<z<1.2 over 8 deg^2. We construct samples over a wide range of IRAC flux limits (SWIRE to GOODS depth) and X-ray flux limits (10 ks to 2 Ms). We compare IR-AGN samples defined using the IRAC color selection of Stern et al. and Donley et al. with X-ray detected AGN samples. For roughly similar depth IR and X-ray surveys, we find that ~75% of IR-AGN are identified as X-ray AGN. This fraction increases to ~90% when comparing against the deepest X-ray data, indicating that only ~10% of IR-selected AGN may be heavily obscured. The IR-AGN selection proposed by Stern et al. suffers from contamination by star-forming galaxies at various redshifts when using deeper IR data, though the selection technique works well for shallow IR data. While similar overall, the IR-AGN samples preferentially contain more luminous AGN, while the X-ray AGN samples preferentially contain lower specific accretion rate AGN, where the host galaxy light dominates at IR wavelengths. The host galaxy populations of the IR and X-ray AGN samples have similar restframe colors and stellar masses; both selections identify AGN in blue, star-forming and red, quiescent galaxies.

We present an observationally motivated model to connect the AGN and galaxy populations at 0.2<z<1.0 and predict the AGN X-ray luminosity function (XLF). We start with measurements of the stellar mass function of galaxies (from the Prism Multi-object Survey) and populate galaxies with AGNs using models for the probability of a galaxy hosting an AGN as a function of specific accretion rate. Our model is based on measurements indicating that the specific accretion rate distribution is a universal function across a wide range of host stellar mass with slope gamma_1 = -0.65 and an overall normalization that evolves with redshift. We test several simple assumptions to extend this model to high specific accretion rates (beyond the measurements) and compare the predictions for the XLF with the observed data. We find good agreement with a model that allows for a break in the specific accretion rate distribution at a point corresponding to the Eddington limit, a steep power-law tail to super-Eddington ratios with slope gamma_2 = -2.1 +0.3 -0.5, and a scatter of 0.38 dex in the scaling between black hole and host stellar mass. Our results show that samples of low luminosity AGNs are dominated by moderately massive galaxies (M* ~ 10^{10-11} M_sun) growing with a wide range of accretion rates due to the shape of the galaxy stellar mass function rather than a preference for AGN activity at a particular stellar mass. Luminous AGNs may be a severely skewed population with elevated black hole masses relative to their host galaxies and in rare phases of rapid accretion.

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

We report observations of the CO(2-1) emission of SDSSJ1506+54, a compact (r_e~135pc) starburst galaxy at z=0.6. SDSSJ1506+54 appears to be forming stars close to the limit allowed by stellar radiation pressure feedback models: the measured L_IR/L'_CO 1500 is one of the highest measured for any galaxy. With its compact optical morphology but extended low surface brightness envelope, post-starburst spectral features, high infrared luminosity (L_IR>10^12.5 L_Sun), low gas fraction (M_H2/M_stars~15%), and short gas depletion time (tens of Myr), we speculate that this is a feedback- limited central starburst episode at the conclusion of a major merger. Taken as such, SDSSJ1504+54 epitomizes the brief closing stage of a classic model of galaxy growth: we are witnessing a key component of spheroid formation during what we term a 'redline' starburst.

Local luminous infrared (IR) galaxies (LIRGs) have both high star formation rates (SFR) and a high AGN (Seyfert and AGN/starburst composite) incidence. Therefore, they are ideal candidates to explore the co-evolution of black hole (BH) growth and star formation (SF) activity, not necessarily associated with major mergers. Here, we use Spitzer/IRS spectroscopy of a complete volume-limited sample of local LIRGs (distances of <78Mpc). We estimate typical BH masses of 3x10^7 M_sun using [NeIII]15.56micron and optical [OIII]5007A gas velocity dispersions and literature stellar velocity dispersions. We find that in a large fraction of local LIRGs the current SFR is taking place not only in the inner nuclear ~1.5kpc region, as estimated from the nuclear 11.3micron PAH luminosities, but also in the host galaxy. We next use the ratios between the SFRs and BH accretion rates (BHAR) to study whether the SF activity and BH growth are contemporaneous in local LIRGs. On average, local LIRGs have SFR to BHAR ratios higher than those of optically selected Seyferts of similar AGN luminosities. However, the majority of the IR-bright galaxies in the RSA Seyfert sample behave like local LIRGs. Moreover, the AGN incidence tends to be higher in local LIRGs with the lowest SFRs. All this suggests that in local LIRGs there is a distinct IR-bright star forming phase taking place prior to the bulk of the current BH growth (i.e., AGN phase). The latter is reflected first as a composite and then as a Seyfert, and later as a non-LIRG optically identified Seyfert nucleus with moderate SF in its host galaxy.

We present optical spectroscopy of a sample of 38 post-starburst quasars (PSQs) at z ~ 0.3, 29 of which have morphological classifications based on Hubble Space Telescope imaging. These broad-lined active galactic nuclei (AGNs) possess the spectral signatures of massive intermediate-aged stellar populations making them potentially useful for studying connections between nuclear activity and host galaxy evolution. We model the spectra in order to determine the ages and masses of the host stellar populations, and the black hole masses and Eddington fractions of the AGNs. Our model components include an instantaneous starburst, a power-law, and emission lines. We find the PSQs have MBH ~ 10^8 Msun accreting at a few percent of Eddington luminosity and host ~ 10^10.5 Msun stellar populations which are several hundred Myr to a few Gyr old. We investigate relationships among these derived properties, spectral properties, and morphologies. We find that PSQs hosted in spiral galaxies have significantly weaker AGN luminosities, older starburst ages, and narrow emission-line ratios diagnostic of ongoing star-formation when compared to their early-type counterparts. We conclude that the early-type PSQs are likely the result of major mergers and were likely luminous infrared galaxies in the past, while spiral PSQs with more complex star-formation histories are triggered by less dramatic events (e.g., harassment, bars). We provide diagnostics to distinguish the early-type and spiral hosts when high spatial resolution imaging is not available.

We present the discovery of compact, obscured star formation in galaxies at z ~ 0.6 that exhibit >1000 km/s outflows. Using optical morphologies from the Hubble Space Telescope and infrared photometry from the Wide-field Infrared Survey Explorer, we estimate star formation rate (SFR) surface densities that approach Sigma_SFR ~ 3000 Msun/yr/kpc^2, comparable to the Eddington limit from radiation pressure on dust grains. We argue that feedback associated with a compact starburst in the form of radiation pressure from massive stars and ram pressure from supernovae and stellar winds is sufficient to produce the high-velocity outflows we observe, without the need to invoke feedback from an active galactic nucleus.

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

(Abridged) We report on the X-ray and multiwavelength properties of 11 radio-quiet quasars with weak or no emission lines identified by the Sloan Digital Sky Survey (SDSS) with redshift z=0.4-2.5. The distribution of relative X-ray brightness for our low-redshift weak-line quasar (WLQ) candidates is significantly different from that of typical radio-quiet quasars, having an excess of X-ray weak sources, but it is consistent with that of high-redshift WLQs. The X-ray weak sources generally show similar UV emission-line properties to those of the X-ray weak quasar PHL 1811; they may belong to the notable class of PHL 1811 analogs. The average X-ray spectrum of these sources is somewhat harder than that of typical radio-quiet quasars. Several other low-redshift WLQ candidates have normal ratios of X-ray-to-optical/UV flux, and their average X-ray spectral properties are also similar to those of typical radio-quiet quasars. The X-ray weak and X-ray normal WLQ candidates may belong to the same subset of quasars having high-ionization "shielding gas" covering most of the wind-dominated broad emission-line region, but be viewed at different inclinations. The mid-infrared-to-X-ray spectral energy distributions (SEDs) of these sources are generally consistent with those of typical SDSS quasars, showing that they are not likely to be BL Lac objects with relativistically boosted continua and diluted emission lines. However, one source in our X-ray observed sample is remarkably strong in X-rays, indicating that a small fraction of low-redshift WLQ candidates may actually be BL Lacs residing in the radio-faint tail of the BL Lac population. We also investigate universal selection criteria for WLQs over a wide range of redshift, finding that it is not possible to select WLQ candidates in a fully consistent way using different prominent emission lines as a function of redshift.

We present estimates of black hole accretion rates and nuclear, extended, and total star-formation rates for a complete sample of Seyfert galaxies. Using data from the Spitzer Space Telescope, we measure the active galactic nucleus (AGN) luminosity using the [O IV] 25.89 micron emission line and the star-forming luminosity using the 11.3 micron aromatic feature and extended 24 micron continuum emission. We find that black hole growth is strongly correlated with nuclear (r<1 kpc) star formation, but only weakly correlated with extended (r>1 kpc) star formation in the host galaxy. In particular, the nuclear star-formation rate (SFR) traced by the 11.3 micron aromatic feature follows a relationship with the black hole accretion rate (BHAR) of the form SFR\proptoBHAR^0.8, with an observed scatter of 0.5 dex. This SFR-BHAR relationship persists when additional star formation in physically matched r=1 kpc apertures is included, taking the form SFR\proptoBHAR^0.6. However, the relationship becomes almost indiscernible when total SFRs are considered. This suggests a physical connection between the gas on sub-kpc and sub-pc scales in local Seyfert galaxies that is not related to external processes in the host galaxy. It also suggests that the observed scaling between star formation and black hole growth for samples of AGNs will depend on whether the star formation is dominated by a nuclear or extended component. We estimate the integrated black hole and bulge growth that occurs in these galaxies and find that an AGN duty cycle of 5-10% would maintain the ratio between black hole and bulge masses seen in the local universe.

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

We present Spitzer Space Telescope photometry of 18 Sloan Digital Sky Survey (SDSS) quasars at 2.7 <= z <= 5.9 which have weak or undetectable high-ionization emission lines in their rest-frame ultraviolet (UV) spectra (hereafter weak-lined quasars, or WLQs). The Spitzer data are combined with SDSS spectra and ground-based, near-infrared (IR) photometry of these sources to produce a large inventory of spectral energy distributions (SEDs) of WLQs across the rest-frame ~0.1-5 mum spectral band. The SEDs of our sources are inconsistent with those of BL Lacertae objects which are dominated by synchrotron emission due to a jet aligned close to our line-of-sight, but are consistent with the SED of ordinary quasars with similar luminosities and redshifts that exhibit a near-to-mid-IR 'bump', characteristic of hot dust emission. This indicates that broad emission lines in WLQs are intrinsically weak, rather than suffering continuum dilution from a jet, and that such sources cannot be selected efficiently from traditional photometric surveys.

We present images of 29 post-starburst quasars (PSQs) from a Hubble Space Telescope (\emph{HST}) Advanced Camera for Surveys (ACS) Wide Field Channel Snapshot program. These broad-lined active galactic nuclei (AGN) possess the spectral signatures of massive ($M_{burst} \sim 10^{10} M_{\odot}$), moderate-aged stellar populations (hundreds of Myrs). Thus, their composite nature provides insight into the AGN-starburst connection. We measure quasar-to-host galaxy light contributions via semi-automated two-dimensional light profile fits of PSF-subtracted images. We examine the host morphologies, as well as, model the separate bulge and disk components. The \emph{HST}/ACS-F606W images reveal an equal number of spiral (13/29) and early-type (13/29) hosts, with the remaining three hosts having indeterminate classifications. AGNs hosted by early-type galaxies have on average greater luminosity than those hosted by spiral galaxies. Disturbances, such as tidal tails, shells, star-forming knots, and asymmetries are seen as signposts of interaction/merger activity. Disturbances such as these were found in 17 of the 29 objects and are evenly distributed among early-type and spiral galaxies. Two of these systems are clearly merging with their companions. Compared to other AGN of similar luminosity and redshift these PSQs have a higher fraction of early-type hosts and disturbances. Our most luminous objects with disturbed early-type host galaxies appear to be consistent with merger products. Thus, these luminous disturbed galaxies may represent a phase in an evolutionary scenario for merger driven activity and of hierarchical galaxy evolution. Our less luminous objects appear to be consistent with Seyfert galaxies not requiring triggering by major mergers. Many of these Seyferts are barred spiral galaxies.

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

We used Spitzer/IRS spectroscopic data on 426 galaxies including quasars, Seyferts, LINER and HII galaxies to investigate the relationship among the mid-IR emission lines. There is a tight linear correlation between the [Ne V]14.3 um and 24.3 um (97.1 eV) and the [O IV]25.9 um (54.9 eV) high-ionization emission lines. The correlation also holds for these high-ionization emission lines and the [Ne III]15.56 um (41 eV) emission line, although only for active galaxies. We used these correlations to calculate the [Ne III] excess due to star formation in Seyfert galaxies. We also estimated the [O IV] luminosity due to star formation in active galaxies and determined that it dominates the [O IV] emission only if the contribution of the active nucleus to the total luminosity is below 5%. We find that the AGN dominates the [O IV] emission in most Seyfert galaxies, whereas star-formation adequately explains the observed [O IV] emission in optically classified HII galaxies. Finally we computed photoionization models to determine the physical conditions of the narrow line region where these high-ionization lines originate. The estimated ionization parameter range is -2.8 < log U < -2.5 and the total hydrogen column density range is 20 < log nH (cm-2) < 21.

We present Spitzer measurements of the aromatic (also known as PAH) features for 35 Seyfert galaxies from the revised Shapley-Ames sample and find that the relative strengths of the features differ significantly from those observed in star-forming galaxies. Specifically, the features at 6.2, 7.7, and 8.6 micron are suppressed relative to the 11.3 micron feature in Seyferts. Furthermore, we find an anti-correlation between the L(7.7 micron)/L(11.3 micron) ratio and the strength of the rotational H2 (molecular hydrogen) emission, which traces shocked gas. This suggests that shocks suppress the short-wavelength features by modifying the structure of the aromatic molecules or destroying the smallest grains. Most Seyfert nuclei fall on the relationship between aromatic emission and [Ne II] emission for star-forming galaxies, indicating that aromatic-based estimates of the star-formation rate in AGN host galaxies are generally reasonable. For the outliers from this relationship, which have small L(7.7 micron)/L(11.3 micron) ratios and strong H2 emission, the 11.3 micron feature still provides a valid measure of the star-formation rate.

We present Gemini-North K-band spectra of two representative members of the class of high-redshift quasars with exceptionally weak rest-frame ultraviolet emission lines (WLQs), SDSS J114153.34+021924.3 at z=3.55 and SDSS J123743.08+630144.9 at z=3.49. In both sources we detect an unusually weak broad H_beta line and we place tight upper limits on the strengths of their [O III] lines. Virial, H_beta-based black-hole mass determinations indicate normalized accretion rates of L/L_Edd=0.4 for these sources, which is well within the range observed for typical quasars with similar luminosities and redshifts. We also present high-quality XMM-Newton imaging spectroscopy of SDSS J114153.34+021924.3 and find a hard-X-ray photon index of Gamma=1.91^{+0.24}_{-0.22} which supports the virial L/L_Edd determination in this source. Our results suggest that the weakness of the broad-emission lines in WLQs is not a consequence of an extreme continuum-emission source but instead due to abnormal broad-emission line region properties.

We present radio and X-ray observations, as well as optical light curves, for a subset of 26 BL Lac candidates from the Sloan Digital Sky Survey (SDSS) lacking strong radio emission and with z<2.2. Half of these 26 objects are shown to be stars, galaxies, or absorbed quasars. We conclude that the other 13 objects are Active Galactic Nuclei (AGN) with abnormally weak emission features; ten of those 13 are definitively radio-quiet, and, for those with available optical light curves, their level of optical flux variability is consistent with radio-quiet quasars. We cannot exclude the possibility that some of these 13 AGN lie on the extremely radio-faint tail of the BL Lac distribution, but our study generally supports the notion that all BL Lac objects are radio-loud. These radio-quiet AGN appear to have intrinsically weak or absent broad emission line regions, and, based on their X-ray properties, we argue that some are low-redshift analogs to weak line quasars (WLQs). SDSS BL Lac searches are so far the only systematic surveys of the SDSS database capable of recovering such exotic low-redshift WLQs. There are 71 more z<2.2 radio-quiet BL Lac candidates already identified in the SDSS not considered here, and many of those might be best unified with WLQs as well. Future studies combining low- and high-redshift WLQ samples will yield new insight on our understanding of the structure and formation of AGN broad emission line regions.

Type 2 AGNs with intrinsically weak broad emission lines (BELs) would be exceptions to the unified model. After examining a number of proposed candidates critically, we find that the sample is contaminated significantly by objects with BELs of strengths indicating that they actually contain intermediate-type AGNs, plus a few Compton-thick sources as revealed by extremely low ratios of X-ray to nuclear IR luminosities. We develop quantitative metrics that show two (NGC 3147 and NGC 4594) of the remaining candidates to have BELs 2-3 orders of magnitude weaker than those of typical type-1 AGNs. Several more galaxies remain as candidates to have anomalously weak BELs, but this status cannot be confirmed with the existing information. Although the parent sample is poorly defined, the two confirmed objects are well under 1% of its total number of members, showing that the absence of a BEL is possible, but very uncommon in AGN. We evaluate these two objects in detail using multi-wavelength measurements. They have little X-ray extinction with N_H < 10^21 cm^{-2}. Their IR spectra show strong silicate emission (NGC 4594) or weak aromatic features on a generally power law continuum with a suggestion of silicates in emission (NGC 3147). No polarized BEL is detected in NGC 3147. These results indicate that the two unobscured type-2 objects have circumnuclear tori that are approximately face-on. Combined with their X-ray and optical/UV properties, this behavior implies that we have an unobscured view of the nuclei and thus that they have intrinsically weak BELs. We compare their properties with those of the other less-extreme candidates. We then compare the distributions of bolometric luminosities and accretion rates of these objects with theoretical models that predict weak BELs.

We present a sample of 723 optically selected BL Lac candidates from the SDSS DR7 spectroscopic database encompassing 8250 deg^2 of sky; our sample constitutes one of the largest uniform BL Lac samples yet derived. Each BL Lac candidate has a high-quality SDSS spectrum from which we determine spectroscopic redshifts for ~60% of the objects. Redshift lower limits are estimated for the remaining objects utilizing the lack of host galaxy flux contamination in their optical spectra; we find that objects lacking spectroscopic redshifts are likely at systematically higher redshifts. Approximately 80% of our BL Lac candidates match to a radio source in FIRST/NVSS, and ~40% match to a ROSAT X-ray source. The homogeneous multiwavelength coverage allows subdivision of the sample into 637 radio-loud BL Lac candidates and 86 weak-featured radio-quiet objects. The radio-loud objects broadly support the standard paradigm unifying BL Lac objects with beamed radio galaxies. We propose that the majority of the radio-quiet objects may be lower-redshift (z<2.2) analogs to high-redshift weak line quasars (i.e., AGN with unusually anemic broad emission line regions). These would constitute the largest sample of such objects, being of similar size and complementary in redshift to the samples of high-redshift weak line quasars previously discovered by the SDSS. However, some fraction of the weak-featured radio-quiet objects may instead populate a rare and extreme radio-weak tail of the much larger radio-loud BL Lac population. Serendipitous discoveries of unusual white dwarfs, high-redshift weak line quasars, and broad absorption line quasars with extreme continuum dropoffs blueward of rest-frame 2800 Angstroms are also briefly described.

We present Spitzer IRS observations of a complete sample of 57 SDSS type-1 quasars at z~1. Aromatic features at 6.2 and/or 7.7 um are detected in about half of the sample and show profiles similar to those seen in normal galaxies at both low- and high-redshift, indicating a star-formation origin for the features. Based on the ratio of aromatic to star-formation IR (SFIR) luminosities for normal star-forming galaxies at z~1, we have constructed the SFIR luminosity function (LF) of z~1 quasars. As we found earlier for low-redshift PG quasars, these z~1 quasars show a flatter SFIR LF than do z~1 field galaxies, implying the quasar host galaxy population has on average a higher SFR than the field galaxies do. As measured from their SFIR LF, individual quasar hosts have on average LIRG-level SFRs, which mainly arise in the circumnuclear regions. By comparing with similar measurements of low-redshift PG quasars, we find that the comoving SFIR luminosity density in quasar hosts shows a much larger increase with redshift than that in field galaxies. The behavior is consistent with pure density evolution since the average SFR and the average SFR/BH-accretion-rate in quasar hosts show little evolution with redshift. For individual quasars, we have found a correlation between the aromatic-based SFR and the luminosity of the nuclear radiation, consistent with predictions of some theoretical models. We propose that type 1 quasars reside in a distinct galaxy population that shows elliptical morphology but that harbors a significant fraction of intermediate-age stars and is experiencing intense circumnuclear star formation.

The [O IV] 25.89 micron line has been shown to be an accurate indicator of active galactic nucleus (AGN) intrinsic luminosity in that it correlates well with hard (10-200 keV) X-ray emission. We present measurements of [O IV] for 89 Seyfert galaxies from the unbiased Revised Shapley-Ames (RSA) sample. The [O IV] luminosity distributions of obscured and unobscured Seyferts are indistinguishable, indicating that their intrinsic AGN luminosities are quite similar and that the RSA sample is well suited for tests of the unified model. In addition, we analyze several commonly used proxies for AGN luminosity, including [O III] 5007 A, 6 cm radio, and 2-10 keV X-ray emission. We find that the radio luminosity distributions of obscured and unobscured AGNs show no significant difference, indicating that radio luminosity is a useful isotropic luminosity indicator. However, the observed [O III] and 2-10 keV luminosities are systematically smaller for obscured Seyferts, indicating that they are not emitted isotropically.

We identify a sample of 74 high-redshift quasars (z>3) with weak emission lines from the Fifth Data Release of the Sloan Digital Sky Survey and present infrared, optical, and radio observations of a subsample of four objects at z>4. These weak emission-line quasars (WLQs) constitute a prominent tail of the Lya+NV equivalent width distribution, and we compare them to quasars with more typical emission-line properties and to low-redshift active galactic nuclei with weak/absent emission lines, namely BL Lac objects. We find that WLQs exhibit hot (T~1000 K) thermal dust emission and have rest-frame 0.1-5 micron spectral energy distributions that are quite similar to those of normal quasars. The variability, polarization, and radio properties of WLQs are also different from those of BL Lacs, making continuum boosting by a relativistic jet an unlikely physical interpretation. The most probable scenario for WLQs involves broad-line region properties that are physically distinct from those of normal quasars.

We compare [O IV] 25.89 micron emission line luminosities with very hard (10-200 keV) X-rays from Swift, Integral, and BeppoSAX for a complete sample of 89 Seyferts from the Revised Shapley-Ames sample. Using Seyfert 1s, we calibrate [O IV] as a measure of AGN intrinsic luminosity, for particular use in high-obscuration environments. With this calibration, we measure the average decrement in 14-195 keV X-ray to [O IV] luminosity ratio for Seyfert 2s compared to type 1s. We find a decrement of 3.1 +- 0.8 for Seyfert 2s, and a decrement of 5.0 +- 2.7 for known Compton-thick Seyfert 2s. These decrements imply column densities of approximately log N(H)=24.6 and 24.7 cm^-2, respectively. Thus, we infer that the average Seyfert 2 is more highly obscured and intrinsically more luminous than would be inferred even from the very hard X-rays. We demonstrate two applications of the hard X-ray to [O IV] ratio. We measure a column density for the extremely obscured NGC 1068 of log N(H)=25.3-25.4 cm^-2. Finally, by comparing [O IV] luminosities to total infrared luminosities for twelve bright ultraluminous infrared galaxies, we find that four have substantial AGN contributions.

We present Chandra observations of nine high-redshift quasars (z=2.7-5.9) discovered by the Sloan Digital Sky Survey with weak or undetectable high-ionization emission lines in their UV spectra (WLQs). Adding archival X-ray observations of six additional sources of this class has enabled us to place the strongest constraints yet on the X-ray properties of this remarkable class of AGNs. Although our data cannot rule out the possibility that the emission lines are overwhelmed by a relativistically boosted continuum, as manifested by BL Lac objects, we find that WLQs are considerably weaker in the X-ray and radio bands than the majority of BL Lacs found at much lower redshifts. If WLQs are high-redshift BL Lacs, then it is difficult to explain the lack of a large parent population of X-ray and radio bright weak-lined sources at high redshift. We also consider the possibility that WLQs are quasars with extreme properties, and in particular that the emission lines are suppressed by high accretion rates. Using joint spectral fitting of the X-ray spectra of 11 WLQs we find that the mean photon index in the hard X-ray band is consistent with those observed in typical radio-quiet AGNs with no hint of an unusually steep hard-X-ray spectrum. This result poses a challenge to the hypothesis that WLQs have extremely high accretion rates, and we discuss additional observations required to test this idea.

[abridged] We present multi-epoch high-resolution optical spectroscopy, UV/radio/X-ray imaging, and archival Hubble and Spitzer observations of an intermediate luminosity optical transient recently discovered in the nearby galaxy NGC300. We find that the transient (NGC300 OT2008-1) has a peak absolute magnitude of M_bol~-11.8 mag, intermediate between novae and supernovae, and similar to the recent events M85 OT2006-1 and SN2008S. Our high-resolution spectra, the first for this event, are dominated by intermediate velocity (~200-1000 km/s) hydrogen Balmer lines and CaII emission and absorption lines that point to a complex circumstellar environment, reminiscent of the yellow hypergiant IRC+10420. In particular, we detect broad CaII H&K absorption with an asymmetric red wing extending to ~1000 km/s, indicative of gas infall onto a massive and relatively compact star (blue supergiant or Wolf-Rayet star); an extended red supergiant progenitor is unlikely. The origin of the inflowing gas may be a previous ejection from the progenitor or the wind of a massive binary companion. The low luminosity, intermediate velocities, and overall similarity to a known eruptive star indicate that the event did not result in a complete disruption of the progenitor. We identify the progenitor in archival Spitzer observations, with deep upper limits from Hubble data. The spectral energy distribution points to a dust-enshrouded star with a luminosity of about 6x10^4 L_sun, indicative of a ~10-20 M_sun progenitor (or binary system). This conclusion is in good agreement with our interpretation of the outburst and circumstellar properties. The lack of significant extinction in the transient spectrum indicates that the dust surrounding the progenitor was cleared by the outburst.

Numerical simulations suggest that active galactic nuclei (AGNs) play an important role in the formation of early-type galaxies by expelling gas and dust in powerful galactic winds and quenching star formation. However, the existence of AGN feedback capable of halting galaxy-wide star formation has yet to be observationally confirmed. To investigate this question, we have obtained spectra of 14 post-starburst galaxies at z~0.6 to search for evidence of galactic winds. In 10/14 galaxies we detect Mg II 2796,2803 absorption lines which are blueshifted by 490 - 2020 km/s with respect to the stars. The median blueshift is 1140 km/s. We hypothesize that the outflowing gas represents a fossil galactic wind launched near the peak of the galaxy's activity, a few 100 Myr ago. The velocities we measure are intermediate between those of luminous starbursts and broad absorption line quasars, which suggests that feedback from an AGN may have played a role in expelling cool gas and shutting down star formation.

Observations of candidate BL Lacertae objects spectroscopically selected from the Sloan Digital Sky Survey (SDSS) reveal a large fraction with high polarization (P > 3%). This result confirms that synchrotron radiation makes an important contribution to the observed optical continuum for most objects in the sample. The SDSS sample can be divided into separate categories, with objects of undetermined redshift generally having the highest optical polarization. Polarization as high as 23% and the lack of spectral features suggests that the synchrotron continuum completely dominates the spectra of these sources. The mean polarization levels observed for objects having measured redshifts is much lower, with the maximum observed polarization for this group being ~10%. The lower polarizations of these objects are reminiscent of the less spectacular polarization levels shown by BL Lac objects discovered in X-ray surveys. We find no SDSS BL Lac candidates at z > 1 with P > 3%, calling their classification as BL Lac objects into question. In addition, the existence of radio-quiet BL Lac objects is not verified since none of 10 potentially radio-weak BL Lac candidates observed are highly polarized. Regardless of whether the high-redshift and radio-weak objects are included in this optical sample, the overall levels of polarization observed are intermediate between those seen for X-ray and radio-selected BL Lac objects.

We present early photometric and spectroscopic data on the afterglow of GRB 060218 and report the evolution of the underlying supernova 2006aj. Our data span a time-range of 4 days to 22 days after the GRB and clearly establish that SN 2006aj is a fast-evolving broad-lined Type Ic SN with an extremely short rise-time (~ 10 days) and a large optical luminosity (M_V = -18.7 mag). The SN properties are deduced well since the GRB afterglow does not contribute a significant amount to the total light output. The spectra show broad lines indicative of large expansion velocities, but are better matched by those of SN 2002ap and SN 1997ef (that are not associated with a GRB) than those of the proto-typical GRB-related SN 1998bw. We refine the redshift estimate to z = 0.0335 +/- 0.00007. The host-galaxy is a low-metallicity dwarf galaxy (with M_V \~ -16.0 mag), similar to host-galaxies of other GRB-associated SNe.

We present new Chandra observations of 21 z>4 quasars, including 11 sources at z>5. These observations double the number of X-ray detected quasars at z>5, allowing investigation of the X-ray spectral properties of a substantial sample of quasars at the dawn of the modern Universe. By jointly fitting the spectra of 15 z>5 radio-quiet quasars (RQQs), including sources from the Chandra archive, with a total of 185 photons, we find a mean X-ray power-law photon index of Gamma=1.95^{+0.30}_{-0.26}, and a mean neutral intrinsic absorption column density of N_H<~6x10^{22} cm^{-2}. These results show that quasar X-ray spectral properties have not evolved up to the highest observable redshifts. We also find that the mean optical-X-ray spectral slope (alpha_ox) of optically-selected z>5 RQQs, excluding broad absorption line quasars, is alpha_ox=-1.69+/-0.03, which is consistent with the value predicted from the observed relationship between alpha_ox and ultraviolet luminosity. Four of the sources in our sample are members of the rare class of weak emission-line quasars, and we detect two of them in X-rays. We discuss the implications our X-ray observations have for the nature of these mysterious sources and, in particular, whether their weak-line spectra are a consequence of continuum boosting or a deficit of high-ionization line emitting gas.

We present the discovery of seven quasars at z>5.7, selected from ~2000 deg^2 of multicolor imaging data of the Sloan Digital Sky Survey (SDSS). The new quasars have redshifts z from 5.79 to 6.13. Five are selected as part of a complete flux-limited sample in the SDSS Northern Galactic Cap; two have larger photometric errors and are not part of the complete sample. One of the new quasars, SDSS J1335+3533 (z=5.93), exhibits no emission lines; the 3-sigma limit on the rest-frame equivalent width of Ly alpha+NV line is 5 A. It is the highest redshift lineless quasar known, and could be a gravitational lensed galaxy, a BL Lac object or a new type of quasar. Two new z>6 quasars, SDSS 1250+3130 (z=6.13) and SDSS J1137+3549 (z=6.01), show deep Gunn-Peterson absorption gaps in Ly alpha. These gaps are narrower the complete Gunn-Peterson absorption troughs observed among quasars at z>6.2 and do not have complete Ly beta absorption.