BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jorge Peñarrubia

Jorge Peñarrubia appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COgr-qcastro-ph.SRhep-th

Catalog footprint

What is connected

31works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2026arXiv

HSTPROMO Internal Proper Motion Kinematics of Dwarf Spheroidal Galaxies: II. Velocity Anisotropy and Dark Matter Cusp Slope of Sculptor

We analyze three epochs of HST imaging over 20 years for the Sculptor dwarf spheroidal galaxy, measuring precise proper motions for 119 stars and combining them with 1760 existing line-of-sight velocities. This catalog yields the first radially-resolved 3D velocity dispersion profiles for Sculptor. We confirm mild oblate rotation, with major-axis velocities reaching $\sim 2$ km s$^{-1}$ beyond 20.0 arcmin. Using a methodology similar to that in the first paper in this series, we solve the Jeans equations in oblate axisymmetric geometry to infer the galaxy's mass profile. Our modeling reveals a significant degeneracy due to the unknown galaxy inclination, which is overlooked under spherical symmetry assumptions. This degeneracy allows acceptable fits across a range of dark matter profiles, from cuspy to cored. While we do not directly constrain the inclination with our Jeans models, higher-order line-of-sight velocity moments provide useful additional constraints: comparisons with scalefree models from de Bruijne et al. (1996) favor highly flattened (more face-on) configurations. Adopting an inclination well consistent with these comparisons ($i = 57.1$ degrees), we find, alongside radial velocity anisotropy, a dark matter density slope of $Γ_{\rm dark} = 0.29^{+0.31}_{-0.41}$ within the radial extent of the 3D velocity data, ruling out a cusp with $Γ_{\rm dark} \leq -1$ at 99.8% confidence. This confidence increases for lower inclinations and decreases drastically for nearly edge-on configurations. The results qualitatively agree with $Λ$CDM, SIDM, and Fuzzy DM scenarios that predict core formation, while our specific measurements provide quantitative constraints on the prescriptions of feedback, cross sections, or particle masses required by these models, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Signatures of dark subhalos in dwarf spheroidal galaxies: I. Fluctuations in surface density

Dark matter (DM) subhalos offer critical tests of cosmological models through their abundance and properties, yet most remain undetectable due to their lack of stars. We investigate whether their presence leaves measurable imprints on the projected stellar density fields of dwarf spheroidal galaxies (dSphs). Building on literature $N$-body experiments, we show that subhalo interactions induce subtle out-of-equilibrium fluctuations appearing as density corrugations. In a CDM framework, these fluctuations are dominated by the most massive subhalos in the host halo. We develop a Fourier-based framework to quantify these features, identifying characteristic peaks in the spatial frequency spectrum that are well described by Voigt profiles. The peak parameters are sensitive to both the subhalo mass function and the number of stellar tracers. For the configurations tested, $N_{\star} \sim 10^5$ stars suffice to detect subhalo populations with $M_{\rm subhalo} \lesssim 10^6~\mathrm{M}_{\odot}$, while larger masses produce stronger and more complex signatures. We assess the feasibility of this technique by analyzing Gaia and HST data: in this context, the Fornax dwarf shows residual low-frequency structures resembling those in our controlled subhalo experiments, making it an interesting case for follow-up. Prospectively, wide-field surveys such as Euclid, the Nancy Grace Roman Space Telescope, and the Vera C. Rubin Observatory are expected to deliver stellar samples of $N_{\star} \sim 10^5$ per dwarf, offering compelling prospects for probing subhalo imprints. Our results introduce a novel pathway to constrain the subhalo mass function in dSphs, and motivate follow-up work that incorporates alternative DM models and additional dynamical perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effects on the local dark matter distribution due to the Large Magellanic Cloud

We study the local dark matter distribution in two models for the Milky Way (MW)-Large Magellanic Cloud (LMC) interaction. The effect of the LMC on the local dark matter distribution is dependent on the evolution of the MW-LMC system, such that a static model is insufficient to accurately model the dark matter velocity distribution in the solar neighbourhood. An evolved model boosts local LMC dark matter particle velocities by nearly 50%, to a median value of $\approx750$km/s. MW dark matter particles also experience a velocity boost, which we identify as being caused by reflex motion owing to the infall of the LMC. We study the implications of LMC particles in the solar neighbourhood for dark matter detection experiments. Specifically, the directionality of LMC particles is distinguishable from the MW particles, with a difference in the apparent origin centroid location between the MW and LMC particles of $26\pm6 ^\circ$. This unique identifier, along with their high velocities, can be utilised by directional detectors to search for dark matter particles originating in the LMC.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Forward and Back: Kinematics of the Palomar 5 Tidal Tails

