We present the results of N-body models of the Milky Way and Large Magellanic Cloud system and study the kinematic reflex motion in the stellar halo owing to the barycentre displacement of the disc. In agreement with previous studies, we find that the Milky Way disc may be moving at 40 km/s relative to the barycentre prior to the Large Magellanic Cloud infall. The resulting reflex motion is visible in tangential velocities of the stellar halo as a simple dipole. The signal is strongest for stars with long dynamical times, identifiable in position-velocity data as stars with large apocentres, whose dynamical memory is still well-represented by the unperturbed Milky Way potential. The signal varies across the sky depending on the stellar tracer and may be up to the same magnitude as the velocity of the disc centre-of-mass, making reflex motion a source of bias for Milky Way potential determinations based on the modeling of stellar streams and/or smooth halo tracers such as blue horizontal branch or RR Lyrae stars.

The clumpiness of dark matter on sub-kpc scales is highly sensitive to the tidal evolution and survival of subhaloes. In agreement with previous studies, we show that N-body realisations of cold dark matter subhaloes with centrally-divergent density cusps form artificial constant-density cores on the scale of the resolution limit of the simulation. These density cores drive the artificial tidal disruption of subhaloes. We run controlled simulations of the tidal evolution of a single subhalo where we repeatedly reconstruct the density cusp, preventing artificial disruption. This allows us to follow the evolution of the subhalo for arbitrarily large fractions of tidally stripped mass. Based on this numerical evidence in combination with simple dynamical arguments, we argue that cuspy dark matter subhaloes cannot be completely disrupted by smooth tidal fields. Modelling stars as collisionless tracers of the underlying potential, we furthermore study the tidal evolution of Milky Way dwarf spheroidal galaxies. Using a model of the Tucana III dwarf as an example, we show that tides can strip dwarf galaxies down to sub-solar luminosities. The remnant 'micro-galaxies' would appear as co-moving groups of metal-poor, low-mass stars of similar age, embedded in sub-kpc dark matter subhaloes.

This paper presents $N$-body and stochastic models that describe the motion of tracer particles in a potential that contains a large population of extended substructures. Fluctuations of the gravitational field induce a random walk of orbital velocities that is fully specified by drift and diffusion coefficients. In the impulse and local approximations the coefficients are computed analytically from the number density, mass, size and relative velocity of substructures without arbitrary cuts in forces or impact parameters. The resulting Coulomb logarithm attains a well-defined geometrical meaning, $\ln(\Lambda)=\ln (D/c)$, where $D/c$ is the ratio between the average separation and the individual size of substructures. Direct-force and Monte-Carlo $N$-body experiments show excellent agreement with the theory if substructures are sufficiently extended ($c/D\gtrsim 10^{-3}$) and not spatially overlapping ($c/D\lesssim 10^{-1}$). However, close encounters with point-like objects ($c/D\ll 10^{-3}$) induce a heavy-tailed, non-Gaussian distribution of high-energy impulses that cannot be described with Brownian statistics. In the point-mass limit ($c/D\approx 0$) the median Coulomb logarithm measured from $N$-body models deviates from the theoretical relation, converging towards a maximum value $\langle \ln(\Lambda)\rangle \approx 8.2$ independently of the mass and relative velocity of nearby substructures.

Large galaxies grow through the accumulation of dwarf galaxies. In principle it is possible to trace this growth history using the properties of a galaxy's stellar halo. Previous investigations of the galaxy M31 (Andromeda) have shown that outside a radius of 25 kpc the population of halo globular clusters is rotating in alignment with the stellar disk, as are more centrally located clusters. The M31 halo also contains coherent stellar substructures, along with a smoothly distributed stellar component. Many of the globular clusters outside 25 kpc are associated with the most prominent substructures, while others are part of the smooth halo. Here we report a new analysis of the kinematics of these globular clusters. We find that the two distinct populations are rotating with perpendicular orientations. The rotation axis for the population associated with the smooth halo is aligned with the rotation axis for the plane of dwarf galaxies that encircles M31. We interpret these separate cluster populations as arising from two major accretion epochs, likely separated by billions of years. Stellar substructures from the first epoch are gone, but those from the more recent second epoch still remain.

We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS survey and enables us to confirm that Andromeda XXVII's heliocentric distance is 827 $\pm$ 47 kpc and spectroscopic metallicity is -2.1$^{+0.4}_{-0.5}$. We also re-derive Andromeda XXVII's kinematic properties, measuring a systemic velocity = -526.1$^{+10.0}_{-11.0}${\kms} and a velocity dispersion that we find to be non-Gaussian but for which we derive a formal value of 27.0$^{+2.2}_{-3.9}${\kms}. In the upper segment of the North West Stream we measure mean values for the metallicity = -1.8$\pm$0.4, systemic velocity = -519.4 $\pm$4.0{\kms} and velocity dispersion = 10.0$\pm$4.0{\kms}. We also detect a velocity gradient of 1.7$\pm$0.3 {\kms} kpc$^{-1}$ on an infall trajectory towards M31. With a similar gradient, acting in the same direction, in the lower segment we suggest that the North West Stream is not a single structure. As the properties of the upper segment of the North West Stream and Andromeda XXVII are consistent within ~90\% confidence limits, it is likely that the two are related and plausible that Andromeda XXVII is the progenitor of this stream.

We propose that the Ultra-Diffuse Galaxy (UDG) population represents a set of satellite galaxies born in $\sim10^{10}-10^{11}$ M$_\odot$ halos, similar to field dwarfs, which suffer a dramatic reduction in surface brightness due to tidal stripping and heating. This scenario is observationally motivated by the radial alignment of UDGs in Coma as well as the significant dependence of UDG abundance on cluster mass. As a test of this formation scenario, we apply a semi-analytic model describing the change in stellar mass and half-light radius of dwarf satellites, occupying either cored or cuspy halos, to cluster subhalos in the Illustris-dark simulation. Key to this model are results from simulations which indicate that galaxies in cored dark-matter halos expand significantly in response to tidal stripping and heating, whereas galaxies in cuspy halos experience limited size evolution. Our analysis indicates that a population of tidally-stripped dwarf galaxies, residing in cored halos (like those hosting low-surface brightness field dwarfs), is able to reproduce the observed sizes and stellar masses of UDGs in clusters remarkably well.

