Showing 1–50 of 97 results for author: Giersz, M

MOCCA: Global properties of tidally filling and underfilling globular star clusters with multiple stellar populations

Authors: Arkadiusz Hypki, Enrico Vesperini, Mirek Giersz, Jongsuk Hong, Abbas Askar, Magdalena Otulakowska-Hypka, Lucas Hellstrom, Grzegorz Wiktorowicz

Abstract: We explore the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that present-day properties (core and half-light radii, ratio of the number of second-generation (SG) stars to the total number of stars, NSG/NTOT) of these models cover the observed values o… ▽ More We explore the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that present-day properties (core and half-light radii, ratio of the number of second-generation (SG) stars to the total number of stars, NSG/NTOT) of these models cover the observed values of these quantities for Milky Way GCs. Starting with a relatively small value of the SG fraction (NSG/NTOT ~ 0.25) and a SG system concentrated in the inner regions of the cluster, we find, in agreement with previous studies, that systems in which the first-generation (FG) is initially tidally filling or slightly tidally underfilling best reproduce the observed ratios of NSG/NTOT and have values of the core and half-light radii typical of those of many Galactic globular clusters. Models in which the FG is initially tidally underfilling retain values of NSG/NTOT close to their initial values. These simulations expand previous investigations and serve to further constrain the viable range of initial parameters and better understand their influence on present-day GC properties. The results of this investigation also provide the basis for our future survey aimed at building specific models to reproduce the observed trends (or lack thereof) between the properties of multiple stellar populations and other clusters properties. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: Accepted to Astronomy and Astrophysics (A&A)

Probing Intracluster Dynamics and Evolution of Globular Clusters through Cataclysmic Variable Populations

Authors: Kwangmin Oh, Jongsuk Hong, C. Y. Hui, Sangin Kim, Mirek Giersz

Abstract: Dynamical interactions in globular clusters (GCs) significantly impact the formation and evolution of binary sources, including cataclysmic variables (CVs). This study investigates the connection between dynamical states of GCs and X-ray luminosity ($L_{x}$) distributions of CV populations through both simulations and actual observations. Utilizing a Monte Carlo simulation tool, MOCCA, we categori… ▽ More Dynamical interactions in globular clusters (GCs) significantly impact the formation and evolution of binary sources, including cataclysmic variables (CVs). This study investigates the connection between dynamical states of GCs and X-ray luminosity ($L_{x}$) distributions of CV populations through both simulations and actual observations. Utilizing a Monte Carlo simulation tool, MOCCA, we categorize the simulated GCs into three different evolutionary stages which are referred to as Classes I/II/III. Significant differences are found in the $L_{x}$ distributions of the CVs among these three Classes. In observational aspects, we have analyzed 179 CV candidates in 18 GCs observed by the {\it Chandra} X-ray Observatory. By dividing these GCs into three Families of different dynamical ages, namely Families I/II/III, the $L_{x}$ distributions of the CV candidates also show significant differences among these three Families. Both simulations and observational results suggest that CVs in more dynamically evolved clusters (Class/Family III) exhibit brighter X-ray emission. This highlights the influence of the dynamical status of a GC on the properties of its hosted compact binaries. Similar to blue stragglers, CV populations can serve as tracers of a GC's dynamical history. Our findings provide insights for understanding the interplay between intracluster dynamics and the evolution of compact binaries in GCs. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 11 pages, 9 figures, 4 tables, Accepted for MNRAS

Double white dwarf binary population in MOCCA star clusters -- Comparisons with observations of close and wide binaries

Authors: Lucas Hellström, Mirosław Giersz, Arkadiusz Hypki, Diogo Belloni, Abbas Askar, Grzegorz Wiktorowicz

Abstract: There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct… ▽ More There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct initial populations: one dense and another less dense. Due to the lower density of one population, a large number of objects escape during the early GC evolution, leading to a high mass-loss rate. In this dataset we have analysed three main groups of DWDs, namely in-cluster binaries, escaped binaries, and isolated evolution of primordial binaries. We compared the properties of these groups to observations of close and wide binaries. We find that the number of escaping DWDs is significantly larger than the number of in-cluster binaries and those that form via the isolated evolution of all promiridial binaries in our GC models. This suggests that dynamics play an important role in the formation of DWDs. For close binaries, we found a good agreement in the separations of escaped binaries and isolated binaries, but in-cluster binaries showed slight differences. We could not reproduce the observed extremely low mass WDs due to the limitations of our stellar and binary evolution prescriptions. For wide binaries, we also found a good agreement in the separations and masses, after accounting for observational selection effects. We conclude that, even though the current observational samples of DWDs are extremely biased and incomplete, our results compare reasonably well with observations. △ Less

Submitted 9 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Discovery of a hot post-AGB star in Galactic globular cluster E3

Authors: R. Kumar, A. Moharana, S. Piridi, A. C. Pradhan, K. G. Hełminiak, N. Ikonnikova, A. Dodin, R. Szczerba, M. Giersz, D. K. Ojha, M. R. Samal

Abstract: We report a new hot post-asymptotic giant branch (PAGB) star in the Galactic globular cluster (GC) E3, which is one of the first of the identified PAGB stars in a GC to show a binary signature. The star stands out as the brightest source in E3 in the \mbox{\em Astrosat}/UVIT images. We confirmed its membership with the cluster E3 using Gaia DR3 kinematics and parallax measurements. We supplemented… ▽ More We report a new hot post-asymptotic giant branch (PAGB) star in the Galactic globular cluster (GC) E3, which is one of the first of the identified PAGB stars in a GC to show a binary signature. The star stands out as the brightest source in E3 in the \mbox{\em Astrosat}/UVIT images. We confirmed its membership with the cluster E3 using Gaia DR3 kinematics and parallax measurements. We supplemented the photometric observations with radial velocities (RVs) from high-resolution spectroscopic observations at two epochs and with ground- and space-based photometric observations from 0.13 $μ$m to 22 $μ$m. We find that the RVs vary over $\sim$6 \kms\ between the two epochs. This is an indication of the star being in a binary orbit. A simulation of possible binary systems with the observed RVs suggests a binary period of either 39.12 days or 17.83 days with mass ratio q$\geq$1.0. The [Fe/H] derived using the high-resolution spectra is $\sim -$0.7 dex, which closely matches the cluster metallicity. The spectroscopic and photometric measurements suggest \Teff\ and $\log g$ of the star as 17\,500$\pm$1\,000~K and 2.37$\pm$0.20~dex, respectively. Various PAGB evolutionary tracks on the Hertzsprung--Russell (H-R) diagram suggest a current mass of the star in the range 0.51$-$0.55 \Msun. The star is enriched with C and O abundances, showing similar CNO abundances compared to the other PAGB stars in GCs with the evidence of the third dredge-up on the AGB phase. △ Less

Submitted 29 April, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Published in A&A Letters, Published Version

Journal ref: A&A 685, L6 (2024)

The Dragon-II simulations -- III. Compact binary mergers in clusters with up to 1 million stars: mass, spin, eccentricity, merger rate and pair instability supernovae rate

Authors: Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Francesco Paolo Rizzuto, Mirek Giersz, Thorsten Naab, Peter Berczik

