Stellar Stream Research Articles

Detection of Strong Epicyclic Density Spikes in the GD-1 Stellar Stream: An Absence of Evidence for the Influence of Dark Matter Subhalos?

Abstract The density variations in thin stellar streams may encode important information on the nature of dark matter. For instance, if dark matter aggregates into massive subhalos, these perturbers are expected to scatter stars out of dynamically cold stellar streams, possibly leading to detectable gaps in those structures. Here, we reexamine the density variations in the GD-1 stream, using Gaia Data Release 2 (DR2) astrometry and Panoramic Survey Telescope and Rapid Response System photometry, together with high-precision radial velocities measured with the Canada–France–Hawaii Telescope/ESPaDOnS and Very Large Telescope/Ultraviolet and Visual Echelle Spectrograph instruments and complemented with public radial velocity catalogs. We show that after correcting for projection effects, the density profile exhibits high contrast periodic peaks, separated by 2.64 ± 0.18 kpc. An N-body simulation is presented that reproduces this striking morphology with simple epicyclic motion in a smooth Galactic potential. We also discuss the reliability of measuring density variations using ground-based photometric surveys, and for the particular case of GD-1 we highlight some of the artifacts present in the Gaia DR2 catalog along its track. Massive dark subhalos do not appear to be required to explain the density clumping along GD-1.

Gravitational potential from small-scale clustering in action space: application to Gaia Data Release 2

ABSTRACT Most measurements of mass in astronomy that use kinematics of stars or gas rely on assumptions of equilibrium that are often hard to verify. Instead, we develop a novel idea that uses the clustering in action space, as a probe of the underlying gravitational potential: the correct potential should maximize small-scale clustering in the action space. We provide a first-principle derivation of likelihood using the two-point correlation function in action space, and we test it against simulations of stellar streams. We then apply this method to the second data release of Gaia, and we use it to measure the radial force fraction fh and logarithmic slope α of the dark matter halo profile. We investigate stars within 9–11 kpc and 11.5–15 kpc from the Galactic Centre, and we find (fh, α) = (0.391 ± 0.009, 1.835 ± 0.092) and (0.351 ± 0.012, 1.687 ± 0.079), respectively. We also confirm that the set of parameters that maximize the likelihood function does correspond to the most clustering in the action space. The best-fitting circular velocity curve for the Milky Way potential is consistent with past measurements (although it is ∼5–10 per cent lower than previous methods that use masers or globular clusters). Our work provides a clear demonstration of the full statistical power that lies in the full phase space information, relieving the need for ad hoc assumptions such as virial equilibrium, circular motion or steam-finding algorithms.

Lithium-rotation connection in the newly discovered young stellar stream Psc–Eri (Meingast 1)

Context. As a fragile element, lithium is a sensitive probe of physical processes occurring in stellar interiors. Aims. We aim to investigate the relationship between lithium abundance and rotation rate in low-mass members of the newly discovered 125 Myr-old Psc–Eri stellar stream. Methods. We obtained high-resolution optical spectra and measured the equivalent width of the 607.8 nm LiI line for 40 members of the Psc–Eri stream, whose rotational periods have been previously derived. Results. We show that a tight correlation exists between the lithium content and rotation rate among the late-G to early-K-type stars of the Psc–Eri stream. Fast rotators are systematically Li rich, while slow rotators are Li depleted. This trend mimics that previously reported for the similar age Pleiades cluster. Conclusions. The lithium-rotation connection thus seems to be universal over a restricted effective temperature range for low-mass stars at or close to the zero-age main sequence, and does not depend on environmental conditions.

The kinematical and space structures of IC 2391 open cluster and moving group with Gaia-DR2

The kinematical parameters, spatial shape and structure of the open cluster IC 2391 and the associated stellar stream are studied here using Gaia Data Release 2 (GDR2) astrometry data. The apex positions are determined for the open cluster IC 2391 (data taken from Cantat-Gaudin et al.) and for the kinematical stream’s stars mentioned in Montes et al. employing both convergent point and AD-diagram methods. The values of apex coordinates are: (A, D)CP = (; for cluster) and (, ; stream), and (A0, D0) = (, ; cluster) and (, ; stream). The results are in good agreement with the previously calculated values. The positions of the stars in the disk and the spatial dispersion velocities are determined. The paths of cluster and associated stream are traced in the disk by orbit calculation back in time to their places of formation. A possible genetic relationship between the cluster and stream has been detected. The approximation of the spatial and kinematical shape of the stream and the cluster is made. According to this study, even though currently the cluster and stream seem to have a spatial difference in their locations, they appear to have formed in the same region of the Galactic disk.

