We extend our previous work on the age-chemical abundance structure of the Galactic outer disc to the inner disc (4 < r < 8 kpc) based on the SDSS/APOGEE survey. Different from the outer disc, the inner disc stars exhibit a clear bimodal distribution in the [Mg/Fe]-[Fe/H] plane. While a number of scenarios have been proposed in the literature, it remains challenging to recover this bimodal distribution with theoretical models. To this end, we present a chemical evolution model embedding a complex multi-phase inner disc formation scenario that matches the observed bimodal [Mg/Fe]-[Fe/H] distribution. In this scenario, the formation of the inner disc is dominated by two main starburst episodes 6 Gyr apart with secular, low-level star formation activity in between. In our model, the first starburst occurs at early cosmic times (t~1 Gyr) and the second one 6 Gyr later at a cosmic time of t~7 Gyr. Both these starburst episodes are associated with gas accretion events in our model, and are quenched rapidly. The first starburst leads to the formation of the high-$\alpha$ sequence, and the second starburst leads to the formation of the metal-poor low-$\alpha$ sequence. The metal-rich low-$\alpha$ stars, instead, form during the secular evolution phase between the two bursts. Our model shows that the $\alpha$-dichotomy originates from the rapid suppression of star formation after the first starburst. The two starburst episodes are likely to be responsible for the formation of the geometric thick disc (z >1 kpc), with the old inner thick disc and the young outer thick disc forming during the first and the second starbursts, respectively.

Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups.

Wolf-Rayet (WR) galaxies are a rare population of galaxies that host living high-mass stars during their WR phase (i.e. WR stars) and are thus expected to provide interesting constraints on the stellar Initial Mass Function, massive star formation, stellar evolution models, etc. Spatially resolved spectroscopy should in principle provide a more efficient way of identifying WR galaxies than single-fiber surveys of galactic centers such as SDSS-I & II, as WR stars should be more preferentially found in discs. Using IFU data from the ongoing SDSS-IV MaNGA survey, we have performed a thorough search for WR galaxies. We first identify H II regions in each datacube and carry out full spectral fitting to the stacked spectra. We then visually inspect the residual spectrum of each H II region and identify WR regions that present a significant "blue bump" at 4600-4750 A. The resulting WR catalog includes 267 WR regions of ~500pc (radius) sizes, distributed in 90 galaxies from the current sample of MaNGA (MaNGA Product Launch 7). We find WR regions are exclusively found in galaxies that show bluest colors and highest star formation rates for their mass. Most WR galaxies have late-type morphologies and show relatively large asymmetry in their images, implying that WR regions are more preferentially found in interacting/merging galaxies. We estimate the stellar mass function of WR galaxies and the mass-dependent detection rate. The detection rate of WR galaxies is typically ~2%, with weak dependence on stellar mass. This detection rate is about 40 times higher than previous studies with SDSS single fiber data, and by a factor of 2 lower than the CALIFA-based WR catalog. We make comparisons with SDSS and CALIFA studies, and conclude that different detection rates can be explained mainly by three factors: spatial coverage, spectral signal-to-noise ratio, and redshift ranges of the parent sample.

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

The Open Cluster Chemical Abundances and Mapping (OCCAM) survey aims to constrain key Galactic dynamical and chemical evolution parameters by the construction of a large, comprehensive, uniform, infrared-based spectroscopic data set of hundreds of open clusters. This fourth contribution from the OCCAM survey presents analysis using SDSS/APOGEE DR16 of a sample of 128 open clusters, 71 of which we designate to be "high quality" based on the appearance of their color-magnitude diagram. We find the APOGEE DR16 derived [Fe/H] abundances to be in good agreement with previous high resolution spectroscopic open cluster abundance studies. Using the high quality sample, we measure Galactic abundance gradients in 16 elements, and find evolution of some of the [X/Fe] gradients as a function of age. We find an overall Galactic [Fe/H] vs R_GC gradient of $-0.068 \pm 0.001$ dex kpc$^{-1}$ over the range of $6 <$ R_GC $< 13.9$ kpc; however, we note that this result is sensitive to the distance catalog used, varying as much as 15%. We formally derive the location a break in the [Fe/H] abundance gradient as a free parameter in the gradient fit for the first time. We also measure significant Galactic gradients in O, Mg, S, Ca, Mn, Cr, Cu, Na, Al, and K, some of which are measured for the first time. Our large sample allows us to explore four well-populated age bins to explore the time evolution of gradients for a large number of elements and comment on possible implications for Galactic chemical evolution and radial migration.

Studies of the kinematics and chemical compositions of Galactic globular clusters (GCs) enable the reconstruction of the history of star formation, chemical evolution, and mass assembly of the Galaxy. Using the latest data release (DR16) of the SDSS/APOGEE survey, we identify 3,090 stars associated with 46 GCs. Using a previously defined kinematic association, we break the sample down into eight separate groups and examine how the kinematics-based classification maps into chemical composition space, considering only $\alpha$ (mostly Si and Mg) elements and Fe. Our results show that: (i) The loci of both in situ and accreted subgroups in chemical space match those of their field counterparts; (ii) GCs from different individual accreted subgroups occupy the same locus in chemical space. This could either mean that they share a similar origin or that they are associated with distinct satellites which underwent similar chemical enrichment histories; (iii) The chemical compositions of the GCs associated with the low orbital energy subgroup defined by Massari and collaborators is broadly consistent with an in situ origin. However, at the low metallicity end, the distinction between accreted and in situ populations is blurred; (iv) Regarding the status of GCs whose origin is ambiguous, we conclude the following: the position in Si-Fe plane suggests an in situ}origin for Liller 1 and a likely accreted origin for NGC 5904 and NGC 6388. The case of NGC 288 is unclear, as its orbital properties suggest an accretion origin, its chemical composition suggests it may have formed in situ.

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.

Studies of the ages, abundances, and motions of individual stars in the Milky Way provide one of the best ways to study the evolution of disk galaxies over cosmic time. The formation of the Milky Way's barred inner region in particular is a crucial piece of the puzzle of disk galaxy evolution. Using data from APOGEE and Gaia, we present maps of the kinematics, elemental abundances, and age of the Milky Way bulge and disk that show the barred structure of the inner Milky Way in unprecedented detail. The kinematic maps allow a direct, purely kinematic determination of the bar's pattern speed of 41+/-3 km/s/kpc and of its shape and radial profile. We find the bar's age, metallicity, and abundance ratios to be the same as those of the oldest stars in the disk that are formed in its turbulent beginnings, while stars in the bulge outside of the bar are younger and more metal-rich. This implies that the bar likely formed ~8 Gyr ago, when the decrease in turbulence in the gas disk allowed a thin disk to form that quickly became bar-unstable. The bar's formation therefore stands as a crucial epoch in the evolution of the Milky Way, a picture that is in line with the evolutionary path that emerges from observations of the gas kinematics in external disk galaxies over the last ~10 Gyr.

An improved high-resolution and deep A$_{Ks}$ foreground dust extinction map is presented for the Galactic disk area within $295^{\circ} \lesssim l \lesssim 350^{\circ}$, $-1.0^{\circ} \lesssim b \lesssim +1.0^{\circ}$. At some longitudes the map reaches up to $|b|\sim2.25^{\circ}$, for a total of $\sim$148 deg$^2$. The map was constructed via the Rayleigh-Jeans Color Excess (RJCE) technique based on deep near-infrared (NIR) and mid-infrared (MIR) photometry. The new extinction map features a maximum bin size of 1', and relies on NIR observations from the Two Micron All-Sky Survey (2MASS) and new data from ESO's Vista Variables in the V\'ia L\'actea (VVV) survey, in concert with MIR observations from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE). The VVV photometry penetrates $\sim$4 magnitudes fainter than 2MASS, and provides enhanced sampling of the underlying stellar populations in this heavily obscured region. Consequently, the new results supersede existing RJCE maps tied solely to brighter photometry, revealing a systematic underestimation of extinction in prior work that was based on shallower data. The new high-resolution and large-scale extinction map presented here is readily available to the community through a web query interface.

Stellar ages are a crucial component to studying the evolution of the Milky Way. Using Gaia DR2 distance estimates, it is now possible to estimate stellar ages for a larger volume of evolved stars through isochrone matching. This work presents [M/H]-age and [$\alpha$/M]-age relations derived for different spatial locations in the Milky Way disc. These relations are derived by hierarchically modelling the star formation history of stars within a given chemical abundance bin. For the first time, we directly observe that significant variation is apparent in the [M/H]-age relation as a function of both Galactocentric radius and distance from the disc mid-plane. The [M/H]-age relations support claims that radial migration has a significant effect in the plane of the disc. Using the [M/H] bin with the youngest mean age at each radial zone in the plane of the disc, the present-day metallicity gradient is measured to be $-0.059 \pm 0.010$ dex kpc$^{-1}$, in agreement with Cepheids and young field stars. We find a vertically flared distribution of young stars in the outer disc, confirming predictions of models and previous observations. The mean age of the [M/H]-[$\alpha$/M] distribution of the solar neighborhood suggests that the high-[M/H] stars are not an evolutionary extension of the low-$\alpha$ sequence. Our observational results are important constraints to Galactic simulations and models of chemical evolution.

