---
abstract: 'We present new FLAMES$+$GIRAFFE spectroscopy of 36 member stars in the isolated Local Group dwarf spheroidal galaxy Tucana. We measure a systemic velocity for the system of $v_{\mathrm{Tuc}}=216.7_{-2.8}^{+2.9}$kms$^{-1}$, and a velocity dispersion of $\sigma_{\mathrm{v,Tuc}}=14.4_{-2.3}^{+2.8}$kms$^{-1}$. We also detect a rotation gradient of $\frac{dv_{r}}{d\chi}=7.6^{+4.2}_{-4.3}$ kms$^{-1}$ kpc$^{-1}$, which reduces the systemic velocity to $v_{\mathrm{Tuc}}=215.2_{-2.7}^{+2.8}$ kms$^{-1}$ and the velocity dispersion to $\sigma_{v,\mathrm{Tuc}}=13.3_{-2.3}^{+2.7}$ kms$^{-1}$. We perform Jeans modelling of the density profile of Tucana using the line–of–sight velocities of the member stars. We find that it favours a high central density consistent with ‘pristine’ subhalos in $\Lambda$CDM, and a massive dark matter halo ($\sim$$10^{10}$M$_\odot$) consistent with expectations from abundance matching. Tucana appears to be significantly more centrally dense than other isolated Local Group dwarfs, making it an ideal laboratory for testing dark matter models.'
author:
- |
    Alexandra L. Gregory$^{1}$[^1], Michelle L. M. Collins$^{1}$, Justin I. Read$^{1}$, Michael J. Irwin$^{2}$, Rodrigo A. Ibata$^{3}$, Nicolas F. Martin$^{3, 4}$, Alan W. McConnachie$^{5}$, Daniel R. Weisz$^{6}$\
    $^{1}$Department of Physics, University of Surrey, Guildford, GU2 7XH, Surrey, UK\
    $^{2}$Institute of Astronomy, Madingley Rise, Cambridge CB3 0HA, UK\
    $^{3}$Université de Strasbourg, CNRS, Observatoire Astronomique de Strasbourg, UMR 7550, F-67000 Strasbourg, France\
    $^{4}$Max-Planck-Institut für Astronomy, Königstuhl 17, D-69117, Heidelberg, Germany\
    $^{5}$National Research Council, Herzberg Institute of Astrophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada\
    $^{6}$Department of Astronomy, University of California, Berkeley, CA 94720, USA
bibliography:
- 'tucana.bib'
date: 'Accepted XXX. Received YYY; in original form ZZZ'
title: 'Kinematics of the Tucana Dwarf Galaxy: An Unusually Dense Dwarf in the Local Group'
---

\[firstpage\]

galaxies:dwarf – galaxies:evolution – galaxies:haloes – galaxies:kinematics and dynamics – Local Group

Introduction {#sec:intro}
============

Isolated dwarf galaxies are insulated from tidal effects and interactions, and they therefore provide a unique window onto the role environment plays in galactic evolution. However, whilst there are many detailed studies of Milky Way and Andromeda satellites (e.g. @simon07 [@tollerud12; @collins13]), comparative observations of isolated dwarfs require a more concerted effort due to the greater distances involved. The ACS LCID project has mapped out the star formation histories of several isolated Local Group dwarfs (e.g. @monelli10 [@gallart15; @aparicio16]), while the Solo Survey [@higgs16] is a wide–field photometric survey targeting isolated dwarfs within 3 Mpc of the Milky Way. Detailed studies of such galaxies are necessary to fully understand galaxy formation and evolution, and may provide answers to several small scale issues with the $\Lambda$ cold dark matter (CDM) paradigm (e.g. @bullock17).

Pure dark matter structure formation simulations in $\Lambda$CDM predict that dwarf galaxies should reside in high density halos [@dubinski91; @navarro96], with a central dark matter density $\rho_{\rm{DM}}(150 \rm{pc}) > 10^8$ M$_{\odot}$ kpc$^{-3}$ [@read18a; @read18b]. This has long been known to be inconsistent with observations of gas rich dwarf irregular galaxies (e.g. @moore94 [@read17a]), which has become known as the ‘cusp-core’ problem. Similarly, the inferred masses of most satellite dwarf galaxies within their half–light radii are also found to be inconsistent with this prediction [@read06b; @boylankolchin11; @boylankolchin12]. Simulated dark matter subhalos appear to be too dense to host the observed dwarf satellites of the Milky Way and Andromeda (M31), a problem referred to as ‘Too Big to Fail’.

  ------------------------------------------------ ------------------------------------------
  **RA$^{a}$**                                     22 41 49.6

  \[3pt\] **Declination$^{a}$**                    -64 25 10

  \[3pt\] **Distance from Milky Way$^{b}$**        887$\pm$49 kpc

  \[3pt\] **Core Radius$^{c}$**                    42$\pm$6“\
                                                   & 176$\pm$26 pc\
                                                   **Half–Light Radius$^{c}$** & 66$\pm$12”

                                                   284$\pm$54 pc

  \[3pt\] **Ellipticity$^{c}$**                    0.48$\pm$0.03

  \[3pt\] **Position Angle of Major Axis$^{c}$**   97$^{\circ}\pm$2$^{\circ}$

  \[3pt\] **Luminosity$^{d}$**                     5.5$\times$10$^{5}$ L$_{\odot}$

  \[3pt\] **Stellar Mass$^{e}$**                   3.2$\times$10$^{6}$ M$_{\odot}$
  ------------------------------------------------ ------------------------------------------

  : Key Properties of Tucana. a) @lavery92; b) @bernard09; c) @saviane96; d) @mateo98; e) @hidalgo13[]{data-label="tab:tuc_properties"}

One elegant solution to both of the above problems is the idea of ‘dark matter heating’ (e.g. @navarro96 [@read05; @pontzen12; @pontzen14; @read16]). In this scenario, repeated gas inflow and outflow cause the central gravitational potential of the dwarf galaxy to fluctuate. The dark matter responds to this by migrating outwards, lowering the inner dark matter density. @read18a have recently found an anti–correlation between the amount of star formation and the central dark matter density in a sample of 16 nearby dwarf galaxies, exactly as expected if dark matter migrates slowly outwards as star formation proceeds. However, their sample of dwarf galaxies is small and may suffer from selection effects. For this reason, it is interesting to measure the inner dark matter density of a larger sample of dwarfs, particularly dwarfs with a purely old stellar population that are expected to retain their ‘pristine’ central dark matter density.

In this paper, we investigate an isolated dwarf spheroidal (dSph) galaxy of the Local Group— Tucana— whose star formation shut down long ago. This makes it a particularly clean test for probing the nature of dark matter, since it is less likely to have had its dark matter ‘heated up’ (e.g. @dicintio14 [@onorbe15; @read16; @bermejo18; @read18a; @read18b]). Tucana was rediscovered and proposed as a Local Group member by @lavery92 after appearances in earlier catalogues. It is located 880 kpc from the Milky Way and 1350 kpc from M31 [@castellani96; @fraternali09], making it one of the most isolated galaxies of the Local Group. Observations by @fraternali09 suggest that Tucana is receding from both the Milky Way and the Local Group, and, if bound, has not yet reached apocentre. It is likely to have been isolated for the majority of its history, although tracing the kinematics to higher redshift may imply a possible interaction with the Milky Way around 10 Gyr ago [@fraternali09]. @sales07 suggest that Tucana’s isolation may be the result of a three body ejection mechanism, potentially involving the Milky Way and the Magellanic Clouds. Along with And XVIII, the other isolated dSph in the Local Group is Cetus, located 775 kpc from the Milky Way [@whiting99; @lewis07]. Both these isolated dwarfs lie in the direction of Sculptor, and it has been postulated that they may form part of a bridge between the Local Group and the Sculptor Group [@whiting99; @fraternali09].

Tucana has experienced no recent star formation. A study by @monelli10 as part of the ACS LCID project found that Tucana experienced a strong burst of star formation $\sim13$ Gyr ago, lasting for  1 Gyr, with star formation having stopped completely by $\sim9$ Gyr ago, with the exception of a small intermediate age, low metallicity population interpreted as owing to contamination by blue stragglers. The authors show that the colour magnitude diagram of Tucana exhibits the typical features of an old stellar population, such as a lack of a blue main sequence. @gallart15 confirm that $90\%$ of Tucana’s stars formed more than 10 Gyr ago. @avilavergara16 find that $75\%$ of Tucana’s history has been spent as a ‘closed box’ with no net inflow or outflow of gas. A similar chemical history was also inferred for Cetus [@avilavergara16]. It has been proposed that the purported interaction between Tucana and the Milky Way some 10 Gyr ago could have stripped enough gas to completely shut down star formation in the galaxy [@teyssier12]. In agreement with the lack of ongoing star formation, @oosterloo96 observe that there is no HI emission within the optical boundary of Tucana. @fraternali09 demonstrate that a nearby (on–sky) HI cloud is more likely to be associated with the Magellanic stream. Through a combination of its unusual location and chemical history, Tucana is a highly unique and interesting Local Group object.

![image](velocity_map.pdf){width="0.9\linewidth"}

