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.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

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

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

The Milky Way provides an ideal laboratory to test our understanding of galaxy evolution, owing to our ability to observe our Galaxy over fine scales. However, connecting the Galaxy to the wider galaxy population remains difficult, due to the challenges posed by our internal perspective and to the different observational techniques employed. Here, we present a sample of galaxies identified as Milky Way Analogs (MWAs) on the basis of their stellar masses and bulge-to-total ratios, observed as part of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We analyse the galaxies in terms of their stellar kinematics and populations as well as their ionised gas contents. We find our sample to contain generally young stellar populations in their outskirts. However, we find a wide range of stellar ages in their central regions, and we detect central AGN-like or composite-like activity in roughly half of the sample galaxies, with the other half consisting of galaxies with central star-forming emission or emission consistent with old stars. We measure gradients in gas metallicity and stellar metallicity that are generally flatter in physical units than those measured for the Milky Way; however, we find far better agreement with the Milky Way when scaling gradients by galaxies' disc scale lengths. From this, we argue much of the discrepancy in metallicity gradients to be due to the relative compactness of the Milky Way, with differences in observing perspective also likely to be a factor.

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.

Analyses of data from spectroscopic and astrometric surveys have led to conflicting results concerning the vertical characteristics of the Milky Way. Ages are often used to provide clarity, but typical uncertainties of $>$ 40\,\% restrict the validity of the inferences made. Using the \textit{Kepler} APOKASC sample for context, we explore the global population trends of two K2 campaign fields (3 and 6), which extend further vertically out of the Galactic plane than APOKASC. We analyse the properties of red giant stars utilising three asteroseismic data analysis methods to cross-check and validate detections. The Bayesian inference tool PARAM is used to determine the stellar masses, radii and ages. Evidence of a pronounced red giant branch bump and an [$\alpha$/Fe] dependence on the position of the red clump is observed from the radii distribution of the K2 fields. Two peaks in the age distribution centred at $\sim$5 and and $\sim$12 Gyr are found using a sample with $\sigma_{\rm{age}}$ $<$ 35\,\%. In a comparison with \textit{Kepler}, we find the older peak to be more prominent for K2. This age bimodality is also observed based on a chemical selection of low- ($\leq$ 0.1) and high- ($>$ 0.1) [$\alpha$/Fe] stars. As a function of vertical distance from the Galactic mid-plane ($|Z|$), the age distribution shows a transition from a young to old stellar population with increasing $|Z|$ for the K2 fields. Further coverage of campaign targets with high resolution spectroscopy is required to increase the yield of precise ages achievable with asteroseismology.

Within the standard model of hierarchical galaxy formation in a {\Lambda}CDM Universe, the environment of galaxies is expected to play a key role in driving galaxy formation and evolution. In this paper we investigate whether and how the gas metallicity and the star formation surface density ({\Sigma}_SFR) depend on galaxy environment. To this end we analyse a sample of 1162 local, star-forming galaxies from the galaxy survey Mapping Nearby Galaxies at APO (MaNGA). Generally, both parameters do not show any significant dependence on environment. However, in agreement with previous studies, we find that low-mass satellite galaxies are an exception to this rule. The gas metallicity in these objects increases while their {\Sigma}SFR decreases slightly with environmental density. The present analysis of MaNGA data allows us to extend this to spatially resolved properties. Our study reveals that the gas metallicity gradients of low-mass satellites flatten and their {\Sigma}SFR gradients steepen with increasing environmental density. By extensively exploring a chemical evolution model, we identify two scenarios that are able to explain this pattern: metal-enriched gas accretion or pristine gas inflow with varying accretion timescales. The latter scenario better matches the observed {\Sigma}SFR gradients, and is therefore our preferred solution. In this model, a shorter gas accretion timescale at larger radii is required. This suggests that 'outside-in quenching' governs the star formation processes of low-mass satellite galaxies in dense environments.

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.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

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.

We analyze the stellar initial mass functions (IMF) of a large sample of early type galaxies (ETGs) provided by MaNGA. The large number of IFU spectra of individual galaxies provide high signal-to-noise composite spectra that are essential for constraining IMF and to investigate possible radial gradients of the IMF within individual galaxies. The large sample of ETGs also make it possible to study how the IMF shape depends on various properties of galaxies. We adopt a novel approach to study IMF variations in ETGs, use Bayesian inferences based on full spectrum fitting. The Bayesian method provides a statistically rigorous way to explore potential degeneracy in spectrum fitting, and to distinguish different IMF models with Bayesian evidence. We find that the IMF slope depends systematically on galaxy velocity dispersion, in that galaxies of higher velocity dispersion prefer a more bottom-heavy IMF, but the dependence is almost entirely due to the change of metallicity, $Z$, with velocity dispersion. The IMF shape also depends on stellar age, $A$, but the dependence is completely degenerate with that on metallicity through a combination $AZ^{-1.42}$. Using independent age and metallicity estimates we find that the IMF variation is produced by metallicity instead of age. The IMF near the centers of massive ETGs appears more bottom-heavy than that in the outer parts, while a weak opposite trend is seen for low-mass ETGs. Uncertainties produced by star formation history, dust extinction, $\alpha$-element abundance enhancement and noise in the spectra are tested.

