The spectral analysis and data products in Data Release 16 (DR16; December 2019) from the high-resolution near-infrared APOGEE-2/SDSS-IV survey are described. Compared to the previous APOGEE data release (DR14; July 2017), APOGEE DR16 includes about 200000 new stellar spectra, of which 100000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observatory in Chile. DR16 includes all data taken up to August 2018, including data released in previous data releases. All of the data have been re-reduced and re-analyzed using the latest pipelines, resulting in a total of 473307 spectra of 437445 stars. Changes to the analysis methods for this release include, but are not limited to, the use of MARCS model atmospheres for calculation of the entire main grid of synthetic spectra used in the analysis, a new method for filling "holes" in the grids due to unconverged model atmospheres, and a new scheme for continuum normalization. Abundances of the neutron capture element Ce are included for the first time. A new scheme for estimating uncertainties of the derived quantities using stars with multiple observations has been applied, and calibrated values of surface gravities for dwarf stars are now supplied. Compared to DR14, the radial velocities derived for this release more closely match those in the Gaia DR2 data base, and a clear improvement in the spectral analysis of the coolest giants can be seen. The reduced spectra as well as the result of the analysis can be downloaded using links provided in the SDSS DR16 web page.

We report the stellar atmospheric parameters for 7503 spectra contained in the first release of the MaNGA stellar library (MaStar) in SDSS DR15. The first release of MaStar contains 8646 spectra measured from 3321 unique stars, each covering the wavelength range 3622 \AA\ to 10354 \AA\ with a resolving power of $R \sim$ 1800. In this work, we first determined the basic stellar parameters: effective temperature ($\rm T_{eff}$), surface gravity ($\log g$), and metallicity ($\rm[Fe/H]$), which best fit the data using an empirical interpolator based on the Medium-resolution Isaac Newton Telescope library of empirical spectra (MILES), as implemented by the University of Lyon Spectroscopic analysis Software (Koleva et al. 2008, ULySS) package. While we analyzed all 8646 spectra from the first release of MaStar, since MaStar has a wider parameter-space coverage than MILES, not all of these fits are robust. In addition, not all parameter regions covered by MILES yield robust results, likely due to the non-uniform coverage of the parameter space by MILES. We tested the robustness of the method using the MILES spectra itself and identified a proxy based on the local density of the training set. With this proxy, we identified 7503 MaStar spectra with robust fitting results. They cover the range from 3179K to 20,517K in effective temperature ($\rm T_{eff}$), from 0.40 to 5.0 in surface gravity ($\log g$), and from $-$2.49 to $+$0.73 in metallicity ($\rm[Fe/H]$).

The Large Magellanic Cloud (LMC) is the closest and most studied example of an irregular galaxy. Among its principal defining morphological features, its off-centred bar and single spiral arm stand out, defining a whole family of galaxies known as the Magellanic spirals (Sm). These structures are thought to be triggered by tidal interactions and possibly maintained via gas accretion. However, it is still unknown whether they are long-lived stable structures. In this work, by combining photometry that reaches down to the oldest main sequence turn-off in the colour-magnitude diagrams (CMD, up to a distance of $\sim$4.4 kpc from the LMC centre) from the SMASH survey and CMD fitting techniques, we find compelling evidence supporting the long-term stability of the LMC spiral arm, dating the origin of this structure to more than 2~Gyr ago. The evidence suggests that the close encounter between the LMC and the Small Magellanic Cloud (SMC) that produced the gaseous Magellanic Stream and its Leading Arm (LA) also triggered the formation of the LMC's spiral arm. Given the mass difference between the Clouds and the notable consequences of this interaction, we can speculate that this should have been one of their closest encounters. These results set important constraints on the timing of LMC-SMC collisions, as well as on the physics behind star formation induced by tidal encounters.

We use the first release of the SDSS/MaStar stellar library comprising ~9000, high S/N spectra, to calculate integrated spectra of stellar population models. The models extend over the wavelength range 0.36-1.03 micron and share the same spectral resolution (R~1800) and flux calibration as the SDSS-IV/MaNGA galaxy data. The parameter space covered by the stellar spectra collected thus far allows the calculation of models with ages and chemical composition in the range t>200 Myr, -2 <=[Z/H]<= + 0.35, which will be extended as MaStar proceeds. Notably, the models include spectra for dwarf Main Sequence stars close to the core H-burning limit, as well as spectra for cold, metal-rich giants. Both stellar types are crucial for modelling lambda>0.7 micron absorption spectra. Moreover, a better parameter coverage at low metallicity allows the calculation of models as young as 500 Myr and the full account of the Blue Horizontal Branch phase of old populations. We present models adopting two independent sets of stellar parameters (T_eff, logg, [Z/H]). In a novel approach, their reliability is tested 'on the fly' using the stellar population models themselves. We perform tests with Milky Way and Magellanic Clouds globular clusters, finding that the new models recover their ages and metallicities remarkably well, with systematics as low as a few per cent for homogeneous calibration sets. We also fit a MaNGA galaxy spectrum, finding residuals of the order of a few per cent comparable to the state-of-art models, but now over a wider wavelength range.

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 investigate the Fe, C, N, O, Mg, Al, Si, K, Ca, Ce and Nd abundances of 2283 red giant stars in 31 globular clusters from high-resolution spectra observed in both the northern and southern hemisphere by the SDSS-IV APOGEE-2 survey. This unprecedented homogeneous dataset, largest to date, allows us to discuss the intrinsic Fe spread, the shape and statistics of Al-Mg and N-C anticorrelations as a function of cluster mass, luminosity, age and metallicity for all 31 clusters. We find that the Fe spread does not depend on these parameters within our uncertainties including cluster metallicity, contradicting earlier observations. We do not confirm the metallicity variations previously observed in M22 and NGC 1851. Some clusters show a bimodal Al distribution, while others exhibit a continuous distribution as has been previously reported in the literature. We confirm more than 2 populations in $\omega$ Cen and NGC 6752, and find new ones in M79. We discuss the scatter of Al by implementing a correction to the standard chemical evolution of Al in the Milky Way. After correction, its dependence on cluster mass is increased suggesting that the extent of Al enrichment as a function of mass was suppressed before the correction. We observe a turnover in the Mg-Al anticorrelation at very low Mg in $\omega$ Cen, similar to the pattern previously reported in M15 and M92. $\omega$ Cen may also have a weak K-Mg anticorrelation, and if confirmed, it would be only the third cluster known to show such a pattern.

We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

