We investigate the disc-halo connection in massive (Mstar/Msun>5e10) disc galaxies from the cosmological hydrodynamical simulations EAGLE and IllustrisTNG, and compare it with that inferred from the study of HI rotation curves in nearby massive spirals from the Spitzer Photometry and Accurate Rotation Curves (SPARC) dataset. We find that discrepancies between the the simulated and observed discs arise both on global and on local scales. Globally, the simulated discs inhabit halos that are a factor ~4 (in EAGLE) and ~2 (in IllustrisTNG) more massive than those derived from the rotation curve analysis of the observed dataset. We also use synthetic rotation curves of the simulated discs to demonstrate that the recovery of the halo masses from rotation curves are not systematically biased. We find that the simulations predict dark-matter dominated systems with stellar-to-total enclosed mass ratios that are a factor of 1.5-2 smaller than real galaxies at all radii. This is an alternative manifestation of the `failed feedback problem', since it indicates that simulated halos hosting massive discs have been too inefficient at converting their baryons into stars, possibly due to an overly efficient stellar and/or AGN feedback implementation.

We study the gas kinematics of a sample of six isolated gas-rich low surface brightness galaxies, of the class called ultra-diffuse galaxies (UDGs). These galaxies have recently been shown to be outliers from the baryonic Tully-Fisher relation (BTFR), as they rotate much slower than expected given their baryonic mass, and to have baryon fractions similar to the cosmological mean. By means of a 3D kinematic modelling fitting technique, we show that the HI in our UDGs is distributed in "thin" regularly rotating discs and we determine their rotation velocity and gas velocity dispersion. We revisit the BTFR adding galaxies from other studies. We find a previously unknown trend between the deviation from the BTFR and the disc scale length valid for dwarf galaxies with circular speeds < 45 km/s, with our UDGs being at the extreme end. Based on our findings, we suggest that the high baryon fractions of our UDGs may originate due to the fact that they have experienced weak stellar feedback, likely due to their low star formation rate surface densities, and as a result they did not eject significant amounts of gas out of their discs. At the same time, we find indications that our UDGs may have higher-than-average stellar specific angular momentum, which can explain their large optical scale lengths.

We study the gas kinematics traced by the 21-cm emission of a sample of six HI$-$rich low surface brightness galaxies classified as ultra-diffuse galaxies (UDGs). Using the 3D kinematic modelling code $\mathrm{^{3D}}$Barolo we derive robust circular velocities, revealing a startling feature: HI$-$rich UDGs are clear outliers from the baryonic Tully-Fisher relation, with circular velocities much lower than galaxies with similar baryonic mass. Notably, the baryon fraction of our UDG sample is consistent with the cosmological value: these UDGs are compatible with having no "missing baryons" within their virial radii. Moreover, the gravitational potential provided by the baryons is sufficient to account for the amplitude of the rotation curve out to the outermost measured point, contrary to other galaxies with similar circular velocities. We speculate that any formation scenario for these objects will require very inefficient feedback and a broad diversity in their inner dark matter content.

In $\Lambda$CDM cosmology, to first order, galaxies form out of the cooling of baryons within the virial radius of their dark matter halo. The fractions of mass and angular momentum retained in the baryonic and stellar components of disc galaxies put strong constraints on our understanding of galaxy formation. In this work, we derive the fraction of angular momentum retained in the stellar component of spirals, $f_j$, the global star formation efficiency $f_M$, and the ratio of the asymptotic circular velocity ($V_{\rm flat}$) to the virial velocity $f_V$, and their scatter, by fitting simultaneously the observed stellar mass-velocity (Tully-Fisher), size-mass, and mass-angular momentum (Fall) relations. We compare the goodness of fit of three models: (i) where the logarithm of $f_j$, $f_M$, and $f_V$ vary linearly with the logarithm of the observable $V_{\rm flat}$; (ii) where these values vary as a double power law; and (iii) where these values also vary as a double power law but with a prior imposed on $f_M$ such that it follows the expectations from widely used abundance matching models. We conclude that the scatter in these fractions is particularly small ($\sim 0.07$ dex) and that the linear model is by far statistically preferred to that with abundance matching priors. This indicates that the fundamental galaxy formation parameters are small-scatter single-slope monotonic functions of mass, instead of being complicated non-monotonic functions. This incidentally confirms that the most massive spiral galaxies should have turned nearly all the baryons associated with their haloes into stars. We call this the failed feedback problem.

We analyse the mass density distribution in the centres of galaxies across five orders of magnitude in mass range. Using high-quality spiral galaxy rotation curves and infrared photometry from SPARC, we conduct a systematic study of their central dark matter fraction ($f_{\rm DM}$) and their mass density slope ($\alpha$), within their effective radius. We show that lower-mass spiral galaxies are more dark matter dominated and have more shallow mass density slopes when compared with more massive galaxies, which have density profiles closer to isothermal. Low-mass ($M_{\star} \lesssim 10^{10}\, \rm M_{\odot}$) gas-rich spirals span a wide range of \fdm\ values, but systematically lower than in gas-poor systems of similar mass. With increasing galaxy mass, the values of \fdm\ decrease and the density profiles steepen. In the most massive late-type gas-poor galaxies, a possible flattening of these trends is observed. When comparing these results to massive ($M_{\star} \gtrsim 10^{10}M_{\odot}$) elliptical galaxies from SPIDER and to dwarf ellipticals from SMACKED, these trends result to be inverted. Hence, the values of both $f_{\rm DM}$ and $\alpha$, as a function of $M_{\star}$, exhibit a U-shape trend. At a fixed stellar mass, the mass density profiles in dwarf ellipticals are steeper than in spirals. These trends can be understood by stellar feedback from a more prolonged star formation period in spirals, causing a transformation of the initial steep density cusp to a more shallow profile via differential feedback efficiency by supernovae, and by galaxy mergers or AGN feedback in higher-mass galaxies.

The velocity distribution of stars is a sensitive probe of the gravitational potential of the Galaxy, and hence of its dark matter distribution. In particular, the shape of the dark halo (e.g. spherical, oblate, or prolate) determines velocity correlations, and different halo geometries are expected to result in measurable differences. Here we explore and interpret the correlations in the $(v_R, v_z)$-velocity distribution as a function of position in the Milky Way. We selected a high-quality sample of stars from the Gaia DR2 catalogue and characterised the orientation of the velocity distribution or tilt angle over a radial distance range of $[4-13]~$kpc and up to $3.5~$kpc away from the Galactic plane while taking into account the effects of the measurement errors. We find that the tilt angles change from spherical alignment in the inner Galaxy ($R\sim4~$kpc) towards more cylindrical alignments in the outer Galaxy ($R\sim11~$kpc) when using distances that take a global zero-point offset in the parallax of $-29~\mu$as. However, if the amplitude of this offset is underestimated, then the inferred tilt angles in the outer Galaxy only appear shallower and are intrinsically more consistent with spherical alignment for an offset as large as $-54~\mu$as. We further find that the tilt angles do not seem to strongly vary with Galactic azimuth and that different stellar populations depict similar tilt angles. Therefore we introduce a simple analytic function that describes the trends found over the full radial range. Since the systematic parallax errors in Gaia DR2 depend on celestial position, magnitude, and colour in complex ways, it is not possible to fully correct for them. Therefore it will be particularly important for dynamical modelling of the Milky Way to thoroughly characterise the systematics in astrometry in future Gaia data releases.

Dwarf galaxies are key objects for small-scale cosmological tests like the abundance problems or the planes-of-satellites problem. It is therefore a crucial task to get accurate information for as many nearby dwarf galaxies as possible. Using extremely deep, ground-based $V$ and $i$-band Subaru Suprime Cam photometry with a completeness of $i=27$ mag, we measure the tip of the red giant branch distance for the dwarf galaxy [TT2009] 25. This dwarf resides in the field around the Milky Way-analog NGC 891. By using a Bayesian approach, we measure a distance of $10.28^{+1.17}_{-1.73}$ Mpc, which is consistent with the distance of NGC 891, thus we confirm it as a member of NGC 891. The dwarf galaxy follows the scaling relations defined by the Local Group dwarfs. We do not find an extended stellar halo around [TT2009] 25. In the small field of view of 100 kpc covered by the survey, only one bright dwarf galaxy and the giant stream are apparent. This is comparable to the Milky Way, where one bright dwarfs reside in the same volume, as well as the Sagittarius stream - excluding satellites which are farther away but would be projected in the line-of-sight. It is thus imperative to survey for additional dwarf galaxies in a larger area around NGC 891 to test the abundance of dwarf galaxies and compare it to the number of satellites around the Milky Way.

It is commonly believed that galaxies use, throughout the Hubble time, a very small fraction of the baryons associated to their dark matter halos to form stars. This so-called low "star formation efficiency" $f_\star\equiv M_\star/f_{\rm b} M_{\rm halo}$, where $f_{\rm b}\equiv\Omega_{\rm b}/\Omega_{\rm c}$ is the cosmological baryon fraction, is expected to reach its peak at nearly $L^\ast$ (at efficiency $\approx 20\%$) and decline steeply at lower and higher masses. We have tested this using a sample of nearby star-forming galaxies, from dwarfs ($M_\star\simeq 10^7 M_\odot$) to high-mass spirals ($M_\star\simeq 10^{11} M_\odot$) with HI rotation curves and 3.6$\mu$m photometry. We fit the observed rotation curves with a Bayesian approach by varying three parameters, stellar mass-to-light ratio, halo concentration and mass. We found two surprising results: 1) the star formation efficiency is a monotonically increasing function of $M_\star$ with no sign of a decline at high masses, and 2) the most massive spirals ($M_\star\simeq 1-3 \times 10^{11} M_\odot$) have $f_\star\approx 0.3-1$, i.e. they have turned nearly all the baryons associated to their haloes into stars. These results imply that the most efficient galaxies at forming stars are massive spirals (not $L^\ast$ galaxies), they reach nearly 100% efficiency and thus, once both their cold and hot gas is considered into the baryon budget, they have virtually no missing baryons. Moreover, there is no evidence of mass quenching of the star formation occurring in galaxies up to halo masses of $M_{\rm halo}\approx {\rm a\, few}\times 10^{12} M_\odot$.

A new family of self-consistent DF-based models of stellar systems is explored. The stellar component of the models is described by a distribution function (DF) depending on the action integrals, previously used to model the Fornax dwarf spheroidal galaxy (dSph). The stellar component may cohabit with either a dark halo, also described by a DF, or with a massive central black hole. In all cases we solve for the model's self-consistent potential. Focussing on spherically symmetric models, we show how the stellar observables vary with the anisotropy prescribed by the DF, with the dominance and nature of the dark halo, and with the mass of the black hole. We show that precise fits to the observed surface brightness profiles of four globular clusters can be obtained for a wide range of prescribed velocity anisotropies. We also obtain precise fits to the observed projected densities of four dSphs. Finally, we present a three-component model of the Scupltor dSph with distinct DFs for the red and blue horizontal branch stars and the dark matter halo.

Large spectroscopic surveys of the Milky Way have revealed that a small population of stars in the halo have light element abundances comparable to those found in globular clusters. The favoured explanation for the peculiar abundances of these stars is that they originated inside a globular cluster and were subsequently lost. Using orbit calculations we assess the likelihood that an existing sample of 57 field stars with globular cluster-like CN band strength originated in any of the currently known Milky Way globular clusters. Using Sloan Digital Sky Survey and Gaia data, we determine orbits and integrals of motion of our sample of field stars, and use these values and metallicity to identify likely matches to globular clusters. The pivot hypothesis is that had these stars been stripped from such objects, they would have remained on very similar orbits. We find that ~ 70% of the sample of field stars have orbital properties consistent with the halo of the Milky Way; however, only 20 stars have likely orbital associations with an existing globular cluster. The remaining ~ 30% of the sample have orbits that place them in the outer Galactic disc. No cluster of similar metallicity is known on analogous disc orbits. The orbital properties of the halo stars seem to be compatible with the globular cluster escapee scenario. The stars in the outer disc are particularly surprising and deserve further investigation to establish their nature.

We estimate the mass of the inner ($<20$ kpc) Milky Way and the axis ratio of its dark matter halo using globular clusters as tracers. At the same time, we constrain the phase-space distribution of the globular cluster system. We use the Gaia DR2 catalogue of 75 globular clusters' proper motions and recent measurements of the proper motions of another 20 distant clusters obtained with the Hubble Space Telescope. We describe the globular cluster system with a 2-component distribution function (DF), with a flat, rotating disc and a rounder, more extended halo. While fixing the Milky Way's disc and bulge, we let the mass and shape of the dark matter halo and we fit these two parameters, together with other six describing the DF, with a Bayesian method. We find the mass of the Galaxy within 20 kpc to be $M(<20{\,\rm kpc})=1.91^{+0.18}_{-0.17} \times 10^{11} M_\odot$, of which $M_{\rm DM}(<20{\,\rm kpc})=1.37^{+0.18}_{-0.17}\times 10^{11}M_\odot$ is in dark matter, and the density axis ratio of the dark matter halo to be $q=1.30 \pm 0.25$. This implies a virial mass $M_{\rm via} = 1.3 \pm 0.3 \times 10^{12} M_\odot$. Our analysis rules out oblate ($q<0.8$) and strongly prolate halos ($q>1.9$) with 99\% probability. Our preferred model reproduces well the observed phase-space distribution of globular clusters and has a disc component that closely resembles that of the Galactic thick disc. The halo component follows a power-law density profile $\rho \propto r^{-3.3}$, has a mean rotational velocity of $V_{\rm rot}\simeq -14\,\rm km\,s^{-1}$ at 20 kpc, and has a mildly radially biased velocity distribution ($\beta\simeq 0.2 \pm 0.07$, fairly constant outside 15 kpc). We also find that our distinction between disc and halo clusters resembles, although not fully, the observed distinction in metal-rich ([Fe/H]$>-0.8$) and metal-poor ([Fe/H]$\leq-0.8$) cluster populations.

We investigate the relation between stellar mass and specific stellar angular momentum, or `Fall relation', for a sample of 17 isolated, regularly rotating disc galaxies at z=1. All galaxies have a) rotation curves determined from Halpha emission-line data; b) HST imaging in optical and infrared filters; c) robust determinations of their stellar masses. We use HST images in f814w and f160w filters, roughly corresponding to rest-frames B and I bands, to extract surface brightness profiles for our systems. We robustly bracket the specific angular momentum by assuming that rotation curves beyond the outermost Halpha rotation point stay either flat or follow a Keplerian fall-off. By comparing our measurements with those determined for disc galaxies in the local Universe, we find no evolution in the Fall relation in the redshift range 0<z<1, regardless of the band used and despite the uncertainties in the stellar rotation curves at large radii. This result holds unless stellar masses at z=1 are systematically underestimated by more than 50%. Our findings are compatible with expectations based on a LCDM cosmological framework and support a scenario where both the stellar Tully-Fisher and mass-size relations for spirals do not evolve significantly in this redshift range.