Gravitating systems surrounded by a dynamic sea of substructures experience fluctuations of the local tidal field which inject kinetic energy into the internal motions. This paper uses stochastic calculus techniques to describe `tidal heating' as a random walk of orbital velocities that leads to diffusion in a 4-dimensional energy--angular momentum space. In spherical, static potentials we derive analytical solutions for the Green's propagators directly from the number density and velocity distribution of substructures with known mass & size functions without arbitrary cuts in forces or impact parameters. Furthermore, a Monte-Carlo method is presented, which samples velocity 'kicks' from a probability function and can be used to model orbital scattering in fully generic potentials. For illustration, we follow the evolution of planetary orbits in a clumpy environment. We show that stochastic heating of (mass-less) discs in a Keplerian potential leads to the formation, and subsequent `evaporation' of Oort-like clouds, and derive analytical expressions for the escape rate and the fraction of comets on retrograde orbits as a function of time. Extrapolation of the subhalo mass function of Milky Way-like haloes down to the WIMP free-streaming length suggests that objects in the outer Solar system experience repeated interactions with dark microhaloes on dynamical time-scales.

Mass estimators are a key tool to infer the dark matter content in pressure-supported systems like dwarf spheroidal galaxies (dSphs). We construct an estimator for enclosed masses based on the virial theorem which is insensitive to anisotropy in the velocity dispersion and tailored to yield masses with minimum uncertainty introduced by our ignorance on (i) the shape of the inner halo profile, and (ii) how deeply the stellar component is embedded within the halo: $M(<1.8\,R_\mathrm{h}) \approx 3.5 \times 1.8\,R_\mathrm{h} \langle \sigma_\mathrm{los}^2 \rangle G^{-1}$, where by $R_\mathrm{h}$ we denote the projected half-light radius and by $\langle \sigma_\mathrm{los}^2 \rangle$ the luminosity-averaged squared line-of-sight velocity dispersion. Tests against controlled simulations show that this estimator provides unbiased enclosed masses with an accuracy of $\sim 10$ per cent. This confirms the robustness of similar previously proposed mass estimators. Application to published kinematic data of Milky Way dSphs reveals a tight correlation between enclosed mass and luminosity. Using $N$-body models we show that tidal stripping has little effect on this relation. Comparison against cuspy and cored dark matter haloes extracted from controlled re-simulations of the Aquarius A2 merger tree shows that the high mass densities of ultrafaint galaxies are not compatible with large dark matter cores, and that the (total) halo masses of the classical Milky Way dSphs span a remarkably narrow range ($8 \lesssim \mathrm{log_{10}}\,(M/\mathrm{M_\odot}) \lesssim 10$) at present, showing no clear trend with either galaxy size or luminosity.

We present a new technique to probe the central dark matter (DM) density profile of galaxies that harnesses both the survival and observed properties of star clusters. As a first application, we apply our method to the `ultra-faint' dwarf Eridanus II (Eri II) that has a lone star cluster ~45 pc from its centre. Using a grid of collisional $N$-body simulations, incorporating the effects of stellar evolution, external tides and dynamical friction, we show that a DM core for Eri II naturally reproduces the size and the projected position of its star cluster. By contrast, a dense cusped galaxy requires the cluster to lie implausibly far from the centre of Eri II (>1 kpc), with a high inclination orbit that must be observed at a particular orbital phase. Our results, therefore, favour a dark matter core. This implies that either a cold DM cusp was `heated up' at the centre of Eri II by bursty star formation, or we are seeing an evidence for physics beyond cold DM.

A large population of extended substructures generates a stochastic gravitational field that is fully specified by the function $p({\bf F})$, which defines the probability that a tracer particle experiences a force $\bf F$ within the interval ${\bf F},{\bf F}+ d\bf F$. This paper presents a statistical technique for deriving the spectrum of random fluctuations directly from the number density of substructures with known mass and size functions. Application to the subhalo population found in cold dark matter simulations of Milky Way-sized haloes shows that, while the combined force distribution is governed by the most massive satellites, the fluctuations of the {\it tidal} field are completely dominated by the smallest and most abundant subhaloes. In light of this result we discuss observational experiments that may be sufficiently sensitive to Galactic tidal fluctuations to probe the "dark" low-end of the subhalo mass function and constrain the particle mass of warm and ultra-light axion dark matter models.

It has been suggested that the internal dynamics of dwarf spheroidal galaxies (dSphs) can be used to test whether or not ultralight axions with $m_a\sim 10^{-22}\text{eV}$ are a preferred dark matter candidate. However, comparisons to theoretical predictions tend to be inconclusive for the simple reason that while most cosmological models consider only dark matter, one observes only baryons. Here we use realistic kinematic mock data catalogs of Milky Way dSph's to show that the "mass-anisotropy degeneracy" in the Jeans equations leads to biased bounds on the axion mass in galaxies with unknown dark matter halo profiles. In galaxies with multiple chemodynamical components this bias can be partly removed by modelling the mass enclosed within each subpopulation. However, analysis of the mock data reveals that the least-biased constraints on the axion mass result from fitting the luminosity-averaged velocity dispersion of the individual chemodynamical components directly. Applying our analysis to two dSph's with reported stellar subcomponents, Fornax and Sculptor, and assuming that the halo profile has not been acted on by baryons, yields core radii $r_{c}>1.5$ kpc and $r_c> 1.2$ kpc respectively, and $m_a<0.4\times 10^{-22}\text{eV}$ at 97.5\% confidence. These bounds are in tension with the number of observed satellites derived from simple (but conservative) estimates of the subhalo mass function in Milky Way-like galaxies. We discuss how baryonic feedback might affect our results, and the impact of such a small axion mass on the growth of structures in the Universe.

