No Arabic abstract
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.
We compare distance resolved, absolute proper motions in the Milky Way bar/bulge region to a grid of made-to-measure dynamical models with well defined pattern speeds. The data are obtained by combining the relative VVV Infrared Astrometric Catalog v1 proper motions with the Gaia DR2 absolute reference frame. We undertake a comprehensive analysis of the various errors in our comparison, from both the data and the models, and allow for additional, unknown, contributions by using an outlier-tolerant likelihood function to evaluate the best fitting model. We quantify systematic effects such as the region of data included in the comparison, with or without possible overlap from spiral arms, and the choice of synthetic luminosity function and bar angle used to predict the data from the models. Resulting variations in the best-fit parameters are included in the final error budget. We measure the bar pattern speed to be Omega_b=35.4+-0.9 km/s/kpc and the azimuthal solar velocity to be V_phi_sun= 251.4+-1.7 km/s. These values, when combined with recent measurements of the Galactic rotation curve, yield the distance of corotation, 6.3 < R_(CR) [kpc] < 6.8, the outer Lindblad resonance (OLR), 10.5 < R_(OLR) [kpc] < 11.5, and the higher order, m=4, OLR, 8.5 < R_(OLR_4) [kpc] < 9.0. The measured low pattern speed provides strong evidence for the long-slow bar scenario.
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.
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 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.