We present new dynamical models of dwarf spheroidal galaxies (dSphs) in which both the stellar component and the dark halo are described by analytic distribution functions that depend on the action integrals. In their most general form these distribution functions can represent axisymmetric and possibly rotating stellar systems. Here, as a first application, we model the Fornax dSph, limiting ourselves, for simplicity, to the non rotating, spherical case. The models are compared with state-of-the-art spectroscopic and photometric observations of Fornax, exploiting the knowledge of the line-of-sight velocity distribution of the models and accounting for the foreground contamination from the Milky Way. The model that best fits the structural and kinematic properties of Fornax has a cored dark halo, with core size $r_{\rm c}\simeq1.03$ kpc. The dark-to-luminous mass ratio is $(M_{\rm dm}/M_{\star})|_{R_{\rm eff}}\simeq9.6$ within the effective radius $R_{\rm eff} \simeq 0.62\,$kpc and $(M_{\rm dm}/M_{\star})|_{3 {\rm kpc}} \simeq 144$ within 3 kpc. The stellar velocity distribution is isotropic almost over the full radial range covered by the spectroscopic data and slightly radially anisotropic in the outskirts of the stellar distribution. The dark-matter annihilation $J$-factor and decay $D$-factor are, respectively, $\log_{10}(J$ $[$GeV$^2$ cm$^{-5}])\simeq18.34$ and $\log_{10}(D$ $[$GeV cm$^{-2}])\simeq18.55$, for integration angle $\theta = 0.5^{\circ}$. This cored halo model of Fornax is preferred, with high statistical significance, to both models with a Navarro, Frenk and White dark halo and simple mass-follows-light models.