Abstract: Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave (GW) sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the \textsc{Dragon-II} simulation database, a suite of 19 direct $N$-body m… ▽ More Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave (GW) sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the \textsc{Dragon-II} simulation database, a suite of 19 direct $N$-body models representing dense star clusters with up to $10^6$ stars and $<33\%$ of stars in primordial binaries. We find a substantial population of black hole binary (BBH) mergers, some of them involving an intermediate-mass BH (IMBH), and a handful mergers involving a stellar BH and either a neutron star (NS) or a white dwarf (WD). Primordial binary mergers, $\sim 30\%$ of the whole population, dominate ejected mergers. Dynamical mergers, instead, dominate the population of in-cluster mergers and are systematically heavier than primordial ones. Around $20\%$ of \textsc{Dragon-II} mergers are eccentric in the LISA band and $5\%$ in the LIGO band. We infer a mean cosmic merger rate of $\mathcal{R}\sim 12(4.4)(1.2)$ yr$^{-1}$ Gpc$^3$ for BBHs, NS-BH, and WD-BH binary mergers, respectively, and discuss the prospects for multimessenger detection of WD-BH binaries with LISA. We model the rate of pair-instability supernovae (PISNe) in star clusters and find that surveys with a limiting magnitude $m_{\rm bol}=25$ can detect $\sim 1-15$ yr$^{-1}$ PISNe. Comparing these estimates with future observations could help to pin down the impact of massive star evolution on the mass spectrum of compact stellar objects in star clusters. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 22 pages, 14 figures, 3 tables. Comments welcome. Submitted to MNRAS

The Dragon-II simulations -- II. Formation mechanisms, mass, and spin of intermediate-mass black holes in star clusters with up to 1 million stars

Authors: Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Francesco Paolo Rizzuto, Mirek Giersz, Thorsten Naab, Peter Berczik

Abstract: The processes that govern the formation of intermediate-mass black holes (IMBHs) in dense stellar clusters are still unclear. Here, we discuss the role of stellar mergers, star-BH interactions and accretion, as well as BH binary (BBH) mergers in seeding and growing IMBHs in the \textsc{Dragon-II} simulation database, a suite of 19 direct $N$-body models representing dense clusters with up to… ▽ More The processes that govern the formation of intermediate-mass black holes (IMBHs) in dense stellar clusters are still unclear. Here, we discuss the role of stellar mergers, star-BH interactions and accretion, as well as BH binary (BBH) mergers in seeding and growing IMBHs in the \textsc{Dragon-II} simulation database, a suite of 19 direct $N$-body models representing dense clusters with up to $10^6$ stars. \textsc{Dragon-II} IMBHs have typical masses of $m_{\rm IMBH} = (100-380)$ M$_\odot$ and relatively large spins $χ_{\rm IMBH} > 0.6$. We find a link between the IMBH formation mechanism and the cluster structure. In clusters denser than $3\times 10^5$ M$_\odot$ pc$^{-3}$, the collapse of massive star collision products represents the dominant IMBH formation process, leading to the formation of heavy IMBHs ($m_{\rm IMBH} > 200$ M$_\odot$), possibly slowly rotating, that form over times $<5$ Myr and grow further via stellar accretion and mergers in just $<30$ Myr. BBH mergers are the dominant IMBH formation channel in less dense clusters, for which we find that the looser the cluster, the longer the formation time ($10-300$ Myr) and the larger the IMBH mass, although remaining within $200$ M$_\odot$. Strong dynamical scatterings and relativistic recoil efficiently eject all IMBHs in \textsc{Dragon-II} clusters, suggesting that IMBHs in this type of cluster are unlikely to grow beyond a few $10^2$ M$_\odot$. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 15 pages, 6 figures, 2 tables, 1 appendix. Comments welcome. Submitted to MNRAS

The Dragon-II simulations -- I. Evolution of single and binary compact objects in star clusters with up to 1 million stars

Authors: Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Mirek Giersz, Peter Berczik, Sara Rastello, Giuliano Iorio, Michela Mapelli, Massimiliano Gatto, Eva K. Grebel

Abstract: We present the first results of the \textsc{Dragon-II} simulations, a suite of 19 $N$-body simulations of star clusters with up to $10^6$ stars, with up to $33\%$ of them initially paired in binaries. In this work, we describe the main evolution of the clusters and their compact objects (COs). All \textsc{Dragon-II} clusters form in their centre a black hole (BH) subsystem with a density $10-100$… ▽ More We present the first results of the \textsc{Dragon-II} simulations, a suite of 19 $N$-body simulations of star clusters with up to $10^6$ stars, with up to $33\%$ of them initially paired in binaries. In this work, we describe the main evolution of the clusters and their compact objects (COs). All \textsc{Dragon-II} clusters form in their centre a black hole (BH) subsystem with a density $10-100$ times larger than the stellar density, with the cluster core containing $50-80\%$ of the whole BH population. In all models, the BH average mass steeply decreases as a consequence of BH burning, reaching values $\langle m_{\rm BH}\rangle < 15$ M$_\odot$ within $10-30$ relaxation times. Generally, our clusters retain only BHs lighter than $30$ M$_\odot$ over $30$ relaxation times. Looser clusters retain a higher binary fraction, because in such environments binaries are less likely disrupted by dynamical encounters. We find that BH-main sequence star binaries have properties similar to recently observed systems. Double CO binaries (DCOBs) ejected from the cluster exhibit larger mass ratios and heavier primary masses than ejected binaries hosting a single CO (SCOBs). Ejected SCOBs have BH masses $m_{\rm BH} = 3-20$ M$_\odot$, definitely lower than those in DCOBs ($m_{\rm BH} = 10-100$ M$_\odot$). △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 22 pages, 21 figures, 4 tables. Comments welcome. Submitted to MNRAS

MOCCA-Survey Database: Extra Galactic Globular Clusters. III. The population of black holes in Milky Way and Andromeda-like galaxies

Authors: A. Leveque, M. Giersz, A. Askar, M. Arca-Sedda, A. Olejak

Abstract: In this work, we investigate the black hole (BH) population of globular clusters (GCs) in Milky Way- (MW) and Andromeda- (M31) like galaxies. We combine the population synthesis code MASinGa and the MOCCA-Survey Database I to infer the properties of GCs harbouring a BH subsystem (BHS), an IMBH, or neither of those. We find that the typical number of GCs with a BHS, an IMBH, or none become comparab… ▽ More In this work, we investigate the black hole (BH) population of globular clusters (GCs) in Milky Way- (MW) and Andromeda- (M31) like galaxies. We combine the population synthesis code MASinGa and the MOCCA-Survey Database I to infer the properties of GCs harbouring a BH subsystem (BHS), an IMBH, or neither of those. We find that the typical number of GCs with a BHS, an IMBH, or none become comparable in the galactic outskirts, whilst the inner galactic regions are dominated by GCs without a significant dark component. Our models suggest that GCs harbouring a BHS are slightly heavier and with larger half-mass radii compared to the overall population. We retrieve the properties of binary BHs (BBHs) that have either merged in the last 3 Gyr or survived in their parent cluster until present-day. We find that around 80\% of the merging BBHs form due to dynamical interactions while the remaining originate from evolution of primordial binaries. We infer a merger rate for BBHs in the local Universe of $1.0-23\,\,\rm{yr^{-1}\,Gpc^{-3}}$, depending on the adopted assumptions. We find around 100-240 BBHs survive until present-day and are mostly concentrated in the inner few kpc of the galaxy. We estimate also the number of BHs transported into the galactic nucleus by infalling star clusters, finding around 1,000-3,000 BHs and 100-200 BBHs are transported over a time span of 12 Gyr. This enables us to constrains the total amount of BHs and BBHs binaries lurking in nuclear star cluster, i.e. $N_{BHs}=(1.4-2.2)\times10^4$ and $N_{BBHs}=700-1,100$. △ Less

