The Large Magellanic Cloud (LMC) is the closest and most studied example of an irregular galaxy. Among its principal defining morphological features, its off-centred bar and single spiral arm stand out, defining a whole family of galaxies known as the Magellanic spirals (Sm). These structures are thought to be triggered by tidal interactions and possibly maintained via gas accretion. However, it is still unknown whether they are long-lived stable structures. In this work, by combining photometry that reaches down to the oldest main sequence turn-off in the colour-magnitude diagrams (CMD, up to a distance of $\sim$4.4 kpc from the LMC centre) from the SMASH survey and CMD fitting techniques, we find compelling evidence supporting the long-term stability of the LMC spiral arm, dating the origin of this structure to more than 2~Gyr ago. The evidence suggests that the close encounter between the LMC and the Small Magellanic Cloud (SMC) that produced the gaseous Magellanic Stream and its Leading Arm (LA) also triggered the formation of the LMC's spiral arm. Given the mass difference between the Clouds and the notable consequences of this interaction, we can speculate that this should have been one of their closest encounters. These results set important constraints on the timing of LMC-SMC collisions, as well as on the physics behind star formation induced by tidal encounters.

The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.

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.

Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution ($m_{dm} = 4 \times 10^4 M_{\odot}$ ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself.

Using N-body simulations of the Large Magellanic Cloud (LMC's) passage through the Milky Way (MW), tailored to reproduce observed kinematic properties of both galaxies, we show that the high-speed tail of the Solar Neighborhood dark matter distribution is overwhelmingly of LMC origin. Two populations contribute at high speeds: 1) Particles that were once bound to the LMC, and 2) MW halo particles that have been accelerated owing to the response of the halo to the recent passage of the LMC. These particles reach speeds of 700-900 km/s with respect to the Earth, above the local escape speed of the MW. The high-speed particles follow trajectories similar to the Solar reflex motion, with peak velocities reached in June. For low-mass dark matter, these high-speed particles can dominate the signal in direct-detection experiments, extending the reach of the experiments to lower mass and elastic scattering cross sections even with existing data sets. Our study shows that even non-disrupted MW satellite galaxies can leave a significant dark-matter footprint in the Solar Neighborhood.

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 current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

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.

We investigate the absorption features associated with a gas-rich dwarf galaxy using cosmological hydrodynamics simulations. Our goal is to explore whether the progenitors of the lowest mass dwarf galaxies known to harbor neutral hydrogen today (M_star~10^6 solar mass, M_halo=4x10^9 solar mass) could possibly be detected as Damped Lyman-alpha Absorbers (DLAs) over cosmic time. We trace the evolution of a single dwarf galaxy, pre-selected to contain DLAs, from the era of the first metal-free, so-called Population~III (Pop~III), stars, down to z=0, thus allowing us to study the metal enrichment history of DLAs associated with the simulated galaxy. We find that the progenitors of the simulated dwarf are expected to be seen for most of their evolution as DLAs that are contaminated by normal, Population~II, stars. The time period during which DLAs are only metal-enriched by Pop~III stars, on the other hand, is likely very brief, confined to high redshifts, z~6. The susceptibility of the dwarfs to the external UV radiation background allows them to preserve neutral gas only at the centre (a few ~100 pc). This results in a small probability that the simulated dwarf would be observed as a DLA. This study suggests that DLAs are unlikely to be hosted in the lowest mass dwarfs that can harbor neutral gas (M_halo~ 4x10^9 solar mass), below which neutral gas is unlikely to exist. However, this study does illustrate that, when detected, absorption lines provide a powerful method for probing ISM conditions inside the smallest dwarf galaxies at intermediate to high redshifts.

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.

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.

We present a study of the three-dimensional (3D) structure of the Large Magellanic Cloud (LMC) using ~2.2 million red clump (RC) stars selected from the Survey of the MAgellanic Stellar History. To correct for line-of-sight dust extinction, the intrinsic RC color and magnitude and their radial dependence are carefully measured by using internal nearly dust-free regions. These are then used to construct an accurate 2D reddening map (165 square degrees with ~10 arcmin resolution) of the LMC disk and the 3D spatial distribution of RC stars. An inclined disk model is fit to the 2D distance map yielding a best-fit inclination angle i = 25.86(+0.73,-1.39) degrees with random errors of +\-0.19 degrees and line-of-nodes position angle theta = 149.23(+6.43,-8.35) degrees with random errors of +/-0.49 degrees. These angles vary with galactic radius, indicating that the LMC disk is warped and twisted likely due to the repeated tidal interactions with the Small Magellanic Cloud (SMC). For the first time, our data reveal a significant warp in the southwestern part of the outer disk starting at rho ~ 7 degrees that departs from the defined LMC plane up to ~4 kpc toward the SMC, suggesting that it originated from a strong interaction with the SMC. In addition, the inner disk encompassing the off-centered bar appears to be tilted up to 5-15 degrees relative to the rest of the LMC disk. These findings on the outer warp and the tilted bar are consistent with the predictions from the Besla et al. simulation of a recent direct collision with the SMC.

We explore the stellar structure of the Large Magellanic Cloud (LMC) disk using data from the Survey of the MAgellanic Stellar History (SMASH) and the Dark Energy Survey. We detect a ring-like stellar overdensity in the red clump star count map at a radius of ~6 degrees (~5.2 kpc at the LMC distance) that is continuous over ~270 degrees in position angle and is only limited by the current data coverage. The overdensity shows an amplitude up to 2.5 times higher than that of the underlying smooth disk. This structure might be related to the multiple arms found by de Vaucouleurs. We find that the overdensity shows spatial correlation with intermediate-age star clusters, but not with young (< 1 Gyr) main-sequence stars, indicating the stellar populations associated with the overdensity are intermediate in age or older. Our findings on the LMC overdensity can be explained by either of two distinct formation mechanisms of a ring-like overdensity: (1) the overdensity formed out of an asymmetric one-armed spiral wrapping around the LMC main body, which is induced by repeated encounters with the Small Magellanic Cloud (SMC) over the last Gyr, or (2) the overdensity formed very recently as a tidal response to a direct collision with the SMC. Although the measured properties of the overdensity alone cannot distinguish between the two candidate scenarios, the consistency with both scenarios suggests that the ring-like overdensity is likely a product of tidal interaction with the SMC, but not with the Milky Way halo.

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

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 the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material.

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.

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 N-body simulations of a Sagittarius-like dwarf spheroidal galaxy (Sgr) that follow its orbit about the Milky Way (MW) since its first crossing of the Galaxy's virial radius to the present day. As Sgr orbits around the MW, it excites vertical oscillations, corrugating and flaring the Galactic stellar disc. These responses can be understood by a two-phase picture in which the interaction is first dominated by torques from the wake excited by Sgr in the MW dark halo before transitioning to tides from Sgr's direct impact on the disc at late times. We show for the first time that a massive Sgr model simultaneously reproduces the locations and motions of arc-like over densities, such as the Monoceros Ring and the Triangulum Andromeda stellar clouds, that have been observed at the extremities of the disc, while also satisfying the solar neighbourhood constraints on the vertical structure and streaming motions of the disc. In additional simulations, we include the Large Magellanic Cloud (LMC) self consistently with Sgr. The LMC introduces coupling through constructive and destructive interference, but no new corrugations. In our models, the excitation of the current structure of the outer disk can be traced to interactions as far back as 6-7 Gyr ago (corresponding to $z\leq1$). Given the apparently quiescent accretion history of the MW over this timescale, this places Sgr as the main culprit behind the vertical oscillations of the disc and the last major accretion event for the Galaxy with the capacity to modulate its chemodynamical structure.

We present the discovery of a very faint stellar system, SMASH 1, that is potentially a satellite of the Large Magellanic Cloud. Found within the Survey of the MAgellanic Stellar History (SMASH), SMASH 1 is a compact ($r_h = 9.1^{+5.9}_{-3.4}$ pc) and very low luminosity (M_V = -1.0 +/- 0.9, $L_V=10^{2.3 +/- 0.4}$ Lsun) stellar system that is revealed by its sparsely populated main sequence and a handful of red-giant-branch candidate member stars. The photometric properties of these stars are compatible with a metal-poor ([Fe/H]=-2.2) and old (13 Gyr) isochrone located at a distance modulus of ~18.8, i.e. a distance of ~57 kpc. Situated at 11.3$^\circ$ from the LMC in projection, its 3-dimensional distance from the Cloud is ~13 kpc, consistent with a connection to the LMC, whose tidal radius is at least 16 kpc. Although the nature of SMASH 1 remains uncertain, its compactness favors it being a stellar cluster and hence dark-matter free. If this is the case, its dynamical tidal radius is only <19 pc at this distance from the LMC, and smaller than the system's extent on the sky. Its low luminosity and apparent high ellipticity ($\epsilon=0.62^{+0.17}_{-0.21}$) with its major axis pointing toward the LMC may well be the tell-tale sign of its imminent tidal demise.

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 investigate the star formation history and chemical evolution of isolated analogues of Local Group (LG) ultra faint dwarf galaxies (UFDs; stellar mass range of 10^2 solar mass < M_star <10^5 solar mass) and gas rich, low mass dwarfs (Leo P analogs; stellar mass range of 10^5 solar mass < M_star <10^6 solar mass). We perform a suite of cosmological hydrodynamic zoom-in simulations to follow their evolution from the era of the first generation of stars down to z=0. We confirm that reionization, combined with supernova (SN) feedback, is primarily responsible for the truncated star formation in UFDs. Specifically, haloes with a virial mass of M_vir < 2 x 10^9 solar mass form> 90\% of stars prior to reionization. Our work further demonstrates the importance of Population~III (Pop~III) stars, with their intrinsically high $\rm [C/Fe]$ yields, and the associated external metal-enrichment, in producing low-metallicity stars ($\rm [Fe/H]\lesssim-4$) and carbon-enhanced metal-poor (CEMP) stars. We find that UFDs are composite systems, assembled from multiple progenitor haloes, some of which hosted only Population~II (Pop~II) stars formed in environments externally enriched by SNe in neighboring haloes, naturally producing, extremely low-metallicity Pop~II stars. We illustrate how the simulated chemical enrichment may be used to constrain the star formation histories (SFHs) of true observed UFDs. We find that Leo P analogs can form in haloes with M_vir ~ 4 x 10^9 solar mass (z=0). Such systems are less affected by reionization and continue to form stars until z=0, causing higher metallicity tails. Finally, we predict the existence of extremely low-metallicity stars in LG UFD galaxies that preserve the pure chemical signatures of Pop~III nucleosynthesis.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

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 report the discovery of a stellar over-density 8$^{\circ}$ north of the center of the Small Magellanic Cloud (Small Magellanic Cloud Northern Over-Density; SMCNOD) using data from the first two years of the Dark Energy Survey (DES) and the first year of the MAGellanic SatelLITEs Survey (MagLiteS). The SMCNOD is indistinguishable in age, metallicity and distance from the nearby SMC stars, being primarly composed of intermediate-age stars (6 Gyr, Z=0.001), with a small fraction of young stars (1 Gyr, Z=0.01). The SMCNOD has an elongated shape with an ellipticity of 0.6 and a size of $\sim$ 6x2 deg. It has an absolute magnitude of $M_V \cong$ -7.7, $r_h = 2.1$ kpc, and $\mu_V(r<r_h)$ = 31.2 mag arcsec$^{-2}$. We estimate a stellar mass of $\sim 10^5$ $M_{\odot}$, following a Kroupa mass function. The SMCNOD was probably removed from the SMC disk by tidal stripping, since it is located near the head of the Magellanic Stream, and the literature indicates likely recent LMC-SMC encounters. This scenario is supported by the lack of significant HI gas. Other potential scenarios for the SMCNOD origin are a transient over-density within the SMC tidal radius or a primordial SMC satellite in advanced stage of disruption.

Space observatories like the Hubble Space Telescope and Gaia are providing unprecedented 6D phase space information of satellite galaxies. Such measurements can shed light on the structure and assembly history of the Local Group, but improved statistical methods are needed to use them efficiently. Here we illustrate such a method using analogs of the Local Group's two most massive satellite galaxies, the Large Magellanic Cloud (LMC) and Triangulum (M33), from the Illustris dark-matter-only cosmological simulation. We use a Bayesian inference scheme combining measurements of positions, velocities, and specific orbital angular momenta ( j ) of the LMC/M33 with importance sampling of their simulated analogs to compute posterior estimates of the Milky Way (MW) and Andromeda's (M31) halo masses. We conclude the resulting host halo mass is more susceptible to bias when using measurements of the current position and velocity of satellites, especially when satellites are at short-lived phases of their orbits (i.e. at pericentre). Instead, the j value of a satellite is well-conserved over time and provides a more reliable constraint on host mass. The inferred virial mass of the MW (M31) using j of the LMC (M33) is $\rm M_{vir, MW}=1.02^{+0.77}_{-0.55}\times10^{12}\; M_{\odot}$ ($\rm M_{vir, M31}=1.37^{+1.39}_{-0.75}\times10^{12}\; M_{\odot}$). Choosing simulated analogs whose j values are consistent with the conventional picture of a previous (< 3 Gyr ago), close encounter (< 100 kpc) of M33 about M31 results in a very low virial mass for M31 ($\rm\sim\!10^{12}\; M_{\odot}$). This supports the new scenario put forth in Patel et al. (2017), wherein M33 is on its first passage about M31 or on a long period orbit. We conclude that this Bayesian inference scheme, utilising satellite j, is a promising method to reduce the current factor of two spread in the mass range of the MW and M31 moving forward.

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.

The Milky Way (MW) and M31 both harbor massive satellite galaxies, the Large Magellanic Cloud (LMC) and M33, which may comprise up to 10 per cent of their host's total mass. Massive satellites can change the orbital barycentre of the host-satellite system by tens of kiloparsecs and are cosmologically expected to harbor dwarf satellite galaxies of their own. Assessing the impact of these effects depends crucially on the orbital histories of the LMC and M33. Here, we revisit the dynamics of the MW-LMC system and present the first detailed analysis of the M31-M33 system utilizing high precision proper motions and statistics from the dark matter-only Illustris cosmological simulation. With the latest Hubble Space Telescope proper motion measurements of M31, we reliably constrain M33's interaction history with its host. In particular, like the LMC, M33 is either on its first passage (t_{inf} < 2 Gyr ago) or if M31 is massive (>=2x10^12 Msun), it is on a long period orbit of about 6 Gyr. Cosmological analogs of the LMC and M33 identified in Illustris support this picture and provide further insight about their host masses. We conclude that, cosmologically, massive satellites like the LMC and M33 are likely completing their first orbits about their hosts. We also find that the orbital energies of such analogs prefer a MW halo mass ~1.5x10^12 Msun and an M31 halo mass >=1.5x10^12 Msun. Despite conventional wisdom, we conclude it is highly improbable that M33 made a close (< 100 kpc) approach to M31 recently ( t_{peri} < 3 Gyr ago). Such orbits are rare (< 1 per cent) within the 4$\sigma$ error space allowed by observations. This conclusion cannot be explained by perturbative effects through four body encounters between the MW, M31, M33, and the LMC. This surprising result implies that we must search for a new explanation for M33's strongly warped gas and stellar discs.

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 N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position and velocity with respect to the Galactic Center. We find that the LMC can produce strong disturbances - warping of the stellar and gaseous discs - in the Galaxy, without violating constraints from the phase-space distribution of stars in the Solar Neighbourhood. All models correctly reproduce the phases of the warp and its anti-symmetrical shape about the disc's mid-plane. If the warp is due to the LMC alone, then the largest mass model is favoured ($2.5\times10^{11}\,\rm{M_{\odot}}$). Still, some quantitative discrepancies remain, including deficits in height of $\Delta Z=0.7 \, \rm{kpc}$ at $R=22 \,\rm{kpc}$ and $\Delta Z=0.7\, \rm{kpc}$ at $R=16 \,\rm{kpc}$. This suggests that even higher infall masses for the LMC's halo are allowed by the data. A comparison with the vertical perturbations induced by a heavy Sagittarius dSph model ($10^{11}\,\rm{M_{\odot}}$) suggest that positive interference with the LMC warp is expected at $R=16 \, \rm{kpc}$. We conclude that the vertical structure of the Galactic disc beyond the Solar Neighbourhood may jointly be shaped by its most massive satellites. As such, the current structure of the MW suggests we are seeing the process of disc heating by satellite interactions in action.

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 report the discovery of one RR Lyrae star in the ultra--faint satellite galaxy Hydra II based on time series photometry in the g, r and i bands obtained with the Dark Energy Camera at Cerro Tololo Interamerican Observatory, Chile. The RR Lyrae star has a mean magnitude of $i = 21.30\pm 0.04$ which translates to a heliocentric distance of $151\pm 8$ kpc for Hydra II; this value is $\sim 13\%$ larger than the estimate from the discovery paper based on the average magnitude of several blue horizontal branch star candidates. The new distance implies a slightly larger half-light radius of $76^{+12}_{-10}$ pc and a brighter absolute magnitude of $M_V = -5.1 \pm 0.3$, which keeps this object within the realm of the dwarf galaxies. The pulsational properties of the RR Lyrae star ($P=0.645$ d, $\Delta g = 0.68$ mag) suggest Hydra II may be a member of the intermediate Oosterhoff or Oosterhoff II group. A comparison with other RR Lyrae stars in ultra--faint systems indicates similar pulsational properties among them, which are different to those found among halo field stars and those in the largest of the Milky Way satellites. We also report the discovery of 31 additional short period variables in the field of view (RR Lyrae, SX Phe, eclipsing binaries, and a likely anomalous cepheid). However, given their magnitudes and large angular separation from Hydra II, they must be field stars not related to Hydra II.

We present deep optical images of the Large and Small Magellanic Clouds (LMC and SMC) using a low cost telephoto lens with a wide field of view to explore stellar substructure in the outskirts of the stellar disk of the LMC (r < 10 degrees from the center). These data have higher resolution than existing star count maps, and highlight the existence of stellar arcs and multiple spiral arms in the northern periphery, with no comparable counterparts in the South. We compare these data to detailed simulations of the LMC disk outskirts, following interactions with its low mass companion, the SMC. We consider interaction in isolation and with the inclusion of the Milky Way tidal field. The simulations are used to assess the origin of the northern structures, including also the low density stellar arc recently identified in the DES data by Mackey et al. 2015 at ~ 15 degrees. We conclude that repeated close interactions with the SMC are primarily responsible for the asymmetric stellar structures seen in the periphery of the LMC. The orientation and density of these arcs can be used to constrain the LMC's interaction history with and impact parameter of the SMC. More generally, we find that such asymmetric structures should be ubiquitous about pairs of dwarfs and can persist for 1-2 Gyr even after the secondary merges entirely with the primary. As such, the lack of a companion around a Magellanic Irregular does not disprove the hypothesis that their asymmetric structures are driven by dwarf-dwarf interactions.

This proceeding overviews our current understanding of the orbital history and mass of the Large and Small Magellanic Clouds. Specifically I will argue that the Clouds are on their first infall about our Milky Way and that their total masses are necessarily ~10 times larger than traditionally estimated. This conclusion is based on the recently revised HST proper motions of the Clouds and arguments concerning the binary status of the LMC-SMC pair and their baryon fractions

This paper explores the effect of the LMC on the mass estimates obtained from the timing argument. We show that accounting for the presence of the LMC systematically lowers the Local Group mass ($M_{\rm LG}$) derived from the relative motion of the Milky Way--Andromeda pair. Motivated by this result we apply a Bayesian technique devised by Pe\~narrubia et al. (2014) to simultaneously fit (i) distances and velocities of galaxies within 3~Mpc and (ii) the relative motion between the Milky Way and Andromeda derived from HST observations, with the LMC mass ($M_{\rm LMC}$) as a free parameter. Our analysis returns a Local Group mass $M_{\rm LG}=2.64^{+0.42}_{-0.38}\times 10^{12}M_\odot$ at a 68\% confidence level. The masses of the Milky Way, $M_{\rm MW}=1.04_{-0.23}^{+0.26}\times 10^{12}M_\odot$, and Andromeda, $M_{\rm M31}=1.33_{-0.33}^{+0.39}\times 10^{12}M_\odot$, are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find a (total) LMC mass $M_{\rm LMC}=0.25_{-0.08}^{+0.09}\times 10^{12}M_\odot$, which is indicative of an extended dark matter halo and supports the scenario where this galaxy is just past its first pericentric approach. Consequently, these results suggest that the LMC may induce significant perturbations on the Galactic potential.

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), implying a well-constrained pericentric passage about the Milky Way (MW) ~ 50 Myr ago. In this scenario, the LMC's gaseous disk has recently experienced maximal ram pressure stripping, suggesting the current extent of its HI disk directly probes the medium in which it is moving. From the observed stellar and HI distributions of the system we find evidence of a truncated gas profile along the windward ``leading edge' of the LMC disk, despite a far more extended stellar component. We explore the implications of this ram pressure stripping signature, using both analytic prescriptions and full 3-dimensional hydrodynamic simulations of the LMC. Our simulations subject the system to a headwind whose velocity components correspond directly to the recent orbital history of the LMC. We vary the density of this headwind, using a variety of sampled parameters for a Beta-profile for a theoretical MW circumgalactic medium (CGM), comparing the resulting HI morphology directly to observations of the LMC HI and stellar components. This model can match the radial extent of the LMC's leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is n(R = 48.2 +/- 5 kpc) = 1.1 (+.44/-.45) x 1e-4 cm^-3. The implied pericentric density proves insensitive to both the broader CGM structure and temperature profile, thus providing a model-independent constraint on the local gas density. This result imposes an important constraint on the density profile of the MW's CGM, and thus the total baryon content of the MW. From our work, assuming a Beta-profile valid to ~ Rvir, we infer a total diffuse CGM mass M(300 kpc) = 2.6 +/- 1.4 x 1e10 Msun or approximately 15% of a 1e12 Msun MW's baryonic mass budget.

We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.

Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC's orbital period and apocentric distance derived assuming a fixed Milky Way are significantly shortened for models where the Milky Way is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the Milky Way center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km/s. Furthermore, we show that the gravitational pull of the LMC and response of the Milky Way are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the Milky Way around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the Milky Way and Local Group.

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

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 introduce a novel theoretical model to explain the long-standing puzzle of the nature of the microlensing events reported towards the Large Magellanic Cloud (LMC) by the MACHO and OGLE collaborations. We propose that a population of tidally stripped stars from the Small Magellanic Cloud (SMC) located ~4-10 kpc behind a lensing population of LMC disk stars can naturally explain the observed event durations (17-71 days), event frequencies and spatial distribution of the reported events. Differences in the event frequencies reported by the OGLE (~0.33 /yr) and MACHO (~1.75 /yr) surveys appear to be naturally accounted for by their different detection efficiencies and sensitivity to faint sources. The presented models of the Magellanic System were constructed without prior consideration of the microlensing implications. These results favor a scenario for the interaction history of the Magellanic Clouds, wherein the Clouds are on their first infall towards the Milky Way and the SMC has recently collided with the LMC 100-300 Myr ago, leading to a large number of faint sources distributed non-uniformly behind the LMC disk. In contrast to self-lensing models, microlensing events are also expected to occur in fields off the LMC's stellar bar since the stellar debris is not expected to be concentrated in the bar region. This scenario leads to a number of observational tests: the sources are low-metallicity SMC stars, they exhibit high velocities relative to LMC disk stars that may be detectable via proper motion studies, and, most notably, there should exist a stellar counterpart to the gaseous Magellanic Stream and Bridge with a V-band surface brightness of > 32 mag/arcsec^2. In particular, the stellar Bridge should contain enough RR Lyrae stars to be detected by the ongoing OGLE survey of this region.

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

We examine Milky Way-Magellanic Cloud systems selected from the Millennium-II Simulation in order to place the orbits of the Magellanic Clouds in a cosmological context. Our analysis shows that satellites massive enough to be LMC analogs are typically accreted at late times. Moreover, those that are accreted at early times and survive to the present have orbital properties that are discrepant with those observed for the LMC. The high velocity of the LMC, coupled with the dearth of unbound orbits seen in the simulation, argues that the mass of the MW's halo is unlikely to be less than 2 x 10^12 Msun. This conclusion is further supported by statistics of halos hosting satellites with masses, velocities, and separations comparable to those of the LMC. We further show that: (1) LMC and SMC-mass objects are not particularly uncommon in MW-mass halos; (2) the apparently high angular momentum of the LMC is not cosmologically unusual; and (3) it is rare for a MW halo to host a LMC-SMC binary system at z=0, but high speed binary pairs accreted at late times are possible. Based on these results, we conclude that the LMC was accreted within the past four Gyr and is currently making its first pericentric passage about the MW.

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.

The outer regions of galactic disks have received increased attention since ultraviolet observations with GALEX demonstrated that nearly 30% of galaxies have UV emission beyond their optical extents, indicating star formation activity. These galaxies have been termed extended UV (XUV) disks. Here, we address whether these observations contradict the gas surface density threshold for star formation inferred from Halpha radial profiles of galaxies. We run smoothed particle hydrodynamics simulations of isolated disk galaxies with fiducial star formation prescriptions and show that over-densities owing to the presence of spiral structure can induce star formation in extended gas disks. For direct comparison with observations, we use the 3-D radiative transfer code Sunrise to create simulated FUV and K_s band images. We find that galaxies classified as Type I XUV disks are a natural consequence of spiral patterns, but we are unable to reproduce Type II XUV disks. We also compare our results to studies of the Kennicutt-Schmidt relation in outer disks.

Dwarf spheroidal galaxies are the most dark matter dominated systems in the nearby Universe and their origin is one of the outstanding puzzles of how galaxies form. Dwarf spheroidals are poor in gas and stars, making them unusually faint, and those known as ultra-faint dwarfs have by far the lowest measured stellar content of any galaxy. Previous theories require that dwarf spheroidals orbit near giant galaxies like the Milky Way, but some dwarfs have been observed in the outskirts of the Local Group. Here we report simulations of encounters between dwarf disk galaxies and somewhat larger objects. We find that the encounters excite a process, which we term ``resonant stripping'', that can transform them into dwarf spheroidals. This effect is distinct from other mechanisms proposed to form dwarf spheroidals, including mergers, galaxy-galaxy harassment, or tidal and ram pressure stripping, because it is driven by gravitational resonances. It may account for the observed properties of dwarf spheroidals in the Local Group, including their morphologies and kinematics. Resonant stripping predicts that dwarf spheroidals should form through encounters, leaving detectable long stellar streams and tails.

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

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.

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.

Previous models of galactic disk heating in interactions invoke restrictive assumptions not necessarily valid in modern LCDM contexts: that satellites and orbits are rigid and circular, with slow decay over many orbital times from dynamical friction. This leads to a linear scaling of disk heating with satellite mass: disk heights and velocity dispersions scale ~M_sat/M_disk. In turn, observed disk thicknesses present strong constraints on merger histories: the implication for the Milky Way is that <5% of its mass could come from mergers since z~2, in conflict with cosmological predictions. More realistically, satellites merge on nearly radial orbits, and once near the disk, resonant interactions efficiently remove angular momentum while tidal effects strip mass, leading to rapid merger/destruction in a couple of free-fall plunges. Under these conditions the proper heating efficiency is non-linear in mass ratio, ~(M_sat/M_disk)^2. We derive the scaling of disk scale heights and velocity dispersions as a function of mass ratio and disk gas content in this regime, and show this accurately describes the results of simulations with proper 'live' halos and disks. Under realistic circumstances, disk heating in minor mergers is suppressed by an order of magnitude relative to expectations of previous models. We show that the Milky Way disk could have absorbed ~5-10 1:10 mass-ratio mergers since z=2, in agreement with cosmological models. These distinctions lead to dramatic differences in which mass ratios are most important for disk heating and in the isophotal shapes of disk+bulge systems.

We present a scenario for the production of dynamically cold rings around the Milky Way via a high-eccentricity, flyby encounter. These initial conditions are more cosmologically motivated than those considered in previous works. We find that the encounters we examine generically produce a series of nearly dynamically cold ring-like features on low-eccentricity orbits that persist over timescales of ~2-4 Gyr via the tidal response of the primary galaxy to the close passage of the satellite. Moreover, they are both qualitatively and quantitatively similar to the distribution, kinematics, and stellar population of the Monoceros ring. Therefore, we find that a high eccentricity flyby by a satellite galaxy represents a cosmologically appealing scenario for forming kinematically distinct ring-like features around the Milky Way.

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.

Narrow dust rings observed around some young stars (e.g., HR 4796A) need to be confined. We present a possible explanation for the formation and confinement of such rings in optically thin circumstellar disks, without invoking shepherding planets. If an enhancement of dust grains (e.g., due to a catastrophic collision) occurs somewhere in the disk, photoelectric emission from the grains can heat the gas to temperatures well above that of the dust. The gas orbits with super(sub)-Keplerian speeds inward (outward) of the associated pressure maximum. This tends to concentrate the grains into a narrow region. The rise in dust density leads to further heating and a stronger concentration of grains. A narrow dust ring forms as a result of this instability. We show that this mechanism not only operates around early-type stars that have high UV fluxes, but also around stars with spectral types as late as K. This implies that this process is generic and may have occurred during the lifetime of each circumstellar disk. We examine the stringent upper-limit on the H2 column density in the HR 4796A disk and find it to be compatible with the presence of a significant amount of hydrogen gas in the disk. We also compute the OI and CII infrared line fluxes expected from various debris disks and show that these will be easily detectable by the upcoming Herschel mission. Herschel will be instrumental in detecting and characterizing gas in these disks.