We analyse a set of cosmological magneto-hydrodynamic simulations of the formation of Milky Way-mass galaxies identified to have a prominent radially anisotropic stellar halo component similar to the so-called "Gaia Sausage" found in the Gaia data. We examine the effects of the progenitor of the Sausage (the Gaia-Enceladus-Sausage, GES) on the formation of major galactic components analogous to the Galactic thick disc and inner stellar halo. We find that the GES merger is likely to have been gas-rich and contribute 10-50$\%$ of gas to a merger-induced centrally concentrated starburst that results in the rapid formation of a compact, rotationally supported thick disc that occupies the typical chemical thick disc region of chemical abundance space. We find evidence that gas-rich mergers heated the proto-disc of the Galaxy, scattering stars onto less-circular orbits such that their rotation velocity and metallicity positively correlate, thus contributing an additional component that connects the Galactic thick disc to the inner stellar halo. We demonstrate that the level of kinematic heating of the proto-galaxy correlates with the kinematic state of the population before the merger, the progenitor mass and orbital eccentricity of the merger. Furthermore, we show that the mass and time of the merger can be accurately inferred from local stars on counter-rotating orbits.

Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$\sigma$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

We use the astrometric and photometric data from Gaia Data Release 2 and line-of-sight velocities from various other surveys to study the 3d structure and kinematics of the Sagittarius dwarf galaxy. The combination of photometric and astrometric data makes it possible to obtain a very clean separation of Sgr member stars from the Milky Way foreground; our final catalogue contains ~2.6e5 candidate members with magnitudes G<18, more than half of them being red clump stars. We construct and analyze maps of the mean proper motion and its dispersion over the region ~30x12 degrees, which show a number of interesting features. The intrinsic 3d density distribution (orientation, thickness) is strongly constrained by kinematics; we find that the remnant is a prolate structure with the major axis pointing at 45deg from the orbital velocity and extending up to ~5 kpc, where it transitions into the stream. We perform a large suite of N-body simulations of a disrupting Sgr galaxy as it orbits the Milky Way over the past 2.5 Gyr, which are tailored to reproduce the observed properties of the remnant (not the stream). The richness of available constraints means that only a narrow range of parameters produce a final state consistent with observations. The total mass of the remnant is ~4e8 Msun, of which roughly a quarter resides in stars. The galaxy is significantly out of equilibrium, and even its central density is below the limit required to withstand tidal forces. We conclude that the Sgr galaxy is terminally wounded, and will be fully disrupted over the next Gyr.

Most binaries are undetected. Astrometric reductions of a system using the assumption that the object moves like a single point mass can be biased by unresolved binary stars. The discrepancy between the centre of mass of the system (which moves like a point mass) and the centre of light (which is what we observe) introduces additional motion. We explore the extent to which binary systems affect single object models fit to astrometric data. This tells us how observations are diluted by binaries and which systems cause the largest discrepancies - but also allows us to make inferences about the binarity of populations based on observed astrometric error. By examining a sample of mock observations, we show that binaries with periods close to one year can mimic parallax and thus bias distance measurements, whilst long period binaries can introduce significant apparent proper motion. Whilst these changes can soak up some of the error introduced by the binary, the total deviation from the best fitting model can be translated into a lower limit on the on-sky separation of the pair. Throughout we link these predictions to data from the Gaia satellite, whilst leaving the conclusions generalizable to other surveys.

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

Milky Way dwarf satellites are unique objects that encode the early structure formation and therefore represent a window into the high redshift Universe. So far, their study was conducted using electromagnetic waves only. The future Laser Interpreter Space Antenna (LISA) has the potential to reveal Milky Way satellites in gravitational waves emitted by double white dwarf (DWD) binaries. We investigate gravitational wave (GW) signals detectable by LISA as a possible tool for the identification and characterisation of the Milky Way satellites. We use the binary population synthesis technique to model the population of DWDs in dwarf satellites and we assess the impact on the number of LISA detections when making changes to the total stellar mass, distance, star formation history and metallicity of satellites. We calibrate predictions for the known Milky Way satellites on their observed properties. We find that DWDs emitting at frequencies $\gtrsim 3\,$mHz can be detected in Milky Way satellites at large galactocentric distances. The number of these high frequency DWDs per satellite primarily depends on its mass, distance, age and star formation history, and only mildly depends on the other assumptions regarding their evolution such as metallicity. We find that dwarf galaxies with $M_\star>10^6\,$M$_{\odot}$ can host detectable LISA sources with a number of detections that scales linearly with the satellite's mass. We forecast that out of the known satellites, Sagittarius, Fornax, Sculptor and the Magellanic Clouds can be detected with LISA. As an all-sky survey that does not suffer from contamination and dust extinction, LISA will provide observations of the Milky Way and dwarf satellites galaxies valuable for Galactic archaeology and near-field cosmology.

Using the data from the Sloan Digital Sky Survey and the Gaia satellite, we assemble a pure sample of $\sim$3000 Blue Horizontal Branch (BHB) stars with 7-D information, including positions, velocities and metallicities. We demonstrate that, as traced with BHBs, the Milky Way's stellar halo is largely unmixed and can not be well represented with a conventional Gaussian velocity distribution. A single-component model fails because the inner portions of the halo are swamped with metal-rich tidal debris from an ancient, head-on collision, known as the "Gaia Sausage". Motivated by the data, we build a flexible mixture model which allows us to track the evolution of the halo make-up across a wide range of radii. It is built from two components, one representing the radially anisotropic Sausage stars with their lobed velocity distribution, the other representing a more metal-poor and more isotropic component built up from minor mergers. We show that inside 25 kpc the "Sausage" contributes at least 50 % of the Galactic halo. The fraction of "Sausage" stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.

Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically-luminous standard candles and their distances can be inferred with small relative precision ($\sim3\%$). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, $d = 20.6 \pm 0.2~\textrm{kpc}$. We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work.

We report the discovery of a young (${\rm age} \sim 130~{\rm Myr}$), low-mass ($M \sim 1200~{\rm M}_\odot$), metal-poor ($[{\rm Fe}/{\rm H}] \sim -1.1$) stellar association at a heliocentric distance $D \approx 29~{\rm kpc}$, placing it far into the Milky Way halo. At its present Galactocentric position $(R, z) \sim (23, 15)~{\rm kpc}$, the association is (on the sky) near the leading arm of the gas stream emanating from the Magellanic cloud system, but is located $\approx 60^\circ$ from the Large Magellanic Cloud (LMC) center on the other side of the Milky Way disk. If we assume that the cluster is co-located with HI gas in the stream, we directly measure the distance to the leading arm of the Magellanic stream. The measured distance is inconsistent with stream predictions from models of the LMC/SMC interaction and infall into the Milky Way that do not account for ram pressure and gas interaction with Milky Way disk. The estimated age of the cluster is consistent with the time of last passage of the leading arm gas through the Galactic midplane, and we therefore speculate that this star-formation event was triggered by its the last disk midplane passage. Most details of this idea remain a puzzle: the Magellanic stream has low column density, the Milky Way disk at this large radius has low gas density, and the relative velocity of the leading arm gas and Milky Way gas is large. However it formed, the discovery of a young stellar cluster in the Milky Way halo presents an interesting opportunity for study. This cluster was discovered with Gaia astrometry and photometry alone, but folow-up DECam photometry was crucial for measuring its properties.

Recent measurements suggest that the Large Magellanic Cloud (LMC) may weigh as much as 25\% of the Milky Way. In this work we explore how such a large satellite affects mass estimates of the Milky Way based on equilibrium modelling of the stellar halo or other tracers. In particular, we show that if the LMC is ignored, the Milky Way mass is overestimated by as much as 50\%. This bias is due to the bulk motion in the outskirts of the Galaxy's halo and can be, at least in part, accounted for with a simple modification to the equilibrium modelling. Finally, we show that the LMC has a substantial effect on the orbit Leo I which acts to increase its present day speed relative to the Milky Way. We estimate that accounting for a $1.5\times10^{11} M_\odot$ LMC would lower the inferred Milky Way mass to $\sim10^{12} M_\odot$.

Data from the Gaia satellite show that the solar neighbourhood of the Milky Way's stellar halo is imprinted with substructure from several accretion events. Evidence of these events is found in "the Shards", stars clustering with high significance in both action space and metallicity. Stars in the Shards share a common origin, likely as ancient satellite galaxies of the Milky Way, so will be embedded in dark matter (DM) counterparts. These "Dark Shards" contain two substantial streams (S1 and S2), as well as several retrograde, prograde and lower energy substructures. The retrograde stream S1 has a very high Earth-frame speed of $\sim 550$ km s$^{-1}$ while S2 moves on a prograde, but highly polar orbit and enhances the peak of the speed distribution at around $300$ km s$^{-1}$. The presence of the Dark Shards locally leads to modifications of many to the fundamental properties of experimental DM signals. The~S2 stream in particular gives rise to an array of effects in searches for axions and in the time dependence of nuclear recoils: shifting the peak day, inducing nonsinusoidal distortions, and increasing the importance of the gravitational focusing of DM by the Sun. Dark Shards additionally bring new features for directional signals, while also enhancing the DM flux towards Cygnus.

We use the extensive $Gaia$ Data Release 2 set of Long Period Variables to select a sample of Oxygen-rich Miras throughout the Milky Way disk and bulge for study. Exploiting the relation between Mira pulsation period and stellar age/chemistry, we slice the stellar density of the Galactic disk and bulge as a function of period. We find the morphology of both components evolves as a function of stellar age/chemistry with the stellar disk being stubby at old ages, becoming progressively thinner and more radially extended at younger stellar ages, consistent with the picture of inside-out and upside-down formation of the Milky Way's disk. We see evidence of a perturbed disk, with large-scale stellar over-densities visible both in and away from the stellar plane. We find the bulge is well modelled by a triaxial boxy distribution with an axis ratio of $\sim [1:0.4:0.3]$. The oldest of the Miras ($\sim$ 9-10 Gyr) show little bar-like morphology, whilst the younger stars appear inclined at a viewing angle of $\sim 21^{\circ}$ to the Sun-Galactic Centre line. This suggests that bar formation and buckling took place 8-9 Gyr ago, with the older Miras being hot enough to avoid being trapped by the growing bar. We find the youngest Miras to exhibit a strong peanut morphology, bearing the characteristic X-shape of an inclined bar structure.

We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the star in the Galactic frame is $1755\pm50$ km/s. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A* with a velocity of $\sim 1800$ km/s and travelled for $4.8$ Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion $V_{y,\odot}= 246.1\pm 5.3$ km/s or position $R_0=8.12\pm 0.23$ kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disk of young stars at the Galactic centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time.

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.

We measure the total stellar halo luminosity using red giant branch (RGB) stars selected from Gaia data release 2. Using slices in magnitude, colour and location on the sky, we decompose RGB stars belonging to the disc and halo by fitting 2-dimensional Gaussians to the Galactic proper motion distributions. The number counts of RGB stars are converted to total stellar halo luminosity using a suite of isochrones weighted by age and metallicity, and by applying a volume correction based on the stellar halo density profile. Our method is tested and calibrated using Galaxia and N-body models. We find a total luminosity (out to 100 kpc) of L_halo = 7.9 +/- 2.0 x 10^8 L_Sun excluding Sgr, and L_halo = 9.4 +/- 2.4 x 10^8 L_Sun including Sgr. These values are appropriate for our adopted stellar halo density profile and metallicity distribution, but additional systematics related to these assumptions are quantified and discussed. Assuming a stellar mass-to-light ratio appropriate for a Kroupa initial mass function (M*/L = 1.5), we estimate a stellar halo mass of M*_halo = 1.4 +/- 0.4 x 10^9 M_Sun. This mass is larger than previous estimates in the literature, but is in good agreement with the emerging picture that the (inner) stellar halo is dominated by one massive dwarf progenitor. Finally, we argue that the combination of a ~10^9 M_Sun mass and an average metallicity of <[Fe/H]> ~ -1.5 for the Galactic halo points to an ancient (~10 Gyr) merger event.

Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation.

Ultrafaint dwarf galaxies ($M_\star\lesssim10^{5-6}\,{\rm M}_\odot$) are relics of an early phase of galaxy formation. They contain some of the oldest and most metal-poor stars in the Universe which likely formed before the epoch of hydrogen reionisation. These galaxies are so faint that they can only be detected as satellites of the Milky Way. They are so small that they are just barely resolved in current cosmological hydrodynamics simulations. Here we combine very high resolution cosmological $N$-body simulations with a semi-analytic model of galaxy formation to study the demographics and spatial distribution of ultrafaint satellites in Milky Way-mass haloes. We show that the abundance of these galaxies is correlated with the assembly history of the host halo: at fixed mass, haloes assembled earlier contain, on average, more ultrafaint satellites today than haloes assembled later. We identify simulated galactic haloes that experience an ancient Gaia-Enceladus-Sausage-like and a recent LMC-like accretion event and find that the former occurs in 33% of the sample and the latter in 9%. Only 3% experience both events and these are especially rich in ultrafaint satellites, most acquired during the ancient accretion event. Our models predict that the radial distribution of satellites is more centrally concentrated in early-forming haloes. Accounting for the depletion of satellites by tidal interactions with the central disc, we find a very good match to the observed radial distribution of satellites in the Milky Way over the entire radial range. This agreement is mainly due to the ability of our model to track 'orphan' galaxies after their subhaloes fall below the resolution limit of the simulation.

We present two catalogues of active galactic nucleus (AGN) candidates selected from the latest data of two all-sky surveys -- Data Release 2 (DR2) of the \emph{Gaia} mission and the unWISE catalogue of the \emph{Wide-field Infrared Survey Explorer} (\emph{WISE}). We train a random forest classifier to predict the probability of each source in the \emph{Gaia}-unWISE joint sample being an AGN, $P_{\rm RF}$, based on \emph{Gaia} astrometric and photometric measurements and unWISE photometry. The two catalogues, which we designate C75 and R85, are constructed by applying different $P_{\rm RF}$ threshold cuts to achieve an overall completeness of 75\% ($\approx$90\% at \emph{Gaia} $G\leq20$ mag) and reliability of 85\% respectively. The C75 (R85) catalogue contains 2,734,464 (2,182,193) AGN candidates across the effective 36,000 deg$^2$ sky, of which $\approx$0.91 (0.52) million are new discoveries. Photometric redshifts of the AGN candidates are derived by a random forest regressor using \emph{Gaia} and \emph{WISE} magnitudes and colours. The estimated overall photometric redshift accuracy is 0.11. Cross-matching the AGN candidates with a sample of known bright cluster galaxies, we identify a high-probability strongly-lensed AGN candidate system, SDSS\,J1326$+$4806, with a large image separation of 21\farcs06. All the AGN candidates in our catalogues will have $\sim$5-year long light curves from \emph{Gaia} by the end of the mission, and thus will be a great resource for AGN variability studies. Our AGN catalogues will also be helpful in AGN target selections for future spectroscopic surveys, especially ones in the southern hemisphere. The C75 catalogue can be downloaded at https://www.ast.cam.ac.uk/~ypshu/AGN_Catalogues.html.

The Gaia Radial Velocity Spectrometer (RVS) provides a sample of 7,224,631 stars with full six-dimensional phase space information. Bayesian distances of these stars are available from the catalogue of Sch\"onrich et al. (2019). We exploit this to map out the behaviour of the velocity ellipsoid within 5 kpc of the Sun. We find that the tilt of the disc-dominated RVS sample is accurately described by the relation $\alpha = (0.952 \pm 0.007)\arctan (|z|/R)$, where ($R,z$) are cylindrical polar coordinates. This corresponds to velocity ellipsoids close to spherical alignment (for which the normalising constant would be unity) and pointing towards the Galactic centre. Flattening of the tilt of the velocity ellipsoids is enhanced close to the plane and Galactic centre, whilst at high elevations far from the Galactic center the population is consistent with exact spherical alignment. Using the LAMOST catalogue cross-matched with Gaia DR2, we construct thin disc and halo samples of reasonable purity based on metallicity. We find that the tilt of thin disc stars straddles $\alpha = (0.909-1.038)\arctan (|z|/R)$, and of halo stars straddles $\alpha = (0.927-1.063)\arctan (|z|/R)$. We caution against the use of reciprocal parallax for distances in studies of the tilt, as this can lead to serious artefacts.

We use the distribution of accreted stars in SDSS-Gaia DR2 to demonstrate that a non-trivial fraction of the dark matter halo within Galactocentric radii of 7.5-10 kpc and $|z| > 2.5$ kpc is in substructure, and thus may not be in equilibrium. Using a mixture likelihood analysis, we separate the contributions of an old, isotropic stellar halo and a younger anisotropic population. The latter dominates and is uniform within the region studied. It can be explained as the tidal debris of a disrupted massive satellite on a highly radial orbit, and is consistent with mounting evidence from recent studies. Simulations that track the tidal debris from such mergers find that the dark matter traces the kinematics of its stellar counterpart. If so, our results indicate that a component of the nearby dark matter halo that is sourced by luminous satellites is in kinematic substructure referred to as debris flow. These results challenge the Standard Halo Model, which is discrepant with the distribution recovered from the stellar data, and have important ramifications for the interpretation of direct detection experiments.

Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc.

We cross-match high-precision astrometric data from Gaia DR2 with accurate multi-band photometry from the Dark Energy Survey (DES) DR1 to confidently measure proper motions for nine stellar streams in the DES footprint: Aliqa Uma, ATLAS, Chenab, Elqui, Indus, Jhelum, Phoenix, Tucana III, and Turranburra. We determine low-confidence proper motion measurements for four additional stellar streams: Ravi, Wambelong, Willka Yaku, and Turbio. We find evidence for a misalignment between stream tracks and the systemic proper motion of streams that may suggest a systematic gravitational influence from the Large Magellanic Cloud. These proper motions, when combined with radial velocity measurements, will allow for detailed orbit modeling which can be used to constrain properties of the LMC and its affect on nearby streams, as well as global properties of the Milky Way's gravitational potential.

We study the orbits of ultra-faint dwarf galaxies in the combined presence of the Milky Way and LMC and we find 6 dwarfs which were likely accreted with the LMC (Car 2, Car 3, Hor 1, Hyi 1, Phe 2, Ret 2), in addition to the SMC, representing strong evidence of dwarf galaxy group infall. This procedure depends on the gravitational pull of the LMC, thus allowing us to place a lower bound on the Cloud's mass of $M_{\rm LMC} > 1.24\times10^{11} M_\odot$. This mass estimate is validated by applying the technique to a cosmological zoom-in simulation of a Milky Way-like galaxy with an LMC analogue where we find that while this lower bound may be overestimated, it will improve in the future with smaller observational errors. We apply this technique to dwarf galaxies lacking radial velocities and find that Eri 3 has a broad range of radial velocities for which it has a significant chance ($> 0.4$) of having being bound to the Cloud. We study the non-Magellanic classical satellites and find that Fornax has an appreciable probability of being an LMC satellite if the LMC is sufficiently massive. In addition, we explore how the orbits of the Milky Way satellites change in the presence of the LMC and find a significant change for several objects. Finally, we find that the LMC satellites are slightly smaller than the Milky Way satellites at a fixed luminosity, possibly due to the different tidal environments they have experienced.

We report the discovery of a Milky Way satellite in the constellation of Antlia. The Antlia 2 dwarf galaxy is located behind the Galactic disc at a latitude of $b\sim 11^{\circ}$ and spans $1.26$ degrees, which corresponds to $\sim2.9$ kpc at its distance of 130 kpc. While similar in spatial extent to the Large Magellanic Cloud, Antlia 2 is orders of magnitude fainter at $M_V=-9$ mag, making it by far the lowest surface brightness system known (at $\sim31.9$ mag/arcsec$^2$), $\sim100$ times more diffuse than the so-called ultra diffuse galaxies. The satellite was identified using a combination of astrometry, photometry and variability data from \textit{Gaia} Data Release 2, and its nature confirmed with deep archival DECam imaging, which revealed a conspicuous BHB signal. We have also obtained follow-up spectroscopy using AAOmega on the AAT, identifying 159 member stars, and we used them to measure the dwarf's systemic velocity, $290.9\pm0.5$km/s, its velocity dispersion, $5.7\pm1.1$ km/s, and mean metallicity, [Fe/H]$=-1.4$. From these properties we conclude that Antlia 2 inhabits one of the least dense Dark Matter (DM) halos probed to date. Dynamical modelling and tidal-disruption simulations suggest that a combination of a cored DM profile and strong tidal stripping may explain the observed properties of this satellite. The origin of this core may be consistent with aggressive feedback, or may even require alternatives to cold dark matter (such as ultra-light bosons).

The Gaia Sausage is the major accretion event that built the stellar halo of the Milky Way galaxy. Here, we provide dynamical and chemical evidence for a second substantial accretion episode, distinct from the Gaia Sausage. The Sequoia Event provided the bulk of the high energy retrograde stars in the stellar halo, as well as the recently discovered globular cluster FSR 1758. There are up to 6 further globular clusters, including $\omega$~Centauri, as well as many of the retrograde substructures in Myeong et al. (2018), associated with the progenitor dwarf galaxy, named the Sequoia. The stellar mass in the Sequoia galaxy is $\sim 5 \times 10^{7} M_\odot$, whilst the total mass is $\sim 10^{10} M_\odot$, as judged from abundance matching or from the total sum of the globular cluster mass. Although clearly less massive than the Sausage, the Sequoia has a distinct chemo-dynamical signature. The strongly retrograde Sequoia stars have a typical eccentricity of $\sim0.6$, whereas the Sausage stars have no clear net rotation and move on predominantly radial orbits. On average, the Sequoia stars have lower metallicity by $\sim 0.3$ dex and higher abundance ratios as compared to the Sausage. We conjecture that the Sausage and the Sequoia galaxies may have been associated and accreted at a comparable epoch.

We extend predictive microlensing event searches using the Vista Variables in the Via Lactea survey and the second Gaia data release. We identify two events with maxima in 2019 that require urgent follow-up. First, we predict that the nearby M2 dwarf L 338-152 will align with a background source with a closest approach of $35^{+35}_{-23}$ mas on 2019 November $16^{+28}_{-27}$ d. This will cause a peak astrometric shift and photometric amplification of the background source of $2.7^{+3.5}_{-1.5}$ mas and $5.6^{+143.2}_{-5.2}$ mmag respectively. This event should be astrometrically detectable by both the Hubble Space Telescope (HST) and the Spectro-Polarimetric High-contrast Exoplanet Research instrument on the Very Large Telescope. Secondly, we predict the likely K dwarf NLTT 45128 will lens a background source with a closest approach of $105.3^{+12.2}_{-11.7}$ mas on 2019 September $26^{+15}_{-15}$ d. This will produce a peak astrometric shift of $0.329^{+0.065}_{-0.059}$ mas. NLTT 45128 is only 3.6 magnitudes brighter than the background source which makes it an excellent candidate for follow-up with HST. Characterisation of these signals will allow direct gravitational masses to be inferred for both L 338-152 and NLTT 45128 with an estimated precision of $\sim9$ and $\sim13$ per cent respectively.

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion.

Using RR Lyrae stars in the Gaia Data Release 2 and Pan-STARRS1 we study the properties of the Pisces Overdensity, a diffuse sub-structure in the outer halo of the Milky Way. We show that along the line of sight, Pisces appears as a broad and long plume of stars stretching from 40 to 110 kpc with a steep distance gradient. On the sky Pisces's elongated shape is aligned with the Magellanic Stream. Using follow-up VLT FORS2 spectroscopy, we have measured the velocity distribution of the Pisces candidate member stars and have shown it to be as broad as that of the Galactic halo but offset to negative velocities. Using a suite of numerical simulations, we demonstrate that the structure has many properties in common with the predicted behaviour of the Magellanic wake, i.e. the Galactic halo overdensity induced by the in-fall of the Magellanic Clouds.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The goal of this survey is to study the formation and evolution of the Milky Way halo to deduce its assembly history and the 3D distribution of mass in the Milky Way. The combination of multi-band photometry, Gaia proper motion and parallax data, and radial velocities and the metallicity and elemental abundances obtained from low-resolution spectra of halo giants with 4MOST, will yield an unprecedented characterisation of the Milky Way halo and its interface with the thick disc. The survey will produce a volume- and magnitude-limited complete sample of giant stars in the halo. It will cover at least 10,000 square degrees of high Galactic latitude, and measure line-of-sight velocities with a precision of 1-2 km/s as well as metallicities to within 0.2 dex.

We model the fastest moving (v_tot > 300 km/s) local (D < 3 kpc) halo stars using cosmological simulations and 6-dimensional Gaia data. Our approach is to use our knowledge of the assembly history and phase-space distribution of halo stars to constrain the form of the high velocity tail of the stellar halo. Using simple analytical models and cosmological simulations, we find that the shape of the high velocity tail is strongly dependent on the velocity anisotropy and number density profile of the halo stars --- highly eccentric orbits and/or shallow density profiles have more extended high velocity tails. The halo stars in the solar vicinity are known to have a strongly radial velocity anisotropy, and it has recently been shown the origin of these highly eccentric orbits is the early accretion of a massive (M_star ~ 10^9 M_Sun) dwarf satellite. We use this knowledge to construct a prior on the shape of the high velocity tail. Moreover, we use the simulations to define an appropriate outer boundary of 2r_200, beyond which stars can escape. After applying our methodology to the Gaia data, we find a local (r_0=8.3 kpc) escape speed of v_esc(r_0) = 528(+24,-25) km/s. We use our measurement of the escape velocity to estimate the total Milky Way mass, and dark halo concentration: M_200,tot = 1.00(+0.31,-0.24) x 10^12 M_Sun, c_200 = 10.9(+4.4,-3.3). Our estimated mass agrees with recent results in the literature that seem to be converging on a Milky Way mass of M_200,tot ~ 10^12 M_Sun.