A Panoramic Landscape of the Sagittarius Stream in Gaia DR2 Revealed with the STREAMFINDER Spyglass

Abstract We present the first full six-dimensional panoramic portrait of the Sagittarius stream, obtained by searching for wide stellar streams in the Gaia DR2 data set with the STREAMFINDER algorithm. We use the kinematic behavior of the sample to devise a selection of Gaia RR Lyrae, providing excellent distance measurements along the stream. The proper motion data are complemented with radial velocities from public surveys. We find that the global morphological and kinematic properties of the Sagittarius stream are still reasonably well reproduced by the simple Law & Majewski model (LM10), although the model overestimates the leading arm and trailing arm distances by up to ∼15%. The sample newly reveals the leading arm of the Sagittarius stream as it passes into very crowded regions of the Galactic disk toward the Galactic anticenter direction. Fortuitously, this part of the stream is almost exactly at the diametrically opposite location from the Galactic center to the progenitor, which should allow an assessment of the influence of dynamical friction and self-gravity in a way that is nearly independent of the underlying Galactic potential model.

An extended Pal 5 stream in Gaia DR2

ABSTRACT We present the results of a detailed search for members of the Pal 5 tidal tail system in Gaia Data Release 2 (DR2). Tidal tails provide a sensitive method for measuring the current and past gravitational potential of their host galaxy as well as for testing predictions for the abundance of dark matter subhaloes. The Pal 5 globular cluster and its associated tails are an excellent candidate for such analysis; however, only ∼23° of arc are currently known, with in particular the leading tail much shorter than the trailing. Using Gaia DR2 and its precise astrometry, we extend the known extent of the Pal 5 tail to ∼30°, 7° of which are newly detected along the leading arm. The detected leading and trailing arms are symmetric in length and remain near constant width. This detection constrains proposed models in which the Galactic bar truncates Pal 5’s leading arm. Follow-up spectroscopic observations are necessary to verify the candidate stream stars are consistent with the known tidal tails. If confirmed, this Pal 5 stream extension opens up new possibilities to constrain the Galactic potential.

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

ABSTRACT 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 $\, {\rm km~s^{-1}}$ 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 modelling of stellar streams and/or smooth halo tracers such as blue horizontal branch or RR Lyrae stars.

The M31/M33 tidal interaction: a hydrodynamic simulation of the extended gas distribution

ABSTRACT We revisit the orbital history of the Triangulum galaxy (M33) around the Andromeda galaxy (M31) in view of the recent Gaia Data Release 2 proper motion measurements for both Local Group galaxies. Earlier studies consider highly idealized dynamical friction, but neglect the effects of dynamical mass loss. We show the latter process to be important using mutually consistent orbit integration and N-body simulations. Following this approach, we find an orbital solution that brings these galaxies to within ∼50 kpc of each other in the past, ∼6.5 Gyr ago. We explore the implications of their interaction using an N-body/hydrodynamical simulation with a focus on the origin of two prominent features: (1) M31’s Giant Stellar Stream; and (2) the M31–M33 H i filament. We find that the tidal interaction does not produce a structure reminiscent of the stellar stream that survives up to the present day. In contrast, the M31–M33 H i filament is likely a fossil structure dating back to the time of the ancient encounter between these galaxies. Similarly, the observed outer disc warp in M33 may well be a relic of this past event. Our model suggests the presence of a tidally induced gas envelope around these galaxies, and the existence of a diffuse gas stream, the ‘Triangulum stream’, stretching for tens of kpc from M33 away from M31. We anticipate upcoming observations with the recently commissioned, Five-hundred-metre Aperture Spherical radio Telescope that will target the putative stream in its first years of operation.

Elemental Abundances in M31: A Comparative Analysis of Alpha and Iron Element Abundances in the the Outer Disk, Giant Stellar Stream, and Inner Halo of M31

