J. Astrophys. Astr. (2021) 42:12 Ó Indian Academy of Sciences

https://doi.org/10.1007/s12036-021-09732-4Sadhana(0123456789().,-volV)FT3](0123456789().,-volV)

LETTER

Is the binding energy of related to their core black hole mass?

C. SIVARAM1 and KENATH ARUN2,*

1Indian Institute of Astrophysics, Bangalore 560 034, India. 2Christ Junior College, Bangalore 560 029, India. E-mail: [email protected]

MS received 28 February 2020; accepted 4 January 2021

Abstract. Most of the large galaxies host a supermassive black hole, but their origin is still not well understood. In this paper, we look at a possible connection between the gravitational binding energies of large galaxies, etc., and the masses of their central black holes. Using this relation (between the gravitational binding energy of the host structure and the black hole energy), we argue why globular clusters are unlikely to harbour large black holes and why dwarf galaxies, if they have to host black holes, should have observed mass-to-light ratios of *100.

Keywords. Binding energy—black holes—dwarf galaxies.

It is by now well established that most large galaxies correlations (Silk & Rees 1998). Also, there is no (spiral, elliptical, etc.) host a supermassive black hole understanding as to whether these relations evolve (SMBH) in their centre (Kaper et al. 2001). Again with the evolution of the (and the central BH) active galactic nuclei, quasars, etc., are powered by or are set at the time of their formation (Di Matteo the gravitational energy of matter accreting onto these et al. 2005; Volonteri & Natarajan 2009). Again, central SMBHs (Melia 2007). These black hole many other large stellar conglomerations such as (BH) masses are typically of several millions of globular clusters or dwarf galaxies also do not seem to solar masses and can be as large as a few billion solar host massive BHs in their core regions. masses, such as in the case of M87. The origin of these Here we propose a possible connection between the SMBHs at the galactic cores is still an enigma. Again gravitational binding energies of large galaxies, etc., the BH masses have been related to the galactic bulge and the masses of their central BHs. The idea is mass, suggesting a common origin (Zel’dovich & similar to what happens in the core collapse of mas- Novikov 1971; Shapiro & Teukolsky 1983). sive . The gravitational binding energy released in Various studies have pointed out the correlation stellar collapse is carried away mainly by neutrinos (as between the masses of the SMBHs at the centre of in the case of SN1987A). The total energy carried galaxies and the different characteristics of their host away is just the gravitational binding energy of the galaxies (Kormendy & Richstone 1995), such as the remnant neutron (Burrows 1990). spheroid luminosity and spheroid mass (Magorrian In the case of SN1987A, the estimated 3 Â 1053 erg et al. 1998), the spheroid velocity dispersion (Gultekin carried away by neutrinos just corresponds to the et al. 2009), and the kinetic energy of random motions binding energy of a 1:4 M neutron star (Bahcall 1989). (Feoli & Mancini 2009). The Kormendy relation More massive stars would release more binding energy, connects the estimated mass of the BH to the K-band and this would correspond to the formation of a BH. The luminosity and the velocity dispersion of the bulge of gravitational binding energy released in the formation the host galaxy measured outside the BH sphere of of a BH would correspond to its rest energy, i.e., to the influence (Kormendy 1977). mass of the BH (for the same reason, a supernova would However, the underlying physics of these various not give rise to a white dwarf as a remnant, as the observations is not yet fully understood. It is thought binding energy of the white dwarf is two or more orders that feedback processes are responsible for such smaller than the energy released in the explosion). 12 Page 2 of 4 J. Astrophys. Astr. (2021) 42:12

By analogy, we argue that whatever processes led The binding energy of the cluster is as follows: to the formation of the BHs at the galaxy cores was GM2 connected with the formation of the galaxy and should ðÞBE ¼ cluster  1051 erg, ð7Þ cluster R thus be related to the gravitational binding energy of cluster the final structure which forms. which is lesser than the binding energy of even neu- We note that in the case of galactic structures, the tron stars. binding energy of the parent galaxy is comparable In the case of some globular clusters such as the with that of the BH hosted by the galaxy. The binding Mayall II (G1) in the , there is 4 energy of a typical galaxy such as the Milky Way is evidence of a 2 Â 10 M BH. There are also studies given by (Meylan et al. 2001) that have suggested that the 2 G1 could in fact be the core of a dwarf GMspiral ðÞBE ¼  1061 erg elliptical galaxy. The BH in G1 could have been spiral R ÀÁspiral formed at an earlier epoch when the cluster was more 12 Mspiral  10 M ; Rspiral  10 kpc : ð1Þ closely packed or due to the merger of few clusters. X-ray and radio emissions from Mayall II appear to be The galaxy harbours a  106 M BH. The energy consistent with an intermediate-mass BH. The energy (mass) of the BH is as follows: 3 2 57 of the 10 M BH is E ¼ MBHc  10 erg. In order 2 2 GMBH 61 for the cluster to have a binding energy of this order, ðÞBE BH¼ MBHc ¼  10 erg: ð2Þ RBH its size should be