Since first noticed by Shapley in 1939, a faint object coincident with the Fornax dwarf spheroidal has long been discussed as a possible sixth globular cluster system. However, debate has continued over whether this overdensity is a statistical artifact or a blended galaxy group. In this Letter we demonstrate, using deep DECam imaging data, that this object is well resolved into stars and is a bona fide star cluster. The stellar overdensity of this cluster is statistically significant at the level of ~ 6 - 6.7 sigma in several different photometric catalogs including Gaia. Therefore, it is highly unlikely to be caused by random fluctuation. We show that Fornax 6 is a star cluster with a peculiarly low surface brightness and irregular shape, which may indicate a strong tidal influence from its host galaxy. The Hess diagram of Fornax 6 is largely consistent with that of Fornax field stars, but it appears to be slightly bluer. However, it is still likely more metal-rich than most of the globular clusters in the system. Faint clusters like Fornax 6 that orbit and potentially get disrupted in the centers of dwarf galaxies can prove crucial for constraining the dark matter distribution in Milky Way satellites.

We explore the origin of the observed decline in [O/Fe] (and [Mg/Fe]) with Galactocentric distance for high-metallicity stars ([Fe/H] > -1.1), based on a sample of halo stars selected within the Apache Point Observatory Galactic Evolution Experiment (APOGEE) fourteenth data release (DR14). We also analyse the characteristics of the [$\alpha$/Fe] distributions in the inner-halo regions inferred from two zoom-in Milky Way mass-sized galaxies that are taken as case studies. One of them qualitatively reproduces the observed trend to have higher fraction of $\alpha$-rich star for decreasing galactocentric distance; the other exhibits the opposite trend. We find that stars with [Fe/H] > -1.1 located in the range [15 - 30] kpc are consistent with formation in two starbursts, with maxima separated by about ~ 1 Gyr. We explore the contributions of stellar populations with different origin to the [$\alpha$/Fe] gradients detected in stars with [Fe/H] > -1.1. Our analysis reveals that the simulated halo that best matches the observed chemical trends is characterised by an accretion history involving low to intermediate-mass satellite galaxies with a short and intense burst of star formation, and contributions from a more massive satellite with dynamical masses about ~ 10$^{10}$M$_{\odot}$, distributing low-[$\alpha$/Fe] stars at intermediate radius.

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber.

We present the results of a systematic Milky Way satellite search performed across an array of publicly available wide-area photometric surveys. Our aim is to complement previous searches by widening the parameter space covered. Specifically, we focus on objects smaller than $1'$ and include old, young, metal poor and metal rich stellar population masks. As a result we find 9 new likely genuine stellar systems in data from GAIA, DES, and Pan-STARRS, which were picked from the candidate list because of conspicuous counterparts in the cut-out images. The presented systems are all very compact ($r_h<1'$) and faint ($M_V\gtrsim-3$), and are associated either with the Galactic disk, or the Magellanic Clouds. While most of the stellar systems look like Open Clusters, their exact classification is, as of today, unclear. With these discoveries, we extend the parameter space occupied by star clusters to sizes and luminosities previously unexplored and demonstrate that rather than two distinct classes of Globular and Open clusters, there appears to be a continuity of objects, unmarked by a clear decision boundary.

We carried out a spectroscopic follow-up program of the four new stellar stream candidates detected by Belokurov & Koposov (2016) in the outskirts of the Large Magellanic Cloud (LMC) using FORS2 (VLT). The medium-resolution spectra were used to measure the line-of-sight velocities, estimate stellar metallicities and to classify stars into Blue Horizontal Branch (BHB) and Blue Straggler (BS) stars. Using the 4-D phase-space information, we attribute approximately one half of our sample to the Magellanic Clouds, while the rest is part of the Galactic foreground. Only two of the four stream candidates are confirmed kinematically. While it is impossible to estimate the exact levels of MW contamination, the phase-space distribution of the entire sample of our Magellanic stars matches the expected velocity gradient for the LMC halo and extends as far as 33 deg (angular separation) or 29 kpc from the LMC center. Our detections reinforce the idea that the halo of the LMC seems to be larger than previously expected, and its debris can be spread in the sky out to very large separations from the LMC center. Finally, we provide some kinematic evidence that many of the stars analysed here have likely come from the Small Magellanic Cloud.

The recently discovered S1 stream passes through the Solar neighbourhood on a low inclination, counter-rotating orbit. The progenitor of S1 is a dwarf galaxy with a total mass comparable to the present-day Fornax dwarf spheroidal, so the stream is expected to have a significant DM component. We compute the effects of the S1 stream on WIMP and axion detectors as a function of the density of its unmeasured dark component. In WIMP detectors the S1 stream supplies more high energy nuclear recoils so will marginally improve DM detection prospects. We find that even if S1 comprises less than 10% of the local density, multi-ton xenon WIMP detectors can distinguish the S1 stream from the bulk halo in the relatively narrow mass range between 5 and 25 GeV. In directional WIMP detectors such as CYGNUS, S1 increases DM detection prospects more substantially since it enhances the anisotropy of the WIMP signal. Finally, we show that axion haloscopes possess by far the greatest potential sensitivity to the S1 stream. Once the axion mass has been discovered, the distinctive velocity distribution of S1 can easily be extracted from the axion power spectrum.

Iorio et al. (2018) mapped out the Milky Way halo using a sample of RR Lyrae stars drawn from a cross-match of Gaia with 2MASS. We investigate the significant residual in their model which we constrain to lie at Galactocentric radii $12<R<27\;\mathrm{kpc}$ and extend over $2600\;\mathrm{deg}^2$ of the sky. A counterpart of this structure exists in both the Catalina Real Time Survey and the sample of RR Lyrae variables identified in Pan-STARRS by Hernitschek et al. (2016), demonstrating that this structure is not caused by the spatial inhomogeneity of Gaia. The structure is likely the Virgo Stellar Stream and/or Virgo Over-Density. We show the structure is aligned with the Magellanic Stream and suggest that it is either debris from a disrupted dwarf galaxy that was a member of the Vast Polar Structure or that it is SMC debris from a tidal interaction of the SMC and LMC $3\;\mathrm{Gyr}$ ago. If the latter then the sub-structure in Virgo may have a Magellanic origin.

Using the astrometry from the ESA's Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significant components with high radial velocity anisotropy, beta > 0.6. However, only in one third of the hosts do the high-beta stars contribute significantly to the accreted stellar halo overall, similar to what is observed in the Milky Way. For this particular subset we reveal the origin of the dominant stellar halo component with high metallicity, [Fe/H]~-1, and high orbital anisotropy, beta>0.8, by tracing their stars back to the epoch of accretion. It appears that, typically, these stars come from a single dwarf galaxy with a stellar mass of order of 10^9-10^10 Msol that merged around 6-10 Gyr ago, causing a sharp increase in the halo mass. Our study therefore establishes a firm link between the excess of radially anisotropic stellar debris in the Milky Way halo and an ancient head-on collision between the young Milky Way and a massive dwarf galaxy

We use the $Gaia$ DR2 RR Lyrae sample to gain an uninterrupted view of the Galactic stellar halo. We dissect the available volume in slices parallel to the Milky Way's disc to show that within $\sim30$ kpc from the Galactic centre the halo is triaxial, with the longest axis misaligned by $\sim70^{\circ}$ with respect to the Galactic $x$-axis. This anatomical procedure exposes two large diffuse over-densities aligned with the semi-major axis of the halo: the Hercules-Aquila Cloud and the Virgo Over-density. We reveal the kinematics of the entire inner halo by mapping out the amplitudes and directions of the RR Lyrae proper motions. These are then compared to simple models with different anisotropies to demonstrate that the inner halo is dominated by stars on highly eccentric orbits. We interpret the shape of the density and the kinematics of the $Gaia$ DR2 RR Lyrae as evidence in favour of a scenario in which the bulk of the halo was deposited in a single massive merger event.

We report the discovery of a thin stellar stream - which we name the Jet stream - crossing the constellations of Hydra and Pyxis. The discovery was made in data from the SLAMS survey, which comprises deep $g$ and $r$ imaging for a $650$ square degree region above the Galactic disc performed by the CTIO Blanco + DECam. SLAMS photometric catalogues will be made publicly available. The stream is approximately 0.18 degrees wide and 10 degrees long, though it is truncated by the survey footprint. Its colour-magnitude diagram is consistent with an old, metal-poor stellar population at a heliocentric distance of approximately 29 kpc. We corroborate this measurement by identifying a spatially coincident overdensity of likely blue horizontal branch stars at the same distance. There is no obvious candidate for a surviving stream progenitor.

We used Gaia Data Release 2 to search for upcoming photometric microlensing events, identifying two candidates with high amplification. In the case of candidate 1, a spectrum of the lens (l1) confirms it is a usdM3 subdwarf with mass $\approx 0.11 M_\odot$, while the event reaches maximum amplification of $20^{+20}_{-10}$ mmag on November 3rd 2019 ($\pm$1d). For candidate 2, the lens (l2) is a metal-poor M dwarf with mass $\approx 0.38 M_\odot$ derived from spectral energy distribution (SED) fitting, and maximum amplification of $10^{+40}_{-10}$ mmag occurs on June 3rd 2019 ($\pm$4d). This permits a new algorithm for mass inference on the microlens. Given the predicted time, the photometric lightcurve of these events can be densely sampled by ground-based telescopes. The lightcurve is a function of the unknown lens mass, together with 8 other parameters for all of which Gaia provides measurements and uncertainties. Leveraging this prior information on the source and lens provided by Gaia's astrometric solution, and assuming that a ground-based campaign can provide 50 measurements at mmag precision, we show for example that the mass of l1 can be recovered to within 20 per cent (68 per cent confidence limit).

We use astrometry and broad-band photometry from Data Release 2 of the ESA's Gaia mission to map out low surface-brightness features in the stellar density distribution around the Large and Small Magellanic Clouds. The LMC appears to have grown two thin and long stellar streams in its Northern and Southern regions, highly reminiscent of spiral arms. We use computer simulations of the Magellanic Clouds' in-fall to demonstrate that these arms were likely pulled out of the LMC's disc due to the combined influence of the SMC's most recent fly-by and the tidal field of the Milky Way.

We measure the orbital properties of halo stars using 7-dimensional information provided by Gaia and the Sloan Digital Sky Survey. A metal-rich population of stars, present in both local main sequence stars and more distant blue horizontal branch stars, have very radial orbits (eccentricity ~ 0.9) and apocenters that coincide with the stellar halo "break radius" at galactocentric distance r ~ 20 kpc. Previous work has shown that the stellar halo density falls off much more rapidly beyond this break radius. We argue that the correspondence between the apocenters of high metallicity, high eccentricity stars and the broken density profile is caused by the build-up of stars at the apocenter of a common dwarf progenitor. Although the radially biased stars are likely present down to metallicities of [Fe/H] ~ -2 the increasing dominance at higher metallicities suggests a massive dwarf progenitor, which is at least as massive as the Fornax and Sagittarius dwarf galaxies, and is likely the dominant progenitor of the inner stellar halo.

We use a sample of ~350 RR Lyrae stars with radial velocities and Gaia DR2 proper motions to study orbital properties of the Hercules-Aquila Cloud (HAC) and Virgo Over-density (VOD). We demonstrate that both structures are dominated by stars on highly eccentric orbits, with peri-centres around ~1 kpc and apo-centres between 15 and 25 kpc from the Galactic centre. Given that the stars in the HAC and the VOD occupy very similar regions in the space spanned by integrals of motion, we conclude that these diffuse debris clouds are part of the same accretion event. More precisely, these inner halo sub-structures likely represent two complementary not-fully-mixed portions of an ancient massive merger, also known as the "sausage" event.

We use the SDSS-Gaia catalogue to search for substructure in the stellar halo. The sample comprises 62\,133 halo stars with full phase space coordinates and extends out to heliocentric distances of $\sim 10$ kpc. As actions are conserved under slow changes of the potential, they permit identification of groups of stars with a common accretion history. We devise a method to identify halo substructures based on their clustering in action space, using metallicity as a secondary check. This is validated against smooth models and numerical constructed stellar halos from the Aquarius simulations. We identify 21 substructures in the SDSS-Gaia catalogue, including 7 high significance, high energy and retrograde ones. We investigate whether the retrograde substructures may be material stripped off the atypical globular cluster $\omega$~Centauri. Using a simple model of the accretion of the progenitor of the $\omega$~Centauri, we tentatively argue for the possible association of up to 5 of our new substructures (labelled Rg1, Rg3, Rg4, Rg6 and Rg7) with this event. This sets a minimum mass of $5 \times 10^8 M_\odot$ for the progenitor, so as to bring $\omega$~Centauri to its current location in action -- energy space. Our proposal can be tested by high resolution spectroscopy of the candidates to look for the unusual abundance patterns possessed by $\omega$~Centauri stars.

The Gaia Sausage is an elongated structure in velocity space discovered by Belokurov et al. (2018) using the kinematics of metal-rich halo stars. It was created by a massive dwarf galaxy ($\sim 5 \times 10^{10} M_\odot$) on a strongly radial orbit that merged with the Milky Way at a redshift $z\lesssim 3$. We search forthe associated Sausage Globular Clusters by analysing the structure of 91 Milky Way globular clusters (GCs) in action space using the Gaia Data Release 2 catalogue, complemented with Hubble Space Telescope proper motions. There is a characteristic energy $E_{\rm crit}$ which separates the in situ objects, such as the bulge/disc clusters, from the accreted objects, such as the young halo clusters. There are 15 old halo GCs that have $E > E_{\rm crit}$. Eight of the high energy, old halo GCs are strongly clumped in azimuthal and vertical action, yet strung out like beads on a chain at extreme radial action. They are very radially anisotropic ($\beta \sim 0.95$) and move on orbits that are all highly eccentric ($e \gtrsim 0.80$). They also form a track in the age-metallicity plane distinct from the bulk of the Milky Way GCs and compatible with a dwarf spheroidal origin. These properties are consistent with GCs associated with the merger event that gave rise to the Gaia Sausage.

We investigate the utility of the 3-dimensional Two-Point Correlation Function (3D 2PCF) in quantifying substructure in the stellar halo of the Milky Way, particularly as a means of constraining the accretion history of our Galaxy. We use RR Lyrae variable stars from two different surveys as tracers of the structure in the Galactic stellar halo. We compare our measurements of the 3D 2PCF in these datasets to a suite of simulations of the formation of the stellar halo from Bullock and Johnston (2005). While there is some room for interpretation, we find that the amounts of structure to be broadly consistent with the simulations, while appearing smoother than average within the inner halo and at small scales. This suggests a preferred accretion history scenario in which the Milky Way's stellar halo acquired most of its mass about ~8 Gigayears ago, and has been largely quiescent since. Finally, we discuss the prospects of statistical tools such as the 2PCF in the Gaia era of galactic archaeology.

We assemble the largest sample of oxygen rich Miras to date and highlight their importance for age-dating the components of the Galaxy. Using data from the Catalina Rapid Transient Survey and the All Sky Automated Survey for Supernovae, we extract a clean sample of $\sim 2,400$ O-Miras, stretching from the Galactic Bulge to the distant halo. Given that the period of O-Miras correlates with age, this offers a new way of determining age gradients throughout the Galaxy. We use our sample to show (i) disk O-Miras have periods increasing on moving outwards from ~ 3 to 15 kpc, so the outer disk O-Miras are younger than the inner disk, (ii) the transition from younger disk to halo O-Miras occurs at r ~ 15 kpc and is marked by a plummeting in period, (iii) there exists a population of young O-Miras likely kicked from the disk to heights of order of |Z| ~ 10 kpc, (iv) great circle counts of old Miras show strong evidence for distant debris agglomeration associated with the Magellanic Clouds, (v) seven stars in our samples are located at distances between 200 and 500 kpc surpassing all previously established records, and, finally, (vi) O-Miras may be present in the Fornax, Sculptor, Sextans and Leo II Galactic dwarf spheroidals, as well as the distant globular cluster Pal 4. We spotlight the importance of O-Mira in the Era of Gaia as universal chronometers of the Galactic populations.

We present results of the first dynamical stream fits to the recently discovered Tucana III stream. These fits assume a fixed Milky Way potential and give proper motion predictions, which can be tested with the upcoming Gaia Data Release 2. These fits reveal that Tucana III is on an eccentric orbit around the Milky Way and, more interestingly, that Tucana III passed within 15 kpc of the Large Magellanic Cloud (LMC) approximately 75 Myr ago. Given this close passage, we fit the Tucana III stream in the combined presence of the Milky Way and the LMC. We find that the predicted proper motions depend on the assumed mass of the LMC and that the LMC can induce a substantial proper motion perpendicular to the stream track. A detection of this misalignment will directly probe the extent of the LMC's influence on our Galaxy, and has implications for nearly all methods which attempt to constraint the Milky Way potential. Such a measurement will be possible with the upcoming Gaia DR2, allowing for a measurement of the LMC's mass.

We used the Tycho-Gaia Astrometric Solution catalogue, part of the Gaia Data Release 1, to search for candidate astrometric microlensing events expected to occur within the remaining lifetime of the Gaia satellite. Our search yielded one promising candidate. We predict that the nearby DQ type white dwarf LAWD 37 (WD 1142-645) will lens a background star and will reach closest approach on November 11th 2019 ($\pm$ 4 days) with impact parameter $380\pm10$ mas. This will produce an apparent maximum deviation of the source position of $2.8\pm0.1$ mas. In the most propitious circumstance, Gaia will be able to determine the mass of LAWD 37 to $\sim3\%$. This mass determination will provide an independent check on atmospheric models of white dwarfs with helium rich atmospheres, as well as tests of white dwarf mass radius relationships and evolutionary theory.

We use a new panoramic imaging survey, conducted with the Dark Energy Camera, to map the stellar fringes of the Large and Small Magellanic Clouds to extremely low surface brightness V $\gtrsim$ 32 mag arcsec$^{-2}$. Our results starkly illustrate the closely interacting nature of the LMC-SMC pair. We show that the outer LMC disk is strongly distorted, exhibiting an irregular shape, evidence for warping, and significant truncation on the side facing the SMC. Large diffuse stellar substructures are present both to the north and south of the LMC, and in the inter-Cloud region. At least one of these features appears co-spatial with the bridge of RR Lyrae stars that connects the Clouds. The SMC is highly disturbed -- we confirm the presence of tidal tails, as well as a large line-of-sight depth on the side closest to the LMC. Young, intermediate-age, and ancient stellar populations in the SMC exhibit strikingly different spatial distributions. In particular, those with ages $\sim$ 1.5-4 Gyr exhibit a spheroidal distribution with a centroid offset from that of the oldest stars by several degrees towards the LMC. We speculate that the gravitational influence of the LMC may already have been perturbing the gaseous component of the SMC several Gyr ago. With careful modeling, the variety of substructures and tidal distortions evident in the Magellanic periphery should tightly constrain the interaction history of the Clouds.

We report the discovery of a nearby dwarf galaxy in the constellation of Hydrus, between the Large and the Small Magellanic Clouds. Hydrus 1 is a mildy elliptical ultra-faint system with luminosity $M_V\sim$ -4.7 and size $\sim$ 50 pc, located 28 kpc from the Sun and 24 kpc from the LMC. From spectroscopy of $\sim$ 30 member stars, we measure a velocity dispersion of 2.7 km/s and find tentative evidence for a radial velocity gradient consistent with 3 km/s rotation. Hydrus 1's velocity dispersion indicates that the system is dark matter dominated, but its dynamical mass-to-light ratio M/L $\sim$ 66 is significantly smaller than typical for ultra-faint dwarfs at similar luminosity. The kinematics and spatial position of Hydrus~1 make it a very plausible member of the family of satellites brought into the Milky Way by the Magellanic Clouds. While Hydrus 1's proximity and well-measured kinematics make it a promising target for dark matter annihilation searches, we find no evidence for significant gamma-ray emission from Hydrus 1. The new dwarf is a metal-poor galaxy with a mean metallicity [Fe/H]=-2.5 and [Fe/H] spread of 0.4 dex, similar to other systems of similar luminosity. Alpha-abundances of Hyi 1 members indicate that star-formation was extended, lasting between 0.1 and 1 Gyr, with self-enrichment dominated by SN Ia. The dwarf also hosts a highly carbon-enhanced extremely metal-poor star with [Fe/H] $\sim$ -3.2 and [C/Fe] $\sim$ +3.0.

We analyse the structure of the local stellar halo of the Milky Way using $\sim$ 60000 stars with full phase space coordinates extracted from the SDSS--{\it Gaia} catalogue. We display stars in action space as a function of metallicity in a realistic axisymmetric potential for the Milky Way Galaxy. The metal-rich population is more distended towards high radial action $J_R$ as compared to azimuthal or vertical action, $J_\phi$ or $J_z$. It has a mild prograde rotation $(\langle v_\phi \rangle \approx 25$ km s$^{-1}$), is radially anisotropic and highly flattened with axis ratio $q \approx 0.6 - 0.7$. The metal-poor population is more evenly distributed in all three actions. It has larger prograde rotation $(\langle v_\phi \rangle \approx 50$ km s$^{-1}$), a mild radial anisotropy and a roundish morphology ($q\approx 0.9$). We identify two further components of the halo in action space. There is a high energy, retrograde component that is only present in the metal-rich stars. This is suggestive of an origin in a retrograde encounter, possibly the one that created the stripped dwarf galaxy nucleus, $\omega$Centauri. Also visible as a distinct entity in action space is a resonant component, which is flattened and prograde. It extends over a range of metallicities down to [Fe/H] $\approx -3$. It has a net outward radial velocity $\langle v_R \rangle \approx 12$ km s$^{-1}$ within the Solar circle at $|z| <3.5$ kpc. The existence of resonant stars at such extremely low metallicities has not been seen before.

We designed a follow-up program to find the spectroscopic properties of the Hercules-Aquila Cloud (HAC) and test scenarios for its formation. We measured the radial velocities (RVs) of 45 RR Lyrae in the southern portion of the HAC using the facilities at the MDM observatory, producing the first large sample of velocities in the HAC. We found a double-peaked distribution in RVs, skewed slightly to negative velocities. We compared both the morphology of HAC projected onto the plane of the sky and the distribution of velocities in this structure outlined by RR Lyrae and other tracer populations at different distances to N-body simulations. We found that the behaviour is characteristic of an old, well-mixed accretion event with small apo-galactic radius. We cannot yet rule out other formation mechanisms for the HAC. However, if our interpretation is correct, HAC represents just a small portion of a much larger debris structure spread throughout the inner Galaxy whose distinct kinematic structure should be apparent in RV studies along many lines of sight.

We present Magellan/IMACS, Anglo-Australian Telescope/AAOmega+2dF, and Very Large Telescope/GIRAFFE+FLAMES spectroscopy of the CarinaII (Car II) & Carina III (Car III) dwarf galaxy candidates, recently discovered in the Magellanic Satellites Survey (MagLiteS). We identify 18 member stars in Car II, including 2 binaries with variable radial velocities and 2 RR Lyrae stars. The other 14 members have a mean heliocentric velocity $v_{\rm hel} = 477.2 \pm 1.2$ km/s and a velocity dispersion of $\sigma_v = 3.4^{+1.2}_{-0.8}$ km/s. Assuming Car II is in dynamical equilibrium, we derive a total mass within the half-light radius of $1.0^{+0.8}_{-0.4} \times 10^{6} M_\odot$, indicating a mass-to-light ratio of $369^{+309}_{-161} M_\odot/L_\odot$. From equivalent width measurements of the calcium triplet lines of 9 RGB stars, we derive a mean metallicity of [Fe/H] = $-2.44 \pm 0.09$ with dispersion $\sigma_{\rm [Fe/H]} = 0.22 ^{+0.10}_{-0.07}$. Considering both the kinematic and chemical properties, we conclude that Car II is a dark-matter-dominated dwarf galaxy. For Car III, we identify 4 member stars, from which we calculate a systemic velocity of $v_{\rm hel} = 284.6^{+3.4}_{-3.1}$ km/s. The brightest RGB member of Car III has a metallicity of [Fe/H] $= -1.97 \pm 0.12$. Due to the small size of the Car III spectroscopic sample, we cannot conclusively determine its nature. Although these two systems have the smallest known physical separation ($\Delta d\sim10~kpc$) among Local Group satellites, the large difference in their systemic velocities, $\sim200$ km/s, indicates that they are unlikely to be a bound pair. One or both systems are likely associated with the Large Magellanic Cloud (LMC), and may remain LMC satellites today. No statistically significant excess of $\gamma$-rays emission is found at the locations of Car II and Car III in eight years of Fermi-LAT data.

