Search

Observed scaling relations in galaxies between baryons and dark matter global properties are key to shed light on the process of galaxy formation and on the nature of dark matter. Here, we study the scaling relation between the neutral hydrogen (HI) and dark matter mass in isolated rotationally-supported disk galaxies at low redshift. We first show that state-of-the-art galaxy formation simulations predict that the HI-to-dark halo mass ratio decreases with stellar mass for the most massive disk galaxies. We then infer dark matter halo masses from high-quality rotation curve data for isolated disk galaxies in the local Universe, and report on the actual universality of the HI-to-dark halo mass ratio for these observed galaxies. This scaling relation holds for disks spanning a range of 4 orders of magnitude in stellar mass and 3 orders of magnitude in surface brightness. Accounting for the diversity of rotation curve shapes in our observational fits decreases the scatter of the HI-to-dark halo mass ratio while keeping it constant. This finding extends the previously reported discrepancy for the stellar-to-halo mass relation of massive disk galaxies within galaxy formation simulations to the realm of neutral atomic gas. Our result reveals that isolated galaxies with regularly rotating extended HI disks are surprisingly self-similar up to high masses, which hints at mass-independent self-regulation mechanisms that have yet to be fully understood., Comment: 14 pages, 4 figures. Accepted for publication in ApJL

Thin stellar discs on both galactic and nuclear, sub-kpc scales are believed to be fragile structures that would be easily destroyed in major mergers. In turn, this makes the age-dating of their stellar populations a useful diagnostics for the assembly history of galaxies. We aim at carefully exploring the fragility of such stellar discs in intermediate- and low- mass encounters, using high-resolution N-body simulations of galaxy models with structural and kinematic properties tailored to actually observed galaxies. As a first but challenging step, we create a dynamical model of FCC 170, a nearly edge-on galaxy in the Fornax cluster with multiple galactic components and including both a galactic scale and nuclear stellar disc (NSD), using detailed kinematic data from the Multi Unit Spectroscopic Explorer and a novel method for constructing distribution function-based self-consistent galaxy models. We then create N-body realisations of this model and demonstrate that it remains in equilibrium and preserves its properties over many Gyr, when evolved with a sufficiently high particle number. However, the NSD is more prone to numerical heating, which gradually increases its thickness by up to 22 per cent in 10 Gyr even in our highest-resolution runs. Nevertheless, these N-body models can serve as realistic representations of actual galaxies in merger simulations., Comment: 13 pages, 9 figures

Correlations between velocity measurements in disk galaxy rotation curves are usually neglected when fitting dynamical models. Here I show how data correlations can be taken into account in rotation curve decompositions using Gaussian Processes. I find that marginalizing over correlation parameters proves critical to obtain unbiased estimates of the luminous and dark matter distributions in galaxies., Comment: 4 pages, 1 figure. Published in RNAAS. Associated jupyter notebook at https://lposti.github.io/MLPages/gaussian_processes/2022/11/02/gp_rotcurves.html

We present new and archival atomic hydrogen (\hi) observations of \galnum\ of the most massive spiral galaxies in the local Universe ($M_\star>10^{11} \, \mathrm{M}_\odot$). From 3D kinematic modeling of the datacubes, we derive extended \hi\ rotation curves, and from these, we estimate masses of the dark matter halos and specific angular momenta of the discs. We confirm that massive spiral galaxies lie at the upper ends of the Tully-Fisher relation (mass vs velocity, $M \propto V^{4}$) and Fall relation (specific angular momentum vs mass, $j \propto M^{0.6}$), in both stellar and baryonic forms, with no significant deviations from single power laws. We study the connections between baryons and dark matter through the stellar (and baryon)-to-halo ratios of mass $f_\mathrm{M} \equiv M_\star/M_\mathrm{h}$ and specific angular momentum $f_\mathrm{j} \equiv j_\star/j_\mathrm{h}$ and $f_\mathrm{j,bar} \equiv j_\mathrm{bar}/j_\mathrm{h}$. Combining our sample with others from the literature for less massive disc-dominated galaxies, we find that $f_\mathrm{M}$ rises monotonically with $M_\star$ and $M_\mathrm{h}$ (instead of the inverted-U shaped $f_\mathrm{M}$ for spheroid-dominated galaxies), while $f_\mathrm{j}$ and $f_\mathrm{j,bar}$ are essentially constant near unity over four decades in mass. Our results indicate that disc galaxies constitute a self-similar population of objects closely linked to the self-similarity of their dark halos. This picture is reminiscent of early analytical models of galaxy formation wherein discs grow by relatively smooth and gradual inflow, isolated from disruptive events such as major mergers and strong AGN feedback, in contrast to the more chaotic growth of spheroids., Comment: Accepted for publication on MNRAS

