On the requirement for mixing-heating to utilize bubble cosmic rays to heat the intracluster medium

# On the requirement for mixing-heating to utilize bubble cosmic rays to heat the intracluster medium

Noam Soker11affiliation: Department of Physics, Technion – Israel Institute of Technology, Haifa 32000, Israel; soker@physics.technion.ac.il 22affiliation: Guangdong Technion Israel Institute of Technology, Shantou 515069, Guangdong Province, China
###### Abstract

I conduct simple analytical estimates and conclude that mixing by vortices is a more efficient process to transfer the cosmic ray energy of jet-inflated bubbles to the intracluster medium (ICM) than streaming of cosmic rays along magnetic field lines. Jets and the bubbles they inflate transfer heat to the ambient gas in cooling flows in cluster of galaxies and in galaxies. The internal energy of the jet-inflated bubbles is dominate by very hot thermal gas and/or cosmic rays. Cosmic rays that stream along magnetic field lines that connect the bubbles with the ICM heat the ICM as their energy is dissipated there. I find that about half of the cosmic ray energy is dissipated in the bubbles themselves. I also find that the ICM volume that the cosmic ray streaming process heats is only about five times as large as the volume of the bubbles. The outcome of heating by streaming only is that the cosmic rays form a larger bubble filled with very hot thermal gas. Therefore, there is a requirement for a more efficient process to transfer the internal energy of the bubbles to the ICM. I suggest that this process is heating by mixing that operates very well for both cosmic rays and the very hot thermal gas inside the bubbles. This leaves mixing-heating to be the dominant heating process of cooling flows.

## 1. Introduction

Although the old “cooling-flow problem” in clusters and groups of galaxies, and in galaxies has been solved, there is still no consensus on the heating mechanism of the intracluster medium (ICM; for a recent review see Soker 2016). The ICM in many clusters of galaxies has short radiative cooling times. If no heating processes of the ICM exist large amounts of ICM should cool to low temperatures (by ICM I refer also to the interstellar medium in galactic cooling flows). The cooling-flow problem refers to the contradiction between the large mass cooling rates of the ICM that the no-heating assumption implies and the observations of much lower mass cooling rates. The solution is simply to relax the no-heating assumption and to consider the heating of the ICM in all cooling flows. The vast majority of heating models involve jets that the active galactic nucleus (AGN) at the center of the cooling flow launches. The jets and the hot low-density bubbles they inflate heat the ICM via a negative feedback mechanism that prevents a run-away heating or cooling (e.g., Fabian 2012; Farage et al. 2012; McNamara & Nulsen 2012; Gaspari et al. 2013; Pfrommer 2013; Barai et al. 2016; Soker 2016; Bîrzan et al. 2017; Iqbal et al. 2017).

The two parts that compose the negative feedback cycle are the feeding of the AGN by the ICM and the heating of the ICM by the AGN. There are tens of observational and theoretical studies in recent years that support models in which the ICM feed the AGN with cold clumps (e.g., limiting the list to the last 3 years, Choudhury & Sharma 2016; Hamer et al. 2016; Loubser et al. 2016; McNamara et al. 2016; Russell et al. 2016; Tremblay et al. 2016; Yang & Reynolds 2016; David et al. 2017; Donahue et al. 2017; Gaspari et al. 2017; Hogan et al. 2017; Meece et al. 2017; Prasad et al. 2017; Russell et al. 2017a; Voit et al. 2017; Babyk et al. 2018; Gaspari et al. 2018; Ji et al. 2018; Prasad et al. 2018; Pulido et al. 2018; Vantyghem et al. 2018; Voit 2018; Yang et al. 2018). These models can be grouped under the cold feedback mechanism that Pizzolato & Soker (2005) developed to replace the then popular Bondi accretion. In the cold feedback mechanism gas flows inward to form cold clouds and to feed the central AGN. Cold clumps are observed now down to the AGN, e.g., in Perseus (Fujita & Nagai, 2017). This implies that a cooling flow does take place, like the extreme cooling flow in the Phoenix cluster (e.g., Pinto et al. 2018), although in most cases the cooling flow is a moderate one. For that, in this study I use the term cooling flow (rather than other names that were invented later and caused some confusions; Soker 2010).

