We constructed a Hubble Space Telescope (HST) astro-photometric catalog of the central region of the Galactic globular cluster NGC 1261. This catalog, complemented with Gaia DR2 data sampling the external regions, has been used to estimate the structural parameters of the system (i.e., core, half-mass, tidal radii and concentration) from its resolved star density profile. We computed high-precision proper motions thanks to multi-epoch HST data and derived the cluster velocity dispersion profile in the plane of the sky for the innermost region, finding that the system is isotropic. The combination with line-of-sight information collected from spectroscopy in the external regions provided us with the cluster velocity dispersion profile along the entire radial extension. We also measured the absolute proper motion of NGC 1261 using a few background galaxies as a reference. The radial distribution of the Blue Straggler Star population shows that the cluster is in a low/intermediate phase of dynamical evolution.

The Wide-Field InfraRed Space Telescope (WFIRST) will be capable of delivering precise astrometry for faint sources over the enormous field of view of its main camera, the Wide-Field Imager (WFI). This unprecedented combination will be transformative for the many scientific questions that require precise positions, distances, and velocities of stars. We describe the expectations for the astrometric precision of the WFIRST WFI in different scenarios, illustrate how a broad range of science cases will see significant advances with such data, and identify aspects of WFIRST's design where small adjustments could greatly improve its power as an astrometric instrument.

In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

Studying our Galaxy, the Milky Way (MW), gives us a close-up view of the interplay between cosmology, dark matter, and galaxy formation. In the next decade our understanding of the MW's dynamics, stellar populations, and structure will undergo a revolution thanks to planned and proposed astrometric, spectroscopic and photometric surveys, building on recent advances by the Gaia astrometric survey. Together, these new efforts will measure three-dimensional positions and velocities and numerous chemical abundances for stars to the MW's edge and well into the Local Group, leading to a complete multidimensional view of our Galaxy. Studies of the multidimensional Milky Way beyond the Gaia frontier---from the edge of the Galactic disk to the edge of our Galaxy's dark matter halo---will unlock new scientific advances across astrophysics, from constraints on dark matter to insights into galaxy formation.

The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions.

We present the first detailed kinematic analysis of the proper motions (PMs) of stars in the Magellanic Bridge, from both the \textit{Gaia} Data Release 2 catalog and from \textit{Hubble Space Telescope} Advanced Camera for Surveys data. For the \textit{Gaia} data, we identify and select two populations of stars in the Bridge region, young main sequence (MS) and red giant stars. The spatial locations of the stars are compared against the known H {\small I} gas structure, finding a correlation between the MS stars and the H {\small I} gas. In the \textit{Hubble Space Telescope} fields our signal comes mainly from an older MS and turn-off population, and the proper motion baselines range between $\sim 4$ and 13 years. The PMs of these different populations are found to be consistent with each other, as well as across the two telescopes. When the absolute motion of the Small Magellanic Cloud is subtracted out, the residual Bridge motions display a general pattern of pointing away from the Small Magellanic Cloud towards the Large Magellanic Cloud. We compare in detail the kinematics of the stellar samples against numerical simulations of the interactions between the Small and Large Magellanic Clouds, and find general agreement between the kinematics of the observed populations and a simulation in which the Clouds have undergone a recent direct collision.

We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber.

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.

The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is confidently detected, consistent with the known line-of-sight rotation curves: $V_{\rm rot} = -206\pm86$ km s$^{-1}$ (counter-clockwise) for M31, and $V_{\rm rot} = 80\pm52$ km s$^{-1}$ (clockwise) for M33. We measure the center-of-mass PM of each galaxy relative to surrounding background quasars in DR2. This yields that $({\mu}_{\alpha*},{\mu}_{\delta})$ equals $(65 \pm 18 , -57 \pm 15)$ $\mu$as yr$^{-1}$ for M31, and $(31 \pm 19 , -29 \pm 16)$ $\mu$as yr$^{-1}$ for M33. In addition to the listed random errors, each component has an additional residual systematic error of 16 $\mu$as yr$^{-1}$. These results are consistent at 0.8$\sigma$ and 1.0$\sigma$ with the (2 and 3 times higher-accuracy) measurements already available from Hubble Space Telescope (HST) optical imaging and VLBA water maser observations, respectively. This lends confidence that all these measurements are robust. The new results imply that the M31 orbit towards the Milky Way is somewhat less radial than previously inferred, $V_{\rm tan, DR2+HST} = 57^{+35}_{-31}$ km s$^{-1}$, and strengthen arguments that M33 may be on its first infall into M31. The results highlight the future potential of Gaia for PM studies beyond the Milky Way satellite system.

The Orphan Stream is one of the most prominent tidal streams in the Galactic halo. Using data on red giants, RR Lyrae, and horizontal branch stars from Gaia and other surveys, we determine the proper motion of the Orphan Stream over a path of more than 90 degrees on the sky. We also provide updated tracks for the sky position, distance, and radial velocity of the stream. Our tracks in these latter dimensions mostly agree with previous results. However, there are significant corrections to the earlier distance and latitude tracks as the stream approaches the Galactic disk. Stream stars selected with three-dimensional kinematics display a very tight red giant sequence. Concordantly, we find that applying a proper motion cut removes the most metal-rich stars from earlier spectroscopic samples of stream stars, though a significant dispersion remains indicating a dwarf galaxy origin. The deceleration of the stream towards its leading end suggests a circular velocity of ~200 km/s at a galactocentric radius ~30 kpc, consistent with other independent evidence. However, the track of the stream departs significantly from an orbit; the spatial track does not point along the same direction as the velocity vector, and it exhibits a lateral wiggle that is unlikely to match any reasonable orbit. The low metallicity and small dispersion of the stream in the various coordinates point to a progenitor with a relatively low dynamical mass ~$10^8$ Msun.

We present a new measurement of the systemic proper motion of the Small Magellanic Cloud (SMC), based on an expanded set of 30 fields containing background quasars and spanning a $\sim$3 year baseline, using the \textit{Hubble Space Telescope} (\textit{HST}) Wide Field Camera 3. Combining this data with our previous 5 \textit{HST} fields, and an additional 8 measurements from the \textit{Gaia}-Tycho Astrometric Solution Catalog, brings us to a total of 43 SMC fields. We measure a systemic motion of $\mu_{W}$ = $-0.82$ $\pm$ 0.02 (random) $\pm$ 0.10 (systematic) mas yr$^{-1}$ and $\mu_{N}$ = $-1.21$ $\pm$ 0.01 (random) $\pm$ 0.03 (systematic) mas yr$^{-1}$. After subtraction of the systemic motion, we find little evidence for rotation, but find an ordered mean motion radially away from the SMC in the outer regions of the galaxy, indicating that the SMC is in the process of tidal disruption. We model the past interactions of the Clouds with each other based on the measured present-day relative velocity between them of $103 \pm 26$ km s$^{-1}$. We find that in 97\% of our considered cases, the Clouds experienced a direct collision $147 \pm 33$ Myr ago, with a mean impact parameter of $7.5 \pm 2.5$ kpc.

The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_{\phi} \rangle, \sigma_r, \sigma_{\phi}, \sigma_{\theta}$, and the velocity anisotropy $\beta$. Using the full HALO7D sample, we find $\beta=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $\beta$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $\beta$ in each of the individual HALO7D fields. We find that the velocity anisotropy $\beta$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $\beta$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $\beta$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $\beta$ variation patterns; spatially mapping $\beta$ is thus a way forward in characterizing the accretion history of the Galaxy.

