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

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.

We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S.

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 present our extensive observational campaign on the Swift-discovered GRB141121A, al- most ten years after its launch. Our observations covers radio through X-rays, and extends for more than 30 days after discovery. The prompt phase of GRB 141121A lasted 1410 s and, at the derived redshift of z = 1.469, the isotropic energy is E{\gamma},iso = 8.0x10^52 erg. Due to the long prompt duration, GRB141121A falls into the recently discovered class of UL-GRBs. Peculiar features of this burst are a flat early-time optical light curve and a radio-to-X-ray rebrightening around 3 days after the burst. The latter is followed by a steep optical-to-X-ray decay and a much shallower radio fading. We analyze GRB 141121A in the context of the standard forward-reverse shock (FS,RS) scenario and we disentangle the FS and RS contributions. Finally, we comment on the puzzling early-time (t ~3 d) behavior of GRB 141121A, and suggest that its interpretation may require a two-component jet model. Overall, our analysis confirms that the class of UL-GRBs represents our best opportunity to firmly establish the prominent emission mechanisms in action during powerful GRB explosions, and future missions (like SVOM, XTiDE, or ISS-Lobster) will provide many more of such objects.

We present optical and near-infrared (NIR) photometry for three gamma-ray burst supernovae (GRB-SNe): GRB 120729A, GRB 130215A / SN 2013ez and GRB 130831A / SN 2013fu. In the case of GRB 130215A / SN 2013ez, we also present optical spectroscopy at t-t0=16.1 d, which covers rest-frame 3000-6250 Angstroms. Based on Fe II (5169) and Si (II) (6355), our spectrum indicates an unusually low expansion velocity of 4000-6350 km/s, the lowest ever measured for a GRB-SN. Additionally, we determined the brightness and shape of each accompanying SN relative to a template supernova (SN 1998bw), which were used to estimate the amount of nickel produced via nucleosynthesis during each explosion. We find that our derived nickel masses are typical of other GRB-SNe, and greater than those of SNe Ibc that are not associated with GRBs. For GRB 130831A / SN 2013fu, we use our well-sampled R-band light curve (LC) to estimate the amount of ejecta mass and the kinetic energy of the SN, finding that these too are similar to other GRB-SNe. For GRB 130215A, we take advantage of contemporaneous optical/NIR observations to construct an optical/NIR bolometric LC of the afterglow. We fit the bolometric LC with the millisecond magnetar model of Zhang & Meszaros (2001), which considers dipole radiation as a source of energy injection to the forward shock powering the optical/NIR afterglow. Using this model we derive an initial spin period of P=12 ms and a magnetic field of B=1.1 x 10^15 G, which are commensurate with those found for proposed magnetar central engines of other long-duration GRBs.

We investigate Schmidt's conjecture (i.e., that the star formation rate scales in a power-law fashion with the gas density) for four well-studied local molecular clouds (GMCs). Using the Bayesian methodology we show that a local Schmidt scaling relation of the form Sigma*(A_K) = kappa x (A_K)^{beta} (protostars pc^{-2}) exists within (but not between) GMCs. Further we find that the Schmidt scaling law, by itself, does not provide an adequate description of star formation activity in GMCs. Because the total number of protostars produced by a cloud is given by the product of Sigma*(A_K) and S'(> A_K), the differential surface area distribution function, integrated over the entire cloud, the cloud's structure plays a fundamental role in setting the level of its star formation activity. For clouds with similar functional forms of Sigma*(A_K), observed differences in their total SFRs are primarily due to the differences in S'(> A_K) between the clouds. The coupling of Sigma*(A_K) with the measured S'(> A_K) in these clouds also produces a steep jump in the SFR and protostellar production above A_K ~ 0.8 magnitudes. Finally, we show that there is no global Schmidt law that relates the star formation rate and gas mass surface densities between GMCs. Consequently, the observed Kennicutt-Schmidt scaling relation for disk galaxies is likely an artifact of unresolved measurements of GMCs and not a result of any underlying physical law of star formation characterizing the molecular gas.

Despite over 30 years of study, the mass-area relationship within and among clouds is still poorly understood both observationally and theoretically. Modern extinction datasets should have sufficient resolution and dynamic range to characterize this relationship for nearby molecular clouds, although recent papers using extinction data seem to yield different interpretations regarding the nature and universality of this aspect of cloud structure. In this paper we try to unify these various results and interpretations by accounting for the different ways cloud properties are measured and analyzed. We interpret the mass-area relationship in terms of the column density distribution function and its possible variation within and among clouds. We quantitatively characterize regional variations in the column density PDF. We show that structures both within and among clouds possess the same degree of "universality", in that their PDF means do not systematically scale with structure size. Because of this, mass scales linearly with area.

