The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are present. Here we correct for this effect and derive the mass formed in stars by $z=2$ for a sample of 16 Local Group dwarf galaxies. We explore early differences between ${\it fast}$ and ${\it slow}$ dwarfs, and evaluate the impact of internal feedback by supernovae (SN) on the baryonic and dark matter (DM) component of the dwarfs. ${\it Fast}$ dwarfs assembled more stellar mass at early times and have larger amounts of DM within the half-light radius than ${\it slow}$ dwarfs. By imposing that ${\it slow}$ dwarfs cannot have lost their gas by $z=2$, we constrain the maximum coupling efficiency of SN feedback to the gas and to the DM to be $\sim$10%. We find that internal feedback alone appears insufficient to quench the SFH of ${\it fast}$ dwarfs by gas deprivation, in particular for the fainter systems. Nonetheless, SN feedback can core the DM halo density profiles relatively easily, producing cores of the sizes of the half-light radius in ${\it fast}$ dwarfs by $z=2$ with very low efficiencies. Amongst the "classical" Milky Way satellites, we predict that the smallest cores should be found in Draco and Ursa Minor, while Sculptor and Fornax should host the largest ones.

We confirm the reality of the recently discovered Milky Way stellar cluster $\textit{Gaia}$ 1 using spectra acquired with the HERMES and AAOmega spectrographs of the Anglo-Australian Telescope. This cluster had been previously undiscovered due to its close angular proximity to Sirius, the brightest star in the sky at visual wavelengths. Our observations identified 41 cluster members, and yielded an overall metallicity of [Fe/H]$=-0.13\pm0.13$ and barycentric radial velocity of $v_r=58.30\pm0.22$ km/s. These kinematics provide a dynamical mass estimate of $12.9^{+4.6}_{-3.9}\times10^3$ M$_{\odot}$. Isochrone fits to $\textit{Gaia}$, 2MASS, and Pan-STARRS1 photometry indicate that $\textit{Gaia}$ 1 is an intermediate age ($\sim3$ Gyr) stellar cluster. Combining the spatial and kinematic data we calculate $\textit{Gaia}$ 1 has a circular orbit with a radius of about 12~kpc, but with a large out of plane motion: $z_\textrm{max}=1.1^{+0.4}_{-0.3}$ kpc. Clusters with such orbits are unlikely to survive long due to the number of plane passages they would experience.

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.

A central tenet of the current cosmological paradigm is that galaxies grow over time through the accretion of smaller systems. Here, we present new kinematic measurements near the centre of one of the densest pronounced substructures, the South-West Cloud, in the outer halo of our nearest giant neighbour, the Andromeda galaxy. These observations reveal that the kinematic properties of this region of the South-West Cloud are consistent with those of PA-8, a globular cluster previously shown to be co-spatial with the stellar substructure. In this sense the situation is reminiscent of the handful of globular clusters that sit near the heart of the Sagittarius dwarf galaxy, a system that is currently being accreted into the Milky Way, confirming that accretion deposits not only stars but also globular clusters into the halos of large galaxies.

We present radial velocities and [Fe/H] abundances for 340 stars in the Fornax dwarf spheroidal from R~16,000 spectra. The targets have been obtained in the outer parts of the galaxy, a region which has been poorly studied before. Our sample shows a wide range in [Fe/H], between -0.5 and -3.0 dex, in which we detect three subgroups. Removal of stars belonging to the most metal-rich population produces a truncated metallicity distribution function that is identical to Sculptor, indicating that these systems have shared a similar early evolution, only that Fornax experienced a late, intense period of star formation (SF). The derived age-metallicity relation shows a fast increase in [Fe/H] at early ages, after which the enrichment flattens significantly for stars younger than ~8 Gyr. Additionally, the data indicate a strong population of stars around 4 Gyr, followed by a second rapid enrichment in [Fe/H]. A leaky-box chemical enrichment model generally matches the observed relation but does not predict a significant population of young stars nor the strong enrichment at late times. The young population in Fornax may therefore originate from an externally triggered SF event. Our dynamical analysis reveals an increasing velocity dispersion with decreasing [Fe/H] from sigma_sys 7.5 km/s to >14 km/s, indicating an outside-in star formation history in a dark matter dominated halo. The large velocity dispersion at low metallicities is possibly the result of a non-Gaussian velocity distribution amongst stars older than ~8 Gyr. Our sample also includes members from the Fornax GCs H2 and H5. In agreement with past studies we find [Fe/H]=-2.04+-0.04 and a mean radial velocity RV=59.36+-0.31 km/s for H2 and [Fe/H]=-2.02+-0.11 and RV=59.39+-0.44 km/s for H5. Overall, we find large complexity in the chemical and dynamical properties, with signatures that additionally vary with galactocentric distance.

We present a detailed kinematic analysis of the outer halo globular cluster (GC) system of M31. Our basis for this is a set of new spectroscopic observations for 78 clusters lying at projected distances between Rproj ~20-140 kpc from the M31 centre. These are largely drawn from the recent PAndAS globular cluster catalogue; 63 of our targets have no previous velocity data. Via a Bayesian maximum likelihood analysis we find that GCs with Rproj > 30 kpc exhibit coherent rotation around the minor optical axis of M31, in the same direction as more centrally- located GCs, but with a smaller amplitude of 86+/-17 km s-1. There is also evidence that the velocity dispersion of the outer halo GC system decreases as a function of projected distance from the M31 centre, and that this relation can be well described by a power law of index ~ -0.5. The velocity dispersion profile of the outer halo GCs is quite similar to that of the halo stars, at least out to the radius up to which there is available information on the stellar kinematics. We detect and discuss various velocity correlations amongst subgroups of GCs that lie on stellar debris streams in the M31 halo. Many of these subgroups are dynamically cold, exhibiting internal velocity dispersions consistent with zero. Simple Monte Carlo experiments imply that such configurations are unlikely to form by chance, adding weight to the notion that a significant fraction of the outer halo GCs in M31 have been accreted alongside their parent dwarf galaxies. We also estimate the M31 mass within 200 kpc via the Tracer Mass Estimator, finding (1.2 - 1.6) +/- 0.2 10^{12}M_sun. This quantity is subject to additional systematic effects due to various limitations of the data, and assumptions built in into the TME. Finally, we discuss our results in the context of formation scenarios for the M31 halo.

