-
Braving the Storm I: Quantifying Disk-wide Ionized Outflows in the Large Magellanic Cloud with ULLYSES
Authors:
Yong Zheng,
Kirill Tchernyshyov,
Knut Olsen,
Yumi Choi,
Chad Bustard,
Julia Roman-Duval,
Robert Zhu,
Enrico M. Di Teodoro,
Jessica Werk,
Mary Putman,
Anna F. McLeod,
Yakov Faerman,
Raymond C. Simons,
Joshua Peek
Abstract:
The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines observed with the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Proc…
▽ More
The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines observed with the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regression and identify reliable LMC interstellar absorption lines over $v_{\rm helio}=175-375$ km s$^{-1}$. Our analyses show disk-wide, warm, ionized outflows in SiIV and CIV across the LMC with bulk velocities of $|v_{\rm out, bulk}|\sim20-60$ km s$^{-1}$, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities ($Σ_{\rm SFR}$), and the outflows in star-forming regions with higher $Σ_{\rm SFR}$ tend to be more ionized. Considering the ionized outflows from both sides of the LMC as traced by CIV, we conservatively estimate a total outflow rate of $\dot{M}_{\rm out}\gtrsim 0.03~M_\odot {\rm yr}^{-1}$ and a mass loading factor of $η\gtrsim 0.2$. We compare the LMC's outflows with those detected in starburst galaxies, and find a universal scaling relation of $|v_{\rm out, bulk}|\propto Σ_{\rm SFR}^{0.22}$ over a range of star-forming conditions ($Σ_{\rm SFR}\sim10^{-2}-10^{2}~M_\odot {\rm yr}^{-1} {\rm kpc}^{-2}$). Lastly, we find that the outflows are co-rotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Deep Learning Cosmic Ray Transport from Density Maps of Simulated, Turbulent Gas
Authors:
Chad Bustard,
John Wu
Abstract:
The coarse-grained propagation of Galactic cosmic rays (CRs) is traditionally constrained by phenomenological models of Milky Way CR propagation fit to a variety of direct and indirect observables; however, constraining the fine-grained transport of CRs along individual magnetic field lines -- for instance, diffusive vs streaming transport models -- is an unsolved challenge. Leveraging a recent tr…
▽ More
The coarse-grained propagation of Galactic cosmic rays (CRs) is traditionally constrained by phenomenological models of Milky Way CR propagation fit to a variety of direct and indirect observables; however, constraining the fine-grained transport of CRs along individual magnetic field lines -- for instance, diffusive vs streaming transport models -- is an unsolved challenge. Leveraging a recent training set of magnetohydrodynamic turbulent box simulations, with CRs spanning a range of transport parameters, we use convolutional neural networks (CNNs) trained solely on gas density maps to classify CR transport regimes. We find that even relatively simple CNNs can quite effectively classify density slices to corresponding CR transport parameters, distinguishing between streaming and diffusive transport, as well as magnitude of diffusivity, with class accuracies between $92\%$ and $99\%$. As we show, the transport-dependent imprints that CRs leave on the gas are not all tied to the resulting density power spectra: classification accuracies are still high even when image spectra are flattened ($85\%$ to $98\%$ accuracy), highlighting CR transport-dependent changes to turbulent phase information. We interpret our results with saliency maps and image modifications, and we discuss physical insights and future applications.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
The Impact of Cosmic Rays on Thermal and Hydrostatic Stability in Galactic Halos
Authors:
Tsun Hin Navin Tsung,
S. Peng Oh,
Chad Bustard
Abstract:
We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As m…
▽ More
We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As mass dropout from thermal instability causes the background density and hence plasma $β\equiv P_g/P_B$ to fall, the CGM becomes globally unstable. At the cool disk to hot halo interface, a sharp drop in density boosts Alfven speeds and CR gradients, driving a transition from diffusive to streaming transport. CR forces and heating strengthen, while countervailing gravitational forces and radiative cooling weaken, resulting in a loss of both hydrostatic and thermal equilibrium. In lower $β$ halos, CR heating drives a hot, single-phase diffuse wind with velocities $v \propto (t_\mathrm{heat}/t_\mathrm{ff})^{-1}$, which exceeds the escape velocity when $t_\mathrm{heat}/t_\mathrm{ff} \lesssim 0.4$. In higher $β$ halos, CR forces drive multi-phase winds with cool, dense fountain flows and significant turbulence. These flows are CR dominated due to "trapping" of CRs by weak transverse B-fields, and have the highest mass loading factors. Thus, local thermal instability can result in winds or fountain flows where either the heat or momentum input of CRs dominates.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
Cosmic Ray Drag and Damping of Compressive Turbulence
Authors:
Chad Bustard,
S. Peng Oh
Abstract:
While it is well-known that cosmic rays (CRs) can gain energy from turbulence via second order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time $t_{\rm damp} \sim ρv^{2}/\dot{E}_{\rm CR} \propto E_{\rm CR}^{-1}$ becomes compa…
▽ More
While it is well-known that cosmic rays (CRs) can gain energy from turbulence via second order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time $t_{\rm damp} \sim ρv^{2}/\dot{E}_{\rm CR} \propto E_{\rm CR}^{-1}$ becomes comparable to the turbulent cascade time. Magnetohydrodynamic (MHD) simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating at $\dot{E}_{\rm CR} \sim \tildeε$, where $\tildeε$ is the turbulent energy input rate. In that case, almost all the energy in large scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge-Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering of $E > 300$ GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large scale motions, with kinetic energy reduced by up to to an order of magnitude in realistic $E_{\rm CR} \sim E_{\rm g}$ scenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e. $\tildeε \gg ρv^{3}/L$.
△ Less
Submitted 14 August, 2023; v1 submitted 10 January, 2023;
originally announced January 2023.
-
Turbulent Reacceleration of Streaming Cosmic Rays
Authors:
Chad Bustard,
S. Peng Oh
Abstract:
Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as non-resonant reacceleration. It is often invoked to explain observed ratios of primary to secondary CRs at $\sim \rm GeV$ energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic…
▽ More
Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as non-resonant reacceleration. It is often invoked to explain observed ratios of primary to secondary CRs at $\sim \rm GeV$ energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasma $β$. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-$β$ environments like the interstellar medium (ISM) but remain unchanged in high-$β$ environments like the intracluster medium (ICM). We also quantify the streaming energy loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution-dependent (hence unconverged in large scale simulations) and heavily suppressed -- by an order of magnitude -- compared to the isotropic loss rate $v_{A} \cdot \nabla P_{\rm CR} / P_{\rm CR} \sim v_{A}/L_{0}$, due to misalignment between the mean field and isotropic CR gradients. Counterintuitively, and unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength over $β\sim 1-100$ and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper (Bustard and Oh, in prep), we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent power spectrum.
△ Less
Submitted 19 January, 2023; v1 submitted 3 August, 2022;
originally announced August 2022.
-
Radiative Turbulent Mixing Layers and the Survival of Magellanic Debris
Authors:
Chad Bustard,
Max Gronke
Abstract:
The Magellanic Stream is sculpted by its infall through the Milky Way's circumgalactic medium, but the rates and directions of mass, momentum, and energy exchange through the Stream-halo interface are relative unknowns critical for determining the origin and fate of the Stream. Complementary to large-scale simulations of LMC-SMC interactions, we apply new insights derived from idealized, high-reso…
▽ More
The Magellanic Stream is sculpted by its infall through the Milky Way's circumgalactic medium, but the rates and directions of mass, momentum, and energy exchange through the Stream-halo interface are relative unknowns critical for determining the origin and fate of the Stream. Complementary to large-scale simulations of LMC-SMC interactions, we apply new insights derived from idealized, high-resolution "cloud-crushing" and radiative turbulent mixing layer simulations to the Leading Arm and Trailing Stream. Contrary to classical expectations of fast cloud breakup, we predict that the Leading Arm and much of the Trailing Stream should be surviving infall and even gaining mass due to strong radiative cooling. Provided a sufficiently supersonic tidal swing-out from the Clouds, the present-day Leading Arm could be a series of high-density clumps in the cooling tail behind the progenitor cloud. We back up our analytic framework with a suite of converged wind-tunnel simulations, finding that previous results on cloud survival and mass growth can be extended to high Mach number ($\mathcal{M}$) flows with a modified drag time $t_{drag} \propto 1 + \mathcal{M}$ and longer growth time. We also simulate the Trailing Stream; we find that the growth time is long ($\sim$ Gyrs) compared to the infall time, and approximate H$α$ emission is low on average ($\sim$few mR) but can be up to tens of mR in bright spots. Our findings also have broader extragalactic implications for e.g. galactic winds, which we discuss.
△ Less
Submitted 21 June, 2022; v1 submitted 18 August, 2021;
originally announced August 2021.
-
Mapping the Supernovae Driven Winds of the Large Magellanic Cloud in H$α$ Emission I
Authors:
Drew A. Ciampa,
Kathleen A. Barger,
Nicolas Lehner,
Madeline Horn,
Michael Hernandez,
L. Matthew Haffner,
Brianna Smart,
Chad Bustard,
Sam Barber,
Henry Boot
Abstract:
We present the first spectroscopically resolved \ha\ emission map of the Large Magellanic Cloud's (LMC) galactic wind. By combining new Wisconsin H-alpha Mapper (WHAM) observations ($I_{\rm Hα}\gtrsim10~{\rm mR}$) with existing \hicm\ emission observations, we have (1) mapped the LMC's near-side galactic wind over a local standard of rest (LSR) velocity range of…
▽ More
We present the first spectroscopically resolved \ha\ emission map of the Large Magellanic Cloud's (LMC) galactic wind. By combining new Wisconsin H-alpha Mapper (WHAM) observations ($I_{\rm Hα}\gtrsim10~{\rm mR}$) with existing \hicm\ emission observations, we have (1) mapped the LMC's near-side galactic wind over a local standard of rest (LSR) velocity range of $+50\le\rm v_{LSR}\le+250~{\rm km}~{\rm s}^{-1}$, (2) determined its morphology and extent, and (3) estimated its mass, outflow rate, and mass-loading factor. We observe \ha\ emission from this wind to typically 1-degree off the LMC's \hi\ disk. Kinematically, we find that the diffuse gas in the warm-ionized phase of this wind persists at both low ($\lesssim100~{\rm km}~{\rm s}^{-1}$) and high ($\gtrsim100~{\rm km}~{\rm s}^{-1}$) velocities, relative to the LMC's \hi\ disk. Furthermore, we find that the high-velocity component spatially aligns with the most intense star-forming region, 30~Doradus. We, therefore, conclude that this high-velocity material traces an active outflow. We estimate the mass of the warm ($T_e\approx10^4~\rm K$) ionized phase of the near-side LMC outflow to be $\log{\left(M_{\rm ionized}/M_\odot\right)=7.51\pm0.15}$ for the combined low and high velocity components. Assuming an ionization fraction of 75\% and that the wind is symmetrical about the LMC disk, we estimate that its total (neutral and ionized) mass is $\log{\left(M_{\rm total}/M_\odot\right)=7.93}$, its mass-flow rate is $\dot{M}_{\rm outflow}\approx1.43~M_\odot~\rm yr^{-1}$, and its mass-loading factor is $η\approx4.54$. Our average mass-loading factor results are roughly a factor of 2.5 larger than previous \ha\ imaging and UV~absorption line studies, suggesting that those studies are missing nearly half the gas in the outflows.
△ Less
Submitted 28 January, 2021;
originally announced January 2021.
-
Cosmic Ray Transport, Energy Loss, and Influence in the Multiphase Interstellar Medium
Authors:
Chad Bustard,
Ellen G. Zweibel
Abstract:
The bulk propagation speed of GeV-energy cosmic rays is limited by frequent scattering off hydromagnetic waves. Most galaxy evolution simulations that account for this confinement assume the gas is fully ionized and cosmic rays are well-coupled to Alfvén waves; however, multiphase density inhomogeneities, frequently under-resolved in galaxy evolution simulations, induce cosmic ray collisions and i…
▽ More
The bulk propagation speed of GeV-energy cosmic rays is limited by frequent scattering off hydromagnetic waves. Most galaxy evolution simulations that account for this confinement assume the gas is fully ionized and cosmic rays are well-coupled to Alfvén waves; however, multiphase density inhomogeneities, frequently under-resolved in galaxy evolution simulations, induce cosmic ray collisions and ionization-dependent transport driven by cosmic ray decoupling and elevated streaming speeds in partially neutral gas. How do cosmic rays navigate and influence such a medium, and can we constrain this transport with observations? In this paper, we simulate cosmic ray fronts impinging upon idealized, partially neutral clouds and lognormally-distributed clumps, with and without ionization-dependent transport. With these high-resolution simulations, we identify cloud interfaces as crucial regions where cosmic ray fronts can develop a stair-step pressure gradient sufficient to collisionlessly generate waves, overcome ion-neutral damping, and exert a force on the cloud. We find that the acceleration of cold clouds is hindered by only a factor of a few when ionization-dependent transport is included, with additional dependencies on magnetic field strength and cloud dimensionality. We also probe how cosmic rays sample the background gas and quantify collisional losses. Hadronic gamma-ray emission maps are qualitatively different when ionization-dependent transport is included, but the overall luminosity varies by only a small factor, as the short cosmic ray residence times in cold clouds are offset by the higher densities that cosmic rays sample.
△ Less
Submitted 22 April, 2021; v1 submitted 11 December, 2020;
originally announced December 2020.
-
The Magellanic Corona and the formation of the Magellanic Stream
Authors:
Scott Lucchini,
Elena D'Onghia,
Andrew J Fox,
Chad Bustard,
Joss Bland-Hawthorn,
Ellen Zweibel
Abstract:
The dominant gaseous structure in the Galactic halo is the Magellanic Stream, an extended network of neutral and ionized filaments surrounding the Large and Small Magellanic Clouds (LMC/SMC), the two most massive satellite galaxies of the Milky Way. Recent observations indicate that the Clouds are on their first passage around our Galaxy, the Stream is made up of gas stripped from both the LMC and…
▽ More
The dominant gaseous structure in the Galactic halo is the Magellanic Stream, an extended network of neutral and ionized filaments surrounding the Large and Small Magellanic Clouds (LMC/SMC), the two most massive satellite galaxies of the Milky Way. Recent observations indicate that the Clouds are on their first passage around our Galaxy, the Stream is made up of gas stripped from both the LMC and the SMC, and the majority of this gas is ionized. While it has long been suspected that tidal forces and ram-pressure stripping contributed to the Stream's formation, a full understanding of its origins has defied modelers for decades. Several recent developments, including the discovery of dwarf galaxies associated with the Magellanic Group, the high mass of the LMC, the detection of highly ionized gas toward stars in the LMC and the predictions of cosmological simulations all support the existence of a halo of warm ionized gas around the LMC at a temperature of $\sim5\times10^{5}\;\mathrm{K}$. Here we show that by including this "Magellanic Corona" in hydrodynamic simulations of the Magellanic Clouds falling onto the Galaxy, we can simultaneously reproduce the Stream and its Leading Arm. Our simulations explain the Stream's filamentary structure, spatial extent, radial velocity gradient, and total ionized gas mass. We predict that the Magellanic Corona will be unambiguously observable via high-ionization absorption lines in the ultraviolet spectra of background quasars lying near the LMC. This prediction is directly testable with the Cosmic Origins Spectrograph on the Hubble Space Telescope.
△ Less
Submitted 9 September, 2020;
originally announced September 2020.
-
Cosmic Ray Driven Outflows from the Large Magellanic Cloud: Contributions to the LMC Filament
Authors:
Chad Bustard,
Ellen G. Zweibel,
Elena D'Onghia,
J. S. Gallagher III,
Ryan Farber
Abstract:
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 r…
▽ More
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.
△ Less
Submitted 27 April, 2020; v1 submitted 5 November, 2019;
originally announced November 2019.
-
The Role of the Parker Instability in Structuring the Interstellar Medium
Authors:
Evan Heintz,
Chad Bustard,
Ellen Zweibel
Abstract:
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 (…
▽ More
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 ($γ_c = 0$), locked to the thermal gas ($γ_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.
△ Less
Submitted 18 March, 2020; v1 submitted 8 October, 2019;
originally announced October 2019.
-
The Propagation of Cosmic Rays from the Galactic Wind Termination Shock: Back to the Galaxy?
Authors:
Lukas Merten,
Chad Bustard,
Ellen G. Zweibel,
Julia Becker Tjus
Abstract:
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 con…
▽ More
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.
△ Less
Submitted 23 April, 2018; v1 submitted 22 March, 2018;
originally announced March 2018.
-
The Fate of Supernova-Heated Gas in Star-Forming Regions of the LMC: Lessons for Galaxy Formation?
Authors:
Chad Bustard,
Stephen A. Pardy,
Elena D'Onghia,
Ellen G. Zweibel,
J. S. Gallagher III
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 July, 2018; v1 submitted 20 February, 2018;
originally announced February 2018.
-
Cosmic Ray Acceleration by a Versatile Family of Galactic Wind Termination Shocks
Authors:
Chad Bustard,
Ellen G. Zweibel,
Cory Cotter
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 December, 2016; v1 submitted 20 October, 2016;
originally announced October 2016.
-
A Versatile Family of Galactic Wind Models
Authors:
Chad Bustard,
Ellen G. Zweibel,
Elena D'Onghia
Abstract:
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-l…
▽ More
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.
△ Less
Submitted 27 January, 2016; v1 submitted 23 September, 2015;
originally announced September 2015.
