No Arabic abstract
We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 kms kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 kms kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.
Based on Gaia Early Data Release 3 (EDR3), we estimate the proper motions for 46 dwarf spheroidal galaxies (dSphs) of the Milky Way. The uncertainties in proper motions, determined by combining both statistical and systematic errors, are smaller by a factor 2.5, when compared with Gaia Data Release 2. We have derived orbits in four Milky Way potential models that are consistent with the MW rotation curve, with total mass ranging from $2.8times10^{11}$ $M_{odot}$ to $15times10^{11}$ $M_{odot}$. Although the type of orbit (ellipse or hyperbola) are very dependent on the potential model, the pericenter values are firmly determined, largely independent of the adopted MW mass model. By analyzing the orbital phases, we found that the dSphs are highly concentrated close to their pericenter, rather than to their apocenter as expected from Keplers law. This may challenge the fact that most dSphs are Milky Way satellites, or alternatively indicates an unexpected large number of undiscovered dSphs lying very close to their apocenters. Between half and two thirds of the satellites have orbital poles that indicate them to orbit along the Vast Polar Structure (VPOS), with the vast majority of these co-orbiting in a common direction also shared by the Magellanic Clouds, which is indicative of a real structure of dSphs.
The second data release from the Gaia mission (DR2) provides a comprehensive and unprecedented picture of the motions of astronomical sources in the plane of the sky, extending from the solar neighborhood to the outer reaches of the Milky Way. I present proper motion measurements based on Gaia DR2 for 17 ultra-faint dwarf galaxies within 100 kpc of the Milky Way. I compile the spectroscopically-confirmed member stars in each dwarf bright enough for Gaia astrometry from the literature, producing member samples ranging from 2 stars in Triangulum II to 68 stars in Bootes I. From the spectroscopic member catalogs I estimate the proper motion of each system. I find good agreement with the proper motions derived by the Gaia collaboration for Bootes I and Leo I. The tangential velocities for 14 of the 17 dwarfs are determined to better than 50 km/s, more than doubling the sample of such measurements for Milky Way satellite galaxies. The orbital pericenters are well-constrained, with a median value of 38 kpc. Only one satellite, Tucana III, is on an orbit passing within 15 kpc of the Galactic center, suggesting that the remaining ultra-faint dwarfs are unlikely to have experienced severe tidal stripping. As a group, the ultra-faint dwarfs are on high-velocity, eccentric, retrograde trajectories, with nearly all of them having space motions exceeding 370 km/s. In a low-mass (M_vir = 0.8 x 10^12 M_sun) Milky Way potential, eight out of the 17 galaxies lack well-defined apocenters and appear likely to be on their first infall, indicating that the Milky Way mass may be larger than previously estimated or that many of the ultra-faint dwarfs are associated with the Magellanic Clouds. The median eccentricity of the ultra-faint dwarf orbits is 0.79, similar to the values seen in numerical simulations, but distinct from the rounder orbits of the more luminous dwarf spheroidals.
We have derived absolute proper motions of the entire Galactic bulge region from VIRAC and Gaia. We present these as both integrated on-sky maps and, after isolating standard candle red clump (RC) stars, as a function of distance using RC magnitude as a proxy. These data provide a new global, 3-dimensional view of the Milky Way barred bulge kinematics. We find a gradient in the mean longitudinal proper motion, $<mu_l^star>$, between the different sides of the bar, which is sensitive to the bar pattern speed. The split RC has distinct proper motions and is colder than other stars at similar distance. The proper motion correlation map has a quadrupole pattern in all magnitude slices showing no evidence for a separate, more axisymmetric inner bulge component. The line-of-sight integrated kinematic maps show a high central velocity dispersion surrounded by a more asymmetric dispersion profile. $sigma_{mu_l} / sigma_{mu_b}$ is smallest, $sim1.1$, near the minor axis and reaches $sim1.4$ near the disc plane. The integrated $<mu_b>$ pattern signals a superposition of bar rotation and internal streaming motion, with the near part shrinking in latitude and the far part expanding. To understand and interpret these remarkable data, we compare to a made-to-measure barred dynamical model, folding in the VIRAC selection function to construct mock maps. We find that our model of the barred bulge, with a pattern speed of 37.5 $mathrm{km , s^{-1} , kpc^{-1}}$, is able to reproduce all observed features impressively well. Dynamical models like this will be key to unlocking the full potential of these data.
Observations of low-mass satellite galaxies in the nearby Universe point towards a strong dichotomy in their star-forming properties relative to systems with similar mass in the field. Specifically, satellite galaxies are preferentially gas poor and no longer forming stars, while their field counterparts are largely gas rich and actively forming stars. Much of the recent work to understand this dichotomy has been statistical in nature, determining not just that environmental processes are most likely responsible for quenching these low-mass systems but also that they must operate very quickly after infall onto the host system, with quenching timescales $lesssim 2~ {rm Gyr}$ at ${M}_{star} lesssim 10^{8}~{rm M}_{odot}$. This work utilizes the newly-available $Gaia$ DR2 proper motion measurements along with the Phat ELVIS suite of high-resolution, cosmological, zoom-in simulations to study low-mass satellite quenching around the Milky Way on an object-by-object basis. We derive constraints on the infall times for $37$ of the known low-mass satellite galaxies of the Milky Way, finding that $gtrsim~70%$ of the `classical satellites of the Milky Way are consistent with the very short quenching timescales inferred from the total population in previous works. The remaining classical Milky Way satellites have quenching timescales noticeably longer, with $tau_{rm quench} sim 6 - 8~{rm Gyr}$, highlighting how detailed orbital modeling is likely necessary to understand the specifics of environmental quenching for individual satellite galaxies. Additionally, we find that the $6$ ultra-faint dwarf galaxies with publicly available $HST$-based star-formation histories are all consistent with having their star formation shut down prior to infall onto the Milky Way -- which, combined with their very early quenching times, strongly favors quenching driven by reionization.
We estimate the mass of the Milky Way (MW) within 21.1 kpc using the kinematics of halo globular clusters (GCs) determined by Gaia. The second Gaia data release (DR2) contained a catalogue of absolute proper motions (PMs) for a set of Galactic GCs and satellite galaxies measured using Gaia DR2 data. We select from the catalogue only halo GCs, identifying a total of 34 GCs spanning $2.0 < r < 21.1$ kpc, and use their 3D kinematics to estimate the anisotropy over this range to be $beta = 0.46^{+0.15}_{-0.19}$, in good agreement, though slightly lower than, a recent estimate for a sample of halo GCs using HST PM measurements further out in the halo. We then use the Gaia kinematics to estimate the mass of the MW inside the outermost GC to be $M(< 21.1 mathrm{kpc}) = 0.21^{+0.04}_{-0.03} 10^{12} mathrm{M_odot}$, which corresponds to a circular velocity of $v_mathrm{circ}(21.1 mathrm{kpc}) = 206^{+19}_{-16}$ km/s. The implied virial mass is $M_mathrm{virial} = 1.28^{+0.97}_{-0.48} 10^{12} mathrm{M_odot}$. The error bars encompass the uncertainties on the anisotropy and on the density profile of the MW dark halo, and the scatter inherent in the mass estimator we use. We get improved estimates when we combine the Gaia and HST samples to provide kinematics for 46 GCs out to 39.5 kpc: $beta = 0.52^{+0.11}_{-0.14}$, $M(< 39.5 mathrm{kpc}) = 0.42^{+0.07}_{-0.06} 10^{12} mathrm{M_odot}$, and $M_mathrm{virial} = 1.54^{+0.75}_{-0.44} 10^{12} mathrm{M_odot}$. We show that these results are robust to potential substructure in the halo GC distribution. While a wide range of MW virial masses have been advocated in the literature, from below $10^{12} mathrm{M_odot}$ to above $2 times 10^{12}mathrm{M_odot}$, these new data imply that an intermediate mass is most likely.