On the other hand, there is no agreement on the processes that contribute the most to the heating of the ICM by the jets that the AGN launches and the bubbles that they inflate. I find it useful to distinguish between heating processes where the jets and jet-inflated bubbles do work on the ICM and between heating by energy transport. In the first class of processes the jets and the jet-inflated bubbles do work on the ICM by exciting sound waves (e.g., Fabian et al. 2006, 2017; Tang & Churazov 2018), by driving shocks (e.g., Forman et al. 2007; Randall et al. 2015; Guo et al. 2018), by powering turbulence (e.g., De Young 2010; Gaspari et al. 2014; Zhuravleva et al. 2014, 2017), and/or by uplifting gas from inner regions (e.g., Gendron-Marsolais et al. 2017). However, there are studies that show that these processes do not heat the ICM efficiently enough, despite the fact that these processes themselves take place to some degree. Although turbulence is observed in the ICM (e.g., Zhuravleva et al. 2014, 2015; Anderson & Sunyaev 2016; Arévalo et al. 2016; Hitomi Collaboration et al. 2016; Hofmann et al. 2016; Hitomi Collaboration et al. 2017), there are studies that question turbulent heating (e.g., Falceta-Gonçalves et al. 2010; Reynolds et al. 2015; Hitomi Collaboration et al. 2016; Hillel & Soker 2017a; Bambic et al. 2018). Sound waves also are unable to supply the entire heating (e.g. Fujita & Suzuki 2005).Heating by shocks has some problems as well (e.g., Soker et al. 2016). Observations show that bubbles can uplift cooler gas (e.g,. Russell et al. 2017b; Su et al. 2017; Gendron-Marsolais et al. 2017), and some studies simulate the uplifting process (e.g., Guo et al. 2018; Churazov et al. 2001), but Hillel & Soker (2018) claim that this cannot be the main heating process of the ICM.

In the second class of heating processes the energy from the hot jet-inflated bubbles is carried into the ICM. One possibility, that is the main subject of the present study, is that cosmic rays that are accelerated within the jet-inflated bubbles, stream into the ICM and heat it (e.g. Fujita et al. 2013; Fujita & Ohira 2013; Pfrommer 2013). The second possibility is heating by mixing that works as follows. As the jets propagate through the ICM and inflate bubbles they form many vortices inside the bubble and on the boundary between the ICM and the bubbles. These vortices mix hot bubble gas into the ICM (e.g., Brüggen & Kaiser 2002; Brüggen et al. 2009; Gilkis & Soker 2012; Hillel & Soker 2014; Yang & Reynolds 2016).

Finally, some studies suggest that two or more processes out of the different process can work together, e.g., thermal conduction and cosmic rays (e.g., Guo & Oh 2008), mixing of cosmic rays that were accelerated inside jet-inflated bubbles with the the ICM, (Pfrommer, 2013), and heating by turbulence together with turbulent-mixing (e.g. Banerjee & Sharma 2014).

In a series of papers (e.g., Hillel & Soker 2017a, 2018 for most recent papers) we argued that although the inflation of bubbles also excite sound waves, shocks, and ICM turbulence, heating by mixing is much more efficient than the heating mechanisms where the jets and bubbles do work on the ICM. The heating by mixing can work for bubbles that are filled with very hot thermal gas or with cosmic rays, or a combination of the two. We did not compare the heating by mixing with the heating process where cosmic rays from the bubbles stream, instead of being mixed by vortices, into the ICM and heat it.

In light of some recent interest in cosmic ray heating by streaming (e.g., Ruszkowski et al. 2017; Ehlert et al. 2018; Jiang & Oh 2018; Ruszkowski et al. 2018), and the possibility that in many cases cosmic ray energy is the main energy content of bubbles (e,g., Abdulla et al. (2018) for MS 0735.6+7421), in this study I examine some of the properties of this process. In section 2 I consider only cosmic ray streaming, without mixing by vortices, and study the partition of cosmic ray energy between heating the bubbles themselves and heating the ICM, and in section 3 I estimate the volume of the ICM that can be heated by streaming cosmic rays. I discuss and summarize my results in section 4 where I consider mixing of cosmic rays by vortices to be much more efficient than streaming.

