We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin.

Galactic winds and fountains driven by supernova-heated gas play an integral role in re-distributing gas in galaxies, depositing metals in the circumgalactic medium (CGM), and quenching star formation. The interplay between these outflows and ram pressure stripping due to the galaxy's motion through an ambient medium may enhance these effects by converting fountain flows into expelled gas. In this paper, we present controlled, 3D simulations of ram pressure stripping combined with thermally driven, local outflows from clustered supernovae in an isolated disk galaxy modeled on the Large Magellanic Cloud (LMC), a dwarf satellite of the Milky Way on its first infall. Observational evidence of local outflows emanating from supergiant shells in the LMC and a trailing filament of HI gas originating from these regions - with no obvious Leading Arm counterpart - may represent a perfect example of this process. Our simulations present a proof-of-concept that ram pressure can convert fountain flows into expelled gas. We find that fountains launched near the peak star formation time of the LMC can comprise part of the LMC filament in the Trailing Stream, but with lower column densities than observed. Larger, more numerous outflows from the LMC may be possible and may contribute more mass, but higher inertia gas will lengthen the timescale for this gas to be swept away by ram pressure. Given the high resolution observations, increased knowledge of star formation histories, and growing evidence of multiphase, ionized outflows, the LMC is an ideal test-bed for future wind models.

We have developed a novel technique based on a clustering algorithm which searches for kinematically- and chemically-clustered stars in the APOGEE DR12 Cannon data. As compared to classical chemical tagging, the kinematic information included in our methodology allows us to identify stars that are members of known globular clusters with greater confidence. We apply our algorithm to the entire APOGEE catalog of 150,615 stars whose chemical abundances are derived by the Cannon. Our methodology found anti-correlations between the elements Al and Mg, Na and O, and C and N previously identified in the optical spectra in globular clusters, even though we omit these elements in our algorithm. Our algorithm identifies globular clusters without a priori knowledge of their locations in the sky. Thus, not only does this technique promise to discover new globular clusters, but it also allows us to identify candidate streams of kinematically- and chemically-clustered stars in the Milky Way.

The Magellanic Stream and Leading Arm of HI that stretches from the Large and Small Magellanic Clouds (LMC and SMC) and over 200 degrees of the Southern sky is thought to be formed from multiple encounters between the LMC and SMC. In this scenario, most of the gas in the Stream and Leading Arm is stripped from the SMC, yet recent observations have shown a bifurcation of the Trailing Arm that reveals LMC origins for some of the gas. Absorption measurements in the Stream also reveal an order of magnitude more gas than in current tidal models. We present hydrodynamical simulations of the multiple encounters between the LMC and SMC at their first pass around the Milky Way, assuming that the Clouds were more extended and gas rich in the past. Our models create filamentary structures of gas in the Trailing Stream from both the LMC and SMC. While the SMC trailing filament matches the observed Stream location, the LMC filament is offset. In addition, the total observed mass of the Stream in these models is underestimated of a factor of four when the ionized component is accounted for. Our results suggest that there should also be gas stripped from both the LMC and SMC in the Leading Arm, mirroring the bifurcation in the Trailing Stream. This prediction is consistent with recent measurements of spatial variation in chemical abundances in the Leading Arm, which show that gas from multiple sources is present, although the nature is still uncertain.