By applying spectroscopic decomposition methods to a sample of MaNGA early-type galaxies, we separate out spatially and kinematically distinct stellar populations, allowing us to explore the similarities and differences between galaxy bulges and discs, and how they affect the global properties of the galaxy. We find that the components have interesting variations in their stellar populations, and display different kinematics. Bulges tend to be consistently more metal rich than their disc counterparts, and while the ages of both components are comparable, there is an interesting tail of younger, more metal poor discs. Bulges and discs follow their own distinct kinematic relationships, both on the plane of the stellar spin parameter, lambda_R, and ellipticity, and in the relation between stellar mass and specific angular momentum, j, with the location of the galaxy as a whole on these planes being determined by how much bulge and disc it contains. As a check of the physical significance of the kinematic decompositions, we also dynamically model the individual galaxy components within the global potential of the galaxy. The resulting components exhibit kinematic parameters consistent with those from the spectroscopic decomposition, and though the dynamical modelling suffers from some degeneracies, the bulges and discs display systematically different intrinsic dynamical properties. This work demonstrates the value in considering the individual components of galaxies rather than treating them as a single entity, which neglects information that may be crucial in understanding where, when and how galaxies evolve into the systems we see today.

The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey provides precise chemical abundances of 18 chemical elements for $\sim$ 176,000 red giant stars distributed over much of the Milky Way Galaxy (MW), and includes observations of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). The APOGEE chemical abundance patterns of Sgr have revealed that it is chemically distinct from the MW in most chemical elements. We employ a \emph{k}-means clustering algorithm to 6-dimensional chemical space defined by [(C+N)/Fe], [O/Fe], [Mg/Fe], [Al/Fe], [Mn/Fe], and [Ni/Fe] to identify 62 MW stars in the APOGEE sample that have Sgr-like chemical abundances. Of the 62 stars, 35 have \emph{Gaia} kinematics and positions consistent with those predicted by \emph{N}-body simulations of the Sgr stream, and are likely stars that have been stripped from Sgr during the last two pericenter passages ($<$ 2 Gyr ago). Another 20 of the 62 stars exhibit chemical abundances indistinguishable from the Sgr stream stars, but are on highly eccentric orbits with median $r_{\rm apo} \sim $ 25 kpc. These stars are likely the `accreted' halo population thought to be the result of a separate merger with the MW 8-11 Gyr ago. We also find one hypervelocity star candidate. We conclude that Sgr was enriched to [Fe/H] $\sim$ -0.2 before its most recent pericenter passage. If the `accreted halo' population is from one major accretion event, then this progenitor galaxy was enriched to at least [Fe/H] $\sim$ -0.6, and had a similar star formation history to Sgr before merging.

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 map the trends of elemental abundance ratios across the Galactic disk, spanning R = 3-15 kpc and midplane distance |Z|= 0-2 kpc, for 15 elements in a sample of 20,485 stars measured by the SDSS/APOGEE survey (O, Na, Mg, Al, Si, P, S, K, Ca, V, Cr, Mn, Fe, Co, Ni). Adopting Mg rather than Fe as our reference element, and separating stars into two populations based on [Fe/Mg], we find that the median trends of [X/Mg] vs. [Mg/H] in each population are nearly independent of location in the Galaxy. The full multi-element cartography can be summarized by combining these nearly universal median sequences with our measured metallicity distribution functions and the relative proportions of the low-[Fe/Mg] (high-alpha) and high-[Fe/Mg] (low-alpha) populations, which depend strongly on R and |Z|. We interpret the median sequences with a semi-empirical "2-process" model that describes both the ratio of core collapse and Type Ia supernova contributions to each element and the metallicity dependence of the supernova yields. These observationally inferred trends can provide strong tests of supernova nucleosynthesis calculations. Our results lead to a relatively simple picture of abundance ratio variations in the Milky Way, in which the trends at any location can be described as the sum of two components with relative contributions that change systematically and smoothly across the Galaxy. Deviations from this picture and future extensions to other elements can provide further insights into the physics of stellar nucleosynthesis and unusual events in the Galaxy's history.

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.