According to LCDM theory, hierarchical evolution occurs on all mass scales, implying that satellites of the Milky Way should also have companions. The recent discovery of ultra-faint dwarf galaxy candidates in close proximity to the Magellanic Clouds provides an opportunity to test this theory. We present proper motion (PM) measurements for 13 of the 32 new dwarf galaxy candidates using Gaia data release 2. All 13 also have radial velocity measurements. We compare the measured 3D velocities of these dwarfs to those expected at the corresponding distance and location for the debris of an LMC analog in a cosmological numerical simulation. We conclude that 4 of these galaxies (Hor1, Car2, Car3 and Hyi1) have come in with the Magellanic Clouds, constituting the first confirmation of the type of satellite infall predicted by LCDM. Ret2, Tuc2 and Gru1 have velocity components that are not consistent within 3 sigma of our predictions and are therefore less favorable. Hya2 and Dra2 could be associated with the LMC and merit further attention. We rule out Tuc3, Cra2, Tri2 and Aqu2 as potential members. Of the dwarfs without measured PMs, 5 of them are deemed unlikely on the basis of their positions and distances alone as being too far from the orbital plane expected for LMC debris (Eri2, Ind2, Cet2, Cet3 and Vir1). For the remaining sample, we use the simulation to predict PMs and radial velocities, finding that Phx2 has an overdensity of stars in DR2 consistent with this PM prediction.

