We performed a series of high-resolution $N$-body simulations to examine whether dark matter candidates in the form of primordial black holes (PBHs) can solve the cusp-core problem in low-mass dwarf galaxies. If some fraction of the dark matter in low-mass dwarf galaxies consists of PBHs and the rest is cold dark matter, dynamical heating of the cold dark matter by the PBHs induces a cusp-to-core transition in the total dark matter profile. The mechanism works for PBHs in the 25-100 M$_{\sun}$ mass window, consistent with the LIGO detections, but requires a lower limit on the PBH mass fraction of 1$\%$ of the total dwarf galaxy dark matter content. The cusp-to-core transition time-scale is between 1 and 8 Gyr. This time-scale is also a constant multiple of the relaxation time between cold dark matter particles and PBHs, which depends on the mass, the mass fraction and the scale radius of the initial density profile of PBHs. We conclude that dark matter cores occur naturally in halos comprised of cold dark matter and PBHs, without the need to invoke baryonic processes.

In this study, the gravitational octree code originally optimized for the Fermi, Kepler, and Maxwell GPU architectures is adapted to the Volta architecture. The Volta architecture introduces independent thread scheduling requiring either the insertion of the explicit synchronizations at appropriate locations or the enforcement of the same implicit synchronizations as do the Pascal or earlier architectures by specifying \texttt{-gencode arch=compute\_60,code=sm\_70}. The performance measurements on Tesla V100, the current flagship GPU by NVIDIA, revealed that the $N$-body simulations of the Andromeda galaxy model with $2^{23} = 8388608$ particles took $3.8 \times 10^{-2}$~s or $3.3 \times 10^{-2}$~s per step for each case. Tesla V100 achieves a 1.4 to 2.2-fold acceleration in comparison with Tesla P100, the flagship GPU in the previous generation. The observed speed-up of 2.2 is greater than 1.5, which is the ratio of the theoretical peak performance of the two GPUs. The independence of the units for integer operations from those for floating-point number operations enables the overlapped execution of integer and floating-point number operations. It hides the execution time of the integer operations leading to the speed-up rate above the theoretical peak performance ratio. Tesla V100 can execute $N$-body simulation with up to $25 \times 2^{20} = 26214400$ particles, and it took $2.0 \times 10^{-1}$~s per step. It corresponds to $3.5$~TFlop/s, which is 22\% of the single-precision theoretical peak performance.

