Supermassive black holes (SMBHs) are found in most galactic nuclei. A significant fraction of these nuclei also contain a nuclear stellar cluster (NSC) surrounding the SMBH. In this paper, we consider the idea that the NSC forms first, from the merger of several stellar clusters that may contain intermediate-mass black holes (IMBHs). These IMBHs can subsequently grow in the NSC and form an SMBH. We carry out $N$-body simulations of the simultaneous merger of three stellar clusters to form an NSC. We investigate the outcome of the simulated runs containing zero, one, two and three IMBHs. We find that IMBHs can efficiently sink to the centre of the merged cluster. If multiple merging clusters contain an IMBH, we find that an IMBH binary is likely to form and subsequently merge by gravitational wave emission. We show that these mergers are catalyzed by dynamical interactions with surrounding stars, which systematically harden the binary and increase its orbital eccentricity. The seed SMBH will be ejected from the NSC by the recoil kick produced when two IMBHs merge, if their mass ratio $q\gtrsim 0.15$. We find that the other channel, namely, dynamical interactions between an IMBH binary and a third single IMBH, is rarely effective in removing the seed. If the seed is ejected then no SMBH will form in the NSC. This is a natural pathway to explain those galactic nuclei that contain an NSC but apparently lack an SMBH, such as M33. However, if an IMBH is retained then it can seed the growth of an SMBH through gas accretion and tidal disruption of stars.

We study the formation rate of binary black hole mergers formed through gravitational-wave emission between unbound, single black holes in globular clusters. While the formation of these binaries in very dense systems such as galactic nuclei has been well studied, we show here that this process can operate in lower-density stellar systems as well, forming binaries at a rate similar to other proposed pathways for creating eccentric mergers. Recent advances in post-Newtonian cluster dynamics indicate that a large fraction of dynamically-assembled binary black holes merge inside their host clusters during weak and strong binary-single and binary-binary interactions, and that these systems may retain measurable eccentricities as they travel through the LIGO and LISA sensitivity bands. Using an analytic approach to modeling binary black holes from globular clusters, we show that the formation of merging binaries from previously unbound black holes can operate at a similar rate to mergers forming during strong binary encounters, and that these binaries inhabit a unique region of the gravitational-wave frequency space which can be identified by proposed deci-Hertz space-based detectors.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

We utilize multi-epoch MUSE spectroscopy to study binaries in the core of NGC 3201. Our sample consists of 3553 stars with 54883 spectra in total comprising 3200 main-sequence stars up to 4 magnitudes below the turn-off. Each star in our sample has between 3 and 63 (with a median of 14) reliable radial velocity (RV) measurements within five years of observations. We introduce a statistical method to determine the probability of a star showing RV variations based on the whole inhomogeneous RV sample. Using HST photometry and an advanced dynamical MOCCA simulation of this specific GC we overcome observational biases that previous spectroscopic studies had to deal with. This allows us to infer a binary frequency in the MUSE FoV and enables us to deduce the underlying true binary frequency of (6.75+-0.72) % in NGC 3201. The comparison of the MUSE observations with the MOCCA simulation suggests a significant fraction of primordial binaries. We can also confirm a radial increase of the binary fraction towards the GC centre due to mass segregation. We discovered that in our sample at least (57.5+-7.9) % of blue straggler stars (BSS) are in a binary system. For the first time in a study of GCs, we were able to fit Keplerian orbits to a significant sample of 95 binaries. We present the binary system properties of eleven BSS and show evidence that two BSS formation scenarios, the mass transfer in binary (or triple) star systems and the coalescence due to binary-binary interactions, are present in our data. We also describe the binary and spectroscopic properties of four sub-subgiant (or red straggler) stars. Furthermore, we discovered two new black hole (BH) candidates with minimum masses (Msini) of (7.68+-0.50) M_sun, (4.4+-2.8) M_sun, and refine the minimum mass estimate on the already published BH to (4.53+-0.21) M_sun. These BHs are consistent with an extensive BH subsystem hosted by NGC 3201.