There is long tradition extending more than a century on the identification of moving groups in the Solar neighborhood. However, with the advent of large kinematic surveys, and especially of the upcoming Gaia data releases there is a need for more sophisticated and automated substructure finders. We analyze the TGASxRAVE dataset to identify moving groups in the Galactic disk. These groups of stars may then be used to map dynamical and star formation processes in the vicinity of the Sun. We use the ROCKSTAR algorithm, a "friends-of-friends"-like substructure finder in 6D phase-space, and analyze the Hertzsprung-Russell diagrams of the groups identified. We find 125 moving groups within 300 pc of the Sun, containing on average 50 stars, and with 3D velocity dispersions smaller than 10 km/s. Most of these groups were previously unknown. Our photometric analysis allows us to isolate a subsample of 30 statistically significant groups likely composed of stars that were born together.

Extended stellar haloes are a natural by-product of the hierarchical formation of massive galaxies. If merging is a non-negligible factor in the growth of our Galaxy, evidence of such events should be encoded in its stellar halo. Reliable identification of genuine halo stars is a challenging task however. The 1st Gaia data release contains the positions, parallaxes and proper motions for over 2 million stars, mostly in the Solar neighbourhood. Gaia DR2 will enlarge this sample to over 1.5 billion stars, the brightest ~5 million of which will have a full phase-space information. Our aim is to develop a machine learning model to reliably identify halo stars, even when their full phase-space information is not available. We use the Gradient Boosted Trees algorithm to build a supervised halo star classifier. The classifier is trained on a sample extracted from the Gaia Universe Model Snapshot, convolved with the errors of TGAS, as well as with the expected uncertainties of the upcoming Gaia DR2. We also trained our classifier on the cross-match between the TGAS and RAVE catalogues, where the halo stars are labelled in an entirely model independent way. We then use this model to identify halo stars in TGAS. When full phase- space information is available and for Gaia DR2-like uncertainties, our classifier is able to recover 90% of the halo stars with at most 30% distance errors, in a completely unseen test set, and with negligible levels of contamination. When line-of-sight velocity is not available, we recover ~60% of such halo stars, with less than 10% contamination. When applied to the TGAS data, our classifier detects 337 high confidence RGB halo stars. Although small, this number is consistent with the expectation from models given the data uncertainties. The large parallax errors are the biggest limitation to identify a larger number of halo stars in all the cases studied.

