The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of $N$-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting a SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50% of Milky Way sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna -- LISA -- which will detect SMBH coalescence in the $10^4-10^7$$M_{\rm \odot}$ mass range.

We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo's orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo's orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 percent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modeling pericentric passages of subhaloes and their associated satellite galaxies.

The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.

Recent observations revealed a mysterious ultra diffuse galaxy, NGC1052-DF2, in the group of a large elliptical galaxy, NGC1052. Compared to expectations from abundance matching models, the dark matter mass contained in NGC1052-DF2 is smaller by a factor of $\sim 400$. We utilize controlled $N$-body simulations of the tidal interaction between NGC1052 and a smaller satellite galaxy, that we suppose as the progenitor of NGC1052-DF2, to test if tidal stripping can explain dark-matter deficiency at such levels. We find that when assuming a tightly bound and quite radial orbit and cored density structure for the dark halo of the satellite, our simulations reproduce well both the mass profile and the effective radius inferred from the observations of NGC1052-DF2. Orbital parameters are in the tail, but still consistent with measurements of their distributions from cosmological simulations. Such strongly dark-matter deficient galaxies, in our scenario, are thus expected to be relatively rare in groups and clusters, and not present in the field.

While feedback from Active Galactic Nuclei (AGN) is an important heating source in the centre of galaxy clusters, it is still unclear how the feedback energy is injected into the intracluster medium (ICM) and what role different numerical approaches play. Here, we compare four hydrodynamical schemes in idealized simulations of a rising bubble inflated by AGN feedback in a hot stratified ICM: (traditional) smoothed particle hydrodynamics (TSPH), a pressure flavour of SPH (PSPH), a meshless finite mass (MFM) scheme, as well as an Eulerian code with adaptive mesh refinement. In the absence of magnetic fields, the bubble is Kelvin-Helmholtz unstable on short enough time scales to dissolve it fully in the ICM, which is captured by MFM and RAMSES simulations, while in the TSPH simulation the bubble survives. When the ICM is turbulent, mixing of the bubble with the ICM is accelerated. This occurs if the numerical scheme can capture the instabilities well. The differences in the evolution of the bubble has a surprisingly small influence on the thermal structure of the ICM. However, in the simulations with MFM and RAMSES the bubble disruption leads to turbulent stirring of the ICM which is suppressed in SPH. In the latter the thermal energy remains trapped in the bubble and is transported to large radii. We discuss if the choice of hydrodynamical schemes can lead to systematic differences in the outcomes of cosmological simulations.

To gain understanding of the complicated, non-linear and numerical processes associated with the tidal evolution of dark matter subhaloes in numerical simulation, we perform a large suite of idealized simulations that follow individual N-body subhaloes in a fixed, analytical host halo potential. By varying both physical and numerical parameters, we investigate under what conditions the subhaloes undergo disruption. We confirm the conclusions from our more analytical assessment in van den Bosch et al. that most disruption is numerical in origin; as long as a subhalo is resolved with sufficient mass and force resolution, a bound remnant survives. This implies that state-of-the-art cosmological simulations still suffer from significant overmerging. We demonstrate that this is mainly due to inadequate force-softening, which causes excessive mass loss and artificial tidal disruption. In addition, we show that subhaloes in N-body simulations are susceptible to a runaway instability triggered by the amplification of discreteness noise in the presence of a tidal field. These two processes conspire to put serious limitations on the reliability of dark matter substructure in state-of-the-art cosmological simulations. We present two criteria that can be used to assess whether individual subhaloes in cosmological simulations are reliable or not, and advocate that subhaloes that satisfy either of these two criteria be discarded from further analysis. We discuss the potential implications of this work for several areas in astrophysics.

Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption. In this paper we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e., artificial). We show that, contrary to naive expectation, subhaloes that experience a tidal shock $\Delta E$ that exceeds the subhalo's binding energy, $|E_{\rm b}|$, do not undergo disruption, even when $\Delta E/|E_{\rm b}|$ is as large as 100. Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy. In addition, we show that tidal heating due to high-speed (impulsive) encounters with other subhaloes (`harassment'), is negligible compared to the tidal effects due to the host halo. Hence, we conclude that, in the absence of baryonic processes, the complete, physical disruption of CDM substructure is extremely rare, and that most disruption in numerical simulations therefore must be artificial. We discuss various processes that have been associated with numerical overmerging, and conclude that inadequate force-softening is the most likely culprit.

