BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Tomoaki Ishiyama

Tomoaki Ishiyama contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.GAastro-ph.SRastro-ph.IMDistributed, Parallel, and Cluster Computing

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 published item(s)

preprint2026arXiv

Rapid emergence of overmassive black holes in the early Universe

The origin of supermassive black holes (SMBHs) remains a long-standing problem in astrophysics. Recent JWST observations reveal an unexpectedly abundant population of overmassive black holes at z>4-6, where the BH masses lie far above local scaling relations and not reproduced by current cosmological models. How such overmassive black holes form and rapidly grow within young galaxies has remained unclear. Here we present fully cosmological radiation-hydrodynamic simulations that, for the first time, self-consistently follow the birth, early growth, and emergent observable signatures of SMBHs in proto-cluster environments. We find that heavy seeds of order $10^6 M_\text{sun}$ naturally form, exceeding typical theoretical expectations by an order of magnitude. These seeds rapidly develop dense, optically thick disks whose strong electron scattering produces broad H$α$ emission comparable to that seen in little red dots (LRDs). Sustained super-Eddington accretion then drives fast growth to $\sim 3 \times 10^7 ~M_\text{sun}$ by $z \sim 8$. These results provide a unified physical scenario in which LRDs correspond to a short-lived, enshrouded phase of heavy-seed formation, naturally evolving into the overmassive quasars detected by JWST and ultimately the progenitors of today's SMBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

The Origin of Spin-Alignment of Dark Matter Subhalos

Subhalo spin is essential for modeling galaxy formation and controlling systematic uncertainties in intrinsic alignment (IA) studies. However, the physical mechanisms governing subhalo spin acquisition within the tidal environments of host halos remain poorly understood. In this work, we investigate the alignment between subhalo and host halo spins using the high-resolution cosmological $N$-body simulation, Shin-Uchuu. We find that the spin alignment between subhalos and host halos becomes increasingly pronounced toward the central regions. Our analysis reveals that subhalos typically acquire spin in the same direction as their orbital angular momentum. Since the orbital angular momentum of most subhalos is aligned with the host halo spin, an overall alignment between subhalo and host spins emerges. When classified by orbital orientation, however, subhalo spins in the inner regions are found to be oriented perpendicularly or anti-parallel to the host spin for polar and retrograde orbits, respectively. These results provide strong evidence that subhalo spins are acquired through torques exerted by the tidal field of the host halo. Furthermore, we demonstrate that the mass ratio and the radial distance from the host center are the primary parameters governing subhalo spin alignment, while the dependence on the accretion redshift is found to be negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Optimizing the Gravitational Tree Algorithm for Many-Core Processors

Gravitational $N$-body simulations calculate numerous interactions between particles. The tree algorithm reduces these calculations by constructing a hierarchical oct-tree structure and approximating gravitational forces on particles. Over the last three decades, the tree algorithm has been extensively used in large-scale simulations, and its parallelization in distributed memory environments has been well studied. However, recent supercomputers are equipped with many CPU cores per node, and optimizations of the tree construction in shared memory environments are becoming crucial. We propose a novel tree construction method in contrast to the conventional top-down approach. It first creates all leaf cells without traversing the tree and then constructs the remaining cells by a bottom-up approach. We evaluated the performance of our novel method on the supercomputer Fugaku and an Intel machine. On a single thread, our method accelerates one of the most time-consuming processes of the conventional tree construction method by a factor of above 3.0 on Fugaku and 2.2 on the Intel machine. Furthermore, as the number of threads increases, our parallel tree construction time reduces considerably. Compared to the conventional sequential tree construction method, we achieve a speedup of over 45 on 48 threads of Fugaku and more than 56 on 112 threads of the Intel machine. In stark contrast to the conventional method, the tree construction with our method no longer constitutes a bottleneck in the tree algorithm, even when using many threads.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Constraining Self-Interacting Dark Matter with Dwarf Spheroidal Galaxies and High-resolution Cosmological $N$-body Simulations

We study the density structures of dark matter subhalos for both cold dark matter and self-interacting dark matter models using high-resolution cosmological $N$-body simulations. We quantify subhalo&#39;s central density at 150 pc from the center of each subhalo at the classical dwarf spheroidal and ultrafaint dwarf scales found in Milky-Way sized halos. By comparing them with observations, we find that the self-interacting scattering cross-section of $σ/m<3\ \rm{cm^{2}g^{-1}}$ is favored. Due to the combination of hosts&#39; tide and self-interactions, the central density of subhalos with small pericenter shows a noticeable difference between the cold and the self-interacting models, indicating that Milky-Way satellites with small pericenter are ideal sites to further constrain the nature of dark matter by future large spectroscopic surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the effect of field variance on the H$α$ luminosity function with the New Numerical Galaxy Catalogue ($ν^2$GC)