We investigate the dissolution process for dynamically evolving star clusters embedded in an external tidal field by exploring the MOCCA Survey Database I, with focus on the presence and evolution of a stellar-mass black hole subsystem. We argue that the presence of a black hole subsystem can lead to the dissolution of tidally filling star clusters and this can be regarded as a third type of cluster dissolution mechanism (in addition to well-known mechanisms connected with strong mass loss due to stellar evolution and mass loss connected with the relaxation process). This third process is characterized by abrupt cluster dissolution connected with the loss of dynamical equilibrium. The abrupt dissolution is powered by strong energy generation from a stellar-mass black hole subsystem accompanied by tidal stripping. Additionally, we argue that such a mechanism should also work for even tidally under-filling clusters with top-heavy initial mass function. Observationally, star clusters which undergo dissolution powered by the third mechanism would look as a 'dark cluster' i.e. composed of stellar mass black holes surrounded by an expanding halo of luminous stars (Banerjee & Kroupa 2011), and they should be different from 'dark clusters' harbouring intermediate mass black holes as discussed by Askar et al. (2017a). An additional observational consequence of an operation of the third dissolution mechanism should be a larger than expected abundance of free floating black holes in the Galactic halo.

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.

Supermassive black holes are found in most galactic nuclei. A large fraction of these nuclei also contain a nuclear stellar cluster surrounding the black hole. Here we consider the idea that the nuclear stellar cluster formed first and that the supermassive black hole grew later. In particular we consider the merger of three stellar clusters to form a nuclear stellar cluster, where some of these clusters contain a single intermediate-mass black hole (IMBH). In the cases where multiple clusters contain IMBHs, we discuss whether the black holes are likely to merge and whether such mergers are likely to result in the ejection of the merged black hole from the nuclear stellar cluster. In some cases, no supermassive black hole will form as any merger product is not retained. This is a natural pathway to explain those galactic nuclei that contain a nuclear stellar cluster but apparently lack a supermassive black hole; M33 being a nearby example. Alternatively, if an IMBH merger product is retained within the nuclear stellar cluster, it may subsequently grow, e.g. via the tidal disruption of stars, to form a supermassive black hole.

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

We have been investigating populations of cataclysmic variables (CVs) in a set of more than 300 globular cluster (GC) models evolved with the MOCCA code. One of the main questions we have intended to answer is whether most CVs in GCs are dynamically formed or not. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in GCs than dynamical formation of CVs. In particular, we found that, on average, the detectable CV population is predominantly composed of CVs formed via a typical common envelope phase (> 70 per cent). However, core-collapsed models tend to have higher fractions of bright CVs than non-core-collapsed ones, which suggests then that the formation of CVs is indeed slightly favoured through strong dynamical interactions in core-collapsed GCs, due to the high stellar densities in their cores.

In this paper we explore the interplay between intermediate-mass black holes (IMBHs) and their nursing globular clusters (GCs), taking advantage of over 2000 Monte Carlo GC models. We find that the average density of IMBHs sphere of influence can be uniquely connected to the host GCs luminosity and half-light radius via a fundamental plane. We propose a statistical approach to systematically identify potential Galactic GCs harbouring either an IMBH or a massive subsystem comprised of stellar BHs. Our models show that the IMBH is often bound to a stellar companion or a stellar BH, which can lead to tidal disruption events or to low-frequency gravitational waves. We show that GCs orbiting close to the Galactic Centre have a larger probability to witness IMBH formation during their early evolution. These low-orbit GCs can deliver several IMBHs into the galaxy innermost regions, with potential impact on both electromagnetic and GW emission. We discuss potential connections between IMBHs and SMBHs inhabiting galactic nuclei, exploring the possibility that in some cases they share similar formation pathways.

We briefly describe and discuss the set-up of the project MOCCA Survey Database I. The database contains more than 2000 Monte Carlo models of evolution of real star cluster performed with the MOCCA code. Then, we very briefly discuss results of analysis of the database regarding the following projects: formation of intermediate mass black holes, abrupt cluster dissolution harboring black hole subsystems, retention fraction of black hole - black hole mergers, and tidal disruption events with intermediate mass black holes.