We show that hard encounters in the central regions of globular clusters embedded in dark matter (DM) haloes necessarily lead to the formation of gravitationally-bound stellar envelopes that extend far beyond the nominal tidal radius of the system. Using statistical arguments and numerical techniques we derive the equilibrium distribution function of stars ejected from the centre of a non-divergent spherical potential. Independently of the velocity distribution with which stars are ejected, GC envelopes have density profiles that approach asymptotically $\rho\sim r^{-4}$ at large distances and become isothermal towards the centre. Adding a DM halo component leaves two clear-cut observational signatures: (i) a flattening, or slightly increase of the projected velocity dispersion profile at large distances, and (ii) an outer surface density profile that is systematically shallower than in models with no dark matter.

The parameterised post-Newtonian (PPN) formalism has enabled stringent tests of static weak-field gravity in a theory-independent manner. Here we incorporate screening mechanisms of modified gravity theories into the framework by introducing an effective gravitational coupling and defining the PPN parameters as functions of position. To determine these functions we develop a general method for efficiently performing the post-Newtonian expansion in screened regimes. For illustration, we derive all the PPN functions for a cubic galileon and a chameleon model. We also analyse the Shapiro time delay effect for these two models and find no deviations from General Relativity insofar as the signal path and the perturbing mass reside in a screened region of space.

We use 12 cosmological $N$-body simulations of Local Group systems (the Apostle models) to inspect the relation between the virial mass of the main haloes ($M_{\rm vir,1}$ and $M_{\rm vir,2}$), the mass derived from the relative motion of the halo pair ($M_{\rm tim}$), and that inferred from the local Hubble flow ($M_{\rm lhf}$). We show that within the Spherical Collapse Model (SCM), the correspondence between the three mass estimates is exact, i.e. $M_{\rm lhf}=M_{\rm tim}=M_{\rm vir,1}+M_{\rm vir,2}$. However, comparison with Apostle simulations reveals that, contrary to what the SCM states, a relatively large fraction of the mass that perturbs the local Hubble flow and drives the relative trajectory of the main galaxies is not contained within $R_{\rm vir}$, and that the amount of "extra-virial" mass tends to increase in galaxies with a slow accretion rate. In contrast, modelling the peculiar velocities around the Local Group returns an unbiased constraint on the virial mass ratio of the main galaxy pair. Adopting the outer halo profile found in $N$-body simulations, which scales as $\rho\sim R^{-4}$ at $R\gtrsim R_{\rm vir}$, indicates that the galaxy masses perturbing the local Hubble flow roughly correspond to the asymptotically-convergent (total) masses of the individual haloes. We show that estimates of $M_{\rm vir}$ based on the dynamics of tracers at $R\gg R_{\rm vir}$ require a priori information on the internal matter distribution and the growth rate of the main galaxies, both of which are typically difficult to quantify.

We use high-resolution $N$-body simulations to study the effect of a galactic disc on the dynamical evolution of dark matter substructures with orbits and structural parameters extracted from the Aquarius A-2 merger tree (Springel et al. 2008). Satellites are modelled as equilibrium $N$-body realizations of generalized Hernquist profiles with $2\times10^6$ particles and injected in the analytical evolving host potential at $z_\mathrm{infall}$, defined by the peak of their mass evolution. We select all substructures with $M_{200}(z_\mathrm{infall})\geq 10^8\,\mathrm{M_\odot}$ and first pericentric distances $r_p<r_{200}$. Motivated by observations of Milky Way dwarf spheroidal galaxies, we also explore satellite models with cored dark matter profiles with a fixed core size $r_c=0.8\,a_s$ where $a_s$ is the Hernquist scale radius. We find that models with cuspy satellites have twice as many surviving substructures at $z=0$ than their cored counterparts, and four times as many if we only consider those on orbits with $r_p\lesssim0.1\,r_{200}$. For a given profile, adding an evolving disc potential reduces the number of surviving substructures further by a factor of $\lesssim2$ for satellites on orbits that penetrate the disc ($r_p\lesssim 20\,\mathrm{kpc}$). For large $r_p$, where tidal forces and the effect of the disc become negligible, the number of satellites per pericentre bin converges to similar values for all four models.

The Horndeski action is the most general scalar-tensor theory with at most second-order derivatives in the equations of motion, thus evading Ostrogradsky instabilities and making it of interest when modifying gravity at large scales. To pass local tests of gravity, these modifications predominantly rely on nonlinear screening mechanisms that recover Einstein's Theory of General Relativity in regions of high density. We derive a set of conditions on the four free functions of the Horndeski action that examine whether a specific model embedded in the action possesses an Einstein gravity limit or not. For this purpose, we develop a new and surprisingly simple scaling method that identifies dominant terms in the equations of motion by considering formal limits of the couplings that enter through the new terms in the modified action. This enables us to find regimes where nonlinear terms dominate and Einstein's field equations are recovered to leading order. Together with an efficient approximation of the scalar field profile, one can then further evaluate whether these limits can be attributed to a genuine screening effect. For illustration, we apply the analysis to both a cubic galileon and a chameleon model as well as to Brans-Dicke theory. Finally, we emphasise that the scaling method also provides a natural approach for performing post-Newtonian expansions in screened regimes.