We present alpha-element abundances of Mg, Si, and Ti for a large sample of field stars in two outer fields of the Fornax dwarf spheroidal galaxy (dSph), obtained with VLT/GIRAFFE (R~16,000). Due to the large fraction of metal-poor stars in our sample, we are able to follow the alpha-element evolution from [Fe/H]=-2.5 continuously to [Fe/H]=-0.7 dex. For the first time we are able to resolve the turnover from the Type II supernovae (SNe) dominated, alpha-enhanced plateau down to subsolar [alpha/Fe] values due to the onset of SNe Ia, and thus to trace the chemical enrichment efficiency of the galaxy. Our data support the general concept of an alpha-enhanced plateau at early epochs, followed by a well-defined "knee", caused by the onset of SNe Ia, and finally a second plateau with sub-solar [alpha/Fe] values. We find the position of this knee to be at [Fe/H]=-1.9 and therefore significantly more metal-poor than expected from comparison with other dSphs and standard evolutionary models. Surprisingly, this value is rather comparable to the knee in Sculptor, a dSph about 10 times less luminous than Fornax. Using chemical evolution models, we find that both the position of the knee as well as the subsequent plateau at sub-solar level can hardly be explained unless the galaxy experienced several discrete star formation events with a drastic variation in star formation efficiency, while a uniform star formation can be ruled out. One possible evolutionary scenario is that Fornax experienced one or several major accretion events from gas-rich systems in the past, so that its current stellar mass is not indicative of the chemical evolution environment at ancient times. If Fornax is the product of several smaller building blocks, this may also have implications of the understanding on the formation process of dSphs in general.

We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into stream-like populations and a smooth halo component. Fitting a three-dimensional halo model reveals that the most metal-poor populations ([Fe/H]<-1.7) are distributed approximately spherically (slightly prolate with ellipticity c/a=1.09+/-0.03), with only a relatively small fraction (42%) residing in discernible stream-like structures. The sphericity of the ancient smooth component strongly hints that the dark matter halo is also approximately spherical. More metal-rich populations contain higher fractions of stars in streams (86% for [Fe/H]>-0.6). The space density of the smooth metal-poor component has a global power-law slope of -3.08+/-0.07, and a non-parametric fit shows that the slope remains nearly constant from 30kpc to 300kpc. The total stellar mass in the halo at distances beyond 2 degrees is 1.1x10^10 Solar masses, while that of the smooth component is 3x10^9 Solar masses. Extrapolating into the inner galaxy, the total stellar mass of the smooth halo is plausibly 8x10^9 Solar masses. We detect a substantial metallicity gradient, which declines from [Fe/H]=-0.7 at R=30kpc to [Fe/H]=-1.5 at R=150kpc for the full sample, with the smooth halo being 0.2dex more metal poor than the full sample at each radius. While qualitatively in-line with expectations from cosmological simulations, these observations are of great importance as they provide a prototype template that such simulations must now be able to reproduce in quantitative detail.

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

We present results from a spectroscopic survey of the dwarf spheroidal And XXII and the two extended clusters EC1 and EC2. These three objects are candidate satellites of the Triangulum galaxy, M33, which itself is likely a satellite of M31. We use the DEep Imaging Multi-Object Spectrograph mounted on the Keck-II telescope to derive radial velocities for candidate member stars of these objects and thereby identify the stars that are most likely actual members. Eleven most probable stellar members (of 13 candidates) are found for AndXXII. We obtain an upper limit of sigma_v < 6.0 km s-1 for the velocity dispersion of AndXXII, [Fe/H] ~ -1.6 for its metallicity, and 255pc for the Plummer radius of its projected density profile. We construct a colour magnitude diagram for AndXXII and identify both the red giant branch and the horizontal branch. The position of the latter is used to derive a heliocentric distance to And XXII of 853 pm 26 kpc. The combination of the radial velocity, distance, and angular position of AndXXII indicates that it is a strong candidate for being the first known satellite of M33 and one of the very few examples of a galactic satellite of a satellite. N-body simulations imply that this conclusion is unchanged even if M31 and M33 had a strong encounter in the past few Gyr. We test the hypothesis that the extended clusters highlight tidally stripped galaxies by searching for an excess cloud of halo-like stars in their vicinity. We find such a cloud for the case of EC1 but not EC2. The three objects imply a dynamical mass for M33 that is consistent with previous estimates.

We present the results of a spectroscopic survey performed in the outskirts of the globular cluster NGC1851 with VIMOS@VLT. The radial velocities of 107 stars in a region between 12' and 33' around the cluster have been derived. We clearly identify the cluster stellar population over the entire field of view, indicating the presence of a significant fraction of stars outside the tidal radius predicted by King models. We also find tentative evidence of a cold (sigma_v< 20 km/s) peak in the distribution of velocities at v_r~180 km/s constituted mainly by Main Sequence stars whose location in the color-magnitude diagram is compatible with a stream at a similar distance of this cluster. If confirmed, this evidence would strongly support the extra-Galactic origin of this feature.

The structure of the Sagittarius stream in the Southern Galactic hemisphere is analysed with the Sloan Digital Sky Survey Data Release 8. Parallel to the Sagittarius tidal track, but ~ 10deg away, there is another fainter and more metal-poor stream. We provide evidence that the two streams follow similar distance gradients but have distinct morphological properties and stellar populations. The brighter stream is broader, contains more metal-rich stars and has a richer colour-magnitude diagram with multiple turn-offs and a prominent red clump as compared to the fainter stream. Based on the structural properties and the stellar population mix, the stream configuration is similar to the Northern "bifurcation". In the region of the South Galactic Cap, there is overlapping tidal debris from the Cetus Stream, which crosses the Sagittarius stream. Using both photometric and spectroscopic data, we show that the blue straggler population belongs mainly to Sagittarius and the blue horizontal branch stars belong mainly to the Cetus stream in this confused location in the halo.

We develop, implement and characterise an enhanced data reduction approach which delivers precise, accurate, radial velocities from moderate resolution spectroscopy with the fibre-fed VLT/FLAMES+GIRAFFE facility. This facility, with appropriate care, delivers radial velocities adequate to resolve the intrinsic velocity dispersions of the very faint dSph dwarf galaxies. Importantly, repeated measurements let us reliably calibrate our individual velocity errors ($0.2 \leq \delta_V\leq 5$ km s$^{-1}$) and directly detect stars with variable radial velocities. We show, by application to the Bootes-1 dwarf spheroidal, that the intrinsic velocity dispersion of this system is significantly below 6.5\,km/s reported by previous studies. Our data favor a two-population model of Bootes-1, consisting of a majority `cold' stellar component, with velocity dispersion $2.4^{+0.9}_{-0.5}$\,km/s, and a minority `hot' stellar component, with velocity dispersion $\sim 9$\,km/s, although we can not completely rule out a single component distribution with velocity dispersion $4.6^{0.8}_{-0.6}$\,km/s. We speculate this complex velocity distribution actually reflects the distribution of velocity anisotropy in Bootes-1, which is a measure of its formation processes.

We present the results of deep imaging obtained at the CFHT with MegaCam in the Anticenter direction at two different heights above the Galactic disk. We detect the presence of the Monoceros ring in both fields as a conspicuous and narrow Main Sequence feature which dominates star counts over a large portion of the color-magnitude diagram down to g'~24. The comparison of the morphology and density of this feature with a large variety of Galactic models excludes the possibility that it can be due to a flare of the Galactic disk, supporting an extra-Galactic origin for this ring-like structure.

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.