We present elemental abundance analysis of high-resolution spectra for five giant stars, deriving Fe, Mg, Al, C, N, O, Si and Ce abundances, and spatially located within the innermost regions of the bulge globular cluster NGC 6522, based on H-band spectra taken with the multi-object APOGEE-north spectrograph from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. Of the five cluster candidates, two previously unremarked stars are confirmed to have second-generation (SG) abundance patterns, with the basic pattern of depletion in C and Mg simultaneous with enrichment in N and Al as seen in other SG globular cluster populations at similar metallicity. } In agreement with the most recent optical studies, the NGC 6522 stars analyzed exhibit (when available) only mild overabundances of the s-process element Ce, contradicting the idea of the NGC 6522 stars being formed from gas enriched by spinstars and indicating that other stellar sources such as massive AGB stars could be the primary intra-cluster medium polluters. The peculiar abundance signature of SG stars have been observed in our data, confirming the presence of multiple generations of stars in NGC 6522.

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

We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Galaxy's chemical and kinematical history. It was part of the Sloan Digital Sky Survey III (SDSS-III) from 2011 -- 2014 using the 2.5 m Sloan Foundation Telescope at Apache Point Observatory, New Mexico. The APOGEE-2 survey is now using the spectrograph as part of SDSS-IV, as well as a second spectrograph, a close copy of the first, operating at the 2.5 m du Pont Telescope at Las Campanas Observatory in Chile. Although several fiber-fed, multi-object, high resolution spectrographs have been built for visual wavelength spectroscopy, the APOGEE spectrograph is one of the first such instruments built for observations in the near-infrared. The instrument's successful development was enabled by several key innovations, including a "gang connector" to allow simultaneous connections of 300 fibers; hermetically sealed feedthroughs to allow fibers to pass through the cryostat wall continuously; the first cryogenically deployed mosaic volume phase holographic grating; and a large refractive camera that includes mono-crystalline silicon and fused silica elements with diameters as large as ~ 400 mm. This paper contains a comprehensive description of all aspects of the instrument including the fiber system, optics and opto-mechanics, detector arrays, mechanics and cryogenics, instrument control, calibration system, optical performance and stability, lessons learned, and design changes for the second instrument.

Kepler-730 is a planetary system hosting a statistically validated hot Jupiter in a 6.49-day orbit and an additional transiting candidate in a 2.85-day orbit. We use spectroscopic radial velocities from the APOGEE-2N instrument, Robo-AO contrast curves, and Gaia distance estimates to statistically validate the planetary nature of the additional Earth-sized candidate. We perform astrophysical false positive probability calculations for the candidate using the available Kepler data and bolster the statistical validation by using radial velocity data to exclude a family of possible binary star solutions. Using a radius estimate for the primary star derived from stellar models, we compute radii of $1.100^{+0.047}_{-0.050}\ R_{Jup}$ and $0.140\pm0.012\ R_{Jup}$ ($1.57\pm0.13\ R_{\oplus}$) for Kepler-730b and Kepler-730c, respectively. Kepler-730 is only the second compact system hosting a hot Jupiter with an inner, transiting planet.

In this letter we present evidence of coherent outward motion of a sample of young stars ($t<$30 Myr) in the Perseus Arm, whose apparent origin is located in the vicinity of the W3/W4/W5 complex. Using astrometric and photometric data from the Gaia DR2 catalog of an 8$^\circ$ radius field centered near W4, we selected a sample of young, intermediate to high-mass star candidates. The sample is limited to sources with parallax uncertainties below 20% and Bayesian distance estimates within 1800 and 3100 pc. The selection includes embedded stellar populations as well as young open clusters. Projected velocities derived from perspective-corrected proper motions clearly suggest that the young star population emerge from the Perseus arm, with a possible convergence zone near W3/W4/W5 region, tracing a front that expands away at a rate of about $15~{\rm km~s}^{-1}~{\rm kpc}^{-1}$.

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

We present a search for new members of the 300 km s$^{-1}$ stream (300S) near the dwarf galaxy Segue 1 using wide-field survey data. We identify 11 previously unknown bright stream members in the APOGEE-2 and SEGUE-1 and 2 spectroscopic surveys. Based on the spatial distribution of the high-velocity stars, we confirm for the first time that this kinematic structure is associated with a 24$^{\circ}$-long stream seen in SDSS and Pan-STARRS imaging data. The 300S stars display a metallicity range of $-2.17 < {\rm [Fe/H]} < -1.24$, with an intrinsic dispersion of 0.21$_{-0.09}^{+0.12}$ dex. They also have chemical abundance patterns similar to those of Local Group dwarf galaxies, as well as that of the Milky Way halo. Using the open-source code galpy to model the orbit of the stream, we find that the progenitor of the stream passed perigalacticon about 70 Myr ago, with a closest approach to the Galactic Center of about 4.1 kpc. Using Pan-STARRS DR1 data, we obtain an integrated stream luminosity of $4 \times 10^3$ L$_{\odot}$. We conclude that the progenitor of the stream was a dwarf galaxy that is probably similar to the satellites that were accreted to build the present-day Milky Way halo.

IC 166 is an intermediate-age open cluster ($\sim 1$ Gyr) which lies in the transition zone of the metallicity gradient in the outer disc. Its location, combined with our very limited knowledge of its salient features, make it an interesting object of study. We present the first high-resolution spectroscopic and precise kinematical analysis of IC 166, which lies in the outer disc with $R_{GC} \sim 12.7$ kpc. High resolution \textit{H}-band spectra were analyzed using observations from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We made use of the Brussels Automatic Stellar Parameter (BACCHUS) code to provide chemical abundances based on a line-by-line approach for up to eight chemical elements (Mg, Si, Ca, Ti, Al, K, Mn and Fe). The $\alpha-$element (Mg, Si, Ca and whenever available Ti) abundances, and their trends with Fe abundances have been analysed for a total of 13 high-likelihood cluster members. No significant abundance scatter was found in any of the chemical species studied. Combining the positional, heliocentric distance, and kinematic information we derive, for the first time, the probable orbit of IC 166 within a Galactic model including a rotating boxy bar, and found that it is likely that IC 166 formed in the Galactic disc, supporting its nature as an unremarkable Galactic open cluster with an orbit bound to the Galactic plane.

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.

Bars are common galactic structures in the local universe that play an important role in the secular evolution of galaxies, including the Milky Way. In particular, the velocity distribution of individual stars in our galaxy is useful to shed light on stellar dynamics, and provides information complementary to that inferred from the integrated light of external galaxies. However, since a wide variety of models reproduce the distribution of velocity and the velocity dispersion observed in the Milky Way, we look for signatures of the bar on higher-order moments of the line-of-sight velocity ($V_{los}$) distribution. We make use of two different numerical simulations --one that has developed a bar and one that remains nearly axisymmetric-- to compare them with observations in the latest APOGEE data release (SDSS DR14). This comparison reveals three interesting structures that support the notion that the Milky Way is a barred galaxy. A high skewness region found at positive longitudes constrains the orientation angle of the bar, and is incompatible with the orientation of the bar at $\ell=0^\circ$ proposed in previous studies. We also analyse the $V_{los}$ distributions in three regions, and introduce the Hellinger distance to quantify the differences among them. Our results show a strong non-Gaussian distribution both in the data and in the barred model, confirming the qualitative conclusions drawn from the velocity maps. In contrast to earlier work, we conclude it is possible to infer the presence of the bar from the kurtosis distribution.

The nature of the Triangulum-Andromeda (TriAnd) system has been debated since the discovery of this distant, low-latitude Milky Way (MW) overdensity more than a decade ago. Explanations for its origin are either as a halo substructure from the disruption of a dwarf galaxy or a distant extension of the Galactic disk. We test these hypotheses using chemical abundances of a dozen TriAnd members from the Sloan Digital Sky Survey's 14th Data Release of Apache Point Observatory Galactic Evolution Experiment (APOGEE) data to compare to APOGEE abundances of stars with similar metallicity from both the Sagittarius (Sgr) dSph, and the outer MW disk. We find that TriAnd stars are chemically distinct from Sgr across a variety of elements, (C+N), Mg, K, Ca, Mn, and Ni, with a separation in [X/Fe] of about 0.1 to 0.4 dex depending on the element. Instead, the TriAnd stars, with a median metallicity of about -0.8, exhibit chemical abundance ratios similar to those of the lowest metallicity ([Fe/H] ~ -0.7) stars in the outer Galactic disk, and are consistent with expectations of extrapolated chemical gradients in the outer disk of the MW. These results suggest that TriAnd is associated with the MW disk, and, therefore, that the disk extends to this overdensity --- i.e., past a Galactocentric radius of 24 kpc --- albeit vertically perturbed about 7 kpc below the nominal disk midplane in this region of the Galaxy.

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.

