Search

The circumgalactic medium (CGM) is the largest baryon reservoir around galaxies, but its extent, mass, and temperature distribution remain uncertain. We propose that cold gas in the CGM resides primarily in $\sim 100 \hbox{--} 10^4$ cloud complexes (CCs), each containing a mist of tiny cold cloudlets dispersed in a warm/hot medium ($\sim 10^5 \hbox{--} 10^6$~K). Modeling CCs as uniform and misty simplifies the calculation of observables like ion absorption columns compared to resolving tiny individual cloudlets. Using Monte Carlo simulations, we explore how CC properties affect the observed spread in column densities. A power-law distribution of CCs ($dN_{\rm CC}/dR \propto R^{-1.2}$) reproduces MgII equivalent width (EW) and column density distribution trends with impact parameter ($R_\perp$). We show that the area-averaged MgII column density, combined with the covering fraction of CCs, provides a robust proxy for estimating the cold CGM mass, independent of other model parameters. Modeling individual CCs demonstrates that turbulent broadening blends cloudlet absorption lines, allowing CCs to approximate the observational effects of their constituent cloudlets analytically. Direct simulations of cloudlets within multiple CCs confirm the computational challenges of fully resolving the mist-like structure. Comparing modeled MgII absorption with observations, we estimate the cold CGM mass of Milky Way-like galaxies to be $\sim 10^{10} \, M_\odot$, about $10\%$ of the total CGM mass. This work provides a practical framework for connecting CGM models with observations, shedding light on the cold gas distribution in galaxy halos and its role in the galactic baryon cycle., Comment: submitted to MNRAS. 22 pages, 2 tables, 18 figures including Appendix

