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.

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

We describe the DR14 APOGEE-TGAS catalogue, a new SDSS value-added catalogue that provides precise astrophysical parameters, chemical abundances, astro-spectro-photometric distances and extinctions, as well as orbital parameters for $\sim 30,000$ APOGEE-TGAS stars, among them $\sim5,000$ high-quality giant stars within 1 kpc.

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

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

We study the properties of the central spheroids located within 10 kpc of the centre of mass of Milky Way mass-sized galaxies simulated in a cosmological context. The simulated central regions are dominated by stars older than 10 Gyr, mostly formed in situ, with a contribution of ~30 per cent from accreted stars. These stars formed in well-defined starbursts, although accreted stars exhibit sharper and earlier ones. The fraction of accreted stars increases with galactocentric distance, so that at a radius of ~8-10 kpc a fraction of ~40 per cent, on average, are detected. Accreted stars are slightly younger, lower metallicity, and more $\alpha$-enhanced than in situ stars. A significant fraction of old stars in the central regions come from a few ($2-3$) massive satellites ($\sim 10^{10}{\rm M}_\odot$). The bulge components receive larger contributions of accreted stars formed in dwarfs smaller than $\sim 10^{9.5}{\rm M}_\odot$. The difference between the distributions of ages and metallicities of old stars is thus linked to the accretion histories -- those central regions with a larger fraction of accreted stars are those with contributions from more massive satellites. The kinematical properties of in situ and accreted stars are consistent with the latter being supported by their velocity dispersions, while the former exhibit clear signatures of rotational support. Our simulations demonstrate a range of characteristics, with some systems exhibiting a co-existing bar and spheroid in their central regions, resembling in some respect the central region of the Milky Way.

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.

The galaxy formation process in the $\Lambda$-Cold Dark Matter scenario can be constrained from the analysis of stars in the Milky Way's halo system. We examine the variation of chemical abundances in distant halo stars observed by the Apache Point Galactic Evolution Experiment (APOGEE), as a function of distance from the Galactic center ($r$) and iron abundance ([M/H]), in the range 5 $\lesssim r \lesssim$ 30 kpc and $-2.5 <$ [M/H] $<$ 0.0. We perform a statistical analysis of the abundance ratios derived by the APOGEE pipeline (ASPCAP) and distances calculated by several approaches. Our analysis reveals signatures of a different chemical enrichment between the inner and outer regions of the halo, with a transition at about 15 kpc. The derived metallicity distribution function exhibits two peaks, at [M/H] $\sim -1.5$ and $\sim -2.1$, consistent with previously reported halo metallicity distributions. We obtain a difference of $\sim 0.1$ dex for $\alpha$-element-to-iron ratios for stars at $r > 15$ kpc and [M/H] $> -1.1$ (larger in the case of O, Mg and S) with respect to the nearest halo stars. This result confirms previous claims for low-$\alpha$ stars found at larger distances. Chemical differences in elements with other nucleosynthetic origins (Ni, K, Na, and Al) are also detected. C and N do not provide reliable information about the interstellar medium from which stars formed because our sample comprises RGB and AGB stars and can experience mixing of material to their surfaces.

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.

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

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

The Kapteyn moving group has been postulated as tidal debris from $\omega$ Centauri. If true, members of the group should show some of the chemical abundance patterns known for stars in the cluster. We present an optical and near-infrared high-resolution, high-S/N spectroscopic study of 14 stars of the Kapteyn group, plus 10 additional stars (the $\omega$ Cen-group) that, while not listed as members of the Kapteyn group as originally defined, have been nevertheless associated dynamically with $\omega$ Centauri. Abundances for Na, O, Mg, Al, Ca and Ba were derived from the optical spectra, while the strength of the chromospheric He I 10830 {\AA} line is studied as a possible helium abundance indicator. The resulting Na-O and Mg-Al patterns for stars of the combined Kapteyn and $\omega$ Cen-group samples do not resemble those of $\omega$ Centauri, and are not different from those of field stars of the Galactic halo. The distribution of equivalent widths of the He I 10830 {\AA} line is consistent with that found among non-active field stars. Therefore, no evidence is found for second-generation stars within our samples, which most likely rules out a globular-cluster origin. Moreover, no hint of the unique Ba-overabundance at the metal-rich end, well-established for $\omega$ Centauri stars, is seen among stars of the combined samples. Because this specific Ba pattern is present in $\omega$ Centauri irrespective of stellar generation, this would rule out the possibility that our entire sample might be composed of only first generation stars from the cluster. Finally, for the stars of the Kapteyn group, the possibility of an origin in the hypothetical $\omega$ Centauri's parent galaxy is disfavored by the different run of $\alpha$-elements with metallicity between our targets and stars from present-day dwarf galaxies.