Machine learning is a powerful technique, becoming increasingly popular in astrophysics. In this paper, we apply machine learning to more than a thousand globular cluster (GC) models simulated as part of the 'MOCCA-Survey Database I' project in order to correlate present-day observable properties with the presence of a subsystem of stellar mass black holes (BHs). The machine learning model is then applied to available observed parameters for Galactic GCs to identify which of them that are most likely to be hosting a sizeable number of BHs and reveal insights into what properties lead to the formation of BH subsystems. With our machine learning model, we were able to shortlist 21 Galactic GCs that are most likely to contain a BH subsystem. We show that the clusters shortlisted by the machine learning classifier include those in which BH candidates have been observed (M22, M10 and NGC 3201) and that our results line up well with independent simulations and previous studies that manually compared simulated GC models with observed properties of Galactic GCs. These results can be useful for observers searching for elusive stellar mass BH candidates in GCs and further our understanding of the role BHs play in GC evolution. In addition, we have released an online tool that allows one to get predictions from our model after they input observable properties.

Machine learning is a powerful technique, becoming increasingly popular in astrophysics. In this paper, we apply machine learning to more than a thousand globular cluster (GC) models simulated as part of the 'MOCCA-Survey Database I' project in order to correlate present-day observable properties with the presence of a subsystem of stellar mass black holes (BHs). The machine learning model is then applied to available observed parameters for Galactic GCs to identify which of them that are most likely to be hosting a sizeable number of BHs and reveal insights into what properties lead to the formation of BH subsystems. With our machine learning model, we were able to shortlist 21 Galactic GCs that are most likely to contain a BH subsystem. We show that the clusters shortlisted by the machine learning classifier include those in which BH candidates have been observed (M22, M10 and NGC 3201) and that our results line up well with independent simulations and previous studies that manually compared simulated GC models with observed properties of Galactic GCs. These results can be useful for observers searching for elusive stellar mass BH candidates in GCs and further our understanding of the role BHs play in GC evolution. In addition, we have released an online tool that allows one to get predictions from our model after they input observable properties.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We investigate here populations of cataclysmic variables (CVs) in a set of 288 globular cluster (GC) models evolved with the MOCCA code. This is by far the largest sample of GC models ever analysed with respect to CVs. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in GCs than dynamical formation of CVs, and that dynamically formed CVs and CVs formed under no/weak influence of dynamics have similar white dwarf mass distributions. In addition, we found that, on average, the detectable CV population is predominantly composed of CVs formed via typical common envelope phase (CEP) ($\gtrsim70$ per cent), that only $\approx2-4$ per cent of all CVs in a GC is likely to be detectable, and that core-collapsed models tend to have higher fractions of bright CVs than non-core-collapsed ones. We also consistently show, for the first time, that the properties of bright and faint CVs can be understood by means of the pre-CV and CV formation rates, their properties at their formation times and cluster half-mass relaxation times. Finally, we show that models following the initial binary population proposed by Kroupa and set with low CEP efficiency better reproduce the observed amount of CVs and CV candidates in NGC 6397, NGC 6752 and 47 Tuc. To progress with comparisons, the essential next step is to properly characterize the candidates as CVs (e.g. by obtaining orbital periods and mass ratios).

