The dwarf satellite galaxies of the Milky Way Crater II and Antlia II have uncommonly low dynamical mass densities, due to their large size and low velocity dispersion. Previous work have failed to identify formation scenarios within the $\Lambda$CDM framework and have invoked cored dark matter haloes, processed by tides. I show that the tidal evolution of $\Lambda$CDM NFW haloes is richer than previously recognised: tidal heating causes the innermost regions of haloes that fall short of the mass-concentration relation to expand significantly, resulting in the formation of giant, kinematically cold satellites like Crater II and Antlia II. Furthermore, while the satellite is reaching apocenter, extra-tidal material can cause an even more inflated appearance. When present, as likely for the larger Antlia II, nominally unbound material can be recognised thanks to its somewhat hotter kinematics and line-of-sight velocity gradient. Contrary to other formation scenarios, Crater II and Antlia II may well have experienced very little mass loss, as in fact hinted by their observed metallicity. If indeed a satellite of NGC1052, tidal evolution of a low-concentration halo may similarly have led to the formation of NGC1052-DF2.

We use Monte Carlo $\Lambda$CDM assembly histories, minor-merger N-body simulations and empirical relations between halo mass and the globular cluster (GC) abundance to study the kinematical properties of halo GCs in massive galaxies, with $M_{vir}(z=0)=10^{13.5}M_\odot$. While the accreted stellar halo is dominated by the contributions of massive satellites, we show that satellites with low virial mass (i.e. low satellite-to-host virial mass ratio, VMRs) are important contributors to the population of accreted GCs. The relative contribution of accretion events with low VMRs is highest for the halo population of blue GCs and gradually decreases for red GCs and accreted stars. As a consequence of the reduced efficiency of dynamical friction on minor mergers, our populations of cosmologically accreted blue GCs are systematically more spatially extended and have higher velocity dispersions than for red GCs, in agreement with observations. For the same reason, assembly histories including a higher fraction of minor mergers result in more spatially extended GC populations. GC line-of-sight velocity distributions featuring negative values of the kurtosis $\kappa$, as recently observed, are ubiquitous in our models. Therefore, $\kappa<0$ is not at odds with an accretion scenario, and in fact a fingerprint of the important contribution of minor mergers. However, our populations of accreted GCs remain mostly radially biased, with profiles of the anisotropy parameter $\beta$ that are mildly radial in the center ($\beta(r<10~{\rm kpc})\sim0.2$) and strongly radially anisotropic at large galactocentric distances ($\beta(r>30~{\rm kpc})\gtrsim0.6$), for both red and blue populations.

Shells are low surface brightness tidal debris that appear as interleaved caustics with large opening angles, often situated on both sides of the galaxy center. In this paper, we study the incidence and formation processes of shell galaxies in the cosmological gravity+hydrodynamics Illustris simulation. We identify shells at redshift z=0 using stellar surface density maps, and we use stellar history catalogs to trace the birth, trajectory and progenitors of each individual star particle contributing to the tidal feature. Out of a sample of the 220 most massive galaxies in Illustris ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_{\odot}$), $18\%\pm3\%$ of the galaxies exhibit shells. This fraction increases with increasing mass cut: higher mass galaxies are more likely to have stellar shells. Furthermore, the fraction of massive galaxies that exhibit shells decreases with increasing redshift. We find that shell galaxies observed at redshift $z=0$ form preferentially through relatively major mergers ($\gtrsim$1:10 in stellar mass ratio). Progenitors are accreted on low angular momentum orbits, in a preferred time-window between $\sim$4 and 8 Gyrs ago. Our study indicates that, due to dynamical friction, more massive satellites are allowed to probe a wider range of impact parameters at accretion time, while small companions need almost purely radial infall trajectories in order to produce shells. We also find a number of special cases, as a consequence of the additional complexity introduced by the cosmological setting. These include galaxies with multiple shell-forming progenitors, satellite-of-satellites also forming shells, or satellites that fail to produce shells due to multiple major mergers happening in quick succession.

