Showing papers by "Julio F. Navarro published in 2021"

The asymptotic tidal remnants of cold dark matter subhaloes

Abstract: We use N-body simulations to study the evolution of cuspy cold dark matter (CDM) halos in the gravitational potential of a massive host. Tidal mass losses reshape CDM halos, leaving behind bound remnants whose characteristic densities are set by the mean density of the host at the pericentre of their respective orbit. The evolution to the final bound remnant state is essentially complete after ~5 orbits for nearly circular orbits, while reaching the same remnant requires ~25 and ~40 orbits for eccentric orbits with 1:5 and 1:20 pericentre-to-apocentre ratios, respectively. The density profile of tidal remnants is fully specified by the fraction of mass lost, and approaches an exponentially-truncated Navarro-Frenk-White profile in the case of heavy mass loss. Resolving tidal remnants requires excellent numerical resolution; poorly resolved subhalos have systematically lower characteristic densities and are more easily disrupted. Even simulations with excellent spatial and time resolution fail when the final remnant is resolved with fewer than 3000 particles. We derive a simple empirical model that describes the evolution of the mass and the density profile of the tidal remnant applicable to a wide range of orbital eccentricities and pericentric distances. Applied to the Milky Way, our results suggest that $10^8$ - $10^{10}$ solar mass halos accreted ~10 Gyrs ago on 1:10 orbits with pericentric distance ~10 kpc should have been stripped to 0.1 - 1 per cent of their original mass. This implies that estimates of the survival and structure of such halos (the possible hosts of ultra-faint Milky Way satellites) based on direct cosmological simulations may be subject to substantial revision.

Revisiting the tension between fast bars and the ΛCDM paradigm

Abstract: The pattern speed with which galactic bars rotate is intimately linked to the amount of dark matter in the inner regions of their host galaxies. In particular, dark matter haloes act to slow down bars via torques exerted through dynamical friction. Observational studies of barred galaxies tend to find that bars rotate fast, while hydrodynamical cosmological simulations of galaxy formation and evolution in the Lambda cold dark matter (ΛCDM) framework have previously found that bars slow down excessively. This has led to a growing tension between fast bars and the ΛCDM cosmological paradigm. In this study we revisit this issue, using the Auriga suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations of galaxy formation and evolution in the ΛCDM framework, finding that bars remain fast down to z = 0. In Auriga, bars form in galaxies that have higher stellar-to-dark matter ratios and are more baryon-dominated than in previous cosmological simulations; this suggests that in order for bars to remain fast, massive spiral galaxies must lie above the commonly used abundance matching relation. While this reduces the aforementioned tension between the rotation speed of bars and ΛCDM, it accentuates the recently reported discrepancy between the dynamically inferred stellar-to-dark matter ratios of massive spirals and those inferred from abundance matching. Our results highlight the potential of using bar dynamics to constrain models of galaxy formation and evolution.Key words: galaxies: evolution / galaxies: kinematics and dynamics / galaxies: spiral / galaxies: general / dark matter

The Pristine survey – XII. Gemini-GRACES chemo-dynamical study of newly discovered extremely metal-poor stars in the Galaxy

Abstract: High-resolution optical spectra of 30 metal-poor stars selected from the Pristine survey are presented, based on observations taken with the Gemini Observatory GRACES spectrograph. Stellar parameters Teff and log g are determined using a Gaia DR2 colour–temperature calibration and surface gravity from the Stefan–Boltzmann equation. GRACES spectra are used to determine chemical abundances (or upper limits) for 20 elements (Li, O, Na, Mg, K, Ca, Ti, Sc, Cr, Mn, Fe, Ni, Cu, Zn, Y, Zr, Ba, La, Nd, Eu). These stars are confirmed to be metal-poor ([Fe/H] < −2.5), with higher precision than from earlier medium-resolution analyses. The chemistry for most targets is similar to other extremely metal-poor stars in the Galactic halo. Three stars near [Fe/H] = −3.0 have unusually low Ca and high Mg, suggestive of contributions from few SN II where alpha-element formation through hydrostatic nucleosynthesis was more efficient. Three new carbon-enhanced metal-poor (CEMP) stars are also identified (two CEMP-s and one potential CEMP-no star) when our chemical abundances are combined with carbon from previous medium-resolution analyses. The GRACES spectra also provide precision radial velocities (σRV ≤ 0.2 km s−1) for dynamical orbit calculations with the Gaia DR2 proper motions. Most of our targets are dynamically associated with the Galactic halo; however, five stars with [Fe/H] < −3 have planar-like orbits, including one retrograde star. Another five stars are dynamically consistent with the Gaia-Sequoia accretion event; three have typical halo [α/Fe] ratios for their metallicities, whereas two are [Mg/Fe]-deficient, and one is a new CEMP-s candidate. These results are discussed in terms of the formation and early chemical evolution of the Galaxy.