Submitted 20 January, 2023; v1 submitted 4 September, 2022; originally announced September 2022.

Comments: 18 pages,15 figures

Forming short period sub-stellar companions in 47 Tucanae -- II. Analytic expressions for the orbital evolution of planets in dense environments

Authors: Andrew J. Winter, Cathie J. Clarke, Giovanni Rosotti, Mirek Giersz

Abstract: Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularis… ▽ More Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularising tides in the planet. We present an analytic framework for quantifying the role of dynamical encounters in high density environments during HEM. We show that encounters can enhance or suppress HEM, depending on the local stellar density and the initial semi-major axis $a_0$. For moderate densities, external perturbations can excite large eccentricities that allow a planet to circularise over the stellar lifetime. At extremely high densities, these perturbations can instead result in tidal disruption of the planet, thus yielding no HJ. This may explain the apparent excess of HJs in M67 compared with their local field star abundance versus their apparent deficit in 47 Tuc. Applying our analytic framework, we demonstrate that for an initial massive planet population similar to the field, the expected HJ occurrence rate in 47 Tuc is $f_\mathrm{HJ}=2.2\times 10^{-3}$, which remains consistent with present constraints. Future large (sample sizes $\gtrsim 10^5$) or sensitive transit surveys of stars in globular clusters are required to refute the hypothesis that the initial planet population is similar to the solar neighbourhood average. Non-detection in such surveys would have broad consequences for planet formation theory, implying planet formation rates in globular clusters must be suppressed across a wide range of $a_0$. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS - 28 pages, 17 figures

MOCCA-Survey Database: Extra Galactic Globular Clusters. II. Milky Way and Andromeda

Authors: A. Leveque, M. Giersz, M. Arca-Sedda, A. Askar

Abstract: A comprehensive study of the co-evolution of globular cluster systems (GCS) in galaxies requires the ability to model both the large scale dynamics (0.01 - 10 kpc) regulating their orbital evolution, and the small scale dynamics (sub-pc - AU) regulating the internal dynamics of each globular cluster (GC). In this work we present a novel method that combine semi-analytic models of GCS with fully se… ▽ More A comprehensive study of the co-evolution of globular cluster systems (GCS) in galaxies requires the ability to model both the large scale dynamics (0.01 - 10 kpc) regulating their orbital evolution, and the small scale dynamics (sub-pc - AU) regulating the internal dynamics of each globular cluster (GC). In this work we present a novel method that combine semi-analytic models of GCS with fully self-consistent Monte Carlo models to simultaneously evolve large GCSs. We use the population synthesis code MASinGa and the MOCCA-Survey Database I to create synthetic GC populations aimed at representing the observed features of GCs in the Milky Way (MW) and Andromeda (M31). Our procedure enables us to recover the spatial and mass distribution of GCs in such galaxies, and to constrain the amount of mass that GCs left either in the halo as dispersed debris, or in the galactic centre, where they can contribute to the formation of a nuclear star cluster (NSC) and can bring stellar and possibly intermediate mass black holes there. The final masses reported by our simulations are of a few order of magnitudes smaller than the observed values. These differences show that mass build-up of a NSC and central BHs in galaxies like MW and M31 cannot be solely explained by the infalling GC scenario. This build-up is likely to depend on the interplay between interactions and mergers of infalling GCs and gas. The latter can contribute to both in-situ star formation in the NSC and growth of the central BH. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 17 pages, 10 figures

A Monte Carlo study of early gas expulsion and evolution of star clusters: new simulations with the MOCCA code in the AMUSE framework

Authors: A. Leveque, M. Giersz, S. Banerjee, E. Vesperini, J. Hong, S. Portegies Zwart

Abstract: We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in… ▽ More We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 13 pages, 6 figures

MOCCA: Dynamics and evolution of binary stars of multiple stellar populations in tidally filling and underfilling globular star clusters

Authors: Arkadiusz Hypki, Mirek Giersz, Jongsuk Hong, Agostino Leveque, Abbas Askar, Diogo Belloni, Magdalena Otulakowska-Hypka

Abstract: We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we f… ▽ More We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we focus on the number ratio of FG and SG binaries, its spatial variation, and the way their abundances are affected by various cluster initial properties. We find that present-day SG stars are more abundant in clusters that were initially tidally filling. Conversely, FG stars stay more abundant in clusters that were initially tidally underfilling. We find that the ratio between binary fractions is not affected by the way we calculate these fractions (e.g. only main-sequence binaries (MS) or observational binaries, i.e. MS stars $> 0.4 M_{\odot}$ mass ratios $> 0.5$). This implies that the MS stars themselves are a very good proxy for probing entire populations of FG and SG. We also discuss how it relates to the observations of Milky Way GCs. We show that \MOCCA models are able to reproduce the observed range of SG fractions for Milky Way GCs for which we know these fractions. We show how the SG fractions depend on the initial conditions and provide some constraints for the initial conditions to have more numerous FG or SG stars at the Hubble time. △ Less

Submitted 3 October, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Acceped to MNRAS

The impact of stellar evolution on rotating star clusters: the gravothermal-gravogyro catastrophe and the formation of a bar of black holes

Authors: A. W. H. Kamlah, R. Spurzem, P. Berczik, M. Arca Sedda, F. Flammini Dotti, N. Neumayer, X. Pang, Q. Shu, A. Tanikawa, M. Giersz

Abstract: We present results from a suite of eight direct N-body simulations, performed with \textsc{Nbody6++GPU}, representing realistic models of rotating star clusters with up to $1.1\times 10^5$ stars. Our models feature primordial (hard) binaries, a continuous mass spectrum, differential rotation, and tidal mass loss induced by the overall gravitational field of the host galaxy. We explore the impact o… ▽ More We present results from a suite of eight direct N-body simulations, performed with \textsc{Nbody6++GPU}, representing realistic models of rotating star clusters with up to $1.1\times 10^5$ stars. Our models feature primordial (hard) binaries, a continuous mass spectrum, differential rotation, and tidal mass loss induced by the overall gravitational field of the host galaxy. We explore the impact of rotation and stellar evolution on the star cluster dynamics. In all runs for rotating star clusters we detect a previously predicted mechanism: an initial phase of violent relaxation followed by the so-called gravogyro catastrophe. We find that the gravogyro catastrophe reaches a finite amplitude, which depends in strength on the level of the bulk rotation, and then levels off. After this phase the angular momentum is transferred from high-mass to low-mass particles in the cluster (both stars and compact objects). Simultaneously, the system becomes gravothermally unstable and collapses, thus undergoing the so-called gravothermal-gravogyro catastrophe. Comparing models with and without stellar evolution, we find an interesting difference. When stellar evolution is not taken into account, the whole process proceeds at a faster pace. The population of heavy objects tend to form a triaxial structure that rotates in the cluster centre. When stellar evolution is taken into account, we find that such a {\it rotating bar} is populated by stellar black holes and their progenitors. The triaxial structure becomes axisymmetric over time, but we also find that the models without stellar evolution suffer repeated gravogyro catastrophes as sufficient angular momentum and mass are removed by the tidal field. △ Less

