With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: i.) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); ii.) Magellanic satellites that were recently captured by the MCs $<$ 1 Gyr ago (Reticulum 2, Phoenix 2); iii.) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and iv.) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC's gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic satellites.

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

Understanding the evolutionary history of the Magellanic Clouds requires an in-depth exploration and characterization of the stellar content in their outer regions, which ultimately are key to tracing the epochs and nature of past interactions. We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at ~1.9 degr from the center of the SMC in the north-east direction.The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2. The shell is composed of a young population with an age ~ 150 Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute g-band magnitude of the shell is M_{g,shell} = -10.78+/- 0.02, with a surface brightness of mu_{g,shell} = 25.81+/- 0.01 mag/arcsec^{2}. We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, ~150 Myr ago.

The Hubble Space Telescope has produced astonishing science over the past thirty years. Hubble's productivity can continue to soar for years to come provided some worn out components get upgraded. While powerful new ground-based and space telescopes are expected to come online over the next decade, none of them will have the UV capabilities that make Hubble a unique observatory. Without Hubble, progress in UV and blue optical astrophysics will be halted. Observations at these wavelengths are key for a range of unresolved astrophysics questions, ranging from the characterization of solar system planets to understanding interaction of galaxies with the intergalactic medium and the formation history of the universe. Hubble will remain our only source of high-angular resolution UV imaging and high-sensitivity UV spectroscopy for the next two decades, offering the ability for continued unique science and maximizing the science return from complementary observatories. Therefore, we recommend that NASA, ESA, and the private sector study the scientific merit, technical feasibility, and risk of a new servicing mission to Hubble to boost its orbit, fix aging components, and expand its instrumentation. Doing so would: 1) keep Hubble on its path to reach its unmet full potential, 2) extend the mission's lifetime past the next decade, which will maximize the synergy of Hubble with other upcoming facilities, and 3) enable and enhance the continuation of scientific discoveries in UV and optical astrophysics.