We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{\rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{\rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{\rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $\beta = 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (<39.5 \rm{kpc}) = 0.61^{+0.18}_{-0.12} \times 10^{12}$ M$_{\odot}$, and the virial mass $M_\rm{vir} = 2.05^{+0.97}_{-0.79} \times 10^{12}$ M$_{\odot}$, are consistent with, but on the high side among recent mass estimates in the literature.

We present a catalog of relative proper motions for 368,787 stars in the 30 Doradus region of the Large Magellanic Cloud (LMC), based on a dedicated two-epoch survey with the Hubble Space Telescope (HST) and supplemented with proper motions from our pilot archival study. We demonstrate that a relatively short epoch difference of 3 years is sufficient to reach a $\sim$0.1 mas yr$^{-1}$ level of precision or better. A number of stars have relative proper motions exceeding a 3-sigma error threshold, representing a mixture of Milky Way denizens and 17 potential LMC runaway stars. Based upon 183 VFTS OB-stars with the best proper motions, we conclude that none of them move faster than $\sim$0.3 mas yr$^{-1}$ in each coordinate -- equivalent to $\sim$70 km s$^{-1}$. Among the remaining 351 VFTS stars with less accurate proper motions, only one candidate OB runaway can be identified. We rule out any OB star in our sample moving at a tangential velocity exceeding $\sim$120 km s$^{-1}$. The most significant result of this study is finding 10 stars over wide range of masses, which appear to be ejected from the massive star cluster R136 in the tangential plane to angular distances from $35^{\prime\prime}$ out to $407^{\prime\prime}$, equivalent to 8-98 pc. The tangential velocities of these runaways appear to be correlated with apparent magnitude, indicating a possible dependence on the stellar mass.

Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic potentials. We find that debris at the apocenters must be dynamically young, in the sense of being stripped off in the last two pericentric passages, while the Sagittarius dwarf is currently experiencing a third passage. The ratio of apocenters is sensitive to both dynamical friction and the outer slope of the Galactic rotation curve. These dependences can be understood with simple regularities connecting the apocentric radii, circular velocities, and orbital period of the progenitor. The effect of dynamical friction on the stream can be constrained using substructure within the leading apocenter. Our ensembles of models are not intended as statistically proper fits to the stream. Nevertheless, out of the range of models we consider, we consistently find the mass within 100 kpc to be $\sim 7 \times 10^{11} \, M_{\odot}$, with a nearly flat rotation curve between 50 and 100 kpc. This points to a more extended Galactic halo than assumed in some current models. As in previous work, we find prolate or triaxial halos ease agreement with the track of the leading stream. We display the behavior of our models in various observational spaces and characterize the substructure expected within the stream. In particular, the young trailing stream visible near trailing apocenter should exhibit a tight trend of velocity with distance separate from the older debris, and we suggest that this will serve as an especially useful probe of the outer Galactic potential.

High precision proper motion (PM) measurements are available for approximately 20% of all known dwarf satellite galaxies of the Milky Way (MW). Here we extend the Bayesian framework of Patel et al. (2017b) to include all MW satellites with measured 6D phase space information and apply it with the Illustris-Dark simulation to constrain the MW's mass. Using the properties of each MW satellite individually, we find that the scatter among mass estimates is reduced when the magnitude of specific orbital angular momentum (j) is adopted rather than their combined instantaneous positions and velocities. We also find that high j satellites (i.e. Leo II) constrain the upper limits for the MW's mass and low j satellites rather than the highest speed satellites (i.e. Leo I and LMC), set the lower mass limits. When j of all classical satellites is used to simultaneously estimate the MW's mass, we conclude the halo mass is $0.85^{+0.23}_{-0.26}\times 10^{12} \rm \, M_{\odot}$ (including Sagittarius dSph) and $0.96^{+0.29}_{-0.28} \times 10^{12} \rm \, M_{\odot}$ (excluding Sagittarius dSph), cautioning that low j satellites on decaying orbits like Sagittarius dSph may bias the distribution. These estimates markedly reduce the current factor of two spread in the mass range of the MW. We also find a well-defined relationship between host halo mass and satellite j distribution, which yields the prediction that upcoming PMs for ultra-faint dwarfs should reveal j within $ 5\times10^3 - 10^4$ kpc km s$^{-1}$. This is a promising method to significantly constrain the cosmologically expected mass range for the MW and eventually M31 as more satellite PMs become available.

Metal-poor globular clusters (GCs) exhibit intriguing Al-Mg anti-correlations and possible Si-Al correlations, which are important clues to decipher the multiple-population phenomenon. NGC 5053 is one of the most metal-poor GCs in the nearby Universe, and has been suggested to be associated with the Sagittarius (Sgr) dwarf galaxy, due to its similarity in location and radial velocity with one of the Sgr arms. In this work, we simulate the orbit of NGC 5053, and argue against a physical connection between Sgr and NGC 5053. On the other hand, the Mg, Al, and Si spectral lines, which are difficult to detect in the optical spectra of NGC 5053 stars, have been detected in the near-infrared APOGEE spectra. We use three different sets of stellar parameters and codes to derive the Mg, Al, and Si abundances. Regardless of which method is adopted, we see a large Al variation, and a substantial Si spread. Along with NGC 5053, metal-poor GCs exhibit different Mg, Al, and Si variations. Moreover, NGC 5053 has the lowest cluster mass among the GCs that have been identified to exhibit an observable Si spread until now.

Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 \mu m}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals.

We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu_W,mu_N)_Scl = (-0.0296+/-0.0209,-0.1358 +/-0.0214) mas/yr. The implied Galactocentric velocity vectors for Draco and Sculptor have radial and tangential components: (V_rad,V_tan)_Dra = (-88.6,161.4) +/- (4.4,5.6) km/s; and (V_rad,V_tan)_Scl = (72.6,200.2) +/- (1.3,10.8) km/s. We study the detailed orbital history of both Draco and Sculptor via numerical orbit integrations. Orbital periods of Draco and Sculptor are found to be 1-2 and 2-5 Gyrs, respectively, accounting for uncertainties in the MW mass. We also study the influence of the Large Magellanic Cloud (LMC) on the orbits of Draco and Sculptor. Overall, the inclusion of the LMC increases the scatter in the orbital results. Based on our calculations, Draco shows a rather wide range of orbital parameters depending on the MW mass and inclusion/exclusion of the LMC, but Sculptor's orbit is very well constrained with its most recent pericentric approach to the MW being 0.3-0.4 Gyr ago. Our new PMs imply that the orbital trajectories of both Draco and Sculptor are confined within the Disk of Satellites (DoS), better so than implied by earlier PM measurements, and likely rule out the possibility that these two galaxies were accreted together as part of a tightly bound group.

We present a proper motion measurement for the halo globular cluster Pyxis, using HST/ACS data as the first epoch, and GeMS/GSAOI Adaptive Optics data as the second, separated by a baseline of about 5 years. This is both the first measurement of the proper motion of Pyxis and the first calibration and use of Multi-Conjugate Adaptive Optics data to measure an absolute proper motion for a faint, distant halo object. Consequently, we present our analysis of the Adaptive Optics data in detail. We obtain a proper motion of mu_alpha cos(delta)=1.09+/-0.31 mas/yr and mu_delta=0.68+/-0.29 mas/yr. From the proper motion and the line-of-sight velocity we find the orbit of Pyxis is rather eccentric with its apocenter at more than 100 kpc and its pericenter at about 30 kpc. We also investigate two literature-proposed associations for Pyxis with the recently discovered ATLAS stream and the Magellanic system. Combining our measurements with dynamical modeling and cosmological numerical simulations we find it unlikely Pyxis is associated with either system. We examine other Milky Way satellites for possible association using the orbit, eccentricity, metallicity, and age as constraints and find no likely matches in satellites down to the mass of Leo II. We propose that Pyxis probably originated in an unknown galaxy, which today is fully disrupted. Assuming that Pyxis is bound and not on a first approach, we derive a 68% lower limit on the mass of the Milky Way of 0.95*10^12 M_sun.

The Hubble Space Telescope recently celebrated 25 years of operation. Some of the first images of extragalactic optical jets were taken by HST in the mid-1990s; with time baselines on the order of 20 years and state-of-the-art astrometry techniques, we are now able to reach accuracies in proper-motion measurements on the order of a tenth of a milliarcsecond per year. We present the results of a recent HST program to measure the kiloparsec-scale proper motions of eleven nearby optical jets with Hubble, the first sample of its kind. When paired with VLBI proper-motion measurements on the parsec scale, we are now able to map the full velocity profile of these jets from near the black hole to the final deceleration as they extend out into and beyond the host galaxy. We see convincing evidence that weak-flavor jets (i.e., FR Is) have a slowly increasing jet speed up to 100 pc from the core, where superluminal components are first seen.

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields along the Orphan Stream. We determine absolute PMs of several individual stars per target field using established techniques that utilize distant background galaxies to define a stationary reference frame. Five Orphan Stream stars are identified in one of the four fields based on combined color-magnitude and PM information. The average PM is consistent with the existing model of the Orphan stream by Newberg et al. In addition to the Orphan stream stars, we detect stars that likely belong to other stellar streams. To identify which stellar streams these stars belong to, we examine the 2-d bulk motion of each group of stars on the sky by subtracting the PM contribution of the solar motion (which is a function of position on the sky and distance) from the observed PMs, and comparing the vector of net motion with the spatial extent of known stellar streams. By doing this, we identify candidate stars in the Sagittarius and Lethe streams, and a newly-found stellar stream at a distance of ~17 kpc, which we tentatively name the "Parallel stream". Together with our Sagittarius stream study (Sohn et al., 2015, ApJ, 803, 56), this work demonstrates that even in the Gaia era, HST will continue to be advantageous in measuring PMs of old stellar populations on a star-by-star basis, especially for distances beyond ~10 kpc.

We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions. The proper motions were measured using long baseline (5-7 years) multi-epoch HST/ACS photometry, and the LOS velocities were extracted from deep (5-6 hour integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic $(l,b,$ LOS) coordinate system are $\langle v^2_l \rangle^{1/2} = 138^{+43}_{-26}$ km/s, $\langle v^2_b \rangle^{1/2} = 88^{+28}_{-17}$ km/s, and $\langle v^2_{\rm{LOS}} \rangle^{1/2} = 91^{+27}_{-14}$ km/s. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter $\beta$ using 3D kinematics outside of the solar neighborhood. We find $\beta=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $\beta$ measurements by 2$\sigma$ ($\beta_{SN} \sim 0.5-0.7$). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to $\beta=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1\sigma$. The potential decrease in $\beta$ with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at $r \sim 25$ kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.

The radio galaxy 3C 273 hosts one of the nearest and best-studied powerful quasar jets. Having been imaged repeatedly by the Hubble Space Telescope (HST) over the past twenty years, it was chosen for an HST program to measure proper motions in the kiloparsec-scale resolved jets of nearby radio-loud active galaxies. The jet in 3C 273 is highly relativistic on sub-parsec scales, with apparent proper motions up to 15$c$ observed by VLBI (Lister et al., 2013). In contrast, we find that the kpc-scale knots are compatible with being stationary, with a mean speed of $-$0.2$\pm$0.5$c$ over the whole jet. Assuming the knots are packets of moving plasma, an upper limit of 1c implies a bulk Lorentz factor $\Gamma<$2.9. This suggests that the jet has either decelerated significantly by the time it reaches the kpc scale, or that the knots in the jet are standing shock features. The second scenario is incompatible with the inverse Compton off the Cosmic Microwave Background (IC/CMB) model for the X-ray emission of these knots, which requires the knots to be in motion, but IC/CMB is also disfavored in the first scenario due to energetic considerations, in agreement with the recent finding of Meyer & Georganopoulos (2014) which ruled out the IC/CMB model for the X-ray emission of 3C 273 via gamma-ray upper limits.

We have imaged with HST's WFC3/UVIS the central 2.7$\times$2.7 arcmin$^2$ region of the giant elliptical galaxy M 87, using the ultraviolet filter F275W. In combination with archival ACS/WFC data taken through the F606W and F814W filters, covering the same field, we have constructed integrated-light UV-optical colors and magnitudes for 1460 objects, most of which are believed to be globular clusters belonging to M 87. The purpose was to ascertain whether the multiple-populations syndrome, ubiquitous among Galactic globular clusters (GCs), exists also among the M 87 family of clusters. To achieve this goal, we sought those GCs with exceptionally blue UV-to-optical colors, because helium-enriched sub-populations produce a horizontal-branch morphology that is well populated at high effective temperature. For comparison, integrated, synthetic UV$-$optical and purely optical colors and magnitudes have been constructed for 45 Galactic GCs, starting from individual-star photometry obtained with the same instruments and the same filters. We identify a small group of M 87 clusters exhibiting a radial UV$-$optical color gradient, representing our best candidate GCs hosting multiple populations with extreme helium content. We also find that the central spatial distribution of the bluer GCs is flattened in a direction parallel to the jet, while the distribution of redder GCs is more spherical. We release to the astronomical community our photometric catalog in F275W, F606W and F814W bands and the high-quality image stacks in the same bands.

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields spanning 200 degrees along the Sagittarius (Sgr) stream: one trailing arm field, one field near the Sgr dwarf spheroidal tidal radius, and two leading arm fields. We determine absolute PMs of dozens of individual stars per field, using established techniques that use distant background galaxies as stationary reference frame. Stream stars are identified based on combined color-magnitude diagram and PM information. The results are broadly consistent with the few existing PM measurements for the Sgr galaxy and the trailing arm. However, our new results provide the highest PM accuracy for the stream to date, the first PM measurements for the leading arm, and the first PM measurements for individual stream stars; we also serendipitously determine the PM of the globular cluster NGC~6652. In the trailing-arm field, the individual PMs allow us to kinematically separate trailing-arm stars from leading-arm stars that are 360 degrees further ahead in their orbit. Also, in three of our fields we find indications that two distinct kinematical components may exist within the same arm and wrap of the stream. Qualitative comparison of the HST data to the predictions of the Law & Majewski and Penarrubia et al. N-body models show that the PM measurements closely follow the predicted trend with Sgr longitude. This provides a successful consistency check on the PM measurements, as well as on these N-body approaches (which were not tailored to fit any PM data).

High-precision astrometry throughout the Local Group is a unique capability of the Hubble Space Telescope (HST), with potential for transformative science, including constraining the nature of dark matter, probing the epoch of reionization, and understanding key physics of galaxy evolution. While Gaia will provide unparalleled astrometric precision for bright stars in the inner halo of the Milky Way, HST is the only current mission capable of measuring accurate proper motions for systems at greater distances (> 80 kpc), which represents the vast majority of galaxies in the Local Group. The next generation of proper-motion measurements will require long time baselines, spanning many years to decades and possibly multiple telescopes, combining HST with the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). However, the current HST allocation process is not conducive to such multi-cycle/multi-mission science, which will bear fruit primarily over many years. We propose an HST astrometry initiative to enable long-time-baseline, multi-mission science, which we suggest could be used to provide comprehensive kinematic measurements of all dwarf galaxies and high surface-density stellar streams in the Local Group with HST's Advanced Camera for Surveys (ACS) or Wide Field Camera 3 (WFC3). Such an initiative not only would produce forefront scientific results within the next 5 years of HST's life, but also would serve as a critical anchor point for future missions to obtain unprecedented astrometric accuracy, ensuring that HST leaves a unique and lasting legacy for decades to come.

We present deep $UBVI_C$ photometric data of the young open cluster Westerlund 2. An abnormal reddening law of $R_{V,cl}=4.14\pm0.08$ was found for the highly reddened early-type members ($E(B-V)\geq 1.45$), whereas a fairly normal reddening law of $R_{V,fg}=3.33\pm0.03$ was confirmed for the foreground early-type stars ($E(B-V)_{fg}<1.05$). The distance modulus was determined from zero-age main-sequence (ZAMS) fitting to the reddening-corrected colour-magnitude diagram of the early-type members to be $V_0-M_V=13.9\pm0.14$ (random error) $_{-0.1}^{+0.4}$ (the upper limit of systematic error) mag ($d = 6.0 \pm 0.4 _{-0.3}^{+1.2}$ kpc). To obtain the initial mass function, pre-main-sequence (PMS) stars were selected by identifying the optical counterparts of Chandra X-ray sources and mid-infrared emission stars from the Spitzer GLIMPSE source catalog. The initial mass function shows a shallow slope of $\Gamma=-1.1 \pm 0.1$ down to $\log m = 0.7$. The total mass of Westerlund 2 is estimated to be at least 7,400 $M_{\odot}$. The age of Westerlund 2 from the main-sequence turn-on and PMS stars is estimated to be $\lesssim$ 1.5 Myr. We confirmed the existence of a clump of PMS stars located $\sim1$ arcmin north of the core of Westerlund 2, but we could not find any clear evidence for an age difference between the core and the northern clump.

This is the first in a series of three papers describing a project with the Cosmic Origins Spectrograph on the Hubble Space Telescope to measure abundances of the neutral interstellar medium (ISM) in a sample of 9 nearby star-forming galaxies. The goal is to assess the (in)homogeneities of the multiphase ISM in galaxies where the bulk of metals can be hidden in the neutral phase, yet the metallicity is inferred from the ionized gas in the HII regions. The sample, spanning a wide range in physical properties, is to date the best suited to investigate the metallicity behavior of the neutral gas at redshift z=0. ISM absorption lines were detected against the far-ultraviolet spectra of the brightest star-forming region(s) within each galaxy. Here we report on the observations, data reduction, and analysis of these spectra. Column densities were measured by a multi-component line-profile fitting technique, and neutral-gas abundances were obtained for a wide range of elements. Several caveats were considered including line saturation, ionization corrections, and dust depletion. Ionization effects were quantified with `ad-hoc' CLOUDY models reproducing the complex photoionization structure of the ionized and neutral gas surrounding the UV-bright sources. An `average spectrum of a redshift z=0 star-forming galaxy' was obtained from the average column densities of unsaturated profiles of neutral-gas species. This template can be used as a powerful tool for studies of the neutral ISM at both low and high redshift.

The Hubble Space Telescope (HST) has proven to be uniquely suited for the measurement of proper motions (PMs) of stars and galaxies in the nearby Universe. Here we summarize the main results and ongoing studies of the HSTPROMO collaboration, which over the past decade has executed some two dozen observational and theoretical HST projects on this topic. This is continuing to revolutionize our dynamical understanding of many objects, including: globular clusters; young star clusters; stars and stellar streams in the Milky Way halo; Local Group galaxies, including dwarf satellite galaxies, the Magellanic Clouds, and the Andromeda galaxy; and AGN Black Hole Jets.

We report new results from an HST archival program to study proper motions in the optical jet of the nearby radio galaxy M87. Using over 13 years of archival imaging, we reach accuracies below 0.1c in measuring the apparent velocities of individual knots in the jet. We confirm previous findings of speeds up to 4.5c in the inner 6" of the jet, and report new speeds for optical components in the outer part of the jet. We find evidence of significant motion transverse to the jet axis on the order of 0.6c in the inner jet features, and superluminal velocities parallel and transverse to the jet in the outer knot components, with an apparent ordering of velocity vectors possibly consistent with a helical jet pattern. Previous results suggested a global deceleration over the length of the jet in the form of decreasing maximum speeds of knot components from HST-1 outward, but our results suggest that superluminal speeds persist out to knot C, with large differentials in very nearby features all along the jet. We find significant apparent accelerations in directions parallel and transverse to the jet axis, along with evidence for stationary features in knots D, E, and I. These results are expected to place important constraints on detailed models of kpc-scale relativistic jets.

We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations.

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M_MW). Despite Leo I's large Galacto-centric space velocity (200 km/s) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M_MW > 10^{12} M_sun at 95% confidence for a variety of Bayesian priors on M_MW. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M_MW would increase by 30%. Imposing a mass weighted LCDM prior, we find a median Milky Way virial mass of M_MW=1.6 x 10^{12} M_sun, with a 90% confidence interval of [1.0-2.4] x 10^{12} M_sun. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

