Search

In this study, we investigate the contribution of millisecond pulsars (MSPs) to the gamma-ray excess observed in the Galactic Center by analyzing data from high-resolution direct N-body simulations of six globular clusters (GCs) that experience close encounters with the nuclear star cluster. Using the {\phi}-GPU code, we tracked the orbits of individual neutron stars (NSs) formed during the simulations, assuming a fraction of these NSs evolve into MSPs. Our model includes state-of-the-art single stellar evolution code including prescription for neutron star formation. We estimated the gamma-ray flux from these MSPs, considering known values for their gamma-ray emission. Our results show that MSPs originating from the six modeled GCs contribute a small but non-negligible fraction of the observed gamma-ray flux. This finding suggests that the actual gamma-ray flux from MSPs could be much higher when considering the entire population of GCs, potentially significantly contributing to the gamma-ray excess. This study highlights the importance of considering MSPs in the Galactic Center, originating from nearby globular clusters, as a potential source of the observed gamma-ray excess. Future work will involve more sophisticated simulations incorporating binary stellar evolution and comparing the fraction of MSPs in observed GCs to refine our models and improve the accuracy of our estimates., Comment: in Russian language

Intermediate Mass Black Hole (IMBH) mergers with masses $10^4 - 10^6$ $M_{\odot}$ are expected to produce gravitational waves (GWs) detectable by the Laser Interferometer Space Antenna (LISA) with high signal to noise ratios out to redshift 20. IMBH mergers are expected to take place within dwarf galaxies, however, the dynamics, timescales, and effect on their hosts are largely unexplored. In a previous study, we examined how IMBHs would pair and merge within nucleated dwarf galaxies. IMBHs in nucleated hosts evolve very efficiently, forming a binary system and coalescing within a few hundred million years. Although the fraction of dwarf galaxies ($10^7$ M$_{\odot} \leq$ $M_{\star} \leq 10^{10}$ M$_{\odot}$) hosting nuclear star clusters is between 60-100\%, this fraction drops to 20-70\% for lower-mass dwarfs ($M_{\star}\approx 10^7$ M$_{\odot}$), with the largest drop in low-density environments. Here, we extend our previous study by performing direct $N-$body simulations to explore the dynamics and evolution of IMBHs within non-nucleated dwarf galaxies, under the assumption that IMBHs exist within these dwarfs. To our surprise, none of IMBHs in our simulation suite merge within a Hubble time, despite many attaining high eccentricities $e \sim 0.7-0.95$. We conclude that extremely low stellar density environments in the centers of non-nucleated dwarfs do not provide an ample supply of stars to interact with IMBHs binary resulting in its stalling, in spite of triaxiality and high eccentricity, common means to drive a binary to coalescence. Our findings underline the importance of considering all detailed host properties to predict IMBH merger rates for LISA., Comment: Submitted to ApJ

Membership analysis is an important tool for studying star clusters. There are various approaches to membership determination, including supervised and unsupervised machine learning (ML) methods. We perform membership analysis using the supervised machine learning approach. We train and test our ML models on two sets of star cluster data: snapshots from $N$-body simulations and 21 different clusters from the Gaia Data Release 3 data. We explore five different ML models: Random Forest (RF), Decision Trees, Support Vector Machines, Feed-Forward Neural Networks, and K-Nearest Neighbors. We find that all models produce similar results, with RF showing slightly better accuracy. We find that a balance of classes in datasets is optional for successful learning. The classification accuracy depends strongly on the astrometric parameters. The addition of photometric parameters does not improve performance. We do not find a strong correlation between the classification accuracy and clusters' age, mass, and half-mass radius. At the same time, models trained on clusters with a larger number of members generally produce better results., Comment: Accepted for publication in A&A. 10 pages, 12 figures