MaNGA provides the opportunity to make precise spatially resolved measurements of the IMF slope in galaxies owing to its unique combination of spatial resolution, wavelength coverage and sample size. We derive radial gradients in age, element abundances and IMF slope analysing optical and near-infrared absorption features from stacked spectra out to the half-light radius of 366 early-type galaxies with masses $9.9 - 10.8\;\log M/M_{\odot}$. We find flat gradients in age and [$\alpha$/Fe] ratio, as well as negative gradients in metallicity, consistent with the literature. We further derive significant negative gradients in the [Na/Fe] ratio with galaxy centres being well enhanced in Na abundance by up to 0.5 dex. Finally, we find a gradient in IMF slope with a bottom-heavy IMF in the centre (typical mass excess factor of 1.5) and a Milky Way-type IMF at the half-light radius. This pattern is mass-dependent with the lowest mass galaxies in our sample featuring only a shallow gradient around a Milky Way IMF. Our results imply the local IMF-$\sigma$ relation within galaxies to be even steeper than the global relation and hint towards the local metallicity being the dominating factor behind the IMF variations. We also employ different stellar population models in our analysis and show that a radial IMF gradient is found independently of the stellar population model used. A similar analysis of the Wing-Ford band provides inconsistent results and further evidence of the difficulty in measuring and modelling this particular feature.

In the age of high-resolution spectroscopic stellar surveys of the Milky Way, the number of stars with detailed abundances of multiple elements is rapidly increasing. These elemental abundances are directly influenced by the evolutionary history of the Galaxy, but this can be difficult to interpret without an absolute timeline of the abundance enrichment. We present age-abundance trends for [M/H], [{\alpha}/M], and 17 individual elements using a sample of 721 solar neighbourhood Hipparcos red giant stars observed by APOGEE. These age trends are determined through a Bayesian hierarchical modelling method presented by Feuillet et al. (2016). We confirm that the [{\alpha}/M]- age relation in the solar neighbourhood is steep and relatively narrow (0.20 dex age dispersion), as are the [O/M]- and [Mg/M]-age relations. The age trend of [C/N] is steep and smooth, consistent with stellar evolution. The [M/H]-age relation has a mean age dispersion of 0.28 dex and a complex overall structure. The oldest stars in our sample are those with the lowest and highest metallicities, while the youngest stars are those with solar metallicity. These results provide strong constraints on theoretical models of Galactic chemical evolution (GCE). We compare them to the predictions of one-zone GCE mod- els and multi-zone mixtures, both analytic and numerical. These comparisons support the hypothesis that the solar neighbourhood is composed of stars born at a range of Galactocentric radii, and that the most metal-rich stars likely migrated from a region with earlier and more rapid star formation such as the inner Galaxy.

In our previous work, we found that only two scenarios are capable of reproducing the observed integrated mass-metallicity relations for the gas and stellar components of local star-forming galaxies simultaneously. One scenario invokes a time-dependent metal outflow loading factor with stronger outflows at early times. The other scenario uses a time-dependent IMF slope with a steeper IMF at early times. In this work, we extend our study to investigate the radial profile of gas and stellar metallicity in local star-forming galaxies using spatially resolved spectroscopic data from the SDSS-IV MaNGA survey. We find that most galaxies show negative gradients in both gas and stellar metallicity with steeper gradients in stellar metallicity. The stellar metallicity gradient tend to be mass dependent with steeper gradients in more massive galaxies while no clear mass dependence is found for the gas metallicity gradient. Then we compare the observations with the predictions from a chemical evolution model of the radial profiles of gas and stellar metallicities. We confirm that the two scenarios proposed in our previous work are also required to explain the metallicity gradients. Based on these two scenarios we successfully reproduce the radial profiles of gas metallicity, stellar metallicity, stellar mass surface density, and star formation rate (SFR) surface density simultaneously. The origin of the negative gradient in stellar metallicity turns out to be driven by either radially dependent metal outflow or IMF slope. In contrast, the radial dependence of the gas metallicity is less constrained because of the degeneracy in model parameters.

In this work we make use of DR14 APOGEE spectroscopic data to study a sample of 92 known OB stars. We developed a near-infrared semi-empirical spectral classification method that was successfully used in case of four new exemplars, previously classified as later B-type stars. Our results agree well with those determined independently from ECHELLE optical spectra, being in line with the spectral types derived from the "canonical" MK blue optical system. This confirms that the APOGEE spectrograph can also be used as a powerful tool in surveys aiming to unveil and study large number of moderately and highly obscured OB stars still hidden in the Galaxy.

Understanding the formation and evolution of our Galaxy requires accurate distances, ages and chemistry for large populations of field stars. Here we present several updates to our spectro-photometric distance code, that can now also be used to estimate ages, masses, and extinctions for individual stars. Given a set of measured spectro-photometric parameters, we calculate the posterior probability distribution over a given grid of stellar evolutionary models, using flexible Galactic stellar-population priors. The code (called {\tt StarHorse}) can acommodate different observational datasets, prior options, partially missing data, and the inclusion of parallax information into the estimated probabilities. We validate the code using a variety of simulated stars as well as real stars with parameters determined from asteroseismology, eclipsing binaries, and isochrone fits to star clusters. Our main goal in this validation process is to test the applicability of the code to field stars with known {\it Gaia}-like parallaxes. The typical internal precision (obtained from realistic simulations of an APOGEE+Gaia-like sample) are $\simeq 8\%$ in distance, $\simeq 20\%$ in age,$\simeq 6\ %$ in mass, and $\simeq 0.04$ mag in $A_V$. The median external precision (derived from comparisons with earlier work for real stars) varies with the sample used, but lies in the range of $\simeq [0,2]\%$ for distances, $\simeq [12,31]\%$ for ages, $\simeq [4,12]\%$ for masses, and $\simeq 0.07$ mag for $A_V$. We provide StarHorse distances and extinctions for the APOGEE DR14, RAVE DR5, GES DR3 and GALAH DR1 catalogues.

Metal-poor globular clusters (GCs) exhibit intriguing Al-Mg anti-correlations and possible Si-Al correlations, which are important clues to decipher the multiple-population phenomenon. NGC 5053 is one of the most metal-poor GCs in the nearby Universe, and has been suggested to be associated with the Sagittarius (Sgr) dwarf galaxy, due to its similarity in location and radial velocity with one of the Sgr arms. In this work, we simulate the orbit of NGC 5053, and argue against a physical connection between Sgr and NGC 5053. On the other hand, the Mg, Al, and Si spectral lines, which are difficult to detect in the optical spectra of NGC 5053 stars, have been detected in the near-infrared APOGEE spectra. We use three different sets of stellar parameters and codes to derive the Mg, Al, and Si abundances. Regardless of which method is adopted, we see a large Al variation, and a substantial Si spread. Along with NGC 5053, metal-poor GCs exhibit different Mg, Al, and Si variations. Moreover, NGC 5053 has the lowest cluster mass among the GCs that have been identified to exhibit an observable Si spread until now.

The Hercules stream is a group of co-moving stars in the Solar neighbourhood, which can potentially be explained as a signature of either the outer Lindblad resonance (OLR) of a fast Galactic bar or the corotation resonance of a slower bar. In either case, the feature should be present over a large area of the disc. With the recent commissioning of the APOGEE-2 Southern spectrograph we can search for the Hercules stream at $(l,b)=(270^\circ,0)$, a direction in which the Hercules stream, if caused by the bar's OLR, would be strong enough to be detected using only the line-of-sight velocities. We clearly detect a narrow, Hercules-like feature in the data that can be traced from the solar neighbourhood to a distance of about 4 kpc. The detected feature matches well the line-of-sight velocity distribution from the fast-bar (OLR) model. Confronting the data with a model where the Hercules stream is caused by the corotation resonance of a slower bar leads to a poorer match, as the corotation model does not predict clearly separated modes, possibly because the slow-bar model is too hot.

We carried out an optical polarimetric study in the direction of the RCW95 star forming region in order to probe the sky-projected magnetic field structure by using the distribution of linear polarization segments which seem to be well aligned with the more extended cloud component. A mean polarization angle of $\theta=49.8^o\pm7.7^o$ was derived. Through the spectral dependence analysis of polarization it was possible to obtain the total-to-selective extinction ratio ($R_V$) by fitting the Serkowski function, resulting in a mean value of $R_V=2.93\pm0.47$. The foreground polarization component was estimated and is in agreement with previous studies in this direction of the Galaxy. Further, near-infrared images from Vista Variables in the Via L\'actea (VVV) survey were collected to improve the study of the stellar population associated with the HII region. The Automated Stellar Cluster Analysis (ASteCA) algorithm was employed to derive structural parameters for two clusters in the region, and a set of PAdova and TRieste Stellar Evolution Code (PARSEC) isochrones was superimposed on the decontaminated colour-magnitude diagrams (CMDs) to estimate an age of about 3 Myr for both clusters. Finally, from the near-infrared photometry study combined with spectra obtained with the Ohio State Infrared Imager and Spectrometer (OSIRIS) mounted at the Southern Astrophysics Research Telescope (SOAR) we derived the spectral classification of the main ionizing sources in the clusters associated with IRAS 15408$-$5356 and IRAS 15412$-$5359, both objects classified as O4 V stars.