We present a novel method to constrain the mass of ultra-light bosons as the dark matter using stellar streams formed by disrupting Globular Clusters in the Milky Way. The turbulent density field of Fuzzy Dark Matter (FDM) haloes results in perturbations and dynamical heating of thin streams. Using numerical simulations based on an effective model, we explore the magnitude of this phenomenon and show that this is observable for the range of axion masses $m_a$ that is interesting for alleviating the `small-scale problems' of $\Lambda$CDM. We derive an analytical model for the thickening of thin stellar streams and obtain an early conservative lower limit for the boson mass of $m_a > 1.5 \times 10^{-22}~$eV, using pre-Gaia literature data for six Milky Way streams and after marginalizing over physical parameters. This demonstrates the great promise for using this novel dynamical method as a fully independent probe of FDM, to complement results based on Lyman-$\alpha$ forest data.

We use the observed abundances of ultra-diffuse galaxies (UDGs) in clusters and groups and \Lambda CDM subhalo mass functions to put constraints on the distribution {of present-day halo masses of satellite} UDGs. If all of the most massive subhaloes in the cluster host a UDG, UDGs occupy all subhaloes with \log M_{sub}/M_\odot\gtrsim11. For a model in which the efficiency of UDG formation is higher around some characteristic halo mass, higher fractions of massive UDGs require larger spreads in the UDG mass distribution. In a cluster with a virial mass of 10^{15}M_\odot, the 90% upper limit for the fraction of UDGs with \log M_{sub}/M_\odot>12 is 7%, occupying 70% of all cluster subhaloes above the same mass. To reproduce the observed abundances, however, the mass distribution of satellite UDGs has to be broad, with >30% having \log M_{sub}/M_\odot<10.9. This strongly supports that UDGs are part of a continuous distribution in which a majority are hosted by low mass haloes. The abundance of satellite UDGs may fall short of the linear relation with the cluster/group mass M_{host} in low-mass hosts, \log M_{host}/M_\odot\sim 12. Characterising these deviations -- or the lack thereof -- will allow for stringent constraints on the UDG virial mass distribution.

We use data from the HST Coma Cluster Treasury program to assess the richness of the Globular Cluster Systems (GCSs) of 54 Coma ultra-diffuse galaxies (UDGs), 18 of which have a half-light radius exceeding 1.5 kpc. We use a hierarchical Bayesian method tested on a large number of mock datasets to account consistently for the high and spatially varying background counts in Coma. These include both background galaxies and intra-cluster GCs (ICGCs), which are disentangled from the population of member GCs in a probabilistic fashion. We find no candidate for a GCS as rich as that of the Milky Way, our sample has GCSs typical of dwarf galaxies. For the standard relation between GCS richness and halo mass, 33 galaxies have a virial mass $M_{vir}\leq10^{11}M_\odot$ at 90% probability. Only three have $M_{vir}>10^{11}M_\odot$ with the same confidence. The mean colour and spread in colour of the UDG GCs are indistinguishable from those of the abundant population of ICGCs. The majority of UDGs in our sample are consistent with the relation between stellar mass and GC richness of 'normal' dwarf galaxies. Nine systems, however, display GCSs that are richer by a factor of 3 or more (at 90% probability). Six of these have sizes $\lesssim1.4$ kpc. Our results imply that the physical mechanisms responsible for the extended size of the UDGs and for the enhanced GC richness of some cluster dwarfs are at most weakly correlated.

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.

Stellar shells are low surface brightness arcs of overdense stellar regions, extending to large galactocentric distances. In a companion study, we identified 39 shell galaxies in a sample of 220 massive ellipticals ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_\odot$) from the Illustris cosmological simulation. We used stellar history catalogs to trace the history of each individual star particle inside the shell substructures, and we found that shells in high-mass galaxies form through mergers with massive satellites (stellar mass ratios $\mu_{\mathrm{stars}}\gtrsim1:10$). Using the same sample of shell galaxies, the current study extends the stellar history catalogs in order to investigate the metallicity of stellar shells around massive galaxies. Our results indicate that outer shells are often times more metal-rich than the surrounding stellar material in a galaxy's halo. For a galaxy with two different satellites forming $z=0$ shells, we find a significant difference in the metallicity of the shells produced by each progenitor. We also find that shell galaxies have higher mass-weighted logarithmic metallicities ([Z/H]) at $2$-$4\,\mathrm{R}_{\mathrm{eff}}$ compared to galaxies without shells. Our results indicate that observations comparing the metallicities of stars in tidal features, such as shells, to the average metallicities in the stellar halo can provide information about the assembly histories of galaxies.

The recent detection of two faint and extended star clusters in the central regions of two Local Group dwarf galaxies, Eridanus II and Andromeda XXV, raises the question of whether clusters with such low densities can survive the tidal field of cold dark matter haloes with central density cusps. Using both analytic arguments and a suite of collisionless N-body simulations, I show that these clusters are extremely fragile and quickly disrupted in the presence of central cusps $\rho\sim r^{-\alpha}$ with $\alpha\gtrsim 0.2$. Furthermore, the scenario in which the clusters where originally more massive and sank to the center of the halo requires extreme fine tuning and does not naturally reproduce the observed systems. In turn, these clusters are long lived in cored haloes, whose central regions are safe shelters for $\alpha\lesssim 0.2$. The only viable scenario for hosts that have preserved their primoridal cusp to the present time is that the clusters formed at rest at the bottom of the potential, which is easily tested by measurement of the clusters proper velocity within the host. This offers means to readily probe the central density profile of two dwarf galaxies as faint as $L_V\sim5\times 10^5 L_\odot$ and $L_V\sim6\times10^4 L_\odot$, in which stellar feedback is unlikely to be effective.

Theory and observations agree that the accreted stellar halos (ASHs) of Milky Way-like galaxies display significant scatter. I take advantage of this stochasticity to invert the link between halo assembly history (HAH) and ASH, using mock ASHs corresponding to 750 $\Lambda$CDM HAHs, sharing a final virial mass of $M_{h}(z=0)=10^{12.25}M_\odot$. Hosts with poor/rich ASHs assemble following orthogonal growth-patterns. Hosts with rich ASHs experience accretion events (AEs) with high virial mass ratios (HVMRs, $M_s/M_h\gtrsim 0.1$) at $0.5\lesssim z_{infall}\lesssim1.5$, in a phase of fast growth. This maximizes the accreted stellar mass under the condition these satellites are disrupted by $z=0$. At similar times, hosts with poor ASHs grow slowly through minor mergers, with only very recent HVMR AEs: this results in a globally more abundant satellite population and in distinctive surviving massive satellites (stellar mass $\log M_{s,*}/M_\odot\gtrsim 9$). Several properties of the Milky Way are in agreement with the predictions of this framework for hosts with poor, concentrated ASHs, including: i) the recent infall of Sagittarius and Magellanic Clouds, ii) the likely higher-than-average concentration of its dark halo, iii) the signatures of fast chemical enrichment of a sizable fraction of its halo stellar populations.

The accreted component of stellar halos is composed of the contributions of several satellites, falling onto their host with their different masses, at different times, on different orbits. This work uses a suite of idealised, collisionless N-body simulations of minor mergers and a particle tagging technique to understand how these different ingredients shape each contribution to the accreted halo, in both density and kinematics. I find that more massive satellites deposit their stars deeper into the gravitational potential of the host, with a clear segregation enforced by dynamical friction. Earlier accretion events contribute more to the inner regions of the halo; more concentrated subhaloes sink deeper through increased dynamical friction. The orbital circularity of the progenitor at infall is only important for low-mass satellites: dynamical friction efficiently radialises the most massive minor mergers erasing the imprint of the infall orbit for satellite-to-host virial mass ratios $\gtrsim1/20$. The kinematics of the stars contributed by each satellite is also ordered with satellite mass: low-mass satellites contribute fast-moving populations, in both ordered rotation and radial velocity dispersion. In turn, contributions by massive satellites have lower velocity dispersion and lose their angular momentum to dynamical friction, resulting in a strong radial anisotropy.

I use a library of controlled minor merger N-body simulations, a particle tagging technique and Monte Carlo generated $\Lambda$CDM accretion histories to study the highly stochastic process of stellar deposition onto the accreted stellar halos (ASHs) of $L_*$ galaxies. I explore the main physical mechanisms that drive the connection between the accretion history and the density profile of the ASH. I find that: i) through dynamical friction, more massive satellites are more effective at delivering their stars deeper into the host; ii) as a consequence, ASHs feature a negative gradient between radius and the local mass-weighed virial satellite-to-host mass ratio; iii) in $L_*$ galaxies, most ASHs feature a density profile that steepens towards sharper logarithmic slopes at increasing radii, though with significant halo-to-halo scatter; iv) the ASHs with the largest total ex-situ mass are such because of the chance accretion of a small number of massive satellites (rather than of a large number of low-mass ones).

As a result of their internal dynamical coherence, thin stellar streams formed by disrupting globular clusters (GCs) can act as detectors of dark matter (DM) substructure in the Galactic halo. Perturbations induced by close flybys amplify into detectable density gaps, providing a probe both of the abundance and of the masses of DM subhaloes. Here, we use N-body simulations to show that the Galactic population of giant molecular clouds (GMCs) can also produce gaps (and clumps) in GC streams, and so may confuse the detection of DM subhaloes. We explore the cases of streams analogous to the observed Palomar 5 and GD1 systems, quantifying the expected incidence of structure caused by GMC perturbations. Deep observations should detect such disturbances regardless of the substructure content of the Milky Way's halo. Detailed modelling will be needed to demonstrate that any detected gaps or clumps were produced by DM subhaloes rather than by molecular clouds.

Recent observations have revealed the existence of an abundant population of faint, low surface brightness (SB) galaxies, which appear to be numerous and ubiquitous in nearby galaxy clusters, including the Virgo, Coma and Fornax clusters. With median stellar masses of dwarf galaxies, these ultra-diffuse galaxies (UDGs) have unexpectedly large sizes, corresponding to a mean SB of $24\lesssim\langle\mu_e\rangle_r\ {\rm mag}^{-1} {\rm arcsec}^2\lesssim27$ within the effective radius. We show that the UDG population represents the tail of galaxies formed in dwarf-sized haloes with higher-than-average angular momentum. By adopting the standard model of disk formation -- in which the size of galaxies is set by the spin of the halo -- we recover both the abundance of UDGs as a function of the host cluster mass and the distribution of sizes within the UDG population. According to this model, UDGs are not failed $L_*$ galaxies, but genuine dwarfs, and their low SB is not uniquely connected to the harsh cluster environment. We therefore expect a correspondingly abundant population of UDGs in the field, with likely different morphologies and colours.

We study the chemodynamical properties of the Carina dwarf spheroidal by combining an intermediate spectroscopic resolution dataset of more than 900 red giant and red clump stars, with high-precision photometry to derive the atmospheric parameters, metallicities and age estimates for our targets. Within the red giant branch population, we find evidence for the presence of three distinct stellar sub-populations with different metallicities, spatial distributions, kinematics and ages. As in the Fornax and Sculptor dwarf spheroidals, the subpopulation with the lowest average metallicity is more extended and kinematically hotter than all other populations. However, we identify an inversion in the parallel ordering of metallicity, kinematics and characteristic length scale in the two most metal rich subpopulations, which therefore do not contribute to a global negative chemical gradient. Contrary to common trends in the chemical properties with radius, the metal richest population is more extended and mildly kinematically hotter than the main component of intermediate metallicity. More investigations are required to ascertain the nature of this inversion, but we comment on the mechanisms that might have caused it.

(abridged) We present a dynamical analysis of the extended stellar stream encircling NGC 1097. Within a statistical framework, we model its surface brightness using mock streams as in Amorisco (2015) and deep imaging data from the CHART32 telescope (Stellar Tidal Stream Survey). We reconstruct the post-infall evolution of the progenitor, which has experienced 3 pericentric passages and lost more than 2 orders of magnitude in mass. At infall, $5.4\pm0.6$ Gyr ago, the progenitor was a disky dwarf with mass of $\log_{10}[m(<3.4\pm1 {\rm kpc})/ M_\odot]=10.35\pm0.25$. We illustrate how the 90$^\circ$ turn in the stream, identifying the `dog leg', is the signature of the progenitor's prograde rotation. Today, the remnant is a nucleated dwarf, with a LOS velocity of $v_{\rm p, los}^{\rm obs}=-30\pm 30$ kms$^{-1}$, and a luminosity of $3.3\times 10^7 L_{V,\odot}$ (Galianni et al. 2010). Our independent analysis predicts $v_{\rm p, los}=-51^{-17}_{+14}$ kms$^{-1}$, and measures $\log_{10}(m/ M_\odot)=7.4^{+0.6}_{-0.8}$, so that the compact nucleus is soon becoming a low-luminosity UCD. We find that NGC 1097 has a mass of $M_{200}=1.8^{+0.5}_{-0.4} \times 10^{12}\; M_{\odot}$, and its concentration $c_{200}=6.7^{+2.4}_{-1.3}$ is in agreement with LCDM. The stream is described almost down to the noise in a spherical host potential, we find this would not be possible if the halo was substantially triaxial at large radii. Its morphology shows that the slope of the total density profile bends from an inner $\gamma(r_{\rm peri})=1.5\pm0.15$. The progenitor's orbit reaches $r_{\rm apo}=150\pm 15$ kpc, more than a half of the virial radius of the host, so that, for the first time on an individual extragalactic halo, we measure the outer density slope, $\gamma(0.6r_{200,c})=3.9\pm0.5$. This demonstrates the promise of the newborn field of detailed, statistical modelling of extragalactic tidal streams.