We study the kinematics and scaling relations of a sample of 43 giant spiral galaxies that have stellar masses exceeding $10^{11}$ $M_\odot$ and optical discs up to 80 kpc in radius. We use a hybrid 3D-1D approach to fit 3D kinematic models to long-slit observations of the H$\alpha$-[NII] emission lines and we obtain robust rotation curves of these massive systems. We find that all galaxies in our sample seem to reach a flat part of the rotation curve within the outermost optical radius. We use the derived kinematics to study the high-mass end of the two most important scaling relations for spiral galaxies: the stellar/baryonic mass Tully-Fisher relation and the Fall (mass-angular momentum) relation. All galaxies in our sample, with the possible exception of the two fastest rotators, lie comfortably on both these scaling relations determined at lower masses, without any evident break or bend at the high-mass regime. When we combine our high-mass sample with lower-mass data from the Spitzer Photometry & Accurate Rotation Curves catalog, we find a slope of $\alpha=4.25\pm0.19$ for the stellar Tully-Fisher relation and a slope of $\gamma=0.64\pm0.11$ for the Fall relation. Our results indicate that most, if not all, of these rare, giant spiral galaxies are scaled up versions of less massive discs and that spiral galaxies are a self-similar population of objects up to the very high-mass end., Comment: 12 pages, 5 figures + Appendix. Accepted for publication in MNRAS

The relations between the specific angular momenta ($j$) and masses ($M$) of galaxies are often used as a benchmark in analytic models and hydrodynamical simulations as they are considered to be amongst the most fundamental scaling relations. Using accurate measurements of the stellar ($j_\ast$), gas ($j_{\rm gas}$), and baryonic ($j_{\rm bar}$) specific angular momenta for a large sample of disc galaxies, we report the discovery of tight correlations between $j$, $M$, and the cold gas fraction of the interstellar medium ($f_{\rm gas}$). At fixed $f_{\rm gas}$, galaxies follow parallel power laws in 2D $(j,M)$ spaces, with gas-rich galaxies having a larger $j_\ast$ and $j_{\rm bar}$ (but a lower $j_{\rm gas}$) than gas-poor ones. The slopes of the relations have a value around 0.7. These new relations are amongst the tightest known scaling laws for galaxies. In particular, the baryonic relation ($j_{\rm bar}-M_{\rm bar}-f_{\rm gas}$), arguably the most fundamental of the three, is followed not only by typical discs but also by galaxies with extreme properties, such as size and gas content, and by galaxies previously claimed to be outliers of the standard 2D $j-M$ relations. The stellar relation ($j_{\ast}-M_{\ast}-f_{\rm gas}$) may be connected to the known $j_\ast-M_\ast-$bulge fraction relation; however, we argue that the $j_{\rm bar}-M_{\rm bar}-f_{\rm gas}$ relation can originate from the radial variation in the star formation efficiency in galaxies, although it is not explained by current disc instability models., Comment: A&A Letters, in press. Data catalogue will be available via CDS and at this link https://unishare.nl/index.php/s/NMQRYfrrDpj8iaj

