Search

Sgr A lobes are a pair of 15-pc-sized bipolar bubbles with co-aligned major axes perpendicular to the Galactic plane found in X-ray and radio observations of the Galactic center (GC). Their elusive origin is a vital ingredient in understanding the ongoing high energy processes at the GC. Here we perform a suite of hydrodynamic simulations to explore the tidal disruption event (TDE) outflow model for the origin of the Sgr A lobes. By following the outflow evolution in the circumnuclear medium, we find that TDE outflows naturally produce bipolar lobes delimited by forward shocks, and the dense postshock shell contributes significantly to the lobe's X-ray emission. Our fiducial run reproduces the morphology, density, temperature, and X-ray surface brightness distribution of the observed Sgr A lobes reasonably well. The current lobe age is ~3300 yr. Our model further predicts that the uplifted wake flow breaks through the ejecta-induced shock front, producing a shock-enclosed head region which is relatively dim in X-rays compared to other lobe regions. Both the dim head region and the predicted limb-brightening feature of the lobes are slightly inconsistent with current X-ray observations, and may be used to further test our model with more sensitive future X-ray observations. While light narrow jets and massive wide winds from TDE events usually do not reproduce the observed oval-shaped morphology of the lobes, the TDE outflow in our fiducial run is massive and yet narrow. Whether it is a jet or wind remains uncertain and future simulations covering both the outflow acceleration region and its pc-scale evolution would be very helpful in determining whether the Sgr A lobes indeed originate from a pair of TDE jets or winds., Comment: 10 pages, 10 figures. Comments welcome

We report the discovery of the 10 kilo-parsec (kpc) scale radio lobes in the Sombrero galaxy (NGC 4594), using data from the Continuum Halos in Nearby Galaxies - an Expanded Very Large Array (VLA) Survey (CHANG-ES) project. We further examine the balance between the magnetic pressure inside the lobes and the thermal pressure of the ambient hot gas. At the radii $r$ of ~(1-10) kpc, the magnetic pressure inside the lobes and the thermal pressure of the ambient hot gas are generally in balance. This implies that the jets could expand into the surroundings at least to r ~ 10 kpc. The feedback from the active galactic nucleus (AGN) jet responsible for the large-scale lobes may help to explain the unusually high X-ray luminosity of this massive quiescent isolated disk galaxy, although more theoretical work is needed to further examine this possibility., Comment: 8 pages, 3 figures, 2 tables, accepted for publication in ApJ

Multi-wavelength observations indicate that the intracluster medium in some galaxy clusters contains cold filaments, while their formation mechanism remains debated. Using hydrodynamic simulations, we show that cold filaments could naturally condense out of hot gaseous wake flows uplifted by the jet-inflated active galactic nucleus (AGN) bubbles. Consistent with observations, the simulated filaments extend to tens of kiloparsecs from the cluster center, with a representative mass of $\rm 10^{8}- 10^{9}\ M_{\odot}$ for a typical AGN outburst energy of $10^{60}~ \rm erg$. They show smooth velocity gradients, stretching typically from inner inflows to outer outflows with velocity dispersions of several hundred kilometers per second. The properties of cold filaments are affected substantially by jet properties. Compared to kinetic-energy-dominated jets, thermal-energy-dominated jets are easier to produce long cold filaments with large masses as observed. AGN jets with an early turn-on time, a low jet base, or a very high power tend to overheat the cluster center, and produce short cold filaments that take a relatively long time to condense out., Comment: 15 pages, 8 figures, revised version, published in ApJ

The Hot Universe Baryon Surveyor (HUBS) is a proposed space-based X-ray telescope for detecting X-ray emissions from the hot gas content in our universe. With its unprecedented spatially-resolved high-resolution spectroscopy and large field of view, the HUBS mission will be uniquely qualified to measure the physical and chemical properties of the hot gas in the interstellar medium, the circumgalactic medium, the intergalactic medium, and the intracluster medium. These measurements will be valuable for two key scientific goals of HUBS, namely to unravel the AGN and stellar feedback physics that governs the formation and evolution of galaxies, and to probe the baryon budget and multi-phase states from galactic to cosmological scales. In addition to these two goals, the HUBS mission will also help us solve some problems in the fields of galaxy clusters, AGNs, diffuse X-ray backgrounds, supernova remnants, and compact objects. This paper discusses the perspective of advancing these fields using the HUBS telescope., Comment: 52 pages, 22 figures. Accepted for publication in Science China: Physics, Mechanics and Astronomy