A kinematic study of Tucana was previously undertaken by @fraternali09, hereafter F09. This study used the FORS2 instrument at the Very Large Telescope (VLT) to measure the radial velocities of red giant branch stars in Tucana, and hence estimate Tucana’s systemic velocity and velocity dispersion, constraining the mass and metallicity of the dSph. The stellar radial velocities were found by cross correlating a template with the Ca triplet of each star. 20 stars were identified as candidate members of Tucana based on their velocity. Using a maximum likelihood method to fit a Gaussian profile to the histogram of velocities, the authors obtain a systemic velocity of $v=193.0\pm$4.9 kms$^{-1}$ and a velocity dispersion of $\sigma_{{v}}=17.4^{+4.5}_{-3.5}$ kms$^{-1}$. From this they calculate a mass–to–light ratio of $M/L=105^{+95}_{-49}$, implying a mass of $M\mathrm{_{half}}$$\approx6\times10^{7}$ M$_{\odot}$. F09 also detect a rotation signature of magnitude $v_{\mathrm{rot}}=16$ kms$^{-1}$ along the major axis. Accounting for this slightly increases the systemic velocity to $v=194.0\pm4.3$ kms$^{-1}$, and reduces the velocity dispersion to $\sigma_{v}=15.8^{+4.1}_{-3.1}$ kms$^{-1}$. This corresponds to a system mass of $M\mathrm{_{half}}$$\approx5\times10^{7}$ M$_{\odot}$. F09 also measured a mean metallicity of \[Fe/H\]=-1.95$\pm$0.15, with a metallicity dispersion of 0.32$\pm$0.06dex.

The aim of our study is to measure the velocity and velocity dispersion of Tucana to a higher accuracy than before. The velocity dispersion measured by F09 is unusually high and may be consistent with the most massive surviving subhalos in pure dark matter $\Lambda$CDM simulations; however the uncertainties on the result are large. We aim to use higher resolution spectroscopy of a larger stellar sample to constrain the velocity dispersion and density profile of Tucana. This paper is structured as follows: Section \[sec:obs\] describes the observations and data reduction process. We detail the calculation of the velocity dispersion, including determination of the radial velocities and errors, membership probabilities of the Tucana candidates, and calculation of the mass of Tucana in section \[sec:veldisp\]. Details of the process for modelling the density profile are given in section \[sec:gravsphere\]. The implications of these results are discussed in section \[sec:disc\], and we conclude in section \[sec:conclusion\]. Key properties of Tucana are listed in Table \[tab:tuc\_properties\]. These values are used throughout this paper unless otherwise stated.

