Based on our current high resolution direct N-body modelling of the Milky Way typical Star Cluster systems dynamical evolution we try to numerically estimate the influence of individual spin values and orientations on gravitational wave (GW) waveforms and observed time-frequency maps during multiple cycles for binary black hole (BBH) mergers. In our up to date N-body dynamical simulations we use the high order relativistic post-Newtonian corrections for the BH binary particles (3.5 post-Newtonian (PN) terms including spin-spin and spin-orbit terms). In the current work, we present the GW waveforms catalogue which covers the large parameter space in mass ratios 0.05 - 0.82 and extreme possible individual spin cases.

The Orion Nebula Cluster (ONC) is one of the nearest open clusters, which we can directly compare to numerical simulations. We performed a simulation of star cluster formation similar to the ONC using our new N-body/smoothed particle hydrodynamics code, ASURA+BRIDGE. We found that the hierarchical formation of star clusters via clump mergers can explain the observed three peaks in the stellar age distribution as well as the dynamically anisotropic structures of the ONC.

We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-calledMulti-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range fields, ensuring that simulation time scales with N log N in the number of particles, as opposed to N2 for direct N-body. The collisional relaxation effects are modelled by computing particle interactions based on a collision operator approach that ensures rigorous conservation of energy and momenta and depends only on particles velocities and cell-based integrated quantities.

The purpose of the present work is a detailed investigation of the dynamical evolution of Collinder 135 and UBC 7 star clusters. We present a set of dynamical numerical simulations using realistic star cluster -body modeling technique with the forward integration of the star-by-star cluster models to the present day, based on best-available 3D coordinates and velocities obtained 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.

Magnetic field plays an important role in star formation and galaxy evolution. Previous studies discussed about the origin of magnetic field and its effect to the environment. With the recent advancement of supercomputers, adding the magnetic field to a cosmological hydrodynamic simulations only become feasible. In this proceeding, we present the results of high-resolution magneto-hydrodynamic simulation with GIZMO and compare our simulation result with the previous literature and the observations.

By means of numerical simulations we study the radial-orbit instability in anisotropic self-gravitating N–body systems under the effect of noise. We find that the presence of additive or multiplicative noise has a different effect on the onset of the instability, depending on the initial value of the orbital anisotropy.

For the past 25 years, we have been considering the Stark effect for neutral helium lines in DB white dwarfs using the standard Stark broadening theory in both the impact regime (in the center of the lines) and the quasi-static regime (in the wings) for the electrons, while neglecting the effect of ions in motion. Although this is probably a good approximation based on previous theoretical work, the transition between the two regimes for the electrons and the contribution of the ions very near the core might be poorly represented. To better represent these particularities, we report here the results of a new series of simulations that treat the local dynamics and interactions of both electrons and ions around a neutral helium atom. From these simulations, we produce new improved line profiles, which we compare with our previous analytical results.

One of the major and widely known small scale problem with the Lambda CDM model of cosmology is the “core-cusp” problem. In this study we investigate whether this problem can be resolved using bar instabilities. We see that all the initial bars are thin (b/a < 0.3) in our simulations and the bar becomes thick ( b /a > 0.3) faster in the high resolution simulations. By increasing the resolution, we mean a larger number of disk particles. The thicker bars in the high resolution simulations transfer less angular momentum to the halo. Hence, we find that in the high resolution simulations it takes around 7 Gyr for the bar to remove inner dark matter cusp which is too long to be meaningful in galaxy evolution timescales. Physically, the reason is that as the resolution increases, the bar buckles faster and becomes thicker much earlier on.

Evolution and the formation of bars in the galactic disks is studied in the context of Modified Gravity (MOG) by using N-body simulations. It is found that changing the value of free parameters of the model can effectively alter the strength of the bar and disk’s stability.

We study the minor mergers of galaxies using simulations. For this we use GADGET2 code. We present results of simulations of minor mergers of disc galaxies of mass ratio 1:10. These simulations consist of collisionless as well as hydrodynamical runs including a gaseous component in the galactic disc of primary galaxy. Our goal is to establish the characteristics of discs obtained after the merger.We observe that the primary galaxy discs are not destroyed after the merger. We take different initial conditions for the primary galaxy varying the gas percentage in disc from 0–40 percentage and study the thickness of the disc after the merger. We generally observe that the thickness of the disc increases after the merger for any gas percentage. We also observe that as the gas percentage increases in the disc of initial primary galaxy, the increase in the thickness keeps decreasing.

We describe an algorithm that can fit the properties of the dwarf galaxy progenitor of a tidal stream, given the properties of that stream. We show that under ideal conditions (the Milky Way potential, the orbit of the dwarf galaxy progenitor, and the functional form of the dwarf galaxy progenitor are known exactly), the density and angular width of stars along the stream can be used to constrain the mass and radial profile of both the stellar and dark matter components of the progenitor dwarf galaxy that was ripped apart to create the stream. Our provisional fit for the parameters of the dwarf galaxy progenitor of the Orphan Stream indicates that it is less massive and has fewer stars than previous works have indicated.