Although true metal-free "Population III" stars have so-far escaped discovery, their nature, and that of their supernovae, is revealed in the chemical products left behind in the next generations of stars. Here we report the detection of an ultra-metal poor star in the Sculptor dwarf spheroidal galaxy, AS0039. With [Fe/H]$_{\rm LTE}=-4.11$, it is the most metal-poor star so far discovered in any external galaxy. Contrary to the majority of Milky Way stars at this metallicity, AS0039 is clearly not enhanced in carbon, with [C/Fe]$_{\rm LTE}=-0.75$ and A(C)=+3.60, making it the lowest detected carbon abundance in any star to date. It furthermore lacks $\alpha$-element uniformity, having extremely low [Mg/Ca]$_{\rm NLTE}=-0.60$ and [Mg/Ti]$_{\rm NLTE}=-0.86$, in stark contrast with the near solar ratios observed in C-normal stars within the Milky Way halo. The unique abundance pattern indicates that AS0039 formed out of material that was predominantly enriched by a $\sim$20$ M_\odot$ progenitor star with an unusually high explosion energy $E=10\times10^{51}$ erg. The star AS0039 is thus one of the first observational evidence for zero-metallicity hypernovae and provides a unique opportunity to investigate the diverse nature of Population III stars., Comment: Accepted in ApJL, 7 pages, 3 figures, 1 table, and an appendix (1 figure and 3 tables)

We derive the stellar-to-halo mass relation (SHMR), namely $f_\star\propto M_\star/M_{\rm h}$ versus $M_\star$ and $M_{\rm h}$, for early-type galaxies from their near-IR luminosities (for $M_\star$) and the position-velocity distributions of their globular cluster systems (for $M_{\rm h}$). Our individual estimates of $M_{\rm h}$ are based on fitting a dynamical model with a distribution function expressed in terms of action-angle variables and imposing a prior on $M_{\rm h}$ from the concentration-mass relation in the standard $\Lambda$CDM cosmology. We find that the SHMR for early-type galaxies declines with mass beyond a peak at $M_\star\sim 5\times 10^{10}M_\odot$ and $M_{\rm h}\sim 10^{12}M_\odot$ (near the mass of the Milky Way). This result is consistent with the standard SHMR derived by abundance matching for the general population of galaxies, and with previous, less robust derivations of the SHMR for early types. However, it contrasts sharply with the monotonically rising SHMR for late types derived from extended HI rotation curves and the same $\Lambda$CDM prior on $M_{\rm h}$ as we adopt for early types. The SHMR for massive galaxies varies more or less continuously, from rising to falling, with decreasing disc fraction and decreasing Hubble type. We also show that the different SHMRs for late and early types are consistent with the similar scaling relations between their stellar velocities and masses (Tully-Fisher and Faber-Jackson relations). Differences in the relations between the stellar and halo virial velocities account for the similarity of the scaling relations. We argue that all these empirical findings are natural consequences of a picture in which galactic discs are built mainly by smooth and gradual inflow, regulated by feedback from young stars, while galactic spheroids are built by a cooperation between merging, black-hole fuelling, and feedback from AGNs., Comment: Accepted for publication in A&A. Figure 3 shows the main result of the paper

Ultra-diffuse galaxies (UDGs) are very low-surface brightness galaxies with large effective radii. Spectroscopic measurements of a few UDGs have revealed a low dark matter content, based on the internal motion of stars or globular clusters (GCs). This is in contrast to the large number of GCs found for these systems, from which it would be expected to correspond to a large dark matter halo mass. Here we present HST+ACS observations for the UDG MATLAS-2019 in the NGC5846 group. Using the F606W and F814W filters, we trace the GC population two magnitudes below the peak of the GC luminosity function (GCLF). Employing Bayesian considerations, we identify 26+-6 GCs associated with the dwarf, yielding a large specific frequency of S_N=58+-14. We use the turnover of the GCLF to derive a distance of 21+-2 Mpc, which is consistent with the NGC5846 group of galaxies. Due to the superior image quality of the HST, we are able to resolve the GCs and measure their sizes, which are consistent with the sizes of GCs around Local Group galaxies. Using the linear relation between the total mass of galaxies and of GCs, we derive a halo mass of 0.9(+-0.2) *10^11 M_solar (M_solar/L_solar>1000). The high abundance of GCs, together with the small uncertainties, make MATLAS-2019 one of the most extreme UDGs, which likely sets an upper limit of the number of GCs for UDGs., Comment: 16 pages, 10 figures, 1 table, accepted for publication in APJ