## 2. Heating the Bubble and the ICM

The basic assumptions of the cosmic ray heating are as follows. The heating per units volumes is (e.g., Ehlert et al. 2018)

 Hcr=|vA⋅→∇Pcr|, (1)

where is the Alfven speed, is the density, and is the cosmic ray pressure. I consider only cosmic ray streaming heating as it dominates over Coulomb and hadronic heating (e.g., Ruszkowski et al. 2017).

Consider an ideal situation where the bubble is filled with cosmic rays that are disconnected from the ICM by magnetic fields. Then in a short time magnetic field lines from the bubble reconnect to the magnetic field lines in the ICM, and the cosmic rays stream to the ICM and heat it. Namely, there is a magnetic flux tube of cross section that connects the bubble to the ICM.

Let as take the geometry as Ehlert et al. (2018) obtain in their simulations, where the cosmic rays leaks from the bottom of the bubbles in a pseudo-cylindrical stream. Namely, the cross section of the cosmic ray stream stays about constant. We can treat then the flow as a one-dimensional flow. The situation before reconnection is depicted schematically in the upper panel of Fig. 1. I take the coordinate to be along the magnetic field lines, with at the boundary between the bubble and the ICM. The initial cosmic ray pressure inside the bubble is while in the ICM it is zero. The magnetic field in the bubble, , is expected to be stronger than that in the ICM, , and the density inside the bubble much smaller than that in the ISM. As such, the Alfven speed inside the bubble is much larger than that in the ICM .

The lower panel of Fig. 1 schematically presents the flow some time after reconnection. The value of the cosmic ray pressure in the boundary between the bubble and the ICM decreases to a value of , and it extends to distance into the ICM, i.e., , while the reduction in cosmic ray pressure extends to a distance of into the bubble.

I assume that the Alfven speed in the ICM is constant, and so is the Alfven speed in the bubble, but the later is larger. The cosmic ray heating rate of the ICM per unit area, where the units area is perpendicular to the magnetic field lines, is

 ˙EI=∫DI0Hcrdx≃vA,I∫DI0∣∣∣dPcrdx∣∣∣dx≃vA,IPcr,h. (2)

In the same manner, the heating rate per unit area inside the bubble is

 ˙Eb=∫0DbHcrdx≃vA,b(Pcr,0−Pcr,h). (3)

I now further assume that in a short time following the reconnection of the bubble and the ICM field lines, the value of the cosmic ray pressure at drops from to a more or less constant value of . The value of will continue to decrease after the information of the dissipation of cosmic ray has reached the other side of the bubble, , where is the radius of the bubble. For a time of I take the value of not to change much. The energy of the cosmic ray that has been dissipated is supplied by cosmic ray energy flux from the bubble , where is the cosmic ray energy density inside the bubble. I calculate the cosmic ray energy that the bubble supplies over time by calculating the decrease of the cosmic ray energy and pressure inside the bubble per unit area. For an approximate linear relation between the cosmic ray pressure inside the bubble and distance into the bubble, as I mark in Fig. 1, over a time this flux has supplied a cosmic ray energy per unit area of

 Ecr,sup≃(4Pcr,0−4Pcr,h)Db2≃2(Pcr,0−Pcr,h)vA,bt, (4)

where I take . The dissipated cosmic ray energy per unit area at the same time is given by

 Ecr,dis=(˙EI+˙Eb)t. (5)

Equating the dissipated energy to the supplied energy and using equations (2) and (3), gives the relation

 (Pcr,0−Pcr,h)vA,b≃Pcr,hvA,I. (6)

Now we can use back equations (2) and (3), and conclude from equation (6) that the cosmic ray dissipation rate inside the bubble is about equal to that in the ICM

 ˙EI≃˙Eb. (7)