We investigate the minor interactions of two disk galaxies with mass ratio of 10:1 in fly-by encounters that do not lead to the merging of the galaxies. In our N-body simulations, we vary only the pericenter distances to see the effect of the fly-by on the bulge of the major galaxy over the course of the trajectory. At different time steps of the evolution, we did two-dimensional fittings of disk, bulge and bar to trace the variation in the sersic index of the bulge. Our results suggest that galaxy bulges can become boxy/disky through flyby interactions of galaxies.

We collected radial velocities of more than 50.000 individual stars in 156 Galactic globular clusters (GGC) and matched them with HST photometry and Gaia DR2 proper motions. This allowed us to derive the GGC’s mean proper motions and space velocities. By fitting a large set of N-body simulations to their velocity dispersion and surface density profiles, combined with new measurements of their internal radially dependent mass functions, we have determined their present-day masses and structural parameters, and for 144 GGCs their internal kinematics. We also derive the initial cluster masses by calculating the cluster orbits backwards in time applying suitable recipes to account for mass-loss and dynamical friction. The new fundamental parameters of GGCs are publicly available via an online database, which will regularly be updated.

Using direct N-body simulations of self-gravitating systems we study the dependence of dynamical chaos on the system size N. We find that the N-body chaos quantified in terms of the largest Lyapunov exponent Λmax decreases with N. The values of its inverse (the so-called Lyapunov time tλ) are found to be smaller than the two-body collisional relaxation time but larger than the typical violent relaxation time, thus suggesting the existence of another collective time scale connected to many-body chaos.

We investigate the dynamics of supermassive black holes (SMBHs) in galactic cores by means of a semi-analytic model based on the Langevin equation, including dynamical friction and stochastic noise accounting for the gravitational interactions with stars. The model is validated against direct N-body simulations of intermediate-mass black holes in stellar clusters where a realistic number of particles is accessible. For the galactic case, we find that the SMBH experiences a Brownian-like motion with a typical displacement from the geometric center of the Galaxy of a few parsecs, for system parameters compatible with M87.

Star clusters are often born as star-cluster systems, which include several stellar clumps. Such star-cluster complexes could have formed from turbulent molecular clouds. Since Gaia Data Release 2 provided us high quality velocity data of individual stars in known star-cluster complexes, we now can compare the velocity structures of the observed star-cluster complexes with simulated ones. We performed a series of N-body simulations for the formation of star-cluster complexes starting from turbulent molecular clouds. We measured the inter-cluster velocity dispersions of our simulated star-cluster complexes and compared them with the Carina region and NGC 2264. We found that the Carina region and NGC 2264 formed from molecular clouds with a mass of ∼4 × 105M⊙ and ∼4 × 104M⊙, respectively. In our simulations, we also found that the maximum cluster mass (Mc,max) in the complex follows ${M_{{\rm{c}},{\rm{max}}}} = 0.{\rm{2}}0M_g^{0.76}$, where Mg is the initial gas mass.

We present a model for hydrodynamic + N-body simulations of star cluster formation and evolution using AMUSE. Our model includes gas dynamics, star formation in regions of dense gas, stellar evolution and a galactic tidal spiral potential, thus incorporating most of the processes that play a role in the evolution of star clusters.

We test our model on initial conditions of two colliding molecular clouds as well as a section of a spiral arm from a previous galaxy simulation.

Three-body interactions of stellar-mass binaries with intermediate mass black holes (IMBHs) in nuclei of globular clusters may produce specific features that may serve as an independent indicator of existence of the IMBHs. By means of direct N-body integrations we follow the dynamical evolution of globular clusters of moderate extension and mass with 50% binary population over a time span of ≍ 0.8 Gyr and compare the cases with and without the primordial binaries as well as with and without the IMBH. We show that (i) presence of the IMBH leads to rapid formation of a density cusp regardless of the initial binary fraction, (ii) binary rich clusters with the IMBH produce high velocity escapers at a rate of ≍ 0.1 Myr−1 and (iii) clusters hosting an IMBH together with high number of binaries form a denser halo of marginally unbound stars than clusters that lack either the IMBH or the binary population.

In this work, we proposed a possible mechanism for the formation S-type planet in close binaries (0.5 au < aB < 3 au). Numerical simulations showed that the maximum capture probability is ∼ 10%, which can be comparable to the tidal capture probability of hot Jupiters in single star systems. The capture probability is related to binary configurations. Furthermore, we find that S-type planets with retrograde orbits can be naturally produced via capture process. These planets on retrograde orbits can help us distinguish in situ formation and post-capture origin for S-type planet in close binaries. The forthcoming missions (PLATO or TESS) will provide the opportunity and feasibility to detect such planets.

We compared the number of lunar craters with diameters greater than 15 km with age less than 1.1 Gyr in the region of the Oceanus Procellarum with the estimates of the number of craters made based on the number of near-Earth objects and on the characteristic times elapsed before collisions of near-Earth objects with the Moon. Our estimates allow the increase of the number of near-Earth objects after a recent catastrophic disruption of a large main-belt asteroid. However, destruction of some old craters and variations in orbital distribution of near-Earth objects with time could allow that the mean number of near-Earth objects during the last billion years could be close to the present value.