Dark matter (DM) haloes forming near the thermal cut-off scale of the density perturbations are unique, since they are the smallest objects and form through monolithic gravitational collapse, while larger haloes contrastingly have experienced mergers. While standard cold dark matter (CDM) simulations readily produce haloes that follow the universal Navarro-Frenk-White (NFW) density profile with an inner slope, $\rho \propto r^{-\alpha}$, with $\alpha=1$, recent simulations have found that when the free-streaming cut-off expected for the CDM model is resolved, the resulting haloes follow nearly power-law density profiles of $\alpha\sim1.5$. In this paper, we study the formation of density cusps in haloes using idealized $N$-body simulations of the collapse of proto-haloes. When the proto-halo profile is initially cored due to particle free-streaming at high redshift, we universally find $\sim r^{-1.5}$ profiles irrespective of the proto-halo profile slope outside the core and large-scale non-spherical perturbations. Quite in contrast, when the proto-halo has a power-law profile, then we obtain profiles compatible with the NFW shape when the density slope of the proto-halo patch is shallower than a critical value, $\alpha_{\rm ini} \sim 0.3$, while the final slope can be steeper for $\alpha_{\rm ini}\ga 0.3$. We further demonstrate that the $r^{-1.5}$ profiles are sensitive to small scale noise, which gradually drives them towards an inner slope of $-1$, where they become resilient to such perturbations. We demonstrate that the $r^{-1.5}$ solutions are in hydrostatic equilibrium, largely consistent with a simple analytic model, and provide arguments that angular momentum appears to determine the inner slope.

Motivated by a recently found interesting property of the dark halo surface density within a radius, $r_{\rm max}$, giving the maximum circular velocity, $V_{\rm max}$, we investigate it for dark halos of the Milky Way's and Andromeda's dwarf satellites based on cosmological simulations. We select and analyze the simulated subhalos associated with Milky Way-sized dark halos and find that the values of their surface densities, $\Sigma_{V_{\rm max}}$, are in good agreement with those for the observed dwarf spheroidal satellites even without employing any fitting procedures. This implies that this surface density would not be largely affected by any baryonic feedbacks and thus universal. Moreover, all subhalos on the small scales of dwarf satellites are expected to obey the relation $\Sigma_{V_{\rm max}}\propto V_{\rm max}$, irrespective of differences in their orbital evolutions, host halo properties, and observed redshifts. Therefore, we find that the universal scaling relation for dark halos on dwarf galaxy mass scales surely exists and provides us important clues to understanding fundamental properties of dark halos. We also investigate orbital and dynamical evolutions of subhalos to understand the origin of this universal dark halo relation and find that most of subhalos evolve generally along the $r_{\rm max}\propto V_{\rm max}$ sequence, even though these subhalos have undergone different histories of mass assembly and tidal stripping. This sequence, therefore, should be the key feature to understand the nature of the universality of $\Sigma_{V_{\rm max}}$.

Primordial dark matter (DM) haloes are the smallest gravitationally bound DM structures from which the first stars, black holes, and galaxies form and grow in the early universe. However, their structures are sensitive to the free streaming scale of DM, which in turn depends on the nature of DM particles. In this work, we test the hypothesis that the slope of the central cusps in primordial DM haloes near the free streaming scale depends on the nature of merging process. By combining and analysing data from a cosmological simulation with the cutoff in the small-scale matter power spectrum as well as a suite of controlled, high-resolution simulations of binary mergers, we find that (1) the primordial DM haloes form preferentially through major mergers in radial orbits; (2) their central DM density profile is more susceptible to a merging process compared to that of galaxy and cluster-size DM haloes; (3) consecutive major mergers drive the central density slope to approach the {\it universal} form characterized by the NFW profile, which is shown to be robust to the impacts of mergers and serves an attractor solution for the density structure of DM haloes. Our work highlights the importance of dynamical processes on the structure formation during the Dark Ages.

We study the dynamical response of extended systems, hosts, to smaller systems, satellites, orbiting around the hosts using extremely high-resolution N-body simulations with up to one billion particles. This situation corresponds to minor mergers which are ubiquitous in the scenario of hierarchical structure formation in the universe. According to Chandrasekhar (1943), satellites create density wakes along the orbit and the wakes cause a deceleration force on satellites, i.e. dynamical friction. This study proposes an analytical model to predict the dynamical response of hosts as reflected in their density distribution and finds not only traditional wakes but also mirror images of over- and underdensities centered on the host. Our controlled N-body simulations with high resolutions verify the predictions of the analytical model. We apply our analytical model to the expected dynamical response of nearby interacting galaxy pairs, the Milky Way - Large Magellanic Cloud system and the M31 - M33 system.

