The Structure of Multiphase Galactic Winds

We present a novel analytic framework to model the steady-state structure of multiphase galactic winds comprised of a hot, volume-filling component and a cold, clumpy component. We first derive general expressions for the structure of the hot phase for arbitrary mass, momentum, and energy source terms. Next, informed by recent simulations, we parameterize the cloud–wind mass transfer rates, which are set by the competition between turbulent mixing and radiative cooling. This enables us to cast the cloud–wind interaction as a source term for the hot phase and thereby simultaneously solve for the evolution of both phases, fully accounting for their bidirectional influence. With this model, we explore the nature of galactic winds over a broad range of conditions. We find that (i) with realistic parameter choices, we naturally produce a hot, low-density wind that transports energy while entraining a significant flux of cold clouds, (ii) mixing dominates the cold cloud acceleration and decelerates the hot wind, (iii) during mixing thermalization of relative kinetic energy provides significant heating, (iv) systems with low hot phase mass loading factors and/or star formation rates can sustain higher initial cold phase mass loading factors, but the clouds are quickly shredded, and (v) systems with large hot phase mass loading factors and/or high star formation rates cannot sustain large initial cold phase mass loading factors, but the clouds tend to grow with distance from the galaxy. Our results highlight the necessity of accounting for the multiphase structure of galactic winds, both physically and observationally, and have important implications for feedback in galactic systems.

The radio continuum perspective on cosmic-ray transport in external galaxies

Radio continuum observations of external galaxies provide us with an excellent outside view on the distribution of cosmic-ray electrons in the disc and halo. In this review, we summarise the current state of what we have learned from modelling such observations with cosmic-ray transport, paying particular attention to the question to what extent we can exploit radio haloes when studying galactic winds. We have developed the user-friendly framework spinnaker to model radio haloes with either pure advection or diffusion, allowing us to study both diffusion coefficients and advection speeds in nearby galaxies. Using these models, we show that we can identify galaxies with winds using both morphology and radio spectral indices of radio haloes. Advective radio haloes are ubiquitous, indicating that already fairly low values of the star formation rate (SFR) surface density ($\Sigma_{ \mathrm{SFR}}$ Σ SFR ) can trigger galactic winds. The advection speeds scale with SFR, $\Sigma_{\mathrm{SFR}}$ Σ SFR , and rotation speed as expected for stellar feedback-driven winds. Accelerating winds are in agreement with our radio spectral index data, but this is sensitive to the magnetic field parametrisation, so that constant wind speeds cannot be ruled out either. The question to what extent cosmic rays can be a driving force behind winds is still an open issue and we discuss only in passing how a simple iso-thermal wind model could fit our data. Nevertheless, the comparison with inferences from observations and theory looks promising with radio continuum offering a complementary view on galactic winds. We finish with a perspective on future observations and challenges lying ahead.

Characterizing mass, momentum, energy and metal outflow rates of multi-phase galactic winds in the FIRE-2 cosmological simulations

We characterize mass, momentum, energy and metal outflow rates of multi-phase galactic winds in a suite of FIRE-2 cosmological ‘zoom-in’ simulations from the Feedback in Realistic Environments (FIRE) project. We analyse simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass haloes, and high-redshift massive haloes. Consistent with previous work, we find that dwarfs eject about 100 times more gas from their interstellar medium (ISM) than they form in stars, while this mass ‘loading factor’ drops below one in massive galaxies. Most of the mass is carried by the hot phase (>105 K) in massive haloes and the warm phase (103 − 105 K) in dwarfs; cold outflows (<103 K) are negligible except in high-redshift dwarfs. Energy, momentum and metal loading factors from the ISM are of order unity in dwarfs and significantly lower in more massive haloes. Hot outflows have 2 − 5 × higher specific energy than needed to escape from the gravitational potential of dwarf haloes; indeed, in dwarfs, the mass, momentum, and metal outflow rates increase with radius whereas energy is roughly conserved, indicating swept up halo gas. Burst-averaged mass loading factors tend to be larger during more powerful star formation episodes and when the inner halo is not virialized, but we see effectively no trend with the dense ISM gas fraction. We discuss how our results can guide future controlled numerical experiments that aim to elucidate the key parameters governing galactic winds and the resulting associated preventative feedback.

Dust entrainment in galactic winds

ABSTRACT Winds driven by stellar feedback are an essential part of the galactic ecosystem and are the main mechanism through which low-mass galaxies regulate their star formation. These winds are generally observed to be multiphase with detections of entrained neutral and molecular gas. They are also thought to enrich the circumgalactic medium around galaxies with metals and dust. This ejected dust encodes information about the integrated star formation and outflow history of the galaxy. Therefore it is important to understand how much dust is entrained and driven out of the disc by galactic winds. Here, we demonstrate that stellar feedback is efficient in driving dust-enriched winds and eject enough material to account for the amount of extraplanar dust observed in nearby galaxies. The amount of ejected dust depends on the sites from where they are launched, with dustier galaxies launching more dust-enriched outflows. Moreover, the outflowing cold and dense gas is significantly more dust enriched than the volume filling hot and tenuous material. These results provide an important new insight into the dynamics, structure, and composition of galactic winds and their role in determining the dust content of the extragalactic gas in galaxies.