Using a large sample of Main Sequence stars with 7-D measurements supplied by Gaia and SDSS, we study the kinematic properties of the local (within ~10 kpc from the Sun) stellar halo. We demonstrate that the halo's velocity ellipsoid evolves strongly with metallicity. At the low [Fe/H] end, the orbital anisotropy (the amount of motion in the radial direction compared to the tangential one) is mildly radial with 0.2<beta<0.4. However, for stars with [Fe/H]>-1.7 we measure extreme values of beta~0.9. Across the metallicity range considered, i.e. -3<[Fe/H]<-1, the stellar halo's spin is minimal, at the level of 20<v_theta (km/s) <30. Using a suite of cosmological zoom-in simulations of halo formation, we deduce that the observed acute anisotropy is inconsistent with the continuous accretion of dwarf satellites. Instead, we argue, the stellar debris in the inner halo were deposited in a major accretion event by a satellite with Mvir>10^10 Msun around the epoch of the Galactic disc formation, i.e. between 8 and 11 Gyr ago. The radical halo anisotropy is the result of the dramatic radialisation of the massive progenitor's orbit, amplified by the action of the growing disc.

We present the first spatially-resolved rest-frame UV study of the gravitationally lensed galaxy, the 'Cosmic Horseshoe' (J1148+1930) at z=2.38. Our gravitational lens model shows that the system is made up of four star-forming regions, each ~4-8 kpc^2 in size, from which we extract four spatially exclusive regional spectra. We study the interstellar and wind absorption lines, along with CIII] doublet emission lines, in each region to investigate any variation in emission/absorption line properties. The mapped CIII] emission shows distinct kinematical structure, with velocity offsets of ~+/-50 km/s between regions suggestive of a merging system, and a variation in equivalent width that indicates a change in ionisation parameter and/or metallicity between the regions. Absorption line velocities reveal a range of outflow strengths, with gas outflowing between -200<v(km/s)<-50 relative to the systemic velocity of that region. Interestingly, the strongest gas outflow appears to emanate from the most diffuse star-forming region. The star-formation rates remain relatively constant (~8-16 M_sol/yr), mostly due to large uncertainties in reddening estimates. As such, the outflows appear to be 'global' rather than 'locally' sourced. We measure electron densities with a range of log(Ne)=3.92-4.36 cm^-3, and point out that such high densities may be common when measured using the CIII] doublet due to its large critical density. Overall, our observations demonstrate that while it is possible to trace variations in large scale gas kinematics, detecting inhomogeneities in physical gas properties and their effects on the outflowing gas may be more difficult. This study provides important lessons for the spatially-resolved rest-frame UV studies expected with future observatories, such as JWST.

Using deep photometric data from CFHT/Megacam, we study the morphology and density of the GD-1 stream, one of the longest and coldest stellar streams in the Milky Way. Our deep data recovers the lower main sequence of the stream with unprecedented quality, clearly separating it from Milky Way foreground and background stars. An analysis of the distance to different parts of the stream shows that GD-1 lies at a heliocentric distance between 8 and 10 kpc, with only a shallow gradient across 45 deg on the sky. Matched filter maps of the stream density show clear density variations, such as deviations from a single orbital track and tentative evidence for stream fanning. We also detect a clear under-density in the middle of the stream track at $\varphi_{1}$=-45 deg surrounded by overdense stream segments on either side. This location is a promising candidate for the elusive missing progenitor of the GD-1 stream. We conclude that the GD-1 stream has clearly been disturbed by interactions with the Milky Way disk or other sub-halos.

We report the discovery of two ultra-faint satellites in the vicinity of the Large Magellanic Cloud (LMC) in data from the Magellanic Satellites Survey (MagLiteS). Situated 18$^{\circ}$ ($\sim 20$ kpc) from the LMC and separated from each other by only $18^\prime$, Carina~II and III form an intriguing pair. By simultaneously modeling the spatial and the color-magnitude stellar distributions, we find that both Carina~II and Carina~III are likely dwarf galaxies, although this is less clear for Carina~III. There are in fact several obvious differences between the two satellites. While both are well described by an old and metal poor population, Carina~II is located at $\sim 36$ kpc from the Sun, with $M_V\sim-4.5$ and $r_h\sim 90$ pc, and it is further confirmed by the discovery of 3 RR Lyrae at the right distance. In contrast, Carina~III is much more elongated, measured to be fainter ($M_V\sim-2.4$), significantly more compact ($r_h\sim30$ pc), and closer to the Sun, at $\sim 28$ kpc, placing it only 8 kpc away from Car~II. Together with several other systems detected by the Dark Energy Camera, Carina~II and III form a strongly anisotropic cloud of satellites in the vicinity of the Magellanic Clouds.

We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of $g > 23.5$ with a relative photometric calibration uncertainty of $< 1 \%$. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.

We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess the significance of halo substructures based on their clustering in velocity space. The two most substantial structures are multiple wraps of a stream which has undergone considerable phase mixing (S1, with 94 members) and a kinematically cold stream (S2, with 61 members). The member stars of S1 have a median position of ($X,Y,Z$) = ($8.12, -0.22, 2.75$) kpc and a median metallicity of [Fe/H] $= -1.78$. The stars of S2 have median coordinates ($X,Y,Z$) = ($8.66, 0.30, 0.77$) kpc and a median metallicity of [Fe/H] $= -1.91$. They lie in velocity space close to some of the stars in the stream reported by Helmi et al. (1999). By modelling, we estimate that both structures had progenitors with virial masses $\approx 10^{10} M_\odot$ and infall times $\gtrsim 9$ Gyr ago. Using abundance matching, these correspond to stellar masses between $10^6$ and $10^7 M_\odot$. These are somewhat larger than the masses inferred through the mass-metallicity relation by factors of 5 to 15. Additionally, we identify two further substructures (S3 and S4 with 55 and 40 members) and two clusters or moving groups (C1 and C2 with 24 and 12) members. In all 6 cases, clustering in kinematics is found to correspond to clustering in both configuration space and metallicity, adding credence to the reliability of our detections.

We exploit the exquisite, deep Hyper Suprime-Cam (HSC) imaging data to probe the Galactic halo out to 200 kpc. Using the ~100 square degree, multi-band photometry of the first HSC Wide survey data release, we identify blue horizontal branch (BHB) stars beyond 50 kpc in the halo. The presence of the Sagittarius (Sgr) stream in the HSC fields produces a notable excess of stars at the apocentre of the leading arm (~50-60 kpc). For fields excluding Sgr, the BHB counts are consistent with a continuation of a -4 power-law from the inner halo. However, we find that the majority of the non-Sgr BHB stars beyond 50 kpc reside in one 27 square degree HSC field called "VVDS". Curiously, this field is located close to the Magellanic plane, and we hypothesize that the excess of stars between 50 and 200 kpc could be associated with distant Magellanic debris. Indeed, without the VVDS, there are very few BHBs in the remaining portions of the Galaxy probed by the HSC. Accordingly, this scarcity of tracers is consistent with a significant decline in stellar density beyond 50 kpc, with a power-law of -4 or steeper.

We explore the connection between galaxies and dark matter halos in the Milky Way (MW) and quantify the implications on properties of the dark matter particle and the phenomenology of low-mass galaxy formation. This is done through a probabilistic comparison of the luminosity function of MW dwarf satellite galaxies to models based on two suites of zoom-in simulations. One suite is dark-matter-only while the other includes a disk component, therefore we can quantify the effect of the MW's baryonic disk on our results. We apply numerous Stellar-Mass-Halo-Mass (SMHM) relations allowing for multiple complexities: scatter, a characteristic break scale, and subhalos which host no galaxy. In contrast to previous works we push the model/data comparison to the faintest dwarfs by modeling observational incompleteness, allowing us to draw three new conclusions. Firstly, we constrain the SMHM relation for $10^2<M_*/M_\odot<10^8$ galaxies, allowing us to bound the peak halo mass of the faintest MW satellite to $M_\mathrm{vir}>2.4\times10^8M_\odot$ ($1\sigma$). Secondly, by translating to a Warm Dark Matter (WDM) cosmology, we bound the thermal relic mass $m_\mathrm{WDM}>2.9$ keV at 95\% confidence, on a par with recent constraints from the Lyman-$\alpha$ forest. Lastly, we find that the observed number of ultra-faint MW dwarfs is in tension with the theoretical prediction that reionisation prevents galaxy formation in almost all $10^8M_\odot$ halos. This can be tested with the next generation of deep imaging surveys. To this end, we predict the likely number of detectable satellite galaxies in the Subaru/HSC survey and the LSST. Confronting these predictions with future observations will be amongst our strongest tests of WDM and the effect reionisation on low-mass systems.

We show that tagging RR Lyrae stars according to their location in the period-amplitude diagram can be used to shed light on the genesis of the Galactic stellar halo. The mixture of RR Lyrae of ab type, separated into classes along the lines suggested by Oosterhoff, displays a strong and coherent evolution with Galactocentric radius. The change in the RR Lyrae composition appears to coincide with the break in the halo's radial density profile at ~25 kpc. Using simple models of the stellar halo, we establish that at least three different types of accretion events are necessary to explain the observed RRab behavior. Given that there exists a correlation between the RRab class fraction and the total stellar content of a dwarf satellite, we hypothesize that the field halo RRab composition is controlled by the mass of the progenitor contributing the bulk of the stellar debris at the given radius. This idea is tested against a suite of cosmological zoom-in simulations of Milky Way-like stellar halo formation. Finally, we study some of the most prominent stellar streams in the Milky Way halo and demonstrate that their RRab class fractions follow the trends established previously.

We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.

We present a comprehensive and precise description of the Sagittarius (Sgr) stellar stream's 3D geometry as traced by its old stellar population. This analysis draws on the sample of ${\sim}44,000$ RR Lyrae (RRab) stars from the Pan-STARRS1 (PS1) 3$\pi$ survey (Hernitschek et al. 2016,Sesar et al. 2017b), which is ${\sim}80\%$ complete and ${\sim}90\%$ pure within 80~kpc, and extends to ${\gtrsim} 120$~kpc with a distance precision of ${\sim} 3\%$. A projection of RR Lyrae stars within $|\tilde{B}|_{\odot}<9^\circ$ of the Sgr stream's orbital plane reveals the morphology of both the leading and the trailing arms at very high contrast, across much of the sky. In particular, the map traces the stream near-contiguously through the distant apocenters. We fit a simple model for the mean distance and line-of-sight depth of the Sgr stream as a function of the orbital plane angle $\tilde{\Lambda}_{\odot}$, along with a power-law background-model for the field stars. This modeling results in estimates of the mean stream distance precise to ${\sim}1\%$ and it resolves the stream's line-of-sight depth. These improved geometric constraints can serve as new constraints for dynamical stream models.

We study the kinematic properties of the Eastern Banded Structure (EBS) and Hydra I overdensity using exquisite proper motions derived from the Sloan Digital Sky Survey (SDSS) and Gaia source catalog. Main sequence turn-off stars in the vicinity of the EBS are identified from SDSS photometry; we use the proper motions and, where applicable, spectroscopic measurements of these stars to probe the kinematics of this apparent stream. We find that the EBS and Hydra I share common kinematic and chemical properties with the nearby Monoceros Ring. In particular, the proper motions of the EBS, like Monoceros, are indicative of prograde rotation (V_phi ~ 180-220 km/s), which is similar to the Galactic thick disc. The kinematic structure of stars in the vicinity of the EBS suggest that it is not a distinct stellar stream, but rather marks the "edge" of the Monoceros Ring. The EBS and Hydra I are the latest substructures to be linked with Monoceros, leaving the Galactic anti-centre a mess of interlinked overdensities which likely share a unified, Galactic disc origin.

We use a new contiguous imaging survey conducted using the Dark Energy Camera to investigate the distribution and properties of young stellar populations in the Magellanic inter-Cloud region. These young stars are strongly spatially clustered, forming a narrow chain of low-mass associations that trace the densest HI gas in the Magellanic Bridge and extend, in projection, from the SMC to the outer disk of the LMC. The associations in our survey footprint have ages $\lesssim 30$ Myr, masses in the range $\sim 100-1200\,{\rm M}_\odot$, and very diffuse structures with half-light radii of up to $\sim 100$ pc. The two most populous are strongly elliptical, and aligned to $\approx 10^{{\rm o}}$ with the axis joining the centres of the LMC and SMC. These observations strongly suggest that the young inter-Cloud populations formed in situ, likely due to the compression of gas stripped during the most recent close LMC-SMC encounter. The associations lie at distances intermediate between the two Clouds, and we find no evidence for a substantial distance gradient across the imaged area. Finally, we identify a vast shell of young stars surrounding a central association, that is spatially coincident with a low column density bubble in the HI distribution. The properties of this structure are consistent with a scenario where stellar winds and supernova explosions from massive stars in the central cluster swept up the ambient gas into a shell, triggering a new burst of star formation. This is a prime location for studying stellar feedback in a relatively isolated environment.

We study the structure of the inner Milky Way using the latest data release of the Vista Variables in Via Lactea (VVV) survey. The VVV is a deep near-infrared, multi-colour photometric survey with a coverage of 300 square degrees towards the Bulge/Bar. We use Red Clump (RC) stars to produce a high-resolution dust map of the VVV's field of view. From de-reddened colour-magnitude diagrams we select Red Giant Branch stars to investigate their 3D density distribution within the central 4 kpc. We demonstrate that our best-fit parametric model of the Bulge density provides a good description of the VVV data, with a median percentage residual of 5$\%$ over the fitted region. The strongest of the otherwise low-level residuals are overdensities associated with a low-latitude structure as well as the so-called X-shape previously identified using the split RC. These additional components contribute only $\sim5\%$ and $\sim7\%$ respectively to the Bulge mass budget. The best-fit Bulge is `boxy' with an axis ratio of [1:0.44:0.31] and is rotated with respect to the Sun-Galactic Centre line by at least $20^{\circ}$. We provide an estimate of the total, full sky, mass of the Bulge of $M_\mathrm{Bulge}^{\mathrm{Chabrier}} = 2.36 \times 10^{10} M_{\odot}$ for a Chabrier initial mass function. We show there exists a strong degeneracy between the viewing angle and the dispersion of the Red Clump absolute magnitude distribution. The value of the latter is strongly dependent on the assumptions made about the intrinsic luminosity function of the Bulge.

We report the discovery of a stellar over-density 8$^{\circ}$ north of the center of the Small Magellanic Cloud (Small Magellanic Cloud Northern Over-Density; SMCNOD) using data from the first two years of the Dark Energy Survey (DES) and the first year of the MAGellanic SatelLITEs Survey (MagLiteS). The SMCNOD is indistinguishable in age, metallicity and distance from the nearby SMC stars, being primarly composed of intermediate-age stars (6 Gyr, Z=0.001), with a small fraction of young stars (1 Gyr, Z=0.01). The SMCNOD has an elongated shape with an ellipticity of 0.6 and a size of $\sim$ 6x2 deg. It has an absolute magnitude of $M_V \cong$ -7.7, $r_h = 2.1$ kpc, and $\mu_V(r<r_h)$ = 31.2 mag arcsec$^{-2}$. We estimate a stellar mass of $\sim 10^5$ $M_{\odot}$, following a Kroupa mass function. The SMCNOD was probably removed from the SMC disk by tidal stripping, since it is located near the head of the Magellanic Stream, and the literature indicates likely recent LMC-SMC encounters. This scenario is supported by the lack of significant HI gas. Other potential scenarios for the SMCNOD origin are a transient over-density within the SMC tidal radius or a primordial SMC satellite in advanced stage of disruption.

We present the first detailed study of the behaviour of the stellar proper motion across the entire Galactic Anti-centre area visible in the Sloan Digital Sky Survey data. We use recalibrated SDSS astrometry in combination with positions from {\it Gaia} DR1 to provide tangential motion measurements with a systematic uncertainty $<$5 kms$^{-1}$ for the Main Sequence stars at the distance of the Monoceros Ring. We demonstrate that Monoceros members rotate around the Galaxy with azimuthal speeds of ~230 kms$^{-1}$, only slightly lower than that of the Sun. Additionally, both vertical and azimuthal components of their motion are shown to vary considerably but gradually as a function of Galactic longitude and latitude. The stellar over-density in the Anti-centre region can be split into two components, the narrow, stream-like ACS and the smooth Ring. According to our analysis, these two structures show very similar but clearly distinct kinematic trends, which can be summarised as follows: the amplitude of the velocity variation in $v_{\phi}$ and $v_z$ in the ACS is higher compared to the Ring, whose velocity gradients appear to be flatter. Currently, no model available can explain the entirety of the data in this area of the sky. However, the new accurate kinematic map introduced here should provide strong constraints on the genesis of the Monoceros Ring and the associated sub-structure.

Only in the Milky Way is it possible to conduct an experiment which uses stellar streams to detect low-mass dark matter subhaloes. In smooth and static host potentials, tidal tails of disrupting satellites appear highly symmetric. However, perturbations from dark subhaloes, as well as from GMCs and the Milky Way bar, can induce density fluctuations that destroy this symmetry. Motivated by the recent release of unprecedentedly deep and wide imaging data around the Pal~5 stellar stream, we develop a new probabilistic, adaptive and non-parametric technique which allows us to bring the cluster's tidal tails into clear focus. Strikingly, we uncover a stream whose density exhibits visible changes on a variety of angular scales. We detect significant bumps and dips, both narrow and broad: two peaks on either side of the progenitor, each only a fraction of a degree across, and two gaps, $\sim2^{\circ}$ and $\sim9^{\circ}$ wide, the latter accompanied by a gargantuan lump of debris. This largest density feature results in a pronounced inter-tail asymmetry which cannot be made consistent with an unperturbed stream according to a suite of simulations we have produced. We conjecture that the sharp peaks around Pal 5 are epicyclic overdensities, while the two dips are consistent with impacts by subhaloes. Assuming an age of 3.4 Gyr for Pal 5, these two gaps would correspond to the characteristic size of gaps created by subhaloes in the mass range of $10^6-10^7 M_\odot$ and $10^7-10^8 M_\odot$ respectively. In addition to dark substructure, we find that the bar of the Milky Way can plausibly produce the asymmetric density seen in Pal 5 and that GMCs could cause the smaller gap.

We combine Gaia data release 1 astrometry with Sloan Digital Sky Survey (SDSS) images taken some ~10-15 years earlier, to measure proper motions of stars in the halo of our Galaxy. The SDSS-Gaia proper motions have typical statistical errors of 2 mas/yr down to r ~ 20 mag, and are robust to variations with magnitude and colour. Armed with this exquisite set of halo proper motions, we identify RR Lyrae, blue horizontal branch (BHB), and K giant stars in the halo, and measure their net rotation with respect to the Galactic disc. We find evidence for a gently rotating prograde signal (< $V_\phi$ > ~ 5-25 km/s) in the halo stars, which shows little variation with Galactocentric radius out to 50 kpc. The average rotation signal for the three populations is < $V_\phi$ > = 14 +/- 2 +/- 10 (syst.) km/s. There is also tentative evidence for a kinematic correlation with metallicity, whereby the metal richer BHB and K giant stars have slightly stronger prograde rotation than the metal poorer stars. Using the Auriga simulation suite we find that the old (T >10 Gyr) stars in the simulated halos exhibit mild prograde rotation, with little dependence on radius or metallicity, in general agreement with the observations. The weak halo rotation suggests that the Milky Way has a minor in situ halo component, and has undergone a relatively quiet accretion history.

We present the results of the very first search for faint Milky Way satellites in the Gaia data. Using stellar positions only, we are able to re-discover objects detected in much deeper data as recently as the last couple of years. While we do not identify new prominent ultra-faint dwarf galaxies, we report the discovery of two new star clusters, Gaia 1 and Gaia 2. Gaia 1 is particularly curious, as it is a massive (2.2$\times$10$^4$ M$_\odot$), large ($\sim$9 pc) and nearby (4.6 kpc) cluster, situated 10' away from the brightest star on the sky, Sirius! Even though this satellite is detected at significance in excess of 10, it was missed by previous sky surveys. We conclude that Gaia possesses powerful and unique capabilities for satellite detection thanks to its unrivalled angular resolution and highly efficient object classification.

We report the discovery of tidal tails around the Galactic globular cluster NGC 7492, based on the Data Release 1 of the Pan-STARRS 1 survey. The tails were detected with a version of the matched filter technique applied to the $(g-r,r)$ and $(g-i,i)$ color-magnitude diagrams. Tidal tails emerging from the cluster extend at least $\sim$3.5 degrees in the North-East to South-East direction, equivalent to $\sim1.5$ kpc in projected length.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

We identify a halo substructure in the Tycho Gaia Astrometric Solution (TGAS) dataset, cross-matched with the RAVE-on data release. After quality cuts, the stars with large radial action ($J_R > 800$ kms$^{-1}$ kpc) are extracted. A subset of these stars is clustered in longitude and velocity and can be selected with further cuts. The 14 stars are centered on $(X,Y,Z) \approx (9.0,-1.0,-0.6)$ kpc and form a coherently moving structure in the halo with median $(v_R,v_\phi,v_z) = (167.33,0.86,-94.85)$ kms$^{-1}$. They are all metal-poor giants with median [Fe/H] $=-0.83$. To guard against the effects of distance errors, we compute spectrophotometric distances for the 8 out of 14 stars where this is possible. We find that 6 of the stars are still comoving. These 6 stars also have a much tighter [Fe/H] distribution $\sim -0.7$ with one exception ([Fe/H] = -2.12). We conclude that the existence of the comoving cluster is stable against changes in distance estimation and conjecture that this is the dissolving remnant of a long-ago accreted globular cluster.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We exploit the first data release from the Gaia mission to identify candidate Mira variables in the outskirts of the Magellanic Clouds. The repeated observations of sources during the initial phase of the Gaia mission is used to identify stars that show signs of variability. This variability information, combined with infrared photometry from 2MASS and WISE, allows us to select a clean sample of giants in the periphery of the LMC. We find evidence for Miras surrounding the LMC out to ~20 deg in all directions, apart from the North-West quadrant. Our sample does not generally follow the gas distribution of the Magellanic system; Miras are notably absent in the gaseous bridge between the LMC and SMC, but they are likely related to the stellar RR Lyrae bridge reported by Belokurov et al. (2016). The stellar stream discovered by Mackey et al. (2016) to the North of the LMC is almost perfectly delineated by our Mira variables, and likely extends further East toward the Galactic plane. The presence of an intermediate-age population in this stream advocates an LMC disc origin. We also find a significant excess of Miras to the East of the LMC; these more diffusely distributed stars are likely stripped SMC stars due to interactions with the LMC. Miras are also identified in regions of the sky away from the Clouds; we locate stars likely associated with known massive substructures, and also find potential associations with stripped SMC debris above the Galactic plane.

We present the results of the spectroscopic and photometric follow-up of two field galaxies that were selected as possible stellar counterparts of local high velocity clouds. Our analysis shows that the two systems are distant (D>20 Mpc) dwarf irregular galaxies unrelated to the local HI clouds. However, the newly derived distance and structural parameters reveal that the two galaxies have luminosities and effective radii very similar to the recently identified Ultra Diffuse Galaxies (UDGs). At odds with classical UDGs, they are remarkably isolated, having no known giant galaxy within ~2.0 Mpc. Moreover, one of them has a very high gas content compared to galaxies of similar stellar mass, with a HI to stellar mass ratio M_HI/M_* ~90, typical of almost-dark dwarfs. Expanding on this finding, we show that extended dwarf irregulars overlap the distribution of UDGs in the M_V vs. log(r_e) plane and that the sequence including dwarf spheroidals, dwarf irregulars and UDGs appears as continuously populated in this plane.

We measure the variation of the escape speed of the Galaxy across a range of $\sim$ 40 kpc in Galactocentric radius. The local escape speed is found to be $521^{+46}_{-30}\,\mathrm{km\,s^{-1}}$, in good agreement with other studies. We find that this has already fallen to $379^{+34}_{-28}\,\mathrm{km\,s^{-1}}$ at a radius of 50 kpc. Through measuring the escape speed and its variation, we obtain constraints on the Galactic potential as a whole. In particular, the gradient in the escape speed with radius suggests that the total mass contained within 50 kpc is $29^{+7}_{-5}\times10^{10}\mathrm{M}_\odot$, implying a relatively light dark halo for the Milky Way. Our method represents a novel way of estimating the mass of the Galaxy, and has very different systematics to more commonly used models of tracers, which are more sensitive to the central parts of the halo velocity distributions. Using our inference on the escape speed, we then investigate the orbits of high-speed Milky Way dwarf galaxies. For each dwarf we consider, we predict small pericenter radii and large orbital eccentricities. This naturally explains the large observed ellipticities of two of the dwarfs, which are likely to have been heavily disrupted as they passed through pericenter.

We take advantage of the deep and wide coverage of the VST ATLAS survey to study the line-of-sight structure of the Sagittarius stellar stream in the Southern hemisphere, only ~40{\deg} away from the progenitor. We use photometrically selected Sub-Giant Branch (SGB) stars to reveal a complex debris morphology of the trailing arm and detect at least two clear peaks in the SGB distance modulus distribution. The separation between the two line-of-sight components is at least 5 kpc at the edge of the VST ATLAS footprint, but appears to change along the stream, which allows us to conclude that these detections correspond to two physically independent stellar structures, rather than a mix of co-distant stellar populations within a single stream. Our discovery of a fork in the Sgr trailing arm is verified using Blue Horizontal Branch stars and our distance measurements are calibrated using RR Lyrae stars from the Catalina Real-Time Transient Survey. Comparing with numerical simulations of the Sgr dwarf disruption, the more distant of the two components in the fork matches perfectly with the track of the trailing debris. However, no obvious counterpart exists in the simulation for the closer line-of-sight component.