The isolated elliptical (IE) NGC 7796 is accompanied by an interesting early-type dwarf galaxy, named NGC7796-DW1. It exhibits a tidal tail, very boxy isophotes, and multiple nuclei or regions (A, B, and C) that are bluer than the bulk population of the galaxy, indicating a younger age. These properties are suggestive of a dwarf-dwarf merger remnant. We use the Multi-Unit Spectroscopic Explorer (MUSE) at the VLT to investigate NGC 7796-DW1. We extract characteristic spectra to which we apply the STARLIGHT population synthesis software to obtain ages and metallicities of the various population components of the galaxy. The galaxy's main body is old and metal-poor. A surprising result is the extended line emission in the galaxy, forming a ring-like structure with a projected diameter of 2.2 kpc. The line ratios fall into the regime of HII-regions, although OB-stellar populations cannot be identified by spectral signatures. Nucleus A is a relatively old (7 Gyr or older) and metal-poor super star cluster, most probably the nucleus of the dwarf, now displaced. The star-forming regions B and C show younger and distinctly more metal-rich components. The emission line ratios of regions B and C indicate an almost solar oxygen abundance, if compared with radiation models of HII regions. NGC7796-DW1 occupies a particular role in the group of transition-type galaxies with respect to its origin and current evolutionary state, being the companion of an IE. The dwarf-dwarf merger scenario is excluded because of the missing metal-rich merger component. A viable alternative is gas accretion from a reservoir of cold, metal-rich gas. NGC7796 has to provide this gas within its X-ray bright halo. As illustrated by NGC7796-DW1, cold accretion may be a general solution to the problem of extended star formation histories in transition dwarf galaxies. (abridged)

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.

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 fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

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 Monoceros Ring (MRi) structure is an apparent stellar overdensity that has been postulated to entirely encircle the Galactic plane and has been variously described as being due to line-of-sight effects of the Galactic warp and flare or of extragalactic origin (via accretion). Despite being intensely scrutinised in the literature for more than a decade, no studies to-date have been able to definitively uncover its origins. Here we use $N$-body simulations and a genetic algorithm to explore the parameter space for the initial position, orbital parameters and, for the first time, the final location of a satellite progenitor. We fit our models to the latest Pan-STARRS data to determine whether an accretion scenario is capable of producing an in- Plane ring-like structure matching the known parameters of the MRi. Our simulations produce streams that closely match the location, proper motion and kinematics of the MRi structure. However, we are not able to reproduce the mass estimates from earlier studies based on Pan-STARRS data. Furthermore, in contrast with earlier studies our best-fit models are those for progenitors on retrograde orbits. If the MRi was produced by satellite accretion, we find that its progenitor has an initial mass upper limit of ~$10^{10}$M$_\odot$ and the remnant is likely located behind the Galactic bulge, making it diffcult to locate observationally. While our models produce realistic MRi-like structures we cannot definitively conclude that the MRi was produced by the accretion of a satellite galaxy.

We find two chemically distinct populations separated relatively cleanly in the [Fe/H] - [Mg/Fe] plane, but also distinguished in other chemical planes, among metal-poor stars (primarily with metallicities [Fe/H] $< -0.9$) observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and analyzed for Data Release 13 (DR13) of the Sloan Digital Sky Survey. These two stellar populations show the most significant differences in their [X/Fe] ratios for the $\alpha$-elements, C+N, Al, and Ni. In addition to these populations having differing chemistry, the low metallicity high-Mg population (which we denote the HMg population) exhibits a significant net Galactic rotation, whereas the low-Mg population (or LMg population) has halo-like kinematics with little to no net rotation. Based on its properties, the origin of the LMg population is likely as an accreted population of stars. The HMg population shows chemistry (and to an extent kinematics) similar to the thick disk, and is likely associated with $\it in$ $\it situ$ formation. The distinction between the LMg and HMg populations mimics the differences between the populations of low- and high-$\alpha$ halo stars found in previous studies, suggesting that these are samples of the same two populations.

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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

Without the interference of a number of events, galaxies may suffer in crowded environments (e.g., stripping, harassment, strangulation); isolated elliptical galaxies provide a control sample for the study of galaxy formation. We present the study of a sample of isolated ellipticals using imaging from a variety of telescopes, focusing on their globular cluster systems as tracers of their stellar halos. Our main findings are: (a) GC color bimodality is common even in the most isolated systems; (b) the specific frequency of GCs is fairly constant with galaxy mass, without showing an increase towards high-mass systems like in the case of cluster ellipticals; (c) on the other hand, the red fraction of GCs follows the same inverted V shape trend with mass as seen in cluster ellipticals; and (d) the stellar halos show low S\'ersic indices which are consistent with a major merger origin.

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.

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

Multiple populations revealed in globular clusters (GCs) are important windows to the formation and evolution of these stellar systems. The metal-rich GCs in the Galactic bulge are an indispensable part of this picture, but the high optical extinction in this region has prevented extensive research. In this work, we use the high resolution near-infrared (NIR) spectroscopic data from APOGEE to study the chemical abundances of NGC 6553, which is one of the most metal-rich bulge GCs. We identify ten red giants as cluster members using their positions, radial velocities, iron abundances, and NIR photometry. Our sample stars show a mean radial velocity of $-0.14\pm5.47$ km s$^{-1}$, and a mean [Fe/H] of $-0.15\pm 0.05$. We clearly separate two populations of stars in C and N in this GC for the first time. NGC 6553 is the most metal-rich GC where the multiple stellar population phenomenon is found until now. Substantial chemical variations are also found in Na, O, and Al. However, the two populations show similar Si, Ca, and iron-peak element abundances. Therefore, we infer that the CNO, NeNa, and MgAl cycles have been activated, but the MgAl cycle is too weak to show its effect on Mg. Type Ia and Type II supernovae do not seem to have significantly polluted the second generation stars. Comparing with other GC studies, NGC 6553 shows similar chemical variations as other relatively metal-rich GCs. We also confront current GC formation theories with our results, and suggest possible avenues for improvement in the models.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We present the stellar surface mass density {\it vs.} gas metallicity ($\Sigma_*-Z$) relation for more than 500,000 spatially-resolved star-forming resolution elements (spaxels) from a sample of 653 disk galaxies included in the SDSS IV MaNGA survey. We find a tight relation between these local properties, with higher metallicities as the surface density increases. This relation extends over three orders of magnitude in the surface mass density and a factor of four in metallicity. We show that this local relationship can simultaneously reproduce two well-known properties of disk galaxies: their global mass-metallicity relationship {\it and} their radial metallicity gradients. We also find that the $\Sigma_* - Z$ relation is largely independent of the galaxy's total stellar mass and specific star-formation rate (sSFR), except at low stellar mass and high sSFR. These results suggest that in the present-day universe local properties play a key role in determining the gas-phase metallicity in typical disk galaxies.

