Search for: All records

We analyse the M31 halo and its substructure within a projected radius of 120 kpc using a combination of Subaru/HSC NB515 and Canada France Hawaii Telescope/MegaCam g and i bands. We succeed in separating M31’s halo stars from foreground contamination with $\sim$90 per cent accuracy by using the surface gravity sensitive NB515 filter. Based on the selected M31 halo stars, we discover three new substructures, which associate with the Giant Southern Stream (GSS) based on their photometric metallicity estimates. We also produce the distance and photometric metallicity estimates for the known substructures. While these quantities for the GSS are reproduced in our study, we find that the north-western stream shows a steeper distance gradient than found in an earlier study, suggesting that it is likely to have formed in an orbit closer to the Milky Way. For two streams in the eastern halo (Stream C and D), we identify distance gradients that had not been resolved. Finally, we investigate the global halo photometric metallicity distribution and surface brightness profile using the NB515-selected halo stars. We find that the surface brightness of the metal-poor and metal-rich halo populations, and the all population can be fitted to a power-law profile with an index of $\alpha =-1.65\pm 0.02$, $-2.82\pm 0.01$, and $-2.44\pm 0.01$, respectively. In contrast to the relative smoothness of the halo profile, its photometric metallicity distribution appears to be spatially non-uniform with non-monotonic trends with radius, suggesting that the halo population had insufficient time to dynamically homogenize the accreted populations.

 

We map the 3D kinematics of the Galactic disc out to 3.5 kpc from the Sun, and within 0.75 kpc from the mid-plane of the Galaxy. To this end, we combine high-quality astrometry from Gaia EDR3, with heliocentric line-of-sight velocities from Gaia DR2, and spectroscopic surveys including APOGEE, GALAH, and LAMOST. We construct an axisymmetric model for the mean velocity field, and subtract this on a star-by-star basis to obtain the residual velocity field in the Galactocentric components (Vϕ, VR, Vz), and Vlos. The velocity residuals are quantified using the power spectrum, and we find that the peak power (A/[km s−1]) in the mid-plane (|z| < 0.25 kpc) is (Aϕ, AR, AZ, Alos) = (4.2,8.5,2.6,4.6), at 0.25 < |z|/[kpc] < 0.5, is (Aϕ, AR, AZ, Alos) = (4.0,7.9,3.6,5.3), and at 0.5 < |z|/[kpc] < 0.75, is (Aϕ, AR, AZ, Alos) = (1.9,6.9,5.2,6.4). Our results provide a sophisticated measurement of the streaming motion in the disc and in the individual components. We find that streaming is most significant in VR, and at all heights (|Z|) probed, but is also non-negligible in other components. Additionally, we find that patterns in velocity field overlap spatially with models for spiral arms in the Galaxy. Our simulations show that phase-mixing of disrupting spiral arms can generate such residuals in the velocity field, where the radial component is dominant, just as in real data. We also find that with time evolution, both the amplitude and physical scale of the residual motion decrease.

 

K2 was a community-driven NASA mission where all targets were proposed through guest observer programmes. Here we provide an overview of one of the largest of these endeavours, the K2 Galactic Archaeology Programme (K2GAP), with about 25 per cent of the observed targets being allocated to this programme. K2GAP provides asteroseismic parameters for about 23 000 giant stars across the Galaxy, which together with spectroscopic stellar parameters can give age and masses of stars. We discuss in detail the target selection procedure and provide a python program that implements the selection function (github.com/sanjibs/k2gap). Broadly speaking, the targets were selected on 2MASS colour J − Ks > 0.5, with finely tuned adjustments for each campaign. We discuss the detection completeness of the asteroseismic parameters νmax and Δν. About 14 per cent of giants were found to miss νmax detections and it was difficult to detect Δν for RC stars. Making use of the selection function, we compare the observed distribution of asteroseismic masses to theoretical predictions. The median asteroseismic mass is higher by about 4 per cent compared to predictions. We provide a selection-function-matched mock catalogue of stars based on a synthetic model of the Galaxy for the community to use in subsequent analyses of the K2GAP data set (physics.usyd.edu.au/k2gap/download/).

 

Abstract We present the first detailed comparison of populations of dwarf galaxy stellar streams in cosmological simulations and the Milky Way. In particular, we compare streams identified around 13 Milky Way analogs in the FIRE-2 simulations to streams observed by the Southern Stellar Stream Spectroscopic Survey ( S 5 ). For an accurate comparison, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to >110 kpc) and surviving only at larger apocenters (≳40 kpc) than those observed in the Milky Way. We find that the population of high-stellar-mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as intact satellites in DES-like observations, since their tidal tails have too low surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low-surface-brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations. 

ABSTRACT In the currently favoured cosmological paradigm galaxies form hierarchically through the accretion of satellites. Since a satellite is less massive than the host, its stars occupy a smaller volume in action space. Actions are conserved when the potential of the host halo changes adiabatically, so stars from an accreted satellite would remain clustered in action space as the host evolves. In this paper, we identify recently disrupted accreted satellites in three Milky Way-like disc galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellites through simulation snapshots. We try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the z = 0 snapshot. Even with completely error-free mock data we find that only 35 per cent (14/39) satellites are well recovered while the rest (25/39) are poorly recovered (i.e. either contaminated or split up). Most (10/14 ∼70 per cent) of the well-recovered satellites have infall times <7.1 Gyr ago and total mass >4 × 108M⊙ (stellar mass more than 1.2 × 106 M⊙, although our upper mass limit is likely to be resolution dependent). Since cosmological simulations predict that stellar haloes include a population of in situ stars, we test our ability to recover satellites when the data include 10–50 per cent in situ contamination. We find that most previously well-recovered satellites stay well recovered even with 50 per cent contamination. With the wealth of 6D phase space data becoming available we expect that cluster analysis in action space will be useful in identifying the majority of recently accreted and moderately massive satellites in the Milky Way. 

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}$.