## MOCCA-SURVEY Database I: Dissolution of tidally filling star clusters harboring black hole subsystem

Mirek Giersz , Abbas Askar , Long Wang 1 2 3 4 5 , , , Arkadiusz Hypki , 6 Agostino Leveque 1 and Rainer Spurzem 7 8 9 , ,

1 Nicolaus Copernicus Astronomical Centre, Polish Academy of Sciences, Warsaw, Poland email: mig@camk.edu.pl

2 Lund Observatory, Department of Astronomy, and Theoretical Physics, Lund University, Box 43, SE-221 00 Lund, Sweden email: askar@astro.lu.se

3 Helmholtz-Institut f¨r Strahlen- und Kernphysik, University of Bonn, Nussallee 14-16, u D-53115 Bonn, Germany

5

4 Argelander Institut F¨r Astronomie, Auf Dem H¨gel 71, 53121, Bonn, Germany u u

RIKEN Center for Computational Science, 7-1-26 Minatojima-minami-machi, Chuo-ku,

Kobe, Hyogo 650-0047, Japan email:

longwang.astro@live.com

6 Astronomical Observatory Institute, Faculty of Physics, A. Mickiewicz University, S/suppress loneczna 36, 60-286 Pozna´, Poland n email: ahypki@amu.edu.pl

7

National Astronomical Observatories and Key Laboratory for Computational Astrophysics, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing 100012, China 8 Astronomisches Rechen-Institut, Zentrum f¨r Astronomie, University of Heidelberg, u Monchhofstrasse 12-14, D-69120 Heidelberg, Germany 9 Kavli Institute for Astronomy and Astrophysics, Peking University, Yi He Yuan Lu 5, HaiDian District, Beijing 100871, China email: spurzem@ari.uni-heidelberg.de

Abstract. We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium. Such cluster dissolution can be regarded as a third type of cluster dissolution mechanism. We additionally argue that such a mechanism should also work for tidally under-filling clusters with a top-heavy initial mass function.

Keywords. stellar dynamics, methods: numerical, globular clusters: general, stars: black holes

## 1. Introduction

There are generally two types of cluster dissolution mechanisms discussed in the literature. One, slow, connected with relaxation-driven mass loss, and second, abrupt, connected with stellar evolution mass loss (e.g., Contenta et al. 2015). Fukushige &amp; Heggie (1995) attributed the final abrupt cluster dissolution to the loss of dynamical equilibrium within the cluster. In such a situation, the cluster will not undergo core collapse. Many authors pointed out additional physical processes (apart from relaxation or stellar evolution), which can influence the way in which a cluster is dissolved. In the 90s, Chernoff

&amp; Weinberg (1990) showed that star cluster structure, and the tidal field of a parent galaxy have a very strong influence on cluster evolution and its dissolution time. This was a seminal study, which later resulted in many more detailed papers discussing the final state of star clusters. It was shown that several other factors can also influence the lifetime of star clusters, such as: crossing time-scale (Whitehead et al. 2013), primordial binary population (Tanikawa &amp; Fukushige 2009), type of the Galactic orbit (Baumgardt &amp; Makino 2003), the form of the Galactic potential and tidal shocking (Gnedin &amp; Ostriker 1997), and the properties of dark remnants retained in star clusters (Banerjee &amp; Kroupa 2011 and Contenta et al. 2015). Inspired by those works, we decided to look at the cluster dissolution process in the models from the MOCCA Survey Database I (Askar et al. 2017), models of real star clusters simulated by the MOCCA code (Giersz et al. 2013 and Hypki &amp; Giersz 2013). We introduced a third mechanism for cluster dissolution, which operates in tidally filling star clusters harboring a black hole subsystem (BHS). The physical processes responsible for the abrupt cluster dissolution are connected with strong energy generation in dynamical interactions between black holes (BH) and BHBH binaries (BBHs) and strong cluster tidal stripping. The details of the third cluster dissolution mechanism are described by Giersz et al. (2019). Here, we will only provide short a summary and point out for the most important ingredients of the process of cluster dissolution.