Color distributions of globular clusters (GCs) in most massive galaxies are bimodal. Assuming linear color-to-metallicity conversions, bimodality is viewed as the presence of merely two GC subsystems with distinct metallicities, which serves as a critical backbone of various galaxy formation theories. Recent studies, however, revealed that the color-metallicity relations (CMRs) often used to derive GC metallicities (e.g., CMRs of g-z, V-I and C-T1) are in fact inflected. Such inflection can create bimodal color distributions if the underlying GC metallicity spread is simply broad as expected from the hierarchical merging paradigm of galaxy formation. In order to test the nonlinear-CMR scenario for GC color bimodality, the u-band photometry is proposed because the u-related CMRs (e.g., CMRs of u-g and u-z) are theoretically predicted to be least inflected and most distinctive among commonly used optical CMRs. Here, we present Hubble Space Telescope (HST)/WFC3 F336W (u-band) photometry of the GC system in M84, a giant elliptical in the Virgo galaxy cluster. Combining the u data with the existing HST ACS/WFC g and z data, we find that the u-z and u-g color distributions are different from the g-z distribution in a very systematic manner and remarkably consistent with our model predictions based on the nonlinear-CMR hypothesis. The results lend further confidence to validity of the nonlinear-CMR scenario as an explanation for GC color bimodality. There are some GC systems showing bimodal spectroscopic metallicity, and in such systems the inflected CMRs often create stronger bimodality in the color domain.

Based on long baseline (5-7 years) multi-epoch HST/ACS photometry, used previously to measure the proper motion of M31, we present the proper motions (PMs) of 13 main-sequence Milky Way halo stars. The sample lies at an average distance of r ~24 kpc from the Galactic center, with a root-mean-square spread of 6 kpc. At this distance, the median PM accuracy is 5 km/s. We devise a maximum likelihood routine to determine the tangential velocity ellipsoid of the stellar halo. The velocity second moments in the directions of the Galactic (l,b) system are < vl^2 >^{1/2} = 123 (+29, -23) km/s, and < vb^2 >^{1/2} = 83 (+24, -16) km/s. We combine these results with the known line-of-sight second moment, < vlos^2 >^{1/2} = 105 \pm 5$ km/s, at this < r > to study the velocity anisotropy of the halo. We find approximate isotropy between the radial and tangential velocity distributions, with anisotropy parameter beta = 0.0 (+0.2, -0.4). Our results suggest that the stellar halo velocity anisotropy out to r ~ 30 kpc is less radially biased than solar neighborhood measurements. This is opposite to what is expected from violent relaxation, and may indicate the presence of a shell-type structure at r ~ 24 kpc. With additional multi-epoch HST data, the method presented here has the ability to measure the transverse kinematics of the halo for more stars, and to larger distances. This can yield new improved constraints on the stellar halo formation mechanism, and the mass of the Milky Way.