We study the dynamical properties of halo stars located in the Solar Neighbourhood. Our goal is to explore how the properties of the halo depend on the selection criteria used for defining a sample of halo stars. Once this is understood we proceed to measure the shape and orientation of the halo's velocity ellipsoid and we use this information to put constraints on the gravitational potential of the Galaxy. We use the recently released Gaia DR1 catalogue cross-matched to the RAVE dataset for our analysis. We develop a dynamical criterion based on the distribution function of stars in various Galactic components, using action integrals to identify halo members, and compare this to metallicity and to kinematically selected samples. With this new method, we find 1156 stars in the Solar Neighbourhood to be likely members of the stellar halo. Our dynamically selected sample consists mainly of distant giants on elongated orbits. Their metallicity distribution is rather broad, with roughly half of the stars having [M/H] $\ge -1$ dex. The use of different selection criteria has an important impact on the characteristics of the velocity distributions obtained. Nonetheless, for our dynamically selected and for the metallicity selected samples, we find the local velocity ellipsoid to be aligned in spherical coordinates in a Galactocentric reference frame: this suggests that the total gravitational potential is rather spherical in the region spanned by the orbits of the halo stars in these samples.

We derive the stellar-to-halo specific angular momentum relation (SHSAMR) of galaxies at $z=0$ by combining i) the standard $\Lambda$CDM tidal torque theory ii) the observed relation between stellar mass and specific angular momentum (Fall relation) and iii) various determinations of the stellar-to-halo mass relation (SHMR). We find that the ratio $f_j = j_\ast/j_{\rm h}$ of the specific angular momentum of stars to that of the dark matter i) varies with mass as a double power-law, ii) it always has a peak in the mass range explored and iii) it is $3-5$ times larger for spirals than for ellipticals. The results have some dependence on the adopted SHMR and we provide fitting formulae in each case. For any choice of the SHMR, the peak of $f_j$ occurs at the same mass where the stellar-to-halo mass ratio $f_\ast = M_\ast/M_{\rm h}$ has a maximum. This is mostly driven by the straightness and tightness of the Fall relation, which requires $f_j$ and $f_\ast$ to be correlated with each other roughly as $f_j\propto f_\ast^{2/3}$, as expected if the outer and more angular momentum rich parts of a halo failed to accrete onto the central galaxy and form stars (biased collapse). We also confirm that the difference in the angular momentum of spirals and ellipticals at a given mass is too large to be ascribed only to different spins of the parent dark-matter haloes (spin bias).