We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.

A cold high-velocity (HV, $\sim$ 200 km/s) peak was first reported in several Galactic bulge fields based on the APOGEE commissioning observations. Both the existence and the nature of the high-velocity peak are still under debate. Here we revisit this feature with the latest APOGEE DR13 data. We find that most of the low latitude bulge fields display a skewed Gaussian distribution with a HV shoulder. However, only 3 out of 53 fields show distinct high-velocity peaks around 200 km/s. The velocity distribution can be well described by Gauss-Hermite polynomials, except the three fields showing clear HV peaks. We find that the correlation between the skewness parameter ($h_{3}$) and the mean velocity ($\bar{v}$), instead of a distinctive HV peak, is a strong indicator of the bar. It was recently suggested that the HV peak is composed of preferentially young stars. We choose three fields showing clear HV peaks to test this hypothesis using the metallicity, [$\alpha$/M] and [C/N] as age proxies. We find that both young and old stars show HV features. The similarity between the chemical abundances of stars in the HV peaks and the main component indicates that they are not systematically different in terms of chemical abundance or age. In contrast, there are clear differences in chemical space between stars in the Sagittarius dwarf and the bulge stars. The strong HV peaks off-plane are still to be explained properly, and could be different in nature.

We report the peculiar chemical abundance patterns of eleven atypical Milky Way (MW) field red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). These atypical giants exhibit strong Al and N enhancements accompanied by C and Mg depletions, strikingly similar to those observed in the so-called second-generation (SG) stars of globular clusters (GCs). Remarkably, we find low-Mg abundances ([Mg/Fe]$<$0.0) together with strong Al and N overabundances in the majority (5/7) of the metal-rich ([Fe/H]$\gtrsim - 1.0$) sample stars, which is at odds with actual observations of SG stars in Galactic CGs of similar metallicities. This chemical pattern is unique and unprecedented among MW stars, posing urgent questions about its origin. These atypical stars could be former SG stars of dissolved GCs formed with intrinsically lower abundances of Mg and enriched Al (subsequently self-polluted by massive AGB stars) or the result of exotic binary systems. We speculate that the stars Mg-deficiency as well as the orbital properties suggest that they could have an extragalactic origin. This discovery should guide future dedicated spectroscopic searches of atypical stellar chemical patterns in our Galaxy; a fundamental step forward to understand the Galactic formation and evolution.

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

We utilize elemental-abundance information for Galactic red giant stars in five open clusters (NGC 7789, NGC 6819, M67, NGC 188, and NGC 6791) from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR13 dataset to age-date the chemical evolution of the high- and low-$\alpha$ element sequences of the Milky Way. Key to this time-stamping is the cluster NGC 6791, whose stellar members have mean abundances that place it in the high-$\alpha$, high-[Fe/H] region of the [$\alpha$/Fe]-[Fe/H] plane. Based on the cluster's age ($\sim 8$ Gyr), Galactocentric radius, and height above the Galactic plane, as well as comparable chemistry reported for APOGEE stars in Baade's Window, we suggest that the two most likely origins for NGC 6791 are as an original part of the thick-disk, or as a former member of the Galactic bulge. Moreover, because NGC 6791 lies at the \textit{high metallicity end} ([Fe/H] $\sim 0.4$) of the high-$\alpha$ sequence, the age of NGC 6791 places a limit on the \textit{youngest age} of stars in the high-metallicity, high-$\alpha$ sequence for the cluster's parent population (i.e., either the bulge or the disk). In a similar way, we can also use the age and chemistry of NGC 188 to set a limit of $\sim 7$ Gyr on the \textit{oldest age} of the low-$\alpha$ sequence of the Milky Way. Therefore, NGC 6791 and NGC 188 are potentially a pair of star clusters that bracket both the timing and the duration of an important transition point in the chemical history of the Milky Way.

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 the detection of a curved magnetic field in the ring-like shell of the bubble N4, derived from near-infrared polarization of reddened diskless stars located behind this bubble. The magnetic field in the shell is curved and parallel to the ring-like shell, and its strength is estimated to be $\sim120\,\mu$G in the plane of the sky. The magnetic field strength in the shell is significantly enhanced compared to the local field strength. We calculate the mass-to-flux ratio for the submillimeter clumps in the shell and find that they are all magnetically subcritical. Our results demonstrate that the magnetic field strengthens as the interstellar medium is compressed into a shell, and suggest that the magnetic field has the potential to hinder star formation triggered by \ion{H}{2} region expansion.

We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.

Baade's window (BW) is one of the most observed Galactic bulge fields in terms of chemical abundances. Due to its low and homogeneous interstellar absorption it is considered as a calibration field for Galactic bulge studies. In the era of large spectroscopic surveys, calibration fields such as BW are necessary to cross calibrate the stellar parameters and individual abundances of the APOGEE survey. We use the APOGEE BW stars to derive their metallicity distribution function (MDF) and individual abundances, for $\alpha$- and iron-peak elements of the APOGEE ASPCAP pipeline (DR13), as well as the age distribution for stars in BW. We determine the MDF of APOGEE stars in BW and find a remarkable agreement with that of the Gaia-ESO survey (GES). Both exhibit a clear bimodal distribution. We also find that the Mg-metallicity planes of both surveys agree well, except for the metal-rich part ([Fe/H] >0.1), where APOGEE finds systematically higher Mg abundances with respect to the GES. The ages based on the [C/N] ratio reveal a bimodal age distribution, with a major old population at 10 Gyr, with a decreasing tail towards younger stars. A comparison between APOGEE estimates and stellar parameters, and those determined by other sources reveals detectable systematic offsets, in particular for spectroscopic surface gravity estimates. In general, we find a good agreement between individual abundances of O, Na, Mg, Al, Si, K, Ca, Cr, Mn, Co, and Ni from APOGEE with that of literature values. We have shown that in general APOGEE data show a good agreement in terms of MDF and individual chemical abundances with respect to literature works. Using the [C/N] ration we found a significant fraction of young stars in BW which is in agreement with the model of Haywood et al. (2016).

We explore to which extent stars within Galactic disk open clusters resemble each other in the high-dimensional space of their photospheric element abundances, and contrast this with pairs of field stars. Our analysis is based on abundances for 20 elements, homogeneously derived from APOGEE spectra (with carefully quantified uncertainties, with a median value of $\sim 0.03$ dex). We consider 90 red giant stars in seven open clusters and find that most stars within a cluster have abundances in most elements that are indistinguishable (in a $\chi^2$-sense) from those of the other members, as expected for stellar birth siblings. An analogous analysis among pairs of $>1000$ field stars shows that highly significant abundance differences in the 20-dimensional space can be established for the vast majority of these pairs, and that the APOGEE-based abundance measurements have high discriminating power. However, pairs of field stars whose abundances are indistinguishable even at 0.03~dex precision exist: $\sim 0.3$ percent of all field star pairs, and $\sim 1.0$ percent of field star pairs at the same (solar) metallicity [Fe/H] = $0 \pm 0.02$. Most of these pairs are presumably not birth siblings from the same cluster, but rather doppelganger. Our analysis implies that 'chemical tagging' in the strict sense, identifying birth siblings for typical disk stars through their abundance similarity alone, will not work with such data. However, our approach shows that abundances have extremely valuable information for probabilistic chemo-orbital modeling and combined with velocities, we have identified new cluster members from the field.

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

We study the internal gradients of stellar population properties within $1.5\;R_{\rm e}$ for a representative sample of 721 galaxies with stellar masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV MaNGA IFU survey. Through the use of our full spectral fitting code FIREFLY, we derive light and mass-weighted stellar population properties and their radial gradients, as well as full star formation and metal enrichment histories. We also quanfify the impact that different stellar population models and full spectral fitting routines have on the derived stellar population properties, and the radial gradient measurements. In our analysis, we find that age gradients tend to be shallow for both early-type and late-type galaxies. {\em Mass-weighted} age gradients of early-types are positive ($\sim 0.09\; {\rm dex}/R_{\rm e}$) pointing to "outside-in" progression of star formation, while late-type galaxies have negative {\em light-weighted} age gradients ($\sim -0.11\; {\rm dex}/R_{\rm e}$), suggesting an "inside-out" formation of discs. We detect negative metallicity gradients in both early and late-type galaxies, but these are significantly steeper in late-types, suggesting that radial dependence of chemical enrichment processes and the effect of gas inflow and metal transport are far more pronounced in discs. Metallicity gradients of both morphological classes correlate with galaxy mass, with negative metallicity gradients becoming steeper with increasing galaxy mass. The correlation with mass is stronger for late-type galaxies, with a slope of $d(\nabla [Z/H])/d(\log M)\sim -0.2\pm 0.05\;$, compared to $d(\nabla [Z/H])/d(\log M)\sim -0.05\pm 0.05\;$ for early-types. This result suggests that the merger history plays a relatively small role in shaping metallicity gradients of galaxies.