We combine high-resolution spectroscopic data from APOGEE-2 Survey Data Release 16 (DR16) with broad-band photometric data from several sources, as well as parallaxes from {\it Gaia} Data Release 2 (DR2). Using the Bayesian isochrone-fitting code {\tt StarHorse}, we derive distances, extinctions and astrophysical parameters for around 388,815 APOGEE stars, achieving typical distance uncertainties of $\sim 6\%$ for APOGEE giants, $\sim 2\%$ for APOGEE dwarfs, as well as extinction uncertainties of $\sim 0.07$ mag when all photometric information is available, and $\sim 0.17$ mag if optical photometry is missing. {\tt StarHorse} uncertainties vary with the input spectroscopic catalogue, with the available photometry, and with the parallax uncertainties. To illustrate the impact of our results, we show that, thanks to {\it Gaia} DR2 and the now larger sky coverage of APOGEE-2 (including APOGEE-South), we obtain an extended map of the Galactic plane, providing unprecedented coverage of the disk close to the Galactic mid-plane ($|Z_{Gal}|<1$ kpc) from the Galactic Centre out to $R_{\rm Gal}\sim 20$ kpc. The improvements in statistics as well as distance and extinction uncertainties unveil the presence of the bar in stellar density, as well as the striking chemical duality in the innermost regions of the disk, now clearly extending to the inner bulge. We complement this paper with distances and extinctions for stars in other public released spectroscopic surveys: 324,999 in GALAH DR2, 4,928,715 in LAMOST DR5, 408,894 in RAVE DR6, and 6,095 in GES DR3

We present a semi-empirical spectral classification scheme for normal B-type stars using near-infrared spectra (1.5-1.7 $\mu$m) from the SDSS APOGEE2-N DR14 database. The main motivation for working with B-type stars is their importance in the evolution of young stellar clusters, however we also take advantage of having a numerous sample (316 stars) of B-type star candidates in APOGEE2-N, for which we also have optical (3600-9100 \angstrom) counterparts from the LAMOST survey. By first obtaining an accurate spectral classification of the sources using the LAMOST DR3 spectra and the canonical spectral classification scheme Gray & Corbally 2009, we found a linear relation between optical spectral types and the equivalent widths of the hydrogen lines of the Brackett series in the APOGEE2-N NIR spectra. This relation extends smoothly from a similar relation for O and early-B stars found by Roman-Lopes et al. (2018). This way, we obtain a catalog of B-type sources with features in both the optical and NIR, and a classification scheme refined down to one spectral sub-class.

Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic "discovery" threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches.

We report the first APOGEE metallicities and alpha-element abundances measured for 3600 red giant stars spanning a large radial range of both the Large (LMC) and Small Magellanic Clouds (SMC), the largest Milky Way dwarf galaxies. Our sample is an order of magnitude larger than that of previous studies, and extends to much larger radial distances. These are the first results presented that make use of the newly installed Southern APOGEE instrument on the du Pont telescope at Las Campanas Observatory. Our unbiased sample of the LMC spans a large range in metallicity, from [Fe/H]=-0.2 to very metal-poor stars with [Fe/H]=-2.5, the most metal-poor Magellanic Clouds (MCs) stars detected to date. The LMC [alpha/Fe]-[Fe/H] distribution is very flat over a large metallicity range, but rises by ~0.1 dex at -1.0<[Fe/H]<-0.5. We interpret this as a sign of the known recent increase in MC star-formation activity, and are able to reproduce the pattern with a chemical evolution model that includes a recent "starburst". At the metal-poor end, we capture the increase of [alpha/Fe] with decreasing [Fe/H], and constrain the "alpha-knee" to [Fe/H]<-2.2 in both MCs, implying a low star-formation efficiency of ~0.01 Gyr^-1. The MC knees are more metal poor than those of less massive Milky Way (MW) dwarf galaxies such as Fornax, Sculptor, or Sagittarius. One possible interpretation is that the MCs formed in a lower-density environment than the MW, a hypothesis that is consistent with the paradigm that the MCs fell into the MW's gravitational potential only recently.

Using 3D positions and kinematics of stars relative to the Sagittarius (Sgr) orbital plane and angular momentum, we identify 166 Sgr stream members observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) that also have Gaia DR2 astrometry. This sample of 63/103 stars in the Sgr trailing/leading arm are combined with an APOGEE sample of 710 members of the Sgr dwarf spheroidal core (385 of them newly presented here) to establish differences of 0.6 dex in median metallicity and 0.1 dex in [$\alpha$/Fe] between our Sgr core and dynamically older stream samples. Mild chemical gradients are found internally along each arm, but these steepen when anchored by core stars. With a model of Sgr tidal disruption providing estimated dynamical ages (i.e., stripping times) for each stream star, we find a mean metallicity gradient of 0.12 +/- 0.03 dex/Gyr for stars stripped from Sgr over time. For the first time, an [$\alpha$/Fe] gradient is also measured within the stream, at 0.02 +/- 0.01 dex/Gyr using magnesium abundances and 0.04 +/- 0.01 dex/Gyr using silicon, which imply that the Sgr progenitor had significant radial abundance gradients. We discuss the magnitude of those inferred gradients and their implication for the nature of the Sgr progenitor within the context of the current family of Milky Way satellite galaxies, and suggest that more sophisticated Sgr models are needed to properly interpret the growing chemodynamical detail we have on the Sgr system.

We present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range of 3,622-10,354A at a resolving power of R~1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the SDSS-IV/APOGEE-2N observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [alpha/Fe], achieved by a sophisticated target selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar population synthesis, and is particularly well-suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of Sloan Digital Sky Survey Data Release 15. We expect the final release of the library to contain more than 10,000 stars.

Understanding the evolutionary history of the Magellanic Clouds requires an in-depth exploration and characterization of the stellar content in their outer regions, which ultimately are key to tracing the epochs and nature of past interactions. We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at ~1.9 degr from the center of the SMC in the north-east direction.The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2. The shell is composed of a young population with an age ~ 150 Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute g-band magnitude of the shell is M_{g,shell} = -10.78+/- 0.02, with a surface brightness of mu_{g,shell} = 25.81+/- 0.01 mag/arcsec^{2}. We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, ~150 Myr ago.

Whether supernovae are a significant source of dust has been a long-standing debate. The large quantities of dust observed in high-redshift galaxies raise a fundamental question as to the origin of dust in the Universe since stars cannot have evolved to the AGB dust-producing phase in high-redshift galaxies. In contrast, supernovae occur within several millions of years after the onset of star formation. This white paper focuses on dust formation in supernova ejecta with US-Extremely Large Telescope (ELT) perspective during the era of JWST and LSST.

Supernovae are among the most powerful and influential explosions in the universe. They are also ideal multi-messenger laboratories to study extreme astrophysics. However, many fundamental properties of supernovae related to their diverse progenitor systems and explosion mechanisms remain poorly constrained. Here we outline how late-time spectroscopic observations obtained during the nebular phase (several months to years after explosion), made possible with the next generation of Extremely Large Telescopes, will facilitate transformational science opportunities and rapidly accelerate the community towards our goal of achieving a complete understanding of supernova explosions. We highlight specific examples of how complementary GMT and TMT instrumentation will enable high fidelity spectroscopy from which the line profiles and luminosities of elements tracing mass loss and ejecta can be used to extract kinematic and chemical information with unprecedented detail, for hundreds of objects. This will provide uniquely powerful constraints on the evolutionary phases stars may experience approaching a supernova explosion; the subsequent explosion dynamics; their nucleosynthesis yields; and the formation of compact objects that may act as central engines.

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $\Sigma_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