We present ultraviolet (UV) integrated colors of 44 Galactic globular clusters (GGCs) observed with the Galaxy Evolution Explorer (GALEX) in both FUV and NUV bands. We find for the first time that GCs associated with the Sagittarius dwarf galaxy have (FUV-V) colors systematically redder than GGCs with the same metallicity. M31 GCs show almost the same UV colors as GGCs, while M87 are systematically bluer. We speculate about the presence of an interesting trend, suggesting that the UV color of GCs may be correlated with the mass of the host galaxy, in the sense that more massive galaxies possess bluer clusters.

We present wide-field multiband photometry of globular cluster (GC) systems in NGC 1399, NGC 1404, and NGC 1387 located at the central region of the Fornax galaxy cluster. Observation was carried out through U, B, V, and I bands, which marks one of the widest and deepest U-band studies on extragalactic GC systems. The present U-band photometry enables us to significantly reduce the contamination by a factor of two for faint sources (Vo~23.5). The main results based on some 2000 GC candidates around NGC 1399, NGC 1404, and NGC 1387 are as follows: (1) the GC system in each galaxy exhibits bimodal color distributions in all colors examined, but the shape of color histograms varies systematically depending on colors; (2) NGC 1399 shows that the mean colors of both blue and red GCs become bluer with increasing galactocentric radius; (3) NGC 1399 shows overabundance of GCs in the directions of NGC 1404 and NGC 1387, indicating their ongoing interactions; and (4) NGC 1399 also exhibits a ~0.5' offset between the center of the inner GC distribution and the galaxy's optical center, suggesting that NGC 1399 is not yet dynamically relaxed and may be undergoing merger events.

Westerlund 1 is the most important starburst cluster in the Galaxy due to its massive star content. We have performed BVIc and JKs photometry to investigate the initial mass function (IMF). By comparing the observed color with the spectral type - intrinsic color relation, we obtain the mean interstellar reddening of <E(B-V)>=4.19+/-0.23 and <E(J-Ks)>=1.70+/-0.21. Due to the heavy extinction toward the cluster, the zero-age main sequence fitting method based on optical photometry proved to be inappropriate for the distance determination, while the near-infrared photometry gave a reliable distance to the cluster -- 3.8 kpc from the empirical relation. Using the recent theoretical stellar evolution models with rotation, the age of the cluster is estimated to be 5.0+/-1.0 Myr. We derived the IMF in the massive part and obtained a fairly shallow slope of {\Gamma} = -0.8 +/- 0.1. The integration of the IMF gave a total mass for the cluster in excess of 5.0 x 10^4 solar mass. The IMF shows a clear radial variation indicating the presence of mass segregation. We also discuss the possible star formation history of Westerlund 1 from the presence of red supergiants and relatively low-luminosity yellow hypergiants.

We present ultraviolet (UV) integrated colors of 44 Galactic globular clusters (GGCs) observed with the Galaxy Evolution Explorer (GALEX) in both FUV and NUV bands. This data-base is the largest homogeneous catalog of UV colors ever published for stellar systems in our Galaxy. The proximity of GGCs makes it possible to resolve many individual stars even with the somewhat low spatial resolution of GALEX. This allows us to determine how the integrated UV colors are driven by hot stellar populations, primarily horizontal branch stars and their progeny. The UV colors are found to be correlated with various parameters commonly used to define the horizontal branch morphology. We also investigate how the UV colors vary with parameters like metallicity, age, helium abundance and concentration. We find for the first time that GCs associated with the Sagittarius dwarf galaxy have (FUV-V) colors systematically redder than GGCs with the same metallicity. Finally, we speculate about the presence of an interesting trend, suggesting that the UV color of GCs may be correlated with the mass of the host galaxy, in the sense that more massive galaxies possess bluer clusters.

We study the future orbital evolution and merging of the MW-M31-M33 system, using a combination of collisionless N-body simulations and semi-analytic orbit integrations. Monte-Carlo simulations are used to explore the consequences of varying the initial phase-space and mass parameters within their observational uncertainties. The observed M31 transverse velocity implies that the MW and M31 will merge t = 5.86 (+1.61-0.72) Gyr from now, after a first pericenter at t = 3.87 (+0.42-0.32) Gyr. M31 may (probability p=41%) make a direct hit with the MW (defined here as a first-pericenter distance less than 25 kpc). Most likely, the MW and M31 will merge first, with M33 settling onto an orbit around them. Alternatively, M33 may make a direct hit with the MW first (p=9%), or M33 may get ejected from the Local Group (p=7%). The MW-M31 merger remnant will resemble an elliptical galaxy. The Sun will most likely (p=85%) end up at larger radius from the center of the MW-M31 merger remnant than its current distance from the MW center, possibly further than 50 kpc (p=10%). The Sun may (p=20%) at some time in the next 10 Gyr find itself moving through M33 (within 10 kpc), but while dynamically still bound to the MW-M31 merger remnant. The arrival and possible collision of M31 (and possibly M33) with the MW is the next major cosmic event affecting the environment of our Sun and solar system that can be predicted with some certainty. (Abridged)

We present the first proper motion measurements for the galaxy M31. We obtained new V-band imaging data with the HST ACS/WFC and WFC3/UVIS of a spheroid field near the minor axis, an outer disk field along the major axis, and a field on the Giant Southern Stream. The data provide 5-7 year time baselines with respect to pre-existing deep first-epoch observations. We measure the positions of thousands of M31 stars and hundreds of compact background galaxies in each field. High accuracy and robustness is achieved by building and fitting a unique template for each individual object. The average proper motion for each field is obtained from the average motion of the M31 stars between the epochs with respect to the background galaxies. For the three fields, the observed proper motions (mu_W,mu_N) are (-0.0458, -0.0376), (-0.0533, -0.0104), and (-0.0179,-0.0357) mas/yr, respectively. The ability to average over large numbers of objects and over the three fields yields a final accuracy of 0.012 mas/yr. The robustness of the proper-motion measurements and uncertainties are supported by the fact that data from different instruments, taken at different times and with different telescope orientations, as well as measurements of different fields, all yield statistically consistent results. Papers II and III explore the implications for our understanding of the history, future, and mass of the Local Group. (Abridged)

We determine the velocity vector of M31 with respect to the Milky Way and use this to constrain the mass of the Local Group, based on HST proper-motion measurements presented in Paper I. We construct N-body models for M31 to correct the measurements for the contributions from stellar motions internal to M31. We also estimate the center-of-mass motion independently, using the kinematics of satellite galaxies of M31 and the Local Group. All estimates are mutually consistent, and imply a weighted average M31 heliocentric transverse velocity of (v_W,v_N) = (-125.2+/-30.8, -73.8+/-28.4) km/s. We correct for the reflex motion of the Sun using the most recent insights into the solar motion within the Milky Way. This implies a radial velocity of M31 with respect to the Milky Way of V_rad = -109.3+/-4.4 km/s, and a tangential velocity V_tan = 17.0 km/s (<34.3 km/s at 1-sigma confidence). Hence, the velocity vector of M31 is statistically consistent with a radial (head-on collision) orbit towards the Milky Way. We revise prior estimates for the Local Group timing mass, including corrections for cosmic bias and scatter. Bayesian combination with other mass estimates yields M_LG = M_MW(vir) + M_M31(vir) = (3.17 +/- 0.57) x 10^12 solar masses. The velocity and mass results imply at 95% confidence that M33 is bound to M31, consistent with expectation from observed tidal deformations. (Abridged)

We present GALEX data for 44 Galactic globular clusters obtained during 3 GALEX observing cycles between 2004 and 2008. This is the largest homogeneous data set on the UV photometric properties of Galactic globular clusters ever collected. The sample selection and photometric analysis are discussed, and color-magnitude diagrams are presented. The blue and intermediate-blue horizontal branch is the dominant feature of the UV color-magnitude diagrams of old Galactic globular clusters. Our sample is large enough to display the remarkable variety of horizontal branch shapes found in old stellar populations. Other stellar types that are obviously detected are blue stragglers and post core-He burning stars. The main features of UV color-magnitude diagrams of Galactic globular clusters are briefly discussed. We establish the locus of post-core He burning stars in the UV color-magnitude diagram and present a catalog of candidate AGB-manqu \'e, post early-AGB, and post-AGB stars within our cluster sample.