Anisotropy of gravitational wave (GW) emission results in a net momentum gained by the black hole (BH) merger product, leading to a recoil velocity up to $\sim10^3\text{ km s}^{-1}$, which may kick it out of a globular cluster (GC). We estimate GW kick retention fractions of merger products assuming different models for BH spin magnitude and orientation (MS0 - random, MS1 - spin as a function of mass and metalicity, MS2 - constant value of $0.5$). We check how they depend on BH-BH merger time and properties of the cluster. We analyze the implications of GW kick retention fractions on intermediate massive BH (IMBH) formation by repeated mergers in a GC. We also calculate final spin of the merger product, and investigate how it correlates with effective spin of the binary. We used data from MOCCA (MOnte Carlo Cluster simulAtor) GC simulations to get a realistic sample of BH-BH mergers, assigned each BH spin value according to a studied model, and calculated recoil velocity and final spin based on most recent theoretical formulas. We discovered that for physically motivated models, GW kick retention fractions are about $30\%$ and display small dependence on assumptions about spin, but are much more prone to cluster properties. In particular, we discovered a strong dependence of GW kick retention fractions on cluster density. We also show that GW kick retention fractions are high in final life stages of the cluster, but low at the beginning. Finally, we derive formulas connecting final spin with effective spin for primordial binaries, and with maximal effective spin for dynamical binaries.

The dense environment of globular clusters (GCs) can facilitate the formation of binary black holes (BBHs), some of which can merge with gravitational waves (GW) within the age of the Universe. We have performed a survey of Monte-Carlo simulations following the dynamical evolution of GCs with different masses, sizes and binary fractions and explored the impact of the host GC properties on the formation of BBH mergers. We find that the number of BBH mergers from GCs is determined by the GC's initial mass, size and primordial binary fraction. We identify two groups of BBH mergers: a primordial group whose formation does not depend on cluster's dynamics and a dynamical group whose formation is driven by the cluster's dynamical evolution. We show how the BBH origin affects the BBH mergers' main properties such as the chirp mass and merging time distributions. We provide analytic expressions for the dependence of the number of BBH mergers from individual GCs on the main cluster's structural properties and the time evolution of the merger rates of these BBHs. These expressions provide an essential ingredient for a general framework allowing to estimate the merger rate density. Using the relations found in our study, we find a local merger rate density of 0.18$-$1.8 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for primordial BBH mergers and 0.6$-$18 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for dynamical BBH mergers, depending on the GC mass and size distributions, initial binary fraction and the number density of GCs in the Universe.

There have been increasing theoretical speculations and observational indications that certain globular clusters (GCs) could contain a sizeable population of stellar mass black holes (BHs). In this paper, we shortlist at least 29 Galactic GCs that could be hosting a subsystem of BHs (BHS). In a companion paper, we analysed results from a wide array of GC models (simulated with the MOCCA code for cluster simulations) that retained few tens to several hundreds of BHs at 12 Gyr and showed that the properties of the BHS in those GCs correlate with the GC's observable properties. Building on those results, we use available observational properties of 140 Galactic GCs to identify 29 GCs that could potentially be harbouring up to a few hundreds of BHs. Utilizing observational properties and theoretical scaling relations, we estimate the density, size and mass of the BHS in these GCs. We also calculate the total number of BHs and the fraction of BHs contained in a binary system for our shortlisted Galactic GCs. Additionally, we mention other Galactic GCs that could also contain significant number of single BHs or BHs in binary systems.

In this paper, we discuss how globular clusters (GCs) structural and observational properties can be used to infer the presence of a black hole system (BHS) inhabiting their inner regions. We propose a novel way to identify the BHS size, defined as the GC radius containing a mass contributed equally from stars and stellar BHs. Using this definition, similar to the well-known concept of "influence radius", we found a "fundamental plane" connecting the BHS typical density with the GC central surface density profile, total luminosity and observational half-mass radius. Our approach allows us to define a unique way to connect the observational GCs parameters with their dark content. Comparing our results with observed Milky Way GCs, we found that many of them likely host, at the present time, as many as several hundreds of BHs. These BHS are characterized by a relatively low typical density, $\rho_\bhs \sim 10-10^5\Ms$ pc$^{-3}$ and composed of relatively massive BHs, with average masses in the range $m_\bhs = 14-22\Ms$. We also show that a similar approach can be used to find Milky Way GCs potentially hosting an intermediate-mass black hole.

The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.