Magellanic satellites in ΛCDM cosmological hydrodynamical simulations of the Local Group

Abstract: We use the APOSTLE ΛCDM cosmological hydrodynamical simulations of the Local Group to study the recent accretion of massive satellites into the halo of Milky Way (MW)-sized galaxies. These systems are selected to be close analogues to the Large Magellanic Cloud (LMC), the most massive satellite of the MW. The simulations allow us to address, in a cosmological context, the impact of the Clouds on the MW, including the contribution of Magellanic satellites to the MW satellite population, and the constraints placed on the Galactic potential by the motion of the LMC. We show that LMC-like satellites are twice more common around Local Group-like primaries than around isolated haloes of similar mass; these satellites come from large turnaround radii and are on highly eccentric orbits whose velocities at first pericentre are comparable with the primary’s escape velocity. This implies VMWesc(50 kpc) ∼ 365 km s−1, a strong constraint on Galactic potential models. LMC analogues contribute about two satellites with M∗>105M⊙⁠, having thus only a mild impact on the luminous satellite population of their hosts. At first pericentre, LMC-associated satellites are close to the LMC in position and velocity, and are distributed along the LMC’s orbital plane. Their orbital angular momenta roughly align with the LMC’s, but, interestingly, they may appear to ‘counter-rotate’ the MW in some cases. These criteria refine earlier estimates of the LMC association of MW satellites: only the SMC, Hydrus1, Car3, Hor1, Tuc4, Ret2, and Phoenix2 are compatible with all criteria. Carina, Grus2, Hor2, and Fornax are less probable associates given their large LMC relative velocity.

The star formation history of Eridanus II: on the role of SNe feedback in the quenching of ultra-faint dwarf galaxies

Abstract: Eridanus II (EriII) is an ultra-faint dwarf (UFD) galaxy (M_V=-7.1) located at a distance close to the Milky Way virial radius. Early shallow color-magnitude diagrams (CMD) indicated that it possibly hosted an intermediate-age or even young stellar population, which is unusual for a galaxy of this mass. In this paper, we present new ACS/HST CMDs reaching the oldest main sequence turnoff with excellent photometric precision, and derive a precise star formation history (SFH) for this galaxy through CMD-fitting. This SFH shows that the bulk of the stellar mass in Eri II formed in an extremely short star formation burst at the earliest possible time. The derived star formation rate profile has a width at half maximum of 500 Myr and reaches a value compatible with null star formation 13 Gyr ago. However, tests with mock stellar populations and with the CMD of the globular cluster M92 indicate that the star formation period could be shorter than 100 Myr. From the quantitative determination of the amount of mass turned into stars in this early star formation burst (~2x10^5 Msun) we infer the number of SNe events and the corresponding energy injected into the interstellar medium. For reasonable estimates of the EriII virial mass and values of the coupling efficiency of the SNe energy, we conclude that EriII could be quenched by SNe feedback alone, thus casting doubts on the need to invoke cosmic reionization as the preferred explanation for the early quenching of old UFD galaxies.

The Pristine Inner Galaxy Survey (PIGS) III: carbon-enhanced metal-poor stars in the bulge