Recent studies found the densities of dark matter (DM) subhaloes which surround nearby dwarf spheroidal galaxies (dSphs) to be significantly lower than those of the most massive subhaloes expected around Milky Way sized galaxies in cosmological simulations, the so called "too-big-to-fail" (TBTF) problem. A caveat of previous work has been that dark substructures were assumed to contain steep density cusps in the center of DM haloes even though the central density structure of DM haloes is still under debate. In this study, we re-examine the TBTF problem for models of DM density structure with cores or shallowed cusps. Our analysis demonstrates that the TBTF problem is alleviated as the logarithmic slope of the central cusp becomes shallower. We find that the TBTF problem is avoided if the central cusps of DM haloes surrounding dSphs are shallower than $r^{-0.6}$.

This study investigates the dynamical response of dark matter (DM) halos to recurrent starbursts in forming less-massive galaxies to solve the core-cusp problem. The gas, which is heated by supernova feedback after a starburst, expands and the star formation then terminates. This expanding gas loses energy by radiative cooling and then falls back toward the galactic center. Subsequently, the starburst is enhanced again. This cycle of expansion and contraction of the interstellar gas leads to a repetitive change in the gravitational potential of the gas. The resonance between DM particles and the density wave excited by the oscillating potential plays a key role in understanding the physical mechanism of the cusp-core transition of DM halos. DM halos effectively gain kinetic energy from the baryon potential through the energy transfer driven by the resonance between the particles and density waves. We determine that the critical condition for the cusp-core transition is such that the oscillation period of the gas potential is approximately the same as the local dynamical time of DM halos. We present the resultant core radius of a DM halo after the cusp-core transition induced by the resonance by using the conventional mass density profile predicted by the cold dark matter models. Moreover, we verify the analytical model by using $N$-body simulations, and the results validate the resonance model.

Observations have revealed interesting universal properties of dark matter (DM) halos especially around low-mass galaxies. Strigari et al. (2008) showed that DM halos have common enclosed masses within 300pc (Strigari relation). Kormendy & Freeman (2004) reported DM halos having almost identical central surface densities (the $\mu_{\rm 0D}$ relation). In addition, there exists a core-cusp problem, a discrepancy of the central density distribution between simulated halos and observations. We investigate whether a scenario where cuspy halos transform into cores by some dynamical processes can also explain their universal structural properties. It is shown that a cusp-to-core transformation model naturally reproduces the $\mu_{\rm 0D}$ relation and that Strigari relation follows from the $\mu_{\rm 0D}$ relation for dwarf galaxies. We also show that the central densities of cored dark halos provide valuable information about their formation redshifts.

A heterogeneous CPU-GPU node is getting popular in HPC clusters. We need to rethink algorithms and optimization techniques for such system depending on the relative performance of CPU vs. GPU. In this paper, we report a performance optimized particle simulation code "OTOO", that is based on the octree method, for heterogenous systems. Main applications of OTOO are astrophysical simulations such as N-body models and the evolution of a violent merger of stars. We propose optimal task split between CPU and GPU where GPU is only used to compute the calculation of the particle force. Also, we describe optimization techniques such as control of the force accuracy, vectorized tree walk, and work partitioning among multiple GPUs. We used OTOO for modeling a merger of two white dwarf stars and found that OTOO is powerful and practical to simulate the fate of the process.

The core-cusp problem remains as one of the unsolved discrepancies between observations and theories predicted by the standard paradigm of cold dark matter (CDM) cosmology. To solve this problem, we perform N-body simulations to study the nonlinear response of CDM halos to the variance of the gravitational potential induced by gas removal from galaxy centers. In this study, we focus on the timescale of the gas ejection, which is strongly correlated with stellar activities, and demonstrate that it is one of the key factors in determining the dynamical response of CDM halos. The results of simulations show that the power-low index of the mass-density profile of the dark matter halo correlated with the timescale of the mass loss, and it is flatter when the mass loss occurs over a short time than when it occurs over a long time. However, it is still larger than typical observational values; in other words, the central cusp remains for any mass loss model in the simulations. Moreover, for the slow mass-loss case, the final density profile of the dark matter halo recovers the universal density profiles predicted by the CDM cosmology. Therefore, mass loss driven by stellar feedback may not be an effective mechanism to flatten the central cusp.