Observations {#sec:obs}
============

Photometric data for Tucana was obtained using the Magellan/ Megacam instrument on the Clay telescope at Las Campanas observatory on 14th November 2012 as part of the Solo (Solitary Local Dwarfs) observing campaign [@higgs16]. Magellan/ Megacam is a 9$\times$4 array of pixel CCDs with a pixel scale 0.08 arcsec pix$^{-1}$. Two pointings were targeted, with 3 exposures were stacked in both the $g$–band and $i$–band for each field. Each $g$–band exposure was 150s, for a total integration time of 450s. In the $i$-band the exposures were 300s each, for a total integration time of 900s. Seeing ranged from 0.68–0.90 arcsec in the $g$–band and 0.55–0.72 in the $i$–band. The photometry was reduced using the Cambridge Astronomical Survey Unit (CASU) following the process described in @richardson11 and @higgs16. Including only point–like sources, defined as those with classifiers $c_{g}$, $c_{i}=-1$ or -2, this catalogue featured 7203 objects.

To obtain accurate spectroscopy of the faint stellar population of Tucana, the 8.2m Very Large Telescope (VLT) in Paranal, Chile was used. Observations were taken using the FLAMES$+$GIRAFFE spectrograph over 6 nights through June, August and September 2015. GIRAFFE is a fibre–fed spectrograph for the visible range 3700–9000Å. The instrument was used in the Medusa mode, allowing observations of up to 132 objects simultaneously, each with an aperture of 1.2" on the sky. Two fibre configurations were used to maximise the number of targets available. Each exposure was 1200s, and 39 exposures were taken in total: 21 in the first configuration for an integration time of 7 hours; and 18 in the second configuration for a 6 hour integration time. The LR8 grating was used, encompassing a 1190Å wavelength band centred on 8817Å which covers the three peaks of the Ca II triplet at around 8500Å. This setup provides a spectral resolution of $R=6500$, capable of resolving velocity dispersions even in very faint dSphs [@koposov11]. Targets were selected using the Magellan photometry. Overall, 165 individual objects were observed, along with 24 sky regions.

The FLAMES$+$GIRAFFE spectrograph is a higher resolution instrument than FORS2 (which was used in F09), with a resolving power of R=6500 as opposed to R=3200. This allows more accurate measurements of stellar radial velocities, potentially reducing the uncertainties by a factor of 2, and thus better constraining the dispersion. F09 observe 23 stars to a S/N suitable for determining velocities, with 20 of these identified as members. We observe a much larger number of stars out to a wider radius than F09, generating a larger catalogue of Tucana members and reducing the uncertainties in the results.

For the reduction, we used the pipeline for the GIRAFFE instrument provided by ESO. The raw science data was reduced using the `giscience` recipe via the graphical interface `Gasgano`. The pipeline provides corrections for detector effects, including dark and bias corrections, then traces the fibre positions and aligns each fibre with the corresponding spectrum. The spectra are output along with a wavelength calibration and descriptors of the fibre setup and observation. Where possible, calibration frames taken alongside the observations were used (standard calibration files were used for PSF$\_$WIDTH and PSF$\_$CENTROID, which provide the width and centre of the fitted fibre profile).

All spectra were normalised to have unity continuum flux by dividing through by the continuum level. For each observation, the individual sky spectra were median combined to generate a master sky spectrum for each pointing. A median average was chosen as it removes any spurious lines present in individual sky spectra from the master. In order to obtain an accurate sky subtraction, both science and sky were shifted to a continuum level of 0, and the master sky for the relevant pointing normalised by scaling the height of the sky line at 8500Å to match the height of the same sky line in the science spectrum. The sky template was then subtracted from the science spectrum. The resulting spectrum was shifted back to a continuum level of one, and a heliocentric correction was applied to the wavelength scale. Finally, all the spectra of a given object were median combined to give one master spectrum for each object observed. The overall data reduction process results in a set of 165 sky subtracted, heliocentric corrected stellar spectra from Tucana.

Determination of the Velocity Dispersion {#sec:veldisp}
========================================

![Histogram of the radial velocities of all objects, except those with a cross correlation failure. The systemic velocity of Tucana as determined by F09 is shown by the green dashed line, with its uncertainty range shaded, as are the fitted velocities of the two peaks representing Tucana and the Milky Way.[]{data-label="fig:histogram"}](rv_hist_rmfailure.pdf){width="\linewidth"}

Radial Velocities and Errors {#sec:radvel}
----------------------------

The line–of–sight velocities of the stars in Tucana were found by cross correlating the spectrum for each object with a template spectrum. The template used was a Gaussian model of the rest Ca II triplet, with peaks at 8499Å, 8543Å and 8663Å, retaining the relative equivalent widths of the lines. The cross correlation function was calculated using the `pyasl.crosscorrRV` function from PyAstronomy[^2]. Many spectra exhibit obscuring noise around the third Ca II triplet line; to ensure an accurate result, we therefore only use the first two lines in the cross correlation. Based on the expected velocity of Tucana (following the results of F09) and the Milky Way (from a Besançon model— see section \[sec:pvel\]), the allowed velocities were limited to the range -100 kms$^{-1}<v<$300 kms$^{-1}$. The velocity where the cross correlation function is maximised is then taken to be the radial velocity of the object.

Errors on the radial velocities were found by following the iterative Monte Carlo process outlined in @tollerud12. 1000 iterations of each object spectrum were generated by adding noise seeded by the variance per pixel, assuming independent, Poisson distributed noise. Each iteration is cross correlated with the template Ca II triplet, and all 1000 velocities are plotted as a histogram. A Gaussian profile is fitted to the primary peak of the histogram, and the best fit mean and standard deviation are taken as the radial velocity and corresponding error for the given object. There are 25 objects for which the low S/N of the spectrum leads to a cross correlation failure, and 9 for which a Gaussian could not be fitted to the profile. These are removed from the dataset to leave 131 successfully reduced objects. The average S/N of these spectra is 8.4/pixel (with a pixel size of 0.2Å /pixel).

![image](isochrone_4.pdf){width="0.9\linewidth"}

Details of the observed objects are listed in table \[tab:allobj\], along with the details of the corresponding F09 object where applicable. The positions of these objects relative to Tucana are shown in Fig. \[fig:radec\], colour coded by radial velocity. The velocity map indicates a concentration of stars with velocities near the systemic velocity of Tucana in the dense centre of the galaxy. The velocities are plotted as a histogram in Fig. \[fig:histogram\], exhibiting well defined Tucana and Milky Way peaks (see section \[sec:pvel\]).

Determination of Membership {#sec:memprob}
---------------------------

To establish the true velocity profile of Tucana, we must robustly determine which objects are members of Tucana. To do so, we use elements of the probabilistic method outlined in @collins13; namely, position on the colour–magnitude diagram, distance from the centre of Tucana, and the velocity of the object. The probability of membership for a given object is defined by

$$\begin{aligned}
P_{i}\propto P_{\mathrm{CMD}} \times P_{\mathrm{dist}} \times P_{\mathrm{vel}}.
\label{eq:probability}\end{aligned}$$

Below we outline the method for determining $P\mathrm{_{CMD}}$, $P\mathrm{_{dist}}$ and $P\mathrm{_{vel}}$, and hence $P_{i}$.

### Membership Probability Based on CMD Position {#sec:pcmd}

The stellar members of Tucana present in our dataset are expected to lie on the red giant branch (RGB) of the colour–magnitude diagram (CMD). By following the method of @tollerud12 and overlaying appropriate isochrones onto the colour–magnitude diagram, the proximity of each object to the isochrone can be measured using equation \[eq:isomem\],

$$\begin{aligned}
P\mathrm{_{CMD}}=\exp\bigg[-\bigg(\frac{\Delta(g-i)^2}{2\sigma_c}-\frac{\Delta(i)^2}{2\sigma_m}\bigg)\bigg],
\label{eq:isomem}\end{aligned}$$

where $g-i$ is the difference in $g$–band and $i$–band magnitudes, $i$ is the $i$–band magnitude, and $\sigma_c$ and $\sigma_m$ are free parameters which take account of distance and photometry factors. @tollerud12’s value of $\sigma_c=0.1$ was used, but we used $\sigma_m=0.1$ instead of $\sigma_m=0.5$, as this gave a better fit to the Tucana CMD. $P\mathrm{_{CMD}}$ serves as a proxy for the probability of membership.

   **Object**   **Velocity (kms$^{-1}$)**   **RA (deg)**   **Dec. (deg)**   ***g***   ***i***   **S/N (pix$^{-1}$)**   **M?**   **F09**   **F09 Vel. (kms$^{-1}$)**
  ------------ --------------------------- -------------- ---------------- --------- --------- ---------------------- -------- --------- ---------------------------
     21138           116.7$\pm$30.0           340.5952        -64.2066       23.05     21.15            13.8             N         –                  –
     51120            240.7$\pm$8.4           340.3689        -64.4116       22.82     21.71            10.4             Y         –                  –
     34818            -63.0$\pm$3.6           340.2907        -64.3046       22.61     21.11            10.6             N         –                  –
      ...                  ...                  ...             ...           ...       ...             ...             ...       ...                ...
     100021          192.7$\pm$13.3           340.5692        -64.3835       22.50     20.40            10.8             N        22           201.6$\pm$11.5
      ...                  ...                  ...             ...           ...       ...             ...             ...       ...                ...

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.)

A CMD was generated for Tucana as shown in the first panel of Fig. \[fig:CMD\], using the Magellan photometry. To reduce the density of sources around the RGB, only those sources within 4 arcmin of the centre of Tucana were used. A fiducial bounding box is marked to highlight the position of the RGB. The objects in our dataset were matched against the photometric data using on–sky position in order to obtain the photometric properties of each object. Given that the photometric data was used for targeting the spectroscopic observations, most objects had a near direct match in RA/ Dec. These data are overlaid onto Tucana’s CMD in the upper right hand plot of \[fig:CMD\]. As expected, the majority of catalogue sources appear to lie directly on the RGB.

Isochrones were taken from the Dartmouth Stellar Evolution database [@dotter08], and overlaid onto this CMD. Visual inspection suggests that the isochrone which best represents the RGB of Tucana is an isochrone of age 11 Gyr, \[$\alpha$/Fe\]=0, and \[Fe/H\]=-1.5 dex, which is also in agreement with previous measurements [@saviane96; @fraternali09]. Note that the isochrone is purely used as a guide for selecting the most probable members, and so a formal fitting procedure is not required. The isochrone is shifted by the distance modulus $m-M=24.7$, to account for Tucana’s distance of $D=887$ kpc. Equation \[eq:isomem\] is used to determine the proximity of each object to this isochrone. The lower left hand plot of Fig. \[fig:CMD\] shows the CMD of photometric sources with the 11 Gyr, \[$\alpha$/Fe\]=0, \[Fe/H\]=-1.5dex isochrone overlaid. The observed objects are plotted, colour–coded by probability based on isochrone proximity.

### Membership Probability Based on Velocity {#sec:pvel}

We establish the probability of membership based on the velocity of a given object using a maximum likelihood method, first laid out in @martin07 and adapted from @collins13. The radial velocities of all objects are plotted as a histogram in Fig. \[fig:histogram\]. This shows a peak at $v$$\approx$0 kms$^{-1}$ comprising stars in the Milky Way, and a second peak at $v$$\approx$200 kms$^{-1}$ representing the Tucana population. Each peak is well represented by a Gaussian function of the form

$$\begin{gathered}
    P\mathrm{_{peak}}=\frac{1}{\sqrt{2\pi}\sqrt{\sigma_{v,\mathrm{peak}}^2+v_{\mathrm{err}, i}^2}} \\ 
    \times\exp\bigg(-\frac{1}{2}\bigg[\frac{v\mathrm{_{peak}}-v_i}{\sqrt{\sigma_{v,\mathrm{peak}}^2+v_{\mathrm{err}, i}^2}}\bigg]^2\bigg),
\label{eq:probmw}\end{gathered}$$

where $v_i$ is the velocity of the given star, $v_{\mathrm{err}}$ is the uncertainty on its velocity, $v_{\mathrm{peak}}$ is the prior on the velocity of the peak, and $\sigma_{v,\mathrm{peak}}$ is the prior on the velocity dispersion. If $P_{\mathrm{Tuc}}$ is the probability of membership of the Tucana peak, and $P_{\mathrm{MW}}$ is the probability of membership of the Milky Way peak, the overall likelihood function becomes

$$\begin{aligned}
    \mathcal L=\sum_{i=1}^{N} \log(\eta\mathrm{_{MW}} P_{\mathrm{MW},i} + \eta\mathrm{_{Tuc}} P_{\mathrm{Tuc},i}),
\label{eq:likelihood}\end{aligned}$$

where $\eta_{\mathrm{peak}}$ is the fraction of the stars which belong to that peak. Using Bayesian techniques, the probability that a star belongs to Tucana based on its velocity is then given by

$$\begin{aligned}
P\mathrm{_{vel}}=\frac{P_{\mathrm{Tuc},i}}{P_{\mathrm{MW},i} + P_{\mathrm{Tuc},i}}.
\label{eq:pvel}\end{aligned}$$

To measure the velocity and dispersion of Tucana, an MCMC routine is used to fit equation \[eq:probmw\] to the Tucana peak. For this the `emcee` python package from @foremanmackey13 is used. The values of the systemic velocity and dispersion of Tucana published in F09 are used as initial values for $v_{\mathrm{Tuc}}$ and $\sigma_{v,\mathrm{Tuc}}$. The velocity of the Milky Way background in this region of the sky is determined using the Besançon model [@robin03], which established that Milky Way contaminants are expected to have velocities in the range -50 kms$^{-1}<v<$50 kms$^{-1}$. Therefore, $v_{\mathrm{MW}}=0$ kms$^{-1}$ is taken as the initial value for the velocity of this peak, and $\sigma_{v,\mathrm{MW}}=50$ kms$^{-1}$ as the initial value of the velocity dispersion. The analysis uses 100 walkers taking 30000 steps, with a burn in of 6000 steps. Stars with a radial velocity $v>$295 kms$^{-1}$ or $v<$-95 kms$^{-1}$ are rejected, as these velocities lie at the limits of the accepted range for the cross correlation, and so are likely to indicate a cross correlation failure. The initial values of $\eta_{\mathrm{MW}}$ and $\eta_{\mathrm{Tuc}}$ are both taken to be 0.5, based on a visual inspection of the histogram shown in Fig. \[fig:histogram\], and are normalised to sum to 1. The process generates posterior values for the velocity, velocity dispersion and fraction of stars in each peak. This establishes a velocity profile for the two peaks, with the velocity, velocity dispersion and fraction of member stars defined. Including all observed objects, the Tucana peak is found to have a velocity of $v_{\mathrm{Tuc}}=215.7_{-3.1}^{+3.0}$ kms$^{-1}$, and a dispersion of $\sigma_{v,\mathrm{Tuc}}=20.9_{-2.5}^{+2.9}$ kms$^{-1}$, while the Milky Way peak has a velocity of $v_{\mathrm{MW}}=27.5_{-9.5}^{+10.7}$ kms$^{-1}$, and a dispersion of $\sigma_{v,\mathrm{MW}}=61.4_{-7.4}^{+9.5}$ kms$^{-1}$. The quoted uncertainties on these values are the 1 sigma uncertainties returned by the MCMC routine These results are used to produce the Gaussian fits plotted in Fig. \[fig:histogram\]. This velocity is offset from the value measured by F09, and the velocity dispersion is significantly higher. However, this velocity dispersion is artificially increased by the inclusion of non–member stars in the Tucana peak.

We insert the results of the MCMC routine, along with the individual object velocities, into equation \[eq:probmw\] to calculate the probability of each object belonging to the Tucana peak and the Milky Way peak respectively. These probabilities are then inserted into equation \[eq:pvel\] to determine the probability of membership of the given object based on velocity.

### Membership Probability Based on Distance {#sec:pdist}

![Plot showing how the measured systemic velocity and velocity dispersion of Tucana vary as a function of $\mu$, a multiplicative weight applied to the half–light radius in P$_{\rm{dist}}$. The dashed line marks the optimum value used in our analysis. The systemic velocity and velocity dispersion are generally stable, even at large $\mu$.[]{data-label="fig:eta_test"}](eta_test2.pdf){width="\linewidth"}

We also introduce a parameter to account for the distance of the object to the centre of Tucana. This factor is introduced to ensure that objects with a high probability of membership based on their velocity and CMD position, but which are located far outside the centre of the galaxy, are not included as definite members. From @tollerud12, the probability that a star is a Tucana member based on its distance from the centre of the galaxy is given by

$$\begin{aligned}
P\mathrm{_{dist}}=\exp\bigg[-\bigg(\frac{\Delta\alpha^{2}+\Delta\delta^{2}}{2\eta r\mathrm{_{half}}^{2}}\bigg)\bigg],
\label{eq:probdist}\end{aligned}$$

where $r\mathrm{_{half}}$ is the half–light radius of Tucana, $\Delta\alpha^{2}$ and $\Delta\delta^{2}$ are the distances from the object to the centre of Tucana in RA and declination respectively, and $\eta=1.5$ is a free parameter. We use the literature value of $r\mathrm{_{half}}$.

This expression assumes that all member stars will lie within the half–light radius, and hence is designed to give extra weighting to objects within $r\mathrm{_{half}}$. In our dataset, most objects are located outside the half–light radius, as shown in Fig. \[fig:radec\], such that equation \[eq:probdist\] assigns a very low probability of membership to the majority of our dataset. Therefore, to increase the weighting of stars outside the half–light radius, we introduce a multiplicative factor $\mu$ to our definition of $r\mathrm{_{half}}$, such that equation \[eq:probdist\] becomes

$$\begin{aligned}
P\mathrm{_{dist}}=\exp\bigg[-\bigg(\frac{\Delta\alpha^{2}+\Delta\delta^{2}}{2\eta (\mu r\mathrm{_{half}})^{2}}\bigg)\bigg].
\label{eq:probdist2}\end{aligned}$$

In doing so, we are effectively assuming that most observed members lie within $\mu r\mathrm{_{half}}$ of the centre of Tucana, thus ensuring a more gradual decline in membership probability with distance. To find the optimum value of $\mu$, we allow it to vary between 1 and 15, and measure the systemic velocity and velocity dispersion of the resulting sample, as shown in Fig. \[fig:eta\_test\]. We find that the systemic velocity and velocity dispersion are reasonably stable, even at large $\mu$, because the kinematics of Tucana are well separated from the Milky Way contaminants. The optimum value is $\mu=7.0$, which we use in equation \[eq:probdist2\] to obtain $P\mathrm{_{dist}}$. Note that this equation is only used to define our membership sample; it is not influenced by and has no effect on the surface brightness profile used in section \[sec:gravsphere\].

The probability of membership based on velocity, $P_{\mathrm{vel}}$, is combined with the probabilities based on CMD position, $P\mathrm{_{CMD}}$, and distance, $P\mathrm{_{dist}}$, using equation \[eq:probability\], to determine the probability of membership of each object. By removing objects with a membership probability $P_{i}<0.15$, we obtain a population of 37 member stars, details of which are provided in Table \[tab:member\_stars\]. We compare this sample to the dataset generated by F09. Our data was matched to the objects listed in Table 3 of F09 using on–sky position. Table 3 consists of 23 observed objects, 20 of which are identified as members of Tucana. Position matching found that 13 of these members were present in both datasets. Given the significant difference between the systemic velocities obtained by each study, we do not combine the datasets (see section \[sec:f09comparison\] for full discussion).

![Two–dimensional and marginalised PDFs for the systemic velocity and velocity dispersion (both in kms$^{-1}$) of the identified Tucana members, assuming a purely dispersion supported system. The dashed lines represent the mean value and 1$\sigma$ uncertainties.[]{data-label="fig:cornerplot2"}](corner_30000_probability_obj37852.pdf){width="\linewidth"}

   **Object ID**   **Velocity (kms$^{-1}$)**   **RA (deg)**   **Declination (deg)**   **Radius (arcsec)**   ***g***   ***i***   **S/N (pix$^{-1}$)**
  --------------- --------------------------- -------------- ----------------------- --------------------- --------- --------- ----------------------
       51120             240.7$\pm$8.4           340.3689           -64.4116                141.96           22.82     21.71            10.4
       48384             242.4$\pm$8.2           340.6717           -64.3961                342.56           22.92     21.12            13.3
       49196            249.1$\pm$12.3           340.4029           -64.4011                109.11           23.12     21.86            13.4
      100005             202.3$\pm$5.5           340.4345           -64.4083                 55.21           22.43     20.81            8.9
       27260            194.1$\pm$12.7           340.4753           -64.2531                600.41           22.83     21.27            9.4
      100017             211.5$\pm$3.6           340.3983           -64.4114                 97.94           22.52     20.77            10.4
      100006            231.2$\pm$12.3           340.4922           -64.4236                 54.61           22.55     21.18            12.7
      100009             221.8$\pm$4.2           340.4961           -64.4117                 65.44           22.58     20.84            12.5
      100007             196.7$\pm$4.5           340.4203           -64.4161                 60.52           22.77     21.05            7.8
       50350             224.5$\pm$4.2           340.4505           -64.4075                 45.56           22.80     21.50            9.6
       40705            220.6$\pm$10.0           340.6201           -64.3452                367.98           23.55     21.93            12.3
       54894            248.7$\pm$15.3           340.6610           -64.4284                316.65           23.04     21.59            10.9
       50951             234.8$\pm$4.3           340.4445           -64.4107                 38.76           22.94     21.93            11.4
       53700             198.5$\pm$5.1           340.4661           -64.4231                 17.28           23.07     21.73            9.0
       51422             208.6$\pm$2.6           340.3645           -64.4132                147.77           23.17     21.75            8.8
      100013            222.4$\pm$20.8           340.5074           -64.4072                 88.65           22.34     20.66            12.9
      100002             215.2$\pm$5.1           340.4645           -64.4238                 17.61           24.64     24.03            9.9
       52007             219.9$\pm$3.3           340.5419           -64.4159                130.64           23.06     21.86            9.5
       30380            212.4$\pm$10.8           340.2958           -64.2742                581.93           23.01     21.51            13.1
      100016             223.0$\pm$6.7           340.4097           -64.4023                 98.24           22.43     20.79            10.1
       48450             224.7$\pm$9.9           340.6969           -64.3965                380.28           22.85     21.83            12.1
      100001            193.3$\pm$13.6           340.4470           -64.4223                 19.97           22.72     20.78            12.5
       39027            235.7$\pm$18.1           340.5467           -64.3336                339.29           23.05     21.11            12.2
      100003             199.1$\pm$8.9           340.4745           -64.4183                 25.74           24.40     23.64            8.1
       45504            246.5$\pm$21.2           340.2180           -64.3769                404.36           23.02     21.27            12.0
       49630             220.1$\pm$5.7           340.5278           -64.4035                122.74           23.09     21.73            10.4
       24956            205.4$\pm$19.5           340.3787           -64.2358                673.64           23.51     21.87            13.1
      100014            224.6$\pm$10.1           340.4582           -64.3949                 89.42           22.10     20.74            13.0
       28435             228.4$\pm$9.2           340.5639           -64.2612                593.80           23.11     21.73            9.2
       41813             193.1$\pm$2.9           340.5327           -64.3526                267.91           22.64     21.14            14.0
       47569             209.1$\pm$3.4           340.3431           -64.3906                207.54           22.50     21.11            13.8
       54691            226.4$\pm$11.6           340.6335           -64.4276                273.76           23.08     21.93            10.6
       49958             201.0$\pm$0.4           340.4230           -64.4053                 75.50           22.95     21.56            10.7
       34523            247.6$\pm$16.3           340.5140           -64.3024                431.24           22.65     21.01            9.9
       52613             214.1$\pm$7.2           340.6170           -64.4186                246.79           23.10     21.67            9.5
       50487            174.5$\pm$11.6           340.4345           -64.4083                 55.21           22.43     20.81            9.0

Systemic Velocity and Velocity Dispersion {#sysvel}
-----------------------------------------

To generate the final Gaussian fit, we applied an adapted MCMC routine to the 37 identified member stars. We adjust the likelihood function to become

$$\begin{aligned}
\mathcal L=\sum_{i=1}^{N} P_{i} \eta\mathrm{_{Tuc}} P_{\mathrm{Tuc},i},
\label{eq:likelihood2}\end{aligned}$$

thus accounting for the probability of membership in the final routine. This returns a final systemic velocity of Tucana of $v_{\mathrm{Tuc}}=218.3_{-3.1}^{+3.2}$ kms$^{-1}$, and a velocity dispersion of $\sigma_{v,\mathrm{Tuc}}=16.6_{-2.6}^{+3.1}$ kms$^{-1}$. The velocity dispersion is now found to be within 1$\sigma$ of the value calculated by F09, but with significantly smaller error bars. The decrease in the size of the uncertainties is consistent with the increase in sample size relative to F09.

Our probabilistic method returns member stars out to $\sim10r_{\rm{half}}$. This is noted as being a particularly large radius at which to find spectroscopically confirmed members. However, it is not without precedent, as @walker09b find spectroscopic members of Leo V beyond $10r_{\rm{half}}$. In addition, we have shown that the systemic velocity and velocity dispersion are stable to our choice of calibration of the distance probability (see Fig. \[fig:eta\_test\]). Given their strong probabilities of membership, we therefore choose to retain the more distant stars in our membership sample.

One identified member, object ID 37852, has a particularly high velocity of $v_i=277.7\pm12.8$ kms$^{-1}$, some 60 kms$^{-1}$ higher than the measured systemic velocity of Tucana, and hence could be considered a potential contaminant in the sample. It has no counterpart in the F09 dataset, so we cannot make a direct comparison of the velocity measurements, although it does not appear to have significant error bars compared to the rest of our sample. The object is located 5.68 arcmin from the centre of Tucana. Given that it lies $>3\sigma_{v}$ from the mean velocity, we choose to remove this object from our membership sample, in order to avoid including any spurious members. Removing object 37852 reduces the velocity dispersion to $\sigma_{v,\mathrm{Tuc}}=14.4^{+2.8}_{-2.3}$ kms$^{-1}$, and the systemic velocity is reduced to $v_{\mathrm{Tuc}}=216.7^{+2.9}_{-2.8}$ kms$^{-1}$. The 1– and 2– dimensional probability distributions for the velocity and velocity dispersion are shown in Fig. \[fig:cornerplot2\]. These values are smaller than the results from the full sample, but fall within the uncertainties. We therefore use this sample of 36 member stars for the remainder of the analysis. We note that there is still a strong probability that object 37852 is a member of Tucana, but choose to be conservative in our membership definition so as not to overstate our results. Details of object 37852 are included in table \[tab:allobj\].

The histogram of radial velocities of the identified Tucana members is shown in Fig. \[fig:velrad\]. The systemic velocity of the peak has increased slightly, and is more than 3$\sigma$ outside of F09’s value. It corresponds to a velocity relative to the Local Group of $v_{\mathrm{LG}}=+95.4$ kms$^{-1}$, confirming that the galaxy is receding from the Local Group.

### Comparison with Previous Study {#sec:f09comparison}

As previously noted, there is a significant offset between the systemic velocities of our FLAMES dataset and that measured by F09 of $\Delta v_{\mathrm{Tuc}}=22.7$ kms$^{-1}$. To investigate this, we tested the wavelength calibration of the spectra, which is provided by the initial data reduction pipeline, by cross correlating the raw spectra (before sky subtraction) with a reference sky template. This sky spectrum, which was produced during observations of the Sculptor dwarf galaxy by @tolstoy04, provides a reliable reference with which we can calibrate our data. Spectra with a velocity shift $>7$ kms$^{-1}$ were rejected as having a significant wavelength offset which could later affect our analysis of the radial velocity of the object. This process returned negligible shifts between the science and template spectra, confirming the validity of the wavelength calibration. We also rigorously tested the heliocentric correction applied to the spectra to ensure this was accurate. In a further test of the data, we randomly selected a sample of 20 stars (the size of the F09 sample) from our 36 members, and measured the systemic velocity of this sample. After 100 iterations, the difference between the sample velocity and the $v=194.0$ kms$^{-1}$ measured by F09 was $\Delta v>10$ kms$^{-1}$ in all cases, $\Delta v>20$ kms$^{-1}$ in 75 cases, and $\Delta v>22.7$ kms$^{-1}$ in 49 cases. The average offset was $\bar{\Delta v}=22.9$ kms$^{-1}$, consistent with the offset measured for the full dataset.

![Top panel: Histogram of all measured velocities with members highlighted in blue, and the systemic velocity of Tucana as measured by this study and F09 marked as dashed lines. Bottom panel: Velocities of all observed objects plotted as a function of radius from the centre of Tucana, with the identified members marked as blue stars. The dashed lines mark 1$\times$, 2$\times$, 5$\times$ and 10$\times$ the half–light radius of Tucana.[]{data-label="fig:velrad"}](velocity_radius_histogram.pdf){width="\linewidth"}

Comparing the velocities of the stars that overlap the two samples (see table \[tab:allobj\]), we find some variation in the measurements but no systematic offset. Excluding object 100020 (21 in F09), which has an offset of 210 kms$^{-1}$ and is clearly indicative of a mismeasurement in one of the samples; in the overlapping objects we measure an average velocity offset of 11.3kms$^{-1}$. This value is similar to the average velocity error of our identified members ($\bar{{v}_{\mathrm{err}}}=9.3$ kms$^{-1}$) but is somewhat smaller than the offset between the measured systemic velocities. We also note that the 3 non members listed in Table 3 of F09 are matched to non members in the FLAMES data. Many of the objects display an offset $>20$ kms$^{-1}$; however all but one (with the exception of object 100020) have velocities consistent within 3–$\sigma$.

Therefore, whilst the offset appears significant, we believe our result for the velocity can be considered reliable, especially given the larger sample size relative to F09. Furthermore, the offset does not affect our measurement of the velocity dispersion nor any conclusions drawn from this.

### Mass of Tucana {#sec:tucmass}

A number of estimators have been defined to determine the virial mass of a dwarf galaxy from its velocity dispersion. These estimators take the general form

$$\begin{aligned}
    M\mathrm{_{est}}(<\lambda\rm{R})=\frac{\mu r\mathrm{_{half}} \sigma_{v,\mathrm{half}}^{2}}{G},
\label{eq:mass}\end{aligned}$$

where $r_{\mathrm{half}}$ is the half–light radius of the galaxy, and $\sigma_{v,\mathrm{half}}$ is the velocity dispersion at that radius. Assuming a flat velocity dispersion profile, @walker09 define $\lambda$=1 and $\mu$=3.5, such that the mass estimate becomes

$$\begin{aligned}
    M(<r_{\rm{half}})=580 r\mathrm{_{half}} \sigma_{v,\mathrm{half}}^{2}
\label{eq:walkermass}\end{aligned}$$

From this, we determine the mass of Tucana to be $M_{\mathrm{half}}=3.4^{+1.5}_{-1.3}\times10^{7}$ M$_{\odot}$, corresponding to a mass–to–light ratio within the half–light radius of $M/L\mathrm{_{half}}$$\approx62^{+27}_{-23}$ M$_{\odot}$/L$_{\odot}$, assuming a luminosity for Tucana of $5.5\times10^{5}$ L$_{\odot}$. This result is slightly lower than the estimated mass quoted in F09, who found a pressure supported total mass of $M\approx5.0\times10^{7}$ M$_{\odot}$. However, this was determined by assuming that mass follows light in the system, and using an estimator to determine the mass to light ratio, from which the mass is calculated. Given that Tucana is a centrally dense galaxy which is likely to be highly dark matter dominated, this method is less accurate than equation \[eq:walkermass\], which derives the mass directly from the measured velocity dispersion. Using the @walker09 method with the velocity dispersion calculated by F09 returns a value of $M_{\mathrm{half}}=4.1^{+2.3}_{-1.8}\times10^{7}$ M$_{\odot}$, which is consistent with our result.

Recently, @errani18 redefined the mass estimator with $\lambda$=1.8 and $\mu$=3.5. This does not require a flat velocity dispersion profile, and so is insensitive to any fluctuations in the profile (see e.g. Fig. \[fig:gravsphere\_result\]). Using this new estimator, we determine a mass for Tucana of $M_{\mathrm{half}}=8.6^{+3.7}_{-3.2}\times10^{7}$ M$_{\odot}$, corresponding to a mass–to–light ratio within the half–light radius of $M/L\mathrm{_{half}}\approx156^{+68}_{-58}$ M$_{\odot}$/L$_{\odot}$. This is consistent with the result of F09, though it is not directly comparable to the Walker estimate as the enclosed radius is larger.

Rotation in Tucana
------------------

![Radial velocities of the 36 identified Tucana members as a function of their projected position along the major axis. The red dashed line highlights the measured velocity gradient of $\frac{dv_{r}}{d\chi}=7.6^{+4.2}_{-4.3}$ kms$^{-1}$ kpc$^{-1}$. The grey dashed line shows the measured systemic velocity of the system.[]{data-label="fig:velmajorax"}](velocity_gradient_majoraxis_new.pdf){width="\linewidth"}

F09 find a rotation signature for Tucana of $v_{\mathrm{rot}}=16$ kms$^{-1}$, by fitting a rotation curve to the radial velocities of their member stars. To search for this in our data, we plot the radial velocities of our member stars as a function of projected distance along the major axis in Fig. \[fig:velmajorax\]. This highlights a possible gradient across the data. We quantify this gradient by utilising another `emcee` routine. @martin10 redefine the Gaussian representing the Tucana peak as

$$\begin{aligned}
P_{\mathrm{Tuc}}=\frac{1}{\sqrt{2\pi}\sqrt{\sigma_{v, \mathrm{Tuc}}^2+v_{\mathrm{err}, i}^2}} \times  \exp\bigg(-\frac{1}{2}\bigg[\frac{\Delta v_{r,i}^{2}}{\sqrt{\sigma_{v, \mathrm{Tuc}}^2+v_{\mathrm{err}, i}^2}}\bigg]^2\bigg),
\label{eq:gradient}\end{aligned}$$

where $\Delta v_{r,i}$ is the difference between the velocity of a star and a velocity gradient $\frac{dv_{r}}{d\chi}$ acting along the angular distance of a star along an axis $y_i$ with position angle $\theta$, as shown in equation \[eq:deltav\].

$$\begin{aligned}
\Delta v_{r,i}=v_{r,i}-\frac{dv_{r}}{d\chi}y_{i}+\bar{v_{r}}
\label{eq:deltav}\end{aligned}$$

$y_i$ can be determined from the RA and Dec. of the star, ($\alpha_{i}$, $\delta_{i}$), and of the centre of Tucana, ($\alpha_{0}$, $\delta_{0}$), using

$$\begin{aligned}
y_i&=&X_{i}\sin{\theta}+Y_{i}\cos{\theta}\\
X_i&=&(\alpha_{i}-\alpha_{0})\cos(\delta_{0})\\
Y_i&=&\delta_{i}-\delta_{0}.
\label{eq:gradeqns}\end{aligned}$$

By replacing equation \[eq:probmw\] with equation \[eq:gradient\] in the MCMC routine, we can generate a fit to the data which accounts for the velocity gradient $\frac{dv_{r}}{d\chi}$ produced by rotation. We introduce flat priors for the new parameter $\frac{dv_{r}}{d\chi}$ such that $-40<\frac{dv_{r}}{d\chi}<40$. If $\theta$ is fixed to match the position angle of the major axis ($\theta=97^{\circ}$), we measure a rotation gradient of $\frac{dv_{r}}{d\chi}=7.6^{+4.2}_{-4.3}$ kms$^{-1}$ kpc$^{-1}$, with the systemic velocity measured as $v_{\mathrm{Tuc}}=215.2_{-2.7}^{+2.8}$ kms$^{-1}$ and the velocity dispersion as $\sigma_{v,\mathrm{Tuc}}=13.3_{-2.3}^{+2.7}$ kms$^{-1}$. This rotation gradient, marked as the red dashed line in Fig. \[fig:velmajorax\], equates to a rotation velocity of $2.2\pm1.2$kms$^{-1}$ at the half light radius. If $\theta$ is allowed to evolve freely (with a flat prior of $0<\theta<\pi$), we determine the best values to be: $v_{\mathrm{Tuc}}=214.9_{-3.2}^{+3.2}$ kms$^{-1}$; $\sigma_{v,\mathrm{Tuc}}=13.5_{-2.3}^{+2.8}$ kms$^{-1}$; $\frac{dv_{r}}{d\chi}=6.1^{+4.6}_{-4.8}$ kms$^{-1}$kpc$^{-1}$; $\theta=86.5^{\circ}$$^{+37.8}_{-35.5}$; amounting to a rotation velocity of $1.7\pm1.3$kms$^{-1}$ at the half light radius.

It therefore appears that there is a small velocity gradient in Tucana due to the presence of rotation, consistent with alignment with the major axis. The mass estimator from @walker09 (equation \[eq:mass\]) assumes a velocity dispersion dominated system. If we recalculate the dynamical mass using the smaller velocity dispersion accounting for rotation, we obtain a mass of $M_{\mathrm{half}}=2.9^{+1.3}_{-0.9}\times10^{7}$ M$_{\odot}$, corresponding to a mass–to–light ratio within the half–light radius of $M/L\mathrm{_{half}}$$\approx53^{+24}_{-17}$ M$_{\odot}$/L$_{\odot}$. These results are within the uncertainty ranges of the non–rotating result, and given that the rotation gradient is small relative to the velocity dispersion, Tucana is still classed as a dispersion dominated galaxy.

Modelling the Dark Matter Density Profile {#sec:gravsphere}
=========================================

The high velocity dispersion of Tucana suggests that it has a high central dark matter density. In the absence of star formation, a steep central cusp is predicted to be present in the density profiles of all galaxies [@navarro96], but thus far the majority of observations of dSphs show a slight preference towards flattened cores (see the discussion in section \[sec:intro\]). The dearth of recent star formation in Tucana make it a strong candidate for hosting a ‘pristine’ cusp [@brook15; @read18a; @bermejo18]. In this section, we perform Jeans modelling of the data for Tucana to estimate its central dark matter density and quantitatively test this idea.

![image](output_Sigstar.pdf){width="\linewidth"}

![image](output_sigLOS.pdf){width="\linewidth"}

The <span style="font-variant:small-caps;">GravSphere</span> Mass Modelling Code
--------------------------------------------------------------------------------

@read17b introduced a new non–parametric Jeans code <span style="font-variant:small-caps;">GravSphere</span>, which returns the density profile $\rho$(r) and velocity anisotropy $\beta$(r) of a system using only line–of–sight velocities (see also @read18a). <span style="font-variant:small-caps;">GravSphere</span> solves the Jeans equation [@jeans1922] for a set of ‘tracers’ (i.e. stellar members) of a spherical mass distribution defined by the radial density profile $\rho$(r) and velocity anisotropy $\beta$(r). The projected spherical Jeans equation is given by @binney82 as

$$\begin{aligned}
\sigma^{2}_{\mathrm{LOS}}(R)=\frac{2}{\Sigma(R)}\int_{R}^{\infty}\bigg(1-\beta\frac{R^{2}}{r^{2}}\bigg)v\sigma^{2}_{r}\frac{r \mathrm{d}r}{\sqrt{r^{2}-R^{2}}},
\label{eq:jeans}\end{aligned}$$

where $\Sigma(R)$ is the surface mass profile at projected radius R, $v(R)$ is the spherically averaged tracer density, and $\beta$(r) is the velocity anisotropy,

$$\begin{aligned}
\beta=1-\frac{\sigma_t^2}{\sigma_r^2}.
\label{eq:beta}\end{aligned}$$

$\sigma_t$ and $\sigma_r$ are the tangential and radial velocity dispersion profiles. $\sigma_r$ is given by [@vandermarel94; @mamon05]

$$\begin{aligned}
\sigma_r^2(r)&=&\frac{1}{v(r) g(r)}\int_r^\infty\frac{G M(\tilde{r}) v(\tilde{r})}{\tilde{r}^2} g(\tilde{r}) \mathrm{d}\tilde{r}\\
g(r)&=&\exp\bigg(2\int \frac{\beta(r)}{r}\mathrm{d}r \bigg).
\label{eq:sigr}\end{aligned}$$

$M(r)$ is the cumulative mass of the galaxy. <span style="font-variant:small-caps;">GravSphere</span> uses a non–parametric model for $M(r)$, consisting of a contribution from all visible matter, and a contribution from dark matter modelled by a series of power laws centred on a set of radial bins. The tracer light profile uses a series sum of Plummer spheres [@plummer11] and so is also non–parametric. The code fits this model to the surface density profile $\Sigma_{*}(R)$ and line–of–sight velocity dispersion profile $\sigma_{\mathrm{LOS}}(R)$ using `emcee` [@foremanmackey13]. To avoid infinities in $\beta$, a symmetrised version of $\beta$ is used in the model. $\tilde{\beta}$ is defined in equation \[eq:betastar\], and describes the distribution of velocities in the system, with $\tilde{\beta}=0$ describing an isotropic velocity distribution, $\tilde{\beta}=1$ a fully radial distribution and $\tilde{\beta}=-1$ a fully tangential distribution [@read06b].

$$\begin{aligned}
\tilde{\beta}=\frac{\sigma_r^2-\sigma_t^2}{\sigma_r^2+\sigma_t^2}=\frac{\beta}{2-\beta}.
\label{eq:betastar}\end{aligned}$$

$\tilde{\beta}$ is given a flat prior of $-1<\tilde{\beta}<1$, to give equal weight to fully radial and fully tangential distributions. <span style="font-variant:small-caps;">GravSphere</span> also uses virial shape parameters to obtain constraints on $\beta$ using only line–of–sight velocities, thus breaking the well–known $\rho$–$\beta$ degeneracy [@binney82; @merrifield90; @read17b]. The code can be used to determine $\rho$(r) and $\beta$(r) for any near–spherical stellar system, such as star clusters, spheroidal or elliptical galaxies and galaxy clusters. <span style="font-variant:small-caps;">GravSphere</span> has been extensively tested on mock data (@read17b, and see below). For a more complete description of <span style="font-variant:small-caps;">GravSphere</span>, we refer the reader to @read17b.

In addition to the ‘free–form’ dark matter model described above, GravSphere can also fit the cosmologically–motivated <span style="font-variant:small-caps;">coreNFW</span> profile. This profile was originally designed to fit simulations of halos that have undergone ‘dark matter heating’ [@read16], but also provides a good fit to dark matter halos in a self–interacting dark matter cosmology [@read18a]. It has the advantage that it fits a dark matter core size parameter that can be connected to a self–interaction cross section for dark matter [@read18a], and a halo virial mass, $M_{200}$, that can be compared with cosmological expectations from abundance matching (e.g. @read18c). We present the results of fitting this <span style="font-variant:small-caps;">coreNFW</span> model to Tucana in appendix \[appendixa\].

<span style="font-variant:small-caps;">GravSphere</span> has been extensively tested on mock data, including mocks that break spherical symmetry, that include foreground contamination and binary stars, and that are being tidally stripped by a larger host galaxy [@read17b; @read18a]. However, all tests to date have focussed on mocks with $>$500 member velocities, which is an order of magnitude more than we have available for Tucana. For this reason, we present additional mock data tests in appendix \[appendixb\]. These are set up to mimic Tucana, modelling the selection function, contamination and sampling similarly to our real Tucana data. We show that we are able to correctly infer the dark matter density at 150pc from the centre of Tucana with kinematics for just 36 member stars, albeit with substantially larger uncertainties than we obtain for data with 500 member stars. Similarly to the findings in @read18a, we find that <span style="font-variant:small-caps;">GravSphere</span> is not able to distinguish cusps from cores with this many member stars. However, as explained in @read18a [@read18b], an inference of the central dark matter density, $\rho_{\rm DM}(150\,{\rm pc})$, is sufficient to constrain interesting dark matter models. For this reason, we focus in this paper on obtaining an estimate of $\rho_{\rm DM}(150\,{\rm pc})$ for Tucana.

Dark Matter Density Profile
---------------------------

![The density profile of Tucana generated by <span style="font-variant:small-caps;">GravSphere</span>. Grey contours show 1– and 2–$\sigma$ uncertainty ranges, and the vertical line marks the half–light radius of Tucana as modelled by <span style="font-variant:small-caps;">GravSphere</span>. Dashed red lines mark NFW profiles of different pre–infall halo masses. The data weakly favour the $10^{10}$ M$_{\odot}$ model, while the central density $\rho_{\rm DM}(150\,{\rm pc}) > 10^8$M$_\odot$kpc$^{-3}$ at better than 2-sigma confidence, consistent with a dark matter cusp.[]{data-label="fig:gravsphere_result"}](tucana_paper.pdf){width="\linewidth"}

In this section, we use <span style="font-variant:small-caps;">GravSphere</span> to estimate the dark matter density profile of Tucana. For this, we require the surface brightness profile, $\Sigma_*(R)$, and velocity dispersion profile, $\sigma_{\rm{LOS}}(R)$, of our member stars. The former was generated from RGB stars in the Magellan/ MegaCam photometry, adjusted to correct for the fact that our observed fields are misaligned with the central coordinates of Tucana. The velocity dispersion profile was modelled using the radial velocities we have obtained with FLAMES+GIRAFFE. To generate the velocity dispersion profile, the probabilities of membership of our stellar sample are summed, resulting in a membership probability weighted number of stars $N_{\rm{eff}}=\sum_{i=1}^{N_{\rm{mem}}}P_{\rm{mem,i}}$. For our 36–strong sample, this equates to $\sim$20 effective members. The data is binned by radius from the centre of Tucana, with each bin containing the same membership weighted number of stars. The $\Sigma_{*}(R)$ and $\sigma_{\mathrm{LOS}}(R)$ profiles are shown in Fig. \[fig:gravsphere\_models\]. Note that the vertical blue line on these plots marks the half–light radius, $r_{\mathrm{half}}=340$ pc, as modelled from the surface brightness profile. This is just outside 1–$\sigma$ larger than the literature value (see table \[tab:tuc\_properties\]).

We model Tucana with <span style="font-variant:small-caps;">GravSphere</span> under the assumption that it is a spherical, non rotating system. Although we detect a small rotation gradient in Tucana, its effect on the velocity dispersion measurement is negligible with respect to the size of the uncertainties, and so it has no effect on the <span style="font-variant:small-caps;">GravSphere</span> model. We use a stellar mass for Tucana of M$_*$=3.2$\times10^6 M_{\odot}$ [@hidalgo13], assuming an error on M$_*$ of 25$\%$. The final density profile is shown in Fig. \[fig:gravsphere\_result\]. The red dashed lines highlight NFW profiles at different masses. These were produced using the concentration–mass relation from @dutton14, but multiplied by 1.4 to account for the fact that subhalos in the Aquarius simulations are found to be systematically more concentrated than field halos [@springel08].

The profile favours a high central density of $\rho_{\mathrm{DM}}$(150 pc)$=5.5_{-2.5}^{+3.2} \times 10^8$ M$_{\odot}$ kpc$^{-3}$, suggesting that Tucana is as dense, if not more so, than Draco [@read18a]. The profile appears consistent with the presence of a pristine cusp within 1$\sigma$. Using the <span style="font-variant:small-caps;">core</span>NFW fit presented in appendix \[appendixa\], we also infer a pre–infall halo mass of $M_{200}=1.37^{+0.49}_{-0.44} \times 10^{10}$ M$_{\odot}$ for Tucana, consistent with the findings of @brook15.

Discussion {#sec:disc}
==========

Comparison to Previous Work
---------------------------

![image](tbtfplot_new.pdf){width="90.00000%"}

Our velocity dispersion of $\sigma_{\mathrm{v,Tuc}}=14.4_{-2.3}^{+2.8}$kms$^{-1}$ is remarkably similar to that measured by F09, with a difference of just $\sim1$ kms$^{-1}$. It should be noted that the $\sigma_{v,\mathrm{Tuc}}=15.8$ kms$^{-1}$ quoted by F09 accounts for a potential rotation in signature Tucana of 16 kms$^{-1}$. Our non–rotational dispersion result is within 1$\sigma$ of their non–rotational value of $\sigma_{v,\mathrm{Tuc}}=17.4^{+4.5}_{-3.5}$ kms$^{-1}$, so is also consistent with the like–for–like result. We detect a slightly smaller rotation signature of $\frac{dv_{r}}{d\chi}=7.6^{+4.2}_{-4.3}$ kms$^{-1}$.

Fig. \[fig:velrad\] plots the velocities of the observed stars as a function of radius from the centre of Tucana, with those identified as members highlighted in blue. All members are shown to lie within $\sim10 r\mathrm{_{half}}$ of Tucana. There is a selection bias in our observation method; by using a fibre spectrograph, we preferentially observe more distant members, as they are less closely packed, compared to those in the dense centre. Indeed, there is a large population of stars with velocities close to the systemic velocity of Tucana outside $10 r\mathrm{_{half}}$. We have defined our probability functions such that these objects are not selected as members due to their large distances.

There is a substantial velocity offset between the systemic velocity quoted by F09 and that measured from our data of around $\Delta v_{\mathrm{Tuc}}$ $\approx23$ kms$^{-1}$. This could be due to a number of factors, such as a wavelength miscalibration or a bias in the cross correlation procedure used. We also note that as a low resolution slit–based spectrograph, FORS2 (as used in F09) is not optimised for velocity measurements and may suffer from variations in the velocity zero point, as discussed in @kacharov17’s study of the Phoenix dwarf galaxy. Without repeat measurements from FLAMES+GIRAFFE, we cannot conclusively identify the origin of the velocity offset. However, we have thoroughly tested our calibration, as described in section \[sec:f09comparison\], and have confidence in our result. It also should be noted that the offset does not influence our velocity dispersion measurement or any ensuing conclusions. Despite the offset, our result supports the conclusion of F09 that the galaxy is receding from the Milky Way, and could continue to do so, unbound from the Local Group.

Tucana: A Massive Failure
-------------------------

Our velocity dispersion result suggests a high central density consistent with Tucana residing within a supposed ‘massive failure’ halo as predicted to exist by @boylankolchin11. This would make Tucana the first known exception to the ‘too–big–to–fail’ problem, whereby simulated subhalos are too centrally dense to host the observed dwarf galaxies [@boylankolchin11; @boylankolchin12]. In @wang12, too–big–to–fail is restated in terms of the galaxies’ maximum circular velocities; all known satellites are observed with $v_{\mathrm{max}}<30$ kms$^{-1}$ (with the exception of the Magellanic Clouds— @jiang15), yet the simulated halos should host galaxies with 40 kms$^{-1}<v_{\mathrm{max}}<60$ kms$^{-1}$. Our measured velocity dispersion corresponds to a circular velocity for Tucana of $v_{\mathrm{circ}}=22.7_{-4.7}^{+5.4}$ kms$^{-1}$. We plot this velocity alongside those of the nine bright Local Group dSphs (those plotted in @boylankolchin12, Fig. 1), plus the isolated dSph Cetus, as a function of half–light radius in Fig. \[fig:tbtf\_plot\]. The circular velocity of Tucana is significantly higher than that of other Local Group dwarf galaxies, indicating that Tucana is much more centrally dense than the typical dSph. Plotting rotation curves for different maximum circular velocities suggests that Tucana resides in a halo with maximum velocity $v_{\mathrm{max}}>40$ kms$^{-1}$. In other words, Tucana appears to reside in a ‘massive failure’ halo. If this is the case, Tucana would be the only dwarf spheroidal galaxy in the Local Group known to reside in such a halo.

![Inner DM density as a function of pre–infall halo mass for a selection of Local Group dwarf galaxies, as modelled in @read18b. Black points shows galaxies whose star formation was truncated more than 6 Gyrs ago; blue points show galaxies truncated less than 3 Gyrs ago; and the purple points were truncated between 3 and 6 Gyr ago. The shaded regions highlight the expected DM density of cusped and cored profiles. Tucana is shown as the red point. In addition to highlighting the correlation between star formation history and the shape of the density profile, this plot shows the high central DM density of Tucana relative to other Local Group dwarf galaxies, a result which appears to be consistent with a central cusp.[]{data-label="fig:dm_heating"}](dmheating_M200.pdf){width="\linewidth"}

The too–big–to–fail problem has been widely viewed as a universal problem with the $\Lambda$CDM model. Originally discovered to be an issue in the Milky Way satellite population, both @collins14 and @tollerud14 find that the problem is also present in M31 satellites. @read17a show that both too–big–to–fail and ‘missing satellites’ are confined to group environments, and so must be the result of galaxy formation physics as opposed to some exotic formulation of the cosmology. Our result shows that there are conditions under which galaxies can retain their central mass and reside in the $\Lambda$CDM predicted halos. Tucana’s isolation from tidal effects and quenched star formation history mean it has been unaffected by the baryonic feedback effects usually invoked to resolve too–big–to–fail. This would make Tucana unique even among isolated dwarfs: @kirby14 find that star formation in other isolated Local Group galaxies has been energetic enough to lower the central density, such that these galaxies are fully consistent with Milky Way/ M31 satellites. However, the shut–down in star formation in Tucana around 10 Gyr ago [@monelli10] likely allowed it to retain its central mass.

A Pristine Cusp in Tucana
-------------------------

Fig. \[fig:dm\_heating\] shows the central DM density for a range of Local Group dwarfs (the gas rich ones were not modelled with GravSphere) as a function of pre–infall halo mass [@read18b]. The figure highlights a strong correlation between star formation history and inner DM density, with galaxies whose star formation shut down long ago possessing central densities consistent with a central cusp, and those with extended star formation histories consistent with a cored profile. As described in @read18b, this supports the hypothesis that dark matter is heated up by extended, bursty star formation, reducing the central density and transforming cusps to cores.

Using the abundance matching technique from @read18c, with a stellar mass of M$_*=3.2\times10^6$ M$_{\odot}$ [@hidalgo13] and the star formation history for Tucana from @monelli10, we estimate a halo mass for Tucana of M$_{200, \rm{abund}}=9.4 \pm 3.6\times10^9$ M$_{\odot}$, in excellent agreement with our <span style="font-variant:small-caps;">GravSphere</span> models (see appendix \[appendixa\]). To be consistent with the results presented in @read18b, we use M$_{200, \rm{abund}}$ and $\rho_{DM}$(150 pc) as inferred using <span style="font-variant:small-caps;">GravSphere</span> to plot Tucana on Fig. \[fig:dm\_heating\] (red data point). The high central density returned by our <span style="font-variant:small-caps;">GravSphere</span> modelling places Tucana above the other dwarf galaxies in this plot, and is consistent with the expectations of a cusped profile within 1$\sigma$. This result is again consistent with the limited period of star formation experienced by Tucana.

The errors on the inner DM density are very large. However, Tucana is so dense (i.e. its velocity dispersion is so high) that we can be confident that it is more dense than WLM, Fornax and the other isolated dIrrs at 150pc at better than 95$\%$ confidence. Tucana is consistent with expectations for the inner density of ‘pristine’ DM halos in $\Lambda$CDM that have undergone no DM heating, as expected from Tucana’s old stellar population. Based on abundance matching and stellar kinematics, Tucana has a total halo mass consistent with WLM ($M_{200}=0.83\pm 0.2 \times 10^{10}$ M$_{\odot}$, @read17a). WLM’s inner rotation curve favours a central dark matter core over a cusp and is substantially less dense than what we find for Tucana in this work, as shown in Fig. \[fig:dm\_heating\]. This is because unlike WLM, Tucana’s star formation ceased after just $\sim$1-2 Gyr –most likely due to ram pressure stripping by the Milky Way (e.g. @teyssier12 [@gatto13]). Our results are consistent with the idea that Tucana is much denser than WLM because it had insufficient star formation to undergo significant DM heating.

Conclusions {#sec:conclusion}
===========

Using the GIRAFFE spectrograph, we have taken high resolution spectra of the Tucana dwarf galaxy, identifying 36 member stars, and used a cross correlation method to measure their radial velocities. We make several key findings:

-   We find the systemic velocity of Tucana to be $v_{\mathrm{Tuc}}=216.7_{-2.8}^{+2.9}$kms$^{-1}$, corresponding to $v_{\mathrm{GSR}}=121.7_{-2.8}^{+2.9}$kms$^{-1}$. This velocity is receding from the Local Group, and is consistent with the conclusion of F09 that Tucana has long been an isolated dwarf galaxy, which may have interacted with the Milky Way some 10 Gyr ago.

-   We measure a rotation gradient across Tucana of $\frac{dv_{r}}{d\chi}=7.6^{+4.2}_{-4.3}$ kms$^{-1}$ kpc$^{-1}$, which equates to a rotation velocity of $2.2\pm1.2$kms$^{-1}$ at the half light radius. The rotation appears to be aligned with the major axis of Tucana. Despite this rotation signature, Tucana is still primarily a dispersion supported system.

-   The velocity dispersion of Tucana is found to be $\sigma_{\mathrm{v,Tuc}}=14.4_{-2.3}^{+2.8}$kms$^{-1}$. This dispersion is consistent with the result of F09, and suggests that Tucana is significantly more centrally dense than other dSphs. Tucana is found to be compatible with high density subhalos predicted by $\Lambda$CDM simulations, and hence becomes the first known exception to the too–big–to–fail problem. This proves that these ‘massive failure’ halos do exist in nature, confirming one of the key predictions of pure dark matter structure formation simulations in $\Lambda$CDM.

-   We use Jeans modelling to estimate the dark matter density profile of Tucana. Although the uncertainties are large, our results favour a high central density ($\rho_{\rm{DM}}$(150 pc)=$5.5\pm3.2\times10^8$ M$_{\odot}$ kpc$^{-3}$) and a halo mass $M_{200}=1.37^{+0.49}_{-0.44} \times 10^{10}$ M$_{\odot}$ consistent with abundance matching. Tucana’s old–age population distinguishes it from other isolated, gas–rich galaxies which are still forming stars today. In models where dark matter is ‘heated’ by baryonic feedback, Tucana is therefore expected to retain a higher central density than other isolated dwarfs. As anticipated from the lack of recent star formation, Tucana is consistent with residing in a ‘pristine’ dark matter halo, unaffected by dark matter heating (see Fig. \[fig:dm\_heating\]). Further spectroscopic follow–up, particularly in the poorly sampled central regions of the galaxy, would be required to confirm the presence of a potential cusp in the density profile.

Acknowledgements {#acknowledgements .unnumbered}
================

The authors would like to thank John Pritchard of ESO Operations Support for his assistance with the wavelength calibration of the spectra. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 095.B-0133(A). This paper includes data gathered with the 6.5m Magellan Telescopes located at Las Campanas Observatory, Chile. This work makes use of PyAstronomy and `emcee`. DRW is supported by a fellowship from the Alfred P. Sloan Foundation. NFM acknowledges that this work has been published under the framework of the IdEx Unistra and benefits from a funding from the state managed by the French National Research Agency as part of the investments for the future program. NFM also acknowledge support by the Programme National Cosmology et Galaxies (PNCG) of CNRS/INSU with INP and IN2P3, co–funded by CEA and CNES.

Modelling Tucana in SIDM {#appendixa}
========================

In this Appendix, we present the results of fitting Tucana with a <span style="font-variant:small-caps;">coreNFW</span> profile (see @read16 [@read18a] for details of the functional form of this profile). This allows us to assess whether our results are sensitive to our choice of mass model for the dark matter halo of Tucana, and has the advantage that one of <span style="font-variant:small-caps;">coreNFW</span> parameters is the halo mass, M$_{200}$. The results, shown in Fig. \[fig:tucana\_sidm\], are consistent with those of the free–form model (Fig. \[fig:gravsphere\_result\]). With <span style="font-variant:small-caps;">coreNFW</span>, we obtain a central density of $\rho_{\mathrm{DM}}$(150 pc)$=6.0_{-2.9}^{+3.7}\times 10^8$ M$_{\odot}$ kpc$^{-3}$ at the 68$\%$ confidence level. As was expected from @read18a, the density profile is systematically shallower in the innermost regions than with the free–form model. This effect is a result of the priors used, but is smaller than the measured uncertainties. The <span style="font-variant:small-caps;">coreNFW</span> model returns a pre–infall halo mass of $M_{200}=1.37^{+0.49}_{-0.44} \times 10^{10}$ M$_{\odot}$. Overall, the results of the two models are fully consistent with each other. This demonstrates that our inference of a high central density in Tucana is not dependent on our choice of mass model and priors.

![The density profile of Tucana, modelled using a <span style="font-variant:small-caps;">coreNFW</span> profile. Grey contours show the 1– and 2–$\sigma$ uncertainty ranges. The results are in excellent agreement with our default ‘free-form’ mass model (see Fig. \[fig:gravsphere\_result\]).[]{data-label="fig:tucana_sidm"}](tucana_paper_sidm.pdf){width="\linewidth"}

Mock Data Testing of <span style="font-variant:small-caps;">GravSphere</span> {#appendixb}
=============================================================================

![image](tucana_paper_mocks.pdf){width="\linewidth"}

In this appendix we outline a series of mock data tests designed to investigate the effect of using <span style="font-variant:small-caps;">GravSphere</span> to model a small stellar sample.

We generate a series of mock data sets designed to reproduce the characteristics of the real kinematic data for Tucana. Each mock generates 2500 stars to reproduce the radial distribution of sources, the distribution of their velocity errors and their CMD positions. We set up the tracer density profile and kinematics for the mock assuming a Plummer light profile for the stars, with a stellar mass of $M_*=0.56\times10^6$ M$_\odot$ and a scale length of $r_P=0.284$ kpc, and a <span style="font-variant:small-caps;">coreNFW</span> profile for the dark matter halo with $M_{200}=10^10$ M$_\odot$ and a concentration parameter chosen to be twice the M$_{200}$–c$_{200}$ relation from @dutton14 (to ensure a high velocity dispersion similar to that found in Tucana). The velocities for the stars were sampled from an isotropic distribution function generated using the AGAMA code [@vasiliev19]. Two mocks were constructed: one designed to represent a cusped galaxy; and one to represent a cored galaxy. The mocks include a foreground contribution mimicking any potential contamination of our real sample. Samples of $\sim$20, $\sim$100 and $\sim$500 effective members were then randomly drawn from the mock datasets. These samples represent the size of the true membership sample—with sampling designed to produce four membership–weighted radial bins as in the real Tucana model— and two levels of increased sampling to illustrate the improvement achieved with more data.

<span style="font-variant:small-caps;">GravSphere</span> was then used to model the density profile of these mock data. The velocity dispersion profile $\sigma_{\rm{LOS}}$ was derived from the mock velocity data. To generate the surface brightness profile $\Sigma_{*}(R)$, mock photometry was produced, simulating the CMD positions and radial distribution of the Magellan/ MegaCam imaging catalogue. The resulting density profiles are shown for both the cusped and cored mocks in Fig. \[fig:gravsphere\_mocks\]. The red lines mark the underlying profile used to seed the mocks. In the cusped case (bottom row), the profile is well recovered by <span style="font-variant:small-caps;">GravSphere</span>, with the constraints tightening as the sample size increases. For the lowest-sampled mock (left panels), the recovered density profile has uncertainties similar to that for the real Tucana data (compare with Fig. \[fig:gravsphere\_result\]). In the cored case (top row), the central density is reasonably well recovered (within 1–$\sigma$), but the density profile is slightly biased towards lower densities than the input profile at large radii, such that the model appears cuspier than the true profile. The slight bias towards cusped profiles is noted and explored further in @read18a. There, it was shown that the bias diminishes with improved sampling to 1000 or 2000 effective members (their Figure B1). This is something that we will improve on in future work. The central density at 150pc– which is of interest in this paper– is well recovered, supporting our use of <span style="font-variant:small-caps;">GravSphere</span> to measure $\rho_{\rm{DM}}$(150pc) for Tucana.

\[lastpage\]

[^1]: E-mail: a.gregory@surrey.ac.uk

[^2]: https://github.com/sczesla/PyAstronomy
