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

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

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.

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

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

Well-calibrated spectropolarimetry studies at resolutions of $R>$10,000 with signal-to-noise ratios (SNRs) better than 0.01\% across individual line profiles, are becoming common with larger aperture telescopes. Spectropolarimetric studies require high SNR observations and are often limited by instrument systematic errors. As an example, fiber-fed spectropolarimeters combined with advanced line-combination algorithms can reach statistical error limits of 0.001\% in measurements of spectral line profiles referenced to the continuum. Calibration of such observations is often required both for cross-talk and for continuum polarization. This is not straightforward since telescope cross-talk errors are rarely less than $\sim$1\%. In solar instruments like the Daniel K. Inouye Solar Telescope (DKIST), much more stringent calibration is required and the telescope optical design contains substantial intrinsic polarization artifacts. This paper describes some generally useful techniques we have applied to the HiVIS spectropolarimeter at the 3.7m AEOS telescope on Haleakala. HiVIS now yields accurate polarized spectral line profiles that are shot-noise limited to 0.01\% SNR levels at our full spectral resolution of 10,000 at spectral sampling of $\sim$100,000. We show line profiles with absolute spectropolarimetric calibration for cross-talk and continuum polarization in a system with polarization cross-talk levels of essentially 100\%. In these data the continuum polarization can be recovered to one percent accuracy because of synchronized charge-shuffling model now working with our CCD detector. These techniques can be applied to other spectropolarimeters on other telescopes for both night and day-time applications such as DKIST, TMT and ELT which have folded non-axially symmetric foci.