Stellar Stream Research Articles

Variation in the stellar mass function along stellar streams

ABSTRACT Stellar streams are the inevitable end product of star cluster evolution, with the properties of a given stream being related to its progenitor. We consider how the dynamical history of a progenitor cluster, as traced by the evolution of its stellar mass function, is reflected in the resultant stream. We generate model streams by evolving star clusters with a range of initial half-mass relaxation times and dissolution times via direct N-body simulations. Stellar streams that dissolve quickly show no variation in the stellar mass function along the stream. Variation is, however, observed along streams with progenitor clusters that dissolve after several relaxation times. The mass function at the edges of a stream is approximately primordial, as it is populated by the first stars to escape the cluster before segregation occurs. Moving inwards the mass function steepens as the intermediate parts of the stream consist of mostly low-mass stars that escaped the cluster after some segregation has occurred. The centre of the stream is then marked by a flatter mass function, as the region is dominated by high-mass stars that quickly segregated to the progenitor cluster’s centre and were the last stars to become unbound. We further find that the maximum slope of the mass function along the stream and the rate at which it decreases with distance from the dissolved progenitor serve as proxies for the dynamical state reached by the progenitor cluster before dissolution; this may be able to be applied to observed streams with near-future observations.

Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by S 5

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 × 1010 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 that the LMC’s effect depends both on the magnitude of the velocity kick imparted to the stream and the direction of this kick.

Via machinae: Searching for stellar streams using unsupervised machine learning

ABSTRACT We develop a new machine learning algorithm, via machinae, to identify cold stellar streams in data from the Gaia telescope. via machinae is based on ANODE, a general method that uses conditional density estimation and sideband interpolation to detect local overdensities in the data in a model agnostic way. By applying ANODE to the positions, proper motions, and photometry of stars observed by Gaia, via machinae obtains a collection of those stars deemed most likely to belong to a stellar stream. We further apply an automated line-finding method based on the Hough transform to search for line-like features in patches of the sky. In this paper, we describe the via machinae algorithm in detail and demonstrate our approach on the prominent stream GD-1. Though some parts of the algorithm are tuned to increase sensitivity to cold streams, the via machinae technique itself does not rely on astrophysical assumptions, such as the potential of the Milky Way or stellar isochrones. This flexibility suggests that it may have further applications in identifying other anomalous structures within the Gaia data set, for example debris flow and globular clusters.

Orbital phase-driven biases in galactic mass constraints from stellar streams

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. With 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).

Destruction of Galaxies as a Cause of the Appearance of Stellar Streams

Destruction of Galaxies as a Cause of the Appearance of Stellar Streams

Signature of a Massive Rotating Metal-poor Star Imprinted in the Phoenix Stellar Stream*

Abstract The Phoenix stellar stream has a low intrinsic dispersion in velocity and metallicity that implies the progenitor was probably a low-mass globular cluster. In this work we use Magellan/Magellan Inamori Kyocera Echelle (MIKE) high-dispersion spectroscopy of eight Phoenix stream red giants to confirm this scenario. In particular, we find negligible intrinsic scatter in metallicity ( σ ( [ Fe II / H ] ) = 0.04 − 0.03 + 0.11 ) and a large peak-to-peak range in [Na/Fe] and [Al/Fe] abundance ratios, consistent with the light element abundance patterns seen in the most metal-poor globular clusters. However, unlike any other globular cluster, we also find an intrinsic spread in [Sr ii/Fe] spanning ∼1 dex, while [Ba ii/Fe] shows nearly no intrinsic spread ( σ ( [ Ba II / H ] ) = 0.03 − 0.02 + 0.10 ). This abundance signature is best interpreted as slow-neutron-capture element production from a massive fast-rotating metal-poor star (15–20 M ⊙, v ini/v crit = 0.4, [Fe/H] = −3.8). The low inferred cluster mass suggests the system would have been unable to retain supernovae ejecta, implying that any massive fast-rotating metal-poor star that enriched the interstellar medium must have formed and evolved before the globular cluster formed. Neutron-capture element production from asymptotic giant branch stars or magneto-rotational instabilities in core-collapse supernovae provide poor fits to the observations. We also report one Phoenix stream star to be a lithium-rich giant (A(Li) = 3.1 ± 0.1). At [Fe/H ] = −2.93; it is among the most metal-poor lithium-rich giants known.