The age and metallicity gradients for a sample of group and cluster galaxies from N-body+hydrodynamical simulation are analyzed in terms of galaxy stellar mass. Dwarf galaxies show null age gradient with a tail of high and positive values for systems in groups and cluster outskirts. Massive systems have generally zero age gradients which turn to positive for the most massive ones. Metallicity gradients are distributed around zero in dwarf galaxies and become more negative with mass; massive galaxies have steeper negative metallicity gradients, but the trend flatten with mass. In particular, fossil groups are characterized by a tighter distribution of both age and metallicity gradients. We find a good agreement with both local observations and independent simulations. The results are also discussed in terms of the central age and metallicity, as well as the total colour, specific star formation and velocity dispersion.

The evolution of marginally bound supercluster-like objects in an accelerating LambdaCDM Universe is followed, by means of cosmological simulations, from the present time to an expansion factor a = 100. The objects are identified on the basis of the binding density criterion introduced by Dunner et al. (2006). superclusters are identified with the ones whose mass M > 10^15 M_sun/h, the most massive one with M ~ 8x10^15 M_sun/h, comparable to the Shapley supercluster. The spatial distribution of the superclusters remains essentially the same after the present epoch, reflecting the halting growth of the Cosmic Web as Lambda gets to dominate the expansion of the Universe. The same trend can be seen in the stagnation of the development of the mass function of virialized haloes and bound objects. The situation is considerably different when looking at the internal evolution, quantified in terms of their shape, compactness and density profile, and substructure in terms of their multiplicity function. We find a continuing evolution from a wide range of triaxial shapes at a = 1 to almost perfect spherical shapes at a = 100. We also find a systematic trend towards a higher concentration. Meanwhile, we see their substructure gradually disappearing, as the surrounding subclumps fall in and merge to form one coherent, virialized system.

An exponentially expanding Universe, possibly governed by a cosmological constant, forces gravitationally bound structures to become more and more isolated, eventually becoming causally disconnected from each other and forming so-called "island universes". This new scenario reformulates the question about which will be the largest structures that will remain gravitationally bound, together with requiring a systematic tool that can be used to recognize the limits and mass of these structures from observational data, namely redshift surveys of galaxies. Here we present a method, based on the spherical collapse model and N-body simulations, by which we can estimate the limits of bound structures as observed in redshift space. The method is based on a theoretical criterion presented in a previous paper that determines the mean density contrast that a spherical shell must have in order to be marginally bound to the massive structure within it. Understanding the kinematics of the system, we translated the real-space limiting conditions of this "critical" shell to redshift space, producing a projected velocity envelope that only depends on the density profile of the structure. From it we created a redshift-space version of the density contrast that we called "density estimator", which can be calibrated from N-body simulations for a reasonable projected velocity envelope template, and used to estimate the limits and mass of a structure only from its redshift-space coordinates.