Inflows at the dusty torus and smaller scales is crucial to investigate the process of supermassive black hole accretion. However, only few cases of inflowing gas at small scales have been reported through redshifted broad absorption lines so far. Here we report 9 redshifted narrow absorption lines (NALs) of $\rm Mg^+$ ions with inflowing speeds of 1071 -- 1979 km/s, which are likely along the directions close to the axes of accretion disks. The quasars showing inflowing Mg II NALs have on average slightly smaller Eddington ratios when compared to the sources with outflow Mg II NALs. The upper limits of locations of the detected NALs are at parsec scale, around the distances of dusty tori to central SMBHs. The one possible origin of these infalling NALs is from dusty tori. Of course, these infalling NALs can also be naturally explained by chaotic cold accretion resulted from the nonlinear interaction of active galactic nucleus (AGN) jets with the interstellar medium, and these cold gaseous blobs may originally precipitate in metal-rich trailing outflows uplifted by AGN jet ejecta. The infalling NALs may thus provide direct evidence for cold gas precipitation and accretion in AGN feedback processes, and provide the direct evidence of inflowing gas along the directions close to quasar jets and at parsec scale. It does not matter whether these infalling NALs are from the dusty tori or the interaction of AGN jets with the ISM, the infalling NALs cannot provide sufficient fuels to power the quasars., Comment: 12 page, 7 figures, Accepted by A&A

The hot halo gas distribution in the inner Milky Way (MW) contains key fossil records of the past energetic feedback processes in the Galactic center. Here we adopt a variety of spherical and disk-like MW halo gas models as initial conditions in a series of simulations to investigate the formation of the Fermi bubbles in the jet-shock scenario. The simulation results are compared directly with relevant X-ray and gamma-ray observations of the Fermi bubbles to constrain the halo gas distribution in the inner Galaxy before the Fermi bubble event. Our best-fit gas density distribution can be described by a power law in radius $n_{\rm e}(r)=0.01(r/1 \text{~kpc})^{-1.5}$ cm$^{-3}$. Our study can not determine if there is an inner density core, which if exists, should be very small with size $r_{c} \lesssim 0.5$ kpc. When extrapolating to large radii $r\sim 50-90$ kpc, our derived density distribution lies appreciably below the recently estimated gas densities from ram-pressure stripping calculations, suggesting that the halo gas density profile either flattens out or has one or more discontinuities within $10 \lesssim r \lesssim 50$ kpc. Some of these discontinuities may be related to the eROSITA bubbles, and our derived gas density profile may correspond to the hot gas distribution in the inner eROSITA bubbles about $5$ Myr ago., Comment: 16 pages, 12 figures. Accepted for publication ApJ

Tidal disruption events (TDEs) occurred in active galactic nuclei (AGNs) are a special class of sources with outstanding scientific significance. TDEs can generate ultrafast winds, which should almost inevitably collide with the preexisting AGN dusty tori. We perform analytical calculations and simulations on the wind-torus interactions and find such a process can generate considerable X-ray afterglow radiation several years or decades later after the TDE outburst. This provides a new origin for the years delayed X-rays in TDEs. The X-ray luminosity can reach 10^{41-42} erg/s, and the light curve characteristics depend on the parameters of winds and tori. We apply the model to two TDE candidates, and provide lower limits on the masses of the disrupted stars, as well as rigorous constraints on the gas densities of tori. Our results suggest that the observations of the time delay, spectral shape, luminosity and the light curve of the X-ray afterglow can be used to constrain the physical parameters of both TDE winds and tori, including the wind velocity, wind density and torus density., Comment: 13 pages, submitted to ApJ