Wrapping around the Milky Way, the Sagittarius stream is the dominant substructure in the halo. Our statistical selection method has allowed us to identify 106 highly likely members of the Sagittarius stream. Spectroscopic analysis of metallicity and kinematics of all members provides us with a new mapping of the Sagittarius stream. We find correspondence between the velocity distribution of stream stars and those computed for a triaxial model of the Milky Way dark matter halo. The Sagittarius trailing arm exhibits a metallicity gradient, ranging from $-0.59$ dex to $-0.97$ dex over 142$^{\circ}$. This is consistent with the scenario of tidal disruption from a progenitor dwarf galaxy that possessed an internal metallicity gradient. We note high metallicity dispersion in the leading arm, causing a lack of detectable gradient and possibly indicating orbital phase mixing. We additionally report on a potential detection of the Sextans dwarf spheroidal in our data.

As tracers of star formation, galaxy assembly and mass distribution, globular clusters have provided important clues to our understanding of early-type galaxies. But their study has been mostly constrained to galaxy groups and clusters where early-type galaxies dominate, leaving the properties of the globular cluster systems (GCSs) of isolated ellipticals as a mostly uncharted territory. We present Gemini-South/GMOS $g'i'$ observations of five isolated elliptical galaxies: NGC 3962, NGC 2865, IC 4889, NGC 2271 and NGC 4240. Photometry of their GCSs reveals clear color bimodality in three of them, remaining inconclusive for the other two. All the studied GCSs are rather poor with a mean specific frequency $S_N\sim 1.5$, independently of the parent galaxy luminosity. Considering also previous work, it is clear that bimodality and especially the presence of a significant, even dominant, population of blue clusters occurs at even the most isolated systems, casting doubts on a possible accreted origin of metal-poor clusters as suggested by some models. Additionally, we discuss the possible existence of ultra-compact dwarfs around the isolated elliptical NGC 3962.

Previous studies have shown that the kinematics of the field elliptical galaxy NGC 7507 do not necessarily require dark matter. This is troubling because, in the context of LCDM cosmologies, all galaxies should have a large dark matter component. We use penalised pixel fitting software to extract velocities and velocity dispersions from GMOS slit mask spectra. Using Jeans and MONDian modelling we produce best fit models to the velocity dispersion. We find that NGC 7507 has a two component stellar halo, with the outer halo and inner haloes counter rotating. The velocity dispersion profile exhibits an increase at ~70" (~7.9 kpc), reminiscent of several other elliptical galaxies. Our best fit models are those under mild anisotropy which include ~100 times less dark matter than predicted by LCDM, although mildly anisotropic models that are completely dark matter free fit almost equally well. Our MONDian models, both isotropic and anisotropic, systematically fail to reproduce the measured velocity dispersions at almost all radii. The counter rotating outer halo implies a merger remnant, as does the increase in velocity dispersion at ~70". From simulations it seems plausible that the merger that caused the increase in velocity dispersion was a spiral-spiral merger. Our Jeans models are completely consistent with a no dark matter scenario, however, some dark matter can be accommodated, although at much lower concentrations that predicted by LCDM simulations. This indicates that NGC 7507 may be a dark matter free elliptical galaxy. Whether NGC 7507 is completely dark matter free or very dark matter poor, this is at odds with predictions from current LCDM cosmological simulations. It may be possible that the observed velocity dispersions could be reproduced if the galaxy is significantly flattened along the line of sight (e.g. due to rotation), however, invoking this flattening is problematic.

We investigate the globular cluster system of the isolated elliptical NGC 7796, present new photometry of the galaxy, and use published kinematical data to constrain the dark matter content. Deep images in B and R, obtained with the VIsible MultiObject Spectrograph (VIMOS) at the VLT, form the data base. We present isotropic and anisotropic Jeans-models and give a morphological description of the companion dwarf galaxy. The globular cluster system has about 2000 members, so it is not as rich as those of giant ellipticals in galaxy clusters with a comparable stellar mass, but richer than many cluster systems of other isolated ellipticals. The colour distribution of GCs is bimodal, which does not necessarily mean a metallicity bimodality. The kinematic literature data are somewhat inconclusive. The velocity dispersion in the inner parts can be reproduced without dark matter under isotropy. Radially anisotropic models need a low stellar mass-to-light ratio, which would contrast with the old age of the galaxy. A MONDian model is supported by X-ray analysis and previous dynamical modelling, but better data are necessary for a confirmation. The dwarf companion galaxy NGC 7796-1 exhibits tidal tails, multiple nuclei, and very boxy isophotes. NGC 7796 is an old, massive isolated elliptical galaxy with no indications of later major star formation events as seen frequently in other isolated ellipticals. Its relatively rich globular cluster system shows that isolation does not always mean a poor cluster system. The properties of the dwarf companion might indicate a dwarf-dwarf merger. (abridged)

