Showing 1–50 of 151 results for author: Makino, J

Statistical Study of Planetesimal-Driven Migration

Authors: Tenri Jinno, Takayuki R. Saitoh, Junichiro Makino

Abstract: Recent exoplanet observations have revealed a diversity of exoplanetary systems, which suggests the ubiquity of radial planetary migration. One powerful known mechanism of planetary migration is planetesimal-driven migration (PDM), which makes planets undergo significant migration through gravitational scattering with planetesimals. Here, we present the results of our high-resolution self-consiste… ▽ More Recent exoplanet observations have revealed a diversity of exoplanetary systems, which suggests the ubiquity of radial planetary migration. One powerful known mechanism of planetary migration is planetesimal-driven migration (PDM), which makes planets undergo significant migration through gravitational scattering with planetesimals. Here, we present the results of our high-resolution self-consistent $N$-body simulations of PDM, in which gravitational interactions among planetesimals, the gas drag, and Type-I migration are all taken into account. Our results show that even small protoplanets can actively migrate through PDM. Moreover, a fair fraction of them migrate outward. This outward migration can give a solution for the ''planet migration problem'' caused by Type-I migration and explain the origin of Jovian planets. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: 25 pages, 13 figures, submitted

The jet and resolved features of the central supermassive black hole of M 87 observed with EHT in 2017 -- Comparison with the GMVA 86 GHz results

Authors: Makoto Miyoshi, Yoshiaki Kato, Junichiro Makino, Masato Tsuboi

Abstract: M 87 is the best target for studying black hole accretion and jet formation. Reanalysis of the EHT public data at 230 GHz shows a core-knots structure at the center and jet features (Miyoshi et al. 2022a). We here compare this with the new results of GMVA at 86 GHz showing a spatially resolved central core (Lu et al. 2023a). There are similarities and differences between the two. At 86 GHz,"two br… ▽ More M 87 is the best target for studying black hole accretion and jet formation. Reanalysis of the EHT public data at 230 GHz shows a core-knots structure at the center and jet features (Miyoshi et al. 2022a). We here compare this with the new results of GMVA at 86 GHz showing a spatially resolved central core (Lu et al. 2023a). There are similarities and differences between the two. At 86 GHz,"two bright regions" are seen on the ring in the core."Core-Knot-Westknot", triple structure in the 230 GHz image shows apparent appearance of two peaks similar to the "two bright regions" when convolved with the GMVA beam. This similarity suggests that both frequencies reveal the same objects in the core area. Protrusions are observed on both the south and north sides of the core at both frequencies, becoming prominent and wing-like at 230 GHz. The 86 GHz image shows a triple ridge jet structure, while the 230 GHz image shows only a bright central ridge with two roots. Both frequencies show a shade between the core and the central ridge. To detect the faint features from the EHT2017 data, we found that the use of all baseline data is essential. Using all including the ultrashort baseline data, revealed the jet and faint structures. Without the ultrashort baselines, these structures were not detectable. The lack of detection of any faint structures other than the ring in the M 87 data by the EHTC is presumably due to the exclusion of ultrashort baselines from their analysis. △ Less

Submitted 31 January, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 17 pages, 4 figures, ApJL, Accepted 27, January, 2024

Surrogate Modeling for Computationally Expensive Simulations of Supernovae in High-Resolution Galaxy Simulations

Authors: Keiya Hirashima, Kana Moriwaki, Michiko S. Fujii, Yutaka Hirai, Takayuki R. Saitoh, Junichiro Makino, Shirley Ho

Abstract: Some stars are known to explode at the end of their lives, called supernovae (SNe). The substantial amount of matter and energy that SNe release provides significant feedback to star formation and gas dynamics in a galaxy. SNe release a substantial amount of matter and energy to the interstellar medium, resulting in significant feedback to star formation and gas dynamics in a galaxy. While such fe… ▽ More Some stars are known to explode at the end of their lives, called supernovae (SNe). The substantial amount of matter and energy that SNe release provides significant feedback to star formation and gas dynamics in a galaxy. SNe release a substantial amount of matter and energy to the interstellar medium, resulting in significant feedback to star formation and gas dynamics in a galaxy. While such feedback has a crucial role in galaxy formation and evolution, in simulations of galaxy formation, it has only been implemented using simple {\it sub-grid models} instead of numerically solving the evolution of gas elements around SNe in detail due to a lack of resolution. We develop a method combining machine learning and Gibbs sampling to predict how a supernova (SN) affects the surrounding gas. The fidelity of our model in the thermal energy and momentum distribution outperforms the low-resolution SN simulations. Our method can replace the SN sub-grid models and help properly simulate un-resolved SN feedback in galaxy formation simulations. We find that employing our new approach reduces the necessary computational cost to $\sim$ 1 percent compared to directly resolving SN feedback. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 11 pages, 9 figures, Accepted for the NeurIPS 2023 AI4Science Workshop

$N$-body simulation of planetary formation through pebble accretion in a radially structured protoplanetary disk

Authors: Tenri Jinno, Takayuki R. Saitoh, Yota Ishigaki, Junichiro Makino

Abstract: In the conventional theory of planet formation, it is assumed that protoplanetary disks are axisymmetric and have a smooth radial profile. However, recent radio observations of protoplanetary disks have revealed that many of them have complex radial structures. In this study, we perform a series of $N$-body simulations to investigate how planets are formed in protoplanetary disks with radial struc… ▽ More In the conventional theory of planet formation, it is assumed that protoplanetary disks are axisymmetric and have a smooth radial profile. However, recent radio observations of protoplanetary disks have revealed that many of them have complex radial structures. In this study, we perform a series of $N$-body simulations to investigate how planets are formed in protoplanetary disks with radial structures. For this purpose, we consider the effect of continuous pebble accretion onto the discontinuity boundary within the terrestrial planet-forming region ($\sim0.6$ AU). We found that protoplanets grow efficiently at the discontinuity boundary, reaching the Earth mass within $\sim10^4$ years. We confirmed that giant collisions of protoplanets occur universally in our model. Moreover, we found that multiple planet-sized bodies form at regular intervals in the vicinity of the discontinuity boundary. These results indicate the possibility of the formation of solar system-like planetary systems in radially structured protoplanetary disks. △ Less

Submitted 10 August, 2023; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: 19 pages, 19 figures, accepted for publication in PASJ (Publications of the Astronomical Society of Japan). The paper is the same as v1, but comments on the paper status added

3D-Spatiotemporal Forecasting the Expansion of Supernova Shells Using Deep Learning toward High-Resolution Galaxy Simulations

Authors: Keiya Hirashima, Kana Moriwaki, Michiko S. Fujii, Yutaka Hirai, Takayuki R. Saitoh, Junichiro Makino

Abstract: Supernova (SN) plays an important role in galaxy formation and evolution. In high-resolution galaxy simulations using massively parallel computing, short integration timesteps for SNe are serious bottlenecks. This is an urgent issue that needs to be resolved for future higher-resolution galaxy simulations. One possible solution would be to use the Hamiltonian splitting method, in which regions req… ▽ More Supernova (SN) plays an important role in galaxy formation and evolution. In high-resolution galaxy simulations using massively parallel computing, short integration timesteps for SNe are serious bottlenecks. This is an urgent issue that needs to be resolved for future higher-resolution galaxy simulations. One possible solution would be to use the Hamiltonian splitting method, in which regions requiring short timesteps are integrated separately from the entire system. To apply this method to the particles affected by SNe in a smoothed-particle hydrodynamics simulation, we need to detect the shape of the shell on and within which such SN-affected particles reside during the subsequent global step in advance. In this paper, we develop a deep learning model, 3D-MIM, to predict a shell expansion after a SN explosion. Trained on turbulent cloud simulations with particle mass $m_{\rm gas}=1$M$_\odot$, the model accurately reproduces the anisotropic shell shape, where densities decrease by over 10 per cent by the explosion. We also demonstrate that the model properly predicts the shell radius in the uniform medium beyond the training dataset of inhomogeneous turbulent clouds. We conclude that our model enables the forecast of the shell and its interior where SN-affected particles will be present. △ Less

Submitted 17 September, 2023; v1 submitted 31 January, 2023; originally announced February 2023.

Comments: 14 pages, 14 figures, 3 tables, accepted for MNRAS

Comments on "Imaging Reanalyses of EHT Data"

Authors: Makoto Miyoshi, Yoshiaki Kato, Jun Makino