According to the latest evidence, the Universe is entering an era of exponential expansion, where gravitationally bound structures will get disconnected from each other, forming isolated `island universes'. In this scenario, we present a theoretical criterion to determine the boundaries of gravitationally bound structures and a physically motivated definition of superclusters as the largest bound structures in the Universe. We use the spherical collapse model self-consistently to obtain an analytical condition for the mean density enclosed by the last bound shell of the structure (2.36 times the critical density in the present Universe, assumed to be flat, with 30 per cent matter and 70 per cent cosmological constant, in agreement with the previous, numerical result of Chiueh and He). $N$-body simulations extended to the future show that this criterion, applied at the present cosmological epoch, defines a sphere that encloses $\approx 99.7$ per cent of the particles that will remain bound to the structure at least until the scale parameter of the Universe is 100 times its present value. On the other hand, $(28\pm 13)$ per cent of the enclosed particles are in fact not bound, so the enclosed mass overestimates the bound mass, in contrast with the previous, less rigorous criterion of, e.~g., Busha and collaborators, which gave a more precise mass estimate. We also verify that the spherical collapse model estimate for the radial infall velocity of a shell enclosing a given mean density gives an accurate prediction for the velocity profile of infalling particles, down to very near the centre of the virialized core.

We use numerical simulations to investigate the orbital characteristics of tidal debris from satellites whose orbits are dragged into the plane of galactic disks by dynamical friction before disruption. We find that these satellites may deposit a significant fraction of their stars into the disk components of a galaxy, and use our results to motivate the search for accretion relicts in samples of metal-poor disk stars in the vicinity of the Sun. Satellites disrupted on very eccentric orbits coplanar with the disk are expected to shed stars in "trails" of distinct orbital energy and angular momentum during each pericentric passage. To an observer located between the pericenter and apocenter of such orbits, these trails would show as distinct groupings of stars with low vertical velocity and a broad, symmetric, often double-peaked distribution of Galactocentric radial velocities. One group of stars with these characteristics stands out in available compilations of nearby metal-poor stars. These stars have specific angular momenta similar to that of the globular cluster omegaCen, long hypothesized to be the nucleus of a dwarf galaxy disrupted by the Milky Way tidal field. In addition to their kindred kinematics, stars in the omegaCen group share distinct chemical abundance characteristics, and trace a well-defined track in the [alpha/Fe] versus [Fe/H] plane, consistent with simple closed-box enrichment models and a protracted star formation history. The dynamical and chemical coherence of this group suggests that it consists of stars that once belonged to the dwarf that brought omegaCen into the Galaxy. The presence of this and other "tidal relicts" in the solar neighbourhood suggest an extra-Galactic origin for the presence of nearby stars with odd kinematics and chemistry.

We examine in detail the dynamical structure of an elliptical galaxy simulated in the Lambda CDM scenario. The morphology of the galaxy evolves dramatically over time in response to the mode and timing of mass accretion; smooth deposition of cooled gas leads to the formation of centrifugally supported disks, whilst major mergers disperse stellar disks into spheroids. These two modes of accretion alternate successively until z~0.6, when the galaxy undergoes one last major (1:2) merger that consumes much of the remaining gas into stars. Little gas cools and accretes subsequently and, as a result, most stars at z=0 are in a spheroidal component that resembles present-day elliptical galaxies. Dynamically, the galaxy is well approximated by an E4 oblate rotator. Boxy isophotes are obtained when the galaxy is seen face-on and Vrot/sigma<<1. On the other hand, disky isophotes are found for inclinations which maximize Vrot/sigma. The line-of-sight velocity distribution is nearly Gaussian at all radii. The sign of the Gauss-Hermite skewness parameter h3 is anti-correlated with the apparent sense of rotation, in agreement with observed trends. The simulated galaxy has much higher effective surface brightness than normal ellipticals of similar luminosity, in a way reminiscent of the less common type of compact ellipticals. Our simulation shows that repeated episodes of dissipational collapse, followed by merger events, lead to stellar spheroids that are only mildly triaxial and of relatively simple kinematic structure. This is in better agreement with observation than earlier models based on dissipationless mergers of stellar disks, and a positive step towards reconciling the observed structure of ellipticals with the hierarchical merging scenario.

We examine the tidal disruption of satellite galaxies in a cosmological simulation of the formation of a disk galaxy in the $\Lambda$CDM scenario. We find that the disruption of satellite galaxies in orbits roughly coplanar with the disk leads naturally to the formation of ring-like stellar structures similar to that recently discovered in the outer disk of the Milky Way. Two interpretations appear plausible in this context. One is that the ring is a transitory, localized radial density enhancement reflecting the apocenter of particles stripped from a satellite during a recent pericentric passage (a ``tidal arc'' reminiscent of the tidal arms seen in disk galaxy mergers). In the second scenario, the ring is analogous to the ``shells'' found around some elliptical galaxies, and would result from a minor merger that took place several Gyr ago. The two interpretations differ in several ways. Tidal arcs are expected to span a limited longitude range; may carry a substantial fraction of the original mass of the satellite; should exhibit a significant velocity gradient with Galactic longitude; and are in all likelihood asymmetric above and below the plane of the disk. Shells, on the other hand, may comprise at most a small fraction of the original mass of the satellite and, due to their more relaxed state, ought to be symmetric above and below the plane, with no discernible velocity gradients across the structure. If confirmed as a tidal feature, the ring discovered by SDSS in the outer Galactic disk would strengthen the view--supported by numerical simulations--that minor mergers have played a significant role building up not only the stellar halo, but also the disk components of the Galaxy.

The stability of rotating isotropic spherical stellar systems is investigated by using N-body simulations. Four spherical models with realistic density profiles are studied: one of them fits the luminosity profile of globular clusters, while the remaining three models provide good approximations to the surface brightness of elliptical galaxies. The phase-space distribution function f(E) of each one of these non-rotating models satisfies the sufficient condition for stability df/dE < 0. Different amounts of rotation are introduced in these models by changing the sign of the z-component of the angular momentum for a given fraction of the particles. Numerical simulations show that all these rotating models are stable to both radial and non-radial perturbations, irrespective of their degree of rotation. These results suggest that rotating isotropic spherical models with realistic density profiles might generally be stable. Furthermore, they show that spherical stellar systems can rotate very rapidly without becoming oblate.

We have investigated the stability of a set of non-rotating anisotropic spherical models with a phase-space distribution function of the Osipkov-Merritt type. The velocity distribution in these models is isotropic near the center and becomes radially anisotropic at large radii. They are special members of the family studied by Dehnen and Tremaine et al. where the mass density has a power-law cusp $\rho\propto r^{-\gamma}$ at small radii and decays as $\rho\propto r^{-4}$ at large radii. The radial-orbit instability of models with $\gamma$ = 0, 1/2, 1, 3/2, and 2, was studied using an N-body code written by one of us and based on the `self-consistent field' method developed by Hernquist and Ostriker. These simulations have allowed us to delineate a boundary in the $(\gamma,r_{a})$-plane that separates the stable from the unstable models. This boundary is given by $2T_{r}/T_{t} = 2.31 \pm 0.27$, for the ratio of the total radial to tangential kinetic energy. We also found that the stability criterion $df/dQ\le 0$, recently raised by Hjorth, gives lower values compared with our numerical results.