We study the gravitational wave (GW) frequency and chirp mass distribution of binary black hole (BBH) mergers assembled through three-body interactions in globular clusters (GCs), when GW emission at the 2.5 post-Newtonian (PN) level is included in the $N$-body equation-of-motion (EOM). From performing $\sim 2.5\times10^{6}$ PN binary-single interactions based on GC data from the `MOCCA-Survey Database I' project, and by the use of analytical methods, we find that $5-10\%$ of all the three-body assembled GC BBH mergers have a GW frequency at formation that is $\gtrsim 10^{-1}$ Hz, implying they enter the LIGO band without having drifted through the LISA band first. If PN terms are not included in the EOM, one finds instead that all BBH mergers drifts through both LISA and LIGO. As the fraction of BBH mergers that only show up in LIGO is expected to be $\sim 0\%$ for standard field binary BBH mergers, future joint measurements with LISA and LIGO can be used to gain insight into the formation of BBH mergers.

We estimate the population of eccentric gravitational wave (GW) binary black hole (BBH) mergers forming during binary-single interactions in globular clusters (GCs), using ~ 800 GC models that were evolved using the MOCCA code for star cluster simulations as part of the MOCCA-Survey Database I project. By re-simulating binary-single interactions (only involving 3 BHs) extracted from this set of GC models using an N-body code that includes GW emission at the 2.5 post-Newtonian level, we find that ~ 10% of all the BBHs assembled in our GC models that merge at present time form during chaotic binary-single interactions, and that about half of this sample have an eccentricity > 0.1 at 10 Hz. We explicitly show that this derived rate of eccentric mergers is ~ 100 times higher than one would find with a purely Newtonian N-body code. Furthermore, we demonstrate that the eccentric fraction can be accurately estimated using a simple analytical formalism when the interacting BHs are of similar mass; a result that serves as the first successful analytical description of eccentric GW mergers forming during three-body interactions in realistic GCs.

We introduce and present results from the COCOA (Cluster simulatiOn Comparison with ObservAtions) code that has been developed to create idealized mock photometric observations using results from numerical simulations of star cluster evolution. COCOA is able to present the output of realistic numerical simulations of star clusters carried out using Monte Carlo or \textit{N}-body codes in a way that is useful for direct comparison with photometric observations. In this paper, we describe the COCOA code and demonstrate its different applications by utilizing globular cluster (GC) models simulated with the MOCCA (MOnte Carlo Cluster simulAtor) code. COCOA is used to synthetically observe these different GC models with optical telescopes, perform PSF photometry and subsequently produce observed colour magnitude diagrams. We also use COCOA to compare the results from synthetic observations of a cluster model that has the same age and metallicity as the Galactic GC NGC 2808 with observations of the same cluster carried out with a 2.2 meter optical telescope. We find that COCOA can effectively simulate realistic observations and recover photometric data. COCOA has numerous scientific applications that maybe be helpful for both theoreticians and observers that work on star clusters. Plans for further improving and developing the code are also discussed in this paper.

Numerical simulations were carried out to study the origin of multiple stellar populations in the intermediate-age clusters NGC 411 and NGC 1806 in the Magellanic Clouds. We performed NBODY6++ simulations based on two different formation scenarios, an ad hoc formation model where second-generation (SG) stars are formed inside a cluster of first-generation (FG) stars using the gas accumulated from the external intergalactic medium and a minor merger model of unequal mass ($M_{SG}$/$M_{FG}$ ~5-10%) clusters with an age difference of a few hundred million years. We compared our results such as the radial profile of the SG-to-FG number ratio with observations on the assumption that the SG stars in the observations are composed of cluster members, and confirmed that both the ad hoc formation and merger scenarios reproduce the observed radial trend of the SG-to-FG number ratio which shows less centrally concentrated SG than FG stars. It is difficult to constrain the formation scenario for the multiple populations by only using the spatial distribution of the SG stars. SG stars originating from the merger scenario show a significant velocity anisotropy and rotational features compared to those from the ad hoc formation scenario. Thus, observations aimed at kinematic properties like velocity anisotropy or rotational velocities for SG stars should be obtained to better understand the formation of the multiple populations in these clusters. This is, however, beyond current instrumentation capabilities.

Colour-magnitude diagrams (CMDs) are powerful tools that might be used to infer stellar properties in globular clusters (GCs), for example, the binary fraction and their mass ratio ($q$) distribution. In the past few years, observations have revealed that q distributions of GC main-sequence (MS) binaries are generally flat, and a distribution characterized by a strong increase towards q $\approx$ 1 is not typical in GCs. In numerical simulations of GC evolution with the initial binary population (IBP) described by Kroupa, synthetic CMD colour distributions exhibit a peak associated with binaries that have q $\approx$ 1. While the Kroupa IBP reproduces binary properties in star-forming regions, clusters and the Galactic field, the peak in the q distribution towards q $\approx$ 1 observed for GC simulations is not consistent with distributions derived from observations. The objective of this paper is to refine and further improve the physical formulation of pre-main-sequence eigenevolution proposed by Kroupa in order to achieve CMD colour distributions of simulated GC models similar to those observed in real GCs, and to get a similarly good agreement with binary properties for late-type binaries in the Galactic field. We present in this paper a modified Kroupa IBP, in which early-type stars follow observational distributions, and late-type stars are generated according to slightly modified pre-main-sequence eigenevolution prescriptions. Our modifications not only lead to a qualitatively good agreement with respect to long-term observations of late-type binaries in the Galactic field, but also resolve the above-mentioned problem related to binary distributions in GC models.

In this third of a series of papers related to cataclysmic variables (CVs) and related objects, we analyse the population of CVs in a set of 12 globular cluster models evolved with the MOCCA Monte Carlo code, for two initial binary populations (IBPs), two choices of common-envelope phase (CEP) parameters, and three different models for the evolution of CVs and the treatment of angular momentum loss. When more realistic models and parameters are considered, we find that present-day cluster CV duty cycles are extremely-low ($\lesssim 0.1$ per cent) which makes their detection during outbursts rather difficult. Additionally, the IBP plays a significant role in shaping the CV population properties, and models that follow the Kroupa IBP are less affected by enhanced angular momentum loss. We also predict from our simulations that CVs formed dynamically in the past few Gyr (massive CVs) correspond to bright CVs (as expected), and that faint CVs formed several Gyr ago (dynamically or not) represent the overwhelming majority. Regarding the CV formation rate, we rule out the notion that it is similar irrespective of the cluster properties. Finally, we discuss the differences in the present-day CV properties related to the IBPs, the initial cluster conditions, the CEP parameters, formation channels, the CV evolution models, and the angular momentum loss treatments.

The detection of intermediate mass black holes (IMBHs) in Galactic globular clusters (GCs) has so far been controversial. In order to characterize the effectiveness of integrated-light spectroscopy through integral field units, we analyze realistic mock data generated from state-of-the-art Monte Carlo simulations of GCs with a central IMBH, considering different setups and conditions varying IMBH mass, cluster distance, and accuracy in determination of the center. The mock observations are modeled with isotropic Jeans models to assess the success rate in identifying the IMBH presence, which we find to be primarily dependent on IMBH mass. However, even for a IMBH of considerable mass (3% of the total GC mass), the analysis does not yield conclusive results in 1 out of 5 cases, because of shot noise due to bright stars close to the IMBH line-of-sight. This stochastic variability in the modeling outcome grows with decreasing BH mass, with approximately 3 failures out of 4 for IMBHs with 0.1% of total GC mass. Finally, we find that our analysis is generally unable to exclude at 68% confidence an IMBH with mass of $10^3~M_\odot$ in snapshots without a central BH. Interestingly, our results are not sensitive to GC distance within 5-20 kpc, nor to mis-identification of the GC center by less than 2'' (<20% of the core radius). These findings highlight the value of ground-based integral field spectroscopy for large GC surveys, where systematic failures can be accounted for, but stress the importance of discrete kinematic measurements that are less affected by stochasticity induced by bright stars.

This is the second in a series of papers associated with cataclysmic variables (CVs) and related objects, formed in a suite of simulations for globular cluster evolution performed with the MOCCA Monte Carlo code. We study the properties of our simulated CV populations throughout the entire cluster evolution. We find that dynamics extends the range of binary CV progenitor properties, causing CV formation from binary progenitors that would otherwise not become CVs. The CV formation rate in our simulations can be separated into two regimes: an initial burst ($\lesssim$ 1 Gyr) connected with the formation of the most massive white dwarfs, followed by a nearly constant formation rate. This result holds for all models regardless of the adopted initial conditions, even when most CVs form dynamically. Given the cluster age-dependence of CV properties, we argue that direct comparisons to observed Galactic field CVs could be misleading, since cluster CVs can be up to 4 times older than their field counterparts. Our results also illustrate that, due mainly to unstable mass transfer, some CVs that form in our simulations are destroyed before the present-day. Finally, some field CVs might have originated from globular clusters, as found in our simulations, although the fraction of such escapers should be small relative to the entire Galactic field CV population.

We describe the dynamical evolution of a unique type of dark star cluster model in which the majority of the cluster mass at Hubble time is dominated by an intermediate-mass black hole (IMBH). We analyzed results from about 2000 star cluster models (Survey Database I) simulated using the Monte-Carlo code MOCCA and identified these dark star cluster models. Taking one of these models, we apply the method of simulating realistic "mock observations" by utilizing the COCOA and SISCO codes to obtain the photometric and kinematic observational properties of the dark star cluster model at 12 Gyr. We find that the perplexing Galactic globular cluster NGC 6535 closely matches the observational photometric and kinematic properties of the dark star cluster model presented in this paper. Based on our analysis and currently observed properties of NGC 6535, we suggest that this globular cluster could potentially harbour an IMBH. If it exists, the presence of this IMBH can be detected robustly with proposed kinematic observations of NGC 6535.

In this first of a series of papers, we utilize results for around two thousand star cluster models simulated using the MOCCA code for star cluster evolution (Survey Database I) to determine the astrophysical properties and local merger rate densities for coalescing binary black holes (BBHs) originating from globular clusters (GCs). We extracted information for all coalescing BBHs that escape the cluster models and subsequently merge within a Hubble time along with BBHs that are retained in our GC models and merge inside the cluster via gravitational wave (GW) emission. By obtaining results from a substantial number of realistic star cluster models that cover different initial parameters, we have an extremely large statistical sample of BBHs with stellar mass and massive stellar BH ($\lesssim 100M_{\odot}$) components that merge within a Hubble time. Using this data, we estimate local merger rate densities for these BBHs originating from GCs to be at least 5.4 ${\rm Gpc}^{-3}\,{\rm yr}^{-1}$

In this paper, which is the first in a series of papers associated with cataclysmic variables and related objects, we introduce the CATUABA code, a numerical machinery written for analysis of the MOCCA simulations, and show some first results by investigating the present-day population of cataclysmic variables in globular clusters. Emphasis was given on their properties and the observational selection effects when observing and detecting them. In this work we analysed in this work six models, including three with Kroupa distributions of the initial binaries. We found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place. We show and explain why this is inconsistent with observational and theoretical results. Our results indicate that the population of cataclysmic variables in globular clusters is, mainly, in the last stage of their evolution and observational selection effects can drastically change the expected number of observed cataclysmic variables. We show that the probability of observing them during the outbursts is extremely small and conclude that the best way of looking for cataclysmic variables in globular clusters is by searching for variabilities during quiescence, instead of during outbursts. For that, one would need a very deep observation which could reach magnitudes $\gtrsim$ 27 mag. Finally, we argue that cataclysmic variables in globular clusters are not necessarily magnetic.

We introduce and present preliminary results from COCOA (Cluster simulatiOn Comparison with ObservAtions) code for a star cluster after 12 Gyrs of evolution simulated using the MOCCA code. The COCOA code is being developed to quickly compare results of numerical simulations of star clusters with observational data. We use COCOA to obtain parameters of the projected cluster model. For comparison, a FITS file of the projected cluster was provided to observers so that they could use their observational methods and techniques to obtain cluster parameters. The results show that the similarity of cluster parameters obtained through numerical simulations and observations depends significantly on the quality of observational data and photometric accuracy.

We discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the influence of an intermediate mass black hole on the host global globular cluster properties. We discuss two scenarios for intermediate mass black hole formation. The SLOW and FAST scenarios. They occur later or earlier in the cluster evolution and require smaller or extremely large central densities, respectively. In our simulations, the formation of intermediate mass black holes is highly stochastic. In general, higher formation probabilities follow from larger cluster concentrations (i.e. central densities). We further discuss possible observational signatures of the presence of intermediate mass black holes in globular clusters that follow from our simulations. These include the spatial and kinematic structure of the host cluster, possible radio, X-ray and gravitational wave emissions due to dynamical collisions or mass-transfer and the creation of hypervelocity main sequence escapers during strong dynamical interactions between binaries and an intermediate mass black hole. All simulations discussed in this paper were performed with the MOCCA Monte Carlo code. MOCCA accurately follows most of the important physical processes that occur during the dynamical evolution of star clusters but, as with other dynamical codes, it approximates the dissipative processes connected with stellar collisions and binary mergers.

We very briefly discuss proposed in the literature possible scenarios for intermediate mass black holes (IMBH) formation in globular clusters. We also discuss the results of the MOCCA simulations of about 2000 models (BigSurvey) regarding the distribution of events connected with electromagnetic and gravitational radiations, namely: mass transfer on IMBH, collisions and mergers with IMBH and mergers with IMBH due to gravitational radiation. The rates of these events are very small, so their observation is very improbable.

In this first investigation of the MOCCA database with respect to cataclysmic variables, we found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common-envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place. Our results also indicate that the population of cataclysmic variables in globular clusters are, mainly, in the last stage of their evolution and observational selection effects can change drastically the expected number and properties of observed cataclysmic variables.

Introducing the DRAGON simulation project, we present direct $N$-body simulations of four massive globular clusters (GCs) with $10^6$ stars and 5$\%$ primordial binaries at a high level of accuracy and realism. The GC evolution is computed with NBODY6++GPU and follows the dynamical and stellar evolution of individual stars and binaries, kicks of neutron stars and black holes, and the effect of a tidal field. We investigate the evolution of the luminous (stellar) and dark (faint stars and stellar remnants) GC components and create mock observations of the simulations (i.e. photometry, color-magnitude diagrams, surface brightness and velocity dispersion profiles). By connecting internal processes to observable features we highlight the formation of a long-lived 'dark' nuclear subsystem made of black holes (BHs), which results in a two-component structure. The inner core is dominated by the BH subsystem and experiences a core collapse phase within the first Gyr. It can be detected in the stellar (luminous) line-of-sight velocity dispersion profiles. The outer extended core - commonly observed in the (luminous) surface brightness profiles - shows no collapse features and is continuously expanding. We demonstrate how a King (1966) model fit to observed clusters might help identify the presence of post core-collapse BH subsystems. For global observables like core and half-mass radii the direct simulations agree well with Monte-Carlo models. Variations in the initial mass function can result in significantly different GC properties (e.g. density distributions) driven by varying amounts of early mass loss and the number of forming BHs.

We will discuss the evolution of star clusters with an large initial binary fraction, up to 95%. The initial binary population is chosen to follow the invariant orbital-parameter distributions suggested by Kroupa (1995). The Monte Carlo MOCCA simulations of star cluster evolution are compared to the observations of Milone et al. (2012) for photometric binaries. It is demonstrated that the observed dependence on cluster mass of both the binary fraction and the ratio of the binary fractions inside and outside of the half mass radius are well recovered by the MOCCA simulations. This is due to a rapid decrease in the initial binary fraction due to the strong density-dependent destruction of wide binaries described by Marks, Kroupa & Oh (2011). We also discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We will discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the dependence of its mass accretion rate on the global cluster properties.