We study the internal radial gradients of stellar population properties within $1.5\;R_{\rm e}$ and analyse the impact of galaxy environment. We use a representative sample of 721 galaxies with masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV survey MaNGA. We split this sample by morphology into early-type and late-type galaxies. Using the full spectral fitting code FIREFLY, we derive the light and mass-weighted stellar population properties age and metallicity, and calculate the gradients of these properties. We use three independent methods to quantify galaxy environment, namely the $N^{th}$ nearest neighbour, the tidal strength parameter $Q$ and distinguish between central and satellite galaxies. In our analysis, we find that early-type galaxies generally exhibit shallow light-weighted age gradients in agreement with the literature and mass-weighted median age gradients tend to be slightly positive. Late-type galaxies, instead, have negative light-weighted age gradients. We detect negative metallicity gradients in both early and late-type galaxies that correlate with galaxy mass, with the gradients being steeper and the correlation with mass being stronger in late-types. We find, however, that stellar population gradients, for both morphological classifications, have no significant correlation with galaxy environment for all three characterisations of environment. Our results suggest that galaxy mass is the main driver of stellar population gradients in both early and late-type galaxies, and any environmental dependence, if present at all, must be very subtle.

We report chemical abundances obtained by SDSS-III/APOGEE for giant stars in five globular clusters located within 2.2 kpc of the Galactic centre. We detect the presence of multiple stellar populations in four of those clusters (NGC 6553, NGC 6528, Terzan 5, and Palomar 6) and find strong evidence for their presence in NGC 6522. All clusters present a significant spread in the abundances of N, C, Na, and Al, with the usual correlations and anti-correlations between various abundances seen in other globular clusters. Our results provide important quantitative constraints on theoretical models for self-enrichment of globular clusters, by testing their predictions for the dependence of yields of elements such as Na, N, C, and Al on metallicity. They also confirm that, under the assumption that field N-rich stars originate from globular cluster destruction, they can be used as tracers of their parental systems in the high- metallicity regime.

We study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $\omega$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

The nature of warm, ionized gas outside of galaxies may illuminate several key galaxy evolutionary processes. A serendipitous observation by the MaNGA survey has revealed a large, asymmetric H$\alpha$ complex with no optical counterpart that extends $\approx8"$ ($\approx6.3$ kpc) beyond the effective radius of a dusty, starbursting galaxy. This H$\alpha$ extension is approximately three times the effective radius of the host galaxy and displays a tail-like morphology. We analyze its gas-phase metallicities, gaseous kinematics, and emission-line ratios, and discuss whether this H$\alpha$ extension could be diffuse ionized gas, a gas accretion event, or something else. We find that this warm, ionized gas structure is most consistent with gas accretion through recycled wind material, which could be an important process that regulates the low-mass end of the galaxy stellar mass function.

By means of the fossil record method implemented through Pipe3D, we reconstruct the global and radial stellar mass growth histories (MGHs) of an unprecedentedly large sample of galaxies, ranging from dwarf to giant objects, from the "Mapping Nearby Galaxies at the Apache Point Observatory" survey. We confirm that the main driver of the global MGHs is mass, with more massive galaxies assembling their masses earlier (downsizing), though for a given mass, the global MGHs segregate by color, specific star formation rate (sSFR), and morphological type. From the inferred radial mean MGHs, we find that at the late evolutionary stages (or for fractions of assembled mass larger than ~ 80%), the innermost regions formed stars on average earlier than the outermost ones (inside-out). At earlier epochs, when the age resolution of the method becomes poor, the mass assembly seems to be spatially homogeneous or even in the outside-in mode, specially for the red/quiescent/early-type galaxies. The innermost MGHs are in general more regular (less scatter around the mean) than the outermost ones. For dwarf and low-mass galaxies, we do not find evidence of an outside-in formation mode; instead their radial MGHs are very diverse most of the time, with periods of outside- in and inside-out modes (or strong radial migration), suggesting this an episodic SF history. Blue/star-forming/late-type galaxies present on average a significantly more pronounced inside-out formation mode than red/quiescent/early-type galaxies, independently of mass. We discuss our results in the light of the processes of galaxy formation, quenching, and radial migration. We discuss also on the uncertainties and biases of the fossil record method and how they could affect our results.

Using kinematic maps from the Sloan Digital Sky Survey (SDSS) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we reveal that the majority of low-mass quenched galaxies exhibit coherent rotation in their stellar kinematics. Our sample includes all 39 quenched low-mass galaxies observed in the first year of MaNGA. The galaxies are selected with $M_{r} > -19.1$, stellar masses $10^{9}$ M$_{\odot} < M_{\star} < 5\times10^{9}$ M$_{\odot}$, EW$_{H\alpha} <2$ \AA, and all have red colours $(u-r)>1.9$. They lie on the size-magnitude and $\sigma$-luminosity relations for previously studied dwarf galaxies. Just six ($15\pm5.7$ per cent) are found to have rotation speeds $v_{e,rot} < 15$ km s$^{-1}$ at $\sim1$ $R_{e}$, and may be dominated by pressure support at all radii. Two galaxies in our sample have kinematically distinct cores in their stellar component, likely the result of accretion. Six contain ionised gas despite not hosting ongoing star formation, and this gas is typically kinematically misaligned from their stellar component. This is the first large-scale Integral Field Unit (IFU) study of low mass galaxies selected without bias against low-density environments. Nevertheless, we find the majority of these galaxies are within $\sim1.5$ Mpc of a bright neighbour ($M_{K} < -23$; or M$_{\star} > 5\times10^{10}$ M$_{\odot}$), supporting the hypothesis that galaxy-galaxy or galaxy-group interactions quench star formation in low-mass galaxies. The local bright galaxy density for our sample is $\rho_{proj} = 8.2\pm2.0$ Mpc$^{-2}$, compared to $\rho_{proj} = 2.1\pm0.4$ Mpc$^{-2}$ for a star forming comparison sample, confirming that the quenched low mass galaxies are preferentially found in higher density environments.

Quiescent galaxies with little or no ongoing star formation dominate the galaxy population above $M_{*}\sim 2 \times 10^{10}~M_{\odot}$, where their numbers have increased by a factor of $\sim25$ since $z\sim2$. Once star formation is initially shut down, perhaps during the quasar phase of rapid accretion onto a supermassive black hole, an unknown mechanism must remove or heat subsequently accreted gas from stellar mass loss or mergers that would otherwise cool to form stars. Energy output from a black hole accreting at a low rate has been proposed, but observational evidence for this in the form of expanding hot gas shells is indirect and limited to radio galaxies at the centers of clusters, which are too rare to explain the vast majority of the quiescent population. Here we report bisymmetric emission features co-aligned with strong ionized gas velocity gradients from which we infer the presence of centrally-driven winds in typical quiescent galaxies that host low-luminosity active nuclei. These galaxies are surprisingly common, accounting for as much as $10\%$ of the population at $M_* \sim 2 \times 10^{10}~ M_{\odot}$. In a prototypical example, we calculate that the energy input from the galaxy's low-level active nucleus is capable of driving the observed wind, which contains sufficient mechanical energy to heat ambient, cooler gas (also detected) and thereby suppress star formation.

We study the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA. Making use of Baldwin-Philips-Terlevich diagnostic diagrams we demonstrate the ubiquitous presence of extended (kpc scale) low ionisation emission-line regions (LIERs) in both star forming and quiescent galaxies. In star forming galaxies LIER emission can be associated with diffuse ionised gas, most evident as extra-planar emission in edge-on systems. In addition, we identify two main classes of galaxies displaying LIER emission: `central LIER' (cLIER) galaxies, where central LIER emission is spatially extended, but accompanied by star formation at larger galactocentric distances, and `extended LIER' (eLIER) galaxies, where LIER emission is extended throughout the whole galaxy. In eLIER and cLIER galaxies, LIER emission is associated with radially flat, low H$\alpha$ equivalent width of line emission ($<$ 3 \AA) and stellar population indices demonstrating the lack of young stellar populations, implying that line emission follows tightly the continuum due to the underlying old stellar population. The H$\alpha$ surface brightness radial profiles are always shallower than $\rm 1/r^{2}$ and the line ratio [OIII]$\lambda$5007/[OII]$\lambda$3727,29 (a tracer of the ionisation parameter of the gas) shows a flat gradient. This combined evidence strongly supports the scenario in which LIER emission is not due to a central point source but to diffuse stellar sources, the most likely candidates being hot, evolved (post-asymptotic giant branch) stars. Shocks are observed to play a significant role in the ionisation of the gas only in rare merging and interacting systems.