Kinematics of Antlia 2 and Crater 2 from the Southern Stellar Stream Spectroscopic Survey (S 5)

Abstract We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S 5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity–metallicity relation and is likely only now experiencing substantial stellar mass loss. Current dynamical models of Antlia 2 require it to have lost over 90% of its stars to tides, in tension with the low stellar mass loss implied by the updated metallicity. Overall, the new kinematic measurements support a tidal disruption scenario for the origin of these large and extended dwarf spheroidal galaxies.

Dancing in the void: hydrodynamical N-body simulations of the extremely metal-poor galaxy DDO 68

ABSTRACT Using hydrodynamical N-body simulations, we show that the observed structure and kinematics of the extremely metal-poor dwarf irregular galaxy DDO 68 is compatible with a merger event with at least two smaller satellite galaxies. We were able to obtain a self-consistent model that simultaneously reproduces several of its observed features, including the very asymmetric and disturbed shape of the stellar component, the overall ${\rm H\, \small {\rm I}}$ distribution and its velocity field, the arc-like stellar structure to the west, and the low surface brightness stellar stream to the north. The model implies the interaction of the main progenitor of DDO 68 with two systems with dynamical masses $7\times 10^8\, \mathrm{ M}_{\odot }$ and almost $10^8\, \mathrm{ M}_{\odot }$ – 1/20 and 1/150 times the dynamical mass of DDO 68, respectively. We show that the merger between DDO 68 and the most massive of its satellites offers a route to explain the large offset of DDO 68 from the mass–metallicity relation. Assuming that the interacting galaxies have metallicities prior to the merger compatible with those of galaxies with similar stellar masses, we provide quantitative evidence that gas mixing alone does not suffice at diluting the gas of the two components; according to our simulations, the ${\rm H\, \small {\rm II}}$ regions observed along the cometary tail trace the low metallicity of the accreted satellite rather than that of DDO 68’s main body. In this case, the mass corresponding to the low metallicity is that of the secondary body and DDO 68 becomes consistent with the mass–metallicity relation.

The Galaxy Progenitors of Stellar Streams around Milky Way–mass Galaxies in the FIRE Cosmological Simulations

Abstract Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)–mass galaxies from the FIRE-2 suite at z = 0, including seven isolated MW-mass systems and six hosts resembling the MW–M31 pair. In total, we identify 106 simulated stellar streams, with no significant differences in the number of streams and the masses of their progenitors between the isolated and paired environments. We resolve simulated streams with stellar masses ranging from ∼5 × 105 up to ∼109 M ⊙, similar to the mass range between the Orphan and Sagittarius streams in the MW. We confirm that present-day simulated satellite galaxies are good proxies for stellar stream progenitors, with similar properties including their stellar mass function, velocity dispersion, [Fe/H] and [α/H] evolution tracks, and orbital distribution with respect to the galactic disk plane. Each progenitor’s lifetime is marked by several important timescales: its infall, star formation quenching, and stream formation times. We show that the ordering of these timescales is different between progenitors with stellar masses higher and lower than ∼2 × 106 M ⊙. Finally, we show that the main factor controlling the rate of phase mixing, and therefore the rate of fading, of tidal streams from satellite galaxies in MW-mass hosts is the nonadiabatic evolution of the host potential. Other factors commonly used to predict phase-mixing timescales, such as progenitor mass and orbital circularity, show virtually no correlation with the number of dynamical times required for a stream to become phase-mixed.

Evidence of a Dwarf Galaxy Stream Populating the Inner Milky Way Halo