We examine the dust geometry and Ly{\alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{\alpha} emitting and starbursting systems with Hubble Space Telescope Ly{\alpha}, H{\alpha}, and H{\beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{\alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{\alpha}/H{\alpha} and H{\alpha}/H{\beta} maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{\alpha} photons only scatter. We show that the dust that affects the escape of Ly{\alpha} is more restricted to the galaxies' central regions, while the larger Ly{\alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{\alpha} scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Ly{\alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{\alpha} and global dust properties.

We present the first robust detection of HI 21 cm emission in the blue compact galaxy Haro 11 using the 100m Robert C. Byrd Green Bank Telescope (GBT). Haro 11 is a luminous blue compact galaxy with emission in both Lyman Alpha and the Lyman continuum. We detect (5.1 $\pm$ 0.7 $\times$10$^8$) M$_{\odot}$ of HI gas at an assumed distance of 88 Mpc, making this galaxy HI deficient compared to other local galaxies with similar optical properties. Given this small HI mass, Haro 11 has an elevated M$_{H2}$/M$_{HI}$ ratio and a very low gas fraction compared to most local galaxies, and contains twice as much mass in ionized hydrogen as in neutral hydrogen. The HI emission has a linewidth of 71 kms$^{-1}$ and is offset 60 kms$^{-1}$ redward of the optical line center. It is undergoing a starburst after a recent merger which has elevated the star formation rate, and will deplete the gas supply in $<$ 0.2 Gyr. Although this starburst has elevated the SFR compared to galaxies with similar HI masses and linewidths, Haro 11 matches a trend of lower gas fractions toward higher star formation rates and is below the general trend of increasing HI mass with increasing luminosity. Taken together, our results paint Haro 11 as a standard low-mass galaxy that is undergoing an unusually efficient star formation episode.

Magellanic spiral galaxies are a class of one-armed systems that often exhibit an offset stellar bar, and are rarely found around massive spiral galaxies. Using a set of N-body and hydrodynamic simulations we consider a dwarf-dwarf galaxy interaction as the driving mechanism for the formation of this peculiar class of systems. We investigate here the relation between the dynamical, stellar and gaseous disk center and the bar. In all our simulations the bar center always coincides with the dynamical center, while the stellar disk becomes highly asymmetric during the encounter causing the photometric center of the Magellanic galaxy disk to become mismatched with both the bar and the dynamical center. The disk asymmetries persist for almost 2 Gyrs, the time that it takes for the disk to be re-centered with the bar, and well after the companion has passed. This explains the nature of the offset bar found in many Magellanic-type galaxies, including the Large Magellanic Cloud (LMC) and NGC 3906. In particular, these results, once applied to the LMC, suggest that the dynamical center should reside in the bar center instead of the HI center as previously assumed, pointing to a variation in the current estimate of the north component of the LMC proper motion.

We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $\alpha$ (H$\alpha$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$\alpha$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$\alpha$ and Lyman $\alpha$ (Ly$\alpha$) images. Only for individual galaxies a causal connection between spatially resolved H$\alpha$ kinematics and Ly$\alpha$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $\sigma_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/\sigma_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/\sigma_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/\sigma_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$\alpha$ properties, we find that dispersion dominated systems show higher Ly$\alpha$ equivalent widths and higher Ly$\alpha$ escape fractions than systems with $v_\mathrm{shear}/\sigma_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$\alpha$ radiation.

We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.

We present new H I imaging and spectroscopy of the 14 UV-selected star-forming galaxies in the Lyman Alpha Reference Sample (LARS), aimed for a detailed study of the processes governing the production, propagation, and escape of Ly$\alpha$ photons. New H I spectroscopy, obtained with the 100m Green Bank Telescope (GBT), robustly detects the H I spectral line in 11 of the 14 observed LARS galaxies (although the profiles of two of the galaxies are likely confused by other sources within the GBT beam); the three highest redshift galaxies are not detected at our current sensitivity limits. The GBT profiles are used to derive fundamental H I line properties of the LARS galaxies. We also present new pilot H I spectral line imaging of 5 of the LARS galaxies obtained with the Karl G. Jansky Very Large Array (VLA). This imaging localizes the H I gas and provides a measurement of the total H I mass in each galaxy. In one system, LARS 03 (UGC 8335 or Arp 238), VLA observations reveal an enormous tidal structure that extends over 160 kpc from the main interacting systems and that contains $>$10$^9$ M$_{\odot}$ of H I. We compare various H I properties with global Ly$\alpha$ quantities derived from HST measurements. The measurements of the Ly$\alpha$ escape fraction are coupled with the new direct measurements of H I mass and significantly disturbed H I velocities. Our robustly detected sample reveals that both total H I mass and linewidth are tentatively correlated with key Ly$\alpha$ tracers. Further, on global scales, these data support a complex coupling between Ly$\alpha$ propagation and the H I properties of the surrounding medium.

We present new HI spectral line imaging of the extremely metal-poor, star-forming dwarf irregular galaxy Leo P. Our HI images probe the global neutral gas properties and the local conditions of the interstellar medium (ISM). The HI morphology is slightly elongated along the optical major-axis. We do not find obvious signatures of interaction or infalling gas at large spatial scales. The neutral gas disk shows obvious rotation, although the velocity dispersion is comparable to the rotation velocity. The rotation amplitude is estimated to be V_c = 15 +/- 5 km/s. Within the HI radius probed by these observations, the mass ratio of gas to stars is roughly 2:1, while the ratio of the total mass to the baryonic mass is ~15:1. We use this information to place Leo P on the baryonic Tully-Fisher relation, testing the baryonic content of cosmic structures in a sparsely populated portion of parameter space that has hitherto been occupied primarily by dwarf spheroidal galaxies. We detect the signature of two temperature components in the neutral ISM of Leo P; the cold and warm components have characteristic velocity widths of 4.2 +/- 0.9 km/s and 10.1 +/- 1.2 km/s, corresponding to kinetic temperatures of ~1100 K and ~6200 K, respectively. The cold HI component is unresolved at a physical resolution of 200 pc. The highest HI surface densities are observed in close physical proximity to the single HII region. A comparison of the neutral gas properties of Leo P with other extremely metal-deficient (XMD) galaxies reveals that Leo P has the lowest neutral gas mass of any known XMD, and that the dynamical mass of Leo P is more than two orders of magnitude smaller than any known XMD with comparable metallicity.

We report upon new results regarding the Lya output of galaxies, derived from the Lyman alpha Reference Sample (LARS), focusing on Hubble Space Telescope imaging. For 14 galaxies we present intensity images in Lya, Halpha, and UV, and maps of Halpha/Hbeta, Lya equivalent width (EW), and Lya/Halpha. We present Lya and UV light profiles and show they are well-fitted by S\'ersic profiles, but Lya profiles show indices systematically lower than those of the UV (n approx 1-2 instead of >~4). This reveals a general lack of the central concentration in Lya that is ubiquitous in the UV. Photometric growth curves increase more slowly for Lya than the FUV, showing that small apertures may underestimate the EW. For most galaxies, however, flux and EW curves flatten by radii ~10 kpc, suggesting that if placed at high-z, only a few of our galaxies would suffer from large flux losses. We compute global properties of the sample in large apertures, and show total luminosities to be independent of all other quantities. Normalized Lya throughput, however, shows significant correlations: escape is found to be higher in galaxies of lower star formation rate, dust content, mass, and several quantities that suggest harder ionizing continuum and lower metallicity. Eight galaxies could be selected as high-z Lya emitters, based upon their luminosity and EW. We discuss the results in the context of high-z Lya and UV samples. A few galaxies have EWs above 50 AA, and one shows f_escLya of 80%; such objects have not previously been reported at low-z.

We present new low-resolution HI spectral line imaging, obtained with the Karl G. Jansky Very Large Array (JVLA), of the star-forming Magellanic irregular galaxy UGCA 105. This nearby (D = 3.39+/-0.25 Mpc), low mass [M_HI=(4.3+/-0.5)x10^8 Solar masses] system harbors a large neutral gas disk (HI radius ~7.2 kpc at the N_HI=10^20 cm^-2 level) that is roughly twice as large as the stellar disk at the B-band R_25 isophote. We explore the neutral gas dynamics of this system, fitting tilted ring models in order to extract a well-sampled rotation curve. The rotation velocity rises in the inner disk, flattens at 72+/-3 km/s, and remains flat to the last measured point of the disk (~7.5 kpc). The dynamical mass of UGCA 105 at this outermost point, (9+/-2)x10^9 Solar masses, is ~10 times as large as the luminous baryonic components (neutral atomic gas and stars). The proximity and favorable inclination (55 degrees) of UGCA 105 make it a promising target for high-resolution studies of both star formation and rotational dynamics in a nearby low-mass galaxy.