We carry out controlled $N$-body simulations that follow the dynamical evolution of binary stars in the dark matter (DM) haloes of ultra-faint dwarf spheroidals (dSphs). We find that wide binaries with semi-major axes $a\gtrsim a_t$ tend to be quickly disrupted by the tidal field of the halo. In smooth potentials the truncation scale, $a_t$, is mainly governed by (i) the mass enclosed within the dwarf half-light radius ($R_h$) and (ii) the slope of the DM halo profile at $R\approx R_h$, and is largely independent of the initial eccentricity distribution of the binary systems and the anisotropy of the stellar orbits about the galactic potential. For the reported velocity dispersion and half-light radius of Segue I, the closest ultra-faint, our models predict $a_t$ values that are a factor 2--3 smaller in cuspy haloes than in any of the cored models considered here. Using mock observations of Segue I we show that measuring the projected two-point correlation function of stellar pairs with sub-arcsecond resolution may provide a useful tool to constrain the amount and distribution of DM in the smallest and most DM-dominated galaxies.

This paper explores the effect of the LMC on the mass estimates obtained from the timing argument. We show that accounting for the presence of the LMC systematically lowers the Local Group mass ($M_{\rm LG}$) derived from the relative motion of the Milky Way--Andromeda pair. Motivated by this result we apply a Bayesian technique devised by Pe\~narrubia et al. (2014) to simultaneously fit (i) distances and velocities of galaxies within 3~Mpc and (ii) the relative motion between the Milky Way and Andromeda derived from HST observations, with the LMC mass ($M_{\rm LMC}$) as a free parameter. Our analysis returns a Local Group mass $M_{\rm LG}=2.64^{+0.42}_{-0.38}\times 10^{12}M_\odot$ at a 68\% confidence level. The masses of the Milky Way, $M_{\rm MW}=1.04_{-0.23}^{+0.26}\times 10^{12}M_\odot$, and Andromeda, $M_{\rm M31}=1.33_{-0.33}^{+0.39}\times 10^{12}M_\odot$, are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find a (total) LMC mass $M_{\rm LMC}=0.25_{-0.08}^{+0.09}\times 10^{12}M_\odot$, which is indicative of an extended dark matter halo and supports the scenario where this galaxy is just past its first pericentric approach. Consequently, these results suggest that the LMC may induce significant perturbations on the Galactic potential.

This paper uses dynamical invariants to describe the evolution of collisionless systems subject to time-dependent gravitational forces without resorting to maximum-entropy probabilities. We show that collisionless relaxation can be viewed as a special type of diffusion process in the integral-of-motion space. In time-varying potentials with a fixed spatial symmetry the diffusion coefficients are closely related to virial quantities, such as the specific moment of inertia, the virial factor and the mean kinetic and potential energy of microcanonical particle ensembles. The non-equilibrium distribution function (DF) is found by convolving the initial DF with the Green function that solves Einstein's equation for freely diffusing particles. Such a convolution also yields a natural solution to the Fokker-Planck equations in the energy space. Our mathematical formalism can be generalized to potentials with a time-varying symmetry, where diffusion extends over multiple dimensions of the integral-of-motion space. The new probability theory is in many ways analogous to stochastic calculus, with two significant differences: (i) the equations of motion that govern the trajectories of particles are fully deterministic, and (ii) the diffusion coefficients can be derived self-consistently from microcanonical phase-space averages without relying on ergodicity assumptions. For illustration we follow the cold collapse of $N$-body models in a time-dependent logarithmic potential. Comparison between the analytical and numerical results shows excellent agreement in regions where the potential evolution does not depart too strongly from the adiabatic regime.

We perform in-depth dynamical modelling of the luminous and dark matter (DM) content of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equations for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM density profiles to be $\gamma = 0.6\pm0.4$, which is consistent with either a DM cusp ($\gamma = 1$) or with a DM core $(\gamma = 0)$. Our findings are consistent with a Salpeter IMF, and marginally consistent with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red GCs and stars. The modelling results are consistent with the virial mass--concentration relation predicted by $\Lambda$CDM simulations. The virial mass of NGC 1407 is $\log$ $M_{\rm vir} = 13.3 \pm 0.2 M_{\odot}$, whereas the stellar mass is $\log M_* = 11.8 \pm 0.1 M_{\odot}$. The overall uncertainties on the mass of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the disagreement between our results and previous X-ray results to the gas not being in hydrostatic equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM dominated, with a dynamical mass-to-light ratio of $M/L=260_{-100} ^{+174} M_{\odot}/L_{\odot, B}$. However, this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale radius.

We use high-resolution N-body simulations to follow the formation and evolution of tidal streams associated to dwarf spheroidal galaxies (dSphs). The dSph models are embedded in dark matter (DM) haloes with either a centrally-divergent 'cusp', or an homogeneous-density 'core'. In agreement with previous studies, we find that as tides strip the galaxy the evolution of the half-light radius and the averaged velocity dispersion follows well-defined tracks that are mainly controlled by the amount of mass lost. Crucially, the evolutionary tracks behave differently depending on the shape of the DM profile: at a fixed remnant mass, dSphs embedded in cored haloes have larger sizes and higher velocity dispersions than their cuspy counterparts. The divergent evolution is particularly pronounced in galaxies whose stellar component is strongly segregated within their DM halo and becomes more disparate as the remnant mass decreases. Our analysis indicates that the DM profile plays an important role in defining the internal dynamics of tidal streams. We find that stellar streams associated to cored DM models have velocity dispersions that lie systematically above their cuspy counterparts. Our results suggest that the dynamics of streams with known dSph progenitors may provide strong constraints on the distribution of DM on the smallest galactic scales.

We run controlled N-body experiments to study the evolution of the dark matter (DM) halo profiles of dwarf galaxies driven by the accretion of DM substructures. Our initial conditions assume that supernova feedback erases the primordial DM cusps of haloes with $10^{9}-10^{10} \rm{M_{\odot}}$ at $z=0$. The orbits and masses of the infalling substructures are borrowed from the {\it Aquarius} simulations. Our experiments show that a fraction of haloes that undergo 1:3 down to 1:30 mergers are susceptible to reform a DM cusp by $z\approx 0$. Cusp regrowth is driven the accretion of DM substructures that are dense enough reach the central regions of the main halo before being tidally disrupted. The infall of substructures with a mass ratio above 1:6 on the mean of the reported mass-concentration relation systematically lead to cusp regrowth. Between 1:6 to 1:8, and 1:8 to 1:30 substructures need to be located one and two-sigma above the mean, respectively. The merging timescales of these dense, low-mass substructures is relatively long $(5-11 \rm{Gyrs})$, which may pose a timescale problem for the longevity of DM cores in dwarf galaxies. These results suggest that a certain level of scatter in the central density slopes of galactic haloes acted-on by feedback is to be expected given the stochastic mass accretion histories of low-mass haloes and the diverse star formation histories observed in the Local Group dwarves.

