Showing 1–50 of 132 results for author: Spurzem, R

Efficiency of black hole formation via collisions in stellar systems: An analysis of data from simulations and observations

Authors: M. C. Vergara, D. R. G. Schleicher, A. Escala, B. Reinoso, F. Flammini Dotti, A. W. H. Kamlah, M. Liempi, N. Hoyer, N. Neumayer, R. Spurzem

Abstract: This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar syst… ▽ More This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar systems leading to the formation of very massive stars that subsequently collapse into black holes. We analyze data from numerical simulations and observations of diverse stellar systems, encompassing various initial conditions, initial mass functions, and evolution scenarios. We computed a critical mass, determined by the interplay of collision time, system age, and initial properties of the star cluster. The efficiency of black hole formation ($ε_{\mathrm{BH}}$) is defined as the ratio of initial stellar mass divided by critical mass. We find that stellar systems with a ratio of initial stellar mass over critical mass above 1 exhibit a high efficiencies in terms of black hole formation, ranging from $30-100\%$. While there is some scatter, potentially attributed to complex system histories and the presence of gas, the results highlight the potential for achieving high efficiencies via a purely collisional channel in black hole formation. In conclusion, this theoretical exploration elucidates the connection between star clusters and black hole formation. The study underscores the significance of UCDs, GCs, and NSCs as environments conducive to the black hole formation scenario via stellar collisions. The defined black hole formation efficiency ($ε_{\mathrm{BH}}$) is shown to be influenced by the ratio of the initial stellar mass to the critical mass. △ Less

Submitted 11 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 16 pages, 5 figures

Hot Jupiter Formation in Dense Star Clusters

Authors: Leonard Benkendorff, Francesco Flammini Dotti, Katja Stock, Maxwell Xu Cai, Rainer Spurzem

Abstract: Hot Jupiters (HJ) are defined as Jupiter-mass exoplanets orbiting around their host star with an orbital period < 10 days. It is assumed that HJ do not form in-situ but ex-situ. Recent discoveries show that star clusters contribute to the formation of HJ. We present direct $N$-body simulations of planetary systems in star clusters and analyze the formation of HJ in them. We combine two direct $N$-… ▽ More Hot Jupiters (HJ) are defined as Jupiter-mass exoplanets orbiting around their host star with an orbital period < 10 days. It is assumed that HJ do not form in-situ but ex-situ. Recent discoveries show that star clusters contribute to the formation of HJ. We present direct $N$-body simulations of planetary systems in star clusters and analyze the formation of HJ in them. We combine two direct $N$-body codes: NBODY6++GPU for the dynamics of dense star clusters with 32 000 and 64 000 stellar members and LonelyPlanets used to follow 200 identical planetary systems around solar mass stars in those star clusters. We use different sets with 3, 4, or 5 planets and with the innermost planet at a semi-major axis of 5 au or 1 au and follow them for 100 Myr in our simulations. The results indicate that HJs are generated with high efficiency in dense star clusters if the innermost planet is already close to the host star at a semi-major axis of 1 au. If the innermost planet is initially beyond a semi-major axis of 5 au, the probability of a potential HJ ranges between $1.5-4.5$ percent. Very dense stellar neighborhoods tend to eject planets rather than forming HJs. A correlation between HJ formation and angular momentum deficit (AMD) is not witnessed. Young Hot Jupiters ($t_{\rm age} < 100$ Myrs) have only been found, in our simulations, in planetary systems with the innermost planet at a semi-major axis of 1 au. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Analysis of Kozai Cycles in Equal-Mass Hierarchical Triple Supermassive Black Hole Mergers in the Presence of a Stellar Cluster

Authors: Wei Hao, M. B. N. Kouwenhoven, Rainer Spurzem, Pau Amaro Seoane, Rosemary A. Mardling, Xiuming Xu

Abstract: Supermassive black holes (SMBHs) play an important role in galaxy evolution. Binary and triple SMBHs can form after galaxy mergers. A third SMBH may accelerate the SMBH merging process, possibly through the Kozai mechanism. We use N -body simulations to analyze oscillations in the orbital elements of hierarchical triple SMBHs with surrounding star clusters in galaxy centers. We find that SMBH trip… ▽ More Supermassive black holes (SMBHs) play an important role in galaxy evolution. Binary and triple SMBHs can form after galaxy mergers. A third SMBH may accelerate the SMBH merging process, possibly through the Kozai mechanism. We use N -body simulations to analyze oscillations in the orbital elements of hierarchical triple SMBHs with surrounding star clusters in galaxy centers. We find that SMBH triples spend only a small fraction of time in the hierarchical merger phase (i.e., a binary SMBH with a distant third SMBH perturber). Most of the time, the enclosed stellar mass within the orbits of the innermost or the outermost SMBH is comparable to the SMBH masses, indicating that the influence of the surrounding stellar population cannot be ignored. We search for Eccentric Kozai-Lidov (EKL) oscillations for which (i) the eccentricity of the inner binary and inclination are both oscillate and are anti-phase or in-phase and (ii) the oscillation period is consistent with EKL timescale. We find that EKL oscillations are short-lived and rare: the triple SMBH spends around 3% of its time in this phase over the ensemble of simulations, reaching around 8% in the best-case scenario. This suggests that the role of the EKL mechanism in accelerating the SMBH merger process may have been overestimated in previous studies. We follow-up with three-body simulations, using initial conditions extracted from the simulation, and the result can to some extent repeat the observed EKL-like oscillations. This comparison provides clues about why those EKL oscillations with perturbing stars are short-lived. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 24 pages, 16 figures, Published: 20 December 2023 on MNRAS https://doi.org/10.1093/mnras/stad3908

Growth of Seed Black Holes in Galactic Nuclei

Authors: Rainer Spurzem, Francesco Rizzuto, Manuel Arca Sedda, Albrecht Kamlah, Peter Berczik, Qi Shu, Ataru Tanikawa, Thorsten Naab

Abstract: The evolution of dense star clusters is followed by direct high-accuracy N-body simulation. The problem is to first order a gravitational N-body problem, but stars evolve due to astrophysics and the more massive ones form black holes or neutron stars as compact remnants at the end of their life. After including updates of stellar evolution of massive stars and for the relativistic treatment of bla… ▽ More The evolution of dense star clusters is followed by direct high-accuracy N-body simulation. The problem is to first order a gravitational N-body problem, but stars evolve due to astrophysics and the more massive ones form black holes or neutron stars as compact remnants at the end of their life. After including updates of stellar evolution of massive stars and for the relativistic treatment of black hole binaries we find the growth of intermediate mass black holes and we show that in star clusters binary black hole mergers in the so-called pair creation supernova (PSN) gap occur easily. Such black hole mergers have been recently observed by the LIGO-Virgo-KAGRA (LVK) collaboration, a network of ground based gravitational wave detectors. △ Less

Submitted 16 July, 2023; originally announced July 2023.

Comments: Article in Proceedings of NIC Symposium 2022, 14 pages, 5 figures

Journal ref: In: NIC Symposium 2022, M.Müller, Ch. Peter, A. Trautmann (Editors), ISBN 978-3-95806-646-5 (Book), DNB Link for this article: https://d-nb.info/1286760135

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