## Wo = 6, tidally filling

Figure 1. Evolution of the fraction of cluster bound mass as a function of time for the MOCCA tidally filling models with King parameter, W 0 , and different number of stars, N , tidal radii, R t , mass FALLBACK ON or OFF. N-body model for N = 100000 is depicted with thick black line. Description of models and lines are given inside the figure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00266v1-figure-7d939963-1545-49c5-8ce0-98402360f17c.jpg)

## 2. Results

Cluster dissolution powered by stellar evolution mass loss coupled with strong cluster expansion and overflow of the tidal radius, R t , is characterized by the abrupt loss of

equilibrium and very fast dissolution. Usually, such clusters do not enter the core collapse phase of evolution. Dynamical evolution of models which are more concentrated is controlled by the relaxation process, not by stellar evolution mass loss. They evolve much more slowly and can enter the core collapse phase, provided that the relaxation was fast enough. The third cluster dissolution mechanism is a mixture of the above two processes - cluster enters the core collapse and post-collapse phases of evolution, and yet undergoes the abrupt final dissolution. The cluster core collapse is connected with strong mass segregation and BHS formation, and is followed, for usually a few Gyrs, by the balanced evolution between the BHS and the rest of the cluster (Breen &amp; Heggie 2013), and then abrupt dissolution.

Fig. 1 shows evolution of clusters which can form BHS (SNe kick velocities are reduced by mass fallback according to Belczynski et al. (2002) prescription) and clusters which do not form BHS because of high SNe kick velocities. Models have different numbers of objects, different tidal radii and binary fractions. Additionally, it shows an N-body model † to check if MOCCA models properly account for the final stages of the cluster evolution. The N-body model shows a very similar shape of abrupt cluster dissolution as the MOCCA models, which clearly confirms that the abrupt dissolution of clusters observed in the MOCCA simulations are real features of cluster evolution.

It is clear from Fig.1 that the abrupt dissolution does not depend on the global cluster properties, like tidal radius, cluster concentration, binary fraction or number of objects. They only determine the speed of the cluster evolution, as the relaxation time scale directly depends on those parameters. Generally, the larger the number of objects or the larger the binary fraction or the smaller the cluster concentration, the longer the abrupt dissolution time. It is worth pointing out that the observed behavior of cluster dissolution happens only for tidally filling models (the ratio between tidal and half-mass radii is typically about 7 for King model with W 0 = 6), models which do not have the room to freely expand. This, together with the fact that BHS is present until the final stages of cluster dissolution suggests that BHS and strong tidal stripping have a crucial role in destroying the balanced cluster evolution and lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium.

As it is clearly shown by Giersz et al. (2019) cluster dissolution is connected with the decoupling of the BHS from the rest of the cluster and the loss of dynamical equilibrium by other objects. Luminous objects become hot and expand. BHS stops to collapse further and slowly starts to expand and disrupts itself like an isolated cluster. This picture further supports the evolution of the escape velocity and the escape rate for tidally filing models. The cluster is constantly tidally stripping, but the internal cluster structure is not strongly changing, so the escape velocity becomes smaller and smaller, and more and more objects can be kicked out from the cluster. The positive feedback between strong energy generation, high flow of escaping stars, decreasing escape velocity and enhanced tidal stripping leads finally to the loss of cluster dynamical equilibrium and abrupt dissolution.

The above described third mechanism of cluster dissolution is universal and should work for all clusters which harbor BHSs to the end of their life, and are exposed to a strong tidal field, regardless of the initial cluster structure or surrounding environment.

† N-body simulation was done by Long Wang with the use of the Nbody6++GPU code (Wang et al. 2015) with initial conditions exactly the same as for the MOCCA model

## 3. Conclusions

Tidally filling clusters with BHS in the center can show abrupt dissolution, provided that the initial King concentration parameter is moderate and BHS can survive until cluster death. The dissolution is controlled by a strong energy generation in dynamical interactions inside the BHS and strong tidal stripping, which leads finally to sudden loss of dynamical equilibrium and decoupling of BHS from the rest of the cluster. Such mechanism should operate also for tidally under-filling clusters with top-heavy initial mass functions.