We report the detection of a pair of degree-long tidal tails associated with the globular cluster Palomar 14, using images obtained at the CFHT. We reveal a power-law departure from a King profile at large distances to the cluster center. The density map constructed with the optimal matched filter technique shows a nearly symmetrical and elongated distribution of stars on both sides of the cluster, forming a S-shape characteristic of mass loss. This evidence may be the telltale signature of tidal stripping in action. This, together with its large Galactocentric distance, imposes strong constraints on its orbit and/or origin: i) it must follow an external orbit confined to the peripheral region of the Galactic halo and/or ii) it formed in a satellite galaxy later accreted by the Milky Way.

[Abridged] Within the hierarchical framework for galaxy formation, minor merging and tidal interactions are expected to shape all large galaxies to the present day. As a consequence, most seemingly normal disk galaxies should be surrounded by spatially extended stellar 'tidal features' of low surface brightness. As part of a pilot survey for such interaction signatures, we have carried out ultra deep, wide field imaging of 8 isolated spiral galaxies in the Local Volume, with data taken at small (D=0.1-0.5m) robotic telescopes that provide exquisite surface brightness sensitivity (mu_V)~28.5$ mag/arcsec^2). This initial observational effort has led to the discovery of six previously undetected extensive (to ~30 kpc) stellar structures in the halos surrounding these galaxies, likely debris from tidally disrupted satellites. In addition, we confirm and clarify several enormous stellar over-densities previously reported in the literature, but never before interpreted as tidal streams. Even this pilot sample of galaxies exhibits strikingly diverse morphological characteristics of these extended stellar features: great circle-like features that resemble the Sagittarius stream surrounding the Milky Way, remote shells and giant clouds of presumed tidal debris far beyond the main stelar body, as well as jet-like features emerging from galactic disks. A qualitative comparison with available simulations set in a Lambda-Cold Dark Matter cosmology shows that the extraordinary variety of stellar morphologies detected in this pilot survey matches that seen in those simulations. The common existence of these tidal features around 'normal' disk galaxies and the morphological match to the simulations constitutes new evidence that these theoretical models also apply to a large number of other Milky Way-mass disk galaxies in the Local Volume.

We use N-body simulations to explore the possibility that the Sagittarius (Sgr) dwarf galaxy was originally a late-type, rotating disc galaxy, rather than a non-rotating, pressure-supported dwarf spheroidal galaxy, as previously thought. We find that bifurcations in the leading tail of the Sgr stream, similar to those detected by the SDSS survey, naturally arise in models where the Sgr disc is misaligned with respect to the orbital plane. Moreover, we show that the internal rotation of the progenitor may strongly alter the location of the leading tail projected on the sky, and thus affect the constraints on the shape of the Milky Way dark matter halo that may be derived from modelling the Sgr stream. Our models provide a clear, easily-tested prediction: although tidal mass stripping removes a large fraction of the original angular momentum in the progenitor dwarf galaxy, the remnant core should still rotate with a velocity amplitude ~20 km/s that could be readily detected in future, wide-field kinematic surveys of the Sgr dwarf.

We use analytical and N-body methods to examine the survival of wide stellar binaries against repeated encounters with dark substructures orbiting in the dark matter haloes of dwarf spheroidal galaxies (dSphs). Our models adopt cosmologically-motivated conditions wherein dSphs are dark-matter dominated systems that form hierarchically and orbit about a host galaxy. Our analytical estimates show that wide binaries are disrupted at a rate that is proportional to the local density of dark substructures averaged over the life-time of the binary population. The fact that external tides can efficiently strip dark substructures from the outskirts of dSphs implies that the present number and distribution of binaries is strongly coupled with the mass evolution of individual galaxies. Yet we show that for the range of dynamical masses and Galactocentric distances spanned by Milky Way dSphs, a truncation in the separation function at a_max <~ 0.1 pc is expected in all these galaxies. An exception may be the Sagittarius dSph, which has lost most of is dark matter envelope to tides and is close to full disruption. Our simulations indicate that at separations larger than a_max the perturbed binary distribution scales as dN/da \propto a^{-2.1} independently of the mass and density of substructures. These results may be used to determine whether the binary separation function found in dwarf galaxies is compatible with the scale-free hierarchical picture that envisions the existence of dark substructures in all galactic haloes. We show that the ACS camera on board of the Hubble telescope may be able to test this prediction in dSphs at heliocentric distances <100 kpc, even if the binary fraction amounts only 10% of the stellar population.

(Abridged) We use N-body simulations to study the effects that a divergent (i.e. "cuspy") dark matter (DM) profile introduces on the tidal evolution of dwarf spheroidal galaxies (dSphs). Our models assume cosmologically-motivated initial conditions where dSphs are DM-dominated systems on eccentric orbits about a host galaxy composed of a dark halo and a baryonic disc. We find that the resilience of dSphs to tidal stripping is extremely sensitive to the halo cuspiness; whereas dwarfs with a cored profile can be easily destroyed by the host disc, those with cusps always retain a bound remnant. For a given halo profile the evolution of the structural parameters as driven by tides is controlled solely by the total amount of mass lost. This information is used to construct a semi-analytic code that simulates the hierarchical build-up of spiral galaxies assuming different halo profiles and disc masses. We find that tidal encounters with discs tend to decrease the average mass of satellites at all galactocentric radii. Interestingly, satellites accreted before re-ionization (z>6), which may be singled out by anomalous metallicity patterns, survive only if haloes are cuspy. We show that the size-mass relation established from Milky Way (MW) dwarfs strongly supports the presence of cusps in the majority of these systems, as cored models systematically underestimate the masses of the known Ultra-Faint dSphs. Our models also indicate that a massive M31 disc may explain why many of its dSphs fall below the size-mass relationship derived from MW dSphs. We use our models to constrain the mass threshold below which star formation is suppressed in DM haloes, finding that luminous satellites must be accreted with masses above 10^8--10^9 M_sol in order to explain the size-mass relation observed in MW dwarfs.