We present results from the largest systematic investigation of broad absorption line (BAL) acceleration to date. We use spectra of 140 quasars from three Sloan Digital Sky Survey programs to search for global velocity offsets in BALs over timescales of ~2.5-5.5 years in the quasar rest frame. We carefully select acceleration candidates by requiring monolithic velocity shifts over the entire BAL trough, avoiding BALs with velocity shifts that might be caused by profile variability. The CIV BALs of two quasars show velocity shifts consistent with the expected signatures of BAL acceleration, and the BAL of one quasar shows a velocity-shift signature of deceleration. In our two acceleration candidates, we see evidence that the magnitude of the acceleration is not constant over time; the magnitudes of the change in acceleration for both acceleration candidates are difficult to produce with a standard disk-wind model or via geometric projection effects. We measure upper limits to acceleration and deceleration for 76 additional BAL troughs and find that the majority of BALs are stable to within about 3% of their mean velocities. The lack of widespread acceleration/deceleration could indicate that the gas producing most BALs is located at large radii from the central black hole and/or is not currently strongly interacting with ambient material within the host galaxy along our line of sight.

In this work, we present a spectroscopic study of very massive stars found outside the center of the massive stellar cluster NGC3603. From the analysis of SOAR spectroscopic data and related optical-NIR photometry, we confirm the existence of several very massive stars in the periphery of NGC 3603. The first group of objects (MTT58, WR42e and RFS7) is compound by three new Galactic exemplars of the OIf*/WN type, all of them with probable initial masses well above 100 Msun and estimated ages of about 1 Myr. Based on Goodman blue-optical spectrum of MTT68, we can confirm the previous finding in the NIR of the only other Galactic exemplar (besides HD93129A) of the O2If* type known to date. Based on its position relative to a set of theoretical isochrons in a Hertzprung-Russel diagram, we concluded that the new O2If* star could be one of the most massive (150 Msun) and luminous (Mv=-7.3) O-star in the Galaxy. Also, another remarkable result is the discovery of a new O2V star (MTT31) that is the first exemplar of the class so far identified in the Milk Way. From its position in the Hertzprung-Russel diagram it is found that this new star probably had an initial mass of 80 Msun, as well as an absolute magnitude Mv=-6.0 corresponding to a luminosity similar to other known O2V stars in the LMC. We also communicate the discovery of a new Galactic O3.5If* star (RFS8) which case is quite intriguing. Indeed, It is located far to the south of the NGC 3603 center, in apparent isolation at a large radial projected linear distance of 62 pc. Its derived luminosity is similar to that of the other O3.5If* (Sh18) found in the NGC 3603's innermost region. The fact that a such high mass star is observed far isolated in the field led us to speculate that perhaps it could have been expelled from the innermost parts of the complex by a close fly-by dynamical encounter with a very massive hard binary system.