The equality of Equations (1) and (2) is suggestive GM2 of what happens in the release of gravitational binding R ¼ cluster  10À4 pc, ð8Þ M c2 energy in stellar collapses leading to the formation of BH 6 a remnant compact object. where Mcluster  10 M . Again for large elliptical galaxies, mass is few As the potential decreases, the kinetic energy of the 14 orders higher, that is Melliptical  10 M , so the cluster increases, hence spreading it out over time. binding energy is as follows: One of the possible models for the formation of an ðÞBE  1063 erg: ð3Þ IMBH is the collapse of a cluster of stars (Miller & elliptical Hamilton 2002). The collapsed core can accrete Thus, galaxies, such as M87 harbour billion solar matter ejected during the formation. If a collection of mass BH, whose binding energy is as follows: a thousand 10 solar mass stars in a volume of a 63 cube collapses, ejecting 30% of its mass, then the total ðÞBE BH  10 erg: ð4Þ 33 mass of the ambient gas, nmP  6 Â 10 kg. If *30% Therefore, we have, in the case of galaxies, the of the mass of each star is ejected, then the accretion 2 2 3 7 binding energy of the galaxy being comparable with rate given by m_ ¼ 16pnmPG M =c  4 Â 10 the BH energy: kg sÀ1. For a constant density, the dynamical friction GM2 (Chandrasekhar 1942) is given as follows: M c2 ¼ galaxy : ð5Þ BH R G2M2q Hence, the fraction of BH mass to the galactic mass Fd ¼ C ; ð9Þ v2 is given by m M GM where G is the gravitational constant, M the mass of BH ¼ galaxy  10À5 À 10À6: ð6Þ M Rc2 the moving object, q the density, and vm the velocity galaxy of the object in the frame in which the surrounding 8 This implies that smaller galaxies of 10 M could matter was initially at rest. C is not a constant but harbour BHs of 100 M . depends on how vm compares with the velocity dis- This is in agreement with the observations of host persion of the surrounding matter. galaxy mass and its central BH mass, as shown in The dynamical friction comes into effect in the Table 1 (Bandara et al. 2009). evolution of the cluster due to the interaction between The result from Equation (6) could also provide a the ambient gas and dust, and the central Intermediate reason as to why globular clusters do not host a BH. Mass Black Hole (IMBH). The initial BH is subjected J. Astrophys. Astr. (2021) 42:12 Page 3 of 4 12

Table 1. Correlation between SMBH mass (MBH) and mass of its host galaxy (Mgalaxy).