I present an organic description of the regimes of collisionless tidal streams and define the orderings between the physical quantities that shape their morphology. Three fundamental dichotomies are identified in the form of dimensionless inequalities. These govern i) the speed of the stream's growth, ii) its internal coherence, iii) its thickness or opening angles. The mechanisms that regulate such main properties are analysed. The slope of the host's density profile influences the speed of the stream's growth, in both length and width, as steeper profiles enhance differential streaming. Internal coherence is the requirement for the appearance of substructure in tidal debris, and I concentrate on the `feathering' typical of GC streams. Overdensities are associated with minima in the relative streaming velocity of the stream members. For streams with high circularity, these are caused by the epicyclic oscillations of stars; however, for highly non-circular progenitor's orbits, substructure is caused by the oscillating differences in energy and actions with which material is shed at different orbital phases of the progenitor. This modulation results in different streaming speeds: the streakline of material shed between two successive apocentric passages is folded along its length, pulled at its centre by the faster streaming of particles released near pericenter, which are therefore more widely scattered. When the stream is coherent enough, this mechanism is potentially capable of generating a bimodal profile in the density distributions of the longer wraps of more massive progenitors, which I dub `bifurcations'. The conditions for internal coherence are explored and I comment on the cases of Palomar 5, Willman 1, the Anticenter and Sagittarius' streams. Analytical methods are accompanied by numerical experiments, performed using a purposely built generative model, also presented here.