We construct a model of H$α$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($ν^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$α$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$α$ luminosity function (LF) at $z=0.4$. The large volume of the $ν^2$GC makes it possible to examine the spatial distribution of HAEs over a region of (411.8 Mpc)$^3$ in the comoving scale. The surface number density of $z=0.4$ HAEs with $L_{\rm Hα} \geq 10^{40}$ erg s$^{-1}$ is 308.9 deg$^{-2}$. We have confirmed that the HAE is a useful tracer for the large-scale structure of the Universe because of their significant overdensity ($>$ 5$σ$) at clusters and the filamentary structures. The H$α$ LFs within a survey area of $\sim$2 deg$^2$ (typical for previous observational studies) show a significant field variance up to $\sim$1 dex. Based on our model, one can estimate the variance on the H$α$ LFs within given survey areas.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-analytic modelling of AGNs: auto-correlation function and halo occupation

The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Abundance and Structure of Subhaloes near the Free Streaming Scale and Their Impact on Indirect Dark Matter Searches

The free streaming motion of dark matter particles imprints a cutoff in the matter power spectrum and set the scale of the smallest dark matter halo. Recent cosmological $N$-body simulations have shown that the central density cusp is much steeper in haloes near the free streaming scale than in more massive haloes. Here, we study the abundance and structure of subhaloes near the free streaming scale at very high redshift using a suite of unprecedentedly large cosmological $N$-body simulations, over a wide range of the host halo mass. The subhalo abundance is suppressed strongly below the free streaming scale, but the ratio between the subhalo mass function in the cutoff and no cutoff simulations is well fitted by a single correction function regardless of the host halo mass and the redshift. In subhaloes, the central slopes are considerably shallower than in field haloes, however, are still steeper than that of the NFW profile. Contrary, the concentrations are significantly larger in subhaloes than haloes and depend on the subhalo mass. We compare two methods to extrapolate the mass-concentration relation of haloes and subhaloes to z=0 and provide a new simple fitting function for subhaloes, based on a suite of large cosmological $N$-body simulations. Finally, we estimate the annihilation boost factor of a Milky-Way sized halo to be between 1.8 and 6.2.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The impact of filamentary accretion of subhaloes on the shape and orientation of haloes

Dark matter haloes are formed through hierarchical mergers of smaller haloes in large-scale cosmic environments, and thus anisotropic subhalo accretion through cosmic filaments have some impacts on halo structures. Recent studies using cosmological simulations have shown that the orientations of haloes correlate with the direction of cosmic filaments, and these correlations significantly depend on the halo mass. Using high-resolution cosmological $N$-body simulations, we quantified the strength of filamentary subhalo accretion for galaxy- and group-sized host haloes ($M_{\rm host}=5\times10^{11-13}M_{\odot}$) by regarding the entry points of subhaloes as filaments and present statistical studies that how the shape and orientation of host haloes at redshift zero correlate with the strength of filamentary subhalo accretion. We confirm previous studies that found the host halo mass dependence of the alignment between orientations of haloes and filaments. We also show that, for the first time, the shape and orientation of haloes weakly correlate with the strength of filamentary subhalo accretion even if the host halo masses are the same. Minor-to-major axis ratios of haloes tend to decrease as their filamentary accretion gets stronger. Haloes with highly anisotropic accretion become more spherical or oblate, while haloes with isotropic accretion become more prolate or triaxial. For haloes with strong filamentary accretion, their major axes are preferentially aligned with the filaments, while their angular momentum vectors tend to be slightly more misaligned.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Discrimination of heavy elements originating from Pop III stars in z = 3 intergalactic medium

We investigate the distribution of metals in the cosmological volume at $z\sim3$, in particular, provided by massive population III (Pop III) stars using a cosmological $N$-body simulation in which a model of Pop III star formation is implemented. Owing to the simulation, we can choose minihaloes where Pop III star formation occurs at $z>10$ and obtain the spatial distribution of the metals at lower-redshifts. To evaluate the amount of heavy elements provided by Pop III stars, we consider metal yield of pair-instability or core-collapse supernovae (SNe) explosions of massive stars. By comparing our results to the Illustris-1 simulation, we find that heavy elements provided by Pop III stars often dominate those from galaxies in low density regions. The median value of the volume averaged metallicity is $Z\sim 10^{-4.5 - -2} Z_{\odot}$ at the regions. Spectroscopic observations with the next generation telescopes are expected to detect the metals imprinted on quasar spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Halo Substructure Boosts to the Signatures of Dark Matter Annihilation

The presence of dark matter substructure will boost the signatures of dark matter annihilation. We review recent progress on estimates of this subhalo boost factor---a ratio of the luminosity from annihilation in the subhalos to that originating the smooth component---based on both numerical $N$-body simulations and semi-analytic modelings. Since subhalos of all the scales, ranging from the Earth mass (as expected, e.g., the supersymmetric neutralino, a prime candidate for cold dark matter) to galaxies or larger, give substantial contribution to the annihilation rate, it is essential to understand subhalo properties over a large dynamic range of more than twenty orders of magnitude in masses. Even though numerical simulations give the most accurate assessment in resolved regimes, extrapolating the subhalo properties down in sub-grid scales comes with great uncertainties---a straightforward extrapolation yields a very large amount of the subhalo boost factor of $\gtrsim$100 for galaxy-size halos. Physically motivated theoretical models based on analytic prescriptions such as the extended Press-Schechter formalism and tidal stripping modeling, which are well tested against the simulation results, predict a more modest boost of order unity for the galaxy-size halos. Giving an accurate assessment of the boost factor is essential for indirect dark matter searches and thus, having models calibrated at large ranges of host masses and redshifts, is strongly urged upon.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of star clusters in a cosmological tidal field

We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gamma-ray Signal from Earth-mass Dark Matter Microhalos

Earth-mass dark matter microhalos with size of ~100 AUs are the first structures formed in the universe, if the dark matter of the Universe are made of neutralino. Here, we report the results of ultra-high-resolution simulations of the formation and evolution of these microhalos. We found that microhalos have the central density cusps of the form $ρ\propto r^{-1.5}$, much steeper than the cusps of larger dark halos. The central regions of these microhalos survive the encounters with stars except in very inner region of the galaxy down to the radius of a few hundreds pcs from the galactic center. The annihilation signals from nearest microhalos are observed as gamma-ray point-sources (radius less than 1&#39;), with unusually large proper motions of ~0.2 degree per year. Their surface brightnesses are ~10% of that of the galactic center. Their S/N ratios might be better if they are far from the galactic plane. Luminosities of subhalos are determined only by their mass, and they are more than one order of magnitude luminous than the estimation by Springel et al. (2008): A boost factor can be larger than 1000. Perturbations to the millisecond pulsars by gravitational attractions of nearby earth-mass microhalos can be detected by the observations of Parkes Pulsar Timing Array (PPTA).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos

We present the results of the &#34;Cosmogrid&#34; cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10^7 Msun. We have analyzed the spherically averaged density profiles of the three most massive halos which are of galaxy group size and contain at least 170 million particles. The slopes of these density profiles become shallower than -1 at the inner most radius. We also find a clear correlation of halo concentration with mass. The mass dependence of the concentration parameter cannot be expressed by a single power law, however a simple model based on the Press-Schechter theory proposed by Navarro et al. gives reasonable agreement with this dependence. The spin parameter does not show a correlation with the halo mass. The probability distribution functions for both concentration and spin are well fitted by the log-normal distribution for halos with the masses larger than ~10^8 Msun. The subhalo abundance depends on the halo mass. Galaxy-sized halos have 50% more subhalos than ~10^{11} Msun halos have.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High Performance Gravitational N-body Simulations on a Planet-wide Distributed Supercomputer

We report on the performance of our cold-dark matter cosmological N-body simulation which was carried out concurrently using supercomputers across the globe. We ran simulations on 60 to 750 cores distributed over a variety of supercomputers in Amsterdam (the Netherlands, Europe), in Tokyo (Japan, Asia), Edinburgh (UK, Europe) and Espoo (Finland, Europe). Regardless the network latency of 0.32 seconds and the communication over 30.000 km of optical network cable we are able to achieve about 87% of the performance compared to an equal number of cores on a single supercomputer. We argue that using widely distributed supercomputers in order to acquire more compute power is technically feasible, and that the largest obstacle is introduced by local scheduling and reservation policies.

0 citations0 reviewsNo code linkedOpen access