Galaxy system MBHðÞM MgalaxyðÞM MBH=Mgalaxy SDSS J0029-0055 2:65  108 1:04  1013 2:5  10À5 SDSS J0252?0039 6:48  107 1:27  1013 5:1  10À6 SDSS J0330-0020 1:88  108 1:51  1013 1:2  10À5 SDSS J0728?3835 1:90  108 1:73  1013 1:1  10À5 SDSS J0737?3216 1:28  109 2:41  1013 5:3  10À5 SDSS J0822?2652 4:67  108 2:00  1013 2:3  10À5 SDSS J0955?0101 1:14  108 1:21  1013 9:4  10À6 SDSS J0959?0410 2:25  108 1:08  1013 2:1  10À5 SDSS J1142?1001 2:81  108 1:68  1013 1:7  10À5 SDSS J1403?0006 1:74  108 1:18  1013 1:5  10À5 SDSS J1451-0239 2:21  108 1:17  1013 1:9  10À5 SDSS J1538?5817 1:06  108 1:17  1013 9:1  10À6 SDSS J2238-0754 1:47  108 1:44  1013 1:0  10À5 SDSS J2303?1422 3:83  108 2:57  1013 1:5  10À5 SDSS J2341?0000 1:44  108 1:87  1013 7:7  10À6 to a gravitational force due to the , which is by about 1 kilo-parsec from the centre. The frequency given by correspondingpffiffiffiffiffiffiffiffiffiffiffiffiffiffi to the oscillation of the BH, x = À13 À1 2 k=MBH  10 s . r /ðÞ¼r 4pGqðÞ¼r 4pGmnðÞ r ; ð10Þ For the M82 system, the number density of the stars 4 À3 where nrðÞis the density distribution of stars and m in the cluster,  10 pc , with a typical mass of 1 M the typical mass of the star. has a velocity dispersion, r  30 km sÀ1. We obtain 12 For a constant density, qðÞ¼r q0, the potential is as the damping time for the system as tdamp  4  10 s, follows: which is of the order of 105 years. With the effects of  1 dynamical friction, the relaxation time for the system /ðÞ¼Àr 2pGq R2 À r2 ; ð11Þ 7 0 3 to form the IMBH is tform  10 years. For a denser core, i.e.,  103=ðÞ0:01 pc 3, the time taken to form the where R is the radius of the cluster. IMBH will be  106 years (Sivaram & Arun 2014). The gravitational force on the BH is given by Similarly, certain dwarf galaxies could also harbour 4 103 M BH, even though their binding energy (with F ¼ÀM r/ðÞ¼Àr pGq M r: ð12Þ g BH 3 0 BH visible mass) suggests that they should not. The The particle inside a homogeneous gravitational binding energy of the is as follows: system performs simple harmonic motion. 2 The equation of motion of the BH in the star cluster GMdwarf 55 ðÞBE dwarf¼  10 erg: ð14Þ is given by Rdwarf The energy of the 103 M BH is of the order of 1057 d2r dr erg. This would suggest that there is a substantial MBH 2 þ kr þ c ¼ 0; ð13Þ dt dt amount of dark matter to make up for the binding where k ¼ð4=3ÞpGq M , c ¼À3pG2mnðÞ r M2 = energy of the dwarf galaxy (Faber & Gallagher 1979; pffiffiffi 0 BH BH ð 2rÞ3 and r the velocity dispersion of stars in the Kolb & Turner 1990): cluster. The BH undergoes a damped oscillation in the M star cluster, with a damping time corresponding to  100: ð15Þ L tdamp ¼ MBH=c. In the case of the M82 system , which harbours an Recent studies (Reines et al. 2020) with Very Large IMBH of mass in the range of 500 M (Pasham et al. Array (VLA) have found 13 massive BHs in dwarf 2014), the BH is not found at the centre but displaced galaxies that are less than a billion light-years from 12 Page 4 of 4 J. Astrophys. Astr. (2021) 42:12

Earth. Their study indicates that roughly half of the Kaper L. et al. 2001, Black Hole in Binaries and Galactic massive BHs in dwarf galaxies are not in the centres Nuclei. Springer, New York of those galaxies. This probably indicates that the Kolb E. W., Turner M. S. 1990, The Early . galaxies likely have merged with others earlier in their Addison-Wesley, Redwood City history, as elaborated in this work, thus maintaining Kormendy J. 1977, Astrophys. J. 218, 333 Kormendy J., Richstone D. 1995, Annu. Rev. Astron. the validity of our argument. Astrophys. 33, 581 Magorrian J. et al. 1998, Astron. J. 115, 2285 Melia F. 2007, The Galactic Supermassive Black Hole. References Princeton University Press, Princeton Meylan G. et al. 2001, Astron. J. 122, 830 Bahcall J. 1989, Neutrino Astrophysics. Cambridge Univer- Miller M. C., Hamilton D. P. 2002, Mon. Not. R. Astron. sity Press, Cambridge Soc. 330, 232 Bandara K., Crampton D., Simard L. 2009, Astrophys. J. Pasham D., Strohmayer T., Mushotzky R. 2014, Nature 513, 74 704, 1135 Reines A. E. et al. 2020, Astrophys. J. 888, 36 Burrows A. 1990, Annu. Rev. Nucl. Part. Sci. 40, 181 Shapiro S., Teukolsky S. 1983, White Dwarf, Neutron Stars Chandrasekhar S. 1942, Introduction to Stellar Dynamics. and Black Holes. Wiley Verlag, Weinheim University of Chicago Press, Chicago Silk J., Rees M. J. 1998, Astron. Astrophys. 331, L1 Di Matteo T., Springel V., Hernquist L. 2005, Nature 433, Sivaram C., Arun K. 2014, Adv. High Energy Phys. Spec. 604 Issue: Black Holes Phys. 14, 924848 Faber S., Gallagher J. 1979, Annu. Rev. Astron. Astrophys. Volonteri M., Natarajan P. 2009, Mon. Not. R. Astron. Soc. 17, 135 400, 1911 Feoli A., Mancini L. 2009, Astrophys. J. 703, 1502 Zel’dovich Y. B., Novikov I. D. 1971, Relativistic Astro- Gultekin K. et al. 2009, Astrophys. J. 698, 198 physics, Part I. University of Chicago Press, Chicago