The recent detection of large column density absorption lines from highly ionized gas in a few directions through the circumgalactic medium (CGM) of the Milky Way (MW) has been puzzling. The inferred temperature from these absorption lines far exceeds the virial temperature of the MW, and the column densities are also too large to be easily explained. In this paper, we propose a novel idea to explain these observations and claim that they may not have originated from the CGM, but from a totally different type of source, namely, stellar ejecta from supernovae (SNe) above the Galactic disk that happen to lie in the line of sight to the background quasars. About $\sim 20\%$ of massive OB stars (progenitors of core-collapse supernovae) are known to be runaway stars that have high ejection velocities near the Galactic plane and can end up exploding as SNe above the Galactic disk. We show that the associated reverse shock in the supernova remnant in the early non-radiative phase can heat the ejecta to temperatures of $\gtrsim 10^7\,{\rm K}$ and can naturally explain the observed high column density of ions in the observed `super-virial' phase along with $\alpha$-enriched super-solar abundance that is typical of core-collapse supernovae. However, SNe from runaway stars has a covering fraction of $\lesssim 0.7 \%$ and thus can only explain the observations along limited sightlines., Comment: 14 pages, accepted for publication in the Astrophysical Journal (ApJ)

The presence of the $\approx 10^6$ K gas in the circumgalactic medium of the Milky Way has been well established. However, the location and the origin of the newly discovered hot gas at `super-virial' temperatures of $\approx 10^7$ K have been puzzling. This hot gas has been detected in both absorption and emission; here we focus on the emitting gas only. We show that both the `virial' and the `super-virial' temperature gas as observed in \emph{emission} occupy disk-like extraplanar regions, in addition to the diffuse virial temperature gas filling the halo of the Milky Way. We perform idealized hydrodynamical simulations to show that the $\approx 10^7$ K emitting gas is likely to be produced by stellar feedback in and around the Galactic disk. We further show that the emitting gas at both super-virial and virial temperatures in the extraplanar regions is metal enriched and is not in hydrostatic equilibrium with the halo but is continuously evolving., Comment: 11 pages, accepted for publication in the Astrophysical Journal (ApJ)

We show that massive young star clusters may be possible candidates that can accelerate Galactic cosmic rays (CRs) in the range of $10^7\hbox{--}10^9$ GeV (between the `knee' and `ankle'). Various plausible scenarios such as acceleration at the wind termination shock (WTS), supernova shocks inside these young star clusters, etc. have been proposed,since it is difficult to accelerate particles up to the $10^7\hbox{--}10^9$ GeV range in the standard paradigm of CR acceleration in supernova remnants. We consider a model for the production of different nuclei in CRs from massive stellar winds using the observed distribution of young star clusters in the Galactic plane. We present a detailed calculation of CR transport in the Galaxy, taking into account the effect of diffusion, interaction losses during propagation, and particle re-acceleration by old supernova remnants to determine the all-particle CR spectrum. Using the maximum energy estimate from the Hillas criterion, we argue that a young massive star cluster can accelerate protons up to a few tens of PeV. Upon comparison with the observed data, our model requires a CR source spectrum with an exponential cutoff of $5\times 10^7 Z$ GeV ($50\,Z$~PeV) from these clusters together with a cosmic-ray injection fraction of $\sim 5\%$ of the wind kinetic energy. We discuss the possibility of achieving these requirements in star clusters, and the associated uncertainties, in the context of considering star clusters as the natural accelerator of the `second component' of Galactic cosmic rays., Comment: 18 pages, 6 figures, accepted for publication in ApJ

Recent observations and simulations reveal that the circumgalactic medium (CGM) surrounding galaxies is multiphase, with the gas temperatures spanning a wide range at most radii, $\sim 10^4\ {\rm K}$ to the virial temperature ($\sim 10^6$ K for Milky Way). Traditional CGM models using simple density profiles are inadequate at reproducing observations that indicate a broad temperature range. Alternatively, a model based on probability distribution functions (PDFs) with parameters motivated by simulations can better match multi-wavelength observations. In this work, we use log-normal distributions, commonly seen in the simulations of the multiphase interstellar and circumgalactic media, to model the multiphase CGM. We generalize the isothermal background model by Faerman et al. 2017 to include more general CGM profiles. We extend the existing probabilistic models from 1D-PDFs in temperature to 2D-PDFs in density-temperature phase space and constrain its parameters using a Milky Way-like {\tt Illustris TNG50-1} halo. We generate various synthetic observables such as column densities of different ions, UV/X-ray spectra, and dispersion and emission measures. X-ray and radio (Fast Radio Burst) observations mainly constrain the hot gas properties. However, interpreting cold/warm phase diagnostics is not straightforward since these phases are patchy, with inherent variability in intercepting these clouds along arbitrary lines of sight. We provide a tabulated comparison of model predictions with observations and plan to expand this into a comprehensive compilation of models and data. Our modeling provides a simple analytic framework that is useful for describing important aspects of the multiphase CGM., Comment: 23 pages, 15 figures, 4 tables; submitted to MNRAS

We consider the dynamics of and emission from growing superbubbles in a stratified interstellar gaseous disc driven by energy release from supernovae explosions in stellar clusters with {masses $M_{cl}= 10^5-1.6\times 10^6~M_\odot$}. Supernovae are spread randomly within a sphere of $r_c=60$ pc, and inject energy episodically with a specific rate $1/130~M_\odot^{-1}$ proportional to the star formation rate (SFR) in the cluster. Models are run for several values of SFR in the range $0.01$ to $0.1~M_\odot$ yr$^{-1}$, with the corresponding average surface energy input rate $\sim 0.04-0.4$ erg cm$^{-2}$ s$^{-1}$. We find that the discrete energy injection by isolated SNe are more efficient in blowing superbubbles: asymptotically they reach heights of up to 3 to 16 kpc for $M_{cl}=10^5-1.6\times 10^5~M_\odot$, correspondingly, and stay filled with a hot and dilute plasma for at least 30 Myr. During this time they emit X-ray, H$\alpha$ and dust infrared emission. X-ray liminosities $L_X\propto {\rm SFR}^{3/5}$ that we derive here are consistent with observations in star-forming galaxies. Even though dust particles of small sizes $a\leq 0.03~\mu$m are sputtered in the interior of bubbles, larger grains still contribute considerably ensuring the bubble luminosity $L_{\rm IR}/{\rm SFR}\sim 5\times 10^7 L_\odot M_\odot^{-1} ~{\rm yr}$. It is shown that the origin of the North Polar Spur in the Milky Way can be connected with activity of a cluster with the stellar mass of $\sim 10^5~M_\odot$ and the ${\rm SFR}\sim 0.1~M_\odot$ yr$^{-1}$ some 25--30 Myr ago. Extended luminous haloes observed in edge-on galaxies (NGC 891 as an example) can be maintained by disc spread stellar clusters of smaller masses $M_\ast \simlt 10^5~M_\odot$., Comment: 13 pages, 18 figures, accepted to MNRAS

One of the most distant galaxies GN-z11 was formed when the Universe was $\le$ 400 Myr old, and it displays a burst-like star formation rate $\sim 25~M_\odot$ yr$^{-1}$ with a metallicity $Z\sim 0.2\pm 0.1Z_\odot$. It resembles $z=2-3$ galaxies (at ``cosmic noon") except for the fact that the measured reddening $E(B-V)=0.01\pm 0.01$ indicates the presence of little or no dust. This marked absence of dust hints towards violent dynamical events that destroy or evacuate dust along with gas out of the galaxy on a relatively short time scale and make it transparent. We apply a 3D numerical model to infer possible physical characteristics of these events. We demonstrate that the energetics of the observed star formation rate is sufficient to tear apart the dusty veil on time scales of $20-25$ Myr. This can explain the apparent lack of evolution of UV luminosity function of galaxies between and $z\ge 10$ and $z\sim 7$, by compensating for the lower galaxy masses at higher redshift by the absence of dust. We show, however, that this is a temporary phenomenon and soon after the last of the supernovae explosions have taken place, the expanding shell shrinks and obscures the galaxy on time scales of $\approx 5-8$ Myr., Comment: accepted to MNRAS, 6 pages, 3 figures

Diffuse radio emission has been detected in a considerable number of galaxy clusters and groups, revealing the presence of pervasive cosmic magnetic fields, and of relativistic particles in the large-scale structure (LSS) of the Universe. Since cluster radio emission is faint and steep spectrum, its observations are largely limited by the instrument sensitivity and frequency of observation, leading to a dearth of information, more so for lower-mass systems. The unprecedented sensitivity of recently commissioned low-frequency radio telescope arrays, aided by the development of advanced calibration and imaging techniques, have helped in achieving unparalleled image quality. At the same time, the development of sophisticated numerical simulations and the availability of supercomputing facilities have paved the way for high-resolution numerical modeling of radio emission, and the structure of the cosmic magnetic fields in LSS, leading to predictions matching the capabilities of observational facilities. In view of these rapidly-evolving scenerio in modeling and observations, in this review, we summarise the role of the new telescope arrays and the development of advanced imaging techniques and discuss the detections of various kinds of cluster radio sources. In particular, we discuss observations of the cosmic web in the form of supercluster filaments, studies of emission in poor clusters and groups of galaxies, and of ultra-steep spectrum sources. We also review the current theoretical understanding of various diffuse cluster radio sources and the associated magnetic field and polarization. As the statistics of detections improve along with our theoretical understanding, we update the source classification schemes based on their intrinsic properties. We conclude by summarising the role of the upgraded GMRT and our expectations from the upcoming Square Kilometre Array (SKA) observatories., Comment: 32 pages, 10 figures, accepted for publication in the Journal of Astrophysics and Astronomy (JoAA) (to appear in the special issue on "Indian participation in the SKA")

We discuss the production of $\gamma$-rays from cosmic rays (CR) in the circumgalactic medium (CGM) of Andromeda (M31) in light of the recent detection of $\gamma$-rays from an annular region of $\sim 5.5-120$ kpc away from the M31 disc. We consider the CRs accelerated as a result of the star-formation in the M31 disk, which are lifted to the CGM by advection due to outflow and CR diffusion. The advection time scale due to bulk flow of gas triggered by star formation activity in the M31 disc is comparable ($\sim$ Gyr) to the diffusion time scale with diffusion coefficient $\ge10^{29}$ cm$^2$ s$^{-1}$ for the propagation of CR protons with energy $\sim 412$ GeV that are responsible for the highest energy photons observed. We show that a leptonic origin of the $\gamma$-rays from cosmic ray (CR) electrons has difficulties, as the inverse Compton time scale ($\sim$Myr) is much lower than advection time scale ($\sim$Gyr) to reach $120$ kpc. Invoking CR electrons accelerated by accretion shocks in the CGM at $\sim100-120$ kpc does not help since it would lead to diffuse X-ray features that are not observed. We, therefore, study the production of $\gamma$-rays via hadronic interaction between CR protons and CGM gas with the help of numerical two-fluid (thermal + CR) hydrodynamical simulation. We find that a combination of these mechanisms, that are related to the star formation processes in M31 in the last $\sim $ Gyr, along with diffusion and hadronic interaction, can explain the observed flux from the CGM of M31., Comment: 10 pages, 5 figures, Accepted for publication in MNRAS on May 19, 2022

We constrain the cosmic ray (CR) population in the circumgalactic medium (CGM) of Milky Way by comparing the observations of absorption lines of OVIII ion with predictions from analytical models of CGM : precipitation (PP) and isothermal (IT) model. For a CGM in hydrostatic equilibrium, the introduction of CR suppresses thermal pressure, and affects the OVIII ion abundance. We explore the allowances given to the ratio of CR pressure to thermal pressure ($\rm{P}_{\rm{CR}}/\rm{P}_{\rm{th}}=\eta$), with varying boundary conditions, CGM mass content, photoionization by extragalactic ultraviolet background and temperature fluctuations. We find that the allowed maximum values of $\eta$ are : $\eta\lesssim10$ in the PP model and $\eta\lesssim6$ in the IT model. We also explore the spatial variation of $\eta$ : rising ($\eta=Ax$) or declining ($\eta=A/x$) with radius, where A is the normalization of the profiles. In particular, the models with declining ratio of CR to thermal pressure fare better than those with rising ratio with suitable temperature fluctuation (larger $\sigma_{\rm ln T}$ for PP and lower for IT). The declining profiles allow $A\lesssim8$ and $A\lesssim10$ in the case of IT and PP models, respectively, thereby accommodating a large value of $\eta \,(\simeq 200)$ in the central region, but not in the outer regions. These limits, combined with the limits derived from $\gamma$-ray and radio background, can be useful for building models of Milky Way CGM including CR population. However, the larger amount of CR can be packed in cold phase which may be one way to circumvent these constraints., Comment: 11 pages, 5 figures, 1 table Accepted for publication in ApJ on Apr 22, 2022

Analysing a large representative sample of local galaxies (8707), we find that the variation in the shape of their dust attenuation curves is driven primarily by their structure, i.e., distribution of stars (and dust) within them. The attenuation curve for spheroid dominated galaxies, as compared to the disc dominated ones, is nearly twice as steep. Both structural types cover distinct ranges of attenuation slope values. Similar findings are reflected in the case of star-forming and passive galaxies. Spheroids and passive galaxies witness minimal attenuation in the optical compared to UV wavelengths underlining the lack of dusty birth-clouds that define complex star-dust geometry. The distinction in the attenuation properties of spheroids and discs is maintained in each stellar mass range emphasising that structure is the primal cause of variation. However, within a structural group, the attenuation curve becomes shallower with both the increase in total stellar mass and optical depth of the galaxy. Overall, with the extinction curve fixed to be the same for all galaxies, the star-dust geometry emerges to be the prime determinant of the variation in their attenuation properties., Comment: Accepted for publication in MNRAS Letters, 5 pages, 3 figures

The presence of dust in the intracluster medium (ICM) has been a long-standing problem that is still awaiting elucidation. Direct observational diagnostics are rather challenging (though not impossible) either because of a sparse distribution of dust in the intracluster space that makes extinction measurements difficult or because of a low surface brightness of infrared emission from dust. Complex indirect approaches are currently available that can overcome uncertainties and provide a reasonable understanding of the basic regulations of the physical state of dust in the ICM. Contrary to the common opinion that the hot ICM does not allow dust to survive and manifest, many sparse observational data either directly point out that dust exists in the intracluster space or its presence is consistent with the data. Highly divergent data in direct evidence and highly uncertain indirect indicators are often connected either with dust fragility in a hot environment, the possible compactness of spatial (clumpy) dust distribution in the ICM, or dynamical features of dust transport. The source of dust is obviously connected with galaxies, and it turns out that in most cases, dust is carried from galaxies into the ICM while being thermally and dynamically shielded against the hostile influence of high-energy ions. In this review, we briefly discuss related issues from observational and theoretical points of view, including the transport of dust into the ICM, and the associated destructive and protective mechanisms and their characteristic time scales.

Active galactic nuclei (AGN) powered by the central Super-Massive Black Holes (SMBHs) play a major role in modifying the thermal properties of the intracluster medium (ICM). In this work, we implement two AGN heating models: (i) by buoyant cavities rising through stratified ICM (effervescent model) and, (ii) by viscous and conductive dissipation of sound waves (acoustic model). Our aim is to determine whether these heating models are consistent with ICM observables and if one is preferred over the other. We assume an initial entropy profile of ICM that is expected from the purely gravitational infall of the gas in the potential of the dark matter halo. We then incorporate heating, radiative cooling, and thermal conduction to study the evolution of ICM over the age of the clusters. Our results are: (i) Both the heating processes can produce comparable thermal profiles of the ICM with some tuning of relevant parameters. (ii) Thermal conduction is crucially important, even at the level of 10\% of the Spitzer values, in transferring the injected energy beyond the central regions, and without which the temperature/entropy profiles are unrealistically high. (iii) The required injected AGN power scales with cluster mass as $M_{\rm vir}^{1.5}$ for both models. (iv) The required AGN luminosity is comparable with the observed radio jet power, reinforcing the idea that AGNs are the dominant heating source in clusters. (v) Finally, we estimate that the fraction of the total AGN luminosity available as the AGN mechanical luminosity at $0.02r_{500}$ is less than 0.05\%., Comment: 13 pages, 5 figures, Accepted in MNRAS

Current and upcoming high angular resolution and multi-frequency experiments are well poised to explore the rich landscape of secondary CMB anisotropies. In this context, we compute for the first time, the power spectrum of CMB fluctuations from a cosmological distribution of evolving lobes of giant radio galaxies. We, also, explicitly take into account the non-thermal electron distribution, which has important implications for the inference of the CMB angular power spectrum. We calculate the mean global non-thermal y-distortion, \ynt. For observationally reasonable distribution of the jet luminosities in the range of $10^{45}-10^{47}$ ergs$^{-1}$, we find \ynt to be less than $10^{-5}$, and hence not violating the COBE limit as previously claimed. Using the unique spectral dependence of the non-thermal SZ, we show that a detection of \ynt can be within reach at the level of $\gtrsim 5\sigma$ from a future PIXIE-like experiment provided we understand the foregrounds precisely. The total non-thermal SZ power spectrum, $C^{NT}_\ell$, from the radio lobes peaks at $\ell \sim 3000$ with an amplitude $\sim 1\%$ of thermal SZ power spectrum from galaxy clusters. A a detection of the $C^{NT}_\ell$, with a PIXIE-like sensitivity experiment, can lead to $\sim 5\sigma$ constraint on the mass dependence of the jet luminosity with the constraint becoming, at least, ten times better for the proposed more ambitious CMB-HD survey. This will, further, lead to the tightest constrain on the central black hole mass -to- host halo mass scaling relations., Comment: Accepted version

We study the effect of cosmic ray (CR) acceleration in the massive compact star cluster Westerlund 1 in light of its recent detection in $\gamma$-rays. Recent observations reveal a $1/r$ radial distribution of the CR energy density. Here we theoretically investigate whether or not this profile can help to distinguish between (1) continuous CR acceleration in the star cluster stellar wind-driven shocks and (2) discrete CR acceleration in multiple supernovae shocks -- which are often debated in the literature. Using idealized two-fluid simulations and exploring different acceleration sites and diffusion coefficients, we obtain the CR energy density profile and luminosity to find the best match for the $\gamma$-ray observations. We find that the inferred CR energy density profiles from observations of $\gamma$-ray luminosity and mass can be much different from the true radial profile. CR acceleration at either the cluster core region or the wind termination shock can explain the observations, if the diffusion coefficient is $\kappa_{\rm cr}\sim 10^{27}$ cm$^2$ s$^{-1}$ and a fraction of $\approx 10\%-20\%$ of the shock power/post-shock pressure is deposited into the CR component. We also study the possibility of discrete supernovae (SN) explosions being responsible for CR acceleration and find that with an injection rate of 1 SN in every $\sim 0.03$ Myr, one can explain the observed $\gamma$-ray profile. This multiple SN scenario is consistent with X-ray observations only if the thermal conductivity is close to the Spitzer value., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

We consider a model of the circumgalactic medium (CGM) in which feedback maintains a constant ratio of cooling time to freefall time throughout the halo, so that the entire CGM is marginally unstable to multiphase condensation. This 'precipitation model' is motivated by observations of multiphase gas in the cores of galaxy clusters and the halos of massive ellipticals. We derive from the model density and temperature profiles for the CGM around galaxies with masses similar to the Milky Way. After taking into consideration the geometrical position of our solar system in the Milky Way, we show that the CGM model is consistent with observed OVI, OVII, and OVIII column densities only if temperature fluctuations with a log-normal dispersion $\sigma_{\ln T} \sim 0.6$-$1.0$ are included. We show that OVI column densities observed around star-forming galaxies require systematically greater values of $\sigma_{\ln T}$ than around passive galaxies, implying a connection between star formation in the disk and the state of the CGM. Photoionization by an extra-galactic UV background does not significantly change these CGM features for galaxies like the Milky Way but has much greater and significant effects on the CGM of lower-mass galaxies., Comment: 12 pages, 10 figures, 1 table, Accepted for publication in MNRAS

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 $\chi^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 $\sim 1$ kpc. The timescales for advection due to outflows and diffusion of cosmic ray electrons are $\gtrsim 5$ and $\sim 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., Comment: 9 Pages, 9 Figures, accepted for publication in MNRAS

Energetic electrons in the lobes of radio galaxies make them potential sources for not only radio and X-rays but also Sunyaev-Zeldovich (SZ) distortions in the cosmic microwave background (CMB) radiation. Previous works have discussed the energetics of radio galaxy lobes, but assuming thermal SZ effect, coming from the non-thermal electron population. We use an improved evolutionary model for radio galaxy lobes to estimate the observed parameters such as the radio luminosity and intensity of SZ-distortions at the redshifts of observation. We, further, quantify the effects of various relevant physical parameters of the radio galaxies, such as the jet power, the time scale over which the jet is active, the evolutionary time scale for the lobe, etc on the observed parameters. For current SZ observations towards galaxy clusters, we find that the non-thermal SZ distortions from radio lobes embedded in galaxy clusters can be non-negligible compared to the amount of thermal SZ distortion from the intra-cluster medium and, hence, can not be neglected. We show that small and young (and preferably residing in a cluster environment) radio galaxies offer better prospects for the detection of the non-thermal SZ signal from these sources. We further discuss the limits on different physical parameters for some sources for which SZ effect has been either detected or upper limits are available. The evolutionary models enable us to obtain limits, previously unavailable, on the low energy cut-off of electron spectrum ($p_{min} \sim 1\hbox{--}2$) in order to explain the recent non-thermal SZ detection. Finally, we discuss how future CMB experiments, which would cover higher frequency bands ($>$400 GHz), may provide clear signatures for non-thermal SZ effect., Comment: Replaced with published version

We demonstrate that young star clusters have a $\gamma$-ray surface brightness comparable to that of the diffuse Galactic emission (DGE), and estimate that their sky coverage in the direction of the inner Galaxy exceeds unity. We therefore suggest that they comprise a significant fraction of the DGE., Comment: 5 pages, 1 figure, Accepted for publication in MNRAS Letters

We present constraints on a simple analytical model for hot diffuse halo gas, derived from a fit spanning two orders of magnitude in halo mass ($M_{500} \sim 10^{12.5}-10^{14.5} M_{\odot}$). The model is motivated by the observed prevalence of a precipitation limit, and its main free parameter is the central ratio of gas cooling timescale to free-fall timescale ($t_{\rm cool}/t_{\rm ff}$). We use integrated X-ray and thermal Sunyaev-Zel'dovich observations of the environments around massive galaxies, galaxy groups and clusters, averaged in halo mass bins, and obtain the best-fitting model parameters. We find $t_{\rm cool}/t_{\rm ff} \sim 50-110$, depending on the model extrapolation beyond the halo virial radius and possibly on biases present in the data-sets used in the fitting analysis. The model adequately describes the entire mass range, except for intermediate mass halos ($M_{500} \sim 10^{13.5} M_{\odot}$) which systematically fall below the model predictions. However, the best fits for $t_{\rm cool}/t_{\rm ff}$ substantially exceed the values typically derived from X-ray observations of individual systems ($t_{\rm cool}/t_{\rm ff} \sim 10-30$). We consider several explanations for those discrepancies, including X-ray selection biases and a potential anti-correlation between X-ray luminosity and the central galaxy's stellar mass., Comment: 13 pages, 5 figures, 1 table. Accepted for publication in MNRAS

We study the interaction of cosmic rays (CRs) with the diffuse circumgalactic gas of Milky Way (MW) galaxy that results in hadronic $\gamma-$ray emission and radio synchrotron emission. We aim to constrain the CR population in our circumgalactic medium (CGM) with the help of observed isotropic $\gamma$-ray background (IGRB), its anisotropy and radio continuum. We modify different models of CGM gas in hydrostatic equilibrium discussed in literature by including a cosmic ray population, parametrized by $\eta \equiv P_{\rm CR}/P_{\rm th}$. For the simplest isothermal model, while the IGRB intensity allows $\eta \lesssim 3$, the anisotropy resulting from the Solar system's off-center position in MW rules out all values of $\eta$. For the precipitation model, in which the cooling of the CGM gas is regulated with an optimum ratio of cooling time to free-fall time,while the observed IGRB intensity allows $\eta \lesssim 230$, the observed anisotropy allows only very large values of $\eta$, of order $\gtrsim 100$. The radio continuum limits $\eta \lesssim 400$ for precipitation model and does not constrain isothermal model, however these constraints are mitigated by synchrotron loss time being comparable to CR diffusion time scales. These bounds are relevant for current numerical simulations that indicate a significant CR population in CGM of galaxies of MW mass., Comment: 8 pages, 5 figures, Accepted for publication in ApJL

Using an idealized set-up, we investigate the dynamical role of cosmic rays (CRs) in the early stages of galactic outflows for galaxies of halo masses $10^{8}$, $10^{11}$ and $10^{12}$ $M_\odot$. The outflow is launched from a central region in the galactic disk where we consider three different constant star formation rates ($0.1$, $1$, and $10$ $M_{\odot}$ yr$^{-1}$) over a dynamical timescale of $50$ Myr. We determine the temperature distribution of the gas and find that CRs can reduce the temperature of the shocked gas, which is consistent with previous results. However, we show that CRs do not have any noticeable effect on the mass loading by the outflow. We find that CRs can reduce the size of the outflow, which contradicts previous claims of efficient dynamical impact of CRs; however, it is consistent with earlier theoretical models of cosmic ray driven blastwave as well as stellar wind. We discuss the dependence of our results on CR injection prescriptions and compare them with earlier studies. We conclude that in the early stages of galactic outflows the dynamical role of CRs is not important., Comment: 19 pages, 2 tables, 19 figures, accepted for publication in MNRAS, Link to time evolution of the simulation performed : https://youtu.be/3-JZAv7vwo8

We consider the size distribution of superbubbles in a star forming galaxy. Previous studies have tried to explain the distribution by using adiabatic self-similar evolution of wind driven bubbles, assuming that bubbles stall when pressure equilibrium is reached. We show, with the help of hydrodynamical numerical simulations, that this assumption is not valid. We also include radiative cooling of shells. In order to take into account non-thermal pressure in the ambient medium, we assume an equivalent higher temperature than implied by thermal pressure alone. Assuming that bubbles stall when the outer shock speed becomes comparable to the ambient sound speed (which includes non-thermal components), we recover the size distribution with a slope of $\sim -2.7$ for typical values of ISM pressure in Milky Way, which is consistent with observations. Our simulations also allow us to follow the evolution of size distribution in the case of different values of non-thermal pressure, and we show that the size distribution steepens with lower pressure, to slopes intermediate between only-growing and only-stalled cases., Comment: 12 pages, 16 figures, Published in MNRAS

Supermassive black holes (BHs) and their host galaxies are interlinked by virtue of feedbacks and are thought to be co-eval across the Hubble time. This relation is highlighted by an approximate proportionality between the BH mass $M_\bullet$ and the mass of a stellar bulge $M_\ast$ of the host galaxy. However, a large spread of the ratio $M_\bullet/M_\ast$ and a considerable excess of BH mass at redshifts $z\sim8$, indicate that the coevolution of central massive BHs and stellar populations in host galaxies may have experienced variations in its intensity. These issues require a robust determination of the relevant masses (BH, stars and gas), which is difficult in the case of distant high-redshift galaxies that are unresolved. In this paper, we seek to identify spectral diagnostics that may tell us about the relative masses of the BH, the gas mass and stellar mass. We consider general features of SEDs of galaxies that harbour growing massive BHs, forming stars and interstellar/circumgalactic gas. We focus on observational manifestations of possible predominances or intermittent variations in evolutionary episodes of growing massive BHs and forming stellar populations. We consider simplified scenarios for star formation and massive BHs growth, and simple models for chemical composition of gas, for dust free gas as well as for gas with dust mass fraction of $1/3$ of the metal content. We argue that wideband multi-frequency observations (X-ray to submillimeter) of the composite emission spectra of growing BH, stellar population and nebular emission of interstellar gas are sufficient to infer their masses., Comment: 15 pages, 13 figures, submitted to MNRAS, comments welcome

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$_{\odot}$ 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 $\sim 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 $\rm{v}_{\rm adv} \propto \Sigma_{\rm SFR}^{0.3}$ which is consistent with earlier observations and analytical estimates., Comment: 12 pages, 12 figures, Under Review with MNRAS

Cosmic ray (CR) sources leave signatures in the isotopic abundances of CRs. Current models of Galactic CRs that consider supernovae (SNe) shocks as the main sites of particle acceleration cannot satisfactorily explain the higher ${\rm ^{22}Ne/^{20}Ne}$ ratio in CRs compared to the interstellar medium. Although stellar winds from massive stars have been invoked, their contribution relative to SNe ejecta has been taken as a free parameter. Here we present a theoretical calculation of the relative contributions of wind termination shocks (WTSs) and SNe shocks in superbubbles, based on the hydrodynamics of winds in clusters, the standard stellar mass function, and stellar evolution theory. We find that the contribution of WTSs towards the total CR production is at least $25\%$, which rises to $\gtrsim 50\%$ for young ($\lesssim 10$ Myr) clusters, and explains the observed $^{22}{\rm Ne}/^{20} {\rm Ne}$ ratio. We argue that since the progenitors of apparently isolated supernovae remnants (SNRs) are born in massive star clusters, both WTS and SNe shocks can be integrated into a combined scenario of CRs being accelerated in massive clusters. This scenario is consistent with the observed ratio of SNRs to $\gamma$-ray bright ($L_\gamma \gtrsim 10^{35}$ erg s$^{-1}$) star clusters, as predicted by star cluster mass function. Moreover, WTSs can accelerate CRs to PeV energies, and solve other longstanding problems of the standard supernova paradigm of CR acceleration., Comment: 19 pages, 18 figures, 1 table; Submitted to MNRAS; Comments welcome!

We have investigated the long term average spectral properties of galactic X-ray binaries in the energy range of 3-20 keV, using long term monitoring data from MAXI -Gas Slit Camera (GSC). These long term average spectra are used to construct separately the composite spectra of galactic High Mass X-ray binaries (HMXBs) and Low Mass X-ray binaries (LMXBs). These composite spectra can be described empirically with piece-wise power-law with three components. X-rays from HMXBs are considered as important contributors to heating and ionization of neutral hydrogen in the intergalactic medium during the Epoch of Reionization. Using the above empirical form of the composite HMXB spectra extrapolated to lower energies as an input, we have studied the impact of these sources on the 21-cm signal using the outputs of N-body simulation and 1D radiative transfer. The heating due to the composite spectrum is less patchy compared to power-law spectrum with a spectral index $\alpha = 1.5$, used in previous studies. The amplitude of the heating peak of large scale power spectrum, when plotted as a function of the redshift, is less for the composite spectrum., Comment: 13 pages, 6 figures. Accepted for publication in MNRAS

Galactic gas outflows are driven by stellar feedback with dominant contribution from supernovae (SN) explosions. The question of whether the energy deposited by SNe initiates a large scale outflow or gas circulation on smaller scales -- between discs and intermediate haloes, depends on SN rate and their distribution in space and time. We consider here gas circulation by disc-wide unclustered SNe with galactic star formation rate in the range from {$\simeq 6\times 10^{-4}$ to $\simeq 6\times 10^{-2}~M_\odot$~yr$^{-1}$~kpc$^{-2}$, corresponding to mid-to-high star formation observed in galaxies. We show that such disc-wide SN explosion regime can form circulation of warm ($T\sim 10^4$ K) and cold ($T<10^3$ K) phases within a few gas scale heights, and elevation of hot {($T>10^5$ K)} gas at higher ($z>1$ kpc) heights. We found that the threshold energy input rate for hot gas outflows with disc-wide supernovae explosions is estimated to be of the order $\sim 4\times 10^{-4}$~erg~s$^{-1}$~cm$^{-2}$. We discuss the observational manifestations of such phenomena in optical and X-ray bands. In particular, we found that for face-on galaxies with SF ($\Sigma_{_{\rm SF}}>0.02~M_\odot$~yr$^{-1}$~kpc$^{-2}$), the line profiles of ions typical for warm gas show a double-peak shape, corresponding to out-of-plane outflows. In the X-ray bands, galaxies with high SF rates ($\Sigma_{_{\rm SF}}>0.006~M_\odot$~yr$^{-1}$~kpc$^{-2}$) can be bright, with a smooth surface brightness in low-energy bands ($0.1\hbox{--}0.3$~keV) and patchy at higher energies ($1.6\hbox{--}8.3$~keV).}, Comment: 11 pages, 11 figures

We consider the synchrotron emission from high energy electrons accelerated in supernova explosions of massive Population III stars in high redshift minihaloes of mass $10^{5\hbox {--}7} \rm M_\odot$. We show the resulting intensity of radio background from this process can be substantial, which could potentially explain the recently reported EDGES result, if not for the associated heating of the IGM by CR protons which are also produced at the same time. The trade-off between the radio background and heating is such that the 21 cm brightness temperature cannot be larger than $\vert \Delta T_{21}\vert \sim 0.25$ K. The radio background and heating are both produced by energetic particles, although one by energetic electrons and the other by energetic protons. The two competing processes, production of radio background and heating of IGM by Pop III supernovae, determine the depth of the trough in the 21 cm brightness temperature which can be observed in future experiments and used as a test of this scenario., Comment: 6 pages, 2 figures, Accepted for publication in MNRAS

We study the heating of the intergalactic medium (IGM) surrounding high redshift star forming galaxies due to cosmic rays (CR). We take into account the diffusion of low energy cosmic rays and study the patchiness of the resulting heating. We discuss the case of IGM heating around a high redshift minihalo ($z\sim 10\hbox{--}20$, M$\sim 10^5\hbox{--}10^7$ M$_\odot$),and put an upper limit on the diffusion coefficient $D\le 1\times 10^{26}$ cm$^2$ s$^{-1}$ for the heating to be inhomogeneous at $z\sim 10$ and $D\le 5\hbox{--}6 \times 10^{26}$ cm$^2$ s$^{-1}$ at $z\sim 20$. For typical values of $D$, our results suggest uniform heating by CR at high redshift, although there are uncertainties in magnetic field and other CR parameters. We also discuss two cases with continuous star formation, one in which the star formation rate (SFR) of a galaxy is high enough to make the IGM in the vicinity photoionized, and another in which the SFR is low enough to keep it neutral but high enough to cause significant heating by cosmic ray protons. In the neutral case (low SFR), we find that the resulting heating can make the gas hotter than the cosmic microwave background (CMB) radiation for $D < 10^{30}$ cm$^2$ s$^{-1}$, within a few kpc of the galaxy, and unlikely to be probed by near future radio observations. In the case of photoionized IGM (high SFR), the resulting heating of the gas in the vicinity of high redshift ($z\sim 4)$ galaxies of mass $\ge 10^{12}$ M$_\odot$ can suppress gas infall into the galaxy. At lower redshifts ($z\sim 0$), an SFR of $\sim 1$ M$_\odot$ yr$^{-1}$ can suppress the infall into galaxies of mass $\le 10^{10}$ M$_\odot$., Comment: 10 pages, 9 figures, Accepted for publication in MNRAS

We use $1$D and $3$D two-fluid cosmic ray (CR) hydrodynamic simulations to investigate the role of CRs in the vicinity of a compact young star cluster. We model a self-gravitating cloud (density profile $\rho \propto r^{-1}$), include important thermal and non-thermal processes, and explore two different CR injection scenarios. We show that if internal shocks in the wind-driving region are the main site for CR acceleration, then the resulting $\gamma$-ray luminosity ($L_{\rm \gamma}$) can reach $\approx 5\%$ of the mechanical luminosity ($L_{\rm w}$), independent of the fraction of wind energy ($\sim 1-20\%$) injected into CRs. In contrast, if the forward/reverse shock of a bubble is the injection site then $L_{\rm \gamma}$ increases linearly with the CR injection fraction, as expected analytically. We find that the X-ray luminosity ($L_{\rm x}$) in the forward/reverse shock injection scenario is $\gtrsim 10^{-3} L_{\rm w}$, which is $\sim 10$ times larger than in the central wind-driving injection case. We predict the corresponding range of the synchrotron radio luminosity. We show how multi-wavelength observations can constrain the CR parameters. Comparing the predicted multi-wavelength luminosities with those of 30 Doradus we identify the reverse shock as the most probable CR injection site, and that thermal conduction is important. We do not find significant dynamical impact of CRs in our models., Comment: 16 pages, 14 figures, 2 tables; submitted to MNRAS, Comments welcome

We study the excess entropy and the corresponding non-gravitational feedback energy ($E_{feedback}$) in the intra-cluster medium (ICM) by considering a sample of 38 galaxy clusters using Chandra X-ray and NRAO VLA Sky Survey (NVSS)/Giant Metrewave Radio Telescope (GMRT) radio observations. We find moderate correlation of the feedback energy and brightest cluster galaxy (BCG) radio luminosity ($L_R$) with the various cluster thermal properties. We show conclusively that the active galactic nucleus (AGN) is more effective in transferring feedback energy to the ICM in less massive clusters. We find that within $0.3r_{500}$, the feedback energy correlates with cluster temperature as $E_{feedback}\propto T_{obs}^{0.98\pm0.37}$. Moreover, for radio detected BCG sample we find that BCG radio luminosity at 1.4 GHz scales with gas mass as $L_R\propto m_{g,obs}^{ 1.76\pm0.71}$ and with X-ray luminosity as $L_R\propto L_{X,obs}^{0.94\pm0.35}$. Finally, we discuss the implications of our results with regard to feedback in clusters., Comment: 6 pages, 5 figures, 1 table, Major revision, Accepted in MNRAS Letters

We use observations of stacked X-ray luminosity and Sunyaev-Zel'dovich (SZ) signal from a cosmological sample of $\sim 80,000$ and $104,000$ massive galaxies, respectively, with $ 10^{12.6}\lesssim M_{500} \lesssim 10^{13} M_{\odot}$ and mean redshift, \={z} $\sim$ 0.1 - 0.14 to constrain the hot Circumgalactic Medium (CGM) density and temperature. The X-ray luminosities constrain the density and hot CGM mass, while the SZ signal helps in breaking the density-temperature degeneracy. We consider a simple power-law density distribution ($n_e \propto r^{-3\beta}$) as well as a hydrostatic hot halo model, with the gas assumed to be isothermal in both cases. The datasets are best described by the mean hot CGM profile $\propto r^{-1.2}$, which is shallower than an NFW profile. For halo virial mass $\sim 10^{12}$ - $10^{13} M_{\odot}$, the hot CGM contains $\sim$ 20 - 30\% of galactic baryonic mass for the power-law model and 4 - 11\% for the hydrostatic halo model, within the virial radii. For the power-law model, the hot CGM profile broadly agrees with observations of the Milky Way. The mean hot CGM mass is comparable to or larger than the mass contained in other phases of the CGM for $L^*$ galaxies., Comment: Accepted for publication in MNRAS

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$ 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\hbox{-}10^5$ K (in warm clumps) and $10^6\hbox{-}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., Comment: 18 pages, 18 figures- Accepted in MNRAS. Comments and suggestions welcomed

We consider various possible scenarios to explain the recent observation of a claimed broad H$\alpha$ absorption in our Galactic halo with peak optical depth $\tau \simeq 0.01$ and equivalent width $W \simeq 0.17 \, \rm \AA$. We show that the absorbed feature cannot arise from the circumgalactic and ISM H$\alpha$ absorption. As the observed absorption feature is quite broad $\Delta\lambda \simeq 30 \, \rm \AA$, we also consider CNO lines that lie close to H$\alpha$ as possible alternatives to explain the feature. We show such lines could also not account for the observed feature. Instead, we suggest that it can arise from diffuse interstellar bands (DIBs) carriers or polyaromatic hydrocarbons (PAHs) absorption. While we identify several such lines close to the H$\alpha$ transition, we are unable to determine the molecule responsible for the observed feature, partly because of selection effects that prevents us from identifying DIBs/PAHs features close to H$\alpha$ using local observations. Deep integration on a few extragalactic sources with high spectral resolution might allow us to distinguish between different possible explanations., Comment: 5 pages, to appear in the Astrophysical Journal Letters

Using analytic methods and $1$-D two-fluid simulations, we study the effect of cosmic rays (CRs) on the dynamics of interstellar superbubbles (ISBs) driven by multiple supernovae (SNe)/stellar winds in OB associations. In addition to CR advection and diffusion, our models include thermal conduction and radiative cooling. We find that CR injection at the reverse shock or within a central wind-driving region can affect the thermal profiles of ISBs and hence their X-ray properties. Even if a small fraction ($10-20\%$) of the total mechanical power is injected into CRs, a significant fraction of the ram pressure at the reverse shock can be transferred to CRs. The energy transfer becomes efficient if (1) the reverse shock gas Mach number exceeds a critical value ($M_{\rm th}\gtrsim 12$) and (2) the CR acceleration time scale $\tau_{\rm acc}\sim \kappa_{\rm cr}/v^2$ is shorter than the dynamical time, where $\kappa_{\rm cr}$ is CR diffusion constant and $v$ is the upstream velocity. We show that CR affected bubbles can exhibit a volume averaged hot gas temperature $1-5\times10^{6}$ K, lower by a factor of $2-10$ than without CRs. Thus CRs can potentially solve the long-standing problem of the observed low ISB temperatures., Comment: 18 pages, 18 figures, 3 tables; submitted to MNRAS; Comments welcome

AGN feedback is regarded as an important non-gravitational process in galaxy clusters, providing useful constraints on large-scale structure formation. It modifies the structure and energetics of the intra-cluster medium (ICM) and hence its understanding is crucially needed in order to use clusters as high precision cosmological probes. In this context, particularly keeping in mind the upcoming high quality radio data expected from radio surveys like SKA with its higher sensitivity, high spatial and spectral resolutions, we review our current understanding of AGN feedback, its cosmological implications and the impact that SKA can have in revolutionizing our understanding of AGN feedback in large-scale structures. Recent developments regarding the AGN outbursts and its possible contribution to excess entropy in the hot atmospheres of groups and clusters, its correlation with the feedback energy in ICM, quenching of cooling flows and the possible connection between cool core clusters and radio mini-halos, are discussed. We describe current major issues regarding modeling of AGN feedback and its impact on the surrounding medium. With regard to the future of AGN feedback studies, we examine the possible breakthroughs that can be expected from SKA observations. In the context of cluster cosmology, for example, we point out the importance of SKA observations for cluster mass calibration by noting that most of $z>1$ clusters discovered by eROSITA X-ray mission can be expected to be followed up through a 1000 hour SKA-1 mid programme. Moreover, approximately $1000$ radio mini halos and $\sim 2500$ radio halos at $z<0.6$ can be potentially detected by SKA1 and SKA2 and used as tracers of galaxy clusters and determination of cluster selection function., Comment: 14 pages, 10 figures, Review article accepted in Journal of Astrophysics and Astronomy (JOAA)

We study the merging and evolution of isolated supernovae (SNe) remnants in a stellar cluster into a collective superbubble, with the help of 3-D hydrodynamic simulations. We particularly focus on the transition stage when the isolated SNe remnants gradually combine to form a superbubble. We find that when the SN rate is high ($\nu_{\rm sn}\sim 10^{-9}$ pc$^{-3}$ yr$^{-1}$), the merging phase lasts for $\sim 10^4$ yr, for $n=1\hbox{--}10$ cm$^{-3}$, and the merging phase lasts for a longer time ($\sim 0.1$ Myr or more) for lower SN rates ($\nu_{\rm sn}\le 10^{-10}$ pc$^{-3}$ yr$^{-1}$). During this transition phase, the growing superbubble is filled with dense and cool fragments of shells and most of the energy is radiated away during this merging process. After passing through the intermediate phase, the superbubble eventually settles on to a new power-law wind asymptote that is smaller than estimated in a continuous wind model. This results in a significant (more than {\it several times}) underestimation of the mechanical luminosity needed to feed the bubble. We determine the X-ray and H$\alpha$ surface brightnesses as functions of time for such merging SNe in a stellar cluster and find that clusters with high SN rate shine predominantly in soft X-rays and H$\alpha$. In particular, a low value of the volume averaged H$\alpha$ to H$\beta$ ratio and its large spread can be a good indicator of the transition phase of merging SNe., Comment: 15 pages, 15 figures, accepted to MNRAS

We estimate the "non-gravitational" entropy-injection profiles, $\Delta K$, and the resultant energy feedback profiles, $\Delta E$, of the intracluster medium for 17 clusters using their Planck SZ and ROSAT X-Ray observations, spanning a large radial range from $0.2r_{500}$ up to $r_{200}$. The feedback profiles are estimated by comparing the observed entropy, at fixed gas mass shells, with theoretical entropy profiles predicted from non-radiative hydrodynamic simulations. We include non-thermal pressure and gas clumping in our analysis. The inclusion of non-thermal pressure and clumping results in changing the estimates for $r_{500}$ and $r_{200}$ by 10\%-20\%. When clumpiness is not considered it leads to an under-estimation of $\Delta K\approx300$ keV cm$^2$ at $r_{500}$ and $\Delta K\approx1100$ keV cm$^2$ at $r_{200}$. On the other hand, neglecting non-thermal pressure results in an over-estimation of $\Delta K\approx 100$ keV cm$^2$ at $r_{500}$ and under-estimation of $\Delta K\approx450$ keV cm$^2$ at $r_{200}$. For the estimated feedback energy, we find that ignoring clumping leads to an under-estimation of energy per particle $\Delta E\approx1$ keV at $r_{500}$ and $\Delta E\approx1.5$ keV at $r_{200}$. Similarly, neglect of the non-thermal pressure results in an over-estimation of $\Delta E\approx0.5$ keV at $r_{500}$ and under-estimation of $\Delta E\approx0.25$ keV at $r_{200}$. We find entropy floor of $\Delta K\approx300$ keV cm$^2$ is ruled out at $\approx3\sigma$ throughout the entire radial range and $\Delta E\approx1$ keV at more than 3$\sigma$ beyond $r_{500}$, strongly constraining ICM pre-heating scenarios. We also demonstrate robustness of results w.r.t sample selection, X-Ray analysis procedures, entropy modeling etc., Comment: 17 pages, 15 figures, 5 tables, Accepted in MNRAS

We outline two concepts to explain Ultra High Energy Cosmic Rays (UHECRs), one based on radio galaxies and their relativistic jets and terminal hot spots, and one based on relativistic Super-Novae (SNe) or Gamma Ray Bursts (GRBs) in starburst galaxies, one matching the arrival direction data in the South (the radio galaxy Cen A) and one in the North (the starburst galaxy M82). Ubiquitous neutrino emission follows accompanied by compact TeV photon emission, detectable more easily if the direction is towards Earth. The ejection of UHECRs is last. We have observed particles up to ZeV, neutrinos up to PeV, photons up to TeV, 30 - 300 Hz GW events, and hope to detect soon of order Hz to mHz GW events. Energy turnover in single low frequency GW events may be of order 10^63 erg. How can we further test these concepts? First of all by associating individual UHECR events, or directional groups of events, with chemical composition in both the Telescope Array (TA) Coll. and the Auger Coll. data. Second by identifying more TeV to PeV neutrinos with recent SMBH mergers. Third by detecting the order < mHz GW events of SMBH binaries, and identifying the galaxies host to the stellar BH mergers and their GW events in the range up to 300 Hz. Fourth by finally detecting the formation of the first generation of SMBHs and their mergers, surely a spectacular discovery., Comment: 19 pages, 1 figure, invited talk in Vulcano Workshop 2016 "Frontier Objects in Astrophysics and Particle Physics" 22nd - 28th, May 2016. The proceedings will be published in electronic from in the Frascati Physics Series (http://www.lnf.infn.it/sis/frascatiseries)

We constrain the origin of Fermi Bubbles using 2D hydrodynamical simulations of both star formation driven and black hole accretion driven wind models. We compare our results with recent observations of OVIII to OVII line ratio within and near Fermi Bubbles. Our results suggest that independent of the driving mechanisms, a low luminosity ($\mathcal{L} \sim 0.7-1\times 10^{41}$ erg s$^{-1}$) energy injection best reproduces the observed line ratio for which the shock temperature is $\approx 3\times 10^6$ K. Assuming the Galactic halo temperature to be $2\times 10^6$K, we estimate the shock velocity to be $\sim 300$ km s$^{-1}$ for a weak shock. The corresponding estimated age of the Fermi bubbles is $\sim 15-25$ Myr. Such an event can be produced either by a star formation rate of $\sim 0.5$ M$_\odot$ yr$^{-1}$ at the Galactic centre or a very low luminosity jet/accretion wind arising from the central black hole. Our analysis rules out any activity that generates an average mechanical luminosity $\gtrsim 10^{41}$ \ergps as a possible origin of the Fermi Bubbles., Comment: 14 pages, 9 figures, accepted version (MNRAS); includes updates on the electron-proton equilibrium time scale and its implications for high energy jets

The X-ray emission from active galactic nucleus (AGN) is a major component of extragalactic X-ray sky. In this paper, we use the X-ray luminosity function (XLF) and halo occupation distribution (HOD) formalism to construct a halo model for the X-ray emission from AGNs. Verifying that the two inputs (XLF and HOD) are in agreement with each other, we compute the auto-correlation power spectrum in the soft X-ray band (0.5-2 keV) due to the AGNs potentially resolved by eROSITA (extended ROentgen Survey with an Imaging Telescope Array) mission and explore the redshift and mass dependence of the power spectrum. Studying the relative contribution of the Poisson and the clustering terms to the total power, we find that at multipoles $l\lesssim 1000$ (i.e. large scales), the clustering term is larger than the Poisson term. We also forecast the potential of X-ray auto-correlation power spectrum and X-ray-lensing cross-correlation power spectrum using eROSITA and eROSITA-LSST (Large Synoptic Survey Telescope) surveys, respectively, to constrain the HOD parameters and their redshift evolution. In addition, we compute the power spectrum of the AGNs lying below the flux resolution limit of eROSITA, which is essential to understand in order to extract the X-ray signal from the hot diffuse gas present in galaxies and clusters., Comment: 14 pages, 13 figures, 3 tables, Accepted for publication in MNRAS

Recent observations have detected molecular outflows in a few nearby starburst nuclei. We discuss the physical processes at work in such an environment in order to outline a scenario that can explain the observed parameters of the phenomenon, such as the molecular mass, speed and size of the outflows. We show that outflows triggered by OB associations, with $N_{OB}\ge 10^5$ (corresponding to a star formation rate (SFR)$\ge 1$ M$_{\odot}$ yr$^{-1}$ in the nuclear region), in a stratified disk with mid-plane density $n_0\sim 200\hbox{--}1000$ cm$^{-3}$ and scale height $z_0\ge 200 (n_0/10^2 \, {\rm cm}^{-3})^{-3/5}$ pc, can form molecules in a cool dense and expanding shell. The associated molecular mass is $\ge 10^7$ M$_\odot$ at a distance of a few hundred pc, with a speed of several tens of km s$^{-1}$. We show that a SFR surface density of $10 \le \Sigma_{SFR} \le 50$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$ favours the production of molecular outflows, consistent with observed values., Comment: Accepted for publication in MNRAS, 17 pages, 18 figures

Energetic winds and radiation from massive star clusters push the surrounding gas and blow superbubbles in the interstellar medium (ISM). Using 1-D hydrodynamic simulations, we study the role of radiation in the dynamics of superbubbles driven by a young star cluster of mass $10^{6}$ M$_{\odot}$. We have considered a realistic time evolution of the mechanical power as well as radiation power of the star cluster, and detailed heating and cooling processes. We find that the ratio of the radiation pressure on the shell (shocked ISM) to the thermal pressure ($\sim10^{7}$ K) of the shocked wind region is almost independent of the ambient density, and it is greater than unity before $\lesssim 1$ Myr. We explore the parameter space of density and dust opacity of the ambient medium, and find that the size of the hot gas ($\sim$ 10$^{7}$ K) cavity is insensitive to the dust opacity ($\sigma_{d}\approx(0.1-1.5)\times 10^{-21}$ cm$^{2}$), but the structure of the photoionized ($\sim10^4$ K) gas depends on it. Most of the radiative losses occur at $\sim10^{4}$ K, with sub-dominant losses at $\lesssim 10^3$ K and $\sim10^{6}-10^{8}$ K. The superbubbles can retain as high as $\sim 10\%$ of its input energy, for an ambient density of $10^{3}\,m{\rm_{H}\,cm^{-3}}$. We discuss the role of ionization parameter and recombination-averaged density in understanding the dominant feedback mechanism. Finally, we compare our results with the observations of 30 Doradus., Comment: 19 pages, 20 figures, 1 table; Accepted for publication in MNRAS; Added references and appendix; Comments welcome

Non-gravitational feedback affects the nature of the intra-cluster medium (ICM). X-ray cooling of the ICM and in situ energy feedback from AGN's and SNe as well as {\it preheating} of the gas at epochs preceding the formation of clusters are proposed mechanisms for such feedback. While cooling and AGN feedbacks are dominant in cluster cores, the signatures of a preheated ICM are expected to be present even at large radii. To estimate the degree of preheating, with minimum confusion from AGN feedback/cooling, we study the excess entropy and non-gravitational energy profiles upto $r_{200}$ for a sample of 17 galaxy clusters using joint data sets of {\it Planck} SZ pressure and {\it ROSAT/PSPC} gas density profiles. The canonical value of preheating entropy floor of $\gtrsim 300$ keV cm$^2$, needed in order to match cluster scalings, is ruled out at $\approx 3\sigma$. We also show that the feedback energy of 1 keV/particle is ruled out at 5.2$\sigma$ beyond $r_{500}$. Our analysis takes both non-thermal pressure and clumping into account which can be important in outer regions. Our results based on the direct probe of the ICM in the outermost regions do not support any significant preheating., Comment: 6 pages, 4 figures, 1 table, Accepted in MNRAS Letters

We investigate the relative importance of two galactic outflow suppression mechanisms : a) Cosmological infall of the intergalactic gas onto the galaxy, and b) the existence of a hot circumgalactic medium (CGM). Considering only radial motion, the infall reduces the speed of outflowing gas and even halts the outflow, depending on the mass and redshift of the galaxy. For star forming galaxies there exists an upper mass limit beyond which outflows are suppressed by the gravitational field of the galaxy. We find that infall can reduce this upper mass limit approximately by a factor of two (independent of the redshift). Massive galaxies ($\gtrsim 10^{12} M_{\odot}$) host large reservoir of hot, diffuse CGM around the central part of the galaxy. The CGM acts as a barrier between the infalling and outflowing gas and provides an additional source of outflow suppression. We find that at low redshifts ($z\lesssim3.5$), the CGM is more effective than the infall in suppressing the outflows. Together, these two processes give a mass range in which galaxies are unable to have effective outflows. We also discuss the impact of outflow suppression on the enrichment history of the galaxy and its environment., Comment: 8 pages, 5 figures, 1 table, accepted for publication in MNRAS

We explore the formation of superbubbles through energy deposition by multiple supernovae (SNe) in a uniform medium. We use total energy conserving, 3-D hydrodynamic simulations to study how SNe correlated in space and time create superbubbles. While isolated SNe fizzle out completely by $\sim 1$ Myr due to radiative losses, for a realistic cluster size it is likely that subsequent SNe go off within the hot/dilute bubble and sustain the shock till the cluster lifetime. For realistic cluster sizes, we find that the bubble remains overpressured only if, for a given $n_{g0}$, $N_{\rm OB}$ is sufficiently large. While most of the input energy is still lost radiatively, superbubbles can retain up to $\sim 5-10\%$ of the input energy in form of kinetic+thermal energy till 10 Myr for ISM density $n_{g0} \approx 1$ cm$^{-3}$. We find that the mechanical efficiency decreases for higher densities ($\eta_{\rm mech} \propto n_{g0}^{-2/3}$). We compare the radii and velocities of simulated supershells with observations and the classical adiabatic model. Our simulations show that the superbubbles retain only $\lesssim 10\%$ of the injected energy, thereby explaining the observed smaller size and slower expansion of supershells. We also confirm that a sufficiently large ($\gtrsim 10^4$) number of SNe is required to go off in order to create a steady wind with a stable termination shock within the superbubble. We show that the mechanical efficiency increases with increasing resolution, and that explicit diffusion is required to obtain converged results., Comment: 20 pages, 2 tables and 21 figures, Accepted and Published in MNRAS, new convergence study presented in comparison to work by Gentry et al. 2016(http://adsabs.harvard.edu/abs/2016arXiv160601242G), Comments are welcome

We study the diffuse X-ray luminosity ($L_X$) of star forming galaxies using 2-D axisymmetric hydrodynamical simulations and analytical considerations of supernovae (SNe) driven galactic outflows. We find that the mass loading of the outflows, a crucial parameter for determining the X-ray luminosity, is constrained by the availability of gas in the central star forming region, and a competition between cooling and expansion. We show that the allowed range of the mass loading factor can explain the observed scaling of $L_X$ with star formation rate (SFR) as $L_X \propto$ SFR$^2$ for SFR $\gtrsim 1$ M$_\odot$yr$^{-1}$, and a flatter relation at low SFRs. We also show that the emission from the hot circumgalactic medium (CGM) in the halo of massive galaxies can explain the sub-linear behaviour of the $L_X-$SFR relation as well as a large scatter in the diffuse X-ray emission for low SFRs ($\lesssim$ few M$_\odot$yr$^{-1}$). Our results point out that galaxies with small SFRs and large diffuse X-ray luminosities are excellent candidates for detection of the elusive CGM., Comment: 5 pages, 4 figures, Accepted for publication in ApJL

We study the cross-correlation of distribution of galaxies, the Sunyaev-Zel'dovich (SZ) and X-ray power spectra of galaxies from current and upcoming surveys and show these to be excellent probes of the nature, i.e. extent, evolution and energetics, of the circumgalactic medium (CGM). The SZ-galaxy cross-power spectrum, especially at large multipoles, depends on the steepness of the pressure profile of the CGM. This property of the SZ signal can, thus, be used to constrain the pressure profile of the CGM. The X-ray cross power spectrum also has a similar shape. However, it is much more sensitive to the underlying density profile. We forecast the detectability of the cross-correlated galaxy distribution, SZ and X-ray signals by combining South Pole Telescope-Dark Energy Survey (SPT-DES) and eROSITA-DES/eROSITA-LSST (extended ROentgen Survey with an Imaging Telescope Array-Large Synoptic Survey Telescope) surveys, respectively. We find that, for the SPT-DES survey, the signal-to-noise ratio (SNR) peaks at high mass and redshift with SNR $\sim 9$ around $M_h\sim 10^{13} h^{-1} M_{\odot}$ and $z\sim 1.5\hbox{--} 2$ for flat density and temperature profiles. The SNR peaks at $\sim 6 (12 )$ for the eROSITA-DES (eROSITA-LSST) surveys. We also perform a Fisher matrix analysis to find the constraint on the gas fraction in the CGM in the presence or absence of an unknown redshift evolution of the gas fraction. Finally, we demonstrate that the cross-correlated SZ-galaxy and X-ray-galaxy power spectrum can be used as powerful probes of the CGM energetics and potentially discriminate between different feedback models recently proposed in the literature; for example, one can distinguish a `no active galactic nuclei feedback' scenario from a CGM energized by `fixed-velocity hot winds' at greater than $3\sigma$., Comment: 14 pages, 10 figures, 4 tables, accepted for publication in MNRAS

Using hydrodynamical simulations, we show for the first time that an episode of star formation in the center of the Milky Way, with a star-formation-rate (SFR) $\sim 0.5$ M$_\odot$ yr$^{-1}$ for $\sim 30$ Myr, can produce bubbles that resemble the Fermi Bubbles (FBs), when viewed from the solar position. The morphology, extent and multi-wavelength observations of FBs, especially X-rays, constrain various physical parameters such as SFR, age, and the circum-galactic medium (CGM) density. We show that the interaction of the CGM with the Galactic wind driven by a star formation in the central region can explain the observed surface brightness and morphological features of X-rays associated with the Fermi Bubbles. Furthermore, assuming that cosmic ray electrons are accelerated {\it in situ} by shocks and/or turbulence, the brightness and morphology of gamma-ray emission and the microwave haze can be explained. The kinematics of the cold and warm clumps in our model also matches with recent observations of absorption lines through the bubbles., Comment: 14 pages, 12 figures, Accepted for publication in MNRAS