We report on ongoing photometric and spectroscopic work on a sample of isolated elliptical galaxies. We investigate their globular cluster systems, and use the kinematics of globular clusters and the integrated galaxy light to constrain their dark halos, which are not found in the cases of NGC 5812 and NGC 7507

Context: There is strong evidence that globular cluster systems (GCSs) of massive galaxies are largely assembled by infall/accretion processes. Therefore, we expect the GCSs of isolated elliptical galaxies to be poor. Alhough not completely isolated, NGC 7507 is a massive field elliptical galaxy with an apparently very low dark matter content. Aims: We determine the richness, the colour distribution, and the structural properties of the GCS of NGC 7507. Methods: We perform wide-field Washington photometry with data obtained with the MOSAIC II camera at the 4m-Blanco telescope, CTIO. Results: The GCS is very poor with S_N ~ 0.6. We identify three subpopulations with peaks at (C-T1) colours of 1.21, 1.42, and 1.72. The bluest population may represent the old, metal-poor component. This interpretation is supported by its shallow density profile. The red population is more concentrated, resembling the galaxy light. The intermediate-colour population is strongly peaked in colour and we interpret this population as the signature of a starburst, whose age depends on the metallicity, but should be quite old, since no signatures of a merger are identifiable. In addition, we find a main sequence in the stellar foreground population, which we attribute to the Sagittarius dwarf tidal stream. Conclusions: The extraordinarily poor GCS of NGC 7507, a massive elliptical galaxy, is an illustration of how important the environmental conditions are for producing rich GCSs.

We present the kinematic results from our ARGOS spectroscopic survey of the Galactic bulge of the Milky Way. Our aim is to understand the formation of the Galactic bulge. We examine the kinematics of about 17,400 stars in the bulge located within 3.5 kpc of the Galactic centre, identified from the 28,000 star ARGOS survey. We aim to determine if the formation of the bulge has been internally driven from disk instabilities as suggested by its boxy shape, or if mergers have played a significant role as expected from Lambda CDM simulations. From our velocity measurements across latitudes b = -5 deg, -7.5 deg and -10 deg we find the bulge to be a cylindrically rotating system that transitions smoothly out into the disk. Within the bulge, we find a kinematically distinct metal-poor population ([Fe/H] < -1.0) that is not rotating cylindrically. The 5% of our stars with [Fe/H] < -1.0 are a slowly rotating spheroidal population, which we believe are stars of the metal weak thick disk and halo which presently lie in the inner Galaxy. The kinematics of the two bulge components that we identified in ARGOS paper III (mean [Fe/H] = -0.25 and [Fe/H] = +0.15, respectively) demonstrate that they are likely to share a common formation origin and are distinct from the more metal poor populations of the thick disk and halo which are colocated inside the bulge. We do not exclude an underlying merger generated bulge component but our results favour bulge formation from instabilities in the early thin disk.

The globular cluster (GC) systems of isolated elliptical galaxies have only recently begun to be studied in detail, and may exhibit morphological connections to the evolutionary histories of their hosts. Here we present the first in a series of wide-field analyses of the GC systems of the isolated ellipticals - Washington C and R photometry of NGC 3585 and NGC 5812 down to ~24 mag. The GC systems are characterised, with each system displaying both the "Universal" blue peak at (C-R)~1.3, and a red peak, but each with differing strengths. The total number of GCs in each system, and their specific frequencies, are estimated. The GC colours and specific frequencies are highly indicative that the host galaxy environment plays a role in shaping its GC system. We produce, and subtract, accurate models of each galaxy, revealing interesting underlying features, including the first definitive evidence that NGC 5812 is interacting with a dwarf companion galaxy. From the galaxy models we also determine surface brightness and colour profiles. Both colour profiles appear quite flat and with (C-R)~1.7 and we discuss the apparent youth of NGC 3585 in the context of this work.

We describe the motivation, field locations and stellar selection for the ARGOS spectroscopic survey of 28,000 stars in the bulge and inner disk of the Milky Way galaxy across latitudes of b = -5 deg to -10 deg. The primary goal of this survey is to constrain the formation processes of the bulge and establish whether it is predominantly a merger or instability remnant. From the spectra (R = 11,000), we have measured radial velocities and determined stellar parameters, including metallicities and [alpha/Fe] ratios. Distances were estimated from the derived stellar parameters and about 14,000 stars are red giants within 3.5 kpc of the Galactic centre. In this paper, we present the observations and analysis methods. Subsequent papers (III and IV) will discuss the stellar metallicity distribution and kinematics of the Galactic bulge and inner disk, and the implications for the formation of the bulge.

We present the metallicity results from the ARGOS spectroscopic survey of the Galactic bulge. Our aim is to understand the formation of the Galactic bulge: did it form via mergers, as expected from Lambda CDM theory, or from disk instabilities, as suggested by its boxy/peanut shape, or both? We have obtained spectra for 28,000 stars at a spectral resolution of R = 11,000. From these spectra, we have determined stellar parameters and distances to an accuracy of < 1.5 kpc. The stars in the inner Galaxy span a large range in [Fe/H], -2.8 < [Fe/H] < +0.6. From the spatial distribution of the red clump stars as a function of [Fe/H] (Ness et al. 2012a), we propose that the stars with [Fe/H] > -0.5 are part of the boxy/peanut bar/bulge. We associate the lower metallicity stars ([Fe/H] < -0.5) with the thick disk, which may be puffed up in the inner region, and with the inner regions of the metal-weak thick disk and inner halo. For the bulge stars with [Fe/H] > -0.5, we find two discrete populations; (i) stars with [Fe/H] ~ -0.25 which provide a roughly constant fraction of the stars in the latitude interval b = -5 deg to -10 deg, and (ii) a kinematically colder, more metal-rich population with mean [Fe/H] ~ +0.15 which is more prominent closer to the plane. The changing ratio of these components with latitude appears as a vertical abundance gradient of the bulge. We attribute both of these bulge components to instability-driven bar/bulge formation from the thin disk. We do not exclude a weak underlying classical merger-generated bulge component, but see no obvious kinematic association of any of our bulge stars with such a classical bulge component. [abridged]