Abstract: The most metal-deficient stars hold important clues about the early build-up and chemical evolution of the Milky Way, and carbon-enhanced metal-poor (CEMP) stars are of special interest. However, little is known about CEMP stars in the Galactic bulge. In this paper, we use the large spectroscopic sample of metal-poor stars from the Pristine Inner Galaxy Survey (PIGS) to identify CEMP stars ([C/Fe] > +0.7) in the bulge region and to derive a CEMP fraction. We identify 96 new CEMP stars in the inner Galaxy, of which 62 are very metal-poor ([Fe/H] < -2.0); this is more than a ten-fold increase compared to the seven previously known bulge CEMP stars. The cumulative fraction of CEMP stars in PIGS is $42^{\,+14\,}_{\,-13} \%$ for stars with [Fe/H] < -3.0, and decreases to $16^{\,+3\,}_{\,-3} \%$ for [Fe/H] < -2.5 and $5.7^{\,+0.6\,}_{\,-0.5} \%$ for [Fe/H] < -2.0. The PIGS inner Galaxy CEMP fraction for [Fe/H] < -3.0 is consistent with the halo fraction found in the literature, but at higher metallicities the PIGS fraction is substantially lower. While this can partly be attributed to a photometric selection bias, such bias is unlikely to fully explain the low CEMP fraction at higher metallicities. Considering the typical carbon excesses and metallicity ranges for halo CEMP-s and CEMP-no stars, our results point to a possible deficiency of both CEMP-s and CEMP-no stars (especially the more metal-rich) in the inner Galaxy. The former is potentially related to a difference in the binary fraction, whereas the latter may be the result of a fast chemical enrichment in the early building blocks of the inner Galaxy.

The Pristine Dwarf-Galaxy survey - III. Revealing the nature of the Milky Way globular cluster Sagittarius II

Abstract: We present a new spectroscopic study of the faint Milky Way satellite Sagittarius II. Using multi-object spectroscopy from the Fibre LargeArray Multi Element Spectrograph, we supplement the dataset of Longeard et al. (2020) with 47 newly observed stars, 19 of which are identified as members of the satellite. These additional member stars are used toput tighter constraints on the dynamics and the metallicity properties of the system. We find a low velocity dispersion of SgrII v = 1.7 +/-0.5 km s-1, in agreement with the dispersion of Milky Way globular clusters of similar luminosity. We confirm the very metal-poor nature of the satellite ([Fe/H]_SgrII = -2.23 +/- 0.07) and find that the metallicity dispersion of Sgr II is not resolved, reaching only 0.20 atthe 95% confidence limit. No star with a metallicity below -2.5 is confidently detected. Therefore, despite the unusually large size of the system (rh = 35.5 +1.4-1.2 pc), we conclude that Sgr II is an old and metal-poor globular cluster of the Milky Way.

Satellites around Milky Way Analogs: Tension in the Number and Fraction of Quiescent Satellites Seen in Observations versus Simulations

Abstract: We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage II release of the Satellites Around Galactic Analogs Survey (SAGA-ii) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star formation indicator for the SAGA-ii sample and derive star formation rates (SFRs) to compare with those from APOSTLE and Auriga. We compare our detection rates from the NUV and FUV bands to the SAGA-ii Hα detections and find that they are broadly consistent with over 85% of observed satellites detected in all three tracers. We apply the same spatial selection criteria used around SAGA-ii hosts to select satellites around the MW-like hosts in APOSTLE and Auriga. We find very good overall agreement in the derived SFRs for the star-forming satellites as well as the number of star-forming satellites per host in observed and simulated samples. However, the number and fraction of quenched satellites in the SAGA-ii sample are significantly lower than those in APOSTLE and Auriga below a stellar mass of M* ∼ 108 M⊙, even when the SAGA-ii incompleteness and interloper corrections are included. This discrepancy is robust with respect to the resolution of the simulations and persists when alternative star formation tracers are employed. We posit that this disagreement is not readily explained by vagaries in the observed or simulated samples considered here, suggesting a genuine discrepancy that may inform the physics of satellite populations around MW analogs.

Globular clusters as tracers of the dark matter content of dwarfs in galaxy clusters