Discoveries of low mass galaxy pairs and groups are increasing. Studies indicate that dwarf galaxy pairs are gas rich in the field and exhibit elevated star formation rates, suggestive of interactions. Lacking are dynamical models of observed dwarf galaxy pairs to disentangle the physical processes regulating their baryon cycles. We present new optical data and the first detailed theoretical model of an observed tidal encounter between two isolated low mass galaxies, NGC 4490 & NGC 4485. This system is an isolated analog of the Magellanic Clouds and is surrounded by a ~50 kpc extended HI envelope. We use hybrid $N$-body and test-particle simulations along with a visualization interface $Identikit$ to simultaneously reproduce the observed present-day morphology and kinematics. Our results demonstrate how repeated encounters between two dwarf galaxies can "park" baryons at very large distances, without the aid of environmental effects. Our best match to the data is an 8:1 mass ratio encounter where a one-armed spiral is induced in the NGC 4490-analog, which we postulate explains the nature of diffuse starlight presented in the new optical data. We predict that the pair will fully merge in ~370 Myr, but that the extended tidal features will continue to evolve and return to the merged remnant over ~5 Gyr. This pre-processing of baryons will affect the efficiency of gas stripping if such dwarf pairs are accreted by a massive host. In contrast, in isolated environments this study demonstrates how dwarf-dwarf interactions can create a long-lived supply of gas to the merger remnant.

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.

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 recommend configuring the LSST coverage of the inner Galactic plane to allow the production of the definitive age/metallicity map of the Galactic bulge from LSST data, matched to external surveys where appropriate. This will allow the formation history of the Galactic bulge to be reconstructed, as well as furnishing a huge legacy dataset to support one of the key LSST science goals (Mapping the Milky Way). We recommend precise multi-color $grizy$~photometry as deeply as spatial crowding will allow (ideally completing early in the project), complemented by a single-filter survey spread over the entire 10-year time baseline. Both strands should cover as broad an area within the Bulge as possible, with image quality sufficient to reach at least the bulge main sequence turn-off in seeing-limited observations. We specify metrics and a figure of merit by which candidate observing strategies could be evaluated with respect to `static' bulge science (proper motions and photometry)

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.

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 use GAIA DR2 proper motions of the RIOTS4 field OB stars in the Small Magellanic Cloud (SMC) to study the kinematics of runaway stars. The data reveal that the SMC Wing has a systemic peculiar motion relative to the SMC Bar of (v_RA, v_Dec) = (62 +/-7, -18+/-5) km/s and relative radial velocity +4.5 +/- 5.0 km/s. This unambiguously demonstrates that these two regions are kinematically distinct: the Wing is moving away from the Bar, and towards the Large Magellanic Cloud with a 3-D velocity of 64 +/- 10 km/s. This is consistent with models for a recent, direct collision between the Clouds. We present transverse velocity distributions for our field OB stars, confirming that unbound runaways comprise on the order of half our sample, possibly more. Using eclipsing binaries and double-lined spectroscopic binaries as tracers of dynamically ejected runaways, and high-mass X-ray binaries (HMXBs) as tracers of runaways accelerated by supernova kicks, we find significant contributions from both populations. The data suggest that HMXBs have lower velocity dispersion relative to dynamically ejected binaries, consistent with the former corresponding to less energetic supernova kicks that failed to unbind the components. Evidence suggests that our fast runaways are dominated by dynamical, rather than supernova, ejections.

A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectroscopically observed in the literature. We derive systemic PMs for 39 galaxies and galaxy candidates out to 420 kpc, and generally find good consistency for the subset with measurements available from other studies. We derive the implied Galactocentric velocities, and calculate orbits in canonical MW halo potentials of "low" ($0.8 \times 10^{12} M_\odot$) and "high" mass ($1.6 \times 10^{12} M_\odot$). Comparison of the distributions of orbital apocenters and 3D velocities to the halo virial radius and escape velocity, respectively, suggests that the satellite kinematics are best explained in the high-mass halo. Tuc III, Crater II, and additional candidates have orbital pericenters small enough to imply significant tidal influences. Relevant to the missing satellite problem, the fact that fewer galaxies are observed to be near apocenter than near pericenter implies that there must be a population of distant dwarf galaxies yet to be discovered. Of the 39 dwarfs: 12 have orbital poles that do not align with the MW plane of satellites (given reasonable assumptions about its intrinsic thickness); 10 have insufficient PM accuracy to establish whether they align; and 17 satellites align, of which 11 are co-orbiting and (somewhat surprisingly, in view of prior knowledge) 6 are counter-orbiting. Group infall might have contributed to this, but no definitive association is found for the members of the Crater-Leo group.

We present a new distance estimation for the Milky Way dwarf spheroidal satellite Sculptor obtained from multi-epoch mid-infrared observations of RR Lyrae stars. The 3.6 {\mu}m observations have been acquired with the Infrared Array Camera on board the Spitzer Space Telescope as part of the SMHASH Program. Mid-infrared light curves for 42 RRL were obtained, from which we measured Sculptor's distance modulus to be {\mu} = 19.60 $\pm$ 0.02 (statistical) $\pm$ 0.04 (photometric) mag (with $\sigma_{sys}=$ = 0.09 mag), using the 3.6 {\mu}m empirical period-luminosity relations derived from the Galactic globular cluster M4, or {\mu} = 19.57 $\pm$ 0.02 (statistical) $\pm$ 0.04 (photometric) mag (with $\sigma_{sys}=$ = 0.11 mag) using empirical relations in the same passband recently derived from the Large Magellanic Cloud globular cluster Reticulum. Both these measurements are in good agreement with values presented in previous works with Sculptor RR Lyrae stars in optical bands, and are also consistent with recent near-infrared RR Lyrae results. Best agreement with the literature is found for the latter modulus which is equivalent to a distance of d = 82 $\pm$ 1 (statistical) $\pm$ 2 (photometric) kpc (with $\sigma_{sys}=$ = 4 kpc). Finally, using a subsample of RR Lyrae stars with spectroscopic metallicities, we demonstrate that these distance estimates are not affected by metallicity effects.

We quantify the frequency of companions of low redshift ($0.013 < z < 0.0252$), dwarf galaxies ($2 \times 10^8$ M$_\odot <$ M$_{*} < 5 \times 10^9$ M$_\odot$) that are isolated from more massive galaxies in SDSS and compare against cosmological expectations using mock observations of the Illustris simulation. Dwarf multiples are defined as 2 or more dwarfs that have angular separations > 55'', projected separations r$_p < 150$ kpc and relative line-of-sight velocities $\Delta V_{\rm LOS} < 150$ km/s. While the mock catalogs predict a factor of 2 more isolated dwarfs than observed in SDSS, the mean number of observed companions per dwarf is $N_c \sim 0.04$, in good agreement with Illustris when accounting for SDSS sensitivity limits. Removing these limits in the mock catalogs predicts $N_c\sim 0.06$ for future surveys (LSST, DESI), which will be complete to M$_* = 2\times 10^8$ M$_\odot$. The 3D separations of mock dwarf multiples reveal a contamination fraction of $\sim$40% in observations from projection effects. Most isolated multiples are pairs; triples are rare and it is cosmologically improbable that bound groups of dwarfs with more than 3 members exist within the parameter range probed in this study. We find that $<$1% of LMC-analogs in the field have an SMC-analog companion. The fraction of dwarf "Major Pairs'' (stellar mass ratio $>$1:4) steadily increases with decreasing Primary stellar mass, whereas the cosmological "Major Merger rate'' (per Gyr) has the opposite behaviour. We conclude that cosmological simulations can be reliably used to constrain the fraction of dwarf mergers across cosmic time.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We present the first proper motion measurement for an ultra-faint dwarf spheroidal galaxy, Segue 1, using SDSS and LBC data as the first and second epochs separated by a baseline of $\sim 10$ years. We obtain a motion of $\mu_{\alpha}\,\cos(\delta) = -0.37\pm0.57$ mas yr$^{-1}$ and $\mu_{\delta} =-3.39\pm0.58$ mas yr$^{-1}$. Combining this with the known line-of-sight velocity, this corresponds to a Galactocentric V$_\mathrm{rad}=84\pm9$ and V$_\mathrm{tan}=164^{+66}_{-55}$ km s$^{-1}$. Applying Milky Way halo masses between 0.8 to 1.6$\times 10^{12}$ M$_\odot$ results in an apocenter at 33.9$^{+21.7}_{-7.4}$ kpc and pericenter at 15.4$^{+10.1}_{-9.0}$ kpc from the Galactic center, indicating Segue~1 is rather tightly bound to the Milky Way. Since neither the orbital pole of Segue 1 nor its distance to the Milky Way is similar to the more massive classical dwarfs, it is very unlikely that Segue 1 was once a satellite of a massive known galaxy. Using cosmological zoom-in simulations of Milky Way-mass galaxies, we identify subhalos on similar orbits as Segue~1, which imply the following orbital properties: a median first infall 8.1$^{+3.6}_{-4.3}$ Gyrs ago, a median of 4 pericentric passages since then and a pericenter of 22.8$^{+4.7}_{-4.8}$ kpc. This is slightly larger than the pericenter derived directly from Segue 1 and Milky Way parameters, because galaxies with a small pericenter are more likely to be destroyed. Of the surviving subhalo analogs only 27\% were previously a satellite of a more massive dwarf galaxy (that is now destroyed), thus Segue 1 is more likely to have been accreted on its own.

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.

Interactions between pairs of isolated dwarf galaxies provide a critical window into low-mass hierarchical, gas-dominated galaxy assembly and the buildup of stellar mass in low-metallicity systems. We present the first VLT/MUSE optical IFU observations of the interacting dwarf pair dm1647+21, selected from the TiNy Titans survey. The H$\alpha$ emission is widespread and corresponds to a total unobscured star formation rate (SFR) of 0.44 M$_{\odot}$ yr$^{-1}$, 2.7 times higher than the SFR inferred from SDSS data. The implied specific SFR (sSFR) for the system is elevated by more than an order of magnitude above non-interacting dwarfs in the same mass range. This increase is dominated by the lower-mass galaxy, which has a sSFR enhancement of $>$ 50. Examining the spatially-resolved maps of classic optical line diagnostics, we find the ISM excitation can be fully explained by star formation. The velocity field of the ionized gas is not consistent with simple rotation. Dynamical simulations indicate that the irregular velocity field and the stellar structure is consistent with the identification of this system as an ongoing interaction between two dwarf galaxies. The widespread, clumpy enhancements in star formation in this system point to important differences in the effect of mergers on dwarf galaxies, compared to massive galaxies: rather than the funneling of gas to the nucleus and giving rise to a nuclear starburst, starbursts in low-mass galaxy mergers may be triggered by large-scale ISM compression, and thus be more distributed.

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.

We constrain the shape of the Milky Way's halo by dynamical modeling of the observed phase-space tracks of the Pal 5 and GD-1 tidal streams. We find that the only information about the potential gleaned from the tracks of these streams are precise measurements of the shape of the gravitational potential---the ratio of vertical to radial acceleration---at the location of the streams, with weaker constraints on the radial and vertical accelerations separately. The latter will improve significantly with precise proper-motion measurements from Gaia. We measure that the overall potential flattening is 0.95 +/- 0.04 at the location of GD-1 ([R,z] ~ [12.5,6.7] kpc) and 0.94 +/- 0.05 at the position of Pal 5 ([R,z] ~ [8.4,16.8] kpc). Combined with constraints on the force field near the Galactic disk, we determine that the axis ratio of the dark-matter halo's density distribution is 1.05 +/- 0.14 within the inner 20 kpc, with a hint that the halo becomes more flattened near the edge of this volume. The halo mass within 20 kpc is 1.1 +/- 0.1 x 10^{11} M_sun. A dark-matter halo this close to spherical is in tension with the predictions from numerical simulations of the formation of dark-matter halos.

The demographics of dwarf galaxy populations have long been in tension with predictions from the Cold Dark Matter (CDM) paradigm. If primordial density fluctuations were scale-free as predicted, dwarf galaxies should themselves host dark matter subhaloes, the most massive of which may have undergone star formation resulting in dwarf galaxy groups. Ensembles of dwarf galaxies are observed as satellites of more massive galaxies, and there is observational and theoretical evidence to suggest that these satellites at z=0 were captured by the massive host halo as a group. However, the evolution of dwarf galaxies is highly susceptible to environment making these satellite groups imperfect probes of CDM in the low mass regime. We have identified one of the clearest examples to date of hierarchical structure formation at low masses: seven isolated, spectroscopically confirmed groups with only dwarf galaxies as members. Each group hosts 3-5 known members, has a baryonic mass of ~4.4 x 10^9 to 2 x 10^10 Msun, and requires a mass-to-light ratio of <100 to be gravitationally bound. Such groups are predicted to be rare theoretically and found to be rare observationally at the current epoch and thus provide a unique window into the possible formation mechanism of more massive, isolated galaxies.

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.

Palomar 5 (Pal 5) is a faint halo globular cluster associated with narrow tidal tails. It is a useful system to understand the process of tidal dissolution, as well as to constrain the potential of the Milky Way. A well-determined orbit for Pal 5 would enable detailed study of these open questions. We present here the first CCD-based proper motion measurement of Pal 5 obtained using SDSS as a first epoch and new LBT/LBC images as a second, giving a baseline of 15 years. We perform relative astrometry, using SDSS as a distortion-free reference, and images of the cluster and also of the Pal 5 stream for the derivation of the distortion correction for LBC. The reference frame is made up of background galaxies. We correct for differential chromatic refraction using relations obtained from SDSS colors as well as from flux-calibrated spectra, finding that the correction relations for stars and for galaxies are different. We obtain mu_alpha=-2.296+/-0.186 mas/yr and mu_delta=-2.257+/-0.181 mas/yr for the proper motion of Pal 5. We use this motion, and the publicly available code galpy, to model the disruption of Pal 5 in different Milky Way models consisting of a bulge, a disk and a spherical dark matter halo. Our fits to the observed stream properties (streak and radial velocity gradient) result in a preference for a relatively large Pal 5 distance of around 24 kpc. A slightly larger absolute proper motion than what we measure also results in better matches but the best solutions need a change in distance. We find that a spherical Milky Way model, with V_0=220 km/s and V_(20 kpc), i.e., approximately at the apocenter of Pal 5, of 218 km/s, can match the data well, at least for our choice of disk and bulge parametrization.

The hierarchical nature of LCDM suggests that the Magellanic Clouds must have been surrounded by a number of satellites before their infall into the Milky Way. Many of those satellites should still be in close proximity to the Clouds, but some could have dispersed ahead/behind the Clouds along their Galactic orbit. Either way, prior association with the Clouds results in strong restrictions on the present-day positions and velocities of candidate Magellanic satellites: they must lie close to the nearly-polar orbital plane of the Magellanic stream, and their distances and radial velocities must follow the latitude dependence expected for a tidal stream with the Clouds at pericenter. We use a cosmological numerical simulation of the disruption of a massive subhalo in a Milky Way-sized LCDM halo to test whether any of the 20 dwarfs recently-discovered in the DES, SMASH, Pan-STARRS, and ATLAS surveys are truly associated with the Clouds. Of the 6 systems with kinematic data, only Hydra II and Hor 1 have distances and radial velocities consistent with a Magellanic origin. Of the remaining dwarfs, six (Hor 2, Eri 3, Ret 3, Tuc 4, Tuc 5, and Phx 2) have positions and distances consistent with a Magellanic origin, but kinematic data are needed to substantiate that possibility. Conclusive evidence for association would require proper motions to constrain the orbital angular momentum direction, which, for true Magellanic satellites, must coincide with that of the Clouds. We use this result to predict radial velocities and proper motions for all new dwarfs. Our results are relatively insensitive to the assumption of first or second pericenter for the Clouds.

In this paper, we introduce the Local Volume TiNy Titans sample (LV-TNT), which is a part of a larger body of work on interacting dwarf galaxies: TNT (Stierwalt et al. 2015). This LV-TNT sample consists of 10 dwarf galaxy pairs in the Local Universe (< 30 Mpc from Milky Way), which span mass ratios of M_(*,1)/M_(*,2) < 20, projected separations < 100 kpc, and pair member masses of log(M_*/M_Sun) < 9.9. All 10 LV-TNT pairs have resolved synthesis maps of their neutral hydrogen, are located in a range of environments and captured at various interaction stages. This enables us to do a comparative study of the diffuse gas in dwarf-dwarf interactions and disentangle the gas lost due to interactions with halos of massive galaxies, from the gas lost due to mutual interaction between the dwarfs. We find that the neutral gas is extended in the interacting pairs when compared to non-paired analogs, indicating that gas is tidally pre-processed. Additionally, we find that the environment can shape the HI distributions in the form of trailing tails and that the gas is not unbound and lost to the surroundings unless the dwarf pair is residing near a massive galaxy. We conclude that a nearby, massive host galaxy is what ultimately prevents the gas from being reaccreted. Dwarf-dwarf interactions thus represent an important part of the baryon cycle of low mass galaxies, enabling the "parking" of gas at large distances to serve as a continual gas supply channel until accretion by a more massive host.

We present in this study a first analysis of the astrometric error budget of absolute astrometry relative to background galaxies using adaptive optics. We use for this analysis multi-conjugated adaptive optics (MCAO) images obtained with GeMS/GSAOI at Gemini South. We find that it is possible to obtain 0.3 mas reference precision in a random field with 1 hour on source using faint background galaxies. Systematic errors are correctable below that level, such that the overall error is approximately 0.4 mas. Because the reference sources are extended, we find it necessary to correct for the dependency of the PSF centroid on the used aperture size, which would otherwise cause an important bias. This effect needs also to be considered for Extremely Large Telescopes (ELTs). When this effect is corrected, ELTs have the potential to measure proper motions of dwarfs galaxies around M31 with 10 km/s accuracy over a baseline of 5 years.

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 introduce TiNy Titans (TNT), the first systematic study of star formation and the subsequent processing of the interstellar medium in interacting dwarf galaxies. Here we present the first results from a multiwavelength observational program based on a sample of 104 dwarf galaxy pairs selected from a range of environments within the SDSS and caught in various stages of interaction. The TNT dwarf pairs span mass ratios of M1/M2 <10, projected separations <50 kpc, and pair member masses of 7< log(M_*/Msun) <9.7. The dwarf-dwarf merger sequence, as defined by TNT, demonstrates conclusively and for the first time that the star formation enhancement observed for massive galaxy pairs also extends to dwarfs. Star formation is enhanced in paired dwarfs in otherwise isolated environments by a factor of 2.3 (+/- 0.7) at pair separations < 50 kpc relative to unpaired analogs. Starbursts, with Ha EQW >100 A, occur in 20% of TNT dwarf pairs, regardless of environment, compared to only 6-8% of matched unpaired dwarfs. Starbursts can be triggered throughout the merger (out to large pair separations) and not just approaching coalescence. Despite their enhanced star formation, most TNT dwarf pairs have similar gas fractions relative to unpaired dwarfs. Thus, there may be significant reservoirs of diffuse, non-starforming gas surrounding the dwarf pairs or the gas consumption timescales may be long in the starburst phase. The only TNT dwarf pairs with low gas fractions (fgas <0.4) and the only dwarfs with Ha EQW <2 A are found near massive galaxy hosts. Thus, dwarf-dwarf interactions are significant drivers of dwarf galaxy evolution, but ultimately environment is responsible for the quenching of star formation. This work is a precursor to an ongoing high resolution HI and optical imaging program to constrain the spatial distribution of star formation and gas throughout the dwarf-dwarf merger sequence.

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 present the first detailed assessment of the large-scale rotation of any galaxy based on full three-dimensional velocity measurements. We do this for the LMC by combining our HST average proper motion (PM) measurements for stars in 22 fields, with existing line-of-sight (LOS) velocity measurements for 6790 individual stars. We interpret these data with a model of circular rotation in a flat disk. The PM and LOS data paint a consistent picture of the LMC rotation and their combination yields several new insights. The PM data imply a stellar dynamical center that coincides with the HI dynamical center, and a rotation curve amplitude consistent with that inferred from LOS velocity studies. The implied disk viewing angles agree with the range of values found in the literature, but continue to indicate variations with stellar population and/or radius. Young (RSG) stars rotate faster than old (RGB/AGB) stars due to asymmetric drift. Outside the central region, the circular velocity is approximately flat at Vcirc = 91.7 +/- 18.8 km/s. This is consistent with the baryonic Tully-Fisher relation, and implies an enclosed mass M(8.7 kpc) = (1.7 +/- 0.7) x 10^10 solar masses. The virial mass is larger and depends on the full extent of the dark halo. The tidal radius is 22.3 +/- 5.2 kpc (24.0 +/- 5.6 degrees). Combination of the PM and LOS data yields kinematic distance estimates for the LMC, but these are not yet competitive with other methods.

We present proper motions for the Large & Small Magellanic Clouds (LMC & SMC) based on three epochs of \textit{Hubble Space Telescope} data, spanning a $\sim 7$ yr baseline, and centered on fields with background QSOs. The first two epochs, the subject of past analyses, were obtained with ACS/HRC, and have been reanalyzed here. The new third epoch with WFC3/UVIS increases the time baseline and provides better control of systematics. The three-epoch data yield proper motion random errors of only 1-2% per field. For the LMC this is sufficient to constrain the internal proper motion dynamics, as will be discussed in a separate paper. Here we focus on the implied center-of-mass proper motions: mu_W(LMC) = -1.910 +/- 0.020 mas/yr, mu_N(LMC) = 0.229 +/- 0.047 mas/yr, and mu_W(SMC) = -0.772 +/- 0.063 mas/yr, mu_N(SMC) = -1.117 +/- 0.061 mas/yr. We combine the results with a revised understanding of the solar motion in the Milky Way to derive Galactocentric velocities: v_{tot,LMC} = 321 +/- 24 km/s and v_{tot,SMC} = 217 +/- 26 km/s. Our proper motion uncertainties are now dominated by limitations in our understanding of the internal kinematics and geometry of the Clouds, and our velocity uncertainties are dominated by distance errors. Orbit calculations for the Clouds around the Milky Way allow a range of orbital periods, depending on the uncertain masses of the Milky Way and LMC. Periods $\lesssim 4$ Gyr are ruled out, which poses a challenge for traditional Magellanic Stream models. First-infall orbits are preferred (as supported by other arguments as well) if one imposes the requirement that the LMC and SMC must have been a bound pair for at least several Gyr.

We present evidence for a new Milky Way stellar tidal stream in the direction of the Andromeda and Triangulum (M31 and M33) galaxies. Using a matched-filter technique, we search the Sloan Digital Sky Survey Data Release 8 by creating stellar density maps which probe the Milky Way halo at distances between 8 and 40 kpc. A visual search of these maps recovers all of the major known stellar streams, as well as a new stream in the direction of M31/M33 which we name the Triangulum stream. The stream spans 0.2 deg by 12 deg on the sky, or 75 pc by 5.5 kpc in physical units with a best fitting distance of 26+/-4 kpc. The width of the stream is consistent with being the tidal remnant of a globular cluster. A color magnitude diagram of the stream region shows an overdensity which, if identified as a main sequence turn-off, corresponds to an old (~12 Gyr) and metal-poor ([Fe/H] \sim -1.0 dex) stellar population. Future kinematic studies of this and similar cold streams will provide tight constraints on the shape of the Galactic gravitational potential.

We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high precision proper motions (Kallivayalil et al. 2006a,b). Given the new velocities, cosmological simulations of structure formation favor a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy (Boylan-Kolchin et al. 2011, Busha et al. 2011). We illustrate here that the observed irregular morphology and internal kinematics of the MCs (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. This picture further supports a first infall scenario (Besla et a. 2007). In particular, we demonstrate that the Magellanic Stream, a band of HI gas trailing behind the MCs 150 degrees across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally explained by a recent direct collision with the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars (de Vaucouleurs & Freeman 1972), the majority of which are not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are important not only for the evolution of dwarf galaxies but also have direct consequences for the buildup of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type halos.

Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW (period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ~150 degrees across the sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW halo gas. We propose an alternative formation mechanism in which material is removed by LMC tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility and generality of this scenario using an N-body/SPH simulation with cosmologically motivated initial conditions constrained by the observations. Under these conditions we demonstrate that it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or interacting in isolation.

We present MMT/Megacam imaging of the Leo~IV dwarf galaxy in order to investigate its structure and star formation history, and to search for signs of association with the recently discovered Leo~V satellite. Based on parameterized fits, we find that Leo~IV is round, with $\epsilon < 0.23$ (at the 68\% confidence limit) and a half-light radius of $r_{h} \simeq 130$ pc. Additionally, we perform a thorough search for extended structures in the plane of the sky and along the line of sight. We derive our surface brightness detection limit by implanting fake structures into our catalog with stellar populations identical to that of Leo~IV. We show that we are sensitive to stream-like structures with surface brightness $\mu_{r}\lesssim29.6$ mag arcsec$^{-2}$, and at this limit, we find no stellar bridge between Leo IV (out to a radius of $\sim$0.5 kpc) and the recently discovered, nearby satellite Leo V. Using the color magnitude fitting package StarFISH, we determine that Leo~IV is consistent with a single age ($\sim$14 Gyr), single metallicity ($[Fe/H]\sim-2.3$) stellar population, although we can not rule out a significant spread in these value. We derive a luminosity of $M_{V}=-5.5\pm0.3$. Studying both the spatial distribution and frequency of Leo~IV's 'blue plume' stars reveals evidence for a young ($\sim$2 Gyr) stellar population which makes up $\sim$2\% of its stellar mass. This sprinkling of star formation, only detectable in this deep study, highlights the need for further imaging of the new Milky Way satellites along with theoretical work on the expected, detailed properties of these possible 'reionization fossils'.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

We study the dynamics of the Magellanic Clouds in a model for the Local Group whose mass is constrained using the timing argument/two-body limit of the action principle. The goal is to evaluate the role of M31 in generating the high angular momentum orbit of the Clouds, a puzzle that has only been exacerbated by the latest $HST$ proper motion measurements. We study the effects of varying the total Local Group mass, the relative mass of the Milky Way and M31, the proper motion of M31, and the proper motion of the LMC on this problem. Over a large part of this parameter-space we find that tides from M31 are insignificant. For a range of LMC proper motions approximately $3\sigma$ higher than the mean and total Local Group mass $> 3.5\times 10^{12} M_\odot$, M31 can provide a significant torque to the LMC orbit. However, if the LMC is bound to the MW, then M31 is found to have negligible effect on its motion and the origin of the high angular momentum of the system remains a puzzle. Finally, we use the timing argument to calculate the total mass of the MW-LMC system based on the assumption that they are encountering each other for the first time, their previous perigalacticon being a Hubble time ago, obtaining $M_{\rm MW} + M_{\rm LMC} = (8.7 \pm 0.8) \times 10^{11} M_\odot$.

In HST Cycles 11 and 13 we obtained two epochs of ACS/HRC data for fields in the Magellanic Clouds centered on background quasars. We used these data to determine the proper motions of the LMC and SMC to better than 5% and 15% respectively. The results had a number of unexpected implications for the Milky Way-LMC-SMC system. The implied three-dimensional velocities were larger than previously believed and close to the escape velocity in a standard 10^12 solar mass Milky Way dark halo, implying that the Clouds may be on their first passage. Also, the relative velocity between the LMC and SMC was larger than expected, leaving open the possibility that the Clouds may not be bound to each other. To further verify and refine our results we requested an additional epoch of data in Cycle 16 which is being executed with WFPC2/PC due to the failure of ACS. We present the results of an ongoing analysis of these WFPC2 data which indicate good consistency with the two-epoch results.

We review our understanding of the kinematics of the LMC and the SMC, and their orbit around the Milky Way. The line-of-sight velocity fields of both the LMC and SMC have been mapped with high accuracy using thousands of discrete traces, as well as HI gas. The LMC is a rotating disk for which the viewing angles have been well-established using various methods. The disk is elliptical in its disk plane. The disk thickness varies depending on the tracer population, with V/sigma ranging from 2-10 from the oldest to the youngest population. For the SMC, the old stellar population resides in a spheroidal distribution with considerable line-of-sight depth and low V/sigma. Young stars and HI gas reside in a more irregular rotating disk. Mass estimates based on the kinematics indicate that each Cloud is embedded in a dark halo. Proper motion measurements with HST show that both galaxies move significantly more rapidly around the Milky Way than previously believed. This indicates that for a canonical 10^12 solar mass Milky Way the Clouds are only passing by us for the first time. Although a higher Milky Way mass yields a bound orbit, this orbit is still very different from what has been previously assumed in models of the Magellanic Stream. Hence, much of our understanding of the history of the Magellanic System and the formation of the Magellanic Stream may need to be revised. The accuracy of the proper motion data is insufficient to say whether or not the LMC and SMC are bound to each other, but bound orbits do exist within the proper motion error ellipse.

The HST proper motion (PM) measurements of the Clouds have severe implications for their interaction history with the Milky Way (MW) and with each other. The Clouds are likely on their first passage about the MW and the SMC's orbit about the LMC is better described as quasi-periodic rather than circular. Binary L/SMC orbits that satisfy observational constraints on their mutual interaction history (e.g. the formation of the Magellanic Bridge during a collision between the Clouds ~300 Myr ago) can be located within 1 sigma of the mean PMs. However, these binary orbits are not co-located with the Magellanic Stream (MS) when projected on the plane of the sky and the line-of-sight velocity gradient along the LMC's orbit is significantly steeper than that along the MS. These combined results ultimately rule out a purely tidal origin for the MS: tides are ineffective without multiple pericentric passages and can neither decrease the velocity gradient nor explain the offset stream in a polar orbit configuration. Alternatively, ram pressure stripping of an extended gaseous disk may naturally explain the deviation. The offset also suggests that observations of the little-explored region between RA 21h and 23h are crucial for characterizing the full extent of the MS.

Recent proper motion measurements of the Large and Small Magellanic Clouds (LMC and SMC, respectively) by Kallivayalil et al (2006a,b) suggest that the 3D velocities of the Clouds are substantially higher (~100 km/s) than previously estimated and now approach the escape velocity of the Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately modeled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds using the new velocities and a LCDM-motivated MW model with virial mass Mvir = 1e12 Msun (e.g. Klypin et al 2002). We conclude that the Clouds are either currently on their first passage about the MW or, if the MW can be accurately modeled by an isothermal sphere to distances >200 kpc (i.e., Mvir > 2e12 Msun), that their orbital period and apogalacticon distance must be a factor of two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss the implications of this orbital analysis for our understanding of the star formation history, the nature of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing the LMC and SMC that extends roughly 100 degrees across the sky. Specifically, as a consequence of the new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and ram pressure stripping models.

We report on the identification of the lens responsible for microlensing event MACHO-LMC-20. As part of a \textit{Spitzer}/IRAC program conducting mid-infrared follow-up of the MACHO Large Magellanic Cloud microlensing fields, we discovered a significant flux excess at the position of the source star for this event. These data, in combination with high resolution near-infrared \textit{Magellan}/PANIC data has allowed us to classify the lens as an early M dwarf in the thick disk of the Milky Way, at a distance of $\sim 2$ kpc. This is only the second microlens to have been identified, the first also being a M dwarf star in the disk. Together, these two events are still consistent with the expected frequency of nearby stars in the Milky Way thin and thick disks acting as lenses.

We present a measurement of the systemic proper motion of the Small Magellanic Cloud (SMC) made using the Advanced Camera for Surveys (ACS) on the \textit{Hubble Space Telescope} (\textit{HST}). We tracked the SMC's motion relative to 4 background QSOs over a baseline of approximately 2 years. The measured proper motion is : $\mu_W = -1.16 \pm 0.18 \masyr, \mu_N = -1.17 \pm 0.18 \masyr$. This is the best measurement yet of the SMC's proper motion. We combine this new result with our prior estimate of the proper motion of the Large Magellanic Cloud (LMC) from the same observing program to investigate the orbital evolution of both Clouds over the past 9 Gyr. The current relative velocity between the Clouds is $105 \pm 42 \kms$. Our investigations of the past orbital motions of the Clouds in a simple model for the dark halo of the Milky Way imply that the Clouds could be unbound from each other. However, our data are also consistent with orbits in which the Clouds have been bound to each other for approximately a Hubble time. Smaller proper motion errors and better understanding of the LMC and SMC masses would be required to constrain their past orbital history and their bound vs. unbound nature unambiguously. The new proper motion measurements should be sufficient to allow the construction of improved models for the origin and properties of the Magellanic Stream. In turn, this will provide new constraints on the properties of the Milky Way dark halo.

We present a measurement of the systemic proper motion of the Large Magellanic Cloud (LMC) from astrometry with the High Resolution Camera (HRC) of the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). We observed LMC fields centered on 21 background QSOs that were discovered from their optical variability in the MACHO database. The QSOs are distributed homogeneously behind the central few degrees of the LMC. With 2 epochs of HRC data and a ~2 year baseline we determine the proper motion of the LMC to better than 5% accuracy: mu_W = -2.03 +/- 0.08 mas/yr; mu_N = 0.44 +/- 0.05 mas/yr. This is the most accurate proper motion measurement for any Milky Way satellite thus far. When combined with HI data from the Magellanic Stream this should provide new constraints on both the mass distribution of the Galactic Halo and models of the Stream.

We have carried out photometry of the microlensing event MACHO-LMC-5 with Spitzer's IRAC ten years after the magnification of the LMC source star was recorded. This event is unique in the annals of gravitational microlensing: the lensing star itself has been observed using HST (once with WFPC2 and twice with ACS/HRC). Since the separation between the source and lens at the epoch of the Spitzer observations was $\sim0.24''$, the two stars cannot be resolved in the Spitzer images. However, the IRAC photometry clearly establishes that the lens is a M5 dwarf star from its infrared excess, which in turn yields a mass of $\sim0.2 M_{\odot}$. This demonstrates the potential of Spitzer to detect the lenses in other gravitational microlensing events.