We explore the stellar structure of the Large Magellanic Cloud (LMC) disk using data from the Survey of the MAgellanic Stellar History (SMASH) and the Dark Energy Survey. We detect a ring-like stellar overdensity in the red clump star count map at a radius of ~6 degrees (~5.2 kpc at the LMC distance) that is continuous over ~270 degrees in position angle and is only limited by the current data coverage. The overdensity shows an amplitude up to 2.5 times higher than that of the underlying smooth disk. This structure might be related to the multiple arms found by de Vaucouleurs. We find that the overdensity shows spatial correlation with intermediate-age star clusters, but not with young (< 1 Gyr) main-sequence stars, indicating the stellar populations associated with the overdensity are intermediate in age or older. Our findings on the LMC overdensity can be explained by either of two distinct formation mechanisms of a ring-like overdensity: (1) the overdensity formed out of an asymmetric one-armed spiral wrapping around the LMC main body, which is induced by repeated encounters with the Small Magellanic Cloud (SMC) over the last Gyr, or (2) the overdensity formed very recently as a tidal response to a direct collision with the SMC. Although the measured properties of the overdensity alone cannot distinguish between the two candidate scenarios, the consistency with both scenarios suggests that the ring-like overdensity is likely a product of tidal interaction with the SMC, but not with the Milky Way halo.

We present a study of the three-dimensional (3D) structure of the Large Magellanic Cloud (LMC) using ~2.2 million red clump (RC) stars selected from the Survey of the MAgellanic Stellar History. To correct for line-of-sight dust extinction, the intrinsic RC color and magnitude and their radial dependence are carefully measured by using internal nearly dust-free regions. These are then used to construct an accurate 2D reddening map (165 square degrees with ~10 arcmin resolution) of the LMC disk and the 3D spatial distribution of RC stars. An inclined disk model is fit to the 2D distance map yielding a best-fit inclination angle i = 25.86(+0.73,-1.39) degrees with random errors of +\-0.19 degrees and line-of-nodes position angle theta = 149.23(+6.43,-8.35) degrees with random errors of +/-0.49 degrees. These angles vary with galactic radius, indicating that the LMC disk is warped and twisted likely due to the repeated tidal interactions with the Small Magellanic Cloud (SMC). For the first time, our data reveal a significant warp in the southwestern part of the outer disk starting at rho ~ 7 degrees that departs from the defined LMC plane up to ~4 kpc toward the SMC, suggesting that it originated from a strong interaction with the SMC. In addition, the inner disk encompassing the off-centered bar appears to be tilted up to 5-15 degrees relative to the rest of the LMC disk. These findings on the outer warp and the tilted bar are consistent with the predictions from the Besla et al. simulation of a recent direct collision with the SMC.

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.

Data and analysis methodology used for the SDSS/APOGEE Data Releases 13 and 14 are described, highlighting differences from the DR12 analysis presented in Holtzman (2015). Some improvement in the handling of telluric absorption and persistence is demonstrated. The derivation and calibration of stellar parameters, chemical abundances, and respective uncertainties are described, along with the ranges over which calibration was performed. Some known issues with the public data related to the calibration of the effective temperatures (DR13), surface gravity (DR13 and DR14), and C and N abundances for dwarfs (DR13 and DR14) are highlighted. We discuss how results from a data-driven technique, The Cannon (Casey 2016), are included in DR14, and compare those with results from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). We describe how using The Cannon in a mode that restricts the abundance analysis of each element to regions of the spectrum with known features from that element leads to Cannon abundances can lead to significantly different results for some elements than when all regions of the spectrum are used to derive abundances.

Data from the SDSS-IV / Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) have been released as part of SDSS Data Releases 13 (DR13) and 14 (DR14). These include high resolution H-band spectra, radial velocities, and derived stellar parameters and abundances. DR13, released in August 2016, contained APOGEE data for roughly 150,000 stars, and DR14, released in August 2017, added about 110,000 more. Stellar parameters and abundances have been derived with an automated pipeline, the APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP). We evaluate the performance of this pipeline by comparing the derived stellar parameters and abundances to those inferred from optical spectra and analysis for several hundred stars. For most elements -- C, Na, Mg, Al, Si, S, Ca, Cr, Mn, Ni -- the DR14 ASPCAP analysis have systematic differences with the comparisons samples of less than 0.05 dex (median), and random differences of less than 0.15 dex (standard deviation). These differences are a combination of the uncertainties in both the comparison samples as well as the ASPCAP-analysis. Compared to the references, magnesium is the most accurate alpha-element derived by ASPCAP, and shows a very clear thin/thick disk separation, while nickel is the most accurate iron-peak element (besides iron).

We undertake a critical and comprehensive comparison of the radial velocities and the main stellar atmospheric parameters for stars in common between the latest data releases from the APOGEE and the LAMOST surveys. There is a total of 42,420 dwarfs/giants stars in common between the APOGEE DR14 - LAMOST DR3 stellar catalogs. A comparison between the RVs shows an offset of 4.54 $\pm$ 0.03 km/s, with a dispersion of 5.8 km/s, in the sense that APOGEE RVs are larger. We observe a small offset in the Teff of about 13 K, with a scatter of 155 K. Small offset in [Fe/H] of about 0.06 dex together with a scatter of 0.13 dex is also observed. We notice that the largest offset between the surveys occurs in the surface gravities. Using only surface gravities in calibrated red giants from APOGEE DR14, where there are 24,074 stars in common, a deviation of 0.14 dex is found with substantial scatter. There are 17,482 red giant stars in common between APOGEE DR14 and those in LAMOST tied to APOGEE DR12 via the Cannon. There is generally good agreement between the two data-sets. However, we find dependencies of the differences of the stellar parameters on effective temperature. For metal-rich stars, a different trend for the [Fe/H] discrepancies is found. Surprisingly, we see no correlation between the internal APOGEE DR14 - DR12 differences in Teff and those in DR14 - LAMOST tied to DR12, where a correlation should be expected since LAMOST has been calibrated to APOGEE DR12. We also find no correlation between the [Fe/H] discrepancies, suggesting that LAMOST/Cannon is not well coupled to the APOGEE DR12 stellar parameters scale. A [Fe/H] dependence between the stellar parameters in APOGEE DR12 and those in DR14 is reported. We find a weak correlation in the differences between APOGEE DR14 - DR12 and LAMOST on DR12 surface gravity for stars hotter than 4800 K and in the log g range between 2.0 and 2.8.

We present an analysis of spectrosopic and astrometric data from APOGEE-2 and Gaia DR2 to identify structures towards the Orion Complex. By applying a hierarchical clustering algorithm to the 6-dimensional stellar data, we identify spatially and/or kinematically distinct groups of young stellar objects with ages ranging from 1 to 12 Myr. We also investigate the star forming history within the Orion Complex, and identify peculiar sub-clusters. With this method we reconstruct the older populations in the regions that are presently largely devoid of molecular gas, such as Orion C (which includes the $\sigma$ Ori cluster), and Orion D (the population that traces Ori OB1a, OB1b, and Orion X). We report on the distances, kinematics, and ages of the groups within the Complex. The Orion D groups is in the process of expanding. On the other hand, Orion B is still in the process of contraction. In $\lambda$ Ori the proper motions are consistent with a radial expansion due to an explosion from a supernova; the traceback age from the expansion exceeds the age of the youngest stars formed near the outer edges of the region, and their formation would have been triggered when they were half-way from the cluster center to their current positions. We also present a comparison between the parallax and proper motion solutions obtained by Gaia DR2, and those obtained towards star-forming regions by Very Long Baseline Array.