Driven by gravity, massive structures like galaxies and clusters of galaxies are believed to grow continuously through hierarchical merging and accretion of smaller systems. Observational evidence of accretion events is provided by the coherent stellar streams crossing the outer haloes of massive galaxies, such as the Milky Way or Andromeda. At similar mass-scales, around $10^{11}$ solar masses in stars, further observational evidence of merging activity is also ample. Mergers of lower-mass galaxies are expected within the hierarchical process of galaxy formation, but have hitherto not been seen for galaxies with less than about $10^9$ solar masses in stars. Here, we report the kinematic detection of a stellar stream in one of the satellite galaxies of Andromeda, the dwarf spheroidal galaxy Andromeda II, which has a mass of only $10^7$ solar masses in stars. The properties of the stream show that we are observing the remnant of a merger between two dwarf galaxies. This had a dramatic influence on the dynamics of the remnant, which is now rotating around its projected major axis. The stellar stream in Andromeda II illustrates the scale-free character of the formation of galaxies, down to the lowest galactic mass scales.

We combine the detailed Star Formation Histories of the Fornax and Sculptor dwarf Spheroidals with the Mass Assembly History of their dark matter (DM) halo progenitors to estimate if the energy deposited by Supernova type II (SNeII) is sufficient to create a substantial DM core. Assuming the efficiency of energy injection of the SNeII into DM particles is $\epsilon_{\rm gc}=0.05$, we find that a single early episode, $z \gtrsim z_{\rm infall}$, that combines the energy of all SNeII due to explode over 0.5 Gyr, is sufficient to create a core of several hundred parsecs in both Sculptor and Fornax. Therefore, our results suggest that it is energetically plausible to form cores in Cold Dark Matter (CDM) halos via early episodic gas outflows triggered by SNeII. Furthermore, based on CDM merger rates and phase-space density considerations, we argue that the probability of a subsequent complete regeneration of the cusp is small for a substantial fraction of dwarf-size haloes.