Abstract: Globular clusters (GCs) are often used to estimate the dark matter content of galaxies, especially dwarfs, where other kinematic tracers are lacking. These estimates typically assume spherical symmetry and dynamical equilibrium, assumptions that may not hold for the sparse GC population of dwarfs in galaxy clusters. We use a catalog of GCs tagged onto the Illustris simulation to study the accuracy of GC-based mass estimates. We focus on galaxies in the stellar mass range 10$^{8} - 10^{11.8}$ M$_{\odot}$ identified in $9$ simulated Virgo-like clusters. Our results indicate that mass estimates are, on average, quite accurate in systems with GC numbers $N_{\rm GC} \geq 10$ and where the uncertainty of individual GC line-of-sight velocities is smaller than the inferred velocity dispersion, $\sigma_{\rm GC}$. In cases where $N_{\rm GC} \leq 10$, however, biases may result depending on how $\sigma_{\rm GC}$ is computed. We provide calibrations that may help alleviate these biases in methods widely used in the literature. As an application, we find a number of dwarfs with $M_{*} \sim 10^{8.5}\, M_{\odot}$ (comparable to the ultradiffuse galaxy DF2, notable for the low $\sigma_{GC}$ of its $10$ GCs) with $\sigma_{\rm GC} \sim 7$ - $15\; \rm km \rm s^{-1}$. These DF2 analogs correspond to relatively massive systems at their infall time ($M_{200} \sim 1$ - $3 \times 10^{11}$ $M_{\odot}$) which have retained only $3$-$17$ GCs and have been stripped of more than 95$\%$ of their dark matter. Our results suggest that extreme tidal mass loss in otherwise normal dwarf galaxies may be a possible formation channel for ultradiffuse objects like DF2.

Pericentric passage-driven star formation in satellite galaxies and their hosts: CLUES from local group simulations

Abstract: Local Group satellite galaxies show a wide diversity of star formation histories (SFHs) whose origin is yet to be fully understood. Using hydrodynamical simulations from the Constrained Local UniversE project, we study the SFHs of satellites of Milky Way-like galaxies in a cosmological context: while in the majority of the cases the accretion onto their host galaxy causes the satellites to lose their gas, with a subsequent suppression in SF, in about 25$\%$ of our sample we observe a clear enhancement of SF after infall. Peaks in SF clearly correlate with the satellite pericentric passage around its host and, in one case, with a satellite-satellite interaction. We identify two key ingredients that result in enhanced SF after infall: galaxies must enter the host's virial radius with a reservoir of cold gas $M_{\rm gas,inf}/M_{\rm vir,inf}\gtrsim 10^{-2}$ and with a minimum pericentric distance $\gtrsim$10 kpc (mean distance $\sim$50 kpc for the full sample), in order to form new stars due to compression of cold gas at pericentric passage. On the other hand, satellites that infall with little gas or whose pericentric distance is too small, have their gas ram-pressure stripped and subsequent SF quenched. The pericentric passage of satellites likewise correlates with SF peaks in their hosts, suggesting that this mechanism induces bursts of SF in satellites and central galaxies alike, in agreement with recent studies of our Galaxy's SFH. Our findings can explain the recently reported multiple stellar populations observed in dwarf galaxies such as Carina and Fornax, and should be taken into account in semi-analytic models of galaxy formation and satellite quenching.

Velocity-dependent J-factors for annihilation radiation from cosmological simulations

Abstract: We determine the dark matter pair-wise relative velocity distribution in a set of Milky Way-like halos in the Auriga and APOSTLE simulations. Focusing on the smooth halo component, the relative velocity distribution is well-described by a Maxwell-Boltzmann distribution over nearly all radii in the halo. We explore the implications for velocity-dependent dark matter annihilation, focusing on four models which scale as different powers of the relative velocity: Sommerfeld, s-wave, p-wave, and d-wave models. We show that the J -factors scale as the moments of the relative velocity distribution, and that the halo-to-halo scatter is largest for d-wave, and smallest for Sommerfeld models. The J -factor is strongly correlated with the dark matter density in the halo, and is very weakly correlated with the velocity dispersion. This implies that if the dark matter density in the Milky Way can be robustly determined, one can accurately predict the dark matter annihilation signal, without the need to identify the dark matter velocity distribution in the Galaxy.

Velocity-dependent J-factors for annihilation radiation from cosmological simulations