Influence of planets on debris disks in star clusters I: the 50 AU Jupiter

Authors: Kai Wu, M. B. N. Kouwenhoven, Rainer Spurzem, Xiaoying Pang

Abstract: Although debris disks may be common in exoplanet systems, only a few systems are known in which debris disks and planets coexist. Planets and the surrounding stellar population can have a significant impact on debris disk evolution. Here we study the dynamical evolution of debris structures around stars embedded in star clusters, aiming to determine how the presence of a planet affects the evoluti… ▽ More Although debris disks may be common in exoplanet systems, only a few systems are known in which debris disks and planets coexist. Planets and the surrounding stellar population can have a significant impact on debris disk evolution. Here we study the dynamical evolution of debris structures around stars embedded in star clusters, aiming to determine how the presence of a planet affects the evolution of such structures. We combine NBODY6++GPU and REBOUND to carry out N-body simulations of planetary systems in star clusters (N=8000; Rh=0.78 pc) for a period of 100 Myr, in which 100 solar-type stars are assigned 200 test particles. Simulations are carried out with and without a Jupiter-mass planet at 50 au. We find that the planet destabilizes test particles and speeds up their evolution. The planet expels most particles in nearby and resonant orbits. Remaining test particles tend to retain small inclinations when the planet is present, and fewer test particles obtain retrograde orbits. Most escaping test particles with speeds smaller than the star cluster's escape speed originate from cold regions of the planetary system or from regions near the planet. We identify three regions within planetary systems in star clusters: (i) the private region of the planet, where few debris particles remain (40 - 60 au), (ii) the reach of the planet, in which particles are affected by the planet (0 - 400 au), and (iii) the territory of the planetary system, most particles outside which will eventually escape (0 - 700 au). △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: 17 pages, 20 figures, and 3 tables. Accepted for publication in MNRAS

Computational Methods for Collisional Stellar Systems

Authors: Rainer Spurzem, Albrecht Kamlah

Abstract: Dense star clusters are spectacular self-gravitating stellar systems in our Galaxy and across the Universe - in many respects. They populate disks and spheroids of galaxies as well as almost every galactic center. In massive elliptical galaxies nuclear clusters harbor supermassive black holes, which might influence the evolution of their host galaxies as a whole. The evolution of dense star cluste… ▽ More Dense star clusters are spectacular self-gravitating stellar systems in our Galaxy and across the Universe - in many respects. They populate disks and spheroids of galaxies as well as almost every galactic center. In massive elliptical galaxies nuclear clusters harbor supermassive black holes, which might influence the evolution of their host galaxies as a whole. The evolution of dense star clusters is not only governed by the aging of their stellar populations and simple Newtonian dynamics. For increasing particle number, unique gravitational effects of collisional many-body systems begin to dominate the early cluster evolution. As a result, stellar densities become so high that stars can interact and collide, stellar evolution and binary stars change the dynamical evolution, black holes can accumulate in their centers and merge with relativistic effects becoming important. Recent high-resolution imaging has revealed even more complex structural properties with respect to stellar populations, binary fractions and compact objects as well as - the still controversial - existence of intermediate mass black holes in clusters of intermediate mass. Dense star clusters therefore are the ideal laboratory for the concomitant study of stellar evolution and Newtonian as well as relativistic dynamics. Not only the formation and disruption of dense star clusters has to be considered but also their galactic environments in terms of initial conditions as well as their impact on galactic evolution. This review deals with the specific computational challenges for modelling dense, gravothermal star clusters. △ Less

Submitted 19 May, 2023; originally announced May 2023.

Comments: 98 pages, 13 figures, invited article for Living Reviews in Computational Astrophysics

Tracing the Evolution of SMBHs and Stellar Objects in Galaxy Mergers: An Multi-mass Direct N-body Model

Authors: Shuo Li, Shiyan Zhong, Peter Berczik, Rainer Spurzem, Xian Chen, F. K. Liu

Abstract: By using direct N-body numerical simulations, we model the dynamical co-evolution of two supermassive black holes (SMBHs) and the surrounding stars in merging galaxies. In order to investigate how different stellar components evolve during the merger, we generate evolved stellar distributions with an initial mass function. Special schemes have also been developed to deal with some rare but interes… ▽ More By using direct N-body numerical simulations, we model the dynamical co-evolution of two supermassive black holes (SMBHs) and the surrounding stars in merging galaxies. In order to investigate how different stellar components evolve during the merger, we generate evolved stellar distributions with an initial mass function. Special schemes have also been developed to deal with some rare but interesting events, such as tidal disruption of main sequence stars, the plunge of low mass stars, white dwarfs, neutron stars and stellar mass black holes, and the partial tidal disruption of red giants or asymptotic giant branch stars. Our results indicate that the formation of a bound supermassive black hole binary (SMBHB) will enhance the capture rates of stellar objects by the SMBHs. Compared to the equal stellar mass model, the multi-mass model tends to result in a higher average mass of disrupted stars. Instead of being tidally disrupted by the SMBH, roughly half of the captured main sequence stars will directly plunge into the SMBH because of their small stellar radius. Giant stars, on the other hand, can be stripped of their envelopes if they are close enough to the SMBH. Though most remnants of the giant stars can survive after the disruption, a small fraction still could plunge into the SMBH quickly or after many orbital periods. Our results also indicate significant mass segregation of compact stars at the beginning of the merger, and then this effect is destroyed as the two SMBHs form a bound binary. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 36 pages, 11 figures, and 1 table. Accepted for publication in ApJ

Revisit the rate of tidal disruption events: the role of the partial tidal disruption event

Authors: Shiyan Zhong, Shuo Li, Peter Berczik, Rainer Spurzem

Abstract: Tidal disruption of stars in dense nuclear star clusters containing supermassive central black holes (SMBH) is modeled by high-accuracy direct N-body simulation. Stars getting too close to the SMBH are tidally disrupted and a tidal disruption event (TDE) happens. TDEs probe properties of SMBH, their accretion disks, and the surrounding nuclear stellar cluster. In this paper we compare rates of ful… ▽ More Tidal disruption of stars in dense nuclear star clusters containing supermassive central black holes (SMBH) is modeled by high-accuracy direct N-body simulation. Stars getting too close to the SMBH are tidally disrupted and a tidal disruption event (TDE) happens. TDEs probe properties of SMBH, their accretion disks, and the surrounding nuclear stellar cluster. In this paper we compare rates of full tidal disruption events (FTDE) with partial tidal disruption events (PTDE). Since a PTDE does not destroy the star, a leftover object emerges; we use the term 'leftover star' for it; two novel effects occur in the simulation: (1) variation of the leftover star's mass and radius, (2) variation of the leftover star's orbital energy. After switching on these two effects in our simulation, the number of FTDEs is reduced by roughly 28%, and the reduction is mostly due to the ejection of the leftover stars from PTDEs coming originally from relatively large distance. The number of PTDEs is about 75% higher than the simple estimation given by Stone et al. (2020), and the enhancement is mainly due to the multiple PTDEs produced by the leftover stars residing in the diffusive regime. We compute the peak mass fallback rate for the PTDEs and FTDEs recorded in the simulation, and find 58% of the PTDEs have peak mass fallback rate exceeding the Eddington limit, and the number of super-Eddington PTDEs is 2.3 times the number of super-Eddington FTDEs. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: 15 pages, 10 figures, accepted for publication in ApJ

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

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Authors: Katja Stock, Dimitri Veras, Maxwell X. Cai, Rainer Spurzem, Simon Portegies Zwart