The tidal tails of Palomar 5 (Pal 5) have been the focus of many spectroscopic studies in an attempt to identify individual stars lying along the stream and characterise their kinematics. The well-studied trailing tail has been explored out to a distance of 15^\text{o} from the cluster centre, while less than four degrees have been examined along the leading tail. In this paper, we present results of a spectroscopic study of two fields along the leading tail that we have observed with the AAOmega spectrograph on the Anglo-Australian telescope. One of these fields lies roughly 7^\text{o} along the leading tail, beyond what has been previously been explored spectroscopically. Combining our measurements of kinematics and line strengths with Pan-STARRS1 photometric data and Gaia EDR3 astrometry, we adopt a probabilistic approach to identify 16 stars with high probability of belonging to the Pal 5 stream. Eight of these stars lie in the outermost field and their sky positions confirm the presence of ``fanning'' in the leading arm. We also revisit previously-published radial velocity studies and incorporate Gaia EDR3 astrometry to remove interloping field stars. With a final sample of 109 {\it bona fide} Pal 5 cluster and tidal stream stars, we characterise the 3D kinematics along the the full extent of the system. We provide this catalogue for future modeling work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Implications of the Milky Way travel velocity for dynamical mass estimates of the Local Group

The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history including the recent pericentric passage of the LMC, and recent work has found that the MW disk is moving with a lower bound "travel velocity" of $\sim 32~{\rm km}~{\rm s}^{-1}$ with respect to the outer stellar halo. Previous Timing Argument measurements attempt to account for this non-equilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10--12\% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Structure and kinematics of tidally limited satellite galaxies in LCDM

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tidally stripped halo stars from the Large Magellanic Cloud in the Galactic North

We examine whether the Large Magellanic Cloud (LMC) is currently losing its stellar halo to Milky Way (MW) tides. We present a live $N$-body model for the ongoing MW-LMC interaction that predicts a prominent stream of stars tidally stripped from the progenitor LMC. We use this model to define a strategy to search for stripped material in kinematic space. Of the available stellar tracers, we conclude that samples of RR Lyrae stars provide the highest density of kinematic tracers at present. Using a sample of RR Lyrae stars with Gaia EDR3 astrometry we show that the LMC stellar halo in the Southern Galactic hemisphere extends at least out to $\sim 30^\circ$ from the galaxy centre. In addition, several leading arm candidates are found in the Northern hemisphere as far above the disc plane as $b=+34^\circ$ (at 68$^\circ$ from the LMC).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two-Point Separation Functions for Modeling Wide Binary Systems in Nearby Dwarf Galaxies

We use a geometric method to derive (two-dimensional) separation functions amongst pairs of objects within populations of specified position function $\mathrm{d} N/\mathrm{d} \vec{R}$. We present analytic solutions for separation functions corresponding to a uniform surface density within a circular field, a Plummer sphere (viewed in projection), and the mixture thereof -- including contributions from binary objects within both sub-populations. These results enable inferences about binary object populations via direct modeling of object position and pair separation data, without resorting to standard estimators of the two-point correlation function. Analyzing mock data sets designed to mimic known dwarf spheroidal galaxies, we demonstrate the ability to recover input properties including the number of wide binary star systems and, in cases where the number of resolved binary pairs is assumed to be $\gtrsim$ a few hundred, characteristic features (e.g., steepening and/or truncation) of their separation function. Combined with forthcoming observational capabilities, this methodology opens a window onto the formation and/or survival of wide binary populations in dwarf galaxies, and offers a novel probe of dark matter substructure on the smallest galactic scales.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Creation/destruction of ultra-wide binaries in tidal streams