Submitted 31 August, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 19 pages, 9 figures, 2 tables, Accepted by MNRAS

Influence of tidal dissipation on outcomes of binary-single encounters between stars and black holes in stellar clusters

Authors: Lucas Hellström, Abbas Askar, Alessandro A. Trani, Mirek Giersz, Ross P. Church, Johan Samsing

Abstract: In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each… ▽ More In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission (using post-Newtonian (PN) corrections) and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH-star mergers by up to 75 per cent but it does not affect the number of BH-BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH-BH mergers during binary-single encounters. However, BH-star mergers are largely unaffected by the inclusion of these terms. △ Less

Submitted 29 September, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 14 pages, 6 figures, 12 tables, accepted for publication in MNRAS

MOCCA-SURVEY Database II -- Properties of Intermediate Mass Black Holes escaping from star clusters

Authors: Konrad Maliszewski, Mirek Giersz, Dorota Gondek-Rosińska, Abbas Askar, Arkadiusz Hypki

Abstract: In this work we investigate properties of intermediate-mass black holes (IMBHs) that escape from star clusters due to dynamical interactions. The studied models were simulated as part of the preliminary second survey carried out using the MOCCA code (MOCCA-SURVEY Database II), which is based on the Monte Carlo N-body method and does not include gravitational wave recoil kick prescriptions of the b… ▽ More In this work we investigate properties of intermediate-mass black holes (IMBHs) that escape from star clusters due to dynamical interactions. The studied models were simulated as part of the preliminary second survey carried out using the MOCCA code (MOCCA-SURVEY Database II), which is based on the Monte Carlo N-body method and does not include gravitational wave recoil kick prescriptions of the binary black hole merger product. We have found that IMBHs are more likely to be formed and ejected in models where both initial central density and central escape velocities have high values. Most of our studied objects escape in a binary with another black hole (BH) as their companion and have masses between $100$ and $140\: M_{\odot}$. Escaping IMBHs tend to build-up mass most effectively through repeated mergers in a binary with BHs due to gravitational wave emission. Binaries play a key role in their ejection from the system as they allow these massive objects to gather energy needed for escape. The binaries in which IMBHs escape tend to have very high binding energy at the time of escape and the last interaction is strong but does not involve a massive intruder. These IMBHs gain energy needed to escape the cluster gradually in successive dynamical interactions. We present specific examples of the history of IMBH formation and escape from star cluster models. We also discuss the observational implications of our findings as well as the potential influence of the gravitational wave recoil kicks on the process. △ Less

Submitted 18 June, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 11 pages, 11 figures, accepted for publish to MNRAS

Forming short period sub-stellar companions in 47 Tuc: I. Dynamical model and brown dwarf tidal capture rates

Authors: Andrew J. Winter, Giovanni P. Rosotti, Cathie Clarke, Mirek Giersz

Abstract: Stars in globular clusters formed and evolved in the most extreme environment: high density and low metallicity. If the formation of stars and planets are at all sensitive to environmental conditions, this should therefore be evident in globular clusters. Observations have indicated that hot Jupiters are at least an order of magnitude less prevalent in the central region of the globular cluster 47… ▽ More Stars in globular clusters formed and evolved in the most extreme environment: high density and low metallicity. If the formation of stars and planets are at all sensitive to environmental conditions, this should therefore be evident in globular clusters. Observations have indicated that hot Jupiters are at least an order of magnitude less prevalent in the central region of the globular cluster 47 Tucanae than in the field. In this work, we explore the claims in the literature for additional consequences for the low mass stellar initial mass function. Tidal capture, the mechanism that produces X-ray binaries in globular clusters, applies also to brown dwarfs (BDs). This process produces tight stellar-BD binaries that would be detectable by transit surveys. Applying a Monte Carlo dynamical evolution model, we compute the overall BD capture rates. We find that the number of captures is lower than previous estimates. Capture efficiency increases steeply with stellar mass, which means that mass segregation reduces capture efficiency as BDs and low mass stars occupy the same regions. The result of this effect is that the current constraints on the short period companion fraction remains marginally consistent with initially equal numbers of BDs and stars. However, our findings suggest that expanding the sample in 47 Tuc or surveying other globular clusters for close sub-stellar companions can yield constraints on the sub-stellar initial mass function in these environments. We estimate the capture rates in other globular clusters and suggest that 47 Tuc remains a promising target for future transit surveys. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 15 pages, 13 figures, accepted for publication in MNRAS

Black hole mergers in compact star clusters and massive black hole formation beyond the mass-gap

Authors: Francesco Paolo Rizzuto, Thorsten Naab, Rainer Spurzem, Manuel Arca-Sedda, Mirek Giersz, Jeremiah Paul Ostriker, Sambaran Banerjee

Abstract: We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears… ▽ More We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears of evolution, each cluster hosts several black hole (BH) merger events which nearly cover the complete mass range of primary and secondary BH masses for current LIGO/Virgo/Kagra gravitational wave detections. The importance of gravitational recoil is estimated statistically. We present several possible formation paths of massive BHs above the assumed lower PISNe mass-gap limit ($45 M_\odot$) into the intermediate-mass BH (IMBH) regime ($> 100 M_\odot$) which include collisions of stars and BHs as well as the direct collapse of stellar merger remnants with low mass cores. The stellar evolution updates result in the early formation of higher mass stellar BHs than for the previous model. The resulting higher collision rates with massive stars support the rapid formation of massive BHs. For models assuming a high accretion efficiency for star-BH mergers, we present a first-generation formation scenario for GW190521-like events, a merger of two BHs in the PISN mass-gap, which is dominated by star-BH mergers. This IMBH formation path is independent of gravitational recoil and therefore conceivable in dense stellar systems with low escape velocities. One simulated cluster even forms an IMBH binary (153$M_\odot$,173$M_\odot$) which is expected to merge within a Hubble time. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Preparing the next gravitational million-body simulations: Evolution of single and binary stars in Nbody6++GPU, MOCCA and McLuster

Authors: A. W. H. Kamlah, A. Leveque, R. Spurzem, M. Arca Sedda, A. Askar, S. Banerjee, P. Berczik, M. Giersz, J. Hurley, D. Belloni, L. Kühmichel, L. Wang