We present results from MMT/Hectochelle spectroscopy of red giant candidate stars along the line of sight to the recently-discovered Galactic satellite Crater 2. Modelling the joint distribution of stellar positions, velocities and metallicities as a mixture of Crater 2 and Galactic foreground populations, we identify 62 members of Crater 2, for which we resolve line-of-sight velocity dispersion 2.7 +/- 0.3 km/s about mean velocity of 87.5 +/- 0.4 km/s. We also resolve a metallicity dispersion 0.22 about a mean of [Fe/H]=-1.98 +/- 0.1 that is 0.28 +/- 0.14 poorer than is estimated from photometry. Despite Crater 2's relatively large size (projected halflight radius R(h)=1 kpc) and intermediate luminosity (M_V =-8), its velocity dispersion is the coldest that has been resolved for any dwarf galaxy. These properties make Crater 2 the most extreme low-density outlier in dynamical as well as structural scaling relations among the Milky Way's dwarf spheroidals. Even so, under assumptions of dynamical equilibrium and negligible contamination by unresolved binary stars, the observed velocity distribution implies a gravitationally dominant dark matter halo, with dynamical mass 4.4 x 10^6 Msun and mass-to-light ratio M/L=53 enclosed within a radius of 1 kpc, where the equivalent circular velocity is 4.3 km/s.

Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.

The search for chemically unevolved galaxies remains prevalent in the nearby Universe, mostly because these systems provide excellent proxies for exploring in detail the physics of high-z systems. The most promising candidates are extremely metal-poor galaxies (XMPs), i.e., galaxies with <1/10 solar metallicity. However, due to the bright emission line based search criteria traditionally used to find XMPs, we may not be sampling the full XMP population. In 2014 we reoriented this search using only morphological properties and uncovered a population of ~150 `blue diffuse dwarf (BDD) galaxies', and published a sub-sample of 12 BDD spectra. Here we present optical spectroscopic observations of a larger sample of 51 BDDs, along with their SDSS photometric properties. With our improved statistics, we use direct-method abundances to confirm that BDDs are chemically unevolved (7.43<12+log(O/H)<8.01), with ~20% of our sample classified as being XMP galaxies, and find they are actively forming stars at rates of 1-33x10^-2 M_sol/yr in HII regions randomly embedded in a blue, low-surface brightness continuum. Stellar masses are calculated from population synthesis models and estimated to be in the range log(M_star/M_sol) ~5-9. Unlike other low-metallicity star-forming galaxies, BDDs are in agreement with the mass-metallicity relation at low masses, suggesting they are not accreting large amounts of pristine gas relative to their stellar mass. BDD galaxies appear to be a population of actively star-forming dwarf irregular (dIrr) galaxies who fall within the class of low-surface brightness dIrr galaxies. Their ongoing star-formation and irregular morphology make them excellent analogues for galaxies in the early Universe.

We present the discovery of stellar tidal tails around the Large and the Small Magellanic Clouds in the Gaia DR1 data. In between the Clouds, their tidal arms are stretched towards each other to form an almost continuous stellar bridge. Our analysis relies on the exquisite quality of the Gaia's photometric catalogue to build detailed star-count maps of the Clouds. We demonstrate that the Gaia DR1 data can be used to detect variable stars across the whole sky, and in particular, RR Lyrae stars in and around the LMC and the SMC. Additionally, we use a combination of Gaia and Gale to follow the distribution of Young Main Sequence stars in the Magellanic System. Viewed by Gaia, the Clouds show unmistakable signs of interaction. Around the LMC, clumps of RR Lyrae are observable as far as ~20 degrees, in agreement with the most recent map of Mira-like stars reported in Deason et al (2016). The SMC's outer stellar density contours show a characteristic S-shape, symptomatic of the on-set of tidal stripping. Beyond several degrees from the center of the dwarf, the Gaia RR Lyrae stars trace the Cloud's trailing arm, extending towards the LMC. This stellar tidal tail mapped with RR Lyrae is not aligned with the gaseous Magellanic Bridge, and is shifted by some ~5 degrees from the Young Main Sequence bridge. We use the offset between the bridges to put constraints on the density of the hot gaseous corona of the Milky Way.

We announce the discovery of the Aquarius~2 dwarf galaxy, a new distant satellite of the Milky Way, detected on the fringes of the VST ATLAS and the SDSS surveys. The object was originally identified as an overdensity of Red Giant Branch stars, but chosen for subsequent follow-up based on the presence of a strong Blue Horizontal Branch, which was also used to measure its distance of $\sim 110$ kpc. Using deeper imaging from the IMACS camera on the 6.5m Baade and spectroscopy with DEIMOS on Keck, we measured the satellite's half-light radius $5.1\pm 0.8$ arcmin, or $\sim 160$ pc at this distance, and its stellar velocity dispersion of $5.4^{+3.4}_{-0.9}$ km s$^{-1}$. With $\mu=30.2$ mag arcsec$^{-2}$ and $M_V=-4.36$, the new satellite lies close to two important detection limits: one in surface brightness; and one in luminosity at a given distance, thereby making Aquarius~2 one of the hardest dwarfs to find.

Ample observational capabilities exist today to detect the small density perturbations that low-mass dark matter subhaloes impart on stellar streams from disrupting Galactic satellites. In anticipation of these observations, we investigate the expected number and size of gaps by combining an analytic prescription for gap evolution on circular orbits with the flux of subhaloes near the stream. We explore the distribution of gap sizes and depths for a typical cold stream around the Milky Way and find that for a given stream age and gap depth, each subhalo mass produces a characteristic gap size. For a stream with an age of a few Gyr, orbiting at a distance of 10-20 kpc from the Galactic center, even modest subhaloes with a mass of $10^6-10^7 M_\odot$ produce gaps with sizes that are on the order of several degrees. We consider the number and distribution of gap sizes created by subhaloes with masses $10^5-10^9 M_\odot$, accounting for the expected depletion of subhaloes by the Milky Way disk, and present predictions for six cold streams around the Milky Way. For Pal 5, we forecast 0.7 gaps with a density depletion of at least 25\% and a typical gap size of $8^\circ$. Thus, there appears to be {\it no tension} between the recent non-detection of density depletions in the Pal 5 tidal tails and $\Lambda$CDM expectations. These predictions can be used to guide the scale of future gap searches.

In aspherical potentials orbital planes continuously evolve. The gravitational torques impel the angular momentum vector to precess, that is to slowly stray around the symmetry axis, and nutate, i.e. swing up and down periodically in the perpendicular direction. This familiar orbital pole motion - if detected and measured - can reveal the shape of the underlying gravitational potential, the quantity only crudely gauged in the Galaxy so far. Here we demonstrate that the debris poles of stellar tidal streams show a very similar straying and swinging behavior, and give analytic expressions to link the amplitude and the frequency of the pole evolution to the flattening of the dark matter distribution. While these results are derived for near-circular orbits, we show they are also valid for eccentric orbits. Most importantly, we explain how the differential orbital plane precession leads to the broadening of the stream and show that streams on polar orbits ought to scatter faster. We provide expressions for the stream width evolution as a function of the axisymmetric potential flattening and the angle from the symmetry plane and prove that our models are in good agreement with streams produced in N-body simulations. Interestingly, the same intuition applies to streams whose progenitors are on short or long-axis loops in a triaxial potential. Finally, we present a compilation of the Galactic cold stream data, and discuss how the simple picture developed here, along with stream modelling, can be used to constrain the symmetry axes and flattening of the Milky Way.

We announce the discovery of the Crater 2 dwarf galaxy, identified in imaging data of the VST ATLAS survey. Given its half-light radius of ~1100 pc, Crater 2 is the fourth largest dwarf in the Milky Way, surpassed only by the LMC, SMC and the Sgr dwarf. With a total luminosity of $M_V\approx-8$, this satellite galaxy is also one of the lowest surface brightness dwarfs. Falling under the nominal detection boundary of 30 mag arcsec$^{-2}$, it compares in nebulosity to the recently discovered Tuc 2 and Tuc IV and UMa II. Crater 2 is located ~120 kpc from the Sun and appears to be aligned in 3-D with the enigmatic globular cluster Crater, the pair of ultra-faint dwarfs Leo IV and Leo V and the classical dwarf Leo II. We argue that such arrangement is probably not accidental and, in fact, can be viewed as the evidence for the accretion of the Crater-Leo group.

We use the SDSS/SEGUE spectroscopic sample of stars in the leading and trailing streams of the Sagittarius (Sgr) to demonstrate the existence of two sub-populations with distinct chemistry and kinematics. The metallicity distribution function (MDF) of the trailing stream is decomposed into two Gaussians describing a metal-rich sub-population with means and dispersions (-0.74, 0.18) dex and a metal-poor with (-1.33, 0.27) dex. The metal-rich sub-population has a velocity dispersion ~ 8 km/s, whilst the metal-poor is nearly twice as hot ~13 km/s. For the leading stream, the MDF is again well-described by a superposition of two Gaussians, though somewhat shifted as compared to the trailing stream. The metal-rich has mean and dispersion (-1.00, 0.34) dex, the metal-poor (-1.39, 0.22) dex. The velocity dispersions are inflated by projection effects to give 15 - 30 km/s for the metal-poor, and 6 - 20 km/s for the metal-rich, depending on longitude. We infer that, like many dwarf spheroidals, the Sgr progenitor possessed a more extended, metal-poor stellar component and less extended, metal-rich one. We study the implications of this result for the progenitor mass by simulating the disruption of the Sgr, represented as King light profiles in dark halos of masses between $10^{10}$ and $10^{11} M_\odot$, in a three-component Milky Way whose halo is a live Truncated Flat potential in the first phase of accretion and a triaxial Law \& Majewski model in the second phase. We show that that the dark halo of the Sgr must have been $\gtrsim 6 \times 10^{10} M_\odot$ to reproduce the run of velocity dispersion with longitude for the metal-rich and metal-poor sub-populations in the tails.

We present an update of a spectroscopic follow-up survey at low-resolution of a large number of RR Lyrae halo overdensity candidates found in the southern sky. The substructure candidates were identified in the RR Lyrae catalog of Torrealba et al. (2015) using Catalina Real-time Transient Survey (CRTS) data. Radial velocities and mean metallicities have been estimated for target stars in almost half of the original overdensities to assess their potential membership to coherent halo features.

We establish the connection between the Magellanic Clouds (MCs) and the dwarf galaxy candidates discovered in the Dark Energy Survey (DES) by building a dynamical model of the MC satellite populations, based on an extensive suite of tailor-made numerical simulations. Our model takes into account the response of the Galaxy to the MCs infall, the dynamical friction experienced by the MCs and the disruption of the MC satellites by their hosts. The simulation suite samples over the uncertainties in the MC's proper motions, the masses of the MW and the Clouds themselves and allows for flexibility in the intrinsic volume density distribution of the MC satellites. As a result, we can accurately reproduce the DES satellites' observed positions and kinematics. Assuming that Milky Way (MW) dwarfs follow the distribution of subhaloes in $\Lambda$CDM, we further demonstrate that, of 14 observed satellites, the MW halo contributes fewer than 4 (8) of these with 68% (95%) confidence and that 7 (12) DES dwarfs have probabilities greater than 0.7 (0.5) of belonging to the LMC. Marginalising over the entire suite, we constrain the total number of the Magellanic satellites at ~70, the mass of the LMC around $10^{11}M_\odot$ and show that the Clouds have likely endured only one Galactic pericentric passage so far. Finally, we give predictions for the line-of-sight velocities and the proper motions of the satellites discovered in the vicinity of the LMC.

We report the discovery of a substantial stellar overdensity in the periphery of the Large Magellanic Cloud (LMC), found using public imaging from the first year of the Dark Energy Survey. The structure appears to emanate from the edge of the outer LMC disk at a radius $\approx 13.5$ degrees due north of its centre, and stretches more than $10$ kpc towards the east. It is roughly $1.5$ kpc wide and has an integrated $V$-band luminosity of at least $M_V = -7.4$. The stellar populations in the feature are indistinguishable from those in the outer LMC disk. We attempt to quantify the geometry of the outer disk using simple planar models, and find that only a disk with mild intrinsic ellipticity can simultaneously explain the observed stellar density on the sky and the azimuthal line-of-sight distance profile. We also see possible non-planar behaviour in the outer disk that may reflect a warp and/or flare, as well as deviations that resemble a ring-like structure between $\sim9-12$ degrees from the LMC centre. Based on all these observations, we conclude that our remote, stream-like feature is likely comprised of material that has been stripped from the outskirts of the LMC disk, although we cannot rule out that it represents a transient overdensity in the disk itself. We conduct a simple $N$-body simulation to show that either type of structure could plausibly arise due to the tidal force of the Milky Way; however we also recognize that a recent close interaction between the LMC and the SMC may be the source of the stripping or perturbation. Finally, we observe evidence for extremely diffuse LMC populations extending to radii of $\sim 18.5$ kpc in the disk plane ($\approx 20$ degrees on the sky), corroborating previous spectroscopic detections at comparable distances.

The shape of the Galactic dark halo can, in principle, be inferred through modelling of stellar tidal streams in the Milky Way halo. The brightest and the longest of these, the Sagittarius stream, reaches out to large Galactocentric distances and hence can deliver the tightest constraints on the Galaxy's potential. In this contribution, we revisit the idea that the Sagittarius Stream was formed from a rotating progenitor. We demonstrate that the angle between the disk's angular momentum and the progenitor's orbital angular momentum does not remain constant throughout the disruption. Instead, it undergoes a dramatic evolution caused, in part, by the changes in the progenitor's moment of inertia tensor. We show that, even in a spherical potential, the streams produced as a result of a disky dwarf disruption appear to be "precessing". Yet, this debris plane evolution is illusory as it is solely caused by the swaying and wobbling of the progenitor's disk. Stream plane precession is therefore not an unambiguous indicator of asphericity of the dark halo.

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1, respectively. Our sample for Tuc 2 is sufficient to resolve an internal velocity dispersion of $8.6_{-2.7}^{+4.4}$ km s$^{-1}$ about a mean of $-129.1_{-3.5}^{+3.5}$ km s$^{-1}$ (solar rest frame), and to estimate a mean metallicity of [Fe/H]= $-2.23_{-0.12}^{+0.18}$. These results place Tuc 2 on chemodynamical scaling relations followed by dwarf galaxies, suggesting a dominant dark matter component with dynamical mass $2.7_{-1.3}^{+3.1}\times 10^6$ $\mathrm{M}_{\odot}$ enclosed within the central $\sim 160$ pc, and dynamical mass-to-light ratio $1900_{-900}^{+2200}$ $\mathrm{M}_{\odot}/L_{V,\odot}$. For Gru 1 we estimate a mean velocity of $-140.5_{-1.6}^{+2.4}$ km s$^{-1}$ and a mean metallicity of [Fe/H]=$-1.42_{-0.42}^{+0.55}$, but our sample does not resolve Gru 1's velocity dispersion. The radial coordinates of Tuc 2 and Gru 1 in Galactic phase space suggest that their orbits are among the most energetic within distance $\leq 300$ kpc. Moreover, their proximity to each other in this space arises naturally if both objects are trailing the Large Magellanic Cloud.

Using Blue Horizontal Branch stars identified in the Dark Energy Survey Year 1 data, we report the detection of an extended and lumpy stellar debris distribution around the Magellanic Clouds. At the heliocentric distance of the Clouds, overdensities of BHBs are seen to reach at least to ~30 degrees, and perhaps as far as ~50 degrees from the LMC. In 3D, the stellar halo is traceable to between 25 and 50 kpc from the LMC. We catalogue the most significant of the stellar sub-structures revealed, and announce the discovery of a number of narrow streams and diffuse debris clouds. Two narrow streams appear approximately aligned with the Magellanic Clouds' proper motion. Moreover, one of these overlaps with the gaseous Magellanic Stream on the sky. Curiously, two diffuse BHB agglomerations seem coincident with several of the recently discovered DES satellites. Given the enormous size and the conspicuous lumpiness of the LMC's stellar halo, we speculate that the dwarf could easily have been more massive than previously had been assumed.

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]

Cold or Warm, the Dark Matter substructure spectrum must extend to objects with masses as low as $10^7 M_\odot$, according to the most recent Lyman-$\alpha$ measurements. Around a Milky Way-like galaxy, more than a thousand of these subhaloes will not be able to form stars but are dense enough to survive even deep down in the potential well of their host. There, within the stellar halo, these dark pellets will bombard tidal streams as they travel around the Galaxy, causing small but recognizable damage to the stream density distribution. The detection and characterization of these stream ruptures will allow us to constrain the details of the subhalo-stream interaction. In this work, for the first time, we will demonstrate how the properties of a subhalo, most importantly its mass and size, can be reliably inferred from the gap it produces in a tidal stream. For a range of realistic observational setups, mimicking e.g. SDSS, DES, Gaia and LSST data, we find that it is possible to measure the {\it complete set} of properties (including the phase-space coordinates during the flyby) of dark perturbers with $M>10^7 M_\odot$, up to a 1d degeneracy between the mass and velocity.

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like halos using the ELVIS suite of N-body simulations. We identify 25 surviving LMC-mass subhalos, and investigate the lower-mass satellites that were associated with these subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z=0 satellite population of MW/M31 halos were in a surviving LMC-group before falling into the MW/M31 halo. This fraction can vary between 1% and 25%, being higher for groups with higher mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group center become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z=0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs are genuine members, then the LMC-group probably fell into the MW very recently, <2 Gyr ago. If the connection with the LMC is established with follow-up velocity measurements, these "satellites of satellites" represent prime candidates to study the affects of group pre-processing on lower mass dwarfs.

We present a search for gamma-ray emission from the direction of the newly discovered dwarf galaxy Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between ~2-10 GeV and is consistent with annihilation of dark matter for particle masses less than a few x 10^2 GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and taking into account trials factors, we detect emission with p-value less than 9.8 x 10^-5 (>3.7 sigma). An alternative, model-independent treatment of background reduces the significance, raising the p-value to 9.7 x 10^-3 (2.3 sigma). Even in this case, however, Reticulum II has the most significant gamma-ray signal of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for annihilation.

There is hope to discover dark matter subhalos free of stars (predicted by the current theory of structure formation) by observing gaps they produce in tidal streams. In fact, this is the most promising technique for dark substructure detection and characterization as such gaps grow with time, magnifying small perturbations into clear signatures observable by ongoing and planned Galaxy surveys. To facilitate such future inference, we develop a comprehensive framework for studies of the growth of the stream density perturbations. Starting with simple assumptions and restricting to streams on circular orbits, we derive analytic formulae that describe the evolution of all gap properties (size, density contrast etc) at all times. We uncover complex, previously unnoticed behavior, with the stream initially forming a density enhancement near the subhalo impact point. Shortly after, a gap forms due to the relative change in period induced by the subhalo's passage. There is an intermediate regime where the gap grows linearly in time. At late times, the particles in the stream overtake each other, forming caustics, and the gap grows like $\sqrt{t}$. In addition to the secular growth, we find that the gap oscillates as it grows due to epicyclic motion. We compare this analytic model to N-body simulations and find an impressive level of agreement. Importantly, when analyzing the observation of a single gap we find a large degeneracy between the subhalo mass, the impact geometry and kinematics, the host potential and the time since flyby.

We present the first detailed quantitative study of the stellar populations of the Sagittarius (Sgr) streams within the Stripe 82 region, using photometric and spectroscopic observations from the Sloan Digital Sky Survey (SDSS). The star formation history (SFH) is determined separately for the bright and faint Sgr streams, to establish whether both components consist of a similar stellar population mix or have a distinct origin. Best fit SFH solutions are characterised by a well-defined, tight sequence in age-metallicity space, indicating that star formation occurred within a well-mixed, homogeneously enriched medium. Star formation rates dropped sharply at an age of ~5-7 Gyr, possibly related to the accretion of Sgr by the MW. Finally, the Sgr sequence displays a change of slope in age-metallicity space at an age between 11-13 Gyr consistent with the Sgr alpha-element knee, indicating that supernovae type Ia started contributing to the abundance pattern ~1-3 Gyr after the start of star formation. Results for both streams are consistent with being drawn from the parent Sgr population mix, but at different epochs. The SFH of the bright stream starts from old, metal-poor populations and extends to a metallicity of [Fe/H]~-0.7, with peaks at ~7 and 11 Gyr. The faint SFH samples the older, more metal-poor part of the Sgr sequence, with a peak at ancient ages and stars mostly with [Fe/H]<-1.3 and age>9 Gyr. Therefore, we argue in favour of a scenario where the faint stream consists of material stripped i) earlier, and ii) from the outskirts of the Sgr dwarf.

Photometry alone is not sufficient to unambiguously distinguish between ultra-faint star clusters and dwarf galaxies because of their overlap in morphological properties. Here we report on VLT/GIRAFFE spectra of candidate member stars in two recently discovered ultra-faint satellites Reticulum 2 and Horologium 1, obtained as part of the ongoing Gaia-ESO Survey. We identify 18 members in Reticulum 2 and 5 in Horologium 1. We find Reticulum 2 to have a velocity dispersion of ~3.22 km/s, implying a M/L ratio of ~ 500. We have inferred stellar parameters for all candidates and we find Reticulum 2 to have a mean metallicity of [Fe/H] = -2.46+/-0.1, with an intrinsic dispersion of ~ 0.29, and is alpha-enhanced to the level of [alpha/Fe]~0.4. We conclude that Reticulum 2 is a dwarf galaxy. We also report on the serendipitous discovery of four stars in a previously unknown stellar substructure near Reticulum 2 with [Fe/H] ~ -2 and V_hel ~ 220 km/s, far from the systemic velocity of Reticulum 2. For Horologium 1 we infer a velocity dispersion of 4.9^{+2.8}_{-0.9} km/s and a consequent M/L ratio of ~ 600, leading us to conclude that Horologium 1 is also a dwarf galaxy. Horologium 1 is slightly more metal-poor than Reticulum 2 [Fe/H] = -2.76 +/- 0.1 and is similarly alpha-enhanced: [alpha/Fe] ~ 0.3. Despite a large error-bar, we also measure a significant spread of metallicities of 0.17 dex which strengthen the evidence that Horologium 1 is indeed a dwarf galaxy. The line-of-sight velocity of Reticulum 2 is offset by some 100 km/s from the prediction of the orbital velocity of the LMC, thus making its association with the Cloud uncertain. However, at the location of Horologium 1, both the backward integrated LMC's orbit and the LMC's halo are predicted to have radial velocities similar to that of the dwarf. Therefore, it is likely that Horologium 1 is or once was a member of the Magellanic Family.

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of 182 stellar targets along the line of sight to the newly-discovered `ultrafaint' object Reticulum 2 (Ret 2). For 38 of these targets, the spectra are sufficient to provide simultaneous estimates of line-of-sight velocity ($v_{\rm los}$, median random error $\delta_{v_{\rm los}}=1.3$ km s$^{-1}$), effective temperature ($T_{\rm eff}$, $\delta_{T_{\rm eff}}=464$ K), surface gravity ($\log g$, $\delta_{\rm logg}=0.54$ dex) and iron abundance ([Fe/H], $\delta_{\mathrm{[Fe/H]}}=0.45$ dex). We use these results to confirm 18 stars as members of Ret 2. From the member sample we estimate a velocity dispersion of $\sigma_{v_{\rm los}}=3.6_{-0.6}^{+0.9}$ km s$^{-1}$ about a mean of $\langle v_{\rm los}\rangle =64.8_{-1.0}^{+1.1}$ km s$^{-1}$ in the solar rest frame ($\sim -90.9$ km s$^{-1}$ in the Galactic rest frame), and a metallicity dispersion of $\sigma_{\rm [Fe/H]}=0.50_{-0.13}^{+0.17}$ dex about a mean of $\langle \mathrm{[Fe/H]} \rangle =-2.67_{-0.34}^{+0.34}$. These estimates marginalize over possible velocity and metallicity gradients, which are consistent with zero. Our results place Ret 2 on chemodynamical scaling relations followed by the Milky Way's dwarf-galactic satellites. Under assumptions of dynamic equilibrium and negligible contamination from binary stars---both of which must be checked with deeper imaging and repeat spectroscopic observations---the estimated velocity dispersion suggests a dynamical mass of $M(R_{\rm h})\approx 5R_{\rm h}\sigma_{v_{\rm los}}^2/(2G)=2.4_{-0.8}^{+1.3}\times 10^5$ $M_{\odot}$ enclosed within projected halflight radius $R_{\rm h}\sim 32$ pc, with mass-to-light ratio $\approx 2M(R_{\rm h})/L_{\rm V}=462_{-157}^{+264}$ in solar units.

We have used the publicly released Dark Energy Survey data to hunt for new satellites of the Milky Way in the Southern hemisphere. Our search yielded a large number of promising candidates. In this paper, we announce the discovery of 9 new unambiguous ultra-faint objects, whose authenticity can be established with the DES data alone. Based on the morphological properties, three of the new satellites are dwarf galaxies, one of which is located at the very outskirts of the Milky Way, at a distance of 380 kpc. The remaining 6 objects have sizes and luminosities comparable to the Segue~1 satellite and can not be classified straightforwardly without follow-up spectroscopic observations. The satellites we have discovered cluster around the LMC and the SMC. We show that such spatial distribution is unlikely under the assumption of isotropy, and, therefore, conclude that at least some of the new satellites must have been associated with the Magellanic Clouds in the past.