Prompted by peculiar spectroscopic variability observed in SDSS/APOGEE $H$-band spectra, we monitored the Be star HD 55606 using optical spectroscopy and found that it is an exotic double-lined spectroscopic binary (SB2) consisting of a Be star and a hot, compact companion that is probably an OB subdwarf (sdOB) star. Motion of the sdOB star is traced by its impact on the strong He~I lines, observed as radial velocity ($V_{r}$) variable, double-peaked emission profiles with narrow central absorption cores. Weak He II 4686 {\AA} absorption associated with the companion star is detected in most spectra. Use of the emission peaks of low-ionization emission lines to trace the Be star $V_{r}$ and the He I lines to trace the companion star $V_{r}$ yields a circular orbital solution with a 93.8-day period and masses of $M_{\rm Be}=6.2$ $M_{\rm \odot}$ and $M_{\rm sdOB}=0.9$ $M_{\rm \odot}$ in the case of $i=80^{\circ}$. HD 55606 exhibits a variety of phase-locked variability, including the development of shell lines twice per orbit. The shell phases coincide with variation in the double emission peak separations, and both forms of variability are likely caused by a two-armed spiral density perturbation in the Be disk. The intensity ratios of the double emission peaks are also phase-locked, possibly indicating heating by the sdOB star of the side of the Be disk facing it. HD 55606 is a new member of the growing sample of Be+sdOB binaries, in which the Be star's rapid rotation and ability to form a disk can be attributed to past mass transfer.

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material.

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 a catalog of stellar properties for a large sample of 6676 evolved stars with APOGEE spectroscopic parameters and \textit{Kepler} asteroseismic data analyzed using five independent techniques. Our data includes evolutionary state, surface gravity, mean density, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them. We employ a new empirical approach for combining asteroseismic measurements from different methods, calibrating the inferred stellar parameters, and estimating uncertainties. With high statistical significance, we find that asteroseismic parameters inferred from the different pipelines have systematic offsets that are not removed by accounting for differences in their solar reference values. We include theoretically motivated corrections to the large frequency spacing ($\Delta \nu$) scaling relation, and we calibrate the zero point of the frequency of maximum power ($\nu_{\rm max}$) relation to be consistent with masses and radii for members of star clusters. For most targets, the parameters returned by different pipelines are in much better agreement than would be expected from the pipeline-predicted random errors, but 22\% of them had at least one method not return a result and a much larger measurement dispersion. This supports the usage of multiple analysis techniques for asteroseismic stellar population studies. The measured dispersion in mass estimates for fundamental calibrators is consistent with our error model, which yields median random and systematic mass uncertainties for RGB stars of order 4\%. Median random and systematic mass uncertainties are at the 9\% and 8\% level respectively for RC stars.

The Orion Star Forming Complex (OSFC) is a central target for the APOGEE-2 Young Cluster Survey. Existing membership catalogs span limited portions of the OSFC, reflecting the difficulty of selecting targets homogeneously across this extended, highly structured region. We have used data from wide field photometric surveys to produce a less biased parent sample of young stellar objects (YSOs) with infrared (IR) excesses indicative of warm circumstellar material or photometric variability at optical wavelengths across the full 420 square degrees extent of the OSFC. When restricted to YSO candidates with H < 12.4, to ensure S/N ~100 for a six visit source, this uniformly selected sample includes 1307 IR excess sources selected using criteria vetted by Koenig & Liesawitz and 990 optical variables identified in the Pan-STARRS1 3$\pi$ survey: 319 sources exhibit both optical variability and evidence of circumstellar disks through IR excess. Objects from this uniformly selected sample received the highest priority for targeting, but required fewer than half of the fibers on each APOGEE-2 plate. We fill the remaining fibers with previously confirmed and new color-magnitude selected candidate OSFC members. Radial velocity measurements from APOGEE-1 and new APOGEE-2 observations taken in the survey's first year indicate that ~90% of the uniformly selected targets have radial velocities consistent with Orion membership.The APOGEE-2 Orion survey will include >1100 bona fide YSOs whose uniform selection function will provide a robust sample for comparative analyses of the stellar populations and properties across all sub-regions of Orion.

Detailed chemical abundance distributions for fourteen elements are derived for eight high-probability stellar members of the solar metallicity old open cluster M67 with an age of $\sim$4 Gyr. The eight stars consist of four pairs, with each pair occupying a distinct phase of stellar evolution: two G-dwarfs, two turnoff stars, two G-subgiants, and two red clump K-giants. The abundance analysis uses near-IR high-resolution spectra ($\lambda$1.5 -- 1.7$\mu$m) from the APOGEE survey and derives abundances for C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe. Our derived stellar parameters and metallicity for 2M08510076+113115 suggest that this star is a solar-twin, exhibiting abundance differences relative to the Sun of $\leq$ 0.04 dex for all elements. Chemical homogeneity is found within each class of stars ($\sim$0.02 dex), while significant abundance variations ($\sim$0.05 -- 0.20 dex) are found across the different evolutionary phases; the turnoff stars typically have the lowest abundances, while the red clump tend to have the largest. Non-LTE corrections to the LTE-derived abundances are unlikely to explain the differences. A detailed comparison of the derived Fe, Mg, Si, and Ca abundances with recently published surface abundances from stellar models that include chemical diffusion, provides a good match between the observed and predicted abundances as a function of stellar mass. Such agreement would indicate the detection of chemical diffusion processes in the stellar members of M67.

In order to study the growth and evolution of circumstellar disks around classical Be stars, we analyze optical time-series photometry from the KELT survey with simultaneous infrared and visible spectroscopy from the APOGEE survey and BeSS database for a sample of 160 Galactic classical Be stars. The systems studied here show variability including transitions from a diskless to a disk-possessing state (and vice versa), and persistent disks that vary in strength, being replenished at either regularly or irregularly occurring intervals. We detect disk-building events (outbursts) in the light curves of 28\% of our sample. Outbursts are more commonly observed in early- (57\%), compared to mid- (27\%) and late-type (8\%) systems. A given system may show anywhere between 0 -- 40 individual outbursts in its light curve, with amplitudes ranging up to $\sim$0.5 mag and event durations between $\sim$2 -- 1000 days. We study how both the photometry and spectroscopy change together during active episodes of disk growth or dissipation, revealing details about the evolution of the circumstellar environment. We demonstrate that photometric activity is linked to changes in the inner disk, and show that, at least in some cases, the disk growth process is asymmetrical. Observational evidence of Be star disks both growing and clearing from the inside out is presented. The duration of disk buildup and dissipation phases are measured for 70 outbursts, and we find that the average outburst takes about twice as long to dissipate as it does to build up in optical photometry. Our analysis hints that dissipation of the inner disk occurs relatively slowly for late-type Be stars.