Abstract: Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of th… ▽ More Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 16 pages, 18 figures, accepted for publication in MNRAS

Stability of planetary systems within the S-star cluster: the Solar system analogues

Authors: N. Davari, R. Capuzzo Dolcetta, R. Spurzem

Abstract: A dynamically relaxed dense cluster comprised of about 40 stars (the so-called S-stars) inhabits the central region of the Galaxy. Their stars revolve around the Sgr A$^*$ massive object. To understand the dynamical evolution of planetary systems in a particular environment like that around Sgr A$^*$, we carry out direct N--body simulations of planetary systems embedded in the S--star cluster. In… ▽ More A dynamically relaxed dense cluster comprised of about 40 stars (the so-called S-stars) inhabits the central region of the Galaxy. Their stars revolve around the Sgr A$^*$ massive object. To understand the dynamical evolution of planetary systems in a particular environment like that around Sgr A$^*$, we carry out direct N--body simulations of planetary systems embedded in the S--star cluster. In this work, we investigate the short-term stability of the planets orbiting around S-stars after their close interactions with the central massive black hole of our galaxy. We find that planetary systems go through encounters with the SMBH and the nearby stars. We determine the frequency and the strength of planetary systems' encounters with the nearby stars as these encounters remarkably increase for systems assigned to S-stars closer to the SMBH. The SMBH severely destabilizes the planetary systems, though we noted that the small oscillations in the mutual eccentricity and inclination of the planetary system could be caused by the planet--planet coupling and the near-resonance effect between the two planets. We obtain estimates of the fraction of survivor planets ($\sim 51\%$), and find that planets stripped from their hosting star are generally captured on close orbits around Sgr A$^*$. We notify while gas giants are tidally disrupted, terrestrial planets do not. We estimate that Sgr A$^*$ flares can be due to the tidal disruption events of \textit{starless} giant planets. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 13 pages, 8 figures

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.

Merging timescale for supermassive black hole binary in interacting galaxy NGC 6240

Authors: Margarita Sobolenko, Peter Berczik, Rainer Spurzem

Abstract: One of the main possible way of creating the supermassive black hole (SMBH) is a so call hierarchical merging scenario. Central SMBHs at the final phase of interacting and coalescing host-galaxies are observed as SMBH binary (SMBHB) candidates at different separations from hundreds of pc to mpc. Today one of the strongest SMBHB candidate is a ULIRG galaxy NGC 6240 which was X-ray spatially and spe… ▽ More One of the main possible way of creating the supermassive black hole (SMBH) is a so call hierarchical merging scenario. Central SMBHs at the final phase of interacting and coalescing host-galaxies are observed as SMBH binary (SMBHB) candidates at different separations from hundreds of pc to mpc. Today one of the strongest SMBHB candidate is a ULIRG galaxy NGC 6240 which was X-ray spatially and spectroscopically resolved by Chandra. Dynamical calculation of central SMBHBs merging in a dense stellar environment allows us to retrace their evolution from kpc to mpc scales. The main goal of our dynamical modeling was to reach the final, gravitational wave (GW) emission regime for the model BHs. We present the direct N-body simulations with up to one million particles and relativistic post-Newtonian corrections for the SMBHs particles up to 3.5 post-Newtonian terms. Generally speaking, the set of initial physical conditions can strongly effect of our merging time estimations. But in a range of our parameters, we did not find any strong correlation between merging time and BHs mass or BH to bulge mass ratios. Varying the model numerical parameters (such as particle number - N) makes our results quite robust and physically more motivated. From our 20 models we found the upper limit of merging time for central SMBHB is less than $\sim$55 Myr. This concrete number are valid only for our set of combination of initial mass ratios. Further detailed research of rare dual/multiple BHs in dense stellar environment (based on observations data) can clarify the dynamical co-evolution of central BHs and their host-galaxies. △ Less

Submitted 30 May, 2021; originally announced May 2021.

Comments: 11 pages, 6 figures, 3 tables, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 652, A134 (2021)

The bound mass of Dehnen models with centrally peaked star formation efficiency

Authors: B. Shukirgaliyev, A. Otebay, M. Sobolenko, M. Ishchenko, O. Borodina, T. Panamarev, S. Myrzakul, M. Kalambay, A. Naurzbayeva, E. Abdikamalov, E. Polyachenko, S. Banerjee, P. Berczik, R. Spurzem, A. Just

Abstract: Understanding the formation of star clusters with low star-formation efficiency (SFE) is very important to know about the star-formation history. In N-body models of star cluster evolution after gas expulsion, the Plummer model with outer power law density profile has been used massively. We study the impact of the density profile slopes on the survivability of the low-SFE star clusters after inst… ▽ More Understanding the formation of star clusters with low star-formation efficiency (SFE) is very important to know about the star-formation history. In N-body models of star cluster evolution after gas expulsion, the Plummer model with outer power law density profile has been used massively. We study the impact of the density profile slopes on the survivability of the low-SFE star clusters after instantaneous gas expulsion. We compare cases when stellar cluster has Plummer profile and Dehnen profiles with cusp of different slopes at the time of formation. We determine the corresponding density profile of the residual gas for a given global SFE, assuming that our model clusters formed with a constant efficiency per free-fall time and hence have shallower density profile of gas than that of stars. We perform direct $N$-body simulations of evolution of clusters initially in virial equilibrium within gas potential after gas removal. We find that the violent relaxation lasts no longer than 20~Myr independently of the density profile power law slopes. Dehnen model clusters survive after violent relaxation with significantly lower SFEs when the global SFE measured within the Jacobi radius or within a half-mass radius. Dehnen $γ=0$ model clusters show similar final bound fraction with the Plummer model clusters if global SFE is measured within 10 scale radii. The final bound fraction increases with $γ$ values for a given global SFE. We conclude that Dehnen clusters better resist the consequences of the violent relaxation followed the instantaneous gas expulsion than the Plummer clusters. Thus the shallower the outer density slope of the low-SFE clusters, the better for their survivability after gas expulsion. Among Dehnen clusters we find that the steeper the inner slope (cusp) the higher the bound mass fraction is retained after violent relaxation for a given global SFE. △ Less