Extremely metal poor (XMP) galaxies are known to be very rare, despite the large numbers of low-mass galaxies predicted by the local galaxy luminosity function. This paper presents a sub-sample of galaxies that were selected via a morphology-based search on SDSS images with the aim of finding these elusive XMP galaxies. By using the recently discovered extremely metal-poor galaxy, Leo P, as a guide, we obtained a collection of faint, blue systems, each with isolated HII regions embedded in a diffuse continuum, that have remained optically undetected until now. Here we show the first results from optical spectroscopic follow-up observations of 12 of ~100 of these blue, diffuse dwarf (BDD) galaxies yielded by our search algorithm. Oxygen abundances were obtained via the direct method for eight galaxies, and found to be in the range 7.45<12+log(O/H)<8.0, with two galaxies being classified as XMPs. All BDDs were found to currently have a young star-forming population (<10 Myr) and relatively high ionisation parameters of their HII regions. Despite their low luminosities (-11< M_B < -18) and low surface brightnesses (~23-25 mag arcsec^-2), the galaxies were found to be actively star-forming, with current star-formation rates between 0.0003 and 0.078 Msol/yr. From our current subsample, BDD galaxies appear to be a population of non-quiescent dwarf irregular (dIrr) galaxies, or the diffuse counterparts to blue compact galaxies (BCDs) and as such may bridge the gap between these two populations. Our search algorithm demonstrates that morphology-based searches are successful in uncovering more diffuse metal-poor star-forming galaxies, which traditional emission-line based searches overlook.

We present a model for producing tidal streams from disrupting progenitors in arbitrary potentials, utilizing the idea that the majority of stars escape from the progenitor's two Lagrange points. The method involves releasing test particles at the Lagrange points as the satellite orbits the host and dynamically evolving them in the potential of both host and progenitor. The method is sufficiently fast to allow large-dimensional parameter exploration using Monte Carlo methods. We provide the first direct modelling of 6-D stream observations -- assuming a stream rather than an orbit -- by applying our methods to GD-1. This is a kinematically cold stream spanning $60^{\circ}$ of the sky and residing in the outer Galaxy $\approx 15$ kpc distant from the centre. We assume the stream moves in a flattened logarithmic potential characterised by an asymptotic circular velocity $v_0$ and a flattening $q$. We recover values of normalisation $v_0$ = $227.2^{+15.6}_{-18.2}$ kms$^{-1}$ and flattening $q$ = $0.91^{+0.04}_{-0.1}$, if the stream is assumed to leading, and $v_0$ = $226.5^{+17.9}_{-17.0}$ kms$^{-1}$, $q$ = $0.90^{+0.05}_{-0.09}$, if it is assumed to be trailing. This can be compared to the values $v_0 = 224 \pm 13$ kms$^{-1}$ and $q= 0.87^{+0.07}_{-0.04}$ obtained by Koposov et al (2010) using the simpler technique of orbit fitting. Although there are differences between stream and orbit fitting, we conclude that orbit fitting can provide accurate results given the current quality of the data, at least for this kinematically cold stream in this logarithmic model of the Galaxy.

We present a census of blue horizontal branch (BHB) and blue straggler (BS) stars belonging to dwarf galaxies and globular clusters, and compare these counts to that of the Milky Way stellar halo. We find, in agreement with earlier studies, that the ratio of BS-to-BHB stars in these satellite populations is dependent on stellar mass. Dwarf galaxies show an increasing BS-to-BHB ratio with luminosity. In contrast, globular clusters display the reverse trend, with N_BS/N_BHB (< 1) decreasing with luminosity. The faintest (L < 10^5 L_Sun) dwarfs have similar numbers of BS and BHB stars (N_BS/N_BHB ~ 1), whereas more massive dwarfs tend to be dominated by BS stars (N_BS/N_BHB ~ 2-40). We find that the BS-to-BHB ratio in the stellar halo is relatively high (N_BS/N_BHB ~ 5-6), and thus inconsistent with the low ratios found in both ultra-faint dwarfs and globular clusters. Our results favour more massive dwarfs as the dominant "building blocks" of the stellar halo, in good agreement with current predictions from LambdaCDM models.

Using the RAdial Velocity Experiment fourth data release (RAVE DR4), and a new metallicity calibration that will be also taken into account in the future RAVE DR5, we investigate the existence and the properties of super-solar metallicity stars ([M/H] > +0.1 dex) in the sample, and in particular in the Solar neighbourhood. We find that RAVE is rich in super-solar metallicity stars, and that the local metallicity distribution function declines remarkably slowly up to +0.4 dex. Our results show that the kinematics and height distributions of the super-solar metallicity stars are identical to those of the [M/H] < 0 thin-disc giants that we presume were locally manufactured. The eccentricities of the super-solar metallicity stars indicate that half of them are on a roughly circular orbit (e < 0.15), so under the assumption that the metallicity of the interstellar medium at a given radius never decreases with time, they must have increased their angular momenta by scattering at corotation resonances of spiral arms from regions far inside the Solar annulus. The likelihood that a star will migrate radially does not seem to decrease significantly with increasing amplitude of vertical oscillations within range of oscillation amplitudes encountered in the disc.

We present the results of a deep, wide-area variability survey in the Southern hemisphere, the first of its kind. As part of the Catalina Sky Surveys, the Siding Spring Survey (SSS) has covered $14,800$ square degrees in the declination range of $-75^{\circ}\leq\delta\leq-15^{\circ}$. To mine the enormous SSS dataset efficiently we have developed two algorithms: Automatic Period Selection (APS) and Automatic Fourier Decomposition (AFD), which aim to sharpen the period estimation and produce robust lightcurve models. Armed with the APS and AFD outputs we classify $10,540$ ab-type RR Lyrae (RRab) stars ($\sim$90% of which are new) across the Southern sky. As well as the positional information we supply photometric metallicities, and unreddened distances. For the RRab stars in the halo, a study of the photometric metallicity distribution reveals a nearly Gaussian shape with a mean metallicity of ${\rm [Fe/H]}=-1.4$ dex and a dispersion of $0.3$ dex. A spatial study of the RRab metallicities shows no significant radial gradient in the first $\sim7$ kpc from the Galaxy center. However, further out, a small negative gradient is clearly present. This is complemented by a very obvious correlation of the mean RR Lyrae metallicity with distance above the Galactic plane, $z$. We have also carried out an initial substructure search using the discovered RRab, and present the properties of the candidates with significance greater than $2 \sigma$. Most prominent among these is a southern extension of the Sagittarius dwarf galaxy's stream system, reaching down to declinations $\sim -40\deg$.

We show that a combination of infrared photometry from WISE and 2MASS surveys can yield highly pure samples of M giant stars. We take advantage of the new WISE$\cap$2MASS M giant selection to trace the Sagittarius trailing tail behind the Galactic disk in the direction of the anti-centre. The M giant candidates selected via broad-band photometry are confirmed spectroscopically using AAOmega on the AAT in 3 fields around the extremity of the Sgr trailing tail in the Southern Galactic hemisphere. We demonstrate that at the Sgr longitude $\tilde \Lambda_{\odot} = 204^{\circ}$, the line-of-sight velocity of the trailing tail starts to deviate from the track of the Law & Majewski (2010) model, confirming the prediction of Belokurov et al. (2014). This discovery serves to substantiate the measurement of low differential orbital precession of the Sgr stream which in turn may imply diminished dark matter content within 100 kpc.

Motivated by recent observations of the Sagittarius stream, we devise a rapid algorithm to generate faithful representations of the centroids of stellar tidal streams formed in a disruption of a progenitor of an arbitrary mass in an arbitrary potential. Our method works by releasing swarms of test particles at the Lagrange points around the satellite and subsequently evolving them in a combined potential of the host and the progenitor. We stress that the action of the progenitor's gravity is crucial to making streams that look almost indistinguishable from the N-body realizations, as indeed ours do. The method is tested on mock stream data in three different Milky Way potentials with increasing complexity, and is shown to deliver unbiased inference on the Galactic mass distribution out to large radii. When applied to the observations of the Sagittarius stream, our model gives a natural explanation of the stream's apocentric distances and the differential orbital precession. We, therefore, provide a new independent measurement of the Galactic mass distribution beyond 50 kpc. The Sagittarius stream model favours a light Milky Way with the mass 4.1 +/- 0.4 x 10^11 M_sun at 100 kpc, which can be extrapolated to 5.6 +/- 1.2 x 10^11 M_sun at 200 kpc. Such a low mass for the Milky Way Galaxy is in good agreement with estimates from the kinematics of halo stars and from the satellite galaxies (once Leo I is removed from the sample). It entirely removes the "Too Big To Fail Problem".

We demonstrate that the asymmetric distribution of M31 satellites cannot be produced by tides from the Milky Way as such effects are too weak. However, loosely bound associations and groups of satellites can fall into larger haloes and give rise to asymmetries. We compute the survival times for such associations. We prove that the survival time is always shortest in Keplerian potentials, and can be ~ 3 times longer in logarithmic potentials. We provide an analytical formula for the dispersal time in terms of the size and velocity dispersion of the infalling structure. We show that, if an association of ~10 dwarfs fell into the M31 halo, its present aspect would be that of an asymmetric disk of satellites. We also discuss the case of cold substructure in the Andromeda II and Ursa Minor dwarfs.

We explore the Triangulum-Andromeda (TriAnd) overdensity in the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) and SEGUE (the Sloan Extension for Galactic Understanding and Exploration) spectroscopic surveys. Milky Way main sequence turn-off stars in the SPLASH survey reveal that the TriAnd overdensity and the recently discovered PAndAS stream (Martin et al. 2014) share a common heliocentric distance (D ~ 20 kpc), position on the sky, and line-of-sight velocity (V_GSR ~ 50 km/s). Similarly, A-type, giant, and main sequence turn-off stars selected from the SEGUE survey in the vicinity of the Segue 2 satellite show that TriAnd is prevalent in these fields, with a velocity and distance similar to Segue 2. The coincidence of the PAndAS stream and Segue 2 satellite in positional and velocity space to TriAnd suggests that these substructures are all associated, and may be a fossil record of group-infall onto the Milky Way halo. In this scenario, the Segue 2 satellite and PAndAS stream are "satellites of satellites", and the large, metal-rich TriAnd overdensity is the remains of the group central.

In the rest frame of the Local Group (LG), the total momentum of the Milky Way (MW) and Andromeda (M31) should balance to zero. We use this fact to constrain new solutions for the solar motion with respect to the LG centre-of-mass, the total mass of the LG, and the individual masses of M31 and the MW. Using the set of remote LG galaxies at $>350$ kpc from the MW and M31, we find that the solar motion has amplitude $V_{\odot}=299\pm 15 {\rm ~km~s^{-1}}$ in a direction pointing toward galactic longitude $l_{\odot}=98.4^{\circ}\pm 3.6^{\circ}$ and galactic latitude $b_{\odot}=-5.9^{\circ}\pm 3.0^{\circ}$. The velocities of M31 and the MW in this rest frame give a direct measurement of their mass ratio, for which we find $\log_{10} (M_{\rm M31}/M_{\rm MW})=0.36 \pm 0.29$. We combine these measurements with the virial theorem to estimate the total mass within the LG as $M_{\rm LG}=(2.5\pm 0.4)\times 10^{12}~{\rm M}_{\odot}$. Our value for $M_{\rm LG}$ is consistent with the sum of literature values for $M_{\rm MW}$ and $M_{\rm M31}$. This suggests that the mass of the LG is almost entirely located within the two largest galaxies rather than being dispersed on larger scales or in a background medium. The outskirts of the LG are seemingly rather empty. Combining our measurement for $M_{\rm LG}$ and the mass ratio, we estimate the individual masses of the MW and M31 to be $M_{\rm MW}=(0.8\pm 0.5)\times 10^{12}~{\rm M}_{\odot}$ and $M_{\rm M31}=(1.7\pm 0.3)\times 10^{12}~{\rm M}_{\odot}$, respectively. Our analysis favours M31 being more massive than the MW by a factor of $\sim$2.3, and the uncertainties allow only a small probability (9.8%) that the MW is more massive. This is consistent with other properties such as the maximum rotational velocities, total stellar content, and numbers of globular clusters and dwarf satellites, which all suggest that $M_{\rm M31}/M_{\rm MW}>1$.

We employ abundances from the Sloan Digital Sky Survey (SDSS) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) to study the alpha-element distribution of the stellar members of the Sagittarius stream. To test the reliability of SDSS/SEGUE abundances for the study of Sagittarius, we select high-likelihood samples tracing the different components of the Milky Way, and recover known literature alpha-element distributions. Using selection criteria based on the spatial position, radial velocity, distance and colours of individual stars, we obtain a robust sample of Sagittarius-stream stars. The alpha-element distribution of the Sagittarius stream forms a narrow sequence at intermediate metallicities with a clear turn-down, consistent with the presence of an alpha-element "knee". This is the first time that the alpha-element knee of the Sagittarius dwarf galaxy has been detected. Fitting a toy model to our data, we determine that the alpha-knee in Sagittarius takes place at [Fe/H]=-1.27pm0.05, only slightly less metal-poor than the knee in the Milky Way. This indicates that a small number of Sagittarius-like galaxies could have contributed significantly to the build-up of the Milky Way's stellar halo system at ancient times.

We explore the properties of lensing galaxies and lensed faint sources at redshifts between 1.5 and 3.0. Our sample consists of 9 carefully selected strongly-lensed galaxies discovered by the CASSOWARY survey in the Sloan Digital Sky Survey (SDSS) data. We find that, despite some limitations of the original SDSS data, the homogeneous sample of lensing systems can provide a useful insight into lens and source properties. We also explore the limitations of using low-resolution data to model and analyse galaxy-galaxy lensing. We derive the relative alignment of mass and light in fitted lens profiles. The range of magnification extends above 5, hence we are able to analyse potentially small and low-mass galaxies at high redshifts. We confirm the likely evolution of the size-luminosity relation for blue star-forming galaxies as a function of redshift