In a $\Lambda$CDM Universe, the specific stellar angular momentum ($j_\ast$) and stellar mass ($M_\ast$) of a galaxy are correlated as a consequence of the scaling existing for dark matter haloes ($j_{\rm h}\propto M_{\rm h}^{2/3}$). The shape of this law is crucial to test galaxy formation models, which are currently discrepant especially at the lowest masses, allowing to constrain fundamental parameters, e.g. the retained fraction of angular momentum. In this study, we accurately determine the empirical $j_\ast-M_\ast$ relation (Fall relation) for 92 nearby spiral galaxies (from S0 to Irr) selected from the Spitzer Photometry and Accurate Rotation Curves (SPARC) sample in the unprecedented mass range $7 \lesssim \log M_\ast/M_\odot \lesssim 11.5$. We significantly improve all previous estimates of the Fall relation by determining $j_\ast$ profiles homogeneously for all galaxies, using extended HI rotation curves, and selecting only galaxies for which a robust $j_\ast$ could be measured (converged $j_\ast(<R)$ radial profile). We find the relation to be well described by a single, unbroken power-law $j_\ast\propto M_\ast^\alpha$ over the entire mass range, with $\alpha=0.55\pm 0.02$ and orthogonal intrinsic scatter of $0.17\pm 0.01$ dex. We finally discuss some implications for galaxy formation models of this fundamental scaling law and, in particular, the fact that it excludes models in which discs of all masses retain the same fraction of the halo angular momentum.