We present a comprehensive catalog of 700 confirmed star clusters in the field of M31 compiled from three major existing catalogs. We detect 418 and 257 star clusters in Galaxy Evolution Explorer (GALEX) near-ultraviolet (NUV) and far-ultraviolet (FUV) imaging, respectively. Our final catalog includes photometry of star clusters in up to 16 passbands ranging from FUV to NIR as well as ancillary information such as reddening, metallicity, and radial velocities. In particular, this is the most extensive and updated catalog of UV integrated photometry for M31 star clusters. Ages and masses of star clusters are derived by fitting the multi-band photometry with model spectral energy distribution (SED); UV photometry enables more accurate age estimation of young clusters. Our catalog includes 182 young clusters with ages less than 1 Gyr. Our estimated ages and masses of young clusters are in good agreement with previously determined values in the literature. The mean age and mass of young clusters are about 300 Myr and 10^4 M_sun, respectively. We found that the compiled [Fe/H] values of young clusters included in our catalog are systematically lower (by more than 1 dex) than those from recent high-quality spectroscopic data and our SED fitting result. We confirm that most of the young clusters kinematics show systematic rotation around the minor axis and association with the thin disk of M31. The young clusters distribution exhibits a distinct peak in the M31 disk around 10 - 12 kpc from the center and follow a spatial distributions similar to other tracers of disk structure such as OB stars, UV star-forming regions, and dust. Some young clusters also show concentration around the ring splitting regions found in the southern part of the M31 disk and most of them have systematically younger (< 100 Myr) ages.

One of the conundrums in extragalactic astronomy is the discrepancy in observed metallicity distribution functions (MDFs) between the two prime stellar components of early-type galaxies-globular clusters (GCs) and halo field stars. This is generally taken as evidence of highly decoupled evolutionary histories between GC systems and their parent galaxies. Here we show, however, that new developments in linking the observed GC colors to their intrinsic metallicities suggest nonlinear color-to-metallicity conversions, which translate observed color distributions into strongly-peaked, unimodal MDFs with broad metal-poor tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved field stars in nearby elliptical galaxies and those produced by chemical evolution models of galaxies. The GC MDF shape, characterized by a sharp peak with a metal-poor tail, indicates a virtually continuous chemical enrichment with a relatively short timescale. The characteristic shape emerges across three orders of magnitude in the host galaxy mass, suggesting a universal process of chemical enrichment among various GC systems. Given that GCs are bluer than field stars within the same galaxy, it is plausible that the chemical enrichment processes of GCs ceased somewhat earlier than that of field stellar population, and if so, GCs preferentially trace the major, vigorous mode of star formation events in galactic formation. We further suggest a possible systematic age difference among GC systems, in that the GC systems in more luminous galaxies are older. This is consistent with the downsizing paradigm of galaxies and supports additionally the similar nature shared by GCs and field stars. Our findings suggest that GC systems and their parent galaxies have shared a more common origin than previously thought, and hence greatly simplify theories of galaxy formation.

The optical color distributions of globular clusters (GCs) in most large elliptical galaxies are bimodal. Based on the assumed linear relationship between GC colors and their metallicities, the bimodality has been taken as evidence of two GC subsystems with different metallicities in each galaxy and led to a number of theories in the context of galaxy formation. More recent observations and modeling of GCs, however, suggest that the color-metallicity relations (CMRs) are inflected, and thus colors likely trace metallicities in a nonlinear manner. The nonlinearity could produce bimodal color distributions from a broad underlying metallicity spread, even if it is unimodal. Despite the far-reaching implications, whether CMRs are nonlinear and whether the nonlinearity indeed causes the color bimodality are still open questions. Given that the spectroscopic refinement of CMRs is still very challenging, we here propose a new photometric technique to probe the possible nonlinear nature of CMRs. In essence, a color distribution of GCs is a "projected" distribution of their metallicities. Since the form of CMRs hinges on which color is used, the shape of color distributions varies depending significantly on the colors. Among other optical colors, the u-band related colors (e.g., u - g and u - z) are theoretically predicted to exhibit significantly less inflected CMRs than other preferred CMRs (e.g., for g - z). As a case study, we performed the HST/WFPC2 archival u-band photometry for the M87 (NGC 4486) GC system with confirmed color bimodality. We show that the u-band color distributions are significantly different from that of g - z, and consistent with our model predictions. With more u-band measurements, this method will support or rule out the nonlinear-CMR scenario for the origin of GC color bimodality with high confidence. [Abridged]

We derive the star formation histories of eight dwarf spheroidal (dSph) Milky Way satellite galaxies from their alpha element abundance patterns. Nearly 3000 stars from our previously published catalog (Paper II) comprise our data set. The average [alpha/Fe] ratios for all dSphs follow roughly the same path with increasing [Fe/H]. We do not observe the predicted knees in the [alpha/Fe] vs. [Fe/H] diagram, corresponding to the metallicity at which Type Ia supernovae begin to explode. Instead, we find that Type Ia supernova ejecta contribute to the abundances of all but the most metal-poor ([Fe/H] < -2.5) stars. We have also developed a chemical evolution model that tracks the star formation rate, Types II and Ia supernova explosions, and supernova feedback. Without metal enhancement in the supernova blowout, massive amounts of gas loss define the history of all dSphs except Fornax, the most luminous in our sample. All six of the best-fit model parameters correlate with dSph luminosity but not with velocity dispersion, half-light radius, or Galactocentric distance.

We present a catalog of Fe, Mg, Si, Ca, and Ti abundances for 2961 red giant stars that are likely members of eight dwarf satellite galaxies of the Milky Way (MW): Sculptor, Fornax, Leo I, Sextans, Leo II, Canes Venatici I, Ursa Minor, and Draco. For the purposes of validating our measurements, we also observed 445 red giants in MW globular clusters and 21 field red giants in the MW halo. The measurements are based on Keck/DEIMOS medium-resolution spectroscopy combined with spectral synthesis. We estimate uncertainties in [Fe/H] by quantifying the dispersion of [Fe/H] measurements in a sample of stars in monometallic globular clusters. We estimate uncertainties in Mg, Si, Ca, and Ti abundances by comparing our medium-resolution spectroscopic measurements to high-resolution spectroscopic abundances of the same stars. For this purpose, our DEIMOS sample included 132 red giants with published high-resolution spectroscopy in globular clusters, the MW halo field, and dwarf galaxies. The standard deviations of the differences in [Fe/H] and [alpha/Fe] (the average of [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe]) between the two samples is 0.15 and 0.16, respectively. This catalog represents the largest sample of multi-element abundances in dwarf galaxies to date. The next papers in this series draw conclusions on the chemical evolution, gas dynamics, and star formation histories from the catalog presented here. The wide range of dwarf galaxy luminosity reveals the dependence of dwarf galaxy chemical evolution on galaxy stellar mass.

We present ultraviolet (UV) color-magnitude relations (CMRs) of early-type dwarf galaxies in the Virgo cluster, based on Galaxy Evolution Explorer (GALEX) UV and Sloan Digital Sky Survey (SDSS) optical imaging data. We find that dwarf lenticular galaxies (dS0s), including peculiar dwarf elliptical galaxies (dEs) with disk substructures and blue centers, show a surprisingly distinct and tight locus separated from that of ordinary dEs, which is not clearly seen in previous CMRs. The dS0s in UV CMRs follow a steeper sequence than dEs and show bluer UV-optical color at a given magnitude. We also find that the UV CMRs of dEs in the outer cluster region are slightly steeper than that of their counterparts in the inner region, due to the existence of faint, blue dEs in the outer region. We explore the observed CMRs with population models of a luminosity-dependent delayed exponential star formation history. We confirm that the feature of delayed star formation of early-type dwarf galaxies in Virgo cluster is strongly correlated with their morphology and environment. The observed CMR of dS0s is well matched by models with relatively long delayed star formation. Our results suggest that dS0s are most likely transitional objects at the stage of subsequent transformation of late-type progenitors to ordinary red dEs in the cluster environment. In any case, UV photometry provides a powerful tool to disentangle the diverse subpopulations of early-type dwarf galaxies and uncover their evolutionary histories.