We report the discovery of two ultra-faint satellites in the vicinity of the Large Magellanic Cloud (LMC) in data from the Magellanic Satellites Survey (MagLiteS). Situated 18$^{\circ}$ ($\sim 20$ kpc) from the LMC and separated from each other by only $18^\prime$, Carina~II and III form an intriguing pair. By simultaneously modeling the spatial and the color-magnitude stellar distributions, we find that both Carina~II and Carina~III are likely dwarf galaxies, although this is less clear for Carina~III. There are in fact several obvious differences between the two satellites. While both are well described by an old and metal poor population, Carina~II is located at $\sim 36$ kpc from the Sun, with $M_V\sim-4.5$ and $r_h\sim 90$ pc, and it is further confirmed by the discovery of 3 RR Lyrae at the right distance. In contrast, Carina~III is much more elongated, measured to be fainter ($M_V\sim-2.4$), significantly more compact ($r_h\sim30$ pc), and closer to the Sun, at $\sim 28$ kpc, placing it only 8 kpc away from Car~II. Together with several other systems detected by the Dark Energy Camera, Carina~II and III form a strongly anisotropic cloud of satellites in the vicinity of the Magellanic Clouds.

We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (\drvm) for APOGEE targets is a strong function of \logg, with main sequence stars showing \drvm\ as high as $\sim$300 \kms, and steadily dropping down to $\sim$30 \kms\ for \logg$\sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of \drvm\ values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high \drvm\ systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The \drvm\ distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$\lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$\gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.

We present the discovery of a very faint stellar system, SMASH 1, that is potentially a satellite of the Large Magellanic Cloud. Found within the Survey of the MAgellanic Stellar History (SMASH), SMASH 1 is a compact ($r_h = 9.1^{+5.9}_{-3.4}$ pc) and very low luminosity (M_V = -1.0 +/- 0.9, $L_V=10^{2.3 +/- 0.4}$ Lsun) stellar system that is revealed by its sparsely populated main sequence and a handful of red-giant-branch candidate member stars. The photometric properties of these stars are compatible with a metal-poor ([Fe/H]=-2.2) and old (13 Gyr) isochrone located at a distance modulus of ~18.8, i.e. a distance of ~57 kpc. Situated at 11.3$^\circ$ from the LMC in projection, its 3-dimensional distance from the Cloud is ~13 kpc, consistent with a connection to the LMC, whose tidal radius is at least 16 kpc. Although the nature of SMASH 1 remains uncertain, its compactness favors it being a stellar cluster and hence dark-matter free. If this is the case, its dynamical tidal radius is only <19 pc at this distance from the LMC, and smaller than the system's extent on the sky. Its low luminosity and apparent high ellipticity ($\epsilon=0.62^{+0.17}_{-0.21}$) with its major axis pointing toward the LMC may well be the tell-tale sign of its imminent tidal demise.

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.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

We report the discovery of a stellar over-density 8$^{\circ}$ north of the center of the Small Magellanic Cloud (Small Magellanic Cloud Northern Over-Density; SMCNOD) using data from the first two years of the Dark Energy Survey (DES) and the first year of the MAGellanic SatelLITEs Survey (MagLiteS). The SMCNOD is indistinguishable in age, metallicity and distance from the nearby SMC stars, being primarly composed of intermediate-age stars (6 Gyr, Z=0.001), with a small fraction of young stars (1 Gyr, Z=0.01). The SMCNOD has an elongated shape with an ellipticity of 0.6 and a size of $\sim$ 6x2 deg. It has an absolute magnitude of $M_V \cong$ -7.7, $r_h = 2.1$ kpc, and $\mu_V(r<r_h)$ = 31.2 mag arcsec$^{-2}$. We estimate a stellar mass of $\sim 10^5$ $M_{\odot}$, following a Kroupa mass function. The SMCNOD was probably removed from the SMC disk by tidal stripping, since it is located near the head of the Magellanic Stream, and the literature indicates likely recent LMC-SMC encounters. This scenario is supported by the lack of significant HI gas. Other potential scenarios for the SMCNOD origin are a transient over-density within the SMC tidal radius or a primordial SMC satellite in advanced stage of disruption.

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.

In the updated APOGEE-Kepler catalog, we have asteroseismic and spectroscopic data for over 3000 first ascent red giants. Given the size and accuracy of this sample, these data offer an unprecedented test of the accuracy of stellar models on the post-main-sequence. When we compare these data to theoretical predictions, we find a metallicity dependent temperature offset with a slope of around 100 K per dex in metallicity. We find that this effect is present in all model grids tested and that theoretical uncertainties in the models, correlated spectroscopic errors, and shifts in the asteroseismic mass scale are insufficient to explain this effect. Stellar models can be brought into agreement with the data if a metallicity dependent convective mixing length is used, with $ \Delta\alpha_{\rm ML, YREC} \sim 0.2$ per dex in metallicity, a trend inconsistent with the predictions of three dimensional stellar convection simulations. If this effect is not taken into account, isochrone ages for red giants from the Gaia data will be off by as much as a factor of 2 even at modest deviations from solar metallicity ([Fe/H]=$-$0.5).

We report a new ultra-faint stellar system found in Dark Energy Camera data from the first observing run of the Magellanic Satellites Survey (MagLiteS). MagLiteS J0644-5953 (Pictor II or Pic II) is a low surface brightness ({\mu} = 28.5 mag arcsec$^{-2}$ within its half-light radius) resolved overdensity of old and metal-poor stars located at a heliocentric distance of 45 kpc. The physical size (r$_{1/2}$ = 46 pc) and low luminosity (Mv = -3.2 mag) of this satellite are consistent with the locus of spectroscopically confirmed ultra-faint galaxies. MagLiteS J0644-5953 (Pic II) is located 11.3 kpc from the Large Magellanic Cloud (LMC), and comparisons with simulation results in the literature suggest that this satellite was likely accreted with the LMC. The close proximity of MagLiteS J0644-5953 (Pic II) to the LMC also makes it the most likely ultra-faint galaxy candidate to still be gravitationally bound to the LMC.

We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.

We present observations of near-infrared 2.12 micro-meter molecular hydrogen outflows emerging from 1.1 mm dust continuum clumps in the North America and Pelican Nebula (NAP) complex selected from the Bolocam Galactic Plane Survey (BGPS). Hundreds of individual shocks powered by over 50 outflows from young stars are identified, indicating that the dusty molecular clumps surrounding the NGC 7000 / IC 5070 / W80 HII region are among the most active sites of on-going star formation in the Solar vicinity. A spectacular X-shaped outflow, MHO 3400, emerges from a young star system embedded in a dense clump more than a parsec from the ionization front associated with the Pelican Nebula (IC 5070). Suspected to be a binary, the source drives a pair of outflows with orientations differing by 80 degrees. Each flow exhibits S-shaped symmetry and multiple shocks indicating a pulsed and precessing jet. The `Gulf of Mexico' located south of the North America Nebula (NGC 7000), contains a dense cluster of molecular hydrogen objects (MHOs), Herbig-Haro (HH) objects, and over 300 YSOs, indicating a recent burst of star formation. The largest outflow detected thus far in the North America and Pelican Nebula complex, the 1.6 parsec long MHO 3417 flow, emerges from a 500 Solar mass BGPS clump and may be powered by a forming massive star. Several prominent outflows such as MHO 3427 appear to be powered by highly embedded YSOs only visible at a wavelength > 70 micro-meters. An `activity index' formed by dividing the number of shocks by the mass of the cloud containing their source stars is used to estimate the relative evolutionary states of Bolocam clumps. Outflows can be used as indicators of the evolutionary state of clumps detected in mm and sub-mm dust continuum surveys.