The 3D motions of stars in small galaxies beyond our own are minute and yet they are crucial for our understanding of the nature of gravity and dark matter. Even for the dwarf galaxy Sculptor which is one of the best studied systems and inferred to be strongly dark matter dominated, there are conflicting reports on its mean motion around the Milky Way and the 3D internal motions of its stars have never been measured. Here we report, based on data from the Gaia space mission and the Hubble Space Telescope, a new precise measurement of Sculptor's mean proper motion. From this we deduce that Sculptor is currently at its closest approach to the Milky Way and moving on an elongated high-inclination orbit that takes it much farther away than previously thought. For the first time we are also able to measure the internal motions of stars in Sculptor. We find $\sigma_{R}=11.5 \pm 4.3$ km/s and $\sigma_{T}=8.5\pm3.2$ km/s along the projected radial and tangential directions, implying that the stars in our sample move preferentially on radial orbits as quantified by the anisotropy parameter, which we find to be $\beta\sim 0.86^{+0.12}_{-0.83}$ at a location beyond the core radius. Taken at face value such a high radial anisotropy requires abandoning conventional models for the mass distribution in Sculptor. Our sample is dominated by metal-rich stars and for these we find $\beta^{MR} \sim 0.95^{+0.04}_{-0.27}$, a value consistent with multi-component models where Sculptor is embedded in a cuspy dark halo as expected for cold dark matter.