We present the first spectroscopic analysis of the faint M31 satellite galaxies, AndXI and AndXIII, and a reanalysis of existing spectroscopic data for two further faint companions, And IX and AndXII. By combining data obtained using the DEIMOS spectrograph mounted on the Keck II telescope with deep photometry from the Suprime-Cam instrument on Subaru, we have calculated global properties for the dwarfs, such as systemic velocities, metallicites and half-light radii.We find each dwarf to be very metal poor ([Fe/H] -2 both photometrically and spectroscopically, from their stacked spectrum), and as such, they continue to follow the luminosity-metallicity relationship established with brighter dwarfs. We are unable to resolve a dispersion for And XI due to small sample size and low S/N, but we set a one sigma upper limit of sigma-v <5 km/s. For And IX, And XII and And XIII we resolve velocity dispersions of v=4.5 (+3.4,-3.2), 2.6(+5.1,-2.6) and 9.7(+8.9,-4.5) km/s, and derive masses within the half light radii of 6.2(+5.3,-5.1)x10^6 Msun, 2.4 (+6.5,-2.4)x10^6 Msun and 1.1(+1.4,-0.7)x10^7 Msun respectively. We discuss each satellite in the context of the Mateo relations for dwarf spheroidal galaxies, and the Universal halo profiles established for Milky Way dwarfs (Walker et al. 2009). For both galaxies, this sees them fall below the Universal halo profiles of Walker et al. (2009). When combined with the findings of McConnachie & Irwin (2006a), which reveal that the M31 satellites are twice as extended (in terms of both half-light and tidal radii) as their Milky Way counterparts, these results suggest that the satellite population of the Andromeda system could inhabit halos that are significantly different from those of the Milky Way in terms of their central densities (abridged).

What is the mass of the progenitor of the Sagittarius (Sgr) dwarf galaxy? Here, we reassemble the stellar debris using SDSS and 2MASS data to find the total luminosity and likely mass. We find that the luminosity is in the range 9.6-13.2 x10^7 solar luminosities or M_V ~ -15.1 - 15.5, with 70% of the light residing in the debris streams. The progenitor is somewhat fainter than the present-day Small Magellanic Cloud, and comparable in brightness to the M31 dwarf spheroidals NGC 147 and NGC 185. Using cosmologically motivated models, we estimate that the mass of Sgr's dark matter halo prior to tidal disruption was ~10^10 solar masses.

We use N-body simulations to study the tidal evolution of globular clusters (GCs) in dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated scenario in which the dSph is approximated by a static NFW halo with a triaxial shape. For a large set of orbits and projection angles we examine the spatial and velocity distribution of stellar debris deposited during the complete disruption of stellar clusters. Our simulations show that such debris appears as shells, isolated clumps and elongated over-densities at low surface brightness (>26 mag/arcsec^2), reminiscent of substructure observed in several MW dSphs. Such features arise from the triaxiality of the galaxy potential and do not dissolve in time. Stellar over-densities reported in several MW dSphs may thus be the telltale evidence of dark matter haloes being triaxial in shape. We explore a number of kinematic signatures that would help to validate (or falsify) this scenario. The mean angular momentum of the cluster debris associated with box and resonant orbits, which are absent in spherical potentials, is null. As a result, we show that the line-of-sight velocity distribution may exhibit a characteristic "double-peak" depending on the oriention of the viewing angle with respect to the progenitor's orbital plane. Kinematic surveys of dSphs may help to detect and identify substructures associated with the disruption of stellar clusters, as well as to address the shape of the dark matter haloes in which dSphs are embedded.

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

We have conducted a deep multi-wavelength analysis (0.15-160 mum) to study the outer region of the nearby galaxy M94. We show that the non-optical data support the idea that the outskirts of this galaxy is not formed by a closed stellar ring (as traditionally claimed in the literature) but by a spiral arm structure. In this sense, M94 is a good example of a Type III (anti-truncated) disk galaxy having a very bright outer disk. The outer disk of this galaxy contains ~23% of the total stellar mass budget of the galaxy and contributes ~10% of the new stars created showing that this region of the galaxy is active. In fact, the specific star formation rate of the outer disk (~0.012 Gyr^{-1}) is a factor of ~2 larger (i.e. the star formation is more efficient per unit stellar mass) than in the inner disk. We have explored different scenarios to explain the enhanced star formation in the outer disk. We find that the inner disk (if considered as an oval distortion) can dynamically create a spiral arm structure in the outer disk which triggers the observed relatively high star formation rate as well as an inner ring similar to what is found in this galaxy.

We use N-body simulations to study the tidal evolution of globular clusters (GCs) in dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated scenario in which the dSph is approximated by a static NFW halo with a triaxial shape. We apply our models to five GCs spanning three orders of magnitude in stellar density and two in mass, chosen to represent the properties exhibited by the five GCs of the Fornax dSph. We show that only the object representing Fornax's least dense GC (F1) can be fully disrupted by Fornax's internal tidal field--the four denser clusters survive even if their orbits decay to the centre of Fornax. For a large set of orbits and projection angles we examine the spatial and velocity distribution of stellar debris deposited during the complete disruption of an F1-like GC. Our simulations show that such debris appears as shells, isolated clumps and elongated over-densities at low surface brightness (>26 mag/arcsec^2), reminiscent of substructure observed in several MW dSphs. Such features arise from the triaxiality of the galaxy potential and do not dissolve in time. The kinematics of the debris depends strongly on the progenitor's orbit. Debris associated with box and resonant orbits does not display stream motions and may appear "colder"/"hotter" than the dSph's field population if the viewing angle is perpendicular/parallel to progenitor's orbital plane. In contrast, debris associated with loop orbits shows a rotational velocity that may be detectable out to few kpc from the galaxy centre. Chemical tagging that can distinguish GC debris from field stars may reveal whether the merger of GCs contributed to the formation of multiple stellar components observed in dSphs.

(Abridged) We use N-body simulations to explore the effects of tidal stripping on the structure of dwarf spheroidal galaxies (dSphs). We model dSphs as King spheres embedded in NFW dark halos that orbit the Galactic potential on eccentric orbits. We find that episodes of stellar mass loss remove the initial cutoff of the bound stellar core. Once equilibrium has been re-established the outer mass profile approaches a power-law in projection that is well described by a simple Plummer model. Tides also result in transient features in the outer density profile. As the system re-equilibrates, an outward-moving "excess" of stars is found at radii where the local crossing time exceeds the time elapsed since pericenter. If the orbit is known, these results provide a simple way to assess whether "breaks" and "bumps" in the profile of dSphs are actually tidal in origin. We apply this to the Sagittarius dwarf and, encouragingly, identify two features in the profile that may be traced to its two last pericentric passages. Applied to Leo I, we find that any tidal break would locate beyond the area surveyed by current data, casting doubt on recent claims of the detection of tidal debris around this galaxy. For Carina, the tidal break should occur at a radius twice farther than observed, suggesting that the outer excess of stars is not tidal in origin unless its orbit is in error. A similar comment applies to Sculptor, whose pericenter appears too large for Galactic tides to be important but whose outer profile, like that of Draco, nonetheless follows closely a Plummer-law. Fornax and Leo II show no sign of a power-law outer profile, suggesting that they have not suffered stellar tidal stripping. Published profiles for other, fainter Milky Way dSph companions do not extend sufficiently far to allow for conclusive assessment.