Abstract: On June 14, 2022, the EHT collaboration (hereafter EHTC) made the web page (https://eventhorizontelescope.org/blog/imaging-reanalyses-eht-data) with the title "Imaging Reanalyses of EHT Data," in which they made comments on our recent Miyoshi et al .2022 published in the Astrophysical Journal. We investigated the EHTC comments and found that all of the five points raised by the EHTC are subjecti… ▽ More On June 14, 2022, the EHT collaboration (hereafter EHTC) made the web page (https://eventhorizontelescope.org/blog/imaging-reanalyses-eht-data) with the title "Imaging Reanalyses of EHT Data," in which they made comments on our recent Miyoshi et al .2022 published in the Astrophysical Journal. We investigated the EHTC comments and found that all of the five points raised by the EHTC are subjective and unsubstantiated claims. Thus they do not prove the correctness of the result of EHTC. Sincerely we hope that the EHTC will publish, not a collection of unsubstantiated claims, but a discussion based on scientific arguments. △ Less

Submitted 1 August, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: Comments on the EHTC description on their web page

Origin of highly $r$-process-enhanced stars in a cosmological zoom-in simulation of a Milky Way-like galaxy

Authors: Yutaka Hirai, Timothy C. Beers, Masashi Chiba, Wako Aoki, Derek Shank, Takayuki R. Saitoh, Takashi Okamoto, Junichiro Makino

Abstract: The $r$-process-enhanced (RPE) stars provide fossil records of the assembly history of the Milky Way and the nucleosynthesis of the heaviest elements. Observations by the $R$-Process Alliance (RPA) and others have confirmed that many RPE stars are associated with chemo-dynamically tagged groups, which likely came from accreted dwarf galaxies of the Milky Way (MW). However, we do not know how RPE s… ▽ More The $r$-process-enhanced (RPE) stars provide fossil records of the assembly history of the Milky Way and the nucleosynthesis of the heaviest elements. Observations by the $R$-Process Alliance (RPA) and others have confirmed that many RPE stars are associated with chemo-dynamically tagged groups, which likely came from accreted dwarf galaxies of the Milky Way (MW). However, we do not know how RPE stars are formed. Here, we present the result of a cosmological zoom-in simulation of an MW-like galaxy with $r$-process enrichment, performed with the highest resolution in both time and mass. Thanks to this advancement, unlike previous simulations, we find that most highly RPE ($r$-II; [Eu/Fe] $> +0.7$) stars are formed in low-mass dwarf galaxies that have been enriched in $r$-process elements for [Fe/H] $\,<-2.5$, while those with higher metallicity are formed in situ, in locally enhanced gas clumps that were not necessarily members of dwarf galaxies. This result suggests that low-mass accreted dwarf galaxies are the main formation site of $r$-II stars with [Fe/H] $\,<-2.5$. We also find that most low-metallicity $r$-II stars exhibit halo-like kinematics. Some $r$-II stars formed in the same halo show low dispersions in [Fe/H] and somewhat larger dispersions of [Eu/Fe], similar to the observations. The fraction of simulated $r$-II stars is commensurate with observations from the RPA, and the distribution of the predicted [Eu/Fe] for halo $r$-II stars matches that observed. These results demonstrate that RPE stars can be valuable probes of the accretion of dwarf galaxies in the early stages of their formation. △ Less

Submitted 13 December, 2022; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: 20 pages, 15 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 517, Issue 4, December 2022, Pages 4856-4874

The jet and resolved features of the central supermassive black hole of M 87 observed with EHT

Authors: Makoto Miyoshi, Yoshiaki Kato, Junichiro Makino

Abstract: We report our independent image reconstruction of the M 87 from the public data of the Event Horizon Telescope Collaborators (EHTC). Our result is different from the image published by the EHTC. Our analysis shows that (a) the structure at 230 GHz is consistent with those of lower frequency VLBI observations, (b) the jet structure is evident at 230 GHz extending from the core to a few mas, though… ▽ More We report our independent image reconstruction of the M 87 from the public data of the Event Horizon Telescope Collaborators (EHTC). Our result is different from the image published by the EHTC. Our analysis shows that (a) the structure at 230 GHz is consistent with those of lower frequency VLBI observations, (b) the jet structure is evident at 230 GHz extending from the core to a few mas, though the intensity rapidly decreases along the axis, and (c) the unresolved core is resolved into bright three features presumably showing an initial jet with a wide opening angle of about 70 deg. The ring-like structures of the EHTC can be created not only from the public data, but also from the simulated data of a point image. Also, the rings are very sensitive to the FOV size. The u-v coverage of EHT lack about 40 micro-asec fringe spacings. Combining with a very narrow FOV, it created the 40 micro-asec ring structure. We conclude that the absence of the jet and the presence of the ring in the EHTC result are both artifacts owing to the narrow FOV setting and the u-v data sampling bias effect of the EHT array. Because the EHTC's simulations only take into account the reproduction of the input image models, and not those of the input noise models, their optimal parameters can enhance the effects of sampling bias and produce artifacts such as the 40 micro-asec ring structure, rather than reproducing the correct image. △ Less

Submitted 19 May, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 68 pages, 34 figures, ApJ accepted 5th May 2022

Particle-particle Particle-tree Code for Planetary System Formation with Individual Cut-off Method: GPLUM

Authors: Yota Ishigaki, Junko Kominami, Junichiro Makino, Masaki Fujimoto, Masaki Iwasawa

Abstract: In a standard theory of the formation of the planets in our Solar System, terrestrial planets and cores of gas giants are formed through accretion of kilometer-sized objects (planetesimals) in a protoplanetary disk. Gravitational $N$-body simulations of a disk system made up of numerous planetesimals are the most direct way to study the accretion process. However, the use of $N$-body simulations h… ▽ More In a standard theory of the formation of the planets in our Solar System, terrestrial planets and cores of gas giants are formed through accretion of kilometer-sized objects (planetesimals) in a protoplanetary disk. Gravitational $N$-body simulations of a disk system made up of numerous planetesimals are the most direct way to study the accretion process. However, the use of $N$-body simulations has been limited to idealized models (e.g. perfect accretion) and/or narrow spatial ranges in the radial direction, due to the limited number of simulation runs and particles available. We have developed new $N$-body simulation code equipped with a particle-particle particle-tree (${\rm P^3T}$) scheme for studying the planetary system formation process: GPLUM. For each particle, GPLUM uses the fourth-order Hermite scheme to calculate gravitational interactions with particles within cut-off radii and the Barnes-Hut tree scheme for particles outside the cut-off radii. In existing implementations, ${\rm P^3T}$ schemes use the same cut-off radius for all particles, making a simulation become slower when the mass range of the planetesimal population becomes wider. We have solved this problem by allowing each particle to have an appropriate cut-off radius depending on its mass, its distance from the central star, and the local velocity dispersion of planetesimals. In addition to achieving a significant speed-up, we have also improved the scalability of the code to reach a good strong-scaling performance up to 1024 cores in the case of $N=10^6$. GPLUM is freely available from https://github.com/YotaIshigaki/GPLUM with MIT license. △ Less

Submitted 16 June, 2021; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 17 pages, 16 figures

Journal ref: Publications of the Astronomical Society of Japan, Volume 73, Issue 3, June 2021, Pages 660-676

PeTar: a high-performance N-body code for modeling massive collisional stellar systems

Authors: Long Wang, Masaki Iwasawa, Keigo Nitadori, Junichiro Makino

Abstract: The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been performed by using the NBODY6++GPU code. Such models took half a year computational time on a GPU based super-computer. In this work, we develop a new N-body code,… ▽ More The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been performed by using the NBODY6++GPU code. Such models took half a year computational time on a GPU based super-computer. In this work, we develop a new N-body code, PeTar, by combining the methods of Barnes-Hut tree, Hermite integrator and slow-down algorithmic regularization (SDAR). The code can accurately handle an arbitrary fraction of multiple systems (e.g. binaries, triples) while keeping a high performance by using the hybrid parallelization methods with MPI, OpenMP, SIMD instructions and GPU. A few benchmarks indicate that PeTar and NBODY6++GPU have a very good agreement on the long-term evolution of the global structure, binary orbits and escapers. On a highly configured GPU desktop computer, the performance of a million-body simulation with all stars in binaries by using PeTar is 11 times faster than that of NBODY6++GPU. Moreover, on the Cray XC50 supercomputer, PeTar well scales when number of cores increase. The ten million-body problem, which covers the region of ultra compact dwarfs and nuclearstar clusters, becomes possible to be solved. △ Less

Submitted 27 July, 2020; v1 submitted 30 June, 2020; originally announced June 2020.

Comments: 20 pages, 17 figures, accepted for MNRAS

Journal ref: MNRAS 497 (2020) 536-555

Further comments on "Individual external dose monitoring of all citizens of Date City by passive dosimeter 5 to 51 months after the Fukushima NPP accident (series): 1." : Inconsistencies in Table 1 2014 Q3 and Figure 4f

Authors: Shin-ichi Kurokawa, Yutaka Hamaoka, Kyo Kageura, Jun Makino, Masaki Oshikawa, Yoh Tanimoto

Abstract: We point out serious inconsistencies of the first paper of the series, written by Makoto Miyazaki and Ryugo Hayano, which discusses the correlation between the personal doses of the citizens of Date City measured by glass badges with the ambient dose rates measured by six airborne surveys. The last of the six airborne survey was made in the period of 2014 Q3 (from October 2014 to December 2014). T… ▽ More We point out serious inconsistencies of the first paper of the series, written by Makoto Miyazaki and Ryugo Hayano, which discusses the correlation between the personal doses of the citizens of Date City measured by glass badges with the ambient dose rates measured by six airborne surveys. The last of the six airborne survey was made in the period of 2014 Q3 (from October 2014 to December 2014). The real number of participants of the period is about 14,500; however, in Table 1 2014 Q3 it is written that the number of participants is 21,080 and in Fig. 4f 21,052. We conclude that the analysis of the paper with respect to Table 1 2014 Q3 and Fig. 4f are done without using real correct data and we cannot obtain any meaningful information from the table and figure. Since the period 2014 Q3 is also included in Fig. 5 of the second paper of the series, it is quite possible that Fig. 5 of the second paper is made on the basis of, at least partially, false data and is not reliable. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 7 pages 3 figures

A slow-down time-transformed symplectic integrator for solving the few-body problem

Authors: Long Wang, Keigo Nitadori, Junichiro Makino

Abstract: An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which combines the time-transformed explicit symplectic integrator (Preto & Tremaine 1999; Mikkola & Tanikawa 1999) and the slow-down method (Mikkola & Aarseth 1996). Th… ▽ More An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which combines the time-transformed explicit symplectic integrator (Preto & Tremaine 1999; Mikkola & Tanikawa 1999) and the slow-down method (Mikkola & Aarseth 1996). The former conserves the Hamiltonian and the angular momentum for a long-term evolution, while the latter significantly reduces the computational cost for a weakly perturbed binary. In this work, the Hamilton equations of this algorithm are analyzed in detail. We mathematically and numerically show that it can correctly reproduce the secular evolution like the orbit averaged method and also well conserve the angular momentum. For a weakly perturbed binary, the method is possible to provide a few order of magnitude faster performance than the classical algorithm. A publicly available code written in the c++ language, SDAR, is available on GitHub (https://github.com/lwang-astro/SDAR). It can be used either as a stand alone tool or a library to be plugged in other $N$-body codes. The high precision of the floating point to 62 digits is also supported. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 14 pages, 13 figures, accepted to MNRAS

Comment on "Individual external dose monitoring of all citizens of Date City by passive dosimeter 5 to 51 months after the Fukushima NPP accident (series): 1."

Authors: Masaki Oshikawa, Yutaka Hamaoka, Kyo Kageura, Shin-ichi Kurokawa, Jun Makino, Yoh Tanimoto

Abstract: We point out numerous inconsistencies and inappropriate statements in M. Miyazaki and H. Hayano, Journal of Radiological Protection 37, 1 (2016). We point out numerous inconsistencies and inappropriate statements in M. Miyazaki and H. Hayano, Journal of Radiological Protection 37, 1 (2016). △ Less

Submitted 31 January, 2020; originally announced January 2020.

Comments: 6 pages, 1 figure

Are long-term $N$-body simulations reliable?

Authors: David M. Hernandez, Sam Hadden, Junichiro Makino

Abstract: $N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about $200$ e-folding times. Using nearby initial conditions, we can cons… ▽ More $N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about $200$ e-folding times. Using nearby initial conditions, we can construct statistics of the long-term phase-space structure and compare to rough estimates of resonant widths of the system. We compared statistics for a wide range of numerical methods, including a Runge--Kutta method, Wisdom--Holman method, symplectic corrector methods, and a method by Laskar and Robutel. "Improving" an integrator did not increase the phase space accuracy, but simply increasing the number of initial conditions did. In fact, the statistics of a higher order symplectic corrector method were inconsistent with the other methods in one test. △ Less

Submitted 6 February, 2020; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 14 pages, 13 figures, 4 tables. Matches accepted MNRAS version

Accelerated FDPS --- Algorithms to Use Accelerators with FDPS

Authors: Masaki Iwasawa, Daisuke Namekata, Keigo Nitadori, Kentaro Nomura, Long Wang, Miyuki Tsubouchi, Junichiro Makino

Abstract: In this paper, we describe the algorithms we implemented in FDPS to make efficient use of accelerator hardware such as GPGPUs. We have developed FDPS to make it possible for many researchers to develop their own high-performance parallel particle-based simulation programs without spending large amount of time for parallelization and performance tuning. The basic idea of FDPS is to provide a high-p… ▽ More In this paper, we describe the algorithms we implemented in FDPS to make efficient use of accelerator hardware such as GPGPUs. We have developed FDPS to make it possible for many researchers to develop their own high-performance parallel particle-based simulation programs without spending large amount of time for parallelization and performance tuning. The basic idea of FDPS is to provide a high-performance implementation of parallel algorithms for particle-based simulations in a "generic" form, so that researchers can define their own particle data structure and interparticle interaction functions and supply them to FDPS. FDPS compiled with user-supplied data type and interaction function provides all necessary functions for parallelization, and using those functions researchers can write their programs as though they are writing simple non-parallel program. It has been possible to use accelerators with FDPS, by writing the interaction function that uses the accelerator. However, the efficiency was limited by the latency and bandwidth of communication between the CPU and the accelerator and also by the mismatch between the available degree of parallelism of the interaction function and that of the hardware parallelism. We have modified the interface of user-provided interaction function so that accelerators are more efficiently used. We also implemented new techniques which reduce the amount of work on the side of CPU and amount of communication between CPU and accelerators. We have measured the performance of N-body simulations on a systems with NVIDIA Volta GPGPU using FDPS and the achieved performance is around 27 \% of the theoretical peak limit. We have constructed a detailed performance model, and found that the current implementation can achieve good performance on systems with much smaller memory and communication bandwidth. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Implementation and Performance of Barnes-Hut N-body algorithm on Extreme-scale Heterogeneous Many-core Architectures

Authors: Masaki Iwasawa, Daisuke Namekata, Ryo Sakamoto, Takashi Nakamura, Yasuyuki Kimura, Keigo Nitadori, Long Wang, Miyuki Tsubouchi, Jun Makino, Zhao Liu, Haohuan Fu, Guangwen Yang

Abstract: In this paper, we report the implementation and measured performance of our extreme-scale global simulation code on Sunway TaihuLight and two PEZY-SC2 systems: Shoubu System B and Gyoukou. The numerical algorithm is the parallel Barnes-Hut tree algorithm, which has been used in many large-scale astrophysical particle-based simulations. Our implementation is based on our FDPS framework. However, th… ▽ More In this paper, we report the implementation and measured performance of our extreme-scale global simulation code on Sunway TaihuLight and two PEZY-SC2 systems: Shoubu System B and Gyoukou. The numerical algorithm is the parallel Barnes-Hut tree algorithm, which has been used in many large-scale astrophysical particle-based simulations. Our implementation is based on our FDPS framework. However, the extremely large numbers of cores of the systems used (10M on TaihuLight and 16M on Gyoukou) and their relatively poor memory and network bandwidth pose new challenges. We describe the new algorithms introduced to achieve high efficiency on machines with low memory bandwidth. The measured performance is 47.9, 10.6 PF, and 1.01PF on TaihuLight, Gyoukou and Shoubu System B (efficiency 40\%, 23.5\% and 35.5\%). The current code is developed for the simulation of planetary rings, but most of the new algorithms are useful for other simulations, and are now available in the FDPS framework. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comment on "Individual external dose monitoring of all citizens of Date City by passive dosimeter 5 to 51 months after the Fukushima NPP accident (series): II"

Authors: Shin-ichi Kurokawa, Yoh Tanimoto, Yutaka Hamaoka, Kyo Kageura, Jun Makino, Masaki Oshikawa

Abstract: This Letter to the Editor points out several inconsistencies found in the paper written by Makoto Miyazaki and Ryugo Hayano published as J. Radiol. Prot. 37(2017) 623-634. This Letter to the Editor had been in the stage of "is ready to accept" from November 17, 2018 to March 23, 2020 when it has been accepted for publication. Onto this Letter to the Editor, an additional Letter to the Editor on th… ▽ More This Letter to the Editor points out several inconsistencies found in the paper written by Makoto Miyazaki and Ryugo Hayano published as J. Radiol. Prot. 37(2017) 623-634. This Letter to the Editor had been in the stage of "is ready to accept" from November 17, 2018 to March 23, 2020 when it has been accepted for publication. Onto this Letter to the Editor, an additional Letter to the Editor on the same paper has been added. The added Letter to the Editor has been provisionally accepted on April 7, 2020. △ Less

Submitted 5 August, 2020; v1 submitted 29 December, 2018; originally announced December 2018.

Hermite integrator for high-order mesh-free schemes

Authors: Satoko Yamamoto, Junichiro Makino

Abstract: In most of mesh-free methods, the calculation of interactions between sample points or particles is the most time consuming. When we use mesh-free methods with high spatial orders, the order of the time integration should also be high. If we use usual Runge-Kutta schemes, we need to perform the interaction calculation multiple times per one time step. One way to reduce the number of interaction ca… ▽ More In most of mesh-free methods, the calculation of interactions between sample points or particles is the most time consuming. When we use mesh-free methods with high spatial orders, the order of the time integration should also be high. If we use usual Runge-Kutta schemes, we need to perform the interaction calculation multiple times per one time step. One way to reduce the number of interaction calculations is to use Hermite schemes, which use the time derivatives of the right hand side of differential equations, since Hermite schemes require smaller number of interaction calculations than RK schemes do to achieve the same order. In this paper, we construct a Hermite scheme for a mesh-free method with high spatial orders. We performed several numerical tests with fourth-order Hermite schemes and Runge-Kutta schemes. We found that, for both of Hermite and Runge-Kutta schemes, the overall error is determined by the error of spatial derivatives, for timesteps smaller than the stability limit. The calculation cost at the timestep size of the stability limit is smaller for Hermite schemes. Therefore, we conclude that Hermite schemes are more efficient than Runge-Kutta schemes and thus useful for high-order mesh-free methods for Lagrangian Hydrodynamics. △ Less

Submitted 20 November, 2018; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: 22 pages, 16 figures, Accepted for publication in Pub. Astron. Soc. J

INO: Interplanetary Network of Optical Lattice Clocks

Authors: Toshikazu Ebisuzaki, Hidetoshi Katori, Jun'ichiro Makino, Atsushi Noda, Hisaaki Shinkai, Toru Tamagawa

Abstract: The new technique of measuring frequency by optical lattice clocks now approaches to the relative precision of $(Δf/f)=O(10^{-18})$. We propose to place such precise clocks in space and to use Doppler tracking method for detecting low-frequency gravitational wave below 1 Hz. Our idea is to locate three spacecrafts at one A.U. distance (say at L1, L4 & L5 of the Sun-Earth orbit), and apply the Dopp… ▽ More The new technique of measuring frequency by optical lattice clocks now approaches to the relative precision of $(Δf/f)=O(10^{-18})$. We propose to place such precise clocks in space and to use Doppler tracking method for detecting low-frequency gravitational wave below 1 Hz. Our idea is to locate three spacecrafts at one A.U. distance (say at L1, L4 & L5 of the Sun-Earth orbit), and apply the Doppler tracking method by communicating "the time" each other. Applying the current available technologies, we obtain the sensitivity for gravitational wave with three or four-order improvement ($h_{\rm n}\sim 10^{-17}$ or $10^{-18}$ level in $10^{-5}$Hz -- $1$ Hz) than that of Cassini spacecraft in 2001. This sensitivity enables us to observe black-hole mergers of their mass greater than $10^5 M_\odot$ in the cosmological scale. Based on the hierarchical growth model of black-holes in galaxies, we estimate the event rate of detection will be 20-50 a year. We nickname "INO" (Interplanetary Network of Optical Lattice Clocks) for this system, named after Tadataka Ino (1745--1818), a Japanese astronomer, cartographer, and geodesist. △ Less

Submitted 11 February, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: 11 pages, 4 figures; added references; to be published in Int. J. Modern Physics D

Fortran interface layer of the framework for developing particle simulator FDPS

Authors: Daisuke Namekata, Masaki Iwasawa, Keigo Nitadori, Ataru Tanikawa, Takayuki Muranushi, Long Wang, Natsuki Hosono, Kentaro Nomura, Junichiro Makino

Abstract: Numerical simulations based on particle methods have been widely used in various fields including astrophysics. To date, simulation softwares have been developed by individual researchers or research groups in each field, with a huge amount of time and effort, even though numerical algorithms used are very similar. To improve the situation, we have developed a framework, called FDPS, which enables… ▽ More Numerical simulations based on particle methods have been widely used in various fields including astrophysics. To date, simulation softwares have been developed by individual researchers or research groups in each field, with a huge amount of time and effort, even though numerical algorithms used are very similar. To improve the situation, we have developed a framework, called FDPS, which enables researchers to easily develop massively parallel particle simulation codes for arbitrary particle methods. Until version 3.0, FDPS have provided API only for C++ programing language. This limitation comes from the fact that FDPS is developed using the template feature in C++, which is essential to support arbitrary data types of particle. However, there are many researchers who use Fortran to develop their codes. Thus, the previous versions of FDPS require such people to invest much time to learn C++. This is inefficient. To cope with this problem, we newly developed a Fortran interface layer in FDPS, which provides API for Fortran. In order to support arbitrary data types of particle in Fortran, we design the Fortran interface layer as follows. Based on a given derived data type in Fortran representing particle, a Python script provided by us automatically generates a library that manipulates the C++ core part of FDPS. This library is seen as a Fortran module providing API of FDPS from the Fortran side and uses C programs internally to interoperate Fortran with C++. In this way, we have overcome several technical issues when emulating `template' in Fortran. By using the Fortran interface, users can develop all parts of their codes in Fortran. We show that the overhead of the Fortran interface part is sufficiently small and a code written in Fortran shows a performance practically identical to the one written in C++. △ Less

Submitted 25 April, 2018; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: 10 pages, 10 figures; accepted for publication in PASJ; a typo in author name is corrected

The Detection Rates of Merging Binary Black Holes Originating from Star Clusters and Their Mass Function

Authors: Michiko Fujii, Ataru Tanikawa, Junichiro Makino

Abstract: Advanced LIGO achieved the first detection of the gravitational wave, which was from a merging binary black hole (BBH). In the near future, more merger events will be observed, and the mass distribution of them will become available. The mass distribution of merger events reflects the evolutionary path of BBHs: dynamical formation in dense star clusters or common envelope evolution from primordial… ▽ More Advanced LIGO achieved the first detection of the gravitational wave, which was from a merging binary black hole (BBH). In the near future, more merger events will be observed, and the mass distribution of them will become available. The mass distribution of merger events reflects the evolutionary path of BBHs: dynamical formation in dense star clusters or common envelope evolution from primordial binaries. In this paper, we estimate the detection rate of merging BBHs which dynamically formed in dense star clusters by combining the results of $N$-body simulations, modeling of globular clusters, and cosmic star-cluster formation history. We estimate that the merger rate density in the local universe within the redshift of 0.1 is 13--57 Gpc$^{-3}$ yr$^{-1}$. We find that the detection rate is 0.23--4.6 per year for the current sensitivity limit and that it would increase to 5.1--99 per year for the designed sensitivity which will be achieved in 2019. The distribution of merger rate density in the local universe as a function of redshifted chirp mass has a peak close to the low-mass end.The chirp mass function of the detected mergers, on the other hand, has a peak at the high-mass end, but is almost flat. This difference is simply because the detection range is larger for more massive BBHs. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 18 pages, 16 figures, PASJ accepted

A Formulation of Consistent Particle Hydrodynamics in Strong Form

Authors: Satoko Yamamoto, Junichiro Makino

Abstract: In fluid dynamical simulations in astrophysics, large deformations are common and surface tracking is sometimes necessary. Smoothed Particle Hydrodynamics (SPH) method has been used in many of such simulations. Recently, however, it has been shown that SPH cannot handle contact discontinuities or free surfaces accurately. There are several reasons for this problem. The first one is that SPH requir… ▽ More In fluid dynamical simulations in astrophysics, large deformations are common and surface tracking is sometimes necessary. Smoothed Particle Hydrodynamics (SPH) method has been used in many of such simulations. Recently, however, it has been shown that SPH cannot handle contact discontinuities or free surfaces accurately. There are several reasons for this problem. The first one is that SPH requires that the density is continuous and differentiable. The second one is that SPH does not have the consistency, and thus the accuracy is zeroth order in space. In addition, we cannot express accurate boundary conditions with SPH. In this paper, we propose a novel, high-order scheme for particle-based hydrodynamics of compress- ible fluid. Our method is based on kernel-weighted high-order fitting polynomial for intensive variables. With this approach, we can construct a scheme which solves all of the three prob- lems described above. For shock capturing, we use a tensor form of von-Neumann-Richtmyer artificial viscosity. We have applied our method to many test problems and obtained excel- lent result. Our method is not conservative, since particles do not have mass or energy, but only their densities. However, because of the Lagrangian nature of our scheme, the violation of the conservation laws turned out to be small. We name this method Consistent Particle Hydrodynamics in Strong Form (CPHSF). △ Less

Submitted 19 January, 2017; originally announced January 2017.

Comments: 42 pages, 34figures

Unconvergence of Very Large Scale GI Simulations

Authors: Natsuki Hosono, Masaki Iwasawa, Ataru Tanikawa, Keigo Nitadori, Takayuki Muranushi, Junichiro Makino

Abstract: The giant impact (GI) is one of the most important hypotheses both in planetary science and geoscience, since it is related to the origin of the Moon and also the initial condition of the Earth. A number of numerical simulations have been done using the smoothed particle hydrodynamics (SPH) method. However, GI hypothesis is currently in a crisis. The "canonical" GI scenario failed to explain the i… ▽ More The giant impact (GI) is one of the most important hypotheses both in planetary science and geoscience, since it is related to the origin of the Moon and also the initial condition of the Earth. A number of numerical simulations have been done using the smoothed particle hydrodynamics (SPH) method. However, GI hypothesis is currently in a crisis. The "canonical" GI scenario failed to explain the identical isotope ratio between the Earth and the Moon. On the other hand, little has been known about the reliability of the result of GI simulations. In this paper, we discuss the effect of the resolution on the results of the GI simulations by varying the number of particles from $3 \times10^3$ to $10^8$. We found that the results does not converge, but shows oscillatory behaviour. We discuss the origin of this oscillatory behaviour. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: Accepted to PASJ, an animation is available at https://vimeo.com/194156367

Implementation and evaluation of data-compression algorithms for irregular-grid iterative methods on the PEZY-SC processor

Authors: Naoki Yoshifuji, Ryo Sakamoto, Keigo Nitadori, Jun Makino

Abstract: Iterative methods on irregular grids have been used widely in all areas of comptational science and engineering for solving partial differential equations with complex geometry. They provide the flexibility to express complex shapes with relatively low computational cost. However, the direction of the evolution of high-performance processors in the last two decades have caused serious degradation… ▽ More Iterative methods on irregular grids have been used widely in all areas of comptational science and engineering for solving partial differential equations with complex geometry. They provide the flexibility to express complex shapes with relatively low computational cost. However, the direction of the evolution of high-performance processors in the last two decades have caused serious degradation of the computational efficiency of iterative methods on irregular grids, because of relatively low memory bandwidth. Data compression can in principle reduce the necessary memory memory bandwidth of iterative methods and thus improve the efficiency. We have implemented several data compression algorithms on the PEZY-SC processor, using the matrix generated for the HPCG benchmark as an example. For the SpMV (Sparse Matrix-Vector multiplication) part of the HPCG benchmark, the best implementation without data compression achieved 11.6Gflops/chip, close to the theoretical limit due to the memory bandwidth. Our implementation with data compression has achieved 32.4Gflops. This is of course rather extreme case, since the grid used in HPCG is geometrically regular and thus its compression efficiency is very high. However, in real applications, it is in many cases possible to make a large part of the grid to have regular geometry, in particular when the resolution is high. Note that we do not need to change the structure of the program, except for the addition of the data compression/decompression subroutines. Thus, we believe the data compression will be very useful way to improve the performance of many applications which rely on the use of irregular grids. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: Talk given at IA3 2016 Sixth Workshop on Irregular Applications: Architectures and Algorithms http://hpc.pnl.gov/IA3/IA3/Program.html

Santa Barbara Cluster Comparison Test with DISPH

Authors: Takayuki Saitoh, Junichiro Makino

Abstract: The Santa Barbara cluster comparison project (Frenk et al. Frenk+1999) revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: Mesh codes gave a core with a constant entropy whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galile… ▽ More The Santa Barbara cluster comparison project (Frenk et al. Frenk+1999) revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: Mesh codes gave a core with a constant entropy whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much "protection" to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density independent formulation of SPH (DISPH), which can handle contact discontinuities accurately, to simulations of a cluster of galaxies, and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations, and is comparable to those obtained with quasi-Lagrangian schemes such as "moving mesh" and "mesh free" schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it. △ Less

Submitted 17 April, 2016; originally announced April 2016.

Comments: 21 pages, 25 figures, accepted for publication in ApJ

The Giant Impact Simulations with Density Independent Smoothed Particle Hydrodynamics

Authors: Natsuki Hosono, Takayuki R. Saitoh, Junichiro Makino, Hidenori Genda, Shigeru Ida

Abstract: At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free s… ▽ More At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity. △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: Accepted for publication in Icarus

A comparison of SPH artificial viscosities and their impact on the Keplerian disk

Authors: Natsuki Hosono, Takayuki R. Saitoh, Junichiro Makino

Abstract: Hydrodynamical simulations of rotating disk play important roles in the field of astrophysical and planetary science. Smoothed Particle Hydrodynamics (SPH) has been widely used for such simulations. It, however, has been known that with SPH, a cold and thin Kepler disk breaks up due to the unwanted angular momentum transfer. Two possible reasons have been suggested for this breaking up of the disk… ▽ More Hydrodynamical simulations of rotating disk play important roles in the field of astrophysical and planetary science. Smoothed Particle Hydrodynamics (SPH) has been widely used for such simulations. It, however, has been known that with SPH, a cold and thin Kepler disk breaks up due to the unwanted angular momentum transfer. Two possible reasons have been suggested for this breaking up of the disk; the artificial viscosity (AV) and the numerical error in the evaluation of pressure gradient in SPH. Which one is dominant has been still unclear. In this paper, we investigate the reason for this rapid breaking up of the disk. We implemented most of popular formulations of AV and switches and measured the angular momentum transfer due to both AV and the error of SPH estimate of pressure gradient. We found that the angular momentum transfer due to AV at the inner edge triggers the breaking up of the disk. We also found that the classical von-Neumann-Richtmyer-Landshoff type AV with a high order estimate for $\nabla \cdot \vec{v}$ can maintain the disk for $\sim 100$ orbits even when used with the standard formulation of SPH. △ Less

Submitted 3 March, 2016; v1 submitted 22 January, 2016; originally announced January 2016.

Comments: Accepted for publication in The Astrophysical Journal. Animations are available at http://v1.jmlab.jp/~hosono/1601.05903/

Global High-resolution $N$-body Simulation of Planet Formation I. Planetesimal Driven Migration

Authors: Junko Kominami, Hiroshi Daisaka, Junichiro Makino, Masaki Fujimoto

Abstract: We investigated whether outward Planetesimal Driven Migration (PDM) takes place or not in simulations when the self gravity of planetesimals is included. We performed $N$-body simulations of planetesimal disks with large width (0.7 - 4AU) which ranges over the ice line. The simulations consisted of two stages. The first stage simulations were carried out to see the runaway growth phase using the p… ▽ More We investigated whether outward Planetesimal Driven Migration (PDM) takes place or not in simulations when the self gravity of planetesimals is included. We performed $N$-body simulations of planetesimal disks with large width (0.7 - 4AU) which ranges over the ice line. The simulations consisted of two stages. The first stage simulations were carried out to see the runaway growth phase using the planetesimals of initially the same mass. The runaway growth took place both at the inner edge of the disk and at the region just outside the ice line. This result was utilized for the initial setup of the second stage simulations in which the runaway bodies just outside the ice line were replaced by the protoplanets with about the isolation mass. In the second stage simulations, the outward migration of the protoplanet was followed by the stopping of the migration due to the increase of the random velocity of the planetesimals. Due to this increase of random velocities, one of the PDM criteria derived in Minton and Levison (2014) was broken. In the current simulations, the effect of the gas disk is not considered. It is likely that the gas disk plays an important role in planetesimal driven migration, and we plan to study its effect in future papers. △ Less

Submitted 14 February, 2017; v1 submitted 19 January, 2016; originally announced January 2016.

Comments: Accepted to ApJ

Implementation and performance of FDPS: A Framework Developing Parallel Particle Simulation Codes

Authors: Masaki Iwasawa, Ataru Tanikawa, Natsuki Hosono, Keigo Nitadori, Takayuki Muranushi, Junichiro Makino

Abstract: We present the basic idea, implementation, measured performance and performance model of FDPS (Framework for developing particle simulators). FDPS is an application-development framework which helps the researchers to develop particle-based simulation programs for large-scale distributed-memory parallel supercomputers. A particle-based simulation program for distributed-memory parallel computers n… ▽ More We present the basic idea, implementation, measured performance and performance model of FDPS (Framework for developing particle simulators). FDPS is an application-development framework which helps the researchers to develop particle-based simulation programs for large-scale distributed-memory parallel supercomputers. A particle-based simulation program for distributed-memory parallel computers needs to perform domain decomposition, redistribution of particles, and gathering of particle information for interaction calculation. Also, even if distributed-memory parallel computers are not used, in order to reduce the amount of computation, algorithms such as Barnes-Hut tree method should be used for long-range interactions. For short-range interactions, some methods to limit the calculation to neighbor particles are necessary. FDPS provides all of these necessary functions for efficient parallel execution of particle-based simulations as "templates", which are independent of the actual data structure of particles and the functional form of the interaction. By using FDPS, researchers can write their programs with the amount of work necessary to write a simple, sequential and unoptimized program of O(N^2) calculation cost, and yet the program, once compiled with FDPS, will run efficiently on large-scale parallel supercomputers. A simple gravitational N-body program can be written in around 120 lines. We report the actual performance of these programs and the performance model. The weak scaling performance is very good, and almost linear speedup was obtained for up to the full system of K computer. The minimum calculation time per timestep is in the range of 30 ms (N=10^7) to 300 ms (N=10^9). These are currently limited by the time for the calculation of the domain decomposition and communication necessary for the interaction calculation. We discuss how we can overcome these bottlenecks. △ Less

Submitted 24 April, 2016; v1 submitted 13 January, 2016; originally announced January 2016.

Comments: 22 pages, 27 figures,accepted for publication in PASJ. The FDPS package is here https://github.com/fdps/fdps

GPU-Enabled Particle-Particle Particle-Tree Scheme for Simulating Dense Stellar Cluster System

Authors: Masaki Iwasawa, Simon Portegies Zwart, Junichiro Makino

Abstract: We describe the implementation and performance of the ${\rm P^3T}$ (Particle-Particle Particle-Tree) scheme for simulating dense stellar systems. In ${\rm P^3T}$, the force experienced by a particle is split into short-range and long-range contributions. Short-range forces are evaluated by direct summation and integrated with the fourth order Hermite predictor-corrector method with the block times… ▽ More We describe the implementation and performance of the ${\rm P^3T}$ (Particle-Particle Particle-Tree) scheme for simulating dense stellar systems. In ${\rm P^3T}$, the force experienced by a particle is split into short-range and long-range contributions. Short-range forces are evaluated by direct summation and integrated with the fourth order Hermite predictor-corrector method with the block timesteps. For long-range forces, we use a combination of the Barnes-Hut tree code and the leapfrog integrator. The tree part of our simulation environment is accelerated using graphical processing units (GPU), whereas the direct summation is carried out on the host CPU. Our code gives excellent performance and accuracy for star cluster simulations with a large number of particles even when the core size of the star cluster is small. △ Less

Submitted 15 June, 2015; originally announced June 2015.

Smoothed Particle Hydrodynamics with Smoothed Pseudo-Density

Authors: Satoko Yamamoto, Takayuki R. Saitoh, Junichiro Makino

Abstract: In this paper, we present a new formulation of smoothed particle hydrodynamics (SPH), which, unlike the standard SPH (SSPH), is well-behaved at the contact discontinuity. The SSPH scheme cannot handle discontinuities in density (e.g. the contact discontinuity and the free surface), because it requires that the density of fluid is positive and continuous everywhere. Thus there is inconsistency in t… ▽ More In this paper, we present a new formulation of smoothed particle hydrodynamics (SPH), which, unlike the standard SPH (SSPH), is well-behaved at the contact discontinuity. The SSPH scheme cannot handle discontinuities in density (e.g. the contact discontinuity and the free surface), because it requires that the density of fluid is positive and continuous everywhere. Thus there is inconsistency in the formulation of the SSPH scheme at discontinuities of the fluid density. To solve this problem, we introduce a new quantity associated with particles and "density" of that quantity. This "density" evolves through the usual continuity equation with an additional artificial diffusion term, in order to guarantee the continuity of "density". We use this "density" or pseudo density, instead of the mass density, to formulate our SPH scheme. We call our new method as SPH with smoothed pseudo-density (SPSPH). We show that our new scheme is physically consistent and can handle discontinuities quite well. △ Less

Submitted 24 January, 2015; originally announced January 2015.

Comments: 32 pages, 13 figures

A development of an accelerator board dedicated for multi-precision arithmetic operations and its application to Feynman loop integrals

Authors: Shinji Motoki, Hiroshi Daisaka, Naohito Nakasato, Tadashi Ishikawa, Fukuko Yuasa, Toshiyuki Fukushige, Atsushi Kawai, Junichiro Makino

Abstract: Higher order corrections in perturbative quantum field theory are required for precise theoretical analysis to investigate new physics beyond the Standard Model. This indicates that we need to evaluate Feynman loop diagram with multi-loop integral which may require multi-precision calculation. We developed a dedicated accelerator system for multi-precision calculation (GRAPE9-MPX). We present perf… ▽ More Higher order corrections in perturbative quantum field theory are required for precise theoretical analysis to investigate new physics beyond the Standard Model. This indicates that we need to evaluate Feynman loop diagram with multi-loop integral which may require multi-precision calculation. We developed a dedicated accelerator system for multi-precision calculation (GRAPE9-MPX). We present performance results of our system for the case of Feynman two-loop box and three-loop selfenergy diagrams with multi-precision. △ Less

Submitted 30 November, 2014; v1 submitted 13 October, 2014; originally announced October 2014.

Comments: 6 pages, 8 figures, 1 table; submitted to the proceedings of the 16th International workshop on Advanced Computing and Analysis Techniques in physics research (ACAT 2014), Prague

Binary Formation in Planetesimal Disks II. Planetesimals with Mass Spectrum

Authors: Junko D. Kominami, Junichiro Makino

Abstract: Many massive objects have been found in the outer region of the Solar system. How they were formed and evolved has not been well understood, although there have been intensive studies on accretion process of terrestrial planets. One of the mysteries is the existence of binary planetesimals with near-equal mass components and highly eccentric orbits. These binary planetesimals are quite different f… ▽ More Many massive objects have been found in the outer region of the Solar system. How they were formed and evolved has not been well understood, although there have been intensive studies on accretion process of terrestrial planets. One of the mysteries is the existence of binary planetesimals with near-equal mass components and highly eccentric orbits. These binary planetesimals are quite different from the satellites observed in the asteroid belt region. The ratio of the Hill radius to the physical radius of the planetesimals is much larger for the outer region of the disk, compared to the inner region of the disk. The Hill radius increases with the semi major axis. Therefore, planetesimals in the outer region can form close and eccentric binaries, while those in the inner region would simply collide. In this paper, we carried out $N$-body simulations in different regions of the disk and studied if binaries form in the outer region of the disk. We found that large planetesimals tend to form binaries. A significant fraction of large planetesimals are components of the binaries. Planetesimals that become the components of binaries eventually collide with a third body, through three-body encounters. Thus, the existence of binaries can enhance the growth rate of planetesimals in the Trans-Neptunian Object (TNO) region. △ Less

Submitted 16 September, 2014; originally announced September 2014.

Comments: Accepted, PASJ

Evolution of star clusters in a cosmological tidal field

Authors: Steven Rieder, Tomoaki Ishiyama, Paul Langelaan, Junichiro Makino, Stephen L. W. McMillan, Simon Portegies Zwart

Abstract: 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 initializ… ▽ More 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. △ Less

Submitted 27 September, 2013; originally announced September 2013.

Comments: 14 pages, 14 figures, accepted for publication in MNRAS

Density Independent Smoothed Particle Hydrodynamics for Non-Ideal Equation of State

Authors: Natsuki Hosono, Takayuki R. Saitoh, Junichiro Makino

Abstract: The smoothed particle hydrodynamics (SPH) method is a useful numerical tool for the study of a variety of astrophysical and planetlogical problems. However, it turned out that the standard SPH algorithm has problems in dealing with hydrodynamical instabilities. This problem is due to the assumption that the local density distribution is differentiable. In order to solve this problem, a new SPH for… ▽ More The smoothed particle hydrodynamics (SPH) method is a useful numerical tool for the study of a variety of astrophysical and planetlogical problems. However, it turned out that the standard SPH algorithm has problems in dealing with hydrodynamical instabilities. This problem is due to the assumption that the local density distribution is differentiable. In order to solve this problem, a new SPH formulation, which does not require the differentiability of the density, have been proposed. This new SPH method improved the treatment of hydrodynamical instabilities. This method, however, is applicable only to the equation of state (EOS) of the ideal gas. In this paper, we describe how to extend the new SPH method to non-ideal EOS. We present the results of various standard numerical tests for non-ideal EOS. Our new method works well for non-ideal EOS. We conclude that our new SPH can handle hydrodynamical instabilities for an arbitrary EOS and that it is an attractive alternative to the standard SPH. △ Less

Submitted 3 July, 2013; originally announced July 2013.

Comments: 15 pages, 8 figures

Flaring up of the Compact Cloud G2 during the Close Encounter with Sgr A*

Authors: Takayuki R. Saitoh, Junichiro Makino, Yoshiharu Asaki, Junichi Baba, Shinya Komugi, Makoto Miyoshi, Tohru Nagao, Masaaki Takahashi, Takaaki Takeda, Masato Tsuboi, Ken-ichi Wakamatsu

Abstract: A compact gas cloud G2 is predicted to reach the pericenter of its orbit around the super massive black hole (SMBH) of our galaxy, Sagittarius A* (Sgr A*). This event will give us a rare opportunity to observe the interaction between SMBH and gas around it. We report the result of the fully three-dimensional simulation of the evolution of G2 during the first pericenter passage. The strong tidal fo… ▽ More A compact gas cloud G2 is predicted to reach the pericenter of its orbit around the super massive black hole (SMBH) of our galaxy, Sagittarius A* (Sgr A*). This event will give us a rare opportunity to observe the interaction between SMBH and gas around it. We report the result of the fully three-dimensional simulation of the evolution of G2 during the first pericenter passage. The strong tidal force by the SMBH stretches the cloud along its orbit, and compresses it strongly in the vertical direction, resulting in the heating up and flaring up of the cloud. The bolometric luminosity will reach the maximum of $\sim100 L_{\odot}$. This flare should be easily observed in the near infrared. △ Less

Submitted 9 August, 2013; v1 submitted 3 December, 2012; originally announced December 2012.

Comments: 8 pages, 4 figures. Accepted for publication in PASJ. Animation movie and stereograms are available at http://v1.jmlab.jp/~saitoh/G2/

4.45 Pflops Astrophysical N-Body Simulation on K computer -- The Gravitational Trillion-Body Problem

Authors: Tomoaki Ishiyama, Keigo Nitadori, Junichiro Makino

Abstract: As an entry for the 2012 Gordon-Bell performance prize, we report performance results of astrophysical N-body simulations of one trillion particles performed on the full system of K computer. This is the first gravitational trillion-body simulation in the world. We describe the scientific motivation, the numerical algorithm, the parallelization strategy, and the performance analysis. Unlike many p… ▽ More As an entry for the 2012 Gordon-Bell performance prize, we report performance results of astrophysical N-body simulations of one trillion particles performed on the full system of K computer. This is the first gravitational trillion-body simulation in the world. We describe the scientific motivation, the numerical algorithm, the parallelization strategy, and the performance analysis. Unlike many previous Gordon-Bell prize winners that used the tree algorithm for astrophysical N-body simulations, we used the hybrid TreePM method, for similar level of accuracy in which the short-range force is calculated by the tree algorithm, and the long-range force is solved by the particle-mesh algorithm. We developed a highly-tuned gravity kernel for short-range forces, and a novel communication algorithm for long-range forces. The average performance on 24576 and 82944 nodes of K computer are 1.53 and 4.45 Pflops, which correspond to 49% and 42% of the peak speed. △ Less

Submitted 13 April, 2015; v1 submitted 19 November, 2012; originally announced November 2012.

Comments: 10 pages, 6 figures, Proceedings of Supercomputing 2012 (http://sc12.supercomputing.org/), Gordon Bell Prize Winner. Additional information is http://www.ccs.tsukuba.ac.jp/CCS/eng/gbp2012

Journal ref: SC '12 Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis, Article No. 5 (2012)

Few-Body Modes of Binary Formation in Core Collapse

Authors: Ataru Tanikawa, Douglas C. Heggie, Piet Hut, Junichiro Makino

Abstract: At the moment of deepest core collapse, a star cluster core contains less than ten stars. This small number makes the traditional treatment of hard binary formation, assuming a homogeneous background density, suspect. In a previous paper, we have found that indeed the conventional wisdom of binary formation, based on three-body encounters, is incorrect. Here we refine that insight, by further diss… ▽ More At the moment of deepest core collapse, a star cluster core contains less than ten stars. This small number makes the traditional treatment of hard binary formation, assuming a homogeneous background density, suspect. In a previous paper, we have found that indeed the conventional wisdom of binary formation, based on three-body encounters, is incorrect. Here we refine that insight, by further dissecting the subsequent steps leading to hard binary formation. For this purpose, we add some analysis tools in order to make the study less subjective. We find that the conventional treatment does remain valid for direct three-body scattering, but fails for resonant three-body scattering. Especially democratic resonance scattering, which forms an important part of the analytical theory of three-body binary formation, takes too much space and time to be approximated as being isolated, in the context of a cluster core around core collapse. We conclude that, while three-body encounters can be analytically approximated as isolated, subsequent strong perturbations typically occur whenever those encounters give rise to democratic resonances. We present analytical estimates postdicting our numerical results. If we only had been a bit more clever, we could have predicted this qualitative behaviour. △ Less

Submitted 20 November, 2013; v1 submitted 20 August, 2012; originally announced August 2012.

Comments: 44 pages, 18 figures, 1 table

A Density Independent Formulation of Smoothed Particle Hydrodynamics

Authors: Takayuki R. Saitoh, Junichiro Makino

Abstract: The standard formulation of the smoothed particle hydrodynamics (SPH) assumes that the local density distribution is differentiable. This assumption is used to derive the spatial derivatives of other quantities. However, this assumption breaks down at the contact discontinuity. At the contact discontinuity, the density of the low-density side is overestimated while that of the high-density side is… ▽ More The standard formulation of the smoothed particle hydrodynamics (SPH) assumes that the local density distribution is differentiable. This assumption is used to derive the spatial derivatives of other quantities. However, this assumption breaks down at the contact discontinuity. At the contact discontinuity, the density of the low-density side is overestimated while that of the high-density side is underestimated. As a result, the pressure of the low (high) density side is over (under) estimated. Thus, unphysical repulsive force appears at the contact discontinuity, resulting in the effective surface tension. This tension suppresses fluid instabilities. In this paper, we present a new formulation of SPH, which does not require the differentiability of density. Instead of the mass density, we adopt the internal energy density (pressure), and its arbitrary function, which are smoothed quantities at the contact discontinuity, as the volume element used for the kernel integration. We call this new formulation density independent SPH (DISPH). It handles the contact discontinuity without numerical problems. The results of standard tests such as the shock tube, Kelvin-Helmholtz and Rayleigh-Taylor instabilities, point like explosion, and blob tests are all very favorable to DISPH. We conclude that DISPH solved most of known difficulties of the standard SPH, without introducing additional numerical diffusion or breaking the exact force symmetry or energy conservation. Our new SPH includes the formulation proposed by Ritchie & Thomas (2001) as a special case. Our formulation can be extended to handle a non-ideal gas easily. △ Less

Submitted 26 February, 2013; v1 submitted 20 February, 2012; originally announced February 2012.

Comments: 24 pages, 21 figures. Movies and high resolution figures are available at http://v1.jmlab.jp/~saitoh/sph/index.html

PSDF: Particle Stream Data Format for N-Body Simulations

Authors: Will M. Farr, Jeff Ames, Piet Hut, Junichiro Makino, Steve McMillan, Takayuki Muranushi, Koichi Nakamura, Keigo Nitadori, Simon Portegies Zwart

Abstract: We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same… ▽ More We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same approach as XML but with a simpler syntax. Together with a specification of PSDF, we provide background and motivation, as well as specific examples in a variety of computer languages. We also offer a web site from which these examples can be retrieved, in order to make it easy to augment existing codes in order to give them the option to produce PSDF output. △ Less

Submitted 9 January, 2012; originally announced January 2012.

Comments: 5 pages; submitted to New Astronomy

Journal ref: New Astronomy 17 (2012), pp. 520--523

Origin of multiple nuclei in ultraluminous infrared galaxies

Authors: Hidenori Matsui, Takayuki R. Saitoh, Junichiro Makino, Keiichi Wada, Kohji Tomisaka, Eiichiro Kokubo, Hiroshi Daisaka, Takashi Okamoto, Naoki Yoshida

Abstract: Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies.… ▽ More Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 12 pages, 9 figures, accepted for publication in ApJ

Unexpected Formation Modes of the First Hard Binary in Core Collapse

Authors: Ataru Tanikawa, Piet Hut, Junichiro Makino

Abstract: The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that both assumptions are often incorrect: 1) quite a few three-body interactions… ▽ More The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that both assumptions are often incorrect: 1) quite a few three-body interactions produce a hard binary from scratch; 2) and in many cases there are more than three bodies directly and simultaneously involved in the production of the first binary. The main reason for the discrepancies is that the core of a star cluster, at the first deep collapse, contains typically only five or so stars. Therefore, the homogeneous background assumption, which still would be reasonable for, say, 25 stars, utterly breaks down. There have been some speculations in this direction, but we demonstrate this result here explicitly, for the first time. △ Less

Submitted 19 July, 2011; originally announced July 2011.

Comments: 15 pages, 11 figures, 1 table, submitted to New Astronomy

Binary Formation in Planetesimal Disks I. Equal Mass Planetesimals

Authors: Junko K. Daisaka, Junichiro Makino, Hiroshi Daisaka

Abstract: As to April 2010, 48 TNO (trans-Neptunian Object) binaries have been found. This is about 6% of known TNOs. However, in previous theoretical studies of planetary formation in the TNO region, the effect of binary formation has been neglected. TNO binaries can be formed through a variety of mechanisms, such as three-body process, dynamical friction on two massive bodies, inelastic collisions between… ▽ More As to April 2010, 48 TNO (trans-Neptunian Object) binaries have been found. This is about 6% of known TNOs. However, in previous theoretical studies of planetary formation in the TNO region, the effect of binary formation has been neglected. TNO binaries can be formed through a variety of mechanisms, such as three-body process, dynamical friction on two massive bodies, inelastic collisions between two bodies etc. Most of these mechanisms become more effective as the distance from the Sun increases. In this paper, we studied three-body process using direct N-body simulations. We systematically changed the distance from the Sun, the number density of planetesimals, and the radius of the planetesimals and studied the effect of the binaries on the collision rate of planetesimals. In the TNO region, binaries are involved in 1/3 - 1/2 of collisions, and the collision rate is increased by about a factor of a few compared to the theoretical estimate for the direct two-body collisions. Thus, it is possible that the binaries formed through three-body process significantly enhance the collision rate and reduce the growth time scale. In the terrestrial planet region, binaries are less important, because the ratio between the Hill radius and physical size of the planetesimals is relatively small. Although the time scale of our simulations is short, they clearly demonstrated that the accretion process in the TNO region is quite different from that in the terrestrial planet region. Simulations which cover longer time scale are required to obtain more accurate estimate for the accretion enhancement. △ Less

Submitted 3 February, 2011; originally announced February 2011.

Particle-Particle Particle-Tree: A Direct-Tree Hybrid Scheme for Collisional N-Body Simulations

Authors: Shoichi Oshino, Yoko Funato, Junichiro Makino

Abstract: In this paper, we present a new hybrid algorithm for the time integration of collisional N-body systems. In this algorithm, gravitational force between two particles is divided into short-range and long-range terms, using a distance-dependent cutoff function. The long-range interaction is calculated using the tree algorithm and integrated with the constant-timestep leapfrog integrator. The short-r… ▽ More In this paper, we present a new hybrid algorithm for the time integration of collisional N-body systems. In this algorithm, gravitational force between two particles is divided into short-range and long-range terms, using a distance-dependent cutoff function. The long-range interaction is calculated using the tree algorithm and integrated with the constant-timestep leapfrog integrator. The short-range term is calculated directly and integrated with the high-order Hermite scheme. We can reduce the calculation cost per orbital period from O(N^2) to O(N log N), without significantly increasing the long-term integration error. The results of our test simulations show that close encounters are integrated accurately. Long-term errors of the total energy shows random-walk behaviour, because it is dominated by the error caused by tree approximation. △ Less

Submitted 28 January, 2011; originally announced January 2011.

Comments: 22 pages, 15 figures

Shock-induced star cluster formation in colliding galaxies

Authors: Takayuki R. Saitoh, Hiroshi Daisaka, Eiichiro Kokubo, Junichiro Makino, Takashi Okamoto, Kohji Tomisaka, Keiichi Wada, Naoki Yoshida

Abstract: We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found betwe… ▽ More We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters. △ Less

Submitted 19 January, 2011; originally announced January 2011.

Comments: 4 pages, 3 figures, Proceedings of IAU Symposium 270, Computational Star Formation

The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos

Authors: Tomoaki Ishiyama, Steven Rieder, Junichiro Makino, Simon Portegies Zwart, Derek Groen, Keigo Nitadori, Cees de Laat, Stephen McMillan, Kei Hiraki, Stefan Harfst

Abstract: We present the results of the "Cosmogrid" 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… ▽ More We present the results of the "Cosmogrid" 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. △ Less

Submitted 8 April, 2013; v1 submitted 10 January, 2011; originally announced January 2011.

Comments: 15 pages, 18 figures, accepted by ApJ

Journal ref: 2013, ApJ, 767, 146

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

Authors: Derek Groen, Simon Portegies Zwart, Tomoaki Ishiyama, Junichiro Makino

Abstract: 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 seco… ▽ More 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. △ Less

Submitted 3 January, 2011; originally announced January 2011.

Comments: 30 pages, 11 figures, accepted by Comp. Science and Discovery

MSC Class: 68M14 (primary); 68M20; 85-08; 85A40 (secondary) ACM Class: C.2.4; C.2.5

Journal ref: Comput. Sci. Disc. 4 (2011) 015001

Eccentric evolution of SMBH binaries

Authors: Masaki Iwasawa, Sangyong An, Tatsushi Matsubayashi, Yoko Funato, Junichiro Makino

Abstract: In recent numerical simulations \citep{matsubayashi07,lockmann08}, it has been found that the eccentricity of supermassive black hole(SMBH) - intermediate black hole(IMBH) binaries grows toward unity through interactions with stellar background. This increase of eccentricity reduces the merging timescale of the binary through the gravitational radiation to the value well below the Hubble Time. It… ▽ More In recent numerical simulations \citep{matsubayashi07,lockmann08}, it has been found that the eccentricity of supermassive black hole(SMBH) - intermediate black hole(IMBH) binaries grows toward unity through interactions with stellar background. This increase of eccentricity reduces the merging timescale of the binary through the gravitational radiation to the value well below the Hubble Time. It also gives the theoretical explanation of the existence of eccentric binary such as that in OJ287 \citep{lehto96, valtonen08}. In self-consistent N-body simulations, this increase of eccentricity is always observed. On the other hand, the result of scattering experiment between SMBH binaries and field stars \citep{quinlan96} indicated no increase of eccentricity. This discrepancy leaves the high eccentricity of the SMBH binaries in $N$-body simulations unexplained. Here we present a stellar-dynamical mechanism that drives the increase of the eccentricity of an SMBH binary with large mass ratio. There are two key processes involved. The first one is the Kozai mechanism under non-axisymmetric potential, which effectively randomizes the angular momenta of surrounding stars. The other is the selective ejection of stars with prograde orbits. Through these two mechanisms, field stars extract the orbital angular momentum of the SMBH binary. Our proposed mechanism causes the increase in the eccentricity of most of SMBH binaries, resulting in the rapid merger through gravitational wave radiation. Our result has given a definite solution to the "last-parsec problem". △ Less

Submitted 17 November, 2010; originally announced November 2010.

Comments: 10 pages, 5 figures

Gamma-ray Signal from Earth-mass Dark Matter Microhalos

Authors: Tomoaki Ishiyama, Junichiro Makino, Toshikazu Ebisuzaki

Abstract: 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… ▽ More 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'), 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). △ Less

Submitted 8 April, 2013; v1 submitted 17 June, 2010; originally announced June 2010.

Comments: 15 pages, 5 figures, accepted by ApJL

Journal ref: 2010, ApJ, 723, L195

The dynamics of spiral arms in pure stellar disks

Authors: M. S. Fujii, J. Baba, T. R. Saitoh, J. Makino, E. Kokubo, K. Wada

Abstract: It has been believed that spirals in pure stellar disks, especially the ones spontaneously formed, decay in several galactic rotations due to the increase of stellar velocity dispersions. Therefore, some cooling mechanism, for example dissipational effects of the interstellar medium, was assumed to be necessary to keep the spiral arms. Here we show that stellar disks can maintain spiral features f… ▽ More It has been believed that spirals in pure stellar disks, especially the ones spontaneously formed, decay in several galactic rotations due to the increase of stellar velocity dispersions. Therefore, some cooling mechanism, for example dissipational effects of the interstellar medium, was assumed to be necessary to keep the spiral arms. Here we show that stellar disks can maintain spiral features for several tens of rotations without the help of cooling, using a series of high-resolution three-dimensional $N$-body simulations of pure stellar disks. We found that if the number of particles is sufficiently large, e.g., $3\times 10^6$, multi-arm spirals developed in an isolated disk can survive for more than 10 Gyrs. We confirmed that there is a self-regulating mechanism that maintains the amplitude of the spiral arms. Spiral arms increase Toomre's $Q$ of the disk, and the heating rate correlates with the squared amplitude of the spirals. Since the amplitude itself is limited by the value of $Q$, this makes the dynamical heating less effective in the later phase of evolution. A simple analytical argument suggests that the heating is caused by gravitational scattering of stars by spiral arms, and that the self-regulating mechanism in pure-stellar disks can effectively maintain spiral arms on a cosmological timescale. In the case of a smaller number of particles, e.g., $3\times 10^5$, spiral arms grow faster in the beginning of the simulation (while $Q$ is small) and they cause a rapid increase of $Q$. As a result, the spiral arms become faint in several Gyrs. △ Less

Submitted 5 March, 2011; v1 submitted 7 June, 2010; originally announced June 2010.

Comments: 18 pages, 19 figures, accepted for ApJ