A simple model is presented of the formation of the spiral the (z,v_z) phase plane of solar-neighbourhood stars that was recently discovered in Gaia data. The key is that the frequency Omega_z at which stars oscillate vertically depends on angular momentum about the z axis in addition to the amplitude of the star's vertical oscillations. Spirals should form in both <v_phi> and <v_R> whenever a massive substructure, such as the Sgr dwarf galaxy, passes through the Galactic plane. The model yields similar spirals to those observed in both <v_phi> and <v_R>. The primary driver is the component of the tidal force that lies in the plane. We investigate the longevity of the spirals and the mass of the substructure, but the approximations inherent in the model make quantitative results unreliable. The work relies heavily on a self-consistent, multi-component model of our Galaxy produced by the AGAMA package for f(J) modelling.

Our Galaxy is thought to have undergone an active evolutionary history dominated by star formation, the accretion of cold gas, and, in particular, mergers up to 10 gigayear ago. The stellar halo reveals rich fossil evidence of these interactions in the form of stellar streams, substructures, and chemically distinct stellar components. The impact of dwarf galaxy mergers on the content and morphology of the Galactic disk is still being explored. Recent studies have identified kinematically distinct stellar substructures and moving groups, which may have extragalactic origin. However, there is mounting evidence that stellar overdensities at the outer disk/halo interface could have been caused by the interaction of a dwarf galaxy with the disk. Here we report detailed spectroscopic analysis of 14 stars drawn from two stellar overdensities, each lying about 5 kiloparsecs above and below the Galactic plane - locations suggestive of association with the stellar halo. However, we find that the chemical compositions of these stars are almost identical, both within and between these groups, and closely match the abundance patterns of the Milky Way disk stars. This study hence provides compelling evidence that these stars originate from the disk and the overdensities they are part of were created by tidal interactions of the disk with passing or merging dwarf galaxies.

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.

