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.

The Parker instability, a Rayleigh-Taylor like instability of thermal gas supported against gravity by magnetic fields and cosmic rays, is thought to be dynamically important for galaxy evolution, possibly promoting molecular cloud formation and the galactic dynamo. In previous work, we examined the effect of three different cosmic ray transport models on the Parker instability: decoupled ($\gamma_c = 0$), locked to the thermal gas ($\gamma_c = 4/3$) and coupled to the gas with streaming by self-confinement. We expand upon that work here by considering radiative cooling, a smooth gravitational potential, and simulations into the nonlinear regime. We determine that cosmic ray transport away from compression points, whether by diffusion or streaming, is the largest driver of the instability. Heating due to cosmic ray streaming is also destabilizing and especially affects the nonlinear regime. While cooling de-pressurizes the dense gas, streaming cosmic rays heat and inflate the diffuse extraplanar gas, greatly modifying the phase structure of the medium. In 3D, we find that the fastest growth favors short wavelength modes in the horizontal direction perpendicular to the background magnetic field; this is imprinted on Faraday rotation measure maps that may be used to detect the Parker instability. The modifications to the Parker instability that we observe in this work have large implications for the structure and evolution of galaxies, and they highlight the major role that cosmic rays play in shaping their environments.

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.

Although several theories for the origin of cosmic rays in the region between the spectral `knee' and `ankle' exist, this problem is still unsolved. A variety of observations suggest that the transition from Galactic to extragalactic sources occurs in this energy range. In this work we examine whether a Galactic wind which eventually forms a termination shock far outside the Galactic plane can contribute as a possible source to the observed flux in the region of interest. Previous work by Bustard et al. (2017) estimated that particles can be accelerated up to energies above the `knee' up to $R_\mathrm{max} = 10^{16}$ eV for parameters drawn from a model of a Milky Way wind (Everett et al. 2017). A remaining question is whether the accelerated cosmic rays can propagate back into the Galaxy. To answer this crucial question, we simulate the propagation of the cosmic rays using the low energy extension of the CRPropa framework, based on the solution of the transport equation via stochastic differential equations. The setup includes all relevant processes, including three-dimensional anisotropic spatial diffusion, advection, and corresponding adiabatic cooling. We find that, assuming realistic parameters for the shock evolution, a possible Galactic termination shock can contribute significantly to the energy budget in the `knee' region and above. We estimate the resulting produced neutrino fluxes and find them to be below measurements from IceCube and limits by KM3NeT.

There are two distinct breaks in the cosmic ray (CR) spectrum: the so-called "knee" around $3 \times 10^{15}$ eV and the so-called "ankle" around $10^{18}$ eV. Diffusive shock acceleration (DSA) at supernova remnant (SNR) shock fronts is thought to accelerate galactic CRs to energies below the knee, while an extragalactic origin is presumed for CRs with energies beyond the ankle. CRs with energies between $3 \times 10^{15}$ and $10^{18}$ eV, which we dub the "shin," have an unknown origin. It has been proposed that DSA at galactic wind termination shocks, rather than at SNR shocks, may accelerate CRs to these energies. This paper uses the galactic wind model of Bustard et al. (2016) to analyze whether galactic wind termination shocks may accelerate CRs to shin energies within a reasonable acceleration time and whether such CRs can subsequently diffuse back to the galaxy. We argue for acceleration times on the order of 100 Myrs rather than a few billion years, as assumed in some previous works, and we discuss prospects for magnetic field amplification at the shock front. Ultimately, we generously assume that the magnetic field is amplified to equipartition. This formalism allows us to obtain analytic formulae, applicable to any wind model, for CR acceleration. Even with generous assumptions, we find that very high wind velocities are required to set up the necessary conditions for acceleration beyond $10^{17}$ eV. We also estimate the luminosities of CRs accelerated by outflow termination shocks, including estimates for the Milky Way wind.

We present a versatile family of model galactic outflows including non-uniform mass and energy source distributions, a gravitational potential from an extended mass source, and radiative losses. The model easily produces steady-state wind solutions for a range of mass-loading factors, energy-loading factors, galaxy mass and galaxy radius. We find that, with radiative losses included, highly mass-loaded winds must be driven at high central temperatures, whereas low mass-loaded winds can be driven at low temperatures just above the peak of the cooling curve, meaning radiative losses can drastically affect the wind solution even for low mass-loading factors. By including radiative losses, we are able to show that subsonic flows can be ignored as a possible mechanism for expelling mass and energy from a galaxy compared to the more efficient transonic solutions. Specifically, the transonic solutions with low mass-loading and high energy-loading are the most efficient. Our model also produces low-temperature, high-velocity winds that could explain the prevalence of low-temperature material in observed outflows. Finally, we show that our model, unlike the well-known Chevalier & Clegg (1985) model, can reproduce the observed linear relationship between wind X-ray luminosity and star formation rate (SFR) over a large range of SFR from $1-1000$ M$_{\odot}$/yr assuming the wind mass-loading factor is higher for low-mass, and hence, low-SFR galaxies. We also constrain the allowed mass-loading factors that can fit the observed X-ray luminosity vs. SFR trend, further suggesting an inverse relationship between mass-loading and SFR as explored in advanced numerical simulations.