I examine some recent works on cosmic ray heating in light of equation (7). Equation (7) implies that one cannot ignore the cosmic ray energy that is transferred to heat the jet-inflated bubbles themselves. Ehlert et al. (2018), for example, turn off the cosmic ray dissipation in regions where the matter from the bubble dominate. I find this unjustified. Jacob & Pfrommer (2017a) and Jacob & Pfrommer (2017b) derive the cosmic ray heating from a steady-state solution to clusters. They do not include the inflation of bubbles self-consistently in their treatment, and therefore I claim that they substantially underestimate the dissipation of cosmic ray inside the bubbles and near the bubble surface. Namely, they substantially underestimate the dissipation of cosmic ray energy near the source of the cosmic rays.

Ruszkowski et al. (2017) present simulations of feedback heating by cosmic rays. Unfortunately they do not present the evolution with time of the velocity flow of the ICM, nor the magnetic field in the regions where they inject the jets. Therefore, I cannot use their results to estimate the role of mixing or cosmic ray heating inside the hot bubbles. But I do note that they assume that a fraction of of the energy of the jets ends as cosmic rays. Equations (7) implies that in that case, at most a fraction of of the jets’ energy will end up as cosmic ray heating of the ICM. But I argue in this paper that most of the cosmic ray energy is carried to the ICM by mixing, and not by streaming.

I neglected above the process by which turbulence inside the bubble can smooth the cosmic ray pressure, and by that reducing the dissipation of cosmic ray energy inside the bubble. However, this does not seem to be a problem to the derivation of equation (7). Firstly, I do think that there is a strong turbulence, including large vortices, inside the bubble, as the vortices play a significant role in the heating by mixing process, which I take to be the most significant heating process. Secondly, turbulence entangles magnetic field lines that reconnect. This process prevents a smooth flow of cosmic rays into the ICM, increasing further the problem that the cosmic rays do not heat a large enough volume of the ICM. I turn to discuss this problem next.

## 3. The heated ICM volume

The front of the cosmic rays advances at a speed of into the ICM. The speed might be higher than this by at most few tens of percents (e.g., Ruszkowski et al. 2017). At the same time the information that cosmic ray energy is lost from the bubble travels at a speed of into the bubble. When this information covers the entire volume of the bubble , at time , the volume of the ICM that was heated is

 VI(tb)≈VBvA,IvA,b. (8)

I calculate the energy of the cosmic ray in the bubble that supplies the dissipated cosmic ray energy from equation (4). I do this by substituting and taking to be the volume of the bubble, where is the cross section of the magnetic flux tube that connects the bubble to the ICM. This energy that the cosmic rays supply is a fraction of the entire cosmic ray energy that the bubble can supply, where is given by

 η=SEcr,sup(tb)4Pcr,0Vb≃Pcr,0−Pcr,h2Pcr,0. (9)

I crudely take the time to extract most of the available cosmic energy from the bubble to be . The ICM volume that is covered with cosmic ray heating by streaming, until the entire bubble cosmic ray energy is consumed, is

 VI,heat≈tb/ηtbVI(tb)≈VBvA,IvA,b2Pcr,0Pcr,0−Pcr,h≈4VB, (10)

where in the second equality I substituted from equations (8) and (9), and in the third equality I substitute equation (6) with the approximation that the average value of is .

The approximate equation (10) together with equation (7) imply that by streaming alone the cosmic ray energy is channeled to thermal energy that heats a region of approximately only five times the initial volume of the bubbles. From then on, we have larger hot bubbles that still need to heat most of the ICM. There must be another much more efficient process to transfer the energy from the hot bubbles to the ICM.

I further clarify one issue related to the simple model that I have been using here. I assumed above that the cosmic rays stream along magnetic filed lines. Cosmic rays can also diffuse between field lines. The diffusion between field lines is a slow process and might make the cosmic ray distribution less steep, resulting in a slower energy dissipation in the ICM. However, for the cosmic ray streaming to be significant in the heating process it must transfer the cosmic ray energy before the bubbles move outward. Let me consider an example.

