BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jorge Peñarrubia

Jorge Peñarrubia contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.CO

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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
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 &#34;universal&#34; dark matter halo. We show that in addition to the power law M\propto r^{1.4}, viable candidates include a cuspy &#34;NFW&#34; 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