Abstract Stellar streams produced from dwarf galaxies provide direct evidence of the hierarchical formation of the Milky Way. Here, we present the first comprehensive study of the LMS-1 stellar stream, that we detect by searching for wide streams in the Gaia EDR3 data set using the STREAMFINDER algorithm. This stream was recently discovered by Yuan et al. We detect LMS-1 as a 60° long stream to the north of the galactic bulge, at a distance of ∼20 kpc from the Sun, together with additional components that suggest that the overall stream is completely wrapped around the inner Galaxy. Using spectroscopic measurements from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the Sloan Digital Sky Survey, and the Apache Point Observatory Galactic Evolution Experiment, we infer that the stream is very metal-poor (〈[Fe/H]〉 = −2.1) with a significant metallicity dispersion (σ [Fe/H] = 0.4), and it possesses a large radial velocity dispersion (σ v = 20 ± 4 km s−1). These estimates together imply that LMS-1 is a dwarf galaxy stream. The orbit of LMS-1 is close to polar, with an inclination of 75° to the galactic plane. Both the orbit and metallicity of LMS-1 are remarkably similar to the globular clusters NGC 5053, NGC 5024, and the stellar stream Indus. These findings make LMS-1 an important contributor to the stellar population of the inner Milky Way halo.

The Effects of Dark Matter and Baryonic Physics on the Milky Way Subhalo Population in the Presence of the Large Magellanic Cloud

Abstract Given recent developments in our understanding of the Large Magellanic Cloud's (LMC) impact on the Milky Way's (MW) dark matter subhalo population, we compare the signatures of dark matter and baryonic physics on subhalos in MW systems with realistic LMC analogs. In particular, we study the effects of self-interacting dark matter (SIDM), warm dark matter (WDM), and the Galactic disk on the peak maximum circular velocity (V peak) function, radial distribution, and spatial distribution of MW and LMC-associated subhalos using cosmological dark matter-only zoom-in simulations of MW+LMC systems. For a fixed abundance of subhalos expected to host dwarf galaxies (V peak ≳ 20 km s−1), SIDM and WDM can produce a similar mass-dependent suppression of the subhalo V peak function, while disk disruption is mass independent. Subhalos in the inner regions of the MW are preferentially disrupted by both self-interactions and the disk, while suppression in WDM is radially independent. The relative abundance of LMC-associated subhalos is not strongly affected by disk disruption or WDM, but is significantly suppressed in SIDM due to self-interactions with the LMC at early times and with the MW during LMC infall at late times, erasing spatial anisotropy in the MW subhalo population. These results provide avenues to distinguish dark matter and baryonic physics by combining properties of the MW and LMC subhalo populations probed by upcoming observations of satellite galaxies and stellar streams.

Novel constraints on the particle nature of dark matter from stellar streams

Tidal streams are highly sensitive to perturbations from passing dark matter (DM) subhalos and thus provide a means of measuring their abundance. In a recent paper, we analyzed the distribution of stars along the GD-1 stream with a combination of data from the Gaia satellite and the Pan-STARRS survey, and we demonstrated that the population of DM subhalos predicted by the cold dark matter (CDM) paradigm are necessary and sufficient to explain the perturbations observed in the linear density of stars. In this paper, we use the measurements of the subhalo mass function (SHMF) from the GD-1 data combined with a similar analysis of the Pal 5 stream to provide novel constraints on alternative DM scenarios that predict a suppression of the SHMF on scales smaller than the mass of dwarf galaxies, marginalizing over uncertainties in the slope and normalization of the unsuppressed SHMF and the susceptibility of DM subhalos in the inner Milky Way to tidal disruption. In particular, we derive a 95% lower limit on the mass of warm dark matter (WDM) thermal relics m WDM > 3.6 keV from streams alone that strengthens to m WDM > 6.2 keV when adding dwarf satellite counts. Similarly, we constrain the axion mass in ultra-light (“fuzzy”) dark matter (FDM) models to be m FDM > 1.4 × 10-21 eV from streams alone or m FDM > 2.2 × 10-21 eV when adding dwarf satellite counts. Because we make use of simple approximate forms of the streams' SHMF measurement, our analysis is easy to replicate with other alternative DM models that lead to a suppression of the SHMF.

Linking nearby stellar streams to more distant halo overdensities