Abstract: We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of… ▽ More We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of 10~Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf natal kicks enabled/disabled. Compared to \texttt{Nbody6++GPU}, the \texttt{MOCCA} models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The \texttt{MOCCA} models produce more black holes (BHs) and helium white dwarfs (WDs), whilst \texttt{Nbody6++GPU} models are characterised by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorise the stellar evolution recipes available in \texttt{Nbody6++GPU}, \texttt{MOCCA} and \texttt{Mcluster} into four levels: the one implemented in previous \texttt{Nbody6++GPU} and \texttt{MOCCA} versions (\texttt{level A}), state-of-the-art prescriptions (\texttt{level B}), some in a testing phase (\texttt{level C}), and those that will be added in future versions of our codes. △ Less

Submitted 21 December, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: 31 pages, 15 figures, 5 tables, Accepted by MNRAS

Breaching the limit: formation of GW190521-like and IMBH mergers in young massive clusters

Authors: Manuel Arca Sedda, Francesco Paolo Rizzuto, Thorsten Naab, Jeremiah Ostriker, Mirek Giersz, Rainer Spurzem

Abstract: The LIGO-Virgo-Kagra collaboration (LVC) discovered recently GW190521, a gravitational wave (GW) source associated with the merger between two black holes (BHs) with mass $66$ M$_\odot$ and $>85$ M$_\odot$. GW190521 represents the first BH binary (BBH) merger with a primary mass falling in the "upper mass-gap" and the first leaving behind a $\sim 150$ M$_\odot$ remnant. So far, the LVC reported th… ▽ More The LIGO-Virgo-Kagra collaboration (LVC) discovered recently GW190521, a gravitational wave (GW) source associated with the merger between two black holes (BHs) with mass $66$ M$_\odot$ and $>85$ M$_\odot$. GW190521 represents the first BH binary (BBH) merger with a primary mass falling in the "upper mass-gap" and the first leaving behind a $\sim 150$ M$_\odot$ remnant. So far, the LVC reported the discovery of four further mergers having a total mass $>100$ M$_\odot$, i.e. in the intermediate-mass black holes (IMBH) mass range. Here, we discuss results from a series of 80 $N$-body simulations of young massive clusters (YMCs) that implement relativistic corrections to follow compact object mergers. We discover the development of a GW190521-like system as the result of a 3rd-generation merger, and four IMBH-BH mergers with total mass $~(300-350)$ M$_\odot$. We show that these IMBH-BH mergers are low-frequency GW sources detectable with LISA and DECIGO out to redshift $z=0.01-0.1$ and $z>100$, and we discuss how their detection could help unravelling IMBH natal spins. For the GW190521 test case, we show that the 3rd-generation merger remnant has a spin and effective spin parameter that matches the $90\%$ credible interval measured for GW190521 better than a simpler double merger and comparably to a single merger. Due to GW recoil kicks, we show that retaining the products of these mergers require birth-sites with escape velocities $\gtrsim 50-100$ km s$^{-1}$, values typically attained in galactic nuclei and massive clusters with steep density profiles. △ Less

Submitted 3 September, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 23 pages, 11 figures, 2 tables, 2 Appendix. ApJ accepted

MOCCA-SURVEY Database -- I. Tidal disruption events of white dwarfs in globular clusters and young massive clusters

Authors: Ataru Tanikawa, Mirek Giersz, Manuel Arca Sedda

Abstract: We exploit more than 1000 star cluster Monte Carlo models from the MOCCA-SURVEY Database I, to infer the local rate density of white dwarf (WD) tidal disruption events (TDEs) in globular clusters (GCs) and young massive clusters (YMCs). We suggest that the WD TDE rate for GCs and YMCs in the local universe is $\sim 90$-$500~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, with 90 % of WD TDEs occurring in GCs. The… ▽ More We exploit more than 1000 star cluster Monte Carlo models from the MOCCA-SURVEY Database I, to infer the local rate density of white dwarf (WD) tidal disruption events (TDEs) in globular clusters (GCs) and young massive clusters (YMCs). We suggest that the WD TDE rate for GCs and YMCs in the local universe is $\sim 90$-$500~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, with 90 % of WD TDEs occurring in GCs. The total WD TDE rate density is $\sim 9$-$50$ times larger than estimated before. Our results show that thermonuclear explosions induced by WD TDEs can be observed at a rate of $\sim 100$-$550~{\rm yr}^{-1}$ by the next generation optical surveys, such as the Legacy Survey of Space & Time by the Vera C. Rubin Observatory. We also find that massive WDs are preferentially disrupted due to mass segregation, and that 20 % of exploding WDs have $\gtrsim 1.0 M_\odot$ despite of the small population of such WDs. △ Less

Submitted 14 July, 2022; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 19 pages, 16 figures, MNRAS accepted

Dynamical Evolution of Multiple-Population Globular Clusters

Authors: E. Vesperini, J. Hong, M. Giersz, A. Hypki

Abstract: We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolu… ▽ More We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolution of the fraction of 2G stars and the link between the evolution of the fraction of 2G stars and various dynamical parameters; 2) the link between the fraction of 2G stars inside the cluster and in the population of escaping stars during a cluster's dynamical evolution; 3) the dynamics of the spatial mixing of the first-generation (1G) and 2G stars and the details of the structural properties of the two populations as they evolve toward mixing; 4) the implications of the initial differences between the spatial distribution of 1G and 2G stars for the evolution of the anisotropy in the velocity distribution and the expected radial profile of the 1G and 2G anisotropy for clusters at different stages of their dynamical history; 5) the variation of the degree of energy equipartition of the 1G and the 2G populations as a function of the distance from the cluster's centre and the cluster's evolutionary phase. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 16 pages, 19 figures. Accepted for publication in MNRAS

MOCCA SURVEY Database I: Binary Black Hole Mergers from Globular Clusters with Intermediate Mass Black Holes

Authors: Jongsuk Hong, Abbas Askar, Mirek Giersz, Arkadiusz Hypki, Suk-Jin Yoon

Abstract: The dynamical formation of black hole binaries in globular clusters that merge due to gravitational waves occurs more frequently in higher stellar density. Meanwhile, the probability to form intermediate mass black holes (IMBHs) also increases with the density. To explore the impact of the formation and growth of IMBHs on the population of stellar mass black hole binaries from globular clusters, w… ▽ More The dynamical formation of black hole binaries in globular clusters that merge due to gravitational waves occurs more frequently in higher stellar density. Meanwhile, the probability to form intermediate mass black holes (IMBHs) also increases with the density. To explore the impact of the formation and growth of IMBHs on the population of stellar mass black hole binaries from globular clusters, we analyze the existing large survey of Monte-Carlo globular cluster simulation data (MOCCA SURVEY Database I). We show that the number of binary black hole mergers agrees with the prediction based on clusters' initial properties when the IMBH mass is not massive enough or the IMBH seed forms at a later time. However, binary black hole formation and subsequent merger events are significantly reduced compared to the prediction when the present-day IMBH mass is more massive than $\sim10^4 \rm M_{\odot}$ or the present-day IMBH mass exceeds about 1 per cent of cluster's initial total mass. By examining the maximum black hole mass in the system at the moment of black hole binary escaping, we find that $\sim$ 90 per cent of the merging binary black holes escape before the formation and growth of the IMBH. Furthermore, large fraction of stellar mass black holes are merged into the IMBH or escape as single black holes from globular clusters in cases of massive IMBHs, which can lead to the significant under-population of binary black holes merging with gravitational waves by a factor of 2 depending on the clusters' initial distributions. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: 9 pages, 8 figures, Accepted for publication in MNRAS