Abstract We measured [Fe/H] and [α/Fe] using spectral synthesis of low-resolution stellar spectroscopy for 70 individual red-giant-branch stars across four fields spanning the outer disk, Giant Stellar Stream (GSS), and inner halo of M31. Fields at M31-centric projected distances of 23 kpc in the halo, 12 kpc in the halo, 22 kpc in the GSS, and 26 kpc in the outer disk are α-enhanced, with [α/Fe] = 0.43, 0.50, 0.41, and 0.58, respectively. The 23 and 12 kpc halo fields are relatively metal-poor, with [Fe/H] = −1.54 and −1.30, whereas the 22 kpc GSS and 26 kpc outer disk fields are relatively metal-rich with [Fe/H] = −0.84 and −0.92, respectively. For fields with substructure, we separated the stellar populations into kinematically hot stellar halo components and kinematically cold components. We did not find any evidence of a radial [α/Fe] gradient along the high surface brightness core of the GSS between ∼17 and 22 kpc. However, we found tentative suggestions of a negative radial [α/Fe] gradient in the stellar halo, which may indicate that different progenitor(s) or formation mechanisms contributed to the build up of the inner versus outer halo. Additionally, the [α/Fe] distribution of the metal-rich ([Fe/H] > −1.5), smooth inner stellar halo (rproj ≲ 26 kpc) is inconsistent with having formed from the disruption of a progenitor(s) similar to present-day M31 satellite galaxies. The 26 kpc outer disk is most likely associated with the extended disk of M31, where its high α-enhancement provides support for an episode of rapid star formation in M31's disk possibly induced by a major merger.

The Density Structure of Simulated Stellar Streams

Abstract Star particles in a set of dense clusters are self-consistently evolved within an LCDM dark matter distribution with an n-body code. The clusters are started on nearly circular orbits in the more massive sub-halos. Each cluster develops a stellar tidal stream, initially within its original sub-halo. When a sub-halo merges into the main halo the early time stream is dispersed as a somewhat chaotic thick stream, roughly the width of the orbit of the cluster in the sub-halo. Once the cluster orbits freely in the main halo the star stream forms a thin stream again, usually resulting in a thin stream surrounded by a wider distribution of star particles lost at earlier times. To examine the role of the lower-mass dark matter sub-halos in the creation of density variations along the thin tidal star streams two realizations of the simulation are run with and without a normal cold dark matter sub-halo population below 4 × 108 . About 70(40)% of thin streams show density variations that are 2(5) times the star count noise level, irrespective of the presence or absence of low-mass sub-halos. A counts-in-cells analysis (related to the two-point correlation function and power spectrum) of the density along nearly 8000° of streams in the two well matched models finds that the full sub-halo population leads to slightly larger, but statistically significant, density fluctuations on scales of 2°–6°.

Variations in the Width, Density, and Direction of the Palomar 5 Tidal Tails

Abstract Stars that escape globular clusters form tidal tails that are predominantly shaped by the global distribution of mass in the Galaxy, but also preserve a historical record of small-scale perturbations. Using deep grz photometry from DECaLS, we present highly probable members of the tidal tails associated with the evaporating globular cluster Palomar 5. These data yield the cleanest view of a stellar stream beyond ≈20 kpc and reveal: (1) a wide, low surface-brightness extension of the leading tail; (2) significant density variations along the stream; and (3) sharp changes in the direction of both the leading and the trailing tail. In the fiducial Milky Way model, a rotating bar perturbs the Palomar 5 tails and can produce streams with similar width and density profiles to those observed. However, the deviations of the stream track in this simple model do not match those observed in the Palomar 5 trailing tail, indicating the need for an additional source of perturbation. These discoveries open up the possibility of measuring the population of perturbers in the Milky Way, including dark-matter subhalos, with an ensemble of stellar streams and deep photometry alone.

Metallicity and α-Element Abundance Gradients along the Sagittarius Stream as Seen by APOGEE

Abstract Using 3D positions and kinematics of stars relative to the Sagittarius (Sgr) orbital plane and angular momentum, we identify 166 Sgr stream members observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) that also have Gaia DR2 astrometry. This sample of 63/103 stars in the Sgr trailing/leading arm is combined with an APOGEE sample of 710 members of the Sgr dwarf spheroidal core (385 of them newly presented here) to establish differences of 0.6 dex in median metallicity and 0.1 dex in [α/Fe] between our Sgr core and dynamically older stream samples. Mild chemical gradients are found internally along each arm, but these steepen when anchored by core stars. With a model of Sgr tidal disruption providing estimated dynamical ages (i.e., stripping times) for each stream star, we find a mean metallicity gradient of 0.12 ± 0.03 dex Gyr−1 for stars stripped from Sgr over time. For the first time, an [α/Fe] gradient is also measured within the stream, at 0.02 ± 0.01 dex Gyr−1 using magnesium abundances and at 0.04 ± 0.01 dex Gyr−1 using silicon, which imply that the Sgr progenitor had significant radial abundance gradients. We discuss the magnitude of those inferred gradients and their implication for the nature of the Sgr progenitor within the context of the current family of Milky Way satellite galaxies, and we suggest that more sophisticated Sgr models are needed to properly interpret the growing chemodynamical detail we have on the Sgr system.