Near the minor axis of the Galactic bulge, at latitudes b < -5 degrees, the red giant clump stars are split into two components along the line of sight. We investigate this split using the three fields from the ARGOS survey that lie on the minor axis at (l,b) = (0,-5), (0,-7.5), (0,-10) degrees. The separation is evident for stars with [Fe/H] > -0.5 in the two higher-latitude fields, but not in the field at b = -5 degrees. Stars with [Fe/H] < -0.5 do not show the split. We compare the spatial distribution and kinematics of the clump stars with predictions from an evolutionary N-body model of a bulge that grew from a disk via bar-related instabilities. The density distribution of the peanut-shaped model is depressed near its minor axis. This produces a bimodal distribution of stars along the line of sight through the bulge near its minor axis, very much as seen in our observations. The observed and modelled kinematics of the two groups of stars are also similar. We conclude that the split red clump of the bulge is probably a generic feature of boxy/peanut bulges that grew from disks, and that the disk from which the bulge grew had relatively few stars with [Fe/H] < -0.5

The Galactic globular cluster 47 Tucanae (47 Tuc) shows a rare increase in its velocity dispersion profile at large radii, indicative of energetic, yet bound, stars at large radii dominating the velocity dispersion and, potentially, of ongoing evaporation. Escaping stars will form tidal tails, as seen with several Galactic globular clusters, however, the tidal tails of 47 Tuc are yet to be uncovered. We model these tails of 47 Tuc using the most accurate input data available, with the specific aim of determining their locations, as well as the densities of the epicyclic overdensities within the tails. The overdensities from our models show an increase of 3-4% above the Galactic background and, therefore, should be easily detectable using matched filtering techniques. We find that the most influential parameter with regard to both the locations and densities of the epicyclic overdensities is the Heliocentric distance to the cluster. Hence, uncovering these tidal features observationally will contribute greatly to the ongoing problem of determining the distance to 47 Tuc, tightly constraining the distance of the cluster independent of other methods. Using our streakline method for determining the locations of the tidal tails and their overdensities, we show how, in principle, the shape and extent of the tidal tails of any Galactic globular cluster can be determined without resorting to computationally expensive N-body simulations.

We derive distance, density and metallicity distribution of the stellar Monoceros Overdensity (MO) in the outer Milky Way, based on deep imaging with the Subaru Telescope. We applied CMD fitting techniques in three stripes at galactic longitudes: l=130 deg, 150 deg, 170 deg; and galactic latitudes: +15 < b [deg] < +25 . The MO appears as a wall of stars at a heliocentric distance of ~ 10.1\pm0.5 kpc across the observed longitude range with no distance change. The MO stars are more metal rich ([Fe/H] ~ -1.0) than the nearby stars at the same latitude. These data are used to test three different models for the origin of the MO: a perturbed disc model, which predicts a significant drop in density adjacent to the MO that is not seen; a basic flared disc model, which can give a good match to the density profile but the MO metallicity implies the disc is too metal rich to source the MO stars; and a tidal stream model, which bracket the distances and densities we derive for the MO, suggesting that a model can be found that would fully fit the MO data. Further data and modeling will be required to confirm or rule out the MO feature as a stream or as a flaring of the disc.

We model the formation of Auriga's Wheel - a recently discovered collisional ring galaxy. Auriga's Wheel has a number of interesting features including a bridge of stars linking the neighbouring elliptical to the ring galaxy, and evidence for components of expansion and rotation within the ring. Using N-body/SPH modelling, we study collisions between an elliptical galaxy and a late-type disk galaxy. A near direct collision, with a mildy inclined disk, is found to reasonably reproduce the general system morphology ~50 Myr following the collision. The collision must have a relatively low velocity (initially ~150 km s^{-1}) in order to form the observed bridge, and simultaneously match the galaxies separation. Our best-match model suggests the total disk galaxy is ~5 times more massive than the elliptical. We find that the velocity of expansion of the ring is sensitive to the mass of the elliptical, while insensitive to the encounter velocity. We evolve our simulation beyond the current epoch to study the future destiny of the galaxy pair. In our model, the nucleus moves further away from the plane of the ring in the direction of the stellar bridge. The nucleus eventually merges with the elliptical galaxy ~100 Myr after the present time. The ring continues to expand for ~200 Myr before collapsing back. The low initial relative velocity of the two galaxies will eventually result in a complete merger.