Context. It has recently been shown that the halo near the Sun contains several kinematic substructures associated with past accretion events. For the more distant halo, there is evidence of large-scale density variations – in the form of stellar clouds or overdensities. Aims. We study the link between the local halo kinematic groups and three of these stellar clouds: the Hercules-Aquila cloud, the Virgo Overdensity, and the Eridanus-Phoenix overdensity. Methods. We perform orbital integrations in a standard Milky Way potential of a local halo sample extracted from Gaia EDR3 with the goal of predicting the location of the merger debris elsewhere in the Galaxy. We specifically focus on the regions occupied by the three stellar clouds and compare their kinematic and distance distributions with those predicted from the orbits of the nearby debris. Results. We find that the local halo substructures have families of orbits that tend to pile up in the regions where the stellar clouds have been found. The distances and velocities of the cloud’s member stars are in good agreement with those predicted from the orbit integrations, particularly for Gaia-Enceladus stars. This is the dominant contributor of all three overdensities, with a minor part stemming from the Helmi streams and to an even smaller extent from Sequoia. The orbital integrations predict no asymmetries in the sky distribution of halo stars, and they pinpoint where additional debris associated with the local halo substructures may be located.

Evidence for Ultra-diffuse Galaxy Formation through Tidal Heating of Normal Dwarfs

Abstract We have followed up two ultra-diffuse galaxies (UDGs), detected adjacent to stellar streams, with Hubble Space Telescope (HST) imaging and H i mapping with the Jansky Very Large Array in order to investigate the possibility that they might have a tidal origin. With the HST F814W and F555W images we measure the globular cluster (GC) counts for NGC 2708-Dw1 and NGC 5631-Dw1 as 2 − 1 + 1 and 5 − 2 + 1 , respectively. NGC 2708-Dw1 is undetected in H i down to a 3σ limit of log (M H i /M ⊙) = 7.3, and there is no apparent H i associated with the nearby stellar stream. There is a 2σ H i feature coincident with NGC 5631-Dw1. However, this emission is blended with a large gaseous tail emanating from NGC 5631 and is not necessarily associated with the UDG. The presence of any GCs and the lack of clear H i connections between the UDGs and their parent galaxies strongly disfavor a tidal dwarf galaxy origin, but cannot entirely rule it out. The GC counts are consistent with those of normal dwarf galaxies, and the most probable formation mechanism is one where these UDGs were born as normal dwarfs and were later tidally stripped and heated. We also identify an overluminous (M V = −11.1) GC candidate in NGC 2708-Dw1, which may be a nuclear star cluster transitioning to an ultra-compact dwarf as the surrounding dwarf galaxy gets stripped of stars.

Flattening of Dark Matter Cusps during Mergers: Model of M31

Abstract We run high-resolution fully GPU N-body simulations to demonstrate that the dark matter (DM) distribution in M31 is well fitted by a core-like profile. Rich observational data especially on the giant stellar stream provides stringent constraints on the initial conditions of our simulations. We demonstrate that the accretion of a satellite on a highly eccentric orbit heats up the central parts of M31, causes an outward migration of DM particles, flattens the central cusp over more than a decade in scale, and generates a new model-independent DM profile that is well fitted by a core. Our results imply that core-like central profiles could be a common feature of massive galaxies that have been initially cuspy but have accreted satellites on nearly radial orbits.

Regular and chaotic orbits in axisymmetric stellar systems

ABSTRACT The gravitational potentials of realistic galaxy models are in general non-integrable, in the sense that they admit orbits that do not have three independent isolating integrals of motion and are therefore chaotic. However, if chaotic orbits are a small minority in a stellar system, it is expected that they have negligible impact on the main dynamical properties of the system. In this paper, we address the question of quantifying the importance of chaotic orbits in a stellar system, focusing, for simplicity, on axisymmetric systems. Chaotic orbits have been found in essentially all (non-Stäckel) axisymmetric gravitational potentials in which they have been looked for. Based on the analysis of the surfaces of section, we add new examples to those in the literature, finding chaotic orbits, as well as resonantly trapped orbits among regular orbits, in Miyamoto-Nagai, flattened logarithmic and shifted Plummer axisymmetric potentials. We define the fractional contributions in mass of chaotic (ξc) and resonantly trapped (ξt) orbits to a stellar system of given distribution function (DF), which are very useful quantities, for instance in the study of the dispersal of stellar streams of galaxy satellites. As a case study, we measure ξc and ξt in two axisymmetric stellar systems obtained by populating flattened logarithmic potentials with the Evans ergodic DF, finding ξc ∼ 10−4 − 10−3 and ξt ∼ 10−2 − 10−1.