Elemental Abundances in M31: The Kinematics and Chemical Evolution of Dwarf Spheroidal Satellite Galaxies*

Abstract We present deep spectroscopy from Keck/DEIMOS of Andromeda I, III, V, VII, and X, all of which are dwarf spheroidal satellites of M31. The sample includes 256 spectroscopic members across all five dSphs. We confirm previous measurements of the velocity dispersions and dynamical masses, and we provide upper limits on bulk rotation. Our measurements confirm that M31 satellites obey the same relation between stellar mass and stellar metallicity as Milky Way (MW) satellites and other dwarf galaxies in the Local Group. The metallicity distributions show trends with stellar mass that are similar to those of MW satellites, including evidence in massive satellites for external influence, like pre-enrichment or gas accretion. We present the first measurements of individual element ratios, like [Si/Fe], in the M31 system, as well as measurements of the average [α/Fe] ratio. The trends of [α/Fe] with [Fe/H] also follow the same galaxy mass–dependent patterns as MW satellites. Less massive galaxies have more steeply declining slopes of [α/Fe] that begin at lower [Fe/H]. Finally, we compare the chemical evolution of M31 satellites to M31's Giant Stellar Stream and smooth halo. The properties of the M31 system support the theoretical prediction that the inner halo is composed primarily of massive galaxies that were accreted early. As a result, the inner halo exhibits higher [Fe/H] and [α/Fe] than surviving satellite galaxies.

The halo mass function in alternative dark matter models

ABSTRACT The claimed detection of large amounts of substructure in lensing flux anomalies, and in Milky Way stellar stream gap statistics, has led to a step change in constraints on simple warm dark matter models. In this study, we compute predictions for the halo mass function both for these simple models and for comprehensive particle physics models of sterile neutrinos and dark acoustic oscillations. We show that the mass function fit of Lovell et al. underestimates the number of haloes less massive than the half-mode mass, $M_\mathrm {hm}$, by a factor of 2, relative to the extended Press–Schechter (EPS) method. The alternative approach of applying EPS to the Viel et al. matter power spectrum fit instead suggests good agreement at $M_\mathrm {hm}$ relative to the comprehensive model matter power spectrum results, although the number of haloes with mass $\rm{\lt} M_\mathrm {hm}$ is still suppressed due to the absence of small-scale power in the fitting function. Overall, we find that the number of dark matter haloes with masses $\rm{\lt} 10^{8}{\, \rm M_\odot }$ predicted by competitive particle physics models is underestimated by a factor of ∼2 when applying popular fitting functions, although careful studies that follow the stripping and destruction of subhaloes will be required in order to draw robust conclusions.

Search for Stellar Streams in The Galactic Halo From Gaia DR2, GALAH DR2, RAVE DR5 and LAMOST DR4 Data

Stellar streams are stars which are trapped in the same potensial caused by dynamical resonance or tidal force. We aim to analyze kinematic sub- structures (streams) in the Galactic halo by V vs √ U2+ 2V 2 planes of Arifyanto & Fuchs. We crossmatched data from Gaia DR2, GALAH DR2, RAVE DR5 and LAMOST DR4 based on positions. We have 3D kinematics and metallicity data of halo stars selected from kinematics criteria from ratio of probability of thick disk (TD) over halo (H) less than 0.01. Substructures are detected by using wavelet transformation and corrected using 15 Monte Carlo simulations. We obtained four kinematic structures on V vs √U2 + 2V 2 plane which two of them are associated to BB17-1 and BB17-2 streams. All the streams had a high probability from the extragalactic origin.

Identifying stellar streams in Gaia DR2 with data mining techniques

ABSTRACT Streams of stars from captured dwarf galaxies and dissolved globular clusters are identifiable through the similarity of their orbital parameters, a fact that remains true long after the streams have dispersed spatially. We calculate the integrals of motion for 31 234 stars, to a distance of 4 kpc from the Sun, which have full and accurate 6D phase space positions in the Gaia DR2 catalogue. We then apply a novel combination of data mining, numerical, and statistical techniques to search for stellar streams. This process returns five high confidence streams (including one which was previously undiscovered), all of which display tight clustering in the integral of motion space. Colour–magnitude diagrams indicate that these streams are relatively simple, old, metal-poor populations. One of these resolved streams shares very similar kinematics and metallicity characteristics with the Gaia-Enceladus dwarf galaxy remnant, but with a slightly younger age. The success of this project demonstrates the usefulness of data mining techniques in exploring large data sets.