Abstract: We determine the dark matter pair-wise relative velocity distribution in a set of Milky Way-like halos in the Auriga and APOSTLE simulations. Focusing on the smooth halo component, the relative velocity distribution is well-described by a Maxwell-Boltzmann distribution over nearly all radii in the halo. We explore the implications for velocity-dependent dark matter annihilation, focusing on four models which scale as different powers of the relative velocity: Sommerfeld, s-wave, p-wave, and d-wave models. We show that the J-factors scale as the moments of the relative velocity distribution, and that the halo-to-halo scatter is largest for d-wave, and smallest for Sommerfeld models. The J-factor is strongly correlated with the dark matter density in the halo, and is very weakly correlated with the velocity dispersion. This implies that if the dark matter density in the Milky Way can be robustly determined, one can accurately predict the dark matter annihilation signal, without the need to identify the dark matter velocity distribution in the Galaxy.

A unified scenario for the origin of spiral and elliptical galaxy structural scaling laws

Abstract: Elliptical (E) and spiral (S) galaxies follow tight, but different, scaling laws linking their stellar masses, radii, and characteristic velocities. Mass and velocity, for example, scale tightly in spirals with little dependence on galaxy radius (the "Tully-Fisher relation"; TFR). On the other hand, ellipticals appear to trace a 2D surface in size-mass-velocity space (the "Fundamental Plane"; FP). Over the years, a number of studies have attempted to understand these empirical relations, usually in terms of variations of the virial theorem for E galaxies and in terms of the scaling relations of dark matter halos for spirals. We use Lambda Cold Dark Matter (LCDM) cosmological hydrodynamical simulations to show that the observed relations of both ellipticals and spirals arise as the result of (i) a tight galaxy mass-dark halo mass relation, and (ii) the self-similar mass profile of CDM halos. In this interpretation, E and S galaxies of given stellar mass inhabit halos of similar mass, and their different scaling laws result from the varying amounts of dark matter enclosed within their luminous radii. This scenario suggests a new galaxy distance indicator applicable to galaxies of all morphologies, and provides simple and intuitive explanations for long-standing puzzles, such as why the TFR is independent of surface brightness, or what causes the "tilt" in the FP. Our results provide strong support for the predictions of LCDM in the strongly non-linear regime, as well as guidance for further improvements to cosmological simulations of galaxy formation.

The Pristine Dwarf-Galaxy survey -- IV. Probing the outskirts of the dwarf galaxy Bo\"otes I

Abstract: We present a new spectroscopic study of the dwarf galaxy Bootes I (Boo I) with data from the Anglo-Australian Telescope and its AAOmega spectrograph together with the Two Degree Field multi-object system. We observed 36 high-probability Boo I stars selected using Gaia Early Data Release 3 proper motions and photometric metallicities from the Pristine survey. Out of those, 29 are found to be Boo I's stars, resulting in an excellent success rate of 80% at finding new members. Our analysis uses a new pipeline developed to estimate radial velocities and equivalent widths of the calcium triplet lines from Gaussian and Voigt line profile fits. The metallicities of 18 members are derived, including 3 extremely metal-poor stars ([Fe/H] < -3.0), which translates into an exceptional success rate of 25% at finding them with the combination of Pristine and Gaia. Using the large spatial extent of our new members that spans up to 4.1 half-light radii and spectroscopy from the literature, we are able to detect a systemic velocity gradient of 0.15+/-0.10 km s-1 arcmin-1 and a small but resolved metallicity gradient of -0.007+/-0.003 dex arcmin-1. Finally, we show that Boo I is more elongated than previously thought with an ellipticity of {\epsilon} = 0.68+/-0.15. Its velocity and metallicity gradients as well as its elongation suggest that Boo I may have been affected by tides, a result supported by direct dynamical modelling.

Uncovering fossils of the distant Milky Way with UNIONS: NGC 5466 and its stellar stream

