skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Substructure at High Speed. I. Inferring the Escape Velocity in the Presence of Kinematic Substructure
Abstract The local escape velocity provides valuable inputs to the mass profile of the galaxy, and requires understanding the tail of the stellar speed distribution. Following Leonard & Tremaine, various works have since modeled the tail of the stellar speed distribution as ( v esc v ) k , wherevescis the escape velocity, andkis the slope of the distribution. In such studies, however, these two parameters were found to be largely degenerate and often a narrow prior is imposed onkin order to constrainvesc. Furthermore, the validity of the power-law form can breakdown in the presence of multiple kinematic substructures or other mis-modeled features in the data. In this paper, we introduce a strategy that for the first time takes into account the presence of kinematic substructure. We model the tail of the velocity distribution as a sum of multiple power laws as a way of introducing a more flexible fitting framework. Using mock data and data from FIRE simulations of Milky Way-like galaxies, we show the robustness of this method in the presence of kinematic structure that is similar to the recently discovered Gaia Sausage. In a companion paper, we present the new measurement of the escape velocity and subsequently the mass of the Milky Way using Gaia eDR3 data.  more » « less
Award ID(s):
2020275
PAR ID:
10363092
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
926
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 188
Size(s):
Article No. 188
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be v esc = 445 8 + 25 km s−1. A fit with a single component to the same data gives v esc = 472 12 + 17 km s−1. Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of c 200 = 19 7 + 11 and a mass of M 200 = 4.6 0.8 + 1.5 × 10 11 M , which is considerably lighter than previous measurements. 
    more » « less
  2. Abstract We present the discovery of Aquarius III, an ultra-faint Milky Way satellite galaxy identified in the second data release of the DECam Local Volume Exploration survey. Based on deeper follow-up imaging with DECam, we find that Aquarius III is a low-luminosity ( M V = 2.5 0.5 + 0.3 ; L V = 850 260 + 380 L ), extended ( r 1 / 2 = 41 8 + 9 pc) stellar system located in the outer halo (D= 85 ± 4 kpc). From medium-resolution Keck/DEIMOS spectroscopy, we identify 11 member stars and measure a mean heliocentric radial velocity of v sys = 13.1 0.9 + 1.0 km s 1 for the system and place an upper limit ofσv< 3.5 km s−1v< 1.6 km s−1) on its velocity dispersion at the 95% (68%) credible level. Based on calcium-triplet metallicities of the six brightest red giant members, we find that Aquarius III is very metal-poor ([Fe/H]= − 2.61 ± 0.21) with a statistically significant metallicity spread ( σ [ Fe / H ] = 0.46 0.14 + 0.26 dex). We interpret this metallicity spread as strong evidence that the system is a dwarf galaxy as opposed to a star cluster. Combining our velocity measurement with Gaia proper motions, we find that Aquarius III is currently situated near its orbital pericenter in the outer halo (rperi= 78 ± 7 kpc) and that it is plausibly on first infall onto the Milky Way. This orbital history likely precludes significant tidal disruption from the Galactic disk, notably unlike other satellites with comparably low velocity dispersion limits in the literature. Thus, if further velocity measurements confirm that its velocity dispersion is truly belowσv≲ 2 km s−1, Aquarius III may serve as a useful laboratory for probing galaxy formation physics in low-mass halos. 
    more » « less
  3. Abstract The 3D radial escape-velocity profile of galaxy clusters has been suggested to be a promising and competitive tool for constraining mass profiles and cosmological parameters in an accelerating universe. However, the observed line-of-sight escape profile is known to be suppressed compared to the underlying 3D radial (or tangential) escape profile. Past work has suggested that velocity anisotropy in the phase-space data is the root cause. Instead, we find that the observed suppression is from the statistical undersampling of the phase spaces and that the 3D radial escape edge can be accurately inferred from projected data. We build an analytical model for this suppression that only requires the number of observed galaxiesNin the phase-space data within the sky-projected range 0.3 ≤r/R200,critical≤ 1. The radially averaged suppression function is an inverse power law Z v = 1 + ( N 0 / N ) λ withN0= 17.818 andλ= 0.362. We test our model withN-body simulations, using dark matter particles, subhalos, and semianalytic galaxies as the phase-space tracers, and find excellent agreement. We also assess the model for systematic biases from cosmology (ΩΛ,H0), cluster mass (M200,critical), and velocity anisotropy (β). We find that varying these parameters over large ranges can impart a maximal additional fractional change in 〈Zv〉 of 2.7%. These systematics are highly subdominant (by at least a factor of 13.7) to the suppression fromN. 
    more » « less
  4. Abstract Ultra-high-energy cosmic rays (UHECRs), particles characterized by energies exceeding 1018eV, are generally believed to be accelerated electromagnetically in high-energy astrophysical sources. One promising mechanism of UHECR acceleration is magnetized turbulence. We demonstrate from first principles, using fully kinetic particle-in-cell simulations, that magnetically dominated turbulence accelerates particles on a short timescale, producing a power-law energy distribution with a rigidity-dependent, sharply defined cutoff well approximated by the form f cut E , E cut = sech ( E / E cut ) 2 . Particle escape from the turbulent accelerating region is energy dependent, withtesc∝E−δandδ∼ 1/3. The resulting particle flux from the accelerator follows dN / dEdt E s sech ( E / E cut ) 2 , withs∼ 2.1. We fit the Pierre Auger Observatory’s spectrum and composition measurements, taking into account particle interactions between acceleration and detection, and show that the turbulence-associated energy cutoff is well supported by the data, with the best-fitting spectral index being s = 2.1 0.13 + 0.06 . Our first-principles results indicate that particle acceleration by magnetically dominated turbulence may constitute the physical mechanism responsible for UHECR acceleration. 
    more » « less
  5. Abstract The genericity of Arnold diffusion in the analytic category is an open problem. In this paper, we study this problem in the followinga prioriunstable Hamiltonian system with a time-periodic perturbation H ε ( p , q , I , φ , t ) = h ( I ) + i = 1 n ± 1 2 p i 2 + V i ( q i ) + ε H 1 ( p , q , I , φ , t ) , where ( p , q ) R n × T n , ( I , φ ) R d × T d withn,d⩾ 1,Viare Morse potentials, andɛis a small non-zero parameter. The unperturbed Hamiltonian is not necessarily convex, and the induced inner dynamics does not need to satisfy a twist condition. Using geometric methods we prove that Arnold diffusion occurs for generic analytic perturbationsH1. Indeed, the set of admissibleH1isCωdense andC3open (a fortiori,Cωopen). Our perturbative technique for the genericity is valid in theCktopology for allk∈ [3, ∞) ∪ {∞,ω}. 
    more » « less