We present optical and near-IR spectroscopic observations of the luminous blue variable SN 2009ip during its remarkable photometric evolution of 2012. The spectra sample three key points in the SN 2009ip lightcurve, corresponding to its initial brightening in August (2012-A) and its dramatic rebrightening in early October (2012-B). Based on line fluxes and velocities measured in our spectra, we find a surprisingly low I(H-alpha)/I(H-beta) ~ 1.3-1.4 in the 2012-B spectra. Such a ratio implies either a rare Case B recombination scenario where H-alpha, but not H-beta, is optically thick, or an extremely high density for the circumstellar material of n_e > 10^13 cm^(-3) The H-alpha line intensity yields a minimum radiating surface area of >~20,000 AU^2 in H-alpha at the peak of SN 2009ip's photometric evolution. Combined with the nature of this object's spectral evolution in 2012, a high circumstellar density and large radiating surface area imply the presence of a thin disk geometry around the central star (and, consequently, a possible binary companion), suggesting that the observed 2012-B rebrightening of SN 2009ip can be attributed to the illumination of the disk's inner rim by fast-moving ejecta produced by the underlying events of 2012-A.

The Bolocam Galactic Plane Survey (BGPS) is a 1.1 mm continuum survey of dense clumps of dust throughout the Galaxy covering 170 square degrees. We present spectroscopic observations using the Heinrich Hertz Submillimeter Telescope of the dense gas tracers, HCO+ and N2H+ 3-2, for all 6194 sources in the Bolocam Galactic Plane Survey v1.0.1 catalog between 7.5 <= l <= 194 degrees. This is the largest targeted spectroscopic survey of dense molecular gas in the Milky Way to date. We find unique velocities for 3126 (50.5%) of the BGPS v1.0.1 sources observed. Strong N2H+ 3-2 emission (T_{mb} > 0.5 K) without HCO+ 3-2 emission does not occur in this catalog. We characterize the properties of the dense molecular gas emission toward the entire sample. HCO+ is very sub-thermally populated and the 3-2 transitions are optically thick toward most BGPS clumps. The median observed line width is 3.3 km/s consistent with supersonic turbulence within BGPS clumps. We find strong correlations between dense molecular gas integrated intensities and 1.1 mm peak flux and the gas kinetic temperature derived from previously published NH3 observations. These intensity correlations are driven by the sensitivity of the 3-2 transitions to excitation conditions rather than by variations in molecular column density or abundance. We identify a subset of 113 sources with stronger N2H+ than HCO+ integrated intensity, but we find no correlations between the N2H+ / HCO+ ratio and 1.1 mm continuum flux density, gas kinetic temperature, or line width. Self-absorbed profiles are rare (1.3%).

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-\sigma RMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

Two optically obscured Wolf-Rayet (WR) stars have been recently discovered by means of their infrared (IR) circumstellar shells, which show signatures of interaction with each other. Following the systematics of the WR star catalogues, these stars obtain the names WR\,120bb and WR\,120bc. In this paper, we present and analyse new near-IR, $J$, $H$, and $K$-band, spectra using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. For that purpose, the atomic data base of the code has been extended in order to include all significant lines in the near-IR bands. The spectra of both stars are classified as WN9h. As their spectra are very similar the parameters that we obtained by the spectral analyses hardly differ. Despite their late spectral subtype, we found relatively high stellar temperatures of 63 kK. The wind composition is dominated by helium, while hydrogen is depleted to 25 per cent by mass. Because of their location in the Scutum-Centaurus arm, WR\,120bb and WR\,120bc appear highly reddened, $A_{K_{\rm s}} \approx 2$ mag. We adopt a common distance of 5.8\,kpc to both stars, which complies with the typical absolute $K$-band magnitude for the WN9h subtype of -6.5 mag, is consistent with their observed extinction based on comparison with other massive stars in the region, and allows for the possibility that their shells are interacting with each other. This leads to luminosities of $\log(L/L_\odot) = 5.66$ and 5.54 for WR\,120bb and WR\,120bc, with large uncertainties due to the adopted distance. The values of the luminosities of WR\,120bb and WR\,120bc imply that the immediate precursors of both stars were red supergiants (RSG). This implies in turn that the circumstellar shells associated with WR\,120bb and WR\,120bc were formed by interaction between the WR wind and the dense material shed during the preceding RSG phase.

The last part of SpS5 dealt with the circumstellar environment. Structures are indeed found around several types of massive stars, such as blue and red supergiants, as well as WRs and LBVs. As shown in the last years, the potential of IR for their study is twofold: first, IR can help discover many previously unknown nebulae, leading to the identification of new massive stars as their progenitors; second, IR can help characterize the nebular features. Current and new IR facilities thus pave the way to a better understanding of the feedback from massive stars.