We present ~47,000 periodic variables found during the analysis of 5.4 million variable star candidates within a 20,000 square degree region covered by the Catalina Surveys Data Release-1 (CSDR1). Combining these variables with type-ab RR Lyrae from our previous work, we produce an on-line catalog containing periods, amplitudes, and classifications for ~61,000 periodic variables. By cross-matching these variables with those from prior surveys, we find that > 90% of the ~8,000 known periodic variables in the survey region are recovered. For these sources we find excellent agreement between our catalog and prior values of luminosity, period and amplitude, as well as classification. We investigate the rate of confusion between objects classified as contact binaries and type-c RR Lyrae (RRc's) based on periods, colours, amplitudes, metalicities, radial velocities and surface gravities. We find that no more than few percent of these variables in these classes are misidentified. By deriving distances for this clean sample of ~5,500 RRc's, we trace the path of the Sagittarius tidal streams within the Galactic halo. Selecting 146 outer-halo RRc's with SDSS radial velocities, we confirm the presence of a coherent halo structure that is inconsistent with current N-body simulations of the Sagittarius tidal stream. We also find numerous long-period variables that are very likely associated within the Sagittarius tidal streams system. Based on the examination of 31,000 contact binary light curves we find evidence for two subgroups exhibiting irregular lightcurves. One subgroup presents significant variations in mean brightness that are likely due to chromospheric activity. The other subgroup shows stable modulations over more than a thousand days and thereby provides evidence that the O'Connell effect is not due to stellar spots.

We report the discovery of a narrow stellar stream crossing the constellations of Sculptor and Fornax in the Southern celestial hemisphere. The portion of the stream detected in the Data Release 1 photometry of the ATLAS survey is at least 12 degrees long, while its width is $\approx$ 0.25 deg. The Color Magnitude Diagram of this halo sub-structure is consistent with a metal-poor [Fe/H] $\lesssim -1.4$ stellar population located at a heliocentric distance of 20 $\pm$ 2 kpc. There are three globular clusters that could tentatively be associated with the stream: NGC 7006, NGC 7078 (M15) and Pyxis, but NGC 7006 and 7078 seem to have proper motions incompatible with the stream orbit.

We use A-type stars selected from Sloan Digital Sky Survey data release 9 photometry to measure the outer slope of the Milky Way stellar halo density profile beyond 50 kpc. A likelihood-based analysis is employed that models the ugr photometry distribution of blue horizontal branch (BHB) and blue straggler (BS) stars. In the magnitude range, 18.5 < g < 20.5, these stellar populations span a heliocentric distance range of: 10 kpc < D_BS < 75 kpc, 40 kpc < D_BHB < 100 kpc. Contributions from contaminants, such as QSOs, and the effect of photometric uncertainties, are also included in our modeling procedure. We find evidence for a very steep outer halo profile, with power-law index alpha ~ 6 beyond Galactocentric radii r=50 kpc, and even steeper slopes favored (alpha ~ 6-10) at larger radii. This result holds true when stars belonging to known overdensities, such as the Sagittarius stream, are included or excluded. We show that, by comparison to numerical simulations, stellar halos with shallower slopes at large distances tend to have more recent accretion activity. Thus, it is likely that the Milky Way has undergone a relatively quiet accretion history over the past several Gyr. Our measurement of the outer stellar halo profile may have important implications for dynamical mass models of the Milky Way, where the tracer density profile is strongly degenerate with total mass-estimates.

We announce the discovery of a new Galactic companion found in data from the ESO VST ATLAS survey, and followed up with deep imaging on the 4m William Herschel Telescope. The satellite is located in the constellation of Crater (the Cup) at a distance of $\sim$ 170 kpc. Its half-light radius is $r_h=30$ pc and its luminosity is $M_V=-5.5$. The bulk of its stellar population is old and metal-poor. We would probably have classified the newly discovered satellite as an extended globular cluster were it not for the presence of a handful of Blue Loop stars and a sparsely populated Red Clump. The existence of the core helium burning population implies that star-formation occurred in Crater perhaps as recently as 400 Myr ago. No globular cluster has ever accomplished the feat of prolonging its star-formation by several Gyrs. Therefore, if our hypothesis that the blue bright stars in Crater are Blue Loop giants is correct, the new satellite should be classified as a dwarf galaxy with unusual properties. Note that only ten degrees to the North of Crater, two ultra-faint galaxies Leo IV and V orbit the Galaxy at approximately the same distance. This hints that all three satellites may once have been closely associated before falling together into the Milky Way halo.

We map the large-scale sub-structure in the Galactic stellar halo using accurate 3D positions of ~14,000 RR Lyrae reported by the Catalina Sky Survey. In the heliocentric distance range of 10-25 kpc, in the region of the sky approximately bounded by 30{\deg} < l < 55{\deg} and -45{\deg} < b < -25{\deg}, there appears to be a strong excess of RRab stars. This overdensity, peaking at 18 kpc, is most likely associated with the so-called Hercules-Aquila Cloud, previously detected using Main Sequence tracers at similar distances in the Sloan Digital Sky Survey data. Our analysis of the period-amplitude distribution of RR Lyrae in this region indicates that the HAC is dominated by the Oosterhoff I type population. By comparing the measured RR Lyrae number density to models of a smooth stellar halo, we estimate the significance of the observed excess and provide an updated estimate of the total luminosity of the Cloud's progenitor.

Using a variety of stellar tracers -- blue horizontal branch stars, main-sequence turn-off stars and red giants -- we follow the path of the Sagittarius (Sgr) stream across the sky in Sloan Digital Sky Survey data. Our study presents new Sgr debris detections, accurate distances and line-of-sight velocities that together help to shed new light on the puzzle of the Sgr tails. For both the leading and the trailing tail, we trace the points of their maximal extent, or apo-centric distances, and find that they lie at $R^L$ = 47.8 $\pm$ 0.5 kpc and $R^T$ = 102.5 $\pm$ 2.5 kpc respectively. The angular difference between the apo-centres is 93.2 $\pm$ 3.5 deg, which is smaller than predicted for logarithmic haloes. Such differential orbital precession can be made consistent with models of the Milky Way in which the dark matter density falls more quickly with radius. However, currently, no existing Sgr disruption simulation can explain the entirety of the observational data. Based on its position and radial velocity, we show that the unusually large globular cluster NGC 2419 can be associated with the Sgr trailing stream. We measure the precession of the orbital plane of the Sgr debris in the Milky Way potential and show that, surprisingly, Sgr debris in the primary (brighter) tails evolves differently to the secondary (fainter) tails, both in the North and the South.

We employ a large sample of 20171 optically-selected groups and clusters at 0.15 < z < 0.4 in the SDSS to investigate how the stacked stellar mass fraction varies across a wide range of total mass, $M_{500}$. Our study improves upon previous observational studies in a number of important ways, including having a much larger sample size, an explicit inclusion of the intracluster light (ICL) component, and a thorough examination of the accuracy of our total mass estimates via comparisons to simulations and weak lensing observations. We find that the stellar mass fraction depends only weakly on total mass and that the contribution of ICL to the total stellar mass fraction is significant (typically 20-40 per cent). Both of these findings are in excellent accordance with the predictions of cosmological simulations. Under the assumption of a Chabrier (Salpeter) IMF, the derived star formation efficiency ($f_{star}$/$f_{b}$, where $f_b=\Omega_b/\Omega_m$) is relatively low at 8 per cent (14 per cent) and is consistent with the global star formation efficiency of semi-analytic models that reproduce the galaxy stellar mass function. When our measured stellar mass fractions are combined with the observed relation between hot gas mass fraction and total mass from X-ray observations, our results imply that galaxy groups have significantly lower baryon fractions than massive clusters. Ejection of gas due to energetic AGN feedback (most likely at high redshift) provides a plausible mechanism for explaining the trends we observe.

We investigate the relationship between the halo mass, M_200, and concentration, c, for a sample of 26 group- and cluster-scale strong gravitational lenses. In contrast with previous results, we find that these systems are only ~ 0.1 dex more over-concentrated than similar-mass halos from dark matter simulations; the concentration of a halo with M_200 = 10^14 M_sun is log c = 0.78\pm0.05, while simulations of halos with this mass at similar redshifts (z ~ 0.4) predict log c ~ 0.56 - 0.71. We also find that we are unable to make informative inference on the slope of the M_200-c relation in spite of our large sample size; we note that the steep slopes found in previous studies tend to follow the slope in the covariance between M_200 and c, indicating that these results may be measuring the scatter in the data rather than the intrinsic signal. Furthermore, we conclude that our inability to constrain the M_200-c slope is due to a limited range of halo masses, as determined by explicitly modelling our halo mass distribution, and we suggest that other studies may be producing biased results by using an incorrect distribution for their halo masses.

We construct simple triaxial generalisations of Navarro-Frenk-White haloes. The models have elementary gravitational potentials, together with a density that is cusped like 1/r at small radii and falls off like 1/r^3 at large radii. The ellipticity varies with radius in a manner that can be tailored to the user's specification. The closed periodic orbits in the planes perpendicular to the short and long axes of the model are well-described by epicyclic theory, and can be used as building blocks for long-lived discs. As an application, we carry out the simulations of thin discs of satellites in triaxial dark halo potentials. This is motivated by the recent claims of an extended, thin disc of satellites around the M31 galaxy with a vertical rms scatter of ~12 kpc and a radial extent of ~ 300 kpc (Ibata et al. 2013). We show that a thin satellite disc can persist over cosmological times if and only if it lies in the planes perpendicular to the long or short axis of a triaxial halo, or in the equatorial or polar planes of a spheroidal halo. In any other orientation, then the disc thickness doubles on ~5 Gyr timescales and so must have been born with an implausibly small vertical scaleheight.

From the archaeological point of view, the local dwarf galaxies are unique objects in which the imprint of the conditions that shaped the early structure formation can be studied today at high resolution. Over the last decade, this new window into the high redshift Universe has started to be exploited using deep wide-field imaging, high resolution spectroscopy and cutting edge N-body and hydro-dynamical simulations. We review the recent advances in the observational studies of the Milky Way dwarf galaxies, with the aim to understand the properties of the population as a whole and to assist an objective comparison between the models and the data.

We exploit the group environment of the CAmbridge Sloan Survey Of Wide ARcs in the skY (CASSOWARY) z=0.3 lens J2158+0257 to measure the group dynamical mass as a complement to the central dynamical and lensing mass constraints. Follow-up spectroscopy of candidate group members is performed using VLT/FORS2. From the resulting N=21 confirmed members we measure the group dynamical mass by calibrating an analytic tracer mass estimator with cosmological simulations. The luminosity weighted line-of-sight velocity dispersion and the Einstein radius of the lens are used as mass probes in the inner regions of the galaxy. Combining these three observational probes allows us to independently constrain the mass and concentration of the dark matter halo, in addition to the total stellar mass of the central galaxy. We find a dark matter halo in remarkably good agreement with simulations (log_10 M_200/M_sun = 14.2 +/- 0.2, c_200 = 4.4 (+1.6, -1.4)) and a stellar mass-to-light ratio which favors a Salpeter initial mass function ((M/L)* = 5.7 +/- 1.2). Our measurement of a normal halo concentration suggests that there is no discrepancy between simulations and observations at the group mass scale. This is in contrast to the cluster mass scale for which a number of studies have claimed over-concentrated halos. While the halo mass is robustly determined, and the halo concentration is not significantly affected by systematics, the resulting stellar mass-to-light ratio is sensitive to the choice of stellar parameters, such as density profile and velocity anisotropy.

We analyse a 10,000 star spectroscopic survey, focused on Galactic thick disk stars typically 2-5 kpc from the Sun, carried out using the AAOmega Spectrograph on the AAT. We develop methods for completeness-correction of the survey based on SDSS photometry, and we derive star distances using an improved isochrone-fitting method with accuracies better than 10%. We determine the large-scale kinematic (Vphi=172 km/s and sigma(phi,r,z) = (49,51,40) km/s), and abundance properties of the thick disk, showing these representative values are a fair description within about 3 kpc of the sun, and in the range 1-3 kpc from the Galactic Plane. We identify a substantial overdensity in lines of sight towards the inner Galaxy, with metallicity [Fe/H] -1 dex, and higher line of sight velocities than the thick disk, localised along the direction (l,b)=(48,-26). This overdensity appears to be towards, but closer than, the known Hercules-Aquila halo overdensity and may be related to the Humphreys-Larsen inner galaxy thick disk asymmetry.

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

Using multi-epoch observations of the Stripe 82 region done by Sloan Digital Sky Survey, we measure precise statistical proper motions of the stars in the Sagittarius stellar stream. The multi-band photometry and SDSS radial velocities allow us to efficiently select Sgr members and thus enhance the proper motion precision to ~0.1 mas/yr. We measure separately the proper motion of a photometrically selected sample of the main sequence turn-off stars, as well as of a spectroscopically selected Sgr giants. The data allow us to determine the proper motion separately for the two Sgr streams in the South found in Koposov et al.(2012). Together with the precise velocities from SDSS, our proper motion provide exquisite constraints of the 3-D motions of the stars in the Sgr streams.

We use the Bullock & Johnston suite of simulations to study the density profiles of L*-type galaxy stellar haloes. Observations of the Milky Way and M31 stellar haloes show contrasting results: the Milky Way has a `broken' profile, where the density falls off more rapidly beyond ~ 25 kpc, while M31 has a smooth profile out to 100 kpc with no obvious break. Simulated stellar haloes, built solely by the accretion of dwarf galaxies, also exhibit this behavior: some haloes have breaks, while others don't. The presence or absence of a break in the stellar halo profile can be related to the accretion history of the galaxy. We find that a break radius is strongly related to the build up of stars at apocentres. We relate these findings to observations, and find that the `break' in the Milky Way density profile is likely associated with a relatively early (~ 6-9 Gyr ago) and massive accretion event. In contrast, the absence of a break in the M31 stellar halo profile suggests that its accreted satellites have a wide range of apocentres. Hence, it is likely that M31 has had a much more prolonged accretion history than the Milky Way.

We present the analysis of 12227 type-ab RR Lyrae found among the 200 million public lightcurves in the Catalina Surveys Data Release 1 (CSDR1). These stars span the largest volume of the Milky Way ever surveyed with RR Lyrae, covering ~20,000 square degrees of the sky (0 < RA < 360, -22 < Dec < 65 deg) to heliocentric distances of up to 60kpc. Each of the RR Lyrae are observed between 60 and 419 times over a six-year period. Using period finding and Fourier fitting techniques we determine periods and apparent magnitudes for each source. We find that the periods at generally accurate to sigma = 0.002% by comparison with 2842 previously known RR Lyrae and 100 RR Lyrae observed in overlapping survey fields. We photometrically calibrate the light curves using 445 Landolt standard stars and show that the resulting magnitudes are accurate to ~0.05 mags using SDSS data for ~1000 blue horizontal branch stars and 7788 of the RR Lyrae. By combining Catalina photometry with SDSS spectroscopy, we analyze the radial velocity and metallicity distributions for > 1500 of the RR Lyrae. Using the accurate distances derived for the RR Lyrae, we show the paths of the Sagittarius tidal streams crossing the sky at heliocentric distances from 20 to 60 kpc. By selecting samples of Galactic halo RR Lyrae, we compare their velocity, metallicity, and distance with predictions from a recent detailed N-body model of the Sagittarius system. We find that there are some significant differences between the distances and structures predicted and our observations.

We build a sample of distant (D > 80 kpc) stellar halo stars with measured radial velocities. Faint ($20 < g <22$) candidate blue horizontal branch (BHB) stars were selected using the deep, but wide, multi-epoch Sloan Digital Sky Survey photometry. Follow-up spectroscopy for these A-type stars was performed using the VLT-FORS2 instrument. We classify stars according to their Balmer line profiles, and find 7 are bona fide BHB stars and 31 are blue stragglers (BS). Owing to the magnitude range of our sample, even the intrinsically fainter BS stars can reach out to D ~ 90 kpc. We complement this sample of A-type stars with intrinsically brighter, intermediate-age, asymptotic giant branch stars. A set of 4 distant cool carbon stars is compiled from the literature and we perform spectroscopic follow-up on a further 4 N-type carbon stars using the WHT-ISIS instrument. Altogether, this provides us with the largest sample to date of individual star tracers out to r ~ 150 kpc. We find that the radial velocity dispersion of these tracers falls rapidly at large distances and is surprisingly cold (sigma_r ~ 50-60 km/s) between 100-150 kpc. Relating the measured radial velocities to the mass of the Milky Way requires knowledge of the (unknown) tracer density profile and anisotropy at these distances. Nonetheless, by assuming the stellar halo stars between 50-150 kpc have a moderate density fall-off (with power-law slope alpha < 5) and are on radial orbits (sigma^2_t/sigma^2_r < 1), we infer that the mass within 150 kpc is less than 10^12 M_solar and suggest it probably lies in the range (5-10) x 10^11 M_solar. We discuss the implications of such a low mass for the Milky Way.

We present VLT/FORS2 spectroscopy of candidate blue horizontal branch (BHB) stars in the vicinity of the Hercules ultrafaint dwarf galaxy. We identify eight convincing Hercules BHB members, and a further five stars with similar systemic velocities to that of Hercules, but ~ 0.5 kpc from the centre of the galaxy along its major axis. It is likely that these stars once belonged to Hercules, but have been tidally stripped and are now unbound. We emphasise the usefulness of looking for any gradient in the systemic velocity of this stretched system, which would further support our interpretation of the origin of its elongated and distended morphology.

We investigate the possibility of constraining the dark energy equation of state by measuring the ratio of Einstein radii in a strong gravitational lens system with two source planes. This quantity is independent of the Hubble parameter and directly measures the growth of angular diameter distances as a function of redshift. We investigate the prospects for a single double source plane system and for a forecast population of systems discovered by re-observing a population of single source lenses already known from a photometrically selected catalogue such as CASSOWARY or from a spectroscopically selected catalogue such as SLACS. We find that constraints comparable to current data-sets (15% uncertainty on the dark equation of state at 68%CL) are possible with a handful of double source plane systems. We also find that the method's degeneracy between Omega_M and w is almost orthogonal to that of CMB and BAO measurements, making this method highly complimentary to current probes.

We use distant Blue Horizontal Branch stars with Galactocentric distances 16 < r/kpc < 48 as kinematic tracers of the Milky Way dark halo. We model the tracer density as an oblate, power-law embedded within a spherical power-law potential. Using a distribution function method, we estimate the overall power-law potential and the velocity anisotropy of the halo tracers. We measure the slope of the potential to be gamma ~ 0.4 and the overall mass within 50 kpc is ~ 4 x 10^11 M_sol. The tracer velocity anisotropy is radially biased with beta ~ 0.5, which is in good agreement with local solar neighbourhood studies. Our results provide an accurate outer circular velocity profile for the Milky Way and suggest a relatively high concentration dark matter halo (c_vir ~ 20).

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 present a new catalogue of 55,121 groups and clusters centred on Luminous Red Galaxies from SDSS DR7 in the redshift range 0.15<z<0.4. We provide halo mass estimates for each of these groups derived from a calibration between the optical richness of bright galaxies (M_r<-20.5) within 1 Mpc, and X-ray-derived mass for a small subset of 129 groups and clusters with X-ray measurements. We derive the mean (stacked) surface number density profiles of galaxies as a function of total halo mass in different mass bins. We find that derived profiles can be well-described by a projected NFW profile with a concentration parameter (<c>~2.6) that is approximately a factor of two lower than that of the dark matter (as predicted by N-body cosmological simulations) and nearly independent of halo mass. Interestingly, in spite of the difference in shape between the galaxy and dark matter radial distributions, both exhibit a high degree of self-similarity. A self-consistent comparison to several recent semi-analytic models of galaxy formation indicates that: (1) beyond ~0.3 r_500 current models are able to reproduce both the shape and normalisation of the satellite profiles; and (2) within ~0.3 r_500 the predicted profiles are sensitive to the details of the satellite-BCG merger timescale calculation. The former is a direct result of the models being tuned to match the global galaxy luminosity function combined with the assumption that the satellite galaxies do not suffer significant tidal stripping, even though their surrounding dark matter haloes can be removed through this process. Combining our results with measurements of the intracluster light should provide a way to inform theoretical models on the efficacy of the tidal stripping and merging processes.

We present the star formation history (SFH) of the faintest known star-forming galaxy, Leo T, based on imaging taken with the Hubble Space Telescope (HST) Wide Field Planetary Camera 2 (WFPC2). The HST/WFPC2 color-magnitude diagram (CMD) of Leo T is exquisitely deep, extending ~ 2 magnitudes below the oldest main sequence turnoff, permitting excellent constraints on star formation at all ages. We use a maximum likelihood CMD fitting technique to measure the SFH of Leo T assuming three different sets of stellar evolution models: Padova (solar-scaled metallicity) and BaSTI (both solar-scaled and alpha-enhanced metallicities). The resulting SFHs are remarkably consistent at all ages, indicating that our derived SFH is robust to the choice of stellar evolution model. From the lifetime SFH of Leo T, we find that 50% of the total stellar mass formed prior to z ~ 1 (7.6 Gyr ago). Subsequent to this epoch, the SFH of Leo T is roughly constant until the most recent ~ 25 Myr, where the SFH shows an abrupt drop. This decrease could be due to a cessation of star formation or stellar initial mass function sampling effects, but we are unable to distinguish between the two scenarios. Overall, our measured SFH is consistent with previously derived SFHs of Leo T. However, the HST-based solution provides improved age resolution and reduced uncertainties at all epochs. The SFH, baryonic gas fraction, and location of Leo T are unlike any of the other recently discovered faint dwarf galaxies in the Local Group, and instead bear strong resemblance to gas-rich dwarf galaxies (irregular or transition), suggesting that gas-rich dwarf galaxies may share common modes of star formation over a large range of stellar mass (~ 10^5-10^9 Msun).

The Sagittarius dwarf galaxy is the archetype of a tidally disrupting system. Both leading and trailing tails can be observed across at least 180 degrees of the sky and measurements of their luminosity density profiles have recently become available. Using numerical simulations, we explore the factors that control the appearance of such profiles. We use two possible models for the Sgr progenitor. The first is a one-component Plummer model, which may represent either a dark matter free progenitor, or one in which pre-existing dark matter has already been largely stripped. The second is a two-component model in which the stars are represented by a Hernquist sphere embedded in a cosmologically modish Navarro-Frenk-White dark halo. Disruption of the models in the Milky Way galaxy provides us with two tellings of the tale of the formation of the Sgr stream. The initial disintegration of the baryons proceeds more slowly for the two-component models because of the protective cocoon of dark matter. Once this has been stripped, though, matters proceed apace. In both cases, the profiles after ~6 pericentric passages provide good matches to the observational data, but the tails are more extended for the two-component models. The leading and trailing tails are symmetric at apocentre or pericentre. At other orbital phases, asymmetries are present, as tails are compressed as they approach apocentre and stretched out as they approach pericentre. There may exist density enhancements corresponding to such pile-ups which may be observable in current survey data. We re-visit the calculation of Niederste-Ostholt et al. (2010) and slightly revise upwards the luminosity of the Sgr progenitor to 9.9-14.4 x 10^7 solar luminosities based on insights from the simulations.

We estimate the masses of elliptical galaxies out to five effective radii using planetary nebulae and globular clusters as tracers. A sample of 15 elliptical galaxies with a broad variation in mass is compiled from the literature. A distribution function-maximum likelihood analysis is used to estimate the overall potential slope, normalisation and velocity anisotropy of the tracers. We assume power-law profiles for the potential and tracer density and a constant velocity anisotropy. The derived potential power-law indices lie in between the isothermal and Keplerian regime and vary with mass: there is tentative evidence that the less massive galaxies have steeper potential profiles than the more massive galaxies. We use stellar mass-to-light ratios appropriate for either a Chabrier/KTG (Kroupa, Tout & Gilmore) or Salpeter initial mass function to disentangle the stellar and dark matter components. The fraction of dark matter within five effective radii increases with mass, in agreement with several other studies. We employ simple models to show that a combination of star formation efficiency and baryon extent are able to account for this trend. These models are in good agreement with both our measurements out to five effective radii and recent SLACS measurements within one effective radii when a universal Chabrier/KTG initial mass function is adopted.

We introduce a new maximum likelihood method to model the density profile of Blue Horizontal Branch and Blue Straggler stars and apply it to the Sloan Digital Sky Survey Data Release 8 (DR8) photometric catalogue. There are a large number (20,000) of these tracers available over an impressive 14,000 square degrees in both Northern and Southern Galactic hemispheres, and they provide a robust measurement of the shape of the Milky Way stellar halo. After masking out stars in the vicinity of the Virgo Overdensity and the Sagittarius stream, the data are consistent with a smooth, oblate stellar halo with a density that follows a broken power-law. The best fitting model has an inner slope = 2.3 and an outer slope = 4.6, together with a break radius occurring at 27 kpc and a constant halo flattening (that is, ratio of minor axis to major axis) of q = 0.6. Although a broken power-law describes the density fall-off most adequately, it is also well fit by an Einasto profile. There is no strong evidence for variations in flattening with radius, or for triaxiality of the stellar halo.

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

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 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 present the first high-resolution images of CSWA 31, a gravitational lens system observed as part of the SLUGS (Sloan Lenses Unravelled by Gemini Studies) program. These systems exhibit complex image structure with the potential to strongly constrain the mass distribution of the massive lens galaxies, as well as the complex morphology of the sources. In this paper, we describe the strategy used to reconstruct the unlensed source profile and the lens galaxy mass profiles. We introduce a prior distribution over multi-wavelength sources that is realistic as a representation of our knowledge about the surface brightness profiles of galaxies and groups of galaxies. To carry out the inference computationally, we use Diffusive Nested Sampling, an efficient variant of Nested Sampling that uses Markov Chain Monte Carlo (MCMC) to sample the complex posterior distributions and compute the normalising constant. We demonstrate the efficacy of this approach with the reconstruction of the group-group gravitational lens system CSWA 31, finding the source to be composed of five merging spiral galaxies magnified by a factor of 13.

We search for signs of rotation in the subsystems of the Milky Way and M31 that are defined by their satellite galaxies, their globular cluster populations, and their BHB stars. A set of simple distribution functions are introduced to describe anisotropic and rotating stellar populations embedded in dark haloes of approximate Navarro-Frenk-White form. The BHB stars in the Milky Way halo exhibit a dichotomy between a prograde rotating, comparatively metal-rich component ([Fe/H] > -2) and a retrograde rotating, comparatively metal-poor ([Fe/H] < -2) component. The prograde metal-rich population may be associated with the accretion of a massive satellite (~ 10^9 solar masses). The metal-poor population may characterise the primordial stellar halo and the net retrograde rotation could then reflect an underestimate in our adopted local standard of rest circular velocity Theta_0. If Theta_0 is ~ 240 km/s then the metal-poor component has no rotation and there is a net prograde rotation signal of ~45 km/s in the metal-rich component. There is reasonable evidence that the Milky Way globular cluster and satellite galaxy systems are rotating with <v_phi> ~50 km/s and <v_phi> ~ 40 km/s respectively. Furthermore, a stronger signal is found for the satellite galaxies when the angular momentum vector of the satellites is inclined with respect to the normal of the disc. The dwarf spheroidal satellites of M31 exhibit prograde rotation relative to the M31 disc with <v_phi> ~ 40 km/s. We postulate that this group of dwarf spheroidals may share a common origin. We also find strong evidence for systemic rotation in the globular clusters of M31 particularly for the most metal-rich.

Using the Optimal Filter Technique applied to Sloan Digital Sky Survey photometry, we have found extended tails stretching about 1 degree (or several tens of half-light radii) from either side of the ultra-faint globular cluster Palomar 1. The tails contain roughly as many stars as does the cluster itself. Using deeper Hubble Space Telescope data, we see that the isophotes twist in a chacteristic S-shape on moving outwards from the cluster centre to the tails. We argue that the main mechanism forming the tails may be relaxation driven evaporation and that Pal 1 may have been accreted from a now disrupted dwarf galaxy ~500 Myr ago.

We present an AAOmega spectroscopic study of red giant stars in Bootes I, which is an ultra-faint dwarf galaxy, and Segue 1, suggested to be either an extremely low-luminosity dwarf galaxy or a star cluster. Our focus is quantifying the mean abundance and abundance dispersion in iron and carbon, and searching for distant radial-velocity members, in these systems. The primary conclusion of our investigation is that the spread of carbon abundance in both Bootes I and Segue 1 is large. For Bootes I, 4 of our 16 velocity members have [C/H] < ~-3.1, while 2 have [C/H] > ~-2.3, suggesting a range of Delta[C/H] ~ 0.8. For Segue 1 there exists a range Delta[C/H] ~ 1.0, including our discovery of a star with [Fe/H] = -3.5 and [C/Fe] = +2.3, which is a radial velocity member at a distance of 4 half-light radii from the system center. The accompanying ranges in iron abundance are Delta[Fe/H] ~ 1.6 for both Bootes I and Segue 1. For [Fe/H] < -3.0, the Galaxy's dwarf galaxy satellites exhibit a dependence of [C/Fe] on [Fe/H] which is very similar to that observed in its halo populations. We find [C/Fe] ~ 0.3 for stars in the dwarf systems that we believe are the counterpart of the Spite et al. (2005) ``unmixed'' giants of the Galactic halo and for which they report [C/Fe] ~ 0.2, and which presumably represents the natal relative abundance of carbon for material with [Fe/H] = -3.0 to -4.0. We confirm the correlation between luminosity and both mean metallicity and abundance dispersion in the Galaxy's dwarf satellites, which extends to at least as faint as Mv = -5. The very low mean metallicity of Segue 1, and the high carbon dispersion in Bootes I, consistent with inhomogeneous chemical evolution in near zero-abundance gas, suggest these ultra-faint systems could be surviving examples of the very first bound systems.

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 describe the application of non-negative matrix factorisation to generate compact reconstructions of quasar spectra from the Sloan Digital Sky Survey (SDSS), with particular reference to broad absorption line quasars (BALQSOs). BAL properties are measured for SiIV lambda1400, CIV lambda1550, AlIII lambda1860 and MgII lambda2800, resulting in a catalogue of 3547 BALQSOs. Two corrections, based on extensive testing of synthetic BALQSO spectra, are applied in order to estimate the intrinsic fraction of CIV BALQSOs. First, the probability of an observed BALQSO spectrum being identified as such by our algorithm is calculated as a function of redshift, signal-to-noise ratio and BAL properties. Second, the different completenesses of the SDSS target selection algorithm for BALQSOs and non-BAL quasars are quantified. Accounting for these selection effects the intrinsic CIV BALQSO fraction is 41+/-5 per cent. Our analysis of the selection effects allows us to measure the dependence of the intrinsic CIV BALQSO fraction on luminosity and redshift. We find a factor of 3.5+/-0.4 decrease in the intrinsic fraction from the highest redshifts, z~4.0, down to z~2.0. The redshift dependence implies that an orientation effect alone is not sufficient to explain the presence of BAL troughs in some but not all quasar spectra. Our results are consistent with the intrinsic BALQSO fraction having no strong luminosity dependence, although with 3-sigma limits on the rate of change of the intrinsic fraction with luminosity of -6.9 and 7.0 per cent dex^-1 we are unable to rule out such a dependence.

We confirm an eighth gravitational lens system in the CASSOWARY catalogue. Exploratory observations with the X-shooter spectrograph on the VLT show the system CSWA5 to consist of at least three images of a blue star-forming galaxy at z = 1.0686, lensed by an apparent foreground group of red galaxies one of which is at z = 0.3877. The lensed galaxy exhibits a rich spectrum with broad interstellar absorption lines and a wealth of nebular emission lines. Preliminary analysis of these features shows the galaxy to be young, with an age of 25-50 Myr. With a star-formation rate of approximately 20 solar masses/yr, the galaxy has already assembled a stellar mass of 3 x 10^9 solar masses and reached half-solar metallicity. Its blue spectral energy distribution and Balmer line ratios suggest negligible internal dust extinction. A more in-depth analysis of the properties of this system is currently hampered by the lack of a viable lensing model. However, it is already clear that CSWA5 shares many of its physical characteristics with the general population of UV-selected galaxies at redshifts z = 1-3, motivating further study of both the source and the foreground mass concentration responsible for the gravitational lensing.

The primary goal of this paper is to provide the evidence that can either prove or falsify the hypothesis that dark matter in the Galactic halo can clump into stellar-mass compact objects. If such objects existed, they would act as lenses to external sources in the Magellanic Clouds, giving rise to an observable effect of microlensing. We present the results of our search for such events, based on the data from the second phase of the OGLE survey (1996-2000) towards the SMC. The data set we used is comprised of 2.1 million monitored sources distributed over an area of 2.4 square degrees. We found only one microlensing event candidate, however its poor quality light curve limited our discussion on the exact distance to the lensing object. Given a single event, taking the blending (crowding of stars) into account for the detection efficiency simulations, and deriving the HST-corrected number of monitored stars, the microlensing optical depth is tau=(1.55+-1.55)10e-7. This result is consistent with the expected SMC self-lensing signal, with no need of introducing dark matter microlenses. Rejecting the unconvincing event leads to the upper limit on the fraction of dark matter in the form of MACHOs to f<20 per cent for deflectors' masses around 0.4 Msun and f<11 per cent for masses between 0.003 and 0.2 Msun (95 per cent confidence limit). Our result indicates that the Milky Way's dark matter is unlikely to be clumpy and form compact objects in the sub-solar-mass range.

We report the discovery of two new Milky Way satellites in the neighboring constellations of Pisces and Pegasus identified in data from the Sloan Digital Sky Survey. Pisces II, an ultra-faint dwarf galaxy lies at the distance of ~180 kpc, some 15 degrees away from the recently detected Pisces I. Segue 3, an ultra-faint star cluster lies at the distance of 16 kpc. We use deep follow-up imaging obtained with the 4-m Mayall telescope at Kitt Peak National Observatory to derive their structural parameters. Pisces II has a half-light radius of ~60 pc, while Segue 3 is twenty times smaller at only 3pc.

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 have used spectra obtained with X-shooter, the triple arm optical-infrared spectrograph recently commissioned on the Very Large Telescope of the European Southern Observatory, to confirm the gravitational lens nature of the CASSOWARY candidate CSWA 20. This system consists of a luminous red galaxy at redshift z(abs) = 0.741, with a very high velocity dispersion sigma(lens) = 500 km/s, which lenses a blue star-forming galaxy at z(em) = 1.433 into four images with mean separation of about 6 arcseconds. The source shares many of its properties with those of UV-selected galaxies z = 2-3: it is forming stars at a rate of 25 solar masses per year, has a metallicity of about 1/4 solar, and shows nebular emission from two components separated by 0.4 arcseconds in the image plane, possibly indicating a merger. It appears that foreground interstellar material within the galaxy has been evacuated from the sight-line along which we observe the starburst, giving an unextinguished view of its stars and H II regions. CSWA 20, with its massive lensing galaxy producing a high magnification of an intrinsically luminous background galaxy, is a promising target for future studies at a variety of wavelengths.

Our aim is to provide as clean and as complete a sample as possible of red giant branch stars that are members of the Hercules dSph galaxy. With this sample we explore the velocity dispersion and the metallicity of the system. Stromgren photometry and multi-fibre spectroscopy are combined to provide information about the evolutionary state of the stars (via the Stromgren c_1 index) and their radial velocities. Based on this information we have selected a clean sample of red giant branch stars, and show that foreground contamination by Milky Way dwarf stars can greatly distort the results. Our final sample consists of 28 red giant branch stars in the Hercules dSph galaxy. Based on these stars we find a mean photometric metallicity of -2.35 dex which is consistent with previous studies. We find evidence for an abundance spread. Using those stars for which we have determined radial velocities we find a systemic velocity of 45.2 km/s with a dispersion of 3.72 km/s, this is lower than values found in the literature. Furthermore we identify the horizontal branch and estimate the mean magnitude of the horizontal branch of the Hercules dSph galaxy to be V_0=21.17, which corresponds to a distance of 147 kpc. We have shown that a proper cleaning of the sample results in a smaller value for the velocity dispersion of the system. This has implications for galaxy properties derived from such velocity dispersions.

We construct a new sample of ~1700 solar neighbourhood halo subdwarfs from the Sloan Digital Sky Survey, selected using a reduced proper motion diagram. Radial velocities come from the SDSS spectra and proper motions from the light-motion curve catalogue of Bramich et al. (2008). Using a photometric parallax relation to estimate distances gives us the full phase-space coordinates. Typical velocity errors are in the range 30-50 km/s. This halo sample is one of the largest constructed to-date and the disc contamination is at a level of < 1 per cent. This enables us to calculate the halo velocity dispersion to excellent accuracy. We find that the velocity dispersion tensor is aligned in spherical polar coordinates and that (sigma_r, sigma_phi, sigma_theta) = (143 \pm 2, 82 \pm 2, 77 \pm 2) km/s. The stellar halo exhibits no net rotation, although the distribution of v_phi shows tentative evidence for asymmetry. The kinematics are consistent with a mildly flattened stellar density falling with distance like r^{-3.75}. Using the full phase-space coordinates, we look for signs of kinematic substructure in the stellar halo. We find evidence for four discrete overdensities localised in angular momentum and suggest that they may be possible accretion remnants. The most prominent is the solar neighbourhood stream previously identified by Helmi et al. (1999), but the remaining three are new. One of these overdensities is potentially associated with a group of four globular clusters (NGC5466, NGC6934, M2 and M13) and raises the possibility that these could have been accreted as part of a much larger progenitor.

We apply the optimal filter technique to Sloan Digital Sky Survey photometry around Segue 1 and find that the outer parts of the cluster are distorted. There is strong evidence for ~ 1 degree elongations of extra-tidal stars, extending both eastwards and southwestwards of the cluster. The extensions have similar differential Hess diagrams to Segue 1. A Kolmogorov-Smirnov test suggests a high probability that both come from the same parent distribution. The location of Segue 1 is close to crossings of the tidal wraps of the Sagittarius stream. By extracting blue horizontal branch stars from Sloan's spectral database, two kinematic features are isolated and identified with different wraps of the Sagittarius stream. We show that Segue 1 is moving with a velocity that is close to one of the wraps. At this location, we estimate that there are enough Sagittarius stars, indistinguishable from Segue 1 stars, to inflate the velocity dispersion and hence the mass-to-light ratio. All the available evidence is consistent with the interpretation that Segue 1 is a star cluster, originally from the Sagittarius galaxy, and now dissolving in the Milky Way.

We present an analysis of the substructure revealed by 407 RR Lyraes in Sloan Digital Sky Survey (SDSS) Stripe 82. Period estimates are determined to high accuracy using a string-length method. A subset of 178 RR Lyraes with spectrally derived metallicities are employed to derive metallicity-period-amplitude relations, which are then used to find metallicities and distances for the entire sample. The RR Lyraes lie between 5 and 115 kpc from the Galactic center. They are divided into subsets of 316 RRab types and 91 RRc types based on their period, colour and metallicity. The density distribution is not smooth, but dominated by clumps and substructure. Samples of 55 and 237 RR Lyraes associated with the Sagittarius Stream and the Hercules-Aquila Cloud respectively are identified. Hence, ~ 70 % of the RR Lyraes in Stripe 82 belong to known substructure. There is a sharp break in the density distribution at Galactocentric radii of 40 kpc, reflecting the fact that the dominant substructure in Stripe 82 - the Hercules-Aquila Cloud and the Sagittarius Stream - lies within 40 kpc. In fact, almost 60 % of all the RR Lyraes in Stripe 82 are associated with the Hercules-Aquila Cloud alone, which emphasises its pre-eminence. Additionally, evidence of a new and distant substructure - the Pisces Overdensity - is found, consisting of 28 faint RR Lyraes centered on Galactic coordinates (80 deg, -55 deg) and with distances of ~ 80 kpc. The total stellar mass in the Pisces Overdensity is ~10000 solar masses and its metallicity is [Fe/H] ~ -1.5.

We present the results from the OGLE-II survey (1996-2000) towards the Large Magellanic Cloud (LMC), which has the aim of detecting the microlensing phenomena caused by dark matter compact objects in the Galactic Halo (Machos). We use high resolution HST images of the OGLE fields and derive the correction for the number of monitored stars in each field. This also yield blending distributions which we use in 'catalogue level' Monte Carlo simulations of the microlensing events in order to calculate the detection efficiency of the events. We detect two candidates for microlensing events in the All Stars Sample, which translates into an optical depth of 0.43+-0.33x 10e-7. If both events were due to Macho the fraction of mass of compact dark matter objects in the Galactic halo would be 8+-6 per cent. This optical depth, however, along with the characteristics of the events, seems to be consistent with the self-lensing scenario, i.e., self-lensing alone is sufficient to explain the observed microlensing signal. Our results indicate a non-detection of Machos lensing towards the LMC with an upper limit on their abundance in the Galactic halo of 19 per cent for M=0.4 Msun and 10 per cent for masses between 0.01 and 0.2 Msun.

We present radial velocity measurements of four wide halo binary candidates from the sample in Chaname & Gould (2004; CG04) which, to date, is the only sample containing a large number of such candidates. The four candidates that we have observed have projected separations >0.1 pc, and include the two widest binaries from the sample, with separations of 0.45 and 1.1 pc. We confirm that three of the four CG04 candidates are genuine, including the one with the largest separation. The fourth candidate, however, is spurious at the 5-sigma level. In the light of these measurements we re-examine the implications for MACHO models of the Galactic halo. Our analysis casts doubt on what MACHO constraints can be drawn from the existing sample of wide halo binaries.

We announce the discovery of a new Milky Way satellite Segue 2 found in the data of the Sloan Extension for Galactic Understanding and Exploration (SEGUE). We followed this up with deeper imaging and spectroscopy on the Multiple Mirror Telescope. From this, we derive a luminosity of M_v = -2.5, a half-light radius of 34 pc and a systemic velocity of -40$ km/s. Our MMT data also provides evidence for a stream around Segue 2 at a similar heliocentric velocity, and the SEGUE data show that it is also present in neighboring fields. We resolve the velocity dispersion of Segue 2 as 3.4 km/s and the possible stream as about 7 km/s. This object shows points of comparison with other recent discoveries, Segue 1, Boo II and Coma. We speculate that all four objects may be representatives of a population of satellites of satellites -- survivors of accretion events that destroyed their larger but less dense parents. They are likely to have formed at redshifts z > 10 and are good candidates for fossils of the reionization epoch.

We present medium-resolution spectra of 16 radial velocity red-giant members of the low-luminosity Bootes I dwarf spheroidal (dSph) galaxy, that have sufficient S/N for abundance determination, based on the strength of the Ca II K line. Assuming [Ca/Fe] ~ +0.3, the abundance range in the sample is Delta [Fe/H] ~ 1.7 dex, with one star having [Fe/H] = -3.4. The dispersion is sigma([Fe/H]) = 0.45 +/- 0.08 -- similar to those of the Galaxy's more luminous dSph systems and Omega Centauri. This suggests that the large mass (greater than approximately 10 million solar masses) normally assumed to foster self-enrichment and the production of chemical abundance spreads was provided by the non-baryonic material in Bootes I.

We present a spectroscopic study of Leo V, a recently discovered satellite of the Milky Way (MW). From stellar spectra obtained with the MMT/Hectochelle spectrograph we identify seven likely members of Leo V. Five cluster near the Leo V center (R < 3 arcmin) and have velocity dispersion 2.4_{-1.4}^{+2.4} km/s. The other two likely members lie near each other but far from the center (R~13 arcmin ~ 700 pc) and inflate the global velocity dispersion to 3.7_{-1.4}^{+2.3} km/s. Assuming the five central members are bound, we obtain a dynamical mass of M=3.3_{-2.5}^{+9.1} x 10^5M_{sun} (M/L_V=75_{-58}^{+230}[M/L_V]_{sun}). From the stacked spectrum of the five central members we estimate a mean metallicity of [Fe/H]=-2.0\pm 0.2 dex. Thus with respect to dwarf spheroidals of similar luminosity, Leo V is slightly less massive and slightly more metal-rich. Since we resolve the central velocity dispersion only marginally, we do not rule out the possibility that Leo V is a diffuse star cluster devoid of dark matter. The wide separation of its two outer members implies Leo V is losing mass; however, its large distance (D ~ 180 kpc) is difficult to reconcile with MW tidal stripping unless the orbit is very radial.

Self-Organizing Map (SOM) is a promising tool for exploring large multi-dimensional data sets. It is quick and convenient to train in an unsupervised fashion and, as an outcome, it produces natural clusters of data patterns. An example of application of SOM to the new OGLE-III data set is presented along with some preliminary results. Once tested on OGLE data, the SOM technique will also be implemented within the Gaia mission's photometry and spectrometry analysis, in particular, in so-called classification-based Science Alerts. SOM will be used as a basis of this system as the changes in brightness and spectral behaviour of a star can be easily and quickly traced on a map trained in advance with simulated and/or real data from other surveys.

Modern spectroscopic databases provide a wealth of information about the physical processes and environments associated with astrophysical populations. Techniques such as blind source separation (BSS), in which sets of spectra are decomposed into a number of components, offer the prospect of identifying the signatures of the underlying physical emission processes. Principle Component Analysis (PCA) has been applied with some success but is severely limited by the inherent orthogonality restriction that the components must satisfy. Non-negative matrix factorisation (NMF) is a relatively new BSS technique that incorporates a non-negativity constraint on its components. In this respect, the resulting components may more closely reflect the physical emission signatures than is the case using PCA. We discuss some of the considerations that must be made when applying NMF and, through its application to the quasar spectra in the Sloan Digital Sky Survey (SDSS) DR6, we show that NMF is a fast method for generating compact and accurate reconstructions of the spectra. The ability to reconstruct spectra accurately has numerous astrophysical applications. Combined with improved SDSS redshifts, we apply NMF to the problem of defining robust continua for quasars that exhibit strong broad absorption line (BAL) systems. The resulting catalogue of SDSS DR6 BAL quasars will be the largest available. Importantly, the NMF approach allows quantitative error estimates to be derived for the Balnicity Indices as a function of key astrophysical and observational parameters, such as the quasar redshifts and the signal-to-noise ratio of the spectra.

We present discovery images, together with follow-up imaging and spectroscopy, of two large separation gravitational lenses found by our survey for wide arcs (the CASSOWARY). The survey exploits the multicolor photometry of the Sloan Digital Sky Survey to find multiple blue components around red galaxies. CASSOWARY~2 (or "the Cheshire Cat") is composed of two massive early-type galaxies at z = 0.426 and 0.432 respectively lensing two background sources, the first a star-forming galaxy at z = 0.97 and the second a high redshift galaxy (z> 1.4). There are at least three images of the former source and probably four or more of the latter, arranged in two giant arcs. The mass enclosed within the larger arc of radius 11 arcsecs is about 33 x 10^{12} solar masses. CASSOWARY~3 comprises an arc of three bright images of a z = 0.725 source, lensed by a foreground elliptical at z = 0.274. The radius of the arc is about 4 arcsecs and the enclosed mass is 2.5 x 10^{12} solar masses. Together with earlier discoveries like the Cosmic Horseshoe and the 8 O'Clock Arc, these new systems, with separations intermediate between the arcsecond separation lenses of typical strong galaxy lensing and the larger separation cluster lenses, probe the very high end of the galaxy mass function.

We report on the abundance analysis of two red giants in the faint Hercules dwarf spheroidal (dSph) galaxy. These stars show a remarkable deficiency in the neutron-capture elements, while the hydrostatic alpha-elements (O, Mg) are strongly enhanced. Our data indicate [Ba/Fe] and [Mg/Fe] abundance ratios of <-2 dex and ~+0.8 dex, respectively, with essentially no detection of other n-capture elements. In contrast to the only other dSph star with similar abundance patterns, Dra 119, which has a very low metallicity at [Fe/H]=-2.95 dex, our objects, at [Fe/H]~-2.0 dex, are only moderately metal poor. The measured ratio of hydrostatic/explosive alpha-elements indicates that high-mass (~35 M_sun) Type II supernovae progenitors are the main, if not only, contributors to the enrichment of this galaxy. This suggests that star formation and chemical enrichment in the ultrafaint dSphs proceeds stochastically and inhomogeneously on small scales, or that the IMF was strongly skewed to high mass stars. The neutron capture deficiencies and the [Co/Fe] and [Cr/Fe] abundance ratios in our stars are similar to those in the extremely low metallicity Galactic halo. This suggests that either our stars are composed mainly of the ejecta from the first, massive, population III stars (but at moderately high [Fe/H]), or that SN ejecta in the Hercules galaxy were diluted with ~30 times less hydrogen than typical for extreme metal-poor stars.

We present a new suite of photometric and spectroscopic data for the faint Bootes II dwarf spheroidal galaxy candidate. Our deep photometry, obtained with the INT/WFC, suggests a distance of 46 kpc and a small half-light radius of 4.0 arcmin (56 pc), consistent with previous estimates. Follow-up spectroscopy obtained with the Gemini/GMOS instrument yielded radial velocities and metallicities. While the majority of our targets covers a broad range in velocities and metallicities, we find five stars which share very similar velocities and metallicities and which are all compatible with the colors and magnitudes of the galaxy's likely red giant branch. We interpret these as a spectroscopic detection of the Bootes II system. These stars have a mean velocity of -117 km/s, a velocity dispersion of (10.5+-7.4) km/s and a mean [Fe/H] of -1.79 dex, with a dispersion of 0.14 dex. At this metallicity, Boo II is not consistent with the stellar-mass-metallicity relation for the more luminous dwarf galaxies. Coupled with our distance estimate, its high negative systemic velocity rules out any physical connection with its projected neighbor, the Bootes I dwarf spheroidal, which has a velocity of ~+100 km/s. The velocity and distance of Bootes II coincide with those of the leading arm of Sagittarius, which passes through this region of the sky, so that it is possible that Bootes II may be a stellar system associated with the Sagittarius stream. Finally, we note that the properties of Bootes II are consistent with being the surviving remnant of a previously larger and more luminous dSph galaxy.

We report the discovery of a new Milky Way dwarf spheroidal galaxy in the constellation of Leo identified in data from the Sloan Digital Sky Survey. Leo V lies at a distance of about 180 kpc, and is separated by about 3 degrees from another recent discovery, Leo IV. We present follow-up imaging from the Isaac Newton Telescope and spectroscopy from the Hectochelle fiber spectrograph at the Multiple Mirror Telescope. Leo V's heliocentric velocity is 173.4 km/s, which is offset by about 40 km/s from that of Leo IV. A simple interpretation of the kinematic data is that both objects may lie on the same stream, though the implied orbit is only modestly eccentric (e = 0.2)

We present new spectroscopic data for twenty six stars in the recently-discovered Canes Venatici I (CVnI) dwarf spheroidal galaxy. We use these data to investigate the recent claim of the presence of two dynamically inconsistent stellar populations in this system (Ibata et al., 2006). We do not find evidence for kinematically distinct populations in our sample and we are able to obtain a mass estimate for CVnI that is consistent with all available data, including previously published data. We discuss possible differences between our sample and the earlier data set and study the general detectability of sub-populations in small kinematic samples. We conclude that in the absence of supporting observational evidence (for example, metallicity gradients), sub-populations in small kinematic samples (typically fewer than 100 stars) should be treated with extreme caution, as their detection depends on multiple parameters and rarely produces a signal at the 3sigma confidence level. It is therefore essential to determine explicitly the statistical significance of any suggested sub-population.

We model the extremely massive and luminous lens galaxy in the Cosmic Horseshoe Einstein ring system, recently discovered in the Sloan Digital Sky Survey. We use the semi-linear method of Warren & Dye (2003), which pixelises the source surface brightness distribution, to invert the Einstein ring for sets of parameterised lens models. Here, the method is refined by exploiting Bayesian inference to optimise adaptive pixelisation of the source plane and to choose between three differently parameterised models: a singular isothermal ellipsoid, a power law model and a NFW profile. The most probable lens model is the power law with a volume mass density that scales as r^(-1.96+/-0.02) and an axis ratio of ~0.8. The mass within the Einstein ring (i.e., within a cylinder with projected distance of ~30 kpc from the centre of the lens galaxy) is (5.02+/-0.09)*10^12 M_solar, and the mass-to-light ratio is ~30. Even though the lens lies in a group of galaxies, the preferred value of the external shear is almost zero. This makes the Cosmic Horseshoe unique amongst large separation lenses, as almost all the deflection comes from a single, very massive galaxy with little boost from the environment.

We quantify the detectability of stellar Milky Way satellites in the Sloan Digital Sky Survey (SDSS) Data Release 5. We show that the effective search volumes for the recently discovered SDSS--satellites depend strongly on their luminosity, with their maximum distance, $D_{max}$, substantially smaller than the Milky Way halo's virial radius. Calculating the maximum accessible volume, $V_{max}$, for all faint detected satellites, allows the calculation of the luminosity function for Milky Way satellite galaxies, accounting quantitatively for their detectability. We find that the number density of satellite galaxies continues to rise towards low luminosities, but may flatten at $M_V \sim -5$; within the uncertainties, the luminosity function can be described by a single power law $dN/dM_{V}= 10 \times 10^{0.1 (M_V+5)}$, spanning luminosities from $M_V=-2$ all the way to the luminosity of the Large Magellanic Cloud. Comparing these results to several semi-analytic galaxy formation models, we find that their predictions differ significantly from the data: either the shape of the luminosity function, or the surface brightness distributions of the models, do not match.

There has been a vast recent improvement in photometric and kinematic data for star clusters, Ultra Compact dwarfs, galactic nuclei, and local dSph galaxies, with Subaru contributing substantially to the photometric studies in particular. These data show that there is a bimodal distribution in half-light radii, with stable star clusters always being smaller than 35pc, while stable galaxies are always larger than 120pc. We extend the previously known observational relationships and interpret them in terms of a more fundamental pair of intrinsic properties of dark matter itself: dark matter forms cored mass distributions, with a core scale length of greater than about 100pc, and always has a maximum central mass density with a narrow range. The dark matter in dSph galaxies appears to be clustered such that there is a mean volume mass density within the stellar distribution which has the very low value of about 0.1$M_{\odot}$ pc$^{-3}$. None of the dSphs displays kinematics which require the presence of an inner cusp, while in two dSphs there is evidence that the density profile is shallow (cored) in the inner regions. The maximum central dark matter density derived is model dependent, but is likely to have a mean value (averaged over a volume of radius 10pc) of about 0.1$M_{\odot}$ pc$^{-3}$, which is 5GeV/c$^2$cm$^{-3}$). Galaxies are embedded in dark matter halos with these properties; smaller systems containing dark matter are not observed.