Intermediate Mass Black Hole Formation in compact Young Massive Star Clusters

Authors: Francesco Paolo Rizzuto, Thorsten Naab, Rainer Spurzem, Mirek Giersz, J. P. Ostriker, N. C. Stone, Long Wang, Peter Berczik, M. Rampp

Abstract: Young dense massive star clusters are a promising environment for the formation of intermediate mass black holes (IMBHs) through collisions. We present a set of 80 simulations carried out with Nbody6++GPU of 10 initial conditions for compact $\sim 7 \times 10^4 M_{\odot}$ star clusters with half-mass radii $R_\mathrm{h} \lesssim 1 pc$, central densities… ▽ More Young dense massive star clusters are a promising environment for the formation of intermediate mass black holes (IMBHs) through collisions. We present a set of 80 simulations carried out with Nbody6++GPU of 10 initial conditions for compact $\sim 7 \times 10^4 M_{\odot}$ star clusters with half-mass radii $R_\mathrm{h} \lesssim 1 pc$, central densities $��_\mathrm{core} \gtrsim 10^5 M_\odot pc^{-3}$, and resolved stellar populations with 10\% primordial binaries. Very massive stars (VMSs) with masses up to $\sim 400 M_\odot$ grow rapidly by binary exchange and three-body scattering events with main sequences stars in hard binaries. Assuming that in VMS - stellar BH collisions all stellar material is accreted onto the BH, IMBHs with masses up to $M_\mathrm{BH} \sim 350 M_\odot$ can form on timescales of $\lesssim 15$ Myr. This process was qualitatively predicted from Monte Carlo MOCCA simulations. Despite the stochastic nature of the process - typically not more than 3/8 cluster realisations show IMBH formation - we find indications for higher formation efficiencies in more compact clusters. Assuming a lower accretion fraction of 0.5 for VMS - BH collisions, IMBHs can also form. The process might not work for accretion fractions as low as 0.1. After formation, the IMBHs can experience occasional mergers with stellar mass BHs in intermediate mass-ratio inspiral events (IMRIs) on a 100 Myr timescale. Realised with more than $10^5$ stars, 10 \% binaries, the assumed stellar evolution model with all relevant evolution processes included and 300 Myr simulation time, our large suite of simulations indicates that IMBHs of several hundred solar masses might form rapidly in massive star clusters right after their birth while they are still compact. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: Submitted

MOCCA Survey Database: Extra Galactic Globular Clusters. I. Method and first results

Authors: Agostino Leveque, Mirosław Giersz, Maurizio Paolillo

Abstract: Over the last few decades, exhaustive surveys of extra Galactic globular clusters (EGGCs) have become feasible. Only recently, kinematical information of globular clusters (GCs) were available through Gaia DR2 spectroscopy and also proper motions. On the other hand, simulations of GCs can provide detailed information about the dynamical evolution of the system. We present a preliminary study of EG… ▽ More Over the last few decades, exhaustive surveys of extra Galactic globular clusters (EGGCs) have become feasible. Only recently, kinematical information of globular clusters (GCs) were available through Gaia DR2 spectroscopy and also proper motions. On the other hand, simulations of GCs can provide detailed information about the dynamical evolution of the system. We present a preliminary study of EGGCs' properties for different dynamical evolutionary stages. We apply this study to 12 Gyr-old GCs simulated as part of the MOCCA Survey Database. Mimicking observational limits, we consider only a subssample of the models in the database, showing that it is possible to represent observed Galactic GCs. In order to distinguish between different dynamical states of EGGCs, at least three structural parameters are necessary. The best distinction is achieved by considering the central parameters, those being observational core radius, central surface brightness, ratio between central and half-mass velocity dispersion, or similarly considering the central color, the central V magnitude and the ratio between central and half-mass radius velocity dispersion, although such properties could be prohibitive with current technologies. A similar but less solid result is obtained considering the average properties at the half-light radius, perhaps accessible presently in the Local Group. Additionally, we mention that the color spread in EGGCs due to internal dynamical models, at fixed metallcity, could be just as important due to the spread in metallicity. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: 16 pages, 11 figures

On the absence of symbiotic stars in globular clusters

Authors: Diogo Belloni, Joanna Mikołajewska, Krystian Iłkiewicz, Matthias R. Schreiber, Mirek Giersz, Liliana E. Rivera Sandoval, Claudia V. Rodrigues

Abstract: Even though plenty of symbiotic stars (SySts) have been found in the Galactic field and nearby galaxies, not a single one has ever been confirmed in a Galactic globular cluster (GC). We investigate the lack of such systems in GCs for the first time by analysing 144 GC models evolved with the MOCCA code, which have different initial properties and are roughly representative of the Galactic GC popul… ▽ More Even though plenty of symbiotic stars (SySts) have been found in the Galactic field and nearby galaxies, not a single one has ever been confirmed in a Galactic globular cluster (GC). We investigate the lack of such systems in GCs for the first time by analysing 144 GC models evolved with the MOCCA code, which have different initial properties and are roughly representative of the Galactic GC population. We focus here on SySts formed through the wind-accretion channel, which can be consistently modelled in binary population synthesis codes. We found that the orbital periods of the majority of such SySts are sufficiently long (${\gtrsim10^3}$ d) so that, for very dense GC models, dynamical interactions play an important role in destroying their progenitors before the present day (${\sim11-12}$ Gyr). In less dense GC models, some SySts are still predicted to exist. However, these systems tend to be located far from the central parts (${\gtrsim70}$ per cent are far beyond the half-light radius) and are sufficiently rare (${\lesssim1}$ per GC per Myr), which makes their identification rather difficult in observational campaigns. We propose that future searches for SySts in GCs should be performed in the outskirts of nearby low-density GCs with sufficiently long half-mass relaxation times and relatively large Galactocentric distances. Finally, we obtained spectra of the candidate proposed in $ω$ Cen (SOPS IV e-94) and showed that this object is most likely not a SySt. △ Less

Submitted 30 June, 2020; v1 submitted 11 April, 2020; originally announced April 2020.

Comments: accepted by MNRAS, 12 pages

MOCCA survey database I: preliminary mock Extra Galactic Globular Cluster observations

Authors: Agostino Leveque, Mirosław Giersz

Abstract: The photometric properties that we could observe for Extra-Galactic Globular Clusters (EGGCs) are the integrated light of the system and for nearby EGGCs it also is possible to measure both half-light radii and the color spatial distribution, e.g. for areas smaller and larger than the half-light radius. No information about the internal dynamical state of the system could be directly obtained from… ▽ More The photometric properties that we could observe for Extra-Galactic Globular Clusters (EGGCs) are the integrated light of the system and for nearby EGGCs it also is possible to measure both half-light radii and the color spatial distribution, e.g. for areas smaller and larger than the half-light radius. No information about the internal dynamical state of the system could be directly obtained from observations. On the other hand, simulations of Globular Clusters (GCs) can provide detailed information about the dynamical evolution of the system. We present a preliminary study of EGGCs' photometric properties for different dynamical evolutionary stages. We apply this study to 12 Gyr old GCs simulated as part of the MOCCA Survey Database. We determine the magnitudes in different bands from their projected snapshots using the Flexible Stellar Population Synthesis (FSPS) code and we measure the half-light radii from the surface brightness. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 4 pages, 3 figures, Proceedings IAU Symposium No. 351, 2019

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

