In this paper, we build from previous work (Bustard et al. 2018) and present simulations of recent (within the past Gyr), magnetized, cosmic ray driven outflows from the Large Magellanic Cloud (LMC), including our first attempts to explicitly use the derived star formation history of the LMC to seed outflow generation. We run a parameter set of simulations for different LMC gas masses and cosmic ray transport treatments, and we make preliminary comparisons to published outflow flux estimates, neutral and ionized hydrogen observations, and Faraday rotation measure maps. We additionally report on the gas mass that becomes unbound from the LMC disk and swept by ram pressure into the Trailing Magellanic Stream. We find that, even for our largest outburst, the mass contribution to the Stream is still quite small, as much of the outflow-turned-halo gas is shielded on the LMC's far-side due to the LMC's primarily face-on infall through the Milky Way halo over the past Gyr. On the LMC's near-side, past outflows have fought an uphill battle against ram pressure, with near-side halo mass being at least a factor of a few smaller than the far-side. Absorption line studies probing only the LMC foreground, then, may be severely underestimating the total mass of the LMC halo formed by outflows.

Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) -high-energy charged particles accelerated in supernovae and young stars - can impact the efficiency of wind driving. The streaming instability limits the speed at which they can escape. However, in the presence of turbulence, the streaming instability is subject to suppression that depends on the magnetization of turbulence given by its Alfv\'en Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfv\'enic streaming of CRs enhances galactic winds, in the present paper we take into account a realistic model of streaming suppression. We perform three-dimensional magnetohydrodynamic simulations of a section of a galactic disk and find that turbulent damping dependent on local magnetization of turbulent interstellar medium (ISM) leads to more spatially extended gas and CR distributions compared to the earlier streaming calculations, and that scale-heights of these distributions increase for stronger turbulence. Our results indicate that the star formation rate increases with the level of turbulence in the ISM. We also find that the instantaneous wind mass loading is sensitive to local streaming physics with the mass loading dropping significantly as the strength of turbulence increases.

The role of cosmic rays generated by supernovae and young stars has very recently begun to receive significant attention in studies of galaxy formation and evolution due to the realization that cosmic rays can efficiently accelerate galactic winds. Microscopic cosmic ray transport processes are fundamental for determining the efficiency of cosmic ray wind driving. Previous studies focused on modeling of cosmic ray transport either via constant diffusion coefficient or via streaming proportional to the Alfven speed. However, in predominantly cold, neutral gas, cosmic rays can propagate faster than in the ionized medium and the effective transport can be substantially larger; i.e., cosmic rays can decouple from the gas. We perform three-dimensional magnetohydrodynamical simulations of patches of galactic disks including the effects of cosmic rays. Our simulations include the decoupling of cosmic rays in the cold, neutral interstellar medium. We find that, compared to the ordinary diffusive cosmic ray transport case, accounting for the decoupling leads to significantly different wind properties such as the gas density and temperature, significantly broader spatial distribution of cosmic rays, and larger wind speed. These results have implications for X-ray, gamma-ray and radio emission, and for the magnetization and pollution of the circumgalactic medium by cosmic rays.

Cosmic ray anisotropy has been observed in a wide energy range and at different angular scales by a variety of experiments over the past decade. However, no comprehensive or satisfactory explanation has been put forth to date. The arrival distribution of cosmic rays at Earth is the convolution of the distribution of their sources and of the effects of geometry and properties of the magnetic field through which particles propagate. It is generally believed that the anisotropy topology at the largest angular scale is adiabatically shaped by diffusion in the structured interstellar magnetic field. On the contrary, the medium- and small-scale angular structure could be an effect of non-diffusive propagation of cosmic rays in perturbed magnetic fields. In particular, a possible explanation of the observed small-scale anisotropy observed at TeV energy scale, may come from the effect of particle scattering in turbulent magnetized plasmas. We perform numerical integration of test particle trajectories in low-$\beta$ compressible magnetohydrodynamic turbulence to study how the cosmic rays arrival direction distribution is perturbed when they stream along the local turbulent magnetic field. We utilize Liouville's theorem for obtaining the anisotropy at Earth and provide the theoretical framework for the application of the theorem in the specific case of cosmic ray arrival distribution. In this work, we discuss the effects on the anisotropy arising from propagation in this inhomogeneous and turbulent interstellar magnetic field.