Abstract: We examine the spatial clustering of blue horizontal branch (BHB) stars from the u-band of the Canada–France Imaging Survey (CFIS, a component of the Ultraviolet Near-Infrared Optical Northern Survey, or UNIONS). All major groupings of stars are associated with previously known satellites, and among these is NGC 5466, a distant (16 kpc) globular cluster. NGC 5466 reportedly possesses a long stellar stream, although no individual members of the stream have previously been identified. Using both BHBs and more numerous red giant branch stars cross-matched to Gaia Data Release 2, we identify extended tidal tails from NGC 5466 that are both spatially and kinematically coherent. Interestingly, we find that this stream does not follow the same path as the previous detection at large distances from the cluster. We trace the stream across 31° of sky and show that it exhibits a very strong distance gradient in the range 10 < Rhelio < 30 kpc. We compare our observations to simple dynamical models of the stream and find that they are able to broadly reproduce the overall path and kinematics. The fact that NGC 5466 is so distant, traces a wide range of Galactic distances, has an identified progenitor, and appears to have recently had an interaction with the Galaxy’s disc makes it a unique test-case for dynamical modelling of the Milky Way.

The Pristine survey XIV: chemical analysis of two ultra-metal-poor stars

Abstract: Elemental abundances of the most metal-poor stars reflect the conditions in the early Galaxy and the properties of the first stars. We present a spectroscopic follow-up of two ultra metal-poor stars ([Fe/H]<-4.0) identified by the survey {\em Pristine}: Pristine 221.8781+9.7844 and Pristine 237.8588+12.5660 (hereafter Pr 221 and Pr 237, respectively). Combining data with earlier observations, we find a radial velocity of -149.25 $\pm$ 0.27 and -3.18 $\pm$ 0.19 km/s for Pr 221 and Pr 237, respectively, with no evidence of variability between 2018 and 2020. From a one-dimensional (1D) local thermodynamic equilibrium (LTE) analysis, we measure [Fe/H]$_{\rm LTE}$=-4.79 $\pm$ 0.14 for Pr 221 and [Fe/H]$_{\rm LTE}$=-4.22 $\pm$ 0.12 for Pr 237, in good agreement with previous studies. Abundances of Li, Na, Mg, Al, Si, Ca, Ti, Fe, and Sr were derived based on the non-LTE (NLTE) line formation calculations. When NLTE effects are included, we measure slightly higher metallicities: [Fe/H]$_{\rm NLTE}$=-4.40 $\pm$ 0.13 and [Fe/H]$_{\rm NLTE}$=-3.93 $\pm$ 0.12, for Pr 221 and Pr 237, respectively. Analysis of the G-band yields [C/Fe]$_{\rm 1D-LTE} \leq$ +2.3 and [C/Fe]$_{\rm 1D-LTE} \leq$ +2.0 for Pr 221 and Pr 237. Both stars belong to the low-carbon band. Upper limits on nitrogen abundances are also derived. Abundances for other elements exhibit good agreement with those of stars with similar parameters. Finally, to get insight into the properties of their progenitors, we compare NLTE abundances to theoretical yields of zero-metallicity supernovae. This suggests that the supernovae progenitors had masses ranging from 10.6 to 14.4 M$_{\odot}$ and low-energy explosions with 0.3-1.2 $\times$ 10$^{51}$ erg.

Satellite mass functions and the faint end of the galaxy mass-halo mass relation in LCDM

Abstract: The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low mass halos are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass ($M_*$)-halo mass ($M_{200}$) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the Hydrogen-cooling limit. In the second model, $M_*$ scales as a steep power-law of $M_{200}$ with no explicit cutoff, as suggested by recent semianalytic work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultra-faint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the $M_*$-$M_{200}$ relation at the faint end and to elucidate the mechanisms that determine which low-mass halos "light up" or remain dark in the LCDM scenario.

The tidal evolution of the Fornax dwarf spheroidal and its globular clusters