This paper uses statistical and $N$-body methods to explore a new mechanism to form binary stars with extremely large separations ($> 0.1\,{\rm pc}$), whose origin is poorly understood. Here, ultra-wide binaries arise via chance entrapment of unrelated stars in tidal streams of disrupting clusters. It is shown that (i) the formation of ultra-wide binaries is not limited to the lifetime of a cluster, but continues after the progenitor is fully disrupted, (ii) the formation rate is proportional to the local phase-space density of the tidal tails, (iii) the semimajor axis distribution scales as $p(a)d a\sim a^{1/2}d a$ at $a\ll D$, where $D$ is the mean interstellar distance, and (vi) the eccentricity distribution is close to thermal, $p(e)d e= 2 e d e$. Owing to their low binding energies, ultra-wide binaries can be disrupted by both the smooth tidal field and passing substructures. The time-scale on which tidal fluctuations dominate over the mean field is inversely proportional to the local density of clumps. Monte-Carlo experiments show that binaries subject to tidal evaporation follow $p(a)d a\sim a^{-1}d a$ at $a\gtrsim a_{\rm peak}$, known as Öpik's law, with a peak semi-major axis that contracts with time as $a_{\rm peak}\sim t^{-3/4}$. In contrast, a smooth Galactic potential introduces a sharp truncation at the tidal radius, $p(a)\sim 0$ at $a\gtrsim r_t$. The scaling relations of young clusters suggest that most ultra-wide binaries arise from the disruption of low-mass systems. Streams of globular clusters may be the birthplace of hundreds of ultra-wide binaries, making them ideal laboratories to probe clumpiness in the Galactic halo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reflex motion in the Milky Way stellar halo resulting from the Large Magellanic Cloud infall

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finding Horndeski theories with Einstein gravity limits

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of a disc on the population of cuspy and cored dark matter substructures in Milky Way-like galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A probability theory for non-equilibrium gravitational systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A timing constraint on the (total) mass of the Large Magellanic Cloud

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ñarrubia 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining the distribution of dark matter in dwarf spheroidal galaxies with stellar tidal streams

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The SLUGGS Survey: Multi-population dynamical modelling of the elliptical galaxy NGC 1407 from stars and globular clusters

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 $γ= 0.6\pm0.4$, which is consistent with either a DM cusp ($γ= 1$) or with a DM core $(γ= 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 $Λ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A dynamical model of the local cosmic expansion

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 $Ω_Λ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How chameleons core dwarfs with cusps

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Serendipitous Discovery of a Thin Stellar Stream near the Galactic Bulge in the Pan-STARRS1 3Pi Survey

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° 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°,+32°), 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Complex Structure of Stars in the Outer Galactic Disk as revealed by Pan-STARRS1

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The PAndAS Field of Streams: stellar structures in the Milky Way halo toward Andromeda and Triangulum

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 Σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Under the sword of Damocles: plausible regeneration of dark matter cusps at 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A kinematic study of the Andromeda dwarf spheroidal system

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical invariants and diffusion of merger substructures in time-dependent gravitational potentials

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Measuring the slopes of mass profiles for dwarf spheroidals in triaxial CDM potentials

We generate stellar distribution functions (DFs) in triaxial haloes in order to examine the reliability of slopes $Γ\equiv Δ{\rm log} M / Δ{\rm log} r$ inferred by applying mass estimators of the form $M\propto R_eσ^2$ (i.e. assuming spherical symmetry, where $R_e$ and $σ$ 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 $σ$ 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 $Γ$, irrespective of the degree of triaxiality or shape of the tracer number density profile.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The masses of Local Group dwarf spheroidal galaxies: The death of the universal mass profile

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A statistical method for measuring the Galactic potential and testing gravity with cold tidal streams

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The coupling between the core/cusp and missing satellite problems

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 Δ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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Method for Measuring (Slopes of) the Mass Profiles of Dwarf Spheroidal Galaxies

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 $Γ\equiv Δ\log M / Δ\log r=2.61_{-0.37}^{+0.43}$ and $Γ=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 $Γ$ (implying that the stated NFW exclusion levels are conservative).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Outer density profiles of 19 Galactic globular clusters from deep and wide-field imaging

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 $-γ$ 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 < γ< 4$ and a group of low-mass clusters ($ \le 10^5 $ M_sun), with steep profiles ($γ> 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Universal Mass Profile for Dwarf Spheroidal Galaxies?

(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 <ρ> ~ 5 Msun/pc^3 (200 GeV/cm^3).

0 citations0 reviewsNo code linkedOpen access