Galaxy scaling laws, such as the Tully-Fisher, mass-size and Fall relations, can provide extremely useful clues on our understanding of galaxy formation in a cosmological context. Some of these relations are extremely tight and well described by one single parameter (mass), despite the theoretical existence of secondary parameters such as spin and concentration, which are believed to impact these relations. In fact, the residuals of these scaling laws appear to be almost uncorrelated with each other, posing significant constraints on models where secondary parameters play an important role. Here, we show that a possible solution is that such secondary parameters are correlated amongst themselves, in a way that removes correlations in observable space. In particular, we focus on how the existence of an anti-correlation between the dark matter halo spin and its concentration -- which is still debated in simulations -- can weaken the correlation of the residuals of the Tully-Fisher and mass-size relations. Interestingly, using simple analytic galaxy formation models, we find that this happens only for a relatively small portion of the parameter space that we explored, which suggests that this idea could be used to derive constraints to galaxy formation models that are still unexplored., Comment: Accepted for publication in A&A

(Abridged) Specific angular momentum is one of the key parameters that control the evolution of galaxies. We derive the baryonic specific angular momentum of disc galaxies and study its relation with the dark matter specific angular momentum. Using a combination of high-quality HI rotation curves and HI/near-IR surface densities, we homogeneously measure the stellar ($j_{\rm *}$) and gas ($j_{\rm gas}$) specific angular momenta for a large sample of local disc galaxies. This allows us to determine the baryonic specific angular momentum ($j_{\rm bar}$) with high accuracy and across a very wide range of masses. The $j_{\ast}-M_\ast$ relation is an unbroken power-law from $7 \lesssim$ log($M_\ast$/$M_\odot) \lesssim 11.5$, with slope $0.54 \pm 0.02$. For the gas component, we find that the $j_{\rm gas}-M_{\rm gas}$ relation is also an unbroken power-law from $6 \lesssim$ log($M_{\rm gas}$/$M_\odot) \lesssim 11$, with a steeper slope of $1.02 \pm 0.04$. Regarding the baryonic relation, our data support a correlation characterized by single power-law with slope $0.60 \pm 0.02$. Our most massive spirals and smallest dwarfs lie along the same $j_{\rm bar}-M_{\rm bar}$ sequence. While the relations are tight and unbroken, we find internal correlations inside them: At fixed $M_\ast$, galaxies with larger $j_\ast$ have larger disc scale lengths, and at fixed $M_{\rm bar}$, gas-poor galaxies have lower $j_{\rm bar}$ than expected. We estimate the retained fraction of baryonic specific angular momentum, finding it constant across our entire mass range with a value of $\sim 0.6$, indicating that the $j_{\rm bar}$ of present-day disc galaxies is comparable to the initial specific angular momentum of their dark matter haloes. These results set important constraints for hydrodynamical simulations and semi-analytical models aiming to reproduce galaxies with realistic specific angular momenta., Comment: Matching the published version of January (2021A&A...647A..76M) The full version of the main table is available at https://www.dropbox.com/s/xizhjqiuahqbdtp/tbsam.dat

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., Comment: Accepted for publication in MNRAS. v2: a few typos have been corrected and a couple of references added

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., Comment: A&A in press (17 pages, 3 appendices)

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., Comment: Accepted for publication by The Astrophysical Journal Letters (ApJL). V2: Acknowledgments have been updated, and a typo has been corrected

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., Comment: 5 pages, 3 figures, 1 table. Accepted for publication in A&A Letters

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., Comment: submitted to MNRAS

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., Comment: 15 pages, 13 figures, A&A

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$., Comment: accepted for publication in A&A

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., Comment: Accepted for publication in A&A

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., Comment: 6 pages, 5 figures, Submitted to A&A

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., Comment: Accepted for publication in Astronomy & Astrophysics. 13 pages, 9 figure, 2 tables

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(

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., Comment: Accepted 2018 July 6; Received 2018 June 1; Submitted in original to MNRAS form 2018 February 5

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., Comment: matches A&A published version

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)., Comment: matches MNRAS published version

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., Comment: 5 pages, 4 figures, accepted for publication by A&A Letters

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., Comment: Accepted by MNRAS

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., Comment: 9 pages, 8 figures, 1 appendix, MNRAS accepted

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., Comment: 16 pages, 8 figures, ApJ in press. Minor changes from previous version

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., Comment: 15 pages, 14 figures, 1 table. Matches the Accepted version from MNRAS

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., Comment: 16 pages, 3 figures, MNRAS accepted. New figures and corrected equations with respect to previous version. Results unchanged