Wide-field deep gri images obtained with the Megacam of the Canada-France-Hawaii Telescope (CFHT) are used to investigate the spatial configuration of stars around five metal-poor globular cluster M15, M30, M53, NGC 5053, and NGC 5466, in a field-of-view ~3 degree. Applying a mask filtering algorithm to the color-magnitude diagrams of the observed stars, we sorted cluster's member star candidates that are used to examine the characteristics of the spatial stellar distribution surrounding the target clusters. The smoothed surface density maps and the overlaid isodensity contours indicate that all of the five metal-poor globular clusters exhibit strong evidence of extratidal overdensity features over their tidal radii, in the form of extended tidal tails around the clusters. The orientations of the observed extratidal features show signatures of tidal tails tracing the clusters' orbits, inferred from their proper motions, and effects of dynamical interactions with the Galaxy. Our findings include detections of a tidal bridge-like feature and an envelope structure around the pair of globular clusters M53 and NGC 5053. The observed radial surface density profiles of target clusters have a deviation from theoretical King models, for which the profiles show a break at 0.5~0.7r_t, extending the overdensity features out to 1.5~2r_t. Both radial surface density profiles for different angular sections and azimuthal number density profiles confirm the overdensity features of tidal tails around the five metal-poor globular clusters. Our results add further observational evidence that the observed metal-poor halo globular clusters originate from an accreted satellite system, indicative of the merging scenario of the formation of the Galactic halo.

There is a growing body of evidence for the presence of multiple stellar populations in some globular clusters, including NGC 1851. For most of these peculiar globular clusters, however, the evidence for the multiple red giant-branches (RGBs) having different heavy elemental abundances as observed in Omega Centauri is hitherto lacking, although spreads in some lighter elements are reported. It is therefore not clear whether they also share the suggested dwarf galaxy origin of Omega Cen or not. Here we show from the CTIO 4m UVI photometry of the globular cluster NGC 1851 that its RGB is clearly split into two in the U - I color. The two distinct RGB populations are also clearly separated in the abundance of heavy elements as traced by Calcium, suggesting that the type II supernovae enrichment is also responsible, in addition to the pollutions of lighter elements by intermediate mass asymptotic giant branch stars or fast-rotating massive stars. The RGB split, however, is not shown in the V - I color, as indicated by previous observations. Our stellar population models show that this and the presence of bimodal horizontal-branch distribution in NGC 1851 can be naturally reproduced if the metal-rich second generation stars are also enhanced in helium.

We explore the age distribution of the globular cluster (GC) system of the nearby elliptical galaxy NGC 5128 using ultraviolet (UV) photometry from Galaxy Evolution Explorer (GALEX) observations, with UV - optical colors used as the age indicator. Most GCs in NGC 5128 follow the general trends of GCs in M31 and Milky Way in UV - optical color-color diagram, which indicates that the majority of GCs in NGC 5128 are old similar to the age range of old GCs in M31 and Milky Way. A large fraction of spectroscopically identified intermediate-age GC (IAGC) candidates with ~ 3-8 Gyr are not detected in the FUV passband. Considering the nature of intermediate-age populations being faint in the far-UV (FUV) passband, we suggest that many of the spectroscopically identified IAGCs may be truly intermediate in age. This is in contrast to the case of M31 where a large fraction of spectroscopically suggested IAGCs are detected in FUV and therefore may not be genuine IAGCs but rather older GCs with developed blue horizontal branch stars. Our UV photometry strengthens the results previously suggesting the presence of GC and stellar subpopulation with intermediate age in NGC 5128. The existence of IAGCs strongly indicates the occurrence of at least one more major star formation episode after a starburst at high redshift.

We present a wide-field (4.5 deg^2) photometric and spectroscopic survey of the Leo I dwarf spheroidal (dSph) galaxy to explore its extended morphology and dynamics. As in previous papers in this series, we take advantage of photometry in the M, T_2, and DDO51 filter system to select LeoI red giant branch star candidates, and, so far, this selection technique has proven 100% reliable in selecting actual Leo I members among more than 100 M < 21.5 Leo I giant candidates having previous or new Keck DEIMOS spectroscopy to a radius >1.3 times the limiting radius of the fitted, central King profile. The two-dimensional distribution of all similarly-selected Leo I giant candidates is well fitted by a central single-component King profile of limiting radius 13.3 arcmin, but many giant stars are found outside this newly derived King limiting radius. The density profile thus shows a break at a major axis radial distance of ~10 arcmin produced by an excess of stars at and beyond the King limiting radius (spectroscopically confirmed to be made of true Leo I members), and primarily along the major axis of the main body of the rather elongated satellite. This spatial configuration, a rather flat velocity dispersion profile and an asymmetric radial velocity (RV) distribution among the Leo I members at large radii together support a picture where Leo I has been tidally disrupted on at least one, but at most two, perigalactic passages of a massive Local Group member. (abridged)

We study the UV and optical properties of 38 massive GCs in M87, imaged using the HST/STIS and WFPC2 instruments. The majority of these GCs appear extremely bright in the FUV - roughly a magnitude brighter than their Galactic counterparts with similar metallicities. The observed FUV flux is several times larger than predictions of canonical old stellar population models. These canonical models, which assume a fiducial helium enrichment parameter, dY/dZ=2, are able to reproduce the observed FUV fluxes only if ages are 3--5 Gyr larger than the WMAP age of the Universe, although the same models fit the UV photometry of Galactic and M31 GCs for ages less than the WMAP age. A similar discrepancy is found between the mass-weighted and UV-luminosity weighted ages of the massive Galactic GC omega Cen, whose CMD (including peculiar features on its well-populated horizontal branch) can be accurately reproduced by invoking a small super-He-rich (dY/dZ ~ 90) stellar component. By comparison to omega Cen, we propose that the majority of M87 GCs in our sample contain strong signatures of similarly minor super-He-rich sub-components. Although we cannot prove or disprove the extreme helium scenario at the moment, we show that the same phenomenon that causes the extended horizontal branch of omega Cen explains the UV brightness of our sample. If this is indeed due to the extreme helium, this study would be the first to find its signatures in extragalactic objects.

We present ultraviolet (UV) photometry of M31 globular clusters (GCs) found in 23 Galaxy Evolution Explorer (GALEX) images covering the entirety of M31. We detect 485 and 273 GCs (and GC candidates) in the near-ultraviolet (NUV; 2267 A) and far-ultraviolet (FUV; 1516 A), respectively. Comparing M31 data with those of Galactic GCs in the UV with the aid of population models, we find that the age ranges of old GCs in M31 and the Galactic halo are similar. Three metal-rich ([Fe/H]>-1) GCs in M31 produce significant FUV flux making their FUV-V colors unusually blue for their metallicities. These are thought to be analogs of the two peculiar Galactic GCs NGC 6388 and NGC 6441 with extended blue HB stars. Based on the models incorporating helium enriched subpopulations in addition to the majority of the population that have a normal helium abundance, we suggest that even small fraction of super-helium-rich subpopulations in GCs can reproduce the observed UV bright metal-rich GCs. Young clusters in M31 show distinct UV and optical properties from GCs in Milky Way. Population models indicate that their typical age is less than ~ 2 Gyrs. A large fraction of young GCs have the kinematics of the thin, rapidly rotating disk component. However, a subset of the old GCs also shares the thin-disk kinematics of the younger clusters. The existence of young GCs on the outskirts of M31 disk suggests the occurrence of a significant recent star formation in the thin-disk of M31. Old thin-disk GCs may set constraints on the epoch of early formation of the M31 thin-disk. We detect 12 (10) intermediate-age GC candidates in NUV (FUV). We suggest that some of spectroscopically identified intermediate-age GCs may not be truly intermediate in age, but rather older GCs that possess developed HB.

