Source author record

Yohei Miki

Yohei Miki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2020arXiv

Cusp-to-core transition in low-mass dwarf galaxies induced by dynamical heating of cold dark matter by primordial black holes

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.

preprint2017arXiv

MAGI: many-component galaxy initialiser

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}.

preprint2016arXiv

Collision tomography: Physical properties of possible progenitors of the Andromeda stellar stream

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.

preprint2016arXiv

Formation of the Andromeda Giant Stream: Asymmetric Structure and Disc Progenitor

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.

preprint2016arXiv

GOTHIC: Gravitational oct-tree code accelerated by hierarchical time step controlling

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.

preprint2014arXiv

Hunting A Wandering Supermassive Black Hole in M31 Halo -- Hermitage of Black Hole

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.

preprint2014arXiv

Puzzling Outer-Density Profile of the Dark Matter Halo in the Andromeda 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 $α$ of the outer density distribution of the dark matter halo. The best-fit parameter is $α=-3.7$, which is steeper than the CDM prediction.

preprint2014arXiv

Relics of Galaxy Merging: Observational Predictions for a Wandering Massive Black Hole and Accompanying Star Cluster in the Halo of M31

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.

preprint2013arXiv

Resolving the outer density profile of dark matter halo in Andromeda galaxy

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.

preprint2012arXiv

Astrophysical Particle Simulations on Heterogeneous CPU-GPU Systems

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.