Mechanical active galactic nucleus (AGN) feedback plays a key role in massive galaxies, galaxy groups and clusters. However, the energy content of AGN jets that mediate this feedback process is still far from clear. Here we present a preliminary study of radial elongations $\tau$ of a large sample of X-ray cavities, which are apparently produced by mechanical AGN feedback. All the cavities in our sample are elongated along the angular (type-I) or jet directions (type-II), or nearly circular (type-III). The observed value of $\tau$ roughly decreases as the cavities rise buoyantly, confirming the same trend found in hydrodynamic simulations. For young cavities, both type-I and II cavities exist, and the latter dominates. Assuming a spheroidal cavity shape, we derive an analytical relation between the intrinsic radial elongation $\bar{\tau}$ and the inclination-angle-dependent value of $\tau$, showing that projection effect makes cavities appear more circular, but does not change type-I cavities into type-II ones, or vice versa. We summarize radial elongations of young cavities in simulations, finding that $\bar{\tau}$ increases with the kinetic fraction of AGN jets. While mild jets always produce type-II cavities, thermal-energy-dominated strong jets produce type-I cavities, and kinetic-energy-dominated strong jets produce type-II cavities. Our results suggest that some AGN jets are strong and dominated by thermal energy (or cosmic rays). However, these jets do not dominate in AGN feedback. If most jets are dominated by non-kinetic energies, they should be mainly mild jets. If most jets are strong, they must be mainly dominated by the kinetic energy., Comment: Slightly revised version accepted for publication in ApJ. 11 pages, 6 figures

We propose a novel method to constrain the Milky Way (MW) mass $M_{\rm vir}$ with its corona temperature observations. For a given corona density profile, one can derive its temperature distribution assuming a generalized equilibrium model with non-thermal pressure support. While the derived temperature profile decreases substantially with radius, the X-ray-emission-weighted average temperature, which depends most sensitively on $M_{\rm vir}$, is quite uniform toward different sight lines, consistent with X-ray observations. For an Navarro-Frenk-White (NFW) total matter distribution, the corona density profile should be cored, and we constrain $M_{\rm vir}=(1.19$ - $2.95) \times 10^{12} M_{\rm sun}$. For a total matter distribution contributed by an NFW dark matter profile and central baryons, the corona density profile should be cuspy and $M_{\rm vir,dm}=(1.34$ - $5.44) \times 10^{12} M_{\rm sun}$. Non-thermal pressure support leads to even higher values of $M_{\rm vir}$, while a lower MW mass may be possible if the corona is accelerating outward. This method is independent of the total corona mass, its metallicity, and temperature at very large radii., Comment: Significantly revised version, accepted for publication in ApJ Letters. Discussions largely expanded. 7 pages, 4 figures

Active galactic nucleus (AGN) jets are believed to be important in solving the cooling flow problem in the intracluster medium (ICM), while the detailed mechanism is still in debate. Here we present a systematic study on the energy coupling efficiency $\eta_{\rm cp}$, the fraction of AGN jet energy transferred to the ICM. We first estimate the values of $\eta_{\rm cp}$ analytically in two extreme cases, which are further confirmed and extended with a parameter study of spherical outbursts in a uniform medium using hydrodynamic simulations. We find that $\eta_{\rm cp}$ increases from $\sim 0.4$ for a weak isobaric injection to $\gtrsim 0.8$ for a powerful point injection. For any given outburst energy, we find two characteristic outburst powers that separate these two extreme cases. We then investigate the energy coupling efficiency of AGN jet outbursts in a realistic ICM with hydrodynamic simulations, finding that jet outbursts are intrinsically different from spherical outbursts. For both powerful and weak jet outbursts, $\eta_{\rm cp}$ is typically around $0.7-0.9$, partly due to the non-spherical nature of jet outbursts, which produce backflows emanating from the hotspots, significantly enhancing the ejecta-ICM interaction. While for powerful outbursts a dominant fraction of the energy transferred from the jet to the ICM is dissipated by shocks, shock dissipation only accounts for $\lesssim 30\%$ of the injected jet energy for weak outbursts. While both powerful and weak outbursts could efficiently heat cooling flows, powerful thermal-energy-dominated jets are most effective in delaying the onset of the central cooling catastrophe., Comment: 20 pages, 12 figures, slightly revised version, Accepted for publication in ApJ