We have used data from the Sloan Digital Sky Survey (SDSS) Data Release 5 to explore the overall structure and substructure of the stellar halo of the Milky Way using about 4 million color-selected main sequence turn-off stars. We fit oblate and triaxial broken power-law models to the data, and found a `best-fit' oblateness of the stellar halo 0.5<c/a<0.8, and halo stellar masses between Galactocentric radii of 1 and 40kpc of (3.7+/-1.2)x10^8 M_sun. The density profile of the stellar halo is approximately r^{-3}; it is possible that the power law slope is shallower inside 20kpc and steeper outside that radius. Yet, we found that all smooth and symmetric models were very poor fits to the distribution of stellar halo stars because the data exhibit a great deal of spatial substructure. We quantified deviations from a smooth oblate/triaxial model using the RMS of the data around the model profile on scales >~100pc, after accounting for the (known) contribution of Poisson uncertainties. The fractional RMS deviation of the actual stellar distribution from any smooth, parameterized halo model is >~40%: hence, the stellar halo is highly structured. We compared the observations with simulations of galactic stellar halos formed entirely from the accretion of satellites in a cosmological context by analysing the simulations in the same way as the data. While the masses, overall profiles, and degree of substructure in the simulated stellar halos show considerable scatter, the properties and degree of substructure in the Milky Way's halo match well the properties of a `typical' stellar halo built exclusively out of the debris from disrupted satellite galaxies. Our results therefore point towards a picture in which an important fraction of the Milky Way's stellar halo has been accreted from satellite galaxies.

We present a new public archive of light-motion curves in Sloan Digital Sky Survey (SDSS) Stripe 82, covering 99 deg in right ascension from RA = 20.7 h to 3.3 h and spanning 2.52 deg in declination from Dec = -1.26 to 1.26 deg, for a total sky area of ~249 sq deg. Stripe 82 has been repeatedly monitored in the u, g, r, i and z bands over a seven-year baseline. Objects are cross-matched between runs, taking into account the effects of any proper motion. The resulting catalogue contains almost 4 million light-motion curves of stellar objects and galaxies. The photometry are recalibrated to correct for varying photometric zeropoints, achieving ~20 mmag and ~30 mmag root-mean-square (RMS) accuracy down to 18 mag in the g, r, i and z bands for point sources and extended sources, respectively. The astrometry are recalibrated to correct for inherent systematic errors in the SDSS astrometric solutions, achieving ~32 mas and ~35 mas RMS accuracy down to 18 mag for point sources and extended sources, respectively. For each light-motion curve, 229 photometric and astrometric quantities are derived and stored in a higher-level catalogue. On the photometric side, these include mean exponential and PSF magnitudes along with uncertainties, RMS scatter, chi^2 per degree of freedom, various magnitude distribution percentiles, object type (stellar or galaxy), and eclipse, Stetson and Vidrih variability indices. On the astrometric side, these quantities include mean positions, proper motions as well as their uncertainties and chi^2 per degree of freedom. The here presented light-motion curve catalogue is complete down to r~21.5 and is at present the deepest large-area photometric and astrometric variability catalogue available.

We report on 2MASS J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Database. Additional photometry from the Sloan Digital Sky Survey yields an 8-passband lightcurve, from which we derive an orbital period of 2.6390157 +/- 0.0000016 days. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses (M_1 = 0.66 +/- 0.03 M_sun; M_2 = 0.62 +/- 0.03 M_sun) and radii (R_1 = 0.64 +/- 0.08 R_sun; R_2 = 0.61 +/- 0.09 R_sun) of the components, which are consistent with empirical mass-radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the lightcurves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of H-alpha emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.

We have tested the application to Sloan Digital Sky Survey data of the software package MATCH, which fits color-magnitude diagrams (CMDs) to estimate stellar population parameters and distances. These tests on a set of six globular clusters show that these techniques recover their known properties. New ways of using the CMD-fitting software enable us to deal with an extended distribution of stars along the line-of-sight, to constrain the overall properties of sparsely populated objects, and to detect the presence of stellar overdensities in wide-area surveys. We then also apply MATCH to CMDs for twelve recently discovered Milky Way satellites to derive in a uniform fashion their distances, ages and metallicities. While the majority of them appear consistent with a single stellar population, CVn I, UMa II, and Leo T exhibit (from SDSS data alone) a more complex history with multiple epochs of star formation.

We investigate a wide range of possible evolutionary histories for the recently discovered Bootes dwarf spheroidal galaxy, a Milky Way satellite. By means of N-body simulations we follow the evolution of possible progenitor galaxies of Bootes for a variety of orbits in the gravitational potential of the Milky Way. The progenitors considered cover the range from dark-matter-free star clusters to massive, dark-matter dominated outcomes of cosmological simulations. For each type of progenitor and orbit we compare the observable properties of the remnant after 10 Gyr with those of Bootes observed today. Our study suggests that the progenitor of Bootes must have been, and remains now, dark matter dominated. In general our models are unable to reproduce the observed high velocity dispersion in Bootes without dark matter. Our models do not support time-dependent tidal effects as a mechanism able to inflate significantly the internal velocity dispersion. As none of our initially spherical models is able to reproduce the elongation of Bootes, our results suggest that the progenitor of Bootes may have had some intrinsic flattening. Although the focus of the present paper is the Bootes dwarf spheroidal, these models may be of general relevance to understanding the structure, stability and dark matter content of all dwarf spheroidal galaxies.

We present Giant Meterwave Radio Telescope (GMRT) and Westerbork ynthesis Radio Telescope (WSRT) observations of the recently discovered Local Group dwarf galaxy, Leo T. The peak HI column density is measured to be 7x10^20 cm^-2, and the total HI mass is 2.8Xx10^5 Msun, based on a distance of 420 kpc. Leo T has both cold (~ 500 K) and warm (~ 6000 K) HI at its core, with a global velocity dispersion of 6.9 km/s, from which we derive a dynamical mass within the HI radius of 3.3x10^6 Msun, and a mass-to-light ratio of greater than 50. We calculate the Jeans mass from the radial profiles of the HI column density and velocity dispersion, and predict that the gas should be globally stable against star formation. This finding is inconsistent with the half light radius of Leo T, which extends to 170 pc, and indicates that local conditions must determine where star formation takes place. Leo T is not only the lowest luminosity galaxy with on-going star formation discovered to date, it is also the most dark matter dominated, gas-rich dwarf in the Local Group.

A 2.5 x 100 degree region along the celestial equator (Stripe 82) has been imaged repeatedly from 1998 to 2005 by the Sloan Digital Sky Survey. A new catalogue of ~4 million light-motion curves, together with over 200 derived statistical quantities, for objects in Stripe 82 brighter than r~21.5 has been constructed by combining these data by Bramich et al. (2007). This catalogue is at present the deepest catalogue of its kind. Extracting the ~130000 objects with highest signal-to-noise ratio proper motions, we build a reduced proper motion diagram to illustrate the scientific promise of the catalogue. In this diagram disk and halo subdwarfs are well-separated from the cool white dwarf sequence. Our sample of 1049 cool white dwarf candidates includes at least 8 and possibly 21 new ultracool H-rich white dwarfs (T_eff < 4000K) and one new ultracool He-rich white dwarf candidate identified from their SDSS optical and UKIDSS infrared photometry. At least 10 new halo white dwarfs are also identified from their kinematics.

We report the discovery of an almost complete Einstein ring of diameter 10" in Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5). Spectroscopic data from the 6m telescope of the Special Astrophysical Observatory reveals that the deflecting galaxy has a line-of-sight velocity dispersion in excess of 400 km/s and a redshift of 0.444, whilst the source is a star-forming galaxy with a redshift of 2.379. From its color and luminosity, we conclude that the lens is an exceptionally massive Luminous Red Galaxy (LRG) with a mass within the Einstein radius of 5 x 10^12 solar masses. This remarkable system provides a laboratory for probing the dark matter distribution in LRGs at distances out to 3 effective radii, and studying the properties of high redshift star-forming galaxies.

