skip to main content

Title: Galactic potential constraints from clustering in action space of combined stellar stream data
ABSTRACT Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighbouring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Palomar 5 (Pal 5), Orphan, and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a Stäckel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1, Pal 5, and Orphan streams whose combined data set yields $M(\lt 20\, \mathrm{kpc}) = 2.96^{+0.25}_{-0.26} \times 10^{11} \, \mathrm{ M}_{\odot}$. When including the Helmi stream in the data set, the mass uncertainty increases to $M(\lt 20\, \mathrm{kpc}) = 3.12^{+3.21}_{-0.46} \times 10^{11} \, \mathrm{M}_{\odot}$.
; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
4170 to 4193
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT One of the most promising tracers of the Galactic potential in the halo region is stellar streams. However, individual stream fits can be limited by systematic biases. To study these individual stream systematics, we fit streams in Milky Way-like galaxies from Feedback In Realistic Environments cosmological galaxy formation simulations with an analytical gravitational potential by maximizing the clustering of stream stars in action space. We show that for coherent streams the quality of the constraints depends on the orbital phase of the observed stream stars, despite the fact that the phase information is discarded in action-clustering methods. Streams on intermediate phases give the most accurate results, whereas pericentre streams can be highly biased. This behaviour is tied to the amount of correlation present between positions and momenta in each stream’s data: weak correlation in pericentre streams prohibits efficient differentiation between potentials, while strong correlation in intermediate streams promotes it. Although simultaneous fitting of multiple streams is generally prescribed as the remedy to combat individual stream biases, we find that combining multiple pericentric streams is not enough to yield a bias-free result. We finally show that adopting the two-component Stäckel model does not fundamentally induce a biased mass estimate. Withmore »our full data set of two multiwrap streams, we recovered the true rotation curve of the simulated galaxy within $12{{\ \rm per\ cent}}$ over the entire range of radii covered by our set of stars (10–176 kpc) and within $6.5{{\ \rm per\ cent}}$ between the 5 and 95 percentile distance range (23–109 kpc).« less

    We report a tentative detection of the circumgalactic medium (CGM) of Wolf–Lundmark–Melotte (WLM), an isolated, low-mass (logM*/M⊙ ≈ 7.6), dwarf irregular galaxy in the Local Group (LG). We analyse an HST/COS archival spectrum of a quasar sightline (PHL2525), which is 45 kpc (0.5 virial radius) from WLM and close to the Magellanic Stream (MS). Along this sightline, two ion absorbers are detected in Si ii, Si iii, Si iv, C ii, and C iv at velocities of ∼−220 km s−1 (Component v-220) and ∼−150 km s−1 (Component v-150). To identify their origins, we study the position–velocity alignment of the components with WLM and the nearby MS. Near the magellanic longitude of PHL2525, the MS-related neutral and ionized gas moves at ≲−190 km s−1, suggesting an MS origin for Component v-220, but not for Component v-150. Because PHL2525 passes near WLM and Component v-150 is close to WLM’s systemic velocity (∼−132 km s−1), it is likely that Component v-150 arises from the galaxy’s CGM. This results in a total Si mass in WLM’s CGM of $M_{\rm Si}^{\rm CGM}\sim (0.2-1.0)\times 10^5~\mathrm{M}_\odot$ using assumption from other COS dwarf studies. Comparing $M_{\rm Si}^{\rm CGM}$ to the total Si mass synthesized in WLM over its lifetime (∼1.3 × 105 M⊙), we find ∼3 per cent is locked in stars,more »∼6 per cent in the ISM, ∼15–77 per cent in the CGM, and the rest (∼14–76 per cent) is likely lost beyond the virial radius. Our finding resonates with other COS dwarf galaxy studies and theoretical predictions that low-mass galaxies can easily lose metals into their CGM due to stellar feedback and shallow gravitational potential.

    « less
  3. Abstract Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey ( S 5 ), Gaia EDR3, and the Dark Energy Survey, to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from ∼14–19 × 10 10 M ⊙ , probed over a range of closest approach times and distances. With the most constraining stream (Orphan–Chenab), we measure an LMC mass of 18.8 − 4.0 + 3.5 × 10 10 M ⊙ , probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC’s influence on structures in the Milky Way. Using this sample of streams, we find that the LMC’s effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on the impulse approximation and we show thatmore »the LMC’s effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick.« less
  4. Abstract

    We analyse stellar streams in action-angle coordinates combined with recent local direct acceleration measurements to provide joint constraints on the potential of our galaxy. Our stream analysis uses the Kullback–Leibler divergence with a likelihood analysis based on the two-point correlation function. We provide joint constraints from pulsar accelerations and stellar streams for local and global parameters that describe the potential of the Milky Way (MW). Our goal is to build an “acceleration ladder,” where direct acceleration measurements that are currently limited in dynamic range are combined with indirect techniques that can access a much larger volume of the MW. To constrain the MW potential with stellar streams, we consider the Palomar 5, Orphan, Nyx, Helmi, and GD1 streams. Of the potential models that we have considered here, the preferred potential for the streams is a two-component Staeckel potential. We also compare the vertical accelerations from stellar streams and pulsar timing, defining the functionf(z)=α1pulsarzΦz, where Φ is the MW potential determined from stellar streams andα1 pulsarzis the vertical acceleration determined from pulsar timing observations. Our analysis indicates that the Oort limit determined from streams is consistently (regardless of the choice of potential) lowermore »than that determined from pulsar timing observations. The calibration we have derived here may be used to correct the estimate of the acceleration from stellar streams.

    « less
  5. ABSTRACT The fraction of stars that are in binaries or triples at the time of stellar death and the fraction of these systems that survive the supernova explosion are crucial constraints for evolution models and predictions for gravitational wave source populations. These fractions are also subject to direct observational determination. Here, we search 10 supernova remnants containing compact objects with proper motions for unbound binaries or triples using Gaia EDR3 and new statistical methods and tests for false positives. We confirm the one known example of an unbound binary, HD 37424 in G180.0−01.7, and find no other examples. Combining this with our previous searches for bound and unbound binaries, and assuming no bias in favour of finding interacting binaries, we find that 72.0 per cent (52.2–86.4 per cent, 90 per cent confidence) of supernova producing neutron stars are not binaries at the time of explosion, 13.9 per cent (5.4–27.2 per cent) produce bound binaries, and 12.5 per cent (2.8–31.3 per cent) produce unbound binaries. With a strong bias in favour of finding interacting binaries, the medians shift to 76.0 per cent were not binaries at death, 9.5 per cent leave bound binaries, and 13.2 per cent leave unbound binaries. Of explosions that do not leave binaries, ${\lt}18.9{{\ \rm per\ cent}}$ can be fully unbound triples. These limits are conservatively for $M\gtmore »5\, \mathrm{M}_\odot$ companions, although the mass limits for some individual systems are significantly stronger. At birth, the progenitor of PSR J0538+2817 was probably a 13–$19\, \mathrm{M}_\odot$ star, and at the time of explosion, it was probably a Roche limited, partially stripped star transferring mass to HD 37424 and then producing a Type IIL or IIb supernova.« less