We investigate the epicyclic motion of stars escaping from star clusters. Using streaklines, we visualise the path of escaping stars and show how epicyclic motion leads to over- and underdensities in tidal tails of star clusters moving on circular and eccentric orbits about a galaxy. Additionally, we investigate the effect of the cluster mass on the tidal tails, by showing that their structure is better matched when the perturbing effect of the cluster mass is included. By adjusting streaklines to results of N-body computations we can accurately and quickly reproduce all observed substructure, especially the streaky features often found in simulations which may be interpreted in observations as multiple tidal tails. Hence, we can rule out tidal shocks as the origin of such substructures. Finally, from the adjusted streakline parameters we can verify that for the star clusters we studied escape mainly happens from the tidal radius of the cluster, given by x_L = (GM/(\Omega^2-\partial^2\Phi/\partial R^2))^{1/3}. We find, however, that there is another limiting radius, the "edge" radius, which gives the smallest radius from which a star can escape during one cluster orbit about the galaxy. For eccentric cluster orbits the edge radius shrinks with increasing orbital eccentricity (for fixed apocentric distance) but is always significantly larger than the respective perigalactic tidal radius. In fact, the edge radii of the clusters we investigated, which are extended and tidally filling, agree well with their (fitted) King radii, which may indicate a fundamental connection between these two quantities.

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

We report the serendipitous discovery of a collision ring galaxy, identified as 2MASX J06470249+4554022, which we have dubbed 'Auriga's Wheel', found in a SUPRIME-CAM frame as part of a larger Milky Way survey. This peculiar class of galaxies is the result of a near head-on collision between typically, a late type and an early type galaxy. Subsequent GMOS-N long-slit spectroscopy has confirmed both the relative proximity of the components of this interacting pair and shown it to have a redshift of 0.111. Analysis of the spectroscopy reveals that the late type galaxy is a LINER class Active Galactic Nuclei while the early type galaxy is also potentially an AGN candidate, this is very uncommon amongst known collision ring galaxies. Preliminary modeling of the ring finds an expansion velocity of ~200 kms^-1 consistent with our observations, making the collision about 50 Myr old. The ring currently has a radius of about 10 kpc and a bridge of stars and gas is also visible connecting the two galaxies.

Galactic globular clusters have been pivotal in our understanding of many astrophysical phenomena. Here we publish the extracted stellar parameters from a recent large spectroscopic survey of ten globular clusters. A brief review of the project is also presented. Stellar parameters have been extracted from individual stellar spectra using both a modified version of the Radial Velocity Experiment (RAVE) pipeline and a pipeline based on the parameter estimation method of RAVE. We publish here all parameters extracted from both pipelines. We calibrate the metallicity and convert this to [Fe/H] for each star and, furthermore, we compare the velocities and velocity dispersions of the Galactic stars in each field to the Besan\c{c}on Galaxy model. We find that the model does not correspond well with the data, indicating that the model is probably of little use for comparisons with pencil beam survey data such as this.

We have conducted a spectroscopic survey of the inner regions of the Sagittarius (Sgr) dwarf galaxy using the AAOmega spectrograph on the Anglo-Australian Telescope. We determine radial velocities for over 1800 Sgr star members in 6 fields that cover an area 18.84 deg^2, with a typical accuracy of ~2 km/s. Motivated by recent numerical models of the Sgr tidal stream that predict a substantial amount of rotation in the dwarf remnant core, we compare the kinematic data against N-body models that simulate the stream progenitor as (i) a pressure-supported, mass-follows-light system, and (ii) a late-type, rotating disc galaxy embedded in an extended dark matter halo. We find that the models with little, or no intrinsic rotation clearly yield a better match to the mean line-of-sight velocity in all surveyed fields, but fail to reproduce the shape of the line-of-sight velocity distribution. This result rules out models wherein the prominent bifurcation observed in the leading tail of the Sgr stream was caused by a transfer from intrinsic angular momentum from the progenitor satellite into the tidal stream. It also implies that the trajectory of the young tidal tails has not been affected by internal rotation in the progenitor system. Our finding indicates that new, more elaborate dynamical models, in which the dark and luminous components are treated independently, are necessary for simultaneously reproducing both the internal kinematics of the Sgr dwarf and the available data for the associated tidal stream.

Globular clusters have proven to be essential to our understanding of many important astrophysical phenomena. Here we analyse spectroscopic observations of ten Halo globular clusters to determine their dark matter content, their tidal heating by the Galactic disc and halo, describe their metallicities and the likelihood that Newtonian dynamics explain their kinematics. We analyse a large number of members in all clusters, allowing us to address all these issues together, and we have included NGC 288 and M30 to overlap with previous studies. We find that any flattening of the velocity dispersion profiles in the outer regions of our clusters can be explained by tidal heating. We also find that all our GCs have M/L_V < 5, therefore, we infer the observed dynamics do not require dark matter, or a modification of gravity. We suggest that the lack of tidal heating signatures in distant clusters indicates the Halo is not triaxial. The isothermal rotations of each cluster are measured, with M4 and NGC 288 exhibiting rotation at a level of 0.9 +/- 0.1 km/s and 0.25 +/- 0.15 km/s, respectively. We also indirectly measure the tidal radius of NGC 6752, determining a more realistic figure for this cluster than current literature values. Lastly, an unresolved and intriguing puzzle is uncovered with regard to the cooling of the outer regions of all ten clusters.