Within the hierarchical framework for galaxy formation, merging and tidal interactions are expected to shape large galaxies to this day. While major mergers are quite rare at present, minor mergers and satellite disruptions - which result in stellar streams - should be common, and are indeed seen in both the Milky Way and the Andromeda Galaxy. As a pilot study, we have carried out ultra-deep, wide-field imaging of some spiral galaxies in the Local Volume, which has revealed external views of such stellar tidal streams at unprecedented detail, with data taken at small robotic telescopes (0.1-0.5-meter) that provide exquisite surface brightness sensitivity. The goal of this project is to undertake the first systematic and comprehensive imaging survey of stellar tidal streams, from a sample of ~50 nearby Milky-Way-like spiral galaxies within 15 Mpc, that features a surface brightness sensitivity of ~ 30 mag/arcsec^2 The survey will result in estimates of the incidence, size/geometry and stellar luminosity/mass distribution of such streams. This will not only put our Milky Way and M31 in context but, for the first time, also provide an extensive statistical basis for comparison with state-of-the-art, self-consistent cosmological simulations of this phenomenon.

We report the discovery of a giant, loop-like stellar structure around the edge-on spiral galaxy NGC 4013. This arcing feature extends 6 arcmin (~26 kpc in projected distance) northeast from the center and 3 arcmin (~=12 kpc) from the disk plane; likely related features are also apparent on the southwest side of the disk, extending to 4 arcmin (~17 kpc). The detection of this low surface-brightness muR= 27.0+0.3-0.2 mag/sqarcsec) structure is independently confirmed in three separate datasets from three different telescopes. Although its true three dimensional geometry is unknown, the sky- projected morphology of this structure displays a match with the theoretical predictions for the edge-on, projected view of a stellar tidal streams of a dwarf satellite moving in a low inclined (~25deg), nearly circular orbit. Using the recent model of the Monoceros tidal stream in the Milky Way by Penarrubia et al. as template, we find that the progenitor system may have been a galaxy with an initial mass 6*10^8 Msun, of which current position and final fate is unknown. According to this simulation, the tidal stream may be approximately ~2.8 Gyr of age. Our results demonstrate that NGC 4013, previously considered a prototypical isolated disk galaxy in spite of having one of the most prominent HI warps detected thus far, may have in fact suffered a recent minor merger. This discovery highlights that undisturbed disks at high surface brightness levels in the optical but warped in HI maps may in fact reveal complex signatures of recent accretion events in deep photometric surveys.

We present an extragalactic perspective of an extended stellar tidal stream wrapping around the edge-on, spiral galaxy NGC 5907. Our deep images reveal for the first time a large scale complex of arcing loops that is an excellent example of how a low-mass satellite accretion can produce an interwoven, rosette-like structure of debris dispersed in the halo of its host galaxy. The existence of this structure, which has probably formed and survived for several Gigayears, confirms that halos of spiral galaxies in the Local Universe may still contain a significant number of galactic fossils from their hierarchical formation. To examine the validity of the external accretion scenario, we present N-body simulations of the tidal disruption of a dwarf galaxy-like system in a disk galaxy plus dark halo potential that demonstrate that most of the observed tidal features observed in NGC 5907 can be explained by a single accretion event. Unfortunately, with no kinematic data and only the projected geometry of the stream as constraint, the parameters of our model are considerably degenerate and, for now, must be considered illustrative only. Interestingly, NGC 5907 has long been considered a prototypical example of a warped spiral in relative isolation. The presence of an extended tidal stream challenges this picture and suggests that the gravitational perturbations induced by the stream progenitor may be the cause for the warp. The detection of an old, complex tidal stream in a nearby galaxy with rather modest instrumentation points to the viability of surveys to find extragalactic tidal substructures around spiral galaxies in the Local Volume (< 15 Mpc) -- with the prospect of obtaining a census with enough statistical significance to be compared with cosmological simulations.

(abridged) We present a spectroscopic analysis of five stellar streams (`A', `B', `Cr', `Cp' and `D') as well as the extended star cluster, EC4, which lies within streamC, all discovered in the halo of M31 from our CFHT/MegaCam survey. These spectroscopic results were initially serendipitous, making use of our existing observations from the DEep Imaging Multi-Object Spectrograph mounted on the Keck II telescope, and thereby emphasizing the ubiquity of tidal streams that account for ~70% of the M31 halo stars in the targeted fields. Subsequent spectroscopy was then procured in streamCr/p and streamD to trace the velocity gradient along the streams. For the cluster EC4, candidate member stars with average [Fe/H]~-1.4 (Fe/H_spec=-1.6), are found at v_{hel}=-285 km/s suggesting it could be related to streamCp. No similarly obvious cold kinematic candidate is found for streamD, although candidates are proposed in both of two spectroscopic pointings along the stream (both at -400 km/s). Spectroscopy near the edge of streamB suggests a likely kinematic detection, while a candidate kinematic detection of streamA is found (plausibly associated to M33 rather than M31). The low dispersion of the streams in kinematics, physical thickness, and metallicity makes it hard to reconcile with a scenario whereby these stream structures as an ensemble are related to the giant southern stream. We conclude that the M31 stellar halo is largely made up of multiple kinematically cold streams.

We present the second data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities (RVs) and stellar atmosphere parameters of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). It is obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region (8,410--8,795 \AA) for southern hemisphere stars in the magnitude range 9<I<12. Following the first data release (Steinmetz et al. 2006) the current release doubles the sample of published RVs, now containing 51,829 RVs for 49,327 individual stars observed on 141 nights between April 11 2003 and March 31 2005. Comparison with external data sets shows that the new data collected since April 3 2004 show a standard deviation of 1.3 km/s, about twice better than for the first data release. For the first time this data release contains values of stellar parameters from 22,407 spectra of 21,121 individual stars. They were derived by a penalized \chi^2 method using an extensive grid of synthetic spectra calculated from the latest version of Kurucz models. From comparison with external data sets, our conservative estimates of errors of the stellar parameters (for a spectrum with S/N=40) are 400 K in temperature, 0.5 dex in gravity, and 0.2 dex in metallicity. We note however that the internal errors estimated from repeat RAVE observations of 822 stars are at least a factor 2 smaller. We demonstrate that the results show no systematic offsets if compared to values derived from photometry or complementary spectroscopic analyses. The data release includes proper motion and photometric measurements. It can be accessed via the RAVE webpage: http://www.rave-survey.org and through CDS.