Authors: Mirek Giersz, Abbas Askar, Long Wang, Arkadiusz Hypki, Agostino Leveque, Rainer Spurzem

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 connec… ▽ More 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. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: 4 pages, 1 figure. Contribution to the proceedings of the IAU Symposium No. 351: "Star Clusters: from the Milky Way to the Early Universe"

Stellar-Mass Black Holes in Globular Clusters: Dynamical Consequences and Observational Signatures

Authors: Abbas Askar, Mirek Giersz, Manuel Arca Sedda, Ammar Askar, Mario Pasquato, Agostino Leveque

Abstract: 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… ▽ More 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. △ Less

Submitted 31 July, 2019; originally announced July 2019.

Comments: 4 pages, 2 figures, 1 table. Contribution to the proceedings of the IAU Symposium No. 351: "Star Clusters: from the Milky Way to the Early Universe"

Are most Cataclysmic Variables in Globular Clusters dynamically formed?

Authors: Diogo Belloni, Mirek Giersz, Liliana E. Rivera Sandoval, Abbas Askar, Pawel Ciecielag

Abstract: 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… ▽ More 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. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: 4 pages, 1 figures, contribution to the proceedings of the IAU Symposium No. 351, 2019, "Star Clusters: from the Milky Way to the Early Universe"

MOCCA-SURVEY Database I. Intermediate mass black holes in Milky Way globular clusters and their connection to supermassive black holes

Authors: Manuel Arca Sedda, Abbas Askar, Mirek Giersz

Abstract: 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 ident… ▽ More 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. △ Less

Submitted 2 May, 2019; originally announced May 2019.

Comments: 19 pages, 1 appendix, 3 tables, 16 figures. MNRAS submitted. Comments welcome

MOCCA Survey Database I. BHs in star clusters

Authors: Mirek Giersz, Abbas Askar, Jakub Klencki, Jakub Morawski

Abstract: 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 subs… ▽ More 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. △ Less

Submitted 7 April, 2019; originally announced April 2019.

Comments: 10 pages, 4 figures. MG15 Proceedings

MOCCA Survey Database I: Dissolution of tidally filling star clusters harbouring BH subsystems

Authors: Mirek Giersz, Abbas Askar, Long Wang, Arkadiusz Hypki, Agostino Leveque, Rainer Spurzem

Abstract: 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 clust… ▽ More 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. △ Less

Submitted 2 August, 2019; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: 14 pages, 14 figures, accepted to MNRAS

In Search of the Thermal Eccentricity Distribution

Authors: Aaron M. Geller, Nathan W. C. Leigh, Mirek Giersz, Kyle Kremer, Frederic A. Rasio

Abstract: About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary… ▽ More About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two. △ Less

Submitted 31 January, 2019; originally announced February 2019.

Comments: 14 pages, 7 figures, accepted for publication in The Astrophysical Journal

Finding Black Holes with Black Boxes -- Using Machine Learning to Identify Globular Clusters with Black Hole Subsystems

Authors: Ammar Askar, Abbas Askar, Mario Pasquato, Mirek Giersz

Abstract: 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… ▽ More 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. △ Less

Submitted 1 March, 2019; v1 submitted 15 November, 2018; originally announced November 2018.

Comments: 20 pages, 9 figures, 7 tables. Accepted for publication in MNRAS. Source code available at https://github.com/ammaraskar/black-holes-black-boxes

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- IV. cataclysmic variables -- properties of bright and faint populations

Authors: Diogo Belloni, Mirek Giersz, Liliana E. Rivera Sandoval, Abbas Askar, Paweł Ciecieląg

Abstract: 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… ▽ More 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). △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 18 pages, 13 figures; accepted for publication in MNRAS

Binary Black Hole Mergers from Globular Clusters: the Impact of Globular Cluster Properties

Authors: Jongsuk Hong, Enrico Vesperini, Abbas Askar, Mirek Giersz, Magdalena Szkudlarek, Tomasz Bulik

Abstract: 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 form… ▽ More 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. △ Less

Submitted 13 August, 2018; originally announced August 2018.

Comments: 12 pages, 12 figures, 2 tables, accepted for publication in MNRAS

No cataclysmic variables missing: higher merger rate brings into agreement observed and predicted space densities

Authors: Diogo Belloni, Matthias R. Schreiber, Mónica Zorotovic, Krystian Iłkiewicz, Jarrod R. Hurley, Mirek Giersz, Felipe Lagos

Abstract: The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $ρ\lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated vers… ▽ More The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $ρ\lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated version of the Binary Stellar Evolution (BSE) code for single and binary star evolution. We find that the recently suggested empirical prescription of consequential angular momentum loss (CAML) brings into agreement predicted and observed space densities of CVs and period bouncers. To progress with our understanding of CV evolution it is crucial to understand the physical mechanism behind empirical CAML. Our changes to the BSE code are also provided in details, which will allow the community to accurately model mass transfer in interacting binaries in which degenerate objects accrete from low-mass main-sequence donor stars. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: 11 pages, 4 figures, 1 table; accepted for publication in MNRAS

Black Hole Mergers from Globular Clusters Observable by LISA and LIGO: Results from post-Newtonian Binary-Single Scatterings

Authors: Johan Samsing, Daniel J. D'Orazio, Abbas Askar, Mirek Giersz

Abstract: 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 Dat… ▽ More 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. △ Less

Submitted 23 February, 2018; originally announced February 2018.

Comments: 7 pages, 2 figures. comments welcome

MOCCA-SURVEY Database I: Galactic Globular Clusters Harbouring a Black Hole Subsystem

Authors: Abbas Askar, Manuel Arca Sedda, Mirek Giersz

Abstract: 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 c… ▽ More 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. △ Less

Submitted 3 May, 2018; v1 submitted 14 February, 2018; originally announced February 2018.

Comments: 13 pages, 5 figures, 3 tables, accepted for publication in MNRAS. Revised manuscript includes minor corrections and a new subsection (2.5)

MOCCA-SURVEY Database I: Assessing GW kick retention fractions for BH-BH mergers in globular clusters

Authors: J. Morawski, M. Giersz, A. Askar, K. Belczynski

Abstract: 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 m… ▽ More 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. △ Less

Submitted 11 September, 2018; v1 submitted 4 February, 2018; originally announced February 2018.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

MOCCA-SURVEY Database I. Unravelling black hole subsystems in globular clusters

Authors: Manuel Arca Sedda, Abbas Askar, Mirek Giersz

Abstract: 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… ▽ More 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, $ρ_\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. △ Less

Submitted 13 April, 2018; v1 submitted 2 January, 2018; originally announced January 2018.