A high-resolution cosmological simulation of a strong gravitational lens

We present a cosmological hydrodynamical simulation of a 1013 M⊙ galaxy group and its environment (out to 10 times the virial radius) carried out using the Eagle model of galaxy formation. Exploiting a novel technique to increase the resolution of the dark matter calculation independently of that of the gas, the simulation resolves dark matter haloes and subhaloes of mass 5 × 106 M⊙. It is therefore useful for studying the abundance and properties of the haloes and subhaloes targeted in strong lensing tests of the cold dark matter model. We estimate the halo and subhalo mass functions and discuss how they are affected both by the inclusion of baryons in the simulation and by the environment. We find that the halo and subhalo mass functions have lower amplitude in the hydrodynamical simulation than in its dark matter only counterpart. This reflects the reduced growth of haloes in the hydrodynamical simulation due to the early loss of gas by reionisation and galactic winds and, additionally, in the case of subhaloes, disruption by enhanced tidal effects within the host halo due to the presence of a massive central galaxy. The distribution of haloes is highly anisotropic reflecting the filamentary character of mass accretion onto the cluster. As a result, there is significant variation in the number of structures with viewing direction. The median number of structures near the centre of the halo, when viewed in projection, is reduced by a factor of two when baryons are included.

Joint Suzaku and Chandra observations of the MKW4 galaxy group out to the virial radius

We present joint Suzaku and Chandra observations of MKW4. With a global temperature of 1.6 keV, MKW4 is one of the smallest galaxy groups that have been mapped in X-rays out to the virial radius. We measure its gas properties from its center to the virial radius in the north, east, and northeast directions. Its entropy profile follows a power-law of ∝r1.1 between R500 and R200 in all directions, as expected from the purely gravitational structure formation model. The well-behaved entropy profiles at the outskirts of MKW4 disfavor the presence of gas clumping or thermal non-equilibrium between ions and electrons in this system. We measure an enclosed baryon fraction of 11% at R200, remarkably smaller than the cosmic baryon fraction of 15%. We note that the enclosed gas fractions at R200 are systematically smaller for groups than for clusters from existing studies in the literature. The low baryon fraction of galaxy groups, such as MKW4, suggests that their shallower gravitational potential well may make them more vulnerable to baryon losses due to AGN feedback or galactic winds. We find that the azimuthal scatter of various gas properties at the outskirts of MKW4 is significantly lower than in other systems, suggesting that MKW4 is a spherically symmetric and highly relaxed system.

A new model for including galactic winds in simulations of galaxy formation – I. Introducing the Physically Evolved Winds (PhEW) model

ABSTRACT The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is oversimplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds, that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a subgrid prescription to model galactic winds more robustly both physically and numerically.

Cosmic ray feedback from supernovae in dwarf galaxies

The regulation of the baryonic content in dwarf galaxies is a long-standing problem. Supernovae (SNe) are supposed to play a key role in forming large-scale galactic winds by removing important amounts of gas from galaxies. SNe are efficient accelerators of non-thermal particles, so-called cosmic rays (CRs), which can substantially modify the dynamics of the gas and conditions to form large-scale galactic winds. We investigate how CR injection by SNe impacts the star formation and the formation of large-scale winds in dwarf galaxies, and whether it can produce galaxy star-formation rates (SFR) and wind properties closer to observations. We ran CR magneto-hydrodynamical simulations of dwarf galaxies at high resolution (9 pc) with the adaptive mesh refinement code RAMSES. Those disc galaxies are embedded in isolated halos of mass of 1010 and 1011 M⊙, and CRs are injected by SNe. We included CR isotropic and anisotropic diffusion with various diffusion coefficients, CR radiative losses, and CR streaming. The injection of CR energy into the interstellar medium smooths out the highest gas densities, which reduces the SFR by a factor of 2–3. Mass outflow rates are significantly greater with CR diffusion, by 2 orders of magnitudes for the higher diffusion coefficients. Without diffusion and streaming, CRs are inefficient at generating winds. CR streaming alone allows for the formation of winds but which are too weak to match observations. The formation of galactic winds strongly depends on the diffusion coefficient: for low coefficients, CR energy stays confined in high density regions where CR energy losses are highest, and higher coefficients, which allow for a more efficient leaking of CRs out of dense gas, produce stronger winds. CR diffusion leads to colder and denser winds than without CRs, and brings outflow rates and mass loading factors much closer to observations.