We report the discovery of two extremely low luminosity globular clusters in the Milky Way Halo. These objects were detected in the Sloan Digital Sky Survey Data Release 5 and confirmed with deeper imaging at the Calar Alto Observatory. The clusters, Koposov 1 and Koposov 2, are located at $\sim 40-50$ kpc and appear to have old stellar populations and luminosities of only $M_V \sim -1$ mag. Their observed sizes of $\sim 3$ pc are well within the expected tidal limit of $\sim$10 pc at that distance. Together with Palomar 1, AM 4 and Whiting 1, these new clusters are the lowest luminosity globulars orbiting the Milky Way, with Koposov 2 the most extreme. Koposov 1 appears to lie close to distant branch of the Sagittarius stream. The half-mass relaxation times of Koposov 1 and 2 are only $\sim 70$ and $\sim 55$ Myr respectively (2 orders of magnitude shorter than the age of the stellar populations), so it would seem that they have undergone drastic mass segregation. Since they do not appear to be very concentrated, their evaporation timescales may be as low as $\sim 0.1 t_{\rm Hubble}$. These discoveries show that the structural parameter space of globular clusters in the Milky Way halo is not yet fully explored. They also add, through their short remaining survival times, significant direct evidence for a once much larger population of globular clusters.

The study of substructure in the stellar halo of the Milky Way has made a lot of progress in recent years, especially with the advent of surveys like the Sloan Digital Sky Survey. Here, we study the newly discovered tidal tails of the Galactic globular cluster NGC 5466. By means of numerical simulations, we reproduce the shape, direction and surface density of the tidal tails, as well as the structural and kinematical properties of the present-day NGC 5466. Although its tails are very extended in SDSS data (> 45 degrees), NGC 5466 is only losing mass slowly at the present epoch and so can survive for probably a further Hubble time. The effects of tides at perigalacticon and disc crossing are the dominant causes of the slow dissolution of NGC 5466, accounting for about 60 % of the mass loss over the course of its evolution. The morphology of the tails provides a constraint on the proper motion -- the observationally determined proper motion has to be refined (within the stated error margins) to match the location of the tidal tails.

We present evidence for a substantial overdensity of stars in the direction of the constellations of Hercules and Aquila. The Cloud is centered at a Galactic longitude of about 40 degrees and extends above and below the Galactic plane by at least 50 degrees. Given its off-centeredness and height, it is unlikely that the Hercules-Aquila Cloud is related to the bulge or thick disk. More likely, this is a new structural component of the Galaxy that passes through the disk. The Cloud stretches about 80 degrees in longitude. Its heliocentric distance lies between 10 and 20 kpc so that the extent of the Cloud in projection is roughly 20 kpc by 15 kpc. It has an absolute magnitude of -13 and its stellar population appears to be comparable to, but somewhat more metal-rich than, M92.

In this Letter, we announce the discovery of a new dwarf galaxy, Leo T, in the Local Group. It was found as a stellar overdensity in the Sloan Digital Sky Survey Data Release 5 (SDSS DR5). The color-magnitude diagram of Leo T shows two well-defined features, which we interpret as a red giant branch and a sequence of young, massive stars. As judged from fits to the color-magnitude diagram, it lies at a distance of about 420 kpc and has an intermediate-age stellar population with a metallicity of [Fe/H]= -1.6, together with a young population of blue stars of age of 200 Myr. There is a compact cloud of neutral hydrogen with mass roughly 10^5 solar masses and radial velocity 35 km/s coincident with the object visible in the HIPASS channel maps. Leo T is the smallest, lowest luminosity galaxy found to date with recent star-formation. It appears to be a transition object similar to, but much lower luminosity than, the Phoenix dwarf.

Prominent in the `Field of Streams' -- the Sloan Digital Sky Survey map of substructure in the Galactic halo -- is an `Orphan Stream' without obvious progenitor. In this numerical study, we show a possible connection between the newly found dwarf satellite Ursa Major II (UMa II) and the Orphan Stream. We provide numerical simulations of the disruption of UMa II that match the observational data on the position, distance and morphology of the Orphan Stream. We predict the radial velocity of UMa II as -100 km/s as well as the existence of strong velocity gradients along the Orphan Stream. The velocity dispersion of UMa II is expected to be high, though this can be caused both by a high dark matter content or by the presence of unbound stars in a disrupted remnant. However, the existence of a gradient in the mean radial velocity across UMa II provides a clear-cut distinction between these possibilities. The simulations support the idea that some of the anomalous, young halo globular clusters like Palomar 1 or Arp 2 or Ruprecht 106 may be physically associated with the Orphan Stream.

We use Sloan Digital Sky Survey Data Release 5 photometry and spectroscopy to study a tidal stream that extends over 50 degrees in the North Galactic Cap. From the analysis of the path of the stream and the colors and magnitudes of its stars, the stream is about 20 kpc away at its nearest detection (the celestial equator). We detect a distance gradient -- the stream is farther away from us at higher declination. The contents of the stream are made up from a predominantly old and metal-poor population that is similar to the globular clusters M13 and M92. The integrated absolute magnitude of the stream stars is estimated to be M_r = -7.5. There istentative evidence for a velocity signature, with the stream moving at -40 km/s at low declinations and +100 km/s at high declinations. The stream lies on the same great circle as Complex A, a roughly linear association of HI high velocity clouds stretching over 30 degrees on the sky, and as Ursa Major II, a recently discovered dwarf spheroidal galaxy. Lying close to the same great circle are a number of anomalous, young and metal-poor globular clusters, including Palomar 1 and Ruprecht 106.

Using neural networks, Belokurov, Evans & Le Du (2003, 2004) showed that 7 out of the 29 microlensing candidates towards the Large Magellanic Cloud (LMC) of the MACHO collaboration are consistent with blended microlensing and added Gaussian noise. They then estimated the microlensing optical depth to the LMC to be between 0.3 x 10^{-7} and 0.5 x 10^{-7}, lower by about a factor of two than that found by the MACHO collaboration. There has been a recent independent claim of a low optical depth to the LMC by the EROS collaboration, who find 0.15 x 10^{-7} from 3 candidates (Tisserand et al 2006). Griest & Thomas (2005) have contested our calculations. Unfortunately, their paper contains a number of scientific misrepresentations of our work. We stand by our application of neural networks to microlensing searches and believe it to be a technique of great promise. Rather, the main cause of the disparity between Griest & Thomas and Belokurov et al. lies in the very different datasets through which these investigators look for microlensing events. We don't exclude the possibility that some of our non-microlensing designations may change with access to clean data. Whilst not everything is understood about the microlensing datasets towards the LMC, the latest downward revisions of the optical depth means that Griest & Thomas' microlensing puzzle is a roughly 1 sigma effect. Efficiency calculations can correct for the effects of false negatives, but they cannot correct for the effects of false positives (variable stars that are mistaken for microlensing). Therefore, the best strategy in a microlensing experiment is to eschew a decision boundary altogether. Rather, each lightcurve should be assigned a probability and the microlensing rate calculated by summing over the probabilities of all such lightcurves.

We present five new satellites of the Milky Way discovered in Sloan Digital Sky Survey (SDSS) imaging data, four of which were followed-up with either the Subaru or the Isaac Newton Telescopes. They include four probable new dwarf galaxies -- one each in the constellations of Coma Berenices, Canes Venatici, Leo and Hercules -- together with one unusually extended globular cluster, Segue 1. We provide distances, absolute magnitudes, half-light radii and color-magnitude diagrams for all five satellites. The morphological features of the color-magnitude diagrams are generally well described by the ridge line of the old, metal-poor globular cluster M92. In the last two years, a total of ten new Milky Way satellites with effective surface brightness mu_v >~ 28 mag/sq. arcsec have been discovered in SDSS data. They are less luminous, more irregular and appear to be more metal-poor than the previously-known nine Milky Way dwarf spheroidals. The relationship between these objects and other populations is discussed. We note that there is a paucity of objects with half-light radii between ~40 pc and ~ 100 pc. We conjecture that this may represent the division between star clusters and dwarf galaxies.

In this Letter, we study a localized stellar overdensity in the constellation of Ursa Major, first identified in Sloan Digital Sky Survey (SDSS) data and subsequently followed up with Subaru imaging. Its color-magnitude diagram (CMD) shows a well-defined sub-giant branch, main sequence and turn-off, from which we estimate a distance of ~30 kpc and a projected size of ~250 x 125 pc. The CMD suggests a composite population with some range in metallicity and/or age. Based on its extent and stellar population, we argue that this is a previously unknown satellite galaxy of the Milky Way, hereby named Ursa Major II (UMa II) after its constellation. Using SDSS data, we find an absolute magnitude of M_V \~ -3.8, which would make it the faintest known satellite galaxy. UMa II's isophotes are irregular and distorted with evidence for multiple concentrations; this suggests that the satellite is in the process of disruption.

In this Letter, we announce the discovery of a new satellite of the Milky Way in the constellation of Bootes at a distance of 60 kpc. It was found in a systematic search for stellar overdensities in the North Galactic Cap using Sloan Digital Sky Survey Data Release 5 (SDSS DR5). The color-magnitude diagram shows a well-defined turn-off, red giant branch, and extended horizontal branch. Its absolute magnitude is -5.8, which makes it one of the faintest galaxies known. The half-light radius is 220 pc. The isodensity contours are elongated and have an irregular shape, suggesting that Boo may be a disrupted dwarf spheroidal galaxy.

The latest Sloan Digital Sky Survey data reveal a prominent bifurcation in the distribution of debris of the Sagittarius dwarf spheroidal (Sgr) beginning at a right ascension of roughly 190 degrees. Two branches of the stream (A and B) persist at roughly the same heliocentric distance over at least 50 degrees of arc. There is also evidence for a more distant structure (C) well behind the A branch. This paper provides the first explanation for the bifurcation. It is caused by the projection of the young leading (A) and old trailing (B) tidal arms of the Sgr, whilst the old leading arm (C) lies well behind A. This explanation is only possible if the halo is close to spherical, as the angular difference between the branches is a measure of the precession of the orbital plane.

We use Sloan Digital Sky Survey (SDSS) Data Release 5 (DR5) u,g,r,i,z photometry to study Milky Way halo substructure in the area around the North Galactic Cap. A simple color cut (g-r < 0.4) reveals the tidal stream of the Sagittarius dwarf spheroidal, as well as a number of other stellar structures in the field. Two branches (A and B) of the Sagittarius stream are clearly visible in an RGB-composite image created from 3 magnitude slices, and there is also evidence for a still more distant wrap behind the A branch. A comparison of these data with numerical models suggests that the shape of the Galactic dark halo is close to spherical.

In this Letter, we announce the discovery of a new dwarf satellite of the Milky Way, located in the constellation Canes Venatici. It was found as a stellar overdensity in the North Galactic Cap using Sloan Digital Sky Survey Data Release 5 (SDSS DR5). The satellite's color-magnitude diagram shows a well-defined red giant branch, as well as a horizontal branch. As judged from the tip of the red giant branch, it lies at a distance of ~220 kpc. Based on the SDSS data, we estimate an absolute magnitude of Mv ~ -7.9, a central surface brightness of mu_0,V ~ 28 mag arcsecond^-2, and a half-light radius of \~ 8.5' (~ 550 pc at the measured distance). The outer regions of Canes Venatici appear extended and distorted. The discovery of such a faint galaxy in proximity to the Milky Way strongly suggests that more such objects remain to be found.

We report the discovery of tidal tails around the high-latitude Galactic globular cluster NGC 5466 in Sloan Digital Sky Survey (SDSS) data. Neural networks are used to reconstruct the probability distribution of cluster stars in u,g,r,i and z space. The tails are clearly visible, once extra-galactic contaminants and field stars have been eliminated. They extend roughly 4 degrees on the sky, corresponding to about 1 kpc in projected length. The orientation of the tails is in good agreement with the cluster's Galactic orbit, as judged from the proper motion data.

The POINT-AGAPE collaboration is carrying out a search for gravitational microlensing toward M31 to reveal galactic dark matter in the form of MACHOs (Massive Astrophysical Compact Halo Objects) in the halos of the Milky Way and M31. A high-threshold analysis of 3 years of data yields 6 bright, short--duration microlensing events, which are confronted to a simulation of the observations and the analysis. The observed signal is much larger than expected from self lensing alone and we conclude, at the 95% confidence level, that at least 20% of the halo mass in the direction of M31 must be in the form of MACHOs if their average mass lies in the range 0.5-1 M$_\odot$. This lower bound drops to 8% for MACHOs with masses $\sim 0.01$ M$_\odot$. In addition, we discuss a likely binary microlensing candidate with caustic crossing. Its location, some 32' away from the centre of M31, supports our conclusion that we are detecting a MACHO signal in the direction of M31.

We report the results of our statistical study of the POINT-AGAPE CNe catalogue. The first Sloan $r'$ and $i'$ analysis of the MMRD relationship in M31 is performed using the resulting POINT-AGAPE CN catalogue. Good fits are produced to the MMRD in the two filters. However, we are unable to verify the $t_{15}$ relationship for either Sloan filter. The subsequent analysis of our automated CN detection pipeline has provided us with the most thorough knowledge of the completeness of a CN survey to-date. In addition, the large field of view of the survey has permitted us to probe the outburst rate well into the galactic disk, unlike previous CCD imaging surveys. Using this analysis we are able to probe the CN distribution of M31 and evaluate the global nova rate. Using models of the galactic surface brightness of M31, we show that the observed CN distribution consists of a separate bulge and disk population. We also show that the M31 bulge CN eruption rate per unit $r'$ flux is more than five times greater than that of the disk. We deduce a global M31 CN rate of $65^{+16}_{-15}$ year$^{-1}$, a value much higher than found by previous surveys. Using the global rate, we derive a M31 bulge rate of $38^{+15}_{-12}$ year$^{-1}$ and a disk rate of $27^{+19}_{-15}$ year$^{-1}$. Given our understanding of the completeness and an analysis of other sources of error, we conclude that the true global nova rate of M31 is at least 50% higher than was previously thought and this has consequent implications for the presumed CN rate in the Milky Way. We deduce a Galactic bulge rate of $14^{+6}_{-5}$ year$^{-1}$, a disk rate of $20^{+14}_{-11}$ year$^{-1}$ and a global Galactic rate of $34^{+15}_{-12}$ year$^{-1}$.

Perhaps as many as 30 parallax microlensing events are known, thanks to the efforts of the MACHO, OGLE, EROS and MOA experiments monitoring the bulge. Using Galactic models, we construct mock catalogues of microlensing light curves towards the bulge, allowing for the uneven sampling and observational error bars of the OGLE-II experiment. The fraction of parallax events with delta chi^2 > 50 in the OGLE-II database is around ~1%, though higher fractions are reported by some other surveys. This is in accord with expectations from standard Galactic models. The fraction of parallax events depends strongly on the Einstein crossing time (t_E), being less than 5% at t_E = 50 days but rising to 50% at t_E > 1 yr. We find that the existence of parallax signatures is essentially controlled by the acceleration of the observer normalised to the projected Einstein radius on the observer plane divided by t_E^2. The properties of the parallax events - time-scales, projected velocities, source and lens locations - in our mock catalogues are analysed. Typically, ~38% of parallax events are caused by a disk star microlensing a bulge source, while \~33% are caused by a disk star microlensing a disk source (of these disk sources, one sixth are at a distance of 5 kpc or less). There is a significant shift in mean time-scale from 32 d for all events to ~130d for our parallax events. There are corresponding shifts for other parameters, such as the lens-source velocity projected onto the observer plane (~1110 km/s for all events versus ~80 km/s for parallax events) and the lens distance (6.7 kpc versus 3.7 kpc). We also assess the performance of parallax mass estimators and investigate whether our mock catalogue can reproduce events with features similar to a number of conjectured `black hole' lens candidates.

The POINT-AGAPE collaboration surveyed M31 with the primary goal of optical detection of microlensing events, yet its data catalogue is also a prime source of lightcurves of variable and transient objects, including classical novae (CNe). A reliable means of identification, combined with a thorough survey of the variable objects in M31, provides an excellent opportunity to locate and study an entire galactic population of CNe. This paper presents a set of 440 neural networks, working in 44 committees, designed specifically to identify fast CNe. The networks are developed using training sets consisting of simulated novae and POINT-AGAPE lightcurves, in a novel variation on K-fold cross-validation. They use the binned, normalised power spectra of the lightcurves as input units. The networks successfully identify 9 of the 13 previously identified M31 CNe within their optimal working range (and 11 out of 13 if the network error bars are taken into account). They provide a catalogue of 19 new candidate fast CNe, of which 4 are strongly favoured.

Many of the science goals of the GAIA mission, especially for bursting or time-varying phenomena like supernovae or microlensing, require an early identification, analysis and release of preliminary data. The alerting on rare and unusual events is the scientific equivalent of the finding of needles in haystacks or the panning for gold-dust in rivers. In modern signal processing, such novelty detection is routinely performed with self-organising maps (SOMs), which are an unsupervised clustering algorithm invented by Kohonen. Here, we describe the application of SOMS to the classification of data provided by large-scale surveys such as GAIA and to the despatching of scientific alerts. We illustrate our ideas by processing the publically available OGLE II dataset towards the Bulge, identifying major classes of variable stars (such as novae, small amplitude red giant variables, eclipsing binaries and so on) and extracting the rare, discrepant lightcurves from which the alerts can be drawn.

This article reviews the life and death of a scientific theory.

An automated search is carried out for microlensing events using a catalogue of 44554 variable superpixel lightcurves derived from our three-year monitoring program of M31. Each step of our candidate selection is objective and reproducible by a computer. Our search is unrestricted, in the sense that it has no explicit timescale cut. So, it must overcome the awkward problem of distinguishing long-timescale microlensing events from long-period stellar variables. The basis of the selection algorithm is the fitting of the superpixel lightcurves to two different theoretical models, using variable star and blended microlensing templates. Only if microlensing is preferred is an event retained as a possible candidate. Further cuts are made with regard to (i) sampling, (ii) goodness of fit of the peak to a Paczynski curve, (iii) consistency of the microlensing hypothesis with the absence of a resolved source, (iv) achromaticity, (v) position in the colour-magnitude diagram and (vi) signal-to-noise ratio. Our results are reported in terms of first-level candidates, which are the most trustworthy, and second-level candidates, which are possible microlensing but have lower signal-to-noise and are more questionable. The pipeline leaves just 3 first-level candidates, all of which have very short full-width half-maximum timescale (<5 days) and 3 second-level candidates, which have timescales of 31, 36 and 51 days respectively. We also show 16 third-level lightcurves, as an illustration of the events that just fail the threshold for designation as microlensing candidates. They are almost certainly mainly variable stars. Two of the 3 first-level candidates correspond to known events (PA 00-S3 and PA 00-S4) already reported by the POINT-AGAPE project. The remaining first-level candidate is new.

Automatic classification of variability is now possible with tools like neural networks. Here, we present two neural networks for the identification of microlensing events -- the first discriminates against variable stars and the second against supernovae. The inputs to the networks include parameters describing the shape and the size of the lightcurve, together with colour of the event. The network computes the posterior probability of microlensing, together with an estimate of the likely error. An algorithm is devised for direct calculation of the microlensing rate from the output of the neural networks. We present a new analysis of the microlensing candidates towards the Large Magellanic Cloud (LMC). The neural networks confirm the microlensing nature of only 7 of the possible 17 events identified by the MACHO experiment. This suggests that earlier estimates of the microlensing optical depth towards the LMC may have been overestimated. A smaller number of events is consistent with the assumption that all the microlensing events are caused by the known stellar populations in the outer Galaxy/LMC.

This note addresses possible applications of the Tikhonov regularization to image reconstruction of gravitational lens systems. Several modifications of the regularization algorithm are discussed. Our illustrative example is the close quadruple gravitational lens QSO 2237+0305 (Einstein Cross). The restored image of the lens is decomposed into two parts - the quasar components and the background galaxy.

GAIA is the ``super-Hipparcos'' satellite scheduled for launch in 2010 by the European Space Agency. It is a scanning satellite that carries out multi-colour, multi-epoch photometry on all objects brighter than 20th magnitude. We conduct detailed simulations of supernovae (SNe) detection by GAIA. Supernovae of each type are chosen according to the observed distributions of absolute magnitudes, and located in nearby galaxies according to the local large-scale structure. Using an extinction model of the Galaxy and the scanning law of the GAIA satellite, we calculate how many SNe are detectable as a function of the phase of the lightcurve. Our study shows that GAIA will report data on 21400 SNe during the five-year mission lifetime, of which 14300 are SNe Ia, 1400 are SNe Ib/c and 5700 are SNe II. Using the simulations, we estimate that the numbers caught before maximum are 6300 SNe Ia, 500 SNe Ib/c and 1700 SNe II. GAIA will issue about 5 SNe alerts a day. The most distant SNe accessible to GAIA are at a redshift of 0.14 and so GAIA will provide a huge sample of local SNe. There will be many examples of the rarer subluminous events, over-luminous events, SNe Ib/c and SNe II-L. SNe rates will be found as a function of galaxy type, as well as extinction and position in the host galaxy. Amongst other applications, there may be about 26 SNe each year for which detection of gravitational waves is possible and about 180 SNe each year for which detection of gamma-rays is possible.

This paper pioneers the use of neural networks to provide a fast and automatic way to classify lightcurves in massive photometric datasets. As an example, we provide a working neural network that can distinguish microlensing lightcurves from other forms of variability, such as eruptive, pulsating, cataclysmic and eclipsing variable stars. The network has five input neurons, a hidden layer of five neurons and one output neuron. The five input variables for the network are extracted by spectral analysis from the lightcurve datapoints and are optimised for the identification of a single, symmetric, microlensing bump. The output of the network is the posterior probability of microlensing. The committee of neural networks successfully passes tests on noisy data taken by the MACHO collaboration. When used to process 5000 lightcurves on a typical tile towards the bulge, the network cleanly identifies the single microlensing event. When fed with a sub-sample of 36 lightcurves identified by the MACHO collaboration as microlensing, the network corroborates this verdict in the case of 27 events, but classifies the remaining 9 events as other forms of variability. For some of these discrepant events, it looks as though there are secondary bumps or the bump is noisy or not properly contained. Neural networks naturally allow for the possibility of novelty detection -- that is, new or unexpected phenomena which we may want to follow up. The advantages of neural networks for microlensing rate calculations, as well as the future developments of massive variability surveys, are both briefly discussed.

GAIA is an astrometric satellite which has been approved by the European Space Agency for launch in about 2010. It will measure the angles between objects in fields that are separated on the sky by about a radian. Data will stream continuously at 1 Mbps from GAIA's three telescopes, providing information on the positions and velocities of the billion or more astrophysical objects brighter than 20th magnitude. The motion of objects across the sky caused by their space motion and their parallactic motion, and the variability of objects in 15 wavebands, will be measured because each object is observed at least 150 times during the 5 year mission lifetime. From the raw time series, a three-dimensional movie of the motions of stars in the Galaxy will be synthesized. The alert despatcher in the GAIA mission will provide forewarning of many kinds of bursting and variable phenomena. We discuss two applications of the GAIA movie show in detail here. First, the astrometric microlensing signal will allow us to take a complete inventory of all objects -- no matter how dark -- in the solar neighbourhood. Second, the catalogue of supernovae detections will be the largest ever taken. It will provide opportunities both for follow-ups in other wavebands and numerous examples of scarce phenomena (like subluminous SNe).

At any instant, there are roughly 1000 ongoing microlensing events to sources brighter than 20th magnitude in the Milky Way Galaxy. Large-scale maps of the microlensing optical depth and the mean timescale are constructed for a number of models of the Galactic bar and disk, incorporating the effects of streaming and spiral structure. Freudenreich's model can reproduce the high optical depths towards the Bulge. It is also in good agreement with the data towards the spiral arms (except for the gamma Norma field). Spiral structure tends to increase the optical depth by about 20% and the mean timescale by about 100%. Different bar morphologies give characteristically different shaped contours, especially at low Galactic latitudes (|b| < 2 degrees). These could be traced out with a K band microlensing survey, consuming about 100 minutes per night on a telescope like VISTA.

GAIA is the ``super-Hipparcos'' survey satellite selected as a Cornerstone 6 mission by the European Space Agency. GAIA can measure microlensing by the small excursions of the light centroid that occur during events. The all-sky source-averaged astrometric microlensing optical depth is about 10^{-5}. Some 25000 sources will have a significant variation of the centroid shift, together with a closest approach, during the lifetime of the mission. A covariance analysis is used to study the propagation of errors and the estimation of parameters from realistic sampling of the GAIA datastream of transits in the along-scan direction during microlensing events. Monte Carlo simulations are used to study the 2500 events for which the mass can be recovered with an error of less than 50 per cent. These high quality events are dominated by disk lenses within a few tens of parsecs and source stars within a few hundred parsecs. We show that the local mass function can be recovered from the high quality sample to good accuracy. GAIA is the first instrument with the capabilities of measuring the mass locally in very faint objects like black holes and very cool white and brown dwarfs. For only 5 per cent of all astrometric events will GAIA record even one photometric datapoint. There is a need for a dedicated telescope that densely samples the Galactic Centre and spiral arms, as this can improve the accuracy of parameter estimation by a factor of about 10.