The presence of a scalar field that couples nonminimally and universally to matter can enhance gravitational forces on cosmological scales while restoring general relativity in the Solar neighborhood. In the intermediate regime, kinematically inferred masses experience an additional radial dependence with respect to the underlying distribution of matter, which is caused by the increment of gravitational forces with increasing distance from the Milky Way center. The same effect can influence the internal kinematics of subhalos and cause cuspy matter distributions to appear core-like. Specializing to the chameleon model as a worked example, we demonstrate this effect by tracing the scalar field from the outskirts of the Milky Way halo to its interior, simultaneously fitting observed velocity dispersions of chemo-dynamically discriminated red giant populations in the Fornax and Sculptor dwarf spheroidals. Whereas in standard gravity these observations suggest that the matter distribution of the dwarfs is cored, we find that in the presence of a chameleon field the assumption of a cuspy Navarro-Frenk-White profile becomes perfectly compatible with the data. Importantly, chameleon models also predict the existence of slopes between two stellar subcomponents that in Newtonian gravity would be interpreted as a depletion of matter in the dwarf center. Hence, an observation of such an apparently pathological scenario may serve as a smoking gun for the presence of a chameleon field or a similar modification of gravity, independent of baryonic feedback effects. In general, measuring the dynamic mass profiles of the Milky Way dwarfs provides stronger constraints than those inferred from the screening scale of the Solar System since these are located at greater distances from the halo center.

We combine the equations of motion that govern the dynamics of galaxies in the local volume with Bayesian techniques in order to fit orbits to published distances and velocities of galaxies within $\sim 3$ Mpc. We find a Local Group (LG) mass $2.3\pm 0.7\times 10^{12}{\rm M}_\odot$ that is consistent with the combined dynamical masses of M31 and the Milky Way, and a mass ratio $0.54^{+0.23}_{-0.17}$ that rules out models where our Galaxy is more massive than M31 with $\sim 95\%$ confidence. The Milky Way's circular velocity at the solar radius is relatively high, $245\pm 23$ km/s, which helps to reconcile the mass derived from the local Hubble flow with the larger value suggested by the `timing argument'. Adopting {\it Planck}'s bounds on $\Omega_\Lambda$ yields a (local) Hubble constant $H_0=67\pm 5$km/s/Mpc which is consistent with the value found on cosmological scales. Restricted N-body experiments show that substructures tend to fall onto the LG along the Milky Way-M31 axis, where the quadrupole attraction is maximum. Tests against mock data indicate that neglecting this effect slightly overestimates the LG mass without biasing the rest of model parameters. We also show that both the time-dependence of the LG potential and the cosmological constant have little impact on the observed local Hubble flow.

We present a panoptic view of the stellar structure in the Galactic disk's outer reaches commonly known as the Monoceros Ring, based on data from Pan-STARRS1. These observations clearly show the large extent of the stellar overdensities on both sides of the Galactic disk, extending between b = -25 and b = +35 degrees and covering over 130 degrees in Galactic longitude. The structure exhibits a complex morphology with both stream-like features and a sharp edge to the structure in both the north and the south. We compare this map to mock observations of two published simulations aimed at explaining such structures in the outer stellar disk, one postulating an origin as a tidal stream and the other demonstrating a scenario where the disk is strongly distorted by the accretion of a satellite. These morphological comparisons of simulations can link formation scenarios to observed structures, such as demonstrating that the distorted-disk model can produce thin density features resembling tidal streams. Although neither model produces perfect agreement with the observations--the tidal stream predicts material at larger distances which is not detected while in the distorted disk model the midplane is warped to an excessive degree--future tuning of the models to accommodate these latest data may yield better agreement.

We report the discovery of a thin stellar stream found in Pan-STARRS1 photometry near the Galactic bulge in the constellation of Ophiuchus. It appears as a coherent structure in the colour-selected stellar density maps produced to search for tidal debris around nearby globular clusters. The stream is exceptionally short and narrow; it is about 2.5{\deg} long and 6' wide in projection. The colour-magnitude diagram of this object, which harbours a blue horizontal-branch, is consistent with an old and relatively metal-poor population ([Fe/H]~-1.3) located 9.5 +/- 0.9 kpc away at (l,b) ~ (5{\deg},+32{\deg}), and 5.0 +/- 1.0 kpc from the Galactic centre. These properties argue for a globular cluster as progenitor. The finding of such a prominent, nearby stream suggests that many streams could await discovery in the more densely populated regions of our Galaxy.

We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of \Sigma_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances.

Using a variety of stellar tracers -- blue horizontal branch stars, main-sequence turn-off stars and red giants -- we follow the path of the Sagittarius (Sgr) stream across the sky in Sloan Digital Sky Survey data. Our study presents new Sgr debris detections, accurate distances and line-of-sight velocities that together help to shed new light on the puzzle of the Sgr tails. For both the leading and the trailing tail, we trace the points of their maximal extent, or apo-centric distances, and find that they lie at $R^L$ = 47.8 $\pm$ 0.5 kpc and $R^T$ = 102.5 $\pm$ 2.5 kpc respectively. The angular difference between the apo-centres is 93.2 $\pm$ 3.5 deg, which is smaller than predicted for logarithmic haloes. Such differential orbital precession can be made consistent with models of the Milky Way in which the dark matter density falls more quickly with radius. However, currently, no existing Sgr disruption simulation can explain the entirety of the observational data. Based on its position and radial velocity, we show that the unusually large globular cluster NGC 2419 can be associated with the Sgr trailing stream. We measure the precession of the orbital plane of the Sgr debris in the Milky Way potential and show that, surprisingly, Sgr debris in the primary (brighter) tails evolves differently to the secondary (fainter) tails, both in the North and the South.

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