The Fermi bubbles are two giant bubbles in gamma rays lying above and below the Galactic center (GC). Despite numerous studies on the bubbles, their origin and emission mechanism remain elusive. Here we use a suite of hydrodynamic simulations to study the scenario where the cosmic rays (CRs) in the bubbles are mainly accelerated at the forward shocks driven by a pair of opposing jets from Sgr A*. We find that an active galactic nucleus (AGN) jet event happened $5-6$ Myr ago can naturally reproduce the bilobular morphology of the bubbles, and the postshock gas temperature in the bubbles is heated to $\sim0.4$ keV, consistent with recent X-ray observations. The forward shocks compress the hot halo gas, and at low latitudes, the compressed gas shows an X-shaped structure, naturally explaining the biconical X-ray structure in the ROSAT 1.5 keV map in both morphology and X-ray surface brightness. CR acceleration is most efficient in the head regions of the bubbles during the first 2 Myrs. The opposing jets release a total energy of $\sim 10^{55}$ erg with an Eddington ratio of $\sim 10^{-3}$, which falls well in the range of the hot accretion flow mode for black holes. Our simulations further show that the forward shocks driven by spherical winds at the GC typically produce bubbles with much wider bases than observed, and could not reproduce the biconical X-ray structure at low latitudes. This suggests that starburst or AGN winds are unlikely the origin of the bubbles in the shock scenario., Comment: Slightly revised version, accepted for publication in ApJ. 15 pages, 15 figures

Theoretical and observational arguments suggest that there is a large amount of hot ($\sim 10^6$ K), diffuse gas residing in the Milky Way's halo, while its total mass and spatial distribution are still unclear. In this work, we present a general model for the gas density distribution in the Galactic halo, and investigate the gas evolution under radiative cooling with a series of 2D hydrodynamic simulations. We find that the mass inflow rate in the developed cooling flow increases with gas metallicity and the total gas mass in the halo. For a fixed halo gas mass, the spatial gas distribution affects the onset time of the cooling catastrophe, which starts earlier when the gas distribution is more centrally-peaked, but does not substantially affect the final mass inflow rate. The gravity from the Galactic bulge and disk affects gas properties in inner regions, but has little effect on the final inflow rate either. We confirm our results by investigating cooling flows in several density models adopted from the literature, including the Navarro-Frenk-White (NFW) model, the cored-NFW model, the Maller & Bullock model, and the $\beta$ model. Typical mass inflow rates in our simulations range from $\sim 5 M_{\odot}$ yr$^{-1}$ to $\sim 60 M_{\odot}$ yr$^{-1}$, and are much higher than the observed star formation rate in our Galaxy, suggesting that stellar and active galactic nucleus feedback processes may play important roles in the evolution of the Milky Way (MW) and MW-type galaxies., Comment: Slightly revised version, accepted for publication in ApJ. 15 pages, 13 figures

This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected., Comment: 7 pages, 3 figures; proceedings of the IAUS342 "Perseus in Sicily: from black hole to cluster outskirts"

Recent Chandra X-ray observations of many galaxy clusters find evidence for hot metal-rich outflows preferentially aligned with the large-scale axes of X-ray cavities with typical outflow masses of around $10^{9} $ - $10^{10} M_{\odot}$. Here we perform a suite of three hydrodynamic simulations to investigate whether AGN jets could drive these metal-rich outflows in a representative cluster. By using both the tracer variable and virtual particle methods, and additionally following the gas metallicity evolution, we show that metal-rich gas initially located in central regions can indeed be uplifted by the AGN bubble to large distances, a phenomenon called Darwin drift in fluid mechanics, and forming a filamentary trailing outflow extending beyond $100$ kpc behind the bubble. The gas entrained in the trailing outflow is entirely outflowing with an average outflow rate of nearly $100M_{\odot}$/yr during the first $100$ Myr, and at later times, a growing lower part flows back towards the cluster center due to gravity. The outflow mass rises up to about $10^{10} M_{\odot}$ with entrained iron mass of about $10^{6} - 10^{7}M_{\odot}$, consistent with observations and predictions from the drift model. By the end of our simulation ($\sim 800$ Myr after the AGN event), several $10^{9}M_{\odot}$ of the uplifted high-metallicity gas still remains at large altitudes, potentially contributing to the enrichment of the bulk ICM and the broadening of central metallicity peaks observed in cool core clusters., Comment: 9 pages, 5 figures, slightly revised version, Accepted for publication in ApJ