We present Pipe3D, an analysis pipeline based on the FIT3D fitting tool, devel- oped to explore the properties of the stellar populations and ionized gas of Integral Field Spectroscopy data. Pipe3D was created to provide with coherent, simple to distribute, and comparable dataproducts, independently of the origin of the data, focused on the data of the most recent IFU surveys (e.g., CALIFA, MaNGA, and SAMI), and the last generation IFS instruments (e.g., MUSE). Along this article we describe the different steps involved in the analysis of the data, illustrating them by showing the dataproducts derived for NGC 2916, observed by CALIFA and P-MaNGA. As a practical use of the pipeline we present the complete set of dataproducts derived for the 200 datacubes that comprises the V500 setup of the CALIFA Data Release 2 (DR2), making them freely available through the network (ftp://ftp.caha.es/CALIFA/dataproducts/DR2/Pipe3D). Finally, we explore the hypothesis that the properties of the stellar populations and ionized gas of galaxies at the effective radius are representative of the overall average ones, finding that this is indeed the case.

Physical properties were derived for the candidate open cluster La Serena 94, recently unveiled by the VVV collaboration. Thanks to the exquisite angular resolution provided by GeMS/GSAOI, we could characterize this system in detail, for the first time, with deep photometry in JHK$_{s}$ - bands. Decontaminated JHK$_{s}$ diagrams reach about 5 mag below the cluster turnoff in H. The locus of red clump giants in the colour - colour diagram, together with an extinction law, was used to obtain an average extinction of $A_V =14.18 \pm 0.71$. The same stars were considered as standard - candles to derive the cluster distance, $8.5 \pm 1.0$ kpc. Isochrones were matched to the cluster colour - magnitude diagrams to determine its age, $\log{t(yr)}=9.12\pm 0.06$, and metallicity, $Z=0.02\pm0.01$. A core radius of $r_{c}=0.51\pm 0.04$ pc was found by fitting King models to the radial density profile. By adding up the visible stellar mass to an extrapolated mass function, the cluster mass was estimated as $M=(2.65\pm0.57) \times 10^3$ M$_{\odot}$, consistent with an integrated magnitude of $M_{K}=-5.82\pm0.16$ and a tidal radius of $r_{t}=17.2\pm2.1$ pc. The overall characteristics of La Serena 94 confirm that it is an old open cluster located in the Crux spiral arm towards the fourth Galactic quadrant and distant $7.30\pm 0.49$ kpc from the Galactic centre. The cluster distorted structure, mass segregation and age indicate that it is a dynamically evolved stellar system.

VISTA variables in the Via Lactea is an ESO Public survey dedicated to scan the bulge and an adjacent portion of the Galactic disk in the fourth quadrant using the VISTA telescope and the near-infrared camera VIRCAM. One of the leading goals of the VVV survey is to contribute to the knowledge of the star cluster population of the Milky Way. To improve the census of the Galactic star clusters, we performed a systematic scan of the JHKs images of the Galactic plane section of the VVV survey. Our detection procedure is based on a combination of superficial density maps and visual inspection of promising features in the NIR images. The material examined are color-composite images corresponding to the DR1 of VVV. We report the discovery of 493 new star cluster candidates. The analysis of the spatial distribution show that the clusters are very concentrated in he Galactic plane, presenting some local maxima around the position of large star-forming complexes, such as G305, RCW 95, and RCW 106. The vast majority of the cluster candidates are quite compact and generally surrounded by bright and/or dark nebulosities. IRAS point sources are associated with 59% of the sample, while 88% are associated with MSX point sources. GLIMPSE 8 mum images of the cluster candidates show a variety of morphologies, with 292 clusters dominated by knotty sources, while 361 clusters show some kind of nebulosity. Spatial cross-correlation with young stellar objects, masers, and extended green-object catalogs suggest that a large sample of the new cluster candidates are extremely young. In particular, 104 star clusters associated to methanol masers are excellent candidates for ongoing massive star formation. Also, there is a special set of sixteen cluster candidates that present clear signspot of star-forming activity having associated simultaneosly dark nebulae, young stellar objects, EGOs, and masers.

The VISTA Variables in the V\'ia L\'actea (VVV) survey is one of six public ESO surveys, and is now in its 4th year of observing. Although far from being complete, the VVV survey has already delivered many results, some directly connected to the intended science goals (detection of variables stars, microlensing events, new star clusters), others concerning more exotic objects, e.g. novae. Now, at the end of the fourth observing period, and comprising roughly 50% of the proposed observations, the actual status of the survey, as well some of the results based on the VVV data, are presented.

Sh 2-29 is a conspicuous star-forming region marked by the presence of massive embedded stars as well as several notable interstellar structures. In this research, our goals were to determine the role of magnetic fields and to study the size distribution of interstellar dust particles within this turbulent environment. We have used a set of optical and near-infrared polarimetric data obtained at OPD/LNA (Brazil) and CTIO (Chile), correlated with extinction maps, 2MASS data and images from DSS and Spitzer. The region's most striking feature is a swept out interstellar cavity whose polarimetric maps indicate that magnetic field lines were dragged outwards, pilling up along its borders. This led to a higher magnetic strength value ($\approx400\,\mu$G) and an abrupt increase in polarization degree, probably due to an enhancement in alignment efficiency. Furthermore, dense cloud fragmentations with peak $A_{V}$ between 20 and 37 mag were probably triggered by its expansion. The presence of $24\,\mu$m point-like sources indicates possible newborn stars inside this dense environment. A statistical analysis of the angular dispersion function revealed areas where field lines are aligned in a well-ordered pattern, seemingly due to compression effects from the HII region expansion. Finally, Serkowski function fits were used to study the ratio of the total-to-selective extinction, reveling a dual population of anomalous grain particles sizes. This trend suggests that both effects of coagulation and fragmentation of interstellar grains are present in the region.

In this letter we communicate the identification of a new Galactic O2If* star (MTT 68) isolated at a projected linear distance of 3 pc from the centre of the star-burst cluster NGC 3603. From its optical photometry I computed a bolometric luminosity M_Bol = -10.7, which corresponds to a total stellar luminosity of 1.5 x 10^{6} L_Sun. It was found an interesting similarity between MTT 68 and the well known multiple system HD 93129. From Hubble Space Telescope F656N images of the NGC 3603 field, it was found that MTT 68 is actually a visual binary system with an angular separation of 0.38 arcsec, which corresponds to a projected (minimum) linear distance of r_(A-B) = 1.4 x 10^{-2} pc. This value is similar to that for the HD 93129A (O2If*) and HD 93129B (O3.5) pair, r_(A-B) = 3.0 x 10^{-2} pc. On the other hand, HD93129A has a third closer companion named HD 93129Ab (O3.5) at only 0.053 arcsec, and taking into account that the X-ray to total stellar luminosity ratio for the MTT 68 system (L_X/L_Bol 10^{-5}) is about two orders of magnitude above the canonical value expected for single stars, I suspect that the MTT 68 system probably hosts another massive companion not resolved by the HST archive images.

The new multi-epoch near-infrared VVV survey (VISTA Variables in the Via Lactea) is sampling 562 sq. deg of the Galactic bulge and adjacent regions of the disk. Accurate astrometry established for the region surveyed allows the VVV data to be merged with overlapping surveys (e.g., GLIMPSE, WISE, 2MASS, etc.), thereby enabling the construction of longer baseline spectral energy distributions for astronomical targets. However, in order to maximize use of the VVV data, a set of transformation equations are required to place the VVV JHKs photometry onto the 2MASS system. The impetus for this work is to develop those transformations via a comparison of 2MASS targets in 152 VVV fields sampling the Galactic disk. The transformation coefficients derived exhibit a reliance on variables such as extinction. The transformed data were subsequently employed to establish a mean reddening law of E_{J-H}/E_{H-Ks}=2.13 +/- 0.04, which is the most precise determination to date and merely emphasizes the pertinence of the VVV data for determining such important parameters.

In this paper I report the discovery of an O2If*/WN6 star probably still partially embedded in its parental cocoon in the star-burst cluster NGC 3603. From the observed size of the associated compact Hii region, it was possible to derive a probable dynamic age of no more than 600,000 years. Using the computed visual extinction value Av ~ 6 magnitudes, an absolute visual magnitude Mv =-5.7 mag is obtained, which for the assumed heliocentric distance of 7.6 kpc results in a bolometric luminosity of ~ 8x10^5 Lsun. Also from the V magnitude and the V-I color of the new star, and previous models for NGC3603's massive star population, we estimate its mass for the binary (O2If*/WN6 + O3If) and the single-star case (O2If*/WN6). In the former, it was found that the initial mass of each component possibly exceeded 80 Msun and 40 Msun, while in the latter MTT 58's initial mass possibly was in excess of 100 Msun.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

In this work I report the discovery of a new Galactic O2 If*/WN6 star, a rare member of the extremely massive hydrogen core-burning group of stars that due its high intrinsic luminosity (close to the Eddington limit), possess an emission-line spectrum at the beginning of their main-sequence evolution, mimicking the spectral appearance of classical WR stars. The new star is named WR42e and is found in isolation at 2.7 arcmin (about 6 pc) from the core of the star-burst cluster NGC3603. From the computed E(B-V) color excess and observed visual magnitude it was possible to estimate its absolute visual magnitude as MV =-6.3 mag, which is a value similar to those obtained by other researchers for stars of similar spectral type both, in the Galaxy and in the Large Magellanic Cloud. Considering the derived absolute visual magnitude, we computed a bolometric stellar luminosity of about 3.2x106 Lsun. Finally, the mass of the new O2If*/WN6 star was estimated by comparing its observed magnitudes and colors with those of other probable NGC3603 cluster members, founding that the WR42e initial mass possibly exceeds 100 Msun.

We report a study of the stellar content of the Near-infrared cluster [DBS2003]\,157 embedded in the extended H\,{\sc ii} region GAL\,331.31-00.34, which is associated with the IRAS source 16085-5138. $JHK$ photometry was carried out in order to identify potential ionizing candidates, and the follow-up NIR spectroscopy allowed the spectral classification of some sources, including two O-type stars. A combination of NIR photometry and spectroscopy data was used to obtain the distance of these two stars, with the method of spectroscopic parallax: IRS\,298 (O6\,{\sc V}, $3.35 \pm 0.61$\,kpc) and IRS\,339 (O9\,{\sc V}, $3.24 \pm 0.56$\,kpc). Adopting the average distance of $3.29 \pm 0.58$\,kpc and comparing the Lyman continuum luminosity of these stars with that required to account for the radio continuum flux of the H\,{\sc ii} region, we conclude that these two stars are the ionizing sources of GAL\,331.31-00.34. Young stellar objects (YSOs) were searched by using our NIR photometry and MIR data from the GLIMPSE survey. The analysis of NIR and MIR colour-colour diagrams resulted in 47 YSO candidates. The GLIMPSE counterpart of IRAS\,16085-5138, which presents IRAS colour indices compatible with an ultra-compact H\,{\sc ii} region, has been identified. The analysis of its spectral energy distribution between 2 and $100\,\mu$m revealed that this source shows a spectral index $\alpha = 3.6$ between 2 and $25\,\mu$m, which is typical of a YSO immersed in a protostellar envelope. Lower limits to the bolometric luminosity and the mass of the embedded protostar have been estimated as $L=7.7\times10^3L_{\sun}$ and $M = 10\,M_{\sun}$, respectively, which corresponds to a B0--B1\,{\sc V} ZAMS star.

The RCW41 star-forming region is embedded within the Vela Molecular Ridge, hosting a massive stellar cluster surrounded by a conspicuous HII region. Understanding the role of interstellar magnetic fields and studying the newborn stellar population is crucial to build a consistent picture of the physical processes acting on this kind of environment. We have carried out a detailed study of the interstellar polarization toward RCW41, with data from an optical and near-infrared polarimetric survey. Additionally, deep near-infrared images from the NTT 3.5m telescope have been used to study the photometric properties of the embedded young stellar cluster, revealing several YSO's candidates. By using a set of pre-main sequence isochrones, a mean cluster age in the range 2.5 - 5.0 million years was determined, and evidence of sequential star formation were revealed. An abrupt decrease in R-band polarization degree is noticed toward the central ionized area, probably due to low grain alignment efficiency caused by the turbulent environment and/or weak intensity of magnetic fields. The distortion of magnetic field lines exhibit a dual behavior, with the mean orientation outside the area approximately following the borders of the star-forming region, and directed radially toward the cluster inside the ionized area, in agreement with simulations of expanding HII regions. The spectral dependence of polarization allowed a meaningful determination of the total-to-selective extinction ratio by fittings of the Serkowski relation. Furthermore, a large rotation of polarization angle as a function of wavelength is detected toward several embedded stars.

The ESO Public Survey VISTA Variables in the Via Lactea (VVV) started in 2010. VVV targets 562 sq. deg in the Galactic bulge and an adjacent plane region and is expected to run for ~5 years. In this paper we describe the progress of the survey observations in the first observing season, the observing strategy and quality of the data obtained. The observations are carried out on the 4-m VISTA telescope in the ZYJHKs filters. In addition to the multi-band imaging the variability monitoring campaign in the Ks filter has started. Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit. The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing. The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHKs filters taken in the 2010 observing season. The typical image quality is ~0.9-1.0". The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHKs images in the disk area and 90% of the JHKs images in the bulge area. The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge. The first season catalogues contain 1.28x10^8 stellar sources in the bulge and 1.68x10^8 in the disk area detected in at least one of the photometric bands. The combined, multi-band catalogues contain more than 1.63x10^8 stellar sources. About 10% of these are double detections due to overlapping adjacent pointings. These overlapping multiple detections are used to characterise the quality of the data. The images in the JHKs bands extend typically ~4 mag deeper than 2MASS. The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions. The astrometry for Ks=15-18 mag has rms ~35-175 mas.

In this work I communicate the detection of a new Galactic Wolf-Rayet star (WR60a) in Centaurus. The H- and K-band spectra of WR60a, show strong carbon near-infrared emission lines, characteristic of Wolf-Rayet stars of the WC5-7 sub-type. Adopting mean absolute magnitude M$_K$ and mean intrinsic ($J-K_S$) and ($H-K_S$) colours, it was found that WR60a suffer a mean visual extinction of 3.8$\pm$1.3 magnitudes, being located at a probable heliocentric distance of 5.2$\pm$0.8 Kpc, which for the related Galactic longitude (l=312) puts this star probably in the Carina-Sagittarius arm at about 5.9 kpc from the Galactic center. I searched for clusters in the vicinity of WR60a, and in principle found no previously known clusters in a search radius region of several tens arc-minutes. The detection of a well isolated WR star induced us to seek for some still unknown cluster, somewhere in the vicinity of WR60a. From inspection of 5.8$\mu$m and 8.0$\mu$m Spitzer/IRAC GLIMPSE images of the region around the new WR star, it was found strong mid-infrared extended emission at about 13.5 arcmin south-west of WR60a. The study of the the H-K$_S$ colour distribution of point sources associated with the extended emission, reveals the presence of a new Galactic cluster candidate probably formed by at least 85 stars.

In this paper we report the identification of two new Galactic O2 If*/WN6 stars (WR20aa and WR20c), in the outskirt of the massive young stellar cluster Westerlund 2. The morphological similarity between the near-infrared spectra of the new stars with that of WR20a and WR21a (two of the most massive binaries known to date) is remarkable, indicating that probably they are also very massive stars. New optical spectroscopic observations of WR20aa suggest an intermediate O2 If*/WN6 spectral type. Based on a mosaic made from the 3.6 microns Spitzer IRAC images of the region including part of the RCW49 complex, we studied the spatial location of the new emission line stars, finding that WR20aa and WR20c are well displaced from the centre of Westerlund 2, being placed at ~ 36 pc (15.7 arcmin) and ~ 58 pc (25.0 arcmin) respectively, for an assumed heliocentric distance of 8 kpc. Also very remarkably, a radius vector connecting both stars would intercept the Westerlund 2 cluster exactly at the place where its stellar density reaches a maximum. We consequently postulate a scenario in which WR20aa and WR20c had a common origin somewhere in the cluster core, being ejected from their birthplace by dynamical interacion with some other very massive objects, perhaps during some earlier stage of the cluster evolution.

We present a map of the spiral structure of the Galaxy, as traced by molecular CS emission associated with IRAS sources which are believed to be compact HII regions. The CS line velocities are used to determine the kinematic distances of the sources, in order to investigate their distribution in the galactic plane. This allows us to use 870 objects to trace the arms, a number larger than that of previous studies based on classical HII regions. The distance ambiguity of the kinematic distances, when it exists, is solved by different procedures, including the latitude distribution and an analysis of the longitude-velocity diagram. The well defined spiral arms are seen to be confined inside the co-rotation radius, as is often the case in spiral galaxies. We identify a square-shaped sub-structure in the CS map with that predicted by stellar orbits at the 4:1 resonance (4 epicycle oscillations in one turn around the galactic center). The sub-structure is found at the expected radius, based on the known pattern rotation speed and epicycle frequency curve. An inner arm presents an end with strong inward curvature and intense star formation that we tentatively associate with the region where this arm surrounds the extremity of the bar, as seen in many barred galaxies. Finally, a new arm with concave curvature is found in the Sagitta to Cepheus region of the sky.

I report the discovery of two new Galactic Wolf-Rayet stars in Circinus via detection of their C, N and He Near-Infrared emission lines, using ESO-NTT-SOFI archival data. The H- and K-band spectra of WR67a and WR67b, indicate that they are Wolf-Rayet stars of WN6h and WC8 sub-types, respectively. WR67a presents a weak-lined spectrum probably reminiscent of young hydrogen rich main-sequence stars such as WR25 in Car OB1 and HD97950 in NGC3603. Indeed, this conclusion is reinforced by the close morphological match of the WR67a H- and K-band spectra with that for WR21a, a known extremely massive binary system. WR67b is probably a non-dusty WC8 Wolf-Rayet star that has a estimated heliocentric distance of 2.7(0.9) kpc, which for its Galactic coordinates, puts the star probably in the near portion of the Scutum-Centaurus arm.

We describe the public ESO near-IR variability survey (VVV) scanning the Milky Way bulge and an adjacent section of the mid-plane where star formation activity is high. The survey will take 1929 hours of observations with the 4-metre VISTA telescope during five years (2010-2014), covering ~10^9 point sources across an area of 520 deg^2, including 33 known globular clusters and ~350 open clusters. The final product will be a deep near-IR atlas in five passbands (0.9-2.5 microns) and a catalogue of more than 10^6 variable point sources. Unlike single-epoch surveys that, in most cases, only produce 2-D maps, the VVV variable star survey will enable the construction of a 3-D map of the surveyed region using well-understood distance indicators such as RR Lyrae stars, and Cepheids. It will yield important information on the ages of the populations. The observations will be combined with data from MACHO, OGLE, EROS, VST, Spitzer, HST, Chandra, INTEGRAL, WISE, Fermi LAT, XMM-Newton, GAIA and ALMA for a complete understanding of the variable sources in the inner Milky Way. This public survey will provide data available to the whole community and therefore will enable further studies of the history of the Milky Way, its globular cluster evolution, and the population census of the Galactic Bulge and center, as well as the investigations of the star forming regions in the disk. The combined variable star catalogues will have important implications for theoretical investigations of pulsation properties of stars.

This work aims to the study of the Sh 2-307 HII region and related stellar population. Near-infrared imaging and spectroscopic observations in the direction of Sh 2-307 were performed using OSIRIS at SOAR Telescope. From J-, H- and K-band spectra of the brightest source in the cluster, we conclude that it has a near-infrared spectra compatible with that taken for O9v-O9.5v stars. Using the derived spectral type and the respective J, H and K-band photometry, we compute a heliocentric distance of 3.2(0.5) kpc, which for R0 = 8 kpc, puts this cluster at more than 10 kpc from the Galactic centre. From the Brg, H2, and continuum narrow-band images we were able to detect both the NIR counterpart of the associated HII region, as well as, the interface between the ionised and the cool molecular gas. Using the 5 Ghz continuum flux density taken from the PMN catalogue and the Brg narrow band image we estimate that the HII region has a mean diameter of 0.94(0.15) pc, mean electron density of 550 cm-3 and an estimated dynamical age of 1.6 Myears years. The large fraction of sources presenting excess emission at 2micron suggests that the stellar population is very young, with many sources still in the pre-main sequence accreting phase. By the use of theoretical pre-main sequence tracks we derived a cluster mean age of about 2.5 Myears, and from the analyses of the fraction of excess emission sources as a function of their spatial distribution we found evidence for an age spread for the embedded pre-main sequence stellar population. Finally, from the study of the spatial distribution of the low-mass sources relative to the main-cluster source and associated photo-dissociation zones, we conclude that the O-type star probably has been triggering the star formation process in the region.

K-band spectra of young stellar candidates in four southern hemisphere clusters have been obtained with the near-infrared spectrograph GNIRS in Gemini South. The clusters are associated with IRAS sources that have colours characteristic of ultracompact HII regions. Spectral types were obtained by comparison of the observed spectra with those of a NIR library; the results include the spectral classification of nine massive stars and seven objects confirmed as background late-type stars. Two of the studied sources have K-band spectra compatible with those characteristic of very hot stars, as inferred from the presence of Civ, Niii, and Nv emission lines at 2.078 micron, 2.116 micron, and 2.100 micron respectively. One of them, I16177 IRS1, has a K-band spectrum similar to that of Cyg OB2 7, an O3If* supergiant star. The nebular K-band spectrum of the associated UC Hii region shows the s-process [Kriii] and [Seiv] high excitation emission lines, previously identified only in planetary nebula. One young stellar object (YSO) was found in each cluster, associated with either the main IRAS source or a nearby resolved MSX component, confirming the results obtained from previous NIR photometric surveys. The distances to the stars were derived from their spectral types and previously determined JHK magnitudes; they agree well with the values obtained from the kinematic method, except in the case of IRAS15408-5356, for which the spectroscopic distance is about a factor two smaller than the kinematic value.

We report NIR broad and narrow band photometric observations in the direction of the IRAS17149-3916 source that reveal the presence of a young cluster of massive stars embedded in an HII region coincident with RCW121. These observations, together with published radio data, MSX and Spitzer images were used to determine some of the physical parameters of the region. We found 96 cluster member candidates in an area of about 1.5 x 2.0 square arcmin, 30% of them showing excess emission in the NIR. IRS 1, the strongest source in the cluster with an estimated spectral type of O5V-O6V ZAMS based on the color-magnitude diagram, is probably the main ionizing source of the HII region detected at radio wavelengths. Using the integrated Brgamma and the 5 GHz flux densities, we derived a mean visual extinction AV=5.49 magnitudes. From the observed size of the Brgamma extended emission, we calculated the emission measure E and the electron density ne, characteristic of compact HII regions.

We report the discovery of a small cluster of massive stars embedded in a NIR nebula in the direction of the IRAS15411-5352 point source, which is related to the alleged planetary nebula W16-185. The majority of the stars present large NIR excess characteristic of young stellar objects and have bright counterparts in the Spitzer IRAC images; the most luminous star (IRS1) is the NIR counterpart of the IRAS source. We found very strong unresolved Brgamma emission at the IRS1 position and more diluted and extended emission across the continuum nebula. From the sizes and electron volume densities we concluded that they represent ultra-compact and compact HII regions, respectively. Comparing the Brgamma emission with the 7 mm free-free emission, we estimated that the visual extinction ranges between 14 and 20 mag. We found that only one star (IRS1) can provide the number of UV photons necessary to ionize the nebula.

We report the discovery of a young massive stellar cluster embedded in an extended HII region, invisible at optical wavelengths where the extinction is $A_V \approx 28$ magnitudes, associated with the IRAS source 16177-5018. $J, H$ and nb$K$ imaging photometry combined with the $K_S$ 2MASS data show the presence of sources with infrared excess emission at 2.2 $\mu$m, concentrated in an area of about one square parsec around a massive young stellar object identified as the IRAS source. This object has a near-mid infrared spectral index betweem 2.2 and 25 $\mu$m $\alpha({\rm IR}) = d {\rm log}(\lambda F_\lambda)/d {\rm log} \lambda $ =4.78, characteristic of compact H II regions, with bolometric luminosity, inferred from the integrated near to far-infrared flux density of $2.8 \times 10^5 L_\odot$, which corresponds to a ZAMS star of about $42 M_\odot$. From the color-magnitude diagram we were able to classify the majority of the cluster members as reddened massive stars earlier than spectral type B5.