Cosmology of Sub-MeV Dark Matter Freeze-In.

Dark matter (DM) could be a relic of freeze-in through a light mediator, where the DM is produced by extremely feeble, IR-dominated processes in the thermal standard model plasma. In the simplest viable models with DM lighter than 1MeV, the DM has a small effective electric charge and is born with a nonthermal phase-space distribution. This DM candidate would cause observable departures from standard cosmological evolution. In this work, we combine data from the cosmic microwave background (CMB), Lyman-α forest, quasar lensing, stellar streams, and MilkyWay satellite abundances to set limits on freeze-in DM masses up to ∼20 keV, with the exact constraint depending on whether the DM thermalizes in its own sector. We perform forecasts for the CMB-S4 experiment, the Hydrogen Epoch of Reionization Array, and the Vera Rubin Observatory, finding that freeze-in DM masses up to ∼80 keV can be explored.

Phenomenological distribution of dark matter halos of galaxy comprising of stellar streams

Using the potential associated with NFW mass profile (Navarro et al. in Mon. Not. R. Astron. Soc. 402:21, 2010) enclosed at scale radius $r_{s}$ , we calculate the partition function of the clustering of dark halos comprising of stellar streams, which are the outcome of accretion of globular clusters or dwarf galaxies onto a massive host galaxy. The thermodynamics of the system of halos of stellar streams may provide crucial information on the properties of galaxy clustering. The distribution function is studied through statistical formalism. The study can give us an idea of the extended structure of galaxies and the contribution of halos to the thermodynamics and the mass content of the host galaxy. In an attempt to fit the theoretical orbital velocity to the galaxy rotation curve for the data available for galaxy $M33$ , we developed another potential that exactly leads to the best fit orbital velocity of the galaxy. Thus a phenomenological theory for the galaxy with a dark halo is established, with a new density profile.

Uncovering fossils of the distant Milky Way with UNIONS: NGC 5466 and its stellar stream

ABSTRACT We examine the spatial clustering of blue horizontal branch (BHB) stars from the u-band of the Canada–France Imaging Survey (CFIS, a component of the Ultraviolet Near-Infrared Optical Northern Survey, or UNIONS). All major groupings of stars are associated with previously known satellites, and among these is NGC 5466, a distant (16 kpc) globular cluster. NGC 5466 reportedly possesses a long stellar stream, although no individual members of the stream have previously been identified. Using both BHBs and more numerous red giant branch stars cross-matched to Gaia Data Release 2, we identify extended tidal tails from NGC 5466 that are both spatially and kinematically coherent. Interestingly, we find that this stream does not follow the same path as the previous detection at large distances from the cluster. We trace the stream across 31° of sky and show that it exhibits a very strong distance gradient in the range 10 < Rhelio < 30 kpc. We compare our observations to simple dynamical models of the stream and find that they are able to broadly reproduce the overall path and kinematics. The fact that NGC 5466 is so distant, traces a wide range of Galactic distances, has an identified progenitor, and appears to have recently had an interaction with the Galaxy’s disc makes it a unique test-case for dynamical modelling of the Milky Way.

Towards constraining warm dark matter with stellar streams through neural simulation-based inference

ABSTRACT A statistical analysis of the observed perturbations in the density of stellar streams can in principle set stringent constraints on the mass function of dark matter subhaloes, which in turn can be used to constrain the mass of the dark matter particle. However, the likelihood of a stellar density with respect to the stream and subhaloes parameters involves solving an intractable inverse problem which rests on the integration of all possible forward realizations implicitly defined by the simulation model. In order to infer the subhalo abundance, previous analyses have relied on Approximate Bayesian Computation (ABC) together with domain-motivated but handcrafted summary statistics. Here, we introduce a likelihood-free Bayesian inference pipeline based on Amortised Approximate Likelihood Ratios (AALR), which automatically learns a mapping between the data and the simulator parameters and obviates the need to handcraft a possibly insufficient summary statistic. We apply the method to the simplified case where stellar streams are only perturbed by dark matter subhaloes, thus neglecting baryonic substructures, and describe several diagnostics that demonstrate the effectiveness of the new method and the statistical quality of the learned estimator.