We have searched for in-falling stellar streams on to the local Milky Way disc in the CORAVEL and RAVE surveys. The CORAVEL survey consists of local dwarf stars (Nordstrom et al. Geneva-Copenhagen survey) and local Famaey et al. giant stars. We select RAVE stars with radial velocities that are sensitive to the Galactic vertical space velocity (Galactic latitude b < -45 deg). Kuiper statistics have been employed to test the symmetry of the Galactic vertical velocity distribution functions in these samples for evidence of a net vertical flow that could be associated with a (tidal?) stream of stars with vertically coherent kinematics. In contrast to the `Field of Streams' found in the outer halo, we find that the local volumes of the solar neighbourhood sampled by the CORAVEL dwarfs (complete within ~3 x 10^-4 kpc^3), CORAVEL giants (complete within ~5 x 10^-2 kpc^3) and RAVE (5-15% complete within ~8 kpc^3) are devoid of any vertically coherent streams containing hundreds of stars. This is sufficiently sensitive to allow our RAVE sample to rule out the passing of the tidal stream of the disrupting Sagittarius (Sgr) dwarf galaxy through the solar neighbourhood. This agrees with the most recent determination of its orbit and dissociates it from the Helmi et al. halo stream. Our constraints on the absence of the Sgr stream near the Sun could prove a useful tool for discriminating between Galactic potential models. The lack of a net vertical flow through the solar neighbourhood in the CORAVEL giants and RAVE samples argues against the Virgo overdensity crossing the disc near the Sun. There are no vertical streams in the CORAVEL giants and RAVE samples with stellar densities >1.6 x 10^4 and 1.5 x 10^3 stars kpc^-3 respectively and therefore no evidence for locally enhanced dark matter.

(Abridged) We use N-body simulations to study the evolution of dwarf spheroidal galaxies (dSphs) driven by galactic tides. We adopt a cosmologically-motivated model where dSphs are approximated by a King model embedded within an NFW halo. We find that these NFW-embedded King models are extraordinarily resilient to tides; the stellar density profile still resembles a King model even after losing more than 99% of the stars. As tides strip the galaxy, the stellar luminosity, velocity dispersion, central surface brightness, and core radius decrease monotonically. Remarkably, we find that the evolution of these parameters is solely controlled by the total amount of mass lost from within the luminous radius. Of all parameters, the core radius is the least affected: after losing 99% of the stars, R_c decreases by just a factor of ~2. Interestingly, tides tend to make dSphs more dark-matter dominated because the tightly bound central dark matter ``cusp'' is more resilient to disruption than the ``cored'' King profile. We examine whether the extremely large M/L ratios of the newly-discovered ultra-faint dSphs might have been caused by tidal stripping of once brighter systems. Although dSph tidal evolutionary tracks parallel the observed scaling relations in the luminosity-radius plane, they predict too steep a change in velocity dispersion compared with the observational estimates hitherto reported in the literature. The ultra-faint dwarfs are thus unlikely to be the tidal remnants of systems like Fornax, Draco, or Sagittarius. Despite spanning four decades in luminosity, dSphs appear to inhabit halos of comparable peak circular velocity, lending support to scenarios that envision dwarf spheroidals as able to form only in halos above a certain mass threshold.

We examine the dynamics of stellar systems embedded within cold dark matter (CDM) halos in order to assess observational constraints on the dark matter content of Local Group dwarf spheroidals (dSphs). Our analysis shows that the total mass within the luminous radius is reasonably well constrained and approximately independent of the luminosity of the dwarf, highlighting the poor correspondence between luminosity and halo mass. This result implies that the average density of dark matter is substantially higher in physically small systems such as Draco and Sculptor than in larger systems such as Fornax. For example, our results imply that Draco formed in a halo 5 times more massive than Fornax's despite being roughly 70 times fainter. Stellar velocity dispersion profiles, sigma_p(R), provide further constraints; flat sigma_p(R) profiles imply that stars are deeply embedded within their cold dark matter halos and so quite resilient to tidal disruption. We estimate that halos would need to lose more than 90% of their original mass before tides begin affecting the kinematics of stars. We estimate that V_max is about 3 times higher than the central velocity dispersion of the stars, which alleviates significantly the CDM ``substructure crisis''. We use these results to interpret the size differences between the M31 and Milky Way (MW) dSph population. Our modeling indicates that this difference should be reflected in their kinematics, and predicts that M31 dwarfs should have velocity dispersions up to a factor of ~ 2 higher than their MW counterparts. This CDM-motivated prediction may be verified with present observational capabilities.

The dwarf spheroidal (dSph) satellites of the Local Group have long been thought to be simple spheroids of stars embedded within extended dark matter halos. Recently, however, evidence for the presence of spatially and kinematically distinct stellar populations has been accumulating. Here, we examine the influence of such components on dynamical models of dwarf galaxies embedded in cold dark matter halos. We begin by constructing a model of Andromeda II, a dSph satellite of M31 which shows evidence for spatially distinct stellar components. We find that the two-component model predicts an overall velocity dispersion profile that remains approximately constant at $\sim 10 - 11$ km s$^{-1}$ out to $\sim 1$ kpc from the center; this is despite wide kinematic and spatial differences between the two individual components. The presence of two components may also help to explain oddities in the velocity dispersion profiles of other dSphs; we show that velocity dispersion profiles which appear to rise from the center outwards before leveling off--such as those of Leo I, Draco, and Fornax--can result from the gradual transition from a dynamically cold, concentrated component to a second, hotter, and more spatially extended one, both in equilibrium within the same dark halo. Dwarf galaxies with two stellar components generally have a leptokurtic line-of-sight velocity distribution which is well described by a double Maxwellian. Interestingly, we find that multiple equilibrium components could also provide a potential alternative origin for ``extra-tidal'' stars (normally ascribed to tidal effects) in situations where corroborating evidence for tides may be lacking.

We present spectroscopic observations of the AndXII dwarf spheroidal galaxy using DEIMOS/Keck-II, showing it to be moving rapidly through the Local Group (-556 km/s heliocentric velocity, -281 km/s relative to Andromeda from the MW), falling into the Local Group from ~115 kpc beyond Andromeda's nucleus. AndXII therefore represents a dwarf galaxy plausibly falling into the Local Group for the first time, and never having experienced a dense galactic environment. From Green Bank Telescope observations, a limit on the H{I} gas mass of <3000 Msun suggests that AndXII's gas could have been removed prior to experiencing the tides of the Local Group galaxies. Orbit models suggest the dwarf is close to the escape velocity of M31 for published mass models. AndXII is our best direct evidence for the late infall of satellite galaxies, a prediction of cosmological simulations.

The recently discovered Virgo stellar over-density, which expands over \~1000deg^2 perpendicularly to the Galactic disk plane (7< Z <15 kpc, R~7 kpc), is the largest clump of tidal debris ever detected in the outer halo and is likely related with the accretion of a nearby dwarf galaxy by the Milky Way. We carry out N-body simulations of the Sagittarius stream to show that this giant stellar over-density is a confirmation of theoretical model predictions for the leading tail of the Sagittarius stream to cross the Milky Way plane in the Solar neighborhood. Radial velocity measurements are needed to confirm this association and to further constrain the shape of the Milky Way dark matter halo through a new generation of theoretical models. If the identification of Virgo over-density and the Sagittarius leading arm is correct, we predict highly negative radial velocities for the stars of Virgo over-density. The detection of this new portion of the Sagittarius tidal stream would represent an excellent target for the on-going and future kinematic surveys and for dark matter direct detection experiments in the proximity of the Sun.

[Abridged] We present first results from a large-area (~80degx20deg), sparsely sampled two-filter (B,R) imaging survey towards the Canis Major stellar over-density, claimed to be a disrupting Milky Way satellite galaxy. Utilizing stellar colour-magnitude diagrams reaching to B ~ 22 mag, we provide a first delineation of its surface density distribution using main sequence stars. Its projected shape is highly elongated, nearly parallel to the Galactic plane, with an axis ratio of >~ 5:1, substantially more so than what Martin et al. originally found. We also provide a first map of a prominent over-density of blue, presumably younger main sequence stars, which appears to have a maximum near [l,b = 240,-7 deg]. The young population is markedly more localized. We estimate an upper limit on the line-of-sight (l.o.s.) depth of the old population based on the main sequence width, obtaining sigma_los < 1.8 +/- 0.3 kpc, at an adopted D_helio = 7.5 +/- 1 kpc. For the young stellar population, we find sigma_los < 1.5 kpc. There are different explanations for the CMa over-density: (a) a partially disrupting dwarf galaxy on a low-latitude orbit, (b) a projection of the warped outer Galactic disk, and (c) a projection of an out-of-plane spiral arm. While the data provide no firm arguments against the less well-defined third scenario, they have clear implications for each of the others: (a) We infer from the strong elongation in longitude, and simulations in the literature, that the over-density is unlikely to be a gravitationally bound system at the present epoch, but may well be just a recently disrupted satellite remnant. (b) Based on modeling, the line-of-sight depth of the over-density in old stars is clearly inconsistent with published locally axi-symmetric descriptions of the warped Galactic disk.

We explore the possibility that extended disks, such as that recently discovered in M31, are the result of a single dwarf (10^9--10^10 Msun) satellite merger. We carry out N-body simulations of dwarf NFW halos with embedded spheriodal stellar components on co-planar, prograde orbits in a M31-like host galaxy. As the orbit decays due to dynamical friction and the system is disrupted, the stellar particles form an extended, exponential disk-like structure that spans the radial range 30--200 kpc once they have relaxed. The disk scale-length correlates with the initial extent of the stellar component within the satellite halo: the more embedded the stars, the smaller the resulting disk scale-length. If the progenitors start on circular orbits, the kinematics of the stars that make up the extended disk have an average rotational motion that is 30--50 km/s lower than the host's circular velocity. For dwarf galaxies moving on highly eccentric orbits (e=0.7), the stellar debris exhibits a much lower rotational velocity. Our results imply that extended galactic disks might be a generic feature of the hierarchical formation of spiral galaxies such as M31 and the Milky Way.

We present the first data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities and stellar atmosphere parameters (temperature, metallicity, surface gravity) of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). The RAVE program started in 2003, obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region ($\lambda\lambda$ 8,410--8,795 \AA) for southern hemisphere stars drawn from the Tycho-2 and SuperCOSMOS catalogs, in the magnitude range 9<I<12. The first data release is described in this paper and contains radial velocities for 24,748 individual stars (25,274 measurements when including re-observations). Those data were obtained on 67 nights between 11 April 2003 to 03 April 2004. The total sky coverage within this data release is $\sim$4,760 square degrees. The average signal to noise ratio of the observed spectra is 29.5, and 80% of the radial velocities have uncertainties better than 3.4 km/s. Combining internal errors and zero-point errors, the mode is found to be 2 km/s. Repeat observations are used to assess the stability of our radial velocity solution, resulting in a variance of 2.8 km/s. We demonstrate that the radial velocities derived for the first data set do not show any systematic trend with color or signal to noise. The RAVE radial velocities are complemented in the data release with proper motions from Starnet 2.0, Tycho-2 and SuperCOSMOS, in addition to photometric data from the major optical and infrared catalogs (Tycho-2, USNO-B, DENIS and 2MASS). The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

We explore whether stellar tidal streams can provide information on the secular, cosmological evolution of the Milky Way's gravitational potential and on the presence of subhalos. We carry out long-term (~t_hubble) N-body simulations of disrupting satellite galaxies in a semi-analytic Galaxy potential where the dark matter halo and the subhalos evolve according to a LCDM cosmogony. All simulations are constrained to end up with the same position and velocity at present. Our simulations account for: (i) the secular evolution of the host halo's mass, size and shape, (ii) the presence of subhalos and (iii) dynamical friction. We find that tidal stream particles respond adiabatically to the Galaxy growth so that, at present, the energy and angular momentum distribution is exclusively determined by the present Galaxy potential. In other words, all present-day observables can only constrain the present mass distribution of the Galaxy independent of its past evolution. We also show that, if the full phase-space distribution of a tidal stream is available, we can accurately determine (i) the present Galaxy's shape and (ii) the amount of mass loss from the stream's progenitor, even if this evolution spanned a cosmologically significant epoch.

A luminosity-temperature relation for clusters of galaxies is derived. The two models used, take into account the angular momentum acquisition by the proto-structures during their expansion and collapse. The first one is a modification of the self-similar model (SSM) while the second one is a modification of the Punctuated Equilibria Model (Cavaliere et al. 1999). In both models the mass-temperature relation (M-T) used is based on the calculations of Del Popolo (2002b). We show that the above models lead, in X-rays, to a luminosity-temperature relation that scales as L \propto T^5, at scale of groups, flattening to L \propto T^3 for rich clusters and converging to L \propto T^2 at higher temperatures. However a fundamental result of our paper is that the non-similarity in the L-T relation, can be explained by a simple model that takes into account the amount of the angular momentum of a proto-structure. This result is in disagreement with the widely accepted idea that the above non-similarity is due to non-gravitating processes as those of heating/cooling.

We have measured the absolute proper motion of the candidate Canis Major dwarf galaxy (CMa) at $(l,b) = (240\arcdeg, -8\arcdeg)$. Likely main-sequence stars in CMa have been selected from a region in the color-magnitude diagram that has very little contamination from known Milky Way components. We obtain $\mu_{l} cos b = -1.47 \pm 0.37$ and $\mu_{b} = -1.07 \pm 0.38$ mas yr$^{-1}$, on the ICRS system via Hipparcos stars. Together with the radial velocity of 109 km s$^{-1}$, and the assumed distance of 8 kpc, these results imply a space motion of $(\Pi, \Theta, W) = (-5 \pm 12, 188 \pm 10, -49 \pm 15)$ km s$^{-1}$. While CMa has in-plane rotation similar to the mean of thick disk stars, it shows significant $(3\sigma)$ motion perpendicular to the disk, and differs even more $(7\sigma)$ from that expected for the Galactic warp. The $W$ velocity lends support to the argument that the CMa overdensity is part of a satellite galaxy remnant.

The Sloan Digital Sky Survey has recently discovered a coherent ring of stars at low galactic latitude that is believed to be the tidal stream of a merging dwarf galaxy in the Galactic plane (named the Monoceros tidal stream). The existence and location of the core of its progenitor galaxy is still controversial. The best candidate is the Canis Major dwarf galaxy, a distinct overdensity of red stars discovered in the 2MASS survey, but also interpreted as the signature of the Galactic warp viewed in projection. In this paper, we report a variety of new observational evidence that supports the notion that CMa is the remnant of a partially disrupted core of a dwarf satellite. The comparison of the orbit derived from our theoretical model for the parent galaxy of this ring-like structure with an accurate determination of CMa orbit leads to the conclusion that this satellite is the best candidate for the progenitor of the Monoceros tidal stream.

We develop a semi-analytical model that determines the evolution of the mass, position and internal structure of dark matter substructures orbiting in dark matter haloes. We apply this model to the case of the Milky Way. We focus in particular on the effects of mass loss, dynamical friction and substructure--substructure interactions, the last of which has previously been ignored in analytic models of substructure evolution. Our semi-analytical treatment reproduces both the spatial distribution of substructures and their mass function as obtained from the most recent N-body cosmological calculations of Gao et al. (2004). We find that, if mass loss is taken into account, the present distribution of substructures is practically insensitive to dynamical friction and scatterings from other substructures. Implementing these phenomena leads to a slight increase (~5%) in the number of substructures at r<0.35 r_vir, whereas their effects on the mass function are negligible. We find that mass loss processes lead to the disruption of substructures before dynamical friction and gravitational scattering can significantly alter their orbits. Our results suggest that the present substructure distribution at r>0.35 r_vir reflects the orbital properties at infall and is, therefore, purely determined by the dark matter environment around the host halo and has not been altered by dynamical evolution.

We investigate the effect of a density gradient on Chandrasekhar's dynamical friction formula based on the method of 2-body encounters in the local approximation. We apply these generalizations to the orbit evolution of satellite galaxies in Dark Matter haloes. We find from the analysis that the main influence occurs through a position-dependent maximum impact parameter in the Coulomb logarithm, which is determined by the local scale-length of the density distribution. We also show that for eccentric orbits the explicit dependence of the Coulomb logarithm on position yields significant differences for the standard homogeneous force. Including the velocity dependence of the Coulomb logarithm yields ambigous results. The orbital fits in the first few periods are further improved, but the deviations at later times are much larger. The additional force induced by the density gradient, the inhomogeneous force, is not antiparallel to the satellite motion and can exceed 10% of the homogeneous friction force in magnitude. However, due to the symmetry properties of the inhomogeneous force, there is a deformation and no secular effect on the orbit at the first order. Therefore the inhomogeneous force can be safely neglected for the orbital evolution of satellite galaxies. For the homogeneous force we compare numerical N-body calculations with semi-analytical orbits to determine quantitatively the accuracy of the generalized formulae of the Coulomb logarithm in the Chandresekhar approach. With the local scale-length as the maximum impact parameter we find a significant improvement of the orbital fits and a better interpretation of the quantitative value of the Coulomb logarithm.

We have compiled an extensive dataset on potential parts of the Monoceros tidal stream and performed an exhaustive survey of dwarf galaxy semi-analytic orbits in order to constrain its orbital properties. The best-fit orbits are subsequently realized as self-consistent N-body simulations in order to reproduce the spatial and velocity distribution of satellite debris. We find that all kinematic and geometric constraints can be fit by a single stream allowing for multiple wraps. The orbital eccentricity and inclination of the progenitor are strongly constrained to be e = 0.10 \pm 0.05 and i = 25 \pm 5 deg. Ten new estimates of proper motions from the Sloan Digital Sky Survey (SDSS) clearly exclude all retrograde orbits. Particles lost by the satellite populate two, nearly concentric rings naturally explaining the detection of stream stars at both 6--8 kpc (Ibata et al. 2003, Newberg et al. 2003) and 12--18 kpc (the Tri/And stream; Rocha-Pinto et al. 2003) from the Sun. We have attempted to predict the present location of the Monoceros stream progenitor using different information: (i) the kinematical and spatial distribution of detections and (ii) the different mean metallicity in the innner and the outer rings. A current location at l ~ 245, b ~ -18., with a distance to the Sun r_s~15 kpc, reproduces the available constraints. The mass of the progenitor has been poorly determined due to the slow orbital decay. Similar fits have been obtained for masses (3-9)10^8 M_sun. Additionally, we have analyzed to possible common origin of the Canis Major dwarf and the Monoceros stream, unfortunately, without a conclusive result.

We carry out a set of self-consistent N-body calculations to investigate how important the velocity anisotropy in non-spherical dark-matter haloes is for dynamical friction. For this purpose we allow satellite galaxies to orbit within flattened and live dark-matter haloes (DMHs) with aspect-ratio q_h=0.6 and compare the resulting orbit evolution with a semi-analytic code. This code solves the equation of motion of satellites orbiting within the same DMH and implements, according to the dynamical friction model, either Chandrasekhar's formula, which does not incorporate the velocity anisotropy, or Binney's description of dynamical friction in anisotropic systems. In the numerical and the two semi-analytic models the satellites are given different initial orbital inclinations and orbital eccentricities. We find that Binney's approach successfully describes the overall satellite decay and orbital inclination decrease for the whole set of orbits, with an averaged discrepancy of less than 4 per cent in orbital radius during the first 3 orbits. If Chandrasekhar's expression is used instead, the discrepancy increases to 20 per cent and no inclination decrease can be observed. Binney's treatment therefore appears to provide a significantly improved treatment of dynamical friction in anisotropic systems.

We carry out a set of self-consistent N-body calculations to compare the decay rates of satellite dwarf galaxies orbiting a disc galaxy embedded in a dark matter halo (DMH). We consider both spherical and oblate axisymmetric DMHs of aspect ratio q_h=0.6. The satellites are given different initial orbital inclinations, orbital periods and mass. The live flattened DMHs with embedded discs and bulges are set-up using a new fast algorithm, MaGalie (Boily, Kroupa and Pe\~{n}arrubia 2001). We find that the range of survival times of satellites within a flattened DMH becomes of the order of 100% larger than the same satellites within a spherical DMH. In the oblate DMH, satellites on polar orbits have the longest survival time, whereas satellites on coplanar prograde orbits are destroyed most rapidly. The orbital plane of a satellite tilts as a result of anisotropic dynamical friction, causing the satellite's orbit to align with the plane of symmetry of the DMH. Polar orbits are not subjected to alignment. Therefore the decay of a satellites in an axisymmetric DMH may provide a natural explanation for the observed lack of satellites within (0-30) degrees of their host galaxy's disc (Holmberg 1969; Zaritsky and Gonz\'alez 1999). The computations furthermore indicate that the evolution of the orbital eccentricity $e$ is highly dependent of its initial value e(t=0) and the DMH's shape. We also discuss some implications of flattened DMHs for satellite debris streams.