In the simulation of Ehlert et al. (2018) at the lower part of the bubble is at and the upper boundary of the bubble is at . This implies that the bubble leaves the center of the cluster, , at a speed of . The cosmic rays stream at about the Alfven speed. Since the magnetic pressure in the ICM is typically less than of the thermal pressure, the Alfven speed is at most 30 per cent the sound speed. Namely, . If the cosmic rays diffusion time is slower than what I assume in the present study, then the cosmic rays will not have time to diffuse out from the bubbles and heat the center of the cluster before the bubbles move out. It is in the center of the cluster where the heating is most needed.

## 4. Discussion and Summary

By means of simple analytical estimates I examined the efficiency by which streaming of cosmic rays along magnetic field lines can transfer energy from the jet-inflated bubbles to the ICM. In section 2 I found that about half of the cosmic ray energy is dissipated in the bubbles themselves. Even before the dissipation of cosmic ray energy starts, some of the kinetic energy of the jets has been converted to thermal energy of very hot gas inside the bubbles. Over all, more than half of the kinetic energy of the jets ends up as very hot thermal gas inside the bubbles. In section 3 I found that by streaming alone the cosmic rays heat a total ICM volume of about five times the volume of the bubbles, and hence this mechanism does not cover a large enough volume of the ICM for an efficient heating to work.

Taking these two conclusions together, I argue that the effect of cosmic ray heating by streaming alone turns a bubble filled with cosmic rays to a larger bubble filled with very hot thermal gas. Therefore, although heating by cosmic ray streaming does take place, I think there is a need for an additional and a more efficient process to carry the energy from the bubbles to the ICM.

I take this more efficient process to be heating by mixing (for references see section 1). The propagation of jets through the ICM and the inflation of the bubbles form many vortices in the ICM and in the bubbles, and these mix the content of the bubbles, i.e., the very hot thermal gas, the cosmic rays, and the magnetic fields, with the ICM. When the bubble content and the ICM are well mixed we expect that magnetic fields from the bubble content and from the ICM reconnect on small scales. Then streaming of cosmic rays on small scale can make the final cosmic ray energy transfer, as much as heat conduction can locally transfer energy from the very hot thermal gas, which is now well mixed with the ICM, to the ICM.

The heating by mixing process that I refer to results from the vortices that the inflation of bubbles form (e.g., Gilkis & Soker 2012; Hillel & Soker 2016), and it is indiscriminate to the content of the bubble, and hence works for cosmic rays as well as for the thermal gas. Indeed, Pfrommer (2013) already mentioned the mixing of cosmic rays with the ICM and referred to it as an important process. Here I argue that it is more significant than heating by cosmic ray streaming, and hence I strengthen our earlier claim that heating by mixing is the main process by which jets-inflated bubbles heat the ICM (see section 1).

The simulation of the feedback process in cooling flows is a very complicated task, e.g., as Martizzi et al. (2018) show in a very recent study. There are many physical and numerical issues to consider. One of the key ingredients that simulations must include is the launching of kinetic jets. Sternberg & Soker (2008) show that injecting energy off-center instead of launching jets might miss key processes, such as the formation of many vortices. It might be that by inserting their jets in a sphere off-center Weinberger et al. (2017) and Ehlert et al. (2018) underestimated the role of heating by mixing (Hillel & Soker, 2017b).

The heating by mixing also has some implications to the mass feeding part of the feedback cycle. Pizzolato & Soker (2005) showed that non-linear perturbations in the ICM are required to form the cooling clumps that feed the AGN, as was also confirmed in following studies (e.g., Gaspari et al. 2018). Pizzolato & Soker (2005) also suggested that the AGN activity form these non-linear perturbations. This is most likely a result of the vortices and the mixing process. Therefore, in an indirect way, the results of the present study that further strengthen the mixing-heating mechanism, add also to the cold feedback mechanism.

I thank Karen Yang and Yutaka Fujita for very enlightening comments. This research was supported by the Israel Science Foundation, and by the E. and J. Bishop Research Fund at the Technion.

## References