This paper explores a mathematical technique for deriving dynamical invariants (i.e. constants of motion) in time-dependent gravitational potentials. The method relies on the construction of a canonical transformation that removes the explicit time-dependence from the Hamiltonian of the system. By referring the phase-space locations of particles to a coordinate frame in which the potential remains `static' the dynamical effects introduced by the time evolution vanish. It follows that dynamical invariants correspond to the integrals of motion for the static potential expressed in the transformed coordinates. The main difficulty of the method reduces to solving the differential equations that define the canonical transformation, which are typically coupled with the equations of motion. We discuss a few examples where both sets of equations can be exactly de-coupled, and cases that require approximations. The construction of dynamical invariants has far-reaching applications. These quantities allow us, for example, to describe the evolution of (statistical) microcanonical ensembles in time-dependent gravitational potentials without relying on ergodicity or probability assumptions. As an illustration, we follow the evolution of dynamical fossils in galaxies that build up mass hierarchically. It is shown that the growth of the host potential tends to efface tidal substructures in the integral-of-motion space through an orbital diffusion process. The inexorable cycle of deposition, and progressive dissolution, of tidal clumps naturally leads to the formation of a `smooth' stellar halo.

We present the first kinematic analysis of the far outer halo globular cluster (GC) population in the Local Group galaxy M31. Our sample contains 53 objects with projected radii of ~20-130 kpc, of which 44 have no previous spectroscopic information. GCs with projected radii >30 kpc are found to exhibit net rotation around the minor axis of M31, in the same sense as the inner GCs, albeit with a smaller amplitude of 79 +/-19 km/s. The rotation-corrected velocity dispersion of the full halo GC sample is 106 +/-12 km/s, which we observe to decrease with increasing projected radius. We find compelling evidence for kinematic-coherence amongst GCs which project on top of halo substructure, including a clear signature of infall for GCs lying along the North-West stream. Using the tracer mass estimator, we estimate the dynamical mass of M31 within 200 kpc to be M_M31 = (1.2-1.5) +/- 0.2 x 10^12 M_sun. This value is highly dependent on the chosen model and assumptions within.

We generate stellar distribution functions (DFs) in triaxial haloes in order to examine the reliability of slopes $\Gamma\equiv \Delta {\rm log} M / \Delta {\rm log} r$ inferred by applying mass estimators of the form $M\propto R_e\sigma^2$ (i.e. assuming spherical symmetry, where $R_e$ and $\sigma$ are luminous effective radius and global velocity dispersion, respectively) to two stellar sub-populations independently tracing the same gravitational potential. The DFs take the form $f(E)$, are dynamically stable, and are generated within triaxial potentials corresponding directly to subhaloes formed in cosmological dark-matter-only simulations of Milky Way and galaxy cluster haloes. Additionally, we consider the effect of different tracer number density profiles (cuspy and cored) on the inferred slopes of mass profiles. For the isotropic DFs considered here, we find that halo triaxiality tends to introduce an anti-correlation between $R_e$ and $\sigma$ when estimated for a variety of viewing angles. The net effect is a negligible contribution to the systematic error associated with the slope of the mass profile, which continues to be dominated by a bias toward greater overestimation of masses for more-concentrated tracer populations. We demonstrate that simple mass estimates for two distinct tracer populations can give reliable (and cosmologically meaningful) lower limits for $\Gamma$, irrespective of the degree of triaxiality or shape of the tracer number density profile.

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 Andromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS spectrographs. Based on their g-i colors (taken with the CFHT MegaCam imager), physical positions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed star. Using this information, the velocity dispersions, central masses and central densities of the dark matter halos are calculated for these objects, and compared with the properties of the Milky Way dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, all of which have large half-light radii (>700 pc) and low velocity dispersions (sigma_v<5 km/s). In addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]_{half}=10.3^{+7.0}_{-6.7}, making it consistent with a simple stellar system with no appreciable dark matter component within its 1 sigma uncertainties. We suggest that the structure of the dark matter halos of these outliers have been significantly altered by tides.

We calculate the energy that baryons must inject in cold dark matter (CDM) haloes in order to remove centrally-divergent DM cusps on scales relevant to observations of dwarf spheroidal galaxies (dSphs). We estimate that the CDM haloes often associated with the Milky Way's dSphs (M_vir/M_\odot \sim 10^{9-10}) require \Delta E/erg \sim 10^{53-55} in order to form cores on scales comparable to the luminous size of these galaxies. While supernova type II (SNeII) explosions can in principle generate this energy, the actual contribution is limited by the low star formation efficiency implied by the abundance of luminous satellites. Considering that CDM's well-known `core/cusp' and `missing satellite' problems place opposing demands on star formation efficiencies, existing observational evidences for large cores in the most luminous dSphs require that CDM models invoke some combination of the following: (i) efficient (of order unity) coupling of SNeII energy into dark matter particles, (ii) star formation histories peaking at unexpectedly high redshifts (z>6), (iii) a top-heavy stellar IMF, and/or (iv) substantial satellite disruption or other stochastic effects to ease the substructure abundance constraints. Our models show that the tension between CDM problems on small scales would increase if cored DM profiles were to be found in fainter dwarves.