Based on the available observational data from the literature, we analysed the dynamics of the NGC 6240 galaxy central supermassive black hole (SMBH) system. For the dynamical modelling of this triple SBMH system, we used the massively parallel and GPU accelerated phi-GPU direct summation N-body code. Following a long-timescale modelling of the triple system, we carried out a very detailed time output analysis of the von Zeipel-Lidov-Kozai (ZLK) oscillations for the black holes. According to our Newtonian simulation results, for all models and randomisations, the bound system from S1+S2 components formed at ~3.6 Myr. The formation of the bound hierarchical triple system S+N occurred at ~18 Myr. Over the course of these Newtonian simulations of the evolution of the triple SMBH system and the surrounding environment in NGC 6240, ZLK oscillations were detected (in most cases) for the binary components. The inclination angle between the orbital angular momentum of binary components aptly coincides with the theoretical calculations of the ZLK mechanism. In our set of randomised 15 Newtonian $N$-body dynamical galaxy models in 13 systems, we were able to detect a ZLK mechanism. In contrast, our extra few-body post-Newtonian runs (for one randomisation case) show it is only for the large inner binary initial eccentricity (in our case >0.9 that we are able to observe the possibility of the inner binary merging, due to the post-Newtonian energy radiation effects. For the lower eccentricity cases, the test runs show no sign of possible merging or any ZLK oscillations in the system., Comment: 14 pages, 15 figures, accepted for publication to the Astronomy and Astrophysics

Context. Based on the Gaia DR3, we reconstructed the orbital evolution of the known Milky Way globular clusters and found that six objects, NGC 6681, NGC 6981, Palomar 6, NGC 6642, HP 1, and NGC 1904, very likely interact closely with the nuclear star cluster. Aims. We study the dynamical evolution of selected Milky Way globular clusters and their interactions with the Galactic centre over cosmological timescales. We examine the global dynamical mass loss of these globular cluster systems, their close interactions with the Galactic centre, and the potential capture of stars by the Milky Way nuclear star cluster. Methods. For the dynamical modelling of the clusters, we used the parallel N-body code phi-GPU, which allows star-by-star simulations of the systems. Our current code also enabled us to follow the stellar evolution of individual particles, including the formation of high-mass remnants. The modelling was carried out in a Milky Way-like, time-variable potential (with a dynamically changing mass and scale length), obtained from the IllustrisTNG-100 database, with a full integration time of eight billion years. Results. Based on extensive numerical modelling and analysis, we estimated the mass loss and the global and inner structures of the selected six clusters. Over an evolution of eight billion years, the clusters lost 80% of their initial mass. We analysed the phase-space evolution of the individual unbound stars NGC 6681, NGC 6642, HP 1, and NGC 1904. We found that only NGC 6642 could potentially have been a source for populating the Milky Way nuclear star cluster in the past., Comment: 26 pages, 16 figures, 5 Tables, 4 Appendixes. Accepted for publication to the Astronomy and astrophysics

In a previous paper using Gaia DR2 data, we demonstrated that the two closely situated open clusters Collinder 135 and UBC 7 might have formed together about 50 Myr ago. In this work, we performed star-by-star dynamical modelling of the evolution of the open clusters Collinder 135 and UBC 7 from their supposed initial state to their present-day state, reproducing observational distributions of members. Modelling of the Collinder 135 and UBC 7 dynamical evolution was done using the high-order parallel N-body code \phi-GPU with up-to-date stellar evolution. Membership and characteristics of the clusters were acquired based on Gaia DR3 data. The comparison of the present-day radial cumulative star count obtained from the N-body simulations with the current observational data gave us full consistency of the model with observational data, especially in the central 8 pc, where 80% of the stars reside. The proper motion velocity components obtained from the N-body simulations of the stars are also quite consistent with the observed distributions and error bars. These results show that our numerical modelling is able to reproduce the open clusters' current complex 6D observed phase-space distributions with a high level of confidence. Thus, the model demonstrates that the hypothesis of a common origin of Collinder 135 and UBC 7 complies with present-day observational data., Comment: Astronomy & Astrophysics, accepted

The present epoch of the Gaia success gives us a possibility to predict the dynamical evolution of our Solar System in the global Galactic framework with high precision. We statistically investigated the total interaction of globular clusters with the Solar System during six billion years of look-back time. We estimated the gravitational influence of globular clusters' flyby onto the Oort cloud system. To perform the realistic orbital dynamical evolution for each individual cluster, we used our own high-order parallel dynamical $N$ body $\varphi$-GPU code that we developed. To reconstruct the orbital trajectories of clusters, we used five external dynamical time variable galactic potentials selected from the IllustrisTNG-100 cosmological database and one static potential. To detect a cluster's close passages near the Solar System, we adopted a simple distance criterion of below 200 pc. To take into account a cluster's measurement errors (based on Gaia DR3), we generated 1000 initial positions and velocity randomizations for each cluster in each potential. We found 35 globular clusters that have had close passages near the Sun in all the six potentials during the whole lifetime of the Solar System. We can conclude that at a relative distance of 50 pc between a GC and the SolS, we obtain on average $\sim 15$\% of the close passage probability over all six billion years, and at $dR=100$ pc, we get on average $\sim 35$\% of the close passage probability over all six billion years. The globular clusters BH_140, UKS_1, and Djorg_1 have a mean minimum relative distance to the Sun of 9, 19, and 17 pc, respectively. We can assume that a globular cluster with close passages near the Sun is not a frequent occurrence but also not an exceptional event in the Solar System's lifetime., Comment: 11 pages, 6 figures, 3 tables, accepted to the Astronomy & Astrophysics

Aims. We carry out the self-consistent dynamic evolution of the orbital structure of Milky Way globular clusters. This allows us to estimate possible and probable close passages and even collisions of the clusters with each other. Methods. We reproduced the orbits of 147 globular clusters in 10 Gyr lookback time using our own high-order N-body parallel dynamic phi-GPU code. The initial conditions (three coordinates and three velocities for the present time) were derived from the Gaia DR3 catalogue. The galaxy is represented by five external potentials from the IllustrisTNG-100, whose masses and sizes of the disk and halo components are similar to the physical values of the Milky Way at present. Results. We present a statistical analysis of the cumulative close passages rate: About ten close passages with relative distances shorter than 50 pc for every billion years for each of the five external potentials. We present the 22 most reliable collision pairs with a good probability. As an example: Terzan 4 versus Terzan 2 (49%), Terzan 4 versus NGC 6624 (44%), Terzan 4 versus Terzan 5 (40%), Terzan 4 versus NGC 6440 (40%), and Terzan 4 versus Liller 1 (42%). The most active globular cluster in the collision sense is Terzan 4, which has 5.65 collision events on average (averaged over all individual 1000 initial condition realisations). Most collisions are located inside the Galactic disk and form two ring-like structures. The first ring-like structure has the highest collision number density at 1 kpc, and the second sturcture has a maximum at 2 kpc. Conclusions. Based on our numerical simulations, we can conclude that the few dozen Milky Way globular clusters probably undergo some close encounters and even possible collisions during their lifetimes, which can significantly affect their individual dynamical evolution and possibly even their stellar content., Comment: 11 pages, 9 figures, 3 tables, accepted for publication in Astronomy and Astrophysics

For a pair of supermassive black holes (SMBHs) in the remnant of a dual galaxy merger, well-known models exist to describe their dynamical evolution until the final coalescence accompanied by the emission of a low-frequency gravitational wave (GW) signal. In this article, we investigate the dynamical evolution of three SMBH triple systems recovered from the ROMULUS25 cosmological simulation to explore common dynamical evolution patterns and assess typical coalescence times. For this purpose, we construct initial conditions from the ROMULUS25 data and perform high-resolution gravitodynamical \N-body simulations. We track the orbital evolution from the galactic inspiral to the formation of hard binaries at sub-parsec separation and use the observed hardening rates to project the time of coalescence. In all cases, the two heaviest black holes form an efficiently hardening binary that merges within fractions of the Hubble time. The lightest SMBH either gets ejected, forms a stable hierarchical triple system with the heavier binary, forms a hardening binary with the previously merged binary's remnant, or remains on a wide galactic orbit. The coalescence times of the lighter black holes are thus significantly longer than for the heavier binary, as they experience lower dynamical friction and stellar hardening rates. We observe the formation of hierarchical triples when the density profile of the galactic nucleus is sufficiently steep., Comment: 11 pages, 6 figures, 4 tables, accepted for publication in A&A

Context: Potentially hazardous asteroids (PHA) in Earth-crossing orbits pose a constant threat to life on Earth. Several mitigation methods have been proposed, and the most feasible technique appears to be the disintegration of the impactor and the generation of a fragment cloud by explosive penetrators at interception. However, mitigation analyses tend to neglect the effect of orbital dynamics on the trajectory of fragments. Aims: We aim to study the effect of orbital dynamics of the impactor's cloud on the number of fragments that hit the Earth, assuming different interception dates. We investigate the effect of self-gravitational cohesion and the axial rotation of the impactor. Methods: We computed the orbits of 10^5 fragments with a high-precision direct N-body integrator of the eighth order, running on GPUs. We considered orbital perturbations from all large bodies in the Solar System and the self-gravity of the cloud fragments. Results: Using a series of numerical experiments, we show that orbital shear causes the fragment cloud to adopt the shape of a triaxial ellipsoid. The shape and alignment of the triaxial ellipsoid are strongly modulated by the cloud's orbital trajectory and, hence, the impact cross-section of the cloud with respect to the Earth. Therefore, the number of fragments hitting the Earth is strongly influenced by the orbit of the impactor and the time of interception. A minimum number of impacts occur for a well-defined orientation of the impactor rotational axis, depending on the date of interception. Conclusions: To minimise the lethal consequences of an PHA's impact, a well-constrained interception timing is necessary. A too-early interception may not be ideal for PHAs in the Apollo or Aten groups. Thus, we find that the best time to intercept PHA is when it is at the pericentre of its orbit., Comment: Accepted for publication in A&A Letters (some refs. update)

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., Comment: Article in Proceedings of NIC Symposium 2022, 14 pages, 5 figures

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., Comment: 22 pages, 14 figures, 3 tables. Comments welcome. Submitted to MNRAS

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 $\chi_{\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$., Comment: 15 pages, 6 figures, 2 tables, 1 appendix. Comments welcome. Submitted to MNRAS

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$)., Comment: 22 pages, 21 figures, 4 tables. Comments welcome. Submitted to MNRAS

Binary black holes (BHs) can be formed dynamically in the centers of star clusters. The high natal kicks for stellar-mass BHs used in previous works made it hard to retain BHs in star clusters. Recent studies of massive star evolution and supernovae (SN) propose kick velocities that are lower due to the fallback of the SN ejecta. We study the impact of these updates by performing $N$-body simulations following instantaneous gas expulsion. For comparison, we simulate two additional model sets with the previous treatment of stars: one with high kicks and another with artificial removal of the kicks. Our model clusters initially consist of about one hundred thousand stars, formed with centrally-peaked efficiency. We find that the updated treatment of stars, due to the fallback-scaled lower natal kicks, allows clusters to retain SN remnants after violent relaxation. The mass contribution of the retained remnants does not exceed a few percent of the total bound cluster mass during the early evolution. For this reason, the first giga year of evolution is not affected significantly by this effect. Nevertheless, during the subsequent long-term evolution, the retained BHs accelerate mass segregation, leading to the faster dissolution of the clusters., Comment: 8 pages, 7 figures, submitted to MNRAS, comments are welcome

Context. Recent observational data show that the Milky Way (MW) galaxy contains about 170 globular clusters (GCs). A fraction of them is likely formed in dwarf galaxies accreted onto the MW in the past, while the remaining of clusters are formed in-situ. Therefore, different parameters, including orbits, of the globular clusters is a valuable tool for studying the Milky Way evolution. However, since the evolution of the 3D mass distribution of the MW is poorly constrained, the orbits of the clusters are usually calculated in static potentials. Aims. In this work, we study the evolution of the GCs in several external potentials, where we aim to quantify the effects of the evolving galaxy potential on the orbits of the GCs. Methods. For the orbits calculation we used five MW-like potentials from IllustrisTNG-100 simulation. The orbits of 159 GCs were integrated using a high-order N-body parallel dynamic code phi-GPU, with initial conditions obtained from recent Gaia DR3 catalogues. Results. We provide a classification of the GCs orbits according to their 3D shapes and association with different components of the MW (disk, halo, bulge). We also found that the globular clusters in the external potentials have roughly similar energy-angular momentum distributions at the present time. However, both total energy and total angular momentum of the GCs are not conserved due to time-varying nature of the potentials. In some extreme cases, the total energy can change up to 40% (18 objects) over the last 5 Gyr of evolution. We found that the in-situ formed GCs are less affected by the evolution of the TNG potentials as compared to the clusters which are likely formed ex-situ. Therefore, our results suggest that time-varying potentials significantly affect the orbits of the GC, thus making it vital for understanding the formation of the MW., Comment: 19 pages, 17 figures, 2 tables, accepted for publication in Astronomy and Astrophysics

Aims. We estimate the dynamical evolution of the Globular Clusters interaction with the Galactic centre that dynamically changed in the past. Methods. We simulated the orbits of 147 globular clusters over 10 Gyr lookback time using the parallel N-body code phi-GPU. For each globular cluster, we generated 1000 sets of initial data with random proper motions and radial velocities based on the observed values. To distinguish globular clusters interacting with the galactic centre, we used the criterion of a relative distance of less than 100 pc. We used four external potentials from the IllustrisTNG-100 database, which were selected for their similarity to the present-day Milky Way, to simulate the structure of the Galaxy at different times. Results. We obtained 3-4 globular cluster interactions per Gyr at distances of less than 50 pc and 5-6 interactions per Gyr at distances of less than 80 pc among the studied 147 globular clusters that had close passages near the Galactic centre. We selected 10 of them for detailed study and found almost 100% probability of interaction with the Galactic centre for six of them. Conclusions. According to our results, the maximum interaction frequency of globular clusters with the Galactic centre in the Milky Way is likely to be a few dozens of passages per Gyr within a central zone of 100 pc. This low frequency may not be sufficient to fully explain the relatively high mass (of order 10^7 Msol) of the nuclear star cluster in the Milky Way, if we consider only the periodic capture of stars from globular clusters during close encounters. Therefore, we must also consider the possibility that some early globular clusters were completely tidally disrupted during interactions with the forming nuclear star cluster and the Galactic centre., Comment: 25 pages, 21 figures, 5 tables. Accepted for publication in Astronomy and Astrophysics

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., Comment: 36 pages, 11 figures, and 1 table. Accepted for publication in ApJ

This study is aimed at investigating the dynamic evolution of the orbits of stellar globular clusters (GCs). To integrate the orbits backward in time, the authors use models of the time-varying potentials derived from cosmological simulations, which are closest to the potential of the Galaxy. This allows for estimating the probability of close passages (collisions) of GCs with respect to each other and the Galactic center (GalC) in the Galaxy undergoing dynamic changes in the past. To reproduce the dynamics of the Galaxy in time, five potentials selected from the IllustrisTNG-100 large-scale cosmological database, which are similar in their characteristics to the current physical parameters of the Milky Way, are used. With these time-varying potentials, we have reproduced the orbital trajectories of 143 GCs 10 Gyr back in time using our original phi-GPU N-body code. Each GC was treated as a single physical particle with the assigned position and velocity of the cluster center from the Gaia DR2 observations. For each of the potentials, 1000 initial conditions were generated with randomized initial velocities of GCs within the errors of the observational data. In this study, we consider close passages to be passages with a relative distance of less than 100 pc and a relative speed of less than 250 km/s. To select clusters that pass at close distances from the GalC, the following criterion is applied based only on the relative distance: it must be less than 100 pc. Applying the above criteria, the authors obtained statistically significant rates of close passages of GCs with respect to each other and to the GalC. It has been determined that GCs during their evolution have approximately 10 intersecting trajectories with each other on the average and approximately 3 to 4 close passages near the GalC in 1 Gyr at a distance of 50 pc for each of the chosen potentials., Comment: 18 pages, 7 figures, 5 tables, accepted in Kinematics and Physics of Celestial Bodies

The main idea of our research is to estimate the physical coalescence time of the double supermassive black hole (SMBH) system in the centre of NGC 6240 based on the X-ray observations from the Chandra space observatory. The spectra of the Northern and Southern nuclei were fitted by spectral models from Sherpa and both presented the narrow component of the Fe K$\alpha$ emission line. It enabled us to apply the spectral model to these lines and to find relative offset $\approx0.02$ keV. The enclosed dynamical mass of the central region of NGC 6240 with radius 1 kpc was estimated $\approx 2.04\times 10^{11} \rm\; M_{\odot}$. These data allowed us to carry on the high resolution direct N-body simulations with Newtonian and post-Newtonian (up to $2.5\mathcal{PN}$ correction) dynamics for this particular double SMBH system. As a result, from our numerical models we approximated the central SMBH binary merging time for the different binary eccentricities. In our numerical parameters range the upper limit for the merging time, even for the very small eccentricities, is still below $\approx70$ Myr. Gravitational waveforms and amplitude-frequency pictures from such events can be detected using Pulsar Timing Array (PTA) projects at the last merging phase., Comment: 12 pages, 7 figures, 6 tables. Accepted for publication in MNRAS, minor typos were corrected

The recent ALMA maps together with observations of H$_2$O maser emission seem to suggest the presence of a counter-rotation in the obscuring torus of NGC 1068. We propose to explain this phenomenon as due to the influence of a wind, considered as radiation pressure, and the effects of torus orientation. In order to test this idea: 1. we make $N$-body simulation of a clumpy torus taking into account mutual forces between particles (clouds); 2. we apply ray-tracing algorithm with the beams from the central engine to choose the clouds in the torus throat that can be under direct influence of the accretion disk emission; 3. we use semi-analytical model to simulate the influence of the asymmetrical radiation pressure (wind) forced on the clouds in the torus throat. An axis of such a wind is tilted with respect to the torus symmetry axis; 4. we orient the torus relative to an observer and again apply ray-tracing algorithm. In this step the beams go from an observer to the optically thick clouds that allows us to take into account the mutual obscuration of clouds; 5. after projecting on the picture plane, we impose a grid on the resulting cloud distribution and find the mean velocity of clouds in each cells to mimic the ALMA observational maps. By choosing the parameters corresponding to NGC 1068 we obtain the model velocity maps that emulate the effect of an apparent counter-rotation and can explain the discovery made by ALMA., Comment: 11 pages, 11 figures, Accepted for publication in MNRAS

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., Comment: 15 pages, 10 figures, accepted for publication in ApJ

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., Comment: 19 pages, 9 figures, 2 tables, Accepted by MNRAS

The purpose of the present work is a detailed investigation of the dynamical evolutionof Collinder 135 and UBC 7 star clusters. We present a set of dynamical numerical simulationsusing realistic star clusterN-body modeling technique with the forward integration of thestar-by-star cluster models to the present day, based on best-available 3D coordinates and velocitiesobtained from the latest Gaia EDR3 data release. We have established that Collinder 135 and UBC 7 are probably a binary star cluster and have common origin. We carried out a full star-by-star N-body simulation of the stellar population of both clusters using the new algorithm of Single Stellar Evolution and performed a comparison of the results obtained in the observational data (like cumulative number counts), which showed a fairly good agreement., Comment: 6 pages, 6 figures, 3 tables

Coalescing supermassive black hole binaries (BHBs) are expected to be the loudest sources of gravitational waves (GWs) in the Universe. Detection rates for ground or space-based detectors based on cosmological simulations and semi-analytic models are highly uncertain. A major difficulty stems from the necessity to model the BHB from the scale of the merger to that of inspiral. Of particular relevance to the GW merger timescale is the binary eccentricity. Here we present a self-consistent numerical study of the eccentricity of BHBs formed in massive gas-free mergers from the early stages of the merger to the hardening phase, followed by a semi-analytical model down to coalescence. We find that the early eccentricity of the unbound black hole pair is largely determined by the initial orbit. It systematically decreases during the dynamical friction phase. The eccentricity at binary formation is affected by stochasticity and noise owing to encounters with stars, but preserves a strong correlation with the initial orbital eccentricity. Binding of the black holes is a phase characterised by strong perturbations, and we present a quantitative definition of the time of binary formation. During hardening the eccentricity increases in minor mergers, unless the binary is approximately circular, but remains largely unchanged in major mergers, in agreement with predictions from semi-analytical models based on isotropic scattering experiments. Coalescence times due to hardening and GW emission in gas-poor non-rotating ellipticals are <~0.5 Gyr for the large initial eccentricities (0.5 < e < 0.9) typical of galaxy mergers in cosmological simulations., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

Using the data from Gaia (ESA) Data Release 2 we performed the orbital calculations of globular clusters (GCs) of the Milky Way. To explore possible collisions between the GCs, using our developed highorder {\phi}-GRAPE code, we integrated (backwards and forward) the orbits of 119 objects with reliable positions and proper motions. In calculations, we adopted a realistic axisymmetric Galactic potential (bulge + disk + halo). Using different impact conditions, we found five pairs of the GCs that likely experienced collisions: Terzan 3 - NGC 6553, Terzan 3 - NGC 6218, Liller 1 - NGC 6522, Djorg 2 - NGC 6552 and NGC 6355 - NGC 6637. We analyzed the GCs interaction rates with the central supermassive black hole. Assuming the maximum 100 pc distance criteria for separation between them we estimated 11 close encounter events. From our numerical simulations, we estimate the close interaction rate as at least one event per Gyr with the impact parameter less than 30 pc; and one event per Myr with the impact parameter less than 60 pc. Our calculations show one very close encounter of NGC 6121 with the central SMBH near 5.5 pc (practically direct collision). Based on the extended literature search for the possible progenitor of our selected 11 GCs, we found that most of them have a Milky Way main bulge origin., Comment: 9 pages, 9 figures, accepted for publication in Reports of National academy of sciences of the Republic of Kazakhstan

Using the data from Gaia (ESA) Data Release 2 we performed the orbital calculations of globular clusters (GCs) of the Milky Way. To explore possible close encounters (or collisions) between the GCs, using our own developed high-order phi-GRAPE code, we integrated backward and forward the orbits of 119 objects with reliable positions and proper motions. In calculations, we adopted a realistic axisymmetric Galactic potential (bulge + disk + halo). Using different impact conditions, we found four pairs of the six GCs that may have experienced an encounter within twice the sum of the half-mass radii (collisions) over the last 5Gyr: Terzan3 - NGC 6553, Terzan 3 - NGC 6218, Liller 1 - NGC 6522 and Djorg 2 - NGC 6553., Comment: Accepted for publication in Advances in Astronomy and Space Physics (http://aasp.kiev.ua). arXiv admin note: substantial text overlap with arXiv:2201.06891

Aims: We study the impact of Post-Newtonian correction terms on the energetic interaction between a gravitational wave (GW)-emitting supermassive black hole (SMBH) binary system and incoming stars via three-body scattering experiments. Methods: We use the AR-chain code to simulate with high accuracy the interactions between stars and an SMBH binary at separations of$\sim 1000\,R_{sch}$. For all of the interactions, we investigate in detail the energy balance of the three-body systems, using both Newtonian and Post-Newtonian expressions for the SMBH binary orbital energy, taking into account the GW emission by the binary. Results: We find that at these separations, purely Newtonian treatment of the binary orbital energy is insufficient to properly account for the SMBH binary orbital evolution. Instead, along with GW emission, even terms in the PN-corrections must be included in order to describe the energy change of the binary during the stellar interaction.es, 2, Comment: 8 pages, 2 figures, 1 table. Submitted to A&A

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., Comment: 11 pages, 6 figures, 3 tables, Accepted for publication in Astronomy & Astrophysics

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 $\gamma=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 $\gamma$ 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., Comment: 14 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics (A&A)

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., Comment: 31 pages, 15 figures, 5 tables, Accepted by MNRAS

We have developed the dynamical model of a clumpy torus in an active galactic nucleus (AGN) and compared to recent ALMA observations. We present $N$-body simulations of a torus in the field of a supermassive black hole (SMBH), made of up to $N=10^5$ gravitationally interacting clouds. As initial conditions, we choose random distributions of the orbital elements of the clouds with a cut-off in the inclination to mimic the presence of wind cones produced at the early AGN stage. When the torus reaches an equilibrium, it has a doughnut shape. We discuss the presence of box orbits. We have then constructed the velocity and velocity dispersion maps using the resulting distributions of the clouds at equilibrium. The effects of torus inclination and cloud sizes are duly analyzed. We discuss the obscuration effects of the clouds using a ray tracing simulation matching the model maps to ALMA resolution. By comparing the model with the observational maps of NGC 1068 we find that the SMBH mass is $M_\text{smbh}=5\times 10^6 M_\odot$ for the range of the torus inclination angles $45^\circ - 60^\circ$. We also construct the velocity dispersion maps for NGC 1326 and NGC 1672. They show that the peaks in the ALMA dispersion maps are related to the emission of the torus throat. Finally, we obtain the temperature distribution maps with parameters that correspond to our model velocity maps for NGC 1068. They show stratification in temperature distribution with the shape of the high temperature region as in the VLTI/MIDI map., Comment: 16 pages, 19 figures, 1 table, accepted for publication in MNRAS

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 ($\beta$, 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 $\beta$. 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-$\beta$ 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., Comment: 28 pages, 10 figures, 1 table. Accepted for publication in ApJ

Aims: We investigate the orbital and phase space properties of loss cone stars that interact strongly with a hard, high-redshift binary supermassive black hole (SMBH) system formed in a cosmological scenario. Methods: We use a novel hybrid integration approach that combines the direct N-body code $\varphi$-GRAPE with ETICS, a collisionless code that employs the self-consistent field method for force calculation. The hybrid approach shows considerable speed-up over direct summation for particle numbers $> 10^6$, while retaining accuracy of direct N-body for a subset of particles. During the SMBH binary evolution we monitor individual stellar interactions with the binary in order to identify stars that noticeably contribute to the SMBH binary hardening. Results: We successfully identify and analyze in detail the properties of stars that extract energy from the binary. We find that the summed energy changes seen in these stars match very well with the overall binary energy change, demonstrating that stellar interactions are the primary drivers of SMBH binary hardening in triaxial, gas-poor systems. The slight triaxiality of our system results in efficient loss cone refilling, avoiding the final parsec problem. We distinguish three different populations of interactions based on their apocenter. We find a clear prevalence of interactions co-rotating with the binary. Nevertheless, retrograde interactions are the most energetic, contributing only slightly less than the prograde population to the overall energy exchange. The most energetic interactions are also likely to result in a change of sign in the angular momentum of the star. We estimate the merger timescale of the binary to be $\approx 20$ $\mathrm{Myr}$, a value larger by a factor of two than the timescale reported in a previous study., Comment: 17 pages, 13 figures, 2 tables; accepted for publication in Astronomy & Astrophysics

Modern galaxy formation theory suggests that the misalignment between stellar and gaseous components usually results from an external gas accretion and/or interaction with other galaxies. The extreme case of the kinematic misalignment is demonstrated by so-called galaxies with counterrotation that possess two distinct components rotating in opposite directions with respect to each other. We provide an in-deep analysis of galaxies with counterrotation from IllustrisTNG100 cosmological simulations. We have found $25$ galaxies with substantial stellar counterrotation in the stellar mass range of $2\times10^{9}-3\times10^{10}$~\Msun. In our sample the stellar counterrotation is a result of an external gas infall happened $\approx 2-8$~Gyr ago. The infall leads to the initial removal of pre-existing gas, which is captured and mixed together with the infalling component. The gas mixture ends up in the counterrotating gaseous disc. We show that $\approx 90\%$ of the stellar counterrotation formed in-situ, in the counterrotating gas. During the early phases of the infall, gas can be found in inclined extended and rather thin disc-like structures, and in some galaxies they are similar to (nearly-)~polar disc or ring-like structures. We discuss a possible link between the gas infall, AGN activity and the formation of misaligned components. In particular, we suggest that the AGN activity does not cause the counterrotation, although it is efficiently triggered by the retrograde gas infall, and it correlates well with the misaligned component appearance. We also find evidence of the stellar disc heating visible as an increase of the vertical-to-radial velocity dispersion ratio above unity in both co- and counterrotating components, which implies the importance of the kinematical misalignment in shaping the velocity ellipsoids in disc galaxies., Comment: 20 pages, 20 figures, accepted for publication in MNRAS

This study presents a comprehensive chemo-dynamical analysis of LAMOST J1109+0754, a bright (V = 12.8), extremely metal-poor (\abund{Fe}{H} = $-3.17$) star, with a strong \textit{r}-process enhancement (\abund{Eu}{Fe} = +0.94 $\pm$ 0.12). Our results are based on the 7-D measurements supplied by $Gaia$ and the chemical composition derived from a high-resolution ($R\sim 110,000$), high signal-to-noise ratio ($S/N \sim 60)$ optical spectrum obtained by the 2.4\,m Automated Planet Finder Telescope at Lick Observatory. We obtain chemical abundances of 31 elements (from lithium to thorium). The abundance ratios (\abund{X}{Fe}) of the light-elements (Z $\leqslant 30$) suggest a massive Population\,III progenitor in the 13.4-29.5\,M$_\odot$ mass range. The heavy-element ($30 <$ Z $\leqslant 90$) abundance pattern of J1109+075 agrees extremely well with the scaled-Solar \textit{r}-process signature. We have developed a novel approach to trace the kinematic history and orbital evolution of J1109+0754 with a c\textbf{O}smologically de\textbf{RI}ved tim\textbf{E}-varyi\textbf{N}g Galactic po\textbf{T}ential (the ORIENT) constructed from snapshots of a simulated Milky-Way analog taken from the \texttt{Illustris-TNG} simulation. The orbital evolution within this Milky Way-like galaxy, along with the chemical-abundance pattern implies that J1109+0754 likely originated in a low-mass dwarf galaxy located $\sim$ 60\,kpc from the center of the Galaxy, which was accreted $\sim$ 6 - 7\,Gyr ago, and that the star now belongs to the outer-halo population.

Given the closeness of the two open clusters Cr 135 and UBC 7 on the sky, we investigate the possibility of the two clusters to be physically related. We aim to recover the present-day stellar membership in the open clusters Collinder 135 and UBC 7 (300 pc from the Sun), to constrain their kinematic parameters, ages and masses, and to restore their primordial phase space configuration. The most reliable cluster members are selected with our traditional method modified for the use of Gaia DR2 data. Numerical simulations use the integration of cluster trajectories backwards in time with our original high order Hermite4 code \PGRAPE. We constrain the age, spatial coordinates and velocities, radii and masses of the clusters. We estimate the actual separation of the cluster centres equal to 24 pc. The orbital integration shows that the clusters were much closer in the past if their current line-of-sight velocities are very similar and the total mass is more than 7 times larger the mass of the determined most reliable members. We conclude that the two clusters Cr 135 and UBC 7 might very well have formed a physial pair, based on the observational evidence as well as numerical simulations. The probability of a chance coincidence is only about $2\%$., Comment: Accepted for publication as a Letter in Astronomy and Astrophysics. 5 pages, 2 figures

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 $\rho_\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., Comment: Submitted

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., Comment: 7 pages, 5 figures, accepted to A&A

We present a set of self-consistent chemo-dynamical simulations of MW-type galaxies formation to study the origin of the bimodality of $\alpha$-elements in stellar populations. We explore how the bimodality is related to the geometrically and kinematically defined stellar discs, gas accretion and radial migration. We find that the two $\alpha$-sequences are formed in quite different physical environments. The high-$\alpha$ sequence is formed early from a burst of star formation (SF) in a turbulent, compact gaseous disc which forms a thick disc. The low-$\alpha$ stellar populations is the result of quiescent SF supported by the slow accretion of enriched gas onto a radially extended thin disc. Stellar feedback-driven outflows during the formation of the thick disc are responsible for the enrichment of the surrounding gaseous halo, which subsequently feeds the disc on a longer time-scale. During the thin disc phase, chemical evolution reaches an equilibrium metallicity and abundance, where the stars pile-up. This equilibrium metallicity decreases towards the outer disc, generating the ridgeline that forms the low-$\alpha$ sequence. We identify a second mechanism capable of creating a low-$\alpha$ sequence in one of our simulations. Rapid shutdown of the SF, provoked by the feedback at the end of the thick disc phase, suppresses the chemical enrichment of the halo gas, which, once accreted onto the star-forming disc, dilutes the ISM at the beginning of the thin disc formation. Both mechanisms can operate in a galaxy, but the former is expected to occur when SF efficiency ceases to be dominated by the formation of the thick disc, while the latter can occur in the inner regions. Being the result of the presence of low and high gas density environments, the bimodality is independent of any particular merger history, suggesting that it could be much more widespread than has been claimed., Comment: 22 pages, 19 figures, accepted for publication in MNRAS

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.

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations, we confirm previous findings and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to a repopulation of the loss cone and early instability saturation., Comment: 7 pages, 4 figures, MNRAS submitted

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations of a toy model, we confirm previous findings for the harmonic dominating potential and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in the spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to an early instability saturation., Comment: 6 pages, 9 figures, MNRAS accepted

We derive the photometric, kinematic, and abundance characteristics of 18 star-forming MaNGA galaxies with fairly regular velocity fields and surface brightness distributions and with a large offset between the measured position angles of the major kinematic and photometric axes, dPA > 20 degree. The aim is to examine if there is any other distinctive characteristic common to these galaxies. We found morphological signs of interaction in some (in 11 out of 18) but not in all galaxies. The observed velocity fields show a large variety; the maps of the isovelocities vary from an hourglass-like appearance to a set of straight lines. The position angles of the major kinematic axes of the stellar and gas rotations are close to each other. The values of the central oxygen abundance, radial abundance gradient, and star formation rate are distributed within the intervals defined by galaxies with small (no) dPA of similar mass. Thus, we do not find any specific characteristic common to all galaxies with large dPA. Instead, the properties of these galaxies are similar to those of galaxies with small (no) dPA. This suggests that either the reason responsible for the large dPA does not influence other characteristics or the galaxies with large dPA do not share a common origin, they can, instead, originate through different channels., Comment: accepted for publication in Astronomy and Astrophysics

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., Comment: 27 pages, 6 figures, accepted by ApJ

We study the evolution of star clusters in the Galactic tidal field starting from their birth in molecular clumps. Our model clusters form according to the local-density-driven cluster formation model in which the stellar density profile is steeper than that of gas. As a result, clusters resist the gas expulsion better than predicted by earlier models. We vary the impact of the Galactic tidal field {\lambda}, considering different Galactocentric distances (3-18 kpc), as well as different cluster sizes. Our model clusters survive the gas expulsion independent of {\lambda}. We investigated the relation between the cluster mass at the onset of secular evolution and their dissolution time. The model clusters formed with a high star-formation efficiency (SFE) follow a tight mass-dependent dissolution relation, in agreement with previous theoretical studies. However, the low-SFE models present a shallower mass-dependent relation than high-SFE clusters, and most dissolve before reaching 1 Gyr (cluster teenage mortality)., Comment: 4 pages, 2 figures. Submitted for publication in the Proceedings of IAU Symposium 351

We derive the circular velocity curves of the gaseous and stellar discs of UGC 4056, a giant barred galaxy with an active galactic nucleus (AGN). We analyse UGC 4056 using the 2D spectroscopy obtained within the framework of the Mapping Nearby Galaxies at APO (MaNGA) survey. Using images and the colour index g-r from the Sloan Digital Sky Survey (SDSS), we determined the tilt of the galaxy, which allows us to conclude that the galaxy rotates clockwise with trailing spiral arms. We found that the gas motion at the central part of the UGC 4056 shows peculiar features. The rotation velocity of the gaseous disc shows a bump within around three kiloparsecs while the rotation velocity of the stellar disc falls smoothly to zero with decreasing galactocentric distance. We demonstrate that the peculiar radial velocities in the central part of the galaxy may be caused by the inflow of the gas towards the nucleus of the galaxy. The unusual motion of the gas takes place at the region with the AGN-like radiation and can be explained by the gas response to the bar potential., Comment: 11 pages, 10 figures, 1 table, accepted for publication in A&A

We study the impact of the tidal field on the survivability of star clusters following instantaneous gas expulsion. Our model clusters are formed with a centrally-peaked star-formation efficiency profile as a result of star-formation taking place with a constant efficiency per free-fall time. We define the impact of the tidal field as the ratio of the cluster half-mass radius to its Jacobi radius immediately after gas expulsion, $\lambda = r_{h}/R_{J}$. We vary $\lambda$ by varying either the Galactocentric distance, or the size (hence volume density) of star clusters. We propose a new method to measure the violent relaxation duration, in which we compare the total mass-loss rate of star clusters with their stellar evolutionary mass-loss rate. That way, we can robustly estimate the bound mass fraction of our model clusters at the end of violent relaxation. The duration of violent relaxation correlates linearly with the Jacobi radius, when considering identical clusters at different Galactocentric distances. In contrast, it is nearly constant for the solar neighbourhood clusters, slightly decreasing with $\lambda$. The violent relaxation does not last longer than 50 Myr in our simulations. Identical model clusters placed at different Galactocentric distances have the same final bound fraction, despite experiencing different impacts of the tidal field. The solar neighbourhood clusters with different densities experience only limited variations of their final bound fraction. In general, we conclude that the cluster survivability after instantaneous gas expulsion, as measured by their bound mass fraction at the end of violent relaxation, $F_{bound}$, is independent of the impact of the tidal field, $\lambda$., Comment: accepted for publication in MNRAS, 8 pages, 5 figures,3 tables

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., Comment: 27 pages, 15 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics

Using a single N-body simulation ($N=0.14\times 10^9$) we explore the formation, evolution and spatial variation of the phase-space spirals similar to those recently discovered by Antoja et al. in the Milky Way disk, with Gaia DR2. For the first time in the literature, we use a self-consistent N-body simulation of an isolated Milky Way-type galaxy to show that the phase-space spirals develop naturally from vertical oscillations driven by the buckling of the stellar bar. We claim that the physical mechanism standing behind the observed incomplete phase-space mixing process can be internal and not necessarily due to the perturbation induced by a massive satellite. In our model, the bending oscillations propagate outwards and produce axisymmetric variations of the mean vertical coordinate and of the vertical velocity component. As a consequence, the phase-space wrapping results in the formation of patterns with various morphology across the disk, depending on the bar orientation, distance to the galactic center and time elapsed since the bar buckling. Once bending waves appear, they are supported for a long time via disk self-gravity. The underlying physical mechanism implies the link between in-plane and vertical motion that leads directly to phase-space structures whose amplitude and shape are in remarkable agreement with those of the phase-space spirals observed in the Milky Way disk. In our isolated galaxy simulation, phase-space spirals are still distinguishable, at the solar neighbourhood, 3 Gyr after the buckling phase. The long-lived character of the phase-space spirals generated by the bar buckling instability cast doubts on the timing argument used so far to get back at the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several Gyrs ago, and not as recent as suggested so far., Comment: 7 pages, 6 figures; accepted for publication by A&A Letters on Dec 14, 2018