We exploit the detection of three distinct stellar subpopulations in the red giant branch of the Fornax dwarf Spheroidal to probe its density distribution. This allows us to resolve directly the evolution with radius of the dark matter mass profile. We find that a cored dark matter halo provides a perfect fit to the data, being consistent with all three stellar populations well within 1-sigma, and for the first time we are able to put constraints on the core size of such a halo. With respect to previous work, we do not strengthen the statistical exclusion of a dark matter cusp in Fornax, but we find that Navarro-Frenk-White haloes would be required to have unrealistically large scale radii in order to be compatible with the data, hence low values of the concentration parameter. We are then forced to conclude that the Fornax dwarf Spheroidal sits within a dark matter halo having a constant density core, with a core size of between 0.6 and 1.8 kpc.

We present compelling evidence for the complexity of the Fornax dwarf spheroidal. By disentangling three different stellar subpopulations in its red giant branch, we are able to study in detail the dependence between kinematics and metallicity. A well-defined ordering in velocity dispersion, spatial concentration, and metallicity is evident in the subpopulations. We also present evidence for a significant misalignment between the angular momentum vectors of the old and intermediate-age populations. According to the HST measurement of Fornax's proper motion, this corresponds to counter-rotation. These ingredients are used to construct a novel evolutionary history of the Fornax dwarf spheroidal, characterized as a late merger of a bound pair.