Abstract: The dark matter content of the Fornax dwarf spheroidal galaxy inferred from its kinematics is substantially lower than expected from LCDM cosmological simulations. We use N-body simulations to examine whether this may be the result of Galactic tides. We find that, despite improved proper motions from the Gaia mission, the pericentric distance of Fornax remains poorly constrained, mainly because its largest velocity component is roughly anti-parallel to the solar motion. Translating Fornax's proper motion into a Galactocentric velocity is thus sensitively dependent on Fornax's assumed distance: the observed distance uncertainty, $\pm 8\%$, implies pericentric distances that vary between $r_{\rm peri}\sim 50$ and $r_{\rm peri}\sim 150$ kpc. Our simulations show that for $r_{\rm peri}$ in the lower range of that estimate, an LCDM subhalo with maximum circular velocity $V_{\rm max}=40$ km s$^{-1}$ (or virial mass $M_{200}\approx 10^{10} M_\odot$, as expected from LCDM) would be tidally stripped to $V_{\rm max} \sim 23$ km s$^{-1}$ over $10$ Gyr. This would reduce the dark mass within the Fornax stellar half-mass radius to about half its initial value, bringing it into agreement with observations. Tidal stripping affects mainly Fornax's dark matter halo; its stellar component is affected little, losing less than $5\%$ of its initial mass in the process. We also explore the effect of Galactic tides on the dynamical friction decay times of Fornax's population of globular clusters (GC) and find little evidence for substantial changes, compared with models run in isolation. A population of GCs with initial orbital radii between $1$ and $2$ kpc is consistent with the present-day spatial distribution of Fornax GCs, despite assuming a cuspy halo. Both the dark matter content and the GC population of Fornax seem consistent with current models of dwarf galaxy formation in LCDM.

Satellites Around Milky Way Analogs: Tension in the Number and Fraction of Quiescent Satellites Seen in Observations Versus Simulations

Abstract: We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APOSTLE and Auriga. We compare our detection rates from the NUV and FUV bands to the SAGA-II H$\alpha$ detections and find that they are broadly consistent with over $85\%$ of observed satellites detected in all three tracers. We apply the same spatial selection criteria used around SAGA-II hosts to select satellites around the MW-like hosts in APOSTLE and Auriga. We find very good overall agreement in the derived SFRs for the star-forming satellites as well as the number of star-forming satellites per host in observed and simulated samples. However, the number and fraction of quenched satellites in the SAGA-II sample are significantly lower than those in APOSTLE and Auriga below a stellar mass of $M_*\sim10^{8}\,M_{\odot}$, even when the SAGA-II incompleteness and interloper corrections are included. This discrepancy is robust with respect to the resolution of the simulations and persists when alternative star-formation tracers are employed. We posit that this disagreement is not readily explained by vagaries in the observed or simulated samples considered here, suggesting a genuine discrepancy that may inform the physics of satellite populations around MW analogs.

Structure and kinematics of tidally limited satellite galaxies in LCDM.

Abstract: We use N-body simulations to model the tidal evolution of dark matter-dominated dwarf spheroidal galaxies embedded in cuspy Navarro-Frenk-White subhalos. Tides gradually peel off stars and dark matter from a subhalo, trimming it down according to their initial binding energy. This process strips preferentially particles with long orbital times, and comes to an end when the remaining bound particles have crossing times shorter than a fraction of the orbital time at pericentre. The properties of the final stellar remnant thus depend on the energy distribution of stars in the progenitor subhalo, which in turn depends on the initial density profile and radial segregation of the initial stellar component. The stellar component may actually be completely dispersed if its energy distribution does not extend all the way to the subhalo potential minimum, although a bound dark remnant may remain. These results imply that 'tidally-limited' galaxies, defined as systems whose stellar components have undergone substantial tidal mass loss, neither converge to a unique structure nor follow a single tidal track, as claimed in earlier work. On the other hand, tidally limited dwarfs do have characteristic sizes and velocity dispersions that trace directly the characteristic radius ($r_{max}$) and circular velocity ($V_{max}$) of the subhalo remnant. This result places strong upper limits on the size of satellites whose unusually low velocity dispersions are often ascribed to tidal effects. In particular, the large size of kinematically-cold 'feeble giant' satellites like Crater 2 or Antlia 2 cannot be explained as due to tidal effects alone in the Lambda Cold Dark Matter scenario.

The Pristine survey XIV: Uncovering the very metal-poor tail of the thin disc

Abstract: We evaluate the rotational velocity of stars observed by the Pristine survey towards the Galactic anticenter, spanning a wide range of metallicities from the extremely metal-poor regime ($\mathrm{[Fe/H]} -1.5$ dex. Furthermore, we detect signatures of a ''heated disc'', the so-called $Splash$, at metallicities higher than $\sim-1.5$ dex. Finally, at $\mathrm{[Fe/H]} < -1.5$ dex our anticenter sample is dominated by a kinematical halo with a net prograde motion.