We present the underlying relations between colour-magnitude diagrams (CMDs) and synthesis models through the use of stellar luminosity distribution func- tions. CMDs studies make a direct use of the stellar luminosity distribution function while, in general, synthesis models only use its mean value, even though high-order moments can also be obtained. We show that the mean, high-order moments and in- tegrated luminosity distribution functions of stellar ensembles are related to the stellar luminosity distribution function, within the formalism of probabilistic synthesis mod- els. More details have been yet presented in Cervin ~ o & Luridiana (2006) and references therein. As a direct application of this formalism, we discuss two key issues. First, in- ferences on the upper mass limit of the initial mass function as a function of the total mass of clusters. Second, we apply extreme value theory to show that that the cluster mass obtained from normalising the IMF between mmax and mup does not provide the cluster mass in the case where only one star in this mass range is present, as assumed in the IGIMF theory. It provides instead the cluster mass with a 60% probability to have a star with mass larger than mmax, and we argue that in light of this result the basic formulation ofthe IGIMF theory must be revised.

We have used moderate resolution, near-infrared spectra from the SpeX spectrograph on the NASA Infrared Telescope facility to characterize the stellar content of Barnard 59 (B59), the most active star-forming core in the Pipe Nebula. Measuring luminosity and temperature sensitive features in the spectra of 20 candidate YSOs, we identified likely background giant stars and measured each star's spectral type, extinction, and NIR continuum excess. We find that B59 is composed of late type (K4-M6) low-mass (0.9--0.1 M_sun) YSOs whose median stellar age is comparable to, if not slightly older than, that of YSOs within the Rho Oph, Taurus, and Chameleon star forming regions. Deriving absolute age estimates from pre-main sequence models computed by D'Antona et al., and accounting only for statistical uncertainties, we measure B59's median stellar age to be 2.6+/-0.8 Myrs. Including potential systematic effects increases the error budget for B59's median (DM98) stellar age to 2.6+4.1/-2.6 Myrs. We also find that the relative age orderings implied by pre-main sequence evolutionary tracks depend on the range of stellar masses sampled, as model isochrones possess significantly different mass dependencies. The maximum likelihood median stellar age we measure for B59, and the region's observed gas properties, suggest that the B59 dense core has been stable against global collapse for roughly 6 dynamical timescales, and is actively forming stars with a star formation efficiency per dynamical time of ~6%. This maximum likelihood value agrees well with recent star formation simulations that incorporate various forms of support against collapse, such as sub-critical magnetic fields, outflows, and radiative feedback from protostellar heating. [abridged]

We present arcsecond resolution mid-infrared and millimeter observations of the center of the young stellar cluster AFGL961 in the Rosette molecular cloud. Within 0.2 pc of each other, we find an early B star embedded in a dense core, a neighboring star of similar luminosity with no millimeter counterpart, a protostar that has cleared out a cavity in the circumcluster envelope, and two massive, dense cores with no infrared counterparts. An outflow emanates from one of these cores, indicating a deeply embedded protostar, but the other is starless, bound, and appears to be collapsing. The diversity of states implies either that protostellar evolution is faster in clusters than in isolation or that clusters form via quasi-static rather than dynamic collapse. The existence of a pre-stellar core at the cluster center shows that that some star formation continues after and in close proximity to massive, ionizing stars.

We explore here the young stellar populations in the Rosette Molecular Cloud (RMC) region with high spatial resolution X-ray images from the Chandra X-ray Observatory, which are effective in locating weak-lined T Tauri stars as well as disk-bearing young stars. A total of 395 X-ray point sources are detected, 299 of which (76%) have an optical or near-infrared (NIR) counterpart identified from deep FLAMINGOS images. From X-ray and mass sensitivity limits, we infer a total population of about 1700 young stars in the survey region. Based on smoothed stellar surface density maps, we investigate the spatial distribution of the X-ray sources and define three distinctive structures and substructures within them. Structures B and C are associated with previously known embedded IR clusters, while structure A is a new X-ray-identified unobscured cluster. A high mass protostar RMCX #89 = IRAS 06306+0437 and its associated sparse cluster is studied. The different subregions are not coeval but do not show a simple spatial-age pattern. Disk fractions vary between subregions and are generally 20% of the total stellar population inferred from the X-ray survey. The data are consistent with speculations that triggered star formation around the HII region is present in the RMC, but do not support a simple sequential triggering process through the cloud interior. While a significant fraction of young stars are located in a distributed population throughout the RMC region, it is not clear they originated in clustered environments.