We introduce the Minimum Entropy Method, a simple statistical technique for constraining the Milky Way gravitational potential and simultaneously testing different gravity theories directly from 6D phase-space surveys and without adopting dynamical models. We demonstrate that orbital energy distributions that are separable (i.e. independent of position) have an associated entropy that increases under wrong assumptions about the gravitational potential and/or gravity theory. Of known objects, `cold' tidal streams from low-mass progenitors follow orbital distributions that most nearly satisfy the condition of separability. Although the orbits of tidally stripped stars are perturbed by the progenitor's self-gravity, systematic variations of the energy distribution can be quantified in terms of the cross-entropy of individual tails, giving further sensitivity to theoretical biases in the host potential. The feasibility of using the Minimum Entropy Method to test a wide range of gravity theories is illustrated by evolving restricted N-body models in a Newtonian potential and examining the changes in entropy introduced by Dirac, MONDian and f(R) gravity modifications.

We build a sample of distant (D > 80 kpc) stellar halo stars with measured radial velocities. Faint ($20 < g <22$) candidate blue horizontal branch (BHB) stars were selected using the deep, but wide, multi-epoch Sloan Digital Sky Survey photometry. Follow-up spectroscopy for these A-type stars was performed using the VLT-FORS2 instrument. We classify stars according to their Balmer line profiles, and find 7 are bona fide BHB stars and 31 are blue stragglers (BS). Owing to the magnitude range of our sample, even the intrinsically fainter BS stars can reach out to D ~ 90 kpc. We complement this sample of A-type stars with intrinsically brighter, intermediate-age, asymptotic giant branch stars. A set of 4 distant cool carbon stars is compiled from the literature and we perform spectroscopic follow-up on a further 4 N-type carbon stars using the WHT-ISIS instrument. Altogether, this provides us with the largest sample to date of individual star tracers out to r ~ 150 kpc. We find that the radial velocity dispersion of these tracers falls rapidly at large distances and is surprisingly cold (sigma_r ~ 50-60 km/s) between 100-150 kpc. Relating the measured radial velocities to the mass of the Milky Way requires knowledge of the (unknown) tracer density profile and anisotropy at these distances. Nonetheless, by assuming the stellar halo stars between 50-150 kpc have a moderate density fall-off (with power-law slope alpha < 5) and are on radial orbits (sigma^2_t/sigma^2_r < 1), we infer that the mass within 150 kpc is less than 10^12 M_solar and suggest it probably lies in the range (5-10) x 10^11 M_solar. We discuss the implications of such a low mass for the Milky Way.

Dwarf galaxies are widely believed to be among the best targets for indirect dark matter searches using high-energy gamma rays; and indeed gamma-ray emission from these objects has long been a subject of detailed study for ground-based atmospheric Cherenkov telescopes. Here, we update current exclusion limits obtained on the closest dwarf, the Sagittarius dwarf galaxy, in light of recent realistic dark matter halo models. The constraints on the velocity-weighted annihilation cross section of the dark matter particle are of a few 10$^{-23}$ cm$^{3}$s$^{-1}$ in the TeV energy range for a 50 h exposure. The limits are extrapolated to the sensitivities of future Cherenkov Telescope Arrays. For 200 h of observation time, the sensitivity at 95% C.L. reaches 10$^{-25}$ cm$^{3}$s$^{-1}$. Possible astrophysical backgrounds from gamma-ray sources dissembled in Sagittarius dwarf are studied. It is shown that with long-enough observation times, gamma-ray background from millisecond pulsars in a globular cluster contained within Sagittarius dwarf may limit the sensitivity to dark matter annihilations.

Using deep photometric data from WFC@INT and WFI@ESO2.2m we measure the outer number density profiles of 19 stellar clusters located in the inner region of the Milky Way halo (within a Galactocentric distance range of 10-30 kpc) in order to assess the impact of internal and external dynamical processes on the spatial distribution of stars. Adopting power-law fitting templates, with index $-\gamma$ in the outer region, we find that the clusters in our sample can be divided in two groups: a group of massive clusters ($ \ge 10^5 $ M_sun) that has relatively flat profiles with $2.5 < \gamma < 4$ and a group of low-mass clusters ($ \le 10^5 $ M_sun), with steep profiles ($\gamma > 4$) and clear signatures of interaction with the Galactic tidal field. We refer to these two groups as 'tidally unaffected' and 'tidally affected', respectively. Our results also show a clear trend between the slope of the outer parts and the half-mass density of these systems, which suggests that the outer density profiles may retain key information on the dominant processes driving the dynamical evolution of Globular Clusters.

We introduce a method for measuring the slopes of mass profiles within dwarf spheroidal (dSph) galaxies directly from stellar spectroscopic data and without adopting a dark matter halo model. Our method combines two recent results: 1) spherically symmetric, equilibrium Jeans models imply that the product of halflight radius and (squared) stellar velocity dispersion provides an estimate of the mass enclosed within the halflight radius of a dSph stellar component, and 2) some dSphs have chemo-dynamically distinct stellar \textit{sub}components that independently trace the same gravitational potential. We devise a statistical method that uses measurements of stellar positions, velocities and spectral indices to distinguish two dSph stellar subcomponents and to estimate their individual halflight radii and velocity dispersions. For a dSph with two detected stellar subcomponents, we obtain estimates of masses enclosed at two discrete points in the same mass profile, immediately defining a slope. Applied to published spectroscopic data, our method distinguishes stellar subcomponents in the Fornax and Sculptor dSphs, for which we measure slopes $\Gamma\equiv \Delta \log M / \Delta \log r=2.61_{-0.37}^{+0.43}$ and $\Gamma=2.95_{-0.39}^{+0.51}$, respectively. These values are consistent with 'cores' of constant density within the central few-hundred parsecs of each galaxy and rule out `cuspy' Navarro-Frenk-White (NFW) profiles ($d\log M/d\log r \leq 2$ at all radii) with significance $\ga 96%$ and $\ga 99%$, respectively. Tests with synthetic data indicate that our method tends systematically to overestimate the mass of the inner stellar subcomponent to a greater degree than that of the outer stellar subcomponent, and therefore to underestimate the slope $\Gamma$ (implying that the stated NFW exclusion levels are conservative).