We introduce the Stony Brook / SMARTS Atlas of (mostly) Southern Novae. This atlas contains both spectra and photometry obtained since 2003. The data archived in this atlas will facilitate systematic studies of the nova phenomenon and correlative studies with other comprehensive data sets. It will also enable detailed investigations of individual objects. In making the data public we hope to engender more interest on the part of the community in the physics of novae. The atlas is on-line at \url{http://www.astro.sunysb.edu/fwalter/SMARTS/NovaAtlas/} .

We report the discovery of two new Galactic candidate luminous blue variable (cLBV) stars via detection of circular shells (typical of known confirmed and cLBVs) and follow-up spectroscopy of their central stars. The shells were detected at 22 um in the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars of the shells conducted with the renewed Southern African Large Telescope (SALT) showed that their spectra are very similar to those of the well-known LBVs P Cygni and AG Car, and the recently discovered cLBV MN112, which implies the LBV classification for these stars as well. The LBV classification of both stars is supported by detection of their significant photometric variability: one of them brightened in the R- and I-bands by 0.68\pm0.10 mag and 0.61\pm0.04 mag, respectively, during the last 13-18 years, while the second one (known as Hen 3-1383) varies its B,V,R,I and K_s brightnesses by \simeq 0.5-0.9 mag on time-scales from 10 days to decades. We also found significant changes in the spectrum of Hen 3-1383 on a timescale of \simeq 3 months, which provides additional support for the LBV classification of this star. Further spectrophotometric monitoring of both stars is required to firmly prove their LBV status. We discuss a connection between the location of massive stars in the field and their fast rotation, and suggest that the LBV activity of the newly discovered cLBVs might be directly related to their possible runaway status.

We have undertaken a near-IR spectral survey of stars associated with compact nebulae recently revealed by the Spitzer and WISE imaging surveys. These circumstellar nebulae, produced by massive evolved stars, display a variety of symmetries and shapes and are often only evident at mid-IR wavelengths. Stars associated with ~50 of these nebulae have been observed. We also obtained recent spectra of previously confirmed (known) luminous blue variables (LBVs) and candidate LBVs (cLBVs). The spectral similarity of the stars observed when compared directly to known LBVs and Wolf-Rayet (WR) stars indicate many are newly identified cLBVs, with a few being newly discovered WR stars, mostly of WN8-9h spectral type. These results suggest that a large population of previously unidentified cLBVs and related transitional stars reside in the Galaxy and confirm that circumstellar nebulae are inherent to most (c)LBVs.

We have undertaken a near-infrared spectral survey of stars associated with compact mid-IR shells recently revealed by the MIPSGAL (24 micron) and GLIMPSE (8 micron) Spitzer surveys, whose morphologies are typical of circumstellar shells produced by massive evolved stars. Through spectral similarity with known Luminous Blue Variable (LBV) and Wolf-Rayet (WR) stars, a large population of candidate LBVs (cLBVs) and a smaller number of new WR stars are being discovered. This significantly increases the Galactic cLBV population and confirms that nebulae are inherent to most (if not all) objects of this class. Keywords - stars: emission-line, Be, stars: mass loss, stars: winds, outflows, stars: Wolf-Rayet

Type Ia supernovae are key tools for measuring distances on a cosmic scale. They are generally thought to be the thermonuclear explosion of an accreting white dwarf in a close binary system. The nature of the mass donor is still uncertain. In the single-degenerate model it is a main-sequence star or an evolved star, whereas in the double-degenerate model it is another white dwarf. We show that the velocity structure of absorbing material along the line of sight to 35 type Ia supernovae tends to be blueshifted. These structures are likely signatures of gas outflows from the supernova progenitor systems. Thus many type Ia supernovae in nearby spiral galaxies may originate in single-degenerate systems.

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

We present the results of observations of 1882 sources in the Bolocam Galactic Plane Survey (BGPS) at 1.1 mm with the 10m Heinrich Hertz Telescope simultaneously in HCO+ J=3-2 and N2H+ J=3-2. We detect 77% of these sources in HCO^+ and 51% in N2H+ at greater than 3$\sigma$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sources. We determine kinematic distances for 529 sources (440 in the first quadrant breaking the distance ambiguity and 89 in the second quadrant) We derive the size, mass, and average density for this subset of clumps. The median size of BGPS clumps is 0.75 pc with a median mass of 330 M$_{\odot}$ (assuming T_{Dust}=20 K). The median HCO+ linewidth is 2.9 km s$^{-1}$ indicating that BGPS clumps are dominated by supersonic turbulence or unresolved kinematic motions. We find no evidence for a size-linewidth relationship for BGPS clumps. We analyze the effects of the assumed dust temperature on the derived clump properties with a Monte Carlo simulation and we find that changing the temperature distribution will change the median source properties (mass, volume-averaged number density, surface density) by factors of a few. The observed differential mass distribution has a power-law slope that is intermediate between that observed for diffuse CO clouds and the stellar IMF. BGPS clumps represent a wide range of objects (from dense cores to more diffuse clumps) and are typically characterized by larger sizes and lower densities than previously published surveys of high-mass star-forming regions. This collection of objects is a less-biased sample of star-forming regions in the Milky Way that likely span a wide range of evolutionary states.

We use a Minimum Spanning Tree algorithm to characterize the spatial distribution of Galactic Far-IR sources and derive their clustering properties. We aim to reveal the spatial imprint of different types of star forming processes, e.g. isolated spontaneous fragmentation of dense molecular clouds, or events of triggered star formation around HII regions, and highlight global properties of star formation in the Galaxy. We plan to exploit the entire Hi-GAL survey of the inner Galactic plane to gather significant statistics on the clustering properties of star forming regions, and to look for possible correlations with source properties such as mass, temperature or evolutionary stage. In this paper we present a pilot study based on the two 2x2 square degree fields centered at longitudes l=30 and l=59 obtained during the Science Demonstration Phase (SDP) of the Herschel mission. We find that over half of the clustered sources are associated with HII regions and infrared dark clouds. Our analysis also reveals a smooth chromatic evolution of the spatial distribution where sources detected at short-wavelengths, likely proto-stars surrounded by warm circumstellar material emitting in the far-infrared, tend to be clustered in dense and compact groups around HII regions while sources detected at long-wavelengths, presumably cold and dusty density enhancements of the ISM emitting in the sub-millimeter, are distributed in larger and looser groups.

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 {\mu}m images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 {\mu}m background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 {\mu}m, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

We present a first study of the star-forming compact dust condensations revealed by Herschel in the two 2 \times 2 \degr Galactic Plane fields centered at [l;b] = [30\degr; 0 \degr] and [l;b] = [59\degr; 0 \degr], respectively, and observed during the Science Demonstration Phase for the Herschel infrared Galactic Plane survey (Hi-GAL) Key-Project. Compact source catalogs extracted for the two fields in the five Hi-GAL bands (70, 160, 250, 350 and 500 $\mu$m) were merged based on simple criteria of positional association and spectral energy distribution (SED) consistency into a final catalog which contains only coherent SEDs with counterparts in at least three adjacent Herschel bands. These final source lists contain 528 entries for the l = 30\degr field, and 444 entries for the l = 59\degr field. The SED coverage has been augmented with ancillary data at 24 $\mu$m and 1.1 mm. SED modeling for the subset of 318 and 101 sources (in the two fields, respectively) for which the distance is known was carried out using both a structured star/disk/envelope radiative transfer model and a simple isothermal grey-body. Global parameters like mass, luminosity, temperature and dust properties have been estimated. The Lbol/Menv ratio spans four orders of magnitudes from values compatible with the pre-protostellar phase to embedded massive zero-age main sequence stars. Sources in the l = 59\degr field have on average lower L/M, possibly outlining an overall earlier evolutionary stage with respect to the sources in the l = 30\degr field. Many of these cores are actively forming high-mass stars, although the estimated core surface densities appear to be an order of magnitude below the 1 g cm$^{-2}$ critical threshold for high-mass star formation.

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

Abriged version for astroph: The young late-type star V1118 Orionis was in outburst from 2005 to 2006. We followed the outburst with optical and near-infrared photometry; the X-ray emission was further probed with observations taken with XMM-Newton and Chandra during and after the outburst. In addition, we obtained mid-infrared photometry and spectroscopy with Spitzer at the peak of the outburst and in the post-outburst phase. The spectral energy distribution of V1118 Ori varied significantly over the course of the outburst. The optical flux showed the largest variations, most likely due to enhanced emission by a hot spot. The latter dominated the optical and near-infrared emission at the peak of the outburst, while the disk emission dominated in the mid-infrared. The X-ray flux correlated with the optical and infrared fluxes, indicating that accretion affected the magnetically active corona and the stellar magnetosphere. The thermal structure of the corona was variable with some indication of a cooling of the coronal temperature in the early phase of the outburst with a gradual return to normal values. Color-color diagrams in the optical and infrared showed variations during the outburst, with no obvious signature of reddening due to circumstellar matter. Using MC realizations of star+disk+hotspot models to fit the SED in ``quiescence'' and at the peak of the outburst, we determined that the mass accretion rate varied from about 2.5E-7 Msun/yr to 1E-6 Msun/yr; in addition the fractional area of the hotspot increased significantly as well. The multi-wavelength study of the V1118 Ori outburst helped us to understand the variations in spectral energy distributions and demonstrated the interplay between the disk and the stellar magnetosphere in a young, strongly accreting star.

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

We present extensive early photometric (ultraviolet through near-infrared) and spectroscopic (optical and near-infrared) data on supernova (SN) 2008D as well as X-ray data analysis on the associated Swift/X-ray transient (XRT) 080109. Our data span a time range of 5 hours before the detection of the X-ray transient to 150 days after its detection, and detailed analysis allowed us to derive constraints on the nature of the SN and its progenitor; throughout we draw comparisons with results presented in the literature and find several key aspects that differ. We show that the X-ray spectrum of XRT 080109 can be fit equally well by an absorbed power law or a superposition of about equal parts of both power law and blackbody. Our data first established that SN 2008D is a spectroscopically normal SN Ib (i.e., showing conspicuous He lines), and show that SN 2008D had a relatively long rise time of 18 days and a modest optical peak luminosity. The early-time light curves of the SN are dominated by a cooling stellar envelope (for \Delta t~0.1- 4 day, most pronounced in the blue bands) followed by 56^Ni decay. We construct a reliable measurement of the bolometric output for this stripped-envelope SN, and, combined with estimates of E_K and M_ej from the literature, estimate the stellar radius R_star of its probable Wolf-Rayet progenitor. According to the model of Waxman et al. and of Chevalier & Fransson, we derive R_star^{W07}= 1.2+/-0.7 R_sun and R_star^{CF08}= 12+/-7 R_sun, respectively; the latter being more in line with typical WN stars. Spectra obtained at 3 and 4 months after maximum light show double-peaked oxygen lines that we associate with departures from spherical symmetry, as has been suggested for the inner ejecta of a number of SN Ib cores.

Gliese 436 is an M dwarf with a mass of 0.45 Msun and hosts the extrasolar planet GL 436b [3, 6, 7, 2], which is currently the least massive transiting planet with a mass of ~23.17 Mearth [10], and the only planet known to transit an M dwarf. GL 436b represents the first transiting detection of the class of extrasolar planets known as "Hot Neptunes" that have masses within a few times that of Neptune's mass (~17 Mearth) and orbital semimajor axis <0.1 AU about the host star. Unlike most other known transiting extrasolar planets, GL 436b has a high eccentricity (e~0.16). This brings to light a new parameter space for habitability zones of extrasolar planets with host star masses much smaller than typical stars of roughly a solar mass. This unique system is an ideal candidate for orbital perturbation and transit-time variation (TTV) studies to detect smaller, possibly Earth-mass planets in the system. In April 2008 we began a long-term intensive campaign to obtain complete high-precision light curves using the Apache Point Observatory's 3.5-meter telescope, NMSU's 1-meter telescope (located at APO), and Sommers Bausch Observatory's 24" telescope. These light curves are being analyzed together, along with amateur and other professional astronomer observations. Results of our analysis are discussed. Continued measurements over the next few years are needed to determine if additional planets reside in the system, and to study the impact of other manifestations on the light curves, such as star spots and active regions.

We present ground-based observations of the transiting Neptune-mass planet Gl 436b obtained with the 3.5-meter telescope at Apache Point Observatory and other supporting telescopes. Included in this is an observed transit in early 2005, over two years before the earliest reported transit detection. We have compiled all available transit data to date and perform a uniform modeling using the JKTEBOP code. We do not detect any transit timing variations of amplitude greater than ~1 minute over the ~3.3 year baseline. We do however find possible evidence for a self-consistent trend of increasing orbital inclination, transit width, and transit depth, which supports the supposition that Gl 436b is being perturbed by another planet of < 12 Earth Masses in a non-resonant orbit.

We present results from a multi-wavelength campaign to monitor the 2005 outburst of the low-mass young star V1118 Ori. Although our campaign covers the X-ray, optical, infrared, and radio regimes, we focus in this Letter on the properties of the X-ray emission in V1118 Ori during the first few months after the optical outburst. Chandra and XMM-Newton detected V1118 Ori at three epochs in early 2005. The X-ray flux and luminosity stayed similar within a factor of two, and at the same level as in a pre-outburst observation in 2002. The hydrogen column density showed no evidence for variation from its modest pre-outburst value of $N_\mathrm{H} \sim 3 \times 10^{21}$ cm$^{-2}$. However, a spectral change occurred from a dominant hot plasma ($\sim 25$ MK) in 2002 and in January 2005 to a cooler plasma ($\sim 8$ MK) in February 2005 and in March 2005. We hypothesize that the hot magnetic loops high in the corona were disrupted by the closing in of the accretion disk due to the increased accretion rate during the outburst, whereas the lower cooler loops were probably less affected and became the dominant coronal component.

We present optical and/or infrared photometry of the Type Ia supernovae SN 1991T, SN 1991bg, SN 1999ek, SN 2001bt, SN 2001cn, SN 2001cz, and SN 2002bo. All but one of these supernovae have decline rate parameters Delta m_15(B) close to the median value of 1.1 for the whole class of Type Ia supernovae. The addition of these supernovae to the relationship between the near-infrared absolute magnitudes and Delta m_15(B) strengthens the previous relationships we have found, in that the maximum light absolute magnitudes are essentially independent of the decline rate parameter. (SN 1991bg, the prototype of the subclass of fast declining Type Ia supernovae, is a special case.) The dispersion in the Hubble diagram in JHK is only ~0.15 mag. The near-infrared properties of Type Ia supernovae continue to be excellent measures of the luminosity distances to the supernova host galaxies, due to the need for only small corrections from the epoch of observation to maximum light, low dispersion in absolute magnitudes at maximum light, and the minimal reddening effects in the near-infrared.

We report on photometric (BVRIJHK) and low dispersion spectroscopic observations of V1647 Ori, the star that drives McNeil's Nebula, between 10 February and 7 May 2004. The star is photometrically variable atop a general decline in brightness of about 0.3-0.4 magnitudes during these 87 days. The spectra are featureless, aside from H-alpha and the Ca II infrared triplet in emission, and a Na I D absorption feature. The Ca II triplet line ratios are typical of young stellar objects. The H-alpha equivalent width may be modulated on a period of about 60 days. The post-outburst extinction appears to be less than 7 mag. The data are suggestive of an FU Orionis-like event, but further monitoring will be needed to definitively characterize the outburst.

We compiled and investigated the infrared/sub-mm/mm SED of the new outburst star IRAS 05436-0007 in quiescent phase. The star is a flat-spectrum source, with an estimated total luminosity of L_bol ~ 5.6 L_sun, typical of low-mass T Tauri stars. The derived circumstellar mass of 0.5 M_sun is rather high among low-mass YSOs. The observed SED differs from the SEDs of typical T Tauri stars and of 4 well-known EXors, and resembles more the SEDs of FU Orionis objects indicating the presence of a circumstellar envelope. IRAS 05436-0007 seems to be a Class II source with an age of approximately 4x10^5 yr. In this evolutionary stage an accretion disk is already fully developed, though a circumstellar envelope may also be present. Observations of the present outburst will provide additional knowledge on the source.