We present the first high spatial resolution Chandra X-ray study of NGC 2244, the 2 Myr old stellar cluster immersed in the Rosette Nebula. Over 900 X-ray sources are detected; 77% have optical or FLAMINGOS near-infrared (NIR) stellar counterparts and are mostly previously uncatalogued young stellar cluster members. All known OB stars with spectral type earlier than B1 are detected and the X-ray selected stellar population is estimated to be nearly complete between 0.5 and 3 Msun. The X-ray luminosity function (XLF) ranges from 29.4<logLx<32 ergs/s in the hard (2-8keV) band. By comparing the NGC 2244 and Orion Nebula Cluster XLFs, we estimate a total population of 2000 stars in NGC 2244. A number of further results emerge from our analysis: The XLF and the associated K-band luminosity function indicate a normal Salpeter initial mass function (IMF) for NGC 2244. This is inconsistent with the top-heavy IMF reported from earlier optical studies that lacked a good census of <4Msun stars. The spatial distribution of X-ray stars is strongly concentrated around the central O5 star, HD 46150. The other early O star, HD 46223, has few companions. The cluster's stellar radial density profile shows two distinctive structures. This double structure, combined with the absence of mass segregation, indicates that this cluster is not in dynamical equilibrium. The spatial distribution of X-ray selected K-excess disk stars and embedded stars is asymmetric with an apparent deficit towards the north. The fraction of X-ray-selected cluster members with K-band excesses caused by inner protoplanetary disks is 6%, slightly lower than the 10% disk fraction estimated from the FLAMINGOS study based on the NIR-selected sample. This is due to the high efficiency of X-ray surveys in locating disk-free T Tauri stars.[Abridged]

We present the results of a deep near-infrared imaging survey of the Rosette Complex. We studied the distribution of young embedded sources using a variation of the Nearest Neighbor Method applied to a carefully selected sample of near-infrared excess (NIRX) stars which trace the latest episode of star formation in the complex. Our analysis confirmed the existence of seven clusters previously detected in the molecular cloud, and identified four more clusters across the complex. We determined that 60% of the young stars in the complex and 86% of the stars within the molecular cloud are contained in clusters, implying that the majority of stars in the Rosette formed in embedded clusters. We compare the sizes, infrared excess fractions and average extinction towards individual clusters to investigate their early evolution and expansion. We found that the average infrared excess fraction of clusters increases as a function of distance from NGC 2244, implying a temporal sequence of star formation across the complex. This sequence appears to be primordial, possibly resulting from the formation and evolution of the molecular cloud and not from the interaction with the HII region.

We present wide-field near-infrared (JHK) images of the young, 2 Myr IC 348 cluster taken with FLAMINGOS. We use these new data to construct an infrared census of sources, which is sensitive enough to detect a 10 Mjup brown dwarf seen through an extinction of Av=7mag. We examine the cluster's structure and relationship to the molecular cloud and construct the cluster's K band luminosity function. Using our model luminosity function algorithm, we derive the cluster's initial mass function throughout the stellar and substellar regimes and find that the IC 348 IMF is very similar to that found for the Trapezium Cluster with both cluster IMFs having a mode between 0.2 - 0.08 Msun. In particular we find that, similar to our results for the Trapezium, brown dwarfs constitute only 1 in 4 of the sources in the IC 348 cluster. We show that a modest secondary peak forms in the substellar IC 348 KLF, corresponding to the same mass range responsible for a similar KLF peak found in the Trapezium. We interpret this KLF peak as either evidence for a corresponding secondary IMF peak at the deuterium burning limit, or as arising from a feature in the substellar mass-luminosity relation that is not predicted by current theoretical models. Lastly, we find that IC 348 displays radial variations of its sub-solar (0.5 - 0.08 Msun) IMF on a parsec scale. Whatever mechanism that is breaking the universality of the IMF on small spatial scales in IC 348 does not appear to be acting upon the brown dwarf population, whose relative size does not vary with distance from the cluster center.