While the satellites of the Milky Way (MW) have been shown to be largely consistent in terms of their mass contained within one half--light radius (M_{half}) with a "Universal" mass profile, a number of M31 satellites are found to be inconsistent with such relations, and seem kinematically colder in their central regions than their MW cousins. In this work, we present new kinematic and updated structural properties for two M31 dSphs, And V and And VI using data from the Keck Low Resolution Imaging Spectrograph (LRIS) and the DEep Imaging Multi-Object Spectrograph (DEIMOS) instruments and the Subaru Suprime-Cam imager. We measure systemic velocities of v_r=-393.1+/-4.2km/s and -344.8+/-2.5km/s, and dispersions of sigma_v=11.5{+5.3}{-4.4}km/s and sigma_v=9.4{+3.2}{-2.4}km/s for And V and And VI respectively, meaning these two objects are consistent with the trends in sigma_v and r_{half} set by their MW counterparts. We also investigate the nature of this scatter about the MW dSph mass profiles for the "Classical" (i.e. M_V<-8) MW and M31 dSphs. When comparing both the "classical" MW and M31 dSphs to the best--fit mass profiles in the size--velocity dispersion plane, we find general scatter in both the positive (i.e. hotter) and negative (i.e. colder) directions from these profiles. However, barring one exception (CVnI) only the M31 dSphs are found to scatter towards a colder regime, and, excepting the And I dSph, only MW objects scatter to hotter dispersions. We also note that the scatter for the combined population is greater than expected from measurement errors alone. We assess this divide in the context of the differing disc-to-halo mass (i.e. stars and baryons to total virial mass) ratios of the two hosts and argue that the underlying mass profiles for dSphs differ from galaxy to galaxy, and are modified by the baryonic component of the host.

(abridged) We apply the Jeans equation to estimate masses for eight of the brightest dSph galaxies. For Fornax we obtain a model-independent constraint on the maximum-circular velocity, Vmax=18_{-3}^{+5} km/s. Although we obtain only lower-limits of Vmax > 10 km/s for the remaining dSphs, we find that in all cases the enclosed mass at the projected half-light radius is well constrained and robust to a wide range of halo models and velocity anisotropies. We derive a simple analytic formula that estimates M(rhalf) accurately with respect to results from the full Jeans analysis. Applying this formula to the entire population of Local Group dSphs, we demonstrate a correlation such that M(rhalf)\propto rhalf^{1.4\pm 0.4}, or in terms of the mean density interior to the half-light radius, <rho> \propto rhalf^{-1.6\pm 0.4}. This relation is driven by the fact that the dSph data exhibit a correlation between global velocity dispersion and half-light radius. We argue that tidal forces are unlikely to have introduced this relation, but tides may have increased the scatter and/or altered the slope. While the data are well described by mass profiles ranging over a factor of < 2 in normalization (Vmax ~ 10-20 km/s), we consider the hypothesis that all dSphs are embedded within a "universal" dark matter halo. We show that in addition to the power law M\propto r^{1.4}, viable candidates include a cuspy "NFW" halo with Vmax ~ 15 km/s and scale radius r_0 ~ 800 pc, as well as a cored halo with Vmax ~ 13 km/s and r_0 ~ 150 pc. Finally, assuming that their measured velocity dispersions accurately reflect their masses, the smallest dSphs now allow us to resolve dSph densities at radii as small as a few tens of pc. At these small scales we find mean densities as large as <\rho> ~ 5 Msun/pc^3 (200 GeV/cm^3).

We present a 3D kinematical analysis of stars located in Kapteyn Selected Area 71 (l = 167.1, b = -34.7), where previously a stellar excess was found (Dinescu et al. 2002, Newberg et al. 2002). Previous findings indicated that the stellar excess has a cold kinematical signature as inferred from proper motions, and was initially associated with debris from the Sagittarius dwarf galaxy (Sgr) -- namely the southern trailing tail. We have obtained radial velocities using the Hydra multiobject spectrograph on the WIYN 3.5 m telescope. Results for 183 proper-motion selected stars indicate that the dominant population in this stellar excess is not debris from Sgr, but rather a population that kinematically belongs to the ring-like stream that is now known as the Monoceros stream. The orbit determined for this population agrees very well with the predictions for the Monoceros stream from Penarrubia et al. (2005). The radial-velocity dispersion of this population is between 20 and 30 km/s, lower than that of the Galactic field. Also, the shape of the radial-velocity distribution shows a sharp cut-off on one side, which is more in line with model predictions of the disruption of a satellite rather than with the distribution of the Galactic field. Despite the fact that we now believe most of the stars in the stellar excess to be part of Monoceros, about ten stars in this stellar excess have highly negative radial velocities, which is a clear indication of their membership to the Sgr trailing tidal tail.

We present the first determination of the proper motion of the neutron star low mass X-ray binary {Cen X-4} measured from relative astrometry of the secondary star using optical images at different epochs. We determine the Galactic space velocity components of the system and find them to be significantly different from the mean values that characterize the kinematics of stars belonging to the halo, and the thin and the thick disc of the Galaxy. The high metallicity of the secondary star of the system rules out a halo origin and indicates that the system probably originated in the Galactic disc. A statistical analysis of the galactocentric motion revealed that this binary moves in a highly eccentric ($e\simeq 0.85\pm0.1$) orbit with an inclination of $\simeq 110^\circ$ to the Galactic plane. The large Galactic space velocity components strongly support that a high natal kick as a result of a supernova explosion could have propelled the system into such an orbit from a birth place in the Galactic disc. The high Li abundance in the secondary, comparable to that of stars in star forming regions and young stellar clusters like the Pleiades, may suggest a relatively recent formation of the system. Following the orbit backwards in time, we found that the system could have been in the inner regions of the Galactic disc $\sim$100--200 Myr ago. The neutron star might have formed at that moment. However, we cannot rule out the possibility that the system formed at a much earlier time if a Li production mechanism exists in this LMXB.