Aims: We present the first measurement of the proper motion and orbit of the very distant and intriguing globular cluster NCG 2419. Methods: We have combined data from HST and Gaia DR1 to derive the relative proper motions of stars in the direction to the cluster. To tie to an absolute reference frame we have used a background galaxy located in the field. Results: We find the absolute proper motion of NGC 2419 to be $(\mu_{\alpha}\cos(\delta)$, $\mu_{\delta}$)=($-0.17\pm0.26,-0.49\pm0.17$) mas/yr. We have integrated the orbit of the cluster in a Galactic potential and found it to oscillate between $\sim$53 kpc and $\sim$98 kpc on a nearly polar orbit. This makes it very likely that NGC 2419 is a former cluster of the Sagittarius dwarf spheroidal galaxy, also because it shares the same sense of rotation around the Milky Way.

We present a new discrete chemo-dynamical axisymmetric modeling technique, which we apply to the dwarf spheroidal galaxy Sculptor. The major improvement over previous Jeans models is that realistic chemical distributions are included directly in the dynamical modelling of the discrete data. This avoids loss of information due to spatial binning and eliminates the need for hard cuts to remove contaminants and to separate stars based on their chemical properties. Using a combined likelihood in position, metallicity and kinematics, we find that our models naturally separate Sculptor stars into a metal-rich and a metal-poor population. Allowing for non-spherical symmetry, our approach provides a central slope of the dark matter density of $\gamma = 0.5 \pm 0.3$. The metal-rich population is nearly isotropic (with $\beta_r^{red} = 0.0\pm0.1$) while the metal-poor population is tangentially anisotropic (with $\beta_r^{blue} = -0.2\pm0.1$) around the half light radius of $0.26$ kpc. A weak internal rotation of the metal-rich population is revealed with $v_{max}/\sigma_0 = 0.15 \pm 0.15$. We run tests using mock data to show that a discrete dataset with $\sim 6000$ stars is required to distinguish between a core ($\gamma = 0$) and cusp ($\gamma = 1$), and to constrain the possible internal rotation to better than $1\,\sigma$ confidence with our model. We conclude that our discrete chemo-dynamical modelling technique provides a flexible and powerful tool to robustly constrain the internal dynamics of multiple populations, and the total mass distribution in a stellar system.

Clusters of galaxies are embedded in halos of optically thin, gravitationally stratified, weakly magnetized plasma at the system's virial temperature. Due to radiative cooling and anisotropic heat conduction, such intracluster medium (ICM) is subject to local instabilities, which are combinations of the thermal, magnetothermal and heat-flux-driven buoyancy instabilities. If the ICM rotates significantly, its stability properties are substantially modified and, in particular, also the magnetorotational instability (MRI) can play an important role. We study simple models of rotating cool-core clusters and we demonstrate that the MRI can be the dominant instability over significant portions of the clusters, with possible implications for the dynamics and evolution of the cool cores. Our results give further motivation for measuring the rotation of the ICM with future X-ray missions such as ASTRO-H and ATHENA.