Effects of chaos on the detectability of stellar streams

ABSTRACT Observations show that stellar streams originating in satellite dwarf galaxies are frequent in the Universe. While such events are predicted by theory, it is not clear how many of the streams that are generated are washed out afterwards to the point in which it is impossible to detect them. Here, we study how these diffusion times are affected by the fact that typical gravitational potentials of the host galaxies can sustain chaotic orbits. We do this by comparing the behaviour of simulated stellar streams that reside in chaotic or non-chaotic regions of the phase space. We find that chaos does reduce the time interval in which streams can be detected. By analysing detectability criteria in configuration and velocity space, we find that the impact of these results on the observations depends on the quality of both the data and the underlying stellar halo model. For all the stellar streams, we obtain a similar upper limit to the detectable mass.

Globular clusters in the stellar stream surrounding the Milky Way analogue NGC 5907

ABSTRACT We study the globular clusters (GCs) in the spiral galaxy NGC 5907 well-known for its spectacular stellar stream – to better understand its origin. Using wide-field Subaru/Suprime-Cam gri images and deep Keck/DEIMOS multi-object spectroscopy, we identify and obtain the kinematics of several GCs superimposed on the stellar stream and the galaxy disc. We estimate the total number of GCs in NGC 5907 to be 154 ± 44, with a specific frequency of 0.73 ± 0.21. Our analysis also reveals a significant, new population of young star cluster candidates found mostly along the outskirts of the stellar disc. Using the properties of the stream GCs, we estimate that the disrupted galaxy has a stellar mass similar to the Sagittarius dwarf galaxy accreted by the Milky Way, i.e. $\sim 10^8~\rm M_\odot$.

Hunting ghosts: the iconic stellar stream(s) around NGC 5907 under scrutiny

Stellar streams are regarded as crucial objects for testing galaxy formation models because their morphology traces the underlying potentials and their occurrence tracks the assembly history of the galaxies. The existence of one of the most iconic stellar streams, the double loop around NGC 5907, has recently been questioned by new observations with the Dragonfly telescope. This new work only finds parts of the stream, even though a 1σ surface brightness limit of 30.3 mag arcsec−2 in the g band was reached. Using 7.2 h of luminance L-band imaging with the Milanković 1.4-m telescope, we have reobserved the putative double-loop part to confirm or reject this assessment. We do not find signs of the double loop, but see only a single knee-shaped stellar stream. Comparing our results to the data by the Dragonfly team, we find the same features. Our observations reach a 1σ surface brightness limit of 29.7 mag arcsec−2 in the g band. These findings emphasize the need for independent confirmation of detections of very low-surface brightness features.

Probing the merger history of red early-type galaxies with their faint stellar substructures

Several detailed observations, such as those carried out at the Canada-France-Hawaii Telescope (CFHT), have revealed prominent Low Surface Brightness (LSB) fine structures that lead to a change in the apparent morphology of galaxies. Previous photometry surveys have developed observational techniques which make use of the diffuse light detected in the external regions of galaxies. In these studies, the outer perturbations have been identified and classified. These include tidal tails, stellar streams, and shells. These structures serve as tracers for interacting events and merging events and retain some memory of the mass assembly of galaxies. Cosmological numerical simulations are required to estimate their visibility timescale, among other properties, in order to reconstruct the merger history of galaxies. In the present work, we analyze a hydrodynamical cosmological simulation to build up a comprehensive interpretation of the properties of fine structures. We present a census of several types of LSB fine structures compiled using a visual inspection of individual snapshots at various points in time. We reconstruct the evolution of the number of fine structures detected around an early-type galaxy and we compare it with the merger history of the galaxy. We find that most fine structures are associated with major and intermediate mass merger events. Their survival timescale ranges between 0.7 and 4 Gyr. Shells and streams remain visible for a longer time, while tidal tails have a shorter lifetime. These estimates for the survival time of collisional debris provide clues for the interpretation of the shape and frequency of fine structures observed in deep images with regard to their mass assembly. We find that the detectability of stellar streams is most sensitive at the surface brightness limit, demonstrating greater visibility at the deepest surface brightness level used in our simulation. We see between two and three times more streams based on a surface brightness cut of 33 mag arcsec−2 than with 29 mag arcsec−2. We find that the detection of shells is strongly dependent upon the projection angle.