We present JHK-band near-infrared photometry of star clusters in the dwarf irregular galaxy IC 5152. After excluding possible foreground stars, a number of candidate star clusters are identified in the near-infrared images of IC 5152, which include young populations. Especially, five young star clusters are identified in the (J-H, H-K) two color diagram and the total extinction values toward these clusters are estimated to be A_V =2 - 6 from the comparison with the theoretical values given by the Leitherer et al. (1999)'s theoretical star cluster model.

A new large-area Washington M,T_2+DDO51 filter survey of more than 10 deg^2 around the Carina dSph galaxy reveals a spectroscopically confirmed power law radial density "break" population of Carina giant stars extending several degrees beyond the central King profile. Magellan telescope MIKE spectroscopy establishes the existence of Carina stars to at least 4.5 times its central King limiting radius, r_lim and primarily along Carina's major axis. To keep these stars bound to the dSph would require a global Carina mass-to-light ratio of M/L > 6,300 M/L_sun. The MIKE velocities, supplemented with ~950 additional Carina field velocities from archived VLT+GIRAFFE spectra with r<=r_lim, demonstrate a nearly constant Carina velocity dispersion to just beyond r = r_lim, and both a rising velocity dispersion and a velocity shear at still larger radii. Together, the observational evidence suggests that the discovered extended Carina population represents tidal debris from the dSph. Of 65 giant candidates at large angular radii from the Carina center for which MIKE spectra have been obtained 94% are associated either with Carina or a second, newly discovered diffuse, but strongly radial velocity-coherent (velocity dispersion of 9.8 km s^-1), foreground halo system. The fifteen stars in this second, retrograde velocity population have (1) a mean metallicity ~1 dex higher than that of Carina, and (2) colors and magnitudes consistent with the red clump of the Large Magellanic Cloud (LMC). Additional spectroscopy of giant star candidates in fields linking Carina and the LMC shows a smooth velocity gradient between the LMC and the retrograde Carina moving group. We conclude that we have found Magellanic stars almost twice as far (22 deg) from the LMC center than previously known.

We present BV CCD photometry for the open clusters Czernik 24 and Czernik 27. These clusters have never been studied before, and we provide, for the first time, the cluster parameters; reddening, distance, metallicity and age. Czernik 24 is an old open cluster with age 1.8 +/- 0.2 Gyr, metallicity [Fe/H]=-0.41 +/- 0.15 dex, distance modulus (m-M)_0 = 13.1 +/- 0.3 mag (d=4.1 +/- 0.5 kpc), and reddening E(B-V) = 0.54 +/- 0.12 mag. The parameters for Czernik 27 are estimated to be age = 0.63 +/- 0.07 Gyr, [Fe/H]= -0.02 +/- 0.10 dex, (m-M)_0 = 13.8 +/- 0.2 mag (d=5.8 +/- 0.5 kpc), and E(B-V) = 0.15 +/- 0.05 mag. The metallicity and distance values for Czernik 24 are consistent with the relation between the metallicity and the Galactocentric distance of other old open clusters. We find the metallicity gradient of 51 old open clusters including Czernik 24 to be Delta [Fe/H]/Delta R_gc= -0.064 +/- 0.009 dex/kpc.

We analyze velocity dispersion profiles for the Draco and Ursa Minor (UMi) dwarf spheroidal (dSph) galaxies based on published and new Keck HIRES spectra for stars in the outer UMi field. Washington+DDO51 filter photometric catalogs provide additional leverage on membership of individual stars, and beyond 0.5 King limiting radii (R_lim) identify bona fide dSph members up to 4.5 times more efficiently than simple color-magnitude diagram selections. Previously reported ``cold populations'' R_lim are not obvious in the data and appear only with particular binning; more or less constant and platykurtic dispersion profiles are characteristic of these dSphs to large radii. We report discovery of UMi stars to at least 2.7 R_lim (i.e.,210 arcmin or 4 kpc). Even with conservative assumptions, a UMi mass of M > 4.9 x 10^8 M_(sun) is required to bind these stars, implying an unlikely global mass-to-light ratio of M/L > 900 (M/L)_(sun). We conclude that we have found stars tidally stripped from UMi.

We report the discovery of three new star clusters in the halo of the Local Group dwarf irregular galaxy NGC 6822. These clusters were found in the deep images taken with the MegaPrime at the CFHT, covering a total field of 2 deg x 2 deg. The most remote cluster is found to be located as far as 79 arcmin away from the center of NGC 6822. This distance is several times larger than the size of the region in NGC 6822 where star clusters were previously found. Morphological structures of the clusters and color-magnitude diagrams of the resolved stars in the clusters show that at least two of these clusters are proabably old globular clusters.

We have obtained HST/STIS far- and near-UV photometry of globular clusters in four fields in the gE galaxy M87. To a limit of m(FUV) = 25 we detect a total of 66 globular clusters (GCs) in common with the deep HST optical-band study of Kundu et al. (1999). Despite strong overlap in V- and I-band properties, the M87 GCs have UV/optical properties that are distinct from clusters in the Milky Way and in M31. M87 clusters, especially metal-poor ones, produce larger hot HB populations than do Milky Way analogues. Cluster mass is probably not a factor in these distinctions. The most metal-rich M87 GCs in our sample are near Z_sun and overlap the local E galaxy sample in estimated Mg_2 line indices. Nonetheless, the clusters produce much more UV light at a given Mg_2, being up to 1 mag bluer than any gE galaxy in (FUV-V) color. The M87 GCs do not appear to represent a transition between Milky Way-type clusters and E galaxies. The differences are in the correct sense if the clusters are significantly older than the E galaxies. Comparisons with Galactic open clusters indicate that the hot stars lie on the extreme horizontal branch, rather than being blue stragglers, and that the EHB becomes well populated for ages > 5 Gyr. We find that 43 of our UV detections have no optical-band counterparts. Most appear to be UV-bright background galaxies, seen through M87. Eleven NUV variable sources detected at only one epoch in the central field are probably classical novae. [Abridged]

We report the results of a photometric search for giant stars associated with the cores of four high velocity clouds (HVCs) -- two of which are compact HVCs -- using the Las Campanas Du Pont 2.5 meter and Cerro Tololo Blanco 4 meter telescopes in combination with a system of filters (Washington M, T_2 + DDO51) useful for identifying low surface gravity, evolved stars. Identical observations of nearby control fields provide a measure of the ``giant star'' background. Our data reach M_0=22 for three of the HVCs and M_0=21.25 for the fourth, depths that allow the detection of any giant stars within 600 kpc. Although we identify a number of faint late-type giant star candidates, we find neither a coherent red giant branch structure nor a clear excess of giant candidate counts in any HVC. This indicates that the giant candidates are probably not related to the HVCs and are more likely to be either random Milky Way giant stars or one of several classes of potential survey contaminants. Echelle spectroscopy of the brightest giant candidates in one HVC and its control field reveal radial velocities representative of the canonical Galactic stellar populations. In addition to these null results, no evidence of any young HVC stellar populations -- represented by blue main sequence stars -- is found, a result consistent with previous searches. Our methodology, specifically designed to find faint diffuse stellar populations, places the tightest upper limit yet on the total stellar mass of HVCs of a few 10^5 M_{\sun}.

We are exploring the extended stellar distributions of Galactic satellite galaxies and globular clusters. For seven objects studied thus far, the observed profile departs from a King function at large r, revealing a ``break population'' of stars. In our sample, the relative density of the ``break'' correlates to the inferred M/L of these objects. We discuss opposing hypotheses for this trend: (1) Higher M/L objects harbor more extended dark matter halos that support secondary, bound, stellar ``halos''. (2) The extended populations around dwarf spheroidals (and some clusters) consist of unbound, extratidal debris from their parent objects, which are undergoing various degrees of tidal disruption. In this scenario, higher M/L ratios reflect higher degrees of virial non-equilibrium in the parent objects, thus invalidating a precept underlying the use of core radial velocities to obtain masses.