Comments: 14 pages, 12 figures, resubmitted to MNRAS after requested revision

MOCCA-SURVEY Database I: Eccentric Black Hole Mergers During Binary-Single Interactions In Globular Clusters

Authors: Johan Samsing, Abbas Askar, Mirek Giersz

Abstract: 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… ▽ More 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. △ Less

Submitted 27 January, 2018; v1 submitted 17 December, 2017; originally announced December 2017.

Comments: 5 pages, 1 figure, revised manuscript has been accepted for publication in ApJ

The origin of the first neutron star -- neutron star merger

Authors: K. Belczynski, A. Askar, M. Arca-Sedda, M. Chruslinska, M. Donnari, M. Giersz, M. Benacquista, R. Spurzem, D. Jin, G. Wiktorowicz, D. Belloni

Abstract: 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 o… ▽ More 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. △ Less

Submitted 28 March, 2018; v1 submitted 2 December, 2017; originally announced December 2017.

Comments: A&A: accepted

Dynamical equivalence, the origin of the Galactic field stellar and binary population, and the initial radius--mass relation of embedded clusters

Authors: Diogo Belloni, Pavel Kroupa, Helio J. Rocha-Pinto, Mirek Giersz

Abstract: In order to allow a better understanding of the origin of Galactic field populations, dynamical equivalence of stellar-dynamical systems has been postulated by Kroupa and Belloni et al. to allow mapping of solutions of the initial conditions of embedded clusters such that they yield, after a period of dynamical processing, the Galactic field population. Dynamically equivalent systems are defined t… ▽ More In order to allow a better understanding of the origin of Galactic field populations, dynamical equivalence of stellar-dynamical systems has been postulated by Kroupa and Belloni et al. to allow mapping of solutions of the initial conditions of embedded clusters such that they yield, after a period of dynamical processing, the Galactic field population. Dynamically equivalent systems are defined to initially and finally have the same distribution functions of periods, mass ratios and eccentricities of binary stars. Here we search for dynamically equivalent clusters using the {\sc mocca} code. The simulations confirm that dynamically equivalent solutions indeed exist. The result is that the solution space is next to identical to the radius--mass relation of Marks \& Kroupa, $\left( r_h/{\rm pc} \right)= 0.1^{+0.07}_{-0.04}\, \left( M_{\rm ecl}/{\rm M}_\odot \right)^{0.13\pm0.04}$. This relation is in good agreement with the observed density of molecular cloud clumps. According to the solutions, the time-scale to reach dynamical equivalence is about 0.5~Myr which is, interestingly, consistent with the lifetime of ultra-compact HII regions and the time-scale needed for gas expulsion to be active in observed very young clusters as based on their dynamical modelling. △ Less

Submitted 2 December, 2017; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 6 pages, 2 figures; accepted for publication in MNRAS, minor changes made to the text after proofs

The dynamical origin of multiple populations in intermediate-age clusters in the Magellanic Clouds

Authors: Jongsuk Hong, Richard de Grijs, Abbas Askar, Peter Berczik, Chengyuan Li, Long Wang, Licai Deng, M. B. N. Kouwenhoven, Mirek Giersz, Rainer Spurzem

Abstract: 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… ▽ More 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. △ Less

Submitted 28 July, 2017; originally announced July 2017.

Comments: 11 pages, 8 figures, 1 table; accepted for publication in MNRAS

On the initial binary population for star cluster simulations

Authors: Diogo Belloni, Abbas Askar, Mirek Giersz, Pavel Kroupa, Helio J. Rocha-Pinto

Abstract: 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$… ▽ More 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. △ Less

Submitted 13 July, 2017; originally announced July 2017.

Comments: 18 pages, 10 figures, 4 tables; accepted for publication in MNRAS

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

Authors: K. Belczynski, J. Klencki, C. E. Fields, A. Olejak, E. Berti, G. Meynet, C. L. Fryer, D. E. Holz, R. O'Shaughnessy, D. A. Brown, T. Bulik, S. C. Leung, K. Nomoto, P. Madau, R. Hirschi, E. Kaiser, S. Jones, S. Mondal, M. Chruslinska, P. Drozda, D. Gerosa, Z. Doctor, M. Giersz, S. Ekstrom, C. Georgy , et al. (9 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 6 March, 2020; v1 submitted 21 June, 2017; originally announced June 2017.

Comments: A&A accepted: revised comparison with LIGO/Virgo (41 pages)

Journal ref: A&A 636, A104 (2020)

Mass Evaporation Rate of Globular Clusters in a Strong Tidal Field

Authors: Juan P. Madrid, Nathan W. C. Leigh, Jarrod R. Hurley, Mirek Giersz

Abstract: The mass evaporation rate of globular clusters evolving in a strong Galactic tidal field is derived through the analysis of large, multi-mass $N$-body simulations. For comparison, we also study the same evaporation rates using MOCCA Monte Carlo models for globular cluster evolution. Our results show that the mass evaporation rate is a dynamical value; that is, far from a constant single number fou… ▽ More The mass evaporation rate of globular clusters evolving in a strong Galactic tidal field is derived through the analysis of large, multi-mass $N$-body simulations. For comparison, we also study the same evaporation rates using MOCCA Monte Carlo models for globular cluster evolution. Our results show that the mass evaporation rate is a dynamical value; that is, far from a constant single number found in earlier analytical work and commonly used in the literature. Moreover, the evaporation rate derived with these simulations is higher than values previously published. These models also show that the value of the mass evaporation rate depends on the strength of the tidal field. We give an analytical estimate of the mass evaporation rate as a function of time and galactocentric distance xi(R_{GC},t). Upon extrapolating this formula to smaller R_{GC} values, our results provide tentative evidence for a very high xi value at small R_{GC}. Our results suggest that the corresponding mass loss in the inner Galactic potential could be high and it should be accounted for when star clusters pass within it. This has direct relevance to nuclear cluster formation/growth via the infall of globular clusters through dynamical friction. As an illustrative example, we estimate how the evaporation rate increases for a ~10^5 solar masses globular cluster that decays through dynamical friction into the galactic centre. We discuss the findings of this work in relation to the formation of nuclear star clusters by inspiraling globular clusters. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: MNRAS in press

COCOA Code for Creating Mock Observations of Star Cluster Models

Authors: Abbas Askar, Mirek Giersz, Wojciech Pych, Emanuele Dalessandro

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 9 January, 2018; v1 submitted 27 March, 2017; originally announced March 2017.

Comments: 18 pages, 12 figures, accepted for publication in MNRAS. Revised manuscript has a new title, better quality figures and many other improvements. COCOA can be downloaded from: https://github.com/abs2k12/COCOA (comments are welcome)

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- III. Cataclysmic variables -- Implications of model assumptions

Authors: Diogo Belloni, Mónica Zorotovic, Matthias R. Schreiber, Nathan W. C. Leigh, Mirek Giersz, Abbas Askar

Abstract: 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 los… ▽ More 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. △ Less

Submitted 6 March, 2017; originally announced March 2017.

Comments: 18 pages, 7 figures, 2 tables; accepted for publication in MNRAS