We develop a method to extract the shape information of line profiles from discrete kinematic data. The Gauss-Hermite expansion, which is widely used to describe the line of sight velocity distributions extracted from absorption spectra of elliptical galaxies, is not readily applicable to samples of discrete stellar velocity measurements, accompanied by individual measurement errors and probabilities of membership. We introduce two parameter families of probability distributions describing symmetric and asymmetric distortions of the line profiles from Gaussianity. These are used as the basis of a maximum likelihood estimator to quantify the shape of the line profiles. Tests show that the method outperforms a Gauss-Hermite expansion for discrete data, with a lower limit for the relative gain of approx 2 for sample sizes N approx 800. To ensure that our methods can give reliable descriptions of the shape, we develop an efficient test to assess the statistical quality of the obtained fit. As an application, we turn our attention to the discrete velocity datasets of the dwarf spheroidals of the Milky Way. In Sculptor, Carina and Sextans the symmetric deviations are consistent with velocity distributions more peaked than Gaussian. In Fornax, instead, there is an evolution in the symmetric deviations of the line profile from a peakier to more flat-topped distribution on moving outwards. These results suggest a radially biased orbital structure for the outer parts of Sculptor, Carina and Sextans. On the other hand, tangential anisotropy is favoured in Fornax. This is all consistent with a picture in which Fornax may have had a different evolutionary history to Sculptor, Carina and Sextans.

A number of dwarf spheroidal (dSph) galaxies are known to contain a more extended, metal-poor population with a flattish velocity dispersion profile, and a more concentrated, metal-rich population with a velocity dispersion declining with radius. The two populations can be modelled with Michie-King distribution functions (DFs) in the isothermal and in the sharply truncated limits, respectively. We argue that the truncation of the metal-rich population can be traced back to the spatial distribution of the star forming gas. Suppose D is the exponent of the first non-constant term in the Taylor expansion of the total potential at the center (D=1 for NFW halos, D=2 for cored halos). Then, we show that the ratio of the half-light radii of the populations raised to the D/2 power must be smaller than the ratio of the line-of-sight velocity dispersions. Specialising to the case of the Sculptor dSph, we develop a technique to fit simultaneously both populations with Michie-King DFs. This enables us to determine the mass profile of the Sculptor dSph with unprecedented accuracy in the radial range between 0.2 and 1.2 kpc. We show that cored halo models are preferred over cusped halo models, with a likelihood ratio test rejecting NFW models at any significance level higher than 0.05%. Even more worryingly, the best-fit NFW models require concentrations with c<20, which is not in the cosmologically preferred range for dwarf galaxies. We conclude that the kinematics of multiple populations in dSphs provides a substantial new challenge for theories of galaxy formation, with the weight of available evidence strongly against dark matter cusps at the centre.