From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (LTE, 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature, and, in particular, for LTE vs. non-LTE (NLTE) on metallicity of the stars. Here we analyse the [Mg/Fe] and [Fe/H] plane of a sample of 326 stars, comparing LTE and NLTE results obtained using 1D hydrostatic models and averaged <3D> models. We show that compared to the <3D>NLTE benchmark, all other three methods display increasing biases towards lower metallicities, resulting in false trends of [Mg/Fe] against [Fe/H], which have profound implications for interpretations by chemical evolution models. In our best <3D> NLTE model, the halo and disc stars show a clearer behaviour in the [Mg/Fe] - [Fe/H] plane, from the knee in abundance space down to the lowest metallicities. Our sample has a large fraction of thick disc stars and this population extends down to at least [Fe/H] ~ -1.6 dex, further than previously proven. The thick disc stars display a constant [Mg/Fe] ~ 0.3 dex, with a small intrinsic dispersion in [Mg/Fe] that suggests that a fast SN Ia channel is not relevant for the disc formation. The halo stars reach higher [Mg/Fe] ratios and display a net trend of [Mg/Fe] at low metallicities, paired with a large dispersion in [Mg/Fe]. These indicate the diverse origin of halo stars from accreted low-mass systems to stochastic/inhomogeneous chemical evolution in the Galactic halo.

Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known. Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.

Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.

The existence of a vertical age gradient in the Milky Way disc has been indirectly known for long. Here, we measure it directly for the first time with seismic ages, using red giants observed by Kepler. We use Stroemgren photometry to gauge the selection function of asteroseismic targets, and derive colour and magnitude limits where giants with measured oscillations are representative of the underlying population in the field. Limits in the 2MASS system are also derived. We lay out a method to assess and correct for target selection effects independent of Galaxy models. We find that low mass, i.e. old red giants dominate at increasing Galactic heights, whereas closer to the Galactic plane they exhibit a wide range of ages and metallicities. Parametrizing this as a vertical gradient returns approximately 4 Gyr/kpc for the disc we probe, although with a large dispersion of ages at all heights. The ages of stars show a smooth distribution over the last 10 Gyr, consistent with a mostly quiescent evolution for the Milky Way disc since a redshift of about 2. We also find a flat age-metallicity relation for disc stars. Finally, we show how to use secondary clump stars to estimate the present-day intrinsic metallicity spread, and suggest using their number count as a new proxy for tracing the ageing of the disc. This work highlights the power of asteroseismology for Galactic studies; however, we also emphasize the need for better constraints on stellar mass-loss, which is a major source of systematic age uncertainties in red giant stars.

Using G dwarfs from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey, we have determined a vertical metallicity gradient over a large volume of the Milky Way's disk, and examined how this gradient varies for different [a/Fe] subsamples. This sample contains over 40,000 stars with low-resolution spectroscopy over 144 lines of sight. We employ the SEGUE Stellar Parameter Pipeline (SSPP) to obtain estimates of effective temperature, surface gravity, [Fe/H], and [a/Fe] for each star and extract multiple volume-complete subsamples of approximately 1000 stars each. Based on the survey's consistent target-selection algorithm, we adjust each subsample to determine an unbiased picture of the disk in [Fe/H] and [a/Fe]; consequently, each individual star represents the properties of many. The SEGUE sample allows us to constrain the vertical metallicity gradient for a large number of stars over a significant volume of the disk, between ~0.3 and 1.6 kpc from the Galactic plane, and examine the in situ structure, in contrast to previous analyses which are more limited in scope. This work does not pre-suppose a disk structure, whether composed of a single complex population or a distinct thin and thick disk component. The metallicity gradient is -0.243 +0.039 -0.053 dex/kpc for the sample as a whole, which we compare to various literature results. Each [a/Fe] subsample dominates at a different range of heights above the plane of the Galaxy, which is exhibited in the gradient found in the sample as a whole. Stars over a limited range in [a/Fe] show little change in median [Fe/H] with height. If we associate [a/Fe] with age, our consistent vertical metallicity gradients with [a/Fe] suggest that stars formed in different epochs exhibit comparable vertical structure, implying similar star-formation processes and evolution.

Asteroseismology has the capability of precisely determining stellar properties which would otherwise be inaccessible, such as radii, masses and thus ages of stars. When coupling this information with classical determinations of stellar parameters, such as metallicities, effective temperatures and angular diameters, powerful new diagnostics for Galactic studies can be obtained. The ongoing Stroemgren survey for Asteroseismology and Galactic Archaeology (SAGA) has the goal of transforming the Kepler field into a new benchmark for Galactic studies, similarly to the solar neighborhood. Here we present first results from a stripe centred at Galactic longitude 74deg and covering latitude from about 8 to 20deg, which includes almost 1000 K-giants with seismic information and the benchmark open cluster NGC 6819. We describe the coupling of classical and seismic parameters, the accuracy as well as the caveats of the derived effective temperatures, metallicities, distances, surface gravities, masses, and radii. Confidence in the achieved precision is corroborated by the detection of the first and secondary clump in a population of field stars with a ratio of 2 to 1, and by the negligible scatter in the seismic distances among NGC 6819 member stars. An assessment of the reliability of stellar parameters in the Kepler Input Catalogue is also performed, and the impact of our results for population studies in the Milky Way is discussed, along with the importance of an all-sky Stroemgren survey.

We re-examine recent claims of observational evidence for a dual Galactic halo in SEGUE/SDSS data, and trace them back to improper error treatment and neglect of selection effects. In particular, the detection of a vertical abundance gradient in the halo can be explained as a metallicity bias in distance. A similar bias, and the impact of disk contamination, affect the sample of blue horizontal branch stars. These examples highlight why non-volume complete samples require forward-modelling from theoretical models or extensive bias-corrections. We also show how observational uncertainties produce the specific non-Gaussianity in the observed azimuthal velocity distribution of halo stars, which can be erroneously identified as two Gaussian components. A single kinematic component yields an excellent fit to the observed data, when we model the measurement process including distance uncertainties. Furthermore, we show that sample differences in proper motion space are the direct consequence of kinematic cuts, and are enhanced when distance estimates are less accurate. Thus, their presence is neither a proof for a separate population, nor a measure of reliability for the applied distances. We conclude that currently there is no evidence from SEGUE/SDSS that would favour a dual Galactic halo over a single halo full of substructure.

We present the metallicity distribution function (MDF) for 24,270 G and 16,847 K dwarfs at distances from 0.2 to 2.3 kpc from the Galactic plane, based on spectroscopy from the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. This stellar sample is significantly larger in both number and volume than previous spectroscopic analyses, which were limited to the solar vicinity, making it ideal for comparison with local volume-limited samples and Galactic models. For the first time, we have corrected the MDF for the various observational biases introduced by the SEGUE target selection strategy. The SEGUE sample is particularly notable for K dwarfs, which are too faint to examine spectroscopically far from the solar neighborhood. The MDF of both spectral types becomes more metal-poor with increasing |Z|, which reflects the transition from a sample with small [alpha/Fe] values at small heights to one with enhanced [alpha/Fe] above 1 kpc. Comparison of our SEGUE distributions to those of two different Milky Way models reveals that both are more metal-rich than our observed distributions at all heights above the plane. Our unbiased observations of G and K dwarfs provide valuable constraints over the |Z|-height range of the Milky Way disk for chemical and dynamical Galaxy evolution models, previously only calibrated to the solar neighborhood, with particular utility for thin- and thick-disk formation models.

Commissioning observations with the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III, have produced radial velocities (RVs) for ~4700 K/M-giant stars in the Milky Way bulge. These high-resolution (R \sim 22,500), high-S/N (>100 per resolution element), near-infrared (1.51-1.70 um; NIR) spectra provide accurate RVs (epsilon_v~0.2 km/s) for the sample of stars in 18 Galactic bulge fields spanning -1<l<20 deg, |b|<20 deg, and dec>-32 deg. This represents the largest NIR high-resolution spectroscopic sample of giant stars ever assembled in this region of the Galaxy. A cold (sigma_v~30 km/s), high-velocity peak (V_GSR \sim +200 km/s) is found to comprise a significant fraction (~10%) of stars in many of these fields. These high RVs have not been detected in previous MW surveys and are not expected for a simple, circularly rotating disk. Preliminary distance estimates rule out an origin from the background Sagittarius tidal stream or a new stream in the MW disk. Comparison to various Galactic models suggests that these high RVs are best explained by stars in orbits of the Galactic bar potential, although some observational features remain unexplained.

We present a study of the extent to which the Sellwood & Binney radial migration of stars is affected by their vertical motion about the midplane. We use both controlled simulations in which only a single spiral mode is excited, as well as slightly more realistic cases with multiple spiral patterns and a bar. We find that rms angular momentum changes are reduced by vertical motion, but rather gradually, and the maximum changes are almost as large for thick disc stars as for those in a thin disc. We find that particles in simulations in which a bar forms suffer slightly larger angular momentum changes than in comparable cases with no bar, but the cumulative effect of multiple spiral events still dominates. We have determined that vertical action, and not vertical energy, is conserved on average during radial migration.

We present a new method for detecting and correcting systematic errors in the distances to stars when both proper motions and line-of-sight velocities are available. The method, which is applicable for samples of 200 or more stars that have a significant extension on the sky, exploits correlations between the measured U, V and W velocity components that are introduced by distance errors. We deliver a formalism to describe and interpret the specific imprints of distance errors including spurious velocity correlations and shifts of mean motion in a sample. We take into account correlations introduced by measurement errors, Galactic rotation and changes in the orientation of the velocity ellipsoid with position in the Galaxy. Tests on pseudodata show that the method is more robust and sensitive than traditional approaches to this problem. We investigate approaches to characterising the probability distribution of distance errors, in addition to the mean distance error, which is the main theme of the paper. Stars with the most overestimated distances bias our estimate of the overall distance scale, leading to the corrected distances being slightly too small. We give a formula that can be used to correct for this effect. We apply the method to samples of stars from the SEGUE survey, exploring optimal gravity cuts, sample contamination, and correcting the used distance relations.

This note attempts to interpret some of the recent findings about a downtrend in the mean azimuthal velocity of low [alpha/Fe] thin disc stars with increasing metallicity. The presence of such a trend was predicted in the model of Schoenrich & Binney (2009), albeit with a slightly steeper slope. We show that in a simple picture a Galactic disc without mixing in angular momenta would display an exceedingly steep trend, while in the case of complete mixing of all stars the trend has to vanish. The difference between model and observational data can hence be interpreted as the consequence of the radial abundance gradient in the model being too high resulting in an underestimate of the migration strength. We shortly discuss the value of eccentricity distributions in constraining structure and history of the Galactic disc.

In a disc galaxy the distribution of azimuthal components of velocity is very skew. In the past this skewness has been modelled by superposed Gaussians. We use dynamical arguments to derive an analytic formula that can be fitted to observed velocity distributions, and validate it by fits to the velocities derived from a dynamically rigorous model, and to a sample of local stars with accurate space velocities. Our formula is much easier to use than a full distribution function. It has fewer parameters than a multi-Gaussian fit, and the best-fitting model parameters give insight into the underlying disc dynamics. In particular, once the azimuthal velocities of a sample have been successfully fitted, the apparatus provides a prediction for the corresponding distribution of radial velocities vR . An effective formula like ours is invaluable when fitting to data for stars at some distance from the Sun because it enables one to make proper allowance for the errors in distance and proper motion when determining the underlying disc kinematics. The derivation of our formula elucidates the way the horizontal and vertical motions are closely intertwined, and makes it evident that no stellar population can have a scale height and vertical velocity dispersions that are simultaneously independent of radius. We show that the oscillation of a star perpendicular to the Galactic plane modifies the effective potential in which the star moves radially in such a way that the more vertical energy a star has, the larger is the mean radius of its orbit.

We present a re-analysis of the Geneva-Copenhagen survey, which benefits from the infrared flux method to improve the accuracy of the derived stellar effective temperatures and uses the latter to build a consistent and improved metallicity scale. Metallicities are calibrated on high-resolution spectroscopy and checked against four open clusters and a moving group, showing excellent consistency. The new temperature and metallicity scales provide a better match to theoretical isochrones, which are used for a Bayesian analysis of stellar ages. With respect to previous analyses, our stars are on average 100 K hotter and 0.1 dex more metal rich, which shift the peak of the metallicity distribution function around the solar value. From Stromgren photometry we are able to derive for the first time a proxy for alpha elements, which enables us to perform a tentative dissection of the chemical thin and thick disc. We find evidence for the latter being composed of an old, mildly but systematically alpha-enhanced population that extends to super solar metallicities, in agreement with spectroscopic studies. Our revision offers the largest existing kinematically unbiased sample of the solar neighbourhood that contains full information on kinematics, metallicities, and ages and thus provides better constraints on the physical processes relevant in the build-up of the Milky Way disc, enabling a better understanding of the Sun in a Galactic context.

We examine the kinematics of the Galactic halo based on SDSS/SEGUE data by Carollo et al. (2007, 2010). We find that their claims of a counter-rotating halo are the result of substantial biases in distance estimates (of order 50%): the claimed retrograde component, which makes up only a tiny fraction of the entire sample, prone to contaminations, is identified as the tail of distance overestimates. The strong overestimates also result in a lift in the vertical velocity component, which explains the large altitudes those objects were claimed to reach. Errors are worst for the lowest metallicity stars, which explains the metal-poor nature of the artificial component. We also argue that measurement errors were not properly accounted for and that the use of Gaussian fitting on intrinsically non-Gaussian Galactic components invokes the identification of components that are distorted or even artificial. Our evaluation of the data leads to a revision of the estimated velocity ellipsoids and does not yield any reliable evidence for a counterrotating halo component. If a distinct counterrotating halo component exists it must be far weaker than claimed by Carollo et al. Finally we note that their revised analysis presented in Beers et al. (2011) does not alleviate our main concerns.

We re-examine the stellar kinematics of the Solar neighbourhood in terms of the velocity of the Sun with respect to the local standard of rest. We show that the classical determination of its component V_sun in the direction of Galactic rotation via Stroemberg's relation is undermined by the metallicity gradient in the disc, which introduces a correlation between the colour of a group of stars and the radial gradients of its properties. Comparing the local stellar kinematics to a chemodynamical model which accounts for these effects, we obtain (U,V,W)_sun = (11.1 +/- 0.74, 12.24 +/- 0.47, 7.25 +/-0.37) km/s, with additional systematic uncertainties of ~ (1,2,0.5) km/s. In particular, V_sun is 7 km/s larger than previously estimated. The new values of solar motion are extremely insensitive to the metallicity gradient within the disc.

We examine the chemical and dynamical structure in the solar neighbourhood of a model Galaxy that is the endpoint of a simulation of the chemical evolution of the Milky Way in the presence of radial mixing of stars and gas. Although the simulation's star-formation rate declines monotonically from its unique peak and no merger or tidal event ever takes place, the model replicates all known properties of a thick disc, as well as matching special features of the local stellar population such as a metal-poor extension of the thin disc that has high rotational velocity. We divide the disc by chemistry and relate this dissection to observationally more convenient kinematic selection criteria. We conclude that the observed chemistry of the Galactic disc does not provide convincing evidence for a violent origin of the thick disc, as has been widely claimed.

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and alpha elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva-Copenhagen survey (GCS) of the solar neighbourhood, and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age-metallicity plane although this plane was not used in the fitting process. Although this model's star-formation rate is monotonic declining, its disc naturally splits into an alpha-enhanced thick disc and a normal thin disc. In particular the model's distribution of stars in the ([O/Fe],[Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high sigma_z from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of sigma_z with age than predicted by theory. The metallicity gradient in the ISM is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cutoff in the star-formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.

The phase of Super-Soft-Source (SSS) emission of the sixth recorded outburst of the recurrent nova RS Oph was observed twice with Chandra and once with XMM-Newton. The observations were taken on days 39.7, 54.0, and 66.9 after outburst. We confirm a 35-sec period on day 54.0 and found that it originates from the SSS emission and not from the shock. We discus the bound-free absorption by neutral elements in the line of sight, resonance absorption lines plus self-absorbed emission line components, collisionally excited emission lines from the shock, He-like intersystem lines, and spectral changes during an episode of high-amplitude variability. We find a decrease of the oxygen K-shell absorption edge that can be explained by photoionization of oxygen. The absorption component has average velocities of -1286+-267 km/s on day 39.7 and of -771+-65 km/s on day 66.9. The wavelengths of the emission line components are consistent with their rest wavelengths as confirmed by measurements of non-self absorbed He-like intersystem lines. We have evidence that these lines originate from the shock rather than the outer layers of the outflow and may be photoexcited in addition to collisional excitations. We found collisionally excited emission lines that are fading at wavelengths shorter than 15A that originate from the radiatively cooling shock. On day 39.5 we find a systematic blue shift of -526+-114 km/s from these lines. We found anomalous He-like f/i ratios which indicates either high densities or significant UV radiation near the plasma where the emission lines are formed. During the phase of strong variability the spectral hardness light curve overlies the total light curve when shifted by 1000sec. This can be explained by photoionization of neutral oxygen in the line of sight if the densities of order 10^{10}-10^{11} cm^{-3}.