Galaxy clusters host a large reservoir of diffuse plasma with radially-varying temperature profiles. The efficiency of thermal conduction in the intracluster medium (ICM) is complicated by the existence of turbulence and magnetic fields, and has received a lot of attention in the literature. Previous studies suggest that the magnetothermal instability developed in outer regions of galaxy clusters would drive magnetic field lines preferentially radial, resulting in efficient conduction along the radial direction. Using a series of spherically-symmetric simulations, here we investigate the impact of thermal conduction on the observed temperature distributions in outer regions of three massive clusters, and find that thermal conduction substantially modifies the ICM temperature profile. Within 3 Gyr, the gas temperature at a representative radius of $0.3r_{500}$ typically decreases by ~10 - 20% and the average temperature slope between $0.3r_{500}$ and $r_{500}$ drops by ~ 30 - 40%, indicating that the observed ICM would not stay in a long-term equilibrium state in the presence of thermal conduction. However, X-ray observations show that the outer regions of massive clusters have remarkably similar radially-declining temperature profiles, suggesting that they should be quite stable. Our study thus suggests that the effective conductivity along the radial direction must be suppressed below the Spitzer value by a factor of 10 or more, unless additional heating sources offset conductive cooling and maintain the observed temperature distributions. Our study provides a smoking-gun evidence for the suppression of parallel conduction along magnetic field lines in low-collisionality plasmas by kinetic mirror or whistler instabilities., Comment: Slightly revised version, accepted for publication in ApJ. 11 pages, 7 figures

We investigate the effects of AGN feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. The inner boundary of the simulation is chosen so that the Bondi radius is resolved. Compared to previous works, the two accretion modes, namely hot and cold, which correspond to different accretion rates and have different radiation and wind outputs, are carefully discriminated and the feedback effects by radiation and wind in each mode are taken into account. The most updated AGN physics, including the descriptions of radiation and wind from the hot accretion flows and wind from cold accretion disks, are adopted. Physical processes like star formation, Type Ia and Type II supernovae are taken into account. We study the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. We compare our simulation results with observations and find general consistency. Comparisons with previous simulation works find significant differences, indicating the importance of AGN physics. The respective roles of radiation and wind feedbacks are examined and it is found that they are different for different problems of interest such as AGN luminosity and star formation. We find that it is hard to neglect any of them, so we suggest to use the names of "cold feedback mode" and "hot feedback mode" to replace the currently used ones., Comment: 25 pages, 14 figures, accepted for publication to ApJ (Revised to match version published in ApJ)

In this paper, we report the peculiar HI morphology of the cluster spiral galaxy NGC 6145, which has a 150 kpc HI filament on one side that is nearly parallel to its major axis. This filament is made up of several HI clouds and the diffuse HI gas between them, with no optical counterparts. We compare its HI distribution with other one-sided HI distributions in the literature, and find that the overall HI distribution is very different from the typical tidal and ram-pressure stripped HI shape, and its morphology is inconsistent with being a pure accretion event. Only about 30% of the total HI gas is anchored on the stellar disk, while most of HI gas forms the filament in the west. At a projected distance of 122 kpc, we find a massive elliptical companion (NGC 6146) with extended radio emission, whose axis points to an HI gap in NGC 6145. The velocity of the HI filament shows an overall light-of- sight motion of 80 to 180 km/s with respect to NGC 6145. Using the long-slit spectra of NGC 6145 along its major stellar axis, we find that some outer regions show enhanced star formation, while in contrast, almost no star formation activities are found in its center (less than 2 kpc). Pure accretion, tidal or ram-pressure stripping is not likely to produce the observed HI filament. An alternative explanation is the jet-stripping from NGC 6146, although direct evidence for a jet-cold gas interaction has not been found., Comment: 12 pages, 6 figures; Accepted for publication in AJ