We present an approach to the design of distribution functions that depend on the phase-space coordinates through the action integrals. The approach makes it easy to construct a dynamical model of a given stellar component. We illustrate the approach by deriving distribution functions that self-consistently generate several popular stellar systems, including the Hernquist, Jaffe, and Navarro, Frenk and White models. We focus on non-rotating spherical systems, but extension to flattened and rotating systems is trivial. Our distribution functions are easily added to each other and to previously published distribution functions for discs to create self-consistent multi-component galaxies. The models this approach makes possible should prove valuable both for the interpretation of observational data and for exploring the non-equilibrium dynamics of galaxies via N-body simulation.

A long-standing question is whether radiative cooling can lead to local condensations of cold gas in the hot atmospheres of galaxies and galaxy clusters. We address this problem by studying the nature of local instabilities in rotating, stratified, weakly magnetized, optically thin plasmas in the presence of radiative cooling and anisotropic thermal conduction. For both axisymmetric and non-axisymmetric linear perturbations we provide the general equations that can be applied locally to specific systems to establish whether they are unstable and, in case of instability, to determine the kind of evolution (monotonically growing or over-stable) and the growth rates of unstable modes. We present results for models of rotating plasmas representative of Milky Way-like galaxy coronae and cool-cores of galaxy clusters. It is shown that the unstable modes arise from a combination of thermal, magnetothermal, magnetorotational and heat-flux-driven buoyancy instabilities. Local condensation of cold clouds tends to be hampered in cluster cool cores, while it is possible under certain conditions in rotating galactic coronae. If the magnetic field is sufficiently weak the magnetorotational instability is dominant even in these pressure supported systems.

Early-type galaxies (ETGs) are observed to be more compact, on average, at $z \gtrsim 2$ than at $z\simeq 0$, at fixed stellar mass. Recent observational works suggest that such size evolution could reflect the similar evolution of the host dark matter halo density as a function of the time of galaxy quenching. We explore this hypothesis by studying the distribution of halo central velocity dispersion ($\sigma_0$) and half-mass radius ($r_{\rm h}$) as functions of halo mass $M$ and redshift $z$, in a cosmological $\Lambda$-CDM $N$-body simulation. In the range $0\lesssim z\lesssim 2.5$, we find $\sigma_0\propto M^{0.31-0.37}$ and $r_{\rm h}\propto M^{0.28-0.32}$, close to the values expected for homologous virialized systems. At fixed $M$ in the range $10^{11} M_\odot \lesssim M\lesssim 5.5 \times 10^{14} M_\odot$ we find $\sigma_0\propto(1+z)^{0.35}$ and $r_{\rm h}\propto(1+z)^{-0.7}$. We show that such evolution of the halo scaling laws is driven by individual haloes growing in mass following the evolutionary tracks $\sigma_0\propto M^{0.2}$ and $r_{\rm h}\propto M^{0.6}$, consistent with simple dissipationless merging models in which the encounter orbital energy is accounted for. We compare the $N$-body data with ETGs observed at $0\lesssim z\lesssim3$ by populating the haloes with a stellar component under simple but justified assumptions: the resulting galaxies evolve consistently with the observed ETGs up to $z \simeq 2$, but the model has difficulty reproducing the fast evolution observed at $z\gtrsim 2$. We conclude that a substantial fraction of the size evolution of ETGs can be ascribed to a systematic dependence on redshift of the dark matter haloes structural properties.

The thermal stability of a weakly magnetized, rotating, stratified, optically thin plasma is studied by means of linear-perturbation analysis. We derive dispersion relations and criteria for stability against axisymmetric perturbations that generalize previous results on either non-rotating or unmagnetized fluids. The implications for the hot atmospheres of galaxies and galaxy clusters are discussed.