Just before the abrupt dissolution the cluster will look like a 'dark cluster' as described by Banerjee &amp; Kroupa (2011), which is a different kind of 'dark cluster', which harbors an IMBH with mass comparable to the cluster total mass (Askar et al. 2017).

The fast dissolution of massive, tidally filling clusters can have strong influence on the estimated rate of BH-BH mergers from GCs. If clusters are born tidally filling and close to the Galactic center, then we can expect a lot of free floating BHs in the Galactic Bulge.

## Acknowledgements

MG and AL were partially supported by the Polish National Science Center (NCN) through the grant UMO-2016/23/B/ST9/02732. AA is currently supported by the Carl Tryggers Foundation for Scientific Research through the grant CTS 17:113. LW in this work was supported by the funding from Alexander von Humboldt Foundation. AH was supported by Polish National Science Center grant 2016/20/S/ST9/00162. RS has been supported by the National Natural Science Foundation of China (NSFC), grant 11673032. We acknowledge support from the Chinese Academy of Sciences through the Silk Road Project at the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), through the 'Qianren' special foreign experts programme of China, and through the Sino-German collaboration program GZ1284. This work benefited from support by the International Space Science Institute (ISSI), Bern, Switzerland, through its International Team programme ref. no. 393 The Evolution of Rich Stellar Populations &amp; BH Binaries (2017-18).

## References

Askar, A., Szkudlarek, M., Gondek-Rosi´ska, D., Giersz, M., &amp; Bulik, T. 2017, n MNRAS 464, L36

Banerjee, S. &amp; Kroupa, P. 2011, ApJ 741, L12

Baumgardt, H. &amp; Makino, J. 2003, MNRAS 340, 227

Belczynski, K., Kalogera, V. &amp; Bulik, T. 2002, ApJ 572, 407

Breen, P. G. &amp; Heggie, D. C. 2013, MNRAS 432, 2779

Chernoff, D. F. &amp; Weinberg, M. D. 1990, ApJ 351, 121

Contenta, F., Varri, A. L., &amp; Heggie, D. C. 2015, MNRAS 449, L100

Fukushige, T., &amp; Heggie, D. C. 1995, MNRAS 276, 206

Giersz, M., Heggie, D. C., Hurley, J. R. &amp; Hypki, A. 2013, MNRAS 431, 2184

Giersz, M., Askar, A., Wang, L., Hypki, A., Laveque, A., &amp; Spurzem, R. 2019, MNRAS 487, 2417

Gnedin O. Y. &amp; Ostriker, J. P. 1997, ApJ 474, 223

Hypki, A. &amp; Giersz, M. 2013, MNRAS 429, 1221

Tanikawa, A. &amp; Fukushige T. 2009, PASJ 61, 721

Wang, L., Spurzem, R., Aarseth, S., Nitadori, K., Berczik, P., Kouwenhoven, M. B. N. &amp; Naab, T. 2015, MNRAS 450, 4070

Whitehead, A. J., McMillan, S. L. W., Vesperini, E., &amp; Portegies Zwart, S. 2013, ApJ 778, 118
Mirek Giersz Abbas Askar Long Wang Arkadiusz Hypki Agostino Leveque Rainer Spurzem
MOCCA-SURVEY Database I: Dissolution of tidally filling star clusters harboring black hole subsystem
We investigate the dissolution process of star clusters embedded in an
external tidal field and harboring a subsystem of stellar-mass black hole. For
this purpose we analyzed the MOCCA models of real star clusters contained in
the Mocca Survey Database I. We showed that the presence of a stellar-mass
black hole subsystem in tidally filling star cluster can lead to abrupt cluster
dissolution connected with the loss of cluster dynamical equilibrium. Such
cluster dissolution can be regarded as a third type of cluster dissolution
mechanism. We additionally argue that such a mechanism should also work for
tidally under-filling clusters with a top-heavy initial mass function.