Providing initial conditions is an essential procedure for numerical simulations of galaxies. The initial conditions for idealised individual galaxies in $N$-body simulations should resemble observed galaxies and be dynamically stable for time scales much longer than their characteristic dynamical times. However, generating a galaxy model ab initio as a system in dynamical equilibrium is a difficult task, since a galaxy contains several components, including a bulge, disc, and halo. Moreover, it is desirable that the initial-condition generator be fast and easy to use. We have now developed an initial-condition generator for galactic $N$-body simulations that satisfies these requirements. The developed generator adopts a distribution-function-based method, and it supports various kinds of density models, including custom-tabulated inputs and the presence of more than one disc. We tested the dynamical stability of systems generated by our code, representing early- and late-type galaxies, with $N=$~2,097,152 and 8,388,608 particles, respectively, and we found that the model galaxies maintain their initial distributions for at least 1~Gyr. The execution times required to generate the two models were $8.5$ and $221.7$ seconds, respectively, which is negligible compared to typical execution times for $N$-body simulations. The code is provided as open-source software and is publicly and freely available at \url{https://bitbucket.org/ymiki/magi}.

We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream.

We examine the nature, possible orbits and physical properties of the progenitor of the North-western stellar stream (NWS) in the halo of the Andromeda galaxy (M31). The progenitor is assumed to be an accreting dwarf galaxy with globular clusters (GCs). It is, in general, difficult to determine the progenitor's orbit precisely because of many necessary parameters. Recently, Veljanoski et al. 2014 reported five GCs whose positions and radial velocities suggest an association with the stream. We use this data to constrain the orbital motions of the progenitor using test-particle simulations. Our simulations split the orbit solutions into two branches according to whether the stream ends up in the foreground or in the background of M31. Upcoming observations that will determine the distance to the NWS will be able to reject one of the two branches. In either case, the solutions require that the pericentric radius of any possible orbit be over 2 kpc. We estimate the efficiency of the tidal disruption and confirm the consistency with the assumption for the progenitor being a dwarf galaxy. The progenitor requires the mass $\ga 2\times10^6 M_{\sun}$ and half-light radius $\ga 30$ pc. In addition, $N$-body simulations successfully reproduce the basic observed features of the NWS and the GCs' line-of-sight velocities.

The tree method is a widely implemented algorithm for collisionless $N$-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate $N$-body simulations; however, it is infrequently implemented and its potential remains untested in GPU implementations. We have developed a Gravitational Oct-Tree code accelerated by HIerarchical time step Controlling named \texttt{GOTHIC}, which adopts both the tree method and the hierarchical time step. The code adopts some adaptive optimizations by monitoring the execution time of each function on-the-fly and minimizes the time-to-solution by balancing the measured time of multiple functions. Results of performance measurements with realistic particle distribution performed on NVIDIA Tesla M2090, K20X, and GeForce GTX TITAN X, which are representative GPUs of the Fermi, Kepler, and Maxwell generation of GPUs, show that the hierarchical time step achieves a speedup by a factor of around 3--5 times compared to the shared time step. The measured elapsed time per step of \texttt{GOTHIC} is 0.30~s or 0.44~s on GTX TITAN X when the particle distribution represents the Andromeda galaxy or the NFW sphere, respectively, with $2^{24} =$~16,777,216 particles. The averaged performance of the code corresponds to 10--30\% of the theoretical single precision peak performance of the GPU.

We focus on the evidence of a past minor merger discovered in the halo of the Andromeda galaxy (M31). Previous N-body studies have enjoyed moderate success in producing the observed giant stellar stream (GSS) and stellar shells in M31's halo. The observed distribution of stars in the halo of M31 shows an asymmetric surface brightness profile across the GSS; however, the effect of the morphology of the progenitor galaxy on the internal structure of the GSS requires further investigation in theoretical studies. To investigate the physical connection between the characteristic surface brightness in the GSS and the morphology of the progenitor dwarf galaxy, we systematically vary the thickness, rotation velocity and initial inclination of the disc dwarf galaxy in N-body simulations. The formation of the observed structures appears to be dominated by the progenitor's rotation. Besides reproducing the observed GSS and two shells in detail, we predict additional structures for further observations. We predict the detectability of the progenitor's stellar core in the phase-space density distribution, azimuthal metallicity gradient of the western shell-like structure and an additional extended shell in the north-western direction that may constrain the properties of the progenitor galaxy.

To unveil a progenitor of the Andromeda Giant Stellar Stream, we investigate the interaction between an accreting satellite galaxy and the Andromeda Galaxy using an $N$-body simulation. A comprehensive parameter study with 247 models is performed by varying size and mass distribution of the progenitor dwarf galaxy. We show that the binding energy of the progenitor is the crucial parameter in reproducing the Andromeda Giant Stellar Stream and the shell-like structures surrounding the Andromeda Galaxy. As a result of the simulations, the progenitor must satisfy a simple scaling relation between the core radius, the total mass and the tidal radius. Using this relation, we successfully constrain the physical properties of the progenitors to have mass ranging from $5\times10^8 M_\odot$ to $5\times10^9 M_\odot$ and central surface density around $10^3\, M_\odot\, \mathrm{pc}^{-2}$. A detailed comparison between our result and the nearby observed galaxies indicates that possible progenitors of the Andromeda Giant Stellar Stream include a dwarf elliptical galaxy, a dwarf irregular galaxy, and a small spiral galaxy.

The cold dark matter (CDM) cosmology, which is the standard theory of the structure formation in the universe, predicts that the outer density profile of dark matter halos decreases with the cube of distance from the center. However, so far not much effort has examined this hypothesis. In the halo of the Andromeda galaxy (M31), large-scale stellar structures detected by the recent observations provide a potentially suitable window to investigate the mass--density distribution of the dark matter halo. We explore the density structure of the dark matter halo in M31 using an N-body simulation of the interaction between an accreting satellite galaxy and M31. To reproduce the Andromeda Giant Southern Stream and the stellar shells at the east and west sides of M31, we find the sufficient condition for the power-law index $\alpha$ of the outer density distribution of the dark matter halo. The best-fit parameter is $\alpha=-3.7$, which is steeper than the CDM prediction.

Galaxies and massive black holes (BHs) presumably grow via galactic merging events and subsequent BH coalescence. As a case study, we investigate the merging event between the Andromeda galaxy (M31) and a satellite galaxy. We compute the expected observational appearance of the massive BH that was at the center of the satellite galaxy prior to the merger, and is currently wandering in the M31 halo. We demonstrate that a radiatively inefficient accretion flow with a bolometric luminosity of a few tens of solar luminosities develops when Hoyle-Lyttleton accretion onto the BH is assumed. We compute the associated broadband spectrum and show that the radio band (observable with EVLA, ALMA and SKA) is the best frequency range to detect the emission. We also evaluate the mass and the luminosity of the stars bound by the wandering BH and find that such a star cluster is sufficiently luminous that it could correspond to one of the star clusters found by the PAndAS survey. The discovery of a relic massive BH wandering in a galactic halo will provide a direct means to investigate in detail the coevolution of galaxies and BHs. It also means a new population of BHs (off-center massive BHs), and offers targets for clean BH imaging that avoids strong interstellar scattering in the center of galaxies.

In the hierarchical structure formation scenario, galaxies enlarge through multiple merging events with less massive galaxies. In addition, the Magorrian relation indicates that almost all galaxies are occupied by a central supermassive black hole (SMBH) of mass $10^{-3}$ of its spheroidal component. Consequently, SMBHs are expected to wander in the halos of their host galaxies following a galaxy collision, although evidence of this activity is currently lacking. We investigate a current plausible location of an SMBH wandering in the halo of the Andromeda galaxy (M31). According to theoretical studies of $N$-body simulations, some of the many substructures in the M31 halo are remnants of a minor merger occurring about 1 Gyr ago. First, to evaluate the possible parameter space of the infalling orbit of the progenitor, we perform numerous parameter studies using a Graphics Processing Unit (GPU) cluster. To reduce uncertainties in the predicted position of the expected SMBH, we then calculate the time evolution of the SMBH in the progenitor dwarf galaxy from $N$-body simulations using the plausible parameter sets. Our results show that the SMBH lies within the halo ($\sim$20--50 kpc from the M31 center), closer to the Milky Way than the M31 disk. Furthermore, the predicted current positions of the SMBH were restricted to an observational field of $0\degr.6 \times 0\degr.7$ in the northeast region of the M31 halo. We also discuss the origin of the infalling orbit of the satellite galaxy and its relationships with the recently discovered vast thin disk plane of satellite galaxies around M31.

Large-scale faint structure detected by the recent observations in the halo of the Andromeda galaxy (M31) provides an attractive window to explore the structure of outer cold dark matter (CDM) halo in M31. Using an N-body simulation of the interaction between an accreting satellite galaxy and M31, we investigate the mass density profile of the CDM halo. We find the sufficient condition of the outer density profile of CDM halo in M31 to reproduce the Andromeda giant stream and the shells at the east and west sides of M31. The result indicates that the density profile of the outer dark matter halo of M31 is a steeper than the prediction of the theory of the structure formation based on the CDM model.

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.