• Abdulla et al. (2018) Abdulla, Z., Carlstrom, J. E., Mantz, A. B., et al. 2018, arXiv:1806.05050
• Anderson & Sunyaev (2016) Anderson, M., E., & Sunyaev, R. 2016, MNRAS, 459, 2806
• Arévalo et al. (2016) Arévalo, P., Churazov, E., Zhuravleva, I., Forman, W. R., & Jones, C. 2016, ApJ, 818, 14
• Babyk et al. (2018) Babyk, I. V., McNamara, B. R., Nulsen, P. E. J., Russell, H. R., Vantyghem, A. N., Hogan, M. T., & Pulido, F. A. 2018, ApJ, 862, 39
• Bambic et al. (2018) Bambic, C. J., Morsony, B., & Reynolds, C. S. 2018, arXiv:1801.06233
• Banerjee & Sharma (2014) Banerjee, N., & Sharma, P. 2014, MNRAS, 443, 687
• Barai et al. (2016) Barai, P., Murante, G., Borgani, S., Gaspari, Massimo., Granato, G. L., Monaco, P., & Ragone-Figueroa. Cinthia 2016, MNRAS, 461, 1548
• Bîrzan et al. (2017) Bîrzan, L., Rafferty, D. A., Brüggen, M., & Intema, H. T. 2017, MNRAS, 471, 1766
• Brüggen & Kaiser (2002) Brüggen, M., & Kaiser, C. R. 2002, Nature, 418, 301
• Brüggen et al. (2009) Brüggen, M., Scannapieco, E., & Heinz, S. 2009, MNRAS, 395, 2210
• Choudhury & Sharma (2016) Choudhury, P. P., & Sharma, P. 2016, MNRAS, 457, 2554
• Churazov et al. (2001) Churazov, E., Brüggen, M., Kaiser, C. R., Böhringer, H., & Forman, W. 2001, ApJ, 554, 261
• David et al. (2017) David, L. P., Vrtilek, J., O’Sullivan, E., Jones, C., Forman, W., & Sun, M. 2017, ApJ, 842, 84
• De Young (2010) De Young, D. S. 2010, ApJ, 710, 743
• Donahue et al. (2017) Donahue, M., Connor, T., Voit, G. M., & Postman, M. 2017, ApJ, 835, 216
• Ehlert et al. (2018) Ehlert, K., Weinberger, R., Pfrommer, C., Pakmor, R., & Springel, V. 2018, arXiv:1806.05679
• Fabian (2012) Fabian, A. C. 2012, ARA&A, 50, 455
• Fabian et al. (2006) Fabian, A. C., Sanders, J. S., Taylor, G. B., Allen, S. W., Crawford, C. S., Johnstone, R. M., & Iwasawa, K. 2006, MNRAS, 366, 417
• Fabian et al. (2017) Fabian, A. C., Walker, S. A., Russell, H. R., Pinto, C., Sanders, J. S., & Reynolds, C. S. 2017, MNRAS, 464, L1
• Falceta-Gonçalves et al. (2010) Falceta-Gonçalves, D., de Gouveia Dal Pino, E. M., Gallagher, J. S., & Lazarian, A. 2010, ApJ, 708, L57
• Farage et al. (2012) Farage, C. L., McGregor, P. J., & Dopita, M. A. 2012, ApJ, 747, 28
• Forman et al. (2007) Forman, W., Jones, C., Churazov, E., et al. 2007, ApJ, 665, 1057
• Fujita et al. (2013) Fujita, Y., Kimura, S., & Ohira, Y. 2013, MNRAS, 432, 1434
• Fujita & Nagai (2017) Fujita, Y., & Nagai, H. 2017, MNRAS, 465, L94
• Fujita & Ohira (2013) Fujita, Y., & Ohira, Y. 2013, MNRAS, 428, 599
• Fujita & Suzuki (2005) Fujita, Y., & Suzuki, T. K. 2005, ApJ, 630, L1
• Gaspari et al. (2013) Gaspari, M., Brighenti, F., & Ruszkowski, M. 2013, Astronomische Nachrichten, 334, 394
• Gaspari et al. (2014) Gaspari, M., Churazov, E., Nagai, D., Lau, E. T., & Zhuravleva, I. 2014, A&A, 569, A67
• Gaspari et al. (2018) Gaspari, M., McDonald, M., Hamer, S. L., et al. 2018, ApJ, 854, 167
• Gaspari et al. (2017) Gaspari, M., Temi, P., & Brighenti, F. 2017, MNRAS, 466, 677
• Gendron-Marsolais et al. (2017) Gendron-Marsolais, M., Kraft, R. P., Bogdan, A., et al. 2017, ApJ, 848, 26
• Gilkis & Soker (2012) Gilkis, A., & Soker, N. 2012, MNRAS, 427, 1482
• Guo et al. (2018) Guo, F., Duan, X., & Yuan, Y.-F. 2018, MNRAS, 473, 1332
• Guo & Oh (2008) Guo, F., & Oh, S. P. 2008, MNRAS, 384, 251
• Hamer et al. (2016) Hamer, S. L., Edge, A. C., Swinbank, A. M., et al. 2016, MNRAS, 460, 1758
• Hillel & Soker (2014) Hillel, S., & Soker, N. 2014, MNRAS, 445, 4161
• Hillel & Soker (2016) Hillel, S., & Soker, N. 2016, MNRAS, 455, 2139
• Hillel & Soker (2017a) Hillel, S., & Soker, N. 2017a, MNRAS, 466, L39
• Hillel & Soker (2017b) Hillel, S., & Soker, N. 2017b, ApJ, 845, 91
• Hillel & Soker (2018) Hillel, S., & Soker, N. 2018, Research in Astronomy and Astrophysics, 18, 081
• Hitomi Collaboration et al. (2016) Hitomi Collaboration 2016, Nature, 535, 117
• Hitomi Collaboration et al. (2017) Hitomi Collaboration, Aharonian, F., Akamatsu, H., et al. 2017, arXiv:1711.00240
• Hofmann et al. (2016) Hofmann, F., Sanders, J. S., Nandra, K., Clerc, N., & Gaspari, M. 2016, A&A, 585, A130
• Hogan et al. (2017) Hogan, M. T., McNamara, B. R., Pulido, F., et al. 2017, ApJ, 851, 66
• Iqbal et al. (2017) Iqbal, A., Kale, R., Majumdar, S., Nath, B. B., Pandge, M., Sharma, P., Malik, M. A., & Raychaudhury, S. 2017, Journal of Astrophysics and Astronomy, 38, #68
• Jacob & Pfrommer (2017a) Jacob, S., & Pfrommer, C. 2017a, MNRAS, 467, 1449
• Jacob & Pfrommer (2017b) Jacob, S., & Pfrommer, C. 2017b, MNRAS, 467, 1478
• Ji et al. (2018) Ji, S., Oh, S. P., & McCourt, M. 2018, MNRAS, 476, 852
• Jiang & Oh (2018) Jiang, Y.-F., & Oh, S. P. 2018, ApJ, 854, 5
• Loubser et al. (2016) Loubser, S. I., Babul, A., Hoekstra, H., Mahdavi, A., Donahue, M., Bildfell, C., & Voit, G. M. 2016, MNRAS, 456, 1565
• Martizzi et al. (2018) Martizzi, D., Quataert, E., Faucher-Giguere, C.-A., & Fielding, D. 2018, arXiv:1805.06461
• McNamara & Nulsen (2012) McNamara, B. R., & Nulsen, P. E. J. 2012, New Journal of Physics, 14, 055023
• McNamara et al. (2016) McNamara, B. R., Russell, H. R., Nulsen, P. E. J., Hogan, M. T., Fabian, A. C., Pulido, F., & Edge, A. C. 2016, ApJ, 830, 79
• Meece et al. (2017) Meece, G. R., Voit, G. M., & O’Shea, B. W. 2017, ApJ, 841, 133
• Pfrommer (2013) Pfrommer, C. 2013, ApJ, 779, 10
• Pinto et al. (2018) Pinto, C., Bambic, C. J., Sanders, J. S., Fabian, A. C., McDonald, M., Russell, H. R.; Liu, H., &Reynolds, C. S. 2018, arXiv:1808.02872
• Pizzolato & Soker (2005) Pizzolato, F., & Soker, N. 2005, ApJ, 632, 821
• Pizzolato & Soker (2010) Pizzolato, F., & Soker, N. 2010, MNRAS, 408, 961
• Prasad et al. (2017) Prasad, D., Sharma, P., & Babul, A. 2017, MNRAS, 471, 1531
• Prasad et al. (2018) Prasad, D., Sharma, P., & Babul, A. 2018, arXiv:1801.04282
• Pulido et al. (2018) Pulido, F. A., McNamara, B. R., Edge, A. C., et al. 2018, ApJ, 853, 177
• Randall et al. (2015) Randall, S. W., Nulsen, P. E. J., Jones, C., et al. 2015, ApJ, 805, 112
• Reynolds et al. (2015) Reynolds, C. S., Balbus, S. A., & Schekochihin, A. A. 2015, ApJ, 815, 41
• Russell et al. (2017a) Russell, H. R., McDonald, M., McNamara, B. R., et al. 2017a, ApJ, 836, 130
• Russell et al. (2016) Russell, H. R., McNamara, B. R., Fabian, A. C., et al. 2016, MNRAS, 458, 3134
• Russell et al. (2017b) Russell, H. R., McNamara, B. R., Fabian, A. C., et al. 2017b, MNRAS, 472, 4024
• Ruszkowski et al. (2017) Ruszkowski, M., Yang, H.-Y. K., & Reynolds, C. S. 2017, ApJ, 844, 13
• Ruszkowski et al. (2018) Ruszkowski, M., Yang, H.-Y. K., & Reynolds, C. S. 2018, ApJ, 858, 64
• Soker (2010) Soker, N. 2010, arXiv:1007.2249
• Soker (2016) Soker, N. 2016, New Astronomy Reviews, 75, 1
• Soker et al. (2016) Soker, N., Hillel, S., & Sternberg, A. 2016, Research in Astronomy and Astrophysics, 16, 015
• Sternberg & Soker (2008) Sternberg, A., & Soker, N. 2008, MNRAS, 389, L13
• Su et al. (2017) Su, Y., Nulsen, P. E. J., Kraft, R. P., Forman, W. R., Jones, C., Irwin, J. A., Randall, S. W., & Churazov, E. 2017, ApJ, 847, 94
• Sutherland & Dopita (1993) Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253
• Tang & Churazov (2018) Tang, X., & Churazov, E. 2018, arXiv:1801.03034
• Tremblay et al. (2016) Tremblay, G. R., Oonk, J. B. R., Combes, F., et al. 2016, Nature, 534, 218
• Vantyghem et al. (2018) Vantyghem, A. N., McNamara, B. R., Russell, H. R., et al. 2018, arXiv:1807.04331
• Voit (2018) Voit, G. M. 2018, arXiv:1803.06036
• Voit et al. (2017) Voit, G. M., Meece, G., Li, Y., O’Shea, B. W., Bryan, G. L., & Donahue, M. 2017, ApJ, 845, 80
• Weinberger et al. (2017) Weinberger, R., Ehlert, K., Pfrommer, C., Pakmor, R., & Springel, V. 2017, MNRAS, 470, 45303
• Yang & Reynolds (2016) Yang, H.-Y. K., & Reynolds, C. S. 2016, ApJ, 829, 90
• Yang et al. (2018) Yang, L., Tozzi, P., Yu, H., Lusso, E., Gaspari, M., Gilli, R., Nardini, E., & Risaliti, G. 2018, ApJ, 859, 65
• Zhuravleva et al. (2017) Zhuravleva, I., Allen, S. W., Mantz, A. B., & Werner, N. 2017, arXiv:1707.02304
• Zhuravleva et al. (2015) Zhuravleva, I., Churazov, E., Arevalo, P., et al. 2015, MNRAS, 450, 4184
• Zhuravleva et al. (2014) Zhuravleva, I., Churazov, E., Schekochihin, A. A., et al. 2014, Nature, 515, 85
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters