Galactic outflows produced by stellar feedback are known to be multiphase in nature. Observations and simulations indicate that the material within several kiloparsecs of galactic disk midplanes consists of warm clouds embedded within a hot wind. A theoretical understanding of the outflow phenomenon, including both winds and fountain flows, requires study of the interactions among thermal phases. We develop a method to quantify these interactions via measurements of mass, momentum, and energy flux exchanges using temporally and spatially averaged quantities and conservation laws. We apply this method to a star-forming interstellar medium simulation based on the TIGRESS framework, for solar neighborhood conditions. To evaluate the extent of interactions among the phases, we examine the validity of the “ballistic model,” which predicts the trajectories of the warm phase (5050 K < T < 2 × 104 K) treated as non-interacting clouds. This model is successful at intermediate vertical velocities ( ), but at higher velocities, we observe an excess in simulated warm outflow compared to the ballistic model. This discrepancy cannot be fully accounted for by cooling of high-velocity, intermediate-temperature (2 × 104 K < T < 5 × 105 K) gas. We examine the fluxes of mass, momentum, and energy and conclude that the warm phase gains mass via cooling of the intermediate phase and momentum from the hot (T > 5 × 105 K) phase. The large energy flux from the hot outflow, transferred to the warm and intermediate phases, is quickly radiated away. A simple interaction model implies an effective warm cloud size in the fountain flow of a few 100 pc, showing that warm–hot flux exchange mainly involves a few large clouds rather than many small ones.

We study the effects of mass and energy injection due to OB associations spread across the rotating disc of a Milky Way-type galaxy, with the help of 3D hydrodynamic simulations. We compare the resulting X-ray emission with that produced from the injection of mass and energy from a central region. We find that the predicted X-ray image shows a filamentary structure that arises even in the absence of disc gas inhomogeneity. This structure stems from warm clumps made of disc material being lifted by the injected gas. We show that as much as half of the total X-ray emission comes from regions surrounding warm clumps that are made of a mix of disk and injected gas. This scenario has the potential to explain the origin of the observed extra-planar X-ray emission around star forming galaxies and can be used to understand the observed sublinear relation between the L X , the total X-ray luminosity, and SFR. We quantify the mass contained in these 'bow-shock' regions. We also show that the top-most region of the outer shock above the central area emits harder X-rays than the rest. Further, we find that the mass distribution in different temperature ranges is bimodal, peaking at 10 4 -10 5 K (in warm clumps) and 10 6 -10 7 K (X-ray emitting gas). The mass loading factor is found to decrease with increasing SFR, consistent with previous theoretical estimates and simulations.

We study the synchrotron radio emission from extra-planar regions of star forming galaxies. We use ideal magneto-hydrodynamical (MHD) simulations of a rotating Milky Way-type disk galaxy with distributed star formation sites for three star formation rates (SFRs) (0.3, 3, 30 M yr −1 ). From our simulations, we see emergence of galactic-scale magnetised outflows, carrying gas from the disk. We compare the morphology of the outflowing gas with hydrodynamic (HD) simulations. We look at the spatial distribution of magnetic field in the outflows. Assuming that a certain fraction of gas energy density is converted into cosmic ray energy density, and using information about the magnetic field, we obtain synchrotron emissivity throughout the simulation domain. We generate the surface brightness maps at a frequency of 1.4 GHz. The outflows are more extended in the vertical direction than radial and hence have an oblate shape. We further find that the matter right behind the outer shock, shines brighter in these maps than that above or below. To understand whether this feature can be observed, we produce vertical intensity profiles. We convolve the vertical intensity profile with the typical beam sizes of radio telescopes, for a galaxy located at 10 Mpc (similar to NGC 891) in order to estimate the radio scale height to compare with observations. We find that for our SFRs this feature will lie below the RMS noise limit of instruments. The radio scale height is found to be ∼ 300 − 1200 pc , depending on the resolution of the telescope. We relate the advection speed of the outer shock with the surface density of star formation as v adv ∝ Σ 0.3 SFR , which is consistent with earlier observations and analytical estimates.

We present low frequency observations at 315 and 745 MHz from the upgraded Giant Metrewave Radio Telescope (uGMRT) of the edge-on, nearby galaxy NGC 4631. We compare the observed surface brightness profiles along the minor axis of the galaxy with those obtained from hydrodynamical simulations of galactic outflows. These are 3D simulations that replicate star-formation in a Milky-Way mass galaxy and follow magnetized outflows emerging from the disk. We detect a plateau-like feature in the observed emission at a height of 2 − 3 kpc from the mid-plane of the galaxy, in qualitative agreement with that expected from simulations. This feature is believed to be due to the compression of magnetic fields behind the outer shocks of galactic outflows. We model the observed surface brightness profiles by assuming an exponential as well as a Gaussian fitting model. Using χ2 statistics, we find that the exponential model fits the profiles better and we use it to determine the scale heights. We estimate the scale height for the synchrotron radio emission to be ∼1 kpc. The timescales for advection due to outflows and diffusion of cosmic ray electrons are ≳ 5 and ∼160 Myr, respectively. Because advection acts on a timescale much shorter than diffusion, we conclude that in NGC 4631 advection, rather than diffusion, plays the dominant role in the formation of radio halo. The spectral index image with regions of flatter radio spectral index in the halo appears to indicate possible effects of gas outflow from the plane of the galaxy.

The hot component of the circum-galactic medium (CGM) around star forming galaxies is detected as diffuse X-ray emission. The X-ray spectra from the CGM depend on the temperature and metallicity of the emitting plasma, providing important information about the feeding and feedback of the galaxy. The observed spectra are commonly fitted using simple 1-Temperature (1-T) or 2-T models. However, the actual temperature distribution of the gas can be complex because of the interaction between galactic outflows and halo gas. Here we demonstrate this by analysing 3-D hydrodynamical simulations of the CGM with a realistic outflow model. We investigate the physical properties of the simulated hot CGM, which shows a broad distribution in density, temperature, and metallicity. By constructing and fitting the simulated spectra, we show that, while the 1-T and 2-T models are able to fit the synthesized spectra reasonably well, the inferred temperature(s) does not bear much physical meaning. Instead, we propose a log-normal distribution as a more physical model. The log-normal model better fits the simulated spectra while reproducing the gas temperature distribution. We also show that when the star formation rate is high, the spectra inside the bi-conical outflows are distinct from that outside, as outflows are generally hotter and more metal-enriched. Finally, we produce mock spectra for future missions with the eV-level spectral resolution, such as Athena, Lynx, HUBS and XRISM.

No abstract

Interstellar (ISM) and circumgalactic mediums (CGM) around galaxies are linked to several physical processes that drive galaxy evolution. For example, the X-ray emission from the CGM gas around ellipticals has been linked to the AGN feedback occurring in the host. Upcoming telescopes, such as HUBS with ∼2 eV resolution, can provide us with deep insights about the hot gas properties of such galaxies thus constrain these processes. In this project, we discuss X-ray emission of the ISM and CGM of elliptical galaxies simulated using MACER code. We generate X-ray emission data from the MACER simulations with various feedback models and produce mock observations for an instrument with high spectral resolution, which is a necessary step of selecting sources for the future observations with planned mission such as HUBS. More importantly, we establish connections between the physics of AGN and stellar feedback with the emission spectra from the ISM and CGM to investigate the possibility of using observations to constrain feedback models. We fit the X-ray spectra from these simulations with standard fitting procedures and compare the retrieved physical properties with their counterparts from the simulations to understand whether the future high-resolution observations can reliably reveal the properties of the gas in the galaxies.