The globular cluster 47 Tucanae is well studied but it has many characteristics that are unexplained, including a significant rise in the velocity dispersion profile at large radii, indicating the exciting possibility of two distinct kinematic populations. In this Letter we employ a Bayesian approach to the analysis of the largest available spectral dataset of 47 Tucanae to determine whether this apparently two-component population is real. Assuming the two models were equally likely before taking the data into account, we find that the evidence favours the two-component population model by a factor of ~3x10^7. Several possible explanations for this result are explored, namely the evaporation of low-mass stars, a hierarchical merger, extant remnants of two initially segregated populations, and multiple star formation epochs. We find the most compelling explanation for the two-component velocity distribution is that 47 Tuc formed as two separate populations arising from the same proto-cluster cloud which merged <7.3 +/- 1.5 Gyr ago. This may also explain the extreme rotation, low mass-to-light ratio and mixed stellar populations of this cluster.

Globular clusters are an important test bed for Newtonian gravity in the weak-acceleration regime, which is vital to our understanding of the nature of the gravitational interaction. Recent claims have been made that the velocity dispersion profiles of globular clusters flatten out at large radii, despite an apparent paucity of dark matter in such objects, indicating the need for a modification of gravitational theories. We continue our investigation of this claim, with the largest spectral samples ever obtained of 47 Tucanae and M55. Furthermore, this large sample allows for an accurate metallicity calibration based on the equivalent widths of the calcium triplet lines and K band magnitude of the Tip of the Red Giant Branch. Assuming an isothermal distribution, the rotations of each cluster are also measured with both clusters exhibiting clear rotation signatures. The global velocity dispersions of NGC 121 and Kron 3, two globular clusters in the Small Magellanic Cloud, are also calculated. By applying a simple dynamical model to the velocity dispersion profiles of 47 Tuc and M55, we calculate their mass-to-light profiles, total masses and central velocity dispersions. We find no statistically significant flattening of the velocity dispersion at large radii for M55, and a marked increase in the profile of 47 Tuc for radii greater than approximately half the tidal radius. We interpret this increase as an evaporation signature, indicating that 47 Tuc is undergoing, or has undergone, core-collapse, but find no requirement for dark matter or a modification of gravitational theories in either cluster.

Testing Newtonian gravity in the weak-acceleration regime is vital to our understanding of the nature of the gravitational interaction. It has recently been claimed that the velocity dispersion profiles of several globular clusters flatten out at large radii, reminiscent of galaxy rotation curves, even though globular clusters are thought to contain little or no dark matter. We investigate this claim, using AAOmega observations of four globular clusters, namely M22, M30, M53 and M68. M30, one such cluster that has had this claim made for its velocity dispersion, was included for comparison with previous studies. We find no statistically significant flattening of the velocity dispersion at large radii for any of our target clusters and therefore we infer the observed dynamics do not require that globular clusters are dark matter dominated, or a modification of gravity. Furthermore, by applying a simple dynamical model we determine the radial mass-to-light profiles for each cluster. The isothermal rotations of each cluster are also measured, with M22 exhibiting clear rotation, M68 possible rotation and M30 and M53 lacking any rotation, within the uncertainties.

This paper concludes a systematic search for evidence of the Monoceros Ring and Canis Major dwarf galaxy around the Galactic Plane. Presented here are the results for the Galactic longitude range of l = (280 - 025)deg. Testing the claim that the Monoceros Ring encircles the entire Galaxy, this survey attempts to document the position of the Monoceros Ring with increasing Galactic longitude. Additionally, with the discovery of the purported Canis Major dwarf galaxy, searching for more evidence of its interaction is imperative to tracing its path through the Galaxy and understanding its role in the evolution of the Milky Way. Two new detections of the Monoceros Ring have been found at (l, b) = (280,+15)deg and (300,+10)deg. Interestingly, in general there seem to be more detections above the Plane than below it; in this survey around 2/3 of the firm Monoceros Ring detections are in the North. This coincides with the Northern detections appearing to be qualitatively denser and broader than their Southern counterparts. The maximum of the Galactic Warp in the South is also probed in this survey. It is found that these fields do not resemble those in the Canis Major region suggesting that the Warp does not change the shape of the CMD as is witnessed around Canis Major. The origins and morphology of the Monoceros Ring is still elusive primarily due to its enormous extent on the sky. Continued probing of the Galactic Outer Disc is needed before a consensus can be reached on its nature.

We present the results of an AAT wide-field camera survey of the stars in the Monoceros Ring and purported Canis Major overdensity in the Galactic longitudes of {\it l} = (193 - 276)$^\circ$. Current numerical simulations suggest that both of these structures are the result of a single on-going accretion event, although an alternative solution is that the warped and flared disc of the Galaxy can explain the origin of both of these structures. Our results show that, with regards the Monoceros Ring, the warped and flared disc is unable to reproduce the locations and strengths of the detections observed around the Galaxy. This supports a non-Galactic origin for this structure. We report 8 new detections and 2 tentative detections of the Monoceros Ring in this survey. The exact nature of the Canis Major overdensity is still unresolved although this survey provides evidence that invoking the Galactic warp is not a sufficient solution when compared with observation. Several fields in this survey are highly inconsistent with the current Galactic disc models that include a warp and flare, to such an extent that explaining their origins with these structures is problematic. (Abridged)