This paper introduces new phase-space models of dwarf spheroidal galaxies (dSphs). The stellar component has an isotropic, lowered isothermal (or King) distribution function. A physical basis for the isotropization of stellar velocities is given by tidal stirring, whilst the isothermality of the distribution function guarantees the observed flatness of the velocity dispersion profile in the inner parts. Our models reproduce the data on the half-light radius and line of sight central velocity dispersion of the dSphs. We show that different dark halo profiles -- whether cored or cusped -- lead to very similar mass estimates within one particular radius, namely 1.7 half-light radii. Deviations between mass measures due to different density profiles are substantially smaller than the uncertainties propagated by the observational errors. We produce a mass measure for each of the Milky Way dSphs and find that the two most massive are the most luminous, namely Sagittarius (~ 2.8 x 10^8 solar masses) and Fornax (~ 1.3 x 10^8 solar masses). The least massive of the Milky Way satellites are Willman 1 (~ 4 x 10^5 solar masses) and Segue 1 (~ 6 x 10^5 solar masses).

Disk-halo decompositions of galaxy rotation curves are generally performed in a parametric way. We construct self-consistent models of nonspherical isothermal halos embedding a zero-thickness disk, by assuming that the halo distribution function is a Maxwellian. The method developed here can be used to study other physically-based choices for the halo distribution function and the case of a disk accompanied by a bulge. In a preliminary investigation we note the existence of a fine tuning between the scalelengths R_{\Omega} and h, respectively characterizing the rise of the rotation curve and the luminosity profile of the disk, which surprisingly applies to both high surface brightness and low surface brightness galaxies. This empirical correlation identifies a much stronger conspiracy than the one required by the smoothness and flatness of the rotation curve (disk-halo conspiracy). The self-consistent models are characterized by smooth and flat rotation curves for very different disk-to-halo mass ratios, hence suggesting that conspiracy is not as dramatic as often imagined. For a typical rotation curve, with asymptotically flat rotation curve at V_{\infty} (the precise value of which can also be treated as a free parameter), and a typical density profile of the disk, self-consistent models are characterized by two dimensionless parameters, which correspond to the dimensional scales (the disk mass-to-light ratio M/L and the halo central density) of standard disk-halo decompositions. We show that if the rotation curve is decomposed by means of our self-consistent models, the disk-halo degeneracy is removed and typical rotation curves are fitted by models that are below the maximum-disk prescription. Similar results are obtained from a study of NGC 3198. Finally, we quantify the flattening of the spheroidal halo, which is significant, especially on the scale of the visible disk.

Context: Several spiral galaxies, as beautifully exhibited by the case of NGC 6946, display a prominent large-scale spiral structure in their gaseous outer disk. Such structure is often thought to pose a dynamical puzzle, because grand-design spiral structure is traditionally interpreted as the result of density waves carried mostly in the stellar disk. Aims. Here we argue that the outer spiral arms in the cold gas outside the bright optical disk actually have a natural interpretation as the manifestation of the mechanism that excites grand-design spiral structure in the main, star-dominated body of the disk: the excitation is driven by angular momentum transport to the outer regions, through trailing density waves outside the corotation circle that can penetrate beyond the Outer Lindblad Resonance in the gaseous component of the disk. Methods: Because of conservation of the density wave action, these outgoing waves are likely to become more prominent in the outer disk and eventually reach non-linear amplitudes. To calculate the desired amplitude profiles, we make use of the theory of dispersive waves. Results: If the conditions beyond the optical radius allow for an approximate treatment in terms of a linear theory, we show that fitting the observed amplitude profiles leads to a quantitative test on the density of the disk material and thus on the dark matter distribution in the outer parts of the galaxy. Conclusions: This study is thus of interest to the general problem of the disk-halo decomposition of rotation curves.