The cooling flow problem is one of the central problems in galaxy clusters, and active galactic nucleus (AGN) feedback is considered to play a key role in offsetting cooling. However, how AGN jets heat and suppress cooling flows remains highly debated. Using an idealized simulation of a cool-core cluster, we study the development of central cooling catastrophe and how a subsequent powerful AGN jet event averts cooling flows, with a focus on complex gasdynamical processes involved. We find that the jet drives a bow shock, which reverses cooling inflows and overheats inner cool core regions. The shocked gas moves outward in a rarefaction wave, which rarefies the dense core and adiabatically transports a significant fraction of heated energy to outer regions. As the rarefaction wave propagates away, inflows resume in the cluster core, but a trailing outflow is uplifted by the AGN bubble, preventing gas accumulation and catastrophic cooling in central regions. Inflows and trailing outflows constitute meridional circulations in the cluster core. At later times, trailing outflows fall back to the cluster centre, triggering central cooling catastrophe and potentially a new generation of AGN feedback. We thus envisage a picture of cool cluster cores going through cycles of cooling-induced contraction and AGN-induced expansion. This picture naturally predicts an anti-correlation between the gas fraction (or X-ray luminosity) of cool cores and the central gas entropy, which may be tested by X-ray observations., Comment: Slightly revised version, accepted for publication in MNRAS. 14 pages, 10 figures

The nature and origin of the Fermi bubbles detected in the inner Galaxy remain elusive. In this paper, we briefly discuss some recent theoretical and observational developments, with a focus on the AGN jet model. Analogous to radio lobes observed in massive galaxies, the Fermi bubbles could be naturally produced by a pair of opposing jets emanating nearly along the Galaxy's rotation axis from the Galactic center. Our two-fluid hydrodynamic simulations reproduce quite well the bubble location and shape, and interface instabilities at the bubble surface could be effectively suppressed by shear viscosity. We briefly comment on some potential issues related to our model, which may lead to future progress., Comment: 4 pages, 1 figure. Proceedings of IAU Symposium 322 (The Multi-Messenger Astrophysics of the Galactic Centre)

Radio-mode AGN feedback plays a key role in the evolution of galaxy groups and clusters. Its physical origin lies in the kpc-scale interaction of AGN jets with the intracluster medium. Large-scale jet simulations often initiate light internally-supersonic jets with density contrast $0.01<\eta<1$. Here we argue for the first time for the importance of very-light ($\eta<0.01$) internally-subsonic jets. We investigated the shapes of young X-ray cavities produced in a suite of hydrodynamic simulations, and found that bottom-wide cavities are always produced by internally-subsonic jets, while internally-supersonic jets inflate cylindrical, center-wide, or top-wide cavities. We found examples of real cavities with shapes analogous to those inflated in our simulations by internally-subsonic and internally-supersonic jets, suggesting a dichotomy of AGN jets according to their internal Mach numbers. We further studied the long-term cavity evolution, and found that old cavities resulted from light jets spread along the jet direction, while those produced by very light jets are significantly elongated along the perpendicular direction. The northwestern ghost cavity in Perseus is pancake-shaped, providing tentative evidence for the existence of very light jets. Our simulations show that very-light internally-subsonic jets decelerate faster and rise much slower in the ICM than light internally-supersonic jets, possibly depositing a larger fraction of jet energy to cluster cores and alleviating the problem of low coupling efficiencies found previously. The internal Mach number points to the jet's energy content, and internally-subsonic jets are energetically dominated by non-kinetic energy, such as thermal energy, cosmic rays, or magnetic fields., Comment: Slightly revised version accepted for publication in ApJ. 8 pages, 5 figures