Submitted 23 July, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 14 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 654, A53 (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

The long-term evolution of main-sequence binaries in DRAGON simulations

Authors: Qi Shu, Xiaoying Pang, Francesco Flammini Dotti, M. B. N. Kouwenhoven, Manuel Arca Sedda, Rainer Spurzem

Abstract: We present a comprehensive investigation of main-sequence (MS) binaries in the DRAGON simulations, which are the first one-million particles direct $N$-body simulations of globular clusters. We analyse the orbital parameters of the binary samples in two of the DRAGON simulations, D1-R7- IMF93 and D2-R7-IMF01, focusing on their secular evolution and correlations up to 12 Gyr. These two models have… ▽ More We present a comprehensive investigation of main-sequence (MS) binaries in the DRAGON simulations, which are the first one-million particles direct $N$-body simulations of globular clusters. We analyse the orbital parameters of the binary samples in two of the DRAGON simulations, D1-R7- IMF93 and D2-R7-IMF01, focusing on their secular evolution and correlations up to 12 Gyr. These two models have different initial stellar mass functions: Kroupa 1993 (D1-R7-IMF93) and Kroupa 2001 (D2-R7-IMF01); and different initial mass ratio distributions: random paring (D1-R7-IMF93) and a power-law (D1-R7-IMF93). In general, the mass ratio of a population of binaries increases over time due to stellar evolution, which is less significant in D2-R7-IMF01. In D1-R7-IMF93, primordial binaries with mass ratio $q \approx$ 0.2 are most common, and the frequency linearly declines with increasing $q$ at all times. Dynamical binaries of both models have higher eccentricities and larger semi-major axes than primordial binaries. They are preferentially located in the inner part of the star cluster. Secular evolution of binary orbital parameters does not depend on the initial mass-ratio distribution, but is sensitive to the initial binary distribution of the system. At t = 12 Gyr, the binary fraction decreases radially outwards, and mass segregation is present. A color difference of 0.1 mag in $F330W-F814W$ and 0.2 mag in $NUV-y$ between the core and the outskirts of both clusters is seen, which is a reflection of the binary radial distribution and the mass segregation in the cluster. The complete set of data for primordial and dynamical binary systems at all snapshot intervals is made publicly available. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 40 pages, 22 figures, accepted for publication in ApJS on date 29 November 2020

Exploring the origin of stars on bound and unbound orbits causing tidal disruption events

Authors: Shiyan Zhong, Kimitake Hayasaki, Shuo Li, Peter Berczik, Rainer Spurzem

Abstract: Tidal disruption events (TDEs) provide a clue to the properties of a central supermassive black hole (SMBH) and an accretion disk around it, and to the stellar density and velocity distributions in the nuclear star cluster surrounding the SMBH. Deviations of TDE light curves from the standard occurring at a parabolic encounter with the SMBH depends on whether the stellar orbit is hyperbolic or ecc… ▽ More Tidal disruption events (TDEs) provide a clue to the properties of a central supermassive black hole (SMBH) and an accretion disk around it, and to the stellar density and velocity distributions in the nuclear star cluster surrounding the SMBH. Deviations of TDE light curves from the standard occurring at a parabolic encounter with the SMBH depends on whether the stellar orbit is hyperbolic or eccentric (Hayasaki et al. 2018) and the penetration factor ($β$, tidal disruption radius to orbital pericenter ratio). We study the orbital parameters of bound and unbound stars being tidally disrupted by comparison of direct $N$-body simulation data with an analytical model. Starting from the classical steady-state Fokker-Planck model of Cohn & Kulsrud (1978), we develop an analytical model of the number density distribution of those stars as a function of orbital eccentricity ($e$) and $β$. To do so fittings of the density and velocity distribution of the nuclear star cluster and of the energy distribution of tidally disrupted stars are required and obtained from $N$-body data. We confirm that most of the stars causing TDEs in a spherical nuclear star cluster originate from the full loss-cone region of phase space, derive analytical boundaries in eccentricity-$β$ space, and find them confirmed by $N$-body data. Since our limiting eccentricities are much smaller than critical eccentricities for full accretion or full escape of stellar debris, we conclude that those stars are only very marginally eccentric or hyperbolic, close to parabolic. △ Less

Submitted 1 December, 2023; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 28 pages, 10 figures, 1 table. Accepted for publication in ApJ

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

Merging of unequal mass binary black holes in non-axisymmetric galactic nuclei

Authors: Peter Berczik, Manuel Arca Sedda, Margaryta Sobolenko, Marina Ishchenko, Olexander Sobodar, Rainer Spurzem

Abstract: In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem, varying both the total mass and the mass ratio of the SMBH binary (SMBHB). Simulations were carried out with the $\varphi-$GPU $N$-body code, which enabled us t… ▽ More In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem, varying both the total mass and the mass ratio of the SMBH binary (SMBHB). Simulations were carried out with the $\varphi-$GPU $N$-body code, which enabled us to fully exploit supercomputers equipped with graphic processing units (GPUs). As a model for the galactic nuclei, we adopted initial axisymmetric, rotating models, aimed at reproducing the properties of a galactic nucleus emerging from a galaxy merger event, containing two SMBHs which were unbound initially. We found no 'final-parsec problem', as our SMBHs tend to pair and shrink without showing significant signs of stalling. This confirms earlier results and extends them to large particle numbers and rotating systems. We find that the SMBHB hardening depends on the binary-reduced mass ratio via a single parameter function. Our results suggest that, at a fixed value for the SMBHB primary mass, the merger time of highly asymmetric binaries is up to four order of magnitudes smaller than the equal-mass binaries. This can significantly affect the population of SMBHs potentially detectable as gravitational wave sources. △ Less

Submitted 26 July, 2022; v1 submitted 10 August, 2020; originally announced August 2020.

Comments: 7 pages, 5 figures, accepted to A&A

Journal ref: A&A 665, A86 (2022)

Planetary Systems in a Star Cluster II: intermediate-mass black holes and planetary systems

Authors: Francesco Flammini Dotti, M. B. N. Kouwenhoven, Qi Shu, Wei Hao, Rainer Spurzem

Abstract: Most stars form in dense stellar environments. It is speculated that some dense star clusters may host intermediate-mass black holes (IMBHs), which may have formed from runaway collisions between high-mass stars, or from the mergers of less massive black holes. Here, we numerically explore the evolution of populations of planets in star clusters with an IMBH. We study the dynamical evolution of si… ▽ More Most stars form in dense stellar environments. It is speculated that some dense star clusters may host intermediate-mass black holes (IMBHs), which may have formed from runaway collisions between high-mass stars, or from the mergers of less massive black holes. Here, we numerically explore the evolution of populations of planets in star clusters with an IMBH. We study the dynamical evolution of single-planet systems and free-floating planets, over a period of 100~Myr, in star clusters without an IMBH, and in clusters with a central IMBH of mass $100~M_\odot$ or $200~M_\odot$. In the central region ($r\lesssim 0.2$~pc), the IMBH's tidal influence on planetary systems is typically 10~times stronger than the average neighbour star. For a star cluster with a $200M_\odot$ IMBH, the region in which the IMBH's influence is stronger within the virial radius ($\sim 1$~pc). The IMBH quenches mass segregation, and the stars in the core tend to move towards intermediate regions. The ejection rate of both stars and planets is higher when an IMBH is present. The rate at which planets are expelled from their host star rate is higher for clusters with higher IMBH masses, for $t<0.5 t_{rh}$, while remains mostly constant while the star cluster fills its Roche lobe, similar to a star cluster without an IMBH. The disruption rate of planetary systems is higher in initially denser clusters, and for wider planetary orbits, but this rate is substantially enhanced by the presence of a central IMBH. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 17 Pages, 39 Figures, accepted for pubblication on Monthly Notices of the Royal Astronomical Society on date 14 July 2020

On the survival of resonant and non-resonant planetary systems in star clusters

Authors: Katja Stock, Maxwell X. Cai, Rainer Spurzem, M. B. N. Kouwenhoven, Simon Portegies Zwart

Abstract: Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing… ▽ More Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing $N = 8k, 16k, 32k$, and $64k$ stars. The planetary systems, which consist of the four Solar system giant planets Jupiter, Saturn, Uranus, and Neptune, are initialized in different orbital configurations, to study the effect of the system architecture on the dynamical evolution of the entire planetary system, and on the escape rate of the individual planets. We find that the current orbital parameters of the Solar system giants (with initially circular orbits, as well as with present-day eccentricities) and a slightly more compact configuration, have a high resilience against stellar perturbations. A configuration with initial mean-motion resonances of 3:2, 3:2, and 5:4 between the planets, which is inspired by the Nice model, and for which the two outermost planets are usually ejected within the first $10^5$ yr, is in many cases stabilized due to the removal of the resonances by external stellar perturbation and by the rapid ejection of at least one planet. Assigning all planets the same mass of 1 Jovian mass almost equalizes the survival fractions. Our simulations reproduce the broad diversity amongst observed exoplanet systems. We find not only many very wide and/or eccentric orbits, but also a significant number of (stable) retrograde orbits. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 21 pages, 21 figures, accepted by MNRAS

Mass loss from massive globular clusters in tidal fields

Authors: Yohai Meiron, Jeremy J. Webb, Jongsuk Hong, Peter Berczik, Rainer Spurzem, Raymond G. Carlberg

Abstract: Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the… ▽ More Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for tens of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles. △ Less

Submitted 24 March, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

Linking the formation and fate of exo-Kuiper belts within solar system analogues

Authors: Dimitri Veras, Katja Reichert, Francesco Flammini Dotti, Maxwell X. Cai, Alexander J. Mustill, Andrew Shannon, Catriona H. McDonald, Simon Portegies Zwart, M. B. N. Kouwenhoven, Rainer Spurzem

Abstract: Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially em… ▽ More Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS

Direct n-body simulations of tidal disruption rate evolution in unequal-mass galaxy mergers

Authors: Shuo Li, Peter Berczik, Xian Chen, F. K. Liu, Rainer Spurzem, Yu Qiu

Abstract: The hierarchical galaxy formation model predicts supermassive black hole binaries (SMBHBs) in galactic nuclei. Due to the gas poor environment and the limited spatial resolution in observations they may hide in the center of many a galaxy. However, a close encounter of a star with one of the supermassive black holes (SMBHs) may tidally disrupt it to produce a tidal disruption event (TDE) and tempo… ▽ More The hierarchical galaxy formation model predicts supermassive black hole binaries (SMBHBs) in galactic nuclei. Due to the gas poor environment and the limited spatial resolution in observations they may hide in the center of many a galaxy. However, a close encounter of a star with one of the supermassive black holes (SMBHs) may tidally disrupt it to produce a tidal disruption event (TDE) and temporarily light up the SMBH. In a previous work, we investigated with direct N-BODY simulations the evolution of TDE rates of SMBHB systems in galaxy mergers of equal mass. In this work we extend to unequal mass mergers. Our results show that, when two SMBHs are far away from each other, the TDE rate of each host galaxy is similar as in an isolated galaxy. As the two galaxies and their SMBHs separation shrinks, the TDE rate is increasing gradually until it reaches a maximum shortly after the two SMBHs become bound. In this stage, the averaged TDE rate can be enhanced by several times to an order of magnitude relative to isolated single galaxies. Our simulations show that the dependence of the TDE accretion rate on the mass ratio in this stage can be well fitted by power law relations for both SMBHs. After the bound SMBHB forms, the TDE rate decreases with its further evolution. We also find that in minor mergers TDEs of the secondary SMBH during and after the bound binary formation are mainly contributed by stars from the other galaxy. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 27 pages, 6 figures, accepted by ApJ

Planetary systems in a star cluster I: the Solar system scenario

Authors: Francesco Flammini Dotti, M. B. N. Kouwenhoven, Maxwell Xu Cai, Rainer Spurzem

Abstract: Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planeta… ▽ More Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planetary systems escape from the star cluster within 50 Myr; those with low escape speeds often remain intact during and after the escape process. While most planetary systems inside the star cluster remain intact, a subset is strongly perturbed during the first 50 Myr. Over the course of time, 0.3 % - 5.3 % of the planets escape, sometimes up to tens of millions of years after a stellar encounter occurred. Survival rates are highest for Jupiter, while Uranus and Neptune have the highest escape rates. Unless directly affected by a stellar encounter itself, Jupiter frequently serves as a barrier that protects the terrestrial planets from perturbations in the outer planetary system. In low-density environments, Jupiter provides protection from perturbations in the outer planetary system, while in high-density environments, direct perturbations of Jupiter by neighbouring stars is disruptive to habitable-zone planets. The diversity amongst planetary systems that is present in the star clusters at 50 Myr, and amongst the escaping planetary systems, is high, which contributes to explaining the high diversity of observed exoplanet systems in star clusters and in the Galactic field △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 20 pages, 46 figures, accepted for pubblication on Monthly Notices of the Royal Astronomical Society on date 13 August 2019

Interaction of Stars Hosting Planets with Sgr A* Black hole

Authors: Nazanin Davari, Roberto Capuzzo-Dolcetta, Rainer Spurzem

Abstract: We present some preliminary results of our ongoing project about planetary systems around S-stars in the vicinity of Sgr A* black hole. Since S-stars might have migrated in the Galactic Centre (GC) from elsewhere, they probably still keep their planetary systems throughout their voyage. In this work, we study the destiny of their putative planetary systems after close interaction with the central… ▽ More We present some preliminary results of our ongoing project about planetary systems around S-stars in the vicinity of Sgr A* black hole. Since S-stars might have migrated in the Galactic Centre (GC) from elsewhere, they probably still keep their planetary systems throughout their voyage. In this work, we study the destiny of their putative planetary systems after close interaction with the central black hole of our galaxy. △ Less

Submitted 17 August, 2019; originally announced August 2019.

Comments: 4 pages, 3 figures. To appear in the Proceedings of the 351 IAU Symposium "Star Clusters: From the Milky Way to the Early Universe", 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"

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

Planets In Young Massive Clusters: On the survivability of planets in young massive clusters and its implication of planet orbital architectures in globular clusters

Authors: Maxwell Xu Cai, S. Portegies Zwart, M. B. N. Kouwenhoven, Rainer Spurzem

Abstract: As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (compar… ▽ More As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (comparable to Westerlund-1), our simulations show that most wide-orbit planets ($a\geq 20$~au) will be ejected within a timescale of 10 Myr. Interestingly, more than $70\%$ of planets with $a<5$~au survive in the 100 Myr simulations. Ignoring planet-planet scattering and tidal damping, the survivability at $t$ Myr as a function of initial semi-major axis $a_0$ in au in such a YMC can be described as $f_{\rm surv}(a_0, t)=-0.33 \log_{10}(a_0) \left(1 - e^{-0.0482t} \right) + 1$. Upon ejection, about $28.8\%$ of free-floating planets (FFPs) have sufficient speeds to escape from the host cluster at a crossing timescale. The other FFPs will remain bound to the cluster potential, but the subsequent dynamical evolution of the stellar system can result in the delayed ejection of FFPs from the host cluster. Although a full investigation of planet population in GCs requires extending the simulations to multi-Gyr, our results suggest that wide-orbit planets and free-floating planets are unlikely to be found in GCs. △ Less

Submitted 9 September, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: 12 pages, 11 figures, 1 table, accepted by MNRAS

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

Authors: Sambaran Banerjee, Krzysztof Belczynski, Christopher L. Fryer, Peter Berczik, Jarrod R. Hurley, Rainer Spurzem, Long Wang

Abstract: The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We… ▽ More The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster. △ Less

Submitted 9 April, 2020; v1 submitted 20 February, 2019; originally announced February 2019.

Comments: 27 pages, 15 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 639, A41 (2020)

Planetary Systems in Star Clusters: the dynamical evolution and survival

Authors: Francesco Flammini Dotti, Maxwell Xu Cai, Rainer Spurzem, M. B. N. Kouwenhoven

Abstract: Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary co… ▽ More Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary companions. To understand the origin and dynamical evolution of such exoplanet systems, it is necessary to carefully study the effect of their environments. Here, we combine theoretical estimates with state-of-the-art numerical simulations of evolving planetary systems similar to our own solar system in different star cluster environments. We combine the REBOUND planetary system evolution code, and the NBODY6++GPU star cluster evolution code, integrated in the AMUSE multi-physics environment. With our study we can constrain the effect of external perturbations of different environments on the planets and debris structures of a wide variety of planetary systems, which may play a key role for the habitability of exoplanets in the Universe. △ Less

Submitted 30 November, 2018; originally announced November 2018.

Comments: 2 pages, 2 figures, Proceedings for Symposia IAU 345 "Origins: from the Protosun to the First Steps of Life" at the XXX general assembly in Vienna

Direct N-body simulation of the Galactic centre

Authors: Taras Panamarev, Andreas Just, Rainer Spurzem, Peter Berczik, Long Wang, Manuel Arca Sedda

Abstract: We study the dynamics and evolution of the Milky Way nuclear star cluster performing a high resolution direct one-million-body simulation. Focusing on the interactions between such stellar system and the central supermassive black hole, we find that different stellar components adapt their overall distribution differently. After 5 Gyr, stellar mass black holes are characterized by a spatial distri… ▽ More We study the dynamics and evolution of the Milky Way nuclear star cluster performing a high resolution direct one-million-body simulation. Focusing on the interactions between such stellar system and the central supermassive black hole, we find that different stellar components adapt their overall distribution differently. After 5 Gyr, stellar mass black holes are characterized by a spatial distribution with power-slope $-1.75$, fully consistent with the prediction of Bahcall-Wolf pioneering work. Using the vast amount of data available, we infer the rate for tidal disruption events, being $4 \times 10^{-6}$ per yr, and estimate the number of objects that emit gravitational waves during the phases preceding the accretion onto the super-massive black hole, $\sim 270$ per Gyr. We show that some of these sources could form extreme mass-ratio inspirals. We follow the evolution of binary stars population, showing that the initial binary fraction of $5\%$ drops down to $2.5\%$ inside the inner parsec. Also, we explored the possible formation of binary systems containing a compact object, discussing the implications for millisecond pulsars formation and the development of Ia Supernovae. △ Less

Submitted 18 January, 2019; v1 submitted 6 May, 2018; originally announced May 2018.

Comments: Accepted by MNRAS after a major revision

Classification of Tidal Disruption Events Based on Stellar Orbital Properties

Authors: Kimitake Hayasaki, Shiyan Zhong, Shuo Li, Peter Berczik, Rainer Spurzem

Abstract: We study the rates of tidal disruption of stars by intermediate-mass to supermassive black holes on bound to unbound orbits by using high-accuracy direct N-body experiments. The approaching stars from the star cluster to the black hole can take three types of orbit: eccentric, parabolic, and hyperbolic orbits. Since the mass fallback rate shows a different variability depending on these orbital ty… ▽ More We study the rates of tidal disruption of stars by intermediate-mass to supermassive black holes on bound to unbound orbits by using high-accuracy direct N-body experiments. The approaching stars from the star cluster to the black hole can take three types of orbit: eccentric, parabolic, and hyperbolic orbits. Since the mass fallback rate shows a different variability depending on these orbital types, we can classify tidal disruption events (TDEs) into three main categories: eccentric, parabolic, and hyperbolic TDEs. Respective TDEs are characterized by two critical values of the orbital eccentricity: the lower critical eccentricity is the one below which the stars on eccentric orbits cause the finite, intense accretion, and the higher critical eccentricity above which the stars on hyperbolic orbits cause no accretion. Moreover, we find that the parabolic TDEs are divided into three subclasses: precisely parabolic, marginally eccentric, and marginally hyperbolic TDEs. We analytically derive that the mass fallback rate of the marginally eccentric TDEs can be flatter and slightly higher than the standard fallback rate proportional to $t^{-5/3}$, whereas it can be flatter and lower for the marginally hyperbolic TDEs. We confirm by N-body experiments that only few eccentric, precisely parabolic, and hyperbolic TDEs can occur in a spherical stellar system with a single intermediate-mass to supermassive black hole. A substantial fraction of the stars approaching to the black hole would cause the marginally eccentric or marginally hyperbolic TDEs. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 37 pages, 2 tables, 7 figures, accepted for publication in ApJ

Star-disc interaction in galactic nuclei: formation of a central stellar disc

Authors: Taras Panamarev, Bekdaulet Shukirgaliyev, Yohai Meiron, Peter Berczik, Andreas Just, Rainer Spurzem, Chingis Omarov, Emmanuil Vilkoviskij

Abstract: We perform high resolution direct $N$-body simulations to study the effect of an accretion disc on stellar dynamics in an active galactic nucleus (AGN). We show that the interaction of the nuclear stellar cluster (NSC) with the gaseous disc (AD) leads to formation of a stellar disc in the central part of the NSC. The accretion of stars from the stellar disc onto the super-massive black hole is bal… ▽ More We perform high resolution direct $N$-body simulations to study the effect of an accretion disc on stellar dynamics in an active galactic nucleus (AGN). We show that the interaction of the nuclear stellar cluster (NSC) with the gaseous disc (AD) leads to formation of a stellar disc in the central part of the NSC. The accretion of stars from the stellar disc onto the super-massive black hole is balanced by the capture of stars from the NSC into the stellar disc, yielding a stationary density profile. We derive the migration time through the AD to be 3\% of the half-mass relaxation time of the NSC. The mass and size of the stellar disc are 0.7\% of the mass and 5\% of the influence radius of the super-massive black hole. An AD lifetime shorter than the migration time would result in a less massive nuclear stellar disc. The detection of such a stellar disc could point to past activity of the hosting galactic nucleus. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: submitted to MNRAS after revision

Supermassive black holes coalescence mediated by massive perturbers: implications for gravitational waves emission and nuclear cluster formation

Authors: Manuel Arca Sedda, Peter Berczik, Roberto Capuzzo-Dolcetta, Giacomo Fragione, Margaryta Sobolenko, Rainer Spurzem

Abstract: A large fraction of galactic nuclei is expected to host supermassive black hole binaries (BHB), likely formed during the early phase of galaxies assembly and merging. In this paper, we use a large set of state-of-art numerical models to investigate the interplay between a BHB and a massive star cluster (GCs) driven toward the galactic centre by dynamical friction. Varying the BHB mass and mass rat… ▽ More A large fraction of galactic nuclei is expected to host supermassive black hole binaries (BHB), likely formed during the early phase of galaxies assembly and merging. In this paper, we use a large set of state-of-art numerical models to investigate the interplay between a BHB and a massive star cluster (GCs) driven toward the galactic centre by dynamical friction. Varying the BHB mass and mass ratio and the GC orbit, we show that the reciprocal feedback exerted between GCs and the BHB shapes their global properties. We show that, at GC-to-BHB mass ratios above 0.1, the GC affects notably the BHB orbital evolution, possibly boosting its coalescence. This effect is maximized if the GC moves on a retrograde orbit, and for a non-equal mass BHB. We show that the GC debris dispersed around the galactic nucleus can lead to the formation of a nuclear cluster, depending on the BHB tidal field, and that the distribution of compact remnants resulting from the GC disruption can carry information about the BHB orbital properties. We find that red giant stars delivered by the spiralling GC can be disrupted at a rate of $\simeq (0.7-7)\times 10^{-7}$ yr$^{-1}$ for BHB masses $\sim 10^7{\rm M}_\odot$, while tens to hundreds of stars can be possibly observed in the galactic halo as high-velocity stars, with velocities up to $\sim 2000$ km s$^{-1}$, depending on the BHB orbital properties. △ Less

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

Comments: 26 pages, 27 figures, 2 tables. Resubmitted to MNRAS after minor revision

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

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

Stability of Multiplanetary Systems in Star Clusters

Authors: Maxwell Xu Cai, M. B. N. Kouwenhoven, Simon F. Portegies Zwart, Rainer Spurzem

Abstract: Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems ar… ▽ More Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems are assigned to solar-type stars with $\sim 1 M_{\odot}$ and evolved for 50~Myr. We found that following the depletion of protoplanetary disks, external perturbations and planet-planet interactions are two driving mechanisms responsible for the destabilization of planetary systems. The planet survival rate varies from $\sim 95\%$ in the N=2k cluster to $\sim 60\%$ in the N=32k cluster, which suggests that most planetary systems can indeed survive in low-mass clusters, except in the central regions. We also find that planet ejections through stellar encounters are cumulative processes, as only $\sim 3\%$ of encounters are strong enough to excite the eccentricity by $Δe \geq 0.5$. Short-period planets can be perturbed through orbit crossings with long-period planets. When taking into account planet-planet interactions, the planet ejection rate nearly doubles, and therefore multiplicity contributes to the vulnerability of planetary systems. In each ensemble, $\sim 0.2\%$ of planetary orbits become retrograde due to random directions of stellar encounters. Our results predict that young low-mass star clusters are promising sites for next-generation planet surveys, yet low planet detection rates are expected in dense globular clusters such as 47 Tuc. Nevertheless, planets in denser stellar environments are likely to have shorter orbital periods, which enhances their detectability. △ Less

Submitted 12 June, 2017; originally announced June 2017.

Comments: 19 pages, 13 figures, 4 tables, accepted for publication in MNRAS

Planetary systems in star clusters

Authors: M. B. N. Kouwenhoven, Qi Shu, Maxwell Xu Cai, Rainer Spurzem

Abstract: Thousands of confirmed and candidate exoplanets have been identified in recent years. Consequently, theoretical research on the formation and dynamical evolution of planetary systems has seen a boost, and the processes of planet-planet scattering, secular evolution, and interaction between planets and gas/debris disks have been well-studied. Almost all of this work has focused on the formation and… ▽ More Thousands of confirmed and candidate exoplanets have been identified in recent years. Consequently, theoretical research on the formation and dynamical evolution of planetary systems has seen a boost, and the processes of planet-planet scattering, secular evolution, and interaction between planets and gas/debris disks have been well-studied. Almost all of this work has focused on the formation and evolution of isolated planetary systems, and neglect the effect of external influences, such as the gravitational interaction with neighbouring stars. Most stars, however, form in clustered environments that either quickly disperse, or evolve into open clusters. Under these conditions, young planetary systems experience frequent close encounters with other stars, at least during the first 1-10 Myr, which affects planets orbiting at any period range, as well as their debris structures. △ Less

Submitted 4 September, 2016; originally announced September 2016.

Comments: Proceedings of Cosmic-Lab: Star Clusters as Cosmic Laboratories for Astrophysics, Dynamics, and Fundamental Physics - MODEST16

Star-disc interaction in galactic nuclei: orbits and rates of accreted stars

Authors: Gareth F. Kennedy, Yohai Meiron, Bekdaulet Shukirgaliyev, Taras Panamarev, Peter Berczik, Andreas Just, Rainer Spurzem

Abstract: We examine the effect of an accretion disc on the orbits of stars in the central star cluster surrounding a central massive black hole by performing a suite of 39 high-accuracy direct N-body simulations using state-of-the art software and accelerator hardware, with particle numbers up to 128k. The primary focus is on the accretion rate of stars by the black hole (equivalent to their tidal disrupti… ▽ More We examine the effect of an accretion disc on the orbits of stars in the central star cluster surrounding a central massive black hole by performing a suite of 39 high-accuracy direct N-body simulations using state-of-the art software and accelerator hardware, with particle numbers up to 128k. The primary focus is on the accretion rate of stars by the black hole (equivalent to their tidal disruption rate for black holes in the small to medium mass range) and the eccentricity distribution of these stars. Our simulations vary not only the particle number, but disc model (two models examined), spatial resolution at the centre (characterised by the numerical accretion radius) and softening length. The large parameter range and physically realistic modelling allow us for the first time to confidently extrapolate these results to real galactic centres. While in a real galactic centre both particle number and accretion radius differ by a few orders of magnitude from our models, which are constrained by numerical capability, we find that the stellar accretion rate converges for models with N > 32k. The eccentricity distribution of accreted stars, however, does not converge. We find that there are two competing effects at work when improving the resolution: larger particle number leads to a smaller fraction of stars accreted on nearly-circular orbits, while higher spatial resolution increases this fraction. We scale our simulations to some nearby galaxies and find that the expected boost in stellar accretion (or tidal disruption, which could be observed as X-ray flares) in the presence of a gas disc is about a factor of 10. Even with this boost, the accretion of mass from stars is still a factor of ~ 100 slower than the accretion of gas from the disc. Thus, it seems accretion of stars is not a major contributor to black hole mass growth. △ Less

Submitted 18 April, 2016; originally announced April 2016.

Comments: Accepted for publication in MNRAS

The DRAGON simulations: globular cluster evolution with a million stars

Authors: Long Wang, Rainer Spurzem, Sverre Aarseth, Mirek Giersz, Abbas Askar, Peter Berczik, Thorsten Naab, Riko Schadow, M. B. N. Kouwenhoven

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

Submitted 1 February, 2016; originally announced February 2016.

Comments: 18 pages, 8 figures, MNRAS accepted

Journal ref: MNRAS 458, 1450-1465 (2016)

Million-Body Star Cluster Simulations: Comparisons between Monte Carlo and Direct $N$-body

Authors: Carl L. Rodriguez, Meagan Morscher, Long Wang, Sourav Chatterjee, Frederic A. Rasio, Rainer Spurzem

Abstract: We present the first detailed comparison between million-body globular cluster simulations computed with a Hénon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of… ▽ More We present the first detailed comparison between million-body globular cluster simulations computed with a Hénon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find excellent agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (> 1000) are retained for 12 Gyr. Thus, the very accurate direct $N$-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-$N$ dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct $N$-body, is capable of producing a very accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes. △ Less

Submitted 16 January, 2016; originally announced January 2016.

Comments: 9 pages, 7 figures, submitted to MNRAS

GalevNB: a conversion from N-BODY simulations to observations

Authors: Xiaoying Pang, Christoph Olczak, Difeng Guo, Rainer Spurzem, Ralf Kotulla

Abstract: We present $GalevNB$ (Galev for $N$-body simulations), an utility that converts fundamental stellar properties of $N$-body simulations into observational properties using the $GALEV$ (GAlaxy EVolutionary synthesis models) package, and thus allowing direct comparisons between observations and $N$-body simulations. It works by converting fundamental stellar properties, such as stellar mass, temperat… ▽ More We present $GalevNB$ (Galev for $N$-body simulations), an utility that converts fundamental stellar properties of $N$-body simulations into observational properties using the $GALEV$ (GAlaxy EVolutionary synthesis models) package, and thus allowing direct comparisons between observations and $N$-body simulations. It works by converting fundamental stellar properties, such as stellar mass, temperature, luminosity and metallicity into observational magnitudes for a variety of filters of mainstream instruments/telescopes, such as HST, ESO, SDSS, 2MASS, etc.), and into spectra that spans from far-UV (90 $\rm Å$) to near-IR (160 $\rm μ$m). As an application, we use $GalevNB$ to investigate the secular evolution of spectral energy distribution (SED) and color-magnitude diagram (CMD) of a simulated star cluster over a few hundred million years. With the results given by $GalevNB$ we discover an UV-excess in the SED of the cluster over the whole simulation time. We also identify four candidates that contribute to the FUV peak, core helium burning stars, thermal pulsing asymptotic giant branch (TPAGB) stars, white dwarfs and naked helium stars. △ Less

Submitted 19 September, 2015; originally announced September 2015.

Comments: 13 pages, 8 figures, accepted for publication in Research in Astronomy and Astrophysics (RAA)

Boosted Tidal Disruption by Massive Black Hole Binaries During Galaxy Mergers FROM The View of N-Body Simulation

Authors: Shuo Li, F. K. Liu, Peter Berczik, Rainer Spurzem

Abstract: Supermassive black hole binaries (SMBHBs) are productions of the hierarchical galaxy formation model. There are many close connections between central SMBH and its host galaxy because the former plays very important roles on the formation and evolution of a galaxy. For this reason, the evolution of SMBHBs in merging galaxies is an essential problem. Since there are many discussions about SMBHB evo… ▽ More Supermassive black hole binaries (SMBHBs) are productions of the hierarchical galaxy formation model. There are many close connections between central SMBH and its host galaxy because the former plays very important roles on the formation and evolution of a galaxy. For this reason, the evolution of SMBHBs in merging galaxies is an essential problem. Since there are many discussions about SMBHB evolution in gas rich environment, we focus on the quiescent galaxy, using tidal disruption as a diagnostic tool. Our study is based on a series of numerical large particle number direct N-body simulations for dry major mergers. According to the simulation results, the evolution can be divided into three phases. In phase I, the tidal disruption rate for two well separated SMBHs in merging system has similar level to single SMBH in isolate galaxy. After two SMBHs getting close enough to form a bound binary in phase II, the disruption rate can be enhanced for ~ 2 order of magnitudes within a short time. This "boosted" disruption stage finishes after the SMBHB evolving to compact binary system in phase III, corresponding to a drop back of disruption rate to a level of a few times higher than that in Phase I. How to correctly extrapolate our N-body simulation results to reality, and implications of our results to observations, are discussed too. △ Less

Submitted 20 December, 2016; v1 submitted 1 September, 2015; originally announced September 2015.

Comments: 16 pages, Accepted for publication in ApJ

Supermassive Black Holes in Galactic Nuclei with Tidal Disruption of Stars: Paper II - Axisymmetric Nuclei

Authors: Shiyan Zhong, Peter Berczik, Rainer Spurzem

Abstract: Tidal Disruption of stars by supermassive central black holes from dense rotating star clusters is modelled by high-accuracy direct N-body simulation. As in a previous paper on spherical star clusters we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, now according to several classes of orbits which only occur in axisymmetric systems (short… ▽ More Tidal Disruption of stars by supermassive central black holes from dense rotating star clusters is modelled by high-accuracy direct N-body simulation. As in a previous paper on spherical star clusters we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, now according to several classes of orbits which only occur in axisymmetric systems (short axis tube and saucer). Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss-cone in phase space which is raised from the fact that total angular momentum $\bf J$ is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals interesting features. Inside the influence radius of the supermassive black hole the angular distribution of disrupted stars has a conspicuous bimodal structure with a local minimum near the equatorial plane. Outside the influence radius this dependence is weak. We show that the bimodal structure of orbital parameters can be explained by the presence of two families of regular orbits, namely short axis tube and saucer orbits. Also the consequences of our results for the loss cone in axisymmetric galactic nuclei are presented. △ Less

Submitted 12 August, 2015; originally announced August 2015.

Comments: 14 pages, 16 figures, accepted by ApJ

Performance analysis of parallel gravitational $N$-body codes on large GPU cluster

Authors: Siyi Huang, Rainer Spurzem, Peter Berczik

Abstract: We compare the performance of two very different parallel gravitational $N$-body codes for astrophysical simulations on large GPU clusters, both pioneer in their own fields as well as in certain mutual scales - NBODY6++ and Bonsai. We carry out the benchmark of the two codes by analyzing their performance, accuracy and efficiency through the modeling of structure decomposition and timing measureme… ▽ More We compare the performance of two very different parallel gravitational $N$-body codes for astrophysical simulations on large GPU clusters, both pioneer in their own fields as well as in certain mutual scales - NBODY6++ and Bonsai. We carry out the benchmark of the two codes by analyzing their performance, accuracy and efficiency through the modeling of structure decomposition and timing measurements. We find that both codes are heavily optimized to leverage the computational potential of GPUs as their performance has approached half of the maximum single precision performance of the underlying GPU cards. With such performance we predict that a speed-up of $200-300$ can be achieved when up to 1k processors and GPUs are employed simultaneously. We discuss the quantitative information about comparisons of two codes, finding that in the same cases Bonsai adopts larger time steps as well as relative energy errors than NBODY6++, typically ranging from $10-50$ times larger, depending on the chosen parameters of the codes. While the two codes are built for different astrophysical applications, in specified conditions they may overlap in performance at certain physical scale, and thus allowing the user to choose from either one with finetuned parameters accordingly. △ Less

Submitted 11 August, 2015; originally announced August 2015.

Comments: 15 pages, 7 figures, 3 tables, accepted for publication in Research in Astronomy and Astrophysics (RAA)